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:2211.12480 · Physical review. E

Thermodynamics of fermions at any temperature based on parametrized partition function

Yunuo Xiong, Hongwei Xiong

Published May 25, 2023• 23 citations
arXivHTMLPDF

Abstract

In this work, we study the recently developed parametrized partition function formulation and show how we can infer the thermodynamic properties of fermions based on numerical simulation of bosons and distinguishable particles at various temperatures. In particular, we show that in the three-dimensional space defined by energy, temperature, and the parameter characterizing parametrized partition function, we can map the energies of bosons and distinguishable particles to fermionic energies through constant-energy contours. We apply this idea to both noninteracting and interacting Fermi systems and show it is possible to infer the fermionic energies at all temperatures, thus providing a practical and efficient approach to obtain thermodynamic properties of Fermi systems with numerical simulation. As an example, we present energies and heat capacities for 10 noninteracting fermions and 10 interacting fermions and show good agreement with the analytical result for the noninteracting case.

DOI
10.1103/physreve.107.055308

Related papers

On the Thermodynamics of Fermions at any Temperature based on Parametrized Partition Function

Hongwei Xiong, Yunuo Xiong

In this work we study the recently developed parametrized partition function formulation and show how we can infer the thermodynamic properties of fermions based on numerical simulation of bosons and distinguishable particles at various temperatures. In particular, we show that in the three dimensional space defined by energy, temperature and the parameter characterizing parametrized partition function, we can map the energies of bosons and distinguishable particles to fermionic energies through constant-energy contours. We apply this idea to both noninteracting and interacting Fermi systems and show it is possible to infer the fermionic energies at all temperatures, thus providing a practical and efficient approach to obtain thermodynamic properties of large fermion systems with numerical simulation. As an example, we present energies for up to 50 noninteracting fermions and up to 20 interacting fermions at all temperatures and show good agreement with the analytical result for noninteracting case.

•Nov 22, 2022•1 citation

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

Numerical simulation of two-component attractive Fermi gases based on parametrized partition function

Yunuo Xiong, Hongwei Xiong

The zero-temperature and finite-temperature thermodynamics of two-component Fermi gases with finite-range attractive interaction suffer from fermion sign problem, which seems like an insurmountable problem in exact numerical simulations. In a recent work, we find a reliable method to simulate the thermodynamic properties of single-component Fermi gases for both noninteracting and repulsively interacting cases based on the method of parametrized partition function and the $ξ_E$ curve of constant energy. In the present work, this method is generalized to two-component Fermi gases with finite-range attractive interaction, which shows clearly that our method has good chance to apply to various Fermi systems. From the simulated heat capacity, we find a peak at the temperature below the Fermi temperature which implies the pairing of fermions with different spin. At high temperature, the simulated heat capacity approaches the classical value. The reasonable result in this work validates the application of our method to attractive cases, which implies a wide range of applications, from nuclear physics, BCS-BEC crossover, superconductivity, to neutron star, etc..

arXiv:2212.00580•Nov 30, 2022•0 citations

Many-body formalism for fermions: The partition function

D. K. Watson

The partition function, a fundamental tenet in statistical thermodynamics, contains in principle all thermodynamic information about a system. It encapsulates both microscopic information through the quantum energy levels and statistical information from the partitioning of the particles among the available energy levels. For identical particles, this statistical accounting is complicated by the symmetry requirements of the allowed quantum states. In particular, for Fermi systems, the enforcement of the Pauli principle is typically a numerically demanding task, responsible for much of the cost of the calculations. The interplay of these three elements---the structure of the many-body spectrum, the statistical partitioning of the $N$ particles among the available levels, and the enforcement of the Pauli principle---drives the behavior of mesoscopic and macroscopic Fermi systems. In this paper, we develop an approach for the determination of the partition function, a numerically difficult task, for systems of strongly interacting identical fermions and apply it to a model system of harmonically confined, harmonically interacting fermions. This approach uses a recently introduced many-body method that is an extension of the symmetry-invariant perturbation method (SPT) originally developed for bosons. It uses group theory and graphical techniques to avoid the heavy computational demands of conventional many-body methods which typically scale exponentially with the number of particles. The SPT application of the Pauli principle is trivial to implement since it is done ``on paper'' by imposing restrictions on the normal-mode quantum numbers at first order in the perturbation. The method is applied through first order and represents an extension of the SPT method to excited states. Our method of determining the partition function and various thermodynamic quantities is accurate and efficient and has the potential to yield interesting insight into the role played by the Pauli principle and the influence of large degeneracies on the emergence of the thermodynamic behavior of large-$N$ systems.

Physical Review A•Sep 1, 2017•7 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

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
In this work, we study the recently developed parametrized partition function formulation and show how we can infer the thermodynamic properties of fermions based on numerical simulation of bosons and distinguishable particles at various temperatures. In particular, we show that in the three-dimensional space defined by energy, temperature, and the parameter characterizing parametrized partition function, we can map the energies of bosons and distinguishable particles to fermionic energies through constant-energy contours. We apply this idea to both noninteracting and interacting Fermi systems and show it is possible to infer the fermionic energies at all temperatures, thus providing a practical and efficient approach to obtain thermodynamic properties of Fermi systems with numerical simulation. As an example, we present energies and heat capacities for 10 noninteracting fermions and 10 interacting fermions and show good agreement with the analytical result for the noninteracting case.
In this work we study the recently developed parametrized partition function formulation and show how we can infer the thermodynamic properties of fermions based on numerical simulation of bosons and distinguishable particles at various temperatures. In particular, we show that in the three dimensional space defined by energy, temperature and the parameter characterizing parametrized partition function, we can map the energies of bosons and distinguishable particles to fermionic energies through constant-energy contours. We apply this idea to both noninteracting and interacting Fermi systems and show it is possible to infer the fermionic energies at all temperatures, thus providing a practical and efficient approach to obtain thermodynamic properties of large fermion systems with numerical simulation. As an example, we present energies for up to 50 noninteracting fermions and up to 20 interacting fermions at all temperatures and show good agreement with the analytical result for noninteracting case.
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 zero-temperature and finite-temperature thermodynamics of two-component Fermi gases with finite-range attractive interaction suffer from fermion sign problem, which seems like an insurmountable problem in exact numerical simulations. In a recent work, we find a reliable method to simulate the thermodynamic properties of single-component Fermi gases for both noninteracting and repulsively interacting cases based on the method of parametrized partition function and the ξEξ_EξE​ curve of constant energy. In the present work, this method is generalized to two-component Fermi gases with finite-range attractive interaction, which shows clearly that our method has good chance to apply to various Fermi systems. From the simulated heat capacity, we find a peak at the temperature below the Fermi temperature which implies the pairing of fermions with different spin. At high temperature, the simulated heat capacity approaches the classical value. The reasonable result in this work validates the application of our method to attractive cases, which implies a wide range of applications, from nuclear physics, BCS-BEC crossover, superconductivity, to neutron star, etc..
The partition function, a fundamental tenet in statistical thermodynamics, contains in principle all thermodynamic information about a system. It encapsulates both microscopic information through the quantum energy levels and statistical information from the partitioning of the particles among the available energy levels. For identical particles, this statistical accounting is complicated by the symmetry requirements of the allowed quantum states. In particular, for Fermi systems, the enforcement of the Pauli principle is typically a numerically demanding task, responsible for much of the cost of the calculations. The interplay of these three elements---the structure of the many-body spectrum, the statistical partitioning of the NNN particles among the available levels, and the enforcement of the Pauli principle---drives the behavior of mesoscopic and macroscopic Fermi systems. In this paper, we develop an approach for the determination of the partition function, a numerically difficult task, for systems of strongly interacting identical fermions and apply it to a model system of harmonically confined, harmonically interacting fermions. This approach uses a recently introduced many-body method that is an extension of the symmetry-invariant perturbation method (SPT) originally developed for bosons. It uses group theory and graphical techniques to avoid the heavy computational demands of conventional many-body methods which typically scale exponentially with the number of particles. The SPT application of the Pauli principle is trivial to implement since it is done ``on paper'' by imposing restrictions on the normal-mode quantum numbers at first order in the perturbation. The method is applied through first order and represents an extension of the SPT method to excited states. Our method of determining the partition function and various thermodynamic quantities is accurate and efficient and has the potential to yield interesting insight into the role played by the Pauli principle and the influence of large degeneracies on the emergence of the thermodynamic behavior of large-NNN systems.
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.
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.