Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userImproving dynamical properties of metropolized discretizations of overdamped Langevin dynamics
105
0
0.0
(
0
)
Authors
Max Fathi
Ask ChatGPT about the research
No Arabic abstract
The discretization of overdamped Langevin dynamics, through schemes such as the Euler-Maruyama method, can be corrected by some acceptance/rejection rule, based on a Metropolis-Hastings criterion for instance. In this case, the invariant measure sampled by the Markov chain is exactly the Boltzmann-Gibbs measure. However, rejections perturb the dynamical consistency of the resulting numerical method with the reference dynamics. We present in this work some modifications of the standard correction of discretizations of overdamped Langevin dynamics on compact spaces by a Metropolis-Hastings procedure, which allow us to either improve the strong order of the numerical method, or to decrease the bias in the estimation of transport coefficients characterizing the effective dynamical behavior of the dynamics. For the latter approach, we rely on modified numerical schemes together with a Barker rule for the acceptance/rejection criterion.
rate research
Read More
We consider in this work the numerical computation of transport coefficients for Brownian dynamics. We investigate the discretization error arising when simulating the dynamics with the Smart MC algorithm (also known as Metropolis-adjusted Langevin algorithm). We prove that the error is of order one in the time step, when using either the Green-Kubo or the Einstein formula to estimate the transport coefficients. We illustrate our results with numerical simulations.
When an external field drives a colloidal system out of equilibrium, the ensuing colloidal response can be very complex and obtaining a detailed physical understanding often requires case-by-case considerations. In order to facilitate systematic analysis, here we present a general iterative scheme for the determination of the unique external force field that yields a prescribed inhomogeneous stationary or time-dependent flow in an overdamped Brownian many-body system. The computer simulation method is based on the exact one-body force balance equation and allows to specifically tailor both gradient and rotational velocity contributions, as well as to freely control the one-body density distribution. Hence compressibility of the flow field can be fully adjusted. The practical convergence to a unique external force field demonstrates the existence of a functional map from both velocity and density to external force field, as predicted by the power functional variational framework. In equilibrium, the method allows to find the conservative force field that generates a prescribed target density profile, and hence implements the Mermin-Evans classical density functional map from density distribution to external potential. The conceptual tools developed here enable one to gain detailed physical insight into complex flow behaviour, as we demonstrate in prototypical situations.
We show how to derive a simple integrator for the Langevin equation and illustrate how it is possible to check the accuracy of the obtained distribution on the fly, using the concept of effective energy introduced in a recent paper [J. Chem. Phys. 126, 014101 (2007)]. Our integrator leads to correct sampling also in the difficult high-friction limit. We also show how these ideas can be applied in practical simulations, using a Lennard-Jones crystal as a paradigmatic case.
In this work, we consider non-linear corrections to the Langevin effective theory of a heavy quark moving through a strongly coupled CFT plasma. In AdS/CFT, this system can be identified with that of a string stretched between the boundary and the horizon of an asymptotically AdS black-brane solution. We compute the Feynman-Vernon influence phase for the heavy quark by evaluating the Nambu-Goto action on a doubled string configuration. This configuration is the linearised solution of the string motion in the doubled black-brane geometry which has been proposed as the holographic dual of a thermal Schwinger-Keldysh contour of the CFT. Our expression for the influence phase passes non-trivial consistency conditions arising from the underlying unitarity and thermality of the bath. The local effective theory obeys the recently proposed non-linear fluctuation dissipation theorem relating the non-Gaussianity of thermal noise to the thermal jitter in the damping constant. This furnishes a non-trivial check for the validity of these relations derived in the weak coupling regime.
We study the exponential convergence to the stationary state for nonequilibrium Langevin dynamics, by a perturbative approach based on hypocoercive techniques developed for equilibrium Langevin dynamics. The Hamiltonian and overdamped limits (corresponding respectively to frictions going to zero or infinity) are carefully investigated. In particular, the maximal magnitude of admissible perturbations are quantified as a function of the friction. Numerical results based on a Galerkin discretization of the generator of the dynamics confirm the theoretical lower bounds on the spectral gap.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
