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:2501.16323 · Physical review. D/Physical review. D.

Efficient evaluation of real-time path integrals

Job Feldbrugge, Joshua Y. L. Jones

Published Apr 11, 2025• 3 citations
arXivHTMLPDF

Abstract

The Feynman path integral has revolutionized modern approaches to quantum physics. Although the path integral formalism has proven very successful and spawned several approximation schemes, the direct evaluation of real-time path integrals is still extremely expensive and numerically delicate due to its high-dimensional and oscillatory nature. We propose an efficient method for the numerical evaluation of the real-time world-line path integral for theories where the potential is dominated by a quadratic at infinity. This is done by rewriting the high-dimensional oscillatory integral in terms of a series of low-dimensional oscillatory integrals, that we efficiently evaluate with Picard-Lefschetz theory or approximate with the eikonal approximation. Subsequently, these integrals are stitched together with a series of fast Fourier transformations to recover the lattice regularized Feynman path integral. Our method directly applies to problems in quantum mechanics, the word-line quantization of quantum field theory, and quantum gravity.

DOI
10.1103/physrevd.111.083524

Related papers

On the numerical evaluation of real-time path integrals: Double exponential integration and the Maslov correction

Rosenfelder, R.

Ooura's double exponential integration formula for Fourier transforms is applied to the oscillatory integrals occuring in the path-integral description of real-time Quantum Mechanics. Due to an inherent, implicit regularization multi-dimensional Gauss-Fresnel integrals are obtained numerically with high precision but modest number of function calls. In addition, the Maslov correction for the harmonic oscillator is evaluated numerically with an increasing number of time slices in the path integral thereby clearly demonstrating that the real-time propagator acquires an additional phase $ - π/2 $ each time the particle passes through a focal point. However, in the vicinity of these singularities an overall small damping factor is required. Prospects of evaluating scattering amplitudes of finite-range potentials by direct numerical evaluation of a real-time path integral are discussed.

arXiv:2105.02880•May 6, 2021•1 citation

Mathematical aspects of the Feynman path integral

Daniel W. Boutros

Path integration methods are of crucial importance to quantum mechanics and quantum field theory. There are multiple ways the path integral can be constructed, one method uses the link between the Fokker-Planck and the Langevin equations, as is covered in this thesis. This is related to the standard derivation of the path integral in physics, in which the time is discretised. We study the mathematical problems related to the Feynman path integral, in particular the impossibility of a Lebesgue-type measure on the space of paths. It is discussed how oscillatory integrals can be used to have a well-defined `integral' on an infinite dimensional space. This formalism is subsequently applied to both the harmonic and the anharmonic oscillator. We prove an infinite dimensional oscillatory integral exists and obtain a convergent expression for both these cases, under conditions. Examples of these conditions are the initial wavefunctions being Schwartz functions as well as conditions on the endtime, angular frequency and coupling constant. Finally, a possible approach to establish an oscillatory integral for the hydrogen atom is discussed. It is proven that the result is no longer independent of the sequence of projection operators, which is a key step towards a rigorous path integral for this system.

•Jan 1, 2020•0 citations

Existence of real time quantum path integrals

Job Feldbrugge, Neil Turok

Many interesting physical theories have analytic classical actions. We show how Feynman’s path integral may be defined non-perturbatively, for such theories, without a Wick rotation to imaginary time. We start by introducing a class of smooth regulators which render interference integrals absolutely convergent and thus unambiguous. The analyticity of the regulators allows us to use Cauchy’s theorem to deform the integration domain onto a set of relevant, complex “thimbles” (or generalized steepest descent contours) each associated with a classical saddle. The regulator can then be removed to obtain an exact, non-perturbative representation. We show why the usual method of gradient flow, used to identify relevant saddles and steepest descent “thimbles” for finite-dimensional oscillatory integrals, fails in the infinite-dimensional case. For the troublesome high frequency modes, we replace it with a method we call “eigenflow” which we employ to identify the infinite-dimensional, complex “eigenthimble” over which the real time path integral is absolutely convergent. We then bound the path integral over high frequency modes by the corresponding Wiener measure for a free particle. Using the dominated convergence theorem we infer that the interacting path integral defines a good measure. While the real time path integral is more intricate than its Euclidean counterpart, it is superior in several respects. It seems particularly well-suited to theories such as quantum gravity where the classical theory is well developed but the Euclidean path integral does not exist.

arXiv:2207.12798•Annals of Physics•Apr 20, 2023•21 citations

Scattering using real-time path integrals

Wayne Nicholas Polyzou, Ekaterina S. Nathanson

Background: Path integrals are a powerful tool for solving problems in quantum theory that are not amenable to a treatment by perturbation theory. Most path integral computations require an analytic continuation to imaginary time. While imaginary time treatments of scattering are possible, imaginary time is not a natural framework for treating scattering problems. More importantly, quantum algorithms for calculating path integrals require real-time evolution.Purpose: We test a recently introduced method for performing direct calculations of scattering observables using real-time path integrals in order to understand the challenges facing real-time path integral calculations of scattering observables.Method: The computations are based on a new interpretation of the path integral as the expectation value of a potential functional on cylinder sets of continuous paths with respect to a complex probability distribution. The method can in principle be applied to arbitrary short-range potentials.Results: The method is applied to compute matrix elements of M\o{}ller wave operators applied to narrow wave packets. These are used to calculate half-shell sharp-momentum transition matrix elements for one-dimensional potential scattering. The calculations for half-shell transition operator matrix elements converge to the numerical solution of the Lippmann-Schwinger equation.Conclusions: This work presents a proof in principle that scattering observables can be computed using real-time Feynman path integrals. While the computational method is not efficient, it can be improved. It provides a laboratory for studying quantum computational algorithms that are applicable to scattering problems.

arXiv:1712.00046•Physical Review C•Jun 4, 2020•2 citations

Feynman-Kac path integral expansion around the upside-down oscillator

Zbigniew Haba

We discuss path integrals for quantum mechanics with a potential which is a perturbation of the upside-down oscillator. We express the path integral (in the real time) by the Wiener measure. We obtain the Feynman integral for perturbations which are the Fourier-Laplace transforms of a complex measure and for polynomials of the fotm $x^{4n}$ and $x^{4n+2}$ (where $n$ is a natural number). We extend the method to quantum field theory (QFT) with complex scaled spatial coordinates ${\bf x}\rightarrow i{\bf x}$. We show that such a complex extension of the path integral (in the real time) allows a rigorous path integral treatment of a large class of potentials including the ones unbounded from below.

arXiv:2305.11989•May 19, 2023•0 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
The Feynman path integral has revolutionized modern approaches to quantum physics. Although the path integral formalism has proven very successful and spawned several approximation schemes, the direct evaluation of real-time path integrals is still extremely expensive and numerically delicate due to its high-dimensional and oscillatory nature. We propose an efficient method for the numerical evaluation of the real-time world-line path integral for theories where the potential is dominated by a quadratic at infinity. This is done by rewriting the high-dimensional oscillatory integral in terms of a series of low-dimensional oscillatory integrals, that we efficiently evaluate with Picard-Lefschetz theory or approximate with the eikonal approximation. Subsequently, these integrals are stitched together with a series of fast Fourier transformations to recover the lattice regularized Feynman path integral. Our method directly applies to problems in quantum mechanics, the word-line quantization of quantum field theory, and quantum gravity.
Ooura's double exponential integration formula for Fourier transforms is applied to the oscillatory integrals occuring in the path-integral description of real-time Quantum Mechanics. Due to an inherent, implicit regularization multi-dimensional Gauss-Fresnel integrals are obtained numerically with high precision but modest number of function calls. In addition, the Maslov correction for the harmonic oscillator is evaluated numerically with an increasing number of time slices in the path integral thereby clearly demonstrating that the real-time propagator acquires an additional phase $ - π/2 $ each time the particle passes through a focal point. However, in the vicinity of these singularities an overall small damping factor is required. Prospects of evaluating scattering amplitudes of finite-range potentials by direct numerical evaluation of a real-time path integral are discussed.
Path integration methods are of crucial importance to quantum mechanics and quantum field theory. There are multiple ways the path integral can be constructed, one method uses the link between the Fokker-Planck and the Langevin equations, as is covered in this thesis. This is related to the standard derivation of the path integral in physics, in which the time is discretised. We study the mathematical problems related to the Feynman path integral, in particular the impossibility of a Lebesgue-type measure on the space of paths. It is discussed how oscillatory integrals can be used to have a well-defined `integral' on an infinite dimensional space. This formalism is subsequently applied to both the harmonic and the anharmonic oscillator. We prove an infinite dimensional oscillatory integral exists and obtain a convergent expression for both these cases, under conditions. Examples of these conditions are the initial wavefunctions being Schwartz functions as well as conditions on the endtime, angular frequency and coupling constant. Finally, a possible approach to establish an oscillatory integral for the hydrogen atom is discussed. It is proven that the result is no longer independent of the sequence of projection operators, which is a key step towards a rigorous path integral for this system.
Many interesting physical theories have analytic classical actions. We show how Feynman’s path integral may be defined non-perturbatively, for such theories, without a Wick rotation to imaginary time. We start by introducing a class of smooth regulators which render interference integrals absolutely convergent and thus unambiguous. The analyticity of the regulators allows us to use Cauchy’s theorem to deform the integration domain onto a set of relevant, complex “thimbles” (or generalized steepest descent contours) each associated with a classical saddle. The regulator can then be removed to obtain an exact, non-perturbative representation. We show why the usual method of gradient flow, used to identify relevant saddles and steepest descent “thimbles” for finite-dimensional oscillatory integrals, fails in the infinite-dimensional case. For the troublesome high frequency modes, we replace it with a method we call “eigenflow” which we employ to identify the infinite-dimensional, complex “eigenthimble” over which the real time path integral is absolutely convergent. We then bound the path integral over high frequency modes by the corresponding Wiener measure for a free particle. Using the dominated convergence theorem we infer that the interacting path integral defines a good measure. While the real time path integral is more intricate than its Euclidean counterpart, it is superior in several respects. It seems particularly well-suited to theories such as quantum gravity where the classical theory is well developed but the Euclidean path integral does not exist.
Background: Path integrals are a powerful tool for solving problems in quantum theory that are not amenable to a treatment by perturbation theory. Most path integral computations require an analytic continuation to imaginary time. While imaginary time treatments of scattering are possible, imaginary time is not a natural framework for treating scattering problems. More importantly, quantum algorithms for calculating path integrals require real-time evolution.Purpose: We test a recently introduced method for performing direct calculations of scattering observables using real-time path integrals in order to understand the challenges facing real-time path integral calculations of scattering observables.Method: The computations are based on a new interpretation of the path integral as the expectation value of a potential functional on cylinder sets of continuous paths with respect to a complex probability distribution. The method can in principle be applied to arbitrary short-range potentials.Results: The method is applied to compute matrix elements of M\o{}ller wave operators applied to narrow wave packets. These are used to calculate half-shell sharp-momentum transition matrix elements for one-dimensional potential scattering. The calculations for half-shell transition operator matrix elements converge to the numerical solution of the Lippmann-Schwinger equation.Conclusions: This work presents a proof in principle that scattering observables can be computed using real-time Feynman path integrals. While the computational method is not efficient, it can be improved. It provides a laboratory for studying quantum computational algorithms that are applicable to scattering problems.
We discuss path integrals for quantum mechanics with a potential which is a perturbation of the upside-down oscillator. We express the path integral (in the real time) by the Wiener measure. We obtain the Feynman integral for perturbations which are the Fourier-Laplace transforms of a complex measure and for polynomials of the fotm x4nx^{4n}x4n and x4n+2x^{4n+2}x4n+2 (where nnn is a natural number). We extend the method to quantum field theory (QFT) with complex scaled spatial coordinates x→ix{\bf x}\rightarrow i{\bf x}x→ix. We show that such a complex extension of the path integral (in the real time) allows a rigorous path integral treatment of a large class of potentials including the ones unbounded from below.
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.