arxiv: 2604.21710 · v1 · submitted 2026-04-23 · 🌌 astro-ph.HE

Recognition: unknown

Delving into the depths of NGC 3783 with XRISM: V. Broad-band modeling of ionized outflows

Keqin Zhao , Jelle S. Kaastra , Liyi Gu , Missagh Mehdipour , Megan E. Eckart , Keigo Fukumura , Matteo Guainazzi , Chen Li

show 2 more authors

Christos Panagiotou Matilde Signorini

Authors on Pith no claims yet

Pith reviewed 2026-05-09 20:47 UTC · model grok-4.3

classification 🌌 astro-ph.HE

keywords warm absorbersNGC 3783X-ray outflowsSeyfert galaxiesphotoionization modelingionized windslong-term variabilityXRISM

0 comments

The pith

The warm absorbers in NGC 3783 show significant structural and dynamical evolution over the past 24 years.

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

The paper compares new 2024 X-ray spectra of the Seyfert galaxy NGC 3783 from XRISM and XMM-Newton against earlier observations from 2000-2001. It models the multiphase ionized outflows using photoionization calculations to extract ionization, column density, velocity, and turbulence for each component. The total absorbing column has risen by a factor of 1.5, one dominant component has tripled its column while keeping similar ionization, and the set of components now contains fewer low-ionization high-velocity features. These shifts indicate that the outflow is not static but has been replenished and restructured on decade timescales. The result bears on how material is launched and sustained near the central black hole.

Core claim

Eight warm absorber components are identified in the 2024 data, spanning log ξ from 1.08 to 3.38 and outflow velocities from 480 to 1230 km s^{-1}. Their column densities and turbulent velocities are broadly consistent with the 2000-2001 epoch, yet the earlier data showed more low-ionization high-velocity components. The total column density is now 1.5 times higher, the dominant unresolved transition array component has increased its column by a factor of three at similar ionization, and the absorber population therefore requires fresh material. The WAs in NGC 3783 have undergone significant structural and dynamical evolution over the past 24 years.

What carries the argument

Joint photoionization modeling of the XMM-Newton/RGS and XRISM/Resolve spectra with the pion code to parameterize eight absorption components and compare their properties to the 2000-2001 epoch.

If this is right

The outflow requires ongoing replenishment by fresh gas from the accretion flow or disk.
The dominant component can increase its column density by a factor of three while preserving ionization, implying changes in density, location, or covering fraction.
Fewer low-ionization high-velocity components today suggests that some earlier structures have disappeared or weakened.
The overall absorber population is dynamic on 24-year timescales rather than fixed.

Where Pith is reading between the lines

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

If decade-scale changes prove common in other Seyfert galaxies, static single-epoch models of warm absorbers will need to be replaced by time-dependent wind calculations.
Repeated high-resolution observations spaced by years could map how individual components appear, strengthen, or fade, directly constraining wind-launching radii and mass-loss rates.
The observed column increase may correlate with measurable changes in the central accretion rate or disk state that future multi-wavelength monitoring could test.

Load-bearing premise

The eight absorption components identified in the 2024 spectra correspond to physically comparable structures in the 2000-2001 data, so that differences reflect real changes rather than variations in data quality or modeling choices.

What would settle it

Re-fitting the 2000-2001 spectra with the identical pion modeling setup, component selection criteria, and data-quality thresholds used for the 2024 data, then obtaining the same eight components with matching parameters, would show that no evolution occurred.

Figures

Figures reproduced from arXiv: 2604.21710 by Chen Li, Christos Panagiotou, Jelle S. Kaastra, Keigo Fukumura, Keqin Zhao, Liyi Gu, Matilde Signorini, Matteo Guainazzi, Megan E. Eckart, Missagh Mehdipour.

**Figure 1.** Figure 1: XMM-Newton/RGS spectrum of NGC 3783 with our best-fit model. The top panel shows a blow-up of the spectrum near the UTA feature. The strongest absorption features are labeled. Our best-fit model ( [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗

**Figure 2.** Figure 2: XRISM/Resolve spectrum of NGC 3783 with our best-fit model. The top panel show the high ionization WA features. For clarity of display the spectrum in this figure is additionally binned up. The strongest emission and absorption features are labeled. Our best-fit model ( [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

**Figure 3.** Figure 3: Relations between the parameters of the eight outflow components derived from the XRISM Resolve spectrum of NGC 3783. The top panel displays the column density NH as functions of the ionization parameter (log ξ). The middle and bottom panels show the turbulent velocity (σv) and outflow velocity (vout) as functions of log ξ, respectively. The component label ( [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗

**Figure 4.** Figure 4: Curve of thermal stability. The SED is from Mehdipour et al. (2017) in the 2001 unobscured state. The orange triangle and purple square corresponds to the 2000–2001 and 2024 WAs, respectively. The name of the WAs in 2000–2001 and 2004 are marked by the same colour as the data points. 4.3. The variation of the UTA component The Fe M-shell 2p–3d unresolved transition array (UTA) between 16 and 17 Å, togethe… view at source ↗

**Figure 5.** Figure 5: Upper panel: The RGS spectra of 2001 and 2024 with their bestfit model in the UTA energy band. Lower panel: Transmission model of the UTA in 2001 (Comp. 5) and 2024 (Comp. B3). 5. Conclusions We have analysed joint time-averaged XMM-Newton/RGS and XRISM/Resolve spectra of NGC 3783 from the 2024 observations, revealing the rich complexities of eight warm absorber (WA) components. The main results are summ… view at source ↗

**Figure 6.** Figure 6: Column densities of Fe and O in the UTA component B3. The upper panel shows the column density of Fe viii to Fe xvii calculated by pion (green) and slab (orange) model. The lower panel shows the column density of O vi to O viii calculated by pion (green) and slab (orange) model. WA observed in 2000–2001 undergoing failed-wind dynamics and transverse structural evolution. Acknowledgements. K. Zhao acknowled… view at source ↗

read the original abstract

The Seyfert 1 galaxy NGC 3783 hosts a multiphase warm absorber (WA) that has been extensively studied in the X-ray band. High-resolution spectra from 2000-2001 revealed a complex outflow with multiple ionization and velocity components. Two decades later, new XMM-Newton and XRISM observations allow us to investigate the long-term evolution of these outflows. We perform joint spectral modeling of the XMM-Newton/RGS and XRISM/Resolve time-averaged spectra using the pion photoionization code within SPEX. We derive the ionization parameter, column density, turbulent velocity, and outflow velocity for each absorption component, and investigate their thermal stability and Absorption Measure Distribution (AMD) to characterize the physical and dynamical properties of the WA in NGC 3783 in 2024. We compare these results with the 2000-2001 epoch to assess long-term variability, stability, and possible changes in the absorber population. We identify eight WA components spanning log $\xi =$ 1.08-3.38 and outflow velocities of 480-1230 km s$^{-1}$. The ranges of column densities and turbulent velocities remain broadly consistent with the WAs from 2000-2001, but the earlier data contained more low-ionization, high-velocity components. The total column density in 2024 is 1.5 times larger than in 2000-2001, requiring replenishment by fresh material. The dominant Unresolved Transition Array (UTA) absorber (Comp. B3) has increased its column density by a factor of three while maintaining a similar ionization parameter. The WAs in NGC 3783 have undergone significant structural and dynamical evolution over the past 24 years.

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.

Desk Editor's Note private letter to a colleague

New XRISM spectra add a 2024 epoch for NGC 3783's warm absorber with a reported 1.5x column increase, but the evolution claim depends on untested assumptions about component correspondence across epochs.

read the letter

The main thing to know is that this paper models fresh XRISM and XMM-Newton spectra of NGC 3783 from 2024 using pion in SPEX, fits eight absorption components, and directly compares the results to the 2000-2001 data. They find broadly overlapping ranges for column densities and velocities but note a higher total column now, loss of some low-ionization high-velocity parts, and one dominant component tripling in column while keeping similar ionization; they interpret this as requiring fresh material and overall structural change over 24 years. They also look at the absorption measure distribution and thermal stability for context. This is a straightforward observational extension that puts new high-resolution data on a well-studied source and attempts the epoch comparison in one place. The methods are standard and the parameter ranges being consistent is a useful check. The soft spot is the leap to significant evolution. The abstract does not describe explicit criteria for matching the eight 2024 components to the earlier ones, nor does it show tests like whether the old spectra can be adequately described by the new component set or vice versa. Without that, or reported fit statistics and uncertainties, it is difficult to separate real changes from differences in resolution, signal-to-noise, or modeling choices. The stress-test point about assuming physical identity holds up from the summary. This is mainly for researchers who track detailed variability in individual AGN outflows and need the latest numbers on this particular object. It adds a concrete recent data point rather than a broad new method. I would send it to peer review so referees can check the actual fits, figures, and any cross-epoch validation in the full text.

Referee Report

3 major / 3 minor

Summary. The manuscript reports joint spectral modeling of 2024 XMM-Newton/RGS and XRISM/Resolve time-averaged spectra of NGC 3783 using the pion photoionization code in SPEX. It derives ionization parameter, column density, turbulent velocity, and outflow velocity for eight warm absorber components spanning log ξ = 1.08–3.38 and velocities 480–1230 km s⁻¹, analyzes their thermal stability and AMD, and compares these to 2000–2001 epoch results to conclude that the outflows have undergone significant structural and dynamical evolution, including a 1.5× increase in total column density, loss of some low-ionization high-velocity components, and a factor-of-three column increase in the dominant UTA component (B3) while maintaining similar ionization.

Significance. If the cross-epoch component correspondence holds, the work would offer rare observational constraints on AGN warm-absorber variability over ~24-year baselines, including evidence for material replenishment. The high-resolution XRISM data and use of a public photoionization code for joint fitting are strengths that enable better separation of velocity and ionization components than prior lower-resolution observations. The AMD and stability analysis adds physical context. The significance is reduced by the lack of explicit validation for the component-matching assumption that underpins the evolution claim.

major comments (3)

[Abstract and comparison section] Abstract and comparison section: The headline claim that the WAs have undergone significant structural and dynamical evolution (including 1.5× total column increase and population changes) rests on the untested assumption that the eight 2024 components map to the same physical structures as those identified in 2000–2001. No explicit matching criteria (e.g., velocity or ionization proximity thresholds) or cross-validation tests (such as refitting the archival spectra with the 2024 component set or vice versa and reporting ΔC-stat) are described; without these, differences could arise from improved resolution/S/N or modeling choices rather than time evolution.
[Results section on fitted parameters] Results section on fitted parameters: The abstract and text report parameter ranges but provide no fit statistics (C-stat or χ²/dof), 1σ error bars on log ξ, N_H, v_out, or v_turb for the eight components, or discussion of degeneracies among them. This undermines quantitative claims such as the factor-of-three column increase in Comp. B3 and the overall 1.5× total N_H rise, as it is impossible to assess whether the reported changes exceed uncertainties or are driven by component splitting.
[Section on total column density and replenishment] Section on total column density and replenishment: The inference that fresh material is required because the 2024 total column is 1.5 times larger assumes direct additivity and comparability of the summed N_H values across epochs. The manuscript does not address whether the increase could result from the choice to fit eight components (versus fewer in prior work) or from unmodeled continuum or emission features that XRISM resolves but earlier data did not.

minor comments (3)

[Abstract] The abstract states that ranges are 'broadly consistent' with 2000–2001 but does not include a side-by-side table of parameters with uncertainties or quantify the overlap in velocity/ionization space.
[Results] Component labeling (e.g., 'Comp. B3' for the UTA absorber) should explicitly reference the nomenclature used in the 2000–2001 papers to facilitate direct comparison.
[Methods] The manuscript should state the number of free parameters per component and any tied parameters in the joint fit to allow reproducibility assessment.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their constructive and detailed report. We address each major comment below, providing the strongest honest defense of the manuscript while agreeing to revisions where the points identify genuine gaps in presentation or justification.

read point-by-point responses

Referee: [Abstract and comparison section] Abstract and comparison section: The headline claim that the WAs have undergone significant structural and dynamical evolution (including 1.5× total column increase and population changes) rests on the untested assumption that the eight 2024 components map to the same physical structures as those identified in 2000–2001. No explicit matching criteria (e.g., velocity or ionization proximity thresholds) or cross-validation tests (such as refitting the archival spectra with the 2024 component set or vice versa and reporting ΔC-stat) are described; without these, differences could arise from improved resolution/S/N or modeling choices rather than time evolution.

Authors: We agree that the physical correspondence between epochs is central to the evolution interpretation and that the manuscript would benefit from greater transparency on this point. Components were matched by aligning their locations in the log ξ–v_out plane with the ranges and specific values reported for the 2000–2001 epoch in the literature; the dominant UTA component (B3) in particular shows nearly identical ionization while its column density triples, and several low-ionization high-velocity components present earlier are absent, which is difficult to attribute solely to resolution or modeling differences. To strengthen the section we will add explicit proximity thresholds used for matching and a short discussion of why resolution/S/N effects alone are unlikely to produce the observed pattern of changes. A full cross-validation refit of the archival spectra with the 2024 component set is a valuable suggestion but requires re-reduction of the older data; we note this as a planned follow-up rather than part of the current revision. revision: partial
Referee: [Results section on fitted parameters] Results section on fitted parameters: The abstract and text report parameter ranges but provide no fit statistics (C-stat or χ²/dof), 1σ error bars on log ξ, N_H, v_out, or v_turb for the eight components, or discussion of degeneracies among them. This undermines quantitative claims such as the factor-of-three column increase in Comp. B3 and the overall 1.5× total N_H rise, as it is impossible to assess whether the reported changes exceed uncertainties or are driven by component splitting.

Authors: We accept that the results section would be clearer with a consolidated presentation of the fit quality and uncertainties. The joint fit was performed with the pion model in SPEX and the best-fit values are stated in the text, but we will add a dedicated table listing all eight components with their best-fit log ξ, N_H, v_out, v_turb, the associated 1σ uncertainties, the overall C-stat/dof of the joint XMM-Newton + XRISM fit, and a brief paragraph on the main degeneracies (primarily between N_H and log ξ in the intermediate-ionization components). This will allow readers to evaluate whether the reported factor-of-three change in component B3 and the 1.5× total-column increase exceed the uncertainties. revision: yes
Referee: [Section on total column density and replenishment] Section on total column density and replenishment: The inference that fresh material is required because the 2024 total column is 1.5 times larger assumes direct additivity and comparability of the summed N_H values across epochs. The manuscript does not address whether the increase could result from the choice to fit eight components (versus fewer in prior work) or from unmodeled continuum or emission features that XRISM resolves but earlier data did not.

Authors: The total column is obtained by summing the individual component columns, which is the conventional procedure for multi-phase warm absorbers and was applied consistently with the modeling approach used in the earlier epoch. The increase is driven primarily by the tripling of the well-constrained dominant component rather than by the addition of extra components. We will revise the section to (i) state the number of components required for an acceptable fit in each epoch, (ii) note that the continuum and emission modeling were cross-checked against the higher-resolution XRISM data, and (iii) explicitly discuss why the observed column increase is unlikely to be an artifact of component splitting or unresolved features. revision: yes

Circularity Check

0 steps flagged

No significant circularity in observational spectral fitting and epoch comparison

full rationale

The paper conducts joint spectral modeling of XMM-Newton/RGS and XRISM/Resolve data using the public pion code in SPEX to extract ionization parameters, column densities, velocities, and AMD for eight WA components. These fitted values are compared directly to independent 2000-2001 observations without any derivation step that reduces a claimed prediction or result back to the input assumptions or self-citations by construction. Component population differences and total column increase are presented as empirical findings from the fits, not as outputs forced by prior self-referential equations or ansatzes. The analysis remains self-contained against external benchmarks and public codes.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 0 invented entities

The central claim depends on multiple fitted parameters for eight discrete absorption components and on the assumption that the pion photoionization model correctly decomposes the observed spectrum into physically distinct zones.

free parameters (3)

log ξ for each of eight components
Fitted to reproduce observed absorption line strengths and ionization states.
column density N_H for each component
Fitted per component; total reported as 1.5 times higher than prior epoch.
outflow velocity and turbulent velocity
Derived from line centroids and widths for each component.

axioms (2)

domain assumption The warm absorber gas is in photoionization equilibrium
Invoked by use of the pion code to model the spectra.
domain assumption The spectrum can be decomposed into a small number of discrete, uniform absorption components
Standard modeling choice when fitting multiple ionization and velocity zones.

pith-pipeline@v0.9.0 · 5667 in / 1498 out tokens · 71354 ms · 2026-05-09T20:47:50.847692+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

45 extracted references · 1 canonical work pages

[1]

P ., Ró˙ za´nska, A., Hryniewicz, K., Czerny, B., & Behar, E

Adhikari, T. P ., Ró˙ za´nska, A., Hryniewicz, K., Czerny, B., & Behar, E. 2019, ApJ, 881, 78

2019
[2]

2009, ApJ, 703, 1346

Behar, E. 2009, ApJ, 703, 1346

2009
[3]

P ., Blustin, A

Behar, E., Rasmussen, A. P ., Blustin, A. J., et al. 2003, ApJ, 598, 232

2003
[4]

Behar, E., Sako, M., & Kahn, S. M. 2001, ApJ, 563, 497

2001
[5]

C., Williams, P

Bentz, M. C., Williams, P . R., Street, R., et al. 2021, ApJ, 920, 112

2021
[6]

J., Branduardi-Raymont, G., Behar, E., et al

Blustin, A. J., Branduardi-Raymont, G., Behar, E., et al. 2002, A&A, 392, 453

2002
[7]

J., Page, M

Blustin, A. J., Page, M. J., Fuerst, S. V ., Branduardi-Raymont, G., & Ashton, C. E. 2005, A&A, 431, 111 den Herder, J. W., Brinkman, A. C., Kahn, S. M., et al. 2001, A&A, 365, L7

2005
[8]

Fabian, A. C. 2012, ARA&A, 50, 455

2012
[9]

R., Crenshaw, D

Gabel, J. R., Crenshaw, D. M., Kraemer, S. B., et al. 2003, ApJ, 595, 120

2003
[10]

& S ˛ adowski, A

Gaspari, M. & S ˛ adowski, A. 2017, ApJ, 837, 149

2017
[11]

2025, A&A, 704, A146

Gu, L., Fukumura, K., Kaastra, J., et al. 2025, A&A, 704, A146

2025
[12]

2023, A&A, 679, A43

Gu, L., Kaastra, J., Rogantini, D., et al. 2023, A&A, 679, A43

2023
[13]

2012, ApJ, 756, L8

Gupta, A., Mathur, S., Krongold, Y ., Nicastro, F., & Galeazzi, M. 2012, ApJ, 756, L8

2012
[14]

2007, ApJ, 663, 799

Holczer, T., Behar, E., & Kaspi, S. 2007, ApJ, 663, 799

2007
[15]

L., Awaki, H., et al

Ishisaki, Y ., Kelley, R. L., Awaki, H., et al. 2025, Journal of Astronomical Tele- scopes, Instruments, and Systems, 11, 042023

2025
[16]

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. Y amashita & T. Watanabe, 411–414

1996
[17]

Abbott, T

Kaastra, J. S., Raassen, A. J. J., de Plaa, J., & Gu, L. 2024, SPEX X-ray spectral ﬁtting package v3.08.01, Zenodo, https://doi.org/10.5281/zenodo. 12771915

work page doi:10.5281/zenodo 2024
[18]

N., George, I

Kaspi, S., Brandt, W. N., George, I. M., et al. 2002, ApJ, 574, 643

2002
[19]

N., Netzer, H., et al

Kaspi, S., Brandt, W. N., Netzer, H., et al. 2001, ApJ, 554, 216

2001
[20]

L., Ishisaki, Y ., Costantini, E., et al

Kelley, R. L., Ishisaki, Y ., Costantini, E., et al. 2025, Journal of Astronomical

2025
[21]

& Pounds, K

King, A. & Pounds, K. 2015, ARA&A, 53, 115

2015
[22]

B., George, I

Kraemer, S. B., George, I. M., Crenshaw, D. M., et al. 2005, ApJ, 633, 693

2005
[23]

S., et al

Krongold, Y ., Nicastro, F., Brickhouse, N. S., et al. 2003, ApJ, 597, 832

2003
[24]

2007, ApJ, 659, 1022

Krongold, Y ., Nicastro, F., Elvis, M., et al. 2007, ApJ, 659, 1022

2007
[25]

S., Gu, L., & Mehdipour, M

Li, C., Kaastra, J. S., Gu, L., & Mehdipour, M. 2023, A&A, 680, A44

2023
[26]

S., Gu, L., et al

Li, C., Kaastra, J. S., Gu, L., et al. 2026, A&A, 706, A255

2026
[27]

S., Gu, L., et al

Li, C., Kaastra, J. S., Gu, L., et al. 2025, A&A, 694, A302

2025
[28]

Lodders, K., Palme, H., & Gail, H. P . 2009, Landolt Börnstein, 4B, 712

2009
[29]

S., et al

Mao, J., Mehdipour, M., Kaastra, J. S., et al. 2019, A&A, 621, A99

2019
[30]

L., Mathur, S., Das, S., Krongold, Y ., & Gupta, A

McClain, R. L., Mathur, S., Das, S., Krongold, Y ., & Gupta, A. 2024, MNRAS, 527, 5093

2024
[31]

S., Eckart, M

Mehdipour, M., Kaastra, J. S., Eckart, M. E., et al. 2025, A&A, 699, A228

2025
[32]

S., & Kallman, T

Mehdipour, M., Kaastra, J. S., & Kallman, T. 2016, A&A, 596, A65

2016
[33]

S., Kriss, G

Mehdipour, M., Kaastra, J. S., Kriss, G. A., et al. 2017, A&A, 607, A28

2017
[34]

M., Kaastra, J

Miller, J. M., Kaastra, J. S., Miller, M. C., et al. 2015, Nature, 526, 542

2015
[35]

M., Lockman, F

Murphy, E. M., Lockman, F. J., Laor, A., & Elvis, M. 1996, ApJS, 105, 369

1996
[36]

2013, The Physics and Evolution of Active Galactic Nuclei

Netzer, H. 2013, The Physics and Evolution of Active Galactic Nuclei

2013
[37]

2003, ApJ, 599, 933

Netzer, H., Kaspi, S., Behar, E., et al. 2003, ApJ, 599, 933

2003
[38]

& Done, C

Noda, H. & Done, C. 2018, MNRAS, 480, 3898

2018
[39]

& Kallman, T

Proga, D. & Kallman, T. R. 2004, ApJ, 616, 688

2004
[40]

R., Kallman, T

Sadaula, D. R., Kallman, T. R., & Laha, S. 2025, ApJ, 992, 182

2025
[41]

M., Behar, E., et al

Sako, M., Kahn, S. M., Behar, E., et al. 2001, A&A, 365, L168

2001
[42]

2014, MNRAS, 445, 3011 Strüder, L., Briel, U., Dennerl, K., et al

Stern, J., Behar, E., Laor, A., Baskin, A., & Holczer, T. 2014, MNRAS, 445, 3011 Strüder, L., Briel, U., Dennerl, K., et al. 2001, A&A, 365, L18

2014
[43]

1998, A&AS, 130, 333

Theureau, G., Bottinelli, L., Coudreau-Durand, N., et al. 1998, A&AS, 130, 333

1998
[44]

Turner, M. J. L., Abbey, A., Arnaud, M., et al. 2001, A&A, 365, L27

2001
[45]

J., Reeves, J

Turner, T. J., Reeves, J. N., Kraemer, S. B., & Miller, L. 2008, A&A, 483, 161 Xrism Collaboration, Audard, M., Awaki, H., et al. 2025, A&A, 702, A147 Article number, page 9

2008