| 1 |
gezelter |
785 |
// This is the forcefield file for the Clay Force Field (CLAYFF) |
| 2 |
|
|
// |
| 3 |
|
|
// Details can be found in the following article: |
| 4 |
|
|
// |
| 5 |
|
|
// "Molecular Models of Hydroxide, Oxyhydroxid, and Clay Phases and |
| 6 |
|
|
// the Development of a General Force Field" by Randall T. Cygan, |
| 7 |
|
|
// Jian-Jie Liang, and Andrey G. Kalinichev, J. Phys. Chem. B 108, |
| 8 |
|
|
// pp. 1255-1266 (2004). |
| 9 |
|
|
// |
| 10 |
|
|
// This file defines the following atom types: |
| 11 |
|
|
// |
| 12 |
|
|
// h* water hydrogen |
| 13 |
|
|
// ho hydroxyl hydrogen |
| 14 |
|
|
// o* water oxygen |
| 15 |
|
|
// oh hydroxyl oxygen |
| 16 |
|
|
// ob bridging oxygen |
| 17 |
|
|
// obos bridging oxygen with octahedral substitution |
| 18 |
|
|
// obts bridging oxygen with tetrahedral substitution |
| 19 |
|
|
// obss bridging oxygen with double substitution |
| 20 |
|
|
// ohs hydroxyl oxygen with substitution |
| 21 |
|
|
// st tetrahedral silicon |
| 22 |
|
|
// ao octahedral aluminum |
| 23 |
|
|
// at tetrahedral aluminum |
| 24 |
|
|
// mgo octahedral magnesium |
| 25 |
|
|
// mgh hydroxide magnesium |
| 26 |
|
|
// cao octahedral calcium |
| 27 |
|
|
// cah hydroxide calcium |
| 28 |
|
|
// feo octahedral iron |
| 29 |
|
|
// lio octahedral lithium |
| 30 |
|
|
// Na aqueous sodium ion |
| 31 |
|
|
// K aqueous potassium ion |
| 32 |
|
|
// Cs aqueous cesium ion |
| 33 |
|
|
// Ca aqueous calcium ion |
| 34 |
|
|
// Ba aqueous barium ion |
| 35 |
|
|
// Cl aqueous chloride ion |
| 36 |
|
|
|
| 37 |
|
|
begin Options |
| 38 |
|
|
Name = "CLAYFF" |
| 39 |
|
|
vdWtype = "Lennard-Jones" |
| 40 |
|
|
DistanceMixingRule = "Arithmetic" |
| 41 |
|
|
DistanceType = "Rmin" |
| 42 |
|
|
EnergyMixingRule = "Geometric" |
| 43 |
|
|
vdW-14-scale = 0.0 |
| 44 |
|
|
electrostatic-14-scale = 0.0 |
| 45 |
|
|
dielectric = 1.0 |
| 46 |
|
|
end Options |
| 47 |
|
|
|
| 48 |
|
|
begin AtomTypes |
| 49 |
|
|
//Name mass (amu) |
| 50 |
|
|
h* 1.0079 |
| 51 |
|
|
ho 1.0079 |
| 52 |
|
|
o* 15.9994 |
| 53 |
|
|
oh 15.9994 |
| 54 |
|
|
ob 15.9994 |
| 55 |
|
|
obos 15.9994 |
| 56 |
|
|
obts 15.9994 |
| 57 |
|
|
obss 15.9994 |
| 58 |
|
|
ohs 15.9994 |
| 59 |
|
|
st 28.0855 |
| 60 |
|
|
ao 26.981538 |
| 61 |
|
|
at 26.981538 |
| 62 |
|
|
mgo 24.3050 |
| 63 |
|
|
mgh 24.3050 |
| 64 |
|
|
cao 40.078 |
| 65 |
|
|
cah 40.078 |
| 66 |
|
|
feo 55.845 |
| 67 |
|
|
lio 6.941 |
| 68 |
|
|
Na 22.98977 |
| 69 |
|
|
K 39.0983 |
| 70 |
|
|
Cs 132.90545 |
| 71 |
|
|
Ca 40.078 |
| 72 |
|
|
Ba 137.327 |
| 73 |
|
|
Cl 35.453 |
| 74 |
|
|
end AtomTypes |
| 75 |
|
|
|
| 76 |
|
|
begin LennardJonesAtomTypes |
| 77 |
|
|
//Name epsilon R0 |
| 78 |
|
|
o* 0.1554 3.5532 |
| 79 |
|
|
oh 0.1554 3.5532 |
| 80 |
|
|
ob 0.1554 3.5532 |
| 81 |
|
|
obos 0.1554 3.5532 |
| 82 |
|
|
obts 0.1554 3.5532 |
| 83 |
|
|
obss 0.1554 3.5532 |
| 84 |
|
|
ohs 0.1554 3.5532 |
| 85 |
|
|
st 1.8405e-6 3.7064 |
| 86 |
|
|
ao 1.3298e-6 4.7943 |
| 87 |
|
|
at 1.8405e-6 3.7064 |
| 88 |
|
|
mgo 9.0298e-7 5.9090 |
| 89 |
|
|
mgh 9.0298e-7 5.9090 |
| 90 |
|
|
cao 5.0298e-6 6.2484 |
| 91 |
|
|
cah 5.0298e-6 6.2428 |
| 92 |
|
|
feo 9.0298e-6 5.5070 |
| 93 |
|
|
lio 9.0298e-6 4.7257 |
| 94 |
|
|
Na 0.1301 2.6378 |
| 95 |
|
|
K 0.1000 3.7423 |
| 96 |
|
|
Cs 0.1000 4.3002 |
| 97 |
|
|
Ca 0.1000 3.2237 |
| 98 |
|
|
Ba 0.0470 4.2840 |
| 99 |
|
|
Cl 0.1001 4.9388 |
| 100 |
|
|
end LennardJonesAtomTypes |
| 101 |
|
|
|
| 102 |
|
|
begin ChargeAtomTypes |
| 103 |
|
|
// Name charge |
| 104 |
|
|
h* 0.4100 |
| 105 |
|
|
ho 0.4250 |
| 106 |
|
|
o* -0.8200 |
| 107 |
|
|
oh -0.9500 |
| 108 |
|
|
ob -1.0500 |
| 109 |
|
|
obos -1.1808 |
| 110 |
|
|
obts -1.1688 |
| 111 |
|
|
obss -1.2996 |
| 112 |
|
|
ohs -1.0808 |
| 113 |
|
|
st 2.1000 |
| 114 |
|
|
ao 1.5750 |
| 115 |
|
|
at 1.5750 |
| 116 |
|
|
mgo 1.3600 |
| 117 |
|
|
mgh 1.0500 |
| 118 |
|
|
cao 1.3600 |
| 119 |
|
|
cah 1.0500 |
| 120 |
|
|
feo 1.5750 |
| 121 |
|
|
lio 0.5250 |
| 122 |
|
|
Na 1.0 |
| 123 |
|
|
K 1.0 |
| 124 |
|
|
Cs 1.0 |
| 125 |
|
|
Ca 2.0 |
| 126 |
|
|
Ba 2.0 |
| 127 |
|
|
Cl -1.0 |
| 128 |
|
|
end ChargeAtomTypes |
| 129 |
|
|
|
| 130 |
|
|
begin BondTypes |
| 131 |
|
|
|
| 132 |
|
|
//Atom1 Atom2 Harmonic b0 Kb (kcal/mol) |
| 133 |
|
|
//V_Harmonic = 0.5*Kb(b- bo)^2 |
| 134 |
|
|
|
| 135 |
|
|
o* h* Harmonic 1.0000 277.06745 |
| 136 |
|
|
oh ho Harmonic 1.0000 277.06745 |
| 137 |
|
|
ohs ho Harmonic 1.0000 277.06745 |
| 138 |
|
|
|
| 139 |
|
|
end BondTypes |
| 140 |
|
|
|
| 141 |
|
|
begin BendTypes |
| 142 |
|
|
|
| 143 |
|
|
//Harmonic |
| 144 |
|
|
//Atom1 Atom2 Atom3 Harmonic Theta0 Ktheta |
| 145 |
|
|
//V_Harmonic = 0.5*Ktheta(Theta - Theta0)^2 |
| 146 |
|
|
//Ktheta: kcal/mole/rad**2 |
| 147 |
|
|
//Theta0: degrees |
| 148 |
|
|
|
| 149 |
|
|
h* o* h* Harmonic 109.47 22.8848 |
| 150 |
|
|
ao oh ho Harmonic 109.47 15.0 |
| 151 |
|
|
ao ohs ho Harmonic 109.47 15.0 |
| 152 |
|
|
at oh ho Harmonic 109.47 15.0 |
| 153 |
|
|
at ohs ho Harmonic 109.47 15.0 |
| 154 |
|
|
mgh oh ho Harmonic 109.47 15.0 |
| 155 |
|
|
mgh ohs ho Harmonic 109.47 15.0 |
| 156 |
|
|
cah oh ho Harmonic 109.47 15.0 |
| 157 |
|
|
cah ohs ho Harmonic 109.47 15.0 |
| 158 |
|
|
feo oh ho Harmonic 109.47 15.0 |
| 159 |
|
|
feo ohs ho Harmonic 109.47 15.0 |
| 160 |
|
|
lio oh ho Harmonic 109.47 15.0 |
| 161 |
|
|
lio ohs ho Harmonic 109.47 15.0 |
| 162 |
|
|
|
| 163 |
|
|
end BendTypes |
