Related papers: An Efficient Integrator Scheme for Sampling the (Q…
We present a new method of conducting molecular dynamics simulation in isothermal-isobaric ensemble based on Langevin equations of motion. The stochastic coupling to all particle and cell degrees of freedoms is introduced in a correct way,…
We show a unified second-order scheme for constructing simple, robust and accurate algorithms for typical thermostats for configurational sampling for the canonical ensemble. When Langevin dynamics is used, the scheme leads to the BAOAB…
The authors present a new molecular dynamics algorithm for sampling the isothermal-isobaric ensemble. In this approach the velocities of all particles and volume degrees of freedom are rescaled by a properly chosen random factor. The…
A methodology to sample the isothermal-isobaric ensemble using Langevin dynamics is proposed, which combines novel features of geometric integrators for the equations of motion. By employing the Trotter expansion, the methodology generates…
We present a new and improved method for simultaneous control of temperature and pressure in molecular dynamics simulations with periodic boundary conditions. The thermostat-barostat equations are build on our previously developed…
We present generalized-ensemble algorithms for isobaric-isothermal molecular simulations. In addition to the multibaric-multithermal algorithm and replica-exchange method for the isobaric-isothermal ensemble, which have already been…
Heterogeneous multiscale methods (HMM) combine molecular accuracy of particle-based simulations with the computational efficiency of continuum descriptions to model flow in soft matter liquids. In these schemes, molecular simulations…
A recently proposed method for computer simulations in the isothermal-isobaric (NPT) ensemble, based on Langevin-type equations of motion for the particle coordinates and the ``piston'' degree of freedom, is re-derived by straightforward…
Molecular dynamics with the stochastic process provides a convenient way to compute structural and thermodynamic properties of chemical, biological, and materials systems. It is demonstrated that the virtual dynamics case that we proposed…
In this work we introduce a geometric integrator for molecular dynamics simulations of physical systems in the canonical ensemble. In particular, we consider the equations arising from the so-called density dynamics algorithm with any…
It is often difficult to quantitatively determine if a new molecular simulation algorithm or software properly implements sampling of the desired thermodynamic ensemble. We present some simple statistical analysis procedures to allow…
An extended Hamiltonian approach to conduct isothermal-isobaric molecular dynamics simulations with full cell flexibility is presented. The components of the metric tensor are used as the fictitious degrees of freedom for the cell, thus…
Molecules like water have vibrational modes with a zero-point energy well above room temperature. As a consequence, classical molecular dynamics simulations of their liquids largely underestimate the energy of modes with a higher zero-point…
We present a systematic framework for shortening and automating molecular dynamics equilibration through improved position initialization methods and uncertainty quantification analysis, using the Yukawa one-component plasma as an exemplar…
Sampling from flat energy or density distributions has proven useful in equilibrating complex systems with large energy barriers. Several thermostats and barostats are presented to sample these flat distributions by molecular dynamics.…
The integration time step is a critical determinant of performance in molecular dynamics simulations, governing the trade-off between speed and fidelity. Although 2 fs remains the standard in atomistic biomolecular simulations, the push for…
We propose a new generalized-ensemble algorithm, which we refer to as the multibaric-multithermal Monte Carlo method. The multibaric-multithermal Monte Carlo simulations perform random walks widely both in volume space and in potential…
We present a general molecular-dynamics simulation scheme, based on the Nose' thermostat, for sampling according to arbitrary phase space distributions. We formulate numerical methods based on both Nose'-Hoover and Nose'-Poincare'…
We present a novel approach to investigate the long-time stochastic dynamics of multi-dimensional classical systems, in contact with a heat-bath. When the potential energy landscape is rugged, the kinetics displays a decoupling of short and…
We introduce a novel simple algorithm for thermostatting path integral molecular dynamics (PIMD) with the Langevin equation. The staging transformation of path integral beads is employed for demonstration. The optimum friction coefficients…