2
MPQC: Massively Parallel Quantum Chemistry
5
Machine: i686-pc-linux-gnu
7
Start Time: Sun Jan 9 18:38:42 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-4.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 = 9.79516
74
Minimum orthogonalization residual = 4.51364e-05
75
The number of electrons in the projected density = 17.979
83
matrixkit = <ReplSCMatrixKit>
84
filename = basis2_ph3scfpc4cs
85
restart_file = basis2_ph3scfpc4cs.ckpt
97
SCF::compute: energy accuracy = 1.0000000e-08
99
integral intermediate storage = 13451275 bytes
100
integral cache = 17854885 bytes
101
nuclear repulsion energy = 18.1371373021
104
iter 1 energy = -342.2198836133 delta = 2.70830e-02
106
iter 2 energy = -342.4750476404 delta = 8.04035e-03
108
iter 3 energy = -342.4851334287 delta = 9.76994e-04
110
iter 4 energy = -342.4859827659 delta = 2.43201e-04
112
iter 5 energy = -342.4862192867 delta = 1.71273e-04
114
iter 6 energy = -342.4862323834 delta = 4.06361e-05
116
iter 7 energy = -342.4862327690 delta = 6.47574e-06
118
iter 8 energy = -342.4862327953 delta = 2.22426e-06
120
iter 9 energy = -342.4862327961 delta = 2.77310e-07
122
iter 10 energy = -342.4862327963 delta = 1.08951e-07
124
iter 11 energy = -342.4862327963 delta = 3.93533e-08
126
iter 12 energy = -342.4862327963 delta = 1.03108e-08
128
HOMO is 7 A' = -0.366357
129
LUMO is 8 A' = 0.055070
131
total scf energy = -342.4862327963
133
SCF::compute: gradient accuracy = 1.0000000e-06
136
1 P 0.0010271654 -0.0392671030 -0.0000000000
137
2 H 0.0047328489 0.0128781374 -0.0084118631
138
3 H 0.0047328489 0.0128781374 0.0084118631
139
4 H -0.0104928631 0.0135108282 -0.0000000000
140
Value of the MolecularEnergy: -342.4862327963
143
Gradient of the MolecularEnergy:
150
value_accuracy = 3.995639e-09 (1.000000e-08) (computed)
151
gradient_accuracy = 3.995639e-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) ne(H)
208
1 P 0.197727 5.487073 9.281834 0.031466 0.000660 0.000911 0.000329
209
2 H -0.065789 1.055820 0.008046 0.001485 0.000435 0.000003
210
3 H -0.065789 1.055820 0.008046 0.001485 0.000435 0.000003
211
4 H -0.066149 1.056093 0.008109 0.001505 0.000439 0.000003
215
density_reset_frequency = 10
216
level_shift = 0.000000
219
charge = 0.0000000000
223
The following keywords in "basis2_ph3scfpc4cs.in" were ignored:
224
mpqc:mole:guess_wavefunction:multiplicity
225
mpqc:mole:multiplicity
228
mpqc: 3078.58 3078.55
230
calc: 3076.04 3076.01
231
compute gradient: 820.41 820.41
233
one electron gradient: 3.66 3.66
234
overlap gradient: 0.79 0.80
235
two electron gradient: 815.96 815.95
236
contribution: 806.60 806.59
237
start thread: 806.56 806.56
238
stop thread: 0.00 0.00
240
vector: 2255.63 2255.60
244
fock: 2252.50 2252.48
246
ao_gmat: 2249.83 2249.81
247
start thread: 2249.83 2249.81
248
stop thread: 0.00 0.00
250
local data: 0.32 0.31
262
start thread: 0.02 0.01
263
stop thread: 0.00 0.00
265
local data: 0.00 0.00
270
End Time: Sun Jan 9 19:30:01 2005