{"id":"8b4376c6-a9c5-4c05-ada1-a0b52725cb8f","arxiv_id":"2501.15676","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Spherical accretion of a multi-species relativistic fluid with variable gamma produces an acoustic black hole metric whose surface gravity grows with flow energy and proton fraction.","lead":"Ripples in matter falling into a black hole can behave like waves in a curved spacetime, creating a sound horizon inside the flow. This paper computes that analogue gravity for a realistic mixture of electrons, positrons, and protons and shows how the sound horizon's strength changes with flow energy and composition.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The surface-gravity formula Eq. (83) is not derived from the paper's own acoustic metric; its substitutions conflate the pseudo-Newtonian background with the emergent metric, so the quantitative kappa(E_c, xi) results in Figs. 5-6 are unsupported.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: Eq. (83) is quoted from earlier work and simplified with substitutions that mix the pseudo-Newtonian background potential with the acoustic metric. My independent reading confirms this is the most serious issue. The central claim that a sonic horizon exists for the multi-species, variable-gamma flow is supported by a standard derivation: the linear perturbation equations (48)-(50) produce the Unruh-type acoustic metric (58), and the Carter-Penrose construction in Sec. VI A follows the established procedure. That part of the argument is independent of the questionable formula. However, the paper's headline quantitative results, shown in Figs. 5 and 6, are exactly the surface-gravity curves obtained from Eq. (83). If Eq. (83) is not the surface gravity of the acoustic metric, those curves have no demonstrated validity. The paper gives no derivation of Eq. (83) from its own metric, and the text's replacements are internally inconsistent: in the acoustic metric, g_rr is not -1 and g_tt vanishes at the sonic point, while sqrt(1+phi) is not a component of that metric. A first-principles recalculation of kappa from the acoustic metric would settle whether Eq. (83) is numerically close or off by a significant factor. Given that the existence claim survives, the appropriate disposition is conditional acceptance pending this verification, which matches the reader's verdict.","tokens_in":15553,"tokens_out":14148,"duration_ms":125973,"concrete_test":"Compute the acoustic surface gravity from first principles for the PW potential at E = 1.0007, xi = 1, using the standard Killing-horizon definition kappa^2 = -1/2 (nabla_mu chi_nu)(nabla^mu chi^nu) with chi = partial_t and the analytic acoustic metric (58) constructed from the numerically obtained rho0, u0, c_s0; compare this value with Eq. (83). If the two differ by more than a factor of 2, or if the first-principles expression involves g_tt or a non-unit conformal factor, Fig. 5 is quantitatively invalid; repeat at two xi values to test the claimed monotonic rise in Fig. 6.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quantitative central result, the acoustic surface gravity as a function of E_c and composition xi (Figs. 5-6), rests entirely on Eq. (83). That equation is borrowed from Refs. [9,38] and then evaluated by setting sqrt(psi^mu psi_mu) = sqrt(1+phi) and sqrt(-g_rr) = 1. These substitutions are not compatible with the acoustic metric derived in Sec. VI. For the line element (58), g_rr = rho0/c_s0, not -1, and g_tt = (rho0/c_s0)(c_s0^2 - u0^2), which vanishes at the sonic point where u0 = -c_s0. Hence the Killing-norm factor in Eq. (82) should vanish at the horizon, and no sqrt(1+phi) term arises from the acoustic metric itself. The factor introduced instead comes from the pseudo-Newtonian background potential, which played no role in constructing the acoustic spacetime. Thus Eq. (83) has no demonstrated connection to the acoustic metric built in Sec. VI, and the numerical kappa(E_c, xi) curves are unsupported. The existence of the sonic horizon and the Carter-Penrose causal structure are independent of this formula and appear robust, but the paper's new quantitative predictions depend on this unjustified identification.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper models spherically symmetric, inviscid, irrotational accretion onto a non-rotating black hole in a pseudo-Newtonian framework, using the four standard pseudo-potentials of Paczynski-Wiita, Artemova-Bjornsson-Novikov (two versions), and Nowak-Wagoner, with the Ryu-Chattopadhyay multi-species relativistic equation of state (electron-positron-proton composition xi, radially varying adiabatic index). It constructs stationary transonic solutions, classifies the sonic (critical) points as saddle-type fixed points via a dynamical-system analysis, and derives the linearized perturbation equation for the velocity potential, which yields the standard Unruh-Visser acoustic metric. Carter-Penrose diagrams are constructed and identify the sonic point as a null acoustic horizon. The acoustic surface gravity is then evaluated with a formula quoted from earlier work (Eq. (83)) as a function of the specific energy E_c and of xi for all four potentials (Figs. 5-6). Standing-wave (globally subsonic) and WKB traveling-wave (transonic) analyses are used to argue stability of the stationary solutions.","tokens_in":15873,"tokens_out":46735,"duration_ms":403488,"significance":"If the quantitative surface-gravity results were fully supported, the paper's main value would be to extend analogue-gravity treatments of spherical accretion from polytropic/isothermal flows to a multi-species relativistic equation of state with variable Gamma, with a comparison of the sonic-point structure across the four pseudo-Newtonian potentials. The positive content is real: the transonic solutions and critical-point conditions for this EoS are derived consistently and in enough detail to reproduce the phase portraits; the acoustic wave equation (49) and the metric (57)-(58) follow the standard Unruh-Visser framework; and the causal-structure identification of the sonic point as an acoustic Killing horizon is independent of the disputed surface-gravity formula. The main weakness is that the headline quantitative predictions, the kappa(E_c, xi) curves, are evaluated through a quoted formula whose connection to the metric derived in Sec. 6 is not demonstrated; the stability claim additionally relies on a WKB argument that is not uniform at the sonic point.","major_comments":[{"comment":"The quantitative results of the paper — the acoustic surface gravity as a function of E_c and xi shown in Figs. 5 and 6 — rest entirely on Eq. (83), quoted from Refs. [9,38] and evaluated with the substitutions sqrt(psi^mu psi_mu) = sqrt(1+phi) and g_rr = -1. These substitutions are not compatible with the acoustic metric derived in Sec. 6. For the line element (58), g_tt = (rho0/c_s0)(c_s0^2 - u0^2), g_tr = (rho0/c_s0)u0, and g_rr = -rho0/c_s0; hence sqrt(-g_rr) is a conformal factor, not 1, and the Killing-vector norm sqrt(psi^mu psi_mu) = sqrt(|g_tt|) vanishes at the sonic point u0^2 = c_s0^2, so the surface gravity must come from the standard limiting procedure for Killing horizons. The factor sqrt(1+phi) is the pseudo-Newtonian background potential from the Bernoulli integral (18)-(29); it does not appear anywhere in (58), so substituting it for the Killing-norm factor conflates the background gravitational model with the emergent acoustic spacetime. I note that the manuscript's own Carter-Penrose construction uses the quantity kappa = (u0' - c_s0')|_{rc} (Eq. (71)), which arises naturally from the leading-order expansion (65) of the metric function, whereas Eq. (83) adds the unjustified factor sqrt(1+phi)/(1-c_s^2). The authors should derive the surface gravity directly from the metric (58), recompute Figs. 5-6, and either justify or drop the extra factors; the existence of the acoustic horizon and the causal-structure results do not depend on this formula.","section":"§7, Eq. (83); Figs. 5–6"},{"comment":"The traveling-wave stability argument is incomplete at the sonic point. The WKB phase integrals in Eq. (95), k_{-1} = i∫dr/(u0 ∓ c_s0), have a logarithmically divergent integrand at rc because one of the two denominators u0 ± c_s0 vanishes for infall (u0 = -c_s0), so the ordering check (97)-(99), performed only for large r, does not control the expansion near the sonic point; the standard connection-formula treatment of the singular point of Eq. (90) is missing. Since Sec. 8 concludes that 'in both cases ... the stability of the steady state solution is ensured', the stated stability claim is not established by the material presented, although the linearized wave equation (49) itself is not in question.","section":"§8B, Eqs. (90)–(99)"}],"minor_comments":[{"comment":"The statement that the perturbed mass accretion rate Sigma-tilde(r,t) obeys the wave equation (53) with the coefficient matrix (54) is asserted without derivation, and the matrix (54), proportional to (u0/Sigma0), is not the same as f^mu-nu in Eq. (50), proportional to (rho0/c_s0^2); the two differ by a non-constant factor unless additional identities are supplied. Since the acoustic metric is derived from the psi-tilde equation (49), this step should be proven, attributed, or removed.","section":"§5, Eqs. (52)–(54)"},{"comment":"The symbol psi is used for the velocity potential from Eq. (46) onward and again for the Killing vector in Eq. (81), and the potential Phi in Eq. (82) becomes phi in Eq. (83) without comment; introducing xi^mu for the Killing vector and a single symbol for the pseudo-potential would remove the ambiguity.","section":"§5–§7, notation"},{"comment":"Eqs. (29) and (34) are identical and both appear in the text; one of them should be deleted.","section":"§3–§4"},{"comment":"The caption of Fig. 4 refers to a 'green, dashed region' and a 'pink, solid region', but the figure as printed is a monochrome line drawing without such labels; the regions should be marked directly in the figure, and the flow parameters (E_c, xi, potential) used for the diagram should be stated.","section":"§6A, Fig. 4"},{"comment":"The derivation of A_± and of the null coordinates drops the conformal factor rho0/c_s0 of the metric (58) without comment; a single sentence noting that null geodesics (and hence the causal structure) are unchanged by the conformal factor would make this step explicit.","section":"§6A, Eqs. (59)–(77)"},{"comment":"The numerical scheme behind Figs. 1-3 and 5-6 is not described: the integration domain, the outer-boundary value of Theta (equivalently, how E_c fixes the initial data), and the treatment of the critical-point boundary conditions are not stated, which hampers reproducibility.","section":"§2–§3, numerics"},{"comment":"The matrix entries in Eqs. (41) and (B7) contain rendering inconsistencies (Theta appears as theta in the (1,2) entries), and the final trace and determinant expressions in Eqs. (43) and (B9) should be re-checked against these matrices before publication.","section":"Eqs. (41), (B7)"},{"comment":"Refs. [8] and [32] are the same paper (T. K. Das, Class. Quantum Grav. 21, 5253 (2004)) and are cited twice under two numbers.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is an incremental extension of the authors' earlier analogue-accretion program (Refs. [8,9,13,38]) to the Ryu-Chattopadhyay equation of state; the novelty is modest but likely within scope for a specialist journal. The referee's recommendation is driven by a correctness issue, not novelty: Eq. (83) is cited to earlier papers of the same group, and because it is not re-derived for the metric (58), the paper's quantitative claims are not self-contained. I would ask the authors to derive the surface gravity from Eq. (58) and recompute Figs. 5-6, or to weaken the quantitative claims accordingly; verification of the sign conventions in Eqs. (65) and (71) would be a useful by-product of that derivation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know before you spend time on it. The paper's core—applying the standard acoustic-metric formalism to spherically symmetric Bondi accretion with the Ryu–Chattopadhyay multi-species relativistic EoS—is a solid, incremental step and the derivation in Secs. V–VI holds up. But the surface-gravity results in Figs. 5–6, which are the paper's main quantitative claims, rest on a borrowed formula (Eq. 83) that is not connected to the acoustic metric the paper itself constructs. Those κ(E_c, ξ) curves should be treated as unsupported unless the authors re-derive them.\n\nWhat's good: the EoS with variable adiabatic index is genuinely new to this analogue-gravity setting. The critical-point analysis for the four pseudo-Newtonian potentials is standard and appears consistent; the acoustic metric (Eq. 58) is the usual Unruh metric for a nonrelativistic fluid, and the Carter–Penrose diagram correctly shows the sonic horizon as a null surface. The perturbation equations in Sec. V follow the familiar route and the stability conclusions are not obviously wrong.\n\nThe soft spot is real and it's in Sec. VII. Eq. (82) is quoted from Ref. [9] without derivation, and Eq. (83) is obtained by substituting sqrt(ψψ)=sqrt(1+φ) and g_rr=-1. But the acoustic metric has g_rr = -ρ0/c_s0 and g_tt = (ρ0/c_s0)(c_s0^2-u0^2), which vanishes at the sonic point. So the Killing-norm factor in Eq. (82) should vanish there, and the sqrt(1+φ) term has no origin in the acoustic metric. This looks like an unjustified conflation of the pseudo-Newtonian background potential with the emergent spacetime. The existence of the acoustic horizon and the causal structure are unaffected, but the quantitative κ values are not.\n\nIf I were refereeing, I'd ask for a self-contained derivation of κ from Eq. (58), perhaps using the standard Killing-horizon definition adapted to the conformally related acoustic metric. The standing/traveling-wave section adds little and can be trimmed. There are also minor sign and notation slips (e.g., the metric prefactor in Eq. 57 vs. Eq. 58) that need cleanup.\n\nWho this is for: researchers working on analogue gravity in astrophysical accretion, who will find the EoS extension useful. It's not a breakthrough, but it's a legitimate incremental contribution. I'd send it to peer review with a request for major revision on the surface-gravity part.\n\nRecommendation: engage with it; fix the κ derivation and it's publishable.","headline":"Solid incremental extension of acoustic-metric accretion to a multi-species EoS, but the surface-gravity curves rest on a formula that doesn't follow from the paper's own metric.","tokens_in":16350,"tokens_out":8903,"would_cite":true,"duration_ms":72733,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Linear perturbations of spherically symmetric multi-component relativistic accretion onto a Schwarzschild black hole generate a black-hole-like acoustic spacetime whose sound horizon traps acoustic perturbations, and the associated…","keywords":["analogue gravity","acoustic metric","black hole accretion","multi-species equation of state","transonic flow","sonic horizon","surface gravity","Carter-Penrose diagram"],"falsifier":"Launching a sound pulse just inside the sonic radius in a numerical integration of the linear perturbation equations would settle the trapping claim: if any part of the pulse crosses the sonic radius outward, the acoustic horizon is not an event horizon. Recomputing $\\kappa$ at the sonic point with the exact Schwarzschild metric instead of the pseudo-Newtonian potentials would settle the surface-gravity formula: if the resulting band of $\\kappa$ values differs significantly from the four-potential results, Eq. (83) is the weak point.","tokens_in":15404,"feed_emoji":"🕳️","tokens_out":12275,"duration_ms":107130,"temperature":0.7,"pith_summary":"The paper tries to establish that the ordinary radial infall of a multi-component relativistic gas onto a non-rotating black hole carries a hidden curved spacetime within it. When the steady transonic accretion solution is slightly disturbed, the disturbance obeys the same wave equation as a massless scalar field moving on an effective metric, the acoustic metric, whose null surfaces are sonic horizons. Carter-Penrose diagrams built from the stationary solutions show that this sonic point is a genuine causal horizon that acoustic perturbations cannot cross. The paper then computes the acoustic surface gravity at that horizon and finds it increases with the conserved specific energy of the flow and with the proton fraction of the accreting gas. On this view, real astrophysical accretion flows become analogue-gravity systems in which horizon physics is an emergent, fluid-based phenomenon rather than a purely gravitational one.","feed_headline":"Acoustic horizon emerges in multi-species black hole accretion","feed_subtitle":"Sound cannot cross the flow's sonic horizon; the analogue surface gravity grows with energy and proton fraction.","key_machinery":"The central object is the acoustic metric, obtained by matching the linearized perturbation equation $\\partial_\\mu(f^{\\mu\\nu}\\partial_\\nu\\tilde{\\psi})=0$ to the curved-spacetime scalar wave equation $\\partial_\\mu(\\sqrt{-g}\\,g^{\\mu\\nu}\\partial_\\nu\\Phi)=0$, which fixes $g_{\\mu\\nu}$ up to a conformal factor and yields the line element $ds^2 = -(\\rho_0/c_{s0})[-c_{s0}^2 dt^2 + (dx^i-u_0^i\\,dt)\\delta_{ij}(dx^j-u_0^j\\,dt)]$, a Painlevé-Gullstrand form. Its null structure is charted with Carter-Penrose coordinates, and the acoustic surface gravity is extracted from the Killing-vector formula $\\kappa = \\left|\\sqrt{\\psi^\\mu\\psi_\\mu/(-g_{rr})}\\,(1-c_s^2)^{-1}\\,(du/dr-dc_s/dr)\\right|$ at $r_c$, simplified to Eq. (83). The multi-species equation of state enters through the variable polytropic index $N$ and adiabatic index $\\Gamma$, which make the sound speed, density, and critical-point data functions of temperature and composition. This machinery converts the stability question for the steady flow into a statement about the causal geometry of an emergent spacetime.","core_discovery":"The paper claims that, for spherically symmetric accretion described by any of the four post-Newtonian pseudo-Schwarzschild potentials and by a relativistic multi-species equation of state with a radially varying adiabatic index, the linear stability analysis of the stationary transonic solution yields an acoustic metric of Painlevé-Gullstrand form. The sonic point is a saddle-type critical point that becomes a null hypersurface of this acoustic spacetime, so the sound horizon behaves like an event horizon for acoustic perturbations. The acoustic surface gravity, evaluated from Eq. (83), rises monotonically with the specific energy $E_c$ at the critical point and with the proton fraction $\\xi$. The paper presents this as a classical analogue-gravity model of black-hole accretion with a more realistic thermodynamics than the polytropic or isothermal equations of state used in earlier work.","pith_inferences":["If the surface-gravity identification survives contact with full general relativity, the analogue Hawking temperature $T_H=\\kappa/2\\pi$ for each pseudo-potential could be compared with the actual Schwarzschild Hawking temperature, giving a quantitative measure of how much horizon thermodynamics is emergent in accretion flows.","The monotonic rise of $\\kappa$ with proton fraction $\\xi$ suggests that acoustic surface gravity could in principle serve as a composition diagnostic for accreting gas, though detecting an analogue temperature directly is far beyond current instruments.","Applying the same perturbation machinery to pseudo-Kerr potentials with angular momentum should produce multiple sonic points and a shock-induced acoustic white hole; testing that prediction would extend the same method without changing its core mechanism.","Replacing the pseudo-Newtonian potentials by the exact Schwarzschild metric in Eq. (83) would provide a direct check of whether the four-potential results bracket the true relativistic surface gravity; the paper does not perform this comparison."],"forward_implications":["For all four pseudo-Newtonian potentials considered, the stationary transonic solutions are stable and yield one saddle-type sonic point outside the horizon, so the acoustic-horizon result does not depend on which potential model is chosen.","The Carter-Penrose construction shows the sound horizon is a null hypersurface of the acoustic metric, so acoustic perturbations created inside the horizon cannot escape to infinity.","The acoustic surface gravity $\\kappa$ increases with the conserved specific energy $E_c$ and with the proton fraction $\\xi$, so hotter or more proton-rich flows possess a larger analogue Hawking temperature.","Standing-wave analysis for subsonic flows and WKB traveling-wave analysis for supersonic black-hole flows both give non-divergent perturbation amplitudes, supporting linear stability of the steady states.","In spherical symmetry only one acoustic horizon forms; obtaining an acoustic white hole would require axisymmetric flow with angular momentum, as the paper notes in its conclusion."],"supporting_citations":[{"why":"Establishes the original acoustic metric by showing sound perturbations propagate as scalar fields on an effective curved spacetime.","marker":"[2]"},{"why":"Formalizes the acoustic metric construction for irrotational, inviscid fluids, which the paper follows.","marker":"[3]"},{"why":"Provides the acoustic surface gravity expression used in Eq. (82).","marker":"[9]"},{"why":"Supplies the simplified surface gravity formula of Eq. (83) and its pseudo-Newtonian usage.","marker":"[38]"},{"why":"Introduces the multi-species relativistic equation of state with variable adiabatic index.","marker":"[21]"},{"why":"Provides the specific form of the equation of state used for $f(\\Theta,\\xi)$ and the polytropic index $N$.","marker":"[22]"},{"why":"Supplies the Carter-Penrose diagram construction for analogue spacetimes, used to identify sonic horizons.","marker":"[24]"},{"why":"Defines the PW pseudo-Newtonian potential used for one of the four accretion models.","marker":"[26]"},{"why":"Defines the ABN1 and ABN2 pseudo-Newtonian potentials used for two of the accretion models.","marker":"[27]"},{"why":"Defines the NW pseudo-Newtonian potential used for the remaining accretion model.","marker":"[28]"}],"fun_headline_variants":["Emergent acoustic horizon from multi-species accretion","Analogue gravity: sound horizon emerges in black hole accretion","Sonic horizon emerges from multi-component relativistic flow","Acoustic surface gravity rises with energy and proton fraction"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a surface-gravity formula borrowed from earlier analogue-accretion work, simplified by identifying the pseudo-Newtonian potential with part of the acoustic metric and setting a metric component to one, remains valid for the relativistic multi-species equation of state; if that identification fails, the specific $\\kappa(E_c,\\xi)$ values shown in the figures are unsupported, although the existence of the acoustic horizon itself would not be at stake.","fun_headline_variants_meta":{"raw":{"variants":["Emergent acoustic horizon from multi-species accretion","Analogue gravity: sound horizon emerges in black hole accretion","Sonic horizon emerges from multi-component relativistic flow","Acoustic surface gravity rises with energy and proton fraction"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000766,"raw_usage":{"total_tokens":3355,"prompt_tokens":865,"completion_tokens":2490,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":481,"completion_tokens_details":{"reasoning_tokens":2428}},"tokens_in":481,"tokens_out":2490,"duration_ms":16476,"temperature":1.0,"reasoning_tokens":2428,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T14:03:34.145928+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Launching a sound pulse just inside the sonic radius in a numerical integration of the linear perturbation equations would settle the trapping claim: if any part of the pulse crosses the sonic radius outward, the acoustic horizon is not an event horizon. Recomputing $\\kappa$ at the sonic point with the exact Schwarzschild metric instead of the pseudo-Newtonian potentials would settle the surface-gravity formula: if the resulting band of $\\kappa$ values differs significantly from the four-potential results, Eq. (83) is the weak point.","supporting_citations":[{"cited_title":"acoustic metric","cited_arxiv_id":null,"evidence_quote":"Establishes the original acoustic metric by showing sound perturbations propagate as scalar fields on an effective curved spacetime."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Formalizes the acoustic metric construction for irrotational, inviscid fluids, which the paper follows."},{"cited_title":"Barcel´ o, S","cited_arxiv_id":null,"evidence_quote":"Provides the acoustic surface gravity expression used in Eq. (82)."},{"cited_title":"Chandrasekhar, An introduction to the study of stellar structure, Vol","cited_arxiv_id":null,"evidence_quote":"Introduces the multi-species relativistic equation of state with variable adiabatic index."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the specific form of the equation of state used for $f(\\Theta,\\xi)$ and the polytropic index $N$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Carter-Penrose diagram construction for analogue spacetimes, used to identify sonic horizons."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the PW pseudo-Newtonian potential used for one of the four accretion models."},{"cited_title":"Chattopadhyay and D","cited_arxiv_id":null,"evidence_quote":"Defines the ABN1 and ABN2 pseudo-Newtonian potentials used for two of the accretion models."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the NW pseudo-Newtonian potential used for the remaining accretion model."}],"review_version":1}