1
******************************************************************************
2
tstart called on tool.chemistry.gatech.edu
3
Thu Jul 3 18:14:34 2003
10
LABEL = 6-31G** SCF H2O
15
-Geometry before Center-of-Mass shift (a.u.):
17
------------ ----------------- ----------------- -----------------
18
OXYGEN 0.000000000000 0.000000000000 0.000000000000
19
HYDROGEN 0.000000000000 0.000000000000 1.889725988579
20
HYDROGEN 1.840547697619 0.000000000000 -0.428308866009
23
-Rotational constants (cm-1) :
24
A = 24.35463 B = 13.63610 C = 8.74166
25
It is an asymmetric top.
27
-Geometry after Center-of-Mass shift and reorientation (a.u.):
29
------------ ----------------- ----------------- -----------------
30
OXYGEN 0.000000000000 -0.131510165098 0.000000000000
31
HYDROGEN 1.479940996653 1.043580877476 0.000000000000
32
HYDROGEN -1.479940996653 1.043580877476 0.000000000000
35
-SYMMETRY INFORMATION:
36
Computational point group is C2v
37
Number of irr. rep. = 4
39
Number of unique atoms = 2
44
-Basis set on unique center 1:
45
( (S ( 5484.67166000 0.00183107)
46
( 825.23494600 0.01395017)
47
( 188.04695800 0.06844508)
48
( 52.96450000 0.23271434)
49
( 16.89757040 0.47019290)
50
( 5.79963534 0.35852085) )
51
(S ( 15.53961625 -0.11077755)
52
( 3.59993359 -0.14802626)
53
( 1.01376175 1.13076701) )
54
(S ( 0.27000582 1.00000000) )
55
(P ( 15.53961625 0.07087427)
56
( 3.59993359 0.33975284)
57
( 1.01376175 0.72715858) )
58
(P ( 0.27000582 1.00000000) )
59
(D ( 0.80000000 1.00000000) )
62
-Basis set on unique center 2:
63
( (S ( 18.73113696 0.03349460)
64
( 2.82539437 0.23472695)
65
( 0.64012169 0.81375733) )
66
(S ( 0.16127776 1.00000000) )
67
(P ( 1.10000000 1.00000000) )
71
-BASIS SET INFORMATION:
72
Total number of shells = 12
73
Number of primitives = 20
85
-Unique atoms in the canonical coordinate system (a.u.):
87
------------ ----------------- ----------------- -----------------
88
OXYGEN 0.000000000000 0.000000000000 -0.131510165098
89
HYDROGEN 0.000000000000 1.479940996653 1.043580877476
92
-Geometry in the canonical coordinate system (a.u.):
94
------------ ----------------- ----------------- -----------------
95
OXYGEN 0.000000000000 0.000000000000 -0.131510165098
96
HYDROGEN 0.000000000000 1.479940996653 1.043580877476
97
HYDROGEN 0.000000000000 -1.479940996653 1.043580877476
100
-Geometry in the canonical coordinate system (Angstrom):
102
------------ ----------------- ----------------- -----------------
103
OXYGEN 0.000000000000 0.000000000000 -0.069592187382
104
HYDROGEN 0.000000000000 0.783151105291 0.552239257852
105
HYDROGEN 0.000000000000 -0.783151105291 0.552239257852
108
-Geometry in the reference coordinate system (a.u.):
110
------------ ----------------- ----------------- -----------------
111
OXYGEN 0.000000000000 0.000000000000 -0.131510165098
112
HYDROGEN 0.000000000000 1.479940996653 1.043580877476
113
HYDROGEN 0.000000000000 -1.479940996653 1.043580877476
116
--------------------------------------------------------------------------
118
Nuclear Repulsion Energy (a.u.) = 8.804687290999
120
-The Interatomic Distances in angstroms:
125
2 1.0000000 0.0000000
126
3 1.0000000 1.5663022 0.0000000
129
******************************************************************************
130
tstop called on tool.chemistry.gatech.edu
131
Thu Jul 3 18:14:34 2003
133
user time = 0.10 seconds = 0.00 minutes
134
system time = 0.01 seconds = 0.00 minutes
135
total time = 0 seconds = 0.00 minutes
136
******************************************************************************
137
tstart called on tool.chemistry.gatech.edu
138
Thu Jul 3 18:14:36 2003
140
--------------------------------------------
141
CINTS: An integrals program written in C
142
Justin T. Fermann and Edward F. Valeev
143
--------------------------------------------
148
Integral tolerance = 1e-15
149
Max. memory to use = 2500000 double words
150
Number of threads = 1
151
LIBINT's real type length = 64 bit
153
-CALCULATION CONSTANTS:
154
Label = 6-31G** SCF H2O
156
Number of atomic orbitals = 25
157
Number of symmetry orbitals = 25
158
Maximum AM in the basis = 2
160
-SYMMETRY INFORMATION;
161
Computational point group = C2v
163
Wrote 13617 two-electron integrals to IWL file 33
165
******************************************************************************
166
tstop called on tool.chemistry.gatech.edu
167
Thu Jul 3 18:14:36 2003
169
user time = 0.05 seconds = 0.00 minutes
170
system time = 0.01 seconds = 0.00 minutes
171
total time = 0 seconds = 0.00 minutes
172
******************************************************************************
173
tstart called on tool.chemistry.gatech.edu
174
Thu Jul 3 18:14:36 2003
177
------------------------------------------
179
CSCF3.0: An SCF program written in C
181
Written by too many people to mention here
183
------------------------------------------
186
label = 6-31G** SCF H2O
197
nuclear repulsion energy 8.8046872909987
198
first run, so defaulting to core-hamiltonian guess
200
level shift = 0.100000
201
diis scale factor = 1.000000
202
iterations before extrapolation = 0
203
6 error matrices will be kept
205
keeping integrals in 114320 bytes of core
207
The lowest eigenvalue of the overlap matrix was 2.424690e-02
209
Using core guess to determine occupations
212
Symmetry block: A1 A2 B1 B2
216
reading integrals in the IWL format from files 33,35,36,37
217
wrote 6516 integrals to file92
219
iter total energy delta E delta P diiser
220
1 -68.3548843162 7.715957e+01 0.000000e+00 0.000000e+00
221
2 -71.3777628421 3.022879e+00 1.130296e-01 1.175783e+00
222
3 -75.8581629361 4.480400e+00 1.037941e-01 8.860959e-01
223
4 -75.9956701999 1.375073e-01 5.906119e-03 2.044394e-01
224
5 -76.0169750721 2.130487e-02 1.931266e-03 1.024212e-01
225
6 -76.0172760736 3.010015e-04 2.842273e-04 9.429688e-03
226
7 -76.0172962933 2.021965e-05 7.939114e-05 2.474006e-03
227
8 -76.0172965444 2.510913e-07 9.894937e-06 1.856559e-04
228
9 -76.0172965714 2.696443e-08 3.411769e-06 8.255086e-05
229
10 -76.0172965718 4.839507e-10 3.849065e-07 1.190427e-05
230
11 -76.0172965718 5.911716e-12 4.827812e-08 9.566678e-07
232
Orbital energies (a.u.):
234
Doubly occupied orbitals
235
1A1 -20.568998 2A1 -1.320609 1B2 -0.678724
236
3A1 -0.563930 1B1 -0.495003
240
4A1 0.202466 2B2 0.292719 3B2 0.981164
241
5A1 1.056318 6A1 1.129180 2B1 1.168639
242
4B2 1.294986 7A1 1.413596 1A2 1.802875
243
8A1 1.806917 3B1 1.918872 9A1 2.513249
244
5B2 2.537822 6B2 2.713578 2A2 2.921065
245
4B1 2.947486 10A1 3.291755 11A1 3.620731
246
7B2 3.874917 12A1 4.077741
249
SCF total energy = -76.017296571840
250
kinetic energy = 75.657176480828
251
nuc. attr. energy = -198.046606262313
252
elec. rep. energy = 46.372133209644
253
potential energy = -151.674473052669
254
virial theorem = 1.995262655905
255
wavefunction norm = 1.000000000000
256
******************************************************************************
257
tstop called on tool.chemistry.gatech.edu
258
Thu Jul 3 18:14:36 2003
260
user time = 0.01 seconds = 0.00 minutes
261
system time = 0.01 seconds = 0.00 minutes
262
total time = 0 seconds = 0.00 minutes
263
******************************************************************************
264
tstart called on tool.chemistry.gatech.edu
265
Thu Jul 3 18:14:42 2003
267
**********************************************
269
* A simple property program *
270
* by a big TOOL fan *
271
**********************************************
274
TASKS to be performed :
275
$One-particle density will be computed from the eigenvector in checkpoint.
276
$Electric dipole, quadrupole, and octopole moments will be computed.
277
$Reference point for the electric multipole moments calculation is
279
$Reference point for the electric angular momentum calculation is
280
at (0.000000 0.000000 0.000000)
281
$Electrostatic properties at the nuclei will be evaluated.
282
will be evaluated over a rectangular 40x40 grid.
284
Title : '6-31G** SCF H2O '
288
# of molecular orbitals = 25
289
# of basis functions = 25
290
# of atomic orbitals = 25
293
# of unique shells = 12
297
List of GRID PARAMETERS :
298
GRID_ORIGIN = ( -4.00000 -4.00000 -4.00000 )
299
GRID_STEP_X = ( 0.20513 0.00000 0.00000 )
300
GRID_STEP_Y = ( 0.00000 0.20513 0.00000 )
301
GRID_STEP_Z = ( 0.00000 0.00000 0.20513 )
306
--------------------------------------------------------------
307
** Mulliken population analysis of one-particle density **
308
--------------------------------------------------------------
310
-Gross orbital populations :
313
------ ---- --- -----------
341
-Atomic bond populations :
345
1 7.0949687 0.2932585 0.2932585
346
2 0.2932585 0.3216917 -0.0234583
347
3 0.2932585 -0.0234583 0.3216917
350
-Gross atomic populations and net charges :
352
Center Atomic Population Net Charge
353
------ ----------------- ----------
359
--------------------------------------------------------------
360
*** Electric multipole moments ***
361
--------------------------------------------------------------
363
CAUTION : The system has non-vanishing dipole moment, therefore
364
quadrupole and higher moments depend on the reference point.
366
-Coordinates of the reference point (a.u.) :
368
-------------------- -------------------- --------------------
369
0.0000000000 -0.0000000000 -0.0000000000
371
-Electric dipole moment (expectation values) :
373
mu(X) = 0.00000 D = 0.00000000e+00 C*m = 0.00000000 a.u.
374
mu(Y) = -0.00000 D = -9.41286816e-45 C*m = -0.00000000 a.u.
375
mu(Z) = 2.24145 D = 7.47667986e-30 C*m = 0.88185471 a.u.
376
|mu| = 2.24145 D = 7.47667986e-30 C*m = 0.88185471 a.u.
378
-Components of electric quadrupole moment (expectation values) (a.u.) :
380
Q(XX) = -1.76092233 Q(YY) = 1.75857095 Q(ZZ) = 0.00235138
381
Q(XY) = 0.00000000 Q(XZ) = 0.00000000 Q(YZ) = 0.09791257
383
-Principal values (a.u.) and axis of electric quadrupole moment :
385
Q1 = -1.76092233 V1 = ( 1.00000000 0.00000000 0.00000000)
386
Q2 = -0.00309057 V2 = ( 0.00000000 -0.05549403 0.99845902)
387
Q3 = 1.76401290 V3 = ( 0.00000000 -0.99845902 -0.05549403)
389
-Components of electric octopole moment (expectation values) (a.u.) :
391
O(XXX) = 0.00000000 O(XXY) = 0.00000000 O(XXZ) = -1.39735155
392
O(YYY) = -0.00000000 O(XYY) = 0.00000000 O(YYZ) = 3.34015457
393
O(ZZZ) = -1.94280302 O(XZZ) = 0.00000000 O(YZZ) = 0.00000000
397
--------------------------------------------------------------
398
*** Electronic angular momentum ***
399
--------------------------------------------------------------
401
-Coordinates of the reference point (a.u.) :
403
-------------------- -------------------- --------------------
404
0.0000000000 0.0000000000 0.0000000000
406
-Electronic angular momentum (expectation values) :
411
Lx^2 = 0.00000000 a.u.
412
Ly^2 = 0.00000000 a.u.
413
Lz^2 = 0.00000000 a.u.
414
L^2 = 0.00000000 a.u.
417
--------------------------------------------------------------
418
*** Electrostatic properties at atomic centers ***
419
--------------------------------------------------------------
421
-Coordinates of atomic centers (a.u.):
423
--- ------ -------------------- -------------------- --------------------
424
1 8 0.0000000000 0.0000000000 -0.1315101651
425
2 1 0.0000000000 1.4799409967 1.0435808775
426
3 1 0.0000000000 -1.4799409967 1.0435808775
429
-Electrostatic potential and electric field (a.u.) :
432
------ ------------ ------------ ------------ ------------
433
1 -22.30375315 0.00000000 -0.00000000 -0.11900142
434
2 -1.00360324 0.00000000 -0.03191577 -0.03315459
435
3 -1.00360324 0.00000000 0.03191577 -0.03315459
438
-Electric field gradient (regular form) (a.u.):
441
------ -------------------- -------------------- --------------------
442
1 -1214.68596459 -1218.42231060 -1217.05544729
443
2 -1.34576666 -1.73327541 -1.60376479
444
3 -1.34576666 -1.73327541 -1.60376479
447
------ -------------------- -------------------- --------------------
448
1 0.00000000 0.00000000 -0.00000000
449
2 0.00000000 0.00000000 1.48137886
450
3 0.00000000 0.00000000 -1.48137886
453
-Electric field gradient (traceless tensor form) (a.u.):
455
Center XX - RR/3 YY - RR/3 ZZ - RR/3
456
------ -------------------- -------------------- --------------------
457
1 2.03527623 -1.70106977 -0.33420646
458
2 0.21516896 -0.17233979 -0.04282917
459
3 0.21516896 -0.17233979 -0.04282917
462
------ -------------------- -------------------- --------------------
463
1 0.00000000 0.00000000 -0.00000000
464
2 0.00000000 0.00000000 1.48137886
465
3 0.00000000 0.00000000 -1.48137886
468
-Electron density (a.u.):
471
------ --------------------
477
--------------------------------------------------------------
478
*** Evaluating properties over a rectangular 3D grid ***
479
--------------------------------------------------------------
481
-Coordinates of the lower left and upper right corners of
484
---- -------------------- -------------------- --------------------
485
LL -4.0000000000 -4.0000000000 -4.0000000000
486
UR 4.0000000000 4.0000000000 4.0000000000
487
--------------------------------------------------------------
488
*** Miscellaneous properties ***
489
--------------------------------------------------------------
491
-Relativistic MVD one-electron corrections to the energy (a.u.):
493
Mass-velocity (p^4) term : -0.24538837
494
One-electron Darwin term : 0.19443853
495
Total one-electron MVD terms : -0.05094984
498
NOTE : Spatial extents are computed with respect to the same reference point
499
as multipole moments.
501
-Electronic spatial extents (a.u.) :
503
<X^2> = 5.3981 <Y^2> = 7.4322 <Z^2> = 6.5390
506
-Orbital spatial extents of MOs in checkpoint (a.u.) :
508
MO # Symm <X^2> <Y^2> <Z^2> <R^2>
509
------ ---- ----------- ----------- ----------- -----------
510
2 A1 0.5041 0.7356 0.6173 0.0705
511
3 A1 0.4626 0.6907 1.4658 1.8570
512
4 A1 1.7784 5.1551 3.9671 2.6191
513
5 A1 0.9059 2.9915 1.9627 10.9006
514
6 A1 1.1214 1.3495 2.7994 5.8601
515
7 A1 1.6975 1.7126 1.8390 5.2703
516
8 A1 0.6146 1.1631 1.2997 5.2492
517
9 A1 0.8139 1.4535 0.9838 3.0774
518
10 A1 0.3299 2.1249 1.5866 3.2513
519
11 A1 0.6217 2.3871 1.4543 4.0413
520
12 A1 0.8329 1.2143 1.2301 4.4631
521
13 A1 0.8293 1.4930 0.6534 3.2773
522
14 A2 0.7979 1.9431 1.0933 2.9757
523
15 A2 1.2677 0.4420 0.4414 3.8343
524
16 B1 2.3912 0.8023 0.8151 2.1510
525
17 B1 0.9226 0.8071 1.0795 4.0087
526
18 B1 0.9660 1.9422 1.2285 2.8091
527
19 B1 0.4469 1.8301 0.7100 4.1367
528
20 B2 1.5154 6.2984 3.3961 2.9870
529
21 B2 0.8695 3.4650 1.8228 11.2099
530
22 B2 0.9173 3.0304 1.5641 6.1573
531
23 B2 0.2846 2.5910 1.7005 5.5118
532
24 B2 0.4221 2.0860 1.2580 4.5761
533
25 B2 0.4500 2.2811 1.4525 3.7661
534
26 B2 0.0000 0.0000 0.0000 4.1835
536
******************************************************************************
537
tstop called on tool.chemistry.gatech.edu
538
Thu Jul 3 18:14:43 2003
540
user time = 1.03 seconds = 0.02 minutes
541
system time = 0.01 seconds = 0.00 minutes
542
total time = 1 seconds = 0.02 minutes