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:2212.00580

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

Yunuo Xiong, Hongwei Xiong

Published Nov 30, 2022• 0 citations
arXivHTMLPDF

Abstract

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

DOI
10.48550/arxiv.2212.00580

Related papers

Critical Temperature and Thermodynamics of Attractive Fermions at Unitarity

Evgeni Burovski, Nikolay V. Prokof'ev, Boris V. Svistunov +1 more

The unitarity regime of the BCS-BEC crossover can be realized by diluting a system of two-component lattice fermions with an on-site attractive interaction. We perform a systematic-error-free finite-temperature simulation of this system by diagrammatic determinant Monte Carlo method. The critical temperature in units of Fermi energy is found to be T(C)/epsilonF=0.152(7). We also report the behavior of the thermodynamic functions, and discuss the issues of thermometry of ultracold Fermi gases.

arXiv:cond-mat/0602224•Physical Review Letters•Apr 26, 2006•253 citations

Thermodynamic and hydrodynamic behaviour of interacting Fermi gases

Olga Goulko

Fermionic matter is ubiquitous in nature, from the electrons in metals and semiconductors or the neutrons in the inner crust of neutron stars, to gases of fermionic atoms, like 40K or 6Li that can be created and studied under laboratory conditions. It is especially interesting to study these systems at very low temperatures, where we enter the world of quantum mechanical phenomena. Due to the Fermi-Dirac statistics, a dilute system of spin-polarised fermions exhibits no interactions and can be viewed as an ideal Fermi gas. However, interactions play a crucial role for fermions of several spin species. This thesis addresses several questions concerning interacting Fermi gases, in particular the transition between the normal and the superfluid phase and dynamical properties at higher temperatures. First we will look at the unitary Fermi gas: a two-component system of fermions interacting with divergent scattering length. This system is particularly interesting as it exhibits universal behaviour. Due to the strong interactions perturbation theory is inapplicable and no exact theoretical description is available. I will describe the Determinant Diagrammatic Monte Carlo algorithm with which the unitary Fermi gas can be studied from first principles. This algorithm fails in the presence of a spin imbalance (unequal number of particles in the two components) due to a sign problem. I will show how to apply reweighting techniques to generalise the algorithm to the imbalanced case, and present results for the critical temperature and other thermodynamic observables at the critical point, namely the chemical potential, the energy per particle and the contact density. These are the first numerical results for the imbalanced unitary Fermi gas at finite temperature. I will also show how temperatures beyond the critical point can be accessed and present results for the equation of state and the temperature dependence of the contact density. At sufficiently high temperatures a semiclassical description captures all relevant physical features of the system. The dynamics of an interacting Fermi gas can then be studied via a numerical simulation of the Boltzmann equation. I will describe such a numerical setup and apply it to study the collision of two spin-polarised fermionic clouds. When the two components are separated in an elongated harmonic trap and then released, they collide and for sufficiently strong interactions can bounce off each other several times. I will discuss the different types of the qualitative behaviour, show how they can be interpreted in terms of the equilibrium properties of the system, and explain how they relate to the coupling between different excitation modes. I will also demonstrate how transport coefficients, for instance the spin drag, can be extracted from the numerical data.

•Jan 10, 2012•2 citations

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

Thermodynamics of fermions at any temperature based on parametrized partition function

Yunuo Xiong, Hongwei 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 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.

arXiv:2211.12480•Physical review. E•May 25, 2023•23 citations

A complex Langevin approach to ultracold fermions

J.E. Drut, Lukas Rammelmüller, J. Braun

The theoretical treatment of Fermi systems consisting of particles with unequal masses is challenging. Even in one spatial dimension analytic solutions are limited to special configurations and numerical progress with Monte Carlo simulations is hindered by the sign-problem. To circumvent this issue, we exploit the Complex Langevin approach and study one-dimensional mass-imbalanced two-component Fermi gases with attractive and repulsive interactions. We find perfect agreement with results obtained by other methods in a range of parameter space. Promisingly, our approach is not limited to the specific model presented here and can easily be extended to finite spin polarization and, most notably, can also be applied in higher dimensions.

•Nov 9, 2022•0 citations

Pairing and molecule formation along the BCS-BEC crossover for finite range potentials

E Neri, Santiago F. Caballero-Benítez, Vı́ctor Romero-Rochı́n +1 more

Abstract We address the thermodynamics of the zero temperature BCS-BEC crossover of a balanced two component mixture of a finite range interacting fermions within a mean field approach. For the analysis we consider two wide types of finite range potentials describing the interaction between different Fermi species: purely attractive, as the square well, the exponential and the Yukawa potentials, and others having both a repulsive wall at the origin as well as a finite attractive tail, as Van der Waals and dipolar potentials. We study the energy gap, the chemical potential and the pressure as functions of the scattering length for all the interacting potentials considered. The system spatial structure is further analyzed via the the pair functions associated to each potential. We find that for those with pure attractive contributions the crossover can be recognized either, when the chemical potential changes its sign or when the scattering length diverges. In contrast, the pair functions linked to potentials with both repulsive and attractive terms, show wide crossover regions without sharp distinctions at the sign change of the chemical potential nor at unitarity.

Physica Scripta•Nov 12, 2019•6 citations
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..
Fermionic matter is ubiquitous in nature, from the electrons in metals and semiconductors or the neutrons in the inner crust of neutron stars, to gases of fermionic atoms, like 40K or 6Li that can be created and studied under laboratory conditions. It is especially interesting to study these systems at very low temperatures, where we enter the world of quantum mechanical phenomena. Due to the Fermi-Dirac statistics, a dilute system of spin-polarised fermions exhibits no interactions and can be viewed as an ideal Fermi gas. However, interactions play a crucial role for fermions of several spin species. This thesis addresses several questions concerning interacting Fermi gases, in particular the transition between the normal and the superfluid phase and dynamical properties at higher temperatures. First we will look at the unitary Fermi gas: a two-component system of fermions interacting with divergent scattering length. This system is particularly interesting as it exhibits universal behaviour. Due to the strong interactions perturbation theory is inapplicable and no exact theoretical description is available. I will describe the Determinant Diagrammatic Monte Carlo algorithm with which the unitary Fermi gas can be studied from first principles. This algorithm fails in the presence of a spin imbalance (unequal number of particles in the two components) due to a sign problem. I will show how to apply reweighting techniques to generalise the algorithm to the imbalanced case, and present results for the critical temperature and other thermodynamic observables at the critical point, namely the chemical potential, the energy per particle and the contact density. These are the first numerical results for the imbalanced unitary Fermi gas at finite temperature. I will also show how temperatures beyond the critical point can be accessed and present results for the equation of state and the temperature dependence of the contact density. At sufficiently high temperatures a semiclassical description captures all relevant physical features of the system. The dynamics of an interacting Fermi gas can then be studied via a numerical simulation of the Boltzmann equation. I will describe such a numerical setup and apply it to study the collision of two spin-polarised fermionic clouds. When the two components are separated in an elongated harmonic trap and then released, they collide and for sufficiently strong interactions can bounce off each other several times. I will discuss the different types of the qualitative behaviour, show how they can be interpreted in terms of the equilibrium properties of the system, and explain how they relate to the coupling between different excitation modes. I will also demonstrate how transport coefficients, for instance the spin drag, can be extracted from the numerical data.
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.
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.
The theoretical treatment of Fermi systems consisting of particles with unequal masses is challenging. Even in one spatial dimension analytic solutions are limited to special configurations and numerical progress with Monte Carlo simulations is hindered by the sign-problem. To circumvent this issue, we exploit the Complex Langevin approach and study one-dimensional mass-imbalanced two-component Fermi gases with attractive and repulsive interactions. We find perfect agreement with results obtained by other methods in a range of parameter space. Promisingly, our approach is not limited to the specific model presented here and can easily be extended to finite spin polarization and, most notably, can also be applied in higher dimensions.
Abstract We address the thermodynamics of the zero temperature BCS-BEC crossover of a balanced two component mixture of a finite range interacting fermions within a mean field approach. For the analysis we consider two wide types of finite range potentials describing the interaction between different Fermi species: purely attractive, as the square well, the exponential and the Yukawa potentials, and others having both a repulsive wall at the origin as well as a finite attractive tail, as Van der Waals and dipolar potentials. We study the energy gap, the chemical potential and the pressure as functions of the scattering length for all the interacting potentials considered. The system spatial structure is further analyzed via the the pair functions associated to each potential. We find that for those with pure attractive contributions the crossover can be recognized either, when the chemical potential changes its sign or when the scattering length diverges. In contrast, the pair functions linked to potentials with both repulsive and attractive terms, show wide crossover regions without sharp distinctions at the sign change of the chemical potential nor at unitarity.