1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
|
#
# orbital and geometry optimizations of a carbon monoxide
# molecule by the LDA and the cluster method
#
#
# File Name
#
System.CurrrentDirectory ./ # default=./
System.Name CO
level.of.stdout 1 # default=1 (1-3)
level.of.fileout 1 # default=1 (0-2)
#
# Definition of Atomic Species
#
Species.Number 2
<Definition.of.Atomic.Species
O O6.0-s2>1p2>1 O_PBE13
C C6.0-s2>1p2>1 C_PBE13
Definition.of.Atomic.Species>
#
# Atoms
#
Atoms.Number 2
Atoms.SpeciesAndCoordinates.Unit Ang # Ang|AU
<Atoms.SpeciesAndCoordinates
1 O 0.000 0.000 0.000 3.0 3.0
2 C 1.600 0.000 0.000 2.0 2.0
Atoms.SpeciesAndCoordinates>
Atoms.UnitVectors.Unit Ang # Ang|AU
#<Atoms.UnitVectors
# 10.0 0.0 0.0
# 0.0 10.0 0.0
# 0.0 0.0 10.0
#Atoms.UnitVectors>
#
# SCF or Electronic System
#
scf.XcType GGA-PBE # LDA|LSDA-CA|LSDA-PW|GGA-PBE
scf.SpinPolarization off # On|Off|NC
scf.ElectronicTemperature 500.0 # default=300 (K)
scf.energycutoff 150.0 # default=150 (Ry)
scf.maxIter 50 # default=40
scf.EigenvalueSolver cluster # DC|GDC|Cluster|Band
scf.Kgrid 1 1 1 # means n1 x n2 x n3
scf.Mixing.Type rmm-diis # Simple|Rmm-Diis|Gr-Pulay|Kerker|Rmm-Diisk
scf.Init.Mixing.Weight 0.010 # default=0.30
scf.Min.Mixing.Weight 0.001 # default=0.001
scf.Max.Mixing.Weight 0.300 # default=0.40
scf.Mixing.History 6 # default=5
scf.Mixing.StartPulay 6 # default=6
scf.criterion 1.0e-9 # default=1.0e-6 (Hartree)
#
# Orbital Optimization
#
orbitalOpt.Method atoms # Off|Unrestricted|Restricted
orbitalOpt.InitCoes Symmetrical # Symmetrical|Free
orbitalOpt.initPrefactor 0.1 # default=0.1
orbitalOpt.scf.maxIter 40 # default=12
orbitalOpt.Opt.maxIter 5 # default=5
orbitalOpt.per.MDIter 2 # default=1000000
orbitalOpt.criterion 1.0e-4 # default=1.0e-4 (Hartree/borh)
#
# output of contracted orbitals
#
CntOrb.fileout off # on|off, default=off
Num.CntOrb.Atoms 1 # default=1
<Atoms.Cont.Orbitals
1
Atoms.Cont.Orbitals>
#
# MD or Geometry Optimization
#
MD.Type opt # Nomd|Opt|NVE|NVT_VS|NVT_NH
MD.maxIter 3 # default=1
MD.TimeStep 1.0 # default=0.5 (fs)
MD.Opt.criterion 1.0e-4 # default=1.0e-4 (Hartree/bohr)
|
