pith. machine review for the scientific record. sign in

arxiv: 2605.13990 · v1 · submitted 2026-05-13 · 🌌 astro-ph.HE · astro-ph.GA

Recognition: no theorem link

Infrared Line Diagnostics Fail to Constrain Sgr A*'s UV Output

Authors on Pith no claims yet

Pith reviewed 2026-05-15 02:15 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.GA
keywords Sgr A*mid-infrared linesflaresphotoionizationGalactic CenterUV fluxvariability
0
0 comments X

The pith

Infrared emission lines near Sgr A* remain steady despite flares, preventing constraints on the black hole's instantaneous UV output.

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

The paper tests whether mid-infrared emission lines can reveal the ultraviolet output of Sgr A* during its X-ray flares by taking time-resolved JWST spectroscopy of the central 0.3 arcsec region. No statistically significant changes appear in lines such as [Fe II] 5.34 micron, [Ne II] 12.813 micron, [Fe II] 17.936 micron, or [S III] 18.713 micron. CLOUDY photoionization models for flare UV luminosities from 10^32 to 10^39 erg s^-1 predict clear responses in a dense medium, yet the observed constancy matches expectations once light-crossing times of 0.1-10 days and recombination/cooling times longer than flare durations are included. The lines therefore register only the time-averaged radiation field.

Core claim

The absence of detectable variability in mid-infrared emission lines, despite expectations from photoionization models, arises because the spatially extended line-emitting gas has light-crossing timescales of 0.1-10 days and recombination/cooling timescales much longer than individual flare durations; the resulting emission is continuum-dominated, intrinsically weak, and velocity-broadened to ~10^3 km s^-1, so infrared line diagnostics cannot constrain the instantaneous UV flux of Sgr A*.

What carries the argument

The mismatch between short flare durations and the long light-crossing plus recombination timescales of the extended line-emitting gas, which averages any transient response into a steady-state signal.

If this is right

  • The predicted line emission is continuum-dominated and broadened by velocities of order 10^3 km s^-1, which reduces contrast and limits detectability.
  • Extending the search to higher-ionization mid-infrared or near-infrared lines yields no gain in sensitivity.
  • Infrared line diagnostics can only trace the steady-state radiation field and are unsuitable for measuring flare-by-flare UV output.

Where Pith is reading between the lines

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

  • Steady-state photoionization calculations could still be used to bound the time-averaged UV luminosity over weeks or longer.
  • Similar averaging arguments may apply to line diagnostics around other accreting black holes with short variability timescales.
  • Direct constraints on Sgr A* UV output will require wavelengths or techniques that respond on the flare timescale itself.

Load-bearing premise

The line-emitting gas is distributed over scales where light takes days to cross and where recombination and cooling times greatly exceed flare lengths.

What would settle it

Detection of statistically significant changes in any of the listed mid-infrared lines that are synchronized with known Sgr A* flare times on hourly-to-daily scales would contradict the central claim.

Figures

Figures reproduced from arXiv: 2605.13990 by Daryl Haggard, Giacomo Principe, Howard A. Smith, Joey Neilsen, Joseph L. Hora, Joseph M. Michail, Laurent Loinard, Mayura Balakrishnan, Nadeen B. Sabha, Nicole M. Ford, Sebastiano D. von Fellenberg, Sera Markoff, Shuo Zhang, Zach Sumners.

Figure 1
Figure 1. Figure 1: JWST/MIRI Channel 1 5.6 micron image of the Galactic Center field. Prominent sources are labeled in black. The position of Sgr A⋆ is marked by the yellow “X”, and the white circle indicates the 0. ′′3 radius aperture used for spectral extraction. The color scale shows the calibrated surface brightness in units of MJy sr−1 . et al. 2025a) and staring-mode aperture corrections (see Appendix A of J. M. Michai… view at source ↗
Figure 2
Figure 2. Figure 2: Representative lightcurves showing the absence of variability in both [S III] (18.713 µm) and [Ne II] (12.813 µm), shown in blue, in response to the Sgr A⋆ flare described in S. D. von Fellenberg et al. (2025b) and J. M. Michail et al. (2026). The corresponding continuum lightcurve is plotted in orange. Uncertainties on the integrated line flux are estimated via a Monte Carlo procedure in which the spectru… view at source ↗
Figure 3
Figure 3. Figure 3: Comparison of CLOUDY models to the observed JWST/MIRI spectrum of the central 0. ′′3 around Sgr A⋆ . The crimson and gray curve show the MIRI MRS spectra, dereddened and observed, respectively, during the flare window and the dashed crimson curve shows the local continuum-subtracted RMS. Colored curves show the predicted flare-peak spectra for flare luminosities Lflare = 1032–1039 erg s−1 . The black curve… view at source ↗
Figure 4
Figure 4. Figure 4: Predicted time-dependent line fluxes for representative high-ionization MIR (left) and NIR (right) transitions at log Lflare = 39, computed with CLOUDY. The flare is injected for the first hour of the simulation. Fluxes are shown as equivalent Gaussian peak flux densities assuming a FWHM of 1000 km s−1 . In all cases, the emission is extremely faint (≲ 10−5 mJy) and exhibits only weak temporal evolution fo… view at source ↗
Figure 5
Figure 5. Figure 5: Recombination time (blue points), cooling time (orange squares; computed assuming a representative cooling coefficient of 10−22 erg cm3 s −1 ), and density (green diamonds) as a function of the radius at which each ion reaches peak emissivity. Ion labels are shown above the recombination-time points. Even the shortest recombination timescale is of order ∼ 8 days, far exceeding the characteristic variabilit… view at source ↗
read the original abstract

Sgr A*, the 4 x 10^6 solar-mass supermassive black hole at the Galactic Center, exhibits frequent flaring with X-ray luminosities of L_X ~ 10^35--10^36 erg s^-1, while its ultraviolet (UV) emission remains unconstrained due to extreme extinction (A_V ~ 30 mag). We use JWST/MIRI time-resolved spectroscopy of the central Galactic Center's 0.3 arcsec region to search for mid-infrared emission-line variability driven by Sgr A* flares, comparing the results to CLOUDY photoionization models spanning flare luminosities of L_UV = 10^32--10^39 erg s^-1 in a dense medium. We detect no statistically significant variability in any mid-infrared line, including [Fe II] 5.34 micron, [Ne II] 12.813 micron, [Fe II] 17.936 micron, and [S III] 18.713 micron. Despite expectations of a flare-driven response, we show that the lack of variability is consistent with the physical conditions in the spatially extended line-emitting gas, where light-crossing timescales of ~0.1--10 days and recombination and cooling timescales much longer than the flare timescale suppress any observable response to individual flares. We further find that the predicted emission is continuum dominated and that even the brightest lines are intrinsically weak and broadened by velocities of order 10^3 km s^-1, reducing their contrast against the continuum and limiting their detectability. Extending the analysis to higher-ionization mid-infrared and near-infrared lines does not improve sensitivity. These results demonstrate that infrared emission lines trace a steady-state radiation field rather than individual flaring events, and therefore infrared line diagnostics cannot be used to constrain the instantaneous UV flux of Sgr A*.

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

0 major / 2 minor

Summary. The paper reports a non-detection of statistically significant variability in mid-infrared emission lines ([Fe II] 5.34 μm, [Ne II] 12.813 μm, [Fe II] 17.936 μm, [S III] 18.713 μm) from JWST/MIRI time-resolved spectroscopy of the central 0.3 arcsec region around Sgr A*. It compares this null result to CLOUDY photoionization models for L_UV spanning 10^32–10^39 erg s^{-1} and attributes the lack of response to light-crossing timescales of ~0.1–10 days across the extended gas plus recombination/cooling timescales exceeding flare durations, concluding that the lines trace a steady-state radiation field and cannot constrain the instantaneous UV flux of Sgr A*.

Significance. If the result holds, the work is significant for establishing that infrared line diagnostics are insensitive to individual flares from Sgr A* owing to the physical conditions in the spatially extended line-emitting gas. The combination of the reported null detection, explicit model comparison, and timescale argument provides a clear, falsifiable explanation for why such lines are continuum-dominated and velocity-broadened, limiting their utility for transient UV constraints. This is a useful observational limit for Galactic Center studies.

minor comments (2)
  1. The abstract states that lines are 'broadened by velocities of order 10^3 km s^{-1}'; the main text should explicitly state the adopted line profile (Gaussian or otherwise) and how this broadening is incorporated into the CLOUDY predictions for contrast against the continuum.
  2. The manuscript would benefit from a brief table summarizing the exact flare durations and luminosities used in the CLOUDY grid, to allow direct comparison with the observed non-variability limits.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive assessment of our work and for recommending acceptance. We are pleased that the combination of the null detection, model comparison, and timescale arguments is viewed as providing a clear explanation for the lack of observable line variability.

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper reports a direct observational null result from JWST/MIRI time-resolved spectroscopy showing no statistically significant variability in mid-IR lines. This is compared to independent CLOUDY photoionization models run over a range of L_UV values. The explanation for the null result invokes standard physical timescales (light-crossing time across 0.3 arcsec and recombination/cooling times) derived from atomic physics and geometry, without any parameter fitting to the target data that would make the conclusion self-referential. No self-citations, ansatzes, or renamings reduce the central claim to its own inputs by construction; the derivation remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The analysis rests on standard photoionization modeling and domain assumptions about gas geometry and timescales in the Galactic Center; no new free parameters or entities are introduced.

axioms (2)
  • domain assumption CLOUDY photoionization models accurately predict line responses to UV illumination in dense gas
    Invoked to interpret the expected flare-driven variability across L_UV = 10^32-10^39 erg s^-1.
  • domain assumption Light-crossing, recombination, and cooling timescales in the spatially extended line-emitting gas exceed flare durations
    Used to explain the absence of observable response.

pith-pipeline@v0.9.0 · 5699 in / 1302 out tokens · 49084 ms · 2026-05-15T02:15:28.959995+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

61 extracted references · 61 canonical work pages · 2 internal anchors

  1. [1]

    P., & Satyapal, S

    Abel, N. P., & Satyapal, S. 2008, ApJ, 678, 686, doi: 10.1086/529013 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33, doi: 10.1051/0004-6361/201322068 Astropy Collaboration, Price-Whelan, A. M., Sip˝ ocz, B. M., et al. 2018, AJ, 156, 123, doi: 10.3847/1538-3881/aabc4f Astropy Collaboration, Price-Whelan, A. M., Lim, P...

  2. [2]

    2024, ApJ, 974, 98, doi: 10.3847/1538-4357/ad6c08

    Balakrishnan, M., Corrales, L., Markoff, S., et al. 2024, ApJ, 974, 98, doi: 10.3847/1538-4357/ad6c08

  3. [3]

    E., Matthews, K., Neugebauer, G., & Willner, S

    Becklin, E. E., Matthews, K., Neugebauer, G., & Willner, S. P. 1978, ApJ, 219, 121, doi: 10.1086/155761

  4. [4]

    E., & Neugebauer, G

    Becklin, E. E., & Neugebauer, G. 1968, ApJ, 151, 145, doi: 10.1086/149425

  5. [5]

    S., & Wills, B

    Bennert, N., Falcke, H., Schulz, H., Wilson, A. S., & Wills, B. J. 2002, ApJL, 574, L105, doi: 10.1086/342420

  6. [6]

    C., & Katz, S

    Bentz, M. C., & Katz, S. 2015, PASP, 127, 67, doi: 10.1086/679601

  7. [7]

    C., Denney, K

    Bentz, M. C., Denney, K. D., Grier, C. J., et al. 2013, ApJ, 767, 149, doi: 10.1088/0004-637X/767/2/149

  8. [8]

    D., & McKee, C

    Blandford, R. D., & McKee, C. F. 1982, ApJ, 255, 419, doi: 10.1086/159843 Bouffard, ´E., Haggard, D., Nowak, M. A., et al. 2019, ApJ, 884, 148, doi: 10.3847/1538-4357/ab4266

  9. [9]

    2019, ApJ, 871, 161, doi: 10.3847/1538-4357/aaf71f

    Boyce, H., Haggard, D., Witzel, G., et al. 2019, ApJ, 871, 161, doi: 10.3847/1538-4357/aaf71f

  10. [10]

    2022, ApJ, 931, 7, doi: 10.3847/1538-4357/ac6104

    Boyce, H., Haggard, D., Witzel, G., et al. 2022, ApJ, 931, 7, doi: 10.3847/1538-4357/ac6104

  11. [11]

    2013, A&A, 558, A32, doi: 10.1051/0004-6361/201321667

    Clavel, M., Terrier, R., Goldwurm, A., et al. 2013, A&A, 558, A32, doi: 10.1051/0004-6361/201321667

  12. [12]

    R., Mon, B., Haggard, D., et al

    Corrales, L. R., Mon, B., Haggard, D., et al. 2017, ApJ, 839, 76, doi: 10.3847/1538-4357/aa68dd 13

  13. [13]

    Dempsey, R., & Zakamska, N. L. 2018, MNRAS, 477, 4615, doi: 10.1093/mnras/sty941

  14. [14]

    2013, ChiantiPy: Python package for the CHIANTI atomic database,, Astrophysics Source Code Library, record ascl:1308.017 http://ascl.net/1308.017

    Dere, K. 2013, ChiantiPy: Python package for the CHIANTI atomic database,, Astrophysics Source Code Library, record ascl:1308.017 http://ascl.net/1308.017

  15. [15]

    2026, Jdaviz, v4.5.1 Zenodo, doi: 10.5281/zenodo.18894838

    Developers, J., Averbukh, J., Bradley, L., et al. 2026, Jdaviz, v4.5.1 Zenodo, doi: 10.5281/zenodo.18894838

  16. [16]

    K., Ciurlo, A., Morris, M

    Dinh, C. K., Ciurlo, A., Morris, M. R., et al. 2024, AJ, 167, 41, doi: 10.3847/1538-3881/ad10a5

  17. [17]

    K., et al

    Do, T., Witzel, G., Gautam, A. K., et al. 2019, ApJL, 882, L27, doi: 10.3847/2041-8213/ab38c3

  18. [18]

    2009, ApJ, 698, 676, doi: 10.1088/0004-637X/698/1/676

    Dodds-Eden, K., Porquet, D., Trap, G., et al. 2009, ApJ, 698, 676, doi: 10.1088/0004-637X/698/1/676

  19. [19]

    Drappeau, S., Dibi, S., Dexter, J., Markoff, S., & Fragile, P. C. 2013, MNRAS, 431, 2872, doi: 10.1093/mnras/stt388

  20. [20]

    WHAT IS SGR A*? The Starved Black Hole in the Center of the Milky Way

    Falcke, H. 1996, in IAU Symposium, Vol. 169, Unsolved Problems of the Milky Way, ed. L. Blitz & P. J. Teuben, 169, doi: 10.48550/arXiv.astro-ph/9411065

  21. [21]

    M., Balakrishnan, M., von Fellenberg, S

    Ford, N. M., Balakrishnan, M., von Fellenberg, S. D., et al. 2026, ApJ

  22. [22]

    M., Salim, S., Weinberg, N

    Ghez, A. M., Salim, S., Weinberg, N. N., et al. 2008, ApJ, 689, 1044, doi: 10.1086/592738

  23. [23]

    M., Widmann, F., et al

    Gillessen, S., Plewa, P. M., Widmann, F., et al. 2019, ApJ, 871, 126, doi: 10.3847/1538-4357/aaf4f8 GRAVITY Collaboration, Abuter, R., Amorim, A., et al. 2019, A&A, 625, L10, doi: 10.1051/0004-6361/201935656 GRAVITY Collaboration, Abuter, R., Amorim, A., et al. 2021, A&A, 654, A22, doi: 10.1051/0004-6361/202140981

  24. [24]

    M., van Hoof, P

    Gunasekera, C. M., van Hoof, P. A. M., Dehghanian, M., et al. 2025, RMxAA, 61, 120, doi: 10.22201/ia.01851101p.2025.61.03.01

  25. [25]

    2019, ApJ, 886, 96, doi: 10.3847/1538-4357/ab4a7f 18

    Haggard, D., Nynka, M., Mon, B., et al. 2019, ApJ, 886, 96, doi: 10.3847/1538-4357/ab4a7f

  26. [26]

    2023, Spacetime & Spectra: Joint Chandra/JWST/EHT Observations of Sgr A*,, JWST Proposal

    Haggard, D., Balokovic, M., Bower, G., et al. 2023, Spacetime & Spectra: Joint Chandra/JWST/EHT Observations of Sgr A*,, JWST Proposal. Cycle 2, ID. #4572

  27. [27]

    J., Morgan, C

    Hainline, L. J., Morgan, C. W., MacLeod, C. L., et al. 2013, ApJ, 774, 69, doi: 10.1088/0004-637X/774/1/69

  28. [28]

    L., Haggard, D., Witzel, G., et al

    Hora, J. L., Haggard, D., Witzel, G., et al. 2023, Sgr A* as Particle Accelerator: What Drives the Black Hole’s Variable IR and X-ray Emission?,, JWST Proposal. Cycle 2, ID. #3324

  29. [29]

    L., Haggard, D., von Fellenberg, S., et al

    Hora, J. L., Haggard, D., von Fellenberg, S., et al. 2025, A Joint Mid-IR and X-ray Investigation of the Physics Driving Sgr A*’s Flares,, JWST Proposal. Cycle 4, ID. #7532

  30. [30]

    1996, PASJ, 48, 249, doi: 10.1093/pasj/48.2.249

    Koyama, K., Maeda, Y., Sonobe, T., et al. 1996, PASJ, 48, 249, doi: 10.1093/pasj/48.2.249

  31. [31]

    R., Do, T., Ghez, A

    Lu, J. R., Do, T., Ghez, A. M., et al. 2013, ApJ, 764, 155, doi: 10.1088/0004-637X/764/2/155

  32. [32]

    Markoff, S., Falcke, H., Yuan, F., & Biermann, P. L. 2001, A&A, 379, L13, doi: 10.1051/0004-6361:20011346

  33. [33]

    P., Moran, J

    Marrone, D. P., Moran, J. M., Zhao, J.-H., & Rao, R. 2006, in Journal of Physics Conference Series, Vol. 54, Journal of Physics Conference Series, ed. R. Sch¨ odel, G. C

  34. [34]

    Bower, M. P. Muno, S. Nayakshin, & T. Ott (IOP), 354–362, doi: 10.1088/1742-6596/54/1/056

  35. [35]

    P., Baganoff, F

    Marrone, D. P., Baganoff, F. K., Morris, M. R., et al. 2008, ApJ, 682, 373, doi: 10.1086/588806

  36. [36]

    J., et al

    Martins, F., Genzel, R., Hillier, D. J., et al. 2007, A&A, 468, 233, doi: 10.1051/0004-6361:20066688

  37. [37]

    2001, ARA&A, 39, 309, doi: 10.1146/annurev.astro.39.1.309

    Melia, F., & Falcke, H. 2001, ARA&A, 39, 309, doi: 10.1146/annurev.astro.39.1.309

  38. [38]

    M., von Fellenberg, S

    Michail, J. M., von Fellenberg, S. D., Keating, G. K., et al. 2026, ApJ, 997, 282, doi: 10.3847/1538-4357/ae25ef

  39. [39]

    Nahar, S. N. 1995, ApJS, 101, 423, doi: 10.1086/192248

  40. [40]

    E., Popham, R

    Narayan, R., Mahadevan, R., Grindlay, J. E., Popham, R. G., & Gammie, C. 1998, ApJ, 492, 554, doi: 10.1086/305070

  41. [41]

    , eprint =

    Narayan, R., & Yi, I. 1994, ApJL, 428, L13, doi: 10.1086/187381

  42. [42]

    1998, ApJS, 114, 269, doi: 10.1086/313069

    Nayakshin, S., & Melia, F. 1998, ApJS, 114, 269, doi: 10.1086/313069

  43. [43]

    A., Gammie, C., et al

    Neilsen, J., Nowak, M. A., Gammie, C., et al. 2013, ApJ, 774, 42, doi: 10.1088/0004-637X/774/1/42

  44. [44]

    2006, ApJ, 643, 1011, doi: 10.1086/503273

    Paumard, T., Genzel, R., Martins, F., et al. 2006, ApJ, 643, 1011, doi: 10.1086/503273

  45. [45]

    Peterson, B. M. 1997, An Introduction to Active Galactic Nuclei (Cambridge: Cambridge University Press)

  46. [46]

    M., Ferrarese, L., Gilbert, K

    Peterson, B. M., Ferrarese, L., Gilbert, K. M., et al. 2004, ApJ, 613, 682, doi: 10.1086/423269

  47. [47]

    2010, ApJ, 714, 732, doi: 10.1088/0004-637X/714/1/732

    Ponti, G., Terrier, R., Goldwurm, A., Belanger, G., & Trap, G. 2010, ApJ, 714, 732, doi: 10.1088/0004-637X/714/1/732

  48. [48]

    R., et al

    Ponti, G., De Marco, B., Morris, M. R., et al. 2015, MNRAS, 454, 1525, doi: 10.1093/mnras/stv1537

  49. [49]

    2003, A&A, 407, L17, doi: 10.1051/0004-6361:20030983

    Porquet, D., Predehl, P., Aschenbach, B., et al. 2003, A&A, 407, L17, doi: 10.1051/0004-6361:20030983

  50. [50]

    2002, ApJ, 575, 855, doi: 10.1086/341425

    Quataert, E. 2002, ApJ, 575, 855, doi: 10.1086/341425

  51. [51]

    M., Quataert, E., & Stone, J

    Ressler, S. M., Quataert, E., & Stone, J. M. 2018, MNRAS, 478, 3544, doi: 10.1093/mnras/sty1146

  52. [52]

    H., & Lebofsky, M

    Rieke, G. H., & Lebofsky, M. J. 1985, ApJ, 288, 618, doi: 10.1086/162827

  53. [53]

    A., et al

    Schirmer, M., Malhotra, S., Levenson, N. A., et al. 2016, MNRAS, 463, 1554, doi: 10.1093/mnras/stw1819

  54. [54]

    R., Donley, J

    Schmitt, H. R., Donley, J. L., Antonucci, R. R. J., et al. 2003, ApJ, 597, 768, doi: 10.1086/381224 14 Sch¨ odel, R., Morris, M. R., Muzic, K., et al. 2011, A&A, 532, A83, doi: 10.1051/0004-6361/201116994

  55. [55]

    1997, ApJL, 490, L77, doi: 10.1086/311010

    Serabyn, E., Carlstrom, J., Lay, O., et al. 1997, ApJL, 490, L77, doi: 10.1086/311010

  56. [56]

    R., & Loch, S

    Badnell, N. R., & Loch, S. D. 2007, in American Institute of Physics Conference Series, Vol. 901, Atomic and Molecular Data and Their Applications, ed. E. Roueff (AIP), 239–248, doi: 10.1063/1.2727374

  57. [57]

    M., Haggard, D., et al

    Sumners, Z., Ford, N. M., Haggard, D., et al. 2026, ApJ

  58. [58]

    E., Zakamska, N

    Sun, A.-L., Greene, J. E., Zakamska, N. L., et al. 2018, MNRAS, 480, 2302, doi: 10.1093/mnras/sty1394 The Event Horizon Telescope Collaboration. 2023, arXiv e-prints, arXiv:2311.08679, doi: 10.48550/arXiv.2311.08679 von Fellenberg, S. D., Michail, J. M., Willner, S. P., et al. 2025a, ApJ, 995, 215, doi: 10.3847/1538-4357/ae182a von Fellenberg, S. D., Royc...

  59. [59]

    D., Nowak, M

    Wang, Q. D., Nowak, M. A., Markoff, S. B., et al. 2013, Science, 341, 981, doi: 10.1126/science.1240755

  60. [60]

    2015, PASP, 127, 646, doi: 10.1086/682281

    Wells, M., Pel, J.-W., Glasse, A., et al. 2015, PASP, 127, 646, doi: 10.1086/682281

  61. [61]

    2003, ApJ, 598, 301, doi: 10.1086/378716

    Yuan, F., Quataert, E., & Narayan, R. 2003, ApJ, 598, 301, doi: 10.1086/378716