{"id":"ba1495ad-8036-4a4b-8cd4-e0a93d3e574c","arxiv_id":"2411.16825","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":4,"one_line_summary":"Disorder in the form of a repulsive random potential raises the predicted collapse threshold of a 3D attractive BEC, but only for one realization, and it does not actually prevent collapse.","lead":"A variational calculation of three-dimensional attractive Bose-Einstein condensates in a repulsive random potential predicts that the critical atom number for collapse increases. But the result comes from a single disorder realization, and the paper's own numerics show shrinking condensates at higher interaction strengths, so the headline claim that disorder prevents collapse is not supported.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reported 13% increase in ηcr rests on one unspecified disorder realization; since ⟨Udis⟩=0, the linear disorder term in Eq. (6) averages to zero, so the stabilization is a fluctuation, not a statistical property of random potentials.","rationale":"The reader's weakest_assumption identified the absence of ensemble averaging over disorder realizations as the load-bearing flaw, and my independent read reaches the same conclusion. The paper's own definition of the random potential has ⟨Udis⟩=0, and the variational energy depends linearly on Udis, so the averaged energy is that of the clean system. The reported stabilization must therefore come from a single atypical configuration, not from the statistical ensemble. I also note a closely related variational limitation: the Gaussian ansatz is pinned at the trap center and cannot shift to avoid repulsive impurities, which systematically inflates the disorder energy barrier and hence ηcr. This does not change the verdict: the central claim as stated is unsupported and the paper needs either ensemble averaging, a softened claim, or a demonstration of typicality. The proposed concrete test directly quantifies how much of the 13% effect survives averaging and center relaxation; if the effect vanishes, the rejection is confirmed. This is a technical, good-faith check of the manuscript's weakest point and does not rely on any external consensus or ad hominem assessment.","tokens_in":9210,"tokens_out":5317,"duration_ms":59250,"concrete_test":"Repeat the variational calculation of Sec. III A for at least 100 independent random realizations with stated parameters (U0=0.1, σ=0.2, S=300), drawing impurity positions uniformly in a box of size L=30 and computing ηcr for each by locating the saddle-node of E(q)/N in Eq. (6). Then repeat for a typical realization with an extended Gaussian ansatz φ∝exp[-(r-r_c)^2/(2q^2)] allowing center displacement r_c. If the mean ηcr minus 0.671 is within one standard deviation, the reported 13% increase is a realization artifact.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quantitative central claim (ηcr≈0.8 vs clean 0.671 in Fig. 2) is obtained by minimizing Eq. (6) for a single fixed set of impurity positions {r_i}. The paper never specifies those positions, never performs an ensemble average, and reports no variance. This is structurally decisive: Sec. II defines the disorder with ⟨Udis(r)⟩=0, and the disorder contribution to the energy functional (4) is linear in Udis. Therefore the ensemble average of the variational energy is the clean-BEC energy (up to a q-independent constant), so no q-dependent disorder term survives averaging at first order. The 13% shift can only arise from a particular configuration with one or more repulsive impurities near the origin, where the fixed-centered Gaussian ansatz is forced to overlap them. The ansatz (5) contains no center-of-mass displacement parameter, so it cannot lower the energy by shifting away from impurities, which further overestimates any stabilization. In addition, the numerical density at η=1 (Fig. 6) shows decreasing width and increasing amplitude, i.e. collapse-like behavior, and the statement that 'the number of atoms is continuously decreasing' contradicts the conserved norm used in Eqs. (8) and (11). Without a disorder average, or at least a demonstration that the chosen realization is typical, the abstract's claim that disorder prevents collapse is unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":9497,"tokens_out":5608,"duration_ms":52091,"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":[{"comment":"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.","section":"III A, Eq. (6), Fig. 2"},{"comment":"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.","section":"IV, Fig. 6"},{"comment":"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.","section":"III A, Eq. (5)"},{"comment":"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.","section":"IV, numerical method"}],"minor_comments":[{"comment":"The title contains an extraneous space in the word 'condensa tes'.","section":"Title"},{"comment":"The horizontal-axis label is rendered as 'U/OverTilde 0'; it should be typeset as \\tilde{U}_0.","section":"Fig. 2(b)"},{"comment":"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.","section":"Eqs. (6) and (10)"},{"comment":"Reference [28] appears in the bibliography but does not seem to be cited in the main text.","section":"References"}],"recommendation":"reject","confidential_remarks":"The paper's main novelty is the claim that disorder prevents collapse of an attractive BEC. As written, that claim rests on a single unspecified disorder realization, while the ensemble-averaged variational energy contains no disorder contribution at first order. The numerical results at eta=1 also show collapse-like behavior. I do not see a path to the stated conclusion within the current manuscript; a substantially different study, including a proper disorder average and a treatment of the condensate centroid, would be needed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper has a legitimate new calculation, but the central claim that disorder raises the critical interaction strength and prevents collapse does not survive contact with its own model. The reported 13% increase in eta_cr comes from a single fixed set of impurity positions that are never specified, never averaged, and the disorder term in the energy is linear in a zero-mean field, so ensemble averaging would wash it out. The headline effect is likely an artifact of one particular configuration.\n\nWhat is actually new: a Gaussian variational treatment of a harmonically trapped 3D attractive BEC in a speckle-type potential, with closed-form expressions for the energy, effective potential, and breathing-mode frequency. The derivation is long but straightforward, and the clean-limit eta_cr ~ 0.671 checks out. The paper also includes split-step numerical density profiles, which is a reasonable sanity check. I credit them for the honest clean-limit comparison.\n\nThe soft spots are real and load-bearing. No ensemble average is performed anywhere. The model defines Udis with <Udis>=0, yet Eq. (6) is minimized for a single realization. Since the disorder contribution to the energy is linear in Udis, its q-dependence vanishes after averaging over impurity positions; the shift in eta_cr can only come from a configuration with an impurity near the origin. The Gaussian ansatz has no center-of-mass shift, so the trial function is forced to sit on the impurity, further exaggerating any stabilizing effect. The r_i positions are never reported, so the calculation is not reproducible. Second, the paper's own Fig. 6 shows the width shrinking and amplitude growing at eta=1, the classic collapse signature, and the text says this 'indicates that the number of atoms is continuously decreasing,' which directly contradicts the conserved norm in Eqs. (8) and (11). That is not a minor slip; it suggests a misreading of their own numerics. The abstract's 'prevent collapse' is too strong given that the effect, even if real, is a 13% shift for one configuration, and above the new threshold collapse still occurs.\n\nWho this is for: someone working on variational methods for disordered BECs might extract the useful algebra, but the physical conclusion should not be cited. The paper would need proper disorder averaging, a softened claim, and a fix of the atom-number statement before it is referee-ready. In its current form I would not engage it.\n\nRecommendation: reject for peer review; the main result is not supported as stated.","headline":"A careful variational calculation undermined by one unspecified disorder realization and an overclaimed collapse-prevention conclusion.","tokens_in":10023,"tokens_out":3448,"would_cite":false,"duration_ms":34135,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q55","82B44"],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["Bose-Einstein condensate","collapse","disorder","random potential","Gross-Pitaevskii equation","variational method","breathing mode","speckle potential"],"falsifier":"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.","tokens_in":8961,"feed_emoji":"⚛️","tokens_out":2735,"duration_ms":24906,"temperature":0.7,"pith_summary":"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.","feed_headline":"Random potential raises BEC collapse threshold by 13%","feed_subtitle":"A speckle disorder shifts the critical atom number for an attractive condensate from about 0.67 to 0.8.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Provides the experimental observation of collapse of an attractive BEC that the paper uses as the clean-system baseline.","marker":"[3]"},{"why":"Supplies the direct observation of collapse dynamics and the $^7$Li parameters used to estimate the critical atom number.","marker":"[5]"},{"why":"Supplies the speckle potential model used to generate the disorder potential in the paper.","marker":"[30]"},{"why":"Provides the split-step Fourier spectral method used for numerical integration of the Gross-Pitaevskii equation.","marker":"[36]"},{"why":"Supplies the linearization procedure for the condensate width used to derive the breathing-mode frequencies.","marker":"[37]"},{"why":"Provides prior experimental and numerical evidence on attractive BECs in disorder, giving context for the disorder-dependent width.","marker":"[27]"}],"fun_headline_variants":["Random potential delays BEC collapse by 13%","Disorder gives attractive BEC more room before collapse","Speckle disorder raises BEC collapse threshold","Attractive BEC: disorder prevents collapse","Disorder stabilizes attractive BEC against collapse"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Random potential delays BEC collapse by 13%","Disorder gives attractive BEC more room before collapse","Speckle disorder raises BEC collapse threshold","Attractive BEC: disorder prevents collapse","Disorder stabilizes attractive BEC against collapse"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000553,"raw_usage":{"total_tokens":2569,"prompt_tokens":812,"completion_tokens":1757,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":428,"completion_tokens_details":{"reasoning_tokens":1699}},"tokens_in":428,"tokens_out":1757,"duration_ms":12348,"temperature":1.0,"reasoning_tokens":1699,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:52:09.138888+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Andrews, N.J.van Druten, D","cited_arxiv_id":null,"evidence_quote":"Provides the experimental observation of collapse of an attractive BEC that the paper uses as the clean-system baseline."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the direct observation of collapse dynamics and the $^7$Li parameters used to estimate the critical atom number."},{"cited_title":"Stellin, M","cited_arxiv_id":null,"evidence_quote":"Supplies the speckle potential model used to generate the disorder potential in the paper."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the split-step Fourier spectral method used for numerical integration of the Gross-Pitaevskii equation."},{"cited_title":"Cheng and S","cited_arxiv_id":null,"evidence_quote":"Supplies the linearization procedure for the condensate width used to derive the breathing-mode frequencies."},{"cited_title":"Akkermans et al","cited_arxiv_id":null,"evidence_quote":"Provides prior experimental and numerical evidence on attractive BECs in disorder, giving context for the disorder-dependent width."}],"review_version":1}