REVIEW 3 major objections 5 minor 33 references
A failed wind candidate in NGC 3783 from the 2001 year campaign with Chandra/HETGS
T0 review · 3 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read The paper identifies a low-ionization warm-absorber component in NGC 3783 as a failed wind candidate, bound within 0.27 pc of the black hole.
desk verdict A careful reanalysis of the 2001 Chandra campaign that solidifies the Fe UTA variability and adds a tpho-based density/distance estimate for one warm absorber component, but the failed-wind conclusion is not yet supported because the outflow velocity is never stated. 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 load-bearing machinery is the time-dependent photoionization model tpho, which follows the ionization, heating, and cooling of each warm-absorber component in response to an input ionizing light curve built from daily-binned RXTE/PCA monitoring plus the Chandra epochs. The key observable is the ratio spectrum of low-flux over high-flux Chandra/MEG data in the 15–18 Å band, where Fe UTA and O vii edge absorption dominate; the amplitude of the ratio variation as a function of density, set by the recombination timescale, converts the observed variability into a density measurement. The escape criterion $n_{\rm H,upp} = L_{\rm ion} v_{\rm out}^4/[4(GM_{\rm BH})^2 \xi]$ then turns that density limit into the failed-wind conclusion.
What would settle it
A future high-cadence X-ray campaign on NGC 3783 that resolves the Fe UTA and O vii edge during a flux change and measures a recombination lag implying a density below $10^{12.3}\ \mathrm{m}^{-3}$ would falsify the lower limit; so would demonstrating that the observed ratio spectra can be reproduced with a density near $10^{10}\ \mathrm{m}^{-3}$ plus a varying covering fraction.
Extended reading notes
Core claim
The central claim is that component 5 of the NGC 3783 warm absorber, with $\log \xi = 1.65$ and column density $N_{\rm H} = 0.5\times 10^{26}\ \mathrm{m}^{-2}$, has a hydrogen density higher than $10^{12.3}\ \mathrm{m}^{-3}$ at $3\sigma$ confidence, so it lies within 0.27 pc of the central black hole. The paper reaches this by comparing the low-flux to high-flux Chandra/MEG ratio spectra and showing that the Fe UTA and O vii edge variations require the gas to recombine on the observed timescales; lower-density gas has recombination timescales too long to produce the changes. Comparing this density limit to the maximum density for an escaping wind, $n_{\rm H,upp} = L_{\rm ion} v_{\rm out}^4/[4(GM_{\rm BH})^2 \xi]$, shows that the component's velocity is below the escape velocity. The paper concludes that this low-ionization component is a failed wind candidate, gas launched outward but unable to escape the supermassive black hole, and therefore remains bound to the system.
Load-bearing premise
The density and distance limits hold only if the 2001 spectral changes are caused entirely by the gas's ionization state responding to the monitored continuum, and if neither the amount of gas blocking the light, the gas column, nor the shape of the ionizing spectrum changed during the campaign.
Editorial extensions
If this is right
- The low-ionization warm absorber in NGC 3783 is located within 0.27 pc of the black hole, much closer than the earlier 6 pc equilibrium-based estimate.
- Because this gas is denser than the maximum for an escaping wind, it should be counted as bound material that may fall back rather than as part of the AGN's mass and energy output.
- Fe UTA variability is a practical density diagnostic for warm absorbers, since iron's short recombination timescale lets it respond to continuum changes that oxygen cannot follow.
- Photoionization equilibrium is insufficient: the time-dependent model improves the fit at the $3.1\sigma$ level, so equilibrium-based density estimates for variable warm absorbers are biased.
- Oxygen recombination, with its longer timescale, is the more sensitive probe for higher-density plasma in future high-resolution X-ray observations.
Reading between the lines
- The paper leaves implicit that a population of failed winds would lower the net outflow rates of Seyfert galaxies, since only gas above the escape velocity actually leaves the system.
- A stalled wind near 0.27 pc is a natural physical ingredient for the high-ionization broad-line region and for the obscuring material, a connection the paper mentions but does not develop.
- A testable extension is to apply the same ratio-spectrum and tpho analysis to other bright Seyferts with archival multi-epoch gratings; similar Fe UTA variations would turn this single-object failed-wind candidate into a general phenomenon.
- If the covering fraction of component 5 changed during the 2001 campaign, the current density limit could be an artifact of that change rather than of recombination timing; future observations that monitor the covering fraction independently would settle this.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reanalyzes five Chandra/HETGS observations of NGC 3783 from 2001, constructs ratio spectra between low- and high-flux states, detects >10-sigma variability in the Fe UTA and O vii edge region, and uses the tpho time-dependent photoionization model to attribute this variability to component 5 of the warm absorber (log xi = 1.65). A chi-square scan over density yields a 3-sigma lower limit n_H > 10^12.3 m^-3, corresponding to a distance r < 0.27 pc, and the authors argue that this density exceeds the maximum density for an escaping wind, making component 5 a failed wind candidate.
Significance. If the density constraint holds, the paper provides one of the first direct variability-based density measurements placing a low-ionization warm absorber on sub-parsec scales, with consequences for wind-launching and AGN feedback models. The analysis is careful in several respects: the ratio-spectrum method is transparent, the >10-sigma Fe UTA variability is credible, the tpho-versus-pion comparison is a useful diagnostic, and the paper explicitly tests the soft-excess sensitivity in Fig. 10. The density lower limit is an empirical product of a chi-square fit rather than a quantity fixed by construction. However, the failed-wind conclusion is not self-contained because it depends on an unstated outflow velocity, and the density constraint depends on adopting the group's own nine-component Mao et al. (2019) model for the 2001 epoch. These issues do not undermine the density measurement itself, but they currently limit the strength of the central astrophysical claim.
major comments (3)
- [Section 5, Eq. (2)] The failed-wind criterion is not evaluated with an explicit outflow velocity. Eq. (2) defines n_H,upp in terms of v_out, but the text does not quote a v_out for component 5 and instead refers to Figure 8 of Li et al. (2023). Using the paper's own values (L_ion = 6.36 x 10^36 W, M_BH = 2.82 x 10^7 M_sun, log xi = 1.65), Eq. (2) gives n_H,upp ~ 2.5 x 10^12 m^-3 for v_out = 1000 km/s, which already exceeds the reported 3-sigma lower limit of 2.0 x 10^12 m^-3. Thus the claim n_H > n_H,upp holds only for v_out below roughly 950 km/s. The authors must state the adopted v_out and justify it, or the 'failed wind' conclusion is unsupported as presented.
- [Section 3, paragraph 1; Section 4, Fig. 9] The density constraint assumes that the nine-component equilibrium warm-absorber model of Mao et al. (2019), specifically component 5 with log xi = 1.65 and N_H = 0.5 x 10^26 m^-2, is valid for the 2001 epoch and that the observed ratio variability is entirely an ionization response to the RXTE-monitored continuum. A change in covering fraction, column density, or SED shape, or a different absorber decomposition, would alter the derived density and distance limits. The paper should quantify the sensitivity of the 3-sigma lower limit to these assumptions, or at minimum show that the ratio-spectrum shape discriminates against column-density and covering-fraction variability.
- [Section 4, paragraph beginning 'Because the sampling...' and Fig. 9] ObsID 2094 is excluded from the combined chi-square fit because of sparse RXTE sampling, yet the text also states that for this observation the tpho model ratio 'aligns well with the pion model.' The paper should state explicitly that the 10^12.3 m^-3 lower limit is derived from three low/high-flux comparisons only (ObsIDs 2090, 2091, 2092 versus 2093), and clarify whether including 2094 at a lower statistical weight would change the confidence level. This is a transparency issue that directly affects how readers interpret Fig. 9.
minor comments (5)
- [Section 4, paragraph before Fig. 9] The sentence 'we get reliable predictions for densities up to 10 m^-3' is missing an exponent; from the context it should presumably read 'up to 10^13 m^-3', and the following sentence contains the same typo.
- [Section 4, Fig. 9] The reported minimum chi-square value of 128.76 is not accompanied by the number of degrees of freedom or the delta-chi-square value used to define the 3-sigma threshold; please report these so that the statistical statement is reproducible.
- [Section 4, paragraph on the O vii RRC] The assumption that the O vii RRC does not vary intrinsically is reasonable but should be flagged as a model assumption in the conclusions; the current text asserts it without a quantitative test of how a varying RRC would affect the density limit.
- [Introduction, paragraph 3] The phrase 'As a another approach' should be corrected to 'As another approach'.
- [Eq. (2)] The symbols in Eq. (2) are not all defined in the text; in particular, the units of n_H,upp and the meaning of G (Newton's constant versus gravitational parameter) should be stated explicitly for clarity.
Circularity Check
No definitional circularity: the density lower limit is an empirical chi-square fit, but the failed-wind conclusion is imported from a same-group figure with an unstated outflow velocity.
-
self citation load bearing
[Section 5, Discussion, paragraph following Eq. (2)]
"The lower limit density nH of component 5 from our present tpho modelling is larger than nH,upp (see Figure 8 of Li et al. 2023), which indicates that the velocity of component 5 remains below the escape velocity. Therefore, this velocity discrepancy suggests that component 5 may be a failed wind, unable to reach velocities sufficient to escape the gravitational pull of the SMBH, thus remaining bound to the system."
The entire evidence for the failed-wind interpretation is the inequality n_H > n_H,upp, but n_H,upp in Eq. (2) depends on v_out^4, and this paper nowhere states v_out for component 5 or reproduces the comparison. The assertion is instead referred to Figure 8 of Li et al. (2023), a paper sharing the first author and several co-authors. Thus the central interpretive claim is imported from a same-group figure rather than derived self-containedly from the present measurement; without the cited figure's unstated v_out, the inequality cannot be checked. This is load-bearing self-citation, though not a definitional circularity because the lower-limit density itself comes from an empirical chi-square fit to observed ratio spectra.
full rationale
The central density constraint (n_H > 10^12.3 m^-3 at 3 sigma) is obtained by a chi-square scan of tpho model ratios against observed Chandra/MEG ratio spectra (Figs. 8 and 9), so it is an empirical fit, not a quantity defined by construction. Eq. (1) is the standard ionization parameter; r < 0.27 pc follows algebraically from the fitted density, assumed xi, and L_ion. I find no step in which a 'prediction' is equivalent to an input or in which a fitted parameter is renamed as a prediction. The nine-component decomposition (Mao et al. 2019), the SED (Mehdipour et al. 2017), and the tpho code (Rogantini et al. 2022) are prior work by the same group; using them as initial conditions is normal practice and not circular in itself, provided those prior fits are data-based. The one load-bearing self-reference is the failed-wind comparison: the paper states n_H > n_H,upp 'see Figure 8 of Li et al. (2023)', but Eq. (2) depends on v_out^4 and no v_out for component 5 is given or derived here. That makes the central interpretation hinge on a same-group figure rather than a self-contained calculation; it is a missing-support and correctness risk and contributes to a moderate score, but it is not a definitional circularity because the density measurement itself is independent.
Assumptions & free parameters
free parameters (5)
- Density of warm absorber component 5 (nH) =
10^13.4 m^-3 (best fit); lower limit 10^12.3 m^-3 at 3 sigma
- Local continuum scaling factor per epoch =
2.21, 2.30, 1.87, 1.66
- Component 5 ionization parameter log xi =
1.65
- Component 5 column density N_H =
0.5 x 10^26 m^-2
- Ionizing luminosity Lion =
6.36 x 10^36 W
assumptions (6)
- domain assumption The RXTE 2-60 keV light curve tracks the 1-1000 Ryd ionizing continuum with constant SED shape.
- domain assumption The nine-component warm absorber model of Mao et al. (2019), including component 5 with log xi = 1.65 and N_H = 0.5e26 m^-2, describes the 2001 epoch.
- domain assumption Spectral changes are pure ionization and recombination response at fixed column density and covering factor, with no change in covering fraction, column density, or velocity.
- ad hoc to paper The O vii radiative recombination continuum does not vary intrinsically between epochs.
- domain assumption The tpho code and SPEX atomic data correctly compute the time-dependent ionization and transmission.
- domain assumption The failed-wind criterion and the comparison density nH,upp from Li et al. (2023) apply to component 5 without restating v_out.
Cite this review
Pith. "Pith review of A failed wind candidate in NGC 3783 from the 2001 year campaign with Chandra/HETGS." pith.science (2026). https://pith.science/paper/RFMLWBRT
@misc{pith2026250116880,
author = {Pith},
title = {Pith review of: A failed wind candidate in NGC 3783 from the 2001 year campaign with Chandra/HETGS},
year = {2026},
howpublished = {\url{https://pith.science/paper/RFMLWBRT}},
note = {Machine review of arXiv:2501.16880}
}
abstract
We reanalyze the Chandra/HETGS observations of NGC 3783 from the campaign in the year 2001, identifying significant spectral variations in the Fe unresolved transition array (UTA) over timescales of weeks to months. These changes correlate with a $1.4-2$ fold increase in the ionizing continuum and exceed $10 \, \sigma$ significance. The variations primarily originate from a low-ionization state ($\rm log \xi = 1.65$) component of the warm absorber. Time-dependent photoionization modelling confirms the sensitivity of this low-ionization component to continuum variations within the Fe UTA band. Local fitting indicates a lower density limit of $>10^{12.3} \, \rm m^{-3}$ at $3 \, \sigma$ statistical uncertainty, with the component located within $0.27 \, \rm pc$. Our findings suggest that this low-ionization component is a potential failed wind candidate.
Figures
Figures from the paper (6 more)
Reference graph
Works this paper leans on
- [1]
-
[2]
Ballhausen, R., Kallman, T. R., Gu, L., & Paerels, F. 2023, ApJ, 956, 65
work page 2023
-
[3]
2018, in Society of Photo- Optical Instrumentation Engineers (SPIE) Conference Series, V ol
Barret, D., Lam Trong, T., den Herder, J.-W., et al. 2018, in Society of Photo- Optical Instrumentation Engineers (SPIE) Conference Series, V ol. 10699, Space Telescopes and Instrumentation 2018: Ultraviolet to Gamma Ray, ed. J.-W. A. den Herder, S. Nikzad, & K. Nakazawa, 106991G
work page 2018
-
[4]
P., Blustin, A
Behar, E., Rasmussen, A. P., Blustin, A. J., et al. 2003, ApJ, 598, 232
2003
-
[5]
Behar, E., Sako, M., & Kahn, S. M. 2001, ApJ, 563, 497
2001
-
[6]
C., Williams, P
Bentz, M. C., Williams, P. R., Street, R., et al. 2021, ApJ, 920, 112
2021
-
[7]
J., Branduardi-Raymont, G., Behar, E., et al
Blustin, A. J., Branduardi-Raymont, G., Behar, E., et al. 2002, A&A, 392, 453
2002
- [8]
Show all 33 references
-
[9]
2020, MNRAS, 498, 4150
Davies, R., Baron, D., Shimizu, T., et al. 2020, MNRAS, 498, 4150
2020
-
[10]
I., Burtscher, L., Rosario, D., et al
Davies, R. I., Burtscher, L., Rosario, D., et al. 2015, ApJ, 806, 127
2015
-
[11]
R., Kraemer, S
Gabel, J. R., Kraemer, S. B., Crenshaw, D. M., et al. 2005, ApJ, 631, 741
2005
-
[12]
C., Miller, J
Gallo, L. C., Miller, J. M., & Costantini, E. 2023, arXiv e-prints, arXiv:2302.10930
2023 arXiv
-
[13]
& Proga, D
Giustini, M. & Proga, D. 2021, in IAU Symposium, V ol. 356, Nuclear Activity in Galaxies Across Cosmic Time, ed. M. Povi ´c, P. Marziani, J. Masegosa, H. Netzer, S. H. Negu, & S. B. Tessema, 82–86 GRA VITY Collaboration, Amorim, A., Bauböck, M., et al. 2021, A&A, 648, A117
2021
-
[14]
2023, A&A, 679, A43 Juráˇnová, A., Costantini, E., & Uttley, P
Gu, L., Kaastra, J., Rogantini, D., et al. 2023, A&A, 679, A43 Juráˇnová, A., Costantini, E., & Uttley, P. 2022, MNRAS, 510, 4225
2023
-
[15]
S., Detmers, R
Kaastra, J. S., Detmers, R. G., Mehdipour, M., et al. 2012, A&A, 539, A117
2012
-
[16]
S., Mehdipour, M., Behar, E., et al
Kaastra, J. S., Mehdipour, M., Behar, E., et al. 2018, A&A, 619, A112
2018
-
[17]
S., Mewe, R., & Nieuwenhuijzen, H
Kaastra, J. S., Mewe, R., & Nieuwenhuijzen, H. 1996, in UV and X-ray Spec- troscopy of Astrophysical and Laboratory Plasmas, ed. K. Yamashita & T. Watanabe, 411–414
1996
-
[18]
S., Raassen, A
Kaastra, J. S., Raassen, A. J. J., de Plaa, J., & Gu, L. 2024, SPEX X-ray spectral fitting package, (3.08.01). Zenodo
2024
-
[19]
N., George, I
Kaspi, S., Brandt, W. N., George, I. M., et al. 2002, ApJ, 574, 643
2002
-
[20]
H., McKee, C
Krolik, J. H., McKee, C. F., & Tarter, C. B. 1981, ApJ, 249, 422
1981
-
[21]
S., et al
Krongold, Y ., Nicastro, F., Brickhouse, N. S., et al. 2003, ApJ, 597, 832
2003
-
[22]
S., Elvis, M., & Mathur, S
Krongold, Y ., Nicastro, F., Brickhouse, N. S., Elvis, M., & Mathur, S. 2005, ApJ, 622, 842
2005
-
[23]
S., Gu, L., & Mehdipour, M
Li, C., Kaastra, J. S., Gu, L., & Mehdipour, M. 2023, A&A, 680, A44
2023
-
[24]
S., Mehdipour, M., et al
Mao, J., Kaastra, J. S., Mehdipour, M., et al. 2017, A&A, 607, A100
2017
-
[25]
S., et al
Mao, J., Mehdipour, M., Kaastra, J. S., et al. 2019, A&A, 621, A99
2019
-
[26]
S., Kriss, G
Mehdipour, M., Kaastra, J. S., Kriss, G. A., et al. 2017, A&A, 607, A28
2017
-
[27]
2003, ApJ, 599, 933
Netzer, H., Kaspi, S., Behar, E., et al. 2003, ApJ, 599, 933
2003
-
[28]
2022, ApJ, 940, 122
Rogantini, D., Mehdipour, M., Kaastra, J., et al. 2022, ApJ, 940, 122
2022
-
[29]
R., Bautista, M
Sadaula, D. R., Bautista, M. A., García, J. A., & Kallman, T. R. 2023, ApJ, 946, 93
2023
-
[30]
V ., Uttley, P., & Costantini, E
Silva, C. V ., Uttley, P., & Costantini, E. 2016, A&A, 596, A79
2016
-
[31]
B., Tucker, W
Tarter, C. B., Tucker, W. H., & Salpeter, E. E. 1969, ApJ, 156, 943
1969
-
[32]
2020, in Society of Photo-Optical In- strumentation Engineers (SPIE) Conference Series, V ol
Tashiro, M., Maejima, H., Toda, K., et al. 2020, in Society of Photo-Optical In- strumentation Engineers (SPIE) Conference Series, V ol. 11444, Space Tele- scopes and Instrumentation 2020: Ultraviolet to Gamma Ray, ed. J.-W. A. den
2020
-
[33]
1998, A&AS, 130, 333 Article number, page 10 of 10
Theureau, G., Bottinelli, L., Coudreau-Durand, N., et al. 1998, A&AS, 130, 333 Article number, page 10 of 10
1998
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.