[molpro-user] CASSCF calculation

Jeremy Merritt merritjm at unc.edu
Mon Oct 31 16:52:20 GMT 2005


Dear molpro users.
I am new to multi-reference calculations so I hoped I could get some help on a problem.

I am trying to calculate the 3 long range PES is Cs symmetry (1A' ,2A', 1A") for the Gallium (2P) + HCN interaction neglecting spin-orbit.  The ground electronic state of Ga-HCN (van der waals complex) should be doublet sigma, where as HCN-Ga should be doublet pi.  I have confirmed this with single reference calculations.  So I would like to do a CASSCF(1,3) calculation with all three states weighted equally.  I would like to take this to higher levels (MRCI, or MRACPF etc..) but my question starts with the CASSCF.

Here is my input:
  gprint,civectors
  mH=1.00794
  mC=12.0
  mN = 14.01
  rCH = 1.06
  rCN = 1.17
  R = 3.5
  theta = 0.0
  rCM = ((mN*rCN)-(mH*rCH))/(mH+mC+mN)
  
   geometry= {Z,ANGSTROM,
    C,
    X, C, 1.0,
    N, C, rCN, X, 90,
    H, C, rCH, X, 90, N, 180,
    X2, C, rCM, X, 90, N, 0,
    Ga, X2, R, N, theta, X, 0,
  }
 basis =cc-pVDZ
  i=0;
  Do theta=0.0,180.0,180.0,
  Do R=6.5,6.5,-0.2,
  i=i+1
  T(i)=theta
  Rs(i)=R
  rhf
  wf,45,2,1;
  multi,
  occ,19,6
  closed,17,5
  wf,45,1,1;STATE,2;wf,45,2,1;
  
  Ea(i)=energy(1)
  Eb(i)=energy(2)
  Ec(i)=energy(3)

the pertinent output:
        T    RS         EA             EB             EC            
      0.0   6.5   -2016.070571   -2016.069385   -2016.070571   
  180.0   6.5   -2016.071635   -2016.069526   -2016.069526  
so one can see at 0 degrees (HCN-Ga) Ea & Ec are degenerate so they are the two Pi states and Eb is the sigma state.  Also the Pi is lower than the sigma which is correct.
but at 180 degrees (Ga-HCN) now Eb & Ec are degenerate, so they are the Pi states.  Ea is the sigma state and it is predicted to be lower in energy which is correct.
For simple linear geometries it is easy to pick out the switching of the labeling, but when the states cross at nonlinear geometries i don't know how to maintain a consistent labeling.

I've tried looking at the CI vectors:

(HCN-Ga)
 CI vector for state symmetry 1
 ------------------------------

 0+ 0      1.0000000   0.0000000
 +0 0      0.0000000   1.0000000

 TOTAL ENERGIES                     -2016.07057096 -2016.06938539

 CI vector for state symmetry 2
 ------------------------------

 00 +      1.0000000

 TOTAL ENERGIES                     -2016.07057096

(Ga-HCN)

 CI vector for state symmetry 1
 ------------------------------

 0+ 0      1.0000000   0.0000000
 +0 0      0.0000000   1.0000000

 TOTAL ENERGIES                     -2016.07163470 -2016.06952643

 CI vector for state symmetry 2
 ------------------------------

 00 +      1.0000000

 TOTAL ENERGIES                     -2016.06952643

but they don't tell me much.

Any help would be appreciated.
___________________________
Jeremy Merritt
Graduate Student, Roger Miller lab
University of North Carolina
Chapel Hill

Voice:919-962-1579
Fax:919-843-6041
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.molpro.net/pipermail/molpro-user/attachments/20051031/8cc9de90/attachment.html>


More information about the Molpro-user mailing list