Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userTime step rescaling recovers continuous-time dynamical properties for discrete-time Langevin integration of nonequilibrium systems
220
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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 properties. However, for stochastic equations of motion (e.g., Langevin dynamics), there is still broad disagreement over which integration algorithms are most appropriate. While multiple desiderata have been proposed throughout the literature, consensus on which criteria are important is absent, and no published integration scheme satisfies all desiderata simultaneously. Additional nontrivial complications stem from simulating systems driven out of equilibrium using existing stochastic integration schemes in conjunction with recently-developed nonequilibrium fluctuation theorems. Here, we examine a family of discrete time integration schemes for Langevin dynamics, assessing how each member satisfies a variety of desiderata that have been enumerated in prior efforts to construct suitable Langevin integrators. We show that the incorporation of a novel time step rescaling in the deterministic updates of position and velocity can correct a number of dynamical defects in these integrators. Finally, we identify a particular splitting that has essentially universally appropriate properties for the simulation of Langevin dynamics for molecular systems in equilibrium, nonequilibrium, and path sampling contexts.
rate research
Read More
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 sampling features, namely the Boltzmann distribution and Einstein diffusion, in linear systems. This revision, which is based on the introduction of time scaling along with flexibility of a discrete-time velocity attenuation parameter, provides a direct link between the ABO splitting formalism and the GJ methods. This link brings about the conclusion that any GJ method has at least weak second order accuracy in the applied time step. It further helps identify a novel half-step velocity, which simultaneously produces both correct kinetic statistics and correct transport measures for any of the statistically sound GJ methods. Explicit algorithmic expressions are given for the integration of the new half-step velocity into the GJ set of methods. Numerical simulations, including quantum-based molecular dynamics (QMD) using the QMD suite LATTE, highlight the discussed properties of the algorithms as well as exhibit the direct application of robust, time step independent stochastic integrators to quantum-based molecular dynamics.
We construct concrete examples of time operators for both continuous and discrete-time homogeneous quantum walks, and we determine their deficiency indices and spectra. For a discrete-time quantum walk, the time operator can be self-adjoint if the time evolution operator has a non-zero winding number. In this case, its spectrum becomes a discrete set of real numbers.
We present results on the ballistic and diffusive behavior of the Langevin dynamics in a periodic potential that is driven away from equilibrium by a space-time periodic driving force, extending some of the results obtained by Collet and Martinez. In the hyperbolic scaling, a nontrivial average velocity can be observed even if the external forcing vanishes in average. More surprisingly, an average velocity in the direction opposite to the forcing may develop at the linear response level -- a phenomenon called negative mobility. The diffusive limit of the non-equilibrium Langevin dynamics is also studied using the general methodology of central limit theorems for additive functionals of Markov processes. To apply this methodology, which is based on the study of appropriate Poisson equations, we extend recent results on pointwise estimates of the resolvent of the generator associated with the Langevin dynamics. Our theoretical results are illustrated by numerical simulations of a two-dimensional system.
Imaginary-time time-dependent Density functional theory (it-TDDFT) has been proposed as an alternative method for obtaining the ground state within density functional theory (DFT) which avoids some of the difficulties with convergence encountered by the self-consistent-field (SCF) iterative method. It-TDDFT was previously applied to clusters of atoms where it was demonstrated to converge in select cases where SCF had difficulty with convergence. In the present work we implement it-TDDFT propagation for {it periodic systems} by modifying the Quantum ESPRESSO package, which uses a plane-wave basis with multiple $boldsymbol{k}$ points, and has the options of non-collinear and DFT+U calculations using ultra-soft or norm-conserving pseudo potentials. We demonstrate that our implementation of it-TDDFT propagation with multiple $boldsymbol{k}$ points is correct for DFT+U non-collinear calculations and for DFT+U calculations with ultra-soft pseudo potentials. Our implementation of it-TDDFT propagation converges to the exact SCF energy (up to the decimal guaranteed by double precision) in all but one case where it converged to a slightly lower value than SCF, suggesting a useful alternative for systems where SCF has difficulty to reach the Kohn-Sham ground state. In addition, we demonstrate that rapid convergence can be achieved if we use adaptive-size imaginary-time-steps for different kinetic-energy plane-waves.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
