Pith. sign in

REVIEW 4 major objections 4 minor 38 references

Collapsing dynamics of attractive Bose-Einstein condensates in random potentials

T0 review · 4 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper claims that a random potential can delay the collapse of an attractive Bose-Einstein condensate, raising the critical interaction strength by about 13 percent at the parameters studied.

desk verdict A careful variational calculation undermined by one unspecified disorder realization and an overclaimed collapse-prevention conclusion. read the letter →

arxiv 2411.16825 v1 pith:CR2FRJVD submitted 2024-11-25 cond-mat.quant-gas cond-mat.dis-nn

Madjda Kamel , Abdelaali Boudjemaa This is my paper · ORCID
classification cond-mat.quant-gascond-mat.dis-nn MSC 35Q5582B44
keywords Bose-EinsteincondensatecollapsedisorderrandompotentialGross-Pitaevskiiequationvariationalmethodbreathingmodespeckle
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper asks whether a random potential can change the collapse threshold of a trapped Bose-Einstein condensate with attractive interactions. Solving the Gross-Pitaevskii equation with a variational Gaussian ansatz, the authors report that a repulsive speckle disorder raises the critical interaction strength from about 0.67 to about 0.8 at the parameters considered, a 13 percent increase in the critical atom number. They also derive the breathing-mode frequency and show numerically that moderate disorder creates density modulations while very strong disorder fragments the condensate. The intended upshot is that disorder acts as a new mechanism for delaying or preventing the collapse of attractive condensates.

What carries the argument

The variational Gaussian ansatz for the condensate wavefunction, with width $q$ as the static variational parameter, is inserted into the Gross-Pitaevskii energy functional. The resulting effective potential $U_{\rm eff}(q)$ contains the trap, kinetic, attraction, and disorder terms; collapse occurs when the local minimum of $U_{\rm eff}$ disappears, which defines $\eta_{\rm cr}$. The breathing-mode frequency is the curvature of $U_{\rm eff}$ at the equilibrium width. The disorder enters as a sum over $S$ repulsive Gaussian impurities of width $\sigma$, so the entire argument rides on that single-realization sum.

What would settle it

Compute the critical interaction strength for many independent random realizations of the speckle potential at fixed parameters $\tilde U_0=0.1$ and $\sigma=0.2$; if the spread of $\eta_{\rm cr}$ across realizations is comparable to or larger than the 13 percent shift, the stabilization is not a property of disorder itself.

Watch

Extended reading notes

Core claim

The central claim is that the interplay of attractive interactions and a random potential shifts the stability boundary of a three-dimensional harmonically trapped Bose-Einstein condensate: the critical interaction parameter $\eta_{\rm cr}=Na/l$ grows from about $0.67$ in the clean case to about $0.8$ for disorder strength $\tilde U_0=0.1$ and correlation length $\sigma=0.2$, meaning roughly 13 percent more atoms can be held before collapse. The stabilization is traced to the disorder contribution to the energy landscape: for condensate width $q<1$ the random potential deepens the local energy maximum and shallows the minimum, so the metastable state survives to larger attraction. The same energy analysis yields the breathing-mode frequency, which increases with disorder strength and drops to zero at the collapse threshold.

Load-bearing premise

The paper treats one fixed arrangement of impurities as representative of disorder, so the reported stabilization may depend on where the impurities sit rather than being a general property of random potentials.

Editorial extensions

If this is right

  • The critical atom number for collapse is not fixed by the trap and scattering length alone; a repulsive random potential can raise it, so experiments with speckle disorder should see a delayed collapse.
  • The breathing-mode frequency increases with disorder strength at fixed interaction strength, giving a measurable dynamical signature of the stabilizing effect.
  • Moderate disorder induces density modulations in an otherwise smooth condensate, implying that the disorder imprints spatial structure even when the condensate remains phase coherent.
  • Sufficiently strong disorder ($\tilde U_0 \gtrsim 2$) fragments the condensate into several peaks, so the stabilizing regime is limited to moderate disorder.
  • Near the collapse threshold the density modulations are significantly reduced, indicating that the collapse dynamics themselves are altered by the disorder potential.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The reported 13 percent shift is computed for one fixed realization of impurity positions; a proper disorder average over many realizations is likely needed before the stabilization can be called a property of random potentials rather than of a particular configuration.
  • Because the effect arises from repulsive impurities sitting near the cloud center, an experimental test could deposit controlled single impurities at different radial positions and measure how the collapse threshold moves.
  • The same variational machinery could be extended to repulsive interactions, where disorder and interactions compete differently, to check whether the stabilizing role of disorder is specific to the attractive side.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. The manuscript studies a three-dimensional harmonically trapped Bose-Einstein condensate with attractive interactions in a speckle-type random potential. Using a Gaussian variational ansatz, the authors derive closed-form expressions for the energy (Eq. (6)), the effective potential for the width dynamics (Eq. (17)), and the breathing-mode frequency (Eq. (19)). They also solve the stationary and time-dependent Gross-Pitaevskii equation numerically with a split-step Fourier method. The central quantitative result is that the critical interaction parameter increases from approximately 0.671 in the clean case to approximately 0.8 for disorder strength U0=0.1 and correlation length sigma=0.2, presented as evidence that disorder prevents collapse. The paper additionally reports disorder-induced density modulations and a decrease of the breathing frequency as the interaction strength approaches its critical value.

Significance. If the central claim were established, the result would be notable: a repulsive random potential would stabilize an attractive BEC beyond the clean-trap limit, with analytic predictions for the critical atom number, width, and breathing frequency that could be compared with future speckle experiments. The variational calculation is internally consistent and transparent, and all predictions follow from substitution of the Gaussian ansatz into the Gross-Pitaevskii equation without fitting to target outputs. However, the quantitative stabilization result is not established as a property of disorder: it is computed for a single unspecified realization of impurity positions, and because the energy is linear in the zero-mean disorder potential, the ensemble-averaged variational energy reduces to the clean result. The numerical density at eta=1 in Fig. 6 shows the signature of collapse rather than its prevention. These issues make the central claim, as stated in the abstract, unsupported.

major comments (4)
  1. [III A, Eq. (6), Fig. 2] The reported increase of eta_cr from about 0.671 to about 0.8 is obtained from Eq. (6) for one fixed set of impurity positions r_i, but the manuscript never specifies those positions and never performs an ensemble average over disorder realizations. Since Section II defines the disorder with mean zero and the energy functional (4) is linear in U_dis, the ensemble average of the variational energy is the clean-BEC energy up to a q-independent constant; no q-dependent disorder contribution survives averaging at first order. The 13% shift is therefore a configuration-specific fluctuation rather than a statistical property of a random potential. The authors should compute the disorder-averaged eta_cr with an appropriate measure of variance, or demonstrate that the chosen realization is typical; without this, the abstract's claim that disorder prevents collapse is unsupported.
  2. [IV, Fig. 6] Figure 6 and the surrounding text describe the density at eta=1, which is beyond the clean critical value eta_cr approximately 0.671. The text states that the condensate width decreases while its amplitude increases with disorder, which is the standard signature of collapse, not of collapse prevention. The following sentence, stating that the number of atoms is continuously decreasing, contradicts the conserved norm enforced by Eq. (8) and Eq. (11), where dN/dt=0. This internal inconsistency must be resolved; as written, the numerical evidence points in the opposite direction from the abstract's conclusion.
  3. [III A, Eq. (5)] The Gaussian ansatz (5) has no center-of-mass displacement parameter. In a disordered potential the condensate can lower its energy by shifting away from repulsive impurities, so a fixed-centered Gaussian is forced to overlap any impurities located near the origin and will overestimate the energetic cost of that particular disorder configuration, and hence the apparent stabilization. A variational treatment of a disordered potential should include a centroid parameter or justify explicitly why it can be omitted.
  4. [IV, numerical method] The numerical speckle potential is generated by mapping random numbers into the interval [0,L], but the manuscript does not report the specific realization, the random seed, or any check of realization-to-realization variability for the quantities shown in Figs. 5 and 6. The numerical results are therefore not reproducible, and they cannot support a general claim about disorder without a demonstration that the chosen realization is representative.
minor comments (4)
  1. [Title] The title contains an extraneous space in the word 'condensa tes'.
  2. [Fig. 2(b)] The horizontal-axis label is rendered as 'U/OverTilde 0'; it should be typeset as \tilde{U}_0.
  3. [Eqs. (6) and (10)] Equations (6) and (10) are identical in content; the duplication is unnecessary and could be confusing, since one appears in the equilibrium context and the other in the Lagrangian context.
  4. [References] Reference [28] appears in the bibliography but does not seem to be cited in the main text.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the variational derivation is self-contained and contains no fitted-input predictions.

full rationale

The paper's central results are obtained by inserting the Gaussian ansatz (5) into the GPE energy functional (4), producing the closed-form energy (6); eta_cr is then found by locating the saddle point where the local minimum disappears. Neither the disorder parameters (U0, sigma, r_i) nor the clean-limit eta_cr=0.671 are fitted to the reported output eta_cr ~ 0.8; the clean value is reproduced at U0=0 as a consistency check. The breathing frequency (19) follows from a standard linearization of the equation of motion (16), with no target data used in the derivation. Self-citations ([19]-[22], [37]) appear only as contextual references for disordered BECs or for the standard perturbation expansion q=q0+delta-q, and none of these citations supplies the load-bearing premise or forbids alternatives. The concern that a single realization of the speckle pattern may not be statistically representative, and the internal inconsistency in the claim that the atom number decreases despite norm conservation (Eqs. (8), (11)), are correctness/statistical-validity issues, not circularity; they do not make the derivation equivalent to its inputs by construction.

Assumptions & free parameters 4 free parameters · 4 assumptions · 0 invented entities

The central quantitative claim rests on several hand-picked model parameters (sigma, U0, S, and one realization of impurity positions) and on two ad hoc modeling assumptions: the Gaussian ansatz and the Gaussian-impurity speckle model. None of these are derived, and the paper provides no sensitivity analysis or ensemble averaging, so the ledger is heavy for the strength of the conclusion.

free parameters (4)
  • sigma (disorder correlation length) = 0.2
    Chosen by hand; the key quantitative claim (eta_cr = 0.8) depends on this value, and the paper shows eta_cr varies strongly with sigma.
  • U0 (disorder strength) = 0.1 (also 0.2 and 2 in plots)
    Chosen by hand; eta_cr increases linearly with U0, so the reported 13% increase is tied to this specific choice.
  • S (number of Gaussian impurities) = 300
    Chosen for the numerical speckle generation; no convergence check with S is provided, and the variational calculation does not state S.
  • impurity positions r_i = not specified, random in [0,L] with L=30
    The variational energy in Eq. (6) depends on the specific impurity positions, but no seed or distribution is given, and a single realization is used without ensemble averaging.
assumptions (4)
  • domain assumption The Gross-Pitaevskii equation is an adequate mean-field description of the 3D attractive BEC in a random potential, including near collapse.
    Invoked in Eq. (1); no beyond-mean-field corrections are considered, but quantum fluctuations become important near collapse, so this is a nontrivial assumption.
  • ad hoc to paper A single-Gaussian ansatz with width q suffices to describe the ground state and collapse threshold.
    Assumed in Eq. (5); the numerics show fragmentation into multiple peaks for strong disorder, so the ansatz breaks down in regimes the paper discusses, making the variational eta_cr unreliable at large U0.
  • ad hoc to paper The speckle potential is faithfully modeled by a sum of positive Gaussian impurities.
    Eq. (2) models speckle as U0 times a sum of Gaussians; real speckle has Rayleigh amplitude statistics and continuous spatial correlations, so this is a simplification that is not validated in the paper.
  • ad hoc to paper One fixed disorder realization (specific r_i) represents the ensemble of random potentials.
    No averaging over disorder realizations is performed for the variational energy, yet the conclusion 'disorder prevents collapse' is stated for the statistical concept of disorder.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Collapsing dynamics of attractive Bose-Einstein condensates in random potentials." pith.science (2026). https://pith.science/paper/CR2FRJVD

@misc{pith2026241116825,
  author       = {Pith},
  title        = {Pith review of: Collapsing dynamics of attractive Bose-Einstein condensates in random potentials},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/CR2FRJVD}},
  note         = {Machine review of arXiv:2411.16825}
}
read the original abstract

We study the stationary and dynamical properties of three-dimensional trapped Bose-Einstein condensates with attractive interactions subjected to a random potential. To this end, a variational method is applied to solve the underlying Gross-Pitaevskii equation. We derive analytical predictions for the energy, the equilibrium width, and evolution laws of the condensate parameter. The breathing mode oscillations frequency of the condensate has been also calculated in terms of the gas and disorder parameters. We analyze in addition the dynamics of collapse from the Gaussian approximation. Surprisingly, we find that the intriguing interplay of the attractive interaction and disorder effects leads to prevent collapse of the condensate.

Figures

Figures reproduced from arXiv: 2411.16825 by the authors.

Figure 1
Figure 1. FIG. 1: Energy from equation (6) as a function of the width [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: (a) Critical interaction strength [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: Time evolution of the condensate width from Eq. (16) fo [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: FIG. 4: Breathing mode oscillation frequency as a function o [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5: (a) Density profiles of a disordered attractive BEC fo [PITH_FULL_IMAGE:figures/full_fig_p007_5.png]
Figure 6
Figure 6. Figure 6: ). The condensate width decreases while its amplitude increases with the disorder indicating that the number of atoms is continuously decreasing. V. CONCLUSION In this paper we systematically studied the equilibrium and dynamical properties of attractive BECs confined …

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

38 extracted references · 38 canonical work pages

  1. [1]

    The dependence of ω on η is shown in Fig

    7/ 2 S∑ i=1 [ U1 + U2 ] , (19) where U1 =  e − r2 i q2 0 +σ 2 ( 2q2 0σr 2 i q2 0 + σ 2 + σ 3 )  2 − erfc   ri σ √ σ 2 q2 0 + 1     + 2q0σ 2rie− r2 i σ 2 √π √ q2 0 + σ 2   ( 12q2 0 − 3σ 2) − 12q0rie − r2 i q2 0 +σ 2 ( 2q2 0 σr 2 i q2 0 +σ 2 + σ 3 )     q0ri  2 − erfc   ri σ √ σ 2 q2 0 +1     + σ √ q2 0 +σ 2e − r2 i σ 2 ( σ 2 q2 0 +...

    work page

  2. [2]

    M.H.Anderson, J. R. Ensher, M.R. Matthews, C.E.Wieman, E.A. Cornell, Science 269, 198(1995)

    work page 1995

  3. [3]

    Andrews, N.J.van Druten, D

    K.B.Davis, M.O.Mewes, M.R. Andrews, N.J.van Druten, D. S.Durfee, D.M.Kurn, W.Ketterle, Phys.Rev.Lett. 75, 3969(1995)

    work page 1995

  4. [4]

    Bradley, C

    C.C. Bradley, C. A. Sackett, J. J. Tollett, and R.G. Hulet , Phys. Rev. Lett. 75, 1687 (1995)

    work page 1995

  5. [5]

    R. J. Dodd, M. Edwards, C. J. Williams, C. W. Clark, M. J. Ho lland, P. A. Ruprecht and K. Burnett, Phys. Rev. A 54, 661 (1996)

    work page 1996

  6. [6]

    J. M. Gerton, D. Strekalov, I. Prodan, and R. G. Hulet, Nat ure 408, 692 (2000)

    work page 2000

  7. [7]

    Aspect and M

    A. Aspect and M. Inguscio, Phys. Today. 62, 30 (2009)

    work page 2009

  8. [8]

    Modugno, Rep

    G. Modugno, Rep. Prog. Phys. 73 , 102401 (2010). 9

    work page 2010

Show all 38 references
  1. [9]

    C. A. Sackett, H. T. C. Stoof, and R. G. Hulet, Phys. Rev. Le tt. 80, 2031 (1998)

    1998

  2. [10]

    Mardonov, E

    Sh. Mardonov, E. Ya. Sherman, J. G. Muga, H-W. Wang, Y. Ban , X. Chen, Phys. Rev. A 91, 043604 (2015)

    2015

  3. [11]

    Jing Li, B. A. Malomed, Wenliang Xi Chen, E. Ya Sherman, C ommunications in Nonlinear Science and Numerical Simu- lation, 82, 105045 (2020)

    2020

  4. [12]

    M.P.A.Fisher,P.B.Weichman, G.Grinstein, and D.S.Fi sher,Phys.Rev.B 40, 546 (1989)

    1989

  5. [13]

    Crooker, B.Hebral, E.N

    B.C. Crooker, B.Hebral, E.N. Smith,Y.Takano, and J.D. Reppy,Phys.Rev.Lett.51,666 (1983)

    1983

  6. [14]

    Low Temp.Phys

    J.D.Reppy, J. Low Temp.Phys. 87, 205 (1992)

    1992

  7. [15]

    Bogoliubov, J

    N.N. Bogoliubov, J. Phys. USSR 11, 2332 (1947)

    1947

  8. [16]

    Huang and H.-F

    K. Huang and H.-F. Meng, Phys. Rev. Lett. 69, 644 (1992)

    1992

  9. [17]

    Giorgini, L

    S. Giorgini, L. Pitaevskii, and S. Stringari, Phys. Rev . B 49, 12938 (1994)

    1994

  10. [18]

    G. M. Falco, J. Phys. B: At. Mol. Opt. Phys. 42, 215303 (2009)

    2009

  11. [19]

    Lima, Aristeu RP and Pelster, Axel, Phys.Rev.A 84, 041604 (2011)

    2011

  12. [20]

    Boudjemˆ aa, Phys

    A. Boudjemˆ aa, Phys. Rev. A 91, 053619 (2015)

    2015

  13. [21]

    Boudjemˆ aa, Eur

    A. Boudjemˆ aa, Eur. Phys. J. B 92, 145 (2019)

    2019

  14. [22]

    Boudjemˆ aa and K

    A. Boudjemˆ aa and K. Abbas, Phys. Rev. A 102, 023325 (2020)

    2020

  15. [23]

    Yadjour and A

    A. Yadjour and A. Boudjemˆ aa, Physica B. 654, 414683 (2023)

    2023

  16. [24]

    Pasienski, D

    M. Pasienski, D. McKay, M. White, and B. DeMarco, Nature Phys. 6, 677 (2010)

    2010

  17. [25]

    Deissler, M

    B. Deissler, M. Zaccanti, G. Roati, C. D’Errico, M. Fatt ori,M. Modugno, G. Modugno and M. Inguscio, Nature Phys. 6, 354 (2010)

    2010

  18. [26]

    Nagler, M

    B. Nagler, M. Radonji´ c, S. Barbosa, J. Koch, A. Pelster , A. Widera, New. J. Phys. 22, 033021 (2020)

    2020

  19. [27]

    Akkermans et al

    E. Akkermans et al. , J. Phys. B: At. Mol. Opt. Phys. 41, 045302 (2008)

    2008

  20. [28]

    Yong P. Chen, J. Hitchcock, D. Dries, M. Junker, C. Welfo rd, S.E. Pollack, T.A. Corcovilos, R.G. Hulet, Physica D 238, 1321 (2009)

    2009

  21. [29]

    B. Min, T. Li, M. Rosenkranz, and W. Bao, Phys. Rev. A 86, 053612 (2012)

    2012

  22. [30]

    Stellin, M

    F. Stellin, M. Filoche, and F. Dias, Phys. Rev. A 107, 043306 (2023)

    2023

  23. [31]

    Sanchez-Palencia, D

    L. Sanchez-Palencia, D. Cl´ ement, P. Lugan, P. Bouyer a nd A. Aspect, New J. Phys. 10, 045019 (2008)

    2008

  24. [32]

    S. K. Adhikari, Phys. Lett. A 265 , 91 (2000)

    2000

  25. [33]

    W. Bao, D. Jaksch, P. A. Markowich, J. Comput. Phys. 187, 318 (2003)

    2003

  26. [34]

    N Arnold, G

    D. N Arnold, G. David, M. Filoche, D. Jerison, S. Mayboro da, Communications in Partial Differential Equations, 44 1186 (2019)

    2019

  27. [35]

    Altmann, P

    R. Altmann, P. Henning, D. Peterseim, Numerische Mathe matik, 148, 575 (2021)

    2021

  28. [36]

    S. K. Turitsyn, I. Gabitov, E. W. Laedke, V. K. Mezent-se v, S. L. Musher, E. G. Shapiro, T. Schafer, and K. H. Spatschek , Opt. Commun. 151, 117 (1998)

    1998

  29. [37]

    Cheng and S

    Y. Cheng and S. K. Adhikari, Phys. Rev. A 82, 013631 (2010)

    2010

  30. [38]

    Abbas and A

    K. Abbas and A. Boudjemaa , Phys. Rev. A, 107, 033306 (2023)

    2023

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.