1676
1805
Scaling displacements by 1.000000
1678
1807
Internal Coordinate Update in Ang or Rad, aJ/Ang or aJ/Rad
1679
Value Force Displacement New Value
1680
1 1.1453317228 -0.0000439810 -0.0000020456 1.1453296772
1681
2 1.1453317228 -0.0000439810 -0.0000020456 1.1453296772
1682
3 3.1416112536 -0.0000224923 -0.0000185478 3.1415927057
1683
4 3.1415926527 -0.0000000000 -0.0000000000 3.1415926527
1684
MAX force: 0.0000439810 RMS force: 0.0000330702
1686
Back-transformation to cartesian coordinates...
1687
Iter RMS Delta(dx) RMS Delta(dq)
1688
2 0.000003731693 0.000002061476
1689
3 0.000000801340 0.000000458243
1690
4 0.000000178129 0.000000101862
1691
5 0.000000039596 0.000000022643
1692
Convergence to displaced geometry took 5 iterations.
1694
New Cartesian Geometry in a.u.
1695
6.0 0.0000000005 0.0000000000 -0.0056082387
1696
8.0 -0.0000000005 2.1643592566 -0.0056083442
1697
8.0 -0.0000000005 -2.1643592566 -0.0056083442
1699
Geometry written to chkpt
1701
******** OPTKING execution completed ********
1703
******************************************************************************
1704
tstart called on diadem.giga.net
1705
Fri Apr 23 14:43:12 2004
1707
--------------------------------------------
1708
CINTS: An integrals program written in C
1709
Justin T. Fermann and Edward F. Valeev
1710
--------------------------------------------
1715
Integral tolerance = 1e-15
1716
Max. memory to use = 2500000 double words
1717
Number of threads = 1
1718
LIBINT's real type length = 64 bit
1720
-CALCULATION CONSTANTS:
1721
Label = DZP SCF optimization of CO2
1723
Number of atomic orbitals = 48
1724
Number of symmetry orbitals = 48
1725
Maximum AM in the basis = 2
1727
-SYMMETRY INFORMATION;
1728
Computational point group = C2v
1729
Number of irreps = 4
1731
Wrote 167593 two-electron integrals to IWL file 33
1733
******************************************************************************
1734
tstop called on diadem.giga.net
1735
Fri Apr 23 14:43:13 2004
1737
user time = 0.26 seconds = 0.00 minutes
1738
system time = 0.06 seconds = 0.00 minutes
1739
total time = 1 seconds = 0.02 minutes
1740
******************************************************************************
1741
tstart called on diadem.giga.net
1742
Fri Apr 23 14:43:13 2004
1745
------------------------------------------
1747
CSCF3.0: An SCF program written in C
1749
Written by too many people to mention here
1751
------------------------------------------
1753
I think the multiplicity is 1.
1754
If this is wrong, please specify the MULTP keyword
1756
label = DZP SCF optimization of CO2
1767
nuclear repulsion energy 59.1399046218142
1769
using old vector from file30 as initial guess
1770
energy from old vector: -187.67624072
1808
Value Force Displacement New Value
1809
1 1.14518558 0.00311209 0.00014383 1.14532941
1810
2 1.14518558 0.00311209 0.00014383 1.14532941
1811
3 3.14131351 0.00027628 0.00028152 3.14159502
1812
MAX force: 0.0031120932 RMS force: 0.0025460153
1814
Back-transformation to cartesian coordinates...
1815
Iter RMS Delta(dx) RMS Delta(dq)
1816
2 0.000080959665 0.000002716281
1817
3 0.000000739350 0.000002736809
1818
4 0.000000744892 0.000002736620
1819
5 0.000000744829 0.000002736593
1820
6 0.000000744842 0.000002736647
1821
7 0.000000744824 0.000002736458
1822
8 0.000000744767 0.000002736539
1823
9 0.000000744851 0.000002736836
1824
10 0.000000744909 0.000002736647
1825
11 0.000000744824 0.000002736458
1826
12 0.000000744769 0.000002736485
1827
13 0.000000744801 0.000002736620
1828
14 0.000000744836 0.000002736620
1829
15 0.000000744843 0.000002736647
1830
16 0.000000744844 0.000002736620
1831
17 0.000000744841 0.000002736620
1832
18 0.000000744827 0.000002736485
1833
19 0.000000744761 0.000002736404
1834
20 0.000000744782 0.000002736620
1835
21 0.000000744829 0.000002736593
1836
22 0.000000744842 0.000002736647
1837
23 0.000000744826 0.000002736458
1838
24 0.000000744764 0.000002736431
1839
25 0.000000744772 0.000002736485
1840
26 0.000000744769 0.000002736431
1841
27 0.000000744780 0.000002736593
1842
28 0.000000744834 0.000002736620
1843
29 0.000000744827 0.000002736485
1844
30 0.000000744759 0.000002736404
1845
31 0.000000744785 0.000002736647
1846
32 0.000000744840 0.000002736593
1847
33 0.000000744830 0.000002736620
1848
34 0.000000744838 0.000002736620
1849
35 0.000000744836 0.000002736620
1850
36 0.000000744844 0.000002736647
1851
37 0.000000744842 0.000002736620
1852
38 0.000000744847 0.000002736647
1853
39 0.000000744831 0.000002736485
1854
40 0.000000744770 0.000002736431
1855
41 0.000000744777 0.000002736593
1856
42 0.000000744844 0.000002736647
1857
43 0.000000744818 0.000002736431
1858
44 0.000000744765 0.000002736458
1859
45 0.000000744770 0.000002736566
1860
46 0.000000744865 0.000002736836
1861
47 0.000000744887 0.000002736593
1862
48 0.000000744843 0.000002736647
1863
49 0.000000744822 0.000002736458
1864
50 0.000000744775 0.000002736512
1865
51 0.000000744804 0.000002736620
1866
52 0.000000744851 0.000002736647
1867
53 0.000000744818 0.000002736431
1868
54 0.000000744763 0.000002736512
1869
55 0.000000744819 0.000002736647
1870
56 0.000000744826 0.000002736458
1871
57 0.000000744761 0.000002736431
1872
58 0.000000744780 0.000002736593
1873
59 0.000000744835 0.000002736620
1874
60 0.000000744823 0.000002736458
1875
Could not converge new geometry in 60 iterations.
1876
Scaling back displacements by half
1878
Back-transformation to cartesian coordinates...
1879
Iter RMS Delta(dx) RMS Delta(dq)
1880
2 0.000040479832 0.000000005604
1881
3 0.000000001928 0.000000000001
1882
Convergence to displaced geometry took 3 iterations.
1884
New Cartesian Geometry in a.u.
1885
6.0 0.0000000000 0.0000000000 -0.0001089004
1886
8.0 0.0000000000 2.1642228527 0.0000408506
1887
8.0 0.0000000000 -2.1642228527 0.0000408506
1889
Geometry written to chkpt
1891
******** OPTKING execution completed ********
1893
******************************************************************************
1894
tstart called on mac4704.local
1895
Fri Apr 6 12:38:47 2007
1897
--------------------------------------------
1898
CINTS: An integrals program written in C
1899
Justin T. Fermann and Edward F. Valeev
1900
--------------------------------------------
1905
Integral tolerance = 1e-15
1906
Max. memory to use = 2500000 double words
1907
Number of threads = 1
1908
LIBINT's real type length = 64 bit
1910
-CALCULATION CONSTANTS:
1911
Label = DZP RHF optimization of CO2
1913
Number of atomic orbitals = 48
1914
Number of symmetry orbitals = 48
1915
Maximum AM in the basis = 2
1917
-SYMMETRY INFORMATION;
1918
Computational point group = C2v
1919
Number of irreps = 4
1921
Wrote 167644 two-electron integrals to IWL file 33
1923
******************************************************************************
1924
tstop called on mac4704.local
1925
Fri Apr 6 12:38:47 2007
1927
user time = 0.52 seconds = 0.01 minutes
1928
system time = 0.05 seconds = 0.00 minutes
1929
total time = 0 seconds = 0.00 minutes
1930
******************************************************************************
1931
tstart called on mac4704.local
1932
Fri Apr 6 12:38:47 2007
1935
------------------------------------------
1937
CSCF3.0: An SCF program written in C
1939
Written by too many people to mention here
1941
------------------------------------------
1943
I think the multiplicity is 1.
1944
If this is wrong, please specify the MULTP keyword
1946
label = DZP RHF optimization of CO2
1957
nuclear repulsion energy 59.1436319104510
1959
using old vector from file30 as initial guess
1960
energy from old vector: -187.67624061
1962
level shift = 0.100000
1963
diis scale factor = 1.000000
1964
iterations before extrapolation = 0
1965
6 error matrices will be kept
1967
keeping integrals in 1245120 bytes of core
1969
The lowest eigenvalue of the overlap matrix was 1.281386e-02
1972
Reading Occupations from checkpoint file.
1974
Symmetry block: A1 A2 B1 B2
1978
reading integrals in the IWL format from files 33,35,36,37
1979
wrote 70135 integrals to file92
1981
iter total energy delta E delta P diiser
1982
1 -187.6762406820 2.468199e+02 0.000000e+00 0.000000e+00
1983
2 -187.6762406953 1.329704e-08 6.078150e-07 3.077546e-05
1984
3 -187.6762406964 1.150056e-09 1.686229e-07 8.577323e-06
1985
4 -187.6762406966 1.656417e-10 6.321361e-08 4.500182e-06
1986
5 -187.6762406966 2.268052e-11 2.396259e-08 1.456743e-06
1987
6 -187.6762406966 2.359002e-12 8.713075e-09 6.004135e-07
1988
7 -187.6762406966 -5.684342e-14 1.550982e-09 7.938565e-08
1989
8 -187.6762406966 -5.684342e-14 2.316499e-10 1.129897e-08
1990
9 -187.6762406966 8.526513e-14 4.244947e-11 2.012364e-09
1992
Correcting phases of orbitals.
1994
Orbital energies (a.u.):
1996
Doubly occupied orbitals
1997
1B2 -20.652083 1A1 -20.652055 2A1 -11.475157
1998
3A1 -1.544207 2B2 -1.488412 4A1 -0.797263
1999
3B2 -0.744022 1B1 -0.722845 5A1 -0.722845
2000
1A2 -0.543866 4B2 -0.543866
2004
6A1 0.213541 2B1 0.213541 7A1 0.228688
2005
5B2 0.353917 3B1 0.402386 8A1 0.402386
2006
9A1 0.697792 10A1 0.752135 6B2 0.785502
2007
2A2 0.867901 7B2 0.867901 11A1 1.056136
2008
4B1 1.056136 12A1 1.387513 8B2 1.390386
2009
9B2 1.425296 5B1 1.589010 13A1 1.589010
2010
3A2 1.706669 10B2 1.706669 11B2 2.122289
2011
4A2 2.122289 6B1 2.397320 14A1 2.397320
2012
15A1 2.555001 7B1 2.555001 16A1 2.728471
2013
12B2 3.150815 5A2 3.359434 13B2 3.359434
2014
17A1 3.370904 18A1 4.359292 14B2 4.525370
2015
19A1 4.527670 20A1 27.370740 21A1 45.390971
2019
SCF total energy = -187.676240696621
2020
kinetic energy = 187.478066901920
2021
nuc. attr. energy = -561.344567444238
2022
elec. rep. energy = 186.190259845697
2023
potential energy = -375.154307598541
2024
virial theorem = 1.998944065621
2025
wavefunction norm = 1.000000000000
2026
******************************************************************************
2027
tstop called on mac4704.local
2028
Fri Apr 6 12:38:47 2007
2030
user time = 0.10 seconds = 0.00 minutes
2031
system time = 0.03 seconds = 0.00 minutes
2032
total time = 0 seconds = 0.00 minutes
2033
******************************************************************************
2034
tstart called on mac4704.local
2035
Fri Apr 6 12:38:48 2007
2037
--------------------------------------------
2038
CINTS: An integrals program written in C
2039
Justin T. Fermann and Edward F. Valeev
2040
--------------------------------------------
2045
Integral tolerance = 1e-15
2046
Max. memory to use = 2500000 double words
2047
Number of threads = 1
2048
LIBINT's real type length = 64 bit
2050
-CALCULATION CONSTANTS:
2051
Label = DZP RHF optimization of CO2
2053
Number of atomic orbitals = 48
2054
Number of symmetry orbitals = 48
2055
Maximum AM in the basis = 2
2057
-SYMMETRY INFORMATION;
2058
Computational point group = C2v
2059
Number of irreps = 4
2060
Rotational invariance condition satisfied.
2061
|X cross Grad| = 0.000000000000 (it is the accuracy of the computed forces)
2065
-SCF forces in the reference frame (a.u.):
2067
------ ----------------- ----------------- -----------------
2068
1 0.000000000000 0.000000000000 -0.000029001071
2069
2 0.000000000000 -0.000189180058 0.000014500535
2070
3 0.000000000000 0.000189180058 0.000014500535
2072
******************************************************************************
2073
tstop called on mac4704.local
2074
Fri Apr 6 12:38:49 2007
2076
user time = 2.58 seconds = 0.04 minutes
2077
system time = 0.03 seconds = 0.00 minutes
2078
total time = 1 seconds = 0.02 minutes
2080
------------------------------------------------------
2081
OPTKING: for internal coordinate optimizations
2082
------------------------------------------------------
2084
Cartesian geometry and possibly gradient in a.u. with masses
2085
6.0 12.00000000 0.0000000000 0.0000000000 -0.0001089004
2086
8.0 15.99491462 0.0000000000 2.1642228527 0.0000408506
2087
8.0 15.99491462 0.0000000000 -2.1642228527 0.0000408506
2088
0.0000000000 0.0000000000 -0.0000290011
2089
0.0000000000 -0.0001891801 0.0000145005
2090
0.0000000000 0.0001891801 0.0000145005
2091
Searching for geometrical constraints...none found.
2093
Simple Internal Coordinates and Values
2095
(1 1 2) (1.14525750)
2096
(2 1 3) (1.14525750)
2098
(3 2 1 3) (179.99207096)
2100
** Taking normal optimization step. **
2102
Current SCF energy before step -187.6762406966
2104
Taking geometry step number 7
2106
BuB^t Determinant: 5.007489e-03
2108
Force Constants read from PSIF_OPTKING
2110
Performing BFGS Hessian update with previous 6 gradient(s).
2112
Scaling displacements by 1.000000
2114
Internal Coordinate Update in Ang or Rad, aJ/Ang or aJ/Rad
2115
Value Force Displacement New Value
2116
1 1.14525750 0.00155859 0.00007214 1.14532964
2117
2 1.14525750 0.00155859 0.00007214 1.14532964
2118
3 3.14145427 0.00013694 0.00013851 3.14159277
2119
MAX force: 0.0015585935 RMS force: 0.0012750399
2121
Back-transformation to cartesian coordinates...
2122
Iter RMS Delta(dx) RMS Delta(dq)
2123
2 0.000040377441 0.000000131404
2124
3 0.000000035176 0.000000134835
2125
Convergence to displaced geometry took 3 iterations.
2127
New Cartesian Geometry in a.u.
2128
6.0 0.0000000000 0.0000000000 -0.0000000892
2129
8.0 0.0000000000 2.1643591876 0.0000000335
2130
8.0 0.0000000000 -2.1643591876 0.0000000335
2132
Geometry written to chkpt
2134
******** OPTKING execution completed ********
2136
******************************************************************************
2137
tstart called on mac4704.local
2138
Fri Apr 6 12:38:49 2007
2140
--------------------------------------------
2141
CINTS: An integrals program written in C
2142
Justin T. Fermann and Edward F. Valeev
2143
--------------------------------------------
2148
Integral tolerance = 1e-15
2149
Max. memory to use = 2500000 double words
2150
Number of threads = 1
2151
LIBINT's real type length = 64 bit
2153
-CALCULATION CONSTANTS:
2154
Label = DZP RHF optimization of CO2
2156
Number of atomic orbitals = 48
2157
Number of symmetry orbitals = 48
2158
Maximum AM in the basis = 2
2160
-SYMMETRY INFORMATION;
2161
Computational point group = C2v
2162
Number of irreps = 4
2164
Wrote 167600 two-electron integrals to IWL file 33
2166
******************************************************************************
2167
tstop called on mac4704.local
2168
Fri Apr 6 12:38:50 2007
2170
user time = 0.50 seconds = 0.01 minutes
2171
system time = 0.05 seconds = 0.00 minutes
2172
total time = 1 seconds = 0.02 minutes
2173
******************************************************************************
2174
tstart called on mac4704.local
2175
Fri Apr 6 12:38:50 2007
2178
------------------------------------------
2180
CSCF3.0: An SCF program written in C
2182
Written by too many people to mention here
2184
------------------------------------------
2186
I think the multiplicity is 1.
2187
If this is wrong, please specify the MULTP keyword
2189
label = DZP RHF optimization of CO2
2200
nuclear repulsion energy 59.1399065056142
2202
using old vector from file30 as initial guess
2203
energy from old vector: -187.67624070
1772
2205
level shift = 0.100000
1773
2206
diis scale factor = 1.000000
1808
2239
Doubly occupied orbitals
1809
2240
1B2 -20.652100 1A1 -20.652071 2A1 -11.475190
1810
2241
3A1 -1.544166 2B2 -1.488371 4A1 -0.797272
1811
3B2 -0.744011 5A1 -0.722813 1B1 -0.722813
2242
3B2 -0.744011 1B1 -0.722813 5A1 -0.722813
1812
2243
1A2 -0.543862 4B2 -0.543862
1815
2246
Unoccupied orbitals
1816
6A1 0.213489 2B1 0.213489 7A1 0.228667
1817
5B2 0.353915 3B1 0.402385 8A1 0.402385
1818
9A1 0.697811 10A1 0.752030 6B2 0.785494
1819
2A2 0.867905 7B2 0.867905 4B1 1.056120
1820
11A1 1.056120 12A1 1.387490 8B2 1.390174
1821
9B2 1.425283 5B1 1.589027 13A1 1.589027
1822
3A2 1.706585 10B2 1.706585 11B2 2.122296
1823
4A2 2.122296 6B1 2.397253 14A1 2.397253
1824
15A1 2.554882 7B1 2.554882 16A1 2.728424
1825
12B2 3.150711 5A2 3.359295 13B2 3.359295
1826
17A1 3.370721 18A1 4.359272 14B2 4.525336
1827
19A1 4.527430 20A1 27.370691 21A1 45.390918
1831
SCF total energy = -187.676240724599
1832
kinetic energy = 187.477351703086
1833
nuc. attr. energy = -561.336580985421
1834
elec. rep. energy = 186.182988557736
1835
potential energy = -375.153592427685
1836
virial theorem = 1.998940254660
1837
wavefunction norm = 1.000000000000
1838
******************************************************************************
1839
tstop called on diadem.giga.net
1840
Fri Apr 23 14:43:13 2004
1842
user time = 0.05 seconds = 0.00 minutes
1843
system time = 0.03 seconds = 0.00 minutes
1844
total time = 0 seconds = 0.00 minutes
1845
******************************************************************************
1846
tstart called on diadem.giga.net
1847
Fri Apr 23 14:43:13 2004
1849
--------------------------------------------
1850
CINTS: An integrals program written in C
1851
Justin T. Fermann and Edward F. Valeev
1852
--------------------------------------------
1857
Integral tolerance = 1e-15
1858
Max. memory to use = 2500000 double words
1859
Number of threads = 1
1860
LIBINT's real type length = 64 bit
1862
-CALCULATION CONSTANTS:
1863
Label = DZP SCF optimization of CO2
1865
Number of atomic orbitals = 48
1866
Number of symmetry orbitals = 48
1867
Maximum AM in the basis = 2
1869
-SYMMETRY INFORMATION;
1870
Computational point group = C2v
1871
Number of irreps = 4
1872
Rotational invariance condition satisfied.
1873
|X cross Grad| = 0.000000000000 (it is the accuracy of the computed forces)
1877
-SCF forces in the reference frame (a.u.):
1879
------ ----------------- ----------------- -----------------
1880
1 0.000000000000 0.000000000000 0.000000020436
1881
2 0.000000000000 -0.000000020999 -0.000000010218
1882
3 -0.000000000000 0.000000020999 -0.000000010218
1884
******************************************************************************
1885
tstop called on diadem.giga.net
1886
Fri Apr 23 14:43:15 2004
1888
user time = 1.85 seconds = 0.03 minutes
1889
system time = 0.06 seconds = 0.00 minutes
1890
total time = 2 seconds = 0.03 minutes
1892
------------------------------------------------------
1893
OPTKING: for internal coordinate optimizations
1894
------------------------------------------------------
1896
Cartesian geometry and possibly gradient in a.u. with masses
1897
6.0 12.00000000 0.0000000005 0.0000000000 -0.0056082387
1898
8.0 15.99491462 -0.0000000005 2.1643592566 -0.0056083442
1899
8.0 15.99491462 -0.0000000005 -2.1643592566 -0.0056083442
1900
0.0000000000 0.0000000000 0.0000000204
1901
0.0000000000 -0.0000000210 -0.0000000102
1902
-0.0000000000 0.0000000210 -0.0000000102
1904
Simple Internal Coordinates and Values
1906
(1 1 2) (1.14532968)
1907
(2 1 3) (1.14532968)
1909
(3 2 1 3 1) (180.00000558)
1910
(4 2 1 3 2) (179.99999995)
1912
** Taking normal optimization step. **
1914
Current SCF energy before step -187.6762407246
1916
Taking geometry step number 7
1918
BuB^t Determinant: 6.221432e-04
1920
Force Constants read from PSIF_OPTKING
1922
Performing BFGS Hessian update with previous 6 gradient(s).
1924
Scaling displacements by 1.000000
1926
Internal Coordinate Update in Ang or Rad, aJ/Ang or aJ/Rad
1927
Value Force Displacement New Value
1928
1 1.1453296772 0.0000001730 0.0000000079 1.1453296851
1929
2 1.1453296772 0.0000001730 0.0000000079 1.1453296851
1930
3 3.1415927510 -0.0000001178 -0.0000000970 3.1415926541
1931
4 3.1415926527 0.0000000000 0.0000000000 3.1415926527
1932
MAX force: 0.0000001730 RMS force: 0.0000001358
1934
Back-transformation to cartesian coordinates...
1935
Iter RMS Delta(dx) RMS Delta(dq)
1936
2 0.000000019212 0.000000010779
1937
Convergence to displaced geometry took 2 iterations.
1939
New Cartesian Geometry in a.u.
1940
6.0 0.0000000005 0.0000000000 -0.0056083120
1941
8.0 -0.0000000005 2.1643592715 -0.0056082892
1942
8.0 -0.0000000005 -2.1643592715 -0.0056082892
1944
Geometry written to chkpt
1946
******** OPTKING execution completed ********
1948
******************************************************************************
1949
tstart called on diadem.giga.net
1950
Fri Apr 23 14:43:15 2004
1952
--------------------------------------------
1953
CINTS: An integrals program written in C
1954
Justin T. Fermann and Edward F. Valeev
1955
--------------------------------------------
1960
Integral tolerance = 1e-15
1961
Max. memory to use = 2500000 double words
1962
Number of threads = 1
1963
LIBINT's real type length = 64 bit
1965
-CALCULATION CONSTANTS:
1966
Label = DZP SCF optimization of CO2
1968
Number of atomic orbitals = 48
1969
Number of symmetry orbitals = 48
1970
Maximum AM in the basis = 2
1972
-SYMMETRY INFORMATION;
1973
Computational point group = C2v
1974
Number of irreps = 4
1976
Wrote 167437 two-electron integrals to IWL file 33
1978
******************************************************************************
1979
tstop called on diadem.giga.net
1980
Fri Apr 23 14:43:15 2004
1982
user time = 0.32 seconds = 0.01 minutes
1983
system time = 0.05 seconds = 0.00 minutes
1984
total time = 0 seconds = 0.00 minutes
1985
******************************************************************************
1986
tstart called on diadem.giga.net
1987
Fri Apr 23 14:43:15 2004
1990
------------------------------------------
1992
CSCF3.0: An SCF program written in C
1994
Written by too many people to mention here
1996
------------------------------------------
1998
I think the multiplicity is 1.
1999
If this is wrong, please specify the MULTP keyword
2001
label = DZP SCF optimization of CO2
2012
nuclear repulsion energy 59.1399042141966
2014
using old vector from file30 as initial guess
2015
energy from old vector: -187.67624072
2017
level shift = 0.100000
2018
diis scale factor = 1.000000
2019
iterations before extrapolation = 0
2020
6 error matrices will be kept
2022
keeping integrals in 1245120 bytes of core
2024
The lowest eigenvalue of the overlap matrix was 1.281601e-02
2027
Reading Occupations from file30
2029
Symmetry block: A1 A2 B1 B2
2033
reading integrals in the IWL format from files 33,35,36,37
2034
wrote 70090 integrals to file92
2036
iter total energy delta E delta P diiser
2037
1 -187.6762407246 2.468161e+02 0.000000e+00 0.000000e+00
2038
2 -187.6762407246 0.000000e+00 3.362789e-10 2.299502e-08
2039
3 -187.6762407246 0.000000e+00 1.055303e-10 7.227134e-09
2040
4 -187.6762407246 2.842171e-14 2.321214e-11 1.372050e-09
2041
5 -187.6762407246 -2.842171e-14 6.401981e-12 3.762157e-10
2042
6 -187.6762407246 2.842171e-14 2.108778e-12 1.086384e-10
2043
7 -187.6762407246 0.000000e+00 5.051006e-13 2.249493e-11
2044
8 -187.6762407246 -2.842171e-14 5.484130e-14 3.068118e-12
2046
Correcting phases of orbitals.
2048
Orbital energies (a.u.):
2050
Doubly occupied orbitals
2051
1B2 -20.652100 1A1 -20.652071 2A1 -11.475190
2052
3A1 -1.544166 2B2 -1.488371 4A1 -0.797272
2053
3B2 -0.744011 5A1 -0.722813 1B1 -0.722813
2054
4B2 -0.543862 1A2 -0.543862
2058
6A1 0.213489 2B1 0.213489 7A1 0.228667
2059
5B2 0.353915 3B1 0.402385 8A1 0.402385
2060
9A1 0.697811 10A1 0.752030 6B2 0.785494
2061
2A2 0.867905 7B2 0.867905 11A1 1.056120
2062
4B1 1.056120 12A1 1.387490 8B2 1.390174
2063
9B2 1.425283 5B1 1.589027 13A1 1.589027
2064
10B2 1.706585 3A2 1.706585 11B2 2.122296
2065
4A2 2.122296 6B1 2.397253 14A1 2.397253
2066
7B1 2.554882 15A1 2.554882 16A1 2.728424
2067
12B2 3.150711 13B2 3.359295 5A2 3.359295
2068
17A1 3.370721 18A1 4.359272 14B2 4.525336
2069
19A1 4.527430 20A1 27.370691 21A1 45.390918
2073
SCF total energy = -187.676240724599
2074
kinetic energy = 187.477351624874
2075
nuc. attr. energy = -561.336580112015
2076
elec. rep. energy = 186.182987762542
2077
potential energy = -375.153592349473
2078
virial theorem = 1.998940254243
2079
wavefunction norm = 1.000000000000
2080
******************************************************************************
2081
tstop called on diadem.giga.net
2082
Fri Apr 23 14:43:15 2004
2084
user time = 0.06 seconds = 0.00 minutes
2085
system time = 0.01 seconds = 0.00 minutes
2086
total time = 0 seconds = 0.00 minutes
2087
******************************************************************************
2088
tstart called on diadem.giga.net
2089
Fri Apr 23 14:43:15 2004
2091
--------------------------------------------
2092
CINTS: An integrals program written in C
2093
Justin T. Fermann and Edward F. Valeev
2094
--------------------------------------------
2099
Integral tolerance = 1e-15
2100
Max. memory to use = 2500000 double words
2101
Number of threads = 1
2102
LIBINT's real type length = 64 bit
2104
-CALCULATION CONSTANTS:
2105
Label = DZP SCF optimization of CO2
2107
Number of atomic orbitals = 48
2108
Number of symmetry orbitals = 48
2109
Maximum AM in the basis = 2
2111
-SYMMETRY INFORMATION;
2112
Computational point group = C2v
2113
Number of irreps = 4
2114
Rotational invariance condition satisfied.
2115
|X cross Grad| = 0.000000000000 (it is the accuracy of the computed forces)
2119
-SCF forces in the reference frame (a.u.):
2121
------ ----------------- ----------------- -----------------
2122
1 0.000000000000 0.000000000000 -0.000000004422
2123
2 -0.000000000000 -0.000000000315 0.000000002211
2124
3 0.000000000000 0.000000000315 0.000000002211
2126
******************************************************************************
2127
tstop called on diadem.giga.net
2128
Fri Apr 23 14:43:18 2004
2130
user time = 1.95 seconds = 0.03 minutes
2131
system time = 0.35 seconds = 0.01 minutes
2132
total time = 3 seconds = 0.05 minutes
2134
------------------------------------------------------
2135
OPTKING: for internal coordinate optimizations
2136
------------------------------------------------------
2138
Cartesian geometry and possibly gradient in a.u. with masses
2139
6.0 12.00000000 0.0000000005 0.0000000000 -0.0056083120
2140
8.0 15.99491462 -0.0000000005 2.1643592715 -0.0056082892
2141
8.0 15.99491462 -0.0000000005 -2.1643592715 -0.0056082892
2142
0.0000000000 0.0000000000 -0.0000000044
2143
-0.0000000000 -0.0000000003 0.0000000022
2144
0.0000000000 0.0000000003 0.0000000022
2146
Simple Internal Coordinates and Values
2148
(1 1 2) (1.14532969)
2149
(2 1 3) (1.14532969)
2151
(3 2 1 3 1) (179.99999879)
2152
(4 2 1 3 2) (179.99999995)
2154
** Taking normal optimization step. **
2156
Current SCF energy before step -187.6762407246
2158
Taking geometry step number 8
2160
BuB^t Determinant: 6.221432e-04
2162
Force Constants read from PSIF_OPTKING
2164
Performing BFGS Hessian update with previous 6 gradient(s).
2166
Scaling displacements by 1.000000
2168
Internal Coordinate Update in Ang or Rad, aJ/Ang or aJ/Rad
2169
Value Force Displacement New Value
2170
1 1.1453296851 0.0000000026 0.0000000001 1.1453296852
2171
2 1.1453296851 0.0000000026 0.0000000001 1.1453296852
2172
3 3.1415926325 0.0000000255 0.0000000211 3.1415926536
2173
4 3.1415926527 0.0000000000 -0.0000000000 3.1415926527
2174
MAX force: 0.0000000255 RMS force: 0.0000000129
2176
Back-transformation to cartesian coordinates...
2177
Iter RMS Delta(dx) RMS Delta(dq)
2178
2 0.000000004099 0.000000002344
2179
Convergence to displaced geometry took 2 iterations.
2181
New Cartesian Geometry in a.u.
2182
6.0 0.0000000005 0.0000000000 -0.0056082961
2183
8.0 -0.0000000005 2.1643592717 -0.0056083011
2184
8.0 -0.0000000005 -2.1643592717 -0.0056083011
2186
Geometry written to chkpt
2188
******** OPTKING execution completed ********
2190
******************************************************************************
2191
tstart called on diadem.giga.net
2192
Fri Apr 23 14:43:18 2004
2194
--------------------------------------------
2195
CINTS: An integrals program written in C
2196
Justin T. Fermann and Edward F. Valeev
2197
--------------------------------------------
2202
Integral tolerance = 1e-15
2203
Max. memory to use = 2500000 double words
2204
Number of threads = 1
2205
LIBINT's real type length = 64 bit
2207
-CALCULATION CONSTANTS:
2208
Label = DZP SCF optimization of CO2
2210
Number of atomic orbitals = 48
2211
Number of symmetry orbitals = 48
2212
Maximum AM in the basis = 2
2214
-SYMMETRY INFORMATION;
2215
Computational point group = C2v
2216
Number of irreps = 4
2218
Wrote 167346 two-electron integrals to IWL file 33
2220
******************************************************************************
2221
tstop called on diadem.giga.net
2222
Fri Apr 23 14:43:18 2004
2224
user time = 0.28 seconds = 0.00 minutes
2225
system time = 0.12 seconds = 0.00 minutes
2226
total time = 0 seconds = 0.00 minutes
2227
******************************************************************************
2228
tstart called on diadem.giga.net
2229
Fri Apr 23 14:43:18 2004
2232
------------------------------------------
2234
CSCF3.0: An SCF program written in C
2236
Written by too many people to mention here
2238
------------------------------------------
2240
I think the multiplicity is 1.
2241
If this is wrong, please specify the MULTP keyword
2243
label = DZP SCF optimization of CO2
2254
nuclear repulsion energy 59.1399042079730
2256
using old vector from file30 as initial guess
2257
energy from old vector: -187.67624072
2259
level shift = 0.100000
2260
diis scale factor = 1.000000
2261
iterations before extrapolation = 0
2262
6 error matrices will be kept
2264
keeping integrals in 1245120 bytes of core
2266
The lowest eigenvalue of the overlap matrix was 1.281601e-02
2269
Reading Occupations from file30
2271
Symmetry block: A1 A2 B1 B2
2275
reading integrals in the IWL format from files 33,35,36,37
2276
wrote 70011 integrals to file92
2278
iter total energy delta E delta P diiser
2279
1 -187.6762407246 2.468161e+02 0.000000e+00 0.000000e+00
2280
2 -187.6762407246 0.000000e+00 7.229661e-11 5.000456e-09
2281
3 -187.6762407246 0.000000e+00 2.290666e-11 1.571603e-09
2282
4 -187.6762407246 0.000000e+00 5.370639e-12 2.890352e-10
2283
5 -187.6762407246 0.000000e+00 9.069270e-13 4.747209e-11
2284
6 -187.6762407246 0.000000e+00 1.248614e-13 6.623637e-12
2285
7 -187.6762407246 0.000000e+00 3.383132e-14 1.558912e-12
2287
Correcting phases of orbitals.
2289
Orbital energies (a.u.):
2291
Doubly occupied orbitals
2292
1B2 -20.652100 1A1 -20.652071 2A1 -11.475190
2293
3A1 -1.544166 2B2 -1.488371 4A1 -0.797272
2294
3B2 -0.744011 1B1 -0.722813 5A1 -0.722813
2295
4B2 -0.543862 1A2 -0.543862
2299
6A1 0.213489 2B1 0.213489 7A1 0.228667
2300
5B2 0.353915 3B1 0.402385 8A1 0.402385
2301
9A1 0.697811 10A1 0.752030 6B2 0.785494
2302
2A2 0.867905 7B2 0.867905 4B1 1.056120
2303
11A1 1.056120 12A1 1.387490 8B2 1.390174
2304
9B2 1.425283 5B1 1.589027 13A1 1.589027
2247
2B1 0.213489 6A1 0.213489 7A1 0.228667
2248
5B2 0.353915 8A1 0.402385 3B1 0.402385
2249
9A1 0.697811 10A1 0.752030 6B2 0.785494
2250
2A2 0.867905 7B2 0.867905 4B1 1.056120
2251
11A1 1.056120 12A1 1.387490 8B2 1.390174
2252
9B2 1.425283 13A1 1.589027 5B1 1.589027
2305
2253
3A2 1.706585 10B2 1.706585 4A2 2.122296
2306
2254
11B2 2.122296 14A1 2.397253 6B1 2.397253
2307
7B1 2.554882 15A1 2.554882 16A1 2.728424
2308
12B2 3.150711 13B2 3.359295 5A2 3.359295
2309
17A1 3.370721 18A1 4.359272 14B2 4.525336
2310
19A1 4.527430 20A1 27.370691 21A1 45.390918
2314
SCF total energy = -187.676240724599
2315
kinetic energy = 187.477351623680
2316
nuc. attr. energy = -561.336580098677
2317
elec. rep. energy = 186.182987750399
2318
potential energy = -375.153592348278
2319
virial theorem = 1.998940254237
2320
wavefunction norm = 1.000000000000
2321
******************************************************************************
2322
tstop called on diadem.giga.net
2323
Fri Apr 23 14:43:18 2004
2325
user time = 0.06 seconds = 0.00 minutes
2326
system time = 0.01 seconds = 0.00 minutes
2327
total time = 0 seconds = 0.00 minutes
2328
******************************************************************************
2329
tstart called on diadem.giga.net
2330
Fri Apr 23 14:43:18 2004
2332
--------------------------------------------
2333
CINTS: An integrals program written in C
2334
Justin T. Fermann and Edward F. Valeev
2335
--------------------------------------------
2340
Integral tolerance = 1e-15
2341
Max. memory to use = 2500000 double words
2342
Number of threads = 1
2343
LIBINT's real type length = 64 bit
2345
-CALCULATION CONSTANTS:
2346
Label = DZP SCF optimization of CO2
2348
Number of atomic orbitals = 48
2349
Number of symmetry orbitals = 48
2350
Maximum AM in the basis = 2
2352
-SYMMETRY INFORMATION;
2353
Computational point group = C2v
2354
Number of irreps = 4
2355
Rotational invariance condition satisfied.
2356
|X cross Grad| = 0.000000000000 (it is the accuracy of the computed forces)
2360
-SCF forces in the reference frame (a.u.):
2362
------ ----------------- ----------------- -----------------
2363
1 0.000000000000 0.000000000000 0.000000000983
2364
2 0.000000000000 0.000000000000 -0.000000000492
2365
3 0.000000000000 -0.000000000000 -0.000000000492
2367
******************************************************************************
2368
tstop called on diadem.giga.net
2369
Fri Apr 23 14:43:21 2004
2371
user time = 1.91 seconds = 0.03 minutes
2372
system time = 0.37 seconds = 0.01 minutes
2373
total time = 3 seconds = 0.05 minutes
2375
------------------------------------------------------
2376
OPTKING: for internal coordinate optimizations
2377
------------------------------------------------------
2379
Cartesian geometry and possibly gradient in a.u. with masses
2380
6.0 12.00000000 0.0000000005 0.0000000000 -0.0056082961
2381
8.0 15.99491462 -0.0000000005 2.1643592717 -0.0056083011
2382
8.0 15.99491462 -0.0000000005 -2.1643592717 -0.0056083011
2383
0.0000000000 0.0000000000 0.0000000010
2384
0.0000000000 0.0000000000 -0.0000000005
2385
0.0000000000 -0.0000000000 -0.0000000005
2387
Simple Internal Coordinates and Values
2389
(1 1 2) (1.14532969)
2390
(2 1 3) (1.14532969)
2392
(3 2 1 3 1) (180.00000027)
2393
(4 2 1 3 2) (179.99999995)
2395
** Taking normal optimization step. **
2397
Current SCF energy before step -187.6762407246
2399
Taking geometry step number 9
2401
BuB^t Determinant: 6.221432e-04
2403
Force Constants read from PSIF_OPTKING
2405
Performing BFGS Hessian update with previous 6 gradient(s).
2407
Scaling displacements by 1.000000
2409
Internal Coordinate Update in Ang or Rad, aJ/Ang or aJ/Rad
2410
Value Force Displacement New Value
2411
1 1.1453296852 -0.0000000000 -0.0000000000 1.1453296852
2412
2 1.1453296852 -0.0000000000 -0.0000000000 1.1453296852
2413
3 3.1415926583 -0.0000000057 -0.0000000047 3.1415926536
2414
4 3.1415926527 -0.0000000000 0.0000000000 3.1415926527
2415
MAX force: 0.0000000057 RMS force: 0.0000000028
2417
Back-transformation to cartesian coordinates...
2418
Iter RMS Delta(dx) RMS Delta(dq)
2419
2 0.000000000911 0.000000000521
2420
Convergence to displaced geometry took 2 iterations.
2422
New Cartesian Geometry in a.u.
2423
6.0 0.0000000005 0.0000000000 -0.0056082996
2424
8.0 -0.0000000005 2.1643592717 -0.0056082985
2425
8.0 -0.0000000005 -2.1643592717 -0.0056082985
2427
Geometry written to chkpt
2429
******** OPTKING execution completed ********
2431
******************************************************************************
2432
tstart called on diadem.giga.net
2433
Fri Apr 23 14:43:21 2004
2435
--------------------------------------------
2436
CINTS: An integrals program written in C
2437
Justin T. Fermann and Edward F. Valeev
2438
--------------------------------------------
2443
Integral tolerance = 1e-15
2444
Max. memory to use = 2500000 double words
2445
Number of threads = 1
2446
LIBINT's real type length = 64 bit
2448
-CALCULATION CONSTANTS:
2449
Label = DZP SCF optimization of CO2
2451
Number of atomic orbitals = 48
2452
Number of symmetry orbitals = 48
2453
Maximum AM in the basis = 2
2455
-SYMMETRY INFORMATION;
2456
Computational point group = C2v
2457
Number of irreps = 4
2459
Wrote 167265 two-electron integrals to IWL file 33
2461
******************************************************************************
2462
tstop called on diadem.giga.net
2463
Fri Apr 23 14:43:21 2004
2465
user time = 0.30 seconds = 0.01 minutes
2466
system time = 0.06 seconds = 0.00 minutes
2467
total time = 0 seconds = 0.00 minutes
2468
******************************************************************************
2469
tstart called on diadem.giga.net
2470
Fri Apr 23 14:43:21 2004
2473
------------------------------------------
2475
CSCF3.0: An SCF program written in C
2477
Written by too many people to mention here
2479
------------------------------------------
2481
I think the multiplicity is 1.
2482
If this is wrong, please specify the MULTP keyword
2484
label = DZP SCF optimization of CO2
2495
nuclear repulsion energy 59.1399042079830
2497
using old vector from file30 as initial guess
2498
energy from old vector: -187.67624072
2500
level shift = 0.100000
2501
diis scale factor = 1.000000
2502
iterations before extrapolation = 0
2503
6 error matrices will be kept
2505
keeping integrals in 1245120 bytes of core
2507
The lowest eigenvalue of the overlap matrix was 1.281601e-02
2510
Reading Occupations from file30
2512
Symmetry block: A1 A2 B1 B2
2516
reading integrals in the IWL format from files 33,35,36,37
2517
wrote 69718 integrals to file92
2519
iter total energy delta E delta P diiser
2520
1 -187.6762407246 2.468161e+02 0.000000e+00 0.000000e+00
2521
2 -187.6762407246 0.000000e+00 1.606901e-11 1.111415e-09
2522
3 -187.6762407246 -2.842171e-14 5.091003e-12 3.492892e-10
2523
4 -187.6762407246 2.842171e-14 1.194281e-12 6.421326e-11
2524
5 -187.6762407246 0.000000e+00 2.074014e-13 1.048875e-11
2525
6 -187.6762407246 -2.842171e-14 2.062039e-14 1.016882e-12
2527
Correcting phases of orbitals.
2529
Orbital energies (a.u.):
2531
Doubly occupied orbitals
2532
1B2 -20.652100 1A1 -20.652071 2A1 -11.475190
2533
3A1 -1.544166 2B2 -1.488371 4A1 -0.797272
2534
3B2 -0.744011 1B1 -0.722813 5A1 -0.722813
2535
4B2 -0.543862 1A2 -0.543862
2539
2B1 0.213489 6A1 0.213489 7A1 0.228667
2540
5B2 0.353915 8A1 0.402385 3B1 0.402385
2541
9A1 0.697811 10A1 0.752030 6B2 0.785494
2542
7B2 0.867905 2A2 0.867905 4B1 1.056120
2543
11A1 1.056120 12A1 1.387490 8B2 1.390174
2544
9B2 1.425283 5B1 1.589027 13A1 1.589027
2545
10B2 1.706585 3A2 1.706585 4A2 2.122296
2546
11B2 2.122296 6B1 2.397253 14A1 2.397253
2547
15A1 2.554882 7B1 2.554882 16A1 2.728424
2548
12B2 3.150711 13B2 3.359295 5A2 3.359295
2549
17A1 3.370721 18A1 4.359272 14B2 4.525336
2550
19A1 4.527430 20A1 27.370691 21A1 45.390918
2554
SCF total energy = -187.676240724599
2555
kinetic energy = 187.477351623680
2556
nuc. attr. energy = -561.336580098696
2557
elec. rep. energy = 186.182987750418
2558
potential energy = -375.153592348279
2559
virial theorem = 1.998940254237
2560
wavefunction norm = 1.000000000000
2561
******************************************************************************
2562
tstop called on diadem.giga.net
2563
Fri Apr 23 14:43:21 2004
2565
user time = 0.03 seconds = 0.00 minutes
2566
system time = 0.03 seconds = 0.00 minutes
2567
total time = 0 seconds = 0.00 minutes
2568
******************************************************************************
2569
tstart called on diadem.giga.net
2570
Fri Apr 23 14:43:21 2004
2572
--------------------------------------------
2573
CINTS: An integrals program written in C
2574
Justin T. Fermann and Edward F. Valeev
2575
--------------------------------------------
2580
Integral tolerance = 1e-15
2581
Max. memory to use = 2500000 double words
2582
Number of threads = 1
2583
LIBINT's real type length = 64 bit
2585
-CALCULATION CONSTANTS:
2586
Label = DZP SCF optimization of CO2
2588
Number of atomic orbitals = 48
2589
Number of symmetry orbitals = 48
2590
Maximum AM in the basis = 2
2592
-SYMMETRY INFORMATION;
2593
Computational point group = C2v
2594
Number of irreps = 4
2595
Rotational invariance condition satisfied.
2596
|X cross Grad| = 0.000000000000 (it is the accuracy of the computed forces)
2600
-SCF forces in the reference frame (a.u.):
2602
------ ----------------- ----------------- -----------------
2603
1 0.000000000000 0.000000000000 -0.000000000219
2604
2 0.000000000000 -0.000000000000 0.000000000109
2605
3 0.000000000000 0.000000000000 0.000000000109
2607
******************************************************************************
2608
tstop called on diadem.giga.net
2609
Fri Apr 23 14:43:23 2004
2611
user time = 1.96 seconds = 0.03 minutes
2612
system time = 0.36 seconds = 0.01 minutes
2613
total time = 2 seconds = 0.03 minutes
2615
------------------------------------------------------
2616
OPTKING: for internal coordinate optimizations
2617
------------------------------------------------------
2619
Cartesian geometry and possibly gradient in a.u. with masses
2620
6.0 12.00000000 0.0000000005 0.0000000000 -0.0056082996
2621
8.0 15.99491462 -0.0000000005 2.1643592717 -0.0056082985
2622
8.0 15.99491462 -0.0000000005 -2.1643592717 -0.0056082985
2623
0.0000000000 0.0000000000 -0.0000000002
2624
0.0000000000 -0.0000000000 0.0000000001
2625
0.0000000000 0.0000000000 0.0000000001
2627
Simple Internal Coordinates and Values
2629
(1 1 2) (1.14532969)
2630
(2 1 3) (1.14532969)
2632
(3 2 1 3 1) (179.99999994)
2633
(4 2 1 3 2) (179.99999995)
2635
** Taking normal optimization step. **
2637
Current SCF energy before step -187.6762407246
2639
Taking geometry step number 10
2641
BuB^t Determinant: 6.221432e-04
2643
Force Constants read from PSIF_OPTKING
2645
Performing BFGS Hessian update with previous 6 gradient(s).
2647
Scaling displacements by 1.000000
2649
Internal Coordinate Update in Ang or Rad, aJ/Ang or aJ/Rad
2650
Value Force Displacement New Value
2651
1 1.1453296852 0.0000000000 0.0000000000 1.1453296852
2652
2 1.1453296852 0.0000000000 0.0000000000 1.1453296852
2653
3 3.1415926525 0.0000000013 0.0000000010 3.1415926536
2654
4 3.1415926527 0.0000000000 0.0000000000 3.1415926527
2655
MAX force: 0.0000000013 RMS force: 0.0000000006
2657
Back-transformation to cartesian coordinates...
2658
Iter RMS Delta(dx) RMS Delta(dq)
2659
2 0.000000000203 0.000000000116
2660
Convergence to displaced geometry took 2 iterations.
2662
New Cartesian Geometry in a.u.
2663
6.0 0.0000000005 0.0000000000 -0.0056082988
2664
8.0 -0.0000000005 2.1643592717 -0.0056082991
2665
8.0 -0.0000000005 -2.1643592717 -0.0056082991
2667
Geometry written to chkpt
2669
******** OPTKING execution completed ********
2671
******************************************************************************
2672
tstart called on diadem.giga.net
2673
Fri Apr 23 14:43:24 2004
2675
--------------------------------------------
2676
CINTS: An integrals program written in C
2677
Justin T. Fermann and Edward F. Valeev
2678
--------------------------------------------
2683
Integral tolerance = 1e-15
2684
Max. memory to use = 2500000 double words
2685
Number of threads = 1
2686
LIBINT's real type length = 64 bit
2688
-CALCULATION CONSTANTS:
2689
Label = DZP SCF optimization of CO2
2691
Number of atomic orbitals = 48
2692
Number of symmetry orbitals = 48
2693
Maximum AM in the basis = 2
2695
-SYMMETRY INFORMATION;
2696
Computational point group = C2v
2697
Number of irreps = 4
2699
Wrote 166897 two-electron integrals to IWL file 33
2701
******************************************************************************
2702
tstop called on diadem.giga.net
2703
Fri Apr 23 14:43:24 2004
2705
user time = 0.31 seconds = 0.01 minutes
2706
system time = 0.06 seconds = 0.00 minutes
2707
total time = 0 seconds = 0.00 minutes
2708
******************************************************************************
2709
tstart called on diadem.giga.net
2710
Fri Apr 23 14:43:24 2004
2713
------------------------------------------
2715
CSCF3.0: An SCF program written in C
2717
Written by too many people to mention here
2719
------------------------------------------
2721
I think the multiplicity is 1.
2722
If this is wrong, please specify the MULTP keyword
2724
label = DZP SCF optimization of CO2
2735
nuclear repulsion energy 59.1399042079767
2737
using old vector from file30 as initial guess
2738
energy from old vector: -187.67624072
2740
level shift = 0.100000
2741
diis scale factor = 1.000000
2742
iterations before extrapolation = 0
2743
6 error matrices will be kept
2745
keeping integrals in 1245120 bytes of core
2747
The lowest eigenvalue of the overlap matrix was 1.281601e-02
2750
Reading Occupations from file30
2752
Symmetry block: A1 A2 B1 B2
2756
reading integrals in the IWL format from files 33,35,36,37
2757
wrote 68774 integrals to file92
2759
iter total energy delta E delta P diiser
2760
1 -187.6762407246 2.468161e+02 0.000000e+00 0.000000e+00
2761
2 -187.6762407246 -2.842171e-14 3.575474e-12 2.472650e-10
2762
3 -187.6762407246 0.000000e+00 1.132538e-12 7.770369e-11
2763
4 -187.6762407246 0.000000e+00 2.655928e-13 1.428010e-11
2764
5 -187.6762407246 2.842171e-14 4.609926e-14 2.334339e-12
2766
Correcting phases of orbitals.
2768
Orbital energies (a.u.):
2770
Doubly occupied orbitals
2771
1B2 -20.652100 1A1 -20.652071 2A1 -11.475190
2772
3A1 -1.544166 2B2 -1.488371 4A1 -0.797272
2773
3B2 -0.744011 1B1 -0.722813 5A1 -0.722813
2774
4B2 -0.543862 1A2 -0.543862
2778
6A1 0.213489 2B1 0.213489 7A1 0.228667
2779
5B2 0.353915 3B1 0.402385 8A1 0.402385
2780
9A1 0.697811 10A1 0.752030 6B2 0.785494
2781
2A2 0.867905 7B2 0.867905 11A1 1.056120
2782
4B1 1.056120 12A1 1.387490 8B2 1.390174
2783
9B2 1.425283 13A1 1.589027 5B1 1.589027
2784
3A2 1.706585 10B2 1.706585 11B2 2.122296
2785
4A2 2.122296 14A1 2.397253 6B1 2.397253
2786
7B1 2.554882 15A1 2.554882 16A1 2.728424
2787
12B2 3.150711 5A2 3.359295 13B2 3.359295
2788
17A1 3.370721 18A1 4.359272 14B2 4.525336
2789
19A1 4.527430 20A1 27.370691 21A1 45.390918
2793
SCF total energy = -187.676240724599
2794
kinetic energy = 187.477351623679
2795
nuc. attr. energy = -561.336580098683
2796
elec. rep. energy = 186.182987750406
2797
potential energy = -375.153592348278
2798
virial theorem = 1.998940254237
2799
wavefunction norm = 1.000000000000
2800
******************************************************************************
2801
tstop called on diadem.giga.net
2802
Fri Apr 23 14:43:24 2004
2804
user time = 0.01 seconds = 0.00 minutes
2805
system time = 0.04 seconds = 0.00 minutes
2806
total time = 0 seconds = 0.00 minutes
2807
******************************************************************************
2808
tstart called on diadem.giga.net
2809
Fri Apr 23 14:43:24 2004
2811
--------------------------------------------
2812
CINTS: An integrals program written in C
2813
Justin T. Fermann and Edward F. Valeev
2814
--------------------------------------------
2819
Integral tolerance = 1e-15
2820
Max. memory to use = 2500000 double words
2821
Number of threads = 1
2822
LIBINT's real type length = 64 bit
2824
-CALCULATION CONSTANTS:
2825
Label = DZP SCF optimization of CO2
2827
Number of atomic orbitals = 48
2828
Number of symmetry orbitals = 48
2829
Maximum AM in the basis = 2
2831
-SYMMETRY INFORMATION;
2832
Computational point group = C2v
2833
Number of irreps = 4
2834
Rotational invariance condition satisfied.
2835
|X cross Grad| = 0.000000000000 (it is the accuracy of the computed forces)
2839
-SCF forces in the reference frame (a.u.):
2841
------ ----------------- ----------------- -----------------
2842
1 0.000000000000 0.000000000000 0.000000000049
2843
2 0.000000000000 -0.000000000000 -0.000000000024
2844
3 0.000000000000 0.000000000000 -0.000000000024
2846
******************************************************************************
2847
tstop called on diadem.giga.net
2848
Fri Apr 23 14:43:26 2004
2850
user time = 1.94 seconds = 0.03 minutes
2851
system time = 0.45 seconds = 0.01 minutes
2852
total time = 2 seconds = 0.03 minutes
2854
------------------------------------------------------
2855
OPTKING: for internal coordinate optimizations
2856
------------------------------------------------------
2858
Cartesian geometry and possibly gradient in a.u. with masses
2859
6.0 12.00000000 0.0000000005 0.0000000000 -0.0056082988
2860
8.0 15.99491462 -0.0000000005 2.1643592717 -0.0056082991
2861
8.0 15.99491462 -0.0000000005 -2.1643592717 -0.0056082991
2862
0.0000000000 0.0000000000 0.0000000000
2863
0.0000000000 -0.0000000000 -0.0000000000
2864
0.0000000000 0.0000000000 -0.0000000000
2866
Simple Internal Coordinates and Values
2868
(1 1 2) (1.14532969)
2869
(2 1 3) (1.14532969)
2871
(3 2 1 3 1) (180.00000001)
2872
(4 2 1 3 2) (179.99999995)
2874
** Taking normal optimization step. **
2876
Current SCF energy before step -187.6762407246
2878
Taking geometry step number 11
2880
BuB^t Determinant: 6.221432e-04
2882
Force Constants read from PSIF_OPTKING
2884
Performing BFGS Hessian update with previous 6 gradient(s).
2886
Scaling displacements by 1.000000
2888
Internal Coordinate Update in Ang or Rad, aJ/Ang or aJ/Rad
2889
Value Force Displacement New Value
2890
1 1.1453296852 0.0000000000 0.0000000000 1.1453296852
2891
2 1.1453296852 0.0000000000 0.0000000000 1.1453296852
2892
3 3.1415926538 -0.0000000003 -0.0000000002 3.1415926536
2893
4 3.1415926527 0.0000000000 0.0000000000 3.1415926527
2894
MAX force: 0.0000000003 RMS force: 0.0000000001
2896
Back-transformation to cartesian coordinates...
2897
Iter RMS Delta(dx) RMS Delta(dq)
2898
2 0.000000000045 0.000000000026
2899
Convergence to displaced geometry took 2 iterations.
2901
New Cartesian Geometry in a.u.
2902
6.0 0.0000000005 0.0000000000 -0.0056082990
2903
8.0 -0.0000000005 2.1643592717 -0.0056082989
2904
8.0 -0.0000000005 -2.1643592717 -0.0056082989
2906
Geometry written to chkpt
2908
******** OPTKING execution completed ********
2910
******************************************************************************
2911
tstart called on diadem.giga.net
2912
Fri Apr 23 14:43:27 2004
2914
--------------------------------------------
2915
CINTS: An integrals program written in C
2916
Justin T. Fermann and Edward F. Valeev
2917
--------------------------------------------
2922
Integral tolerance = 1e-15
2923
Max. memory to use = 2500000 double words
2924
Number of threads = 1
2925
LIBINT's real type length = 64 bit
2927
-CALCULATION CONSTANTS:
2928
Label = DZP SCF optimization of CO2
2930
Number of atomic orbitals = 48
2931
Number of symmetry orbitals = 48
2932
Maximum AM in the basis = 2
2934
-SYMMETRY INFORMATION;
2935
Computational point group = C2v
2936
Number of irreps = 4
2938
Wrote 166112 two-electron integrals to IWL file 33
2940
******************************************************************************
2941
tstop called on diadem.giga.net
2942
Fri Apr 23 14:43:27 2004
2944
user time = 0.28 seconds = 0.00 minutes
2945
system time = 0.14 seconds = 0.00 minutes
2946
total time = 0 seconds = 0.00 minutes
2947
******************************************************************************
2948
tstart called on diadem.giga.net
2949
Fri Apr 23 14:43:27 2004
2952
------------------------------------------
2954
CSCF3.0: An SCF program written in C
2956
Written by too many people to mention here
2958
------------------------------------------
2960
I think the multiplicity is 1.
2961
If this is wrong, please specify the MULTP keyword
2963
label = DZP SCF optimization of CO2
2974
nuclear repulsion energy 59.1399042079729
2976
using old vector from file30 as initial guess
2977
energy from old vector: -187.67624072
2979
level shift = 0.100000
2980
diis scale factor = 1.000000
2981
iterations before extrapolation = 0
2982
6 error matrices will be kept
2984
keeping integrals in 1245120 bytes of core
2986
The lowest eigenvalue of the overlap matrix was 1.281601e-02
2989
Reading Occupations from file30
2991
Symmetry block: A1 A2 B1 B2
2995
reading integrals in the IWL format from files 33,35,36,37
2996
wrote 66540 integrals to file92
2998
iter total energy delta E delta P diiser
2999
1 -187.6762407246 2.468161e+02 0.000000e+00 0.000000e+00
3000
2 -187.6762407246 -2.842171e-14 7.978108e-13 5.511428e-11
3001
3 -187.6762407246 0.000000e+00 2.521137e-13 1.731483e-11
3002
4 -187.6762407246 2.842171e-14 5.901059e-14 3.172221e-12
3004
Correcting phases of orbitals.
3006
Orbital energies (a.u.):
3008
Doubly occupied orbitals
3009
1B2 -20.652100 1A1 -20.652071 2A1 -11.475190
3010
3A1 -1.544166 2B2 -1.488371 4A1 -0.797272
3011
3B2 -0.744011 5A1 -0.722813 1B1 -0.722813
3012
1A2 -0.543862 4B2 -0.543862
3016
2B1 0.213489 6A1 0.213489 7A1 0.228667
3017
5B2 0.353915 3B1 0.402385 8A1 0.402385
3018
9A1 0.697811 10A1 0.752030 6B2 0.785494
3019
2A2 0.867905 7B2 0.867905 4B1 1.056120
3020
11A1 1.056120 12A1 1.387490 8B2 1.390174
3021
9B2 1.425283 13A1 1.589027 5B1 1.589027
3022
10B2 1.706585 3A2 1.706585 4A2 2.122296
3023
11B2 2.122296 6B1 2.397253 14A1 2.397253
3024
15A1 2.554882 7B1 2.554882 16A1 2.728424
3025
12B2 3.150711 13B2 3.359295 5A2 3.359295
3026
17A1 3.370721 18A1 4.359272 14B2 4.525336
3027
19A1 4.527430 20A1 27.370691 21A1 45.390918
3031
SCF total energy = -187.676240724599
3032
kinetic energy = 187.477351623680
3033
nuc. attr. energy = -561.336580098678
3034
elec. rep. energy = 186.182987750399
3035
potential energy = -375.153592348279
3036
virial theorem = 1.998940254237
3037
wavefunction norm = 1.000000000000
3038
******************************************************************************
3039
tstop called on diadem.giga.net
3040
Fri Apr 23 14:43:27 2004
3042
user time = 0.03 seconds = 0.00 minutes
3043
system time = 0.03 seconds = 0.00 minutes
3044
total time = 0 seconds = 0.00 minutes
3045
******************************************************************************
3046
tstart called on diadem.giga.net
3047
Fri Apr 23 14:43:27 2004
3049
--------------------------------------------
3050
CINTS: An integrals program written in C
3051
Justin T. Fermann and Edward F. Valeev
3052
--------------------------------------------
3057
Integral tolerance = 1e-15
3058
Max. memory to use = 2500000 double words
3059
Number of threads = 1
3060
LIBINT's real type length = 64 bit
3062
-CALCULATION CONSTANTS:
3063
Label = DZP SCF optimization of CO2
3065
Number of atomic orbitals = 48
3066
Number of symmetry orbitals = 48
3067
Maximum AM in the basis = 2
3069
-SYMMETRY INFORMATION;
3070
Computational point group = C2v
3071
Number of irreps = 4
3072
Rotational invariance condition satisfied.
3073
|X cross Grad| = 0.000000000000 (it is the accuracy of the computed forces)
3077
-SCF forces in the reference frame (a.u.):
3079
------ ----------------- ----------------- -----------------
3080
1 0.000000000000 0.000000000000 -0.000000000011
3081
2 -0.000000000000 0.000000000000 0.000000000006
3082
3 0.000000000000 -0.000000000000 0.000000000006
3084
******************************************************************************
3085
tstop called on diadem.giga.net
3086
Fri Apr 23 14:43:30 2004
3088
user time = 2.02 seconds = 0.03 minutes
3089
system time = 0.42 seconds = 0.01 minutes
3090
total time = 3 seconds = 0.05 minutes
3092
------------------------------------------------------
3093
OPTKING: for internal coordinate optimizations
3094
------------------------------------------------------
3096
Cartesian geometry and possibly gradient in a.u. with masses
3097
6.0 12.00000000 0.0000000005 0.0000000000 -0.0056082990
3098
8.0 15.99491462 -0.0000000005 2.1643592717 -0.0056082989
3099
8.0 15.99491462 -0.0000000005 -2.1643592717 -0.0056082989
3100
0.0000000000 0.0000000000 -0.0000000000
3101
-0.0000000000 0.0000000000 0.0000000000
3102
0.0000000000 -0.0000000000 0.0000000000
3104
Simple Internal Coordinates and Values
3106
(1 1 2) (1.14532969)
3107
(2 1 3) (1.14532969)
3109
(3 2 1 3 1) (180.00000000)
3110
(4 2 1 3 2) (179.99999995)
3112
** Taking normal optimization step. **
3114
Current SCF energy before step -187.6762407246
3116
Taking geometry step number 12
3118
BuB^t Determinant: 6.221432e-04
3120
Force Constants read from PSIF_OPTKING
3122
Performing BFGS Hessian update with previous 6 gradient(s).
3124
Scaling displacements by 1.000000
3126
Internal Coordinate Update in Ang or Rad, aJ/Ang or aJ/Rad
3127
Value Force Displacement New Value
3128
1 1.1453296852 -0.0000000000 -0.0000000000 1.1453296852
3129
2 1.1453296852 -0.0000000000 -0.0000000000 1.1453296852
3130
3 3.1415926535 0.0000000001 0.0000000001 3.1415926536
3131
4 3.1415926527 -0.0000000000 -0.0000000000 3.1415926527
3132
MAX force: 0.0000000001 RMS force: 0.0000000000
3134
MAX force is < 1.0e-10. Optimization is complete.
2255
7B1 2.554882 15A1 2.554882 16A1 2.728425
2256
12B2 3.150712 5A2 3.359295 13B2 3.359295
2257
17A1 3.370721 18A1 4.359272 14B2 4.525336
2258
19A1 4.527430 20A1 27.370691 21A1 45.390918
2262
SCF total energy = -187.676240724599
2263
kinetic energy = 187.477352063203
2264
nuc. attr. energy = -561.336585019827
2265
elec. rep. energy = 186.182992232026
2266
potential energy = -375.153592787801
2267
virial theorem = 1.998940256579
2268
wavefunction norm = 1.000000000000
2269
******************************************************************************
2270
tstop called on mac4704.local
2271
Fri Apr 6 12:38:50 2007
2273
user time = 0.10 seconds = 0.00 minutes
2274
system time = 0.03 seconds = 0.00 minutes
2275
total time = 0 seconds = 0.00 minutes
2276
******************************************************************************
2277
tstart called on mac4704.local
2278
Fri Apr 6 12:38:50 2007
2280
--------------------------------------------
2281
CINTS: An integrals program written in C
2282
Justin T. Fermann and Edward F. Valeev
2283
--------------------------------------------
2288
Integral tolerance = 1e-15
2289
Max. memory to use = 2500000 double words
2290
Number of threads = 1
2291
LIBINT's real type length = 64 bit
2293
-CALCULATION CONSTANTS:
2294
Label = DZP RHF optimization of CO2
2296
Number of atomic orbitals = 48
2297
Number of symmetry orbitals = 48
2298
Maximum AM in the basis = 2
2300
-SYMMETRY INFORMATION;
2301
Computational point group = C2v
2302
Number of irreps = 4
2303
Rotational invariance condition satisfied.
2304
|X cross Grad| = 0.000000000000 (it is the accuracy of the computed forces)
2308
-SCF forces in the reference frame (a.u.):
2310
------ ----------------- ----------------- -----------------
2311
1 0.000000000000 0.000000000000 -0.000000024128
2312
2 0.000000000000 -0.000000116407 0.000000012064
2313
3 0.000000000000 0.000000116407 0.000000012064
2315
******************************************************************************
2316
tstop called on mac4704.local
2317
Fri Apr 6 12:38:51 2007
2319
user time = 2.58 seconds = 0.04 minutes
2320
system time = 0.03 seconds = 0.00 minutes
2321
total time = 1 seconds = 0.02 minutes
2323
------------------------------------------------------
2324
OPTKING: for internal coordinate optimizations
2325
------------------------------------------------------
2327
Cartesian geometry and possibly gradient in a.u. with masses
2328
6.0 12.00000000 0.0000000000 0.0000000000 -0.0000000892
2329
8.0 15.99491462 0.0000000000 2.1643591876 0.0000000335
2330
8.0 15.99491462 0.0000000000 -2.1643591876 0.0000000335
2331
0.0000000000 0.0000000000 -0.0000000241
2332
0.0000000000 -0.0000001164 0.0000000121
2333
0.0000000000 0.0000001164 0.0000000121
2334
Searching for geometrical constraints...none found.
2336
Simple Internal Coordinates and Values
2338
(1 1 2) (1.14532964)
2339
(2 1 3) (1.14532964)
2341
(3 2 1 3) (179.99999344)
2343
** Taking normal optimization step. **
2345
Current SCF energy before step -187.6762407246
2347
Taking geometry step number 8
2349
BuB^t Determinant: 4.909957e-03
2351
Force Constants read from PSIF_OPTKING
2353
Performing BFGS Hessian update with previous 6 gradient(s).
2355
Scaling displacements by 1.000000
2357
Internal Coordinate Update in Ang or Rad, aJ/Ang or aJ/Rad
2358
Value Force Displacement New Value
2359
1 1.14532964 0.00000096 0.00000004 1.14532969
2360
2 1.14532964 0.00000096 0.00000004 1.14532969
2361
3 3.14159254 0.00000011 0.00000012 3.14159265
2362
MAX force: 0.0000009590 RMS force: 0.0000007859
2364
MAX force is < 1.0e-05. Optimization is complete.
3135
2365
Final SCF energy is -187.6762407246
3136
2366
The Optimized geometry in a.u.
3137
( C 0.0000000005 0.0000000000 -0.0056082990 )
3138
( O -0.0000000005 2.1643592717 -0.0056082989 )
3139
( O -0.0000000005 -2.1643592717 -0.0056082989 )
2367
( C 0.0000000000 0.0000000000 -0.0000000892 )
2368
( O 0.0000000000 2.1643591876 0.0000000335 )
2369
( O 0.0000000000 -2.1643591876 0.0000000335 )
3141
2381
******** OPTKING execution completed ********