19.01.2015 Views

MOLPRO

MOLPRO

MOLPRO

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

30 LOCAL METHODS FOR EXCITED STATES 227<br />

The densities saved with DENSAVE can be written into a cube file. Corresponding commands<br />

in CUBE module 34.7 are<br />

DENSITY,rec.ifil,cc2<br />

Creates a cube file from the ground state density saved on record<br />

rec.ifil.<br />

DENSITY,rec.ifil,full,state=st.1 Creates a cube file from the excited state density of state<br />

st.1.<br />

DENSITY,rec.ifil,diff,state=st.1 Creates a cube file from the difference density (full -<br />

cc2) for state st.1.<br />

30.5 Examples<br />

Example 1:<br />

***,The five lowest singlet excited states of water molecule<br />

memory,64,m<br />

gdirect<br />

symmetry,nosym;orient,noorient<br />

geometry={<br />

o ,, 0.000, 0.000, 0.119<br />

h1,, 1.423, 0.000, -0.947<br />

h2,, -1.423, 0.000, -0.947 }<br />

basis={<br />

default,vdz<br />

set,mp2fit<br />

default,vdz/mp2fit<br />

set,jkfit<br />

default,vdz/jkfit }<br />

hf<br />

{lt-df-lcc2<br />

!ground state CC2<br />

eom,-6.1<br />

!five lowest states<br />

eomprint,popul=-1,loceom=-1 } !minimize the output<br />

The excitation energies (in eV) stand after the calculation in array OMEGAF S for singlet states<br />

and in OMEGAF T for triplet states.<br />

Example 2:<br />

***,The five lowest triplet excited states and properties<br />

memory,64,m<br />

gdirect<br />

symmetry,nosym;orient,noorient<br />

geometry={<br />

o ,, 0.000, 0.000, 0.119<br />

h1,, 1.423, 0.000, -0.947<br />

h2,, -1.423, 0.000, -0.947 }<br />

basis={<br />

default,vdz<br />

set,mp2fit<br />

default,vdz/mp2fit<br />

set,jkfit

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

Saved successfully!

Ooh no, something went wrong!