From c936d176fc9bea328b097d73e13aeb153a207d4f Mon Sep 17 00:00:00 2001
From: sjplimp These are pair styles contributed by users, which can be used if
@@ -426,8 +426,8 @@ angle potentials. Click on the style itself for a full description:
These are angle styles contributed by users, which can be used if
@@ -460,7 +460,7 @@ description:
dpd/tstat dsmc eam eam/opt eam/alloy eam/alloy/opt eam/fs eam/fs/opt
-eim gauss gayberne gayberne/gpu
-gran/hertz/history gran/hooke gran/hooke/history lj/charmm/coul/charmm
-lj/charmm/coul/charmm/implicit lj/charmm/coul/long lj/charmm/coul/long/opt lj/class2
-lj/class2/coul/cut lj/class2/coul/long lj/cut lj/cut/gpu
-lj/cut/opt lj/cut/coul/cut lj/cut/coul/debye lj/cut/coul/long
-lj/cut/coul/long/tip4p lj/expand lj/gromacs lj/gromacs/coul/gromacs
-lj/smooth lj96/cut lubricate meam
-morse morse/opt peri/lps peri/pmb
-reax resquared soft sw
-table tersoff tersoff/zbl yukawa yukawa/colloid
+
+gran/hertz/history gran/hooke gran/hooke/history hbond/dreiding/lj
+hbond/dreiding/morse lj/charmm/coul/charmm lj/charmm/coul/charmm/implicit lj/charmm/coul/long
+lj/charmm/coul/long/opt lj/class2 lj/class2/coul/cut lj/class2/coul/long
+lj/cut lj/cut/gpu lj/cut/opt lj/cut/coul/cut
+lj/cut/coul/debye lj/cut/coul/long lj/cut/coul/long/tip4p lj/expand
+lj/gromacs lj/gromacs/coul/gromacs lj/smooth lj96/cut
+lubricate meam morse morse/opt
+peri/lps peri/pmb reax resquared
+soft sw table tersoff tersoff/zbl yukawa yukawa/colloid
-none hybrid charmm class2
-cosine cosine/delta cosine/squared harmonic table
+
+cosine cosine/delta cosine/periodic cosine/squared harmonic table
diff --git a/doc/Section_commands.txt b/doc/Section_commands.txt
index 9dfb33cff7..53997e819d 100644
--- a/doc/Section_commands.txt
+++ b/doc/Section_commands.txt
@@ -579,6 +579,8 @@ potentials. Click on the style itself for a full description:
"gran/hertz/history"_pair_gran.html,
"gran/hooke"_pair_gran.html,
"gran/hooke/history"_pair_gran.html,
+"hbond/dreiding/lj"_pair_hbond_dreiding.html,
+"hbond/dreiding/morse"_pair_hbond_dreiding.html,
"lj/charmm/coul/charmm"_pair_charmm.html,
"lj/charmm/coul/charmm/implicit"_pair_charmm.html,
"lj/charmm/coul/long"_pair_charmm.html,
@@ -657,6 +659,7 @@ angle potentials. Click on the style itself for a full description:
"class2"_angle_class2.html,
"cosine"_angle_cosine.html,
"cosine/delta"_angle_cosine_delta.html,
+"cosine/periodic"_angle_cosine_periodic.html,
"cosine/squared"_angle_cosine_squared.html,
"harmonic"_angle_harmonic.html,
"table"_angle_table.html :tb(c=4,ea=c,w=100)
@@ -695,7 +698,8 @@ description:
"hybrid"_improper_hybrid.html,
"class2"_improper_class2.html,
"cvff"_improper_cvff.html,
-"harmonic"_improper_harmonic.html :tb(c=4,ea=c,w=100)
+"harmonic"_improper_harmonic.html,
+"umbrella"_improper_umbrella.html :tb(c=4,ea=c,w=100)
:line
diff --git a/doc/Section_howto.html b/doc/Section_howto.html
index 6fc10741c4..0924aab3ff 100644
--- a/doc/Section_howto.html
+++ b/doc/Section_howto.html
@@ -16,7 +16,7 @@ certain kinds of LAMMPS simulations.
The example input scripts included in the LAMMPS distribution and highlighted in this section also show how to @@ -167,19 +166,18 @@ the same as in 3d.
There are many different ways to compute forces in the CHARMM -and AMBER molecular dynamics codes, only some of which are -available as options in LAMMPS. A force field has 2 parts: the -formulas that define it and the coefficients used for a particular -system. Here we only discuss formulas implemented in LAMMPS. Setting +
A force field has 2 parts: the formulas that define it and the +coefficients used for a particular system. Here we only discuss +formulas implemented in LAMMPS that correspond to formulas commonly +used in the CHARMM, AMBER, and DREIDING force fields. Setting coefficients is done in the input data file via the read_data command or in the input script with commands like pair_coeff or -bond_coeff. See this section for -additional tools that can use CHARMM or AMBER to assign force field -coefficients and convert their output into LAMMPS input. +bond_coeff. See this section +for additional tools that can use CHARMM or AMBER to assign force +field coefficients and convert their output into LAMMPS input.
See (MacKerell) for a description of the CHARMM force field. See (Cornell) for a description of the AMBER force @@ -193,16 +191,56 @@ field. with common options in CHARMM or AMBER. See each command's documentation for the formula it computes.
-DREIDING is a generic force field developed by the Goddard +group at Caltech and is useful for +predicting structures and dynamics of organic, biological and +main-group inorganic molecules. The philosophy in DREIDING is to use +general force constants and geometry parameters based on simple +hybridization considerations, rather than individual force constants +and geometric parameters that depend on the particular combinations of +atoms involved in the bond, angle, or torsion terms. DREIDING has an +explicit hydrogen bond term to describe +interactions involving a hydrogen atom (H___A) on very electronegative +atoms (N, O, F). +
+See (Mayo) for a description of the DREIDING force field +
+These style choices compute force field formulas that are consistent +with the DREIDING force field. See each command's +documentation for the formula it computes. +
+Thermodynamic integration is a widely used method to compute free -energies from atomistic simulations. LAMMPS can be used to run -thermodynamic integration calculations using the methods discussed in -this section and the fix adapt command. Currently, -it is capable of the transformations essential for computing melting -points using the pseudo-supercritical path method developed by Eike -and Maginn. -
-See the examples/TI directory for more information and sample files -that compute a melting point using the techniques described in this -section. That directory has its own README file. See also the paper -by Jayaraman for an example of using this implementation -of thermodynamic integration in LAMMPS to compute melting points of -alkali nitrate salts, using the steps outlined here. -
-In this method, three intermediate "pseudo-supercritical" states are -accessed in the transformation between the liquid and solid -phases. These pseudo-states are a weak liquid, a dense weak liquid, -and an ordered weak phase. The transformation between the liquid and -solid states can also be driven uisng the fix adapt -command. -
-For the transformation from the liquid to the weak phase, the -intermolecular interactions need to be weakened. Appropriate scale -factors, computed by variables you define, and applied to pair styles -by fix adapt, can be used to do this, as in the -example scripts. The compute ti command can -accumulate the value of dU/dlambda. See Jayaraman and -Maginn for more information about calculating a -free energy from dU/dlambda. -
-IMPORTANT NOTE: The pair styles that fix adapt can scale on-the-fly -are listed on the fix adapt doc page. interaction scaling -is desired. If a pair style is not on that list, it is generally -quite easy to add an extract() method to the pair style, to enable fix -adapt to rescale it. -
-Step 2 is the transformation of the simulation box density from the -liquid phase to that of the equilibrated crystal. The parameters for -box1 and box2 should be obtained from equilibrated NPT simulations of -the liquid and crystal phases and used in a fix -deform command to change the box size and/or shape. -It also advisable to use fix adapt on the pair styles -to prevent overlaps which may occur during the box transformation. -
-In step 3, the dense, weak system is transformed to an ordered state, -which has the same ordering as in the equilibrated crystal. Ordering -is achieved by introducing an attractive potential between atoms and -lattice sites. These lattice sites can be calculated as the mean -positions of the atoms in an equilibrium simulation of the -crystal. The pair/gauss command can be used to -introduce an attractive Gaussian potential between the atoms and their -corresponding lattice sites. The prefactor of the Gaussian pair -potential can be scaled by fix adapt to turn on the -attractions. Again, the quantity dU/dlambda can be tracked via the -compute ti command. -
-Step 4 is the transformation of the ordered state to the final -crystal. In this step, the intermolecular interactions are scaled -back to full strength, while the Gaussian tethers are removed, all via -fix adapt. -
-(MacKerell) MacKerell, Bashford, Bellott, Dunbrack, Evanseck, Field, Fischer, Gao, Guo, Ha, et al, J Phys Chem, 102, 3586 (1998).
+ + +(Mayo) Mayo, Olfason, Goddard III, J Phys Chem, 94, 8897-8909 +(1990). +
(Jorgensen) Jorgensen, Chandrasekhar, Madura, Impey, Klein, J Chem @@ -1711,19 +1688,4 @@ Phys, 79, 926 (1983).
(Shinoda) Shinoda, Shiga, and Mikami, Phys Rev B, 69, 134103 (2004).
- - -(Eike and Maginn) Eike and Maginn, J Chem Phys, 124, -164503 (2006). -
- - -(Jayaraman and Maginn) Jayaraman and Maginn, Journal of Chemical Physics, -127, 214504 (2007). -
- - -(Jayaraman) Jayaraman, Thompson, von Lilienfeld and Maginn, Industrial -and Engineering Chemistry Research, 49, 559-571 (2010). -