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:1610.08964

Analytically continued physical states in the path-integral: a sign-problem-free Quantum Monte Carlo simulation of Bell states dynamics

Evgeny Polyakov, Alexey N. Rubtsov

Published Oct 25, 2016• 0 citations
arXivHTMLPDF

Abstract

The derivation of path integrals is reconsidered. It is shown that the expression for the discretized action is not unique, and the path integration domain can be deformed so that at least Gaussian path integrals become probabillistic. This leads to a practical algorithm of sign-problem-free Monte Carlo sampling from the Gaussian path integrals. Moreover, the dynamical influence of Gaussian quantum system (the bath) on any other quantum system can be exactly represented as interaction with classical non-Markovian noise. We discuss the relation of these findings to the Bell's theorem and the Feynman's conjecture on the exponential complexity of the classical simulation of quantum systems. In Feynman's path integral we have quasiprobability distributions for trajectories, and in analitycally continued path integrals we have probability distributions for quasitrajectories.

DOI
10.48550/arxiv.1610.08964

Related papers

Phase space path integrals in Monte Carlo quantum dynamics

Stavros Caratzoulas, Philip Pechukas

In this article we present two phase space path integrals in terms of coherent states. The first one is derived in a standard fashion but using a nonstandard resolution of the identity in terms of coherent states with different width parameters. The second path integral emerges from a novel phase space representation in terms of coherent states distributed on n-dimensional manifolds embedded in the 2n-dimensional phase space of an n-degree-of-freedom system. These states are shown to form locally complete basis sets since we show that fairly smooth and localized functions can be expanded in terms of them in a unique way. In this representation the time evolution operator can be cast in the form of a phase space path integral. Both path integrals can be evaluated by straightforward implementation of Monte Carlo methods. In both cases the probability amplitude between two phase points turns out to be proportional to the average of the phase, eih∫(p dq−H dt), over a Markov process of phase space paths consisting of classical trajectories interrupted at intervals by Gaussian ‘‘quantum jumps.’’ The numerical evaluation of these expressions through importance sampling is demonstrated.

The Journal of Chemical Physics•Apr 22, 1996•10 citations

Exact Classical Effective Potential

John H. Samson

A quantum spin system can be modelled by an equivalent classical system, with an effective Hamiltonian obtained by integrating all non-zero frequency modes out of the path integral. The effective Hamiltonian H_eff(S_i) derived from the coherent-state integral is highly singular: the quasiprobability density exp(-beta H_eff), a Wigner function, imposes quantisation through derivatives of delta functions. This quasiprobability is the distribution of the time-averaged lower symbol of the spin in the coherent-state integral. We relate the quantum Monte Carlo minus-sign problem to the non-positivity of this quasiprobability, both analytically and by Monte Carlo integration.

arXiv:quant-ph/9905089•May 26, 1999•0 citations

Path-integral formulation of bosonic Markovian open quantum dynamics with Monte Carlo stochastic trajectories using the Glauber–Sudarshan P, Wigner, and Husimi Q functions and hybrids

Toma Yoneya, Kazuya Fujimoto, Yuki Kawaguchi

The Monte Carlo trajectory sampling of stochastic differential equations based on the quasiprobability distribution functions, such as the Glauber–Sudarshan P, Wigner, and Husimi Q functions, enables us to investigate bosonic open quantum many-body dynamics described by the Gorini–Kossakowski–Sudarshan–Lindblad (GKSL) equation. In this method, the Monte Carlo samplings for the initial distribution and stochastic noises incorporate quantum fluctuations, and thus, we can go beyond the mean-field approximation. However, description using stochastic differential equations is possible only when the corresponding Fokker–Planck equation has a positive-semidefinite diffusion matrix. In this work, we analytically derive the stochastic differential equations for arbitrary Hamiltonian and jump operators based on the path-integral formula, independently of the derivation of the Fokker–Planck equation. In the course of the derivation, we formulate the path-integral representation of the GKSL equation by using the -ordered quasiprobability distribution function, which systematically describes the aforementioned quasiprobability distribution functions by changing the real parameter . The essential point of this derivation is that we employ the Hubbard–Stratonovich transformation in the path integral, and its application is not always feasible. We find that the feasible condition of the Hubbard–Stratonovich transformation is identical to the positive-semidefiniteness condition of the diffusion matrix in the Fokker–Planck equation. In the benchmark calculations, we confirm that the Monte Carlo simulations of the obtained stochastic differential equations well reproduce the exact dynamics of physical quantities and non-equal time correlation functions of numerically solvable models, including the Bose–Hubbard model. This work clarifies the applicability of the approximation and gives systematic and simplified procedures to obtain the stochastic differential equations to be numerically solved.

arXiv:2508.01991•Annals of Physics•Dec 16, 2025•0 citations

Exact classical effective potentials

John H. Samson

A quantum spin system can be modelled by an equivalent classical system, with an effective Hamiltonian obtained by integrating all non-zero frequency modes out of the path integral. The effective Hamiltonian Heff({Si}) derived from the coherent-state integral is highly singular: the quasiprobability density exp(−βHeff), a Wigner function, imposes quantisation through derivatives of delta functions. This quasiprobability is the distribution of the timeaveraged lower symbol of the spin in the coherent-state integral. We relate the quantum Monte Carlo minus-sign problem to the non-positivity of this quasiprobability, both analytically and by Monte Carlo integration. 1

•May 26, 1999•0 citations

The Feynman path integral

Mark E. Tuckerman

Abstract The focus of Chapter 12 is the Feynman path integral formulations of quantum mechanics and quantum statistical mechanics. The connection between real-time quantum propagator and the density matrix is established via the Wick rotation, and a Trotter expansion is used to derive discrete path integral expressions for both. The continuous limit of the discrete path integral is derived to yield the functional integral formulation of the path integral. The construction of estimators and their averages for computing quantum observables is discussed. Next, the path integral is extended to incorporate Boltzmann, Fermi-Dirac, and Bose-Einstein statistics, and a useful recursive formula for the partition function is derived that accelerates convergence by building up N-body partition functions from those of smaller systems. Quantum free energies are discussed, and algorithms for the evaluation of path integrals using molecular dynamics or Monte Carlo are presented. The chapter concludes with a discussion of higher-order factorizations and ring-polymer contraction methods to reduce the computational overhead of path-integral calculations

•Aug 1, 2023•0 citations

Rejection-free quantum Monte Carlo in continuous time from transition path sampling

Luke Causer, Konstantinos Sfairopoulos, Jamie F. Mair +1 more

Continuous-time quantum Monte Carlo refers to a class of algorithms designed to sample the thermal distribution of a quantum Hamiltonian through exact expansions of the Boltzmann exponential in terms of stochastic trajectories which are periodic in imaginary time. Here, we show that for (sign-problem-free) quantum many-body systems with discrete degrees of freedom—such as spins on a lattice—this sampling can be done in a rejection-free manner using transition path sampling (TPS). The key idea is to converge the trajectory ensemble through updates where one individual degree of freedom is modified across all time while the remaining unaltered ones provide a time-dependent background. The ensuing single-body dynamics provides a way to generate trajectory updates exactly, allowing one to obtain the target ensemble efficiently via rejection-free TPS. We demonstrate our method on the transverse field Ising model in one and two dimensions, and on the quantum triangular plaquette (or Newman-Moore) model. We show that despite large autocorrelation times, our method is able to efficiently recover the respective quantum phase transition of each model. We also discuss the connection to rare event sampling in continuous-time Markov dynamics. Published by the American Physical Society 2024

arXiv:2305.08935•Physical review. B./Physical review. B•Jan 16, 2024•2 citations
The derivation of path integrals is reconsidered. It is shown that the expression for the discretized action is not unique, and the path integration domain can be deformed so that at least Gaussian path integrals become probabillistic. This leads to a practical algorithm of sign-problem-free Monte Carlo sampling from the Gaussian path integrals. Moreover, the dynamical influence of Gaussian quantum system (the bath) on any other quantum system can be exactly represented as interaction with classical non-Markovian noise. We discuss the relation of these findings to the Bell's theorem and the Feynman's conjecture on the exponential complexity of the classical simulation of quantum systems. In Feynman's path integral we have quasiprobability distributions for trajectories, and in analitycally continued path integrals we have probability distributions for quasitrajectories.
In this article we present two phase space path integrals in terms of coherent states. The first one is derived in a standard fashion but using a nonstandard resolution of the identity in terms of coherent states with different width parameters. The second path integral emerges from a novel phase space representation in terms of coherent states distributed on n-dimensional manifolds embedded in the 2n-dimensional phase space of an n-degree-of-freedom system. These states are shown to form locally complete basis sets since we show that fairly smooth and localized functions can be expanded in terms of them in a unique way. In this representation the time evolution operator can be cast in the form of a phase space path integral. Both path integrals can be evaluated by straightforward implementation of Monte Carlo methods. In both cases the probability amplitude between two phase points turns out to be proportional to the average of the phase, eih∫(p dq−H dt), over a Markov process of phase space paths consisting of classical trajectories interrupted at intervals by Gaussian ‘‘quantum jumps.’’ The numerical evaluation of these expressions through importance sampling is demonstrated.
A quantum spin system can be modelled by an equivalent classical system, with an effective Hamiltonian obtained by integrating all non-zero frequency modes out of the path integral. The effective Hamiltonian H_eff(S_i) derived from the coherent-state integral is highly singular: the quasiprobability density exp(-beta H_eff), a Wigner function, imposes quantisation through derivatives of delta functions. This quasiprobability is the distribution of the time-averaged lower symbol of the spin in the coherent-state integral. We relate the quantum Monte Carlo minus-sign problem to the non-positivity of this quasiprobability, both analytically and by Monte Carlo integration.
The Monte Carlo trajectory sampling of stochastic differential equations based on the quasiprobability distribution functions, such as the Glauber–Sudarshan P, Wigner, and Husimi Q functions, enables us to investigate bosonic open quantum many-body dynamics described by the Gorini–Kossakowski–Sudarshan–Lindblad (GKSL) equation. In this method, the Monte Carlo samplings for the initial distribution and stochastic noises incorporate quantum fluctuations, and thus, we can go beyond the mean-field approximation. However, description using stochastic differential equations is possible only when the corresponding Fokker–Planck equation has a positive-semidefinite diffusion matrix. In this work, we analytically derive the stochastic differential equations for arbitrary Hamiltonian and jump operators based on the path-integral formula, independently of the derivation of the Fokker–Planck equation. In the course of the derivation, we formulate the path-integral representation of the GKSL equation by using the -ordered quasiprobability distribution function, which systematically describes the aforementioned quasiprobability distribution functions by changing the real parameter . The essential point of this derivation is that we employ the Hubbard–Stratonovich transformation in the path integral, and its application is not always feasible. We find that the feasible condition of the Hubbard–Stratonovich transformation is identical to the positive-semidefiniteness condition of the diffusion matrix in the Fokker–Planck equation. In the benchmark calculations, we confirm that the Monte Carlo simulations of the obtained stochastic differential equations well reproduce the exact dynamics of physical quantities and non-equal time correlation functions of numerically solvable models, including the Bose–Hubbard model. This work clarifies the applicability of the approximation and gives systematic and simplified procedures to obtain the stochastic differential equations to be numerically solved.
A quantum spin system can be modelled by an equivalent classical system, with an effective Hamiltonian obtained by integrating all non-zero frequency modes out of the path integral. The effective Hamiltonian Heff({Si}) derived from the coherent-state integral is highly singular: the quasiprobability density exp(−βHeff), a Wigner function, imposes quantisation through derivatives of delta functions. This quasiprobability is the distribution of the timeaveraged lower symbol of the spin in the coherent-state integral. We relate the quantum Monte Carlo minus-sign problem to the non-positivity of this quasiprobability, both analytically and by Monte Carlo integration. 1
Abstract The focus of Chapter 12 is the Feynman path integral formulations of quantum mechanics and quantum statistical mechanics. The connection between real-time quantum propagator and the density matrix is established via the Wick rotation, and a Trotter expansion is used to derive discrete path integral expressions for both. The continuous limit of the discrete path integral is derived to yield the functional integral formulation of the path integral. The construction of estimators and their averages for computing quantum observables is discussed. Next, the path integral is extended to incorporate Boltzmann, Fermi-Dirac, and Bose-Einstein statistics, and a useful recursive formula for the partition function is derived that accelerates convergence by building up N-body partition functions from those of smaller systems. Quantum free energies are discussed, and algorithms for the evaluation of path integrals using molecular dynamics or Monte Carlo are presented. The chapter concludes with a discussion of higher-order factorizations and ring-polymer contraction methods to reduce the computational overhead of path-integral calculations
Continuous-time quantum Monte Carlo refers to a class of algorithms designed to sample the thermal distribution of a quantum Hamiltonian through exact expansions of the Boltzmann exponential in terms of stochastic trajectories which are periodic in imaginary time. Here, we show that for (sign-problem-free) quantum many-body systems with discrete degrees of freedom—such as spins on a lattice—this sampling can be done in a rejection-free manner using transition path sampling (TPS). The key idea is to converge the trajectory ensemble through updates where one individual degree of freedom is modified across all time while the remaining unaltered ones provide a time-dependent background. The ensuing single-body dynamics provides a way to generate trajectory updates exactly, allowing one to obtain the target ensemble efficiently via rejection-free TPS. We demonstrate our method on the transverse field Ising model in one and two dimensions, and on the quantum triangular plaquette (or Newman-Moore) model. We show that despite large autocorrelation times, our method is able to efficiently recover the respective quantum phase transition of each model. We also discuss the connection to rare event sampling in continuous-time Markov dynamics. Published by the American Physical Society 2024