{"id":"a692fca4-0e71-4c8e-935b-13c150a7e772","arxiv_id":"1908.07356","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"The maximum spherical accretion rate onto a central black hole in a Hernquist galaxy increases with the galaxy's characteristic linear density and decreases with radiation, approaching the classical Bondi rate when the outer boundary is far beyond the semi-Bondi radius.","lead":"This paper recalculates Bondi accretion onto a central black hole using boundary conditions at a finite radius instead of infinity, then adds the gravity and radiation of the host elliptical galaxy. It finds that the maximum accretion rate rises with the galaxy's linear density and falls with luminosity, and approaches the classic Bondi value when the outer boundary is far away.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The §4 result that λ_c decreases with l and l_g depends on the optically-thin Thomson radiation force of Eq. 29; if absorption or multiple scattering matters, the (1-l)(1-lg) factorization driving Figs. 5-6 breaks down.","rationale":"The reader's weakest assumption matches my own: the Eq. 29 radiation force is the main bridge from a pure Bondi calculation to the paper's astrophysical conclusions about how λ_c depends on galaxy density and luminosity. I agree with the conditional verdict: the algebra and figures are internally consistent under the stated assumptions, but those assumptions are strong and are explicitly flagged by the authors as future work. The finite-radius results of Sections 2-3 and the transonic-saddle discussion of Section 5 are largely independent of this issue, so even if the radiation-force model were wrong, the first part of the paper would stand. A radiative-transfer check can settle whether the concern actually lands. The reader's other concerns—missing quantitative comparison with Korol et al. and Ciotti & Pellegrini, the unproven saddle-type conclusion, and mechanical issues—are secondary and do not change my assessment. No new concern moves the verdict.","tokens_in":13365,"tokens_out":27204,"duration_ms":270946,"concrete_test":"Take the transonic density and temperature profiles from Eqs. 8 and 11 for the λ_c solutions underlying Fig. 6 (γ=1.4, R_f→∞, several m_g/R_g, l, l_g). With a physical normalization (e.g., ρ_f=10^-23 g cm^-3, c_sf=300 km/s), compute the Thomson and absorption optical depths from r_f to r_c using standard ionized-gas opacities. If τ_es+τ_abs ≥ 0.1, re-solve the momentum equation (Eq. 22) with a frequency-dependent or Monte Carlo radiation force including absorption and re-emission; if the resulting λ_c(m_g/R_g,l,l_g) differs from Fig. 6 by more than about 10%, the central claim of §4 is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 4's key result—that λ_c increases with m_g/R_g and decreases with l and l_g—is obtained by inserting f_rad = G M_BH [l/r² + m_g l_g/(r+r_g)²] (Eq. 29) into the momentum equation and absorbing luminosity into effective potential factors (1-l) and (1-lg) in g(R) (Eq. 32). This force law is exact only for an optically thin, Thomson-scattering-dominated, spherically symmetric medium with constant central luminosity and a transparent Hernquist stellar radiation field. In real galactic nuclei, dust absorption, line driving, and multiple electron scatterings alter both the magnitude and radial dependence of f_rad; the simple factorization into (1-l) and (1-lg) then no longer holds, so the monotonic and linear dependencies of λ_c on l, l_g, and m_g/R_g shown in Fig. 6 are not guaranteed. The manuscript itself lists the optically thin restriction as a limitation and defers absorption to future work (§6), so the concern is acknowledged but unresolved.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper revisits steady, spherically symmetric polytropic Bondi accretion with boundary conditions imposed at a finite outer radius r_f rather than at infinity. It derives the modified Bernoulli equation and the critical accretion parameter λ_c in dimensionless form (Eqs 12-19), showing numerically that λ_c approaches the classical Bondi value when r_f is 2-3 orders of magnitude larger than GM/c_sf^2 and that the kinetic energy at r_f is negligible for r_f > about 2 GM/c_sf^2. It then specializes to a Hernquist galaxy, adding the galactic gravitational potential and, in the optically thin Thomson-scattering approximation, the radiation force of the central source and the galaxy (Eqs 20-32). Parameter scans indicate that λ_c increases roughly linearly with the dimensionless galaxy linear density m_g/R_g and decreases with both Eddington-scaled luminosities l and l_g (Figs 3-6). The final section argues from Mach-number plots that the transonic critical point is a saddle.","tokens_in":13629,"tokens_out":16717,"duration_ms":171180,"significance":"If the central claims hold, the paper provides a simple, self-contained extension of Bondi's formula to finite, luminous galactic environments. The dimensionless formulation is clean, the algebra is internally consistent, and the recovery of the classical Bondi limit for large r_f is a reassuring external benchmark. The paper contains no fitted parameters and the parameter scans cover the relevant regime, so the qualitative scaling of λ_c with m_g/R_g, l, and l_g is useful for interpreting X-ray observations of elliptical galaxies. However, the dependence on l and l_g is only as robust as the inverse-square Thomson radiation force assumed in Eq. (29), and the paper's quantitative statements would benefit from a local stability analysis and numerical error estimates.","major_comments":[{"comment":"The 'exact' solution λ_c1 is not actually defined. The displayed algebraic equation involves two unknowns, λ_c1 and R, and is only the sonic-point condition u=1; it is not accompanied by the critical-point condition dg/dR=0 that would determine R_c. Please state the complete system, for example R_c from dg/dR=0 with A = -1/R_f + λ_c1^2/(2 R_f^4), together with λ_c1^{2a} = g(R_c,λ_c1)/f(1), and describe the numerical solution method. Without this, the comparison in Fig. 1 that is used to justify neglecting V_f^2/2 is not reproducible.","section":"§2, equation after Eq. (19), Fig. 1"},{"comment":"The conclusion that the transonic critical point is of saddle type is not established by the plotted Mach-number curves. Saddle, node, and center are distinguished by the eigenvalues of the linearized autonomous system at the critical point, or equivalently by the signs of the two slopes of the integral curves in the (R,u) plane. The observation that 'the slope of u's curve completely changes and becomes positive at R=R_c' is not a local stability analysis. Please either add the eigenvalue calculation or soften the claim to 'consistent with a saddle'.","section":"§5, Eq. (33), Figs 7-8"},{"comment":"The central claim that λ_c decreases monotonically with l and l_g is derived solely from the optically thin Thomson radiation force, where the luminosity enters through exact inverse-square terms and factorizes into (1-l) and (1-l_g). The manuscript lists absorption as future work in §6, but the abstract and summary present the dependence without this caveat. In a real galactic nucleus, dust absorption, line driving, or multiple scattering would modify both the magnitude and the radial dependence of f_rad, so the factorization in Eq. (32) and the linear trends in Fig. 6 are not guaranteed. Please state the optically thin, Thomson-dominated regime as a scope condition in the abstract and conclusion, and if possible give a rough column-density or dust-to-gas criterion for where the approximation breaks down.","section":"§4, Eq. (29), Fig. 6, §6"}],"minor_comments":[{"comment":"The abstract contains the typo 'kinitic' for 'kinetic', and the text has several other typographical errors ('accrection', 'centeral') that should be corrected.","section":"Abstract"},{"comment":"The caption and text say Figure 2 shows R_c, but the rendered ordinate label appears to be λ_c; please verify that every panel in Figs 2 and 3 is labeled with the quantity actually plotted.","section":"Figure 2"},{"comment":"The symbol L is used both for the central source luminosity and for the total galaxy luminosity in the same expression; please use distinct notation, for example L_BH and L_gal.","section":"Eq. (25)"},{"comment":"The statement that 'there are two solutions for u per each λ < λ_c' is imprecise, because for a fixed radius there may be zero, one, or two solutions depending on R; please rephrase as a statement about the topology of the solution curves.","section":"§5"},{"comment":"The paper does not report numerical convergence or error bars for the parameter scans; please state the root-finding tolerance, grid resolution, or other numerical accuracy used to generate the figures.","section":"Figs 1-6"},{"comment":"The critical-radius formulas are valid only for γ < 5/3; the paper treats γ = 1.4 and 1.6, but this domain of validity should be stated explicitly.","section":"Eqs. (15) and (17)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of MNRAS, and I found no indication of duplicate publication or unsupported self-citation. The main risks are the overstatement of the luminosity-dependence result beyond the optically thin regime and the incomplete specification of the 'exact' solution and the saddle-type proof; these are fixable with the revisions requested above. I therefore recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper before going deeper: it's a sound but modest extension of Bondi's spherical accretion problem to a finite outer radius in a Hernquist galaxy, including radiation pressure from the central source and the galaxy. The central derivation is internally consistent and reduces to the classical Bondi value in the right limit. If you work on SMBH accretion estimates from X-ray data, the finite-radius correction here is worth knowing, but it's not a breakthrough.\n\nThe genuinely new pieces are the general polytropic index (earlier work by Korol et al. 2016 and Ciotti & Pellegrini 2017, 2018 was isothermal), and the separate treatment of central and galactic luminosity in the radiation force. The paper also makes a clean quantitative statement: once the outer radius R_f is about 2-3 orders of magnitude larger than GM/c_sf^2, the critical accretion parameter λ_c is essentially the Bondi value. That is useful and testable.\n\nThe soft spots are in proportion to the claim. First, the comparison with Korol et al. and Ciotti & Pellegrini is qualitative. The reader is told these exist but never shown a direct numerical comparison, so it's hard to see exactly what the polytropic generality adds in practice. Second, the saddle-type conclusion in Section 5 is asserted from visual inspection of Mach-number plots. No local stability analysis is given, so that conclusion is under-supported. Third——and this is the most physically fragile piece——the radiation force is assumed to be optically thin Thomson scattering, which makes the (1-l)(1-lg) factorization in Eq. (32) exact. If absorption or multiple scattering matters, the monotonic trends in Fig. 6 are not guaranteed. The authors acknowledge this in Section 6 and defer it to future work, so it's not a hidden flaw, but it does bound the applicability of the main result.\n\nPresentation issues: the abstract in the version I have is truncated (\"e.g.:\"), and some figure captions are garbled. Minor, but worth fixing before publication.\n\nOverall, this paper deserves a serious referee. It's a coherent, self-contained derivation that extends known results in a meaningful way. A specialist in spherical accretion or galaxy-scale feedback will get value; a general reader can skip. My recommendation: send it to review, but require the authors to add a quantitative comparison with the related isothermal papers and to either prove the saddle claim or soften it.","headline":"A sound, modest extension of Bondi accretion to finite radii and Hernquist galaxy radiation; the algebra holds, but the saddle claim and radiation-force assumptions need sharper treatment.","tokens_in":14158,"tokens_out":3292,"would_cite":false,"duration_ms":33516,"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":"Galaxy gravity boosts Bondi accretion rate","keywords":["Bondi accretion","finite radius boundary","elliptical galaxies","massive black holes","Hernquist galaxy","radiation pressure","Thomson scattering","accretion rate"],"falsifier":"Re-solve the sonic-point condition retaining the boundary kinetic term $V_f^2/2$ in the Bernoulli integral for a grid of $R_f$ and $m_g/R_g$ values; if the resulting $\\lambda_c$ differs from the values produced by Eq. (32) by more than the paper's quoted few-percent accuracy for $R_f \\gtrsim 2$, or if $\\lambda_c$ fails to approach $0.625$ ($\\gamma=1.4$) and $0.367$ ($\\gamma=1.6$) as $R_f$ reaches about $10^3$, the central finite-radius claim is falsified.","tokens_in":13165,"feed_emoji":"🌌","tokens_out":16932,"duration_ms":149968,"temperature":0.7,"pith_summary":"This paper asks how spherical accretion onto a central massive black hole changes when the inflow starts not at infinity but at the outer edge of a real galactic gas reservoir, and when the gravity and light of the host elliptical galaxy also push on the gas. It establishes two things. First, in the classical Bondi problem with a finite outer boundary, the kinetic energy of the gas at that boundary is negligible once the boundary is more than about two semi-Bondi radii $GM/c_{sf}^2$, and the maximum dimensionless accretion rate $\\lambda_c$ returns to Bondi's classical values ($0.625$ for $\\gamma=1.4$, $0.367$ for $\\gamma=1.6$) once the boundary is two to three orders of magnitude beyond that radius. Second, for a Hernquist galaxy the galaxy's gravity raises $\\lambda_c$ roughly linearly in the dimensionless linear density $m_g/R_g$, while radiation from the central source and from the galaxy, parameterized by $l$ and $l_g$, lowers it. Getting these scalings right matters because real observations of inflowing gas in elliptical galaxies sit at finite radii, so using them to estimate black hole accretion rates depends on correcting for boundary and galactic radiation effects.","feed_headline":"Galaxy gravity boosts Bondi accretion rate","feed_subtitle":"Finite radii restore the classic Bondi rate; a denser galaxy raises it, radiation lowers it.","key_machinery":"The machinery is the finite-radius Bernoulli integral in dimensionless form. The paper writes the steady spherical flow equations in units where $r = R r_e$, $v = V c_{sf}$, $c_s = C_s c_{sf}$, and $\\rho = Z \\rho_f$, with $r_e = GM/c_{sf}^2$; mass conservation becomes $R^2 V Z = \\lambda$, where $\\lambda = \\dot M / (4\\pi G^2 M^2 \\rho_f/c_{sf}^3)$ is the dimensionless accretion rate. The Bernoulli equation separates into $f(u)$, unchanged from Bondi, and $g(R)$, which encodes all finite-boundary and galactic terms; the maximum of $g(R)$, at the critical radius where the Mach number is one, determines $\\lambda_c$. For the Hernquist galaxy the function is $g(R) = R^{4(\\gamma-1)/(\\gamma+1)}[(1-l)/R + m_g(1-l_g)/(R+R_g) - (1-l)/R_f - m_g(1-l_g)/(R_f+R_g) + V_f^2/2 + 1/(\\gamma-1)]$; its derivative gives the critical radius $R_c$ and its value at $R_c$ gives $\\lambda_c$ (Eqs. 18 and 32). The paper also checks the topology of the transonic solutions and finds saddle-type critical points.","core_discovery":"On its own terms, the paper claims that Bondi's steady spherical accretion solution is robust to replacing the boundary at infinity with a finite radius $r_f$, provided $r_f$ is large enough. The Bernoulli integral retained at the finite boundary shows that the boundary kinetic term $V_f^2/2$ is negligible for $R_f = r_f/(GM/c_{sf}^2) \\gtrsim 2$, so the approximate critical accretion parameter $\\lambda_{c0}$ and the exact $\\lambda_{c1}$ agree to within a few percent. At large $R_f$, $\\lambda_c$ asymptotes to the classical Bondi constants $0.625$ ($\\gamma=1.4$) and $0.367$ ($\\gamma=1.6$). In the Hernquist galaxy case, the galaxy's potential adds a term $m_g/(R+R_g)$ to the Bernoulli function, and optically thin Thomson-scattered radiation from the central source and the stellar population adds repulsive terms $(1-l)/R$ and $m_g(1-l_g)/(R+R_g)$; maximizing the resulting function $g(R)$ (Eq. 32) gives a critical radius and a critical $\\lambda_c$ that increase with the galaxy's linear density $m_g/R_g$ and decrease with both luminosity parameters, with a saddle-type transonic solution at the critical point.","pith_inferences":["The near-linear dependence of $\\lambda_c$ on $m_g/R_g$ hints at a ready correction factor of the form $1 + \\alpha\\, m_g/R_g$ for Bondi-based black hole accretion estimates in ellipticals, with $\\alpha$ set by $\\gamma$, $l$, and $l_g$; fitting that factor to observed galaxies would be a direct test of the model.","Because the radiation terms appear as $(1-l)$ and $m_g(1-l_g)$ multiplying the two gravitational terms, the same equations imply an Eddington-scaled threshold near $l=1$ or $l_g=1$ at which the effective gravity reverses and outflow rather than accretion should set in, a regime the paper does not explore.","The optically thin, Thomson-only treatment is best suited to hot, diffuse coronal gas; at smaller radii where absorption and line driving matter, the derived scalings would need to be replaced, so for bright systems the present $\\lambda_c$ likely acts as an upper bound rather than a full prediction."],"forward_implications":["For observed elliptical galaxies, Bondi estimates of black hole accretion rates made from gas density and temperature at a finite radius remain reliable whenever that radius is more than about $10^2$-$10^3 GM_{\\mathrm{BH}}/c_{sf}^2$; closer in, boundary corrections should be applied.","In a Hernquist host, a denser, more compact stellar distribution (larger $m_g/R_g$) raises the maximum possible accretion rate, so galaxy mass profiles push the inferred rate above the classical Bondi value.","Radiative feedback suppresses accretion: higher Eddington-scaled central luminosity $l$ and galactic luminosity $l_g$ both decrease $\\lambda_c$, with the central luminosity having the stronger effect.","The critical point is always a saddle, so the analytic maximum-accretion solution is the physically selected transonic solution and can serve as the reference state for time-dependent or multi-dimensional work."],"supporting_citations":[{"why":"Defines the classical spherical accretion model and the critical accretion rate that the paper generalizes to a finite outer boundary.","marker":"Bondi (1952)"},{"why":"Supplies the galaxy potential and the deprojected luminosity profile used for the galactic gravity and radiation terms.","marker":"Hernquist (1990)"},{"why":"Provides the constant-luminosity assumption that lets the central radiation force be written as a simple inverse-square term.","marker":"Cassinelli & Castor (1973)"},{"why":"Provides the luminosity-density formula j(r) that the paper integrates to obtain the galaxy's luminosity L_g(r).","marker":"Dehnen (1993)"},{"why":"Together with Dehnen, gives the three-dimensional radiation energy density model for the galaxy's diffuse light.","marker":"Tremaine et al. (1994)"},{"why":"Earlier finite-radius Bondi application to early-type galaxies whose deviations from Bondi values this work extends with galactic gravity and radiation.","marker":"Korol et al. (2016)"},{"why":"Establishes the analytical isothermal accretion framework in Hernquist and Jaffe galaxies against which the present critical parameters are derived.","marker":"Ciotti & Pellegrini (2017)"},{"why":"Shows that radiation from the central source reduces the spherical accretion rate, the effect the present luminosity parameters quantify.","marker":"Fukue (2001)"},{"why":"Supplies the classification of critical points used to identify the transonic solution as saddle-type.","marker":"Kato et al. (2008)"}],"fun_headline_variants":["Finite radius restores classic Bondi accretion","Galaxy gravity boosts Bondi accretion rate","Radiation damps Bondi accretion in ellipticals","Bondi accretion robust to finite boundary","Hernquist galaxy gravity raises accretion limit"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that every gas parcel is optically thin with opacity set only by Thomson scattering and that both the central black-hole luminosity and the galaxy luminosity are constant with radius, so the radiation force acts as pure inverse-square repulsion; if absorption, multiple scattering, or radiation-driven outflows become significant, the derived dependence of $\\lambda_c$ on $m_g/R_g$, $l$, and $l_g$ no longer follows.","fun_headline_variants_meta":{"raw":{"variants":["Finite radius restores classic Bondi accretion","Galaxy gravity boosts Bondi accretion rate","Radiation damps Bondi accretion in ellipticals","Bondi accretion robust to finite boundary","Hernquist galaxy gravity raises accretion limit"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000565,"raw_usage":{"total_tokens":2725,"prompt_tokens":1040,"completion_tokens":1685,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":656,"completion_tokens_details":{"reasoning_tokens":1616}},"tokens_in":656,"tokens_out":1685,"duration_ms":12264,"temperature":1.0,"reasoning_tokens":1616,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:35:41.938806+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-solve the sonic-point condition retaining the boundary kinetic term $V_f^2/2$ in the Bernoulli integral for a grid of $R_f$ and $m_g/R_g$ values; if the resulting $\\lambda_c$ differs from the values produced by Eq. (32) by more than the paper's quoted few-percent accuracy for $R_f \\gtrsim 2$, or if $\\lambda_c$ fails to approach $0.625$ ($\\gamma=1.4$) and $0.367$ ($\\gamma=1.6$) as $R_f$ reaches about $10^3$, the central finite-radius claim is falsified.","supporting_citations":[{"cited_title":"Kyoto Univ","cited_arxiv_id":null,"evidence_quote":"Supplies the galaxy potential and the deprojected luminosity profile used for the galactic gravity and radiation terms."}],"review_version":1}