|
|
|
Sample directory: adf/e_QMMM_Surface/
This is an example of a Ziegler-Natta type catalytic system: a TiCl complex embedded in a MgCl surface with two organic substrates also attached to the surface. To make the computation faster, the QMMM approach is applied. The QM part includes only the active site and a piece of the MgCl surface.
The computation is formally a geometry optimization, but to keep the sample doable in a reasonable time the sample performs only one geometry update step. In the optimization, all of the MgCl surface atoms are frozen.
The standard force field has been modified to accommodate this calculation. The modified force field file is part of the sample run script. In this modified file, bonds are defined between Mg-Cl atoms in the MM connection table. This results in some torsions where the atoms are collinear. To rectify this problem, the torsional potentials for these atoms are set to potential type '0' (no potential).
There are no capping atoms mediating the bonds between the QM and MM regions because the boundary goes through the MgCl surface, which is ionically bound.
*************************************************************************
SYBYL/TRIPOS FORCE FIELD FILE FOR ADF QM/MM
MODIFIED WITH UFF1.01 FOR Si Mg Ti Cl
L. Petitjean 15.11.1999
*************************************************************************
(Most of the contents of the modified force field file is omitted here. You quickly get the difference with the standard sybyl force field file in the ADF database by running a UNIX diff on the two files.
eor
TITLE
noprint sfo,frag,functions
! keywords for calculation methods and optimization
XC
LDA VWN
GGA BLYP
END
GEOMETRY selected cartesian
ITERATIONS 1
HESSUPD BFGS
CONVERGE e=1e-4 grad=1e-3 rad=1e-2
STEP rad=0.15
DIIS ndiis=5 ok=10.0 cyc=3
END
The "Iterations 1" subkey specification in the Geometry block specifies that only one step in the optimization is carried out.
SCF
ITERATIONS 250
CONVERGE 1E-6 1E-6
MIXING 0.2
DIIS N=10 OK=0.5 cyc=5 CX=5.0 BFAC=0
END
! keywords for molecule specification
CHARGE 0 0
ATOMS cartesian
1 Mg x1 y1 z1
(all other atoms in the Atoms block omitted here)
GEOVAR
x1=.00000 F
y1=.00000 F
z1=.00000 F
x2=.00000 F
y2=1.72129 F
z2=1.82068 F
x3=.00000 F
y3=.00000 F
z3=-3.64100 F
x4=.00000 F
y4=-1.72130 F
z4=-1.82068 F
x5=.00000 F
y5=1.72130 F
z5=-1.82032 F
x6=.00000 F
y6=1.72130 F
z6=-5.46132 F
x7=2.53903
y7=.03004
Many of the coordinates have a "F" after their initial value specification under Geovar, indicating that these coordinates will be kept frozen during optimization.
The remaining initial value specifications are omitted here.
QMMM
OPTIMIZE
MAX_STEPS 3000
MAX_GRADIENT 0.01
METHOD BFGS
PRINT_CYCLES 100
SUBEND
FORCE_FIELD_FILE ´pwd´/champ_de_force.ff
The local file "champ_de_force.ff" is used as force field file. Of course, this is the file we've just set up in the run script.
WARNING_LEVEL=1
ELSTAT_COUPLING_MODEL=1
MM_CONNECTION_TABLE
Contents of the MM_CONNECTION_TABLE block is omitted.
CHARGES
1 .957
2 -.608
Initial charges are specified for (all) the atoms. Whether or not the charges on the QM (and LI) atoms are used depends on the type of electrostatic coupling between the QM and MM system. See the QMMM manual for details.
SUBEND
END
FRAGMENTS
Ti t21.Ti.2p.IV
Cl t21.Cl.2p.III
Mg t21.Mg.2p.III
C t21.C.1s.III
H t21.H.III
END
END INPUT
eor