1
******************************************************************************
2
PSI3 started on augustus.chemistry.gatech.edu at Wed Mar 12 18:13:09 2008
4
-----------------------------------------------------------------------
5
PSI3: An Open-Source Ab Initio Electronic Structure Package
8
T. D. Crawford, C. D. Sherrill, E. F. Valeev, J. T. Fermann, R. A. King,
9
M. L. Leininger, S. T. Brown, C. L. Janssen, E. T. Seidl, J. P. Kenny,
10
and W. D. Allen, J. Comput. Chem. 28, 1610-1616 (2007)
11
-----------------------------------------------------------------------
13
PSI3 will perform a ROHF DETCI energy computation.
15
The following programs will be executed:
24
******************************************************************************
25
tstart called on augustus.chemistry.gatech.edu
26
Wed Mar 12 18:13:09 2008
32
LABEL = 6-31G** CISD H2O
37
Parsed basis sets from /theoryfs/ds/home/ashley/psi3/lib/pbasis.dat
38
Coordinates after reading z-matrices
42
1 0.0000000 0.0000000 0.0000000
43
2 0.0000000 0.0000000 1.8897260
44
3 1.8405477 0.0000000 -0.4283089
46
-Geometry before Center-of-Mass shift (a.u.):
48
------------ ----------------- ----------------- -----------------
49
OXYGEN 0.000000000000 0.000000000000 0.000000000000
50
HYDROGEN 0.000000000000 0.000000000000 1.889725988579
51
HYDROGEN 1.840547697619 0.000000000000 -0.428308866009
54
-Rotational constants (cm-1) :
55
A = 24.35463 B = 13.63610 C = 8.74166
56
It is an asymmetric top.
58
-Geometry after Center-of-Mass shift and reorientation (a.u.):
60
------------ ----------------- ----------------- -----------------
61
OXYGEN 0.000000000000 -0.131510165113 0.000000000000
62
HYDROGEN 1.479940996653 1.043580877461 -0.000000000000
63
HYDROGEN -1.479940996653 1.043580877461 -0.000000000000
66
-SYMMETRY INFORMATION:
67
Computational point group is C2v
68
Number of irr. rep. = 4
70
Number of unique atoms = 2
75
-Basis set on unique center 1:
76
( (S ( 5484.67166000 0.00183107)
77
( 825.23494600 0.01395017)
78
( 188.04695800 0.06844508)
79
( 52.96450000 0.23271434)
80
( 16.89757040 0.47019290)
81
( 5.79963534 0.35852085) )
82
(S ( 15.53961625 -0.11077755)
83
( 3.59993359 -0.14802626)
84
( 1.01376175 1.13076701) )
85
(S ( 0.27000582 1.00000000) )
86
(P ( 15.53961625 0.07087427)
87
( 3.59993359 0.33975284)
88
( 1.01376175 0.72715858) )
89
(P ( 0.27000582 1.00000000) )
90
(D ( 0.80000000 1.00000000) )
93
-Basis set on unique center 2:
94
( (S ( 18.73113696 0.03349460)
95
( 2.82539437 0.23472695)
96
( 0.64012169 0.81375733) )
97
(S ( 0.16127776 1.00000000) )
98
(P ( 1.10000000 1.00000000) )
102
-BASIS SET INFORMATION:
103
Total number of shells = 12
104
Number of primitives = 20
116
-Unique atoms in the canonical coordinate system (a.u.):
118
------------ ----------------- ----------------- -----------------
119
OXYGEN 0.000000000000 0.000000000000 -0.131510165113
120
HYDROGEN -0.000000000000 1.479940996653 1.043580877461
123
-Geometry in the canonical coordinate system (a.u.):
125
------------ ----------------- ----------------- -----------------
126
OXYGEN 0.000000000000 0.000000000000 -0.131510165113
127
HYDROGEN -0.000000000000 1.479940996653 1.043580877461
128
HYDROGEN -0.000000000000 -1.479940996653 1.043580877461
131
-Geometry in the canonical coordinate system (Angstrom):
133
------------ ----------------- ----------------- -----------------
134
OXYGEN 0.000000000000 0.000000000000 -0.069592187390
135
HYDROGEN -0.000000000000 0.783151105291 0.552239257844
136
HYDROGEN -0.000000000000 -0.783151105291 0.552239257844
139
-Geometry in the reference coordinate system (a.u.):
141
------------ ----------------- ----------------- -----------------
142
OXYGEN 0.000000000000 0.000000000000 -0.131510165113
143
HYDROGEN -0.000000000000 1.479940996653 1.043580877461
144
HYDROGEN -0.000000000000 -1.479940996653 1.043580877461
147
--------------------------------------------------------------------------
149
Nuclear Repulsion Energy (a.u.) = 8.804687290999
151
-The Interatomic Distances in angstroms:
156
2 1.0000000 0.0000000
157
3 1.0000000 1.5663022 0.0000000
159
Note: To print *all* bond angles, out-of-plane
160
angles, and torsion angles set print = 3
163
******************************************************************************
164
tstop called on augustus.chemistry.gatech.edu
165
Wed Mar 12 18:13:09 2008
167
user time = 0.03 seconds = 0.00 minutes
168
system time = 0.01 seconds = 0.00 minutes
169
total time = 0 seconds = 0.00 minutes
170
******************************************************************************
171
tstart called on augustus.chemistry.gatech.edu
172
Wed Mar 12 18:13:09 2008
174
--------------------------------------------
175
CINTS: An integrals program written in C
176
Justin T. Fermann and Edward F. Valeev
177
--------------------------------------------
182
Integral tolerance = 1e-15
183
Max. memory to use = 2500000 double words
184
Number of threads = 1
185
LIBINT's real type length = 64 bit
187
-CALCULATION CONSTANTS:
188
Label = 6-31G** CISD H2O
190
Number of atomic orbitals = 25
191
Number of symmetry orbitals = 25
192
Maximum AM in the basis = 2
194
-SYMMETRY INFORMATION;
195
Computational point group = C2v
198
Wrote 13617 two-electron integrals to IWL file 33
200
******************************************************************************
201
tstop called on augustus.chemistry.gatech.edu
202
Wed Mar 12 18:13:09 2008
204
user time = 0.02 seconds = 0.00 minutes
205
system time = 0.00 seconds = 0.00 minutes
206
total time = 0 seconds = 0.00 minutes
207
******************************************************************************
208
tstart called on augustus.chemistry.gatech.edu
209
Wed Mar 12 18:13:09 2008
212
------------------------------------------
214
CSCF3.0: An SCF program written in C
216
Written by too many people to mention here
218
------------------------------------------
220
Cannot check consistency of the multiplicity
222
and number of electrons, double check
225
label = 6-31G** CISD H2O
236
nuclear repulsion energy 8.8046872909987
237
first run, so defaulting to core-hamiltonian guess
239
level shift = 1.000000
241
level shifting will stop after 10 cycles
242
diis scale factor = 1.020000
243
iterations before extrapolation = 0
244
4 error matrices will be kept
246
keeping integrals in 171480 bytes of core
248
The lowest eigenvalue of the overlap matrix was 2.336125e-02
251
Using DOCC and SOCC to
252
determine occupations
255
Symmetry block: A1 A2 B1 B2
260
open-shell energy coeffs
261
open shell pair alpha beta
262
1 1 0.000000 -1.000000
263
reading integrals in the IWL format from files 33,35,36,37
265
47 integrals written to file92 in 1 buffers
266
6469 integrals written to file93 in 1 buffers
267
wrote 0 integrals to file92
269
iter total energy delta E delta P diiser
270
1 -70.1800495418 7.898474e+01 0.000000e+00 0.000000e+00
271
2 -75.5626439347 5.382594e+00 4.364931e-02 1.077742e+00
272
3 -75.6058710292 4.322709e-02 1.922631e-03 1.338049e-01
273
4 -75.6172580803 1.138705e-02 9.382708e-04 5.074445e-02
274
5 -75.6207436397 3.485559e-03 6.040278e-04 2.590645e-02
275
6 -75.6212620053 5.183656e-04 2.472676e-04 9.514286e-03
276
7 -75.6213289882 6.698285e-05 9.727534e-05 3.294634e-03
277
8 -75.6213393189 1.033070e-05 4.556339e-05 1.838791e-03
278
9 -75.6213414982 2.179271e-06 2.431360e-05 1.002525e-03
279
10 -75.6213420534 5.552553e-07 1.324908e-05 5.256503e-04
280
11 -75.6213421977 1.443179e-07 7.008908e-06 2.648454e-04
281
12 -75.6213422349 3.720626e-08 3.651811e-06 1.311746e-04
282
13 -75.6213422445 9.555990e-09 1.896921e-06 6.372376e-05
283
14 -75.6213422470 2.532772e-09 9.974065e-07 3.070339e-05
284
15 -75.6213422477 6.992735e-10 5.335607e-07 1.470574e-05
285
16 -75.6213422479 1.988951e-10 2.888512e-07 6.973979e-06
286
17 -75.6213422480 5.773870e-11 1.578632e-07 3.861872e-06
287
18 -75.6213422480 1.678302e-11 8.665503e-08 2.225833e-06
288
19 -75.6213422480 4.888534e-12 4.750228e-08 1.255708e-06
289
20 -75.6213422480 1.364242e-12 2.592910e-08 6.967700e-07
290
21 -75.6213422480 3.836931e-13 1.408333e-08 3.820852e-07
291
22 -75.6213422480 1.563194e-13 7.618907e-09 2.078137e-07
292
23 -75.6213422480 2.842171e-14 4.108980e-09 1.123135e-07
293
24 -75.6213422480 5.684342e-14 2.207817e-09 6.031696e-08
294
25 -75.6213422480 0.000000e+00 1.179750e-09 3.215727e-08
295
26 -75.6213422480 0.000000e+00 6.259136e-10 1.700924e-08
296
27 -75.6213422480 -1.421085e-14 3.296388e-10 8.931934e-09
297
28 -75.6213422480 1.421085e-14 1.725884e-10 4.665238e-09
298
29 -75.6213422480 -5.684342e-14 9.005764e-11 2.428719e-09
300
ci_typ is DETCI so mo vector will be rotated
302
Orbital energies (a.u.):
304
Doubly occupied orbitals
305
1A1 -21.122454 2A1 -1.807728 1B2 -1.168198
308
Singly occupied orbitals
312
4A1 -0.131097 2B2 -0.048500 3B2 0.586122
313
5A1 0.641232 6A1 0.720655 2B1 0.727665
314
4B2 0.870547 7A1 0.987059 1A2 1.346739
315
8A1 1.356305 3B1 1.458174 9A1 2.046744
316
5B2 2.120283 6B2 2.262511 2A2 2.470165
317
4B1 2.502866 10A1 2.856186 11A1 3.171210
318
7B2 3.431043 12A1 3.610004
321
* SCF total energy = -75.621342248017
322
kinetic energy = 75.322101857972
323
nuc. attr. energy = -192.747706958106
324
elec. rep. energy = 41.804262852116
325
potential energy = -150.943444105989
326
virial theorem = 1.996042910888
327
wavefunction norm = 1.000000000000
328
******************************************************************************
329
tstop called on augustus.chemistry.gatech.edu
330
Wed Mar 12 18:13:09 2008
332
user time = 0.01 seconds = 0.00 minutes
333
system time = 0.00 seconds = 0.00 minutes
334
total time = 0 seconds = 0.00 minutes
335
******************************************************************************
336
tstart called on augustus.chemistry.gatech.edu
337
Wed Mar 12 18:13:09 2008
340
**************************************************
341
* TRANSQT2: Program to transform integrals from *
342
* the SO basis to the MO basis. *
344
* Daniel, David, & Justin *
345
**************************************************
350
Wave function = DETCI
356
Delete TEI File = Yes
357
Memory (Mbytes) = 256.0
365
Number of active MOs = 24
367
Label # SOs # FZDC # DOCC # SOCC # VIRT # FZVR
368
----- ----- ------ ------ ------ ------ ------
374
Nuclear Rep. energy (chkpt) = 8.80468729099867
375
SCF energy (chkpt) = -75.62134224801731
377
Presorting SO-basis two-electron integrals.
378
Sorting File: SO Ints (pq,rs) nbuckets = 1
379
Frozen-core energy = -61.217651345039968
380
Starting first half-transformation.
381
Sorting half-transformed integrals.
382
Starting second half-transformation.
383
Two-electron integral transformation complete.
384
******************************************************************************
385
tstop called on augustus.chemistry.gatech.edu
386
Wed Mar 12 18:13:10 2008
388
user time = 0.02 seconds = 0.00 minutes
389
system time = 0.01 seconds = 0.00 minutes
390
total time = 1 seconds = 0.02 minutes
391
******************************************************************************
392
tstart called on augustus.chemistry.gatech.edu
393
Wed Mar 12 18:13:10 2008
396
*******************************************************
402
*******************************************************
408
EX LVL = 2 H0 BLOCKSIZE = 400
409
VAL EX LVL = 0 H0 GUESS SIZE= 400
410
H0COUPLINGSIZE= 0 H0 COUPLING = no
411
NPRINT = 20 MAX DET = 10000
412
MAXITER = 12 FREEZE CORE = yes
413
NUM ROOTS = 1 ICORE = 2
416
E CONV = 6 MIXED4 = yes
417
OEI FILE = 35 R4S = no
418
OEI ERASE = no REPL OTF = no
419
TEI FILE = 72 DIAG METHOD = SEM
420
PRECONDITIONER= DAVIDSON UPDATE = DAVIDSON
421
S = 0.500000 Ms0 = no
422
TEI ERASE = no MAXNVECT = 13
423
RESTART = no RESTART VECS = 0
424
GUESS VECTOR = H0BLOCK OPENTYPE = HIGHSPIN
426
COLLAPSE SIZE = 1 HD AVE = EVANGELISTI
427
LSE = no LSE ITER = 0
428
HD OTF = yes NO DFILE = no
429
MPN = no MPN SCHMIDT = no
431
WIGNER = no ZERO BLOCKS = no
432
PERT Z = 1.0000 ROOT = 0
433
PTHREADS = no NTHREADS = 1
434
EXPORT VECTOR = no NUM EXPORT = 0
435
FILTER_GUESS = no SF_RESTRICT = no
436
OPDM = no TRANS DENSITY= no
441
STATE AVERAGE = 1(1.00)
442
STATE AVERAGE = 1(1.00)
446
ORBS IN CI = 24 NUM ALP EXPL = 4
447
FROZEN CORE = 1 NUM BET = 4
448
RESTR CORE = 0 NUM BET EXPL = 3
450
RAS1 LVL = 3 A RAS3 MAX = 2
451
RAS1 MIN = 7 B RAS3 MAX = 2
452
A RAS1 LVL = 3 RAS4 LVL = 24
453
A RAS1 MIN = 3 A RAS4 MAX = 0
454
A RAS1 MAX = 5 B RAS4 MAX = 0
455
B RAS1 LVL = 3 RAS4 MAX = 0
456
B RAS1 MIN = 2 A RAS34 MAX = 2
457
B RAS1 MAX = 4 B RAS34 MAX = 2
458
RAS3 LVL = 4 RAS34 MAX = 2
463
FROZEN DOCC = 1 0 0 0
464
FROZEN UOCC = 0 0 0 0
469
*******************************************************
472
There are 1221 alpha strings
473
There are 884 beta strings
474
CI space contains 13 blocks
476
CI space requires 3857 determinants
478
Check SCF Energy from 1- and 2-electron integrals
480
SCF Energy (ref): -75.6213422480
481
Nuclear repulsion energy: 8.8046872910
482
One-electron energy: -32.0199322620
483
Two-electron energy: 9.9231725036
484
Frozen core energy: -61.2176513450
485
Total electronic energy: -83.3144111035
486
Total SCF energy: -74.5097238125
488
*** Calculated Energy Differs from SCF Energy in CHKPT ! ***
490
CI vector/subblock length = 1396
492
*** H0 Block Eigenvalue = -75.70469197
494
Find the roots by the Simultaneous Expansion Method (Block Davidson Method)
495
Energy convergence = 1e-06
496
RMS CI vector convergence = 0.0001
498
Using 1 initial trial vectors
499
Iter 0 Root 1 = -75.704691967 Delta_E -2.329E+01 Delta_C 6.333E-01
500
Iter 1 Root 1 = -75.783964002 Delta_E -7.927E-02 Delta_C 6.266E-02
501
Iter 2 Root 1 = -75.784975400 Delta_E -1.011E-03 Delta_C 1.703E-02
502
Iter 3 Root 1 = -75.785060973 Delta_E -8.557E-05 Delta_C 4.425E-03
503
Iter 4 Root 1 = -75.785066807 Delta_E -5.835E-06 Delta_C 1.017E-03
504
Iter 5 Root 1 = -75.785067181 Delta_E -3.736E-07 Delta_C 2.587E-04
505
Iter 6 Root 1 = -75.785067203 Delta_E -2.241E-08 Delta_C 6.476E-05 c
507
* ROOT 1 CI total energy = -75.7850672033928
510
The 20 most important determinants
512
1 0.977151 ( 302, 641) 2A1 X 3A1 X 1B1 A 1B2 X
513
2 -0.037409 ( 324, 667) 2A1 X 3A1 X 1B1 A 2B2 X
514
3 0.024331 ( 307, 646) 2A1 X 1B1 A 1B2 X 6A1 X
515
4 0.024282 ( 936, 239) 2A1 A 3A1 B 1B1 X 1B2 X 3B1 A
516
5 -0.024185 ( 923, 239) 2A1 A 3A1 X 1B1 B 1B2 X 4A1 A
517
6 0.024039 ( 932, 239) 2A1 A 3A1 X 1B1 X 1B2 B 1A2 A
518
* 7 -0.023658 ( 326, 669) 2A1 X 3A1 X 1B1 A 4B2 X
519
8 0.023298 ( 321, 667) 2A1 X 3A1 X 1B2 A 2B1 A 2B2 B
520
9 0.023250 ( 326, 667) 2A1 X 3A1 X 1B1 A 2B2 B 4B2 A
521
10 -0.023146 ( 325, 668) 2A1 X 3A1 X 1B1 A 3B2 X
522
11 -0.022949 ( 321, 669) 2A1 X 3A1 X 1B2 A 2B1 A 4B2 B
523
12 -0.022588 ( 7, 440) 2A1 X 3A1 X 1B1 B 1B2 A 2B2 A
524
13 0.022379 ( 324, 669) 2A1 X 3A1 X 1B1 A 2B2 A 4B2 B
525
* 14 -0.022277 ( 928, 239) 2A1 A 3A1 X 1B1 B 1B2 X 9A1 A
526
15 0.021797 ( 305, 644) 2A1 X 1B1 A 1B2 X 5A1 X
527
16 -0.021693 ( 321, 646) 2A1 X 3A1 A 1B2 X 6A1 B 2B1 A
528
17 0.020706 ( 924, 239) 2A1 A 3A1 X 1B1 B 1B2 X 5A1 A
529
18 0.020347 ( 305, 667) 2A1 X 3A1 B 1B1 A 1B2 A 5A1 A 2B2 B
530
* 19 0.019749 ( 311, 650) 2A1 X 1B1 A 1B2 X 8A1 X
531
20 -0.019717 ( 324, 644) 2A1 X 3A1 A 1B1 A 1B2 B 5A1 B 2B2 A
536
Total Time (s) %Time %Relative
537
-----------------------------------------------------
547
"A good bug is a dead bug"
551
"I didn't write FORTRAN. That's the problem."
555
******************************************************************************
556
tstop called on augustus.chemistry.gatech.edu
557
Wed Mar 12 18:13:11 2008
559
user time = 1.30 seconds = 0.02 minutes
560
system time = 0.03 seconds = 0.00 minutes
561
total time = 1 seconds = 0.02 minutes
563
--------------------------
564
PSI3 Computation Completed
565
--------------------------
567
PSI3 stopped on augustus.chemistry.gatech.edu at Wed Mar 12 18:13:11 2008
569
Total PSI3 wall time 2 seconds = 0.03 minutes
570
******************************************************************************
b'\\ No newline at end of file'