turbomole

Report
TURBOMOLE
Program Package for ab initio Electronic
Structure Calculations
Adem Tekin
Yrd. Doc. Dr.
15.06.2012 Istanbul
Attention!!
Input and lsf scripts can be downloaded from:
http://www.be.itu.edu.tr/~adem.tekin/turbomole/filelist.txt
http://www.be.itu.edu.tr/~adem.tekin/turbomole/water.xyz
http://www.be.itu.edu.tr/~adem.tekin/molpro/filelist.txt
http://www.be.itu.edu.tr/~adem.tekin/molpro/ex1.com
Turbomole
15.06.2012
Activate turbomole and molpro in your shell
— Go to your home directory
— Open your bash shell script: vi .bashrc
— Add the following lines for TURBOMOLE:
export TURBODIR=/RS/progs/TURBOMOLE
PATH=$PATH:$TURBODIR/scripts
PATH=$PATH:$TURBODIR/bin/`sysname`
export PARA_ARCH=MPI
— Add the following lines for MOLPRO:
export INTEL_LICENSE_FILE=/RS/progs/intel/licenses/RS/progs/intel/impi/3.1/bin64/mpivars.sh
export molprodir=/RS/progs/molpro/mpp/2009/1-20/x86_64
PATH=$PATH:$molprodir/bin
export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:/RS/progs/molpro/mpp/2009/120/x86_64/lib
export PARA_ARCH=MPI
— Add the following lines for VESTA:
export PATH=$PATH:/RS/progs/VESTA-x86_64
Turbomole
15.06.2012
— Log
out and log in again then check: export
$TURBODIR
export $molprodir
Outline
— Introduction
— How to create input (define tool)
— Single point calculations at the Hartree-Fock, DFT and RIDFT levels
— Geometry optimization of water dimer at the RIMP2 level
— Geometry optimization of water dimer at the SCS-MP2 level
— Geometry optimization of water dimer at the CP-SCS-MP2 level
Turbomole
15.06.2012
Turbomole
— Developed by Reinhart Ahlrichs (University of Karlsruhe) 1989-2007
— Since 2007 TURBOMOLE GmbH founded by Ahlrichs is responsible
— for more information please visit to:
— http://www.turbomole.com/
— http://www.cosmologic.de
— A forum site is available under http://www.turbo-forum.com
Turbomole
15.06.2012
Turbomole usage philosophy
The usage of TURBOMOLE has been adapted to the way a UNIX user is working:
— Command line driven
— Many different programs, each one specialized for methods and/or properties
— Scripts are used to combine the functionalities and manage workflows
— Input can be changed by text editors
— Output processed by standard UNIX tools (editors, grep, awk, ...)
Turbomole
15.06.2012
Turbomole usage philosophy
Turbomole
15.06.2012
Turbomole modules
Turbomole
15.06.2012
Turbomole modules
DEFINE
interactive input generator which creates the input file control.
UFF
performs a geometry optimization at a force field level.
DSCF
for semi-direct SCF (self-consistent field) and DFT calculations.
GRAD
require a successful DSCF run and calculates the gradient of the energy with respect
to nuclear coordinates.
RIDFT and
RDGRAD
perform DFT calculations – as DSCF and DGRAD – within the RI-J approximation.
MPGRAD
requires a well converged SCF run and performs closed-shell RHF or UHF calculations
yielding single point MP2 energies and, if desired, the corresponding gradient.
RIMP2
calculates MP2 energies and gradients for RHF and UHF wavefunctions significantly
more efficient then MPGRAD.
RICC2
calculates electronic excitation energies, transition moments and properties of excited
states at the CIS, CIS(D), ADC(2) and CC2 level using a closed shell RHF or a UHF
SCF reference function. 15.06.2012
Turbomole
Turbomole modules
RELAX
requires a gradient run – by GRAD, RDGRAD, RIMP2 or MPGRAD and proposes a
new structure based on the gradient and the approximated force constants.
STATPT
performs a structure optimization using the “Trust Radius Image Minimization”
algorithm. It can be used to find minima or transition structure.
FROG
executes one molecular dynamics step.
AOFORCE
require a well converged SCF or DFT run and performs an analytic calculation
of force constants, vibrational frequencies and IR intensities.
ESCF
requires a well converged SCF or DFT run and calculates time dependent and dielectric
properties.
EGRAD
computes gradients and first-order properties of excited states.
MPSHIFT
computes NMR chemical shieldingns for all atoms of the molecules at the SCF, DFT
or MP2 level.
FREEH
calculates thermodynamic functions from molecular data in a control file: an
AOFORCE or a NUMFORCE run is a necessary prerequisite.
Turbomole
15.06.2012
How to create the input
— Build
define has some limited features to build a structure from scratch, but that
is neither convenient nor intuitive, except for some small cases.
→ Use an external builder (molden, ECCE, Hyperchem, Insight, Arguslab,
MAPS, Accelrys DS Visualizer, ...there are many!) and convert the structure
to xyz format.
— Convert
Use the TURBOMOLE script x2t to convert xyz files to TURBOMOLE
coordinates:
x2t struct.xyz > coord
Turbomole
15.06.2012
Define
define is an interactive input generator which creates the input file control:
—
Supports most basis sets in use, especially the only fully atom optimized consistent basis sets of
SVP, TZV, and QZV quality available for the atoms H–Rn (including lanthanides and actinides
up to Lawrencium for SVP and TZV)
— Determines the molecular symmetry
—
Determines the internal coordinates, allowing efficient geometry optimization
— Allows to perform a geometry optimization at a force field level to pre-optimize the geometry
and to calculate a Cartesian Hessian matrix
—
Sets the keywords necessary for single point calculations and geometry optimizations for a variety
of methods
—
Manipulate geometries of molecules
Turbomole
15.06.2012
Define
There are four major menus in define:
•Molecular Geometry Menu
•Atomic Attribute Definition Menu
•Occupation number & Molecular Orbital Definition Menu
•General Menu
Navigation within define:
1.
To finish the menu and proceed to the next one, enter:
*
2.
To go back to the previous menu, enter:
&
3.
If you are in one of the four main menus, hitting <Enter> will just print out the menu
you are in
4.
To quit define immediately (and prematurely), enter:
Turbomole
qq
15.06.2012
Define
To call define just type define in the shell : $define
— First you will be asked whether you want to define a new system or not
— If not press the ENTER button
Turbomole
15.06.2012
Define
Then you will be asked for a n input title:
— either give a name
— or press the ENTER button to proceed
Turbomole
15.06.2012
Define – molecular geometry menu
Turbomole
15.06.2012
Define – atomic attribute definition menu
Suppose we created the coord file via: x2t struct.xyz > coord
Select only these for
a geometry
optimization based
on cartesian coordinates
a coord
desy
ired
to load the geometry
to detect the symmetry, if you expect a higher symmetry, repeat
with increased tolerance desy 0.1
to use internal redundant coordinates
Then pressing * will lead us to the Atomic Attribute Definition Menu
Turbomole
15.06.2012
Define – atomic attribute definition menu
The default basis set is def-SV(P)
You can assign the same basis set for all atoms in the system by:
b all aug-cc-pVTZ
(Dunning’s basis set)
If you want to use a different basis set for any atomic species (let’s say for hydrogen):
b “h” aug-cc-pVDZ
bl
use to see the current basis set assignments
TURBOMOLE has a new set of basis sets abbreviated by def2:
def2-SV(P), def2-SVP, def2-TZVP, def2-TZVPP, def2-QZVP, def2-QZVPP
bi
(At least
Turbomole
information about improved TURBOMOLE basis sets
you must use def2-TZVP for 15.06.2012
a quantative analysis)
Define – atomic attribute definition menu
SV(P) or def-SV(P)
for routine SCF or DFT. Quality is about 6-31G*.
TZVP or def-TZVP
for accurate SCF or DFT. Quality is slightly better than 6-311G**.
TZVPP or def-TZVPP for MP2 or close to basis set limit SCF or DFT. Comparable
to 6-311G(2df).
QZVP and QZVPP for highly correlated treatments; quadruple zeta + 3d2f1g
or 4d2f1g (beyond Ne), 3p2d1f for H.
Turbomole
15.06.2012
Define – occupation number & orbital
definition menu
This menu serves for the definition of occupation numbers (and thereby the charge of the
molecule) and the start orbitals for the SCF calculations. For well-behaved molecules,
The procedure is very easy. Just choose the option:
eht
to perform an extended Huckel calculation.
Turbomole
15.06.2012
Define – occupation number & orbital
definition menu
— If you do not want them type n
— You can accept the default parameters by typing y, then you will be
asked for the charge of the system (default is 0)
Turbomole
15.06.2012
Define – occupation number & orbital
definition menu
— to accept the current occupation type y, then you will enter to the GENERAL MENU
Turbomole
15.06.2012
Define – general menu
Here, you can choose additional parameters for your calculation.
Turbomole
15.06.2012
Define – general menu - scf
Turbomole
15.06.2012
Define – general menu – mp2
— Do not forget to freeze (frozen core approximation)
— Do not forget to add cbas basis: the basis set selected in atomic attribute definition menu
— You can change memory (default is 60 MB)
— You can change convergence, use 0.1E-07
Turbomole
15.06.2012
Define – general menu – cc2
— CC2 is an approximated coupled cluster singles and doubles (CCSD) model
— The CC2 total energy is of second-order Moller-Plesset perturbation theory (MP2) quality
— The scaling of CC2 is N5 whereas the scaling of CCSD is N6, where N is the number of orbitals
— Similar to MP2 submenu
Turbomole
— Use always ri option
— Spin component scaled (SCS) MP2 is hidden in the ricc2 submenu
15.06.2012
For more info about CC2 check Chem. Phys. Lett. 243 409 (1995)
Define – general menu – cc2 – ricc2
Turbomole
15.06.2012
Define – general menu – dft
— switch on DFT
— select functional
— select gridsize
Type on to switch on DFT
(default functional is b-p)
(default gridsize is m3)
Turbomole
15.06.2012
Define – general menu – dft
Turbomole
15.06.2012
S22 test of interaction energies [JCP 132
(2010) 144104]
Mean Absolute Deviation
Turbomole
15.06.2012
Empirical Dispersion Correction for DFT
Calculations
— DFT functionals fails for intermolecular interaction energy calculations
(They do not result any binding between two monomers)
— This is due to insufficient treatment of electron correlation
(which is so important to describe dispersion forces)
— Stefan Grimme suggested dispersion correction to ordinary DFT energies
— When using the dispersion correction, the total energy is given by
— Edisp is an empirical dispersion correction given by
Turbomole
Nat is the number of atoms
C6ij denotes the dispersion coefficient for atom pair ij
15.06.2012on the selected DF
s6 global scaling factor depends
R is the interatomic distance between i and j
Define – general menu – dft
Turbomole offers grids 1 (coarse) to 7 (nest), and the multiple
grids m3 to m5. The latter employ a coarser grid during SCF
iterations, and grid 3 to grid 5 in the final SCF iteration and the
gradient evaluation. Default is grid m3, for clusters with more
than 50 atoms use m4.
Turbomole
15.06.2012
Define – general menu – RI-J (resolution
of identity)
— In the RI approximation two-electron integrals are approximated by three-center expansions
— RI procedure is recommended for non-hybrid functionals (which speeds up calculations
by a factor of 10 at least)
Turbomole
— activate ri with on
— assign jbas
15.06.2012
— then call jobex script with –ri option (jobex -ri)
Define – general menu – MARI-J
— RI-J calculations can be done even more effciently with the Multipole Accelerated
RI-J (MARI-J ) option, especially for larger molecules.
— Just rely on the defaults.
Turbomole
15.06.2012
Define – general menu – stp – transition
state search
— Stationary points are places on the potential energy surface (PES) with a zero gradient.
— Two types of stationary points are of special importance to chemists. These are minima
(reactants, products, intermediates) and first-order saddle points (transition states)
— The type of stationary point optimization depends on the value of irtvec. By default
irtvec is set to 0, which implies a structure minimization. A value irtvec > 0 implies a
transition state optimization using the eigenvalue-following TRIM (Trust Radius Image
Minimization) algorithm .
— Transition state search is invoked by
jobex –trans
Turbomole
15.06.2012
Vibrational normal modes
— After having performed a single point energy calculation or a geometry optimization,
just start the program aoforce (analytic force constant calculation):
aoforce > force.out
— All details will be written to force.out
— The keyword $vibrational spectrum in the control file contains a list of the modes:
— NumForce enables numerical force constant calculations for all levels of theory with a
gradient implemented. To run NumForce just issue the command
NumForce –ri –level mp2 -central
Turbomole
15.06.2012
Structure optimizations using the jobex script
— The shell script Jobex controls and executes automatic optimizations of molecular
geometry parameters.
— It will cycle through the direct SCF, gradient and force relaxation programs and stop if
either the maximum number of cycles is reached or the convergence criteria are fulfilled.
— Jobex has many options:
-energy integer
-gcart integer
-c integer
-dscf
-grad
-statpt
-relax
-trans
-level level
-ri
-ex
-md
Turbomole
converge total energy up to 10-integer Hartree (default:6)
converge maximum norm of cartesian gradient up to 10-integer a.u. (def:3)
perform up to integer cycles (default:20)
begin with a direct SCF step
begin with a gradient step
begin with a force relaxation step
use the RELAX program for force relaxation
perform transition state search
define the optimization level, level=scf, mp2, cc2 or uff
use RI modules
perform excited state geometry optimization using EGRAD
a molecular dynmaics run
15.06.2012
Turbomole output
— There will be an output written to file job.start.
— The convergence is signalled by the file converged; otherwise, you should find the file
not.converged within your working directory.
— If Jobex finds a file named stop or STOP in the working directory, Jobex will stop after the
present step has terminated. You can create stop by the command touch stop.
— The output of the last complete cycle is written to file job.last, while the output of the
running cycle is collected within the file job.<cycle>, where <cycle> is the index of the
cycle.
— Look at the file gradient: it contains all geometries and gradients at each step of
optimization:
grep cycle gradient
— To see the final geometry call
t2x –c > final.xyz
t2x > final_all.xyz (to see all trajectory)
— At the command line, use dist, bend or tors to get bond lengths, angles, dihedral angles
Turbomole
15.06.2012
Keywords in the control file
— The file control is the input file for TURBOMOLE.
— Keywords start with a ‘$’, e.g. $title.
Turbomole
$title
NO2 c2v UKS SVP
$operating system unix
$symmetry c2v
$coord file=coord
$intdef file=coord
$atoms
n 1 \
basis =n def-SVP
o 2-3 \
basis =o def-SVP
…
$energy file=energy
$grad file=grad
$forceapprox file=force
$lock off
$dft
functional b-p
gridsize m3
$last step define 15.06.2012
$end
Submitting Turbomole jobs on UHEM
— single job submission
#!/bin/bash
#BSUB -m karadeniz_temp1
#BSUB -a intelmpi
#BSUB -J triazine
#BSUB -o triazine.out
#BSUB -e triazine.err
#BSUB -q workshop
#BSUB -P workshop
#BSUB -R "span[ptile=1]"
#BSUB -n 1
# ege or anadolu
# bu kismi degistirmeyin !!!
# isinizi tanimlayacak bir isim
# bu kismi degistirmeyin !!!
# bu kismi degistirmeyin !!!
# kuyruk ismi
# proje ismi
# bu kismi degistirmeyin !!!
# kullanilacak olan islemci sayisi
echo 'Starting time:'
date
jobex -ri -level cc2 -energy 7 -gcart 4 -c 1000
echo 'Ending time:'
date
Turbomole
15.06.2012
Submitting Turbomole jobs on UHEM
— parallel job submission
#!/bin/bash
#BSUB -m karadeniz_temp1
#BSUB -m "b074 b075"
#BSUB -a intelmpi
#BSUB -J triazine
#BSUB -o triazine.out
#BSUB -e triazine.err
#BSUB -q workshop
#BSUB -P workshop
#BSUB -R "span[ptile=8]"
#BSUB -n 16
# submit the job to the b074 and b075 nodes
# bu kismi degistirmeyin !!!
# isinizi tanimlayacak bir isim
# bu kismi degistirmeyin !!!
# bu kismi degistirmeyin !!!
# kuyruk ismi
# proje ismi
# bu kismi degistirmeyin !!!
# kullanilacak olan islemci sayisi
export HOSTS_FILE=/AKDENIZ/HOME005/users/tekin/hostsfile
export PARNODES=16
echo 'Starting time:'
date
jobex -ri -level cc2 -energy 7 -gcart 4 -c 1000
echo 'Ending time:'
date
Turbomole
15.06.2012
check availability of nodes with bhosts command
b074
b074
b074
b074
b074
b074
b074
b074
b075
b075
b075
b075
b075
b075
b075
b075
Turbomole examples
— Single point calculations at the Hartree-Fock, DFT and RIDFT
levels
— Geometry optimization at the RIMP2 level
— Geometry optimization at the SCS-MP2 level
— Geometry optimization at the CP-SCS-MP2 level
Turbomole
15.06.2012
Single point calculations – Hartree Fock, DFT
and RIDFT energy calculation of benzene
1.
2.
3.
4.
Create an empty directory, e.g., hf_def2_tzvp_benzene
Call define
Hit <enter>
Give a title:
benzene, HF calculation
5. Load benzene from the structure library ($TURBODIR/structures)
a !benzene
6. define will scan through the structure library and ask you if you want to take the proposed
molecule
Turbomole
15.06.2012
Single point calculations – Hartree Fock, DFT
and RIDFT energy calculation of benzene
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.
17.
Just say y to accept the structure
Switch on the symmetry detection be entering:
desy
Next step is to provide internal coordinates. (That is not needed for single point calculations.)
ired
Now we are done with this menu and proceed to the next one by entering
*
This is the atomic attributes menu. Choose your basis set:
b all def2-TZVP
Now we are done with this menu and proceed to the next one by entering
*
Determine the starting molecular orbitals with an extended Huckel guess:
eht
Then you will be asked a lot of questions. Just accept the defaults.
Now you are in general menu.
The default method is Hartree-Fock (HF) calculation, and therefore we do not have to
change anything here. Simply enter:
*
To run a HF single point energy calculation,
call
Turbomole
15.06.2012
dscf > dscf.out
Single point calculations – Hartree Fock, DFT
and
RIDFT
energy
calculation
of
benzene
18. After the run, look in the output (dscf.out) and/or the energy file. The energy file should give:
$energy
SCF
1 -230.7848712224
$end
SCFKIN
230.6217054044
SCFPOT
-461.4065766268
19. Look in the output (dscf.out) for the dipole and quadrupole moments and whatever you like
to know.
20. Get more details about the orbitals, the HOMO-LUMO gap and occupation:
eiger
21. Let’s run a DFT calculation just after the HF calculation .
22. Call define and just accept all the defaults, then you will be in general menu.
23. Choose dft
24. Choose on to switch on DFT (Now you have chosen the default functional (b-p) and
gridsize m3)
25. Quit from define with *
26. Run dscf again.
27. Look again to the energy file:
$energy
SCF
1 -230.7848712224
Turbomole
2 -232.3365711535
$end
SCFKIN
230.6217054044
15.06.2012
231.4048783475
SCFPOT
-461.4065766268
-463.7414495010
Single point calculations – Hartree Fock,
DFT and RIDFT energy calculation of
benzene
28. Now, let’s switch on RI approximation following previous steps enter to the general menu
29. enter ri and then on.
30. Run ridft program:
ridft > ridft.out
$energy
SCF
1 -230.7848712224
2 -232.3365711535
3 -232.3368258829
$end
Turbomole
SCFKIN
230.6217054044
231.4048783475
231.4106397547
15.06.2012
SCFPOT
-461.4065766268
-463.7414495010
-463.7474656376
Do not be confused with auxiliary basis
1. Coulomb-Fitting ($jbas). Basis sets declared in this data-group are needed for RI-DFT
calculations, where only the Coulomb part of the Fock operator is approximated. They can
be kept rather small.
2. Coulomb- and Exchange-Fitting ($jkbas). If additionally the exchange part of the Fock
operator needs to be approximated (RI-SCF calculations), a larger fitting basis is needed.
It also can be used for RI-DFT runs (if the accuracy of the $jbas auxiliary basis is not
satisfactory for you).
3. Correlation-Fitting ($cbas). These basis sets are needed for RI-MP2 and RI-CC2
calculations. They contain basis functions with high angular momentum.
Turbomole
15.06.2012
Geometry optimizations – RIMP2, SCS-MP2 and CPSCS-MP2 calculations of water dimer
1. Open a new directory, e.g., rimp2_avdz_water
2. Create water.xyz file. Here are the coordinates:
6
O
H
H
O
H
H
Turbomole
0.0640488861
-0.0275603110
-0.8291328873
-0.0539423286
0.3402097464
0.3402097464
-1.5542001572
-0.6109761929
-1.8994223998
1.4751039841
1.8881478835
1.8881478835
3. Convert xyz file to turbomole coord file:
x2t water.xyz > coord
4. Call define
5. Hit <Enter>
6. Give a title
7. Read the coordinates:
a coord
15.06.2012
8. desy
9. ired
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.7547208224
-0.7547208224
Geometry optimizations – RIMP2, SCS-MP2 and
CP-SCS-MP2 calculations of water dimer
10. Type * to go to the next menu
11. Choose basis set:
b all aug-cc-pVDZ
12. Type * to go to the next menu
13. Activate eht and accept the defaults
14. Now you are in the general menu
15. First select scf
16. Change scf option conv and then type 8 and go back to the general menu by pressing to ENTE
17. Now, select mp2
18. Activate freeze and go back to the mp2 submenu by pressing to *
19. Activate cbas: Automatically the auxiliary basis set will be loaded.
20. Press * to go back to the mp2 submenu.
21. Select the denconv and enter 0.1E-07.
22. Press * to go back to the mp2 submenu.
23. Press * to quit define.
24. Run the jobex script with the following parameters:
jobex -ri -level mp2 -energy 7 -gcart 4 -c 1000
25. Check the convergence by grep cycle gradient
26. After the convergence convert the turbomole geometry to xyz:
t2x –c > rimp2_opt.xyz or t2x > final_all_trajectory.xyz
Geometry optimizations – RIMP2, SCS-MP2 and
CP-SCS-MP2 calculations of water dimer
1. Open a new folder, e.g., scsmp2_avdz_water
2. Copy the rimp2 xyz output here:
cp ../rimp2_avdz_water/rimp2_opt.xyz .
3. Convert the xyz file to turbomole file:
x2t rimp2_opt.xyz > coord
4. Call define
5. a coord
6. desy
7. ired
8. *
9. b all aug-cc-pVDZ
10. *
11. Activate eht and accept the defaults
12. Now you are in the general menu
13. First select scf
14. Change scf option conv and then type 8 and go back to the general menu by pressing to ENTE
15. Now, select cc2
16. Activate freeze and go back to the cc2 submenu by pressing to *
17. Activate cbas: Automatically the auxiliary basis set will be loaded.
18. Press * to go back to the cc2 submenu.
19. Select the denconv and enter 0.1E-07.
Geometry optimizations – RIMP2, SCS-MP2 and
CP-SCS-MP2 calculations of water dimer
20. Press * to go back to the cc2 submenu.
21. Activate ricc2
22. Switch on scs correction by typing scs
23. Choose geoopt model:
geoopt mp2
24. *
25. *
26. *
27. vi control to check the ricc2 entries:
$ricc2
mp2 energy only
geoopt model=mp2
scs
cos= 1.20000
state=(x)
css= 0.33333
28. If there is nothing after mp2, enter “energy only” by hand
29. Run the jobex script with the following parameters:
jobex -ri -level cc2 -energy 7 -gcart 4 -c 1000
30. Check the convergence by grep cycle gradient
31. After the convergence convert the turbomole geometry to xyz:
t2x –c > scsmp2_opt.xyz or t2x > final_all_trajectory.xyz
Geometry optimizations – RIMP2, SCS-MP2 and
CP-SCS-MP2 calculations of water dimer
1. Open a new folder, e.g., cp_scsmp2_avdz_water
2. Copy the scsmp2 xyz output here:
cp ../scsmp2_avdz_water/scsmp2_opt.xyz .
3. Convert the xyz file to turbomole file:
x2t scsmp2_opt.xyz > coord
4. Call define
5. a coord
6. desy
7. ired
8. frag
Turbomole
9. Switch frag on by typing on
15.06.2012
Geometry optimizations – RIMP2, SCS-MP2 and
CP-SCS-MP2 calculations of water dimer
Turbomole
15.06.2012
Two times press <ENTER> button to see the fragment numbers. At most you can define
three fragments.
Geometry optimizations – RIMP2, SCS-MP2 and
CP-SCS-MP2 calculations of water dimer
—First three atoms belong to the first water molecule and the last three atoms belong to the
Second water molecule. You can assign them like below:
Turbomole
15.06.2012
Geometry optimizations – RIMP2, SCS-MP2 and
CP-SCS-MP2 calculations of water dimer
10. Enter to leave the frag menu
11. *
12. b all aug-cc-pVDZ
13. *
14. Activate eht and accept the defaults
15. Now you are in the general menu
16. First select scf
17. Change scf option conv and then type 8 and go back to the general menu by pressing to ENTE
18. Now, select cc2
19. Activate freeze and go back to the cc2 submenu by pressing to *
20. Activate cbas: Automatically the auxiliary basis set will be loaded.
21. Press * to go back to the cc2 submenu.
22. Select the denconv and enter 0.1E-07.
20. Press * to go back to the cc2 submenu.
21. Activate ricc2
22. Switch on scs correction by typing scs
23. Choose geoopt model:
geoopt mp2
24. *
Turbomole
15.06.2012
25. *
26. *
Geometry optimizations – RIMP2, SCS-MP2
and CP-SCS-MP2 calculations of water
dimer
27. vi control to check the ricc2 entries:
$ricc2
mp2 energy only
geoopt model=mp2
scs
cos= 1.20000
state=(x)
css= 0.33333
28. If there is nothing after mp2, enter “energy only” by hand
29. Run the jobex script with the following parameters:
jobbsse -ri -level mp2 –opt
30. Check the convergence by grep cycle gradient
31. After the convergence convert the turbomole geometry to xyz:
t2x –c > cp_scsmp2_opt.xyz or t2x > final_all_trajectory.xyz
Turbomole
15.06.2012
Geometry optimizations – Comparison of RIMP2, SCS-MP2
and CP-SCS-MP2 water dimer geometries
RIMP2
Turbomole
SCS-MP2
15.06.2012
CP-SCS-MP2
COSMO – Dealing with solvation effects
1. Open a new folder, e.g., cosmo_benzene
2. Load benzene from the structure library with
a !benzene
3. desy
4. ired
5. *
6. b all def-TZVP
7. *
8. eht and many times <Enter>
9. Switch on dft and ri
10. Exit from define with *
11. cosmoprep
Turbomole
15.06.2012
http://www.gaussian.com/
g_tech/g_ur/k_scrf.htm
Water: ε=78.3553
Acetonitrile: ε=35.688
Methanol: ε=32.613
Ethanol: ε=24.852
IsoQuinoline: ε=11.00
Quinoline: ε=9.16
Chloroform: ε=4.7113
DiethylEther: ε=4.2400
Dichloromethane: ε=8.93
DiChloroEthane: ε=10.125
CarbonTetraChloride: ε=2.2280
Benzene: ε=2.2706
Toluene: ε=2.3741
COSMO – Dealing with solvation effects
Turbomole
15.06.2012
COSMO – Dealing with solvation effects
12. r all o
if you have atoms
like “mg”
r “mg” b
13. *
14. <Enter>
15. jobex –ri > jobex.out
Turbomole
15.06.2012

similar documents