1
******************************************************************************
2
PSI3 started on augustus.chemistry.gatech.edu at Wed Mar 12 18:13:03 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 RHF 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:03 2008
37
Parsed basis sets from /theoryfs/ds/home/ashley/psi3/lib/pbasis.dat
39
-Geometry before Center-of-Mass shift (a.u.):
41
------------ ----------------- ----------------- -----------------
42
OXYGEN 0.000000000000 0.000000000000 -0.074271925400
43
HYDROGEN 0.000000000000 -1.494958998200 -1.072864037300
44
HYDROGEN 0.000000000000 1.494958998200 -1.072864037300
47
-Rotational constants (cm-1) :
48
A = 33.72471 B = 13.36351 C = 9.57098
49
It is an asymmetric top.
51
-Geometry after Center-of-Mass shift and reorientation (a.u.):
53
------------ ----------------- ----------------- -----------------
54
OXYGEN 0.000000000000 0.111757309654 0.000000000000
55
HYDROGEN -1.494958998200 -0.886834802246 0.000000000000
56
HYDROGEN 1.494958998200 -0.886834802246 0.000000000000
59
-SYMMETRY INFORMATION:
60
Computational point group is C2v
61
Number of irr. rep. = 4
63
Number of unique atoms = 2
68
-Basis set on unique center 1:
69
( (S ( 5484.67166000 0.00183107)
70
( 825.23494600 0.01395017)
71
( 188.04695800 0.06844508)
72
( 52.96450000 0.23271434)
73
( 16.89757040 0.47019290)
74
( 5.79963534 0.35852085) )
75
(S ( 15.53961625 -0.11077755)
76
( 3.59993359 -0.14802626)
77
( 1.01376175 1.13076701) )
78
(S ( 0.27000582 1.00000000) )
79
(P ( 15.53961625 0.07087427)
80
( 3.59993359 0.33975284)
81
( 1.01376175 0.72715858) )
82
(P ( 0.27000582 1.00000000) )
85
-Basis set on unique center 2:
86
( (S ( 18.73113696 0.03349460)
87
( 2.82539437 0.23472695)
88
( 0.64012169 0.81375733) )
89
(S ( 0.16127776 1.00000000) )
93
-BASIS SET INFORMATION:
94
Total number of shells = 9
95
Number of primitives = 18
107
-Unique atoms in the canonical coordinate system (a.u.):
109
------------ ----------------- ----------------- -----------------
110
OXYGEN 0.000000000000 0.000000000000 0.111757309654
111
HYDROGEN 0.000000000000 1.494958998200 -0.886834802246
114
-Geometry in the canonical coordinate system (a.u.):
116
------------ ----------------- ----------------- -----------------
117
OXYGEN 0.000000000000 0.000000000000 0.111757309654
118
HYDROGEN 0.000000000000 -1.494958998200 -0.886834802246
119
HYDROGEN 0.000000000000 1.494958998200 -0.886834802246
122
-Geometry in the canonical coordinate system (Angstrom):
124
------------ ----------------- ----------------- -----------------
125
OXYGEN 0.000000000000 0.000000000000 0.059139425678
126
HYDROGEN 0.000000000000 -0.791098290035 -0.469292800970
127
HYDROGEN 0.000000000000 0.791098290035 -0.469292800970
130
-Geometry in the reference coordinate system (a.u.):
132
------------ ----------------- ----------------- -----------------
133
OXYGEN 0.000000000000 0.000000000000 0.111757309654
134
HYDROGEN 0.000000000000 -1.494958998200 -0.886834802246
135
HYDROGEN 0.000000000000 1.494958998200 -0.886834802246
138
--------------------------------------------------------------------------
140
Nuclear Repulsion Energy (a.u.) = 9.234218520912
142
-The Interatomic Distances in angstroms:
147
2 0.9513554 0.0000000
148
3 0.9513554 1.5821966 0.0000000
150
Note: To print *all* bond angles, out-of-plane
151
angles, and torsion angles set print = 3
154
******************************************************************************
155
tstop called on augustus.chemistry.gatech.edu
156
Wed Mar 12 18:13:03 2008
158
user time = 0.03 seconds = 0.00 minutes
159
system time = 0.00 seconds = 0.00 minutes
160
total time = 0 seconds = 0.00 minutes
161
******************************************************************************
162
tstart called on augustus.chemistry.gatech.edu
163
Wed Mar 12 18:13:03 2008
165
--------------------------------------------
166
CINTS: An integrals program written in C
167
Justin T. Fermann and Edward F. Valeev
168
--------------------------------------------
173
Integral tolerance = 1e-15
174
Max. memory to use = 2500000 double words
175
Number of threads = 1
176
LIBINT's real type length = 64 bit
178
-CALCULATION CONSTANTS:
181
Number of atomic orbitals = 13
182
Number of symmetry orbitals = 13
183
Maximum AM in the basis = 1
185
-SYMMETRY INFORMATION;
186
Computational point group = C2v
189
Wrote 1270 two-electron integrals to IWL file 33
191
******************************************************************************
192
tstop called on augustus.chemistry.gatech.edu
193
Wed Mar 12 18:13:03 2008
195
user time = 0.01 seconds = 0.00 minutes
196
system time = 0.00 seconds = 0.00 minutes
197
total time = 0 seconds = 0.00 minutes
198
******************************************************************************
199
tstart called on augustus.chemistry.gatech.edu
200
Wed Mar 12 18:13:03 2008
203
------------------------------------------
205
CSCF3.0: An SCF program written in C
207
Written by too many people to mention here
209
------------------------------------------
211
I think the multiplicity is 1.
212
If this is wrong, please specify the MULTP keyword
225
nuclear repulsion energy 9.2342185209120
226
first run, so defaulting to core-hamiltonian guess
228
level shift = 0.100000
230
level shifting will stop after 10 cycles
231
diis scale factor = 1.000000
232
iterations before extrapolation = 0
233
6 error matrices will be kept
235
keeping integrals in 13856 bytes of core
237
The lowest eigenvalue of the overlap matrix was 6.914584e-02
239
Using core guess to determine occupations
242
Symmetry block: A1 A2 B1 B2
246
reading integrals in the IWL format from files 33,35,36,37
247
wrote 771 integrals to file92
249
iter total energy delta E delta P diiser
250
1 -69.6481965652 7.888242e+01 0.000000e+00 0.000000e+00
251
2 -71.9332180777 2.285022e+00 3.366880e-01 9.628117e-01
252
3 -75.8929735818 3.959756e+00 2.921836e-01 9.449561e-01
253
4 -75.9780616540 8.508807e-02 2.170364e-02 1.710999e-01
254
5 -75.9848563776 6.794724e-03 7.834688e-03 3.772661e-02
255
6 -75.9851536550 2.972774e-04 6.469860e-04 1.020797e-02
256
7 -75.9853226530 1.689980e-04 1.072542e-03 5.710474e-03
257
8 -75.9853236468 9.937899e-07 1.412447e-04 4.011976e-04
258
9 -75.9853236650 1.820227e-08 9.540334e-06 8.341422e-05
259
10 -75.9853236653 2.444978e-10 1.003736e-06 5.871018e-06
260
11 -75.9853236653 2.899014e-12 1.551904e-07 8.200001e-07
261
12 -75.9853236653 1.421085e-14 1.023704e-08 7.852794e-08
262
13 -75.9853236653 -1.421085e-14 4.142088e-09 2.749802e-08
263
14 -75.9853236653 -1.421085e-14 2.239677e-10 1.858808e-09
264
15 -75.9853236653 0.000000e+00 1.360427e-10 6.735990e-10
265
16 -75.9853236653 -1.421085e-14 5.458203e-11 1.976273e-10
267
Orbital energies (a.u.):
269
Doubly occupied orbitals
270
1A1 -20.553032 2A1 -1.350528 1B2 -0.726863
271
3A1 -0.546143 1B1 -0.497682
275
4A1 0.206667 2B2 0.302606 3B2 1.111487
276
5A1 1.152531 2B1 1.167571 6A1 1.204413
277
4B2 1.389485 7A1 1.673595
280
* SCF total energy = -75.985323665263
281
kinetic energy = 76.005864259816
282
nuc. attr. energy = -199.096452988091
283
elec. rep. energy = 47.105265063012
284
potential energy = -151.991187925079
285
virial theorem = 2.000270323183
286
wavefunction norm = 1.000000000000
287
******************************************************************************
288
tstop called on augustus.chemistry.gatech.edu
289
Wed Mar 12 18:13:03 2008
291
user time = 0.01 seconds = 0.00 minutes
292
system time = 0.00 seconds = 0.00 minutes
293
total time = 0 seconds = 0.00 minutes
294
******************************************************************************
295
tstart called on augustus.chemistry.gatech.edu
296
Wed Mar 12 18:13:03 2008
299
**************************************************
300
* TRANSQT2: Program to transform integrals from *
301
* the SO basis to the MO basis. *
303
* Daniel, David, & Justin *
304
**************************************************
309
Wave function = DETCI
315
Delete TEI File = Yes
316
Memory (Mbytes) = 256.0
324
Number of active MOs = 11
326
Label # SOs # FZDC # DOCC # SOCC # VIRT # FZVR
327
----- ----- ------ ------ ------ ------ ------
333
Nuclear Rep. energy (chkpt) = 9.23421852091198
334
SCF energy (chkpt) = -75.98532366526261
336
Presorting SO-basis two-electron integrals.
337
Sorting File: SO Ints (pq,rs) nbuckets = 1
338
Frozen-core energy = -61.325924584078727
339
Starting first half-transformation.
340
Sorting half-transformed integrals.
341
Starting second half-transformation.
342
Two-electron integral transformation complete.
343
******************************************************************************
344
tstop called on augustus.chemistry.gatech.edu
345
Wed Mar 12 18:13:03 2008
347
user time = 0.01 seconds = 0.00 minutes
348
system time = 0.00 seconds = 0.00 minutes
349
total time = 0 seconds = 0.00 minutes
350
******************************************************************************
351
tstart called on augustus.chemistry.gatech.edu
352
Wed Mar 12 18:13:04 2008
355
*******************************************************
361
*******************************************************
365
Note: Calculation requested is a full CI.
366
Resetting EX_LVL to 8 and turning on all excitations
370
EX LVL = 8 H0 BLOCKSIZE = 400
371
VAL EX LVL = 0 H0 GUESS SIZE= 400
372
H0COUPLINGSIZE= 0 H0 COUPLING = no
373
NPRINT = 20 MAX DET = 10000
374
MAXITER = 12 FREEZE CORE = yes
375
NUM ROOTS = 1 ICORE = 1
378
E CONV = 6 MIXED4 = yes
379
OEI FILE = 35 R4S = no
380
OEI ERASE = no REPL OTF = no
381
TEI FILE = 72 DIAG METHOD = SEM
382
PRECONDITIONER= DAVIDSON UPDATE = DAVIDSON
383
S = 0.000000 Ms0 = yes
384
TEI ERASE = no MAXNVECT = 13
385
RESTART = no RESTART VECS = 0
386
GUESS VECTOR = H0BLOCK OPENTYPE = NONE
388
COLLAPSE SIZE = 1 HD AVE = EVANGELISTI
389
LSE = no LSE ITER = 0
390
HD OTF = yes NO DFILE = no
391
MPN = no MPN SCHMIDT = no
393
WIGNER = no ZERO BLOCKS = no
394
PERT Z = 1.0000 ROOT = 0
395
PTHREADS = no NTHREADS = 1
396
EXPORT VECTOR = no NUM EXPORT = 0
397
FILTER_GUESS = no SF_RESTRICT = no
398
OPDM = no TRANS DENSITY= no
402
EX_TYPE = 1 1 1 1 1 1 1 1
403
STATE AVERAGE = 1(1.00)
404
STATE AVERAGE = 1(1.00)
408
ORBS IN CI = 11 NUM ALP EXPL = 4
409
FROZEN CORE = 1 NUM BET = 5
410
RESTR CORE = 0 NUM BET EXPL = 4
412
RAS1 LVL = 3 A RAS3 MAX = 4
413
RAS1 MIN = 2 B RAS3 MAX = 4
414
A RAS1 LVL = 3 RAS4 LVL = 11
415
A RAS1 MIN = 1 A RAS4 MAX = 0
416
A RAS1 MAX = 5 B RAS4 MAX = 0
417
B RAS1 LVL = 3 RAS4 MAX = 0
418
B RAS1 MIN = 1 A RAS34 MAX = 4
419
B RAS1 MAX = 5 B RAS34 MAX = 4
420
RAS3 LVL = 4 RAS34 MAX = 8
425
FROZEN DOCC = 1 0 0 0
426
FROZEN UOCC = 1 0 0 0
431
*******************************************************
434
There are 330 alpha strings
435
There are 330 beta strings
436
CI space contains 4 blocks
438
CI space requires 27268 determinants
440
Check SCF Energy from 1- and 2-electron integrals
442
SCF Energy (ref): -75.9853236653
443
Nuclear repulsion energy: 9.2342185209
444
One-electron energy: -41.5448033078
445
Two-electron energy: 17.6511857057
446
Frozen core energy: -61.3259245841
447
Total electronic energy: -85.2195421862
448
Total SCF energy: -75.9853236653
450
CI vector/subblock length = 27268
452
*** H0 Block Eigenvalue = -76.08301143
454
Find the roots by the Simultaneous Expansion Method (Block Davidson Method)
455
Energy convergence = 1e-06
456
RMS CI vector convergence = 0.0001
458
Using 1 initial trial vectors
459
Iter 0 Root 1 = -76.083011434 Delta_E -2.399E+01 Delta_C 2.591E-01
460
Iter 1 Root 1 = -76.098854190 Delta_E -1.584E-02 Delta_C 6.103E-02
461
Iter 2 Root 1 = -76.099516700 Delta_E -6.625E-04 Delta_C 2.127E-02
462
Iter 3 Root 1 = -76.099611028 Delta_E -9.433E-05 Delta_C 7.375E-03
463
Iter 4 Root 1 = -76.099621173 Delta_E -1.014E-05 Delta_C 2.029E-03
464
Iter 5 Root 1 = -76.099621944 Delta_E -7.706E-07 Delta_C 6.445E-04
465
Iter 6 Root 1 = -76.099622022 Delta_E -7.833E-08 Delta_C 1.908E-04
466
Iter 7 Root 1 = -76.099622029 Delta_E -7.216E-09 Delta_C 6.363E-05 c
468
* ROOT 1 CI total energy = -76.0996220292430
471
The 20 most important determinants
473
1 -0.980081 ( 82, 82) 2A1 X 3A1 X 1B1 X 1B2 X
474
2 0.052723 ( 92, 92) 2A1 X 3A1 X 1B2 X 2B1 X
475
3 0.046546 ( 88, 88) 2A1 X 1B1 X 1B2 X 6A1 X
476
4 0.039161 ( 102, 102) 2A1 X 3A1 X 1B1 X 2B2 X
477
5 -0.034218 ( 88, 92) 2A1 X 3A1 B 1B1 A 1B2 X 6A1 A 2B1 B
478
6 -0.034218 ( 92, 88) 2A1 X 3A1 A 1B1 B 1B2 X 6A1 B 2B1 A
479
7 0.028966 ( 171, 171) 2A1 X 3A1 X 1B1 X 5A1 X
480
* 8 0.026793 ( 144, 144) 2A1 X 3A1 X 1B1 X 4B2 X
481
9 0.026608 ( 122, 122) 2A1 X 3A1 X 1B1 X 3B2 X
482
10 -0.026043 ( 92, 144) 2A1 X 3A1 X 1B1 B 1B2 A 2B1 A 4B2 B
483
11 -0.026043 ( 144, 92) 2A1 X 3A1 X 1B1 A 1B2 B 2B1 B 4B2 A
484
12 0.024926 ( 92, 102) 2A1 X 3A1 X 1B1 B 1B2 A 2B1 A 2B2 B
485
13 0.024926 ( 102, 92) 2A1 X 3A1 X 1B1 A 1B2 B 2B1 B 2B2 A
486
14 0.024901 ( 82, 98) 2A1 X 3A1 A 1B1 A 1B2 X 6A1 B 2B1 B
487
15 0.024901 ( 98, 82) 2A1 X 3A1 B 1B1 B 1B2 X 6A1 A 2B1 A
488
16 -0.024465 ( 102, 144) 2A1 X 3A1 X 1B1 X 2B2 A 4B2 B
489
17 -0.024465 ( 144, 102) 2A1 X 3A1 X 1B1 X 2B2 B 4B2 A
490
18 0.023528 ( 88, 144) 2A1 X 3A1 B 1B1 X 1B2 A 6A1 A 4B2 B
491
19 0.023528 ( 144, 88) 2A1 X 3A1 A 1B1 X 1B2 B 6A1 B 4B2 A
492
20 -0.023371 ( 88, 102) 2A1 X 3A1 B 1B1 X 1B2 A 6A1 A 2B2 B
497
Total Time (s) %Time %Relative
498
-----------------------------------------------------
508
"A good bug is a dead bug"
512
"I didn't write FORTRAN. That's the problem."
516
******************************************************************************
517
tstop called on augustus.chemistry.gatech.edu
518
Wed Mar 12 18:13:05 2008
520
user time = 1.15 seconds = 0.02 minutes
521
system time = 0.05 seconds = 0.00 minutes
522
total time = 1 seconds = 0.02 minutes
524
--------------------------
525
PSI3 Computation Completed
526
--------------------------
528
PSI3 stopped on augustus.chemistry.gatech.edu at Wed Mar 12 18:13:05 2008
530
Total PSI3 wall time 2 seconds = 0.03 minutes
531
******************************************************************************
b'\\ No newline at end of file'