ﻻ يوجد ملخص باللغة العربية
We expand on the previously published Gr{o}nbech-Jensen Farago (GJF) thermostat, which is a thermodynamically sound variation on the St{o}rmer-Verlet algorithm for simulating discrete-time Langevin equations. The GJF method has been demonstrated to give robust and accurate configurational sampling of the phase space, and its applications to, e.g., Molecular Dynamics is well established. A new definition of the discrete-time velocity variable is proposed based on analytical calculations of the kinetic response of a harmonic oscillator subjected to friction and noise. The new companion velocity to the GJF method is demonstrated to yield correct and time-step-independent kinetic responses for, e.g., kinetic energy, its fluctuations, and Green-Kubo diffusion based on velocity autocorrelations. This observation allows for a new and convenient Leap-Frog algorithm, which efficiently and precisely represents statistical measures of both kinetic and configurational properties at any time step within the stability limit for the harmonic oscillator. We outline the simplicity of the algorithm and demonstrate its attractive time-step-independent features for nonlinear and complex systems through applications to a one-dimensional nonlinear oscillator and three-dimensional Molecular Dynamics.
In light of the recently published complete set of statistically correct GJ methods for discrete-time thermodynamics, we revise the differential operator splitting method for the Langevin equation in order to comply with the basic GJ thermodynamic sa
We show how to derive a simple integrator for the Langevin equation and illustrate how it is possible to check the accuracy of the obtained distribution on the fly, using the concept of effective energy introduced in a recent paper [J. Chem. Phys. 12
When simulating molecular systems using deterministic equations of motion (e.g., Newtonian dynamics), such equations are generally numerically integrated according to a well-developed set of algorithms that share commonly agreed-upon desirable proper
We implement the statistically sound G-JF thermostat for Langevin Dynamics simulations into the ESPREesSo molecular package for large-scale simulations of soft matter systems. The implemented integration method is tested against the integrator curren
A novel formulation of fluid dynamics as a kinetic theory with tailored, on-demand constructed particles removes any restrictions on Mach number and temperature as compared to its predecessors, the lattice Boltzmann methods and their modifications. I