Abstract:Cosmic rays (CRs) can significantly impact dense molecular clouds in galaxies, heating the interstellar medium (ISM) and altering its chemistry, ionization, and thermal properties. Their influence is particularly relevant in environments with high CR rates, such as starburst galaxies with supernova remnants or jets and outflows in active galactic nuclei (AGN). CRs transfer energy to the ionized phase of the ISM far from the ionization source, preventing gas cooling and driving large-scale winds. In this work, we use CLOUDY to explore the effect of CRs on nebular gas, a relatively underexplored area, mainly focused on cold molecular gas. Our models cover a broad range of density ($1 - 10^4\,\rm{cm^{-3}}$), ionization parameter ($-3.5 \leq \log U \leq -1.5$), and CR ionization rate ($10^{-16}\, \rm{s^{-1}} - 10^{-12}\, \rm{s^{-1}}$). These are compared to MUSE observations of two AGN, Centaurus A and NGC 1068, and the starburst NGC 253. We find that CR rates $\gtrsim 10^{-13}\, \rm{s^{-1}}$, typical of AGN and strong starburst galaxies, can significantly alter the thermal structure of the ionized gas by forming a deep secondary low-ionization layer beyond the photoionization-dominated region. This enhances emission from low-ionization transitions, such as [\ion{N}{ii}], [\ion{S}{ii}], and [\ion{O}{i}], affecting line-ratio diagnostics, metallicity, and ionization estimates. Unlike pure photoionization models, AGN simulations with high CR ionization rates reproduce the Seyfert loci in BPT diagrams without requiring super-solar metallicities for the narrow-line region. Additionally, star formation simulations with high CR ionization rates can explain line ratios in the LINER domain. We propose new maximum starburst boundaries for BPT diagrams to distinguish regions dominated by AGN photoionization from those that could be explained by star formation plus high CR ionization rates.

Evidence for Radiation Pressure Compression in the X-ray Narrow Line Region of Seyfert galaxies

Radiation pressure confinement -- V. The predicted free-free absorption and emission in active galactic nuclei

Outflow densities and ionisation mechanisms in the NLRs of the prototypical Seyfert galaxies NGC 1068 and NGC 4151

Accretion properties of X-ray AGN: Evidence for radiation-regulated obscuration with redshift-dependent host galaxy contribution

Radiation Pressure, Absorption and AGN Feedback in the Chandra Deep Fields

Evidence for a dusty warm absorber in NGC 3227 ?

The soft X-ray/NLR connection: a single photoionized medium?

X-ray imaging of the ionisation cones in NGC 5252

Cosmic Ray Ionization of Low-Excitation Lines in Active Galactic Nuclei and Starburst Galaxies

Soft X-ray emission lines from a relativistic accretion disk in MCG -6-30-15 and Mrk 766

Evidence for Relativistic Outflows in Narrow-line Seyfert 1 Galaxies

The Soft X-Ray Properties Of Nearby Low-Luminosity Active Galactic Nuclei And Their Contribution To The Cosmic X-Ray Background

EVIDENCE FOR COSMIC RAY ACCELERATION IN CASSIOPEIA A

Ray-tracing simulations of the Soft X-ray Scattered Emission in obscured Active Galactic Nuclei

Do gas clouds in narrow-line regions of Seyfert galaxies come from their nuclei?

The role of absorption and reflection in the soft X-ray excess of Active Galactic Nuclei : 1. Preliminary results

Radial X-ray profiles of simulated galaxies: Contributions from hot gas and XRBs

Properties of dusty warm absorbers and the case of IRAS 17020+4544

The Shape of the Cosmic X-ray Background: Nuclear Starburst Discs and the Redshift Evolution of AGN Obscuration

Exploring the high-density reflection model for the soft excess in RBS 1124