[molpro-user] How to input singlet Pi state ?(1)
Chuanlu Yang
yangchuanlu at 263.net
Fri Oct 7 04:29:12 BST 2005
molpro-user :
Thanks! But I think what you said is not the key points.
The key problem is how to specify the singlet in molpro ? For Sigma state, singlet means close shell. However, Pi singlet cannot the case. In my impression, molpro usually give triplet when two orbitals are open. Obviously, we must get the two orbitals having opposite spin if we want singlet. In fact, molpro manual say the command ’open’ can specify the ‘+’ or ‘-’ of spin, but I find it no use.Now, I do not know what to do.
Chuanlu Yang , yangchuanlu at 263.net
2005-10-7
Here is the input
Geometry={angstrom;
C
O,C,r
}
basis=avtz
r=1.5
UHF
wf,14,2,0
occ,5,2,1
open,5.1,
open,-1.2
AND HERE IS THE OUTPUT
NR ATOM CHARGE X Y Z
1 C 6.00 0.000000000 0.000000000 -1.619103011
2 O 8.00 0.000000000 0.000000000 1.215485972
Bond lengths in Bohr (Angstrom)
1--2 2.834588983
(1.500000000)
NUCLEAR CHARGE: 14
NUMBER OF PRIMITIVE AOS: 134
NUMBER OF SYMMETRY AOS: 116
NUMBER OF CONTRACTIONS: 92 ( 38A1 + 22B1 + 22B2 + 10A2 )
NUMBER OF CORE ORBITALS: 2 ( 2A1 + 0B1 + 0B2 + 0A2 )
NUMBER OF VALENCE ORBITALS: 8 ( 4A1 + 2B1 + 2B2 + 0A2 )
NUCLEAR REPULSION ENERGY 16.93367197
EXTRA SYMMETRY OF AOS IN SYMMETRY 1: 1 1 1 1 1 1 1 1 1 1 2 1 2 1 2 1 2 1 2 1 1 1 1 1 1 1 1 1 1 2 1 2 1 2 1 2 1 2
EXTRA SYMMETRY OF AOS IN SYMMETRY 2: 1 1 1 1 1 1 1 1 2 1 2 1 1 1 1 1 1 1 1 2 1 2
EXTRA SYMMETRY OF AOS IN SYMMETRY 3: 1 1 1 1 1 1 1 1 2 1 2 1 1 1 1 1 1 1 1 2 1 2
Eigenvalues of metric
1 0.777E-03 0.292E-02 0.658E-02 0.121E-01 0.218E-01 0.436E-01 0.450E-01 0.797E-01
2 0.947E-02 0.368E-01 0.465E-01 0.638E-01 0.153E+00 0.190E+00 0.242E+00 0.303E+00
3 0.947E-02 0.368E-01 0.465E-01 0.638E-01 0.153E+00 0.190E+00 0.242E+00 0.303E+00
4 0.150E+00 0.233E+00 0.339E+00 0.433E+00 0.664E+00 0.942E+00 0.112E+01 0.171E+01
Contracted 2-electron integrals neglected if value below 1.0D-11
AO integral compression algorithm 1 Integral accuracy 1.0D-11
12.583 MB (compressed) written to integral file ( 48.3%)
Node minimum: 5.767 MB, node maximum: 6.816 MB
NUMBER OF SORTED TWO-ELECTRON INTEGRALS: 1202747. BUFFER LENGTH: 32768
NUMBER OF SEGMENTS: 1 SEGMENT LENGTH: 1202747 RECORD LENGTH: 524288
Memory used in sort: 1.76 MW
SORT1 READ 3201232. AND WROTE 925318. INTEGRALS IN 3 RECORDS. CPU TIME: 0.24 SEC, REAL TIME: 0.38 SEC
SORT2 READ 1856399. AND WROTE 2415093. INTEGRALS IN 44 RECORDS. CPU TIME: 0.09 SEC, REAL TIME: 0.96 SEC
Node minimum: 1202747. Node maximum: 1212346. integrals
OPERATOR DM FOR CENTER 0 COORDINATES: 0.000000 0.000000 0.000000
**********************************************************************************************************************************
DATASETS * FILE NREC LENGTH (MB) RECORD NAMES
1 18 3.05 600 500 700 960 900 950 970 1000 1100 1400
BASINP VAR GEOM ABASIS SYMINP ZMAT AOBASIS BASIS S T
1410 1200 1210 1080 1600 129 1650 1700
V H0 H01 AOSYM SMH P2S MOLCAS OPER
PROGRAMS * TOTAL INT
CPU TIMES * 1.82 1.75
REAL TIME * 2.92 SEC
DISK USED * 34.19 MB
GA USED * 0.00 MB (max) 0.00 MB (current)
**********************************************************************************************************************************
SYMMETRIES FROM ORBITAL RESTRICTIONS AND WF CARD NOT CONSISTENT: 2 2 MS2= 0 NOPEN= 2
ERROR EXIT
CURRENT STACK: MAIN
**********************************************************************************************************************************
DATASETS * FILE NREC LENGTH (MB) RECORD NAMES
1 18 3.05 600 500 700 960 900 950 970 1000 1100 1400
BASINP VAR GEOM ABASIS SYMINP ZMAT AOBASIS BASIS S T
1410 1200 1210 1080 1600 129 1650 1700
V H0 H01 AOSYM SMH P2S MOLCAS OPER
PROGRAMS * TOTAL UHF-SCF INT
CPU TIMES * 1.83 0.00 1.75
REAL TIME * 2.92 SEC
DISK USED * 34.19 MB
GA USED * 0.00 MB (max) 0.00 MB (current)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.molpro.net/pipermail/molpro-user/attachments/20051007/5454a981/attachment.html>
More information about the Molpro-user
mailing list