This article was accepted into the corpus but its outbound wikilinks were never NER-processed — typical at the deepest BFS hop or when the run's entity cap was reached. No expansion funnel to show.
| SPC/E | |
|---|---|
| Name | SPC/E |
| Type | Molecular model |
| Developer | B. Hess (originated by H. J. C. Berendsen and collaborators) |
| Introduced | 1987 |
| Applications | Molecular dynamics, Monte Carlo methods, Biomolecular simulation, Liquid state theory |
| Components | Water molecule (three-site rigid) |
SPC/E
SPC/E is a rigid three-site water model widely used in Molecular dynamics and Monte Carlo methods simulations of condensed-phase water and aqueous systems. Designed to reproduce bulk thermodynamic and dielectric properties with modest computational cost, the model became a standard in studies of protein folding, lipid bilayers, electrolyte solutions, and interfacial phenomena. It was introduced as an extension of the SPC family by researchers associated with Erasmus University Rotterdam and refinements from groups including Berendsen group and later popularized by implementations in software such as GROMACS, LAMMPS, and AMBER.
The model represents each water molecule as three interaction sites corresponding to the two hydrogen nuclei and the oxygen nucleus arranged in a rigid geometry constrained by bond lengths and bond angle. Interactions combine pairwise Coulombic electrostatics located on the sites and a Lennard-Jones potential centered on the oxygen site to capture repulsion and dispersion; this is analogous to earlier models like SPC and TIP3P. The geometry typically uses an O–H bond length and H–O–H angle chosen to match experimental gas-phase values or effective condensed-phase geometry used in parametrization. SPC/E employs constraints (e.g., SHAKE or SETTLE algorithms) that are implemented across many molecular simulation packages including CHARMM-compatible engines.
Parameterization of the model adjusted point charges and Lennard-Jones parameters to improve agreement with experimental properties relative to SPC. The effective partial charges on oxygen and hydrogen were tuned alongside the oxygen σ and ε parameters to reproduce densities and radial distribution features of liquid water at ambient conditions. An empirical self-polarization correction is incorporated as an average polarization energy term, introduced to account for many-body polarization effects absent in rigid nonpolarizable models; this correction distinguishes SPC/E from non-corrected variants and aligns it conceptually with polarizable efforts such as POL3 and TIP4P/2005f. SPC/E parameters are commonly combined with force fields for biomolecules developed by groups at University of Groningen and elsewhere.
SPC/E reproduces many bulk thermodynamic quantities with reasonable accuracy at standard conditions: liquid density near 1.0 g·cm−3, cohesion reflected in enthalpy of vaporization when the polarization correction is included, and a liquid structure captured by oxygen–oxygen, oxygen–hydrogen, and hydrogen–hydrogen radial distribution functions. Structural comparisons are routinely made to neutron scattering and X-ray diffraction experiments and to ab initio molecular dynamics results from groups using Density functional theory approximations such as PBE or dispersion-corrected functionals. SPC/E yields pronounced hydrogen-bonding networks manifested in coordination numbers and angular distributions, often compared against data from Infrared spectroscopy and Raman spectroscopy studies.
Transport properties like self-diffusion coefficient and viscosity predicted by SPC/E often deviate from experiment but can be corrected by finite-size scaling and thermostat-barostat choices in simulations. The model typically underestimates viscosity and overestimates diffusion relative to experimental values at ambient temperature; adjustments in thermostating (e.g., Nosé–Hoover chains) and long-range electrostatics methods such as Ewald summation or Particle Mesh Ewald influence reported values. Rotational correlation times and dielectric relaxation computed with SPC/E are compared against dielectric spectroscopy and NMR relaxation measurements, with the model capturing qualitative trends in temperature dependence but showing quantitative discrepancies, motivating further model development like TIP4P/2005.
SPC/E has been employed across a wide range of investigations, including solvation free energies in ion hydration studies, conformational sampling in protein and nucleic acid simulations, membrane hydration around lipid bilayers and micelles, interfacial water at oxide and metal surfaces, and transport in nanoporous materials such as zeolites and carbon nanotubes. Its computational efficiency and availability in major packages (GROMACS, LAMMPS, AMBER, NAMD) make it a default choice for benchmarking and large-scale ensembles, often used in tandem with enhanced-sampling methods like umbrella sampling and metadynamics.
Despite strengths, SPC/E has limitations: the rigid, nonpolarizable representation neglects explicit electronic polarization and nuclear quantum effects that are important for accurate dielectric properties, proton transfer, and isotope effects. Compared to extended models such as TIP4P/2005, TIP5P, and explicit-polarization models like AMOEBA and Drude oscillators, SPC/E is less accurate for phase equilibria (melting point, vapor-liquid coexistence curves), interfacial surface tension, and high-pressure behavior. Its empirical polarization correction improves enthalpy predictions but does not substitute for many-body polarization present in models like MB-pol developed for quantitative spectroscopy and cluster energetics. Users often choose SPC/E for efficiency and historical comparability, while selecting newer models for high-precision studies in physical chemistry and materials science.
Category:Water models