REVIEW 2 major objections 4 minor 30 references
Super-Eddington accretion onto black holes and its application to fallback accretion
T0 review · 2 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read Steady spherical super-Eddington accretion onto a black hole splits into two analytic solution branches, and the predicted fallback luminosity is correspondingly faint.
desk verdict A genuinely useful analytic advance with one load-bearing self-consistency gap in the isothermal branch at high accretion rates; deserves refereeing. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is a single second-order ordinary differential equation, Eq. (15), for $\Psi$, a dimensionless measure of how far the flow departs from free fall, with $\Psi'$ proportional to the diffusive luminosity. Introducing $\chi=\Psi-2\dot{e}=2(\Lambda_a-\Lambda_d)$, twice the difference between the advective and diffusive dimensionless luminosities, and linearizing in the small parameter $\delta=1/\dot{m}$ reduces the problem to a linear second-order ODE, Eq. (18). The sign of $\chi$ near the bifurcation point $x_{\rm bf}\simeq 1-\delta$ selects the branch: $\chi<0$ produces a boundary layer at $x=1$ and the isothermal solutions, while $\chi>0$ produces a boundary layer at $x_{\rm bf}$ and the power-law solutions. Boundary layer theory supplies the analytic approximations, and because the sonic point lies outside the photosphere, the integration constant $\dot{e}$ is fixed by outer boundary conditions rather than by smooth passage through the singularity.
What would settle it
Integrate the full nonlinear equation (15) with the photosphere placed inside the adiabatic sonic radius; if a smooth transonic solution appears that matches neither analytic branch, the two-branch classification fails. Observationally, a confirmed fallback-powered transient whose late-time luminosity follows the spherical $t^{-5/3}$ law would validate the branch, while a brighter source with a shallower decay would favor the turbulent-envelope interpretation.
Extended reading notes
Core claim
For a given dimensionless accretion rate $\dot{m}\gg 1$, steady spherical radiation-pressure-dominated accretion in the free-fall limit admits two types of solutions: flows that become isothermal at large radii and flows whose temperature at infinity vanishes as a power law of radius. The trapping radius, where the advective and diffusive luminosities balance, is $r_{\rm tr}=3\dot{m} r_S$ on the isothermal branch and $r_{\rm tr}=(\dot{m}/2)r_S$ on the power-law branch. Beyond the trapping radius on the isothermal branch the diffusive and advective luminosities nearly cancel, so the rest-frame luminosity at the photosphere is a factor of order $\dot{m}$ smaller than the diffusive luminosity. In the fallback application, a fully ionized inflow gives $L_{\rm obs}\propto \dot{m}^2$ early, $\propto \dot{m}$ in the cooling-dominated regime, and $\propto \dot{m}^{5/3}$ once gas pressure dominates; with recombination, the ionization luminosity is $L_{\rm ion}\propto \dot{m}$ for radiation-dominated fronts and $\propto \dot{m}^{4/3}$ for gas-dominated fronts. At about five years after collapse, the spherical fallback luminosity from red or yellow supergiant progenitors is several orders of magnitude below the luminosities of the NGC 6946 and M31 fading sources, so those sources are instead likely powered by accretion of the turbulent convective envelope.
Load-bearing premise
The classification into two branches rests on the flow staying essentially in free fall out to the photosphere, so the adiabatic sonic point remains outside the solved region and the integration constant $\dot{e}$ is fixed by outer boundary conditions; if the flow becomes subsonic near the photosphere, the two-branch picture would need to be revised.
Editorial extensions
If this is right
- - The trapping radius differs by a factor of six between the branches ($3\dot{m} r_S$ versus $(\dot{m}/2)r_S$), so the radius at which photons first escape depends on which branch the flow selects.
- - The photospheric rest-frame luminosity is suppressed by $1/\dot{m}$ relative to the diffusive luminosity, so Eddington-scaled estimates that neglect this suppression overpredict the light from spherical super-Eddington accretion.
- - For a fully ionized fallback inflow, the predicted light curve steepens from $t^{-10/3}$ (proportional to $\dot{m}^2$) to $t^{-5/3}$ (proportional to $\dot{m}$), and finally to $t^{-25/9}$ (proportional to $\dot{m}^{5/3}$) as gas pressure becomes dominant.
- - With recombination, the ionization luminosity scales as $L_{\rm ion}\propto \dot{m}$ for radiation-pressure-dominated fronts and $L_{\rm ion}\propto \dot{m}^{4/3}$ for gas-pressure-dominated fronts, agreeing with earlier radiation-hydro simulations.
- - Spherical fallback accretion after a failed supernova is far fainter than the fading sources in NGC 6946 and M31 at ages of several years, so those sources likely require accretion of the turbulent convective envelope.
Reading between the lines
- - A testable extension: because the two branches have different trapping radii and temperature profiles, time-dependent or slightly non-spherical implementations of the same equations should reveal which branch the boundary conditions select, allowing simulations to be checked against the analytic scalings.
- - The two ionization scalings imply a discriminating observational test: a confirmed fallback-powered transient whose late-time luminosity decays as $t^{-5/3}$ would be in the radiation-dominated regime, while a $t^{-20/9}$ decay would indicate gas-pressure dominance.
- - If turbulent convective envelope accretion powers the observed fading sources, their late-time luminosities should track the envelope's binding energy and convective properties rather than the black hole mass alone; comparing several failed-supernova candidates would test this.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies steady, spherical, radiation-pressure-dominated super-Eddington accretion onto a Schwarzschild black hole in the free-fall limit. The authors reduce the fluid equations to a second-order ODE for the dimensionless variable Psi, linearize it by assuming small deviations from free fall, and use boundary-layer theory to obtain analytic solutions for two branches: an isothermal branch with trapping radius r_tr = 3 mdot r_S and rest-frame luminosity suppressed by a factor ~mdot^{-1} relative to the diffusive luminosity, and a power-law branch with r_tr = (mdot/2) r_S. These scalings are then applied to fallback accretion following failed supernovae, including a treatment of recombination/ionization fronts, and are used to argue that spherical fallback accretion cannot explain the observed fading sources in NGC 6946 and M31; the authors instead attribute those sources to accretion of the turbulent convective envelope of the supergiant progenitor.
Significance. If the analytic solutions are valid, the paper provides compact, parameter-free scalings for a problem with a long theoretical history, and it connects those scalings to concrete observational predictions. The comparison with the Zampieri et al. (1998) numerical simulations in Fig. 8, the use of realistic MESA progenitor profiles, and the explicit testable prediction that slowly rotating Wolf-Rayet or blue-supergiant collapses should be genuinely faint are clear strengths. The internal agreement between the boundary-layer solution and the numerical integration of the nonlinear ODE in Figs. 4 and 6 is also a valuable check. The paper is generally clearly written and the astrophysical application is timely.
major comments (2)
- [§3 (case z_bf<0) and §6, summary item 1] The claim that the adiabatic sonic point 'does not play a role in the dynamics of the flow inside the photosphere' is not reconciled with the statement in §3 that in the isothermal branch 'the flow may become subsonic either slightly below or above the photosphere.' Using Eq. (11) with f_tau=1, Eq. (9), and the isothermal-branch pressure profile (39), the photospheric Mach number is M_ph^2 ~ 5e12 (M/10 M_sun)^{-1} (T_ph/10^4 K)^{-4} mdot^{-4} up to order-unity factors. For T_ph = 10^4 K and M = 10 M_sun, M_ph drops below unity once mdot is of order 10^3, which is within the range shown in Figs. 1 and 2 and used in the fallback application. In that regime the sonic point lies inside the photosphere-to-horizon domain, so the eigenvalue e_dot is no longer fixed by the outer boundary conditions, and the boundary-layer solution (35)-(36) is not a self-consistent free-fall solution. The paper should either derive and impose an explicit supersonicity criterion and restrict the two-branch classification to the regime where it holds, or redo the eigenvalue analysis with the sonic-point crossing condition; it should also state whether the fallback conclusions in §5 survive in the regime actually realized by the MESA models.
- [§3, Eq. (18)] The linearization of Eq. (15) by discarding the Psi^2, Psi'^2, Psi e_dot, and Psi' e_dot terms is not justified by a controlled asymptotic expansion. The supporting argument that O(e_dot) = Psi is an order-of-magnitude estimate (Appendix B), and the resulting Eq. (18) is validated numerically only for the specific models shown in Figs. 4 and 6 (mdot = 700, M = 10 M_sun for the isothermal branch; mdot = 700 for the power-law branch). The analytic scalings L proportional to mdot^2 and r_tr = 3 mdot r_S used throughout §5 are derived from this linearized equation, but the fallback application uses accretion rates up to mdot ~ 10^12 (Fig. 10). Please provide an explicit estimate of the neglected terms over the full parameter range used in the paper, or restrict the analytic predictions to the region where the linearization is controlled.
minor comments (4)
- [§2.2, Eq. (15)] Equation (15) is presented without derivation; since it is the foundation of the subsequent analysis, a derivation in an appendix (or a clear reference to a previous work) would help the reader verify the coefficients F1 and F2.
- [§5, Eq. (53)] The criterion in Eq. (53) is stated as an order-of-magnitude inequality, and the text concludes that the flow 'will thus likely settle onto the isothermal branch,' while the summary presents this as a definite selection. Please give the quantitative condition for the actual MESA models and state the resulting branch choice explicitly.
- [Fig. 8 caption] Figure 8 would be easier to use if each segment of the luminosity curve were labeled with its governing equation number (Eq. 56, Eq. 57, and the gas-pressure-dominated scaling) directly in the caption, rather than only by the inequalities in the legend.
- [§5.2, Eqs. (60)-(61)] The statement that T_g,ph is 'typically of order T_rec/2' is imported from Zampieri et al. (1998); please state explicitly how this estimate enters the normalization of Eqs. (60) and (61), since those normalizations set the vertical offset of the dashed lines in Fig. 8.
Circularity Check
No circularity found: the two branches, trapping radii, and luminosity scalings are derived from the governing ODE with boundary-condition-determined constants, then benchmarked against numerical and external simulations.
full rationale
The analytic derivation is self-contained: Eq. (18) (the free-fall linearization of Eq. (15)) is solved with boundary-layer expansions, with branch selection determined by the sign of z = chi'/chi near x_bf (Sec. 3) and with integration constants fixed by the photospheric temperature/pressure boundary condition through Eqs. (9) and (11), not by the luminosities being predicted. The trapping radii in Eqs. (29) and (36) are zeros of the resulting chi solutions, and the mdot^2, mdot, and mdot^(4/3) scalings in Secs. 5.1-5.2 follow algebraically from the analytic pressure/velocity profiles, radiation-pressure balance, and ionization/recombination energetics. The amplitudes are checked against exact numerical integrations of the nonlinear ODE (Figs. 4 and 6) and against the external Zampieri et al. (1998) simulations, so the fallback comparisons are not fit to the simulated points. The self-citation to Faran et al. (2019) supplies a recombination-layer optical-depth approximation that is benchmarked against those independent simulations and does not underpin the central two-branch result. Two manuscript caveats are noted: Sec. 3 states that in the isothermal branch the flow 'may become subsonic either slightly below or above the photosphere,' which is a domain-consistency caveat to the free-fall assumption, and the MESA models are cited as 'Antoni et al. in prep'; neither is a circularity because no prediction is defined in terms of the target luminosity and no fitted parameter is renamed as a prediction.
Assumptions & free parameters
free parameters (3)
- C1 (isothermal branch amplitude) =
approximately 1.3e-12 mdot^(1/2) M10 T_inf,4^4, calibrated to the numerical ODE solution
- T_ph (photospheric temperature) =
about 1e4 K, assumed for the fully ionized case
- E_sh (shock/explosion energy) =
1e48 erg, chosen for the fiducial light curve
assumptions (5)
- domain assumption Steady, spherically symmetric accretion with radiation pressure dominating gas pressure (gamma = 4/3) and Thomson-scattering opacity, in a Paczynski-Wiita potential, with photon transport treated by the diffusion approximation and the outer boundary at optical depth unity.
- ad hoc to paper The flow is close enough to free fall that Eq. (15) can be linearized into Eq. (18) by dropping Psi^2, Psi'^2, Psi*edot, and Psi'*edot terms.
- domain assumption The adiabatic sonic point (movable singularity) lies outside the solution domain, so the eigenvalue edot is fixed by outer boundary conditions rather than sonic-point crossing.
- domain assumption Fallback accretion rate follows Chevalier (1989)'s self-similar solution Mdot = (8/9) pi^(5/3) (2GM) rho_0 t_0^(8/3) t^(-5/3), with initial conditions from MESA progenitor models.
- ad hoc to paper The photospheric temperature is pinned near 1e4 K by the balance of compressional heating with atomic cooling, and the recombination/ionization front treatment follows Faran et al. (2019).
Cite this review
Pith. "Pith review of Super-Eddington accretion onto black holes and its application to fallback accretion." pith.science (2026). https://pith.science/paper/4C7RI42M
@misc{pith2026250716893,
author = {Pith},
title = {Pith review of: Super-Eddington accretion onto black holes and its application to fallback accretion},
year = {2026},
howpublished = {\url{https://pith.science/paper/4C7RI42M}},
note = {Machine review of arXiv:2507.16893}
}
read the original abstract
We study the problem of steady-state spherical accretion onto a black hole, in which the internal energy of the flow is governed by radiation and photon diffusion dominates the energy flux at large radii. In the free-fall limit, the fluid equations can admit two types of solutions for a given accretion rate: (1) accretion flows that become isothermal at large radii and (2) solutions in which the temperature at infinity vanishes as a power law of the radius. Using boundary layer theory, we obtain analytic solutions for the two cases and apply our results to fallback accretion onto a black hole following a failed supernova explosion. We give predictions for the observational signature of fallback accretion using realistic progenitor properties from MESA, both for a fully ionized inflow and for the more realistic case in which recombination/ionization take place due to low photospheric temperatures. The observed fading sources coincident with the failed-supernova candidates in NGC 6946 and M31 are too luminous to be powered by spherical accretion onto newly formed black holes; the observed sources are instead likely due to accretion of the turbulent, convective envelope of the supergiant progenitor.
Figures
Figures from the paper (7 more)
Reference graph
Works this paper leans on
-
[1]
Adams, S. M., Kochanek, C. S., Gerke, J. R., Stanek, K. Z., & Dai, X. 2017, MNRAS, 468, 4968, doi: 10.1093/mnras/stx816
-
[2]
2023, MNRAS, 525, 1229, doi: 10.1093/mnras/stad2328
Antoni, A., & Quataert, E. 2023, MNRAS, 525, 1229, doi: 10.1093/mnras/stad2328
-
[3]
Beasor, E. R., Smith, N., & Jencson, J. E. 2025, ApJ, 979, 117, doi: 10.3847/1538-4357/ad8f3f
-
[4]
Begelman, M. C. 1978, MNRAS, 184, 53, doi: 10.1093/mnras/184.1.53
-
[5]
Bender, C. M., & Orszag, S. A. 1978, Advanced Mathematical Methods for Scientists and Engineers 18
work page 1978
-
[6]
Blondin, J. M. 1986, ApJ, 308, 755, doi: 10.1086/164548
-
[7]
1974, ApJ, 189, 147, doi: 10.1086/152780
Buff, J., & McCray, R. 1974, ApJ, 189, 147, doi: 10.1086/152780
doi:10.1086/152780 1974
-
[8]
2021, Nature, 589, 29, doi: 10.1038/s41586-020-03059-w
Burrows, A., & Vartanyan, D. 2021, Nature, 589, 29, doi: 10.1038/s41586-020-03059-w
Show all 30 references
- [9]
-
[10]
Chevalier, R. A. 1989, ApJ, 346, 847, doi: 10.1086/168066
1989 doi
-
[11]
Davies, B., & Beasor, E. R. 2018, MNRAS, 474, 2116, doi: 10.1093/mnras/stx2734
2018 doi
-
[12]
E., et al
De, K., MacLeod, M., Jencson, J. E., et al. 2024, arXiv e-prints, arXiv:2410.14778, doi: 10.48550/arXiv.2410.14778
2024 doi
-
[13]
2016, ApJ, 818, 124, doi: 10.3847/0004-637X/818/2/124
Ugliano, M. 2016, ApJ, 818, 124, doi: 10.3847/0004-637X/818/2/124
2016 doi
-
[14]
2019, ApJ, 879, 20, doi: 10.3847/1538-4357/ab218a Fern´ andez, R., Quataert, E., Kashiyama, K., & Coughlin, E
Faran, T., Goldfriend, T., Nakar, E., & Sari, R. 2019, ApJ, 879, 20, doi: 10.3847/1538-4357/ab218a Fern´ andez, R., Quataert, E., Kashiyama, K., & Coughlin, E. R. 2018, MNRAS, 476, 2366, doi: 10.1093/mnras/sty306
2019 doi
-
[15]
Flammang, R. A. 1982, MNRAS, 199, 833, doi: 10.1093/mnras/199.4.833 —. 1984, MNRAS, 206, 589, doi: 10.1093/mnras/206.3.589
1982 doi
-
[16]
2014, MNRAS, 439, 4011, doi: 10.1093/mnras/stu257
Gilkis, A., & Soker, N. 2014, MNRAS, 439, 4011, doi: 10.1093/mnras/stu257
2014 doi
-
[17]
2021, ApJ, 911, 6, doi: 10.3847/1538-4357/abe59e
Ivanov, M., & Fern´ andez, R. 2021, ApJ, 911, 6, doi: 10.3847/1538-4357/abe59e
2021 doi
-
[18]
S., Neustadt, J
Kochanek, C. S., Neustadt, J. M. M., & Stanek, K. Z. 2024, ApJ, 962, 145, doi: 10.3847/1538-4357/ad18d7
2024 doi
-
[19]
Lovegrove, E., & Woosley, S. E. 2013, ApJ, 769, 109, doi: 10.1088/0004-637X/769/2/109
2013 doi
-
[20]
E., & Zhang, W
Lovegrove, E., Woosley, S. E., & Zhang, W. 2017, ApJ, 845, 103, doi: 10.3847/1538-4357/aa7b7d
2017 doi
-
[21]
Nadezhin, D. K. 1980, Ap&SS, 69, 115, doi: 10.1007/BF00638971
1980 doi
-
[22]
2010, ApJ, 725, 904, doi: 10.1088/0004-637X/725/1/904
Nakar, E., & Sari, R. 2010, ApJ, 725, 904, doi: 10.1088/0004-637X/725/1/904
2010 doi
-
[23]
Quataert, E., Lecoanet, D., & Coughlin, E. R. 2019, MNRAS, 485, L83, doi: 10.1093/mnrasl/slz031
2019 doi
-
[24]
2012, ApJL, 750, L31, doi: 10.1088/2041-8205/750/2/L31
Robertson, B., & Goldreich, P. 2012, ApJL, 750, L31, doi: 10.1088/2041-8205/750/2/L31
2012 doi
-
[25]
2009, A&A, 508, 751, doi: 10.1051/0004-6361/200912495
Vink, J. 2009, A&A, 508, 751, doi: 10.1051/0004-6361/200912495
2009 doi
-
[26]
Smartt, S. J. 2015, PASA, 32, e016, doi: 10.1017/pasa.2015.17
2015 doi
-
[27]
J., Eldridge, J
Smartt, S. J., Eldridge, J. J., Crockett, R. M., & Maund, J. R. 2009, MNRAS, 395, 1409, doi: 10.1111/j.1365-2966.2009.14506.x
2009
-
[28]
Janka, H. T. 2016, ApJ, 821, 38, doi: 10.3847/0004-637X/821/1/38
2016 doi
-
[29]
1984, ApJ, 284, 394, doi: 10.1086/162419
Vitello, P. 1984, ApJ, 284, 394, doi: 10.1086/162419
1984 doi
-
[30]
L., & Wasserman, I
Zampieri, L., Colpi, M., Shapiro, S. L., & Wasserman, I. 1998, ApJ, 505, 876, doi: 10.1086/306192
1998 doi
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.