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.

27 THE MRCC PROGRAM OF M. KALLAY (MRCC) 200<br />

***,h2o excitation energies<br />

memory,8,m<br />

gthresh,energy=1.d-8<br />

geometry={<br />

o;h1,o,r;h2,o,r,h1,theta}<br />

theta=104<br />

r=1 ang<br />

basis=vdz<br />

hf<br />

ii=0<br />

s=2 !number of states in each symmetry<br />

do sym=1,4<br />

!loop over irreps<br />

ccsd;eom,-(s+0.1*sym);$p=molpro;save_energy<br />

mrcc,method=ccsd, symm=sym,nstates=2;$p=mrcc;save_energy<br />

mrcc,method=ccsdt,symm=sym,nstates=2;$p=mrcc;save_energy<br />

s=1<br />

enddo<br />

{table,method,prog,states,e,exc<br />

sort,3}<br />

save_energy={ !procedure to save results in variables<br />

!nogprint,variable<br />

e1=energy(1)<br />

do i=1,#energy<br />

ii=ii+1<br />

e(ii)=energy(i)<br />

method(ii)=program<br />

prog(ii)=p<br />

states(ii)=i+0.1*sym<br />

exc(ii)=(e(ii)-e1)*toev<br />

end do<br />

}<br />

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

This yields<br />

METHOD PROG STATES E EXC<br />

CCSD <strong>MOLPRO</strong> 1.1 -76.23580212 0.000<br />

CCSD MRCC 1.1 -76.23580212 0.000<br />

CCSDT MRCC 1.1 -76.23922746 0.000<br />

CCSD <strong>MOLPRO</strong> 1.2 -76.23580212 0.000<br />

CCSD MRCC 1.2 -76.23580212 0.000<br />

CCSDT MRCC 1.2 -76.23922746 0.000<br />

CCSD <strong>MOLPRO</strong> 1.3 -76.23580212 0.000<br />

CCSD MRCC 1.3 -76.23580212 0.000<br />

CCSDT MRCC 1.3 -76.23922746 0.000<br />

CCSD <strong>MOLPRO</strong> 1.4 -76.23580212 0.000<br />

CCSD MRCC 1.4 -76.23580212 0.000<br />

CCSDT MRCC 1.4 -76.23922746 0.000<br />

CCSD <strong>MOLPRO</strong> 2.1 -75.85033256 10.489<br />

CCSD MRCC 2.1 -75.85033257 10.489<br />

CCSDT MRCC 2.1 -75.85316687 10.505<br />

CCSD <strong>MOLPRO</strong> 2.2 -75.95093334 7.752<br />

CCSD MRCC 2.2 -75.95093335 7.752<br />

CCSDT MRCC 2.2 -75.95299013 7.789

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

Saved successfully!

Ooh no, something went wrong!