We apply the Generalised Thimble approach to the computation of exact path integrals and correlators in real-time quantum field theory. We first investigate the details of the numerical implementation and ways of optimizing the algorithm. We subsequently apply the method to an interacting two-field system in 0+1 dimensions, illustrating the scope for addressing realistic physical processes using real-time Generalised Thimble computations.
We apply the Generalised Thimble approach to the computation of exact path integrals and correlators in real-time quantum field theory. We first investigate the details of the numerical implementation and ways of optimizing the algorithm. We subsequently apply the method to an interacting two-field system in 0+1 dimensions, illustrating the scope for addressing realistic physical processes using real-time Generalised Thimble computations.
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.
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 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.
Relativistic quantum field theory (QFT) is commonly formulated in terms of operators, asymptotic states, and covariant amplitudes, a perspective that tends to obscure the real-time origin of field dynamics and correlations. Here we formulate quantum fields in a real-time Schrödinger-picture framework, in which fields evolve as probability ensembles on the space of field configurations. Within this formulation, the wavefunctional Ψ[φ,t] encodes a first-order, causal ensemble dynamics on configuration space. Interactions appear as couplings between configuration-space directions, while propagators arise as derived correlation structures rather than as fundamental postulates. Entanglement, scattering amplitudes, and conformal field theory correlators emerge as distinct projections of the same underlying ensemble evolution, corresponding to equal-time, asymptotic, and symmetry-organized observables. Standard operator, diagrammatic, and path-integral formulations are recovered as computational representations of this single real-time dynamics. This organization makes explicit the distinction between fundamental dynamical structure and representational tools in QFT, and clarifies the scope within which ensemble-averaged correlators account for quantum fluctuations, while also delineating the level at which questions associated with individual realizations and randomness would arise beyond the correlator-based field-theoretic description.
The quantum Monte Carlo method on asymptotic Lefschetz thimbles is a numerical algorithm devised specifically for alleviation of the sign problem appearing in the simulations of quantum many-body systems. In this method, the sign problem is alleviated by shifting the integration domain for the auxiliary fields, appearing, for example, in the conventional determinant quantum Monte Carlo method, from real space to an appropriate manifold in complex space. Here, we extend this method to quantum spin models with generic two-spin interactions, by using the Hubbard-Stratonovich transformation to decouple the exchange interactions and the Popov-Fedotov transformation to map the quantum spins to complex fermions. As a demonstration, we apply the method to the Kitaev model in a magnetic field whose ground state is predicted to deliver a topological quantum spin liquid with non-Abelian anyonic excitations. To illustrate how the sign problem is alleviated in this method, we visualize the asymptotic Lefschetz thimbles in complex space, together with the saddle points and the zeros of the fermion determinant. We benchmark our method in the low-temperature region in a magnetic field and show that the sign of the action is recovered considerably and unbiased numerical results are obtained with sufficient precision.