{"id":"fb56ae54-48c7-4fe6-9455-537e60555df3","arxiv_id":"1909.00884","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Hydrodynamic simulations show that small axisymmetric density anisotropies cap the accretion rate onto a point mass at the Bondi value and eject the surplus in bipolar outflows.","lead":"This paper shows that when gas falling onto a star or compact object is a little denser at the equator than at the poles, the flow spontaneously splits into equatorial infall and bipolar outflow. The central object's accretion rate stays pinned near the classic Bondi value, so all excess mass is ejected instead of swallowed.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed Bondi-limited plateau is not established as a boundary-independent maximum: with free-velocity outer boundaries the inflow rate self-selects, and the paper reports that larger density contrasts or smaller outer radii fail to reach steady state.","rationale":"Reading in good faith, the paper is careful: it presents the choked-accretion phenomenon as suggested by numerical experiments, reproduces the Bondi solution in the spherical limit, performs a self-convergence test, and explicitly reports the unstable regime and sensitivity to boundary conditions. Those are real supporting elements. Nevertheless, the headline physical content, namely the existence of a universal maximum accretion rate, requires that an excess supply be diverted, and the numerical setup never independently sets the supply. The free-velocity outer boundary lets the flow select the accretion rate that matches a steady solution, while the instability reported for larger delta or smaller R occurs exactly where a genuine super-Bondi supply would have to be rejected. This does not refute the mechanism, but it means the current evidence supports the existence of a particular steady inflow/outflow branch rather than a boundary-independent maximum. The reader's conditional verdict is therefore appropriate; the proposed fixed-flux experiment is the decisive next step.","tokens_in":13260,"tokens_out":8821,"duration_ms":102923,"concrete_test":"Run the gamma = 4/3, delta = 1% configuration with the outer boundary modified so that the total inward mass flux through r = 10 rB is prescribed to be 5 Mdot_B rather than diagnosed: keep rho(theta) as in Eq. (3.2) and impose v_r(theta) = -[5 Mdot_B / (4 pi R^2 rho(theta))] on inflow latitudes, with a characteristic-compatible treatment of the other variables. If the time-averaged accretion rate through Racc = 0.1 rB exceeds Mdot_B by more than 10%, or if no steady state is reached, the choked plateau is not a robust maximum. A companion run with Racc = 0.05 rB would test sink-radius independence.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the inner accretion rate chokes at approximately the Bondi value regardless of the mass supply, with all excess material ejected. The simulations do not actually force a super-Bondi supply. In Section 3.1 the outer boundary at R = 10 rB fixes rho(theta) = rho0 (1 - delta cos^2 theta) and P = rho^gamma / gamma, but the authors 'allow free (as free-outflow condition) evolution on both velocity components.' The radial velocity at R, and hence the diagnosed injection rate Mdot_in, is part of the steady solution rather than an independently imposed control parameter. The trend in Figure 8, with Mdot_in/Mdot_B up to roughly 17, is produced by changing delta, which also changes the pressure gradient that drives the outflow; this is not a controlled experiment in which an external supply rate is varied while the symmetry breaking is held fixed. The paper itself states in Section 4 that 'taking larger values of delta or outer boundaries smaller than rB lead to rather unstable, highly dynamic accretion flows that do not seem to relax to steady state configurations,' and that 'very distinct boundary conditions may very well lead to substantially different solutions.' These are precisely the conditions that would supply more mass and test the claimed maximum. In addition, the free-outflow inner boundary at 0.1 rB removes whatever reaches it, so the sink itself cannot enforce a limiting rate. The claimed universal maximum therefore rests on a single, self-selected steady branch.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper investigates the effect of a small axisymmetric density contrast between the equatorial plane and the poles on steady-state hydrodynamic accretion onto a Newtonian point mass. The authors first present an analytic potential-flow toy model for an incompressible fluid (Section 2) that produces an inflow/outflow geometry with polar stagnation points, and then perform axisymmetric hydrodynamic simulations for adiabatic indices γ = 1, 4/3, and 7/5 using the aztekas code (Section 3). The outer boundary imposes a density profile ρ(θ)=ρ0(1−δ cos²θ) and the corresponding pressure, while velocities evolve freely; the inner boundary is a free-outflow sphere at 0.1 rB. The simulations are validated against the Bondi solution to within 1% and show second-order convergence. For δ = 0.1%, 1%, and 10%, the authors find steady-state flows in which the polar regions become underdense and drive a bipolar outflow, while the mass accretion rate through the inner boundary saturates at approximately the Bondi rate, with the ratio Mdot/Mdot_B within 5% of unity. The injection rate Mdot_in/Mdot_B increases with δ up to ~17, and the ratio of ejected to injected mass approaches ~94% at δ = 10%. The paper interprets this as evidence for a 'choked accretion' mechanism that limits the accretion rate to the Bondi value and ejects all excess mass.","tokens_in":13488,"tokens_out":7747,"duration_ms":79858,"significance":"If the choking phenomenon were shown to be independent of the specific boundary setup, it would constitute a new, purely hydrodynamical mechanism for producing bipolar outflows from quasi-spherical accretion, complementing magneto-rotational jet models. The paper has notable strengths: the numerical code is benchmarked against the analytic Bondi solution to better than 1% for three equations of state, and a self-convergence test demonstrates second-order convergence. The analytic toy model transparently illustrates the inflow/outflow geometry, and the numerical results show a consistent trend across γ and δ. However, the central claim of a maximum accretion rate is only established within a narrow parameter window (δ ≤ 10%, R = 10 rB, steady flow); the paper itself notes that larger contrasts or smaller outer radii produce non-steady solutions. The generality of the 'choked accretion' mechanism therefore remains a hypothesis, not a demonstrated theorem.","major_comments":[{"comment":"The claim that the inner accretion rate chokes at the Bondi value 'regardless of how large the injected mass rate is' is not supported by the present simulations, which cover only the steady-state regime with density contrasts δ ≤ 10% and an outer boundary at R = 10 rB (Section 3.1). The paper explicitly reports that larger δ or outer boundaries smaller than rB lead to 'rather unstable, highly dynamic accretion flows that do not seem to relax to steady state configurations.' Since the choking behavior is demonstrated only for the steady subset, the universal statement in the abstract ('any extra material being ejected') is not established. The authors should either extend the analysis into the non-steady regime (e.g., by time-averaging the accretion rate and outflow) or explicitly restrict the conclusion to the parameter range explored.","section":"Section 4 (and abstract)"},{"comment":"The injection rate Mdot_in is not an independent control parameter in these experiments. At the outer boundary the authors impose only the density profile (3.2) and the corresponding pressure, while allowing free evolution of both velocity components; consequently Mdot_in and Mdot_ej are emergent properties of the steady solution, not externally imposed mass supplies. Figure 8 thus demonstrates a correlation between Mdot_in/Mdot_B and Mdot_ej/Mdot_in as δ is varied, but it does not test the physical prediction that a prescribed super-Bondi mass supply would be diverted into the outflow while the inner rate stays fixed at Mdot_B. A test that controls the supply rate independently (e.g., by imposing a radial velocity profile at R) is needed to validate the 'maximum achievable accretion rate' interpretation.","section":"Section 3.1, Figure 8"}],"minor_comments":[{"comment":"The reported rates are given to two decimal places without error bars; since the key claim is that Mdot/Mdot_B is close to unity (within a few percent), the numerical uncertainty from the convergence study should be quantified so that the reader can assess the significance of the deviations.","section":"Table 1"},{"comment":"The analytic model sets the accretion rate to the Bondi value by hand in Eq. (2.7); the paper should state more explicitly that this is an assumed input for the toy model rather than a derived consequence, to avoid any appearance of circularity in the analytic motivation.","section":"Section 2, Eq. (2.7)"},{"comment":"The phrase 'free (as free-outflow condition) evolution' for the velocity components at the outer boundary is ambiguous; a zero-gradient copied ghost-cell condition permits both inflow and outflow, so the authors should specify the exact implementation and verify that it does not artificially constrain the diagnosed injection rate.","section":"Section 3.1"},{"comment":"Equation (2.5) is introduced in the context of the incompressible analytic model; the numerical computation of Mdot should instead be stated directly as the surface integral of ρ v_r over a sphere, to avoid confusion about which equations are used in the simulations.","section":"Section 3.3"},{"comment":"There are several typographical artifacts, including a stray backslash character in 'i.e.⃗ v= ∇Φ' in Section 2 and the appearance of '/Slash1s' in the axis label of Figure 1; these should be corrected in the final version.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is well written and the numerical results appear credible within the stated parameter range. The main concern is the generality of the choking claim: the paper's own caveats in Section 4 should be reflected in the abstract and conclusions so that the 'maximum accretion rate' is presented as a hypothesis supported by a limited set of steady-state solutions, rather than as an established universal result."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nRead Aguayo-Ortiz et al. The paper is worth your time: it shows numerically that a steady, non-rotating, axisymmetric accretion flow with a mild equatorial over-density develops a bipolar outflow while the net accretion onto the central object stays at roughly the Bondi rate. The choking effect—inner Mdot saturating near Mdot_B while excess supply is ejected—is new relative to the earlier isothermal perturbation work. The simulations are careful: Bondi validation at sub-percent level, second-order convergence, three adiabatic indices, and a consistent morphology across the board.\n\nWhat I find genuinely useful: the idea that spherical symmetry breaking alone can produce an in/outflow geometry, without rotation or magnetic fields, and that the inner flow converges to a quasi-Bondi configuration. The authors are also honest about where the model strains, and they don't oversell the astrophysical applicability—they explicitly note that S is large and that radiation/MHD effects would shrink it.\n\nThe soft spots are real but not fatal. The strongest claim—that there is a maximum achievable accretion rate equal to the Bondi value regardless of supply—is not as controlled as it sounds. The outer boundary fixes density and pressure but lets velocity float, so the injection rate Mdot_in is part of the solution, not an externally imposed control parameter. Varying delta changes both the supply and the driving pressure gradient simultaneously. The authors themselves report that larger density contrasts or outer boundaries smaller than rB go unstable and never reach steady state, and they concede that very different boundary conditions may produce very different solutions. So the \"universal\" part is an extrapolation from a single family of self-selected steady branches. The free-outflow inner boundary also means the sink itself imposes no limiting rate; the plateau is a property of the flow, which is interesting, but the parameter space tested is narrow (delta up to 10%, Mdot_in up to about 17 Mdot_B). Table 1 has no error bars, though the convergence test mitigates that.\n\nThe analytic model in Section 2 is explicitly a toy: alpha is set to the Bondi rate by hand, so it doesn't independently support the choking claim. That's fine as motivation, and the simulations are the real evidence.\n\nNet: this is a solid, well-executed numerical study of a plausible new mechanism. It deserves a serious referee and a thoughtful revision that either narrows the claims or expands the boundary-condition study. I'd bring it to group.","headline":"A clean numerical demonstration of a new accretion-flow morphology, with the strong universal choking claim outrunning the parameter space actually explored.","tokens_in":14072,"tokens_out":2453,"would_cite":true,"duration_ms":23353,"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":"A small equatorial-polar density imbalance makes steady hydrodynamic accretion onto a point mass saturate at the Bondi rate and eject all surplus matter in bipolar outflows.","keywords":["accretion","Bondi accretion","bipolar outflows","axisymmetric hydrodynamics","spherical symmetry breaking","astrophysical jets","numerical simulations"],"falsifier":"Run the same hydrodynamic problem with the inner boundary changed from free outflow to a hard absorbing surface, or with the outer injection sphere moved inward from $10\\,r_\\mathrm{B}$ to a few $r_\\mathrm{B}$, and measure the steady-state accretion rate; a value that departs from $\\dot{M}_B$ by more than a few per cent, or a flow that never reaches a stationary state, would falsify the choking claim.","tokens_in":12998,"feed_emoji":"💨","tokens_out":10261,"duration_ms":96052,"temperature":0.7,"pith_summary":"Spherically symmetric accretion has a canonical steady solution, Bondi accretion, but real infall is never perfectly spherical. This paper asks what happens when the equatorial plane is slightly denser than the poles, and finds that the flow reorganizes into a steady pattern of equatorial infall plus bipolar outflows, while the mass flux actually reaching the central object saturates at essentially the Bondi value. The excess injected material is expelled even for density contrasts as small as $0.1\\%$, and for a $10\\%$ contrast the ejected gas leaves faster than the local escape speed. The result matters because it shows that bipolar ejection can arise from pure hydrodynamics, without rotation or magnetic fields, and it suggests a universal ceiling on steady spherical accretion.","feed_headline":"Even a tiny equatorial density excess caps accretion at the Bondi rate","feed_subtitle":"A slim equator-to-pole density gap makes extra matter fly out along the poles, no magnetism needed.","key_machinery":"The central object is a named mechanism, choked accretion, produced by an equatorial-to-polar density contrast at the outer injection sphere. The analytic backbone is the incompressible, irrotational velocity potential $\\Phi = \\alpha/r\\left[1 + r^3/(4S^3)(3\\cos^2\\theta - 1)\\right]$, a Laplace-equation solution with a monopole term for Bondi inflow and a quadrupole term; its stagnation points at $(S,0)$ and $(S,\\pi)$ mark the transition from equatorial infall to polar outflow. In the numerical simulations, the same morphology emerges: a pressure gradient from the density inhomogeneity redirects streamlines, and because the flow becomes quasi-spherical at small radii, the inner accretion rate locks onto the Bondi value. The quantity $S$, the stagnation-point radius, controls the geometry and shrinks as $\\delta$ or $\\gamma$ grow.","core_discovery":"The paper claims that a steady, non-rotating, purely hydrodynamic accretion flow onto a Newtonian point mass does not simply pass through the spherical Bondi solution when the inflow is axisymmetrically over-dense at the equator. The density contrast generates a pressure gradient that deflects part of the infall, producing a bipolar outflow, while the inner flow re-converges toward the Bondi solution and the total accretion rate across the inner boundary 'chokes' at a value within a few per cent of the Bondi rate. In the simulations this holds for adiabatic indices $\\gamma = 1$, $4/3$, $7/5$ and for equatorial-to-polar contrasts of $0.1$, $1$ and $10$ per cent; the ejected fraction grows with the contrast, reaching $94\\%$ at $\\delta = 10\\%$ with outflow speeds above the local escape velocity. The authors therefore propose the choked accretion mechanism as a hydrodynamical bridge between purely radial accretion and jet-generating disc models.","pith_inferences":["If the choke is a general property of subsonic outer-injection flows, the Bondi rate acts as a strict upper bound on steady non-rotating, non-magnetic accretion; any observed super-Bondi accretion would then be a signature of rotation, magnetic stresses, or time dependence.","A natural numerical test would replace the free-outflow inner boundary with a hard absorbing surface (e.g. a star) or place the outer boundary closer than $10\\,r_\\mathrm{B}$; the paper notes the latter regime becomes unsteady, so the plateau may be restricted to gently perturbed, large-box configurations.","The equatorial-overdense boundary condition is the axisymmetric cousin of the transverse density gradient in Bondi-Hoyle-Lyttleton accretion; choked accretion might be the zero-relative-velocity limit of a larger family of gradient-driven outflows.","Observational extension: if molecular outflows in young stellar objects are choked-accretion driven, their mass-loss rates should track infall rates above a threshold set by the Bondi rate, a correlation that could be searched for in protostellar cores."],"forward_implications":["A density contrast as small as $0.1\\%$ between equator and poles is enough to switch the flow from purely radial Bondi infall to a steady inflow/outflow morphology.","The accretion rate onto the central object saturates at approximately $\\dot{M}_B$ regardless of the injection rate at the outer boundary; extra mass is ejected rather than accreted.","For a $10\\%$ contrast the ejected matter exceeds the local escape speed and carries up to $94\\%$ of the injected mass, so the mechanism can in principle launch unbound outflows.","The same qualitative result is obtained for isothermal ($\\gamma=1$), radiation-dominated ($\\gamma=4/3$) and diatomic-gas ($\\gamma=7/5$) equations of state, suggesting the mechanism is not tied to a specific thermodynamics.","Order-of-magnitude estimates place the mechanism's characteristic scale $S$ within reach of young-stellar-object molecular outflows and long-GRB collapsar conditions, though not of X-ray binaries or AGN as-is."],"supporting_citations":[{"why":"Supplies the spherically symmetric steady accretion solution and reference rate against which the choked accretion rate is measured.","marker":"Bondi (1952)"},{"why":"Perturbative isothermal analysis that first showed a polar density gradient yields inflow/outflow steady states, the result the simulations extend.","marker":"Hernandez et al. (2014)"},{"why":"Provides the incompressible, irrotational potential-flow framework used for the analytic model in Section 2.","marker":"Tejeda (2018)"},{"why":"Supplies the multipole expansion of the velocity potential for axisymmetric accretion adopted in the analytic solution.","marker":"Petrich et al. (1988)"},{"why":"Companion general-relativistic analytic model of choked accretion that motivates the incompressible toy model.","marker":"Tejeda et al. (2019)"},{"why":"Provides the HLL approximate Riemann solver used to compute numerical fluxes in the simulations.","marker":"Harten et al. (1983)"}],"fun_headline_variants":["Slight equator overdensity caps accretion, expels matter poleward","Breaking spherical symmetry chokes Bondi accretion, ignites bipolar outflows","Accretion choked by equator-pole pressure gradient, matter flies out","No magnetism needed: density asymmetry alone makes bipolar outflows","Bondi accretion capped: extra mass ejected along poles"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim that the accretion rate saturates at the Bondi value rests on the particular numerical setup — a free-outflow inner boundary at $0.1\\,r_\\mathrm{B}$ and a cosine-squared density excess injected at $10\\,r_\\mathrm{B}$ — and the paper itself reports that larger contrasts or smaller outer radii become unstable, so the choke could be a property of that setup rather than a general law.","fun_headline_variants_meta":{"raw":{"variants":["Slight equator overdensity caps accretion, expels matter poleward","Breaking spherical symmetry chokes Bondi accretion, ignites bipolar outflows","Accretion choked by equator-pole pressure gradient, matter flies out","No magnetism needed: density asymmetry alone makes bipolar outflows","Bondi accretion capped: extra mass ejected along poles"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000789,"raw_usage":{"total_tokens":3495,"prompt_tokens":975,"completion_tokens":2520,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":591,"completion_tokens_details":{"reasoning_tokens":2431}},"tokens_in":591,"tokens_out":2520,"duration_ms":20462,"temperature":1.0,"reasoning_tokens":2431,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:33:04.167020+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same hydrodynamic problem with the inner boundary changed from free outflow to a hard absorbing surface, or with the outer injection sphere moved inward from $10\\,r_\\mathrm{B}$ to a few $r_\\mathrm{B}$, and measure the steady-state accretion rate; a value that departs from $\\dot{M}_B$ by more than a few per cent, or a flow that never reaches a stationary state, would falsify the choking claim.","supporting_citations":[{"cited_title":"D., Leer B., 1983, SIAM Review, 25","cited_arxiv_id":null,"evidence_quote":"Provides the HLL approximate Riemann solver used to compute numerical fluxes in the simulations."}],"review_version":1}