19.01.2015 Views

MOLPRO

MOLPRO

MOLPRO

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

21 MULTIREFERENCE RAYLEIGH SCHRÖDINGER PERTURBATION THEORY 169<br />

r=[3,4,5,6,7,8,9,10] ang<br />

i=1<br />

geometry={Li<br />

F,1,r(i)}<br />

basis=vtz,F=avtz<br />

hf<br />

do i=1,#r<br />

{multi<br />

closed,3,0,0,0<br />

occ, 5,2,2,0<br />

state,2<br />

canonical,ci}<br />

!Hartree-Fock<br />

!loop over range of bond distances<br />

!SA-CASSCF for 2 states<br />

{rs2,MIX=2,INIT<br />

state,1,1} !single state CASPT2 for reference state 1<br />

e1_caspt2(i)=energy<br />

!unmixed caspt2 energy for ground state<br />

{rs2,MIX=2<br />

state,1,2} !single state CASPT2 for reference state 2<br />

e2_caspt2(i)=energy !unmixed caspt2 energy for excited state<br />

e1_mscaspt2(i)=msenergy(1) !ms-caspt2 energy for ground state<br />

e2_mscaspt2(i)=msenergy(2) !ms-caspt2 energy for excited state<br />

enddo<br />

{table,r,e1_caspt2,e2_caspt2,e1_mscaspt2,e2_mscaspt2<br />

title,SS-SR-CASPT2 for LiF<br />

plot,file=’lif_sr_mscaspt2.plot’<br />

}<br />

http://www.molpro.net/info/current/examples/lif_sr_mscaspt2.com<br />

This produces the plot<br />

-107<br />

SS-SR-CASPT2 for LiF<br />

-107.05<br />

-107.1<br />

E1_CASPT2<br />

E2_CASPT2<br />

E1_MSCASPT2<br />

E2_MSCASPT2<br />

-107.15<br />

-107.2<br />

3 4 5 6 7 8 9 10<br />

R

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!