{"id":"dacbdd62-70a1-42b3-a752-fdf6ad019820","arxiv_id":"1909.01527","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"An equatorial overdense, polar underdense accretion flow onto a Schwarzschild black hole leads to a black-hole accretion rate of about the Michel value, with the excess mass ejected as a bipolar outflow.","lead":"This paper shows that gas falling onto a black hole can choke at a nearly fixed accretion rate, with the excess mass pushed out as two opposite jets by fluid pressure alone. An exact relativistic calculation for a special fluid and computer simulations for ordinary gases both produce this inflow-outflow pattern.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Numerical validation of choked accretion is partly circular: Fig. 14 checks mass conservation, and the analytic outflow is encoded in the l=2 boundary data rather than dynamically emergent.","rationale":"The analytic model is exact under its stated assumptions, so the concern is not an internal inconsistency. However, the paper's central claim that \"any injected flux above that rate is ejected as a bipolar outflow\" is supported by a construction in which the outflow is part of the assumed l=2 velocity potential, rather than an emergent response to overfeeding. The numerical experiments are meant to supply that emergence, but they impose an anisotropic density profile at an open boundary with free velocities, so the measured injection and ejection rates are outputs of the simulation, not controlled inputs. The comparison in Figure 14 uses Eq. (2.52), which is a restatement of steady mass conservation; agreement with it is therefore not evidence for the choked-accretion mechanism. This does not overturn the paper: the analytical solution remains a valid toy model, the benchmark test in Section 3.1 is genuine, and the polytropic runs do show a bipolar morphology for the chosen boundary data. But it shifts the evidential weight from the numerical validation back to the analytic construction, and it reinforces the already conditional astrophysical scope. The reader flagged the open-boundary issue and the lack of resolution studies; the additional point here is that Figure 14 cannot discriminate the model because it is tautological. A boundary-placement and injection-prescription test would settle whether the polytropic outflow is a robust steady-state feature or an artifact of the outer boundary treatment.","tokens_in":27445,"tokens_out":7581,"duration_ms":83630,"concrete_test":"Re-run the gamma=5/3, R=100M, a0=0.6, delta=0.5% polytropic case (Table 5) with two modifications: (1) move the outer boundary to R=200M with the same inner setup; (2) at R=100M replace the free velocity boundary with a pure radial injection of the same total Mdot_in, keeping the equatorial-to-polar density ratio. If the bipolar outflow fraction or Mdot/Mdot_M changes by more than about 10% between these runs, the steady state is controlled by the outer-boundary prescription rather than by choking at the horizon. Additionally, replace the analytic curve in Figure 14 with the identity Mdot_ej/Mdot_in = 1 - Mdot/Mdot_in; the points must fall on it by construction, so only deviations, not agreement, carry information.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central analytic construction is internally sound: Eq. (2.42) follows from spherical-harmonic orthogonality, and Eqs. (2.49)-(2.50) describe a family of steady potential flows in which Mdot_in - Mdot_ej = Mdot. The load-bearing weakness is that this family does not demonstrate that an independently specified excess flux is \"choked\" and redirected. In Eq. (2.13) the l=2 term is part of the assumed velocity potential, and the outer boundary data include the polar outflow through Eq. (2.33); the outflow is prescribed by the boundary condition, not produced by a dynamical bottleneck. The numerical polytropic simulations are intended to relax this, but they impose the density contrast at R and leave both velocity components free (Section 3.2), so Mdot_in is diagnosed, not imposed. Moreover, the headline comparison in Figure 14 uses Eq. (2.52), which is algebraically identical to Mdot_ej/Mdot_in = 1 - Mdot/Mdot_in. This is just steady-state mass conservation; any simulation that conserves mass will land on the \"analytic\" curve regardless of the choking mechanism. The only nontrivial numerical content is that Mdot stays within a factor ~1.7 of Mdot_M (Tables 5-6) and that a bipolar morphology appears for the chosen boundary data. The paper's own Section 4.2 further restricts the astrophysical setting to pair-dominated plasmas with T > 10^11 K, so the mechanism's reach is narrow even if the simulations are accepted.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes a hydrodynamical 'choked accretion' mechanism for launching bipolar outflows from a Schwarzschild black hole. The analytic part constructs a stationary, axisymmetric, irrotational perfect-fluid solution with a stiff equation of state, taking a velocity potential composed of a monopole plus an l=2 quadrupole term (Eq. 2.13). The solution gives a horizon accretion rate Mdot = 16πM^2 α0 ρ0 Γ0 (Eq. 2.42) that is independent of the quadrupole amplitude, and a family of inflow-outflow configurations when the equatorial inflow speed V0 exceeds 6M^2/R^2 (Eq. 2.35). The authors define injection and ejection rates (Eqs. 2.49–2.50) and interpret the balance as a choked-accretion relation. They then perform axisymmetric general-relativistic hydrodynamic simulations with the aztekas code for polytropic equations of state, imposing an equatorial density contrast at the injection sphere while leaving the velocity components free. The simulations produce bipolar outflows with accretion rates of order the Michel rate, and the paper compares several flow diagnostics to the analytic model. Astrophysical applicability is discussed, with the authors noting that standard disk temperatures imply stagnation radii much larger than the inner disk and that the mechanism becomes relevant mainly for electron-positron plasmas at T > 10^11 K (Section 4.2).","tokens_in":27782,"tokens_out":8378,"duration_ms":80876,"significance":"If the mechanism operates as described, it offers a purely hydrodynamical route to bipolar outflows that does not rely on magnetic fields or black-hole spin, complementing the Blandford-Payne and Blandford-Znajek mechanisms. The exact analytic solution is a useful benchmark for relativistic hydrodynamics codes; the paper demonstrates convergence to 0.001% in the stiff-fluid benchmark (Fig. 6) and reports open-source code validation through shock-tube tests (Appendix B). The numerical exploration spans several polytropic indices, density contrasts, and sound speeds, and the qualitative inflow-outflow morphology appears robust across the parameter scan. The main limitations are the restrictive stiff-fluid assumption in the analytic model and the acknowledged absence of rotation, magnetic fields, and radiation transport; the paper's own Section 4.2 narrows the astrophysical applicability to pair-dominated plasmas with T > 10^11 K. These caveats are stated openly, which is a strength, but they should be weighed in the verdict.","major_comments":[{"comment":"The comparison in Fig. 14 is tautological and does not validate the choked-accretion model. Eq. (2.52) is simply Mdot_ej/Mdot_in = 1 − Mdot/Mdot_in, which follows directly from steady-state mass conservation and the definitions of the three rates. Any mass-conserving simulation will land on the analytic curve regardless of whether a choking mechanism is present. The 'very good agreement' claimed in Section 4.1 is therefore not evidence for the model. I recommend replacing Fig. 14 with a nontrivial comparison, such as the measured Mdot versus the Michel rate, or the predicted stagnation-point location as in Fig. 15, which is a genuine test.","section":"Section 4.1, Eq. (2.52), Fig. 14"},{"comment":"The analytic model prescribes the outflow at the boundary rather than showing that an excess flux is dynamically redirected. The l=2 term in Eq. (2.13) is part of the assumed velocity potential, and the polar ejection velocity Vej in Eq. (2.33) is fixed by the same parameter V0 that sets the equatorial inflow. Eqs. (2.49)–(2.50) therefore describe a kinematic decomposition of a given steady-flow family, not an emergent choking of an independently specified injection rate. This does not invalidate the construction as an exact toy model, but the paper should be explicit that the 'choking' claim rests on the numerical simulations of Section 3, where the density contrast is imposed and the velocities are free. Alternatively, the authors could formulate a boundary-value problem in which Mdot_in is fixed and show that the steady solution has the outflow determined by the equations.","section":"Section 2.2, Eqs. (2.13), (2.33), (2.49)–(2.50)"},{"comment":"The numerical evidence for flux saturation is stronger than the paper's presentation suggests, but it is not quantified. For fixed a0, the measured Mdot stays nearly constant while Mdot_in increases by up to a factor of six (e.g., Table 5, a0=0.6: Mdot_in from 52.6 to 318.9 while Mdot remains about 10.8), which is the key 'choked' behavior. The paper should present this explicitly, for example as a plot of Mdot/Mdot_M versus Mdot_in/Mdot_M, and give a quantitative range for the ratio Mdot/Mdot_M (0.98–1.67 across Tables 3–6) rather than stating only that it is 'of the order' of the Michel rate.","section":"Section 3.2, Tables 3–6"}],"minor_comments":[{"comment":"The displayed equation has a typesetting issue: '8πρ 0' should be '8πρ0', and the square-root expression should be checked for readability.","section":"Section 2.5, Eq. (2.51)"},{"comment":"The notation Mdot_M for the Michel rate is easily confused with the derivative of the black-hole mass; consider renaming it to, for example, Mdot_Michel for clarity.","section":"Appendix A, Eq. (A.15)"},{"comment":"The paper's own discussion restricts the astrophysical applicability: for ordinary hydrogen plasmas the stagnation radius S is estimated to exceed 10^3–10^5 M, making the mechanism relevant only for pair-dominated plasmas at T > 10^11 K. This scope restriction should be stated more prominently, ideally in the abstract, so that readers do not overgeneralize the jet-launching claim.","section":"Section 4.2"},{"comment":"The sentence introducing Eq. (3.5) says that the density contrast δ is 'the same' as in Eq. (2.38), which is true at the pole, but the functional forms differ; a brief note that Eq. (3.5) is a first-order parametrization qualitatively matching the analytic profile would remove potential confusion.","section":"Section 3.2, Eq. (3.5)"}],"recommendation":"major_revision","confidential_remarks":"The analytic solution is elegant and the numerical work appears competently executed, but the main validation figure is tautological and the analytic model does not by itself demonstrate an emergent choking mechanism. The paper is likely publishable after reframing the analytic contribution as an exact illustrative family of solutions, replacing Fig. 14 with a nontrivial test, and quantifying the numerical saturation behavior. I therefore recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this paper. First, it contains an exact, closed-form general-relativistic potential-flow solution for a Schwarzschild black hole accreting with an l=2 equatorial density enhancement. That is real and new relative to the Petrich-Shapiro-Teukolsky monopole/dipole work and to the authors' own Newtonian papers. Second, the \"choked accretion\" claim is more modest than the title suggests: the outflow is baked into the chosen boundary data, and the headline numerical comparison (Figure 14) is essentially steady-state mass conservation, not a direct test of a choking mechanism. Worth knowing before you read further.\n\nCredit where due: the analytic part is clean and self-contained. The monopole fixes Mdot = 16 pi M^2 alpha0 rho0 Gamma0 by spherical-harmonic orthogonality; the quadrupole produces the stagnation circle and polar outflow once V0 exceeds 6M^2/R^2. The stiff-fluid benchmark reproduces the analytic accretion rate to 0.001% (Figure 6), and the code validation in the appendices is solid. The polytropic simulations are exploratory but competent: they show a bipolar inflow-outflow morphology even for a 0.1% density contrast, and Mdot stays within a factor ~1.7 of the Michel rate across the parameter space. The paper is honest about the narrow astrophysical window: Section 4.2 explicitly restricts the mechanism to pair-dominated plasmas with T > 10^11 K, and it calls out rotation, magnetic fields, and radiative transport as missing physics.\n\nSoft spots, in rough order of importance. First, the stress-test note is right about Figure 14. Equation (2.52) is just Mej/Min = 1 - Mdot/Min, which any steady, mass-conserving simulation will satisfy. The plot confirms conservation, not the choking mechanism. The real numerical content is that Mdot stays near Mdot_M and the morphology is bipolar. The paper's phrase \"remarkable agreement\" oversells what that figure shows. Second, the polytropic runs lack resolution studies and error bars. The stiff-fluid benchmark has a resolution study; the polytropic tables do not. One or two resolution checks would be cheap insurance. Third, the outer boundary is open and the \"injection rate\" is diagnosed, not imposed. The paper acknowledges this, and it is acceptable for a first exploration, but it means the parameter space is not fully under control. None of this sinks the paper. The exact solution is a useful benchmark, and the mechanism is plausible as a toy model. The overclaiming is mostly in the title and in the reach of Section 4, not in the math.\n\nWho this is for: people working on hydrodynamical jet-launching, Bondi/Michel accretion, and code validation. If I were refereeing, I would ask for a modest revision that rephrases the Figure 14 claim, adds a resolution check for the polytropic runs, and tempers the astrophysical discussion. That is a standard referee job, not a rejection. Send it to review.","headline":"A genuine exact l=2 relativistic accretion solution worth refereeing, but the 'choked accretion' mechanism is largely encoded in the boundary data and the headline numerical comparison is just mass conservation.","tokens_in":28344,"tokens_out":3625,"would_cite":true,"duration_ms":38374,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper shows that a radially infalling flow onto a Schwarzschild black hole, when slightly denser at the equator than at the poles, accretes at a fixed Bondi–Michel rate and redirects all excess matter into a bipolar outflow.","keywords":["choked accretion","Bondi accretion","Michel accretion","black hole physics","relativistic hydrodynamics","astrophysical jets","stiff fluid","Schwarzschild spacetime"],"falsifier":"Run a three-dimensional general-relativistic hydrodynamic simulation in which a rotating Keplerian disk feeds the inner boundary with a realistic angular-momentum profile; if the black hole accretes substantially faster than the Michel rate without producing a persistent polar outflow, or if the outflow disappears when the equatorial density contrast is removed with everything else unchanged, the choked-accretion claim would be falsified.","tokens_in":27216,"feed_emoji":"🕳️","tokens_out":8127,"duration_ms":78805,"temperature":0.7,"pith_summary":"The paper establishes a purely hydrodynamical jet-launching mechanism: when matter falls radially onto a Schwarzschild black hole with an equatorial overdensity and underdense poles, the flow cannot accrete faster than a fixed rate close to the Bondi–Michel value. Any excess flux is choked at a gravitational bottleneck and redirected along the polar axis as an outflow. The authors prove this exactly for an ultrarelativistic stiff fluid and show numerically that the same inflow-outflow morphology appears for general polytropic gases, with the accretion rate staying near the Michel value. If correct, the mechanism offers a way to launch jets without invoking magnetic fields or a spinning black hole.","feed_headline":"Choked accretion turns excess infall into a bipolar outflow","feed_subtitle":"Even a 0.1% equator-to-pole density contrast can flip accretion into a bipolar jet.","key_machinery":"The central object is the quadrupolar velocity potential $\\Phi=-e\\left[t+2M\\ln(1-2M/r)-A(3r^2-6Mr+2M^2)(3\\cos^2\\theta-1)\\right]$, built from the general solution of the relativistic potential-flow wave equation in Schwarzschild spacetime. Regularity of the solution across the event horizon fixes the monopole term and thereby fixes the total accretion rate, while the quadrupole term redirects infalling streamlines toward the poles. The stagnation-point radius $S$ marks the dividing streamline: streamlines with $|\\Psi|<1$ accrete and those with $|\\Psi|>1$ escape, with ejection velocity and equatorial-to-polar density contrast related by explicit algebraic formulas.","core_discovery":"For a stationary, axisymmetric, irrotational flow of a stiff fluid with equation of state $P=K\\rho^2$, the paper constructs an exact solution around a Schwarzschild black hole using a quadrupolar $l=2$ perturbation of the velocity potential. The black hole then accretes at the fixed rate $\\dot M=16\\pi M^2\\alpha_0\\rho_0\\Gamma_0$, which is close to the Michel rate, and the solution contains polar stagnation points that divide the flow into streamlines that fall into the hole and streamlines that escape as a bipolar outflow. The ejection rate satisfies $\\dot M_{\\rm ej}/\\dot M_{\\rm in}=1-1/\\Lambda$, where $\\Lambda$ grows with the injection rate, so the excess injected mass is expelled rather than accreted. Numerical simulations with a polytropic equation of state and a boundary density profile $\\rho(\\theta)=\\rho_0(1-\\delta\\cos^2\\theta)$ recover the same stable inflow-outflow configuration even for a density contrast as small as $\\delta=0.1\\%$, with the accretion rate remaining of order the Michel value.","pith_inferences":["Because the choked rate is essentially the Michel rate, the model implies that time variability in the external mass supply maps almost one-to-one onto outflow variability, so stochastic feeding should produce flickering in the jet's mass flux before any magnetic collimation acts.","The same geometric idea could apply to other compact objects, where the horizon regularity condition would be replaced by a surface boundary condition; the paper confines its claims to black holes.","A full three-dimensional simulation with a realistic rotating disk would test whether the disk itself creates the required equatorial overdensity, turning the injection radius from an imposed boundary into an emergent property.","The relation $\\dot M_{\\rm ej}/\\dot M_{\\rm in}=1-1/\\Lambda$ could be used as an observational diagnostic: measuring the jet-to-accretion mass ratio and the injection rate would locate the stagnation radius and effective density contrast in a real source."],"forward_implications":["If the mass supply crossing the injection sphere exceeds the Bondi–Michel rate, the steady state is not higher accretion but a bipolar outflow, with the black hole still accreting at order $\\dot M_{\\rm M}$.","A very small equator-to-pole density contrast, of order $0.1\\%$, suffices to switch the flow from spherical inflow to an inflow-outflow morphology.","The ejected fraction and ejection velocity increase monotonically with the injection rate and density contrast; in the analytic model the ejection velocity approaches the speed of light as the polar density contrast approaches unity.","The mechanism provides a jet-launching route that does not require magnetic fields or a rotating black hole, so it may operate in systems where magnetic or Kerr-related mechanisms are weak.","For ordinary plasmas the mechanism requires high inflow temperatures, and the paper estimates it is most applicable to relativistic electron-positron pair plasmas with temperatures above about $10^{11}$ K."],"supporting_citations":[{"why":"Supplies the exact general solution of the potential-flow wave equation in Schwarzschild spacetime used to build the analytic choked-accretion model; its horizon-regularity condition fixes the monopole accretion rate.","marker":"Petrich et al. 1988"},{"why":"Defines the spherically symmetric relativistic accretion solution whose mass accretion rate sets the threshold at which the flow chokes in both the analytic and numerical models.","marker":"Michel 1972"},{"why":"Provides the non-relativistic spherical accretion rate that the paper identifies as the Newtonian counterpart of the limiting accretion value.","marker":"Bondi 1952"},{"why":"Proposed the original purely hydrodynamical jet-launching mechanism based on an axisymmetric polar density gradient, which this work revisits in the relativistic regime.","marker":"Hernandez et al. 2014"},{"why":"Presents the non-relativistic counterpart of the analytic model and supplies Newtonian numerical results that the relativistic simulations are compared with.","marker":"Aguayo-Ortiz et al. 2019"},{"why":"Gives a general relativistic treatment of Michel-type accretion for general equations of state, used to identify the stiff-fluid Michel accretion rate.","marker":"Chaverra & Sarbach 2015"},{"why":"Describes the numerical scheme and discretization underpinning the polytropic simulations that extend the analytic result.","marker":"Tejeda & Aguayo-Ortiz 2019"}],"fun_headline_variants":["Slight equator bias flips black hole accretion into jet","Even 0.1% density skew makes black hole eject bipolar flow","Choked inflow: tiny asymmetry turns black hole into jet","Bipolar jet from black hole triggered by tiny density contrast"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The mechanism's astrophysical relevance rests on the assumption that gas arriving at the inner edge of the disk has already lost essentially all its angular momentum, so the inflow is radial and irrotational, and that an equator-to-pole density contrast is established and maintained at the injection boundary.","fun_headline_variants_meta":{"raw":{"variants":["Slight equator bias flips black hole accretion into jet","Even 0.1% density skew makes black hole eject bipolar flow","Choked inflow: tiny asymmetry turns black hole into jet","Bipolar jet from black hole triggered by tiny density contrast"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00071,"raw_usage":{"total_tokens":3252,"prompt_tokens":1057,"completion_tokens":2195,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":673,"completion_tokens_details":{"reasoning_tokens":2125}},"tokens_in":673,"tokens_out":2195,"duration_ms":13931,"temperature":1.0,"reasoning_tokens":2125,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:15:48.273975+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a three-dimensional general-relativistic hydrodynamic simulation in which a rotating Keplerian disk feeds the inner boundary with a realistic angular-momentum profile; if the black hole accretes substantially faster than the Michel rate without producing a persistent polar outflow, or if the outflow disappears when the equatorial density contrast is removed with everything else unchanged, the choked-accretion claim would be falsified.","supporting_citations":[{"cited_title":"I., Shapiro S","cited_arxiv_id":null,"evidence_quote":"Supplies the exact general solution of the potential-flow wave equation in Schwarzschild spacetime used to build the analytic choked-accretion model; its horizon-regularity condition fixes the monopole accretion rate."},{"cited_title":"C., 1972, Astrophysics and Space Science, 15, 153","cited_arxiv_id":null,"evidence_quote":"Defines the spherically symmetric relativistic accretion solution whose mass accretion rate sets the threshold at which the flow chokes in both the analytic and numerical models."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the non-relativistic spherical accretion rate that the paper identifies as the Newtonian counterpart of the limiting accretion value."},{"cited_title":"L., Rodr \\' guez-Mota R","cited_arxiv_id":null,"evidence_quote":"Proposed the original purely hydrodynamical jet-launching mechanism based on an axisymmetric polar density gradient, which this work revisits in the relativistic regime."},{"cited_title":"Choked accretion: from radial infall to bipolar outflows by breaking spherical symmetry","cited_arxiv_id":"1909.00884","evidence_quote":"Presents the non-relativistic counterpart of the analytic model and supplies Newtonian numerical results that the relativistic simulations are compared with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives a general relativistic treatment of Michel-type accretion for general equations of state, used to identify the stiff-fluid Michel accretion rate."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Describes the numerical scheme and discretization underpinning the polytropic simulations that extend the analytic result."}],"review_version":1}