Abstract

This study explores the potential of modern implicit solvers for stochastic partial differential equations in the simulation of real-time complex Langevin dynamics. Not only do these methods offer asymptotic stability, rendering the issue of runaway solution moot, but they also allow us to simulate at comparatively large Langevin time steps, leading to lower computational cost. We compare different ways of regularizing the underlying path integral and estimate the errors introduced due to the finite Langevin time steps. Based on that insight, we implement benchmark (non-)thermal simulations of the quantum anharmonic oscillator on the canonical Schwinger-Keldysh contour of short real-time extent.

Highlights

  • Are variants of reweighing, extrapolation from complex parameters and the reformulation of the system of interest in new degrees of freedom unaffected by a sign problem

  • This study explores the potential of modern implicit solvers for stochastic partial differential equations in the simulation of real-time complex Langevin dynamics

  • We focus in this study solely on the question of stability, returning to the remaining two in future work. (I.e. in order to remain in the parameter range where the complex Langevin method itself is known to converge to the correct results, we limit ourselves to a short real-time extent in this study.)

Read more

Summary

Complex Langevin and SDE solvers

The task at hand is to compute quantum statistical expectation values of an observable O. Note that the delta function in the correlator of the noise encompasses all dimensions of x. This prescription places an independent stochastic process at each space-time point. Due to the complex drift term, the field degrees complexify and we can rewrite the evolution equations instead in terms of the real and imaginary part of the field as φ(x, τL) = φR(x, τL) + iφI (x, τL),. We have here used the standard construction in which the noise term η(x, τL) is real These are the stochastic partial differential evolution equations we will solve in the subsequent sections. The first central task is to find appropriate numerical solvers to accommodate these equations

Numerical schemes
On the issue of large excursions
Formulating and discretizing the model
Regularizing the model
Tilted Schwinger-Keldysh contour
Regularization via an implicit scheme
Finite Langevin time step errors
Stable CL simulations at short real-times
Dynamics in thermal equilibrium
Non-equilibrium dynamics
Summary and outlook

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.