Extending the regime of linear response with synthetic forcings

03/06/2023
by   Renato Spacek, et al.
0

Transport coefficients, such as the mobility, thermal conductivity and shear viscosity, are quantities of prime interest in statistical physics. At the macroscopic level, transport coefficients relate an external forcing of magnitude η, with η≪ 1, acting on the system to an average response expressed through some steady-state flux. In practice, steady-state averages involved in the linear response are computed as time averages over a realization of some stochastic differential equation. Variance reduction techniques are of paramount interest in this context, as the linear response is scaled by a factor of 1/η, leading to large statistical error. One way to limit the increase in the variance is to allow for larger values of η by increasing the range of values of the forcing for which the nonlinear part of the response is sufficiently small. In theory, one can add an extra forcing to the physical perturbation of the system, called synthetic forcing, as long as this extra forcing preserves the invariant measure of the reference system. The aim is to find synthetic perturbations allowing to reduce the nonlinear part of the response as much as possible. We present a mathematical framework for quantifying the quality of synthetic forcings, in the context of linear response theory, and discuss various possible choices for them. Our findings are illustrated with numerical results in low-dimensional systems.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
11/19/2022

Error estimates and variance reduction for nonequilibrium stochastic dynamics

Equilibrium properties in statistical physics are obtained by computing ...
research
10/06/2019

Convergence of the likelihood ratio method for linear response of non-equilibrium stationary states

We consider numerical schemes for computing the linear response of stead...
research
05/14/2023

Fixing the flux: A dual approach to computing transport coefficients

We present a method to compute transport coefficients in molecular dynam...
research
11/30/2021

Martingale product estimators for sensitivity analysis in computational statistical physics

We introduce a new class of estimators for the linear response of steady...
research
10/23/2018

Fast Computation of Steady-State Response for Nonlinear Vibrations of High-Degree-of-Freedom Systems

We discuss an integral equation approach that enables fast computation o...
research
05/13/2020

Nonlinear Fokker-Planck Acceleration for Forward-Peaked Transport Problems in Slab Geometry

This paper introduces a nonlinear acceleration technique that accelerate...
research
06/28/2019

Modeling Response Time Distributions with Generalized Beta Prime

We use Generalized Beta Prime distribution, also known as GB2, for fitti...

Please sign up or login with your details

Forgot password? Click here to reset