Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA symplectic integrator for molecular dynamics on a hypersphere
107
0
0.0
(
0
)
Added by
Jean-Michel Caillol
Publication date
2020
fields
Physics
and research's language is
English
Authors
J.-M. Caillol
Ask ChatGPT about the research
No Arabic abstract
We present a reversible and symplectic algorithm called ROLL, for integrating the equations of motion in molecular dynamics simulations of simple fluids on a hypersphere $mathcal{S}^d$ of arbitrary dimension $d$. It is derived in the framework of geometric algebra and shown to be mathematically equivalent to algorithm RATTLE. An application to molecular dynamics simulation of the one component plasma is briefly discussed.
rate research
Read More
Here we present a program aimed at free-energy calculations in molecular systems. It consists of a series of routines that can be interfaced with the most popular classical molecular dynamics (MD) codes through a simple patching procedure. This leaves the possibility for the user to exploit many different MD engines depending on the system simulated and on the computational resources available. Free-energy calculations can be performed as a function of many collective variables, with a particular focus on biological problems, and using state-of-the-art methods such as metadynamics, umbrella sampling and Jarzynski-equation based steered MD. The present software, written in ANSI-C language, can be easily interfaced with both fortran and C/C++ codes.
Deep learning is transforming many areas in science, and it has great potential in modeling molecular systems. However, unlike the mature deployment of deep learning in computer vision and natural language processing, its development in molecular modeling and simulations is still at an early stage, largely because the inductive biases of molecules are completely different from those of images or texts. Footed on these differences, we first reviewed the limitations of traditional deep learning models from the perspective of molecular physics, and wrapped up some relevant technical advancement at the interface between molecular modeling and deep learning. We do not focus merely on the ever more complex neural network models, instead, we emphasize the theories and ideas behind modern deep learning. We hope that transacting these ideas into molecular modeling will create new opportunities. For this purpose, we summarized several representative applications, ranging from supervised to unsupervised and reinforcement learning, and discussed their connections with the emerging trends in deep learning. Finally, we outlook promising directions which may help address the existing issues in the current framework of deep molecular modeling.
This paper had no abstract originally. A second-order symplectic integration algorithm for guiding center motion is presented. The algorithm is based on the Poincare (mid-point) generating function.
To take into account nuclear quantum effects on the dynamics of atoms, the path integral molecular dynamics (PIMD) method used since 1980s is based on the formalism developed by R. P. Feynman. However, the huge computation time required for the PIMD reduces its range of applicability. Another drawback is the requirement of additional techniques to access time correlation functions (ring polymer MD or centroid MD). We developed an alternative technique based on a quantum thermal bath (QTB) which reduces the computation time by a factor of ~20. The QTB approach consists in a classical Langevin dynamics in which the white noise random force is replaced by a Gaussian random force having the power spectral density given by the quantum fluctuation-dissipation theorem. The method has yielded satisfactory results for weakly anharmonic systems: the quantum harmonic oscillator, the heat capacity of a MgO crystal, and isotope effects in 7 LiH and 7 LiD. Unfortunately, the QTB is subject to the problem of zero-point energy leakage (ZPEL) in highly anharmonic systems, which is inherent in the use of classical mechanics. Indeed, a part of the energy of the high-frequency modes is transferred to the low-frequency modes leading to a wrong energy distribution. We have shown that in order to reduce or even eliminate ZPEL, it is sufficient to increase the value of the frictional coefficient. Another way to solve the ZPEL problem is to combine the QTB and PIMD techniques. It requires the modification of the power spectral density of the random force within the QTB. This combination can also be seen as a way to speed up the PIMD.
Kinetic energy equipartition is a premise for many deterministic and stochastic molecular dynamics methods that aim at sampling a canonical ensemble. While this is expected for real systems, discretization errors introduced by the numerical integration may distort such assumption. Fortunately, backward error analysis allows us to identify the quantity that is actually subject to equipartition. This is related to a shadow Hamiltonian, which coincides with the specified Hamiltonian only when the time-step size approaches zero. This paper deals with discretization effects in a straightforward way. With a small computational overhead, we obtain refine
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
