2
MPQC: Massively Parallel Quantum Chemistry
5
Machine: i686-pc-linux-gnu
7
Start Time: Sun Jan 9 18:48:43 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/pc-3.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 = 8.00023
74
Minimum orthogonalization residual = 0.000279986
75
The number of electrons in the projected density = 17.978
83
matrixkit = <ReplSCMatrixKit>
84
filename = basis2_ph3scfpc3cs
85
restart_file = basis2_ph3scfpc3cs.ckpt
97
SCF::compute: energy accuracy = 1.0000000e-08
99
integral intermediate storage = 3823685 bytes
100
integral cache = 27954539 bytes
101
nuclear repulsion energy = 18.1371373021
104
iter 1 energy = -342.2381446563 delta = 4.17229e-02
106
iter 2 energy = -342.4745806855 delta = 1.03069e-02
108
iter 3 energy = -342.4845846432 delta = 1.73537e-03
110
iter 4 energy = -342.4854988965 delta = 5.02879e-04
112
iter 5 energy = -342.4856719395 delta = 2.76468e-04
114
iter 6 energy = -342.4856765649 delta = 4.02136e-05
116
iter 7 energy = -342.4856769323 delta = 1.22284e-05
118
iter 8 energy = -342.4856769526 delta = 3.62396e-06
120
iter 9 energy = -342.4856769531 delta = 2.57704e-07
122
iter 10 energy = -342.4856769532 delta = 1.38131e-07
124
iter 11 energy = -342.4856769532 delta = 4.21314e-08
126
iter 12 energy = -342.4856769532 delta = 1.41510e-08
128
HOMO is 7 A' = -0.366276
129
LUMO is 8 A' = 0.071631
131
total scf energy = -342.4856769532
133
SCF::compute: gradient accuracy = 1.0000000e-06
136
1 P 0.0010282750 -0.0395754837 0.0000000000
137
2 H 0.0048086348 0.0129804217 -0.0085423536
138
3 H 0.0048086348 0.0129804217 0.0085423536
139
4 H -0.0106455446 0.0136146404 -0.0000000000
140
Value of the MolecularEnergy: -342.4856769532
143
Gradient of the MolecularEnergy:
150
value_accuracy = 3.038721e-09 (1.000000e-08) (computed)
151
gradient_accuracy = 3.038721e-07 (1.000000e-06) (computed)
152
hessian_accuracy = 0.000000e+00 (1.000000e-04)
154
Molecular Coordinates:
155
IntMolecularCoor Parameters:
157
scale_bonds = 1.0000000000
158
scale_bends = 1.0000000000
159
scale_tors = 1.0000000000
160
scale_outs = 1.0000000000
161
symmetry_tolerance = 1.000000e-05
162
simple_tolerance = 1.000000e-03
163
coordinate_tolerance = 1.000000e-07
164
have_fixed_values = 0
165
max_update_steps = 100
166
max_update_disp = 0.500000
167
have_fixed_values = 0
169
Molecular formula: H3P
170
molecule<Molecule>: (
173
{ n atoms geometry }={
174
1 P [ -0.0030062008 0.4698128553 0.0000000000]
175
2 H [ -0.6149106543 -0.1558454669 1.0546274364]
176
3 H [ -0.6149106543 -0.1558454669 -1.0546274364]
177
4 H [ 1.2128275196 -0.1581219416 0.0000000000]
181
30.97376 1.00783 1.00783 1.00783
184
STRE s1 1.37044 1 2 P-H
185
STRE s2 1.37044 1 3 P-H
186
STRE s3 1.36841 1 4 P-H
188
BEND b1 100.62737 2 1 3 H-P-H
189
BEND b2 100.79065 2 1 4 H-P-H
190
BEND b3 100.79065 3 1 4 H-P-H
192
OUT o1 73.05249 2 1 3 4 H-P-H-H
193
OUT o2 -73.05249 3 1 2 4 H-P-H-H
194
OUT o3 72.95148 4 1 2 3 H-P-H-H
196
SymmMolecularCoor Parameters:
197
change_coordinates = no
198
transform_hessian = yes
199
max_kappa2 = 10.000000
206
Natural Population Analysis:
207
n atom charge ne(S) ne(P) ne(D) ne(F) ne(G)
208
1 P 0.197339 5.485716 9.283606 0.031697 0.000997 0.000645
209
2 H -0.065668 1.057752 0.007352 0.000559 0.000005
210
3 H -0.065668 1.057752 0.007352 0.000559 0.000005
211
4 H -0.066002 1.058021 0.007406 0.000570 0.000005
215
density_reset_frequency = 10
216
level_shift = 0.000000
219
charge = 0.0000000000
223
The following keywords in "basis2_ph3scfpc3cs.in" were ignored:
224
mpqc:mole:guess_wavefunction:multiplicity
225
mpqc:mole:multiplicity
231
compute gradient: 85.39 85.39
233
one electron gradient: 0.78 0.78
234
overlap gradient: 0.21 0.21
235
two electron gradient: 84.40 84.40
236
contribution: 82.34 82.34
237
start thread: 82.34 82.33
238
stop thread: 0.00 0.00
240
vector: 217.50 217.49
246
ao_gmat: 216.03 216.02
247
start thread: 216.03 216.02
248
stop thread: 0.00 0.00
250
local data: 0.10 0.10
262
start thread: 0.01 0.01
263
stop thread: 0.00 0.00
265
local data: 0.00 0.00
270
End Time: Sun Jan 9 18:53:47 2005