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

The Journal of Chemical Physics

Understanding the sign problem from an exact path integral Monte Carlo model of interacting harmonic fermions

Siu A. Chin

Published Jun 11, 2026• 0 citations
Publisher

Abstract

This work shows that the recently discovered operator contraction identity for solving the discrete path integral of the harmonic oscillator can be applied equally to fermions in any dimension. This then yields an exactly solvable model for studying the sign problem where the path integral Monte Carlo energy at any time step for any number of fermions is known analytically or can be computed numerically. It is found that the sign problem is primarily a property of the free fermion propagator, but repulsive/attractive pairwise interactions can shift the sign problem to larger/smaller imaginary time but do not make it more severe than the non-interacting case. More surprisingly, one can prove analytically that the first closed-shell state in D dimension, with n = D + 1 fermions, has no sign problem at large imaginary time. Direct numerical simulations confirm that this is also true for higher closed-shell states in two and three dimensions. Fourth-order and newly found variable-bead algorithms are used to compute ground state energies of quantum dots with up to 110 electrons and are compared to results obtained by modern neural networks.

DOI
10.1063/5.0324529

Related papers

Understanding the sign problem from an exact Path Integral Monte Carlo model of interacting harmonic fermions

Sa Chin

This work shows that the recently discovered operator contraction identity for solving the discreet Path Integral of the harmonic oscillator can be applied equally to fermions in any dimension. This then yields an exactly solvable model for studying the sign problem where the Path Integral Monte Carlo energy at any time step for any number of fermions is known analytically, or can be computed numerically. It is found that the sign problem is primarily a property of the free fermion propagator, but repulsive/attractive pairwise interaction can shift the sign problem to larger/smaller imaginary time but does not make it more severe than the non-interacting case. More surprisingly, one can prove analytically that the first closed-shell state in $D$ dimension, with $n=D+1$ fermion, has no sign problem at large imaginary time. Direct numerical simulations confirm that this is also true for higher closed-shell states in two and three-dimension. Fourth-order and newly found variable-bead algorithms are used to compute ground state energies of quantum dots with up to 110 electrons and compared to results obtained by modern neural networks.

arXiv:2601.22559•Jan 30, 2026•0 citations

The Path Integral Monte Carlo Sign Problem Is Not Always NP-Hard: Harmonic Fermions Can Be Solved in Quadratic Time

Aarif Chaudhary, Jonas Valenzuela, Sa Chin

In the Path Integral Monte Carlo (PIMC) simulation of fermions in a harmonic trap, with and without pairwise harmonic interactions, the partition functions for any discrete number of imaginary time slices (or beads) and for any choice of the short-time propagator can be analytically obtained from the contracted determinant form of the propagator. This work shows that the resulting recursion relation can be reformulated in the $λ$-ring language, yielding a closed-form finite-bead partition function in two dimensions in terms of permutation statistics. This closed-form partition function can be evaluated by a special algorithm in $O(n^2)$ time, providing an exact and numerically stable scheme for reproducing the energies of the original (undoable) fermion PIMC simulation for $n=10^4$ or more fermions. This result provides a concrete framework in which the numerical instability of the sign problem is completely bypassed, and serves as a counterexample to the prevailing view that all truly fermionic PIMC sign problems are NP-hard.

arXiv:2609.09071•Sep 8, 2026•0 citations

High-order path-integral Monte Carlo methods for solving quantum dot problems

Siu A. Chin

The conventional second-order path-integral Monte Carlo method is plagued with the sign problem in solving many-fermion systems. This is due to the large number of antisymmetric free-fermion propagators that are needed to extract the ground state wave function at large imaginary time. In this work we show that optimized fourth-order path-integral Monte Carlo methods, which use no more than five free-fermion propagators, can yield accurate quantum dot energies for up to 20 polarized electrons with the use of the Hamiltonian energy estimator.

arXiv:1411.6682•Physical Review E•Mar 11, 2015•42 citations

Simulation of excited states and the sign problem in the path integral Monte Carlo method

Alexander P. Lyubartsev

An approach is presented to compute properties of excited states in path integral Monte Carlo simulations of quantum systems. The approach is based on the introduction of several images of the studied system which have the total wavefunction antisymmetric over permutations of these images, and a simulation of the whole system at low enough temperature. The success of the approach relies, however, on the solution of the sign problem for the corresponding system. In the cases when the sign problem may be resolved (as for fermions in one dimension), properties of excited states can be computed with high precision. This was demonstrated by a simulation of excited states of a single particle in one- and two-dimensional harmonic oscillators. An example of the simulation of excited states of several interacting electrons in one dimension is also given.

Journal of Physics A Mathematical and General•Jul 13, 2005•40 citations

Simple proof that there is no sign problem in path integral Monte Carlo simulations of fermions in one dimension

Siu A. Chin

It is widely known that there is no sign problem in path integral Monte Carlo (PIMC) simulations of fermions in one dimension. As far as the author is aware, there is no direct proof of this in the literature. This work shows that the sign of the N-fermion antisymmetric free propagator is given by the product of all possible pairs of particle separations, or relative displacements. For a nonvanishing closed-loop product of such propagators, as required by PIMC, all relative displacements from adjacent propagators are paired into perfect squares, and therefore the loop product must be positive, but only in one dimension. By comparison, permutation sampling, which does not evaluate the determinant of the antisymmetric propagator exactly, remains plagued by a low-level sign problem, even in one dimension.

arXiv:2407.01618•Physical review. E•Jun 25, 2024•3 citations

Solving fermion problems without solving the sign problem: Symmetry-breaking wave functions from similarity-transformed propagators for solving two-dimensional quantum dots

Siu A. Chin

It is well known that the use of the primitive second-order propagator in path-integral Monte Carlo calculations of many-fermion systems leads to the sign problem. This work will show that by using the similarity-transformed Fokker-Planck propagator, it is possible to solve for the ground state of a large quantum dot, with up to 100 polarized electrons, without solving the sign problem. These similarity-transformed propagators naturally produce rotational symmetry-breaking ground-state wave functions previously used in the study of quantum dots and quantum Hall effects. However, instead of localizing the electrons at positions that minimize the potential energy, this derivation shows that they should be located at positions that maximize the bosonic ground-state wave function. Further improvements in the energy can be obtained by using these as initial wave functions in a ground-state path-integral Monte Carlo calculation with second- and fourth-order propagators.

arXiv:1908.10471•Physical review. E•Apr 13, 2020•3 citations
This work shows that the recently discovered operator contraction identity for solving the discrete path integral of the harmonic oscillator can be applied equally to fermions in any dimension. This then yields an exactly solvable model for studying the sign problem where the path integral Monte Carlo energy at any time step for any number of fermions is known analytically or can be computed numerically. It is found that the sign problem is primarily a property of the free fermion propagator, but repulsive/attractive pairwise interactions can shift the sign problem to larger/smaller imaginary time but do not make it more severe than the non-interacting case. More surprisingly, one can prove analytically that the first closed-shell state in D dimension, with n = D + 1 fermions, has no sign problem at large imaginary time. Direct numerical simulations confirm that this is also true for higher closed-shell states in two and three dimensions. Fourth-order and newly found variable-bead algorithms are used to compute ground state energies of quantum dots with up to 110 electrons and are compared to results obtained by modern neural networks.
This work shows that the recently discovered operator contraction identity for solving the discreet Path Integral of the harmonic oscillator can be applied equally to fermions in any dimension. This then yields an exactly solvable model for studying the sign problem where the Path Integral Monte Carlo energy at any time step for any number of fermions is known analytically, or can be computed numerically. It is found that the sign problem is primarily a property of the free fermion propagator, but repulsive/attractive pairwise interaction can shift the sign problem to larger/smaller imaginary time but does not make it more severe than the non-interacting case. More surprisingly, one can prove analytically that the first closed-shell state in DDD dimension, with n=D+1n=D+1n=D+1 fermion, has no sign problem at large imaginary time. Direct numerical simulations confirm that this is also true for higher closed-shell states in two and three-dimension. Fourth-order and newly found variable-bead algorithms are used to compute ground state energies of quantum dots with up to 110 electrons and compared to results obtained by modern neural networks.
In the Path Integral Monte Carlo (PIMC) simulation of fermions in a harmonic trap, with and without pairwise harmonic interactions, the partition functions for any discrete number of imaginary time slices (or beads) and for any choice of the short-time propagator can be analytically obtained from the contracted determinant form of the propagator. This work shows that the resulting recursion relation can be reformulated in the λλλ-ring language, yielding a closed-form finite-bead partition function in two dimensions in terms of permutation statistics. This closed-form partition function can be evaluated by a special algorithm in O(n2)O(n^2)O(n2) time, providing an exact and numerically stable scheme for reproducing the energies of the original (undoable) fermion PIMC simulation for n=104n=10^4n=104 or more fermions. This result provides a concrete framework in which the numerical instability of the sign problem is completely bypassed, and serves as a counterexample to the prevailing view that all truly fermionic PIMC sign problems are NP-hard.
The conventional second-order path-integral Monte Carlo method is plagued with the sign problem in solving many-fermion systems. This is due to the large number of antisymmetric free-fermion propagators that are needed to extract the ground state wave function at large imaginary time. In this work we show that optimized fourth-order path-integral Monte Carlo methods, which use no more than five free-fermion propagators, can yield accurate quantum dot energies for up to 20 polarized electrons with the use of the Hamiltonian energy estimator.
An approach is presented to compute properties of excited states in path integral Monte Carlo simulations of quantum systems. The approach is based on the introduction of several images of the studied system which have the total wavefunction antisymmetric over permutations of these images, and a simulation of the whole system at low enough temperature. The success of the approach relies, however, on the solution of the sign problem for the corresponding system. In the cases when the sign problem may be resolved (as for fermions in one dimension), properties of excited states can be computed with high precision. This was demonstrated by a simulation of excited states of a single particle in one- and two-dimensional harmonic oscillators. An example of the simulation of excited states of several interacting electrons in one dimension is also given.
It is widely known that there is no sign problem in path integral Monte Carlo (PIMC) simulations of fermions in one dimension. As far as the author is aware, there is no direct proof of this in the literature. This work shows that the sign of the N-fermion antisymmetric free propagator is given by the product of all possible pairs of particle separations, or relative displacements. For a nonvanishing closed-loop product of such propagators, as required by PIMC, all relative displacements from adjacent propagators are paired into perfect squares, and therefore the loop product must be positive, but only in one dimension. By comparison, permutation sampling, which does not evaluate the determinant of the antisymmetric propagator exactly, remains plagued by a low-level sign problem, even in one dimension.
It is well known that the use of the primitive second-order propagator in path-integral Monte Carlo calculations of many-fermion systems leads to the sign problem. This work will show that by using the similarity-transformed Fokker-Planck propagator, it is possible to solve for the ground state of a large quantum dot, with up to 100 polarized electrons, without solving the sign problem. These similarity-transformed propagators naturally produce rotational symmetry-breaking ground-state wave functions previously used in the study of quantum dots and quantum Hall effects. However, instead of localizing the electrons at positions that minimize the potential energy, this derivation shows that they should be located at positions that maximize the bosonic ground-state wave function. Further improvements in the energy can be obtained by using these as initial wave functions in a ground-state path-integral Monte Carlo calculation with second- and fourth-order propagators.