Homec4science
Diffusion lammps (efficient_neuronet) Browse

lammps/tools/moltemplate/common/tip3p_attempts/tip3pcharmm.lt
c1b6654dbc9cefficient_neuronet

Last Change
Raw File

tip3pcharmm.lt
View Options

# file "tip3p_charmm.lt"
#
# H1 H2
# \ /
# O
#
# I think this is the TIP3P water model used by CHARMM (and probably AMBER)
# It is (mostly) based on this paper:
# Jorgensen, Chandrasekhar, Madura, Impey, Klein, J Chem Phys, 79, 926 (1983)
TIP3Pcharmm {
write_once("In Init") {
# -- Default styles (for solo "TIP3Pcharmm" water) --
units real
atom_style full
# I'm not sure exactly which cutoffs distances are traditionally used in
# the TIP3P water model used by CHARMM.
# (See the Price JCP 2004 paper for a review.)
# pair_style hybrid lj/charmm/coul/charmm 7.5 8.0 10.0 10.5
# Try this instead:
pair_style hybrid lj/charmm/coul/charmm 10.0 10.5 10.0 10.5
bond_style hybrid harmonic
angle_style hybrid harmonic
pair_modify mix arithmetic
}
write("Data Atoms") {
$atom:O $mol:. @atom:O -0.834 0.0000000 0.00000 0.000000
$atom:H1 $mol:. @atom:H 0.417 0.756950327 0.00000 0.5858822766
$atom:H2 $mol:. @atom:H 0.417 -0.756950327 0.00000 0.5858822766
}
write_once("Data Masses") {
@atom:O 15.9994
@atom:H 1.008
}
write("Data Bonds") {
$bond:OH1 @bond:OH $atom:O $atom:H1
$bond:OH2 @bond:OH $atom:O $atom:H2
}
write("Data Angles") {
$angle:HOH @angle:HOH $atom:H1 $atom:O $atom:H2
}
write_once("In Settings") {
bond_coeff @bond:OH harmonic 450.0 0.9572
angle_coeff @angle:HOH harmonic 55.0 104.52
#########################################################################
#### There are two choices for for the O-O interactions
#########################################################################
#### O-O nonbonded interactions
# For the 1983 Jorgensen version of TIP3P use:
pair_coeff @atom:O @atom:O lj/charmm/coul/charmm 0.1521 3.1507
# For the 2004 Price & Brooks version of TIP3P use:
# pair_coeff @atom:O @atom:O lj/charmm/coul/charmm 0.102 3.188
#########################################################################
#### There are three choices for for the O-H and H-H interactions
#########################################################################
#### 1) CHARMM uses an arithmetic mixing-rule for the O-H sigma parameter
pair_coeff @atom:H @atom:H lj/charmm/coul/charmm 0.0460 0.4000
pair_coeff @atom:O @atom:H lj/charmm/coul/charmm 0.0836 1.7753 #arithmetic
#########################################################################
#### 2) OPLS-AA uses geometric a mixing-fule for the O-H sigma parameter,
#### If you want to use this, uncomment the following two lines:
# pair_coeff @atom:H @atom:H lj/charmm/coul/charmm 0.0460 0.4000
# pair_coeff @atom:O @atom:H lj/charmm/coul/charmm 0.0836 1.1226 #geometric
#########################################################################
#### 3) The original Jorgensen 1983 parameterization has no OH or HH
# lennard-jones interactions. For this behavior, uncomment these lines:
# pair_coeff @atom:H @atom:H lj/charmm/coul/charmm 0.00 0.4000
# pair_coeff @atom:O @atom:H lj/charmm/coul/charmm 0.00 1.7753
#########################################################################
# Define a group for the tip3p water molecules:
group tip3p type @atom:O @atom:H
# Optional: Constrain the angles and distances.
# (Most implementations use this, but it is optional.)
fix fSHAKE tip3p shake 0.0001 10 100 b @bond:OH a @angle:HOH
# (Remember to "unfix" fSHAKE during minimization.)
}
} # "TIP3Pcharmm" water molecule type
c4science · Help