REVIEW 3 major objections 3 minor 49 references
Measuring Outflow Distances in NGC 5548 Using Absorption-Line Variability Diagnostics
T0 review · 3 major / 3 minor · reviewed 2026-08-01 · deepseek-v4-flash
Pith's one-line read Two ultraviolet outflow components in NGC 5548 are located 0.77 and about 1.7 parsecs from the active nucleus, measured through the timing of absorption-line variability rather than traditional density diagnostics.
desk verdict A useful application of the G1/G2 variability method to NGC 5548 with a genuine cross-check against Arav 2015, but component 1's distance is fragile: the DRW timescale systematic shifts it by ~2x and the robustness section contains an inverted scaling relation. 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 the G1 event and its probability F(G1). A G1 event is a triplet of spectra (with one shared epoch) in which the short-interval pair has stronger continuum variability than the long-interval pair, and the absorption-trough variability is also stronger in the short pair; G2 is the opposite outcome. Because gas with a long recombination timescale cannot fully respond within the short interval, F(G1) declines monotonically as t_r grows. The paper computes the mapping by generating damped-random-walk light curves of the ionizing continuum, matching each observed triplet's time intervals and flux changes, and then averaging the simulated continuum over a boxcar of length t_r
What would settle it
Run a dedicated high-cadence UV monitoring campaign on NGC 5548 that measures the actual cross-correlation lag between the far-UV continuum near 200 Å and the C IV absorption-trough variations; if the observed lag distribution is not consistent with a boxcar of length ~0.1 days (component 1) and ~1.8 days (component 6), the simulation mapping is falsified. Alternatively, directly measure the damping timescale of the ionizing continuum at ~200 Å; a value of 35 days instead of 10 would shift the component-1 recombination timescale from ~0.10 to ~0.37 days and its distance by a factor of ~3.7.
Extended reading notes
Core claim
Using two years of ultraviolet spectra of NGC 5548, the paper tracks variability in the C IV absorption troughs of two outflow components. It defines G1 events—triplets where the shorter time interval shows both stronger continuum variation and stronger absorption-trough variation than the longer interval—and measures the G1 fraction observationally: 86.3% for component 1 and 56.9% for component 6. Simulating the ionizing continuum as a damped random walk and letting the absorber respond by averaging flux over a window of length t_r, the authors build a mapping from G1 fraction to recombination timescale and read off t_r ≈ 10^{-0.98} days (0.10 days) and 10^{0.25} days (1.8 days). Combining
Load-bearing premise
Outflow distances here rest on the premise that absorbing gas responds to ionizing continuum changes by a simple boxcar average over a window of length t_r, and that the ~200 Å continuum is a damped random walk with a 10-day timescale (the geometric mean of 2.7 and 35 days, not directly measured); if either assumption fails, the inferred t_r and distances move by factors of a few.
Editorial extensions
If this is right
- Outflow distances in NGC 5548 can be measured from variability timing alone, giving 0.77 pc for component 1 and about 1.7 pc for component 6, consistent with independent methods.
- The G1-probability technique, previously applied to large quasar samples, works on a single, well-monitored AGN, so it can be used with existing and future intensive monitoring campaigns.
- For component 1, whose distance was previously only an upper limit, the method converts the limit into a direct estimate of about 0.1-day recombination timescale.
- The inferred distances imply that these UV outflows are launched within the host galaxy's immediate vicinity (sub-parsec to few-parsec scale), relevant for feedback energy budgets.
Reading between the lines
- Going beyond the paper: if the boxcar response were replaced with a photoionization response function, the inferred t_r values might shift; mock spectral time series could test this.
- Going beyond the paper: a direct measurement of the ~200 Å damping timescale—rather than the 10-day geometric mean—would settle the main calibration uncertainty; the paper's own tau=3 vs 35 day test moves t_r by up to ~0.4 dex.
- Going beyond the paper: applying the same diagnostics to other intensively monitored AGNs would convert a single-object demonstration into a sample, and could test whether outflow radius tracks Eddington ratio or black hole mass.
- Going beyond the paper: component 6's distance uncertainty reaches down to zero; a longer monitoring baseline or additional troughs would decide whether it is truly at 1.7 pc or merely an upper limit.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a method for measuring the radial distance of AGN outflow components in NGC 5548 using variability of C IV absorption troughs. The authors define G1/G2 events based on whether short- or long-interval spectral pairs show stronger absorption variability when the continuum varies more on the short interval. They calibrate the G1 fraction against the recombination timescale t_r using DRW-simulated light curves, assuming the absorber responds by boxcar-averaging the continuum over a window of length t_r. Applying this to troughs A and I, they infer t_r ≈ 10^{-0.98} d and 10^{0.25} d for components 1 and 6, and, using photoionization modeling and Eqs. (1)-(2), derive distances R = 0.77^{+0.10}_{-0.10} pc and R = 1.72^{+1.74}_{-1.72} pc. The component-1 result is claimed to improve on the earlier upper limit. The paper also tests sensitivity to the DRW damping timescale τ in §5.2 and concludes the distances are robust.
Significance. If the central results hold, the paper would demonstrate that the statistical G1-event technique developed for quasar samples can be applied to an individual, well-studied Seyfert galaxy, providing an independent cross-check against traditional excited-state diagnostics. The use of a forward-model inversion with explicit simulations is a transparent approach, and the authors include a systematic test of the DRW timescale. However, the load-bearing systematic uncertainties are not correctly propagated: the relation between t_r and R is misstated, and the inferred distance for component 1 shifts by a factor of ~1.9 when τ is varied within the plausible range, outside the quoted statistical error bars. The boxcar response assumption is also untested. These issues materially affect the claimed precision and the main conclusions.
major comments (3)
- [§5.2, Eqs. (1)–(2)] The statement 'Since Equations (1) and (2) imply ∆log R = 2 ∆log t_r' is incorrect. From Eq. (2), n_H ∝ R^{-2}, and from Eq. (1), t_r ∝ n_e^{-1} ∝ n_H^{-1}, so t_r ∝ R^2, i.e., ∆log R = 0.5 ∆log t_r. The inverted relation is used to argue that the τ_DRW systematic does not alter the conclusions. In fact, for component 1 the baseline t_r = 10^{-0.98} d and the τ=35 d value t_r = 10^{-0.43} d differ by 0.55 dex, which propagates to ∆log R ≈ 0.28 dex, shifting R from 0.77 to ~1.4 pc, well outside the quoted 1σ interval (0.67–0.87 pc). The error budget is therefore dominated by this systematic, and the paper should either fold it into the quoted uncertainties or report component 1 as an upper limit.
- [§3.2.3] The calibration mapping is built on the assumption that the absorber responds to ionizing continuum changes by averaging the continuum flux over a boxcar window of length t_r. This is an ad-hoc response function. The inferred t_r—and hence R—is a direct function of this choice; a different lag distribution (e.g., exponential, which is the more standard recombination response) would change the G1-vs-t_r calibration curve. The paper does not test alternative response functions. Given that the central claim rests on this mapping, the assumption needs either physical justification or a sensitivity test. Without it, the systematic uncertainty in t_r is larger than the quoted Monte Carlo errors.
- [§4, Table 2] For component 1, the inferred t_r = 10^{-0.98} d lies only 0.02 dex above the lower bound of the simulation grid (log t_r = -1.0). The quoted uncertainty of ±0.01 dex is thus a statistical uncertainty within a truncated calibration range. The paper itself notes that for τ_DRW = 3 d the value falls below the grid boundary. Therefore the component-1 measurement is effectively an upper limit, and reporting R = 0.77 ± 0.10 pc overstates the precision. The authors should either present this component as a limit (consistent with the previous work) or add a systematic error term that covers the τ dependence.
minor comments (3)
- [§3.2.2 and Table 4] The G2 classification criterion says 'Nσ > 3', but Table 4 lists rows where Nσ for the short interval is 0.00 and Nσ for the long interval is 6.74 (e.g., row 2), which is classified as G2. Clarify that the Nσ > 3 threshold applies only to the larger of the two Nσ values, or adjust the wording.
- [§4] The F(G1) uncertainties (±0.2%, ±0.5%) are binomial counting errors, but the events are not independent because many spectral triplets share spectra. This likely underestimates the true statistical error. Consider a bootstrap or block-resampling approach over epochs.
- [§4 vs. §1] The paper states that the inferred t_r for component 6 (10^{0.25} = 1.78 d) is 'broadly consistent' with the previous measurement of 4.83 ± 1.28 d. The difference is about 2.4σ by nominal errors; given the quoted 1σ range (3.55–6.11 d) does not include 1.78 d, 'broadly consistent' is an overstatement. Please quantify the comparison or soften the wording.
Circularity Check
No circular derivation: G1 probability is an external observable mapped to t_r by an explicit forward model, and t_r to R uses independent photoionization relations with external benchmarks.
full rationale
The claimed derivation chain is not circular. Observed F(G1) is measured from real HST/COS absorption-line variability (N_sigma classifications), independent of any assumed t_r. The F(G1)-t_r mapping is generated by injecting t_r into boxcar-smoothed DRW light curves and then inverting the observed fraction; this is parameter estimation from an explicit forward model, not an identity. The t_r to R step uses Arav et al. (2012) recombination-timescale equation and the definition of the ionization parameter with Cloudy photoionization parameters adopted from the authors' prior work; those parameters (U_H, n_CV/n_CIV, N_H) are derived from column-density ratios and the SED, not from R, so the distance is not an input to itself. The result is benchmarked against an external method for component 1 (Arav et al. 2015) and against literature values. Acknowledged limitations in the manuscript — component 1's t_r lying near the grid lower bound (Sec. 3.2.3) and the DRW-timescale systematic (Sec. 5.2) — are calibration-range and model-uncertainty issues, not circularity. Note that the paper's robustness statement 'Delta log R = 2 Delta log t_r' in Sec. 5.2 is inverted; Equations (1)-(2) imply Delta log R = 0.5 Delta log t_r, so the tau_DRW systematic is understated and could shift component 1's distance by about a factor of two. This is a correctness/robustness concern, but it does not make the derivation equivalent to its inputs by construction.
Assumptions & free parameters
free parameters (4)
- f (fractional ionizing continuum change) =
0.1
- τ_DRW at 200 Å =
10 days
- SF∞ at 200 Å =
1.579 mag
- Minimum sampled log10(t_r/days) =
-1.0
assumptions (5)
- domain assumption Absorption-line variability is driven by changes in the incident ionizing continuum
- ad hoc to paper The absorber responds by averaging the ionizing continuum over a boxcar window of length t_r
- ad hoc to paper The 200 Å ionizing continuum follows a damped random walk with τ=10 days and SF∞=1.579 mag
- domain assumption Equation (1) from Arav et al. (2012) describes the recombination timescale under ionization-recombination equilibrium
- domain assumption n_H ≈ 0.83 n_e for a fully ionized H/He plasma with negligible metals
Cite this review
Pith. "Pith review of Measuring Outflow Distances in NGC 5548 Using Absorption-Line Variability Diagnostics." pith.science (2026). https://pith.science/paper/XHGMN5SX
@misc{pith2026260721038,
author = {Pith},
title = {Pith review of: Measuring Outflow Distances in NGC 5548 Using Absorption-Line Variability Diagnostics},
year = {2026},
howpublished = {\url{https://pith.science/paper/XHGMN5SX}},
note = {Machine review of arXiv:2607.21038}
}
abstract
AGN-driven outflows serve as a key channel through which the energetic central engine influences host galaxy evolution. Among the physical properties of outflows, their radial distance from the galactic nucleus is particularly important for assessing AGN feedback. In this study, we investigate the UV outflow components in NGC 5548 by analyzing the variability of C IV absorption troughs in optical spectra obtained through multiple HST observations during 2013 and 2014. We construct a set of variability-based diagnostic events, labeled G1 and G2, which are sensitive to the recombination timescale ($t_r$) of ionized gas. By combining these with mock light curves generated using a damped random walk (DRW) model, we numerically establish a mapping between the G1 event probability and $t_r$. This approach allows us to constrain the radial distances of outflow components 1 and 6, whose absorption variability is primarily driven by changes in the incident ionizing continuum, to be $0.77^{+0.10}_{-0.10}$ and $1.72^{+1.74}_{-1.72}$ pc, respectively. These results are consistent with those obtained using a different method in our previous study, as well as with values reported in the literature.
Figures
Figures from the paper (2 more)
Reference graph
Works this paper leans on
-
[1]
Aller, M. C., & Richstone, D. O. 2007, ApJ, 665, 120, doi: 10.1086/519298
doi:10.1086/519298 2007
-
[2]
2018, ApJ, 857, 60, doi: 10.3847/1538-4357/aab494
Arav, N., Liu, G., Xu, X., et al. 2018, ApJ, 857, 60, doi: 10.3847/1538-4357/aab494
-
[3]
2012, A&A, 544, A33, doi: 10.1051/0004-6361/201118501
Arav, N., Edmonds, D., Borguet, B., et al. 2012, A&A, 544, A33, doi: 10.1051/0004-6361/201118501
-
[4]
Arav, N., Chamberlain, C., Kriss, G. A., et al. 2015, A&A, 577, A37, doi: 10.1051/0004-6361/201425302
-
[5]
Kriss, G. A. 2012, ApJ, 751, 107, doi: 10.1088/0004-637X/751/2/107
-
[6]
Bower, R. G., Benson, A. J., Malbon, R., et al. 2006, MNRAS, 370, 645, doi: 10.1111/j.1365-2966.2006.10519.x
arXiv 2006
-
[8]
Capellupo, D. M., Hamann, F., & Barlow, T. A. 2014, MNRAS, 444, 1893, doi: 10.1093/mnras/stu1502
-
[9]
2016, A&A, 592, A27, doi: 10.1051/0004-6361/201628464
Cappi, M., De Marco, B., Ponti, G., et al. 2016, A&A, 592, A27, doi: 10.1051/0004-6361/201628464
Show all 49 references
-
[10]
2023, RMxAA, 59, 327, doi: 10.22201/ia.01851101p.2023.59.02.12
Chatzikos, M., Bianchi, S., Camilloni, F., et al. 2023, RMxAA, 59, 327, doi: 10.22201/ia.01851101p.2023.59.02.12
2023 doi
-
[11]
2026, The Astrophysical Journal, 997, 245, doi: 10.3847/1538-4357/ae2858
Chen, Y., He, Z., & Liu, G. 2026, The Astrophysical Journal, 997, 245, doi: 10.3847/1538-4357/ae2858
2026 doi
-
[12]
C., et al
Chen, Z., He, Z., Ho, L. C., et al. 2022, Nature Astronomy, 6, 339, doi: 10.1038/s41550-021-01561-3 14
2022 doi
-
[13]
M., Kraemer, S
Crenshaw, D. M., Kraemer, S. B., Gabel, J. R., et al. 2003, ApJ, 594, 116, doi: 10.1086/376792 De Rosa, G., Peterson, B. M., Ely, J., et al. 2015, ApJ, 806, 128, doi: 10.1088/0004-637X/806/1/128 Del Zanna, G., Dere, K. P., Young, P. R., Landi, E., &
2003 doi
-
[14]
Mason, H. E. 2015, A&A, 582, A56, doi: 10.1051/0004-6361/201526827 Dovˇ ciak, M., Papadakis, I. E., Kammoun, E. S., & Zhang, W. 2022, A&A, 661, A135, doi: 10.1051/0004-6361/202142358
2015 doi
-
[15]
S., Kriss, G
Ebrero, J., Kaastra, J. S., Kriss, G. A., et al. 2016, A&A, 587, A129, doi: 10.1051/0004-6361/201527808 Filiz Ak, N., Brandt, W. N., Hall, P. B., et al. 2013, ApJ, 777, 168, doi: 10.1088/0004-637X/777/2/168
2016 doi
-
[16]
R., Korista, K
Goad, M. R., Korista, K. T., De Rosa, G., et al. 2016, ApJ, 824, 11, doi: 10.3847/0004-637X/824/1/11
2016 doi
-
[17]
R., Knigge, C., Korista, K
Goad, M. R., Knigge, C., Korista, K. T., et al. 2019, MNRAS, 486, 5362, doi: 10.1093/mnras/stz1186
2019 doi
-
[18]
2017, ApJ, 847, 132, doi: 10.3847/1538-4357/aa8d71
Guo, H., Wang, J., Cai, Z., & Sun, M. 2017, ApJ, 847, 132, doi: 10.3847/1538-4357/aa8d71
2017 doi
-
[19]
2019, MNRAS, 483, 1808, doi: 10.1093/mnras/sty2900
Hamann, F., Herbst, H., Paris, I., & Capellupo, D. 2019, MNRAS, 483, 1808, doi: 10.1093/mnras/sty2900
2019 doi
-
[20]
2025, A&A, 703, A305, doi: 10.1051/0004-6361/202556191
He, Z., & Wang, T. 2025, A&A, 703, A305, doi: 10.1051/0004-6361/202556191
2025 doi
-
[21]
2019, Nature Astronomy, 3, 265, doi: 10.1038/s41550-018-0669-8
He, Z., Wang, T., Liu, G., et al. 2019, Nature Astronomy, 3, 265, doi: 10.1038/s41550-018-0669-8
2019 doi
-
[22]
2022, Science Advances, 8, eabk3291, doi: 10.1126/sciadv.abk3291
He, Z., Liu, G., Wang, T., et al. 2022, Science Advances, 8, eabk3291, doi: 10.1126/sciadv.abk3291
2022 doi
-
[23]
2024, Science China
He, Z., Chen, Z., Liu, G., et al. 2024, Science China
2024
-
[24]
Physics, Mechanics, and Astronomy, 67, 129512, doi: 10.1007/s11433-024-2475-7
-
[25]
F., Hernquist, L., Cox, T
Hopkins, P. F., Hernquist, L., Cox, T. J., et al. 2006, ApJS, 163, 1, doi: 10.1086/499298
2006 doi
-
[26]
M., et al
Horne, K., De Rosa, G., Peterson, B. M., et al. 2021, ApJ, 907, 76, doi: 10.3847/1538-4357/abce60
2021 doi
- [27]
-
[28]
S., Kriss, G
Kaastra, J. S., Kriss, G. A., Cappi, M., et al. 2014, Science, 345, 64, doi: 10.1126/science.1253787
2014 doi
-
[29]
C., Bechtold, J., & Siemiginowska, A
Kelly, B. C., Bechtold, J., & Siemiginowska, A. 2009, ApJ, 698, 895, doi: 10.1088/0004-637X/698/1/895
2009 doi
-
[30]
2013, A&A, 551, L6, doi: 10.1051/0004-6361/201220923 Koz lowski, S., Kochanek, C
Kollatschny, W., & Zetzl, M. 2013, A&A, 551, L6, doi: 10.1051/0004-6361/201220923 Koz lowski, S., Kochanek, C. S., Udalski, A., et al. 2010, ApJ, 708, 927, doi: 10.1088/0004-637X/708/2/927
2013 doi
-
[31]
A., De Rosa, G., Ely, J., et al
Kriss, G. A., De Rosa, G., Ely, J., et al. 2019, ApJ, 881, 153, doi: 10.3847/1538-4357/ab3049
2019 doi
-
[32]
Ferland, G. J. 2022, MNRAS, 516, 4397, doi: 10.1093/mnras/stac2443
2022 doi
-
[33]
M., Dietrich, M., & Barber, S
Leighly, K. M., Dietrich, M., & Barber, S. 2011, ApJ, 728, 94, doi: 10.1088/0004-637X/728/2/94
2011 doi
-
[34]
T., & Dietrich, M
Richards, G. T., & Dietrich, M. 2018, ApJ, 866, 7, doi: 10.3847/1538-4357/aadee6
2018 doi
-
[35]
B., Leighly, K
Lucy, A. B., Leighly, K. M., Terndrup, D. M., Dietrich, M., & Gallagher, S. C. 2014, ApJ, 783, 58, doi: 10.1088/0004-637X/783/1/58
2014 doi
-
[36]
1969, Nature, 223, 690, doi: 10.1038/223690a0
Lynden-Bell, D. 1969, Nature, 223, 690, doi: 10.1038/223690a0
1969 doi
-
[37]
L., Ivezi´ c,ˇZ., Kochanek, C
MacLeod, C. L., Ivezi´ c,ˇZ., Kochanek, C. S., et al. 2010, ApJ, 721, 1014, doi: 10.1088/0004-637X/721/2/1014
2010 doi
-
[38]
L., Ivezi´ c,ˇZ., Sesar, B., et al
MacLeod, C. L., Ivezi´ c,ˇZ., Sesar, B., et al. 2012, ApJ, 753, 106, doi: 10.1088/0004-637X/753/2/106
2012 doi
-
[39]
S., Mehdipour, M., et al
Mao, J., Kaastra, J. S., Mehdipour, M., et al. 2017, A&A, 607, A100, doi: 10.1051/0004-6361/201731378
2017 doi
-
[40]
G., Schaye, J., Ponman, T
McCarthy, I. G., Schaye, J., Ponman, T. J., et al. 2010, MNRAS, 406, 822, doi: 10.1111/j.1365-2966.2010.16750.x
2010
-
[41]
A., Kaastra, J
Mehdipour, M., Kriss, G. A., Kaastra, J. S., et al. 2024, ApJ, 962, 155, doi: 10.3847/1538-4357/ad1bcb
2024 doi
-
[42]
S., Kriss, G
Mehdipour, M., Kaastra, J. S., Kriss, G. A., et al. 2015, A&A, 575, A22, doi: 10.1051/0004-6361/201425373
2015 doi
-
[43]
M., Barth, A
Pei, L., Fausnaugh, M. M., Barth, A. J., et al. 2017, ApJ, 837, 131, doi: 10.3847/1538-4357/aa5eb1
2017 doi
-
[44]
Rees, M. J. 1984, ARA&A, 22, 471, doi: 10.1146/annurev.aa.22.090184.002351 Sch¨ onell, Jr., A. J., Storchi-Bergmann, T., Riffel, R. A., &
1984
-
[45]
2017, MNRAS, 464, 1771, doi: 10.1093/mnras/stw2263
Riffel, R. 2017, MNRAS, 464, 1771, doi: 10.1093/mnras/stw2263
2017 doi
-
[46]
S., Hopkins, P
Somerville, R. S., Hopkins, P. F., Cox, T. J., Robertson, B. E., & Hernquist, L. 2008, MNRAS, 391, 481, doi: 10.1111/j.1365-2966.2008.13805.x
2008
-
[47]
2005, MNRAS, 361, 776, doi: 10.1111/j.1365-2966.2005.09238.x
Springel, V., Di Matteo, T., & Hernquist, L. 2005, MNRAS, 361, 776, doi: 10.1111/j.1365-2966.2005.09238.x
2005
-
[48]
R., Pancoast, A., Treu, T., et al
Williams, P. R., Pancoast, A., Treu, T., et al. 2020, ApJ, 902, 74, doi: 10.3847/1538-4357/abbad7
2020 doi
-
[49]
S., Peterson, B
Yu, Z., Kochanek, C. S., Peterson, B. M., et al. 2020, MNRAS, 491, 6045, doi: 10.1093/mnras/stz3464
2020 doi
-
[50]
2021, ApJL, 906, L8, doi: 10.3847/2041-8213/abd318
Zhao, Q., He, Z., Liu, G., et al. 2021, ApJL, 906, L8, doi: 10.3847/2041-8213/abd318
2021 doi
Reviewed August 1, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.