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:cond-mat/0101395 · Physical Review A

Bose-Einstein condensates in the large-gas-parameter regime

Adelchi Fabrocini, A. Polls

Published Nov 16, 2001• 38 citations
arXivHTMLPDF

Abstract

Bose-Einstein condensates of ${10}^{4}{}^{85}\mathrm{Rb}$ atoms in a cylindrical trap are studied using a recently proposed approach based on a local-density approximation. Since the existence of a Feshbach resonance allows for widely tuning the scattering length of the atoms, values of the peak gas parameter ${x}_{\mathrm{pk}}$ of the order of ${10}^{\ensuremath{-}2}$ may be attained and the standard analysis based on Gross-Pitaevskii and/or Thomas-Fermi equations may result in being questionable. Energy functionals derived from the correlated basis functions theory and the low-density expansion for a homogeneous hard-spheres gas are used to estimate corrections to the Gross-Pitaevskii equation. The two functionals give similar results, both showing large differences with respect to the Thomas-Fermi and Gross-Pitaevskii ones. The column densities at $z=0$ may differ by as much as $\ensuremath{\sim}30%$ and the half maximum radius by $\ensuremath{\sim}20%.$ The scattering lengths estimated by fitting the half maximum radius may differ by $\ensuremath{\sim}40%$ from those given by a Thomas-Fermi-based analysis of the experimental data.

DOI
10.1103/physreva.64.063610

Related papers

Ground state properties and elementary excitation of85Rb Bose-Einstein condensate: analytical method for modified Thomas-Fermi model with large-gas-parameter

Moumita Gupta, Krishna Rai Dastidar

Bose-Einstein condensation in trapped gases of 85 Rb atoms can now be obtained in the large gas-parameter region in presence of Feshbach resonance and the gas density which is characterized by peak gas parameter x pk can attain the value ~ 10 -2 . For such large gas-parameter, the first correction term to the Bogoliubov energy functional is considered to give modified Gross-Pitaevskii (MGP) equation. In this work analytical expressions for the chemical potential, the density of particles and the excitation frequency of the monopole mode have been derived in Thomas-Fermi(TF) approximated approach to the MGP equation (modified TF approach) following method of higher order iteration for Bose-Einstein condensate (BEC) of atoms in spherically symmetric trap. We have also solved the MGP equation by using Runge-Kutta method for integration of nonlinear differential equation considering spherically symmetric trap. The analytical results for chemical potential and density of particles for BEC of 10 3 and 10 4 85 Rb atoms in a spherically symmetric trap are in very good agreement with the numerical results from MGP equation at low as well as at large scattering length i.e, at large gas-parameter region. The frequency of monopole mode calculated using the density of particles in the ground state (found out in higher order iterative method) are also in very good agreement with the existing numerical results. The difference between the analytical and numerical results decreases as the value of scattering length i.e, gas-parameter is increased.

Journal of Physics Conference Series•Sep 1, 2007•4 citations

Beyond-mean-field results for atomic Bose-Einstein condensates at interaction strengths near Feshbach resonances: A many-body dimensional perturbation-theory calculation

Brett Allen McKinney, Martin Dunn, D. K. Watson

We present semianalytical many-body results for energies and excitation frequencies for an inhomogeneous Bose-Einstein condensate over a wide range of atom numbers $N$ for both small $s$-wave scattering lengths, typical of most laboratory experiments, and large scattering lengths, achieved by tuning through a Feshbach resonance. Our dimensional perturbation treatment includes two-body correlations at all orders and yields analytical results through first order by taking advantage of the high degree of symmetry of the condensate at the zeroth-order limit. Because $N$ remains a parameter in our analytical results, the challenge of calculating energies and excitation frequencies does not rise with the number of condensate atoms. In this proof-of-concept paper the atoms are confined in a spherical trap and are treated as hard spheres. Our many-body calculations compare well to Gross-Pitaevskii results in the weakly interacting regime and depart from the mean-field approximation as the density approaches the strongly interacting regime. The excitation frequencies provide a particularly sensitive test of beyond-mean-field corrections. For example, for $N=2000$ atoms and an experimentally realized large scattering length of $a=0.433{a}_{ho}\phantom{\rule{0.3em}{0ex}}({a}_{ho}=\sqrt{\ensuremath{\hbar}∕m{\ensuremath{\omega}}_{ho}})$ we predict a $75%$ shift from the mean-field breathing mode frequency.

Physical Review A•May 13, 2004•19 citations

Beyond the Gross-Pitaevskii approximation: Local density versus correlated basis approach for trapped bosons

Adelchi Fabrocini, A. Polls

We study the ground state of a system of Bose hard spheres trapped in an isotropic harmonic potential to investigate the effect of the interatomic correlations and the accuracy of the Gross-Pitaevskii equation. We compare a local-density approximation, based on the energy functional derived from the low-density expansion of the energy of the uniform hard-sphere gas, and a correlated wave-function approach, which explicitly introduces the correlations induced by the potential. Both higher-order terms in the low-density expansion, beyond Gross-Pitaevskii, and explicit dynamical correlations have effects of the order of percent when the number of trapped particles becomes similar to that attained in recent experiments $(N\ensuremath{\sim}{10}^{7}).$

arXiv:cond-mat/9901342•Physical Review A•Sep 1, 1999•75 citations

Thomas-Fermi approximation for Bose-Einstein condensates in traps

Peter Schuck, X. Viñas

Thomas-Fermi theory for Bose-Einstein condensates in inhomogeneous traps is revisited. The phase-space distribution function of the condensate in the Thomas-Fermi limit $(\stackrel{\ensuremath{\rightarrow}}{\ensuremath{\Elzxh}}0)$ is ${f}_{0}(\mathit{R},\mathit{p})\ensuremath{\propto}\ensuremath{\delta}(\ensuremath{\mu}\ensuremath{-}{H}_{\mathrm{cl}})$ where ${H}_{\mathrm{cl}}$ is the classical counterpart of the self-consistent Gross-Pitaevskii Hamiltonian. Starting from this distribution function the Thomas-Fermi kinetic energy is calculated for any number of particles. Good agreement between the Gross-Pitaevskii and Thomas-Fermi kinetic energies is found even for low and intermediate particle numbers N. Application of this Thomas-Fermi theory to the attractive case and to the calculation of the frequencies of the monopole and quadrupole excitations in the sum rule approach yields conclusive results as well. The difference with the usual Thomas-Fermi approach to the Bose-Einstein condensates (large-N limit) is discussed in detail.

Physical Review A•Mar 8, 2000•32 citations

Single-particle excitations in a trapped gas of Fermi atoms in the BCS-BEC crossover region. II. Broad Feshbach resonance

Yoji Ohashi, Adam Griffin

We apply the formulation developed in a recent paper [Y. Ohashi and A. Griffin, Phys. Rev. A 72, 013601 (2005)] for single-particle excitations in the BCS-BEC crossover to the case of a broad Feshbach resonance. At $T=0$, we solve the Bogoliubov--de Gennes coupled equations taking into account a Bose condensate of bound states (molecules). In the case of a broad resonance, the density profile $n(r)$, as well as the profile of the superfluid order parameter $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\Delta}}(r)$, are spatially spread out to the Thomas-Fermi radius, even in the crossover region. This order parameter $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\Delta}}(r)$ suppresses the effects of low-energy Andreev bound states on the rf tunneling current. As a result, the peak energy in the rf spectrum is found to occur at an energy equal to the superfluid order parameter $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\Delta}}(r=0)$ at the center of the trap, in contrast to the case of a narrow resonance, and in agreement with recent measurements. The local density approximation is found to give a good approximation for the rf-tunneling spectrum.

arXiv:cond-mat/0508213•Physical Review A•Dec 5, 2005•25 citations

Scaling equations for Bose-Einstein condensate dynamics across all interaction regimes

D. C. Marinica, C. Puertas González, T. Estrampes +5 more

We derive a unified set of scaling equations for Bose-Einstein condensates in time-dependent harmonic traps, connecting the weakly interacting Gaussian regime to the strongly interacting Thomas-Fermi regime. The exact Gross-Pitaevskii ground-state density is taken as a fixed profile in rescaled coordinates, with its dynamics described by three scaling factors corresponding to compression or expansion along the three spatial directions. The resulting equations contain no adjustable parameter, since their coefficients are determined once from the initial state. They recover analytically the Gaussian variational and Thomas-Fermi scaling equations in their respective limits and satisfy an axis-resolved virial theorem. The approach also yields the global phase and the low-lying collective-mode frequencies of the condensate. We benchmark the model against three-dimensional Gross-Pitaevskii simulations and against measured expansion energies. The model remains accurate across all interaction regimes, in strongly anisotropic traps, and along time-dependent sequences including a relatively fast quench, up to the point where the underlying frozen-profile hypothesis, shared by all three scaling approaches, breaks down.

cond-mat.quant-gasCross-listedarXiv:2609.13497•Sep 15, 2026•0 citations
Bose-Einstein condensates of 10485Rb{10}^{4}{}^{85}\mathrm{Rb}10485Rb atoms in a cylindrical trap are studied using a recently proposed approach based on a local-density approximation. Since the existence of a Feshbach resonance allows for widely tuning the scattering length of the atoms, values of the peak gas parameter xpk{x}_{\mathrm{pk}}xpk​ of the order of 10\ensuremath−2{10}^{\ensuremath{-}2}10\ensuremath−2 may be attained and the standard analysis based on Gross-Pitaevskii and/or Thomas-Fermi equations may result in being questionable. Energy functionals derived from the correlated basis functions theory and the low-density expansion for a homogeneous hard-spheres gas are used to estimate corrections to the Gross-Pitaevskii equation. The two functionals give similar results, both showing large differences with respect to the Thomas-Fermi and Gross-Pitaevskii ones. The column densities at z=0z=0z=0 may differ by as much as \ensuremath∼30\ensuremath{\sim}30%\ensuremath∼30 and the half maximum radius by \ensuremath∼20\ensuremath{\sim}20%.\ensuremath∼20 The scattering lengths estimated by fitting the half maximum radius may differ by \ensuremath∼40\ensuremath{\sim}40%\ensuremath∼40 from those given by a Thomas-Fermi-based analysis of the experimental data.
Bose-Einstein condensation in trapped gases of 85 Rb atoms can now be obtained in the large gas-parameter region in presence of Feshbach resonance and the gas density which is characterized by peak gas parameter x pk can attain the value ~ 10 -2 . For such large gas-parameter, the first correction term to the Bogoliubov energy functional is considered to give modified Gross-Pitaevskii (MGP) equation. In this work analytical expressions for the chemical potential, the density of particles and the excitation frequency of the monopole mode have been derived in Thomas-Fermi(TF) approximated approach to the MGP equation (modified TF approach) following method of higher order iteration for Bose-Einstein condensate (BEC) of atoms in spherically symmetric trap. We have also solved the MGP equation by using Runge-Kutta method for integration of nonlinear differential equation considering spherically symmetric trap. The analytical results for chemical potential and density of particles for BEC of 10 3 and 10 4 85 Rb atoms in a spherically symmetric trap are in very good agreement with the numerical results from MGP equation at low as well as at large scattering length i.e, at large gas-parameter region. The frequency of monopole mode calculated using the density of particles in the ground state (found out in higher order iterative method) are also in very good agreement with the existing numerical results. The difference between the analytical and numerical results decreases as the value of scattering length i.e, gas-parameter is increased.
We present semianalytical many-body results for energies and excitation frequencies for an inhomogeneous Bose-Einstein condensate over a wide range of atom numbers NNN for both small sss-wave scattering lengths, typical of most laboratory experiments, and large scattering lengths, achieved by tuning through a Feshbach resonance. Our dimensional perturbation treatment includes two-body correlations at all orders and yields analytical results through first order by taking advantage of the high degree of symmetry of the condensate at the zeroth-order limit. Because NNN remains a parameter in our analytical results, the challenge of calculating energies and excitation frequencies does not rise with the number of condensate atoms. In this proof-of-concept paper the atoms are confined in a spherical trap and are treated as hard spheres. Our many-body calculations compare well to Gross-Pitaevskii results in the weakly interacting regime and depart from the mean-field approximation as the density approaches the strongly interacting regime. The excitation frequencies provide a particularly sensitive test of beyond-mean-field corrections. For example, for N=2000N=2000N=2000 atoms and an experimentally realized large scattering length of $a=0.433{a}_{ho}\phantom{\rule{0.3em}{0ex}}({a}_{ho}=\sqrt{\ensuremath{\hbar}∕m{\ensuremath{\omega}}_{ho}})$ we predict a 7575%75 shift from the mean-field breathing mode frequency.
We study the ground state of a system of Bose hard spheres trapped in an isotropic harmonic potential to investigate the effect of the interatomic correlations and the accuracy of the Gross-Pitaevskii equation. We compare a local-density approximation, based on the energy functional derived from the low-density expansion of the energy of the uniform hard-sphere gas, and a correlated wave-function approach, which explicitly introduces the correlations induced by the potential. Both higher-order terms in the low-density expansion, beyond Gross-Pitaevskii, and explicit dynamical correlations have effects of the order of percent when the number of trapped particles becomes similar to that attained in recent experiments $(N\ensuremath{\sim}{10}^{7}).$
Thomas-Fermi theory for Bose-Einstein condensates in inhomogeneous traps is revisited. The phase-space distribution function of the condensate in the Thomas-Fermi limit $(\stackrel{\ensuremath{\rightarrow}}{\ensuremath{\Elzxh}}0)$ is ${f}_{0}(\mathit{R},\mathit{p})\ensuremath{\propto}\ensuremath{\delta}(\ensuremath{\mu}\ensuremath{-}{H}_{\mathrm{cl}})$ where Hcl{H}_{\mathrm{cl}}Hcl​ is the classical counterpart of the self-consistent Gross-Pitaevskii Hamiltonian. Starting from this distribution function the Thomas-Fermi kinetic energy is calculated for any number of particles. Good agreement between the Gross-Pitaevskii and Thomas-Fermi kinetic energies is found even for low and intermediate particle numbers N. Application of this Thomas-Fermi theory to the attractive case and to the calculation of the frequencies of the monopole and quadrupole excitations in the sum rule approach yields conclusive results as well. The difference with the usual Thomas-Fermi approach to the Bose-Einstein condensates (large-N limit) is discussed in detail.
We apply the formulation developed in a recent paper [Y. Ohashi and A. Griffin, Phys. Rev. A 72, 013601 (2005)] for single-particle excitations in the BCS-BEC crossover to the case of a broad Feshbach resonance. At T=0T=0T=0, we solve the Bogoliubov--de Gennes coupled equations taking into account a Bose condensate of bound states (molecules). In the case of a broad resonance, the density profile n(r)n(r)n(r), as well as the profile of the superfluid order parameter $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\Delta}}(r)$, are spatially spread out to the Thomas-Fermi radius, even in the crossover region. This order parameter $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\Delta}}(r)$ suppresses the effects of low-energy Andreev bound states on the rf tunneling current. As a result, the peak energy in the rf spectrum is found to occur at an energy equal to the superfluid order parameter $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\Delta}}(r=0)$ at the center of the trap, in contrast to the case of a narrow resonance, and in agreement with recent measurements. The local density approximation is found to give a good approximation for the rf-tunneling spectrum.
We derive a unified set of scaling equations for Bose-Einstein condensates in time-dependent harmonic traps, connecting the weakly interacting Gaussian regime to the strongly interacting Thomas-Fermi regime. The exact Gross-Pitaevskii ground-state density is taken as a fixed profile in rescaled coordinates, with its dynamics described by three scaling factors corresponding to compression or expansion along the three spatial directions. The resulting equations contain no adjustable parameter, since their coefficients are determined once from the initial state. They recover analytically the Gaussian variational and Thomas-Fermi scaling equations in their respective limits and satisfy an axis-resolved virial theorem. The approach also yields the global phase and the low-lying collective-mode frequencies of the condensate. We benchmark the model against three-dimensional Gross-Pitaevskii simulations and against measured expansion energies. The model remains accurate across all interaction regimes, in strongly anisotropic traps, and along time-dependent sequences including a relatively fast quench, up to the point where the underlying frozen-profile hypothesis, shared by all three scaling approaches, breaks down.