Monte Carlo study of real time dynamics

TL;DR

The article tackles the real-time sign problem in lattice Monte Carlo by deforming the integration domain into a complex manifold using holomorphic gradient flow, then sampling via the contraction algorithm. The method yields real-time correlators for the quantum anharmonic oscillator that agree with exact diagonalization and can reach longer real-time separations than some alternative approaches, albeit with slow convergence. It demonstrates a generic, albeit computationally demanding, framework for real-time dynamics that could extend to quantum field theory with improved proposals and efficient Jacobian estimation. Overall, this work advances stochastic methods for real-time transport coefficients and opens a path toward more scalable real-time quantum simulations.

Abstract

Monte Carlo studies involving real time dynamics are severely restricted by the sign problem that emerges from highly oscillatory phase of the path integral. In this letter, we present a new method to compute real time quantities on the lattice using the Schwinger-Keldysh formalism via Monte Carlo simulations. The key idea is to deform the path integration domain to a complex manifold where the phase oscillations are mild and the sign problem is manageable. We use the previously introduced "contraction algorithm" to create a Markov chain on this alternative manifold. We substantiate our approach by analyzing the quantum mechanical anharmonic oscillator. Our results are in agreement with the exact ones obtained by diagonalization of the Hamiltonian. The method we introduce is generic and in principle applicable to quantum field theory albeit very slow. We discuss some possible improvements that should speed up the algorithm.

Figures (3)

• Figure 1: The Schwinger-Keldysh contour (left) and its discretized form (right) in the complex time plane. $\Delta t_n$ refers to either $\pm a$ or $-ia$ depending on the location of $n$ on the contour.
• Figure 2: Time ordered (Feynman) correlators $\langle T(\dot x(t) \dot x(t^\prime))\rangle$ and $\langle T(x(t) x(t^\prime))\rangle$. The dotted and solid lines represent the exact results obtained by diagonalizing the Hamiltonian.
• Figure 3: Histogram of $\mathbb{S}_I$ (mod $2\pi$) for the $T_\text{flow}=0$ calculation corresponding to an integration over ${\mathds R}^N$ (in red) and the $T_\text{flow}=0.2$ calculation corresponding to an integration over $\Gamma$ (in blue). It is clear that the modest flow $T_\text{flow}=0.2$ reduces the sign problem significantly.