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.
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.
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.
Real-time path integrals provide a natural formulation of quantum dynamics also in situations where a Hamiltonian description is unavailable or impractical, but their direct evaluation is obstructed by severe phase cancellation from the oscillatory weight.A guiding theme of this thesis is that, after complexifying the integration variables, real-time path integrals are naturally organized by complex saddle points.In favorable cases the original contour can be replaced by steepest-descent manifolds (Lefschetz thimbles) attached to the relevant saddles, providing a controlled saddle-based description of real-time observables.In Part I, we study quantum decoherence in a closed-system real-time path-integral formulation: the total system (subsystem plus environment) evolves unitarily, and the environment is traced out only at the final time to obtain the reduced density matrix.For the Caldeira-Leggett model with Gaussian initial conditions, the discretized effective action becomes a quadratic form, so that the reduced density matrix is characterized by a complex saddle point and the thimble integral reduces to a finite-dimensional Gaussian integral.The complex saddle is obtained by solving a large but sparse linear system, allowing an exact evaluation of the reduced density matrix without relying on the usual master-equation approximations.We then extend this methodology to a three-dimensional cubic lattice toy model by separating center-of-mass and relative coordinates, treating the center-of-mass mode as the subsystem and the relative modes as an effective environment.This enables controlled numerical studies of how the suppression of off-diagonal matrix elements depends on the coupling, inverse temperature, and system size, and allows us to extract the early-time scaling behavior of decoherence diagnostics from the complex-saddle evaluation.In Part II, we investigate spacetime dynamics in the Lorentzian IKKT matrix model.To define a Lorentzian formulation that is not continuously connected to the Euclidean model by the standard contour deformation, we include a Lorentz invariant mass term with coefficient γ, for which the classical equations admit real, smooth, expanding solutions when γ > 0. Since the Lorentz group is noncompact, the path integral can contain a divergent volume from the boost directions; we remove this volume nonperturbatively by implementing a boost "gauge-fixing" that selects the frame minimizing tr(A 0 ) 2 , leading to the condition tr(A 0 A i ) = 0 with the associated Faddeev-Popov determinant.We then apply the complex Langevin method to the resulting complex-weight matrix integral.Time and spatial observables are defined from matrix configurations, with holomorphic time ordering implemented by auxiliary variables, and spatial structure characterized through block observables such as the extent R 2 (t) and the moment-of-inertia tensor.In the bosonic model we find approximately real time and no indication of spontaneous breaking of the spatial symmetry in the i parameter range studied.Including fermions, we control the singular-drift problem induced by the Pfaffian by introducing a fermionic mass regulator m f and an anisotropic mass deformation, which enhances fermionic effects while maintaining numerical stability.In this setup we observe a real, expanding spatial extent and a dynamical dominance of three spatial directions at late times, consistent with a (3 + 1)-dimensional expanding spacetime.We finally discuss the intended limiting procedure-taking the large-N limit first and then approaching the undeformed theory-and comment on future directions, including a comparison with partially supersymmetric deformations that preserve 16 supercharges.
Complex (semi-)classical paths, or instantons, form an integral part of our understanding of quantum physics. Whereas real classical paths describe classically allowed transitions in the real-time Feynman path integral, classically forbidden evolution is captured by complex semi-classical paths or instantons. In this paper, we uncover the rich, intricate nature of complex semi-classical paths and interference in the Feynman propagator of a non-relativistic quantum particle in both a smooth Woods-Saxon and a discontinuous Heaviside step potential. We demonstrate that the complex semi-classical paths are connected to caustics and may cease to exist as naive solutions to the boundary value problem when the semi-classical path encounters singularities of the potential. We generalise complex semi-classical paths to equivalence classes. Using this generalisation, we track the contribution of complex semi-classical paths beyond these singularity crossings and identify the instanton responsible for the quantum reflection. Whereas most complex contributions to the path integral are small, we demonstrate that in some cases the contribution of the complex semi-classical path is unsuppressed and persists into the semi-classical limit. Conversely, the Heaviside step potential includes a real classical path that does not contribute in the semi-classical limit. Finally, we develop methods to detect the presence of complex semi-classical paths from propagation amplitudes. The structure of complex semi-classical paths and methods developed here generalises to a large set of problems in real-time quantum physics.
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^ 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.
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.