Saved in:
Bibliographic Details
Main Author: Toxvaerd, Søren
Format: Preprint
Published: 2024
Subjects:
Online Access:https://arxiv.org/abs/2403.01612
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866916145671438336
author Toxvaerd, Søren
author_facet Toxvaerd, Søren
contents Almost all Molecular Dynamics (MD) simulations are discrete dynamics with Newton's algorithm first published in 1687, and much later by L. Verlet in 1967. Discrete Newtonian dynamics has the same qualities as Newton's classical analytic dynamics. Verlet also published a first-order expression for the instant temperature which is inaccurate but presumably used in most MD simulations. One of the motivations for the present article is to correct this unnecessary inaccuracy in $NVT$ MD dynamics. Another motivation is to derive simple algorithms for the Nosé-Hoover $NVT$ dynamics (NH) with the correct temperature constraint. The simulations with NH discrete Newtonian dynamics show that the NH works excellent for a wide range of the response time $τ$ of the NH thermostat, but NH simulations favor a choice of a short response time, even shorter than the discrete-time increment $δt$ used in MD, to avoid large oscillations of the temperature.
format Preprint
id arxiv_https___arxiv_org_abs_2403_01612
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Discrete Newtonian dynamics with Nosé-Hoover thermostats
Toxvaerd, Søren
Statistical Mechanics
Almost all Molecular Dynamics (MD) simulations are discrete dynamics with Newton's algorithm first published in 1687, and much later by L. Verlet in 1967. Discrete Newtonian dynamics has the same qualities as Newton's classical analytic dynamics. Verlet also published a first-order expression for the instant temperature which is inaccurate but presumably used in most MD simulations. One of the motivations for the present article is to correct this unnecessary inaccuracy in $NVT$ MD dynamics. Another motivation is to derive simple algorithms for the Nosé-Hoover $NVT$ dynamics (NH) with the correct temperature constraint. The simulations with NH discrete Newtonian dynamics show that the NH works excellent for a wide range of the response time $τ$ of the NH thermostat, but NH simulations favor a choice of a short response time, even shorter than the discrete-time increment $δt$ used in MD, to avoid large oscillations of the temperature.
title Discrete Newtonian dynamics with Nosé-Hoover thermostats
topic Statistical Mechanics
url https://arxiv.org/abs/2403.01612