pith. sign in

arxiv: 2604.23123 · v1 · submitted 2026-04-25 · 🌌 astro-ph.GA

The Dependence of the Mean Spectral Energy Distributions on the Accretion Rate for Quasars with z < 0.75 from the Sloan Digital Sky Survey

Yan-Song Ma , Yu-Meng Guan , Jian-Xia Jiang , Shao-Jun Li , Xiang-Wei Ning , Yi Tang , Wei-Hao Bian (NJNU) This is my paper

Pith reviewed 2026-05-08 07:46 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords quasarsspectral energy distributionsaccretion rateEddington ratioType 1 AGNFe II strengthH beta widthSloan Digital Sky Survey
0
0 comments X

The pith

Quasar spectral energy distributions redden with rising accretion rate in UV, NIR and MIR, but blue in the optical.

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

The paper assembles average spectral energy distributions from mid-infrared to ultraviolet for 56,969 quasars at redshifts below 0.75 using Sloan Digital Sky Survey data. It bins these averages by the relative strength of optical iron lines, the width of the H-beta emission line, and directly by the Eddington ratio as a measure of accretion rate. Higher accretion rates produce redder continua in the ultraviolet, near-infrared, and mid-infrared, consistent with stronger dust emission, while the optical continuum becomes harder and bluer. The result establishes that the overall shape of Type 1 active galactic nucleus spectra changes systematically with accretion rate.

Core claim

Mean SEDs binned by eigenvector 1 parameters and by Eddington ratio show that Type 1 AGN continua become redder in the MIR, NIR, and UV with increasing accretion rate, indicating more dust emission, while the optical continuum shows the opposite trend of becoming harder and bluer.

What carries the argument

Mean spectral energy distributions built by binning on R_FeII, H-beta line width, and Eddington ratio to isolate accretion-rate effects.

If this is right

  • Quasars with higher Eddington ratios exhibit increased dust reprocessing in the infrared.
  • The optical emission region responds differently to accretion rate than the UV and infrared regions.
  • Larger H-beta widths correspond to bluer MIR continua, consistent with a larger viewing angle to the torus.
  • Eigenvector 1 parameters organize not only optical spectra but also the full multi-wavelength SED shape.

Where Pith is reading between the lines

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

  • These trends could allow SED shape to serve as an auxiliary indicator of accretion rate in samples lacking full spectroscopic coverage.
  • The optical bluing may reflect changes in the inner accretion disk structure that are not captured by standard thin-disk models.
  • Population studies of black hole growth should incorporate accretion-rate-dependent templates to avoid systematic offsets in luminosity estimates.

Load-bearing premise

Binning by R_FeII, H-beta width, and Eddington ratio cleanly isolates accretion-rate effects without residual selection biases or mismatches across the multi-wavelength data.

What would settle it

Repeating the mean SED construction on an independent quasar sample with uniform multi-band photometry and no dependence on the same optical line measurements would yield no trend with Eddington ratio.

Figures

Figures reproduced from arXiv: 2604.23123 by Jian-Xia Jiang, Shao-Jun Li, Wei-Hao Bian (NJNU), Xiang-Wei Ning, Yan-Song Ma, Yi Tang, Yu-Meng Guan.

Figure 1
Figure 1. Figure 1: Mean SED (black) and median SED (orange)—both constructed with gap-repaired photometry—together with the underlying data points (without gap repair) from 749,436 SDSS DR16 quasars used to derive them. The dotted lines are 1 σ levels. The colors from left to right denote data from WISE W4, W3, W2, and W1, 2MASS KHJ, SDSS zirgu, and GALEX NUV and FUV. The mean SEDs from Richards et al. (2006a) (teal represen… view at source ↗
Figure 2
Figure 2. Figure 2: Top: mean corrected SED (gray) , median SED (red) and data points (without gap repair) for 56,969 quasars at z < 0.75 and with 0 ≤ RFe ii ≤ 3. The colors are the same as view at source ↗
Figure 3
Figure 3. Figure 3: The RFe ii dependence of the mean SEDs for three RFe ii bins, i.e., [0, 0.48], (0.48, 0.87], (0.87, 3] with the same sample sizes of about 18,990 (orange, green, and blue lines, respectively). At the bottom of the figure, the four slopes in the MIR, NIR, optical and UV are shown for the three RFe ii bins. The mean SED with the largest RFe ii bin shows the reddest UV, NIR, MIR, and optical continuum. Black … view at source ↗
Figure 4
Figure 4. Figure 4: Top: the LBol/LEdd dependence of the mean SEDs for three LBol/LEdd bins, i.e., (-3.38,-1.26], (-1.26,-0.76], and (-0.76, 1.03] with the same sample sizes of 18,990 (orange, green, and blue lines, respectively). In the bottom, the four slopes in the MIR, NIR, optical, and UV are shown for the three LBol/LEdd bins. The SED of the bin with the largest LBol/LEdd shows the reddest UV, NIR, and MIR continua and … view at source ↗
Figure 5
Figure 5. Figure 5: The FWHM dependence of the mean SEDs for four FWHMHβ bins, i.e., (0, 3000], (3000,5000], (5000, 8000], and [8000-) (km/s) with sample sizes of 18,659, 20,060, 13,312, 4,938 (orange, green, blue, and red lines, respectively). In the bottom, the four slopes in the MIR, NIR, optical and UV are shown for the four FWHMHβ bins. The SED of the bin with the smaller FWHMHβ shows bluer optical and NIR continua and r… view at source ↗
Figure 6
Figure 6. Figure 6: Top: BCs from the integration limits of 1 µm−2 keV as a function of frequency for the RFe ii-binned SEDs in view at source ↗
Figure 8
Figure 8. Figure 8: Relation between BC5100 and L5100. The blue line shows the luminosity-independent BC that assuming a linear dependence BC5100 = 53 − log(L5100) (Netzer 2013). Grey points show every source in DR16, black dashed curves are density contours, and red diamonds mark the medians. M˙ with rs = 0.04, 0.09, and 0.07, respectively, and with pnull < 10−6 . 4. CONCLUSIONS For a large sample of 56,969 SDSS DR16 quasars… view at source ↗
Figure 9
Figure 9. Figure 9: Comparison of the BC5100 and L5100 relation between DR7 and DR16 high-luminosity samples. The blue line shows BC5100 = 53 − log(L5100) (Netzer 2013). Red diamonds indicate the median values for the DR16 high-- luminosity sample, while green diamonds show the median values for the DR7 sample. The contour lines represent the source density. • For three RFe ii-binned SEDs, it is found that the UV slope depend… view at source ↗
read the original abstract

We construct mean spectral energy distributions (SEDs) for a substantial sample of 56,969 Sloan Digital Sky Survey DR16 quasars with $z < 0.75$, utilizing multiwavelength data from the mid-infrared (MIR) to ultraviolet (UV). These SEDs are built on eigenvector 1 parameters -- the relative optical $\rm Fe~ II$ strength ($R_{\rm Fe~II}$) and the H$\beta$ line width ($\rm H\beta$) -- that capture the principal spectral variance of quasar spectra. From three $R_{\rm Fe~II}$-dependent mean SEDs we find that quasars with a larger $R_{\rm Fe~II}$ exhibit redder UV and optical and redder MIR and near-infrared (NIR) continua, indicating more dust emission. We also split our sample directly into Eddington ratio $L_{\rm Bol} /L_{\rm Edd}$ (or dimensionless accretion rate $\dot{\mathscr{M}}$) bins to construct different mean SEDs and find that the continua become increasingly red with increasing $L_{\rm Bol} /L_{\rm Edd}$ (or $\dot{\mathscr{M}}$) in the MIR, NIR, and UV bands. This demonstrates that the shapes of Type 1 AGN SEDs depend on the accretion rate. However, the optical continuum shows the opposite trend (becoming harder and bluer), indicating the complexity of the optical emission region. From $\rm FWHM_{H\beta}$-dependent mean SEDs we find that quasars with a larger $\rm FWHM_{H\beta}$ show redder optical and NIR continua and bluer UV and MIR continua. The bluer MIR continuum suggests that a larger angle between of the line of sight and the torus plane results in weaker torus emission in the MIR.

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

2 major / 2 minor

Summary. The paper constructs mean multi-wavelength SEDs (MIR to UV) for 56,969 SDSS DR16 quasars at z < 0.75. It bins the sample on eigenvector-1 parameters (R_FeII and Hβ FWHM) and directly on Eddington ratio (L_bol/L_edd or dimensionless accretion rate), reporting that higher R_FeII and higher accretion rate produce redder MIR/NIR/UV continua while the optical continuum becomes bluer with increasing L_bol/L_edd; the FWHM_Hβ split shows redder optical/NIR but bluer UV/MIR. The central claim is that Type 1 AGN SED shapes depend on accretion rate.

Significance. If the trends survive checks for selection and coverage biases, the large-sample observational result would provide direct evidence that accretion rate modulates quasar continuum shapes, particularly dust-related emission, with implications for AGN structure and accretion-disk models. The direct binning by L_bol/L_edd (rather than only eigenvector-1 proxies) is a methodological strength, but the significance is limited by the absence of quantified robustness tests against the noted data-coverage and parameter-overlap issues.

major comments (2)
  1. [Methods (Eddington-ratio calculation and sample binning)] The Eddington ratio used for direct binning is computed from Hβ FWHM and continuum luminosity (plus bolometric corrections), quantities that overlap with the eigenvector-1 parameters and the spectral regions averaged into the SEDs. This creates a risk that the reported redder MIR/NIR/UV trends at high L_bol/L_edd partly reflect the construction of the binning variable rather than an independent accretion-rate effect. A test that recomputes L_bol/L_edd with an independent luminosity indicator (e.g., 5100 Å luminosity only or X-ray) and re-derives the mean SEDs is required to establish that the dependence is not circular.
  2. [Methods (multi-wavelength photometry matching and mean-SED construction)] The abstract and methods provide no description of how incomplete multi-wavelength coverage (WISE MIR, GALEX UV, etc.) is treated when stacking mean SEDs, nor whether the accretion-rate bins are matched in redshift or luminosity. Because detection fractions can correlate with luminosity, dust content, or redshift—all of which may differ systematically across L_bol/L_edd bins—the observed reddening in MIR/NIR/UV could partly arise from differential data availability rather than intrinsic SED changes. Explicit completeness fractions per bin and a luminosity/redshift-matched control sample are needed.
minor comments (2)
  1. [Abstract] The abstract states the sample size but omits any mention of error bars on the mean SEDs, bootstrap uncertainties, or the number of objects contributing to each wavelength bin.
  2. [Results and discussion] The opposite trend in the optical continuum (bluer at higher accretion rate) is noted but receives less quantitative discussion than the MIR/NIR/UV reddening; a short comparison to thin-disk models or prior optical-slope studies would clarify the claimed complexity.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for their careful and constructive review of our manuscript. We address each major comment below and have revised the paper where needed to strengthen the analysis and presentation.

read point-by-point responses

  1. Referee: [Methods (Eddington-ratio calculation and sample binning)] The Eddington ratio used for direct binning is computed from Hβ FWHM and continuum luminosity (plus bolometric corrections), quantities that overlap with the eigenvector-1 parameters and the spectral regions averaged into the SEDs. This creates a risk that the reported redder MIR/NIR/UV trends at high L_bol/L_edd partly reflect the construction of the binning variable rather than an independent accretion-rate effect. A test that recomputes L_bol/L_edd with an independent luminosity indicator (e.g., 5100 Å luminosity only or X-ray) and re-derives the mean SEDs is required to establish that the dependence is not circular.

    Authors: We agree there is parameter overlap, as L_bol/L_edd incorporates Hβ FWHM and 5100 Å continuum luminosity. However, the MIR, NIR, and UV portions of the mean SEDs are built from independent WISE and GALEX photometry not used in the binning. The opposite trend (optical continuum blueing with rising L_bol/L_edd) is inconsistent with a purely circular artifact. We have added a dedicated paragraph in the Methods section explaining this independence and performed a supplementary test recomputing L_bol/L_edd using only the 5100 Å luminosity with a fixed bolometric correction factor; the reported trends remain unchanged. We note that X-ray data are unavailable for the full sample. revision: partial

  2. Referee: [Methods (multi-wavelength photometry matching and mean-SED construction)] The abstract and methods provide no description of how incomplete multi-wavelength coverage (WISE MIR, GALEX UV, etc.) is treated when stacking mean SEDs, nor whether the accretion-rate bins are matched in redshift or luminosity. Because detection fractions can correlate with luminosity, dust content, or redshift—all of which may differ systematically across L_bol/L_edd bins—the observed reddening in MIR/NIR/UV could partly arise from differential data availability rather than intrinsic SED changes. Explicit completeness fractions per bin and a luminosity/redshift-matched control sample are needed.

    Authors: We have revised the Methods section to explicitly describe the construction: mean SEDs are formed band-by-band, including only objects with available photometry in each band, and we now tabulate completeness fractions for every Eddington-ratio bin and wavelength range. We have also constructed a control subsample matched in both redshift and bolometric luminosity across the L_bol/L_edd bins and repeated the stacking; the MIR/NIR/UV reddening and optical blueing trends persist at the same significance. These details and the matched-sample results have been added to the text, and the abstract has been updated to reference the completeness treatment. revision: yes

standing simulated objections not resolved
  • X-ray luminosity is unavailable for the full sample, so a complete test using X-ray-based Eddington ratios cannot be performed.

Circularity Check

0 steps flagged

No significant circularity; purely observational binning of public survey data

full rationale

The paper constructs mean SEDs by averaging multi-wavelength photometry for quasars binned according to eigenvector-1 parameters (R_FeII, Hβ FWHM) and separately by Eddington ratio. No derivation, equation, or central claim reduces by construction to a fitted parameter, self-referential definition, or load-bearing self-citation chain. The analysis relies on direct empirical stacking of external survey data (SDSS, WISE, GALEX, etc.) without any predictive step that is statistically forced by the binning inputs themselves. This is the most common honest outcome for large-sample observational papers.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Work rests on standard domain assumptions in quasar spectroscopy and public survey data; no new entities or heavy free parameters introduced beyond binning choices.

axioms (1)
  • domain assumption Eigenvector 1 parameters (R_FeII and H-beta width) capture the principal spectral variance of quasar spectra
    Invoked to construct the three R_FeII-dependent mean SEDs.

pith-pipeline@v0.9.0 · 5675 in / 1203 out tokens · 40442 ms · 2026-05-08T07:46:48.584010+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

66 extracted references · 66 canonical work pages · 1 internal anchor

  1. [1]

    Obscured at the Core: Evidence for Nuclear Dust in Reddened Type-1 AGN

    thebibliography [1] 20pt to REFERENCES 6pt =0pt \@twocolumntrue 12pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key o...

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2604.13031 2017

  2. [2]

    Antonucci R., 1993, , 31, 473

    work page 1993

  3. [3]

    Azadi M., et al., 2022, arXiv:2204.12697

    work page arXiv 2022

  4. [4]

    H., White R

    Becker R. H., White R. L., & Helfand D. J., 1995, ApJ, 450, 559

    work page 1995

  5. [5]

    C., et al., 2013, , 767, 149

    Bentz M. C., et al., 2013, , 767, 149

    work page 2013

  6. [6]

    H., & Zhao Y

    Bian W. H., & Zhao Y. H., 2003, MNRAS, 343, 164

    work page 2003

  7. [7]

    J., Trumper J., et al., 2016, A&A, 588, A103

    Boller T., Freyberg M. J., Trumper J., et al., 2016, A&A, 588, A103

    work page 2016

  8. [8]

    A., & Green R

    Boroson T. A., & Green R. F., 1992, , 80, 109

    work page 1992

  9. [9]

    A., 2002, ApJ, 565, 78

    Boroson T. A., 2002, ApJ, 565, 78

    work page 2002

  10. [10]

    Q., Liu Y

    Chen Y. Q., Liu Y. S., Bian W. H., 2022, ApJ, 940, 50

    work page 2022

  11. [11]

    et al., 2019, ApJ, 870, 84

    Czerny B. et al., 2019, ApJ, 870, 84

    work page 2019

  12. [12]

    Du P., et al., 2016, , 825, 126

    work page 2016

  13. [13]

    Du P., Wang J.-M., 2019, ApJ, 886, 42

    work page 2019

  14. [14]

    J., McDowell J

    Elvis M., Wilkes B. J., McDowell J. C., et al., 1994, , 95, 1

    work page 1994

  15. [15]

    L., 1999, , 111, 63

    Fitzpatrick E. L., 1999, , 111, 63

    work page 1999

  16. [16]

    Fioc M., Rocca-Volmerange B., 1997, A&A, 326, 950

    work page 1997

  17. [17]

    E., et al

    Fukugita M., Ichikawa T., Gunn J. E., et al. 1996, , 111, 1748

    work page 1996

  18. [18]

    H., Jiang X

    Ge X., Bian W. H., Jiang X. L., et al., 2016, MNRAS, 462, 966

    work page 2016

  19. [19]

    Hagen S., Done C., 2023, , arXiv: 2304.01253

    work page arXiv 2023

  20. [20]

    C., 2008, , 46, 475

    Ho L. C., 2008, , 46, 475

    work page 2008

  21. [21]

    Jin C., Ward M., Done C., Gelbord J., 2012, MNRAS, 420, 1825

    work page 2012

  22. [22]

    Kormendy J., & Ho L. C. 2013, , 51, 511

    work page 2013

  23. [23]

    Kaspi S., et al., 2000, , 533, 631

    work page 2000

  24. [24]

    M., et al., 2013, ApJS, 206, 4

    Krawczyk C. M., et al., 2013, ApJS, 206, 4

    work page 2013

  25. [25]

    Kubota A., Done C., 2018, MNRAS, 480, 1247

    work page 2018

  26. [26]

    Gaia Collaboration, Brown A. G. A., Vallenari A., et al., 2018, A&A, 616, A1

    work page 2018

  27. [27]

    M., Turner T

    George I. M., Turner T. J., Yaqoob T., et al., 2000, ApJ, 531, 52

    work page 2000

  28. [28]

    Glikman E., et al., 2022, , 934, 119

    work page 2022

  29. [29]

    N., et al., 2021, ApJ, 910, 103

    Liu H., Luo B., Brandt W. N., et al., 2021, ApJ, 910, 103

    work page 2021

  30. [30]

    Liu Y. S. , Bian W. H., 2022, , 937, 82

    work page 2022

  31. [31]

    W., Higley A

    Lyke B. W., Higley A. N., McLane J. N., et al., 2020, ApJS, 250, 8

    work page 2020

  32. [32]

    Marconi A., Risaliti G., Gilli R., et al., 2004, MNRAS, 351, 169

    work page 2004

  33. [33]

    Meiksin A., 2006, MNRAS, 365, 807

    work page 2006

  34. [34]

    C., Fanson J., Schiminovich D., et al., 2005, ApJL, 619, L1

    Martin D. C., Fanson J., Schiminovich D., et al., 2005, ApJL, 619, L1

    work page 2005

  35. [35]

    J., Almaini O., et al., 2007, MNRAS, 379, 1599

    Lawrence A., Warren S. J., Almaini O., et al., 2007, MNRAS, 379, 1599

    work page 2007

  36. [36]

    Molina J.,et al., 2023, , arXiv: 2304.01017

    work page arXiv 2023

  37. [37]

    Netzer H., 2013, The Physics and Evolution of Active Galactic Nuclei (Cambridge: Cambridge Univ. Press)

    work page 2013

  38. [38]

    B., Sandage A., 1968, ApJ, 154, 21

    Oke J. B., Sandage A., 1968, ApJ, 154, 21

    work page 1968

  39. [39]

    M., et al., 2004, , 613, 682

    Peterson B. M., et al., 2004, , 613, 682

    work page 2004

  40. [40]

    T., et al., 2002, AJ, 123, 2945

    Richards G. T., et al., 2002, AJ, 123, 2945

    work page 2002

  41. [41]

    T., Lacy M., Storrie-Lombardi L

    Richards G. T., Lacy M., Storrie-Lombardi L. J., et al., 2006a, ApJS, 166, 470

    work page

  42. [42]

    T., Strauss M

    Richards G. T., Strauss M. A., Fan X., et al., 2006b, AJ, 131, 2766

    work page

  43. [43]

    C., et al., 2021, , 649, A102

    Rivera G. C., et al., 2021, , 649, A102

    work page 2021

  44. [44]

    R., Webb N

    Rosen S. R., Webb N. A., Watson M. G., et al., 2016, A&A, 590, A1

    work page 2016

  45. [45]

    Siebenmorgen R., Heymann F., Efstathiou A., 2015, , 583, A120

    work page 2015

  46. [46]

    J., Finkbeiner D

    Schlegel D. J., Finkbeiner D. P., & Davis M., 1998, ApJ, 500, 525

    work page 1998

  47. [47]

    P., et al., 2010, , 139, 2360

    Schneider D. P., et al., 2010, , 139, 2360

    work page 2010

  48. [48]

    I., Sunyaev R

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

    work page 1973

  49. [49]

    Brotherton M

    Shang Z. Brotherton M. S., Wills B. J., et al., 2011, ApJS, 196, 2

    work page 2011

  50. [50]

    2011, , 194, 45

    Shen Y., et al. 2011, , 194, 45

    work page 2011

  51. [51]

    Wu, Q., Shen, Y., 2022, , 263, 2, 42

    work page 2022

  52. [52]

    C., 2014, , 513, 210

    Shen Y., Ho L. C., 2014, , 513, 210

    work page 2014

  53. [53]

    H., Ogle P

    Shi Y., Rieke G. H., Ogle P. M., Su K. Y. L., Balog Z., 2014, ApJS, 214, 23

    work page 2014

  54. [54]

    F., Cutri R

    Skrutskie M. F., Cutri R. M., Stiening R., et al. 2006, AJ, 131, 1163

    work page 2006

  55. [55]

    T., Strateva, I., Brandt, W

    Steffen, A. T., Strateva, I., Brandt, W. N., et al. 2006, AJ, 131, 2826

    work page 2006

  56. [56]

    Stolc M., et al., 2023, arXiv: 2304.03015

    work page arXiv 2023

  57. [57]

    W., Marziani P., Dultzin-Hacyan D

    Sulentic J. W., Marziani P., Dultzin-Hacyan D. 2000, , 38, 521

    work page 2000

  58. [58]

    M., Zhao B

    Yu L. M., Zhao B. X., Bian W. H., et al., 2020, , 491, 5881

    work page 2020

  59. [59]

    E., Richards G

    Vanden Berk D. E., Richards G. T., Bauer A., et al., 2001, AJ, 122, 549

    work page 2001

  60. [60]

    Wang H., Zhou H., Yuan W., Wang T., 2012, ApJL, 751, L23

    work page 2012

  61. [61]

    M., Du P., Li Y.-R., Ho L

    Wang J. M., Du P., Li Y.-R., Ho L. C., Hu C., Bai J.-M., 2014, , 792, L13

    work page 2014

  62. [62]

    L., Eisenhardt P

    Wright E. L., Eisenhardt P. R. M., Mainzer A. K., et al., 2010, AJ, 140, 1868

    work page 2010

  63. [63]

    K., et al., 2005, ApJ, 619, L15

    Wyder T. K., et al., 2005, ApJ, 619, L15

    work page 2005

  64. [64]

    A., Johnson W

    Zdziarski A. A., Johnson W. N. Done C. Smith D., & McNaron-Brown K., 1995, ApJL, 438, L63

    work page 1995

  65. [65]

    H., et al., 2014, , 786, 42

    Zhang S. H., et al., 2014, , 786, 42

    work page 2014

  66. [66]

    C., Shangguan J., 2018, ApJ, 862, 118

    Zhuang M.-Y., Ho L. C., Shangguan J., 2018, ApJ, 862, 118

    work page 2018