Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe Effects of Intrinsic Dynamical Ghost Modes in Discrete-Time Langevin Simulations
152
0
0.0
(
0
)
Added by
Niels Gronbech-Jensen
Publication date
2019
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Using the recently published GJF-2GJ Langevin thermostat, which can produce time-step-independent statistical measures even for large time steps, we analyze and discuss the causes for abrupt deviations in statistical data as the time step is increased for some simulations of nonlinear oscillators. Exemplified by the pendulum, we identify a couple of discrete-time dynamical modes in the purely damped pendulum equation as the cause of the observed discrepancies in statistics. The existence, stability and kinetics of the modes are consistent with the acquired velocity distribution functions from Langevin simulations, and we conclude that the simulation deviations from physical expectations are not due to normal, systematic algorithmic time-step errors, but instead due to the inherent properties of discrete time in nonlinear dynamics.
rate research
Read More
We present the complete set of stochastic Verlet-type algorithms that can provide correct statistical measures for both configurational and kinetic sampling in discrete-time Langevin systems. The approach is a brute-force general representation of the Verlet-algorithm with free parameter coefficients that are determined by requiring correct Boltzmann sampling for linear systems, regardless of time step. The result is a set of statistically correct methods given by one free functional parameter, which can be interpreted as the one-time-step velocity attenuation factor. We define the statistical characteristics of both true on-site $v^n$ and true half-step $u^{n+frac{1}{2}}$ velocities, and use these definitions for each statistically correct Stormer-Verlet method to find a unique associated half-step velocity expression, which yields correct kinetic Maxwell-Boltzmann statistics for linear systems. It is shown that no other similar, statistically correct on-site velocity exists. We further discuss the use and features of finite-difference velocity definitions that are neither true on-site, nor true half-step. The set of methods is written in convenient and conventional stochastic Verlet forms that lend themselves to direct implementation for, e.g., Molecular Dynamics applications. We highlight a few specific examples, and validate the algorithms through comprehensive Langevin simulations of both simple nonlinear oscillators and complex Molecular Dynamics.
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 stochastic thermostat, which has been shown to provide correct statistical configurational sampling for any time step that yields stable trajectories. Here, we extend the method and develop a set of discrete-time equations of motion for both particle dynamics and system volume in order to seek pressure control that is insensitive to the choice of the numerical time step. The resulting method is simple, practical, and efficient. The method is demonstrated through direct numerical simulations of two characteristic model systems - a one dimensional particle chain for which exact statistical results can be obtained and used as benchmarks, and a three dimensional system of Lennard-Jones interacting particles simulated in both solid and liquid phases. The results, which are compared against the method of Kolb & Dunweg, show that the new method behaves according to the objective, namely that acquired statistical averages and fluctuations of configurational measures are accurate and robust against the chosen time step applied to the simulation.
Following the recent realisation that periodically driven quantum matter can support new types of spatiotemporal order, now known as discrete time crystals (DTCs), we consider the stability of this phenomenon. Motivated by its conceptual importance as well as its experimental relevance we consider the effect of coupling to an external environment. We use this to argue, both analytically and numerically, that the DTC in disordered one-dimensional systems is destroyed at long times by any such natural coupling. This holds true even in the case where the coupling is such that the system is prevented from heating up by an external thermal bath.
How much time does it take two molecules to react? If a reaction occurs upon contact, the answer to this question boils down to the classic first-passage time problem: find the random time it takes the two molecules to meet. However, this is not always the case as molecules switch stochastically between reactive and non-reactive states. In such cases, the reaction is said to be ``gated by the internal states of the molecules involved which could have a dramatic influence on kinetics. A unified, continuous-time, approach to gated reactions on networks was presented in [Phys. Rev. Lett. 127, 018301, 2021]. Here, we build on this recent advancement and develop an analogous discrete-time version of the theory. Similar to continuous-time, we employ a renewal approach to show that the gated reaction time can always be expressed in terms of the corresponding ungated first-passage and return times; which yields formulas for the generating function of the gated reaction-time distribution and its corresponding mean and variance. In cases where the mean reaction time diverges, we show that the long-time asymptotics of the gated problem is inherited from its ungated counterpart, where only the pre-factor of the power law tail changes. The discretization of time also gives rise to new phenomena that do not exist in the continuous-time analogue. Crucially, when the internal gating dynamics is in, or out of, phase with the spatial process governing molecular encounters resonance and anti-resonance phenomena emerge. These phenomena are illustrated using two case studies which also serve to show how the general approach presented herein greatly simplifies the analysis of gated reactions.
The emph{ab initio} path integral Monte Carlo (PIMC) approach is one of the most successful methods in quantum many-body theory. A particular strength of this method is its straightforward access to imaginary-time correlation functions (ITCF). For example, the well-known density-density ITCF $F(mathbf{q},tau)$ allows one to estimate the linear response of a given system for all wave vectors $mathbf{q}$ from a single simulation of the unperturbed system. Moreover, it constitutes the basis for the reconstruction of the dynamic structure factor $S(mathbf{q},omega)$ -- a key quantity in state-of-the-art scattering experiments. In this work, we present analogous relations between the nonlinear density response in quadratic and cubic order of the perturbation strength and generalized ITCFs measuring correlations between up to four imaginary-time arguments. As a practical demonstration of our new approach, we carry out simulations of the warm dense electron gas and find excellent agreement with previous PIMC results that had been obtained with substantially larger computational effort. In addition, we give a relation between a cubic ITCF and the triple dynamic structure factor $S(mathbf{q}_1,omega_1;mathbf{q}_2,omega_2)$, which evokes the enticing possibility to study dynamic three-body effects on an emph{ab initio} level.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
