Data driven inference for the repulsive exponent of the Lennard-Jones potential in molecular dynamics simulations
OPEN ACCESS
Author / Producer
Date
2017
Publication Type
Journal Article
ETH Bibliography
yes
Citations
Altmetric
OPEN ACCESS
Data
Rights / License
Abstract
The Lennard-Jones (LJ) potential is a cornerstone of Molecular Dynamics (MD) simulations and among the most widely used computational kernels in science. The LJ potential models atomistic attraction and repulsion with century old prescribed parameters (q = 6, p = 12, respectively), originally related by a factor of two for simplicity of calculations. We propose the inference of the repulsion exponent through Hierarchical Bayesian uncertainty quantification We use experimental data of the radial distribution function and dimer interaction energies from quantum mechanics simulations. We find that the repulsion exponent p ≈ 6.5 provides an excellent fit for the experimental data of liquid argon, for a range of thermodynamic conditions, as well as for saturated argon vapour. Calibration using the quantum simulation data did not provide a good fit in these cases. However, values p ≈ 12.7 obtained by dimer quantum simulations are preferred for the argon gas while lower values are promoted by experimental data. These results show that the proposed LJ 6-p potential applies to a wider range of thermodynamic conditions, than the classical LJ 6-12 potential. We suggest that calibration of the repulsive exponent in the LJ potential widens the range of applicability and accuracy of MD simulations.
Permanent link
Publication status
published
External links
Editor
Book title
Journal / series
Volume
7
Pages / Article No.
16576
Publisher
Nature
Event
Edition / version
Methods
Software
Geographic location
Date collected
Date created
Subject
Organisational unit
03499 - Koumoutsakos, Petros (ehemalig) / Koumoutsakos, Petros (former)
02803 - Collegium Helveticum / Collegium Helveticum
02803 - Collegium Helveticum / Collegium Helveticum
02803 - Collegium Helveticum / Collegium Helveticum