Transcript molpro

MOLPRO
A quantum chemistry package
Adem Tekin
Yrd. Doc. Dr.
15.06.2012 Istanbul
Outline
— Introduction
— Input structure and introductory examples
— Brief introduction to SCS-MP2 and DFT-SAPT
— Applications of SAPT: Acetylene-Benzene dimer
— Towards the AcnBzn aggregates
— Exercise: Calculation of interaction energies for
water dimer at MP2, SCS-MP2, B3LYP-D,
DFT-SAPT (PBE0AC), DFT-SAPT (LPBE0AC) and
CCSD(T) levels
Turbomole
15.06.2012
Molpro
— Developed by H. –J. Werner (University of Stuttgart) and P. J. Knowles
(Cardiff University)
— for more information please visit to:
— http://www.molpro.net
— A user mailing list is available under
http://www.molpro.net/mailman/listinfo/molpro-user/?portal=
user&choice=User+mailing+list
— In molpro, the emphasis is given to highly accurate computations,
with extensive treatment of the electron correlation problem through
the multiconfiguration-reference CI and coupled cluster methods.
Turbomole
15.06.2012
Running molpro
— Molpro is accessed using the molpro command. The syntax is:
molpro [options] [datafile]
— datafile is the input. Prepare your input with .com extension.
— output will be written to a datafile.out file
— -d dir1:dir2:…
where dir1:dir2:… is a list of directories which may be used
for creating scratch files.
— There are more options…
— Molpro can be run in parallel:
-n specifies the number of cpu
-N specifies the number of tasks on the node
Turbomole
15.06.2012
Submitting molpro jobs on UHEM
#!/bin/bash
#BSUB -m karadeniz_temp1
#BSUB -a intelmpi
#BSUB -J triazine_avtz16
#BSUB -o triazine_avtz16.out
#BSUB -e triazine_avtz16.err
#BSUB -q workshop
#BSUB -n 16
#BSUB –P workshop
#BSUB -R "span[ptile=8]"
# bu kismi degistirmeyin !!!
#isinizi tanimlayacak bir isim
# bu kismi degistirmeyin !!!
# bu kismi degistirmeyin !!!
#kuyruk ismi
# kullanilacak olan islemci sayisi
# bu kismi degistirmeyin !!!
echo 'Starting time:'
date
molpro -n 16 -N 8 --no-xml-output
$PWD/triazine_avtz16.com>triazine_avtz16.log
echo 'Ending time:'
date
Turbomole
15.06.2012
Molpro input structure
***,title
memory,4,m
file,1,name.int
file,2,name.wfu
gprint,options
gthresh,options
gdirect[,options]
gexpec,opnames
!title (optional)
!memory specification (optional)
!permanent named integral file (optional)
!permanent named wavefunction file (optional)
!global print options (optional)
!global thresholds (optional)
!global direct (optional)
!global definition of one-electron operators
(optional)
basis=basisname
!basis specification. If not present, cc-pVDZ
!is default
geometry={...}
!geometry specification
var1=value,var2=value,... !setting variables for geometry and/or
wavefunction
{command,options
!program or procedure name
directive,data,option !directives for command (optional)
...
}
!end of command block
--!end of input (optional)
Turbomole
15.06.2012
Molpro input examples
(http://www.be.itu.edu.tr/~adem.tekin/molpro)
! SCF calcultion for water. (ex1.com)
TURBOMOLE usage philosophy
***,h2o
!A title
r=1.85,theta=104 !set geometry parameters
geometry={O;
!z-matrix geometry input
H1,O,r;
H2,O,r,H1,theta}
hf
!closed-shell scf
The default basis is cc-pVDZ abbreviated as VDZ
! SCF calculation for water using vtz (ex2.com)
***,h2o cc-pVTZ basis
!A title
r=1.85,theta=104
!set geometry parameters
geometry={O;
!z-matrix geometry input
H1,O,r;
H2,O,r,H1,theta}
basis=VTZ
!use VTZ basis
hf
!closed-shell scf
Molpro input examples
(http://www.be.itu.edu.tr/~adem.tekin/molpro)
:%s/water-ccsdt-avdz/ex1/g
Turbomole
15.06.2012
Molpro examples
! Geometry optimization at the HF level for water (ex3.com)
***,h2o
!A title
r=1.85,theta=104
!set geometry parameters
geometry={O;
!z-matrix geometry input
H1,O,r;
H2,O,r,H1,theta}
basis=6-31g**
!use Pople basis set
hf
!closed-shell scf
optg
!do scf geometry optimization
! Single point CCSD(T) for water (ex4.com)
Turbomole
***,h2o
!A title
r=1.85,theta=104
!set geometry parameters
geometry={O;
!z-matrix geometry input
H1,O,r;
H2,O,r,H1,theta}
basis=VTZ
!use VTZ basis
hf
!closed-shell
scf
15.06.2012
ccsd(t)
!do ccsd(t) calculation
Electron correlation methods
— Hartree-Fock is a single determinant wave function method where molecular orbitals are varied.
— Configuration Interaction (CI) is a multiple determinant wavefunction method where molecular
orbitals are not varied.
— CI wavefunction is constructed by starting with the HF wavefunction and making new determinants by
promoting electrons from the occupied to unoccupied orbitals.
— Depending on the number of excitations used to make the determinants, CI methods are called CIS,
CISD, CISDT, CISDTQ and full CI.
— Multi-configurational self-consistent field (MCSCF) calculations also use multiple determinants. In
contrast to the CI, orbitals are also optimized.
— An MCSCF calculation in which all combinations of the active space orbitals are included is called a
complete active space self-consistent field (CASSCF) calculation.
Turbomole
15.06.2012
Electron correlation methods
— In a CASSCF wavefunction the occupied orbital space is divided into a set of inactive or closed-shell
orbitals and a set of active orbitals. All inactive orbitals are doubly occupied in each Slater determinant.
On the other hand, the active orbitals have varying occupations, and all possible Slater determinants
(or CSFs) are taken into account which can be generated by distributing the Nact=Nel-2mclosed electrons
in all possible ways among the active orbitals, where mclosed is the number of closed-shell (inactive)
orbitals, and Nel is the total number of electrons.
— It is possible to construct a CI wavefunction starting with an MCSCF calculation rather than starting with
a HF wavefunction. This starting wavefunction is called the reference state. These calculations are
called multi-reference configuration interaction (MRCI).
— Coupled cluster (CC) calculations are similar to CI calculations in that the wavefunction is a linear
combination of many determinants. The means for choosing the determinants in CC is more complex than
CI. There are variants of CC: CCSD, CCSDT and CCSD(T).
Turbomole
15.06.2012
Molpro examples
! Geometry optimization at the HF level for water (ex5.com)
***,h2o
!A title
r=1.85,theta=104
!set geometry parameters
geometry={o;
!z-matrix geometry input
h1,O,r;
h2,O,r,H1,theta}
basis=vtz
!use VTZ basis
hf
!closed-shell scf
ccsd(t)
!do ccsd(t) calculation
casscf
!do casscf calculation
mrci
!do mrci calculation
Turbomole
15.06.2012
Molpro example – tables
! Put the results into a table (ex6.com)
***,h2o
!A title
r=1.85,theta=104
!set geometry parameters
geometry={o;
!z-matrix geometry input
h1,O,r;
h2,O,r,H1,theta}
basis=vtz
!use VTZ basis
hf
!closed-shell scf
e(1)=energy
!save scf energy in variable e(1)
method(1)=program
!save the string ’HF’ in variable method(1)
ccsd(t)
!do ccsd(t) calculation
e(2)=energy
!save ccsd(t) energy in variable e(2)
method(2)=program
!save the string ’CCSD(T)’ in variable method(2)
casscf
!do casscf calculation
e(3)=energy
!save scf energy in variable e(3)
method(3)=program
!save the string ’CASSCF’ in variable method(3)
mrci
!do mrci calculation
e(4)=energy
!save scf energy in variable e(4)
method(4)=program
!save the string ’MRCI’ in variable method(4)
table,method,e
!print a table with results
title,Results for H2O, basis=$basis
!title for the table
Molpro example – Procedures (ex7.com)
***,h2o
proc save_e
if(#i.eq.0) i=0
i=i+1
e(i)=energy
method(i)=program
endproc
!A title
!define procedure save_e
!initialize variable i if it does not exist
!increment i
!save scf energy in variable e(i)
!save the present method in variable method(i)
!end of procedure
r=1.85,theta=104
!set geometry parameters
geometry={o;
!z-matrix geometry input
h1,O,r;
h2,O,r,H1,theta}
basis=vtz
!use VTZ basis
hf
!closed-shell scf
save_e
!call procedure, save results
ccsd(t)
!do ccsd(t) calculation
save_e
!call procedure, save results
casscf
!do casscf calculation
save_e
!call procedure, save results
mrci
!do mrci calculation
save_e
!call procedure, save results
table,method,e
!print a table with results
title,Results for H2O, basis=$basis
!title for the table
Molpro example – do loops (ex8.com)
***,H2O potential
symmetry,x
!use cs symmetry
geometry={
o;
!z-matrix
h1,o,r1(i);
h2,o,r2(i),h1,theta(i) }
basis=vdz
!define basis set
angles=[100,104,110]
!list of angles
distances=[1.6,1.7,1.8,1.9,2.0]
!list of distances
i=0
!initialize a counter
do ith=1,#angles
!loop over all angles H1-O-H2
do ir1=1,#distances
!loop over distances for O-H1
do ir2=1,ir1
!loop over O-H2 distances(r1.ge.r2)
i=i+1
!increment counter
r1(i)=distances(ir1)
!save r1 for this geometry
r2(i)=distances(ir2)
!save r2 for this geometry
theta(i)=angles(ith)
!save theta for this geometry
hf;
!do SCF calculation
escf(i)=energy
!save scf energy for this geometry
ccsd(t);
!do CCSD(T) calculation
eccsd(i)=energc
!save CCSD energy
eccsdt(i)=energy
!save CCSD(T) energy
enddo
!end of do loop ith
enddo
!end of do loop ir1
enddo
!end of do loop ir2
{table,r1,r2,theta,escf,eccsd,eccsdt
!produce a table with results
head, r1,r2,theta,scf,ccsd,ccsd(t)
!modify column headers for table
save,h2o.tab
!save the table in file h2o.tab
title,Results for H2O, basis $basis
!title for table
sort,3,1,2}
!sort table
Program control
— ***,text
! Starting a job
— ---
! Ending a job
— MEMORY,n,scale;
! Allocating dynamic memory
—
—
—
—
! allocates 90 000 words of memory
! allocates 500 000 words of memory
!allocates 2 000 000 words of memory
MEMORY,90000
MEMORY,500,K
MEMORY,2,M
1 word = 8 bytes
— IF($method.eq.’HF’) then
...
ENDIF
— GTHRESH,key1=value1,key2=value2,. . . ! Global thresholds
— ENERGY
— ORBITAL
! Convergence threshold for energy (default 1.d-6)
! Convergence threshold for orbital optimization in
the SCF program (default 1.d-5).
— gthresh,energy=1d-8,orbital=1d-7
Turbomole
15.06.2012
Molpro example – Expectation values – dipole and
quadrupole moments (ex9.com)
***,h2o properties
geometry={o;h1,o,r;h2,o,r,h1,theta}
!Z-matrix geometry input
r=1 ang
!bond length
theta=104
!bond angle
gexpec,dm,sm,qm
compute dipole and quarupole moments
$methods=[hf,multi,ci]
!do hf, casscf, mrci
do i=1,#methods
!loop over methods
$methods(i)
!run energy calculation
e(i)=energy
dip(i)=dmz
!save dipole moment in variable dip
quadxx(i)=qmxx
!save quadrupole momemts
quadyy(i)=qmyy
quadzz(i)=qmzz
smxx(i)=xx
!save second momemts
smyy(i)=yy
smzz(i)=zz
enddo
table,methods,dip,smxx,smyy,smzz
!print table of first and
second moments
table,methods,e,quadxx,quadyy,quadzz !print table of quadrupole
moments
Geometry specification – using xyz
coordinates
geomtyp=xyz
geometry={
3 ! number of atoms
This is an example of geometry input for water
with an XYZ file
O ,0.0000000000,0.0000000000,-0.1302052882
H ,1.4891244004,0.0000000000, 1.0332262019
H,-1.4891244004,0.0000000000, 1.0332262019
}
hf
Turbomole
15.06.2012
Basis set
— BASIS,DEFAULT=VTZ,O=AVTZ,H=VDZ
— BASIS=VTZ,O=AVTZ,H=VDZ
— BASIS
DEFAULT=VTZ
O=AVTZ
H=VDZ
END
Turbomole
15.06.2012
Density fitting
— Density fitting can be used to approximate the integrals in spin restricted
Hartree-Fock (HF), density functional theory (KS), second-order Møller-Plesset
perturbation theory (MP2 and RMP2), explicitly correlated MP2 (MP2-F12), and all levels
of closed-shell local correlation methods (LCC2, LMP2-LMP4, LQCISD(T), LCCSD(T)).
— Density fitting is invoked by adding the prefix DF- to the command name, e.g. DF-HF,
DF-KS, DF-MP2 and so on.
— Symmetry is not implemented for density fitting programs.
— By default, a fitting basis set will be chosen automatically that corresponds to the current
orbital basis set and is appropriate for the method. For instance, if the orbital basis set is
VTZ, the default fitting basis is VTZ/JKFIT for DF-HF or DF-KS, and VTZ/MP2FIT for
DF-MP2.
— Examples:
BASIS=VTZ
DF-HF,DF_BASIS=VQZ
DF-MP2,DF_BASIS=VQZ
Turbomole
!use VTZ orbital basis
!use VQZ/JKFIT fitting basis
!use VQZ/MP2FIT fitting basis
15.06.2012
Calculating the intermolecular interactions
in water dimer
— MP2
— SCS-MP2
— B3LYP-D
— DFT-SAPT
— CCSD(T)
Turbomole
15.06.2012
Spin component scaled (SCS) MP2
— SCS-MP2 is based on a separate scaling of the second-order parallel
(αα+ββ) and antiparallel-spin (αβ) pair correlation energies e(2).
— The SCS-MP2 approximation to the correlation energy Ecorr is given by
— Here, pT and pS are empirical parameters determined to be 1/3 and 6/5.
— If you perform an MP2 calculation, you already have the SCS-MP2 energy.
Turbomole
S. Grimme, J. Chem. Phys., 2003, 118, 9095.
The iterative formulation of RayleighSchrödinger perturbation theory
Turbomole
Intermolecular perturbation theory
,
Turbomole
,
Symmetry adapted perturbation theory
• The total interaction energy :
• Contributions of third and higher orders :
Turbomole
15.06.2012
A more simple way to calculate
intermolecular interaction energy
• Supermolecular approach (SMA)
• Basis set superposition error (BSSE)
• Employ Counterpoise (CP) technique
to avoid BSSE in SMA
• SAPT does not contain BSSE
Turbomole
15.06.2012
The analysis of different types of
intermolecular interactions
Turbomole
15.06.2012
Applications of SAPT: Acetylene - Benzene
T
S
-2.9
-5.8
-10.8
SS
Turbomole
15.06.2012
Comparison of SAPT with supermolecular
calculations: Acetylene - Benzene
Turbomole
Aug-cc-pVTZ
Comparison of SAPT with supermolecular
calculations: Acetylene - Furan
Turbomole
15.06.2012
Aug-cc-pVTZ
The intermolecular PES: extrapolation to
the cbs limit for Acetylene - Benzene
E
corr
X 3 E corr ( X )  ( X  1)3 E corr ( X  1)
( X  1, X ) 
with X  ( D : 2, T : 3, Q : 4...)
3
3
X  ( X  1)
• 42 test geometries (including
global minimum region)
• Only
extrapolated
• Extrapolation with doubleand triple- zeta basis sets
sufficient here
formula of Halkier et al. [CPL 286 (1998) 243]
Turbomole
15.06.2012
Intermolecular potential energy surface of the
Acetylene-Benzene dimer: Choice of geometries
Grid Definition
3.0 ≤ R ≤ 7.0 Å
R  ↑ by 0.5 Å
0 ≤ Θ ≤ 90°
0 ≤ Φ ≤ 30°
0 ≤ θ ≤ 90°
0 ≤ φ ≤ 180°
All angles ↑
by 30 °
Resulting Grid
8 × 4 × 2 × 4 × 7 = 1792
1792 – 1099 = 693
C. A. Hunter et al. J. Am. Chem. Soc. 112 (1990), 5525
J. Gauss et al. J. Phys. Chem. 104 (2000),2865
J. L. M. Martin et al. J. Chem. Phys. 108 (1998), 676
Asymptotic interactions: Electrostatic energy
via multipole expansion and Dispersion
energy
Quadrupole-Quadrupole
interaction
Quadrupole-Octapole
interaction
Octapole-Octapole
interaction
•Quadrupole moments and Dispersion coefficients are calculated only once.
•Rotated for other orientations using Euler`s rotation matrices.
•10.0 ≤ R ≤ 49.0 Å, R increased by 3 Å
•10192 orientations produced in this manner.
•10192 + 693 = 10885 (Total number of points used in fit procedure)
Turbomole
15.06.2012
Fitting the intermolecular PES: Fit of
extrapolated total interaction energy
• Site-site fit model:
Repulsion
Dispersion
&
Induction
Electrostatics
• Damping function:
• Merit function (Chi-Square):
C. Leforestier, A. Tekin, G. Jansen , M. Herman, J. Chem. Phys., 135, 234306, (2011)
Fitting the intermolecular PES: Fit quality I
For 10885 points SD [kJ/mol]: 0.92
-唴
0.04
6.9 (15.5 %)
2.6
Turbomole
For 693 points SD [kJ/mol]: 3.63
-唴
唴
Eint [kJ/mol]
15.06.2012
0.18
6.9 (15.5 %)
2.6
唴
Eint [kJ/mol]
Fitting the intermolecular PES: Fit quality II
Turbomole
15.06.2012
Fitting the intermolecular PES: Fit quality III
Turbomole
R (Å)
4.3
5.7
Θ1
15°
80°
Φ1
15°
Θ2
20°
50°
Φ2
25°
145°
15.06.2012
8.5
Attention!!
Turbomole
15.06.2012
Fitting the intermolecular PES: Comparison
to Experiment
(SS1)
Turbomole
They “confirmed” the existence of 1B employing
HF, MP2 and B3LYP with 6-31G++(d,p) geometry
optimizations.
15.06.2012
Intermolecular potential energy surface of
the Acetylene - Acetylene dimer: Choice
of geometries
Grid Definition
3.0 ≤ R ≤ 7.0 Å
R  ↑ by 0.5 Å
0 ≤ Θ1 ≤ 90°
0 ≤ Θ2 ≤ 180°
0 ≤ Φ2 ≤ 360°
Θ1 ↑ by 15 °
Θ2 and Φ2 ↑ by 30 °
Resulting Grid
9 × 7 × 7 × 13 = 5733
5733 – 4904 = 829
Turbomole
15.06.2012
C. Leforestier, A. Tekin, G. Jansen , M. Herman, J. Chem. Phys., 135, 234306, (2011)
Intermolecular potential energy surface of
the Benzene - Benzene dimer: Choice of
geometries
Grid Definition
3.0 ≤ R ≤ 7.0 Å
R  ↑ by 1 Å
0 ≤ Θ ≤ 90°
0 ≤ Φ ≤ 30°
0 ≤ φ ≤ 360°
0 ≤ θ ≤ 360°
0 ≤ ω ≤ 360°
All angles ↑
by 30 °
Resulting Grid
5 × 4 × 2 × 13 × 13 × 13 = 87880
87880 – 86532 = 1348
Turbomole
15.06.2012
AcnBz cluster geometry optimizations:
Test system Ac2Bz
D
Ring
Turbomole
6h
C
2V
(kcal/mol)
Ring
D6h
C2V
MP2/VDZa
0.106
0.
1.517
MP2/TZVP
0.
0.580 2.165
MP2/VDZ/CP
0.
0.241 1.297
Model
0.
0.241 1.401
15.06.2012
a Fujii et al. J. Phys. Chem. A 108 (2004) 2652.
Structural comparison of Ac2Bz
geometries
B
C
A
A
D
C
E
G
B
F
H
(Å)
AB
AC
DE
DF
DG
DH
MP2/TZVP
2.81 2.61 2.66 2.66 2.69 2.68
MP2/VDZ/CP
2.90 2.82 2.82 2.82 2.91 2.91
Model
2.85 2.72 2.82 2.82 2.85 2.85
Turbomole
15.06.2012
(Å)
AB
AC
MP2/TZVP
4.69
2.73
MP2/VDZ/CP 5.15
2.93
Model
4.96
2.85
Acn (n=3-6) clusters
-10.79
-7.11
-10.73
-14.46
-19.30
-14.82
Turbomole
15.06.2012
(Energies are in mH)
AcnBz (n=3-6) clusters
-13.49
-13.51
-19.11
-24.43
-29.99
-30.52
Turbomole
15.06.2012
(Energies are in mH)
Ac1Bz2 clusters
-10.53
-8.74
A
Angew. Chem. Int. Ed. 47 (2008) 10094
Turbomole
-8.31
(Energies
15.06.2012 are in mH)
Ac2Bz2 clusters
-15.49
-16.25
-15.39
-15.37
-15.04
-15.23
-13.36
Turbomole
-12.97
15.06.2012
-14.97
(Energies are in mH)
Intermolecular interactions [in kJ/mol] in
water dimer
Basis
HF
MP2
SCS-MP2
B3LYP-D
PBE0
LPBE0
CCSD(T)
AVDZ
-15.891390
-18.803673
-17.166490
-21.639948
-17.386616
-18.428262
-18.660504
AVTZ
-15.765523
-19.932165
-18.160662
-21.843949
-18.707867
-19.887984
-20.067011
AVQZ
-15.908823
-20.559949
-18.804303
-22.006363
-19.133093
-20.337871
-20.733232
— Molpro energy results are in hartree. Convert them to kJ/mol by multiplying 2625.5.
— DFT-SAPT results are in mH (mili Hartree). Multiply them by 2.6255.
— HF, MP2, SCS-MP2, CCSD and CCSD(T) results are in CCSD(T) outputs.
— Only the SCS-MP2 interaction energy is not written in the output. To calculate it grep
SCS-MP2 energies in the ccsd(t) outputs:
— grep "SCS-MP2 total energy" water-ccsdt-avdz.out
Turbomole
15.06.2012
ESCS-MP2=(-152.52056354-(-76.25674577-76.25727940)) in hartree
Typical CCSD(T) output
!RHF STATE 1.1 Energy
MP2 total energy
SCS-MP2 total energy
!CCSD(T) total energy
SETTING
SETTING
SETTING
SETTING
SETTING
SETTING
SETTING
Turbomole
E_MB_CP_MP2
E_MB_CP_CCSD
E_MB_CP_CCSDT
E_INT_HF
E_INT_MP2
E_INT_CCSD
E_INT_CCSDT
=
=
=
=
=
=
=
15.06.2012
-76.26175127
-76.26944486
-76.27482665
-0.00605271
-0.00716194
-0.00685369
-0.00710741
AU
AU
AU
AU
AU
AU
AU
Typical CCSD(T) output
cp water-ccsdt-avdz.com water-ccsdt-avtz.com
cp water-ccsdt-avtz.com water-ccsdt-avqz.com
Turbomole
15.06.2012
Typical DFT-SAPT output
SETTING DHF
=
-1.15745486
mH
===========
IMW Results
===========
E1pol
E1exch
E1exch(S2)
E2ind
E2ind-exch
E2disp(unc)
E2disp
E2disp-exch(unc)
E2disp-exch
E1tot
E2tot
E1tot+E2tot
[mH]
-11.45483815
10.21741485
10.15101824
-3.98070553
2.32414481
-4.50262689
-3.17861613
0.79013747
0.60784335
-1.23742331
-4.22733350
-5.46475681
SETTING ETOTAL
Turbomole
(
(
(
(
(
(
(
(
(
-0.11454838E+02)
0.10217415E+02)
0.10151018E+02)
-0.39807055E+01)
0.23241448E+01)
-0.45026269E+01)
-0.31786161E+01)
0.79013747E+00)
0.60784335E+00)
( -0.12374233E+01)
( -0.42273335E+01)
( -0.54647568E+01)
=
[kcal/mol]
-7.1880
6.4115
6.3699
-2.4979
1.4584
-2.8254
-1.9946
0.4958
0.3814
-0.7765
-2.6527
-3.4292
-6.62221167
15.06.2012
mH
[kJ/mol]
-30.0747
26.8258
26.6515
-10.4513
6.1020
-11.8216
-8.3455
2.0745
1.5959
-3.2489
-11.0989
-14.3477
I. Supermolecular CCSD(T) input
(http://www.be.itu.edu.tr/~adem.tekin/mol
pro/water-ccsdt-avdz.com)
memory,550,m
gthresh,energy=1d-8,orbital=1d-7
basis,default=avdz
!nosym,
! IF THERE IS NO SYMMETRY
symmetry,z
angstrom
GEOMTERY={
O1,,
0.1134349,
-1.6972642,
0.0000000
H1,,
-0.0489012,
-0.7388274,
0.0000000
H2,,
-0.7694849,
-2.0866817,
0.0000000
O2,,
-0.0837066,
1.2726931,
0.0000000
H3,,
0.3943289,
1.6250401,
0.7627674
H4,,
0.3943289,
1.6250401,
-0.7627674
} Turbomole
15.06.2012
int
II. Supermolecular CCSD(T) input
! --- dimer --hf
e_dim_hf=energy
ccsd(t)
e_dim_mp2=emp2
e_dim_ccsd=energc
e_dim_ccsdt=energt(1)
Turbomole
! --- monomerA-CP --dummy,O2,H3,H4
hf
e_mA_CP_hf=energy
ccsd(t)
e_mA_CP_mp2=emp2
e_mA_CP_ccsd=energc
e_mA_CP_ccsdt=energt(1)
III. Supermolecular CCSD(T) input
! --- monomerB-CP --dummy,O1,H1,H2
hf
e_mB_CP_hf=energy
ccsd(t)
e_mB_CP_mp2=emp2
e_mB_CP_ccsd=energc
e_mB_CP_ccsdt=energt(1)
e_int_hf = e_dim_hf - e_mA_CP_hf - e_mB_CP_hf
e_int_mp2 = e_dim_mp2 - e_mA_CP_mp2 - e_mB_CP_mp2
e_int_ccsd=e_dim_ccsd-e_mA_CP_ccsd-e_mB_CP_ccsd
e_int_ccsdt=e_dim_ccsdt-e_mA_CP_ccsdt-e_mB_CP_ccsdt
Turbomole
15.06.2012
I. Density Fitting MP2 input
(http://www.be.itu.edu.tr/~adem.tekin/molpro/water-mp2-avdz.com)
memory,200,m
gdirect
gthresh,energy=1d-8,orbital=1d-7
basis={
default,avdz
set,mp2fit
default,avdz/MP2FIT
set,jkfit
default,vtz/jkfit
}
nosym,
angstrom
GEOMTERY={
O1,,
0.1134349,
-1.6972642,
H1,,
-0.0489012,
-0.7388274,
H2,,
-0.7694849,
-2.0866817,
O2,,
-0.0837066,
1.2726931,
H3,,
0.3943289,
1.6250401,
H4,,
0.3943289,
1.6250401,
}
int
0.0000000
0.0000000
0.0000000
0.0000000
0.7627674
-0.7627674
II. Density Fitting MP2 input
! --- dimer --df-hf,basis=jkfit,locfit=0
e_dim_hf=energy
df-mp2,basis=mp2fit,locfit=0
e_dim_mp2=emp2
! --- monomerA-CP --dummy,O2,H3,H4
df-hf,basis=jkfit,locfit=0
e_mA_CP_hf=energy
df-mp2,basis=mp2fit,locfit=0
e_mA_CP_mp2=emp2
! --- monomerB-CP --dummy,O1,H1,H2
df-hf,basis=jkfit,locfit=0
e_mB_CP_hf=energy
df-mp2,basis=mp2fit,locfit=0
e_mB_CP_mp2=emp2
e_int_hf = e_dim_hf - e_mA_CP_hf - e_mB_CP_hf
e_int_mp2 = e_dim_mp2 - e_mA_CP_mp2 - e_mB_CP_mp2
Basis set selection for water-b3lyp-davdz.com
for AVDZ
for AVTZ
for AVQZ
basis={
default,avdz
set,jkfit
default,vtz/jkfit
set,mp2fit
default,avdz/mp2fit}
basis={
default,avtz
set,jkfit
default,vqz/jkfit
set,mp2fit
default,avtz/mp2fit}
basis={
default,avqz
set,jkfit
default,v5z/jkfit
set,mp2fit
default,avqz/mp2fit}
Turbomole
15.06.2012
I. DFT-SAPT (PBE0AC) input
(http://www.be.itu.edu.tr/~adem.tekin/molpro/water-pbe0-avdz.com)
memory,170,m
gdirect
gthresh,energy=1d-8,orbital=1d-8
symmetry,nosym
orient,noorient
GEOMTERY={angstrom;
! water
1,O1,,
0.1134349,
-1.6972642,
2,H1,,
-0.0489012,
-0.7388274,
3,H2,,
-0.7694849,
-2.0866817,
4,O2,,
-0.0837066,
1.2726931,
5,H3,,
0.3943289,
1.6250401,
6,H4,,
0.3943289,
1.6250401,
}
basis={
default,avdz
set,jkfit
default,vtz/jkfit
set,mp2fit
default,avdz/mp2fit}
int
0.0000000
0.0000000
0.0000000
0.0000000
0.7627674
-0.7627674
II. DFT-SAPT (PBE0AC) input
ca=2101.2
cb=2102.2
!--HF------------df-hf,basis=jkfit,locfit=0
edm=energy
dummy,O2,H3,H4
df-hf,basis=jkfit,locfit=0
save,$ca
ema=energy
sapt;monomerA
dummy,O1,H1,H2
{df-hf,basis=jkfit,locfit=0; start,atdens
save,$cb}
emb=energy
sapt;monomerB
Turbomole
15.06.2012
III. DFT-SAPT (PBE0AC) input
{grid;gridthr,1d-6}
{sapt,SAPT_NFRQ_DISP=8;intermol,saptlevel=2,ca=$ca,cb=$cb,icp
ks=1,fitlevel=3
dfit,basis_coul=jkfit,basis_exch=jkfit,basis_mp2=mp2fit,cfit_
scf=3}
dHF=(edm-ema-emb)*1000.-e1pol-e1ex-e2ind-e2exind
!--DFT-----------grid;gridthr,1d-8
ca=2103.2
cb=2104.2
dummy,O2,H3,H4
df-ks,pbe0
asymp,0.13,0.5,40,0.05
dfit,basis=jkfit,locfit=0
save,$ca
sapt;monomerA
Turbomole
! 0.13= Ionization energy - HOMO
15.06.2012
IV. DFT-SAPT (PBE0AC) input
dummy,O1,H1,H2
dfit,basis=jkfit,locfit=0
{df-ks,pbe0; asymp,0.13,0.5,40,0.05; start,atdens
save,$cb}
sapt;monomerB
grid;gridthr,1d-6
{sapt,SAPT_LEVEL=3,SAPT_NFRQ_DISP=8;intermol,ca=$ca,cb=$cb,ic
pks=1,fitlevel=3,nlexfac=0.0
dfit,basis_coul=jkfit,basis_exch=jkfit,basis_mp2=mp2fit,cfit_
scf=3}
Etotal=e1pol+e1ex+e2ind+e2exind+e2disp+e2exdisp+dHF
Turbomole
15.06.2012
I. DFT-SAPT (LPBE0AC) input
(http://www.be.itu.edu.tr/~adem.tekin/molpro/water-lpbe0-avdz.com)
memory,300,m
gdirect
gthresh,energy=1d-8,orbital=1d-8
symmetry,nosym
orient,noorient
GEOMTERY={angstrom;
! water
1,O1,,
0.1134349,
-1.6972642,
2,H1,,
-0.0489012,
-0.7388274,
3,H2,,
-0.7694849,
-2.0866817,
4,O2,,
-0.0837066,
1.2726931,
5,H3,,
0.3943289,
1.6250401,
6,H4,,
0.3943289,
1.6250401,
}
basis={
set,orbital; default,avdz
set,jkfit; default,avdz/jkfit
set,mp2fit; default,avdz/mp2fit
set,dflhf; default,avdz/jkfit
}
0.0000000
0.0000000
0.0000000
0.0000000
0.7627674
-0.7627674
II. DFT-SAPT (LPBE0AC) input
ca=2101.2
cb=2102.2
!--HF------------df-hf,basis=jkfit,locorb=0
edm=energy
dummy,O2,H3,H4
df-hf,basis=jkfit,locorb=0
save,$ca
ema=energy
sapt;monomerA
dummy,O1,H1,H2
{df-hf,basis=jkfit,locorb=0
save,$cb}
emb=energy
sapt;monomerB
Turbomole
15.06.2012
III. DFT-SAPT (LPBE0AC) input
{grid;gridthr,1d-6}
{sapt,SAPT_NFRQ_DISP=8;intermol,saptlevel=2,ca=$ca,cb=$cb,icp
ks=1,fitlevel=3
dfit,basis_coul=jkfit,basis_exch=jkfit,,basis_mp2=mp2fit,cfit
_scf=3}
dHF=(edm-ema-emb)*1000.-e1pol-e1ex-e2ind-e2exind
!--DFT-----------grid;gridthr,1d-8
ca=2103.2
cb=2104.2
shift_pyr1=0.13
shift_pyr2=0.13
!monomer A, perform LPBE0AC calculation
dummy,O2,H3,H4
{df-ks,pbex,pw91c,lhf; dftfac,0.75,1.0,0.25;
asymp,shift_pyr1; save,$ca}
sapt;monomerA
IV. DFT-SAPT (LPBE0AC) input
!monomer B, perform LPBE0AC calculation
dummy,O1,H1,H2
{df-ks,pbex,pw91c,lhf; dftfac,0.75,1.0,0.25; start,atdens;
asymp,shift_pyr2; save,$cb}
sapt;monomerB
grid;gridthr,1d-6
{sapt,SAPT_LEVEL=3,SAPT_NFRQ_DISP=8;intermol,ca=$ca,cb=$cb,ic
pks=1,fitlevel=3,nlexfac=0.0
dfit,basis_coul=jkfit,basis_exch=jkfit,basis_mp2=mp2fit,cfit_
scf=3}
Etotal=e1pol+e1ex+e2ind+e2exind+e2disp+e2exdisp+dHF
Turbomole
15.06.2012
DFT-SAPT energy components
Component
AVDZ
AVTZ
DHF
E1pol
E1exch
E1exch(S2)
E2ind
E2ind-exch
E2disp(unc)
E2disp
E2disp-exch(unc)
E2disp-exch
-3.0389
-30.0747
26.8258
26.6515
-10.4513
6.1020
-11.8216
-8.3455
2.0745
1.5959
-3.0862
-30.0541
26.8072
26.6339
-10.9696
6.5146
-13.3978
-9.8115
2.3636
1.8918
-3.1033
-30.1058
26.7828
26.6098
-10.9495
6.4726
-13.8237
-10.2466
2.4813
2.0167
E1tot
E2tot
E1tot+E2tot
-3.2489
-11.0989
-14.3477
-3.2469
-12.3748
-15.6217
-3.3230
-12.7068
-16.0298
Eint
-17.3866
-18.7079
-19.1330
Turbomole
15.06.2012
AVQZ