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:2205.02257 · Journal of High Energy Physics

Reducing the sign problem with line integrals

Rasmus Larsen

Published Feb 7, 2025• 1 citation
arXivHTMLPDF

Abstract

A bstract We present a novel strategy to strongly reduce the severity of the sign problem, using line integrals along paths of changing imaginary action. Highly oscillating regions along these paths cancel out, decreasing their contributions. As a result, sampling with standard Monte-Carlo techniques becomes possible in cases that otherwise require methods taking advantage of complex analysis, such as Lefschetz-thimbles or Complex Langevin. We lay out how to write down an ordinary differential equation for the line integrals. As an example of its usage, we apply the results to a 1d quantum mechanical anharmonic oscillator with a x 4 potential in real time, finite temperature.

DOI
10.1007/jhep02(2025)041

Related papers

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

Sign problems in path integral formulations of quantum mechanics and quantum statistics

Vladimir S. Filinov, Alexander S. Larkin

Nowadays the term 'sign problem' is used to identify two different problems. The ideas to overcome the first type of the 'sign problem' of strongly oscillating complex valued imtegrand in the Feynman path integrals comes from Picard-Lefschetz theory and a complex version of Morse theory. The main idea is to select Lefschetz thimbles as the cycle approaching the critical point at the path-integration, where the imaginary part of the complex action stays constant. Since the imaginary part of the action is constant on each thimble, the sign problem disappears and the integral can be calculated much more effectively. Here based on the Metropolis -- Hastings algorithm a new method of calculations of the integral of the strongly oscillating integrands has been prosed. Some simple test calculation and comparison with available analytical results have been carried out. The second type of the 'sign problem' arises at studies Fermi systems by path integral approach and is caused by the requirement of antisymmetrization of the real valued matrix elements of the density matrix. An explicit analytical expression for effective pair pseudopotential in phase space has been discussed in Wigner formulation of quantum mechanics. Obtained pseudopotential allow to account for Fermi statistical effects as realizes the Pauli blocking of fermions due to the repulsion between identical fermions, which prevents their occupation of same phase space cell. To test this approach, calculations of the momentum distribution function of the ideal system of Fermi particles have been presented over a wide range of momentum and degeneracy parameter.

arXiv:2003.01549•Mar 2, 2020•0 citations

Fermionic sign problem minimization by constant path integral contour shifts

Christoph Gäntgen, Evan Berkowitz, Thomas Luu +2 more

The path integral formulation of quantum mechanical problems including fermions is often affected by a severe numerical sign problem. We show how such a sign problem can be alleviated by a judiciously chosen constant imaginary offset to the path integral. Such integration contour deformations introduce no additional computational cost to the Hamiltonian Monte Carlo algorithm, while its effective sample size is greatly increased. This makes otherwise unviable simulations efficient for a wide range of parameters. Applying our method to the Hubbard model, we find that the sign problem is significantly reduced. Furthermore, we prove that it vanishes completely for large chemical potentials, a regime where the sign problem is expected to be particularly severe without imaginary offsets. In addition to a numerical analysis of such optimized contour shifts, we analytically compute the shifts corresponding to the leading and next-to-leading order corrections to the action. We find that such simple approximations, free of significant computational cost, suffice in many cases. We present a simulation of ${\mathrm{C}}_{60}$ fullerenes (buckyballs) that are successful over a wide parameter range.

arXiv:2307.06785•Physical review. B./Physical review. B•May 20, 2024•8 citations

Monte Carlo Study of Real Time Dynamics on the Lattice

Andrei Alexandru, Gökçe Başar, Paulo F. Bedaque +2 more

Monte Carlo studies involving real time dynamics are severely restricted by the sign problem that emerges from a highly oscillatory phase of the path integral. In this Letter, we present a new method to compute real time quantities on the lattice using the Schwinger-Keldysh formalism via Monte Carlo simulations. The key idea is to deform the path integration domain to a complex manifold where the phase oscillations are mild and the sign problem is manageable. We use the previously introduced "contraction algorithm" to create a Markov chain on this alternative manifold. We substantiate our approach by analyzing the quantum mechanical anharmonic oscillator. Our results are in agreement with the exact ones obtained by diagonalization of the Hamiltonian. The method we introduce is generic and, in principle, applicable to quantum field theory albeit very slow. We discuss some possible improvements that should speed up the algorithm.

arXiv:1605.08040•Physical Review Letters•Aug 17, 2016•106 citations

A simple approach towards the sign problem using path optimisation

Francis Bursa, Michael Kroyter

A bstract We suggest an approach for simulating theories with a sign problem that relies on optimisation of complex integration contours that are not restricted to lie along Lefschetz thimbles. To that end we consider the toy model of a one-dimensional Bose gas with chemical potential. We identify the main contribution to the sign problem in this case as coming from a nearest neighbour interaction and approximately cancel it by an explicit deformation of the integration contour. We extend the obtained expressions to more general ones, depending on a small set of parameters. We find the optimal values of these parameters on a small lattice and study their range of validity. We also identify precursors for the onset of the sign problem. A fast method of evaluating the Jacobian related to the contour deformation is proposed and its numerical stability is examined. For a particular choice of lattice parameters, we find that our approach increases the lattice size at which the sign problem becomes serious from L ≈ 32 to L ≈ 700. The efficient evaluation of the Jacobian ( O ( L ) for a sweep) results in running times that are of the order of a few minutes on a standard laptop.

arXiv:1805.04941•Journal of High Energy Physics•Dec 1, 2018•33 citations

Lefschetz thimble quantum Monte Carlo for spin systems

T. C. Mooney, Jacob Bringewatt, Neill C. Warrington +1 more

Monte Carlo simulations are useful tools for modeling quantum systems, but in some cases they suffer from a sign problem, leading to an exponential slow down in their convergence to a value. While solving the sign problem is generically NP hard, many techniques exist for mitigating the sign problem in specific cases; in particular, the technique of deforming the Monte Carlo simulation's plane of integration onto Lefschetz thimbles (complex hypersurfaces of stationary phase) has seen significant success in the context of quantum field theories. We extend this methodology to spin systems by utilizing spin coherent state path integrals to reexpress the spin system's partition function in terms of continuous variables. Using some toy systems, we demonstrate its effectiveness at lessening the sign problem in this setting, despite the fact that the initial mapping to spin coherent states introduces its own sign problem. The standard formulation of the spin coherent path integral is known to make use of uncontrolled approximations; despite this, for large spins they are typically considered to yield accurate results, so it is somewhat surprising that our results show significant systematic errors. Therefore, possibly of independent interest, our use of Lefschetz thimbles to overcome the intrinsic sign problem in spin coherent state path integral Monte Carlo enables a novel numerical demonstration of a breakdown in the spin coherent path integral.

arXiv:2110.10699•Physical review. B./Physical review. B•Dec 15, 2022•5 citations
A bstract We present a novel strategy to strongly reduce the severity of the sign problem, using line integrals along paths of changing imaginary action. Highly oscillating regions along these paths cancel out, decreasing their contributions. As a result, sampling with standard Monte-Carlo techniques becomes possible in cases that otherwise require methods taking advantage of complex analysis, such as Lefschetz-thimbles or Complex Langevin. We lay out how to write down an ordinary differential equation for the line integrals. As an example of its usage, we apply the results to a 1d quantum mechanical anharmonic oscillator with a x 4 potential in real time, finite temperature.
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.
Nowadays the term 'sign problem' is used to identify two different problems. The ideas to overcome the first type of the 'sign problem' of strongly oscillating complex valued imtegrand in the Feynman path integrals comes from Picard-Lefschetz theory and a complex version of Morse theory. The main idea is to select Lefschetz thimbles as the cycle approaching the critical point at the path-integration, where the imaginary part of the complex action stays constant. Since the imaginary part of the action is constant on each thimble, the sign problem disappears and the integral can be calculated much more effectively. Here based on the Metropolis -- Hastings algorithm a new method of calculations of the integral of the strongly oscillating integrands has been prosed. Some simple test calculation and comparison with available analytical results have been carried out. The second type of the 'sign problem' arises at studies Fermi systems by path integral approach and is caused by the requirement of antisymmetrization of the real valued matrix elements of the density matrix. An explicit analytical expression for effective pair pseudopotential in phase space has been discussed in Wigner formulation of quantum mechanics. Obtained pseudopotential allow to account for Fermi statistical effects as realizes the Pauli blocking of fermions due to the repulsion between identical fermions, which prevents their occupation of same phase space cell. To test this approach, calculations of the momentum distribution function of the ideal system of Fermi particles have been presented over a wide range of momentum and degeneracy parameter.
The path integral formulation of quantum mechanical problems including fermions is often affected by a severe numerical sign problem. We show how such a sign problem can be alleviated by a judiciously chosen constant imaginary offset to the path integral. Such integration contour deformations introduce no additional computational cost to the Hamiltonian Monte Carlo algorithm, while its effective sample size is greatly increased. This makes otherwise unviable simulations efficient for a wide range of parameters. Applying our method to the Hubbard model, we find that the sign problem is significantly reduced. Furthermore, we prove that it vanishes completely for large chemical potentials, a regime where the sign problem is expected to be particularly severe without imaginary offsets. In addition to a numerical analysis of such optimized contour shifts, we analytically compute the shifts corresponding to the leading and next-to-leading order corrections to the action. We find that such simple approximations, free of significant computational cost, suffice in many cases. We present a simulation of C60{\mathrm{C}}_{60}C60​ fullerenes (buckyballs) that are successful over a wide parameter range.
Monte Carlo studies involving real time dynamics are severely restricted by the sign problem that emerges from a highly oscillatory phase of the path integral. In this Letter, we present a new method to compute real time quantities on the lattice using the Schwinger-Keldysh formalism via Monte Carlo simulations. The key idea is to deform the path integration domain to a complex manifold where the phase oscillations are mild and the sign problem is manageable. We use the previously introduced "contraction algorithm" to create a Markov chain on this alternative manifold. We substantiate our approach by analyzing the quantum mechanical anharmonic oscillator. Our results are in agreement with the exact ones obtained by diagonalization of the Hamiltonian. The method we introduce is generic and, in principle, applicable to quantum field theory albeit very slow. We discuss some possible improvements that should speed up the algorithm.
A bstract We suggest an approach for simulating theories with a sign problem that relies on optimisation of complex integration contours that are not restricted to lie along Lefschetz thimbles. To that end we consider the toy model of a one-dimensional Bose gas with chemical potential. We identify the main contribution to the sign problem in this case as coming from a nearest neighbour interaction and approximately cancel it by an explicit deformation of the integration contour. We extend the obtained expressions to more general ones, depending on a small set of parameters. We find the optimal values of these parameters on a small lattice and study their range of validity. We also identify precursors for the onset of the sign problem. A fast method of evaluating the Jacobian related to the contour deformation is proposed and its numerical stability is examined. For a particular choice of lattice parameters, we find that our approach increases the lattice size at which the sign problem becomes serious from L ≈ 32 to L ≈ 700. The efficient evaluation of the Jacobian ( O ( L ) for a sweep) results in running times that are of the order of a few minutes on a standard laptop.
Monte Carlo simulations are useful tools for modeling quantum systems, but in some cases they suffer from a sign problem, leading to an exponential slow down in their convergence to a value. While solving the sign problem is generically NP hard, many techniques exist for mitigating the sign problem in specific cases; in particular, the technique of deforming the Monte Carlo simulation's plane of integration onto Lefschetz thimbles (complex hypersurfaces of stationary phase) has seen significant success in the context of quantum field theories. We extend this methodology to spin systems by utilizing spin coherent state path integrals to reexpress the spin system's partition function in terms of continuous variables. Using some toy systems, we demonstrate its effectiveness at lessening the sign problem in this setting, despite the fact that the initial mapping to spin coherent states introduces its own sign problem. The standard formulation of the spin coherent path integral is known to make use of uncontrolled approximations; despite this, for large spins they are typically considered to yield accurate results, so it is somewhat surprising that our results show significant systematic errors. Therefore, possibly of independent interest, our use of Lefschetz thimbles to overcome the intrinsic sign problem in spin coherent state path integral Monte Carlo enables a novel numerical demonstration of a breakdown in the spin coherent path integral.