| 1 |
// This is the forcefield file for the Dipolar Unified-atom Force Field (DUFF) |
| 2 |
// |
| 3 |
// The sections are divided into AtomTypes |
| 4 |
// |
| 5 |
// many parameters (but not all) are derived from the TRAPPE force field |
| 6 |
// of Siepmann's group. |
| 7 |
|
| 8 |
|
| 9 |
begin BaseAtomTypes |
| 10 |
//Name mass (amu) |
| 11 |
O_TIP4P_FQ 15.9994 |
| 12 |
H_TIP4P_FQ 1.0079 |
| 13 |
M_TIP4P_FQ 0.0 |
| 14 |
end BaseAtomTypes |
| 15 |
|
| 16 |
begin LennardJonesAtomTypes |
| 17 |
//Name epsilon sigma |
| 18 |
O_TIP4P_FQ 0.2862 3.159 |
| 19 |
end LennardJonesAtomTypes |
| 20 |
|
| 21 |
begin ChargeAtomTypes |
| 22 |
// Name charge |
| 23 |
//O_TIP4P_FQ 0.0 |
| 24 |
//H_TIP4P_FQ 0.0 |
| 25 |
//M_TIP4P_FQ 0.0 |
| 26 |
end ChargeAtomTypes |
| 27 |
|
| 28 |
|
| 29 |
begin FluctuatingChargeAtomTypes |
| 30 |
//Name chi (kcal/mol e) hardness (kcal/(mol*e^2)) slaterN slaterZeta ChargeMass (fs^2 kcal/ e^2 mol) |
| 31 |
//Values obtained from Rick, Stuart, and Berne J Chem Phys, 1994, 101 |
| 32 |
//Since they gave chi as an electronegativity difference, we will assign Oxygen's (the M site |
| 33 |
//actually for TIP4P) electronegativity from Rappe and Goddard III, and setting H to equal |
| 34 |
//the difference |
| 35 |
// Rappe and Goddard, O or M chi: 8.741 eV = 201.5705 kcal/mol |
| 36 |
// Rick, Stuart, and Berne, delta Chi between O and H = 68.49 |
| 37 |
// SlaterZeta has units of Bohr radius |
| 38 |
// Fictitious ChargeMass has units of (fs / e)^2 kcal / mol (note the ps->fs difference between |
| 39 |
// this code and the Rick, Stuart, and Berne paper |
| 40 |
H_TIP4P_FQ 133.0805 353.0 1 0.90 60.0 |
| 41 |
M_TIP4P_FQ 201.5705 371.6 2 1.63 60.0 |
| 42 |
end FluctuatingChargeAtomTypes |
