Bringing discrete-time Langevin splitting methods into agreement with thermodynamics.

In light of the recently published complete set of statistically correct Grønbech-Jensen (GJ) methods for discrete-time thermodynamics, we revise a 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 Los Alamos Transferable Tight-Binding for Energetics, highlight the discussed properties of the algorithms as well as exhibit the direct application of robust, time-step-independent stochastic integrators to QMD.

Dissipation-Dependent Thermal Escape from a Potential Well.

Langevin simulations are conducted to investigate the Josephson escape statistics over a large set of parameter values for damping and temperature. The results are compared to both Kramers and Büttiker-Harris-Landauer (BHL) models, and good agreement is found with the Kramers model for high to moderate damping, while the BHL model provides further good agreement down to lower damping values. However, for extremely low damping, even the BHL model fails to reproduce the progression of the escape statistics. In order to explain this discrepancy, we develop a new model which shows that the bias sweep effectively cools the system below the thermodynamic value as the potential well broadens due to the increasing bias. A simple expression for the temperature is derived, and the model is validated against direct Langevin simulations for extremely low damping values.

The challenge of stochastic Størmer-Verlet thermostats generating correct statistics.

In light of the recently developed complete GJ set of single random variable stochastic, discrete-time Størmer-Verlet algorithms for statistically accurate simulations of Langevin equations [N. Grønbech-Jensen, Mol. Phys. 118, e1662506 (2020)], we investigate two outstanding questions: (1) Are there any algorithmic or statistical benefits from including multiple random variables per time step and (2) are there objective reasons for using one or more methods from the available set of statistically correct algorithms? To address the first question, we assume a general form for the discrete-time equations with two random variables and then follow the systematic, brute-force GJ methodology by enforcing correct thermodynamics in linear systems. It is concluded that correct configurational Boltzmann sampling of a particle in a harmonic potential implies correct configurational free-particle diffusion and that these requirements only can be accomplished if the two random variables per time step are identical. We consequently submit that the GJ set represents all possible stochastic Størmer-Verlet methods that can reproduce time step-independent statistics of linear systems. The second question is thus addressed within the GJ set. Based on numerical simulations of complex molecular systems, as well as on analytic considerations, we analyze apparent friction-induced differences in the stability of the methods. We attribute these differences to an inherent, friction-dependent discrete-time scaling, which depends on the specific method. We suggest that the method with the simplest interpretation of temporal scaling, the GJ-I/GJF-2GJ method, be preferred for statistical applications.

Nonequilibrium transient phenomena in the washboard potential.

Transient properties of the one-dimensional washboard potential are investigated in order to understand observed modulations in the statistics of escape events. Specifically, we analyze the effects of different kinds of initial conditions on the escape distribution obtained by linearly increasing the tilt of the potential. Despite the complexity of the dynamics leading up to the eventual escape, we find that the overall statistics can be interpreted in terms of the system parameters, which offers illuminating perspectives for driven one-dimensional systems with washboard potentials. We choose parameters sets relevant for Josephson junctions, a commonly studied system due to both its applications and its use as a model system in condensed matter physics.