On the validity of complex Langevin method for path integral computations

Abstract

The complex Langevin (CL) method is a classical numerical strategy to alleviate the numerical sign problem in the computation of lattice field theories. Mathematically, it is a simple numerical tool to compute a wide class of high-dimensional and oscillatory integrals. However, it is often observed that the CL method converges but the limiting result is incorrect. The literature has several unclear or even conflicting statements, making the method look mysterious. By an in-depth analysis of a model problem, we reveal the mechanism of how the CL result turns biased as the parameter changes, and it is demonstrated that such a transition is difficult to capture. Our analysis also shows that the method works for any observables only if the probability density function generated by the CL process is localized. To generalize such observations to lattice field theories, we formulate the CL method on general groups using rigorous mathematical languages for the first time, and we demonstrate that such localized probability density function does not exist in the simulation of lattice field theories for general compact groups, which explains the unstable behavior of the CL method. Fortunately, we also find that the gauge cooling technique creates additional velocity that helps confine the samples, so that we can still see localized probability density functions in certain cases, as significantly broadens the application of the CL method. The limitations of gauge cooling are also discussed. In particular, we prove that gauge cooling has no effect for Abelian groups, and we provide an example showing that biased results still exist when gauge cooling is insufficient to confine the probability density function.

arXiv comment28 pages,9 figures

References in corpus (24)

• High density QCD on a Lefschetz thimble?
• The Complex Langevin method: When can it be trusted?
• Taming the dynamical sign problem in real-time evolution of quantum many-body problems
• Simulating full QCD at nonzero density using the complex Langevin equation
• Monte Carlo simulations on the Lefschetz thimble: taming the sign problem
• Gauge cooling in complex Langevin for QCD with heavy quarks
• Controlling complex Langevin dynamics at finite density
• The argument for justification of the complex Langevin method and the condition for correct convergence
• A new approach to the complex-action problem and its application to a nonperturbative study of superstring theory
• Real-Time Dynamics in Quantum Impurity Models with Diagrammatic Monte Carlo
• Complex Langevin analysis of the space-time structure in the Lorentzian type IIB matrix model
• Stability of complex Langevin dynamics in effective models
• Simulating QCD at nonzero baryon density to all orders in the hopping parameter expansion
• Complex Langevin and boundary terms
• Localised distributions and criteria for correctness in complex Langevin dynamics
• Status of Complex Langevin
• Controlling Complex Langevin simulations of lattice models by boundary term analysis
• Does the complex Langevin method give unbiased results?
• Applying Complex Langevin Simulations to Lattice QCD at Finite Density
• Real-time gauge theory simulations from stochastic quantization using optimized updating
• Complex Langevin analysis of the spontaneous symmetry breaking in dimensionally reduced super Yang-Mills models
• Thermodynamics of spin-orbit coupled bosons in two dimensions from complex Langevin
• Complex Langevin: Boundary terms at poles
• Beyond Complex Langevin Equations: positive representation of a class of complex measures

Cited by in corpus (1)

• Complex Langevin study for polarons in an attractively interacting one-dimensional two-component Fermi gas