Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Journal of High Energy Physics

Real-time quantum dynamics, path integrals and the method of thimbles

Zong-Gang Mou, Paul M. Saffin, Anders Tranberg B, Simon Woodward

Published Jun 1, 2019• 34 citations
Publisher

Abstract

Abstract Direct numerical evaluation of the real-time path integral has a well-known sign problem that makes convergence exponentially slow. One promising remedy is to use Picard-Lefschetz theory to flow the domain of the field variables into the complex plane, where the integral is better behaved. By Cauchy’s theorem, the final value of the path integral is unchanged. Previous analyses have considered the case of real scalar fields in thermal equilibrium, employing a closed Schwinger-Keldysh time contour, allowing the evaluation of the full quantum correlation functions. Here we extend the analysis by not requiring a closed time path, instead allowing for an initial density matrix for out-of-equilibrium initial value problems. We are able to explicitly implement Gaussian initial conditions, and by separating the initial time and the later times into a two-step Monte-Carlo sampling, we are able to avoid the phenomenon of multiple thimbles. In fact, there exists one and only one thimble for each sample member of the initial density matrix. We demonstrate the approach through explicitly computing the real-time propagator for an interacting scalar in 0+1 dimensions, and find very good convergence allowing for comparison with perturbation theory and the classical-statistical approximation to real-time dynamics.

DOI
10.1007/jhep06(2019)094

Related papers

Real-time quantum dynamics, path integrals and the method of thimbles

Zong-Gang Mou, Paul M. Saffin, Anders Tranberg B +1 more

Direct numerical evaluation of the real-time path integral has a well-known sign problem that makes convergence exponentially slow. One promising remedy is to use Picard-Lefschetz theory to flow the domain of the field variables into the complex plane, where the integral is better behaved. By Cauchy's theorem, the final value of the path integral is unchanged. Previous analyses have considered the case of real scalar fields in thermal equilibrium, employing a closed Schwinger-Keldysh time contour, allowing the evaluation of the full quantum correlation functions. Here we extend the analysis by not requiring a closed time path, instead allowing for an initial density matrix for out-of-equilibrium initial value problems. We are able to explicitly implement Gaussian initial conditions, and by separating the initial time and the later times into a two-step Monte-Carlo sampling, we are able to avoid the phenomenon of multiple thimbles. In fact, there exists one and only one thimble for each sample member of the initial density matrix. We demonstrate the approach through explicitly computing the real-time propagator for an interacting scalar in 0+1 dimensions, and find very good convergence allowing for comparison with perturbation theory and the classical-statistical approximation to real-time dynamics.

arXiv:1902.09147•Feb 25, 2019•2 citations

Quantum tunnelling, real-time dynamics and Picard-Lefschetz thimbles

Zong-Gang Mou, Paul M. Saffin, Anders Tranberg B

Abstract We follow up the work, where in light of the Picard-Lefschetz thimble approach, we split up the real-time path integral into two parts: the initial density matrix part which can be represented via an ensemble of initial conditions, and the dynamic part of the path integral which corresponds to the integration over field variables at all later times. This turns the path integral into a two-stage problem where, for each initial condition, there exits one and only one critical point and hence a single thimble in the complex space, whose existence and uniqueness are guaranteed by the characteristics of the initial value problem. In this paper, we test the method for a fully quantum mechanical phenomenon, quantum tunnelling in quantum mechanics. We compare the method to solving the Schrödinger equation numerically, and to the classical-statistical approximation, which emerges naturally in a well-defined limit. We find that the Picard-Lefschetz result matches the expectation from quantum mechanics and that, for this application, the classical-statistical approximation does not.

arXiv:1909.02488•Journal of High Energy Physics•Nov 1, 2019•29 citations

Computing real-time quantum path integrals on Sewed, almost-Lefschetz thimbles

Zong-Gang Mou, Paul M. Saffin, Anders Tranberg B

A bstract We present a method to compute real-time path integrals numerically, by Monte-Carlo sampling on near-Lefschetz thimbles. We present a collection of tools based on the Lefschetz thimble methods, which together provide an alternative to existing methods such as the Generalised thimble. These involve a convenient coordinate parameterization of the thimble, direct numerical integration along a radial coordinate into an effective path integral weight and locally deforming the Lefschetz thimble using its Gaussian (non-interacting theory) counterpart in a region about the critical point. We apply this to quantum mechanics, identify possible pitfalls and benefits, and benchmark its efficiency.

Journal of High Energy Physics•Mar 7, 2025•1 citation

Lefschetz thimble quantum Monte Carlo for spin systems

T. C. Mooney, Jacob Bringewatt, Neill C. Warrington +1 more

Monte Carlo simulations are useful tools for modeling quantum systems, but in some cases they suffer from a sign problem, leading to an exponential slow down in their convergence to a value. While solving the sign problem is generically NP hard, many techniques exist for mitigating the sign problem in specific cases; in particular, the technique of deforming the Monte Carlo simulation's plane of integration onto Lefschetz thimbles (complex hypersurfaces of stationary phase) has seen significant success in the context of quantum field theories. We extend this methodology to spin systems by utilizing spin coherent state path integrals to reexpress the spin system's partition function in terms of continuous variables. Using some toy systems, we demonstrate its effectiveness at lessening the sign problem in this setting, despite the fact that the initial mapping to spin coherent states introduces its own sign problem. The standard formulation of the spin coherent path integral is known to make use of uncontrolled approximations; despite this, for large spins they are typically considered to yield accurate results, so it is somewhat surprising that our results show significant systematic errors. Therefore, possibly of independent interest, our use of Lefschetz thimbles to overcome the intrinsic sign problem in spin coherent state path integral Monte Carlo enables a novel numerical demonstration of a breakdown in the spin coherent path integral.

arXiv:2110.10699•Physical review. B./Physical review. B•Dec 15, 2022•5 citations

Quantum decoherence in the Caldeira-Leggett model by the real-time path integral on a computer

Jun Nishimura, Hiromasa Watanabe

A bstract We propose first-principle calculations of an open system based on the real-time path integral formalism treating the environment as well as the system of our interest together on a computer. The sign problem that occurs in applying Monte Carlo methods can be overcome in general by using the so-called Lefschetz thimble method, which has been developed over the past decade. Here we focus on the Caldeira-Leggett model, which is well known, in particular, as a model of quantum decoherence. In this case, the calculation simplifies drastically since the path integral becomes Gaussian for typical initial conditions. The relevant saddle point, which is unique and complex, can be determined by solving a linear equation with a huge but sparse coefficient matrix, and the integration over the Lefschetz thimble can be performed analytically. Thus we obtain, without assumptions or approximations, the reduced density matrix after a long-time evolution, tracing out a large number of harmonic oscillators in the environment. In particular, we confirm the dependence of the decoherence time on the coupling constant and the temperature that has been predicted from the master equation in a certain parameter regime.

arXiv:2503.20699•Journal of High Energy Physics•Sep 24, 2025•2 citations

A new picture of quantum tunneling in the real-time path integral from Lefschetz thimble calculations

Jun Nishimura, K. Sakai, Atis Yosprakob

Abstract It is well known that quantum tunneling can be described by instantons in the imaginary-time path integral formalism. However, its description in the real-time path integral formalism has been elusive. Here we establish a statement that quantum tunneling can be characterized in general by the contribution of complex saddle points, which can be identified by using the Picard-Lefschetz theory. We demonstrate this explicitly by performing Monte Carlo simulations of simple quantum mechanical systems, overcoming the sign problem by the generalized Lefschetz thimble method. We confirm numerically that the contribution of complex saddle points manifests itself in a complex “weak value” of the Hermitian coordinate operator $$ \hat{x} $$ x ̂ evaluated at time t, which is a physical quantity that can be measured by experiments in principle. We also discuss the transition to classical dynamics based on our picture.

arXiv:2307.11199•Journal of High Energy Physics•Sep 18, 2023•21 citations
Abstract Direct numerical evaluation of the real-time path integral has a well-known sign problem that makes convergence exponentially slow. One promising remedy is to use Picard-Lefschetz theory to flow the domain of the field variables into the complex plane, where the integral is better behaved. By Cauchy’s theorem, the final value of the path integral is unchanged. Previous analyses have considered the case of real scalar fields in thermal equilibrium, employing a closed Schwinger-Keldysh time contour, allowing the evaluation of the full quantum correlation functions. Here we extend the analysis by not requiring a closed time path, instead allowing for an initial density matrix for out-of-equilibrium initial value problems. We are able to explicitly implement Gaussian initial conditions, and by separating the initial time and the later times into a two-step Monte-Carlo sampling, we are able to avoid the phenomenon of multiple thimbles. In fact, there exists one and only one thimble for each sample member of the initial density matrix. We demonstrate the approach through explicitly computing the real-time propagator for an interacting scalar in 0+1 dimensions, and find very good convergence allowing for comparison with perturbation theory and the classical-statistical approximation to real-time dynamics.
Direct numerical evaluation of the real-time path integral has a well-known sign problem that makes convergence exponentially slow. One promising remedy is to use Picard-Lefschetz theory to flow the domain of the field variables into the complex plane, where the integral is better behaved. By Cauchy's theorem, the final value of the path integral is unchanged. Previous analyses have considered the case of real scalar fields in thermal equilibrium, employing a closed Schwinger-Keldysh time contour, allowing the evaluation of the full quantum correlation functions. Here we extend the analysis by not requiring a closed time path, instead allowing for an initial density matrix for out-of-equilibrium initial value problems. We are able to explicitly implement Gaussian initial conditions, and by separating the initial time and the later times into a two-step Monte-Carlo sampling, we are able to avoid the phenomenon of multiple thimbles. In fact, there exists one and only one thimble for each sample member of the initial density matrix. We demonstrate the approach through explicitly computing the real-time propagator for an interacting scalar in 0+1 dimensions, and find very good convergence allowing for comparison with perturbation theory and the classical-statistical approximation to real-time dynamics.
Abstract We follow up the work, where in light of the Picard-Lefschetz thimble approach, we split up the real-time path integral into two parts: the initial density matrix part which can be represented via an ensemble of initial conditions, and the dynamic part of the path integral which corresponds to the integration over field variables at all later times. This turns the path integral into a two-stage problem where, for each initial condition, there exits one and only one critical point and hence a single thimble in the complex space, whose existence and uniqueness are guaranteed by the characteristics of the initial value problem. In this paper, we test the method for a fully quantum mechanical phenomenon, quantum tunnelling in quantum mechanics. We compare the method to solving the Schrödinger equation numerically, and to the classical-statistical approximation, which emerges naturally in a well-defined limit. We find that the Picard-Lefschetz result matches the expectation from quantum mechanics and that, for this application, the classical-statistical approximation does not.
A bstract We present a method to compute real-time path integrals numerically, by Monte-Carlo sampling on near-Lefschetz thimbles. We present a collection of tools based on the Lefschetz thimble methods, which together provide an alternative to existing methods such as the Generalised thimble. These involve a convenient coordinate parameterization of the thimble, direct numerical integration along a radial coordinate into an effective path integral weight and locally deforming the Lefschetz thimble using its Gaussian (non-interacting theory) counterpart in a region about the critical point. We apply this to quantum mechanics, identify possible pitfalls and benefits, and benchmark its efficiency.
Monte Carlo simulations are useful tools for modeling quantum systems, but in some cases they suffer from a sign problem, leading to an exponential slow down in their convergence to a value. While solving the sign problem is generically NP hard, many techniques exist for mitigating the sign problem in specific cases; in particular, the technique of deforming the Monte Carlo simulation's plane of integration onto Lefschetz thimbles (complex hypersurfaces of stationary phase) has seen significant success in the context of quantum field theories. We extend this methodology to spin systems by utilizing spin coherent state path integrals to reexpress the spin system's partition function in terms of continuous variables. Using some toy systems, we demonstrate its effectiveness at lessening the sign problem in this setting, despite the fact that the initial mapping to spin coherent states introduces its own sign problem. The standard formulation of the spin coherent path integral is known to make use of uncontrolled approximations; despite this, for large spins they are typically considered to yield accurate results, so it is somewhat surprising that our results show significant systematic errors. Therefore, possibly of independent interest, our use of Lefschetz thimbles to overcome the intrinsic sign problem in spin coherent state path integral Monte Carlo enables a novel numerical demonstration of a breakdown in the spin coherent path integral.
A bstract We propose first-principle calculations of an open system based on the real-time path integral formalism treating the environment as well as the system of our interest together on a computer. The sign problem that occurs in applying Monte Carlo methods can be overcome in general by using the so-called Lefschetz thimble method, which has been developed over the past decade. Here we focus on the Caldeira-Leggett model, which is well known, in particular, as a model of quantum decoherence. In this case, the calculation simplifies drastically since the path integral becomes Gaussian for typical initial conditions. The relevant saddle point, which is unique and complex, can be determined by solving a linear equation with a huge but sparse coefficient matrix, and the integration over the Lefschetz thimble can be performed analytically. Thus we obtain, without assumptions or approximations, the reduced density matrix after a long-time evolution, tracing out a large number of harmonic oscillators in the environment. In particular, we confirm the dependence of the decoherence time on the coupling constant and the temperature that has been predicted from the master equation in a certain parameter regime.
Abstract It is well known that quantum tunneling can be described by instantons in the imaginary-time path integral formalism. However, its description in the real-time path integral formalism has been elusive. Here we establish a statement that quantum tunneling can be characterized in general by the contribution of complex saddle points, which can be identified by using the Picard-Lefschetz theory. We demonstrate this explicitly by performing Monte Carlo simulations of simple quantum mechanical systems, overcoming the sign problem by the generalized Lefschetz thimble method. We confirm numerically that the contribution of complex saddle points manifests itself in a complex “weak value” of the Hermitian coordinate operator x^ \hat{x} x^ x ̂ evaluated at time t, which is a physical quantity that can be measured by experiments in principle. We also discuss the transition to classical dynamics based on our picture.