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:1511.02878 · Physical review. B./Physical review. B

Entanglement and the fermion sign problem in auxiliary field quantum Monte Carlo simulations

Peter Broecker, Simon Trebst

Published Aug 22, 2016• 17 citations
arXivHTMLPDF

Abstract

Quantum Monte Carlo simulations of fermions are hampered by the notorious sign problem whose most striking manifestation is an exponential growth of sampling errors with the number of particles. With the sign problem known to be an NP-hard problem and any generic solution thus highly elusive, the Monte Carlo sampling of interacting many-fermion systems is commonly thought to be restricted to a small class of model systems for which a sign-free basis has been identified. Here we demonstrate that entanglement measures, in particular the so-called R\'enyi entropies, can intrinsically exhibit a certain robustness against the sign problem in auxiliary-field quantum Monte Carlo approaches and possibly allow for the identification of global ground-state properties via their scaling behavior even in the presence of a strong sign problem. We corroborate these findings via numerical simulations of fermionic quantum phase transitions of spinless fermions on the honeycomb lattice at and below half filling.

DOI
10.1103/physrevb.94.075144

Related papers

Numerical stabilization of entanglement computation in auxiliary-field quantum Monte Carlo simulations of interacting many-fermion systems

Peter Broecker, Simon Trebst

In the absence of a fermion sign problem, auxiliary-field (or determinantal) quantum Monte Carlo (DQMC) approaches have long been the numerical method of choice for unbiased, large-scale simulations of interacting many-fermion systems. More recently, the conceptual scope of this approach has been expanded by introducing ingenious schemes to compute entanglement entropies within its framework. On a practical level, these approaches, however, suffer from a variety of numerical instabilities that have largely impeded their applicability. Here we report on a number of algorithmic advances to overcome many of these numerical instabilities and significantly improve the calculation of entanglement measures in the zero-temperature projective DQMC approach, ultimately allowing us to reach similar system sizes as for the computation of conventional observables. We demonstrate the applicability of this improved DQMC approach by providing an entanglement perspective on the quantum phase transition from a magnetically ordered Mott insulator to a band insulator in the bilayer square lattice Hubbard model at half filling.

arXiv:1609.07309•Physical review. E•Dec 19, 2016•19 citations

Hybrid Monte Carlo approach to the entanglement entropy of interacting fermions

Joaquín E. Drut, William J. Porter

The Monte Carlo calculation of R\'enyi entanglement entropies ${S}_{n}^{}$ of interacting fermions suffers from a well-known signal-to-noise problem, even for a large number of situations in which the infamous sign problem is absent. A few methods have been proposed to overcome this issue, such as ensemble switching and the use of auxiliary partition-function ratios. Here, we present an approach that builds on the recently proposed free-fermion decomposition method; it incorporates entanglement in the probability measure in a natural way; it takes advantage of the hybrid Monte Carlo algorithm (an essential tool in lattice quantum chromodynamics and other gauge theories with dynamical fermions); and it does not suffer from noise problems. This method displays no sign problem for the same cases as other approaches and is therefore useful for a wide variety of systems. As a proof of principle, we calculate ${S}_{2}^{}$ for the one-dimensional, half-filled Hubbard model and compare with results from exact diagonalization and the free-fermion decomposition method.

arXiv:1506.06654•Physical Review B•Sep 14, 2015•31 citations

Majorana Positivity and the Fermion Sign Problem of Quantum Monte Carlo Simulations

Z. C. Wei, Congjun Wu, Yi Li +2 more

The sign problem is a major obstacle in quantum Monte Carlo simulations for many-body fermion systems. We examine this problem with a new perspective based on the Majorana reflection positivity and Majorana Kramers positivity. Two sufficient conditions are proven for the absence of the fermion sign problem. Our proof provides a unified description for all the interacting lattice fermion models previously known to be free of the sign problem based on the auxiliary field quantum Monte Carlo method. It also allows us to identify a number of new sign-problem-free interacting fermion models including, but not limited to, lattice fermion models with repulsive interactions but without particle-hole symmetry, and interacting topological insulators with spin-flip terms.

arXiv:1601.01994•Physical Review Letters•Jun 23, 2016•119 citations

Computational Complexity and Fundamental Limitations to Fermionic Quantum Monte Carlo Simulations

Matthias Troyer, U.‐J. Wiese

Quantum Monte Carlo simulations, while being efficient for bosons, suffer from the "negative sign problem" when applied to fermions--causing an exponential increase of the computing time with the number of particles. A polynomial time solution to the sign problem is highly desired since it would provide an unbiased and numerically exact method to simulate correlated quantum systems. Here we show that such a solution is almost certainly unattainable by proving that the sign problem is nondeterministic polynomial (NP) hard, implying that a generic solution of the sign problem would also solve all problems in the complexity class NP in polynomial time.

arXiv:cond-mat/0408370•Physical Review Letters•May 4, 2005•1,250 citations

Sign-Problem-Free Fermionic Quantum Monte Carlo: Developments and Applications

Zi-Xiang Li, Hong Yao

Reliable simulations of correlated quantum systems, including high-temperature superconductors and frustrated magnets, are increasingly desired nowadays to further our understanding of essential features in such systems. Quantum Monte Carlo (QMC) is a unique numerically exact and intrinsically unbiased method to simulate interacting quantum many-body systems. More importantly, when QMC simulations are free from the notorious fermion sign problem, they can reliably simulate interacting quantum models with large system size and low temperature to reveal low-energy physics such as spontaneously broken symmetries and universal quantum critical behaviors. Here, we concisely review recent progress made in developing new sign-problem-free QMC algorithms, including those employing Majorana representation and those utilizing hot-spot physics. We also discuss applications of these novel sign-problem-free QMC algorithms in simulations of various interesting quantum many-body models. Finally, we discuss possible future directions of designing sign-problem-free QMC methods.

arXiv:1805.08219•Annual Review of Condensed Matter Physics•Dec 19, 2018•126 citations

Stabilizer-Rényi Microscopy of Critical Correlations in Interacting Fermions

Jun Qi Fang, Fo-Hong Wang, Xiao Yan Xu

Quantum magic---the resource that separates universal quantum computation from efficiently simulable Clifford circuits---has emerged as a diagnostic of many-body quantum states, yet the stabilizer Rényi entropy (SRE) that quantifies it remains largely inaccessible in interacting fermion systems. Estimating the global SRE requires a sum over exponentially many Majorana strings, which in determinant quantum Monte Carlo can be sampled without a sign problem only for the interacting density matrix obtained after averaging the auxiliary fields: sampling strings and fields simultaneously incurs a sign problem even for sign-problem-free models, so that the known sign-free alternative is a nested Monte Carlo that does not scale. We propose instead the two-point SRE---a practical stabilizer-Rényi correlator built from rank-2 SREs of one- and two-site reduced density matrices, which are fixed exactly by single-particle Green functions and density correlations already measured in standard simulations---and show that it obeys universal finite-size scaling at fermionic quantum critical points. In the one-dimensional half-filled spinless $t$-$V$ chain, the correlator distinguishes algebraic and exponential decay regimes and tracks the inverse-logarithmic finite-size drift characteristic of the Berezinskii--Kosterlitz--Thouless transition. On the honeycomb lattice, sign-problem-free quantum Monte Carlo yields Gross--Neveu--Ising scaling with finite anomalous dimension at zero temperature and two-dimensional Ising collapse at the thermal transition. Our results suggest stabilizer-Rényi microscopy as a spatially resolved, quantitatively universal probe of critical correlations in interacting fermionic matter, on par with conventional order-parameter correlators and accessible to quantum-simulator measurements.

cond-mat.str-elCross-listedarXiv:2601.13314•Jan 21, 2026•0 citations
Quantum Monte Carlo simulations of fermions are hampered by the notorious sign problem whose most striking manifestation is an exponential growth of sampling errors with the number of particles. With the sign problem known to be an NP-hard problem and any generic solution thus highly elusive, the Monte Carlo sampling of interacting many-fermion systems is commonly thought to be restricted to a small class of model systems for which a sign-free basis has been identified. Here we demonstrate that entanglement measures, in particular the so-called R'enyi entropies, can intrinsically exhibit a certain robustness against the sign problem in auxiliary-field quantum Monte Carlo approaches and possibly allow for the identification of global ground-state properties via their scaling behavior even in the presence of a strong sign problem. We corroborate these findings via numerical simulations of fermionic quantum phase transitions of spinless fermions on the honeycomb lattice at and below half filling.
In the absence of a fermion sign problem, auxiliary-field (or determinantal) quantum Monte Carlo (DQMC) approaches have long been the numerical method of choice for unbiased, large-scale simulations of interacting many-fermion systems. More recently, the conceptual scope of this approach has been expanded by introducing ingenious schemes to compute entanglement entropies within its framework. On a practical level, these approaches, however, suffer from a variety of numerical instabilities that have largely impeded their applicability. Here we report on a number of algorithmic advances to overcome many of these numerical instabilities and significantly improve the calculation of entanglement measures in the zero-temperature projective DQMC approach, ultimately allowing us to reach similar system sizes as for the computation of conventional observables. We demonstrate the applicability of this improved DQMC approach by providing an entanglement perspective on the quantum phase transition from a magnetically ordered Mott insulator to a band insulator in the bilayer square lattice Hubbard model at half filling.
The Monte Carlo calculation of R\'enyi entanglement entropies Sn{S}_{n}^{}Sn​ of interacting fermions suffers from a well-known signal-to-noise problem, even for a large number of situations in which the infamous sign problem is absent. A few methods have been proposed to overcome this issue, such as ensemble switching and the use of auxiliary partition-function ratios. Here, we present an approach that builds on the recently proposed free-fermion decomposition method; it incorporates entanglement in the probability measure in a natural way; it takes advantage of the hybrid Monte Carlo algorithm (an essential tool in lattice quantum chromodynamics and other gauge theories with dynamical fermions); and it does not suffer from noise problems. This method displays no sign problem for the same cases as other approaches and is therefore useful for a wide variety of systems. As a proof of principle, we calculate S2{S}_{2}^{}S2​ for the one-dimensional, half-filled Hubbard model and compare with results from exact diagonalization and the free-fermion decomposition method.
The sign problem is a major obstacle in quantum Monte Carlo simulations for many-body fermion systems. We examine this problem with a new perspective based on the Majorana reflection positivity and Majorana Kramers positivity. Two sufficient conditions are proven for the absence of the fermion sign problem. Our proof provides a unified description for all the interacting lattice fermion models previously known to be free of the sign problem based on the auxiliary field quantum Monte Carlo method. It also allows us to identify a number of new sign-problem-free interacting fermion models including, but not limited to, lattice fermion models with repulsive interactions but without particle-hole symmetry, and interacting topological insulators with spin-flip terms.
Quantum Monte Carlo simulations, while being efficient for bosons, suffer from the "negative sign problem" when applied to fermions--causing an exponential increase of the computing time with the number of particles. A polynomial time solution to the sign problem is highly desired since it would provide an unbiased and numerically exact method to simulate correlated quantum systems. Here we show that such a solution is almost certainly unattainable by proving that the sign problem is nondeterministic polynomial (NP) hard, implying that a generic solution of the sign problem would also solve all problems in the complexity class NP in polynomial time.
Reliable simulations of correlated quantum systems, including high-temperature superconductors and frustrated magnets, are increasingly desired nowadays to further our understanding of essential features in such systems. Quantum Monte Carlo (QMC) is a unique numerically exact and intrinsically unbiased method to simulate interacting quantum many-body systems. More importantly, when QMC simulations are free from the notorious fermion sign problem, they can reliably simulate interacting quantum models with large system size and low temperature to reveal low-energy physics such as spontaneously broken symmetries and universal quantum critical behaviors. Here, we concisely review recent progress made in developing new sign-problem-free QMC algorithms, including those employing Majorana representation and those utilizing hot-spot physics. We also discuss applications of these novel sign-problem-free QMC algorithms in simulations of various interesting quantum many-body models. Finally, we discuss possible future directions of designing sign-problem-free QMC methods.
Quantum magic---the resource that separates universal quantum computation from efficiently simulable Clifford circuits---has emerged as a diagnostic of many-body quantum states, yet the stabilizer Rényi entropy (SRE) that quantifies it remains largely inaccessible in interacting fermion systems. Estimating the global SRE requires a sum over exponentially many Majorana strings, which in determinant quantum Monte Carlo can be sampled without a sign problem only for the interacting density matrix obtained after averaging the auxiliary fields: sampling strings and fields simultaneously incurs a sign problem even for sign-problem-free models, so that the known sign-free alternative is a nested Monte Carlo that does not scale. We propose instead the two-point SRE---a practical stabilizer-Rényi correlator built from rank-2 SREs of one- and two-site reduced density matrices, which are fixed exactly by single-particle Green functions and density correlations already measured in standard simulations---and show that it obeys universal finite-size scaling at fermionic quantum critical points. In the one-dimensional half-filled spinless ttt-VVV chain, the correlator distinguishes algebraic and exponential decay regimes and tracks the inverse-logarithmic finite-size drift characteristic of the Berezinskii--Kosterlitz--Thouless transition. On the honeycomb lattice, sign-problem-free quantum Monte Carlo yields Gross--Neveu--Ising scaling with finite anomalous dimension at zero temperature and two-dimensional Ising collapse at the thermal transition. Our results suggest stabilizer-Rényi microscopy as a spatially resolved, quantitatively universal probe of critical correlations in interacting fermionic matter, on par with conventional order-parameter correlators and accessible to quantum-simulator measurements.