pith. sign in

arxiv: 2605.30477 · v1 · pith:NBREEXZMnew · submitted 2026-05-28 · 🌌 astro-ph.GA

Simultaneous modeling of FeII emission in the optical and near-infrared in a prototypical Narrow-Line Seyfert 1 galaxy

Denimara Dias dos Santos , Swayamtrupta Panda , Alberto Rodr\'iguez-Ardila , Murilo Marinello This is my paper

Pith reviewed 2026-06-29 06:13 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords FeII emissionActive galactic nucleiBroad-line regionPhotoionization modelingNarrow-Line Seyfert 1 galaxyLyman-alpha fluorescenceNear-infrared spectroscopy
0
0 comments X

The pith

Optical and near-infrared FeII emission in AGN can be simultaneously reproduced with the same gas density and metallicity in the broad-line region.

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

The paper combines optical and near-infrared spectra of the galaxy IZw1 with photoionization models to find conditions that fit both FeII emission bands at once. Updated atomic data sets allow the models to match the observed line ratios when hydrogen density sits between 10^11 and 10^12 cm^{-3} and metallicity is near solar. Lyman-alpha fluorescence is shown to drive most of the near-infrared emission while collisions matter more in the optical. The same models also link FeII to OI and the CaII triplet, pointing to overlapping zones inside the broad-line region.

Core claim

Optical R_4570 and NIR R_1μm FeII emission can be simultaneously reproduced under consistent conditions, with the best agreement obtained using the Smyth et al. (2019) dataset, for hydrogen densities of 10^{11.0} to 10^{12.0} cm^{-3} and near-solar metallicity. Lyman-α fluorescence dominates the NIR while collisional processes dominate the optical. NIR FeII and OI arise in less dense gas than optical FeII and CaT, supporting a vertical structure in the broad-line region.

What carries the argument

Photoionization modeling with updated FeII atomic datasets that explores a wide parameter space of density, metallicity, and radiation field to match observed line ratios in both wavelength regimes.

If this is right

  • NIR FeII and OI trace less dense gas than optical FeII and the CaII triplet inside the same clouds.
  • Lyman-alpha fluorescence must be included to set the correct physical conditions for NIR FeII.
  • CaT and OI lines can serve as proxies for optical FeII strength under the same conditions.
  • A vertical density gradient in the broad-line region is consistent with the differing line ratios.

Where Pith is reading between the lines

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

  • If the single-zone assumption holds, time-series observations could map how density changes with radius in the broad-line region.
  • The fluorescence dominance in the NIR suggests that UV continuum variability will affect infrared FeII more strongly than optical FeII.
  • Extending the same modeling grid to other Seyfert galaxies would test whether the derived density window is universal.

Load-bearing premise

The FeII-emitting gas can be represented by a single set of physical conditions without detailed treatment of the full geometry or multiple cloud components in the broad-line region.

What would settle it

Spectra of additional AGN in which no single density-metallicity pair reproduces both the optical and near-infrared FeII ratios within the observed uncertainties.

read the original abstract

This work investigates the FeII emission in active galactic nuclei (AGN), combining observational data from optical and near-infrared (NIR) spectra of the prototypical FeII emitter IZw1 with state-of-the-art photoionization modeling. Using updated FeII atomic datasets (Smyth et al. 2019; Tayal & Zatsarinny 2018; Bautista et al. 2015), we explore a wide parameter space to determine the physical conditions of FeII-emitting regions in the broad-line region (BLR). Our results show that optical ($R_{\rm 4570}$) and NIR ($R_{\rm 1\mu m}$) FeII emission can be simultaneously reproduced under consistent conditions, with the best agreement obtained using the Smyth et al. (2019) dataset, for hydrogen densities of $10^{11.0}$ to $10^{12.0}$ cm$^{-3}$ and near-solar metallicity. We quantify, for the first time, the impact of Lyman-$\alpha$ fluorescence on the physical conditions of FeII emission in both regimes, revealing its dominant role in the NIR and, in contrast, highlighting the stronger influence of collisional processes in the optical. Additionally, for the first time, we compare optical and NIR FeII emission simultaneously with OI and the CaII triplet (CaT), reinforcing their connection to similar spatial regions and physical properties, as well as their usefulness as better proxies for optical FeII. Our findings support the idea of a vertical BLR structure, with NIR FeII and OI originating in less dense regions of the cloud than optical FeII and CaT.

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 / 1 minor

Summary. The manuscript investigates FeII emission in the NLS1 galaxy IZw1 by combining optical and NIR spectra with photoionization modeling using three updated FeII atomic datasets. It claims that R_4570 (optical) and R_1μm (NIR) FeII ratios can be simultaneously reproduced under consistent physical conditions (n_H = 10^{11.0}–10^{12.0} cm^{-3}, near-solar metallicity), with the Smyth et al. (2019) dataset providing the best match; Lyman-α fluorescence dominates NIR while collisions dominate optical; and comparisons to OI and CaT support a vertical BLR structure with NIR lines arising in less dense gas.

Significance. If the modeling results hold, the work provides new quantitative constraints on BLR conditions for FeII emission and the first explicit assessment of Lyman-α fluorescence contributions across optical and NIR regimes. The simultaneous treatment of multiple line ratios and atomic datasets adds value for interpreting AGN spectra, though the single-zone framework limits direct applicability to stratified geometries.

major comments (2)
  1. [Abstract] Abstract: The central claim that optical and NIR FeII 'can be simultaneously reproduced under consistent conditions' is load-bearing for the paper's conclusions, yet the modeling uses single-zone photoionization runs; this creates tension with the inference that NIR FeII/OI arise in less dense gas than optical FeII/CaT, as a single set of parameters cannot by construction demonstrate a density gradient without explicit multi-component or stratified geometry tests.
  2. [Abstract] Abstract and modeling description: The parameter-space exploration determines best-fit conditions by matching observed ratios, but the manuscript does not report quantitative fit metrics (e.g., χ² values, residual distributions, or uncertainty ranges on n_H and Z) for the simultaneous optical+NIR fits; without these, it is unclear whether the reported density range uniquely supports the 'consistent conditions' result or is an average that would shift under multi-cloud modeling.
minor comments (1)
  1. [Abstract] The abstract states 'for the first time' for both the fluorescence quantification and the OI/CaT comparison; these novelty claims should be supported by explicit references to prior simultaneous optical-NIR FeII studies in the introduction.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough and constructive review. We address the two major comments below. We agree that quantitative fit metrics should be added and will revise accordingly. On the single-zone modeling, we clarify the scope of our claims without altering the core results.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The central claim that optical and NIR FeII 'can be simultaneously reproduced under consistent conditions' is load-bearing for the paper's conclusions, yet the modeling uses single-zone photoionization runs; this creates tension with the inference that NIR FeII/OI arise in less dense gas than optical FeII/CaT, as a single set of parameters cannot by construction demonstrate a density gradient without explicit multi-component or stratified geometry tests.

    Authors: We thank the referee for this observation. The single-zone calculations demonstrate that there is a region of parameter space (n_H = 10^{11.0}–10^{12.0} cm^{-3}, near-solar Z) where both R_4570 and R_1μm are simultaneously matched to the observed values using the Smyth et al. (2019) dataset. This supports the statement that the two FeII ratios can be reproduced under consistent conditions. The suggestion of a vertical BLR structure with NIR FeII/OI arising in less dense gas is drawn from the separate comparisons to the OI and CaT lines, which indicate that those species are better reproduced at somewhat lower densities than the optical FeII/CaT combination. We do not use the single-zone runs to demonstrate a gradient; the multi-line analysis provides that supporting evidence. We will revise the abstract and relevant discussion sections to make this distinction explicit and remove any potential ambiguity. revision: partial

  2. Referee: [Abstract] Abstract and modeling description: The parameter-space exploration determines best-fit conditions by matching observed ratios, but the manuscript does not report quantitative fit metrics (e.g., χ² values, residual distributions, or uncertainty ranges on n_H and Z) for the simultaneous optical+NIR fits; without these, it is unclear whether the reported density range uniquely supports the 'consistent conditions' result or is an average that would shift under multi-cloud modeling.

    Authors: We agree that the absence of quantitative goodness-of-fit metrics limits the clarity of the results. In the revised manuscript we will add χ² values (or equivalent residuals) for the simultaneous matches to R_4570 and R_1μm across the explored grid, together with the range of n_H and Z that satisfy the observed ratios within a stated tolerance. This will allow readers to assess how unique the reported density interval is and how it might behave under multi-component interpretations. revision: yes

Circularity Check

0 steps flagged

No circularity: parameter exploration against external atomic data yields direct matches to observed ratios

full rationale

The paper performs standard photoionization modeling by varying density, metallicity, and other parameters in codes fed with independent atomic datasets (Smyth et al. 2019 and others). It identifies the subset of those parameters that simultaneously reproduce the measured R_4570 and R_1μm ratios in IZw1. This is ordinary fitting, not a derivation that reduces to its own inputs by construction. No load-bearing self-citation, no fitted quantity renamed as a first-principles prediction, and no ansatz smuggled via prior work by the same authors. The single-zone assumption is an explicit modeling choice whose limitations are acknowledged, not a hidden tautology.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The claim depends on the validity of the atomic data updates and the photoionization modeling framework, with two main fitted parameters for density and metallicity.

free parameters (2)
  • hydrogen density = 10^{11.0} to 10^{12.0} cm^{-3}
    Determined as best fit to reproduce observed FeII ratios
  • metallicity = near-solar
    Selected for best agreement with observations
axioms (1)
  • domain assumption Standard assumptions in photoionization codes regarding atomic processes and radiative transfer
    Invoked when using the modeling to explore parameter space for BLR conditions

pith-pipeline@v0.9.1-grok · 5855 in / 1525 out tokens · 42154 ms · 2026-06-29T06:13:29.137886+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Little Red Dots as Intermediate Mass, Super-Eddington Engines: Insights from Type IIn Supernovae and The 1837-1856 Great Eruption of $\eta$ Carinae

    astro-ph.GA 2026-06 unverdicted novelty 6.0

    LRDs are reinterpreted as intermediate-mass super-Eddington systems with wind-driven pseudo-photospheres that explain their spectra and imply engine masses below 10^5 solar masses rather than overmassive black holes.

Reference graph

Works this paper leans on

5 extracted references · 5 canonical work pages · cited by 1 Pith paper

  1. [1]

    A., Ferland, G

    Baldwin, J. A., Ferland, G. J., Korista, K. T., Hamann, F., & LaCluyz´ e, A. 2004, ApJ, 615, 610, doi: 10.1086/424683 Bautista, M. A., Fivet, V., Ballance, C., et al. 2015, The Astrophysical Journal, 808, 174, doi: 10.1088/0004-637x/808/2/174 Bentz, M. C., Peterson, B. M., Netzer, H., Pogge, R. W., & Vestergaard, M. 2009, The Astrophysical Journal, 697, 1...

    work page doi:10.1086/424683 2004

  2. [2]

    2025, A&A, 697, A23, doi: 10.1051/0004-6361/202452067 Garcia-Rissmann, A., Rodr´ ıguez-Ardila, A., Sigut, T

    https://arxiv.org/abs/1705.10877 Floris, A., Pandey, A., Czerny, B., et al. 2025, A&A, 697, A23, doi: 10.1051/0004-6361/202452067 Garcia-Rissmann, A., Rodr´ ıguez-Ardila, A., Sigut, T. A. A., & Pradhan, A. K. 2012, ApJ, 751, 7, doi: 10.1088/0004-637X/751/1/7 Garnica, K., Negrete, C. A., Marziani, P., et al. 2022, A&A, 667, A105, doi: 10.1051/0004-6361/202...

    work page doi:10.1051/0004-6361/202452067 2025

  3. [3]

    J., & Scott, P

    https://adsabs.harvard.edu/pdf/1986ApJ...305..175G Grevesse, N., Asplund, M., Sauval, A. J., & Scott, P. 2010, Ap&SS, 328, 179, doi: 10.1007/s10509-010-0288-z Grier, C. J., Pancoast, A., Barth, A. J., et al. 2017, The Astrophysical Journal, 849, 146, doi: 10.3847/1538-4357/aa901b Grier, C. J., Peterson, B. M., Pogge, R. W., et al. 2012, Astrophysical Jour...

    work page doi:10.1007/s10509-010-0288-z 2010

  4. [4]

    S., Netzer, H., et al

    https:// ui.adsabs.harvard.edu/abs/1987A&A...184...33J/abstract —. 1991, Astronomy and Astrophysics, Vol. 242, 49 —. 1993, Annales de Physique, 241 Kaspi, S., Smith, P. S., Netzer, H., et al. 2000, ApJ, 533, 631, doi: 10.1086/308704 Kovaˇ cevi´ c, J., Popovi´ c, L. v., & Dimitrijevi´ c, M. S. 2010, The Astrophysical Journal Supplement Series, 189, 15, doi...

    work page doi:10.1086/308704 1991

  5. [5]

    2017, The Astrophysical Journal Supplement Series, 124, 3042 Richardson, C

    http://adsabs.harvard.edu/full/1988ApJ...330..751P Rakshit, S., Stalin, C., Chand, H., & Zhang, X. 2017, The Astrophysical Journal Supplement Series, 124, 3042 Richardson, C. T., Allen, J. T., Baldwin, J. A., Hewett, P. C., & Ferland, G. J. 2014, MNRAS, 437, 2376, doi: 10.1093/mnras/stt2056 Rodr´ ıguez-Ardila, A., Pastoriza, M. G., & Donzelli, C. J. 2000,...

    work page doi:10.1093/mnras/stt2056 2017