2
MPQC: Massively Parallel Quantum Chemistry
5
Machine: i686-pc-linux-gnu
7
Start Time: Sun Jan 9 18:48:10 2005
9
Using ProcMessageGrp for message passing (number of nodes = 1).
10
Using PthreadThreadGrp for threading (number of threads = 1).
11
Using ProcMemoryGrp for distributed shared memory.
12
Total number of processors = 1
14
Using IntegralV3 by default for molecular integrals evaluation
16
Reading file /home/cljanss/src/SC/lib/atominfo.kv.
18
IntCoorGen: generated 9 coordinates.
19
Forming optimization coordinates:
20
SymmMolecularCoor::form_variable_coordinates()
21
expected 6 coordinates
22
found 4 variable coordinates
23
found 0 constant coordinates
24
Reading file /home/cljanss/src/SC/lib/basis/4-31g.kv.
25
Reading file /home/cljanss/src/SC/lib/basis/sto-3g.kv.
27
CLSCF::init: total charge = 0
29
Starting from core Hamiltonian guess
31
Using symmetric orthogonalization.
33
Maximum orthogonalization residual = 1.97637
34
Minimum orthogonalization residual = 0.273929
38
CLSCF::init: total charge = 0
40
Projecting guess wavefunction into the present basis set
42
SCF::compute: energy accuracy = 1.0000000e-06
44
integral intermediate storage = 26045 bytes
45
integral cache = 31972707 bytes
46
nuclear repulsion energy = 18.1371373021
49
iter 1 energy = -338.3388187808 delta = 6.57476e-01
51
iter 2 energy = -338.6241201908 delta = 1.66433e-01
53
iter 3 energy = -338.6296004108 delta = 2.56912e-02
55
iter 4 energy = -338.6301007379 delta = 1.05465e-02
57
iter 5 energy = -338.6301095294 delta = 1.34679e-03
59
iter 6 energy = -338.6301096873 delta = 1.87478e-04
61
iter 7 energy = -338.6301097181 delta = 3.97256e-06
63
HOMO is 7 A' = -0.273200
64
LUMO is 3 A" = 0.524454
66
total scf energy = -338.6301097181
68
Projecting the guess density.
70
The number of electrons in the guess density = 18
71
Using symmetric orthogonalization.
73
Maximum orthogonalization residual = 3.44594
74
Minimum orthogonalization residual = 0.0441062
75
The number of electrons in the projected density = 17.9606
83
matrixkit = <ReplSCMatrixKit>
84
filename = basis2_ph3scf431gcs
85
restart_file = basis2_ph3scf431gcs.ckpt
97
SCF::compute: energy accuracy = 1.0000000e-08
99
integral intermediate storage = 40223 bytes
100
integral cache = 31956737 bytes
101
nuclear repulsion energy = 18.1371373021
104
iter 1 energy = -341.8443756254 delta = 3.80839e-01
106
iter 2 energy = -342.0115679291 delta = 6.00420e-02
108
iter 3 energy = -342.0164966415 delta = 9.24441e-03
110
iter 4 energy = -342.0167609387 delta = 2.80305e-03
112
iter 5 energy = -342.0167761725 delta = 8.19284e-04
114
iter 6 energy = -342.0167762885 delta = 8.57380e-05
116
iter 7 energy = -342.0167762959 delta = 1.57384e-05
118
iter 8 energy = -342.0167762917 delta = 5.60526e-06
120
iter 9 energy = -342.0167762923 delta = 1.05327e-06
122
iter 10 energy = -342.0167762917 delta = 1.46598e-07
124
HOMO is 7 A' = -0.363475
125
LUMO is 3 A" = 0.216161
127
total scf energy = -342.0167762917
129
SCF::compute: gradient accuracy = 1.0000000e-06
132
1 P 0.0009131464 -0.0502167284 0.0000000000
133
2 H 0.0084949756 0.0165349499 -0.0149116581
134
3 H 0.0084949756 0.0165349499 0.0149116581
135
4 H -0.0179030976 0.0171468287 -0.0000000000
136
Value of the MolecularEnergy: -342.0167762917
139
Gradient of the MolecularEnergy:
146
value_accuracy = 9.504686e-09 (1.000000e-08) (computed)
147
gradient_accuracy = 9.504686e-07 (1.000000e-06) (computed)
148
hessian_accuracy = 0.000000e+00 (1.000000e-04)
150
Molecular Coordinates:
151
IntMolecularCoor Parameters:
153
scale_bonds = 1.0000000000
154
scale_bends = 1.0000000000
155
scale_tors = 1.0000000000
156
scale_outs = 1.0000000000
157
symmetry_tolerance = 1.000000e-05
158
simple_tolerance = 1.000000e-03
159
coordinate_tolerance = 1.000000e-07
160
have_fixed_values = 0
161
max_update_steps = 100
162
max_update_disp = 0.500000
163
have_fixed_values = 0
165
Molecular formula: H3P
166
molecule<Molecule>: (
169
{ n atoms geometry }={
170
1 P [ -0.0030062008 0.4698128553 0.0000000000]
171
2 H [ -0.6149106543 -0.1558454669 1.0546274364]
172
3 H [ -0.6149106543 -0.1558454669 -1.0546274364]
173
4 H [ 1.2128275196 -0.1581219416 0.0000000000]
177
30.97376 1.00783 1.00783 1.00783
180
STRE s1 1.37044 1 2 P-H
181
STRE s2 1.37044 1 3 P-H
182
STRE s3 1.36841 1 4 P-H
184
BEND b1 100.62737 2 1 3 H-P-H
185
BEND b2 100.79065 2 1 4 H-P-H
186
BEND b3 100.79065 3 1 4 H-P-H
188
OUT o1 73.05249 2 1 3 4 H-P-H-H
189
OUT o2 -73.05249 3 1 2 4 H-P-H-H
190
OUT o3 72.95148 4 1 2 3 H-P-H-H
192
SymmMolecularCoor Parameters:
193
change_coordinates = no
194
transform_hessian = yes
195
max_kappa2 = 10.000000
202
Natural Population Analysis:
203
n atom charge ne(S) ne(P)
204
1 P -0.004899 5.485784 9.519116
205
2 H 0.001628 0.998372
206
3 H 0.001628 0.998372
207
4 H 0.001644 0.998356
211
density_reset_frequency = 10
212
level_shift = 0.000000
215
charge = 0.0000000000
219
The following keywords in "basis2_ph3scf431gcs.in" were ignored:
220
mpqc:mole:guess_wavefunction:multiplicity
221
mpqc:mole:multiplicity
227
compute gradient: 0.07 0.07
229
one electron gradient: 0.01 0.01
230
overlap gradient: 0.00 0.00
231
two electron gradient: 0.06 0.06
232
contribution: 0.04 0.04
233
start thread: 0.04 0.04
234
stop thread: 0.00 0.00
243
start thread: 0.07 0.04
244
stop thread: 0.00 0.00
246
local data: 0.00 0.00
258
start thread: 0.01 0.01
259
stop thread: 0.00 0.00
261
local data: 0.00 0.00
266
End Time: Sun Jan 9 18:48:10 2005