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

Normal liquid He 3 studied by path-integral Monte Carlo with a parametrized partition function

Tommaso Morresi, Giovanni Garberoglio

Published Jan 31, 2025• 15 citations
arXivHTMLPDF

Abstract

We compute the energy per particle of normal liquid $^{3}\mathrm{He}$ in the temperature range $0.15--2$ K using path-integral Monte Carlo simulations, leveraging a recently proposed method to overcome the sign problem---a long-standing challenge in many-body fermionic simulations. This approach is based on introducing a parameter $\ensuremath{\xi}$ into the partition function, which allows a generalization from bosons ($\ensuremath{\xi}=1$) to fermions ($\ensuremath{\xi}=\ensuremath{-}1$). By simulating systems with $\ensuremath{\xi}\ensuremath{\ge}0$, where the sign problem is absent, one can then extrapolate to the fermionic case at $\ensuremath{\xi}=\ensuremath{-}1$. Guided by an independent-particle model that uncovers nonanalytic behavior due to the superfluid transition, which is moderated by finite-size effects, we develop a tailored extrapolation strategy for liquid $^{3}\mathrm{He}$ that departs from the extrapolation schemes shown to be accurate in those cases where quantum degeneracy effects are weak, and enables accurate results in the presence of Bose-Einstein condensation and superfluidity for $\ensuremath{\xi}>0$. Our approach extends the previously proposed framework and yields energy per particle values in good agreement with experimental data.

DOI
10.1103/physrevb.111.014521

Related papers

Path integral Monte Carlo approach to the structural properties and collective excitations of liquid $$^3{\text {He}}$$ without fixed nodes

Tobias Dornheim, Zhandos A. Moldabekov, Jan Vorberger +1 more

Due to its nature as a strongly correlated quantum liquid, ultracold helium is characterized by the nontrivial interplay of different physical effects. Bosonic [Formula: see text] exhibits superfluidity and Bose-Einstein condensation. Its physical properties have been accurately determined on the basis of ab initio path integral Monte Carlo (PIMC) simulations. In contrast, the corresponding theoretical description of fermionic [Formula: see text] is severely hampered by the notorious fermion sign problem, and previous PIMC results have been derived by introducing the uncontrolled fixed-node approximation. In this work, we present extensive new PIMC simulations of normal liquid [Formula: see text] without any nodal constraints. This allows us to to unambiguously quantify the impact of Fermi statistics and to study the effects of temperature on different physical properties like the static structure factor [Formula: see text], the momentum distribution [Formula: see text], and the static density response function [Formula: see text]. In addition, the dynamic structure factor [Formula: see text] is rigorously reconstructed from imaginary-time PIMC data. From simulations of [Formula: see text], we derived the familiar phonon-maxon-roton dispersion function that is well-known for [Formula: see text] and has been reported previously for two-dimensional [Formula: see text] films (Nature 483:576-579 (2012)). The comparison of our new results for both [Formula: see text] and [Formula: see text] with neutron scattering measurements reveals an excellent agreement between theory and experiment.

Scientific Reports•Jan 13, 2022•40 citations

Fermionic Free Energies from Ab Initio Path Integral Monte Carlo Simulations of Fictitious Identical Particles

Tobias Dornheim, Zhandos A. Moldabekov, Sebastian Schwalbe +2 more

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.

Journal of Chemical Theory and Computation•Jul 24, 2025•12 citations

Path Integral Monte Carlo and Bose-Einstein condensation in quantum fluids and solids

R. Rota

Several microscopic theories point out that Bose-Einstein condensation (BEC), i.e., a macroscopic occupation of the lowest energy single particle state in many-boson systems, may appear also in quantum fluids and solids and that it is at the origin of the phenomenon of superfluidity. Nevertheless, the connection between BEC and superfluidity is still matter of debate, since the experimental evidences indicating a non zero condensate fraction in superfluid helium are only indirect. In the theoretical study of BEC in quantum fluids and solids, perturbative approaches are useless because of the strong correlations between the atoms, arising both from the interatomic potential and from the quantum nature of the system. Microscopic Quantum Monte Carlo simulations provide a reliable description of these systems. In particular, the Path Integral Monte Carlo (PIMC) method is very suitable for this purpose. This method is able to provide exact results for the properties of the quantum system, both at zero and finite temperature, only with the definition of the Hamiltonian and of the symmetry properties of the system, giving an easy picture for superfluidity and BEC in many-boson systems. In this thesis, we apply PIMC methods to the study of several quantum fluids and solids. We describe in detail all the features of PIMC, from the sampling methods to the estimators of the physical properties. We present also the most recent techniques, such as the high-order approximations for the thermal density matrix and the worm algorithm, used in PIMC to provide reliable simulations. We study the liquid phase of condensed 4He, providing unbiased estimations of the one-body density matrix g1(r). We analyze the model for g1(r) used to fit the experimental data, highlighting its merits and its faults. In particular we see that, even if it presents some difficulties in the description of the overall behavior of g1(r), it can provide an accurate estimation of the kinetic energy K and of the condensate fraction n0 of the system. Furthermore, we show that our results for n0 as a function of the pressure are in a good agreement with the most recent experimental results. The study of the solid phase of 4He is the most significant part of this thesis. The recent observation of non classical rotational inertia (NCRI) effects in solid helium has generated big interest in the study of an eventual supersolid phase, characterized at the same time by crystalline order and superfluidity. Nevertheless, until now it has been impossible to give a theoretical model able to describe all the experimental evidences. In this work, we perform PIMC simulations of 4He at high densities, according to different microscopic configurations of the atoms. In commensurate crystals we see that BEC does not appear, our model being able to reproduce the momentum distribution obtained form neutron scattering experiments. In a crystal with vacancies, we have been able to see a transition to a superfluid phase at temperatures in agreement with experimental results if the vacancy concentration is low enough. In amorphous solids, superfluid effects are enhanced but appear at temperatures higher than the experimental estimation for the transition temperature. Finally, we study also metastable disordered configurations in molecular para-hydrogen at low temperature. The aim of this study is to investigate if a Bose liquid other than helium can display superfluidity. Choosing accurately a ¿quantum liquid¿ initial configuration and the dimensions of the simulation box, we have been able to frustrate the formation of the crystal and to calculate the temperature dependence of the superfluid density, showing a transition to a superfluid phase at temperatures close to 1 K.

•Dec 20, 2011•5 citations

On the thermodynamic properties of fictitious identical particles and the application to fermion sign problem

Yunuo Xiong, Hongwei Xiong

By generalizing the recently developed path integral molecular dynamics for identical bosons and fermions, we consider the finite-temperature thermodynamic properties of fictitious identical particles with a real parameter ξ interpolating continuously between bosons (ξ = 1) and fermions (ξ = -1). Through general analysis and numerical experiments, we find that the average energy may have good analytical properties as a function of this real parameter ξ, which provides the chance to calculate the thermodynamical properties of identical fermions by extrapolation with a simple polynomial function after accurately calculating the thermodynamic properties of the fictitious particles for ξ ≥ 0. Using several examples, it is shown that our method can efficiently give accurate energy values for finite-temperature fermionic systems. Our work provides a chance to circumvent the fermion sign problem for some quantum systems.

arXiv:2206.08341•The Journal of Chemical Physics•Aug 8, 2022•43 citations

Fermionic physics fromab initiopath integral Monte Carlo simulations of fictitious identical particles

Tobias Dornheim, Panagiotis B. Tolias, Simon Groth +3 more

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).

arXiv:2308.06071•The Journal of Chemical Physics•Oct 27, 2023•42 citations

Parametrized path integral formulation for large fermion systems

Yunuo, Xiong, Hongwei Xiong

The exchange antisymmetry between identical fermions gives rise to the well known fermion sign problem, in the form of large cancellation between positive and negative contribution to the partition function, making any simulation methods which directly sample this partition function exponentially difficult to converge. In this work, we employ path integral molecular dynamics (PIMD) and build upon the recently discovered fictitious particle model to investigate the fermion sign problem further. We consider the validity and invalidity condition for the method of parametrized path integral formulation of the partition function and extrapolation to circumvent the fermion sign problem. For the valid region of our method, our simulation shows that we may give accurate prediction of the energy for large fermion systems, which is much beyond the capability of the direct sampling in the traditional method. In particular, we find and verify a simple universal relation for high temperature noninteracting particles or strongly repulsive interacting particles at low temperatures.

arXiv:2208.13777•Aug 29, 2022•1 citation
We compute the energy per particle of normal liquid 3He^{3}\mathrm{He}3He in the temperature range 0.15−−20.15--20.15−−2 K using path-integral Monte Carlo simulations, leveraging a recently proposed method to overcome the sign problem---a long-standing challenge in many-body fermionic simulations. This approach is based on introducing a parameter \ensuremathξ\ensuremath{\xi}\ensuremathξ into the partition function, which allows a generalization from bosons (\ensuremathξ=1\ensuremath{\xi}=1\ensuremathξ=1) to fermions (\ensuremathξ=\ensuremath−1\ensuremath{\xi}=\ensuremath{-}1\ensuremathξ=\ensuremath−1). By simulating systems with \ensuremathξ\ensuremath≥0\ensuremath{\xi}\ensuremath{\ge}0\ensuremathξ\ensuremath≥0, where the sign problem is absent, one can then extrapolate to the fermionic case at \ensuremathξ=\ensuremath−1\ensuremath{\xi}=\ensuremath{-}1\ensuremathξ=\ensuremath−1. Guided by an independent-particle model that uncovers nonanalytic behavior due to the superfluid transition, which is moderated by finite-size effects, we develop a tailored extrapolation strategy for liquid 3He^{3}\mathrm{He}3He that departs from the extrapolation schemes shown to be accurate in those cases where quantum degeneracy effects are weak, and enables accurate results in the presence of Bose-Einstein condensation and superfluidity for \ensuremathξ>0\ensuremath{\xi}>0\ensuremathξ>0. Our approach extends the previously proposed framework and yields energy per particle values in good agreement with experimental data.
Path integral Monte Carlo approach to the structural properties and collective excitations of liquid 3He^3{\text {He}}3He without fixed nodes
Due to its nature as a strongly correlated quantum liquid, ultracold helium is characterized by the nontrivial interplay of different physical effects. Bosonic [Formula: see text] exhibits superfluidity and Bose-Einstein condensation. Its physical properties have been accurately determined on the basis of ab initio path integral Monte Carlo (PIMC) simulations. In contrast, the corresponding theoretical description of fermionic [Formula: see text] is severely hampered by the notorious fermion sign problem, and previous PIMC results have been derived by introducing the uncontrolled fixed-node approximation. In this work, we present extensive new PIMC simulations of normal liquid [Formula: see text] without any nodal constraints. This allows us to to unambiguously quantify the impact of Fermi statistics and to study the effects of temperature on different physical properties like the static structure factor [Formula: see text], the momentum distribution [Formula: see text], and the static density response function [Formula: see text]. In addition, the dynamic structure factor [Formula: see text] is rigorously reconstructed from imaginary-time PIMC data. From simulations of [Formula: see text], we derived the familiar phonon-maxon-roton dispersion function that is well-known for [Formula: see text] and has been reported previously for two-dimensional [Formula: see text] films (Nature 483:576-579 (2012)). The comparison of our new results for both [Formula: see text] and [Formula: see text] with neutron scattering measurements reveals an excellent agreement between theory and experiment.
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.
Several microscopic theories point out that Bose-Einstein condensation (BEC), i.e., a macroscopic occupation of the lowest energy single particle state in many-boson systems, may appear also in quantum fluids and solids and that it is at the origin of the phenomenon of superfluidity. Nevertheless, the connection between BEC and superfluidity is still matter of debate, since the experimental evidences indicating a non zero condensate fraction in superfluid helium are only indirect. In the theoretical study of BEC in quantum fluids and solids, perturbative approaches are useless because of the strong correlations between the atoms, arising both from the interatomic potential and from the quantum nature of the system. Microscopic Quantum Monte Carlo simulations provide a reliable description of these systems. In particular, the Path Integral Monte Carlo (PIMC) method is very suitable for this purpose. This method is able to provide exact results for the properties of the quantum system, both at zero and finite temperature, only with the definition of the Hamiltonian and of the symmetry properties of the system, giving an easy picture for superfluidity and BEC in many-boson systems. In this thesis, we apply PIMC methods to the study of several quantum fluids and solids. We describe in detail all the features of PIMC, from the sampling methods to the estimators of the physical properties. We present also the most recent techniques, such as the high-order approximations for the thermal density matrix and the worm algorithm, used in PIMC to provide reliable simulations. We study the liquid phase of condensed 4He, providing unbiased estimations of the one-body density matrix g1(r). We analyze the model for g1(r) used to fit the experimental data, highlighting its merits and its faults. In particular we see that, even if it presents some difficulties in the description of the overall behavior of g1(r), it can provide an accurate estimation of the kinetic energy K and of the condensate fraction n0 of the system. Furthermore, we show that our results for n0 as a function of the pressure are in a good agreement with the most recent experimental results. The study of the solid phase of 4He is the most significant part of this thesis. The recent observation of non classical rotational inertia (NCRI) effects in solid helium has generated big interest in the study of an eventual supersolid phase, characterized at the same time by crystalline order and superfluidity. Nevertheless, until now it has been impossible to give a theoretical model able to describe all the experimental evidences. In this work, we perform PIMC simulations of 4He at high densities, according to different microscopic configurations of the atoms. In commensurate crystals we see that BEC does not appear, our model being able to reproduce the momentum distribution obtained form neutron scattering experiments. In a crystal with vacancies, we have been able to see a transition to a superfluid phase at temperatures in agreement with experimental results if the vacancy concentration is low enough. In amorphous solids, superfluid effects are enhanced but appear at temperatures higher than the experimental estimation for the transition temperature. Finally, we study also metastable disordered configurations in molecular para-hydrogen at low temperature. The aim of this study is to investigate if a Bose liquid other than helium can display superfluidity. Choosing accurately a ¿quantum liquid¿ initial configuration and the dimensions of the simulation box, we have been able to frustrate the formation of the crystal and to calculate the temperature dependence of the superfluid density, showing a transition to a superfluid phase at temperatures close to 1 K.
By generalizing the recently developed path integral molecular dynamics for identical bosons and fermions, we consider the finite-temperature thermodynamic properties of fictitious identical particles with a real parameter ξ interpolating continuously between bosons (ξ = 1) and fermions (ξ = -1). Through general analysis and numerical experiments, we find that the average energy may have good analytical properties as a function of this real parameter ξ, which provides the chance to calculate the thermodynamical properties of identical fermions by extrapolation with a simple polynomial function after accurately calculating the thermodynamic properties of the fictitious particles for ξ ≥ 0. Using several examples, it is shown that our method can efficiently give accurate energy values for finite-temperature fermionic systems. Our work provides a chance to circumvent the fermion sign problem for some quantum systems.
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 exchange antisymmetry between identical fermions gives rise to the well known fermion sign problem, in the form of large cancellation between positive and negative contribution to the partition function, making any simulation methods which directly sample this partition function exponentially difficult to converge. In this work, we employ path integral molecular dynamics (PIMD) and build upon the recently discovered fictitious particle model to investigate the fermion sign problem further. We consider the validity and invalidity condition for the method of parametrized path integral formulation of the partition function and extrapolation to circumvent the fermion sign problem. For the valid region of our method, our simulation shows that we may give accurate prediction of the energy for large fermion systems, which is much beyond the capability of the direct sampling in the traditional method. In particular, we find and verify a simple universal relation for high temperature noninteracting particles or strongly repulsive interacting particles at low temperatures.