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:1706.07096 · Physical Review A

Metastability versus collapse following a quench in attractive Bose-Einstein condensates

Jake Golde, Joanna Ruhl, Maxim Olshanii, Vanja Dunjko, Sumita Datta, Boris A. Malomed

Published May 9, 2018• 6 citations
arXivHTMLPDF

Abstract

We consider a Bose-Einstein condensate (BEC) with attractive two-body interactions in a cigar-shaped trap, initially prepared in its ground state for a given negative scattering length, which is quenched to a larger absolute value of the scattering length. Using the mean-field approximation, we compute numerically, for an experimentally relevant range of aspect ratios and initial strengths of the coupling, two critical values of quench. One corresponds to the weakest attraction strength, the quench to which causes the system to collapse before completing even a single return from the narrow configuration (pericenter) in its breathing cycle. The other is a similar critical point for the occurrence of collapse before completing two returns. In the latter case, we also compute the limiting value, as we keep increasing the strength of the postquench attraction towards its critical value, of the time interval between the first two pericenters. We also use a Gaussian variational model to estimate the critical quenched attraction strength below which the system is stable against the collapse for long times. These time intervals and critical attraction strengths, apart from being fundamental properties of nonlinear dynamics of self-attractive BECs, may provide clues to the design of upcoming experiments that are trying to create robust BEC breathers.

DOI
10.1103/physreva.97.053604

Related papers

Instability of a Bose-Einstein condensate with an attractive interaction

Antonios Eleftheriou, Kerson Huang

We study the stability of a Bose-Einstein condensate of harmonically trapped atoms with negative scattering length, specifically ${}^{7}\mathrm{Li}.$ Our method is to solve the time-dependent nonlinear Schr\"odinger equation numerically. For an isolated condensate, with no gain or loss, we find that the system is stable (apart from quantum tunneling) if the particle number N is less than a critical number ${N}_{c}.$ For $N>{N}_{c},$ the system collapses to high-density clumps in a region near the center of the trap. The time for the onset of collapse is on the order of one trap period. Within numerical uncertainty, the results are consistent with the formation of a ``black hole'' of infinite density fluctuations, as predicted by Ueda and Huang [Phys. Rev. A (to be published)]. We numerically obtain ${N}_{c}\ensuremath{\approx}1251.$ We then include gain-loss mechanisms, i.e., the gain of atoms from a surrounding ``thermal cloud,'' and the loss due to two- and three-body collisions. The number N now oscillates in a steady state, with a period of about 145 trap periods. We obtain ${N}_{c}\ensuremath{\approx}1260$ as the maximum value in the oscillations.

arXiv:cond-mat/9908229•Physical Review A•Mar 2, 2000•55 citations

Dynamics and stability of Bose-Einstein condensates with attractive 1/r interaction

Holger Cartarius, Tomaž Fabčič, Jörg Main +1 more

The time-dependent extended Gross-Pitaevskii equation for Bose-Einstein condensates with attractive $1/r$ interaction is investigated with both a variational approach and numerically exact calculations. We show that these condensates exhibit signatures known from the nonlinear dynamics of autonomous Hamiltonian systems. The two stationary solutions created in a tangent bifurcation at a critical value of the scattering length are identified as elliptical and hyperbolical fixed points, corresponding to stable and unstable stationary states of the condensate. The stable stationary state is surrounded by elliptical islands, corresponding to condensates periodically oscillating in time, whereas condensate wave functions in the unstable region undergo a collapse within finite time. For negative scattering lengths below the tangent bifurcation no stationary solutions exist, i.e., the condensate is always unstable and collapses.

arXiv:0803.1799•Physical Review A•Jul 11, 2008•23 citations

Collapse of a self-gravitating Bose-Einstein condensate with attractive self-interaction

Pierre-Henri Chavanis

We study the collapse of a self-gravitating Bose-Einstein condensate with attractive self-interaction. Equilibrium states in which the gravitational attraction and the attraction due to the self-interaction are counterbalanced by the quantum pressure (Heisenberg's uncertainty principle) exist only below a maximum mass ${M}_{\mathrm{max}}=1.012\ensuremath{\hbar}/\sqrt{Gm|{a}_{s}|}$ where ${a}_{s}<0$ is the scattering length of the bosons and $m$ is their mass [P. H. Chavanis, Phys. Rev. D 84, 043531 (2011)]. For $M>{M}_{\mathrm{max}}$ the system is expected to collapse and form a black hole. We study the collapse dynamics by making a Gaussian ansatz for the wave function and reducing the problem to the study of the motion of a particle in an effective potential. We find that the collapse time scales as $(M/{M}_{\mathrm{max}}\ensuremath{-}1{)}^{\ensuremath{-}1/4}$ for $M\ensuremath{\rightarrow}{M}_{\mathrm{max}}^{+}$ and as ${M}^{\ensuremath{-}1/2}$ for $M\ensuremath{\gg}{M}_{\mathrm{max}}$. Other analytical results are given above and below the critical point corresponding to a saddle-node bifurcation. We apply our results to QCD axions with mass $m={10}^{\ensuremath{-}4}\text{ }\text{ }\mathrm{eV}/{c}^{2}$ and scattering length ${a}_{s}=\ensuremath{-}5.8\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}53}\text{ }\text{ }\mathrm{m}$ for which ${M}_{\mathrm{max}}=6.5\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}14}{M}_{\ensuremath{\bigodot}}$ and $R=3.3\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}4}{R}_{\ensuremath{\bigodot}}$. We confirm our previous claim that bosons with attractive self-interaction, such as QCD axions, may form low mass stars (axion stars or dark matter stars) but cannot form dark matter halos of relevant mass and size. These mini axion stars could be the constituents of dark matter. They can collapse into mini black holes of mass $\ensuremath{\sim}{10}^{\ensuremath{-}14}{M}_{\ensuremath{\bigodot}}$ in a few hours. In that case, dark matter halos would be made of mini black holes. We also apply our results to ultralight axions with mass $m=1.93\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}20}\text{ }\text{ }\mathrm{eV}/{c}^{2}$ and scattering length ${a}_{s}=\ensuremath{-}8.29\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}60}\text{ }\text{ }\mathrm{fm}$ for which ${M}_{\mathrm{max}}=0.39\ifmmode\times\else\texttimes\fi{}1{0}^{6}{M}_{\ensuremath{\bigodot}}$ and $R=33\text{ }\text{ }\mathrm{pc}$. These ultralight axions could cluster into dark matter halos. Axionic dark matter halos with attractive self-interaction can collapse into supermassive black holes of mass $\ensuremath{\sim}1{0}^{6}{M}_{\ensuremath{\bigodot}}$ (similar to those reported at the center of galaxies) in about one million years. We point out the limitations of the Gaussian ansatz to describe the late stages of the collapse dynamics. We also mention the possibility that, instead of forming a black hole, the collapse may be accompanied by a burst or relativistic axions (bosenova) leading to a cycle of collapses and explosions as observed for nongravitational Bose-Einstein condensates with attractive self-interaction.

arXiv:1604.05904•Physical review. D/Physical review. D.•Oct 20, 2016•112 citations

Stability of trapped Bose-Einstein condensates

F. Kh. Abdullaev, A. Gammal, Lauro Tomio +1 more

In three-dimensional trapped Bose-Einstein condensate (BEC), described by the time-dependent Gross-Pitaevskii-Ginzburg equation, we study the effect of initial conditions on stability using a Gaussian variational approach and exact numerical simulations. We also discuss the validity of the criterion for stability suggested by Vakhitov and Kolokolov. The maximum initial chirp (initial focusing defocusing of cloud) that can lead a stable condensate to collapse even before the number of atoms reaches its critical limit is obtained for several specific cases. When we consider two- and three-body nonlinear terms, with negative cubic and positive quintic terms, we have the conditions for the existence of two phases in the condensate. In this case, the magnitude of the oscillations between the two phases are studied considering sufficient large initial chirps. The occurrence of collapse in a BEC with repulsive two-body interaction is also shown to be possible.

arXiv:cond-mat/0102018•Physical Review A•Mar 12, 2001•201 citations

Coupled Bose-Einstein condensate: Collapse for attractive interaction

Sadhan K. Adhikari

We study the collapse in a coupled Bose-Einstein condensate of two types of bosons 1 and 2 under the action of a trap using the time-dependent Gross-Pitaevskii equation. The system may undergo collapse when one, two, or three of the scattering lengths ${a}_{\mathrm{ij}}$ for scattering of boson i with j, $i,j=1,2,$ are negative representing an attractive interaction. Depending on the parameters of the problem a single or both components of the condensate may experience collapse.

arXiv:cond-mat/0011526•Physical Review A•Mar 20, 2001•50 citations

Controlled Collapse of a Bose-Einstein Condensate

Jacob Roberts, N. R. Claussen, Simon L. Cornish +3 more

The point of instability of a Bose-Einstein condensate (BEC) due to attractive interactions was studied. Stable 85Rb BECs were created and then caused to collapse by slowly changing the atom-atom interaction from repulsive to attractive using a Feshbach resonance. At a critical value, an abrupt transition was observed in which atoms were ejected from the condensate. By measuring the onset of this transition as a function of number and attractive interaction strength, we determined the stability condition to be N(absolute value of a) / a(ho) = 0.459+/-0.012+/-0.054, slightly lower than the predicted value of 0.574.

arXiv:cond-mat/0102116•Physical Review Letters•May 7, 2001•439 citations
We consider a Bose-Einstein condensate (BEC) with attractive two-body interactions in a cigar-shaped trap, initially prepared in its ground state for a given negative scattering length, which is quenched to a larger absolute value of the scattering length. Using the mean-field approximation, we compute numerically, for an experimentally relevant range of aspect ratios and initial strengths of the coupling, two critical values of quench. One corresponds to the weakest attraction strength, the quench to which causes the system to collapse before completing even a single return from the narrow configuration (pericenter) in its breathing cycle. The other is a similar critical point for the occurrence of collapse before completing two returns. In the latter case, we also compute the limiting value, as we keep increasing the strength of the postquench attraction towards its critical value, of the time interval between the first two pericenters. We also use a Gaussian variational model to estimate the critical quenched attraction strength below which the system is stable against the collapse for long times. These time intervals and critical attraction strengths, apart from being fundamental properties of nonlinear dynamics of self-attractive BECs, may provide clues to the design of upcoming experiments that are trying to create robust BEC breathers.
We study the stability of a Bose-Einstein condensate of harmonically trapped atoms with negative scattering length, specifically 7Li.{}^{7}\mathrm{Li}.7Li. Our method is to solve the time-dependent nonlinear Schr\"odinger equation numerically. For an isolated condensate, with no gain or loss, we find that the system is stable (apart from quantum tunneling) if the particle number N is less than a critical number Nc.{N}_{c}.Nc​. For N>Nc,N>{N}_{c},N>Nc​, the system collapses to high-density clumps in a region near the center of the trap. The time for the onset of collapse is on the order of one trap period. Within numerical uncertainty, the results are consistent with the formation of a ``black hole'' of infinite density fluctuations, as predicted by Ueda and Huang [Phys. Rev. A (to be published)]. We numerically obtain ${N}_{c}\ensuremath{\approx}1251.$ We then include gain-loss mechanisms, i.e., the gain of atoms from a surrounding ``thermal cloud,'' and the loss due to two- and three-body collisions. The number N now oscillates in a steady state, with a period of about 145 trap periods. We obtain ${N}_{c}\ensuremath{\approx}1260$ as the maximum value in the oscillations.
The time-dependent extended Gross-Pitaevskii equation for Bose-Einstein condensates with attractive 1/r1/r1/r interaction is investigated with both a variational approach and numerically exact calculations. We show that these condensates exhibit signatures known from the nonlinear dynamics of autonomous Hamiltonian systems. The two stationary solutions created in a tangent bifurcation at a critical value of the scattering length are identified as elliptical and hyperbolical fixed points, corresponding to stable and unstable stationary states of the condensate. The stable stationary state is surrounded by elliptical islands, corresponding to condensates periodically oscillating in time, whereas condensate wave functions in the unstable region undergo a collapse within finite time. For negative scattering lengths below the tangent bifurcation no stationary solutions exist, i.e., the condensate is always unstable and collapses.
We study the collapse of a self-gravitating Bose-Einstein condensate with attractive self-interaction. Equilibrium states in which the gravitational attraction and the attraction due to the self-interaction are counterbalanced by the quantum pressure (Heisenberg's uncertainty principle) exist only below a maximum mass ${M}_{\mathrm{max}}=1.012\ensuremath{\hbar}/\sqrt{Gm|{a}_{s}|}$ where as<0{a}_{s}<0as​<0 is the scattering length of the bosons and mmm is their mass [P. H. Chavanis, Phys. Rev. D 84, 043531 (2011)]. For M>MmaxM>{M}_{\mathrm{max}}M>Mmax​ the system is expected to collapse and form a black hole. We study the collapse dynamics by making a Gaussian ansatz for the wave function and reducing the problem to the study of the motion of a particle in an effective potential. We find that the collapse time scales as $(M/{M}_{\mathrm{max}}\ensuremath{-}1{)}^{\ensuremath{-}1/4}$ for $M\ensuremath{\rightarrow}{M}_{\mathrm{max}}^{+}$ and as ${M}^{\ensuremath{-}1/2}$ for $M\ensuremath{\gg}{M}_{\mathrm{max}}$. Other analytical results are given above and below the critical point corresponding to a saddle-node bifurcation. We apply our results to QCD axions with mass $m={10}^{\ensuremath{-}4}\text{ }\text{ }\mathrm{eV}/{c}^{2}$ and scattering length ${a}_{s}=\ensuremath{-}5.8\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}53}\text{ }\text{ }\mathrm{m}$ for which ${M}_{\mathrm{max}}=6.5\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}14}{M}_{\ensuremath{\bigodot}}$ and $R=3.3\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}4}{R}_{\ensuremath{\bigodot}}$. We confirm our previous claim that bosons with attractive self-interaction, such as QCD axions, may form low mass stars (axion stars or dark matter stars) but cannot form dark matter halos of relevant mass and size. These mini axion stars could be the constituents of dark matter. They can collapse into mini black holes of mass $\ensuremath{\sim}{10}^{\ensuremath{-}14}{M}_{\ensuremath{\bigodot}}$ in a few hours. In that case, dark matter halos would be made of mini black holes. We also apply our results to ultralight axions with mass $m=1.93\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}20}\text{ }\text{ }\mathrm{eV}/{c}^{2}$ and scattering length ${a}_{s}=\ensuremath{-}8.29\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}60}\text{ }\text{ }\mathrm{fm}$ for which ${M}_{\mathrm{max}}=0.39\ifmmode\times\else\texttimes\fi{}1{0}^{6}{M}_{\ensuremath{\bigodot}}$ and R=33  pcR=33\text{ }\text{ }\mathrm{pc}R=33  pc. These ultralight axions could cluster into dark matter halos. Axionic dark matter halos with attractive self-interaction can collapse into supermassive black holes of mass $\ensuremath{\sim}1{0}^{6}{M}_{\ensuremath{\bigodot}}$ (similar to those reported at the center of galaxies) in about one million years. We point out the limitations of the Gaussian ansatz to describe the late stages of the collapse dynamics. We also mention the possibility that, instead of forming a black hole, the collapse may be accompanied by a burst or relativistic axions (bosenova) leading to a cycle of collapses and explosions as observed for nongravitational Bose-Einstein condensates with attractive self-interaction.
In three-dimensional trapped Bose-Einstein condensate (BEC), described by the time-dependent Gross-Pitaevskii-Ginzburg equation, we study the effect of initial conditions on stability using a Gaussian variational approach and exact numerical simulations. We also discuss the validity of the criterion for stability suggested by Vakhitov and Kolokolov. The maximum initial chirp (initial focusing defocusing of cloud) that can lead a stable condensate to collapse even before the number of atoms reaches its critical limit is obtained for several specific cases. When we consider two- and three-body nonlinear terms, with negative cubic and positive quintic terms, we have the conditions for the existence of two phases in the condensate. In this case, the magnitude of the oscillations between the two phases are studied considering sufficient large initial chirps. The occurrence of collapse in a BEC with repulsive two-body interaction is also shown to be possible.
We study the collapse in a coupled Bose-Einstein condensate of two types of bosons 1 and 2 under the action of a trap using the time-dependent Gross-Pitaevskii equation. The system may undergo collapse when one, two, or three of the scattering lengths aij{a}_{\mathrm{ij}}aij​ for scattering of boson i with j, i,j=1,2,i,j=1,2,i,j=1,2, are negative representing an attractive interaction. Depending on the parameters of the problem a single or both components of the condensate may experience collapse.
The point of instability of a Bose-Einstein condensate (BEC) due to attractive interactions was studied. Stable 85Rb BECs were created and then caused to collapse by slowly changing the atom-atom interaction from repulsive to attractive using a Feshbach resonance. At a critical value, an abrupt transition was observed in which atoms were ejected from the condensate. By measuring the onset of this transition as a function of number and attractive interaction strength, we determined the stability condition to be N(absolute value of a) / a(ho) = 0.459+/-0.012+/-0.054, slightly lower than the predicted value of 0.574.