Shock-driven heating in the circumnuclear star-forming regions of NGC 7582: Insights from JWST NIRSpec and MIRI/MRS spectroscopy
Pith reviewed 2026-05-18 03:21 UTC · model grok-4.3
The pith
Shock heating from slow C-type shocks explains elevated molecular temperatures in NGC 7582's star-forming regions
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Combined NIRSpec and MIRI/MRS spectroscopy shows a power-law temperature distribution in the rotational H2 lines, with the southern star-forming regions having unexpectedly high temperatures. Fitting Paris-Durham shock models to the H2 lines indicates that a slow C-type shock with velocity of approximately 10 km/s is responsible for the heating, while PDR Toolbox models give n_H ~ 10^5 cm^{-3} and G_0 ~ 10^3.
What carries the argument
Fitting of Paris-Durham C-type shock models to the observed pure rotational and rovibrational H2 emission lines to determine the shock parameters that reproduce the high temperatures.
If this is right
- The southern star-forming regions are characterized by high gas density and moderate UV irradiation.
- Shock heating is consistent with near- and mid-IR tracers and diagnostics across the regions.
- Local starburst activity is the likely driver of the shocks and gas dynamics, though an AGN jet contribution is not ruled out.
Where Pith is reading between the lines
- If similar shock signatures appear in other GATOS sample galaxies, it may indicate a common role for circumnuclear star formation in heating molecular gas.
- Higher spatial resolution observations could distinguish between starburst-driven shocks and possible AGN jet contributions.
- These results suggest that shock models should be included when interpreting high-temperature H2 emission in obscured AGN hosts.
Load-bearing premise
That the Paris-Durham shock models and PDR calculations accurately represent the conditions and that other heating mechanisms contribute negligibly once line ratios are matched.
What would settle it
Detection of line ratios that deviate significantly from those predicted by the slow C-type shock models at 10 km/s, such as stronger high-velocity shock signatures or dominant fluorescent excitation patterns.
read the original abstract
We present combined JWST NIRSpec and MIRI/MRS integral field spectroscopy data of the nuclear and circumnuclear regions of the highly dust obscured Seyfert 2 galaxy NGC 7582, which is part of the sample of AGN in the Galaxy Activity, Torus and Outflow Survey (GATOS). Spatially resolved analysis of the pure rotational H$_2$ lines (S(1)-S(7)) reveals a characteristic power-law temperature distribution in different apertures, with the two prominent southern star-forming regions exhibiting unexpectedly high molecular gas temperatures, comparable to those in the AGN powered nuclear region. We investigate potential heating mechanisms including direct AGN photoionisation, UV fluorescent excitation from young star clusters, and shock excitation. We find that shock heating gives the most plausible explanation, consistent with multiple near- and mid-IR tracers and diagnostics. Using photoionisation models from the PhotoDissociation Region Toolbox, we quantify the ISM conditions in the different regions, determining that the southern star-forming regions have a high density ($n_H \sim 10^{5}$ cm$^{-3}$) and are irradiated by a moderate UV radiation field ($G_0 \sim 10^{3}$ Habing). Fitting a suite of Paris-Durham shock models to the rotational H$_2$ lines, as well as rovibrational 1-0 S(1), 1-0 S(2), and 2-1 S(1) H$_2$ emission lines, we find that a slow ($v_s \sim 10$ km/s) C-type shock is likely responsible for the elevated temperatures. Our analysis loosely favours local starburst activity as the driver of the shocks and circumnuclear gas dynamics in NGC 7582, though the possibility of an AGN jet contribution cannot be excluded.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents JWST NIRSpec and MIRI/MRS integral-field spectroscopy of the nuclear and circumnuclear regions of the dust-obscured Seyfert 2 galaxy NGC 7582. Spatially resolved analysis of pure rotational H2 lines (S(1)–S(7)) shows power-law temperature distributions, with unexpectedly high temperatures in the southern star-forming apertures comparable to the AGN nucleus. PDR Toolbox modeling yields n_H ~ 10^5 cm^{-3} and G0 ~ 10^3 in the southern regions. Fitting Paris-Durham C-type shock grids to both rotational and rovibrational H2 lines (including 1-0 S(1), 1-0 S(2), 2-1 S(1)) leads to the conclusion that a slow (v_s ~ 10 km/s) C-type shock provides the dominant heating, with local starburst activity loosely favored over an AGN jet contribution.
Significance. If the central claim holds, the work supplies a concrete example of how shock heating can produce molecular temperatures in circumnuclear star-forming regions that rival those in the AGN itself, using a combination of JWST spectroscopy, PDR Toolbox diagnostics, and external shock grids. The spatially resolved multi-aperture approach and the explicit inclusion of both rotational and rovibrational lines are strengths that could serve as a template for similar GATOS or other AGN-host studies.
major comments (2)
- [heating mechanisms and model-fitting discussion] The claim that shock heating is the most plausible mechanism (abstract and heating-mechanisms discussion) rests on the assumption that AGN photoionisation and UV fluorescence contribute negligibly once PDR Toolbox n_H and G0 are fixed. No quantitative side-by-side predictions of the S(1)–S(7) ladder or 1-0/2-1 ratios from a photoionisation grid at the same n_H ~ 10^5 cm^{-3} and G0 ~ 10^3 are presented, leaving the uniqueness of the Paris-Durham v_s ~ 10 km/s solution untested. This comparison is load-bearing for the central conclusion.
- [shock-model fitting section] The manuscript reports fits of Paris-Durham shock models to the observed H2 lines but supplies no quantitative goodness-of-fit metrics (e.g., reduced chi-squared, residual statistics, or formal error budgets on v_s). Without these, it is difficult to assess how well the slow C-type solution is preferred over other velocities or alternative heating scenarios.
minor comments (2)
- The phrase 'loosely favours local starburst activity' in the abstract could be replaced by a more precise statement of the quantitative evidence that distinguishes starburst-driven shocks from possible AGN-jet contributions.
- Ensure that all aperture positions and extraction regions are explicitly labeled in the figures showing the temperature distributions and line-ratio maps.
Simulated Author's Rebuttal
We are grateful to the referee for their thorough and insightful comments, which have helped us improve the clarity and robustness of our analysis. We address each major comment below and have revised the manuscript accordingly to include the requested comparisons and quantitative metrics.
read point-by-point responses
-
Referee: The claim that shock heating is the most plausible mechanism (abstract and heating-mechanisms discussion) rests on the assumption that AGN photoionisation and UV fluorescence contribute negligibly once PDR Toolbox n_H and G0 are fixed. No quantitative side-by-side predictions of the S(1)–S(7) ladder or 1-0/2-1 ratios from a photoionisation grid at the same n_H ~ 10^5 cm^{-3} and G0 ~ 10^3 are presented, leaving the uniqueness of the Paris-Durham v_s ~ 10 km/s solution untested. This comparison is load-bearing for the central conclusion.
Authors: We thank the referee for this valuable observation. The PDR Toolbox modeling was used to derive the physical conditions, but we recognize that a direct comparison of predicted line ratios from PDR models versus shock models at the same parameters would strengthen our argument for shock dominance. In the revised manuscript, we have added a new subsection in the discussion of heating mechanisms that presents the expected H2 S(1) to S(7) line ratios and 1-0/2-1 ratios from the PDR models at n_H ≈ 10^5 cm^{-3} and G0 ≈ 10^3. These predictions show significantly lower excitation in the higher rotational lines than observed, consistent with additional heating from shocks. We have also referenced literature on UV fluorescence effects to argue that the observed line ratios do not match fluorescent excitation signatures. This addition supports the uniqueness of the slow C-type shock solution. revision: yes
-
Referee: The manuscript reports fits of Paris-Durham shock models to the observed H2 lines but supplies no quantitative goodness-of-fit metrics (e.g., reduced chi-squared, residual statistics, or formal error budgets on v_s). Without these, it is difficult to assess how well the slow C-type solution is preferred over other velocities or alternative heating scenarios.
Authors: We agree that including quantitative goodness-of-fit statistics is essential for a rigorous assessment of the model fits. In the updated shock-model fitting section, we now report the reduced chi-squared values for the Paris-Durham C-type shock models across a range of velocities. The v_s ≈ 10 km/s model provides the best fit, and we include residual plots and formal uncertainties on v_s derived from the chi-squared contours. Higher velocity models yield poorer fits and systematically overpredict certain lines. These metrics confirm the preference for the slow shock solution over alternatives. revision: yes
Circularity Check
No circularity: external model grids fitted to observations
full rationale
The paper's central conclusion—that slow C-type shocks best explain the elevated H2 temperatures—arises from fitting independent Paris-Durham shock models and PDR Toolbox photoionisation calculations to the observed JWST NIRSpec/MIRI line ratios and power-law temperature distributions. These grids originate from prior external literature and are not derived from or defined by the current dataset. The analysis compares multiple tracers (rotational H2 S(1)–S(7), rovibrational lines) against the model predictions at fixed n_H ~10^5 cm^{-3} and G0 ~10^3, without any self-definitional reduction, fitted-input-as-prediction, or load-bearing self-citation chain. The derivation remains falsifiable against the external benchmarks and contains no ansatz smuggling or renaming of known results.
Axiom & Free-Parameter Ledger
free parameters (1)
- shock velocity =
~10 km/s
axioms (1)
- domain assumption The pure rotational H2 lines (S(1)-S(7)) trace a power-law temperature distribution whose parameters can be directly compared to shock-model predictions.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.