Firefly and PC GAMESS-related discussion club


 
Learn how to ask questions correctly  
 
 
We are NATO-free zone
 



WF identification after MCSCF calculations

Andrey
andrei_ipatov@mail.ru


Hi, collagues

I need to do some MCSCF calculations of electronic spectra of diatomic molecules, for instance, the O2 molecule.

I used the input file like that
! ================================================================    
! Input file for FireFly    
! ================================================================    
$CONTRL SCFTYP=MCSCF RUNTYP=ENERGY EXETYP=RUN $END  
$CONTRL COORD=CART UNITS=ANGS $END  
$BASIS GBASIS=N311 NGAUSS=6 NDFUNC=1 $END  
$MCSCF CISTEP=ALDET SOSCF=1 $END

$DET NSTATE=10 NCORE=4 NELS=8 NACT=6 $END
$DRT GROUP=D4H NPRT=1 SOCI=1 $END  
 
$DATA  
Molecule specification  
Dnh 4  
 
O  8.000000 0.000000 0.000000 -0.6035  
O  8.000000 0.000000 0.000000  0.6035  
$END  
 
$GUESS GUESS=MOREAD NORB=38 $END
$VEC
===
$END

In the output file I have the printout like



CI EIGENVECTORS WILL BE LABELED IN GROUP=C1      
PRINTING CI COEFFICIENTS LARGER THAN  0.050000

STATE   1  ENERGY=      -149.7365263188  S=  1.00  SZ=  0.00  SPACE SYM=A  

ALPHA  | BETA   | COEFFICIENT
--------|--------|------------
111010 | 111100 |   0.6865861
111100 | 111010 |  -0.6865861
101110 | 011110 |  -0.1133390
011110 | 101110 |   0.1133390
110011 | 110101 |  -0.0603978
110101 | 110011 |   0.0603978

STATE   2  ENERGY=      -149.7035781216  S=  0.00  SZ=  0.00  SPACE SYM=A  

ALPHA  | BETA   | COEFFICIENT
--------|--------|------------
111010 | 111010 |   0.6454186
111100 | 111100 |  -0.6454186
111100 | 111010 |   0.2134040
111010 | 111100 |   0.2134040
011110 | 101110 |  -0.1514058
101110 | 011110 |  -0.1514058
110011 | 110011 |  -0.0587975
110101 | 110101 |   0.0587975


The questions are the following :

1) Why the symmetry is labelled as C1 for all the states, while
in the input I have used the D4H ?

2) Is there any option to print out the total angular momentum
projections of the states (and any other properties that
can be used to identify the states in the spectrum) ?

Thanks for any help in advance.


[ Previous ] [ Next ] [ Index ]           Mon Feb 27 '12 4:42pm
[ Reply ] [ Edit ] [ Delete ]           This message read 934 times