A Correlation Method for Deriving UA Intra-Molecular Potentials from AA Molecular Dynamics Simulations: Application to Alkanes

23 Jan 2020  ·  Monasse Bernard, Boussinot Frédéric ·

United-Atom (UA) force fields are usually constructed using a Boltzmann-inverse method based on distributions obtained from Monte-Carlo simulations. A new method of constructing UA force fields from All-Atom (AA) molecular dynamics simulations is proposed. In this method, one determines time correlations between oscillators of the same type: between a CC bond and the connected CH bonds; between a valence angle CCC and all valence angles sharing at least one hydrogen with it; between a dihedral angle CCCC and the dihedral angles sharing the two central carbons with it. In the case of no correlation between lengths or angles, energies of the oscillators are independent. In this case, the AA and UA components of the force fields are identical. When a correlation is total, the UA potential of an object is the sum of the energies of the oscillators coupled to it. Partial correlations are also possible. Several kinematic matching tests of the AA molecule with the corresponding UA molecule are considered: mean periods of vibration of the oscillators; energies associated with presence probabilities; mean end-to-end lengths of the molecules. The method is, as an example, applied to aliphatic molecules, simulated with the OPLS-AA force field.

PDF Abstract
No code implementations yet. Submit your code now

Categories


Chemical Physics Statistical Mechanics