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

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

Vladimir S. Filinov, Alexander S. Larkin

Published Mar 2, 2020• 0 citations
arXivHTMLPDF

Abstract

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.

DOI
10.48550/arxiv.2003.01549

Related papers

Solution of the 'sign problem' in pair approximation

Alexander S. Larkin, Vladimir S. Filinov, Vladimir E. Fortov

The main difficulty for path integral Monte Carlo studies of Fermi systems results from the requirement of antisymmetrization of the density matrix and is known in literature as the 'sign problem'. To overcome this issue the new numerical version of the %{\large {\tiny {\Huge }}} Wigner approach to quantum mechanics for treatment thermodynamic properties of degenerate systems of fermions has been developed. The new path integral representation of quantum Wigner function in the phase space has been obtained for canonical ensemble. Explicit analytical expression of the Wigner function accounting for Fermi statistical effects by effective pair pseudopotential has been proposed. Derived pseudopotential depends on coordinates, momenta and degeneracy parameter of fermions and takes into account Pauli blocking of fermions in phase space. The new quantum Monte-Carlo method for calculations of average values of arbitrary quantum operators has been proposed. To test the developed approach calculations of the momentum distribution function of the degenerate ideal system of Fermi particles has been carried out in a good agreement with analytical Fermi distributions. Generalization of this approach for studies influence of interparticle interaction on momentum distribution functions of strongly coupled Fermi system is in progress.

arXiv:1702.04091•Feb 14, 2017•0 citations

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

Lefschetz-thimble analysis of the sign problem in one-site fermion model

Yuya Tanizaki, Yoshimasa Hidaka, Tomoya Hayata

The Lefschetz-thimble approach to path integrals is applied to a one-site model of electrons, i.e., the one-site Hubbard model. Since the one-site Hubbard model shows a non-analytic behavior at the zero temperature and its path integral expression has the sign problem, this toy model is a good testing ground for an idea or a technique to attack the sign problem. Semiclassical analysis using complex saddle points unveils the significance of interference among multiple Lefschetz thimbles to reproduce the non-analytic behavior by using the path integral. If the number of Lefschetz thimbles is insufficient, we found not only large discrepancies from the exact result, but also thermodynamic instabilities. Analyzing such singular behaviors semiclassically, we propose a criterion to identify the necessary number of Lefschetz thimbles. We argue that this interference of multiple saddle points is a key issue to understand the sign problem of the finite-density quantum chromodynamics.

arXiv:1509.07146•New Journal of Physics•Mar 1, 2016•91 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

Phase Space Path Integral Representation for Wigner Function

Alexander S. Larkin, Vladimir S. Filinov

Quantum interference and exchange statistical effects can affect the momentum distribution functions making them non-Maxwellian.Such effects may be important in studies of kinetic properties of matter at low temperatures and under extreme conditions.In this work we have generalized the path integral representation for Wigner function to strongly coupled three-dimensional quantum system of particles with Boltzmann and Fermi statistics.In suggested approach the explicit expression for Wigner function was obtained in harmonic approximation and Monte Carlo method allowing numerical calculation of Wigner function, distribution functions and average quantum values has been developed.As alternative more accurate single-momentum approach and related Monte Carlo method have been developed to calculation of the distribution functions of degenerate system of interacting fermions.It allows partially overcoming the well-known sign problem for degenerate Fermi systems.

Journal of Applied Mathematics and Physics•Jan 1, 2017•16 citations

Reducing the sign problem with line integrals

Rasmus Larsen

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.

arXiv:2205.02257•Journal of High Energy Physics•Feb 7, 2025•1 citation
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 main difficulty for path integral Monte Carlo studies of Fermi systems results from the requirement of antisymmetrization of the density matrix and is known in literature as the 'sign problem'. To overcome this issue the new numerical version of the %{\large {\tiny {\Huge }}} Wigner approach to quantum mechanics for treatment thermodynamic properties of degenerate systems of fermions has been developed. The new path integral representation of quantum Wigner function in the phase space has been obtained for canonical ensemble. Explicit analytical expression of the Wigner function accounting for Fermi statistical effects by effective pair pseudopotential has been proposed. Derived pseudopotential depends on coordinates, momenta and degeneracy parameter of fermions and takes into account Pauli blocking of fermions in phase space. The new quantum Monte-Carlo method for calculations of average values of arbitrary quantum operators has been proposed. To test the developed approach calculations of the momentum distribution function of the degenerate ideal system of Fermi particles has been carried out in a good agreement with analytical Fermi distributions. Generalization of this approach for studies influence of interparticle interaction on momentum distribution functions of strongly coupled Fermi system is in progress.
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.
The Lefschetz-thimble approach to path integrals is applied to a one-site model of electrons, i.e., the one-site Hubbard model. Since the one-site Hubbard model shows a non-analytic behavior at the zero temperature and its path integral expression has the sign problem, this toy model is a good testing ground for an idea or a technique to attack the sign problem. Semiclassical analysis using complex saddle points unveils the significance of interference among multiple Lefschetz thimbles to reproduce the non-analytic behavior by using the path integral. If the number of Lefschetz thimbles is insufficient, we found not only large discrepancies from the exact result, but also thermodynamic instabilities. Analyzing such singular behaviors semiclassically, we propose a criterion to identify the necessary number of Lefschetz thimbles. We argue that this interference of multiple saddle points is a key issue to understand the sign problem of the finite-density quantum chromodynamics.
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.
Quantum interference and exchange statistical effects can affect the momentum distribution functions making them non-Maxwellian.Such effects may be important in studies of kinetic properties of matter at low temperatures and under extreme conditions.In this work we have generalized the path integral representation for Wigner function to strongly coupled three-dimensional quantum system of particles with Boltzmann and Fermi statistics.In suggested approach the explicit expression for Wigner function was obtained in harmonic approximation and Monte Carlo method allowing numerical calculation of Wigner function, distribution functions and average quantum values has been developed.As alternative more accurate single-momentum approach and related Monte Carlo method have been developed to calculation of the distribution functions of degenerate system of interacting fermions.It allows partially overcoming the well-known sign problem for degenerate Fermi systems.
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.