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

arXiv:1909.02488 · Journal of High Energy Physics

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

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

Published Nov 1, 2019• 29 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1007/jhep11(2019)135

Related papers

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

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

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.

Journal of High Energy Physics•Jun 1, 2019•34 citations

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

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

Quantum tunneling in the real-time path integral by the Lefschetz thimble method

Jun Nishimura

Quantum tunneling is mostly discussed in the Euclidean path integral formalism using instantons. On the other hand, it is difficult to understand quantum tunneling based on the real-time path integral due to its oscillatory nature, which causes the notorious sign problem. We show that recent development of the Lefschetz thimble method enables us to investigate this issue numerically. In particular, we find that quantum tunneling occurs due to complex trajectories, which are actually observable experimentally by using the so-called weak measurement.

arXiv:2308.00345•Corfu Summer Institute School and Workshops on Elementary Particle Physics and Gravity (CORFU)•Nov 8, 2023•1 citation

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 Thimbles and Wigner Functions

Vladimir S. Filinov, Alexander S. Larkin

The major difficulty for the Feynman Path Integral Monte Carlo (PIMC) simulations of the quantum systems of particles is the so called “sign problem”, arising due to the fast oscillations of the path integral integrand depending on the complex-valued action. Our aim is to find universal techniques being able to solve this problem. The new method combines the basic ideas of the Metropolis and Hasting algorithms and is based on the Picard-Lefschetz theory and complex-valued version of Morse theory. The basic idea is to choose the Lefschetz thimbles as manifolds approaching the saddle point of the integrand. On this thimble the imaginary part of the complex-valued action remains constant. As a result the integrand on each thimble does not oscillate, so the “sign problem” disappears and the integral can be calculated much more effectively. The developed approach allows also finding saddle points in the complexified space of path integral integration. Some simple test calculations and comparisons with available analytical results have been carried out.

Journal of Applied Mathematics and Physics•Jan 1, 2020•1 citation
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.
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 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.
Quantum tunneling is mostly discussed in the Euclidean path integral formalism using instantons. On the other hand, it is difficult to understand quantum tunneling based on the real-time path integral due to its oscillatory nature, which causes the notorious sign problem. We show that recent development of the Lefschetz thimble method enables us to investigate this issue numerically. In particular, we find that quantum tunneling occurs due to complex trajectories, which are actually observable experimentally by using the so-called weak measurement.
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.
The major difficulty for the Feynman Path Integral Monte Carlo (PIMC) simulations of the quantum systems of particles is the so called “sign problem”, arising due to the fast oscillations of the path integral integrand depending on the complex-valued action. Our aim is to find universal techniques being able to solve this problem. The new method combines the basic ideas of the Metropolis and Hasting algorithms and is based on the Picard-Lefschetz theory and complex-valued version of Morse theory. The basic idea is to choose the Lefschetz thimbles as manifolds approaching the saddle point of the integrand. On this thimble the imaginary part of the complex-valued action remains constant. As a result the integrand on each thimble does not oscillate, so the “sign problem” disappears and the integral can be calculated much more effectively. The developed approach allows also finding saddle points in the complexified space of path integral integration. Some simple test calculations and comparisons with available analytical results have been carried out.