{"id":"6ab9b2f2-39c3-41ed-afc4-8c9d02490d2a","arxiv_id":"2506.21639","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A Kiselev black hole metric is realized as an acoustic spacetime in Gross-Pitaevskii theory, with scalar quasibound and quasinormal spectra computed for quintessence, dust, and radiation cases.","lead":"Using a Bose-Einstein condensate model, the authors build an acoustic analogue of a Kiselev black hole, a spacetime surrounded by an anisotropic fluid, and compute the oscillation spectra of sound waves in that geometry. A generalist might read this because tabletop acoustic setups could let laboratories probe how dark-energy-like fluids affect black hole horizons.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (15) is not obtained by inserting Eqs. (13)-(14) into Eq. (11); the radial and angular sectors carry incompatible conformal factors, so the claimed GP construction is asserted rather than derived.","rationale":"The reader's weakest assumption concerns the physical realizability of the radial velocity profile. My stress test finds a more immediate algebraic obstruction: the stated flow profile and normalization, when inserted into the paper's own acoustic metric formula, do not produce the advertised Kiselev f(r). This directly affects the claim that the Kiselev geometry is obtained 'by using the Gross-Pitaevskii theory.' The spectral analysis of the metric (15) may still be internally consistent and even novel, but the derivation connecting it to a condensate flow must be repaired or justified. I do not see this as fatal: a modified conformal factor, a corrected flow profile, or an explicit external potential could restore the construction. Therefore the verdict should remain conditional rather than accept or reject. The QBS sign error and the absence of an independent QNM benchmark are additional concerns, but the metric-derivation gap is the most load-bearing because it undermines the central construction itself.","tokens_in":17247,"tokens_out":16469,"duration_ms":190601,"concrete_test":"Symbolically substitute A = c_s^2 = 1/sqrt(3), v_mu v^mu = -1, and v_r v^r = 1 - f(r) into Eq. (11), then compute the invariant G_tt G_rr / G_theta theta^2. The claimed metric (15) gives -1/r^4; the substituted expression from Eq. (11) gives -A/[(A + 1) r^4] ~ -0.366/r^4. Also check whether the angular prefactor sqrt(A)/sqrt(A + 1) can equal sqrt(3A); it cannot for real positive A. Either check settles whether the Kiselev acoustic line element is actually derived from the GP action as written.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The load-bearing step is the transition from the GP-derived acoustic line element (11) to the Kiselev line element (15). Substituting the Minkowski background, the normalization v_mu v^mu = -1, and v_r v^r = D/r - C_w/r^(3w+1) = 1 - f(r) from Eq. (13), with A = c_s^2, Eq. (11) gives G_tt = -P(A + f - 1)/(A + 1), G_rr = P A/(A + f - 1), and G_theta theta = P r^2, where P = sqrt(A)/sqrt(A + 1). For Eq. (15) to follow, these must equal a common factor times (-f, 1/f, r^2). But the invariant G_tt G_rr / G_theta theta^2 equals -A/[(A + 1) r^4], whereas Eq. (15) requires -1/r^4. With A = 1/sqrt(3), this ratio is about -0.366/r^4, not -1/r^4. Equivalently, matching the angular prefactor requires sqrt(A)/sqrt(A + 1) = sqrt(3A), which has no positive real solution. Thus the central metric is not a consequence of the stated GP substitution; it is an independent ansatz. This is more fundamental than the realizability of the flow: even the formal algebra does not connect Eqs. (11)-(14) to Eq. (15). A separate sign error in the dust QBS formula (44), noted by the reader, is secondary to this derivation gap.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes an acoustic analogue of the Kiselev (quintessence) black hole based on the Gross–Pitaevskii action. It claims to derive the line element (15) with f(r) = 1 - D/r + C_w/r^(3w+1) and then computes quasibound-state frequencies with the VBK/Heun method and quasinormal modes with sixth-order WKB for the quintessence, dust, and radiation cases. The stated goal is to provide experimentally testable signatures of an analogue spacetime with an anisotropic-fluid background.","tokens_in":17640,"tokens_out":16646,"duration_ms":164697,"significance":"If the construction were valid, the paper would offer a unified analogue-gravity framework for a family of Kiselev-like acoustic black holes and exact spectral predictions that could be compared with BEC experiments. The manuscript contains self-contained Heun-function solutions, detailed WKB convergence tests, and explicit spectral plots for the proposed line element. However, the central derivation connecting the Gross–Pitaevskii action to the Kiselev metric is algebraically incorrect, and the dust quasibound formula is internally inconsistent with the paper's own stability statements. As a result, the main claims are not currently established.","major_comments":[{"comment":"The line element (15) is not a consequence of substituting Eqs. (13) and (14) into Eq. (11). With the Minkowski background (12), A = c_s^2, and v_mu v^mu = -1, the angular component obtained from Eq. (11) is G_theta theta = c_s (A+1)^{-1/2} r^2, while Eq. (15) requires G_theta theta = sqrt(3A) r^2. Equality would demand 1/sqrt(A+1) = sqrt(3), i.e. A = -2/3, which no positive speed of sound satisfies. Equivalently, the ratio G_tt G_rr / G_theta theta^2 from Eq. (11) is -A/[(A+1) r^4], whereas Eq. (15) gives -1/r^4. Thus the KABH metric is an independent ansatz, not derived from the Gross–Pitaevskii action. This invalidates the central construction claim of the paper.","section":"II, Eqs. (11)–(15)"},{"comment":"The dust QBS formula (44) has a sign problem relative to the stability interpretation. For D > C_m, Eq. (44) gives Im omega_n^{(m)} = (n+1)/(2(D - C_m)) > 0, which Sec. IV itself classifies as unstable for the quintessence case. Yet Sec. IV states that the dust QBSs are stable for -1 <= C_m <= +1 with D = 2 and for 0 <= D <= 2 with C_m = -1, both parameter regions satisfying D > C_m. Conversely, for C_m = +1 and 0 < D < 1, Eq. (44) gives negative imaginary parts, but Sec. IV calls the system unstable. The sign in Eq. (44) must be reversed to be consistent with the stability conclusions and with Fig. 2.","section":"III A 2 and IV"},{"comment":"The radial velocity profile in Eq. (13) and the normalization in Eq. (14) are imposed ad hoc so that the effective metric matches the Kiselev form, but the paper does not demonstrate that the Gross–Pitaevskii equation with any trapping potential or boundary condition admits a stationary, irrotational, and real flow with this v_r on the exterior of the horizon. Without such a realizability argument, the proposed experimental signatures are not tied to a concrete condensate setup.","section":"II, Eq. (13)"}],"minor_comments":[{"comment":"There are several typographical issues: 'a n experimental setup' in the abstract, 'spherical ly symmetric' in the abstract, and 'analog gravity models' in the introduction should be corrected.","section":"Abstract and I"},{"comment":"The QBS frequencies in Eqs. (34), (44), and (55) do not depend on the angular separation constant lambda; for spherically symmetric black-hole perturbation theory one would normally expect the spectrum to depend on the multipole number. The authors should explain this feature or verify the polynomial conditions used to derive these formulas.","section":"III A"},{"comment":"The sentence 'All QBSs are overdamped' is inconsistent with the positive imaginary parts of the quintessence QBSs and of the dust QBS for D > C_m; for those modes the imaginary part is positive, not damped.","section":"IV"},{"comment":"In Eq. (58), the symbol w^2 should presumably be omega^2 to match the notation used throughout the paper.","section":"III B, Eq. (58)"}],"recommendation":"reject","confidential_remarks":"The manuscript is within the journal's scope, but the central GP-to-Kiselev derivation is algebraically invalid; the authors would need to provide a genuine derivation or substantially reframe the paper as studying perturbations of an ad hoc metric. The sign inconsistency in the dust QBS sector further undermines the spectral claims. I cannot recommend publication in the present form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the thing you should know: the central construction doesn't close. The claimed transition from the GP effective metric (11) to the Kiselev acoustic metric (15) is not an algebraic consequence of the substitutions given. If you plug the Minkowski background, v^r v_r = 1 - f(r), and v^mu v_mu = -1 into (11), the angular component forces a conformal prefactor P = sqrt(c_s^2)/sqrt(c_s^2+1), while the desired metric needs a different prefactor; the tt and rr components then fail to match -f and 1/f. The paper just asserts Eq. (15). This is load-bearing, because the physical claim is that a GP condensate with the specified flow realizes the Kiselev geometry. As written, that is not demonstrated.\n\nWhat the paper does well: if you take (15) as an ansatz, the treatment of scalar perturbations is standard and competently executed. The Heun-function QBS solutions for w = -2/3, 0, 1/3 are explicit, and the WKB QNM analysis is carried out with enough detail to reproduce. The idea of a Kiselev-family analogue metric is a reasonable extension of the catalogue.\n\nSoft spots in proportion: (1) the derivation gap above is fatal to the 'analogue' interpretation; the metric is an independent ansatz that happens to look like Kiselev. It might be salvageable as a toy model if the authors explicitly disclaim the GP derivation, but then the abstract's claim to 'demonstrate ... mimicked by experimental setup' is too strong. (2) The dust QBS formula (44) has a sign that contradicts their own Fig. 2 and the stability statements in Sec. IV. For D > 1 and Cm = +1, (44) gives positive Im(omega) (growing modes) while the text says stable; for Cm = -1 it predicts instability for all D, again contradicting the text. Either the formula or the conclusions are wrong, and the reader can't tell which. (3) Minor: the WKB results would benefit from an independent check, but that is not essential.\n\nThe citation pattern is honest; the relevant acoustic black hole references are present. There is no invented data or fitted constants—the spectra follow from the metric.\n\nBottom line: this deserves referee time because the target problem is sensible and the spectral machinery is useful, but as submitted it is not correct in its load-bearing logic. I would condition acceptance on a fixed derivation or a reframe as an effective metric model, and a corrected dust spectrum.","headline":"The claimed GP-to-Kiselev derivation doesn't close algebraically, so the paper's central analogy is not supported as written, though the spectral calculations for the ansatz metric are competent.","tokens_in":18187,"tokens_out":8318,"would_cite":false,"duration_ms":87711,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A radial flow in a Bose-Einstein condensate can reproduce a Kiselev black hole, with testable sound spectra.","keywords":["acoustic black holes","analogue gravity","Kiselev spacetime","quintessence","Gross-Pitaevskii equation","Bose-Einstein condensate","quasinormal modes","quasibound states"],"falsifier":"Measure the dispersion of density perturbations in a condensate with the prescribed draining flow. If the effective metric is not the Kiselev one, the quasibound frequencies will not be purely imaginary with the values $\\omega_n^{(q)}=-iC_q n(n+2)/(2(n+1))$, $\\omega_n^{(m)}=i(n+1)/(2(D-C_m))$, and $\\omega_n^{(r)}=-i(n+1)/(2D)$; a single real part in the quasibound spectrum would already falsify the prediction.","tokens_in":17052,"feed_emoji":"🔊","tokens_out":8738,"duration_ms":90795,"temperature":0.7,"pith_summary":"The paper sets out to show that a spherically symmetric black hole surrounded by Kiselev's anisotropic fluid can be emulated in the laboratory as an acoustic black hole. Working from the Gross-Pitaevskii action, the authors derive an effective line element with metric function $f(r)=1-D/r+C_\\varpi/r^{3\\varpi+1}$, which interpolates between quintessence, dust, and radiation analogues. On this background they solve the massless Klein-Gordon equation, obtaining exact quasibound state frequencies and WKB quasinormal mode frequencies for the three matter cases. If the analogy is right, these frequencies are the ``sound of quintessence'' and provide experimental targets for analogue gravity.","feed_headline":"Hearing quintessence: a condensate mimics a Kiselev black hole","feed_subtitle":"Sound in a superfluid can mimic a quintessence black hole; the predicted frequencies are lab-testable.","key_machinery":"The load-bearing mechanism is the acoustic metric $G_{\\mu\\nu}$ extracted from the phase-fluctuation wave equation (8). Its angular and radial pieces are fixed by the background metric, while the time-radial block is fixed by the fluid four-velocity $v^\\mu$; the choices (13)-(14) make the combination $\\frac{c_s^2-v_r v^r}{c_s^2-v_\\mu v^\\mu}g_{tt}$ equal to $-f(r)$, exactly reproducing the Kiselev function. The spectral analysis then runs on two tools: the VBK approach, which rewrites the radial Klein-Gordon equation as a Heun equation and yields closed-form quasibound frequencies, and the WKB formula (58), which converts the peak of the effective potential $V(r)$ into quasinormal frequencies.","core_discovery":"The central claim, stated on the authors' own terms, is that the metric describing a Kiselev black hole surrounded by a fluid with equation-of-state parameter $\\varpi$ is an effective geometry for sound in a Bose-Einstein condensate. The construction starts with the Gross-Pitaevskii action (1), uses the Madelung representation $\\phi=\\sqrt{\\rho}e^{i\\theta}$, and identifies the phase fluctuations $\\theta_1$ with a massless scalar field. Choosing the radial four-velocity component $v_r\\sim\\sqrt{D/r-C_\\varpi/r^{3\\varpi+1}}$ and imposing $v_\\mu v^\\mu=-1$ in the critical-temperature limit converts the acoustic metric (11) into the Kiselev line element (15)-(16). The paper then derives the scalar-field spectrum: quasibound frequencies from Heun-function solutions (Eqs. (34), (44), (55)) and quasinormal frequencies from a sixth-order WKB approximation (Figs. 6, 9, 12), and reads stability off the sign of $\\mathrm{Im}\\,\\omega$.","pith_inferences":["A straightforward laboratory check would impose the radial flow (13) and measure density-perturbation spectra; agreement with the predicted purely imaginary frequencies would confirm that the condensate ``hears'' the Kiselev horizon.","The paper leaves open whether the domain of real $v_r$ covers the exterior region for all parameter choices; a reader building the experiment should first check $D/r-C_\\varpi/r^{3\\varpi+1}\\ge 0$ for $r>r_+$, since the analogue metric is only physical where the flow is real.","The same velocity-metric identity suggests a recipe for other static backgrounds: any desired metric function $f(r)$ can be fed into Eq. (16) and translated into a radial flow, yielding acoustic analogues of regular Kiselev black holes or modified-gravity solutions.","Because the quasibound frequencies have zero real part, a measurement of decay rates alone could identify the effective charge parameter $C_\\varpi$ without resolving oscillations."],"forward_implications":["The same line element (15) reduces to previously known acoustic geometries when $C_\\varpi=0$, and to a Reissner-Nordström-type acoustic metric when $\\varpi=1/3$, so the construction is a unified catalogue rather than a single special case.","All three quasibound spectra are purely imaginary, meaning every mode is overdamped; quintessence QBSs have positive imaginary part (unstable), while radiation and most dust QBSs have negative imaginary part (stable).","In the WKB quasinormal spectrum, increasing $|C_q|$ (quintessence) or $C_m$ (dust) makes the modes decay faster and oscillate faster, whereas increasing $C_r$ (radiation) slows the decay while the modes remain stable.","The horizon structure is the boundary-condition anchor: quintessence and radiation give two acoustic horizons, dust gives one, and the quasibound and quasinormal frequencies depend explicitly on those horizon radii."],"supporting_citations":[{"why":"Introduces acoustic black holes in fluid flows, the concept this paper extends to a Kiselev-like metric.","marker":"[28]"},{"why":"Derives acoustic black holes in Bose-Einstein condensates from the Gross-Pitaevskii equation, the platform used here.","marker":"[35]"},{"why":"Defines the Kiselev black hole with surrounding anisotropic fluid, the spacetime being mimicked.","marker":"[41]"},{"why":"Develops the VBK Heun-function method for black-hole resonant frequencies used to obtain quasibound states.","marker":"[76]"},{"why":"Applies the VBK method to Schwarzschild acoustic black holes, providing the template for quasibound boundary conditions.","marker":"[77]"},{"why":"Supplies the Heun solution and polynomial condition used for the quintessence quasibound spectrum.","marker":"[85]"},{"why":"Supplies the confluent Heun solutions used for the dust and radiation quasibound spectra.","marker":"[86]"},{"why":"Provides the automatic higher-order WKB code used to compute quasinormal frequencies.","marker":"[87]"},{"why":"Gives the sixth-order WKB formula (58) used for the quasinormal mode computation.","marker":"[88]"}],"fun_headline_variants":["Sound mimics a quintessence black hole in a BEC","Acoustic Kiselev black hole from a condensate","Lab acoustic hole mimics a Kiselev black hole","Superfluid sound simulates a quintessence black hole","Bose-Einstein condensate sings a Kiselev black hole"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The construction assumes a real condensate can support the stationary, irrotational radial flow $v_r\\approx\\sqrt{D/r-C_\\varpi/r^{3\\varpi+1}}$ with $v_t$ fixed by $v_\\mu v^\\mu=-1$, while remaining in the critical-temperature limit $m^2\\to 0$; if such a flow cannot be produced, the acoustic metric is formal rather than experimental.","fun_headline_variants_meta":{"raw":{"variants":["Sound mimics a quintessence black hole in a BEC","Acoustic Kiselev black hole from a condensate","Lab acoustic hole mimics a Kiselev black hole","Superfluid sound simulates a quintessence black hole","Bose-Einstein condensate sings a Kiselev black hole"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000332,"raw_usage":{"total_tokens":1850,"prompt_tokens":953,"completion_tokens":897,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":569,"completion_tokens_details":{"reasoning_tokens":823}},"tokens_in":569,"tokens_out":897,"duration_ms":8944,"temperature":1.0,"reasoning_tokens":823,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T22:42:19.190156+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the dispersion of density perturbations in a condensate with the prescribed draining flow. If the effective metric is not the Kiselev one, the quasibound frequencies will not be purely imaginary with the values $\\omega_n^{(q)}=-iC_q n(n+2)/(2(n+1))$, $\\omega_n^{(m)}=i(n+1)/(2(D-C_m))$, and $\\omega_n^{(r)}=-i(n+1)/(2D)$; a single real part in the quasibound spectrum would already falsify the prediction.","supporting_citations":[{"cited_title":"On destabilising quasi-normal modes with a radially concentrated perturbation,","cited_arxiv_id":null,"evidence_quote":"Introduces acoustic black holes in fluid flows, the concept this paper extends to a Kiselev-like metric."},{"cited_title":"Experimental black-hole evaporation?,","cited_arxiv_id":null,"evidence_quote":"Derives acoustic black holes in Bose-Einstein condensates from the Gross-Pitaevskii equation, the platform used here."},{"cited_title":"Event horizons and erg oregions in He-3,","cited_arxiv_id":null,"evidence_quote":"Defines the Kiselev black hole with surrounding anisotropic fluid, the spacetime being mimicked."},{"cited_title":"Thermodynamics and remnants of Kiselev black holes in rainbow gravity,","cited_arxiv_id":null,"evidence_quote":"Develops the VBK Heun-function method for black-hole resonant frequencies used to obtain quasibound states."},{"cited_title":"Joule-Thomson Expansion and Optical Behaviour of Reissner-Nordström-Anti-de Sitte r Black Holes in Rastall Gravity Surrounded by a Quintessence Field,","cited_arxiv_id":null,"evidence_quote":"Applies the VBK method to Schwarzschild acoustic black holes, providing the template for quasibound boundary conditions."},{"cited_title":"Eine Verallgemeinerung der Quantenbedin gungen für die Zwecke der Wellen- mechanik,","cited_arxiv_id":null,"evidence_quote":"Supplies the Heun solution and polynomial condition used for the quintessence quasibound spectrum."},{"cited_title":"Wellenmechanik und halbzahlige Quanti sierung,","cited_arxiv_id":null,"evidence_quote":"Supplies the confluent Heun solutions used for the dust and radiation quasibound spectra."},{"cited_title":"La mécanique ondulatoire de Schrödinger ; une méthode générale de résolution par approximations successives,","cited_arxiv_id":null,"evidence_quote":"Provides the automatic higher-order WKB code used to compute quasinormal frequencies."},{"cited_title":"Structure of a quantized vortex in boson sy stems,","cited_arxiv_id":null,"evidence_quote":"Gives the sixth-order WKB formula (58) used for the quasinormal mode computation."}],"review_version":1}