Mie potential

Mie potential
Name	Mie potential
Field	Statistical mechanics, Molecular dynamics, Computational chemistry
Namedafter	Gustav Mie
Relatedconcepts	Lennard-Jones potential, Van der Waals force, Intermolecular force

Mie potential. The Mie potential is a mathematical model used in statistical mechanics and molecular simulation to describe the interaction energy between a pair of neutral atoms or molecules. It generalizes the widely used Lennard-Jones potential by allowing the exponents of the attractive and repulsive terms to be variable, rather than fixed. This flexibility provides a more accurate and tunable description of intermolecular forces for a broader range of materials and conditions, from noble gases to complex polymers.

Definition and mathematical form

The standard mathematical form of the Mie potential between two particles is given by the equation \( U(r) = \frac{C_n}{r^n} - \frac{C_m}{r^m} \), where \( r \) is the distance between particle centers. The positive term represents the repulsive interaction at short ranges, governed by exponent \( n \), while the negative term represents the attractive London dispersion force, governed by exponent \( m \). The constants \( C_n \) and \( C_m \) are related to the depth of the potential well and the equilibrium separation distance. A common reformulation uses parameters for the well depth \( \epsilon \) and the finite distance at which the inter-particle potential is zero, \( \sigma \), leading to \( U(r) = \frac{n}{n-m} \left( \frac{n}{m} \right)^{m/(n-m)} \epsilon \left[ \left( \frac{\sigma}{r} \right)^n - \left( \frac{\sigma}{r} \right)^m \right] \). This form is directly implemented in simulation software like LAMMPS and GROMACS.

Historical background and development

The potential is named after Gustav Mie, a German physicist known for his work on scattering theory and electrodynamics. Mie introduced the general form in a 1903 paper investigating intermolecular potentials, predating the more famous 1924 work of John Lennard-Jones. The development of the potential was part of early 20th-century efforts to derive equations of state, such as the van der Waals equation, from microscopic principles. Its adoption in modern computational chemistry was spurred by the work of researchers like Aneesur Rahman and Loup Verlet, who pioneered molecular dynamics simulations. The quest for more accurate force fields for materials like water and ionic liquids has renewed interest in this classical formulation.

Physical interpretation and parameters

Physically, the repulsive term models the steep increase in energy as electron clouds overlap, a consequence of the Pauli exclusion principle. The attractive term arises primarily from correlated electron fluctuations described by quantum mechanics, known as London dispersion forces. The exponents \( n \) and \( m \) determine the "hardness" of the repulsive core and the range of the attraction, respectively; larger values of \( n \) yield a harder core. Parameters \( \epsilon \) and \( \sigma \) are typically fitted to experimental data such as second virial coefficients, viscosity, or results from ab initio quantum chemistry methods like those developed at the University of Chicago. For argon, a common reference system, typical values are \( n=14, m=7 \).

Applications in molecular simulation

The Mie potential is extensively used in coarse-grained modeling of soft matter systems, including lipid bilayers, block copolymers, and surfactant micelles. Its adjustable exponents allow for the precise tuning of phase behavior and critical phenomena, which is crucial for studying polymer blends and colloidal suspensions. In engineering, it is applied within the SAFT (Statistical Associating Fluid Theory) equation of state to model the thermodynamic properties of complex fluids like refrigerants and hydrocarbon mixtures. Research institutions like the National Institute of Standards and Technology (NIST) utilize it for developing accurate property tables for industrial design.

Comparison with other intermolecular potentials

Compared to the fixed-exponent Lennard-Jones potential (where n=12, m=6), the Mie potential offers greater flexibility, often leading to improved agreement with experimental data for properties like diffusion coefficients and radial distribution functions. The Buckingham potential includes an exponential repulsion, which can be more realistic but computationally expensive. The purely repulsive Weeks-Chandler-Andersen (WCA) potential is a truncated and shifted version used for studying hard-sphere systems. For charged systems, potentials like the Coulomb potential or the Born-Mayer potential are combined with the Mie form to model ionic crystals or molten salts.

Limitations and extensions

A primary limitation is its inability to describe specific chemical interactions such as hydrogen bonding, metal coordination, or covalent bond formation, which require additional terms or different functional forms like the Morse potential. It also assumes pairwise additivity, neglecting many-body effects important in polarizable media like water. Extensions include the Mie \( \lambda \)-form used in the TraPPE force field, and combination with dipolar or quadrupolar terms for modeling carbon dioxide or benzene. Recent work integrates it with machine learning potentials developed by groups like the Fritz Haber Institute to bridge accuracy and computational cost for materials discovery.

Category:Interatomic potentials Category:Statistical mechanics Category:Computational chemistry