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.

42 GEOMETRY OPTIMIZATION (OPTG) 330<br />

***,HF dimer MP2/CP optimization with relaxed monomers<br />

basis=avtz<br />

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

! INITIAL VALUES OF GEOMETRY VARIABLES<br />

RFF= 5.3<br />

R1= 1.76<br />

R2 = 1.75<br />

THETA1 = 7.0<br />

THETA2 = 111<br />

symmetry,x<br />

orient,noorient<br />

geometry={<br />

f1<br />

f2 f1 rff<br />

h1 f1 r1 f2 theta1<br />

h2 f2 r2 f1 theta2 h1 180.}<br />

label:<br />

text, CALCULATION AT LARGE SEPARATION<br />

rff_save=rff<br />

rff=1000<br />

!save current rff distance<br />

!dimer calculation at large separation<br />

text, HF1<br />

dummy,f2,h2;<br />

{hf;accu,16}<br />

mp2;<br />

ehf1inf=energy<br />

forces;<br />

!second hf is now dummy<br />

!scf for first monomer<br />

!mp2 for first monomer<br />

!save mp2 energy in variable<br />

!compute mp2 gradient for first monomer<br />

text, HF2<br />

dummy,f1,h1;<br />

{hf;accu,16}<br />

mp2;<br />

ehf2inf=energy<br />

forces;<br />

add,1<br />

einf=ehf1inf+ehf2inf<br />

rff=rff_save<br />

!first hf is now dummy<br />

!scf for second monomer<br />

!mp2 for second monomer<br />

!save mp2 energy in variable<br />

!compute mp2 gradient for second monomer<br />

!add to previous gradient<br />

!total energy of unrelaxed momomers<br />

!reset HF - HF distance to current value<br />

text, CP calculation for HF1 MONOMER<br />

dummy,f2,h2;<br />

{hf;accu,16}<br />

mp2;<br />

ehf1=energy<br />

forces;<br />

add,-1<br />

!second hf is now dummy<br />

!scf for first monomer<br />

!mp2 for first monomer<br />

!save mp2 energy in variable<br />

!compute mp2 gradient for first monomer<br />

!subtract from previous gradient<br />

text, CP calculation for HF2 MONOMER<br />

dummy,f1,h1;<br />

{hf;accu,16}<br />

mp2;<br />

ehf2=energy<br />

forces;<br />

add,-1<br />

!first hf is now dummy<br />

!scf for second monomer<br />

!mp2 for second monomer<br />

!save mp2 energy in variable<br />

!compute mp2 gradient for first monomer<br />

!subtract from previous gradient<br />

text, DIMER CALCULATION<br />

dummy<br />

!reset dummies<br />

{hf;accu,16}<br />

!scf for dimer<br />

mp2;<br />

!mp2 for dimer

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

Saved successfully!

Ooh no, something went wrong!