REVIEW 3 major objections 6 minor 26 references
Acoustic black holes, white holes, and wormholes in Bose-Einstein condensates in two dimensions
T0 review · 3 major / 6 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read Smooth 2D condensate flows can form acoustic black holes, white holes, and a one-way wormhole, with a calculable Hawking temperature in uniform-density cases.
desk verdict A solid uniform-density analogue black-hole result is paired with a speed-of-sound error that undermines the non-uniform wormhole claim. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the reduced radial equation for the condensate amplitude $\rho(R)$ in a conformally flat two-dimensional metric, Eq. (2.9): $\rho'' + \rho'/R - B^2/(\rho^3 R^2) + f(R)[(1-V(R))\rho - g_0(R)\rho^3]=0$, where $B$ fixes the radial flow. From a solution one computes the flow velocity $v = 2B/(f(R) R \rho^2)$ and the local sound speed $c \simeq \sqrt{n g_0/m}$; an acoustic horizon is a radius where $|v|=c$. A second central identity is the correspondence principle, Eqs. (2.10)-(2.11): two triples $(f,V,g_0)$ share the same stationary solutions when $f(1-V)$ and $f g_0$ match, which lets the funnel-metric wormhole be redrawn as a flat-space variable-coupling black/white hole. The fluctuation analysis uses the hydrodynamic approximation $\nabla\cdot(\rho_0^2 \nabla N_1)\approx 0$ to turn the coupled linearized equations into a massless scalar in the acoustic metric.
What would settle it
Evolve the full linearized density-phase equations (the uniform-density analogues of (3.9)-(3.10)) around the uniform-density funnel solution without dropping $\nabla\cdot(\rho_0^2\nabla N_1)$; if the radial transmission amplitude across the outer horizon differs from $\exp(-2\pi |B|^3 \omega/\sqrt{B^4-R_0^4})$, the derived Hawking temperature (4.8) fails.
Extended reading notes
Core claim
The central claim is that non-singular stationary solutions of the two-dimensional Gross-Pitaevskii equation exist in which the condensate flows from one asymptotically flat region to another through a funnel-like metric $f(R)=1+(R_0/R)^4$, with the flow becoming supersonic near $R_0$; this is an acoustic one-way wormhole, with a black-hole horizon on one side and a white-hole horizon on the other. The same stationary profiles can be reinterpreted in flat space with a radially varying coupling $g_0(R)=1+(R_0/R)^4$ and a tuned external potential, yielding acoustic black/white holes without the wormhole's second asymptotic region. In the hydrodynamic approximation, linearized density and phase fluctuations combine into a single massless scalar wave equation in an acoustic metric of Painlevé-Gullstrand type. For the uniform-density subclass the acoustic horizons are $R_{p,m}=\sqrt{B^2 \pm \sqrt{B^4-R_0^4}}$, and the Hawking temperature of the outer horizon is given by Eq. (4.8).
Load-bearing premise
The acoustic metric and Hawking temperature rest on the hydrodynamic approximation $\nabla\cdot(\rho_0^2 \nabla N_1)\approx 0$ used to combine the linearized density and phase equations; the paper gives no quantitative check of this approximation for the constructed flow profiles, and for the non-uniform wormhole the background density varies significantly.
Editorial extensions
If this is right
- The funnel-metric solution regularizes the earlier singular acoustic black hole: fluid flows through a second asymptotic region instead of accumulating, so no sink or density divergence is required.
- The flat-space analog with coupling $g_0(R)=1+(R_0/R)^4$ and a tuned external potential reproduces the same stationary profiles, giving a route to laboratory experiments without curved surfaces.
- Uniform-density configurations have horizons at $R_{p,m}=\sqrt{B^2\pm\sqrt{B^4-R_0^4}}$ (extremal when $|B|=R_0$), so horizon positions are known exactly.
- Phase fluctuations in the hydrodynamic limit propagate as massless scalars in an acoustic metric of Painlevé-Gullstrand type, allowing a semiclassical Hawking-temperature and tunneling-coefficient computation.
Reading between the lines
- Because Eq. (4.8) depends only on the ratio $R_0/B$, a uniform-density experiment could in principle scan the horizon temperature by controlling the flow parameter $B$ at fixed $R_0$; the paper does not discuss this tuning.
- The correspondence principle is not limited to the funnel profile: any two triples satisfying (2.10)-(2.11) yield identical stationary condensates, so the method may extend to other curved metrics.
- The flat-space realization needs a singular coupling and an unmodeled sink at $R=0$, so a realistic test would require a finite system with a drain or a bounded domain; the paper leaves this as an open practical step.
- The same hydrodynamic-limit acoustic-metric construction could be applied to anisotropic or non-conformally-flat backgrounds in two dimensions, where the metric is not simply $f(R)(dR^2+R^2 d\varphi^2)$.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies stationary solutions of the two-dimensional Gross-Pitaevskii equation with a conformally flat background metric or, equivalently, with position-dependent coupling and external potential. For a funnel-like metric f(R)=1+(R0/R)^4, it constructs non-singular radial flows connecting two asymptotically flat regions and interprets them as acoustic black/white hole and one-way wormhole configurations. It then linearizes density and phase fluctuations under a hydrodynamic approximation to derive acoustic metrics, identifies horizons in closed form for uniform-density solutions, and computes Hawking temperatures via analytic continuation, Euclidean periodicity, and tunneling arguments. The main quantitative results are the temperature formulas (4.8) and (6.8) and their R0→0 limit TH=1/(2π|B|), which reproduces the photon-gas result of Ref. [22].
Significance. If the claims hold, the paper provides a concrete two-dimensional BEC setting for stationary sonic black/white hole and one-way wormhole configurations, with the notable advantage that the uniform-density cases admit closed-form horizon locations and Hawking temperatures. The temperature derivations are internally cross-checked by three methods (analytic continuation, Wick-rotated periodicity, and tunneling coefficient), and the R0→0 limit correctly reduces to a known result. The paper also uses two independent numerical methods, Chebyshev collocation Newton iteration and a physics-informed neural network, which mutually corroborate the stationary profiles. However, the non-uniform wormhole claim currently rests on an incorrect definition of the local speed of sound, and the numerical evidence lacks convergence diagnostics, so the central existence claim needs revision before the paper can be accepted.
major comments (3)
- [§3.2, Eqs. (3.9)-(3.10)] The local speed of sound is defined in Eq. (3.1) as C = sqrt(2ρ(R))/sqrt(f(R)), with ρ the wave-function amplitude and n=ρ^2 the density. The linearized system (3.9)-(3.10), when combined under the hydrodynamic approximation, implies a phonon speed c_s = sqrt(2)ρ0/sqrt(f), and the acoustic metric (3.11) is consistent with that: its g00 vanishes at ρ0^6 f R^2 = 2B^2. The crossing condition C=|V| used in Fig. 3 therefore reduces to ρ0^5 f R^2 = 2B^2, which differs from the true horizon condition by a factor ρ0. For the non-uniform wormhole solution, where ρ0 is less than 1 in the throat region, the plotted crossings do not coincide with the horizons of the derived phonon metric, and it is not established that any true crossing exists for the presented parameters. The same factor-ρ error appears in Section 5, where the variable-coupling speed is written as C=√(2fρ) instead of C=ρ√(2f). Because the Conclusions claim that all configurations exhibit a crossing and hence wormhole configurations, this issue is load-bearing for the central claim of the paper.
- [§3.2, Eqs. (3.9)-(3.10)] The acoustic metric derivation relies on the approximation ∇·(ρ0^2∇N1)≈0, but no quantitative estimate of its validity is provided. For the non-uniform wormhole solution, the background density varies substantially through the throat, and the same approximation is used for the uniform-density cases, where it is justified only at long wavelengths. The manuscript should either estimate the neglected terms on the actual numerical solutions or explicitly state the low-frequency regime in which the acoustic-metric description applies. Without this, the Hawking-temperature predictions (4.8) and (6.8) remain uncontrolled approximations.
- [§3.1, Eqs. (3.5)-(3.7)] The existence of the non-uniform wormhole solution is established numerically, but the numerical evidence is incomplete. The Chebyshev-Newton method is described with the map (3.5) and boundary conditions, yet the number of collocation points N, the Newton tolerance, the maximum residual, and the sensitivity to the mapping parameter A=10 are not reported. The neural-network comparison in Fig. 2 is only qualitative. Please report the convergence in N and the residual norms for the solutions used in Figs. 3-4, so that the reader can judge whether the claimed horizon crossings are reliably resolved.
minor comments (6)
- [§4.3, §6.3, Appendix A] The corrupted glyph sequence “Leftr⮯g⊸tl⮯ne⇒” appears in Eqs. (4.7), (6.7), and Appendix A; these should be replaced by the intended arrows or words.
- [§4.2-4.3] The diagonalized metrics are central to the temperature derivation, but the line-broken equations make the denominators ambiguous; please typeset (4.4) and (4.6) unambiguously, especially the numerator in (4.6).
- [§2] The statement that two sets of functions “have the exact same set of solutions for ρ(R)” should add the obvious qualifications that the integration constant B and the boundary conditions must also match; as written it is slightly too broad.
- [References] Reference [24] is incomplete: “S W Hawking. ‘Quantum gravity and path integrals’. In: Phys. Rev., D; (United States) (Sept. 1978)” lacks volume, issue, and page numbers.
- [§6] The statement that the variable-coupling uniform-density solutions are “perhaps the solutions that are easier to reproduce in an experimental setting” is too strong, given that the required coupling g0(R)=1+(R0/R)^4 is singular at R=0 and the construction requires an unmodeled sink at the center; the text should soften or justify this experimental-readiness claim.
- [§3.2 to §4.1] The transition from the linearized equations (3.9)-(3.10) to the acoustic metric (3.11) is not shown; displaying the resulting wave equation for θ1 and the explicit identification of the Painlevé-Gullstrand form would make the derivation easier to verify.
Circularity Check
No significant circularity: Hawking temperatures and acoustic metrics are derived from the linearized GPE without fitting; the self-citation to [11] is only for numerical methods and comparison.
full rationale
The paper's derivation chain is self-contained. Static backgrounds are obtained by solving the GPE (Eq. 2.9) numerically; no parameter is fitted to any target Hawking temperature. The acoustic metrics (3.11), (4.3), (5.5), and (6.5) are obtained from the linearized density/phase equations (3.9)-(3.10) under the standard hydrodynamic approximation, cited to Refs. [4,5] (not the authors' own work). The Hawking temperatures (4.8) and (6.8) are computed from the resulting metrics by analytic continuation and by Wick-rotated Euclidean time periodicity; the R0 -> 0 reductions reproduce the independent photon-gas result [22] as a consistency check, not as the derivation. The only self-citation, Ref. [11], supplies numerical BVP techniques and singular solutions used for comparison; it is not load-bearing for the central claims. A separate correctness caveat, not a circularity: the 'scaled local speed of sound' in Eq. (3.1) is amplitude-based, and for the non-uniform density wormhole the plotted crossing with |V| differs from the acoustic-metric horizon condition (g00 = 0) by a factor of rho0; this concerns the validity of the horizon identification, but it is not a reduction of the prediction to its inputs. The flat-space implementation also assumes a singular coupling and an unmodeled sink at R = 0, which weakens experimental claims but is not circular.
Assumptions & free parameters
free parameters (2)
- R0 =
1.4827 in Fig. 2/3 for |B| = 1; 0.1 and 10.0 in Fig. 4
- B =
|B| = 1 in sample solutions; |B| = 50 in Fig. 4(b)
assumptions (4)
- domain assumption Gross-Pitaevskii equation is the correct mean-field description of the BEC dynamics
- domain assumption Spatial metric can be taken conformally flat f(R)(dR^2 + R^2 dphi^2)
- domain assumption Hydrodynamic approximation div(rho0^2 grad N1) = 0
- standard math Standard semi-classical methods for Hawking temperature from a metric (analytic continuation, Euclidean periodicity, Bogoliubov coefficients)
Cite this review
Pith. "Pith review of Acoustic black holes, white holes, and wormholes in Bose-Einstein condensates in two dimensions." pith.science (2026). https://pith.science/paper/6UHCIMWR
@misc{pith2026241202727,
author = {Pith},
title = {Pith review of: Acoustic black holes, white holes, and wormholes in Bose-Einstein condensates in two dimensions},
year = {2026},
howpublished = {\url{https://pith.science/paper/6UHCIMWR}},
note = {Machine review of arXiv:2412.02727}
}
read the original abstract
In a previous article, we studied stationary solutions to the dynamics of a Bose-Einstein condensate (BEC) corresponding to acoustic (or Unruh) black/white holes, namely configurations where the flow becomes supersonic creating a horizon for phonons. In this paper, we consider again the Gross-Pitaevskii Equation (GPE) but looking for stationary numerical solutions in the case where the couplings are position dependent in a prescribed manner. Initially we consider a 2D quantum gas in a funnel-like spatial metric. We then reinterpret this solution as a solution in a flat metric but with spatially dependent coupling and external potential. In these solutions the local speed of sound and magnitude of flow velocity cross, indicating the existence of a supersonic region and therefore of sonic analogues of black/white holes and wormholes. We discuss the numerical techniques used. We also study phase (and density) fluctuations in these solutions and derive approximate acoustic metric tensors. For certain external potentials, we find uniform density acoustic black hole configurations and obtain their Hawking temperature.
Figures
Reference graph
Works this paper leans on
-
[22]
Proposal for an analog Schwarzschild black hole in con- densates of light
L. Liao et al. “Proposal for an analog Schwarzschild black hole in con- densates of light”. In: Phys. Rev. A 99 (2 Feb. 2019), p. 023850
work page 2019
-
[1]
Experimental Black-Hole Evaporation?
W. G. Unruh. “Experimental Black-Hole Evaporation?” In: Phys. Rev. Lett. 46 (21 May 1981), pp. 1351–1353
work page 1981
-
[2]
Particle creation by black holes
S. W. Hawking. “Particle creation by black holes”. eng. In: Commun. Math. Phys., v. 43, no. 3, pp. 199-220 43.3 (1975), pp. 199–220. issn: 0010-3616
work page 1975
-
[3]
Carlos Barcel´ o, Stefano Liberati, and Matt Visser. “Analogue Gravity”. In: Living Reviews in Relativity 14.1 (May 2011), p. 3. issn: 1433-8351
work page 2011
-
[4]
Analogue gravity from Bose-Einstein condensates
Carlos Barcel´ o, S Liberati, and Matt Visser. “Analogue gravity from Bose-Einstein condensates”. In: Classical and Quantum Gravity 18.6 (Mar. 2001), p. 1137
work page 2001
-
[5]
Analogue Models of and for Gravity
Matt Visser, Carlos Barcel´ o, and Stefano Liberati. “Analogue Models of and for Gravity”. In: General Relativity and Gravitation 34.10 (Oct. 2002), pp. 1719–1734. issn: 1572-9532. 24
work page 2002
-
[6]
Sonic Analog of Gravitational Black Holes in Bose- Einstein Condensates
L. J. Garay et al. “Sonic Analog of Gravitational Black Holes in Bose- Einstein Condensates”. In: Phys. Rev. Lett. 85 (22 Nov. 2000), pp. 4643– 4647
work page 2000
-
[7]
Caio C. Holanda Ribeiro, Sang-Shin Baak, and Uwe R. Fischer. “Exis- tence of steady-state black hole analogs in finite quasi-one-dimensional Bose-Einstein condensates”. In: Phys. Rev. D 105 (12 June 2022), p. 124066
work page 2022
Show all 26 references
-
[8]
Testing the upper bound on the speed of scrambling with an analogue of Hawking radiation using trapped ions
Zehua Tian et al. “Testing the upper bound on the speed of scrambling with an analogue of Hawking radiation using trapped ions”. In: The European Physical Journal C 82.3 (2022), p. 212. issn: 1434-6052
2022
-
[9]
“Cosmological
Petr O. Fedichev and Uwe R. Fischer. ““Cosmological” quasiparticle production in harmonically trapped superfluid gases”. In: Phys. Rev. A 69 (3 Mar. 2004), p. 033602
2004
-
[10]
Gibbons-Hawking Effect in the Sonic de Sitter Space-Time of an Expanding Bose-Einstein-Condensed Gas
Petr O. Fedichev and Uwe R. Fischer. “Gibbons-Hawking Effect in the Sonic de Sitter Space-Time of an Expanding Bose-Einstein-Condensed Gas”. In: Phys. Rev. Lett. 91 (24 Dec. 2003), p. 240407
2003
-
[11]
Stationary acoustic black hole solutions in Bose-Einstein condensates and their Borel analysis
Sachin Vaidya and Martin Kruczenski. Stationary acoustic black hole solutions in Bose-Einstein condensates and their Borel analysis . 2024. arXiv: 2411.06678 [cond-mat.quant-gas]
2024 arXiv
-
[12]
Acoustic horizon as a phase-slip surface
Sergei Khlebnikov. “Acoustic horizon as a phase-slip surface”. In: Phys. Rev. A 109 (2 Feb. 2024), p. 023305
2024
-
[13]
Observation of self-oscillating supersonic flow across an acoustic horizon in two dimensions
Hikaru Tamura et al. “Observation of self-oscillating supersonic flow across an acoustic horizon in two dimensions”. In: (Apr. 2023). arXiv: 2304.10667 [cond-mat.quant-gas]
2023 arXiv
-
[14]
Nozieres and D
P. Nozieres and D. Pines. The Theory of Quantum Liquids: Superfluid Bose Liquids . Advanced book Classics. Westview Press, pp. 160–161. isbn: 0201408414. 25
-
[15]
Nozieres and D
P. Nozieres and D. Pines. The Theory of Quantum Liquids: Superfluid Bose Liquids . Advanced book Classics. Westview Press, pp. 168–169. isbn: 0201408414
-
[16]
Nozieres and D
P. Nozieres and D. Pines. The Theory of Quantum Liquids: Superfluid Bose Liquids . Advanced book Classics. Westview Press, p. 66. isbn: 0201408414
-
[17]
Chebyshev differentiation matrices for efficient computation of the eigenvalues of fourth-order Sturm–Liouville problems
A.H. Saleh Taher, A. Malek, and S.H. Momeni-Masuleh. “Chebyshev differentiation matrices for efficient computation of the eigenvalues of fourth-order Sturm–Liouville problems”. In: Applied Mathematical Modelling 37.7 (2013), pp. 4634–4642. issn: 0307-904X
2013
-
[18]
The theory of Newton’s method
A. Gal´ antai. “The theory of Newton’s method”. In: Journal of Com- putational and Applied Mathematics 124.1 (2000). Numerical Analysis
2000
-
[19]
T.Y. Na. Computational Methods in Engineering Boundary Value Prob- lems. Mathematics in science and engineering. Academic Press, 1979, p. 85. isbn: 9780125126502
1979
-
[20]
Sonic black holes in dilute Bose-Einstein conden- sates
L. J. Garay et al. “Sonic black holes in dilute Bose-Einstein conden- sates”. In: Phys. Rev. A 63 (2 Jan. 2001), p. 023611
2001
-
[21]
The river model of black holes
Andrew J. S. Hamilton and Jason P. Lisle. “The river model of black holes”. eng. In: American journal of physics 76.6 (2008), pp. 519–532. issn: 0002-9505
2008
-
[23]
TORTOISE COORDINATE TRANSFORMATION ON APPARENT HORIZON OF A DYNAMICAL BLACK HOLE
XIANMING LIU, ZHENG ZHAO, and WENBIAO LIU. “TORTOISE COORDINATE TRANSFORMATION ON APPARENT HORIZON OF A DYNAMICAL BLACK HOLE”. In: International Journal of Modern Physics: Conference Series 12 (2012), pp. 358–367. 26
2012
-
[24]
Quantum gravity and path integrals
S W Hawking. “Quantum gravity and path integrals”. In: Phys. Rev., D; (United States) (Sept. 1978). 27 A Acoustic wormhole double-funnel In section 3 we consider a 2D quantum fluid in a conformally flat metric given by ds2 = f(R)( dR2+ R2dθ2), f (R)= 1+ R4 0 R4 (A.1) To potent...
1978
-
[26]
To find the exact required shape we consider, in cylindrical coordinates a surface parameterized by (R, θ) of the form: r = r(R), θ = θ, z = z(R) (A.2) such that the induced metric on the surface ds2 = dr2+ r2dθ2+ dz2 (A.3) = ( ∂r(R) ∂R ) 2 dR2+ r(R)2dθ2+( ∂z(R) ∂R ) 2 dR2 (A....
-
[2000]
IV: Optimization and Nonlinear Equations, pp
Vol. IV: Optimization and Nonlinear Equations, pp. 25–44. issn: 0377-0427
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.