Pith. sign in

REVIEW 3 major objections 3 minor 29 references

Interaction between gravitational waves and trapped Bose-Einstein condensates

T0 review · 3 major / 3 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 A gravitational wave shifts the quantum phase of a trapped Bose-Einstein condensate by Nξ(t), and interacting atoms make the shift four orders of magnitude larger.

desk verdict Core phase-shift formulas have a clear algebraic error, so the headline numbers are not reliable; the qualitative BEC-GW idea survives but needs a careful re-derivation. read the letter →

arxiv 2411.15874 v1 pith:ILNDQEHO submitted 2024-11-24 cond-mat.quant-gas gr-qc

A. Perodi , L. Salasnich This is my paper · ORCID
classification cond-mat.quant-gasgr-qc
keywords gravitationalwavesBose-EinsteincondensatesquantumfidelityphaseshiftanisotropicharmonictrapGross-Pitaevskiiequationcoherentstatesmany-bodyenhancement
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 establishes that a passing gravitational wave imprints a many-body phase on a trapped Bose-Einstein condensate, with the phase growing in proportion to the atom number N. Working in the transverse-traceless gauge and using a curved-space Gross-Pitaevskii equation, the authors compute the fidelity amplitude of the condensate's coherent state and find F(t)=1-iNξ(t)+O(N²h²). For an anisotropic harmonic trap, they derive explicit formulas for ξ(t) for both non-interacting and interacting atoms, and show that atomic interactions boost the phase by roughly four orders of magnitude, making interacting condensates the more promising detector medium.

What carries the argument

The key object is the curved-space Gross-Pitaevskii equation, Eq. (11), in which the gravitational wave enters only through the kinetic term h(∂²_x-∂²_y)ψ. The fidelity amplitude F(t)=⟨Ψ_cs|U_int(t)|Ψ_cs⟩ is evaluated at first order in Nh by expanding the time-evolution operator, and the phase ξ(t) is computed as the single-particle expectation of (∂²_y-∂²_x) over the flat-space ground state φ_gs. For interacting atoms, the ground state is a variational Gaussian with the widths determined by minimizing the GPE energy functional, following Ref. [23].

What would settle it

Compute the non-relativistic limit of the Klein-Gordon action in the TT gauge keeping all metric perturbations to first order in h, including h_00 and derivative-of-h terms, and compare the resulting single-particle phase with Eq. (23); if additional terms contribute at the same order, the predicted Nξ(t) changes. Experimentally, a tabletop interferometer with an anisotropic trapped condensate could measure the phase as a function of scattering length and check whether the enhancement follows Eq. (28).

Watch

Extended reading notes

Core claim

The central claim is that the fidelity amplitude of a BEC coherent state, evolved under a gravitational wave perturbation, acquires a phase Nξ(t) that is first-order in the wave amplitude h and linearly enhanced by the particle number N. For a condensate in an anisotropic harmonic trap, the non-interacting phase is Nξ(t)=N h0/2 √(ℏ/m)(√ω_x-√ω_y) $e^{{-ℏω_g²/(4mc²ω_z)}}$ sin(ω_g t), while for interacting atoms the variational Gaussian ground state gives a phase four orders of magnitude larger at typical rubidium scattering lengths. The enhancement mechanism is the coherent-state substitution ψ=√N φ_gs, which converts the single-particle overlap into a many-body phase, and the anisotropy of the trap is essential: a symmetric trap yields zero phase.

Load-bearing premise

The whole calculation rests on the curved-space Gross-Pitaevskii equation in which the gravitational wave enters only as h(∂²_x-∂²_y) in the kinetic term, with all derivatives of h and the time-time metric component neglected; if the true non-relativistic limit of the Klein-Gordon action in TT gauge contains additional couplings, the computed phase shifts are not the ones a real condensate would experience.

Editorial extensions

If this is right

  • If the phase shift scales as Nξ(t), then increasing the condensate atom number directly boosts the gravitational-wave signal, partially compensating the tiny h≈10⁻²⁰ amplitude.
  • Interacting condensates, with the s-wave scattering length tuned via Feshbach resonances, produce a phase about 10⁴ times larger than non-interacting ones at fixed N, making them the preferred configuration for future detectors.
  • The phase amplitude is essentially independent of the gravitational-wave frequency and of the trap frequency along the propagation direction, simplifying detector design for broadband sensitivity.
  • A NOON-state superposition of two condensates would convert the many-body phase Nξ(t) into a measurable relative phase, enabling interferometric readout.
  • An anisotropic harmonic trap is a necessary condition; a condensate symmetric in the plane perpendicular to the wave direction experiences no first-order phase shift.

Reading between the lines

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

  • The derivation omits all derivatives of h and the time-time metric component; a more complete non-relativistic reduction of the Klein-Gordon action might introduce additional coupling terms that alter the phase, so the numerical predictions should be tested against a full metric expansion.
  • The interacting result relies on a variational Gaussian ansatz whose accuracy degrades as the nonlinearity grows; in the strongly interacting regime beyond the N a_s/ℓ_H ≪ 1 limit, the four-orders-of-magnitude enhancement may not hold quantitatively.
  • Extending the coherent-state analysis to include phonon excitations or parametric resonances, as studied for other detector schemes, could reveal additional amplification mechanisms beyond the first-order phase studied here.
  • A concrete experimental test would compare the predicted Nξ(t) with interferometric measurements in an anisotropic trap while sweeping the scattering length through a Feshbach resonance, checking both the sign and the magnitude of the interaction enhancement.
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

3 major / 3 minor

Summary. The paper studies the interaction of a classical gravitational wave with a trapped Bose-Einstein condensate. Starting from a non-relativistic Gross-Pitaevskii equation in the metric of a TT-gauge plane wave, the authors write an interaction Hamiltonian proportional to h(∂_x²−∂_y²), then compute the first-order fidelity amplitude F(t)=1−iNξ(t) for a coherent condensate state. For an anisotropic harmonic trap, they derive phase-shift formulas for non-interacting (Eq. (25)) and interacting (Eq. (28)) condensates, concluding that interacting condensates give a phase shift about four orders of magnitude larger and that the amplitude is independent of the gravitational-wave frequency. They also discuss a NOON-state detection scheme.

Significance. The problem is timely and the use of the fidelity amplitude to expose an N-fold many-body phase is a clean way to frame the question. The paper does not introduce fitted parameters or invented entities, and it identifies a concrete experimental knob (trap anisotropy and scattering length). However, the central quantitative results contain algebraic errors: Eq. (23) mis-evaluates the Gaussian integral, Eq. (25) omits the 1/ω_g factor from the time integral, and Eq. (28) is presented without derivation and is not dimensionally consistent. As a result, the claimed frequency independence, the quoted magnitudes, and the four-order-of-magnitude interacting enhancement are not established by the manuscript as written.

major comments (3)
  1. [V A 1, Eq. (23)] For the normalized Gaussian ground state of Eq. (21), the spatial integral in Eq. (19) evaluates to ∫ d³r φ_gs* (∂_y²−∂_x²) φ_gs = 1/2(1/σ_x²−1/σ_y²) = (m/(2ℏ))(ω_x−ω_y). Equation (23) instead displays a prefactor proportional to (σ_y−σ_x)/(√π σ_x σ_y σ_z) with a spurious √π in the denominator, and the expression as written has the wrong physical dimensions. Because Eq. (23) is the immediate input to Eq. (25), this algebraic error propagates into the central result.
  2. [V A 1, Eq. (25)] With the corrected integral, inserting h(t,z)=h0 cos(k_g z−ω_g t) gives ∫ dz |φ_z(z)|² cos(k_g z−ω_g t) = cos(ω_g t) exp(−k_g²σ_z²/4), and ∫_0^t dt′ cos(ω_g t′) = sin(ω_g t)/ω_g. The resulting phase amplitude is h0 (ω_x−ω_y)/(4ω_g) exp(−ℏω_g²/(4mc²ω_z)), not the expression in Eq. (25). The text's statement that the amplitude is independent of ω_g therefore follows from an omitted factor in the time integral, and the numerical magnitude quoted in Fig. 1 is incorrect.
  3. [V A 2, Eq. (28)] Equation (28) is presented without derivation, and as printed its prefactor is not dimensionless (the factor [ℏ^{11/5}/(m²√Γω_z)]^{2/5} carries residual dimensions). Since the 'four orders of magnitude' interacting enhancement and Fig. 2 are based on Eq. (28), the central comparison between interacting and non-interacting condensates is unsupported.
minor comments (3)
  1. [General] There are several typos: 'dilue' should be 'dilute', 'multimessanger' should be 'multimessenger', 'contibution' should be 'contribution', and the text in Sec. V contains an apparent stray glyph before 'away from 1'.
  2. [Sec. IV and V A 2] The validity condition for the mean-field GPE is stated as 'γ ≪ 1' in Sec. IV but later as 'N a_s/ℓ_H ≪ 1' in Sec. V A 2; the manuscript should use a single, dimensionally consistent criterion.
  3. [Figures] The figures referenced in the text are not visible in the manuscript body provided; if they are part of the submission, they should be placed near their first citation so that the claimed magnitudes can be checked against the corrected formulas.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the GW-induced many-body phase is derived by first-order perturbation theory on an explicit ground state; the only self-citation supplies variational widths and is not load-bearing.

full rationale

The central derivation is self-contained. Eq. (18)-(19) follow by expanding the interaction-picture evolution operator U_int = exp(-i/hbar ∫ H_int dt') and using the coherent-state property of Eqs. (7) and (14); the phase ξ is defined directly by the flat-space ground state and the GW perturbation, with no fitted parameter or target quantity entering its definition. The non-interacting result is obtained by inserting the harmonic-oscillator Gaussian (21) into Eq. (19) and performing the integral. Although the printed algebra in Eqs. (23) and (25) appears to mis-evaluate the Gaussian integral and the time integral (the prefactor and the 1/omega_g factor are incorrect), an algebraic error is not circular reasoning: Eq. (19) is not defined in terms of Eq. (25), and the claimed frequency independence is an artifact of the error, not a constructional equivalence. The interacting case imports variational widths sigma* from Ref. [23], a prior paper by coauthor Salasnich, but those widths are obtained by minimizing the GPE energy functional, not by matching the phase shift or the enhancement factor; the citation therefore supplies an independent ground-state input rather than forcing the predicted enhancement. No uniqueness theorem, self-definition, fitted-input-renamed-as-prediction, or ansatz smuggled in as an external fact is present. The central claim is consequently not equivalent to its inputs by construction; any defect is a matter of derivation correctness, not circularity.

Assumptions & free parameters 1 free parameters · 5 assumptions · 0 invented entities

No new entities or fitted constants are introduced; the only material parameter is the scattering length. The derivation leans on standard BEC/QFT machinery and a self-cited variational formula, but the main quantitative output is compromised by apparent algebraic errors in the printed formulas.

free parameters (1)
  • s-wave scattering length a_s = ≈100 r_B for 87Rb (varied in Fig. 2)
    The interacting-BEC phase in Eq. (28) and the claimed enhancement depend on this material input; it is not fitted here but is scanned to produce the main quantitative claim.
assumptions (5)
  • domain assumption The GW is a classical weak-field perturbation in TT gauge with |h|≪1, and the lab-frame Hamiltonian follows from Eq. (11).
    Secs. II and IV; no derivation of the non-relativistic curved-space GPE from the Klein-Gordon action is shown.
  • domain assumption Contact interaction V=γδ(r-r') with γ=4πℏ²a_s/m is valid for the condensate.
    Sec. III, Eq. (8); standard dilute-gas approximation.
  • standard math Coherent-state formalism: the field is in a Glauber coherent state and the fidelity can be expanded to first order in Nh.
    Secs. III and V, Eqs. (14) and (18).
  • domain assumption The ground-state mode for the interacting case is the Gaussian ansatz with variational widths σ_i* from Ref. [23].
    Sec. V A 2, Eq. (26); the strongly-repulsive formulas are imported from prior work of one author and are not re-derived.
  • domain assumption First-order perturbation theory with the unperturbed flat-space ground state is sufficient for the fidelity phase.
    Sec. V, Eqs. (17)-(19); assumes N h0 is small and neglects back-reaction of the GW on the mode.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Interaction between gravitational waves and trapped Bose-Einstein condensates." pith.science (2026). https://pith.science/paper/ILNDQEHO

@misc{pith2026241115874,
  author       = {Pith},
  title        = {Pith review of: Interaction between gravitational waves and trapped Bose-Einstein condensates},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ILNDQEHO}},
  note         = {Machine review of arXiv:2411.15874}
}
read the original abstract

Inspired by recent proposals for detecting gravitational waves by using Bose-Einstein condensates (BECs), we investigate the interplay between these two phenomena. A gravitational wave induces a phase shift in the fidelity amplitude of the many-body quantum state. We study the enhancement of the phase shift in the case of Bose condensates confined by an anisotropic harmonic potential, considering both ideal and interacting BEC.

Figures

Figures reproduced from arXiv: 2411.15874 by the authors.

Figure 1
Figure 1. shows a two-dimensional plot depicting the amplitude of ξ(t), i.e. Nξ(t = π/(2ωg)), as a function of the confinement frequencies νx = ω/(2π) and νy = ωy/(2π) in the xy plane. The plot has been obtained considering a condensate consisting of non-interacting 87Rb atoms and keeping constant the values of ωz = 2π ×150 Hz, h0 = 10−20 and N = 107 ; and shows how the phase amplitude induced by the GW and enhanced by the nu… view at source ↗
Figure 2
Figure 2. shows the amplitude of Nξ(t), for a condensate of interacting (red line) and non-interacting (blue line) 87Rb atoms, as a function of the confinement frequency along y. Assuming a value of h0 = 10−20 for the GW and N = 107 , the amplitude of the interacting BEC is of the order of 10−13, which in this case also confirms the approximation made by expanding the time evolution operator in a series. For typical values of… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

29 extracted references · 28 canonical work pages

  1. [23]

    Fano, Sullo spettro di assorbimento dei gas nobili presso il limite dello spettro d’arco

    U. Fano, Sullo spettro di assorbimento dei gas nobili presso il limite dello spettro d’arco. Il Nuovo Cimento 12, 154 (1935)

    work page 1935

  2. [1]

    We can now insert Eq

    Non-interacting bosons For non-interacting bosons in the condensate, the ground-state wavefunction is given by ϕgs(⃗ r) = 1 (π3σ2xσ2yσ2z ) 1 4 exp − x2 2σ2x + y2 2σ2y + z2 2σ2z , (21) where σi = r ℏ mωi with i = x, y, z (22) The choice of a Gaussian shape trial function is well justified in the limit of weak interacting bosons, due to the fact that the ex...

    work page

  3. [2]

    Ultracold atoms in curved geometries

    Interacting bosons The GPE is not easily solvable analytically, but nevertheless it is possible to make use of approximate methods to evaluate also the interacting case in which γ ̸= 0. Remind that the GPE is reliable if the nonlinear term is sufficiently small. In the case of harmonic confinement this condition means N as/ℓH ≪ 1, where ℓH is the geometri...

    work page 2023

  4. [3]

    Roitberg, Einstein field equations for Bose-Einstein condensates in cosmology

    A. Roitberg, Einstein field equations for Bose-Einstein condensates in cosmology. J. Phys.: Conf. Ser. 1730, 012017 (2021)

    work page 2021

  5. [4]

    Roitberg, Emergence of effective Lorentzian metric from a vortex defect

    A. Roitberg, Emergence of effective Lorentzian metric from a vortex defect. AIP Conference Proceedings 2872, 012009 (2023)

    work page 2023

  6. [5]

    Liberati, Analogue gravity models of emergent gravity: lessons and pitfalls

    S. Liberati, Analogue gravity models of emergent gravity: lessons and pitfalls. J. Phys.: Conf. Ser. 880, 012009 (2017)

    work page 2017

  7. [6]

    B. P. Abbott, R. Abbott, T. D. Abbott, M. R. Abernathy, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, et al., Observation of gravitational waves from a binary black hole merger. Phys. Rev. Lett. 116, 061102 (2016)

    work page 2016

  8. [7]

    B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, et al., GW170814: a three-detector observation of gravitational waves from a binary black hole coalescence. Phys. Rev. Lett. 119, 141101 (2017)

    work page 2017

Show all 29 references
  1. [8]

    Sch¨ utzhold, Interaction of a Bose-Einstein condensate with a gravitational wave

    R. Sch¨ utzhold, Interaction of a Bose-Einstein condensate with a gravitational wave. Phys. Rev. D 98, 105019 (2018)

    2018

  2. [9]

    Sab ´ ın, D

    C. Sab ´ ın, D. E. Bruschi, M. Ahmadi, and I. Fuentes, Phonon creation by gravitational waves. New J. Phys. 16, 085003 (2014)

    2014

  3. [10]

    M. P. G. Robbins, N. Afshordi, A. O. Jamison, and R. B. Mann, Detection of gravitational waves using parametric resonance in Bose-Einstein condensates. Class. Quantum Grav. 39, 175009 (2022)

    2022

  4. [11]

    Gallerati, G

    A. Gallerati, G. Modanese, and G.A. Ummarino, Interaction Between Macroscopic Quantum Systems and Gravity. Frontiers in Phys. 10, 941858 (2022)

    2022

  5. [12]

    Sen and S

    S. Sen and S. Gangopadhyay, Probing the quantum nature of gravity using a Bose-Einstein condensate, Phys. Rev. D 110, 2, 026014 (2024)

    2024

  6. [13]

    Wald, General Relativity (Chicago Univ

    R.M. Wald, General Relativity (Chicago Univ. Press, 1984)

    1984

  7. [14]

    Salasnich, Quantum Physics of Light and Matter (Springer, 2017)

    L. Salasnich, Quantum Physics of Light and Matter (Springer, 2017)

    2017

  8. [15]

    R. J. Glauber, Coherent and incoherent states of the radiation field. Phys. Rev. 1963, 131, 2766

    1963

  9. [16]

    E. C. G. Sudarshan, Equivalence of semiclassical and quantum mechanical descriptions of statistical light beams. Phys. Rev. Lett. 10, 277 (1963)

    1963

  10. [17]

    E. P. Gross, Structure of a quantized vortex in boson systems. Nuovo Cim. (1955-1965) 20, 454 (1961)

    1961

  11. [18]

    E. P. Gross, Hydrodynamics of a superfluid condensate. J. Math. Phys. 4, 195 (1963)

    1963

  12. [19]

    L. P. Pitaevskii, Vortex lines in an imperfect Bose gas. Sov. Phys. JETP 13, 451 (1961)

    1961

  13. [20]

    Fagnocchi, S

    S. Fagnocchi, S. Finazzi, S. Liberati, M. Kormos, and A. Trombettoni, Relativistic Bose-Einstein condensates: a new system for analogue models of gravity, New J. Phys. 12, 9, 095012 (2010)

    2010

  14. [21]

    Feshbach, Unified theory of nuclear reactions

    H. Feshbach, Unified theory of nuclear reactions. Ann. Phys. 5, 357 (1958). 8

    1958

  15. [22]

    Feshbach, A unified theory of nuclear reactions

    H. Feshbach, A unified theory of nuclear reactions. II. Ann. Phys. 19, 287 (1962)

    1962

  16. [24]

    C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, Feshbach resonances in ultracold gases, Rev. Mod. Phys.82, 1225 (2010)

    2010

  17. [25]

    Salasnich, Time-dependent variational approach to Bose-Einstein condensation

    L. Salasnich, Time-dependent variational approach to Bose-Einstein condensation. Int. J. Mod. Phys. B 14, 1 (2000)

    2000

  18. [26]

    Giorgini, L

    S. Giorgini, L. P. Pitaevskii, and S. Stringari, Thermodynamics of a trapped Bose-condensed gas. J. Low Temp. Phys. 109, 309 (1997)

    1997

  19. [27]

    Dalfovo, S

    F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, Theory of Bose-Einstein condensation in trapped gases. Rev. Mod. Phys. 71, 463 (1999)

    1999

  20. [28]

    I. Afek, O. Ambar, and Y. Silberberg, High-NOON States by Mixing Quantum and Classical Light, Science 328, 879 (2010)

    2010

  21. [29]

    Zhang, M

    J. Zhang, M. Um, D. Lv, J.-N. Zhang, L.-M. Duan, and K. Kim, NOON States of Nine Quantized Vibrations in Two Radial Modes of a Trapped Ion, Phys. Rev. Lett. 121, 160502 (2018)

    2018

Pith tools

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