Cosmic-ray impact on optical and mid-infrared emission line diagnostics in NGC 5728
Pith reviewed 2026-05-21 22:16 UTC · model grok-4.3
The pith
Cosmic rays heat inner gas clouds in NGC 5728 and enhance low-ionization optical emission lines.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Cosmic rays serve as a significant feedback mechanism that heats the interstellar medium and modifies its chemical composition in narrow line region clouds. The study finds that cosmic rays are instrumental in heating the inner regions of gas clouds, enhancing emission of low-ionization optical lines. Mid-infrared data reveal that emission lines like [Ar II] and [Ne II] within the JWST MIRI field of view are sensitive to cosmic rays. In contrast, high-ionization lines such as [Ne V] serve as robust tracers of photoionization insensitive to cosmic rays. Mixed optical and mid-infrared diagnostics help distinguish the relative roles of cosmic rays and shocks.
What carries the argument
CLOUDY photoionization models that vary cosmic ray ionization rate, ionization parameter, and initial hydrogen density to reproduce observed line ratios from the narrow line region of NGC 5728.
If this is right
- Low-ionization mid-infrared lines respond to cosmic ray heating while high-ionization lines do not.
- Mixed optical and mid-infrared diagnostics allow separation of cosmic ray and shock contributions.
- The method resolves the degeneracy between metallicity and cosmic ray ionization in emission line analysis.
- High-ionization lines remain useful tracers of photoionization even when cosmic rays are present.
Where Pith is reading between the lines
- The same modeling approach could be tested on other active galaxies to see how cosmic ray influence scales with nuclear activity.
- Including cosmic rays in standard narrow line region models might change inferred gas densities or metallicities in AGN hosts.
- Future multi-wavelength surveys could use these line ratios to map cosmic ray feedback across galaxy centers.
Load-bearing premise
The CLOUDY models with varied cosmic ray ionization rate, ionization parameter, and initial hydrogen density capture the dominant physical processes in the narrow line region clouds without large unmodeled contributions from detailed shock microphysics or geometry effects.
What would settle it
Observing no increase in the ratios of low-ionization lines such as [Ne II] relative to high-ionization lines like [Ne V] when cosmic ray ionization rate is raised in the models would falsify the central claim.
Figures
read the original abstract
Cosmic rays (CRs), from active galactic nuclei (AGN) jets and supernovae (SNe), serve as a significant feedback mechanism influencing emission lines in narrow line region (NLR) clouds. These highly energetic particles, propelled by shocks, heat the interstellar medium (ISM) and modify its chemical composition. This study investigates the role of CRs, particularly in their ability to excite gas and align with observed line ratios across UV and optical diagnostics. We employ CLOUDY to explore CR ionization rate, ionization parameter, and initial hydrogen density effects on optical and mid-infrared (MIR) emission. Our analysis includes high-quality optical data from the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) for NGC 5728, supplemented by infrared observations from the James Webb Space Telescope (JWST). Our previous results indicate that CRs are instrumental in heating the inner regions of gas clouds, enhancing emission of low-ionization optical lines. Mid-infrared data reveal that emission lines like [Ar II] and [Ne II] within the JWST Mid-Infrared Instrument (MIRI) field of view are sensitive to CRs. In contrast, high-ionization lines (for example, [Ne V]) serve as robust tracers of photoionization insensitive to CRs. Moreover, mixed optical and MIR diagnostics offer insight into the relative roles of CRs and shocks, which often produce similar signatures in emission lines. We find that while both mechanisms can elevate certain line ratios, their influence on MIR diagnostics diverges: shocks and CRs affect low-ionization lines differently, allowing for a better understanding when multi-wavelength data are available. Our approach not only helps to resolve the degeneracy between metallicity and CR ionization but also enables the potential differentiation of shocks and CR-driven processes in AGN.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates the effects of cosmic rays on optical and mid-infrared emission line diagnostics in the narrow line region of NGC 5728. Using CLOUDY photoionization models that vary cosmic-ray ionization rate, ionization parameter, and initial hydrogen density, the authors compare results to MUSE optical spectra and JWST MIRI observations. They conclude that cosmic rays heat inner cloud regions and enhance low-ionization optical lines, that MIR lines such as [Ar II] and [Ne II] are sensitive to CRs while high-ionization lines like [Ne V] are robust photoionization tracers, and that combined optical-MIR diagnostics can distinguish CR-driven processes from shocks despite similar optical signatures.
Significance. If the central modeling results hold under more complete physical treatments, the work would offer a practical multi-wavelength strategy for breaking degeneracies in AGN NLR diagnostics between cosmic-ray ionization, shocks, and photoionization. The direct use of high-quality MUSE and JWST data on a specific target provides a concrete test case for these ideas.
major comments (2)
- Abstract: the modeling strategy and qualitative conclusions are outlined, but no quantitative fit statistics, error budgets, or explicit model-data comparison metrics (e.g., predicted vs. observed line ratios with uncertainties) are supplied. This prevents verification that the data support the stated line sensitivities and is load-bearing for the central claims about CR effects and diagnostic utility.
- Modeling approach (as described in the abstract): the CLOUDY grids vary only CR ionization rate, ionization parameter, and initial hydrogen density. The claim that MIR diagnostics diverge between CRs and shocks (allowing differentiation when optical signatures are similar) requires either explicit inclusion of shock microphysics (e.g., via CLOUDY shock options or MAPPINGS-style prescriptions) or a demonstration that unmodeled effects such as magnetic compression, non-equilibrium ionization, or NLR geometry do not alter low-ionization MIR ratios in ways that undermine the separation. This is load-bearing for the differentiation result.
minor comments (1)
- The abstract references 'our previous results' without a specific citation or quantitative description of how the current grids extend or differ from prior work.
Simulated Author's Rebuttal
We thank the referee for their constructive and detailed report. Below we respond point-by-point to the major comments, indicating where revisions will be made to address the concerns while preserving the core scientific approach of the manuscript.
read point-by-point responses
-
Referee: Abstract: the modeling strategy and qualitative conclusions are outlined, but no quantitative fit statistics, error budgets, or explicit model-data comparison metrics (e.g., predicted vs. observed line ratios with uncertainties) are supplied. This prevents verification that the data support the stated line sensitivities and is load-bearing for the central claims about CR effects and diagnostic utility.
Authors: We agree that the abstract is currently qualitative and would benefit from explicit quantitative anchors. The main text already presents direct model-observation comparisons for key lines (e.g., [Ne II], [Ar II], [Ne V]) with associated uncertainties derived from the MUSE and JWST data, as well as goodness-of-fit assessments across the CLOUDY grid. In the revised manuscript we will condense the most diagnostic quantitative results (specific line-ratio offsets and their uncertainties) into the abstract to make the support for the claimed sensitivities immediately verifiable. revision: yes
-
Referee: Modeling approach (as described in the abstract): the CLOUDY grids vary only CR ionization rate, ionization parameter, and initial hydrogen density. The claim that MIR diagnostics diverge between CRs and shocks (allowing differentiation when optical signatures are similar) requires either explicit inclusion of shock microphysics (e.g., via CLOUDY shock options or MAPPINGS-style prescriptions) or a demonstration that unmodeled effects such as magnetic compression, non-equilibrium ionization, or NLR geometry do not alter low-ionization MIR ratios in ways that undermine the separation. This is load-bearing for the differentiation result.
Authors: We acknowledge that the present grids isolate cosmic-ray effects within a pure photoionization framework and do not run dedicated shock models. The stated divergence rests on comparing our CR-enhanced predictions against published shock-model libraries (primarily MAPPINGS) that already show distinct MIR low-ionization behavior under shock heating. To address the concern directly, the revised manuscript will include a new subsection that (i) tabulates the relevant literature shock predictions for the same MIR lines, (ii) discusses the expected impact of unmodeled factors such as NLR geometry and non-equilibrium ionization on the separation, and (iii) presents a limited sensitivity test varying cloud thickness and density profile within CLOUDY. A full self-consistent shock-plus-CR calculation lies outside the current scope but can be noted as future work. revision: partial
Circularity Check
Minor self-citation to prior team results; central claims rest on independent CLOUDY forward modeling matched to external MUSE and JWST observations
specific steps
-
self citation load bearing
[Abstract]
"Our previous results indicate that CRs are instrumental in heating the inner regions of gas clouds, enhancing emission of low-ionization optical lines."
This sentence invokes prior results from the same team to support the role of CRs. While not the central derivation (which relies on new CLOUDY runs vs. MUSE/JWST data), it constitutes a self-citation that is not independently verified within the current manuscript.
full rationale
The paper's core analysis employs CLOUDY grids varying CR ionization rate, ionization parameter, and n_H to interpret observed line ratios from independent MUSE optical and JWST MIRI data on NGC 5728. No self-definitional loops, fitted inputs renamed as predictions, or ansatzes smuggled via citation are present in the provided text. The sole self-reference appears in the abstract as a supporting statement rather than a load-bearing premise for the new multi-wavelength diagnostics. This qualifies as a minor self-citation (score 2) without reducing the derivation to its own inputs. The modeling is compared against external benchmarks, keeping the work self-contained.
Axiom & Free-Parameter Ledger
free parameters (3)
- CR ionization rate
- ionization parameter
- initial hydrogen density
axioms (2)
- domain assumption CLOUDY correctly incorporates cosmic-ray ionization and heating physics
- domain assumption Observed line ratios primarily reflect the modeled CR, density, and ionization-parameter variations
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We employ CLOUDY simulations to explore CR ionization rate, ionization parameter, and initial hydrogen density effects on optical and mid-infrared (MIR) emission.
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
CRs are instrumental in heating the inner regions of gas clouds, enhancing emission of low-ionization optical lines.
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
, " * write output.state after.block = add.period write newline
ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 ...
-
[2]
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....
-
[3]
Aalto , S., Garcia-Burillo , S., Muller , S., et al. 2015, , 574, A85
work page 2015
-
[4]
Abazajian , K. N., Adelman-McCarthy , J. K., Ag \"u eros , M. A., et al. 2009, , 182, 543
work page 2009
-
[5]
Allen , M. G., Groves , B. A., Dopita , M. A., Sutherland , R. S., & Kewley , L. J. 2008, , 178, 20
work page 2008
-
[6]
Alonso Herrero , A., Hermosa Mu \ n oz , L., Labiano , A., et al. 2024, , 690, A95
work page 2024
-
[7]
\'A lvarez-M \'a rquez , J., Labiano , A., Guillard , P., et al. 2023, , 672, A108
work page 2023
- [8]
-
[9]
Asplund , M., Grevesse , N., Sauval , A. J., & Scott , P. 2009, , 47, 481
work page 2009
-
[10]
Bacon , R., Piqueras , L., Conseil , S., Richard , J., & Shepherd , M. 2016, MPDAF: MUSE Python Data Analysis Framework , Astrophysics Source Code Library, record ascl:1611.003
work page 2016
-
[11]
Bacon, Roland , Brinchmann, Jarle , Conseil, Simon , et al. 2023, A&A, 670, A4
work page 2023
-
[12]
Baldi , R. D. 2023, , 31, 3
work page 2023
- [13]
- [14]
-
[15]
Blanc , G. A., Kewley , L., Vogt , F. P. A., & Dopita , M. A. 2015, , 798, 99
work page 2015
- [16]
-
[17]
2024, astropy/photutils: 2.0.2
Bradley, L., Sip o cz, B., Robitaille, T., et al. 2024, astropy/photutils: 2.0.2
work page 2024
-
[18]
2013, Monthly Notices of the Royal Astronomical Society, 432, 2112
Brinchmann, J., Charlot, S., Kauffmann, G., et al. 2013, Monthly Notices of the Royal Astronomical Society, 432, 2112
work page 2013
-
[19]
Brinchmann, J., Charlot, S., White, S. D. M., et al. 2004, Monthly Notices of the Royal Astronomical Society, 351, 1151
work page 2004
- [20]
-
[21]
2001, Publications of the ASP, 113, 1449
Calzetti , D. 2001, Publications of the ASP, 113, 1449
work page 2001
- [22]
-
[23]
Davies , R., Shimizu , T., Pereira-Santaella , M., et al. 2024, , 689, A263
work page 2024
-
[24]
Dopita , M. A. & Sutherland , R. S. 1995, , 455, 468
work page 1995
-
[25]
Dors , O. L. 2021, , 507, 466
work page 2021
-
[26]
L., Valerdi, M., Freitas-Lemes, P., et al
Dors, O. L., Valerdi, M., Freitas-Lemes, P., et al. 2022, Monthly Notices of the Royal Astronomical Society, 514, 5506–5527
work page 2022
-
[27]
P., Satyapal , S., & Marcu , D
Dudik , R. P., Satyapal , S., & Marcu , D. 2009, , 691, 1501
work page 2009
-
[28]
Dudik , R. P., Weingartner , J. C., Satyapal , S., et al. 2007, , 664, 71
work page 2007
- [29]
- [30]
-
[31]
Feltre , A., Gruppioni , C., Marchetti , L., et al. 2023, , 675, A74
work page 2023
-
[32]
J., Chatzikos , M., Guzm \'a n , F., et al
Ferland , G. J., Chatzikos , M., Guzm \'a n , F., et al. 2017, Revista Mexicana de Astronomia y Astrofisica, 53, 385
work page 2017
-
[33]
J., Done, C., Jin, C., Landt, H., & Ward, M
Ferland, G. J., Done, C., Jin, C., Landt, H., & Ward, M. J. 2020, Monthly Notices of the Royal Astronomical Society, 494, 5917
work page 2020
-
[34]
Ferland , G. J., Fabian , A. C., Hatch , N. A., et al. 2009, , 392, 1475
work page 2009
-
[35]
Ferland , G. J., Porter , R. L., van Hoof , P. A. M., et al. 2013, , 49, 137
work page 2013
-
[36]
A., P \'e rez-Montero , E., V \' lchez , J
Fern \'a ndez-Ontiveros , J. A., P \'e rez-Montero , E., V \' lchez , J. M., Amor \' n , R., & Spinoglio , L. 2021, , 652, A23
work page 2021
-
[37]
A., Spinoglio , L., Pereira-Santaella , M., et al
Fern \'a ndez-Ontiveros , J. A., Spinoglio , L., Pereira-Santaella , M., et al. 2016, , 226, 19
work page 2016
-
[38]
Fern \'a ndez-Ontiveros, J. A. et al. 2025, in preparation
work page 2025
- [39]
-
[40]
1998, The Astrophysical Journal, 498, 579
Genzel, R., Lutz, D., Sturm, E., et al. 1998, The Astrophysical Journal, 498, 579
work page 1998
-
[41]
Gonz \'a lez-Alfonso , E., Fischer , J., Bruderer , S., et al. 2018, , 857, 66
work page 2018
-
[42]
Gonz \'a lez-Alfonso , E., Fischer , J., Bruderer , S., et al. 2013, , 550, A25
work page 2013
-
[43]
2024, The Astrophysical Journal, 966, 204
Goold, K., Seth, A., Molina, M., et al. 2024, The Astrophysical Journal, 966, 204
work page 2024
-
[44]
Goold, K. S. et al. 2025, in preparation
work page 2025
-
[45]
Goulding , A. D. & Alexander , D. M. 2009, , 398, 1165
work page 2009
- [46]
- [47]
-
[48]
2023, arXiv e-prints, arXiv:2304.03726
Hatano , S., Ouchi , M., Nakajima , K., et al. 2023, arXiv e-prints, arXiv:2304.03726
- [49]
-
[50]
Hermosa Mu \ n oz , L., Alonso-Herrero , A., Labiano , A., et al. 2025, , 693, A321
work page 2025
- [51]
-
[52]
2017, Monthly Notices of the Royal Astronomical Society, 472, 2468
Hirschmann, M., Charlot, S., Feltre, A., et al. 2017, Monthly Notices of the Royal Astronomical Society, 472, 2468
work page 2017
- [53]
- [54]
-
[55]
Kauffmann , G., Heckman , T. M., Tremonti , C., et al. 2003, , 346, 1055
work page 2003
-
[56]
J., Groves , B., Kauffmann , G., & Heckman , T
Kewley , L. J., Groves , B., Kauffmann , G., & Heckman , T. 2006, , 372, 961
work page 2006
-
[57]
Koutsoumpou , E., Fern \'a ndez-Ontiveros , J. A., Dasyra , K. M., & Spinoglio , L. 2025, , 693, A215
work page 2025
-
[58]
Lumbreras-Calle , A., Fern \'a ndez-Ontiveros , J. A., Infante-Sainz , R., et al. 2024, arXiv e-prints, arXiv:2412.08327
-
[59]
Luridiana, V., Morisset, C., & Shaw, R. A. 2014, Astronomy & Astrophysics, 573, A42
work page 2014
- [60]
-
[61]
McKee , C. F. 1989, , 345, 782
work page 1989
- [62]
- [63]
-
[64]
Newville, M., Stensitzki, T., Allen, D. B., & Ingargiola, A. 2015, LMFIT: Non-Linear Least-Square Minimization and Curve-Fitting for Python
work page 2015
-
[65]
Oliveira , C. B., Krabbe , A. C., Dors , O. L., et al. 2024, , 531, 199
work page 2024
-
[66]
Padovani , M., Galli , D., & Glassgold , A. E. 2009, , 501, 619
work page 2009
-
[67]
Padovani , M., Ivlev , A. V., Galli , D., & Caselli , P. 2018, , 614, A111
work page 2018
- [68]
-
[69]
Pereira-Santaella , M., \'A lvarez-M \'a rquez , J., Garc \' a-Bernete , I., et al. 2022, , 665, L11
work page 2022
-
[70]
M., Alonso-Herrero , A., & Rieke , G
Pereira-Santaella , M., Diamond-Stanic , A. M., Alonso-Herrero , A., & Rieke , G. H. 2010, , 725, 2270
work page 2010
-
[71]
P \'e rez-D \' az , B., P \'e rez-Montero , E., Fern \'a ndez-Ontiveros , J. A., & V \' lchez , J. M. 2022, , 666, A115
work page 2022
- [72]
- [73]
- [74]
-
[75]
P \'e rez-Montero , E., Dors , O. L., V \' lchez , J. M., et al. 2019, , 489, 2652
work page 2019
-
[76]
A., P \'e rez-D \' az , B., et al
P \'e rez-Montero , E., Fern \'a ndez-Ontiveros , J. A., P \'e rez-D \' az , B., et al. 2024, , 684, A40
work page 2024
-
[77]
2021, Monthly Notices of the Royal Astronomical Society, 505, 4289–4309
Pérez-Díaz, B., Masegosa, J., Márquez, I., & Pérez-Montero, E. 2021, Monthly Notices of the Royal Astronomical Society, 505, 4289–4309
work page 2021
-
[78]
T., Simpson , C., Polimera , M
Richardson , C. T., Simpson , C., Polimera , M. S., et al. 2022, , 927, 165
work page 2022
-
[79]
Rybicki , G. B. & Lightman , A. P. 1979, Radiative processes in astrophysics ( John Wiley & Sons )
work page 1979
-
[80]
2007, Monthly Notices of the Royal Astronomical Society, 382, 1415
Schawinski, K., Thomas, D., Sarzi, M., et al. 2007, Monthly Notices of the Royal Astronomical Society, 382, 1415
work page 2007
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.