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.

42 GEOMETRY OPTIMIZATION (OPTG) 332<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 />

!second hf is now dummy<br />

!scf for first monomer<br />

!mp2 for first monomer<br />

!save mp2 energy in variable<br />

text, HF2<br />

dummy,f1,h1;<br />

{hf;accu,16}<br />

mp2;<br />

ehf2inf=energy<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 />

!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 />

!second hf is now dummy<br />

!scf for first monomer<br />

!mp2 for first monomer<br />

!save mp2 energy in variable<br />

text, CP calculation for HF2 MONOMER<br />

dummy,f1,h1;<br />

{hf;accu,16}<br />

mp2;<br />

ehf2=energy<br />

!first hf is now dummy<br />

!scf for second monomer<br />

!mp2 for second monomer<br />

!save mp2 energy in variable<br />

text, DIMER CALCULATION<br />

dummy<br />

!reset dummies<br />

{hf;accu,16}<br />

!scf for dimer<br />

mp2;<br />

!mp2 for dimer<br />

edimer=energy<br />

!save mp2 energy in variable<br />

etot=edimer-ehf2-ehf1+ehf1inf+ehf2inf !total BSSE corrected energy<br />

optg,numerical,variable=etot,gradient=1.d-4,startcmd=label:<br />

!optimize geometry<br />

text, compute optimized monomer energy

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

Saved successfully!

Ooh no, something went wrong!