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

Simulation of an inhomogeneous Fermi gas through the BCS-BEC crossover

Rocío Jáuregui, R. Paredes, L. Rosales-Zárate, G. Toledo Sánchez

Published Mar 4, 2008• 0 citations
arXivHTMLPDF

Abstract

We perform a variational quantum Monte Carlo simulation of the transition from a Bardeen-Cooper-Schrieffer superfluid (BCS) to a Bose-Einstein condensate (BEC) at zero temperature. The model Hamiltonian involves an attractive short range two body interaction and the atoms number $2N =330$ is chosen so that, in the non-interacting limit, the ground state function corresponds to a closed shell configuration. The system is then characterized by the s-wave scattering length $a$ of the two-particle collisions in the gas, which is varied from negative to positive values, and the Fermi wave number $k_F$. Based on an extensive analysis of the s-wave two-body problem, one parameter variational many-body wave functions are proposed to describe the ground state of the interacting Fermi gas from BCS to BEC states. We exploit properties of antisymmetrized many-body functions to develop efficient techniques that permit variational calculations for a large number of particles. It is shown that a virial relation between the energy per particle and the trapping energy is approximately valid for $-0.1<1/k_Fa<3.4$. The influence of the harmonic trap and the interaction potential as exhibited in two-body correlation functions is also analyzed.

DOI
10.48550/arxiv.0803.0559

Related papers

BEC-BCS crossover of a trapped Fermi gas without using the local density approximation

Rocío Jáuregui, R. Paredes, G. Toledo Sánchez

We perform a variational quantum Monte Carlo simulation of an interacting Fermi gas confined in a three dimensional harmonic potential. This gas is considered as the precursor system from which a molecular bosonic gas is formed. Based on the results of two-body calculations for trapped atoms, we propose a family of variational many-body wave functions that takes into account the qualitative different nature of the Bardeen-Cooper-Schieffer and Bose-Einstein condensate regimes as a function of the scattering length. Energies, densities, and correlation functions are calculated and compared with previous results for homogeneous gases. Universality tests at the unitarity limit are performed including the verification of the virial relation and the evaluation of the $\ensuremath{\beta}$ parameter.

arXiv:cond-mat/0608689•Physical Review A•Jul 10, 2007•15 citations

Virial relations for ultracold trapped Fermi gases with finite range interactions through the Bardeen–Cooper–Schrieffer Bose–Einstein-condensate crossover

L. Rosales-Zárate, Rocío Jáuregui

We study the virial relations for ultracold trapped two-component Fermi gases in the case of short finite range interactions. Numerical verifications for such relations are reported through the Bardeen–Cooper–Schrieffer (BCS) Bose–Einstein-condensate (BEC) crossover. As an intermediate step, it is necessary to evaluate the partial derivatives of the many-body energy with respect to the inverse of the scattering length and with respect to the interaction range. Once the binding energy of the formed molecules in the BEC side is subtracted, the corresponding energy derivatives are found to have extreme values at the unitary limit. The value of the derivative with respect to the potential range in that limit is large enough to yield measurable differences between the total energy and twice the trapping energy unless the interacting system is described by extremely short potential ranges. The virial results are used to check the quality of the variational wavefunction involved in the calculations.

New Journal of Physics•Apr 20, 2010•0 citations

Nonlinear Schrödinger equation for a superfluid Fermi gas in the BCS-BEC crossover

Sadhan K. Adhikari

We introduce a quasianalytic nonlinear Schr\"odinger equation with beyond mean-field corrections to describe the dynamics of a zero-temperature dilute superfluid Fermi gas in the crossover from the weak-coupling Bardeen-Cooper-Schrieffer (BCS) regime, where ${k}_{F}|a|⪡1$ with $a$ the $s$-wave scattering length and ${k}_{F}$ the Fermi momentum, through the unitarity limit ${k}_{F}a\ensuremath{\rightarrow}\ifmmode\pm\else\textpm\fi{}\ensuremath{\infty}$ to the Bose-Einstein condensation (BEC) regime where ${k}_{F}a>0$. The energy of our model is parametrized using the known asymptotic behavior in the BCS, BEC, and the unitarity limits and is in excellent agreement with accurate Green's-function Monte Carlo calculations. The model generates good results for frequencies of collective breathing oscillations of a trapped Fermi superfluid.

arXiv:0802.2222•Physical Review A•Apr 9, 2008•75 citations

Quantum Monte Carlo studies of superfluid Fermi gases

Siyuan Chang, Vijay R. Pandharipande, Joseph Carlson +1 more

We report results of quantum Monte Carlo calculations of the ground state of dilute Fermi gases with attractive short-range two-body interactions. The strength of the interaction is varied to study different pairing regimes which are characterized by the product of the $s$-wave scattering length and the Fermi wave vector, $a{k}_{F}$. We report results for the ground-state energy, the pairing gap $\ensuremath{\Delta}$, and the quasiparticle spectrum. In the weak-coupling regime, $1∕a{k}_{F}<\ensuremath{-}1$, we obtain Bardeen-Cooper-Schrieffer (BCS) superfluid and the energy gap $\ensuremath{\Delta}$ is much smaller than the Fermi gas energy ${E}_{\mathrm{FG}}$. When $a>0$, the interaction is strong enough to form bound molecules with energy ${E}_{\mathrm{mol}}$. For $1∕a{k}_{F}\ensuremath{\gtrsim}0.5$, we find that weakly interacting composite bosons are formed in the superfluid gas with $\ensuremath{\Delta}$ and gas energy per particle approaching $\ensuremath{\mid}{E}_{\mathrm{mol}}\ensuremath{\mid}∕2$. In this region, we seem to have Bose-Einstein condensation (BEC) of molecules. The behavior of the energy and the gap in the BCS-to-BEC transition region, $\ensuremath{-}0.5<1∕a{k}_{F}<0.5$, is discussed.

arXiv:physics/0404115•Physical Review A•Oct 5, 2004•221 citations

Collective modes of low-dimensional superfluid Fermi gases in the BCS-BEC crossover: Time-dependent variational analysis

Yu Long Zhou, Guoxiang Huang

We investigate collective modes and free expansions of quasi-one- and quasi-two-dimensional (Q1D and Q2D) ultracold Fermi gases in the crossover from a Bardeen-Cooper-Schrieffer (BCS) superfluid to a Bose-Einstein condensate (BEC). We solve a superfluid order parameter equation valid for the BCS-BEC crossover by employing a time-dependent variational method. We take a trial wave function of hybrid Gaussian-parabolic type, which not only reflects the low-dimensional character of the system but also allows an essentially analytical approach for the problem. We present Q1D and Q2D criteria that are valid in various superfluid regimes and show clearly the relation between the maximum condensed particle number and the parameters of the trapping potential as well as the atom-atom interaction. We demonstrate that, due to the small particle number in Q1D and Q2D condensates, the contribution to oscillating frequencies of collective modes by the quantum pressure in the strong-confinement direction is significant and hence the Thomas-Fermi approximation cannot be used. We also show that the free expansion of Q1D and Q2D superfluid Fermi gases in the strong-confinement direction is much faster than that in the weak-confinement direction.

Physical Review A•Feb 12, 2007•22 citations

QUANTUM MONTE CARLO STUDY OF THE GROUND-STATE PROPERTIES OF A FERMI GAS IN THE BCS-BEC CROSSOVER

S. Giorgini, G. E. Astrakharchik, Jordi Boronat +1 more

The ground-state properties of a two-component Fermi gas with attractive short-range interactions are calculated using the fixed-node diffusion Monte Carlo method. The interaction strength is varied over a wide range by tuning the value a of the s-wave scattering length of the two-body potential. We calculate the energy per particle, the one- and two-body density matrix as a function of the interaction strength. Results for the momentum distribution of the atoms, as obtained from the Fourier transform of the one-body density matrix, are reported as a function of the interaction strength. Off-diagonal long-range order in the system is investigated through the asymptotic behavior of the two-body density matrix. The condensate fraction of pairs is calculated in the unitary limit and on both sides of the BCS-BEC crossover.

International Journal of Modern Physics B•Dec 15, 2006•0 citations
We perform a variational quantum Monte Carlo simulation of the transition from a Bardeen-Cooper-Schrieffer superfluid (BCS) to a Bose-Einstein condensate (BEC) at zero temperature. The model Hamiltonian involves an attractive short range two body interaction and the atoms number 2N=3302N =3302N=330 is chosen so that, in the non-interacting limit, the ground state function corresponds to a closed shell configuration. The system is then characterized by the s-wave scattering length aaa of the two-particle collisions in the gas, which is varied from negative to positive values, and the Fermi wave number kFk_FkF​. Based on an extensive analysis of the s-wave two-body problem, one parameter variational many-body wave functions are proposed to describe the ground state of the interacting Fermi gas from BCS to BEC states. We exploit properties of antisymmetrized many-body functions to develop efficient techniques that permit variational calculations for a large number of particles. It is shown that a virial relation between the energy per particle and the trapping energy is approximately valid for −0.1<1/kFa<3.4-0.1<1/k_Fa<3.4−0.1<1/kF​a<3.4. The influence of the harmonic trap and the interaction potential as exhibited in two-body correlation functions is also analyzed.
We perform a variational quantum Monte Carlo simulation of an interacting Fermi gas confined in a three dimensional harmonic potential. This gas is considered as the precursor system from which a molecular bosonic gas is formed. Based on the results of two-body calculations for trapped atoms, we propose a family of variational many-body wave functions that takes into account the qualitative different nature of the Bardeen-Cooper-Schieffer and Bose-Einstein condensate regimes as a function of the scattering length. Energies, densities, and correlation functions are calculated and compared with previous results for homogeneous gases. Universality tests at the unitarity limit are performed including the verification of the virial relation and the evaluation of the $\ensuremath{\beta}$ parameter.
We study the virial relations for ultracold trapped two-component Fermi gases in the case of short finite range interactions. Numerical verifications for such relations are reported through the Bardeen–Cooper–Schrieffer (BCS) Bose–Einstein-condensate (BEC) crossover. As an intermediate step, it is necessary to evaluate the partial derivatives of the many-body energy with respect to the inverse of the scattering length and with respect to the interaction range. Once the binding energy of the formed molecules in the BEC side is subtracted, the corresponding energy derivatives are found to have extreme values at the unitary limit. The value of the derivative with respect to the potential range in that limit is large enough to yield measurable differences between the total energy and twice the trapping energy unless the interacting system is described by extremely short potential ranges. The virial results are used to check the quality of the variational wavefunction involved in the calculations.
We introduce a quasianalytic nonlinear Schr\"odinger equation with beyond mean-field corrections to describe the dynamics of a zero-temperature dilute superfluid Fermi gas in the crossover from the weak-coupling Bardeen-Cooper-Schrieffer (BCS) regime, where kF∣a∣⪡1{k}_{F}|a|⪡1kF​∣a∣⪡1 with aaa the sss-wave scattering length and kF{k}_{F}kF​ the Fermi momentum, through the unitarity limit ${k}_{F}a\ensuremath{\rightarrow}\ifmmode\pm\else\textpm\fi{}\ensuremath{\infty}$ to the Bose-Einstein condensation (BEC) regime where kFa>0{k}_{F}a>0kF​a>0. The energy of our model is parametrized using the known asymptotic behavior in the BCS, BEC, and the unitarity limits and is in excellent agreement with accurate Green's-function Monte Carlo calculations. The model generates good results for frequencies of collective breathing oscillations of a trapped Fermi superfluid.
We report results of quantum Monte Carlo calculations of the ground state of dilute Fermi gases with attractive short-range two-body interactions. The strength of the interaction is varied to study different pairing regimes which are characterized by the product of the sss-wave scattering length and the Fermi wave vector, akFa{k}_{F}akF​. We report results for the ground-state energy, the pairing gap $\ensuremath{\Delta}$, and the quasiparticle spectrum. In the weak-coupling regime, $1∕a{k}_{F}<\ensuremath{-}1$, we obtain Bardeen-Cooper-Schrieffer (BCS) superfluid and the energy gap $\ensuremath{\Delta}$ is much smaller than the Fermi gas energy EFG{E}_{\mathrm{FG}}EFG​. When a>0a>0a>0, the interaction is strong enough to form bound molecules with energy Emol{E}_{\mathrm{mol}}Emol​. For $1∕a{k}_{F}\ensuremath{\gtrsim}0.5$, we find that weakly interacting composite bosons are formed in the superfluid gas with $\ensuremath{\Delta}$ and gas energy per particle approaching $\ensuremath{\mid}{E}_{\mathrm{mol}}\ensuremath{\mid}∕2$. In this region, we seem to have Bose-Einstein condensation (BEC) of molecules. The behavior of the energy and the gap in the BCS-to-BEC transition region, $\ensuremath{-}0.5<1∕a{k}_{F}<0.5$, is discussed.
We investigate collective modes and free expansions of quasi-one- and quasi-two-dimensional (Q1D and Q2D) ultracold Fermi gases in the crossover from a Bardeen-Cooper-Schrieffer (BCS) superfluid to a Bose-Einstein condensate (BEC). We solve a superfluid order parameter equation valid for the BCS-BEC crossover by employing a time-dependent variational method. We take a trial wave function of hybrid Gaussian-parabolic type, which not only reflects the low-dimensional character of the system but also allows an essentially analytical approach for the problem. We present Q1D and Q2D criteria that are valid in various superfluid regimes and show clearly the relation between the maximum condensed particle number and the parameters of the trapping potential as well as the atom-atom interaction. We demonstrate that, due to the small particle number in Q1D and Q2D condensates, the contribution to oscillating frequencies of collective modes by the quantum pressure in the strong-confinement direction is significant and hence the Thomas-Fermi approximation cannot be used. We also show that the free expansion of Q1D and Q2D superfluid Fermi gases in the strong-confinement direction is much faster than that in the weak-confinement direction.
The ground-state properties of a two-component Fermi gas with attractive short-range interactions are calculated using the fixed-node diffusion Monte Carlo method. The interaction strength is varied over a wide range by tuning the value a of the s-wave scattering length of the two-body potential. We calculate the energy per particle, the one- and two-body density matrix as a function of the interaction strength. Results for the momentum distribution of the atoms, as obtained from the Fourier transform of the one-body density matrix, are reported as a function of the interaction strength. Off-diagonal long-range order in the system is investigated through the asymptotic behavior of the two-body density matrix. The condensate fraction of pairs is calculated in the unitary limit and on both sides of the BCS-BEC crossover.