Get Force Field Chemistry essential facts below. View Videos or join the Force Field Chemistry discussion. Add Force Field Chemistry to your PopFlock.com topic list for future reference or share this resource on social media.
Force Field Chemistry
A force field is used to minimize the bond stretching energy of this ethane molecule.
All-atom force fields provide parameters for every type of atom in a system, including hydrogen, while united-atom interatomic potentials treat the hydrogen and carbon atoms in each methyl group (terminal methyl) and each methylene bridge as one interaction center. Coarse-grained potentials, which are often used in long-time simulations of macromolecules such as proteins, nucleic acids, and multi-component complexes, provide even cruder representations for higher computing efficiency.
Molecular mechanics potential energy function with continuum solvent.
The basic functional form of potential energy in molecular mechanics includes bonded terms for interactions of atoms that are linked by covalent bonds, and nonbonded (also termed noncovalent) terms that describe the long-range electrostatic and van der Waals forces. The specific decomposition of the terms depends on the force field, but a general form for the total energy in an additive force field can be written as
where the components of the covalent and noncovalent contributions are given by the following summations:
The bond and angle terms are usually modeled by quadratic energy functions that do not allow bond breaking. A more realistic description of a covalent bond at higher stretching is provided by the more expensive Morse potential. The functional form for dihedral energy is highly variable. Additional, "improper torsional" terms may be added to enforce the planarity of aromatic rings and other conjugated systems, and "cross-terms" that describe coupling of different internal variables, such as angles and bond lengths. Some force fields also include explicit terms for hydrogen bonds.
The nonbonded terms are most computationally intensive. A popular choice is to limit interactions to pairwise energies. The van der Waals term is usually computed with a Lennard-Jones potential and the electrostatic term with Coulomb's law, although both can be buffered or scaled by a constant factor to account for electronic polarizability.
As it is rare for bonds to deviate significantly from their reference values the Morse potential is seldom employed for molecular mechanics due to it not being efficient to compute. The most simplistic approaches utilize a Hooke's law formula:
Where is the force constant, is the bond length and is the value for the bond length when all other terms in the force field are set to 0. The term is often referred to as the equilibrium bond length which may cause confusion. The equilibrium bond length would be the value adopted in a minimum energy structure with all other terms contributing.
In most cases the effects on accuracy of modelling a bond as a harmonic oscillator are small however a real bond stretching potential is not harmonic and modelling it as such may lead to inaccuracies in predictions of bond lengths to the thousandth of an angstrom. The strong nature of the interactions between atoms yields large values for the force constant (the stronger the bond the higher the value of the force constant).
Though the Hooke's law formula provides a reasonable level of accuracy at bond lengths near the equilibrium distance it is less accurate as one moves away. In order to model the Morse curve better one can employ cubic and higher powers.
In addition to the functional form of the potentials, force fields define a set of parameters for different types of atoms, chemical bonds, dihedral angles and so on. The parameter sets are usually empirical. A force field would include distinct parameters for an oxygen atom in a carbonylfunctional group and in a hydroxyl group. The typical parameter set includes values for atomic mass, van der Waals radius, and partial charge for individual atoms, and equilibrium values of bond lengths, bond angles, and dihedral angles for pairs, triplets, and quadruplets of bonded atoms, and values corresponding to the effective spring constant for each potential. Most current force fields parameters use a fixed-charge model by which each atom is assigned one value for the atomic charge that is not affected by the local electrostatic environment; proposed developments in next-generation force fields incorporate models for polarizability, in which a particle's charge is influenced by electrostatic interactions with its neighbors. For example, polarizability can be approximated by the introduction of induced dipoles; it can also be represented by Drude particles, massless, charge-carrying virtual sites attached by a springlike harmonic oscillator potential to each polarizable atom. The introduction of polarizability into force fields in common use has been inhibited by the high computational expense associated with calculating the local electrostatic field.
Parameter sets and functional forms are defined by interatomic potentials developers to be self-consistent. Because the functional forms of the potential terms vary extensively between even closely related interatomic potentials (or successive versions of the same interatomic potential), the parameters from one interatomic potential function should clearly never be used together with another interatomic potential function.
There are many approaches to parameterisation of a forcefield. The main classical forcefields are typically paramaterized through a complex process often using approximate quantum mechanics as a basis and modifying potentials to most accurately match experimental observables.  More recently a variety of automated tools have been provided to parameterize new force fields and to assist users to develop their own parameter sets for chemistries which are not parameterized to date.  Efforts are underway within the community to provide open source codes and methods such as openMM and openMD.
All interatomic potentials are based on many approximations and derived from different types of experimental data. Thus, they are termed empirical. Most force fields rely on point charges in reproducing the electrostatic potential around molecules. This leads to substantial failures for anisotropic charge distributions. Some existing energy functions do not account for electronic polarization of the environment, an effect that can significantly reduce electrostatic interactions of partial atomic charges. This problem was addressed by developing polarizable force fields or using macroscopic dielectric constant. However, application of one value of dielectric constant is questionable in the highly heterogeneous environments of proteins or biological membranes, and the nature of the dielectric depends on the model used.
All types of van der Waals forces are also strongly environment-dependent, because these forces originate from interactions of induced and "instantaneous" dipoles (see Intermolecular force). The original Fritz London theory of these forces can only be applied in vacuum. A more general theory of van der Waals forces in condensed media was developed by A. D. McLachlan in 1963 (this theory includes the original London's approach as a special case). The McLachlan theory predicts that van der Waals attractions in media are weaker than in vacuum and follow the like dissolves like rule, which means that different types of atoms interact more weakly than identical types of atoms. This is in contrast to combinatorial rules or Slater-Kirkwood equation applied for development of the classical force fields. The combinatorial rules state that interaction energy of two dissimilar atoms (e.g., C...N) is an average of the interaction energies of corresponding identical atom pairs (i.e., C...C and N...N). According to McLachlan theory, the interactions of particles in a media can even be fully repulsive, as observed for liquid helium. The conclusions of McLachlan theory are supported by direct measurements of attraction forces between different materials (Hamaker constant), as explained by Jacob Israelachvili in his book Intermolecular and surface forces. It was concluded that "the interaction between hydrocarbons across water is about 10% of that across vacuum". Such effects are unaccounted in standard molecular mechanics.
Another round of criticism came from practical applications, such as protein structure refinement. It was noted that Critical Assessment of protein Structure Prediction (CASP) participants did not try to refine their models to avoid "a central embarrassment of molecular mechanics, namely that energy minimization or molecular dynamics generally leads to a model that is less like the experimental structure". The force fields have been applied successfully for protein structure refinement in different X-ray crystallography and NMR spectroscopy applications, especially using program XPLOR. However, such refinement is driven mainly by a set of experimental constraints, whereas the interatomic potentials serve merely to remove interatomic hindrances. The results of calculations are practically the same with rigid sphere potentials implemented in program DYANA (calculations from NMR data), or with programs for crystallographic refinement that do not use any energy functions. The deficiencies of the interatomic potentials remain a major bottleneck in homology modeling of proteins. Such situation gave rise to development of alternative empirical scoring functions specifically for ligand docking,protein folding, homology model refinement, computational protein design, and modeling of proteins in membranes.
There is also an opinion that molecular mechanics may operate with energy which is irrelevant to protein folding or ligand binding. The parameters of typical force fields reproduce enthalpy of sublimation, i.e., energy of evaporation of molecular crystals. However, it was recognized that protein folding and ligand binding are thermodynamically very similar to crystallization, or liquid-solid transitions, because all these processes represent freezing of mobile molecules in condensed media. Thus, free energy changes during protein folding or ligand binding are expected to represent a combination of an energy similar to heat of fusion (energy absorbed during melting of molecular crystals), a conformational entropy contribution, and solvation free energy. The heat of fusion is significantly smaller than enthalpy of sublimation. Hence, the potentials describing protein folding or ligand binding must be weaker than potentials in molecular mechanics. Indeed, the energies of H-bonds in proteins are ~ -1.5 kcal/mol when estimated from protein engineering or alpha helix to coil transition data, but the same energies estimated from sublimationenthalpy of molecular crystals were -4 to -6 kcal/mol. The depths of modified Lennard-Jones potentials derived from protein engineering data were also smaller than in typical potential parameters and followed the like dissolves like rule, as predicted by McLachlan theory.
Popular force fields
Different force fields are designed for different purposes. All are implemented in various computer software.
MM2 was developed by Norman Allinger mainly for conformational analysis of hydrocarbons and other small organic molecules. It is designed to reproduce the equilibrium covalent geometry of molecules as precisely as possible. It implements a large set of parameters that is continuously refined and updated for many different classes of organic compounds (MM3 and MM4).
CFF was developed by Arieh Warshel, Lifson and coworkers as a general method for unifying studies of energies, structures and vibration of general molecules and molecular crystals. The CFF program, developed by Levitt and Warshel, is based on the Cartesian representation of all the atoms, and it served as the basis for many subsequent simulation programs.
ECEPP was developed specifically for modeling of peptides and proteins. It uses fixed geometries of amino acid residues to simplify the potential energy surface. Thus, the energy minimization is conducted in the space of protein torsion angles. Both MM2 and ECEPP include potentials for H-bonds and torsion potentials for describing rotations around single bonds. ECEPP/3 was implemented (with some modifications) in Internal Coordinate Mechanics and FANTOM.
AMBER, CHARMM, and GROMOS have been developed mainly for molecular dynamics of macromolecules, although they are also commonly used for energy minimizing. Thus, the coordinates of all atoms are considered as free variables.
Classical force fields
Assisted Model Building and Energy Refinement (AMBER) - widely used for proteins and DNA.
Chemistry at HARvard Molecular Mechanics (CHARMM) - originally developed at Harvard, widely used for both small molecules and macromolecules
CVFF - also used broadly for small molecules and macromolecules.
COSMOS-NMR - hybrid QM/MM force field adapted to a variety of inorganic compounds, organic compounds and biological macromolecules, including semi-empirical calculation of atomic charges and NMR properties. COSMOS-NMR is optimized for NMR based structure elucidation and implemented in COSMOS molecular modelling package.
GROningen MOlecular Simulation (GROMOS) - a force field that comes as part of the GROMOS software, a general-purpose molecular dynamics computer simulation package for the study of biomolecular systems. GROMOS force field A-version has been developed for application to aqueous or apolar solutions of proteins, nucleotides, and sugars. A B-version to simulate gas phase isolated molecules is also available.
Optimized Potential for Liquid Simulations (OPLS, variants include OPLS-AA, OPLS-UA, OPLS-2001, OPLS-2005) - developed by William L. Jorgensen at the Yale University Department of Chemistry.
ECEPP - first force field for polypeptide molecules - developed by F.A. Momany, H.A. Scheraga and colleagues.
QCFF/PI - A general force fields for conjugated molecules.
Universal Force Field (UFF) - A general force field with parameters for the full periodic table up to and including the actinoids, developed at Colorado State University.
Consistent Force Field (CFF) - a family of forcefields adapted to a broad variety of organic compounds, includes force fields for polymers, metals, etc.
CFF/ind and ENZYMIX - The first polarizable force field which has subsequently been used in many applications to biological systems.
DRF90 developed by P. Th. van Duijnen and coworkers.
PIPF - The polarizable intermolecular potential for fluids is an induced point-dipole force field for organic liquids and biopolymers. The molecular polarization is based on Thole's interacting dipole (TID) model and was developed by Jiali Gao  at the University of Minnesota.
Polarizable Force Field (PFF) - developed by Richard A. Friesner and coworkers.
SP-basis Chemical Potential Equalization (CPE) - approach developed by R. Chelli and P. Procacci.
PHAST - polarizable potential developed by Chris Cioce and coworkers.
CHARMM - polarizable force field developed by S. Patel (University of Delaware) and C. L. Brooks III (University of Michigan).
AMBER - polarizable force field developed by Jim Caldwell and coworkers.
CHARMM - polarizable force field based on the classical Drude oscillator developed by A. MacKerell (University of Maryland, Baltimore) and B. Roux (University of Chicago).
Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) - force field developed by Pengyu Ren (University of Texas at Austin) and Jay W. Ponder (Washington University). AMOEBA force field is gradually moving to more physics-rich AMOEBA+. 
ORIENT - procedure developed by Anthony J. Stone (Cambridge University) and coworkers.
Non-Empirical Molecular Orbital (NEMO) - procedure developed by Gunnar Karlström and coworkers at Lund University (Sweden)
Gaussian Electrostatic Model (GEM) - a polarizable force field based on Density Fitting developed by Thomas A. Darden and G. Andrés Cisneros at NIEHS; and Jean-Philip Piquemal at Paris VI University.
Polarizable procedure based on the Kim-Gordon approach developed by Jürg Hutter and coworkers (University of Zürich)
Computer Simulation of Molecular Structure (COSMOS-NMR) - developed by Ulrich Sternberg and coworkers. Hybrid QM/MM force field enables explicit quantum-mechanical calculation of electrostatic properties using localized bond orbitals with fast BPT formalism. Atomic charge fluctuation is possible in each molecular dynamics step.
Atomistic Polarizable Potential for Liquids, Electrolytes, and Polymers(APPLE&P), developed by Oleg Borogin, Dmitry Bedrov and coworkers, which is distributed by Wasatch Molecular Incorporated.
Reactive force fields
ReaxFF - reactive force field (interatomic potential) developed by Adri van Duin, William Goddard and coworkers. It is fast, transferable and is the computational method of choice for atomistic-scale dynamical simulations of chemical reactions. Parallelized ReaxFF allows reactive simulations on >>1,000,000 atoms.
Empirical valence bond (EVB) - this reactive force field, introduced by Warshel and coworkers, is probably the most reliable and physically consistent way to use force fields in modeling chemical reactions in different environments.[according to whom?] The EVB facilitates calculating activation free energies in condensed phases and in enzymes.
Coarse-grained force fields
Virtual atom molecular mechanics (VAMM) - a coarse-grained force field developed by Korkut and Hendrickson for molecular mechanics calculations such as large scale conformational transitions based on the virtual interactions of C-alpha atoms. It is a knowledge based force field and formulated to capture features dependent on secondary structure and on residue-specific contact information in proteins.
MARTINI - a coarse-grained potential developed by Marrink and coworkers at the University of Groningen, initially developed for molecular dynamics simulations of lipids, later extended to various other molecules. The force field applies a mapping of four heavy atoms to one CG interaction site and is parameterized with the aim of reproducing thermodynamic properties.
SIRAH - a coarse-grained force field developed by Pantano and coworkers of the Biomolecular Simulations Group, Institut Pasteur of Montevideo, Uruguay; developed for molecular dynamics of water, DNA and proteins. Free available for AMBER and GROMACS packages.
Dissipative particle dynamics (DPD) - This is a method commonly applied in chemical engineering. It is typically used for studying the hydrodynamics of various simple and complex fluids which require consideration of time and length scales larger than those accessible to classical Molecular dynamics. The potential was originally proposed by Hoogerbrugge and Koelman  with later modifications by Español and Warren  The current state of the art was well documented in a CECAM workshop in 2008. Recently, work has been undertaken to capture some of the chemical subtitles relevant to solutions. This has led to work considering automated parameterisation of the DPD interaction potentials against experimental observables. 
Machine learning models
ANI is a transferable neural network potential, built from atomic environment vectors, and able to provide DFT accuracy in terms of energies. 
FFLUX (originally QCTFF)  A set of trained Kriging models which operate together to provide a molecular force field trained on Atoms in molecules or Quantum chemical topology energy terms including electrostatic, exchange and electron correlation.
TensorMol a mixed model, a Neural network provides a short range potential, whilst more traditional potentials add screened long range terms.
?-ML not a force field method but a model that adds learnt correctional energy terms to approximate and relatively computationally cheap quantum chemical methods in order to provide an accuracy level of a higher order, more computationally expensive quantum chemical model.
SchNet a Neural network utilising continuous-filter convolutional layers, to predict chemical properties and potential energy surfaces.
The set of parameters used to model water or aqueous solutions (basically a force field for water) is called a water model. Water has attracted a great deal of attention due to its unusual properties and its importance as a solvent. Many water models have been proposed; some examples are TIP3P, TIP4P, SPC, flexible simple point charge water model (flexible SPC), ST2, and mW. Other solvents and methods of solvent representation are also applied within computational chemistry and physics some examples are given on page Solvent model. Recently, novel methods for generating water models have been published.
Post-translational modifications and unnatural amino acids
Forcefield_PTM - An AMBER-based forcefield and webtool for modeling common post-translational modifications of amino acids in proteins developed by Chris Floudas and coworkers. It uses the ff03 charge model and has several side-chain torsion corrections parameterized to match the quantum chemical rotational surface.
Forcefield_NCAA - An AMBER-based forcefield and webtool for modeling common non-natural amino acids in proteins in condensed-phase simulations using the ff03 charge model. The charges have been reported to be correlated with hydration free energies of corresponding side-chain analogs.
Ligand Field Molecular Mechanics (LFMM) - functions for the coordination sphere around transition metals based on the angular overlap model (AOM). Implemented in the Molecular Operating Environment (MOE) as DommiMOE and in Tinker
^Leach, Dr Andrew (2001-01-30). Molecular Modelling: Principles and Applications (2 ed.). Harlow: Prentice Hall. ISBN9780582382107.
^Siu; et al. (2012). "Optimization of the OPLS-AA force field for long hydrocarbons". Journal of Chemical Theory and Computation. 8 (4): 1459-1470. doi:10.1021/ct200908r. PMID26596756.
^Aduri; et al. (2007). "AMBER force field parameters for the naturally occurring modified nucleosides in RNA". Journal of Chemical Theory and Computation. 3 (4): 1464-1475. doi:10.1021/ct600329w. PMID26633217.
^kirschner; et al. (2012). "A glycam-based force field for simulations of lipopolysaccharide membranes: parametrization and validation". Journal of Chemical Theory and Computation. 8 (11): 4719-4731. doi:10.1021/ct300534j. PMID26605626.
^Wang; et al. (2014). "Building force fields: An automatic, systematic, and reproducible approach". The Journal of Physical Chemistry Letters. 5 (11): 1885-1891. doi:10.1021/jz500737m. PMID26273869.
^McDonagh J. L.; et al. (2019). "Utilizing machine learning for efficient parameterization of coarse grained molecular force fields". Journal of Chemical Information and Modeling. 59 (10): 4278-4288. doi:10.1021/acs.jcim.9b00646. PMID31549507.
^Kramer, Christian; Spinn, Alexander; Liedl, Klaus R. (2014). "Charge Anisotropy: Where Atomic Multipoles Matter Most". Journal of Chemical Theory and Computation. 10 (10): 4488-4496. doi:10.1021/ct5005565. PMID26588145.
^ abWarshel, Arieh; Sharma, Pankaz K.; Kato, Mitsunori; Parson, William W. (2006). "Modeling electrostatic effects in proteins". Biochimica et Biophysica Acta (BBA) - Proteins and Proteomics. 1764 (11): 1647-1676. doi:10.1016/j.bbapap.2006.08.007. PMID17049320.
^Schutz, Claudia N.; Warshel, Arieh (2001). "What are the dielectric constants of proteins and how to validate electrostatic models?". Proteins: Structure, Function, and Genetics. 44 (4): 400-417. doi:10.1002/prot.1106. PMID11484218.
^ abcdIsraelachvili, J.N. (1992). Intermolecular and surface forces. Academic Press, San Diego.
^ abLomize A.L., Reibarkh M.Y. and Pogozheva I.D. (2002). "Interatomic potentials and solvation parameters from protein engineering data for buried residues". Protein Sci.11 1984-2000.
^Murphy K.P. and Gill S.J. (1991). "Solid model compounds and the thermodynamics of protein unfolding". J. Mol. Biol.222 699-709.
^Shakhnovich, E.I. and Finkelstein, A.V. (1989). "Theory of cooperative transitions in protein molecules. I. Why denaturation of globular proteins is a first-order phase transition". Biopolymers28 1667-1680.
^Graziano, G., Catanzano, F., Del Vecchio, P., Giancola, C., and Barone, G. (1996). "Thermodynamic stability of globular proteins: a reliable model from small molecule studies". Gazetta Chim. Italiana126 559-567.
^Myers J.K. and Pace C.N. (1996). "Hydrogen bonding stabilizes globular proteins". Biophys. J.71 2033-2039.
^Scholtz J.M., Marqusee S., Baldwin R.L., York E.J., Stewart J.M., Santoro M., and Bolen D.W. (1991). "Calorimetric determination of the enthalpy change for the alpha-helix to coil transition of an alanine peptide in water". Proc. Natl. Acad. Sci. USA88 2854-2858.
^Gavezotti A. and Filippini G. (1994). "Geometry of intermolecular X-H...Y (X,Y=N,O) hydrogen bond and the calibration of empirical hydrogen-bond potentials". J. Phys. Chem.98 4831-4837.
^Allinger, Norman L. (1977). "Conformational analysis. 130. MM2. A hydrocarbon interatomic potential utilizing V1 and V2 torsional terms". J. Am. Chem. Soc.99 8127-8134.
^Allinger, N.L., Yuh, Y.H., & Lii, J-H. (1989). "Molecular Mechanics. The MM3 Force Field for Hydrocarbons. 1". J. Am. Chem. Soc.111 8551-8565.
Lii, J-H., & Allinger, N.L. (1989). "Molecular Mechanics. The MM3 Force Field for Hydrocarbons. 2. Vibrational Frequencies and Thermodynamics". J. Am. Chem. Soc.111 8566-8575.
Lii, J-H., & Allinger, N.L. (1989). "Molecular Mechanics. The MM3 Force Field for Hydrocarbons. 3. The van der Waals Potentials and Crystal data for Aliphatic and Aromatic Hydrocarbons". J. Am. Chem. Soc.111 8576-8582.
^Schaumann, T., Braun, W. and Wutrich, K. (1990). "The program FANTOM for energy refinement of polypeptides and proteins using a Newton-Raphson minimizer in torsion angle space". Biopolymers29 679-694.
^Möllhoff, Margit and Sternberg, Ulrich (2001). "Molecular Mechanics with fluctuating atomic charges - a new force field with semi-empirical charge calculation". J. Mol. Model.7 90-102.
^F. A. Momany, R. F. McGuire, A. W. Burgess, and H. A. Scheraga (1975). "Energy parameters in polypeptides. VII." J. Phys. Chem.79 2361.
^Yelena A. Arnautova , Anna Jagielska, and Harold A. Scheraga. (2006). "A New Force Field (ECEPP-05) for Peptides, Proteins, and Organic Molecules" J. Phys. Chem. B110 5025.
^Warshel, A. (1973). "Quantum Mechanical Consistent Force Field (QCFF/PI) Method: Calculations of Energies, Conformations and Vibronic Interactions of Ground and Excited States of Conjugated Molecules" Israel J. Chem.11 709.
^Warshel, A. and Levitt, M. (1974). "QCFF/PI: A Program for the Consistent Force Field Evaluation of Equilibrium Geometries and Vibrational Frequencies of Molecules". QCPE 247, Quantum Chemistry Program Exchange, Indiana University.
^Warshel A. and Levitt M. (1976) "Theoretical Studies of Enzymatic Reactions: Dielectric Electrostatic and Steric Stabilization of the Carbonium Ion in the Reaction of Lysozyme". J. Mol. Biol.103 227-249.
^Gao, J., Habibollahzadeh, D., and Shao, L. (1995) "A Polarizable Intermolecular Potential Functions for Simulations of Liquid Alcohols". J. Phys. Chem.99 16460-16467.
^Xie, W., Pu, J., MacKerell, A. D., Jr., and Gao, J. (2007) "Development of a Polarizable Intermolecular Potential Function (PIPF) for Liquid Amides and Alkanes". J. Chem. Theory Comput.3 1878-1889.
^Chelli, R., Procacci, P. (2002) "A Transferable Polarizable Electrostatic Force Field for Molecular Mechanics based on the Chemical Potential Equalization Principle." J. Chem. Phys. 17 9175-9189
^Cioce, C; McLaughlin, K; Belof, J; Space, B (2013). "A Polarizable and Transferable PHAST N2 Potential for Use in Materials Simulation". J. Chem. Theory Comput. 9 (12): 5550-5557. doi:10.1021/ct400526a. PMID26592288.
^Patel, S., Brooks, C.L., III (2004) "CHARMM fluctuating charge force field for proteins: I parameterization and application to bulk organic liquid simulations." J. Comput. Chem.25 1-16.
^Patel, S., MacKerell, A.D., Jr., Brooks, C.L., III (2004) "CHARMM fluctuating charge force field for proteins: II Protein/solvent properties from molecular dynamics simulations using a nonadditive electrostatic model." J. Comput. Chem.25 1504-1514.
^Anisimov, V.M., Lamoureux, G., Vorobyov, I.V., Huang, N., Roux, B. and MacKerell, A.D., Jr. (2005) "Determination of Electrostatic Parameters for a Polarizable Force Field Based on the Classical Drude Oscillator". J. Chem. Theory Comput.1 153-168.
^Yu, H., Whitfield, T.W., Harder, E., Lamoureux, G., Vorobyov, I., Anisimov, V. M., MacKerell, A.D., Jr., and Roux, B. (2010) "Simulating Monovalent and Divalent Ions in Aqueous Solution Using a Drude Polarizable Force Field". J. Chem. Theory Comput.6 774-786.
^Liu, C., Piquemal, J-P, and Ren, P. (2020) "Implementation of Geometry-Dependent Charge Flux into the Polarizable AMOEBA+ Potential." J. Phys. Chem. Lett.11 419-426 doi.org/10.1021/acs.jpclett.9b03489
^Liu, C., Piquemal, J-P, and Ren, P. (2019) "AMOEBA+ Classical Potential for Modeling Molecular Interactions." J. Chem. Theory Comput.15 4122-4137 doi.org/10.1021/acs.jctc.9b00261
^O. Engkvist, P.-O. Åstrand, and G. Karlström. "Accurate intermolecular potentials obtained from molecular wave functions: Bridging the gap between quantum chemistry and molecular simulations". Chem. Rev., 100:4087-4108, 2000.
^Gresh, N., Cisneros, G. A., Darden, T. A. and Piquemal, J-P. (2007) "Anisotropic, polarizable molecular mechanics studies of inter-, intra-molecular interactions, and ligand-macromolecule complexes. A bottom-up strategy." J. Chem. Theory Comput.3 1960.
^Piquemal, J-P, Cisneros, G.A., Reinhardt, P., Gresh, N. and Darden, T.A. (2006). "Towards a Force Field based on Density Fitting".
J. Chem. Phys.124 104101.
^Cisneros, G.A., Piquemal, J-P and Darden, T.A. (2006). "Generalization of the Gaussian Electrostatic Model: extension to arbitrary angular momentum, distributed multipoles and speedup with reciprocal space methods". J. Chem. Phys.125 184101.
^Sternberg, Ulrich, Koch, Frank-Thomas and Möllhoff, Margit (1994). "New Approach to the Semiempirical Calculation of Atomic Charges for Polypeptides and Large Molecular Systems" J. Comput. Chem.15 524-531.
^McDonagh J. L.; et al. (2019). "Utilizing machine learning for efficient parameterization of coarse grained molecular force fields". Journal of Chemical Information and Modeling. 59 (10): 4278-4288. doi:10.1021/acs.jcim.9b00646. PMID31549507.
^Fletcher; et al. (2016). "Multipolar electrostatic energy prediction for all 20 natural amino acids using Kriging machine learning". Journal of Chemical Theory and Computation. 12 (6): 2742-2751. doi:10.1021/acs.jctc.6b00457. PMID27224739.
^McDonagh; et al. (2017). "Machine learning of dynamic electron correlation energies from topological atoms". Journal of Chemical Theory and Computation. 14 (1): 216-224. doi:10.1021/acs.jctc.7b01157. PMID29211469.
^Khoury, G. A.; Bhatia, N.; Floudas, C. A. (2014). "Hydration free energies calculated using the AMBER ff03 charge model for natural and unnatural amino acids and multiple water models". Computers & Chemical Engineering. 71: 745-752. doi:10.1016/j.compchemeng.2014.07.017.
^Deeth, Robert J. (2001). "The ligand field molecular mechanics model and the stereoelectronic effects of d and s electrons". Coordination Chemistry Reviews. 212 (212): 11-34. doi:10.1016/S0010-8545(00)00354-4.
^Foscato, Marco; Deeth, Robert J.; Jensen, Vidar R. (2015). "Integration of Ligand Field Molecular Mechanics in Tinker". Journal of Chemical Information and Modeling. 55 (6): 1282-1290. doi:10.1021/acs.jcim.5b00098. PMID25970002.
Israelachvili, J. N. (1992). Intermolecular and surface forces. San Diego: Academic Press. ISBN978-0-12-375181-2.
Schlick, T. (2002). Molecular Modeling and Simulation: An Interdisciplinary Guide. Interdisciplinary Applied Mathematics: Mathematical Biology. New York: Springer-Verlag. ISBN978-0-387-95404-2.
Warshel, A. (1991). Computer Modeling of Chemical Reactions in Enzymes and Solutions. New York: John Wiley & Sons. ISBN978-0-471-53395-5.