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.

33 SYMMETRY-ADAPTED INTERMOLECULAR PERTURBATION THEORY 254<br />

33.7.3 DF-DFT-SAPT calculation of the NeAr dimer using the δ(HF) correction<br />

gdirect; gthresh,energy=1.d-8,orbital=1.d-8,grid=1.d-8<br />

symmetry,nosym<br />

orient,noorient<br />

geometry={<br />

Ne,,0.0,0.0,0.0<br />

Ar,,0.0,0.0,6.5}<br />

basis={<br />

set,orbital; default,avtz<br />

!for orbitals<br />

set,jkfit; default,avtz/jkfit !for JK integrals<br />

set,mp2fit; default,avtz/mp2fit !for E2disp/E2exch-disp<br />

set,dflhf; default,avtz/jkfit !for LHF<br />

}<br />

!=========delta(HF) contribution for higher order interaction terms====<br />

ca=2101.2; cb=2102.2 !sapt files<br />

!dimer<br />

{df-hf,basis=jkfit,locorb=0}<br />

edm=energy<br />

!monomer A<br />

dummy,ar<br />

{df-hf,basis=jkfit,locorb=0; save,$ca}<br />

ema=energy; sapt;monomerA<br />

!monomer B<br />

dummy,ne<br />

{df-hf,basis=jkfit,locorb=0; save,$cb}<br />

emb=energy; sapt;monomerB<br />

!interaction contributions<br />

{sapt,SAPT_LEVEL=2;intermol,ca=$ca,cb=$cb,icpks=1,fitlevel=3<br />

dfit,basis_coul=jkfit,basis_exch=jkfit,cfit_scf=3}<br />

!calculate high-order terms by subtracting 1st+2nd order energies<br />

eint_hf=(edm-ema-emb)*1000 mH<br />

delta_hf=eint_hf-e1pol-e1ex-e2ind-e2exind<br />

!=========DFT-SAPT at second order intermol. perturbation theory====<br />

ca=2103.2; cb=2104.2 !sapt files;<br />

!shifts for asymptotic correction to xc potential<br />

eps_homo_pbe0_ar=-0.440936 !HOMO(Ar)/PBE0 functional<br />

eps_homo_pbe0_ne=-0.589207 !HOMO(Ne)/PBE0<br />

ip_ar=0.5792<br />

!exp. ionisation potential<br />

ip_ne=0.7925<br />

!exp. ionisation potential<br />

shift_ar=ip_ar+eps_homo_pbe0_ar !shift for bulk xc potential (Ar)<br />

shift_ne=ip_ne+eps_homo_pbe0_ne !shift for bulk xc potential (Ne)<br />

!monomer A, perform LPBE0AC calculation<br />

dummy,ar<br />

{df-ks,pbex,pw91c,lhf; dftfac,0.75,1.0,0.25; asymp,shift_ne; save,$ca}<br />

sapt;monomerA<br />

!monomer B, perform LPBE0AC calculation<br />

dummy,ne<br />

{df-ks,pbex,pw91c,lhf; dftfac,0.75,1.0,0.25; start,atdens; asymp,shift_ar; save,$cb}<br />

sapt;monomerB<br />

!interaction contributions<br />

{sapt,SAPT_LEVEL=3;intermol,ca=$ca,cb=$cb,icpks=0,fitlevel=3,nlexfac=0.0<br />

dfit,basis_coul=jkfit,basis_exch=jkfit,cfit_scf=3}<br />

!add high-order approximation to obtain the total interaction energy<br />

eint_dftsapt=e12tot+delta_hf<br />

http:

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

Saved successfully!

Ooh no, something went wrong!