High Resolution Image Download MS PowerPoint Slide We combine the recent η-ensemble path integral Monte Carlo approach to the free energy [Dornheim et al. Phys. Rev. B 2025 111, L041114] with a recent fictitious partition function technique based on inserting a continuous variable that interpolates between the bosonic and Fermionic limits [Xiong and Xiong J. Chem. Phys. 2022 157, 094112] to deal with the Fermion sign problem. As a practical example, we apply our setup to the warm, dense, uniform electron gas over a broad range of densities and temperatures. We obtain accurate results for the exchange–correlation free energy down to half the Fermi temperature and find excellent agreement with the state-of-the-art parametrization by Groth et al. [ Phys. Rev. Lett. 2017 119, 135001]. Our work opens up new avenues for the future study of a host of interacting Fermi systems, including warm dense matter, ultracold atoms, and electrons in quantum dots, and for Fermionic free energy calculations with unprecedented system size.
High Resolution Image Download MS PowerPoint Slide We combine the recent η-ensemble path integral Monte Carlo approach to the free energy [Dornheim et al. Phys. Rev. B 2025 111, L041114] with a recent fictitious partition function technique based on inserting a continuous variable that interpolates between the bosonic and Fermionic limits [Xiong and Xiong J. Chem. Phys. 2022 157, 094112] to deal with the Fermion sign problem. As a practical example, we apply our setup to the warm, dense, uniform electron gas over a broad range of densities and temperatures. We obtain accurate results for the exchange–correlation free energy down to half the Fermi temperature and find excellent agreement with the state-of-the-art parametrization by Groth et al. [ Phys. Rev. Lett. 2017 119, 135001]. Our work opens up new avenues for the future study of a host of interacting Fermi systems, including warm dense matter, ultracold atoms, and electrons in quantum dots, and for Fermionic free energy calculations with unprecedented system size.
The fermion sign problem constitutes one of the most fundamental obstacles in quantum many-body theory. Recently, it has been suggested to circumvent the sign problem by carrying out path integral simulations with a fictitious quantum statistics variable ξ, which allows for a smooth interpolation between the bosonic and fermionic limits [Xiong and Xiong, J. Chem. Phys.157, 094112 (2022)]. This ξ-extrapolation method has subsequently been applied to a variety of systems and has facilitated the analysis of an x-ray scattering measurement taken at the National Ignition Facility with unprecedented accuracy [Dornheim et al., Nat. Commun. 16, 5103 (2025)]. Yet, it comes at the cost of performing an additional 10-20 simulations, which, in combination with the required small error bars, can pose a serious practical limitation. Here, we remove this bottleneck by presenting a new reweighting estimator, which allows the study of the full ξ-dependence from a single path integral Monte Carlo (PIMC) simulation. This is demonstrated for various observables of the uniform electron gas and also warm dense beryllium. We expect our study to be useful for future PIMC simulations of Fermi systems, including ultracold atoms, electrons in quantum dots, and warm dense quantum plasmas.
The ab initio path integral Monte Carlo (PIMC) method is one of the most successful methods in statistical physics, quantum chemistry and related fields, but its application to quantum degenerate Fermi systems is severely hampered by an exponential computational bottleneck: the notorious fermion sign problem. Very recently, Xiong and Xiong [J. Chem. Phys. 157, 094112 (2022)] have suggested to partially circumvent the sign problem by carrying out simulations of fictitious systems guided by an interpolating continuous variable ξ ∈ [-1, 1], with the physical Fermi- and Bose-statistics corresponding to ξ = -1 and ξ = 1. It has been proposed that information about the fermionic limit might be obtained by calculations within the bosonic sector ξ > 0 combined with an extrapolation throughout the fermionic sector ξ < 0, essentially bypassing the sign problem. Here, we show how the inclusion of the artificial parameter ξ can be interpreted as an effective penalty on the formation of permutation cycles in the PIMC simulation. We demonstrate that the proposed extrapolation method breaks down for moderate to high quantum degeneracy. Instead, the method constitutes a valuable tool for the description of large Fermi-systems of weak quantum degeneracy. This is demonstrated for electrons in a 2D harmonic trap and for the uniform electron gas (UEG), where we find excellent agreement (∼0.5%) with exact configuration PIMC results in the high-density regime while attaining a speed-up exceeding 11 orders of magnitude. Finally, we extend the idea beyond the energy and analyze the radial density distribution (2D trap), as well as the static structure factor and imaginary-time density-density correlation function (UEG).
The ab initio thermodynamic simulation of correlated Fermi systems is of central importance for many applications, such as warm dense matter, electrons in quantum dots, and ultracold atoms. Unfortunately, path integral Monte Carlo (PIMC) simulations of fermions are severely restricted by the notorious fermion sign problem (FSP). In this paper, we present a hands-on discussion of the FSP and investigate in detail its manifestation with respect to temperature, system size, interaction-strength and -type, and the dimensionality of the system. Moreover, we analyze the probability distribution of fermionic expectation values, which can be non-Gaussian and fat-tailed when the FSP is severe. As a practical application, we consider electrons and dipolar atoms in a harmonic confinement, and the uniform electron gas in the warm dense matter regime. In addition, we provide extensive PIMC data, which can be used as a reference for the development of new methods and as a benchmark for approximations.
The exchange antisymmetry of identical fermions leads to an exponential computational bottleneck in ab initio simulations, known as the fermion sign problem. The thermodynamic approach of fictitious identical particles (Y. Xiong and H. Xiong, J. Chem. Phys. 157, 094112 (2022)) provides an efficient and accurate means to simulate some fermionic systems by overcoming the fermion sign problem. This method has been significantly promoted and used by National Ignition Facilities for the ab initio simulations and is believed to have wide application prospects in warm dense quantum matter (T. Dornheim et al., arXiv: 2402.19113 (2023)). By utilizing the fictitious identical particles in the bosonic regime and constant energy extrapolation method (Y. Xiong and H. Xiong, Phys. Rev. E 107, 055308 (2023); T. Morresi and G. Garberoglio, Phys. Rev. B 111, 014521 (2025)), there are promising results in simulating the energy of strongly quantum degenerate fermionic systems. The previous works mainly concern the energy of Fermi systems or only consider situations of weak quantum degeneracy. In this study, we extend the concept of the constant energy extrapolation method and demonstrate the potential of the constant density extrapolation method to accurately simulate the density distribution of fermionic systems in strongly quantum degenerate conditions. Furthermore, based on the energy derived from the constant energy extrapolation method, we present simulation results for the entropy of fermions.
Recently, fictitious identical particles have provided a promising way to overcome the fermion sign problem and have been used in path integral Monte Carlo to accurately simulate warm dense matter with up to 1000 electrons [T. Dornheim et al., J. Phys. Chem. Lett. 15, 1305 (2024)1948-718510.1021/acs.jpclett.3c03193]. The inclusion of fictitious identical particles in path integral molecular dynamics (PIMD) can provide another way to simulate fermion systems. In a recent paper [Y. M. Y. Feldman et al., J. Chem. Phys. 159, 154107 (2023)0021-960610.1063/5.0173749], Feldman and Hirshberg improved the recursive formula for PIMD of identical bosons, significantly reducing the computational complexity. In this paper, we extend this latest recursive formula for bosons to PIMD of fictitious identical particles to improve the efficiency of simulating fermion systems. We also provide the virial estimator for calculating energy of fictitious identical particles by using the new recursive technique, which can help suppress fluctuations in energy calculation. As an example, we use the quadratic scaling PIMD for fictitious identical particles to study the simulation of hundreds and even thousands of fermions in a two-dimensional periodic potential, in the hope of providing a simulation tool for the large-scale Fermi-Hubbard model and other strongly correlated fermion systems, such as the simulation of ultracold fermionic gases in optical lattices and the extended Fermi-Hubbard model using quantum dots.
Accurate thermodynamic simulations of correlated fermions using path integral Monte Carlo (PIMC) methods are of paramount importance for many applications such as the description of ultracold atoms, electrons in quantum dots, and warm-dense matter. The main obstacle is the fermion sign problem (FSP), which leads to an exponential increase in computation time both with an increase in the system size and with a decrease in the temperature. Very recently, Hirshberg et al. [J. Chem. Phys. 152, 171102 (2020)] have proposed to alleviate the FSP based on the Bogoliubov inequality. In the present work, we extend this approach by adding a parameter that controls the perturbation, allowing for an extrapolation to the exact result. In this way, we can also use thermodynamic integration to obtain an improved estimate of the fermionic energy. As a test system, we choose electrons in 2D and 3D quantum dots and find in some cases a speed-up exceeding 106, as compared to standard PIMC, while retaining a relative accuracy of ∼0.1%. Our approach is quite general and can readily be adapted to other simulation methods.