pith. machine review for the scientific record. sign in

arxiv: 2604.14708 · v1 · submitted 2026-04-16 · 🌌 astro-ph.HE

Recognition: unknown

Influence of winds on shocked magnetized viscous accretion flows around rotating black holes

Authors on Pith no claims yet

Pith reviewed 2026-05-10 10:51 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords accretion disksblack holeswindsstanding shocksviscous flowsmagnetized accretionrotating black holestransonic solutions
0
0 comments X

The pith

Winds in black hole accretion flows reduce disk luminosity and eliminate standing shocks above a critical mass-loss strength.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper models relativistic magnetized viscous accretion onto spinning black holes while letting the accretion rate drop inward as a power law to represent mass lost in winds. Solving the coupled flow equations shows that this mass and angular momentum loss changes the global transonic structure, lowers overall radiative output, and shifts the location, compression, and strength of any standing shocks that form. The authors identify a critical wind parameter p^crit beyond which steady shocks disappear entirely, and they show that higher viscosity or stronger angular-momentum extraction by the wind lowers this critical value. The central object is therefore the coupled interplay between wind-driven mass loss, viscosity, and shock formation in the disk.

Core claim

Incorporating a radially decreasing mass-accretion rate to account for winds yields global transonic solutions in which standing shocks are possible only for wind parameters below a critical value p^crit; above that threshold no steady shocks exist, the disk luminosity drops markedly, and the remaining shock properties (radius, compression ratio, and strength) are strongly modified, with the critical threshold itself decreasing when viscosity or wind angular-momentum extraction increases.

What carries the argument

The power-law radial decline in mass-accretion rate that self-consistently removes mass and angular momentum from the flow while the governing relativistic viscous MHD equations (including toroidal fields and synchrotron cooling) are integrated for global transonic solutions.

If this is right

  • Steady standing shocks exist only for wind strengths below the critical p^crit.
  • Higher viscosity shrinks the range of wind parameters that still permit shocks.
  • Stronger angular-momentum removal by winds further lowers p^crit.
  • Disk luminosity decreases once winds are included at any appreciable level.
  • Shock radius moves outward and compression ratio and strength change as wind strength rises.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same wind prescription could be tested in non-relativistic or weakly magnetized disks to see whether the critical p^crit behavior persists.
  • Observed changes in X-ray luminosity or quasi-periodic oscillation frequencies might be used to constrain the wind parameter p in real sources.
  • If the power-law mass-loss assumption is replaced by a more detailed wind model, the location of p^crit could shift, offering a direct test of the simplified treatment.

Load-bearing premise

Mass loss in winds can be represented by a simple power-law drop in the local accretion rate with radius without needing to solve the wind-launching equations separately.

What would settle it

Detection of a standing shock in an accretion flow whose measured wind mass-loss index exceeds the calculated p^crit for the observed viscosity and black-hole spin would falsify the claimed disappearance of shocks.

Figures

Figures reproduced from arXiv: 2604.14708 by Camelia Jana, Santabrata Das.

Figure 1
Figure 1. Figure 1: Variation of Mach number (M) with radial dis￾tance (x) for different wind parameters (p) for flows passing through the inner critical points (xin). Here, we choose the model parameters as αB = 0.02, ak = 0.0, βin = 100.0, and m˙ in = 0.01. The left panel corresponds to l = 0.8, xin = 5.14, and λin = 3.42, while the right panel represent the case with l = 1.2, xin = 6.30, and λin = 3.46. See text for the de… view at source ↗
Figure 2
Figure 2. Figure 2: Variation of Mach number (M) with radial distance (x) for different wind parameters (p) for flow containing outer critical points (xout). The flow is injected from xedge = 103 with ˙medge = 0.1 and βedge = 104 , while αB = 0.2 and ak = 0.0 are kept fixed. In panel (a), we consider l = 0.8, Eedge = 1.95 × 10−3 , and λedge = 3.60, whereas in panel (b), l = 1.2, Eedge = 1.58 × 10−3 , and λedge = 4.65 are chos… view at source ↗
Figure 3
Figure 3. Figure 3: Radial variation of flow variables for different wind parameters p. Here, we fix l = 0.8, ak = 0.0, and αB = 0.02 and flow is injected from xedge = 1000 with Eedge = 1.95 × 10−3 , λedge = 3.60, ˙medge = 0.1, and βedge = 104 . The solid (blue), dashed (red), and dotted (green) curves correspond to p = 0.0, 0.2, and 0.4, respectively. See text for further details. that mass loss effectively cools down the di… view at source ↗
Figure 5
Figure 5. Figure 5: Variation of the disk luminosity (Ldisk) as a func￾tion of wind parameter (p) for a set of plasma-β param￾eters. Here, the flow is injected from xedge = 103 with Eedge = 1.95 × 10−3 , λedge = 3.60, ˙medge = 0.1 and we set l = 0.8, αB = 0.02 and ak = 0.0. Circles, diamonds, and asterisks connected with solid (blue), dashed (red), and dotted (green) lines correspond to βedge = 102 , 103 , and 104 , respectiv… view at source ↗
Figure 7
Figure 7. Figure 7: Plot of shock-mediated global accretion solution around non-rotating (ak = 0.0; left panels) and rapidly ro￾tating (ak = 0.95; right panels) black holes. Panels (a) and (d) present the variation of the Mach number (M), while panels (b) and (e) display the corresponding density profiles (ρ) as functions of the radial coordinate (x). The colour bars in the upper and lower panels indicate the distribu￾tion of… view at source ↗
Figure 8
Figure 8. Figure 8: Modification of shock-induced accretion solutions in the presence of wind. Here, we fix l = 1.0 and αB = 0.02. In panel (a), a non-rotating (ak = 0.0) black hole is con￾sidered, and the outer boundary is set at xedge = 1000 with input parameters Eedge = 0.00158, λedge = 4.818, βedge = 5 × 104 , and ˙medge = 0.1. The solid (blue), dashed (red), and dotted (green) curves are plotted for wind pa￾rameters p = … view at source ↗
Figure 9
Figure 9. Figure 9 [PITH_FULL_IMAGE:figures/full_fig_p012_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Variation of the critical wind parameter p crit as a function of αB for different values of l. Here, βin = 100, m˙ in = 0.01, and ak = 0.95 are adopted. See text for the details. out this analysis for a rotating black hole with spin ak = 0.95, and choose βin = 100 and ˙min = 0.01 at the inner critical point xin. By varying the flow energy (Ein) and angular momentum (λin) at xin, we determine p crit for gi… view at source ↗
read the original abstract

We study global transonic solution for a relativistic, magnetized, viscous advective accretion flow around a rotating black hole, incorporating the effects of mass and angular momentum loss through winds. Our model considers dominant toroidal magnetic fields with synchrotron radiation as the primary cooling mechanism. To self-consistently model mass loss, the mass accretion rate is prescribed to decrease inward as a power-law with disk radius. With this, we solve the governing equations that describe the accretion flows in presence of winds and obtain the flow structure in terms of the inflow parameters (energy $\mathcal{E}$, angular momentum $\lambda$, plasma-$\beta$, accretion rate $\dot{m}$, and viscosity $\alpha_{\rm B}$), the wind parameters ($p$, governing mass loss; and $l$, governing angular momentum transport by winds), and the black hole spin ($a_{\rm k}$). Our analysis reveals that winds substantially modify the accretion flow leading to a significant decrease in disk luminosity. We specifically identify global solutions that admit standing shocks and find that winds profoundly alter shock properties, such as the shock radius ($x_{\rm s}$), compression ratio ($R$), and shock strength ($S$). Furthermore, we determine the critical wind parameter $p^{\rm crit}$ beyond which steady shock solutions cease to exist. We demonstrate that increased viscosity and strong angular momentum extraction by winds lead to reduce $p^{\rm crit}$. These findings evidently highlight a complex interplay between viscosity and winds in governing the dynamics of shock formation in accretion disks.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 1 minor

Summary. The manuscript develops a model for global transonic solutions of relativistic, magnetized, viscous accretion flows around rotating black holes that incorporates mass and angular momentum loss through winds. The mass accretion rate is modeled to decrease inward following a power-law dependence on radius. The authors solve for the flow structure and identify solutions with standing shocks, showing that winds reduce disk luminosity, modify shock radius, compression ratio, and strength, and that there exists a critical wind parameter p^crit above which steady shocks no longer exist. They further show that increased viscosity and angular momentum extraction by winds lower this critical value.

Significance. If the modeling assumptions hold, this work contributes to understanding the role of winds in regulating accretion dynamics and shock formation in black hole systems. The identification of p^crit and its dependence on viscosity provides a potentially useful diagnostic for when shocks can persist in the presence of outflows. The inclusion of toroidal magnetic fields and synchrotron cooling adds to the realism for applications to observed high-energy phenomena.

major comments (1)
  1. [Model setup and governing equations] The central modeling choice is the prescription that the mass accretion rate decreases inward as a power-law with disk radius (Ṁ(r) ∝ r^p) to self-consistently model mass loss through winds. This functional form is imposed rather than derived from vertical wind-launching conditions such as centrifugal barrier, magnetic torque, or thermal driving. Because the reported global transonic solutions, standing shocks, altered x_s/R/S, reduced luminosity, and the existence and value of p^crit all follow from integrating the conservation equations under this constraint, the critical boundary is tied to the specific power-law assumption. A derivation from first principles or a sensitivity study to alternative Ṁ(r) profiles is required to establish whether p^crit is a robust physical feature.
minor comments (1)
  1. [Abstract and parameter list] The abstract lists inflow parameters including plasma-β but does not specify how β is initialized or evolved in the global solutions; clarifying this in the methods would aid reproducibility.

Simulated Author's Rebuttal

1 responses · 1 unresolved

We thank the referee for their constructive comments, which help clarify the scope and limitations of our modeling approach. We address the major comment point by point below.

read point-by-point responses
  1. Referee: The central modeling choice is the prescription that the mass accretion rate decreases inward as a power-law with disk radius (Ṁ(r) ∝ r^p) to self-consistently model mass loss through winds. This functional form is imposed rather than derived from vertical wind-launching conditions such as centrifugal barrier, magnetic torque, or thermal driving. Because the reported global transonic solutions, standing shocks, altered x_s/R/S, reduced luminosity, and the existence and value of p^crit all follow from integrating the conservation equations under this constraint, the critical boundary is tied to the specific power-law assumption. A derivation from first principles or a sensitivity study to alternative Ṁ(r) profiles is required to establish whether p^crit is a robust physical feature.

    Authors: We acknowledge that the power-law form Ṁ(r) ∝ r^p is a phenomenological prescription rather than derived from explicit vertical wind-launching physics. This choice follows standard practice in the literature on wind-affected accretion flows (e.g., as used in studies of mass-loss in advective disks) to enable self-consistent radial integration of the conservation equations while capturing the net effect of outflows. A first-principles derivation would require resolving the vertical structure and launching mechanisms (centrifugal, magnetic, or thermal), which necessitates at least 2D axisymmetric simulations and is beyond the scope of the present 1D radial model. We have, however, extensively varied p across a range of values to map the existence and properties of shocks, including the identification of p^crit and its dependence on viscosity and angular momentum loss. To strengthen the manuscript, we will add a dedicated paragraph in the model section discussing the rationale for the power-law assumption, its common usage in similar works, and the associated limitations, along with a brief note that future multi-dimensional extensions could test alternative profiles. revision: partial

standing simulated objections not resolved
  • A derivation of the Ṁ(r) profile from first-principles vertical wind-launching conditions (centrifugal barrier, magnetic torque, or thermal driving), as this requires multi-dimensional modeling outside the current 1D framework.

Circularity Check

0 steps flagged

No significant circularity; parametric study within explicit modeling ansatz

full rationale

The paper adopts an explicit power-law prescription for the radial decline of mass accretion rate (Ṁ(r) ∝ r^p) as the modeling choice to represent wind-driven mass loss, then integrates the conservation equations for energy, angular momentum, and magnetic flux under that fixed functional form to obtain transonic solutions, standing shocks, and the boundary value p^crit at which shocks cease to exist. This is a standard parameter-space exploration in which p is an input that is varied to map solution existence; p^crit emerges from the numerical integration and boundary conditions rather than being identical to the prescription by construction. No self-definitional loop, fitted-input-as-prediction, or load-bearing self-citation is present in the derivation chain.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 0 invented entities

Abstract-only review yields limited visibility into assumptions; the power-law mass-loss prescription is an ad-hoc modeling choice required to close the system.

free parameters (3)
  • p
    Wind parameter controlling the radial power-law index of mass loss; central to locating p^crit.
  • l
    Wind parameter controlling angular momentum extraction; affects p^crit.
  • alpha_B
    Viscosity parameter whose increase is shown to lower p^crit.
axioms (2)
  • domain assumption Dominant toroidal magnetic fields with synchrotron radiation as the primary cooling mechanism.
    Stated as the basis for the energy equation.
  • ad hoc to paper Mass accretion rate decreases inward as a power-law with disk radius.
    Explicitly prescribed to incorporate wind mass loss.

pith-pipeline@v0.9.0 · 5567 in / 1467 out tokens · 40020 ms · 2026-05-10T10:51:34.855159+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

64 extracted references · 60 canonical work pages

  1. [1]

    The Merger Rates and Mass Assembly Histories of Dark Matter Haloes in the Two

    Abbassi, S., Ghanbari, J., & Ghasemnezhad, M. 2010, MNRAS, 409, 1113, doi: 10.1111/j.1365-2966.2010.17364.x

  2. [2]
  3. [3]

    2000, A&A, 355, 271, doi: 10.48550/arXiv.astro-ph/0001103

    Belloni, T., Klein-Wolt, M., M´ endez, M., van der Klis, M., & van Paradijs, J. 2000, A&A, 355, 271, doi: 10.48550/arXiv.astro-ph/0001103

  4. [4]

    , keywords =

    Blandford, R. D., & Begelman, M. C. 1999, MNRAS, 303, L1, doi: 10.1046/j.1365-8711.1999.02358.x —. 2004, MNRAS, 349, 68, doi: 10.1111/j.1365-2966.2004.07425.x

  5. [5]

    , keywords =

    Blandford, R. D., & Payne, D. G. 1982, MNRAS, 199, 883, doi: 10.1093/mnras/199.4.883

  6. [6]

    and Wilkinson, M

    Bu, D.-F., Yuan, F., & Xie, F.-G. 2009, MNRAS, 392, 325, doi: 10.1111/j.1365-2966.2008.14047.x

  7. [7]

    Chakrabarti, S., & Titarchuk, L. G. 1995, ApJ, 455, 623, doi: 10.1086/176610

  8. [8]

    Chakrabarti, S. K. 1989, ApJ, 347, 365, doi: 10.1086/168125 —. 1999, A&A, 351, 185, doi: 10.48550/arXiv.astro-ph/9910014

  9. [9]

    2022, ApJ, 941, 30, doi: 10.3847/1538-4357/ac9d97

    Chatterjee, K., & Narayan, R. 2022, ApJ, 941, 30, doi: 10.3847/1538-4357/ac9d97

  10. [10]

    Chattopadhyay, I., & Chakrabarti, S. K. 2000, International Journal of Modern Physics D, 9, 717, doi: 10.1142/S0218271800000670

  11. [11]

    2008 , month = jan, journal =

    Das, S. 2007, MNRAS, 376, 1659, doi: 10.1111/j.1365-2966.2007.11501.x

  12. [12]
  13. [13]

    Das, S., Chattopadhyay, I., & Chakrabarti, S. K. 2001, ApJ, 557, 983, doi: 10.1086/321692

  14. [14]

    , keywords =

    Das, S., & Sarkar, B. 2018, MNRAS, 480, 3446, doi: 10.1093/mnras/sty2071

  15. [15]

    K., Das, S., Maity, D., & Chakrabarti, S

    Dihingia, I. K., Das, S., Maity, D., & Chakrabarti, S. 2018a, Phys. Rev. D, 98, 083004, doi: 10.1103/PhysRevD.98.083004

  16. [16]

    and Das, Santabrata and Maity, Debaprasad and Nandi, Anuj

    Dihingia, I. K., Das, S., Maity, D., & Nandi, A. 2019, MNRAS, 488, 2412, doi: 10.1093/mnras/stz1933

  17. [17]

    K., Das, S., & Mandal, S

    Dihingia, I. K., Das, S., & Mandal, S. 2018b, MNRAS, 475, 2164, doi: 10.1093/mnras/stx3269

  18. [18]

    K., Das, S., Prabhakar, G., & Mandal, S

    Dihingia, I. K., Das, S., Prabhakar, G., & Mandal, S. 2020, MNRAS, 496, 3043, doi: 10.1093/mnras/staa1687

  19. [19]

    2012 , month = apr, journal =

    Faghei, K., & Mollatayefeh, A. 2012, MNRAS, 422, 672, doi: 10.1111/j.1365-2966.2012.20645.x

  20. [20]

    2004, MNRAS, 351, 1379, doi: 10.1111/j.1365-2966.2004.07876.x

    Fender, R. P., Belloni, T. M., & Gallo, E. 2004, MNRAS, 355, 1105, doi: 10.1111/j.1365-2966.2004.08384.x

  21. [21]

    2007, ApJ, 669, 85, doi: 10.1086/521578

    Fukumura, K., & Kazanas, D. 2007, ApJ, 669, 85, doi: 10.1086/521578

  22. [22]

    2004, ApJ, 611, 964, doi: 10.1086/422243

    Fukumura, K., & Tsuruta, S. 2004, ApJ, 611, 964, doi: 10.1086/422243

  23. [23]

    H., & Stone, J

    Hirose, S., Krolik, J. H., & Stone, J. M. 2006, ApJ, 640, 901, doi: 10.1086/499153

  24. [24]

    2015, ApJ, 807, 108, doi: 10.1088/0004-637X/807/1/108

    Iyer, N., Nandi, A., & Mandal, S. 2015, ApJ, 807, 108, doi: 10.1088/0004-637X/807/1/108

  25. [25]

    Jana, C., & Das, S. 2024, J. Cosmology Astropart. Phys., 2024, 075, doi: 10.1088/1475-7516/2024/07/075

  26. [26]

    1993, MNRAS, 260, 317, doi: 10.1093/mnras/260.2.317

    Kato, S., Wu, X.-B., Yang, L.-T., & Yang, Z.-L. 1993, MNRAS, 260, 317, doi: 10.1093/mnras/260.2.317

  27. [27]

    , keywords =

    Knigge, C. 1999, MNRAS, 309, 409, doi: 10.1046/j.1365-8711.1999.02839.x

  28. [28]

    1999, ApJ, 522, 727, doi: 10.1086/307667

    Koide, S., Shibata, K., & Kudoh, T. 1999, ApJ, 522, 727, doi: 10.1086/307667

  29. [29]

    D., & Lifshitz, E

    Landau, L. D., & Lifshitz, E. M. 1959, Fluid mechanics

  30. [30]

    Le, T., & Becker, P. A. 2005, ApJ, 632, 476, doi: 10.1086/432927

  31. [31]

    doi:10.1111/j.1365-2966.2009.15598.x , archivePrefix =

    Li, S.-L., & Cao, X. 2009, MNRAS, 400, 1734, doi: 10.1111/j.1365-2966.2009.15595.x

  32. [32]

    1999, ApJ, 523, 340, doi: 10.1086/307725

    Lu, J.-F., Gu, W.-M., & Yuan, F. 1999, ApJ, 523, 340, doi: 10.1086/307725

  33. [33]

    E., & Matsumoto, R

    Machida, M., Nakamura, K. E., & Matsumoto, R. 2006, PASJ, 58, 193, doi: 10.1093/pasj/58.1.193

  34. [34]

    K., & Nandi, A

    Majumder, S., Das, S., Agrawal, V. K., & Nandi, A. 2023, MNRAS, 526, 2086, doi: 10.1093/mnras/stad2889

  35. [35]

    Matsumoto, R., Kato, S., Fukue, J., & Okazaki, A. T. 1984, PASJ, 36, 71

  36. [36]

    2024, ApJ, 971, 28, doi: 10.3847/1538-4357/ad55cb

    Mitra, S., & Das, S. 2024, ApJ, 971, 28, doi: 10.3847/1538-4357/ad55cb

  37. [37]

    Nandi, A., Debnath, D., Mandal, S., & Chakrabarti, S. K. 2012, A&A, 542, A56, doi: 10.1051/0004-6361/201117844 16

  38. [38]

    2018, Ap&SS, 363, 90, doi: 10.1007/s10509-018-3314-1

    Nandi, A., Mandal, S., Sreehari, H., et al. 2018, Ap&SS, 363, 90, doi: 10.1007/s10509-018-3314-1

  39. [39]

    , eprint =

    Narayan, R., & Yi, I. 1994, ApJ, 428, L13, doi: 10.1086/187381 —. 1995a, ApJ, 452, 710, doi: 10.1086/176343 —. 1995b, ApJ, 444, 231, doi: 10.1086/175599

  40. [40]

    E., & Matsumoto, R

    Oda, H., Machida, M., Nakamura, K. E., & Matsumoto, R. 2007, PASJ, 59, 457, doi: 10.1093/pasj/59.2.457

  41. [41]

    2005, ApJ, 628, 368, doi: 10.1086/430728

    Ohsuga, K., Mori, M., Nakamoto, T., & Mineshige, S. 2005, ApJ, 628, 368, doi: 10.1086/430728

  42. [42]

    1997, MNRAS, 286, 681, doi: 10.1093/mnras/286.3.681

    Peitz, J., & Appl, S. 1997, MNRAS, 286, 681, doi: 10.1093/mnras/286.3.681

  43. [43]

    , keywords =

    Riffert, H., & Herold, H. 1995, ApJ, 450, 508, doi: 10.1086/176161

  44. [44]

    2022, Journal of Astrophysics and Astronomy, 43, 20, doi: 10.1007/s12036-022-09802-1

    Saedi, K., & Faghei, K. 2022, Journal of Astrophysics and Astronomy, 43, 20, doi: 10.1007/s12036-022-09802-1

  45. [45]

    2016, MNRAS, 461, 190, doi: 10.1093/mnras/stw1327

    Sarkar, B., & Das, S. 2016, MNRAS, 461, 190, doi: 10.1093/mnras/stw1327

  46. [46]

    2018, MNRAS, 473, 2415, doi: 10.1093/mnras/stx2505

    Sarkar, B., Das, S., & Mandal, S. 2018, MNRAS, 473, 2415, doi: 10.1093/mnras/stx2505

  47. [47]

    2022, Journal of Astrophysics and Astronomy, 43, 34, doi: 10.1007/s12036-022-09820-z

    Sarkar, S., & Chattopadhyay, I. 2022, Journal of Astrophysics and Astronomy, 43, 34, doi: 10.1007/s12036-022-09820-z

  48. [48]

    Sen, G., Maity, D., & Das, S. 2022, J. Cosmology Astropart. Phys., 2022, 048, doi: 10.1088/1475-7516/2022/08/048

  49. [49]

    2008, Ap&SS, 317, 201, doi: 10.1007/s10509-008-9882-8

    Shadmehri, M. 2008, Ap&SS, 317, 201, doi: 10.1007/s10509-008-9882-8

  50. [50]

    I., & Sunyaev, R

    Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337

  51. [51]

    L., & Teukolsky, S

    Shapiro, S. L., & Teukolsky, S. A. 1983, Black holes, white dwarfs, and neutron stars : the physics of compact objects

  52. [52]

    2024a, Ap&SS, 369, 1, doi: 10.1007/s10509-023-04263-6 —

    Singh, M., & Das, S. 2024a, Ap&SS, 369, 1, doi: 10.1007/s10509-023-04263-6 —. 2024b, arXiv e-prints, arXiv:2408.02256, doi: 10.48550/arXiv.2408.02256

  53. [53]

    Singh, M., Jana, C., & Das, S. 2025, J. Cosmology Astropart. Phys., 2025, 055, doi: 10.1088/1475-7516/2025/05/055 S ,adowski, A., Narayan, R., Tchekhovskoy, A., & Zhu, Y. 2013, MNRAS, 429, 3533, doi: 10.1093/mnras/sts632

  54. [54]

    2001, MNRAS, 322, 231, doi: 10.1046/j.1365-8711.2001.04022.x

    Stone, J. M., & Pringle, J. E. 2001, MNRAS, 322, 461, doi: 10.1046/j.1365-8711.2001.04138.x

  55. [55]

    , keywords =

    Stone, J. M., Pringle, J. E., & Begelman, M. C. 1999, MNRAS, 310, 1002, doi: 10.1046/j.1365-8711.1999.03024.x

  56. [56]

    2006, ApJ, 645, 1408, doi: 10.1086/500107

    Tsuruta, S. 2006, ApJ, 645, 1408, doi: 10.1086/500107

  57. [57]

    2002, ApJ, 572, 950, doi: 10.1086/340380

    Takahashi, M., Rilett, D., Fukumura, K., & Tsuruta, S. 2002, ApJ, 572, 950, doi: 10.1086/340380

  58. [58]

    2010, ApJ, 714, L176, doi: 10.1088/2041-8205/714/1/L176

    Takahashi, M., & Takahashi, R. 2010, ApJ, 714, L176, doi: 10.1088/2041-8205/714/1/L176

  59. [59]

    Tamilan, M., Hayasaki, K., & Suzuki, T. K. 2025, Progress of Theoretical and Experimental Physics, 2025, 023E02, doi: 10.1093/ptep/ptaf014

  60. [60]

    2008, ApJ, 681, 499, doi: 10.1086/588522

    Xie, F.-G., & Yuan, F. 2008, ApJ, 681, 499, doi: 10.1086/588522

  61. [61]

    1997, ApJ, 489, L29, doi: 10.1086/310956

    Xu, G., & Chen, X. 1997, ApJ, 489, L29, doi: 10.1086/310956

  62. [62]

    2012a, ApJ, 761, 130, doi: 10.1088/0004-637X/761/2/130

    Yuan, F., Bu, D., & Wu, M. 2012a, ApJ, 761, 130, doi: 10.1088/0004-637X/761/2/130

  63. [63]

    Hot Accretion Flows Around Black Holes

    Yuan, F., & Narayan, R. 2014, ARA&A, 52, 529, doi: 10.1146/annurev-astro-082812-141003

  64. [64]

    2012b, ApJ, 761, 129, doi: 10.1088/0004-637X/761/2/129

    Yuan, F., Wu, M., & Bu, D. 2012b, ApJ, 761, 129, doi: 10.1088/0004-637X/761/2/129 Zahra Zeraatgari, F., Mosallanezhad, A., Yuan, Y.-F., Bu, D.-F., & Mei, L. 2020, ApJ, 888, 86, doi: 10.3847/1538-4357/ab594f 17 APPENDIX DERIVATION OF FLOW VARIABLE GRADIENTS Using equation (6) in equations (7), (9), (12) and (16), we get Ev dv dx +E cs dCs dx +E λ dλ dx +E ...