# A biconical ionised gas outflow and evidence of positive feedback in NGC 7172 uncovered by MIRI/JWST

L. Hermosa Muñoz, A. Alonso-Herrero, M. Pereira-Santaella, I.

García-Bernete, S. García-Burillo, B. García-Lorenzo, R. Davies, T. Shimizu,

D. Esparza-Arredondo, E. K. S. Hicks, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

L. Hermosa Muñoz, A. Alonso-Herrero, M. Pereira-Santaella, I. García-Bernete, S. García-Burillo, et al. A biconical ionised gas outflow and evidence of positive feedback in NGC 7172 uncovered by MIRI/JWST. Astronomy & Astrophysics - A&A, 2024, 690, 10.1051/0004-6361/202450262. insu-04822430

HAL Id: insu-04822430

https://insu.hal.science/insu-04822430v1

Submitted on 6 Dec 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
# A biconical ionised gas outflow and evidence of positive feedback in NGC 7172 uncovered by MIRI/JWST

L. Hermosa Muñoz¹,*, A. Alonso-Herrero¹, M. Pereira-Santaella², I. García-Bernete³, S. García-Burillo⁴,

B. García-Lorenzo⁵,⁶, R. Davies⁷, T. Shimizu⁷, D. Esparza-Arredondo⁵,⁶, E. K. S. Hicks⁸,⁹,¹⁰, H. Haidar¹¹,
M. Leist⁹, E. López-Rodríguez¹², C. Ramos Almeida⁵,⁶, D. Rosario¹¹, L. Zhang⁹, A. Audibert⁵,⁶,

E. Bellocchi¹³,¹⁴, P. Boorman¹⁵, A. J. Bunker³, F. Combes¹⁶, S. Campbell¹¹, T. Díaz-Santos¹⁷,¹⁸, L. Fuller⁹,
P. Gandhi¹⁹, O. González-Martín²⁰, S. Hönig¹⁹, M. Imanishi²¹,²², T. Izumi²¹,²², A. Labiano²³,
N. A. Levenson²⁴, C. Packham⁹, C. Ricci²⁵,²⁶, D. Rigopoulou³,¹⁸, D. Rouan²⁷, M. Stalevski²⁸,²⁹,

M. Villar-Martín³⁰, and M. J. Ward³¹

(Affiliations can be found after the references)

Received 5 April 2024 / Accepted 21 July 2024

## ABSTRACT

We present observations of the type-2 Seyfert NGC 7172 obtained with the medium-resolution spectrometer (MRS) of the Mid-Infrared Instrument (MIRI) on board the James Webb Space Telescope (JWST). This galaxy hosts one of the lowest ionised gas mass outflow rates ($M_{\text{out}} \sim 0.005 M_{\odot} \text{yr}^{-1}$) in a sample of six active galactic nuclei (AGN) with similar bolometric luminosities (log $L_{\text{bol}} \sim 44 \text{ erg s}^{-1}$) within the Galactic Activity, Torus, and Outflow Survey (GATOS). We aim to understand the properties of the ionised gas outflow and its impact on the host galaxy. We mainly used the ionised gas emission lines from the neon transitions, which cover a broad range of ionisation potentials (IPs) from ~20 eV to ~130 eV. We applied parametric and non-parametric methods to characterise the line emission and kinematics. The low excitation lines (IP < 25 eV, e.g. [Ne II]) trace the rotating disc emission. The high excitation lines (IP > 90 eV, e.g. [Ne V]), which are likely photoionised exclusively by the AGN, are expanding in the direction nearly perpendicular to the disc of the galaxy, with maximum projected velocities in the range of ~350–500 km s⁻¹. In particular, [Ne V] and [Ne VI] lines reveal a biconical ionised gas outflow emerging north-south from the nuclear region, extending at least ~2.5" N and 3.8" S (projected distance of ~450 and 680 pc, respectively). Most of the emission arising in the northern part of the cone was not previously detected due to obscuration. Given the almost face-on orientation of the outflow and the almost edge-on orientation of the galaxy, NGC 7172 may be a case of weak coupling. Nevertheless, we found evidence of positive feedback in two distinct outflowing clumps at projected distances of 3.1" and 4.3" (i.e. ~560 and 780 pc) south-west of the AGN. We estimated a star formation rate in these regions using the [Ne II] and [Ne III] luminosities of $0.08 M_{\odot} \text{yr}^{-1}$, which is ~10% of that found in the circumnuclear ring. The star formation activity might have been triggered by the interaction between the ionised gas outflow and the interstellar medium of the galaxy.

**Key words.** galaxies: active – galaxies: ISM – galaxies: kinematics and dynamics – galaxies: nuclei – galaxies: Seyfert

## 1. Introduction

The evolution of galaxies over cosmic time is believed to be impacted to a certain degree by feedback processes, particularly for those harbouring an active galactic nucleus (AGN) (see e.g. Kormendy & Ho 2013; Harrison et al. 2018; Harrison & Ramos Almeida 2024). These processes may play a crucial role in altering the gas content of galaxies by driving it outwards in so-called outflows, or by accreting new gas from the intergalactic medium in so-called inflows (see Davies et al. 2014; Storchi-Bergmann & Schnorr-Müller 2019). Several works in the literature have explored the important question of the impact of outflows on the ability of their host galaxies to form new stars (see e.g. Cresci et al. 2015a; Harrison 2017; Venturi et al. 2021; Veilleux et al. 2020; Speranza et al. 2024). On the one hand, if the gas is heated or expelled from the galaxy as outflows, these could potentially prevent star formation (SF) (i.e. negative feedback; e.g. Harrison 2017). On the other hand, if the gas cools down due to fragmentation of the outflow or compression of molecular clouds, SF could be triggered, even out of the galaxy plane (i.e. positive feedback; e.g. Silk

2013; Cresci et al. 2015b; Gallagher et al. 2019; Bellocchi et al. 2023). There are reported cases of positive feedback occurring at low (Maiolino et al. 2017; Shin et al. 2019; Perna et al. 2020; Venturi et al. 2023) and high redshifts (e.g. Cresci et al. 2015a; Nesvadba et al. 2020), although it is a relatively new observational phenomena. Thus, the complete picture of the impact of AGN-launched outflows into their host galaxies is still not clear, despite several attempts in the literature (Fiore et al. 2017; Davies et al. 2020).

NGC 7172 is an almost edge-on spiral galaxy with a type-2 Seyfert nucleus (Véron-Cetty & Véron 2006), log $L_{\text{bol}} \sim 44.1 \text{ erg s}^{-1}$ (Davies et al. 2020), which has been widely studied in the literature (see e.g. Sharples et al. 1984; Thean et al. 2000; Smajić et al. 2012; Thomas et al. 2017; Alonso Herrero et al. 2023). It is a nearby object located at a redshift of 0.00868 (NASA/IPAC Extragalactic Database; NED) and a redshift-independent distance of 37 Mpc (1" ~ 180 pc, see Davies et al. 2015). The galaxy hosts a circumnuclear ring with a diameter of approximately 1.1–1.4 kpc, observed with the CO(3–2) transition (345 GHz; Alonso Herrero et al. 2023, from now on AH23). A dust lane crossing the major axis covers emission from the central and northern regions, especially in the optical range. Thomas et al. (2017) detected an ionisation cone in this

<footnote>* Corresponding author; lhermosa@cab.inta-csic.es</footnote>
range with an opening angle of 120° that appears to be extend-
ing out of the galaxy and to be mostly detected to the south in the
optical, with the northern part nearly fully obscured by the host
galaxy. Davies et al. (2020) derived the nuclear properties of the
ionised gas outflow using X-shooter optical data and obtained a
relatively low ionised mass outflow rate, Ṁout ~ 0.005 M⊙ yr⁻<fnref n="1" />.
Alonso Herrero et al. (2023) also detected the ionised gas out-
flow in the near-infrared (near-IR) using the [Si VI] 1.96 µm line
and concluded that it is causing the molecular gas ring to outflow
in the plane of the galaxy (see also Stone et al. 2016). There have
however been no detections of an outflow in X-rays (Igo et al.
2020).

We observed NGC 7172 with the mid-infrared spectrometer (MRS) of the Mid-Infrared Instrument (MIRI; Wells et al. 2015; Rieke et al. 2015; Wright et al. 2015, 2023) on board the *James Webb* Space Telescope (JWST; Gardner et al. 2023). This dataset was observed within the Galactic Activity, Torus, and Outflow Survey (GATOS) collaboration (García-Burillo et al. 2021; Alonso-Herrero et al. 2021), with the aim of understanding the relation of the outflow properties in a sample of six hard X-ray selected nearby AGN with similar luminosities (log L<sub>AGN</sub> ~ 44 erg s<sup>-1</sup>, Davies et al. 2020; Zhang & Packham 2024). Among these galaxies, NGC 7172 shows one of the lowest ionised mass outflow rates (see Davies et al. 2020). In this paper, we present the analysis of the MIRI data on NGC 7172 (Sect. 2) showing the main results (i.e. kinematic, flux, and line ratios maps; Sect. 3). We also discuss the detection of an ionised gas outflow and the possible existence of positive AGN feedback (Sects. 4 and 5).

**2. Data and methodology**

The target was observed with the medium-resolution spec-
trometer (MRS) of MIRI/JWST within the program ID
1670 (P.I.: T. Shimizu and R. Davies). For the data reduc-
tion, we primarily followed the standard MRS pipeline
procedure (e.g. Labiano et al. 2016; Bushouse et al. 2023
and references therein) and the same configuration of the
pipeline stages described in García-Bernete et al. (2022) and
Pereira-Santaella et al. (2022) to reduce the data using the
pipeline release 1.11.4 and the calibration context 1130. Some
hot and cold pixels are not identified by the current pipeline
version, so we added an extra step before creating the data
cubes to mask them. The data reduction and extra steps are
described in detail in García-Bernete et al. (2024a); the nuclear
spectrum of NGC 7172 can be seen in their Fig. 1 (see also
Pereira-Santaella et al. 2022; García-Bernete et al. 2022). The
total wavelength range of the MRS goes from 4.9 µm to 28 µm,
divided into four channels with different resolutions and fields
of view (FoVs; see details in Pereira-Santaella et al. 2022;
García-Bernete et al. 2024a).

In this work, we individually modelled the main ionisa-
tion emission lines in the data cubes using the LMFIT package
(v1.2.2) inside a PYTHON environment (v3.9.12). We fitted each
emission line using a combination of a linear fit of the contin-
uum around each line (~0.03 µm for each side) with at least one
Gaussian component and a maximum of three per line. To select
the adequate number of Gaussians that are sufficient to obtain a
proper modelling to the lines, we used the ε criteria described
in Cazzoli et al. (2018) (see also Hermosa Muñoz et al. 2024).
In short, if the residuals of the continuum in the line after the
Gaussian modelling are larger than three times the residuals of
the continuum out of the line, then an extra component is added.
When using more than one Gaussian, they are ordered to ensure
continuity in the kinematic and flux maps, typically assuming the

primary component to be the narrowest. We imposed the signal-
to-noise (S/N) of the peak of the lines to be >3. After the mod-
elling, we corrected the velocity dispersion from its instrumental
value (Argyriou et al. 2023) using σ = √σ<sup>2</sup><sub>obs</sub> − σ<sup>2</sup><sub>inst</sub>. The maps
with the modelling with a single Gaussian for the neon emission
lines are shown in Fig. 1, and their corresponding continuum
maps are in Fig. B.1.

Besides the parametric modelling, we applied a non-
parametric modelling to the emission lines following a method
similar to Harrison et al. (2014). Contrary to this work, we did
not use previous Gaussian modelling to define the line, but rather
we defined the line as where the flux is higher than three times
the standard deviation of the continuum near the line (from
±700 to ±1000 km s⁻<fnref n="1" />, see example in Fig. A.2). Then we mea-
sured the parameters v02, v10, v50, v90, and v98, defined as the
velocity corresponding to the 2%, 10%, 50%, 90%, and 98% of
the total accumulated flux from the line, without assuming any
prior conditions to characterise the profiles (i.e. parameters from
Gaussian profiles). The v02 and v98 parameters are indicators of
the maximum and minimum velocities of the line, and thus are
typically used in the literature as measurements of the outflow
velocities (see e.g. Speranza et al. 2024). In general, the results
from the parametric and non-parametric modellings are consis-
tent (see Hervella Seoane et al. 2023, for a detailed comparison
between both methods). We used the parametric modelling for
the spaxel-by-spaxel analysis, and the non-parametric method
only for modelling the profiles from integrated spectra in some
regions of interest.

Given that the kinematic results for the lines are consistent
for different ionisation potentials (IPs; see Sect. 3.1), in this work
we focus mostly on a limited number of them tracing gas emis-
sion from low to high ionisation states. We selected the neon
emission lines, in particular [Ne II] 12.81 µm, [Ne III] 15.55 µm,
[Ne V] 14.32 µm (hereafter [Ne V] unless specified otherwise)<fnref n="1" />,
and [Ne VI] 7.65 µm, as they all are bright lines from the same
element with different IPs (see Table 1), and thus and we can
exploit their line ratios to further understand the physical ori-
gin of the gas (estimation of densities, etc.). The kinematic
maps of the mentioned neon lines obtained with a paramet-
ric modelling with a single Gaussian component are shown in
Fig. 1.

We also generated spatially resolved line ratios with the neon
emission lines flux maps (i.e. [Ne II], [Ne III], and [Ne V]) fitted
with one Gaussian (see Sect. 3.2). These lines are all located in
channel 3 (FoV 5.2" × 6.2"), and thus the proximity in wave-
length of these lines ensures that the ratios are not significantly
dependent on the extinction (Tommasin et al. 2008).

**3. Results**

The majority of the emission lines show signatures of two dif-
ferent kinematic components (see examples in Figs. 2 and A.1).
In the central part of the galaxy, the line profiles are symmetric
and can be modelled adequately with a single Gaussian. How-
ever, further out we detected wings in the lines both to the blue
and to the red (see left and right panels of Fig. 2). To illustrate
this, we selected two regions away from the diffraction ring of
the point spread function (PSF, see Fig. B.1) at both sides of the
AGN towards the blueshifted and redshifted parts of the velocity
field for all the emission lines. They are located at ~3" north-east

<footnote n="1">¹ We also detected the [Ne V]24.32 µm, but we will focus on the
[Ne V]14.32 µm as it has better spatial and spectral resolution.</footnote>
Gaussian profile<fnref n="2" />. We estimated the errors accounting for the standard deviation of the continuum near the line (width of ~100 km s⁻¹ per each side, from ±800 km s⁻¹) and the spectral resolution for each line (Argyriou et al. 2023). Additionally, for the emission lines in the integrated nuclear spectrum, we performed a Monte Carlo test with 30 000 iterations to further constrain the uncertainties from the W80 measurements coming from the continuum noise (see more details in Appendix A). The W80 measurements obtained with this technique are consistent with the reported values in Table A.1 in all cases. For the lines with the highest S/N, as [Ne II], [Ne III], or [Ne V] 14 µm, the statistical uncertainties are smaller than 10 km s⁻¹, and up to ~50 km s⁻¹ for lines with the lowest S/N, such as [Fe VIII] (S/N ~ 5, see upper middle panel in Fig. A.1).

The results, in Fig. 3, mainly show that all the emission lines with similar IPs behave similarly. There is a slight trend of increasing W80 at higher IPs (slope 0.4 ± 0.2), although the correlation is moderate (a Pearson coefficient of ρ ~ 0.5 and a *p*-value of ~0.08, see Table A.2 and Appendix A), and the error bars are large (average of ~50 km s⁻¹). If we were only to account for the statistical errors from the continuum with the Monte Carlo approach, we would derive a tighter and statistically significant trend between these two parameters (slope 0.7 ± 0.2; ρ ~ 0.8, and *p*-value < 0.001). There is a similar trend in the selected SW integrated region (not for the NE, see Fig. A.3). A similar correlation was detected with MIRI data for the galaxy NGC 7469, which was interpreted as a stratified outflow (Armus et al. 2023), and previously with Spitzer data in Dasyra et al. (2011), which was interpreted as stratification of the gas in the AGN narrow line region (NLR), with the high excitation lines typically closer to the AGN (Filippenko & Halpern 1984; Ho 2009).

We do not detect any statistically significant trend for both W80 and v98 parameters with the critical density, $n_{crit}$, of the lines (as in Dasyra et al. 2011), except for W80 in the nuclear and the SW regions where there is a moderate correlation (see Appendix A and Table A.2). This explains the differences in W80 measured between [Ne V] at 14.3 and 24.3 µm in Fig. 3 and Fig. A.3 (see also Table A.1), although in the nuclear region they are still consistent within the errors. The critical density is larger for the [Ne V] at 14.3 µm than at 24.3 µm (log $n_{crit}$ = 4.51 cm⁻³ and 3.77 cm⁻³, respectively) (see Pereira-Santaella et al. 2022), and the lines with high critical densities are produced more efficiently close to the AGN; that is also where the motions are faster, and thus the kinematics could be more perturbed. These trends have been reported in other Seyferts using optical emission lines (see e.g. De Robertis & Osterbrock 1986).

We found that the largest values of the outflow velocity (v98 ~ 400-500 km s⁻¹, maximum in absolute value between v02 and v98) are observed in the majority of lines in the NE region (see Table A.1), indicating that this part is moving faster than the gas in the disc or the southern regions. This was expected from the profiles in the left panels of Fig. 2, and it is also seen in the velocity channel maps (see Figs. B.4, B.5, and B.6), and in the parametric modelling (see Sects. 3.3.1 and 4.1). We note that the average velocity of the lines to the NE (SW) are generally blueshifted (redshifted) with respect to the nuclear spectrum (see Table A.1 and Fig. 2), although some nuclear lines have slightly blueshifted velocities as well.

<table><caption>**Table 1.** Flux measurements of the main neon emission lines (R = 0.5").</caption><thead><tr><th>Line</th><th>IP (eV)</th><th>f<sub>AGN</sub></th><th>f<sub>NE</sub> (×10<sup>-16</sup> erg s<sup>-1</sup> cm<sup>-2</sup>)</th><th>f<sub>SW</sub></th><th>f<sub>C1</sub></th></tr></thead><tbody><tr><td>[Ne II]</td><td>21.6</td><td>238.07 ± 0.65</td><td>8.86 ± 0.30</td><td>12.58 ± 0.34</td><td>8.40 ± 0.05</td></tr><tr><td>[Ne III]</td><td>41.0</td><td>377.35 ± 0.64</td><td>7.07 ± 0.34</td><td>12.82 ± 0.14</td><td>5.68 ± 0.09</td></tr><tr><td>[Ne V]</td><td>97.2</td><td>294.25 ± 0.58</td><td>4.19 ± 0.28</td><td>4.20 ± 0.11</td><td>1.94 ± 0.07</td></tr><tr><td>[Ne VI]</td><td>126.0</td><td>40.12 ± 1.76</td><td>4.23 ± 0.30</td><td>2.47 ± 0.19</td><td>—</td></tr></tbody></table>

**Notes.** The regions of interest are described in Sect. 3 and shown in Fig. 1. All the fluxes have been estimated using a single Gaussian modelling. We have not applied an aperture correction.

## 3.2. Morphology of the ionised gas

The low-ionisation emission lines (IP ≤ 25 eV) in this mid-infrared (mid-IR) range, such as [Ne II], [Ar II] 6.98 µm, [S III] 18.7 µm, and [Fe II] 5.34 µm (this line is discussed in García-Bernete & Rigopoulou 2024b), are typically produced by SF processes (Pereira-Santaella et al. 2010), or shocks (Allen et al. 2008). The flux maps for these lines, in particular that of [Ne II] (see top left panel in Fig. 1), show a disc-like emission oriented east-west, coincident with the direction of the disc in the galaxy (PA = 88°, AH23) and with the SF ring (1.1–1.4 kpc in diameter, i.e. ~6–7", AH23). The spatial resolution with ALMA in AH23 (beam ~ 0.07") is better than that of MIRI/MRS for the [Ne II] (~0.25"; Argyriou et al. 2023), and thus the circumnuclear disc structure is not as well resolved as with CO(3–2).

Higher-ionisation emission lines (IP > 90 eV) are most likely produced by AGN photoionisation (Armus et al. 2007; Pereira-Santaella et al. 2010). This is reflected in the flux maps of the [Ne V] and [Ne VI] lines, where no emission is detected in the disc. Rather, it is concentrated in the nucleus and in the north-south direction with an apparent two-sided X-shaped morphology, which is especially prominent in [Ne VI] (see two bottom rows in Fig. 1). We note that the [Ne V] intensity map in particular could have some contamination from the PSF (see Sect. 3.3.2 and Fig. B.1). The bicone is extended to at least ~2.5" N and ~3.8" S (i.e. 450 pc and 680 pc, respectively, in projection) with an opening angle of ~128°, similar to that found through the optical [O III] λ5007 Å emission (Thomas et al. 2017). The northern part of the cone was not detected in the near-IR line [Si VI] (AH23), likely due to the dust lane crossing the galaxy (see Sect. 1), and was only detected at larger distances in optical wavelengths (Thomas et al. 2017). The emission lines with intermediate IPs (25–90 eV), such as [Ne III], [Ar III] 8.99 µm, and [S IV] 10.51 µm (see discussion on this line in García-Bernete & Rigopoulou 2024b), show intermediate properties and morphologies between the low and high excitation lines. This was expected as they are ionised by a mixture of SF processes and AGN photoionisation (Pereira-Santaella et al. 2010; García-Bernete et al. 2017). As is seen in the second panel from the top of Fig. 1, the flux map of [Ne III] is more concentrated in the nuclear region, but still elongated towards both the disc and the bicone direction, with complex kinematics (see Sect. 3.3.1 and Fig. B.5).

The [Ne III]/[Ne II] ratio map (Fig. 4) shows a similar biconical structure extending in the north-south direction with increased line ratios with respect to the disc direction (average in log of ~0.1 vs ~−0.7, respectively; see left panel in Fig. 4). These values are consistent with those of type-1 Seyferts and H II regions, respectively (see Table 3 in Pereira-Santaella et al. 2010 and Fig. 3 in García-Bernete & Rigopoulou 2024b). There

<footnote n="2">² The W80 parameter is not corrected by the instrumental resolution.</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Emission line profiles of the ionised gas obtained from an integrated spectrum (circular aperture with $R \sim 0.5''$) in the nuclear region (middle), NE (left), and SW (right) parts of the cone (see Fig. 1). All the velocities refer to the systemic value computed with $z = 0.00868$ (see Sect. 1) indicated by the dashed, vertical grey lines. Fluxes were normalised to the maximum of each line.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** W80 (top panel) and v98 (bottom panel) measured in the integrated nuclear spectrum ($R \sim 0.5''$) with respect to the emission lines IP. The dashed grey line in the top panel corresponds to a linear fit of the points and the $1\sigma$ uncertainty. The grey squares in the bottom panel indicate the v98 value reported by Davies et al. (2020) for the optical [S II] and [O III] lines with their uncertainties. The error bars in our data are obtained accounting for the standard deviation of the continuum near the line and the spectral resolution.</figcaption></figure>

is a distinct region located ~3" SW from the nucleus with $\log([\text{Ne III}]/[\text{Ne II}]) \sim -0.2$ and $\log([\text{Ne V}]/[\text{Ne II}]) \sim -0.7$, which corresponds to one of the two kinematically distinct

clumps discussed in Sect. 3.3.1 (see also Sect. 4.2). The [Ne V]/[Ne III] ratio map (right panel in Fig. 4) has several enhanced regions ($\log([\text{Ne V}]/[\text{Ne III}]) > -0.2$) mainly in the nuclear region and several knots in the ionisation cone. These values are consistent with the ratios derived for Seyfert galaxies using Spitzer/IRS data (Pereira-Santaella et al. 2010).

### 3.3. Spatially resolved kinematics of the ionised gas

In this section, we present the results for the kinematics of the ionised gas emission lines according to their IPs (i.e. below and above 90 eV).

#### 3.3.1. Low excitation lines

The mean velocity field of [Ne II] (top middle panel in Fig. 1) shows a rotating pattern with maximum velocities ~330 km s⁻¹ in absolute value ($\Delta v = v_{\text{max}} - v_{\text{min}} \sim 530$ km s⁻¹), oriented along the major axis of the disc (PA ~88° ± 1°, measured from north to east using the minimum and maximum velocities), with the receding side to the west. The velocity dispersion map (top right panel in Fig. 1) is not centrally peaked. Rather, the highest values are detected at ~2" from the nucleus to both east and the west ($\sigma \sim 170$ km s⁻¹), as was already seen in the cold molecular gas (AH23). When fitting the [Ne II] line with two Gaussians (see Fig. B.2), the secondary component is mainly located in the direction of the disc, with a rotating pattern oriented as the primary component. This emission detected as a secondary component along the disc could be attributed to either the SF ring in the central parts, the presence of a nuclear bar, or both (see discussion in AH23). Additionally, the [Ne II] (top row in Fig. 1) and also the [S III] intensity maps (see Fig. B.2) show the presence of two distinct clumps, C1 and C2, at a projected distance of ~3.1" and ~4.3" (i.e. ~560 and 780 pc) with enhanced flux, a low velocity dispersion ($\sigma \sim 50$ km s⁻¹), and blueshifted velocities. These are not detected in other low excitation lines because of their smaller FoVs (ch1 3.2" × 3.7" vs. ch3 5.2" × 6.2" and ch4 6.6" × 7.7"). We estimated that the clumps have a size of R~0.5" (i.e. total size of ~180 pc), so they are resolved in the MRS data (PSF FWHM in ch3 varying from 0.56" to 0.65", and 0.95" for ch4-short, see Fig. 1). We focus on the results from the main clump at ~540 pc from the centre (C1) seen with [Ne II], as the C2 is not fully seen as it is located at the limits of the channel 3 FoV, and the [Ne II] line is detected with better spatial and spectral resolution than [S III]. We further discuss these clumps in Sect. 4.2.
<figure><img src="image_5.png" /></figure>

**Fig. 4.** Maps of the [Ne III]/[Ne II], [Ne V]/[Ne II], and [Ne V]/[Ne III] line ratios (from left to right, respectively) represented in a logarithm scale. All the ratios are estimated after modelling the lines with a single Gaussian component (see Fig. 1). The white star indicates the photometric centre of the flux maps.

The ionised gas traced by intermediate IP lines, such as [Ne III] (see second row of Fig. 1), is rotating along a PA (~80° ± 2°) and with a velocity amplitude similar to [Ne II] (Δv ~ 490 km s⁻¹). The velocity dispersion is irregularly distributed and non-centrally peaked, with several regions with high σ values (average σ ~ 105 km s⁻¹). The central, circular-like structure of the velocity dispersion indicates the presence of some contamination from the strong PSF emission (see ch3-long continuum map in Fig. B.1 and Argyriou et al. 2023). From the channel maps (see Fig. B.5), it can be seen that the velocity is complex, showing emission both in the disc and in the ionisation cone direction. The [Ne III] gas is ionised by a mixture of, on the one hand, SF processes, which explains the similarities of the primary component to that of lines such as [Ne II] (see Fig. 1), and, on the other hand, AGN photoionisation, which explains the detection of emission oriented south from the nucleus out of the disc, potentially related to the high-ionisation line emission (see e.g. Pereira-Santaella et al. 2010).

### 3.3.2. High excitation lines

The kinematic maps for the high-ionisation emission lines are complex. The mean velocity field of the [Ne V] line (third row in Fig. 1) fitted with one Gaussian shows mainly blueshifted velocities, whereas the velocity dispersion map is not homogeneous (σ from 50 to 170 km s⁻¹, with a standard deviation of 100 km s⁻¹; see right panel, third row in Fig. 1), with larger values north from the nucleus. The redshifted velocities (<50 km s⁻¹) are mainly located in two regions NW and SW from the nucleus, whereas the most blueshifted velocities are found ~3'' (i.e. ~540 pc) NE (see Fig. 1). The intensity of the lines (see Figs. 1 and B.6) seems to indicate that the east side of the cone (blueshifted) is more illuminated than the west side (redshifted). Moreover, the NE part of the cone is approaching at larger mean-velocities than the south-east part (~350 vs. ~250 km s⁻¹), which suggests that the gas is not expanding completely parallel to the plane of the sky. We further discuss the geometry of the system in Sect. 4.1.

For [Ne VI] (bottom row in Fig. 1), the emission is concentrated in an X-like shape oriented north-south from the centre (opening angle ~90°), with the approaching (receding) side to the east (west) at a maximum mean velocity of ~200 km s⁻¹ (~80 km s⁻¹). The velocity dispersion of the primary component is rather homogeneous (average σ ~ 107 km s⁻¹), with higher values than the secondary component in the same region (σ ~ 130 km s⁻¹ vs. σ ~ 70 km s⁻¹, respectively).

The inner regions of the velocity and σ maps for [Ne V] and [Ne VI] (see Fig. 1) show the diffraction ring of the JWST PSF (see Argyriou et al. 2023, and continuum maps in Fig. B.1), which indicates that there is also an unresolved component in the line emission. Other galaxies in the GATOS sample, such as MCG-05-23-16 (see Esparza-Arredondo et al., in prep.), show unresolved intensity maps for [Ne V] and [Ne VI] lines, which are dominated by the PSF emission. For NGC 7172, the petals from the PSF (see ch2-long and ch3-medium maps in Fig. B.1) are not clearly seen in the intensity maps (see Fig. 1); thus, the extended emission is likely not much affected by it.

We have also modelled these lines with multiple Gaussian components (see Sect. 2), finding that the two-component modelling was the most preferred fit in several regions with similar properties for both lines. The second component is concentrated mainly to the south and NE of the AGN (see the modelling of [Ne V] in Fig. B.3). There are blue and red wings everywhere along the bicone (see channel maps in Fig. B.6), which would indicate that it is optically thin in these lines. Thus, the detection of various velocity components could be due to both the inclination of the cone with respect to the plane of the sky and its wide angle (see also Roy et al. 2021). We describe the full properties of the two-component modelling in Appendix B.

## 4. Discussion

### 4.1. The ionisation cone and the outflow

The MIRI/MRS observations of NGC 7172 have revealed for the first time a prominent biconical ionisation cone and associated outflow, which are particularly prominent for the [Ne V] and [Ne VI] lines (see Figs. 1 and 4). The emission is observed in the north-south direction, reaching radial distances from the AGN of ~2.5''N and ~3.8''S (i.e. projected distances of 450 pc and 680 pc, respectively). The ionisation cone extends further out than MIRI's FoV as seen from the [O III] emission in the southern part (Thomas et al. 2017). Despite the large extent of the outflow, NGC 7172 has one of the lowest ionised gas mass outflow rates ($M_{\text{out}} \sim 0.005 M_{\odot} \text{ yr}^{-1}$) among the galaxies studied by Davies et al. (2020). However, several effects might have affected the optical estimations of this parameter.

The galaxy is crossed by a prominent dust lane (see Fig. 1 in AH23), which is obscuring the northern part of the ionisation cone as well as the nuclear emission even in the near-IR. Estimates of the optical extinction indicate that $A_V > 3$ mag (Davies et al. 2020). In particular, based on the non-detection of
the northern part in the [Si VI]1.96 µm line in AH23, that region must be optically thick at that wavelength, which implies that the extinction at 2 µm should be $A_{2\mu m} > 1$ mag (or $A_V > 9$ mag Wright et al. 2023). This is corroborated by the extinction values of approximately $A_V = 15-20$ mag derived by Smajić et al. (2012) using near-IR hydrogen recombination lines. Moreover, both sides of the cone are similarly bright in [Ne V] and [Ne VI] (see Fig. 1), which indicates that the extinction effects are not dominant for these lines in the mid-IR. This obscuration could have affected the measurements of the $\dot{M}_{\text{out}}$ for the ionised gas (see Sect. 1 and Davies et al. 2020) by a minimum factor of five, plus at least a factor of two for the previous non-detection of the northern cone.

The derived velocities for the outflow are modest ($v_{98}$ up to ~500 km s⁻¹), but consistent with previous determinations in the optical (see Davies et al. 2020, and Fig. 3). These measurements, however, can be affected by the inclination of the cone with respect to our line of sight (see Sect. 3.3.2), so to correct them a cone-outflow geometrical model would be needed. If the outflow axis were perfectly aligned with the plane of the sky, we would expect a symmetrical velocity structure. Instead, the outflow might be inclined such that the part tilted away from us is still within the line of sight, resulting in an overall blueshift in the average velocity maps. The wide aperture angle (~120°) could also contribute to the observed velocity structure. A significant portion of the outflow could have a velocity component towards us, even if the outflow is not significantly inclined. Additionally, the lack of redshifted emission could imply that the rear part of the bicone is obscured by the dust in the galaxy plane, and we are only seeing the front part due to its orientation. The ALMA detected torus of NGC 7172 was proposed to be inclined at least ~67° (AH23), so if the gas is launched from the inner regions in the perpendicular direction to the torus, it would explain the apparent asymmetry in the velocity maps. Assuming this geometry, the $\dot{M}_{\text{out}}$ would be affected by a factor of ~2.4 due to the inclination of the outflow (see Eq. (4) in García-Burillo et al. 2014).

The electron density, $n_e$, also affects the determination of $\dot{M}_{\text{out}}$. We estimated it from the ratios of the two [Ne V] lines using PYNEB (Luridiana et al. 2015). We measured the line fluxes in the nuclear integrated spectrum ($R \sim 1''$), obtaining a ratio of $f_{[Ne V]14}/f_{[Ne V]24} = 1.52$, which is twice the value that was obtained by Spinoglio et al. (2015) with Spitzer data (integrating 4.7'' for [Ne V]₁₄ and 11.1'' for [Ne V]₂₄), but consistent with typical values for Seyfert galaxies (Pereira-Santaella et al. 2010). If we measure the fluxes at larger (smaller) radii, the ratio decreases (increases). Despite the fact that the mid-IR range is less affected by extinction, this galaxy has prominent silicate features in absorption that indicate that the emission could still be obscured (see García-Bernete et al. 2024a; García-Bernete & Rigopoulou 2024b). Following the extinction laws by Chiar & Tielens (2006) and Pei (1992), both [Ne V] lines fall in a minimum of the mid-IR extinction curve with similar extinction levels (see also Fig. 4 in Hernán-Caballero et al. 2020). This means that the fluxes of these lines are equally, and not heavily, impacted by extinction. Therefore, the line ratio does not need to be corrected from extinction. We assumed a temperature of $10^4$ K for the gas, which gives $n_e \sim 3200$ cm⁻³, three times smaller than that estimated by Davies et al. (2020). If we consider a $R = 2''$ instead, the ratio between the fluxes is 1.12 and the corresponding density would be ~1150 cm⁻³, lower but in the same order of magnitude. We note that if we instead consider a range of temperatures from $T \sim 10^3$ to $10^5$ K, the order of magnitude of $n_e$ remains between $10^3$ and $10^4$ cm⁻³ for the

measured ratio, which is one order of magnitude larger than the estimates obtained using the [S II] emission lines, but consistent with the derivations using the ionisation parameter (see Davies et al. 2020). Thus, considering a $R$ of 1'' to directly compare with the optical estimations, this would add a factor of three in the estimation of the ionised mass outflow rate.

Taking into account all the different factors, the mass outflow rate for the ionised gas in NGC 7172 might have been underestimated in Davies et al. (2020) using optical observations. Indeed, with the MIRI/MRS data Zhang & Packham (2024) followed Riffel et al. (2023) to estimate the mass outflow rate using the [Ne V] 14 µm, and derived a $\dot{M}_{\text{out}} \sim 0.03$ M$_\odot$ yr⁻¹ (not corrected by projection effects). This estimation assumed a spherical geometry for the outflow integrating with an aperture of 0.9'', to directly compare to the aperture used to derive the optical value in Davies et al. (2020). Thus, the integrated value is six times larger than the previous measurement. Despite this, NGC 7172 still remains one of the galaxies from the GATOS sample with the lowest ionised mass outflow rate measurements (see detailed estimations of $\dot{M}_{\text{out}}$ with the radius in Zhang & Packham 2024). We emphasise that this galaxy could be an extreme case among the GATOS sample because of the nearly edge-on orientation of the galaxy and the almost face-on orientation of the ionisation cone. It is likely that NGC 7172 is a case of weak or moderate geometrical coupling (see Fig. 4 in Ramos Almeida et al. 2022), where the expansion of the ionised gas outflow also triggers the expansion of the molecular gas producing another outflow in the direction of the disc, which was detected in CO with ALMA data in AH23. Thus, in general, the ionised mass outflow rate of NGC 7172 will be low compared to other galaxies with similar luminosities (see Zhang & Packham 2024).

## 4.2. Evidence of positive feedback

We detected two kinematically distinct clumps exclusively in the [Ne II] (see C1 and C2 red circles in Figs. 1 and B.4) and [S III] maps (see Fig. B.2), located at ~3.1'' and ~4.3'' SW from the AGN, which are also detected in the MIRI narrow-band images (see Fig. 6, Rosario et al. in prep.). They were not seen in AH23 with ALMA data due to the smaller FoV in the north-south direction. The low velocity dispersion, distinct morphology in the [Ne II] flux map (see Sect. 3.2), and the observed line ratios (log([Ne III]/[Ne II]) ~ -0.3; see Fig. 4) suggest that they are SF regions. Moreover, the velocities are blueshifted with respect to the rotating primary component of the gas in the same spatial region (i.e. the primary component; see Sect. 3.3.1 and Fig. B.4). Given that the approaching side of the galaxy is the eastern part, these regions cannot be located in the plane of the galaxy but are likely close to the edge of the ionised gas cone or outflow. Emission from [Ne V] is also detected in the same region (see Fig. B.6), but the flux map shows no distinct features at the clump position as was observed for [Ne II]. It is likely that the interaction between the expanding outflow and the interstellar medium (ISM) in the disc has entrained gas and triggered an SF episode. Using the [Ne II], [Ne III], and [Ne V] lines, we estimated an SF rate (SFR) for C1 of ~0.08 M$_\odot$ yr⁻¹ using the relations by Zhuang et al. (2019). As a comparison, we estimated for an individual region of the circumnuclear ring (1.5'' west from the nucleus, see Fig. 2 in AH23 and upper-left panel in Fig. 1, integrating the fluxes with an aperture of $R \sim 0.5''$) an SFR of ~1 M$_\odot$ yr⁻¹, which is ten times larger than the SFR measured in the clumps. For comparison, in NGC 92 the SFRs for those H II regions detected in the tidal tails range from 0.003 to 0.01 M$_\odot$ yr⁻¹ (Torres-Flores et al. 2009). Cresci et al. (2015a)
report a case of positive feedback with the detection of two SF clumps, for which they derived an SFR of ~0.03 M$_{\odot}$ yr$^{-1}$, which is the same order of magnitude as the clumps in NGC 7172.

By qualitatively comparing the line ratios of these regions with the mid-IR diagnostic diagrams developed by Feltre et al. (2023), the integrated value of the clumps (log([Ne III]/[Ne II]) ~ -0.11 ± 0.02 and log([Ne V]/[Ne II]) ~ -0.62 ± 0.01) is consistent with photoionisation by a mixture of AGN, SF processes, and shocks<fnref n="3" /> (see Fig. 5). As a comparison, the integrated value for the nuclear region lies in the region exclusively excited by the AGN emission, whereas the selected regions of the cone (see Sect. 3) lie in the AGN and shock regions. Given that these regions are located on the edges of the outflow (see Figs. 1, B.6, and B.7), the ratios could be explained as gas photoionised by the AGN expanding from the centre and interacting with the ISM. Consistently with this scenario, the shocks in the clumps could have been produced due to the compression of the ISM gas by the outflow, eventually triggering the SF process. This suggests that positive feedback is taking place in NGC 7172, similarly to what was found by Cresci et al. (2015b) for the nearby Seyfert NGC 5643. However, unlike this galaxy, the radio emission in NGC 7172 is confined to the inner 2" (Thean et al. 2000) and is not clearly connected with the SF clumps. In AH23, the power of the radio emission was estimated and proposed to be capable of producing the molecular gas outflow seen in the disc, so we cannot discard a contribution from that emission. Bessiere & Ramos Almeida (2022) reported a similar case of enhanced recent SF at the edge of the approaching side of the co-planar ionised outflow of Mrk 34.

Fig. 6 shows a schematic representation of the outflow expansion direction (see Sect. 4.1) along with other relevant features in NGC 7172, such as the SF circumnuclear ring (AH23) and the clumps superimposed to a MIRI image in the F1000W filter (after the nuclear PSF subtraction), which is part of another GATOS Cycle 1 program (PID: 2064, PI: D. Rosario; Rosario et al., in prep.). We see in the same direction as clumps C1 and C2 that there is another clump ~5" from the AGN. This is located outside of the FoV of our MIRI/MRS observations. Other clumps are also visible to the NW in the MIRI image at larger projected distances (~7") from the AGN. Other galaxies in the literature are known to have extra-planar H II regions (see e.g. De Robertis & Osterbrock 1986; Mendes de Oliveira et al. 2004; Stein et al. 2017) at even larger distances from the plane of the disc. Although there are various ways of producing these regions, such as ejections from the disc, the majority of them are typically explained through past or ongoing interactions of the galaxies with other companions (see e.g. Mendes de Oliveira et al. 2004; Torres-Flores et al. 2009; Miralles-Caballero et al. 2012), so that the interaction triggers off-plane H II regions. NGC 7172 is located in a group of galaxies (Sharples et al. 1984) and the faint tidal tails seen in the large-scale structure (see Fig. 1 in AH23 and also Sharples et al. 1984) are a clear signature of a past merger event. Thus, we cannot rule out that some of the bright clumps detected in the MIRI image (see Fig. 6) could have been formed during those processes. However, given that several clumps (i.e. C1, C2, and the third clump at 5") are close (at least in projection) to the west edge of the outflow, they could be the product of a positive feedback episode.

<figure><img src="image_6.png" /><figcaption>Fig. 5. Mid-infrared line ratios following Fig. 5 in Feltre et al. (2023). We represent the spatially resolved line ratios for NGC 7172 as white and grey contours (levels from 1 to 10 spaxels). The yellow star represents the ratio estimated in the integrated spectrum for the nuclear region (AGN). The orange and green circles are the NE and SW regions in the cone. The blue circle represents the star-forming clump C1. The green region (i.e. square) represents AGN+shock models with velocities varying from 200 to 1000 km s⁻¹, the red region represents AGN models at ≥1/3 solar metallicity, and the purple hexagons and line represents AGN+SF models with the ionisation parameter increasing from left to right (for more details, see Fig. 5 of Feltre et al. 2023).</figcaption></figure>

## 5. Summary and conclusions

We have presented JWST MIRI/MRS observations of the GATOS galaxy NGC 7172 to study the ionised gas kinematics, with a particular focus on the neon emission lines (i.e. [Ne II], [Ne III], [Ne V], and [Ne VI]), which present a broad range of ionisation stages of the gas (i.e. IP from 20 to 130 eV). We have derived different kinematic properties of the lines through their modelling with parametric (i.e. one and two Gaussians per line) and non-parametric methods. Our results are as follows:

- Low-intermediate ionisation emission lines: The low-ionisation lines (e.g. [Ne II]) mainly trace the gas in the disc of the galaxy, with a secondary kinematic component located along the major axis probably related to the cold molecular gas ring detected with ALMA (AH23). Intermediate ionisation lines (e.g. [Ne III]) have a complex distribution both along the major and minor axis of the galaxy. They trace both gas emission and motions in the disc and part of the ionisation cone, as they are typically ionised by both the AGN emission and SF processes.

- High-ionisation emission lines: These lines (e.g. [Ne V] and [Ne VI]) are extended in the north-south direction from the centre out to the MIRI/MRS FoV edges (see Fig. 1 and 4), and reveal a double-sided ionisation cone that is oriented almost perpendicular to the galactic disc. Previous works had already detected an ionised gas outflow mainly with optical and near-IR data (Thomas et al. 2017; Davies et al. 2020; Alonso Herrero et al. 2023), but they were missing the northern cone, probably due to the obscuration produced by the dust lane along the major axis of NGC 7172.

- Non-parametric kinematical modelling: We found a mild correlation between the W80 parameter and the IP in the nuclear region (ρ ~ 0.5 and p-value ~ 0.08) and the SW region (projected distance of ~360 pc) in the ionisation cone (ρ ~ 0.65 and p-value ~ 0.03). These correlations are

<footnote n="3">³ Note that these values are estimated using a single Gaussian modelling (see Table 1), which means that both the disc and the clump components contribute to these line ratios.</footnote>
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Schematic figure of the main structures detected in NGC 7172 overplotted on the PSF-subtracted MIRI F1000W image from Rosario et al. (in prep.). We indicate the approximate size of the cold molecular gas ring detected with ALMA in red (see AH23), the AGN position with a white star, and the ionised gas outflow direction with yellow arrows (see Sect. 4.1). Red circles indicate the position of several clumps detected in the MIRI image in the direction of the outflow, two of which coincide with the SF clumps seen with the MRS data (see top left panel of Fig. 1, and Sects. 3.3.1 and 4.2). The bright emission crossing the central region from the SE to the NW across the image is a residual diffraction spike not fully eliminated when subtracting the bright PSF (see more details in Rosario et al., in prep.).</figcaption></figure>

typically interpreted as a stratification of the gas from the NLR, such that the high excitation lines are produced closer to the central parts. However, in this case both the [Ne V] and [Ne VI] emission lines are extended up to the limits of MIRI/MRS’s FoV (~3”, i.e. 540 pc in projection).

- Properties of the outflow: The bicone extends out to ~2.5" N (i.e. 450 pc) and ~3.8" S (i.e. 680 pc) with an opening angle of ~120°, although it could potentially reach larger distances beyond MIRI/MRS’s FoV based on its optical morphology (Thomas et al. 2017). The morphology, intensity, and modest outflow velocities (v98 ~ 400 km s⁻¹) of the high-ionisation lines suggest that the expansion of the cone occurring is not completely co-planar with the plane of the sky. With the MIRI data, we are observing both the receding and the approaching sides of the cone. However, its orientation close to the plane of the sky means that the outflow velocities are only lower limits for the real velocities. We estimated that previous reported measurements of $\dot{M}_{\text{out}} \sim 0.005 M_{\odot} \text{yr}^{-1}$ underestimated the ionised mass outflow rate, mainly due to extinction and projection effects. Nevertheless, NGC 7172 is likely a case of weak or moderate geometrical coupling, so it remains as one of the objects with lowest $\dot{M}_{\text{out}}$ in the GATOS sample (see also Zhang & Packham 2024).

- Feedback: We detected two kinematically distinct clumps in the [Ne II] and [S III] maps at projected distances of ~3.1" and ~4.3" (i.e. 560 and 780 pc) SW from the nucleus (see Figs. 1, B.2, and B.3), blueshifted with respect to the galaxy rotation, and thus likely taking part in the outflow. The observed [Ne III]/[Ne II] ratio, however, is typical of SF and shock regions (SFR ~ 0.08 $M_{\odot} \text{yr}^{-1}$). Given their properties

and location close to the SW edge of the cone (at least in projection), we propose that these regions may represent a case of positive feedback produced by the interaction of the expanding ionised gas outflow with the ISM of the galaxy.

This work highlights the relevance of JWST data for studies on this kind, with the capability of detecting, characterising, and quantifying the properties of outflows in nearby AGN and studying their effects on their host galaxies (see e.g. Álvarez-Márquez et al. 2023; Armus et al. 2023; Zhang & Packham 2024; García-Bernete & Rigopoulou 2024b).

**Acknowledgements.** We thank the referee for his/her comments that have helped to improve the manuscript. LHM, AAH and MVM acknowledge financial support by the grant PID2021-124665NB-I00 funded by the Spanish Ministry of Science and Innovation and the State Agency of Research MCIN/AEI/10.13039/501100011033 PID2021-124665NB-I00 and by ERDF A way of making Europe. MPS acknowledges support from grants RYC2021-033094-I and CNS2023-145506 funded by MCIN/AEI/10.13039/501100011033 and the European Union NextGenerationEU/PRTR. IGB and DR acknowledge support from STFC through grants ST/S000488/1 and ST/W000903/1. SGB acknowledges support from the Spanish grant PID2022-1358560NB-I00, funded by MCIN/AEI/10.13039/501100011033/FEDER, EU. BGL and DEA acknowledge support from the Spanish Ministry of Science and Innovation through the Spanish State Research Agency (AEI-MCINN/10.13039/011000011033) through grants "Participation of the Instituto de Astrofísica de Canarias in the development of HARMONI: D1 and Delta-D1 phases" with references PID2019-107010GB-100 and PID2022-140483NB-C21 and the Severo Ochoa Program 2020-2023 (CEX2019-000920-S). DEA also acknowledge financial support from MICINN through the programme Juan de la Cierva. DEA, AA, and CRA acknowledge funding from the State Research Agency (AEI-MCINN) of the Spanish Ministry of Science and Innovation under the grant "Tracking active galactic nuclei feedback from parsec to kiloparsec scales",
with reference PID2022-141105NB-I00. MTL, CP, EKSH, and LZ acknowl-
edge financial support from program JWST-GO-01670 provided by NASA
through a grant from the Space Telescope Science Institute, which is oper-
ated by the Association of Universities for Research in Astronomy, Inc., under
NASA contract NAS 5-03127. EB acknowledges the María Zambrano pro-
gram of the Spanish Ministerio de Universidades funded by the Next Genera-
tion European Union and is also partly supported by grant RTI2018-096188-
B-I00 funded by the Spanish Ministry of Science and Innovation/State Agency
of Research MCIN/AEI/10.13039/501100011033. AJB acknowledges funding
from the “FirstGalaxies” Advanced Grant from the European Research Council
(ERC) under the European Union’s Horizon 2020 research and innovation pro-
gram (Grant agreement No. 789056). TDS acknowledges the research project
was supported by the Hellenic Foundation for Research and Innovation (HFRI)
under the “2nd Call for HFRI Research Projects to support Faculty Members
& Researchers” (Project Number: 3382). OGM acknowledges financial support
from PAPIIT UNAM IN109123 and “Ciencia de Frontera” CONAHCYT CF-
2023-G100. CR acknowledges support from Fondecyt Regular grant 1230345
and ANID BASAL project FB210003. MJW acknowledges support from a Lev-
erhulme Emeritus Fellowship, EM-2021-064. MS acknowledges support by the
Ministry of Science, Technological Development and Innovation of the Repub-
lic of Serbia (MSTDIRS) through contract no. 451-03-66/2024-03/200002 with
the Astronomical Observatory (Belgrade). This work is based on observations
made with the NASA/ESA/CSA James Webb Space Telescope. The data were
obtained from the Mikulski Archive for Space Telescopes at the Space Tele-
scope Science Institute, which is operated by the Association of Universities for
Research in Astronomy, Inc., under NASA contract NAS 5-03127 for JWST.
These observations are associated with program 1670. This research has made
use of the NASA/IPAC Extragalactic Database (NED), which is operated by
the Jet Propulsion Laboratory, California Institute of Technology, under contract
with the National Aeronautics and Space Administration. This work has made
extensive use of Python (v3.9.12), particularly with ASTROPY (v5.3.3, http://
//www.astropy.org; Astropy Collaboration 2013, 2018), LMFIT (v1.2.2), PYNEB
(v1.1.18; Luridiana et al. 2015), MATPLOTLIB (v3.8.0; Hunter 2007), and NUMPY
(v1.26.0; Harris et al. 2020).

References

Allen, M. G., Groves, B. A., Dopita, M. A., Sutherland, R. S., & Kewley, L. J. 2008, ApJS, 178, 20

Alonso-Herrero, A., García-Burillo, S., Hönig, S. F., et al. 2021, A&A, 652, A99

Alonso Herrero, A., García-Burillo, S., Pereira-Santaella, M., et al. 2023, A&A, 675, A88

Álvarez-Márquez, J., Labiano, A., Guillard, P., et al. 2023, A&A, 672, A108

Argyriou, I., Glasse, A., Law, D. R., et al. 2023, A&A, 675, A111

Armus, L., Charmandaris, V., Bernard-Salas, J., et al. 2007, ApJ, 656, 148

Armus, L., Lai, T., U. V., et al. 2023, ApJ, 942, L37

Astropy Collaboration (Robitaille, T. P., et al.) 2013, A&A, 558, A33

Astropy Collaboration (Price-Whelan, A. M., et al.) 2018, AJ, 156, 123

Bellocchi, E., Martín-Pintado, J., Rico-Villas, F., Martín, S., & Jiménez-Sierra, I. 2023, MNRAS, 519, L68

Bessiere, P. S., & Ramos Almeida, C. 2022, MNRAS, 512, L54

Bushouse, H., Eisenhamer, J., Dencheva, N., et al. 2023, https://doi.org/
10.5281/zenodo.8140011

Cazzoli, S., Márquez, I., Masegosa, J., et al. 2018, MNRAS, 480, 1106

Chiar, J. E., & Tielens, A. G. G. M. 2006, ApJ, 637, 774

Cresci, G., Mainieri, V., Brusa, M., et al. 2015a, ApJ, 799, 82

Cresci, G., Marconi, A., Zibetti, S., et al. 2015b, A&A, 582, A63

Dasyra, K. M., Ho, L. C., Netzer, H., et al. 2011, ApJ, 740, 94

Davies, R. I., Maciejewski, W., Hicks, E. K. S., et al. 2014, ApJ, 792, 101

Davies, R. I., Burtscher, L., Rosario, D., et al. 2015, ApJ, 806, 127

Davies, R., Baron, D., Shimizu, T., et al. 2020, MNRAS, 498, 4150

De Robertis, M. M., & Osterbrock, D. E. 1986, ApJ, 301, 727

Feltre, A., Gruppioni, C., Marchetti, L., et al. 2023, A&A, 675, A74

Filippenko, A. V., & Halpern, J. P. 1984, ApJ, 285, 458

Fiore, F., Feruglio, C., Shankar, F., et al. 2017, A&A, 601, A143

Gallagher, R., Maiolino, R., Belfiore, F., et al. 2019, MNRAS, 485, 3409

García-Bernete, I., Ramos Almeida, C., Landt, H., et al. 2017, MNRAS, 469, 110

García-Bernete, I., Rigopoulou, D., Alonso-Herrero, A., et al. 2022, A&A, 666, L5

García-Bernete, I., Alonso-Herrero, A., Rigopoulou, D., et al. 2024a, A&A, 681, L7

García-Bernete, I., Rigopoulou, D., Donnan, F. R., et al. 2024b, arXiv e-prints
[arXiv:2409.05686]

García-Burillo, S., Combes, F., Usero, A., et al. 2014, A&A, 567, A125

García-Burillo, S., Alonso-Herrero, A., Ramos Almeida, C., et al. 2021, A&A, 652, A98

Gardner, J. P., Mather, J. C., Abbott, R., et al. 2023, PASP, 135, 068001P

Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357

Harrison, C. M. 2017, Nat. Astron., 1, 0165

Harrison, C. M., & Ramos Almeida, C. 2024, Galaxies, 12, 17

Harrison, C. M., Alexander, D. M., Mullaney, J. R., & Swinbank, A. M. 2014, MNRAS, 441, 3306

Harrison, C. M., Costa, T., Tadhunter, C. N., et al. 2018, Nat. Astron., 2, 198

Hermosa Muñoz, L., Cazzoli, S., Márquez, I., et al. 2024, A&A, 683, A43

Hernán-Caballero, A., Spoon, H. W. W., Alonso-Herrero, A., et al. 2020, MNRAS, 497, 4614

Hervella Seoane, K., Ramos Almeida, C., Acosta-Pulido, J. A., et al. 2023, A&A, 680, A71

Ho, L. C. 2009, ApJ, 699, 638

Hunter, J. D. 2007, Comput. Sci. Eng., 9, 90

Igo, Z., Parker, M. L., Matzeu, G. A., et al. 2020, MNRAS, 493, 1088

Kormendy, J., & Ho, L. C. 2013, ARA&A, 51, 511

Labiano, A., Azzollini, R., Bailey, J., et al. 2016, SPIE Conf. Ser., 9910, 99102W

Luridiana, V., Morisset, C., & Shaw, R. A. 2015, A&A, 573, A42

Maiolino, R., Russell, H. R., Fabian, A. C., et al. 2017, Nature, 544, 202

Mendes de Oliveira, C., Cypriano, E. S., Sodré, L.Jr, & Balkowski, C. 2004, ApJ, 605, L17

Miralles-Caballero, D., Colina, L., & Arribas, S. 2012, A&A, 538, A61

Nesvadba, N. P. H., Bicknell, G. V., Mukherjee, D., & Wagner, A. Y. 2020, A&A, 639, L13

Pei, Y. C. 1992, ApJ, 395, 130

Pereira-Santaella, M., Diamond-Stanic, A. M., Alonso-Herrero, A., & Rieke, G. H. 2010, ApJ, 725, 2270

Pereira-Santaella, M., Álvarez-Márquez, J., García-Bernete, I., et al. 2022, A&A, 665, L11

Perna, M., Arribas, S., Catalán-Torrecilla, C., et al. 2020, A&A, 643, A139

Ramos Almeida, C., Bischetti, M., García-Burillo, S., et al. 2022, A&A, 658, A155

Rieke, G. H., Wright, G. S., Böker, T., et al. 2015, PASP, 127, 584

Riffel, R. A., Storchi-Bergmann, T., Riffel, R., et al. 2023, MNRAS, 521, 1832

Roy, N., Bundy, K., Nevin, R., et al. 2021, ApJ, 913, 33

Sharples, R. M., Longmore, A. J., Hawarden, T. G., & Carter, D. 1984, MNRAS, 208, 15

Shin, J., Woo, J.-H., Chung, A., et al. 2019, ApJ, 881, 147

Silk, J. 2013, ApJ, 772, 112

Smajić, S., Fischer, S., Zuther, J., & Eckart, A. 2012, A&A, 544, A105

Speranza, G., Ramos Almeida, C., Acosta-Pulido, J. A., et al. 2024, A&A, 681, A63

Spinoglio, L., Pereira-Santaella, M., Dasyra, K. M., et al. 2015, ApJ, 799, 21

Stein, Y., Bomans, D. J., Ferguson, A. M. N., & Dettmar, R. J. 2017, A&A, 605, A5

Stone, M., Veilleux, S., Meléndez, M., et al. 2016, ApJ, 826, 111

Storchi-Bergmann, T., & Schnorr-Müller, A. 2019, Nat. Astron., 3, 48

Thean, A., Pedlar, A., Kukula, M. J., Baum, S. A., & O'Dea, C. P. 2000, MNRAS, 314, 573

Thomas, A. D., Dopita, M. A., Shastri, P., et al. 2017, ApJS, 232, 11

Tommasin, S., Spinoglio, L., Malkan, M. A., et al. 2008, ApJ, 676, 836

Torres-Flores, S., Mendes de Oliveira, C., de Mello, D. F., et al. 2009, A&A, 507, 723

Veilleux, S., Maiolino, R., Bolatto, A. D., & Aalto, S. 2020, A&A Rev., 28, 2

Venturi, G., Cresci, G., Marconi, A., et al. 2021, A&A, 648, A17

Venturi, G., Treister, E., Finlez, C., et al. 2023, A&A, 678, A127

Véron-Cetty, M. P., & Véron, P. 2006, A&A, 455, 773

Wells, M., Pel, J. W., Glasse, A., et al. 2015, PASP, 127, 646

Wright, G. S., Rieke, G. H., Glasse, A., et al. 2023, PASP, 135, 048003

Wright, G. S., Wright, D., Goodson, G. B., et al. 2015, PASP, 127, 595

Zhuang, M.-Y., Ho, L. C., & Shangguan, J. 2019, ApJ, 873, 103

Zhang, L., Packham, C., Hicks, E. K. S., et al. 2024, arXiv e-prints
[arXiv:2409.09771]

<footnote n="1">¹ Centro de Astrobiología (CAB) CSIC-INTA, Camino bajo el Castillo s/n, 28692 Villanueva de la Cañada, Madrid, Spain</footnote>

<footnote n="2">² Instituto de Física Fundamental, CSIC, Calle Serrano 123, 28006 Madrid, Spain</footnote>

<footnote n="3">³ Department of Physics, University of Oxford, Keble Road, Oxford OX1 3RH, UK</footnote>

<footnote n="4">⁴ Observatorio Astronómico Nacional (OAN-IGN) – Observatorio de Madrid, Alfonso XII, 3, 28014 Madrid, Spain</footnote>

<footnote n="5">⁵ Instituto de Astrofísica de Canarias, C/ Vía Láctea s/n, 38205 La Laguna, Tenerife, Spain</footnote>

<footnote n="6">⁶ Departamento de Astrofísica, Universidad de La Laguna, 38205 La Laguna, Tenerife, Spain</footnote>
⁷ Max Planck Institute for Extraterrestrial Physics (MPE), Giessen-bachstr. 1, 85748 Garching, Germany

⁸ Department of Physics and Astronomy, University of Alaska Anchorage, Anchorage, AK 99508-4664, USA

⁹ Department of Physics and Astronomy, The University of Texas at San Antonio, 1 UTSA Circle, San Antonio, TX 78249, USA

¹⁰ Department of Physics, University of Alaska, Fairbanks, AL 99775-5920, USA

¹¹ School of Mathematics, Statistics and Physics, Newcastle University, Newcastle upon Tyne NE1 7RU, UK

¹² Kavli Institute for Particle Astrophysics & Cosmology (KIPAC), Stanford University, Stanford, CA 94305, USA

¹³ Departmento de Física de la Tierra y Astrofísica, Fac. de CC Físicas, Universidad Complutense de Madrid, 28040 Madrid, Spain

¹⁴ Instituto de Física de Partículas y del Cosmos IPARCOS, Fac. CC Físicas, Universidad Complutense de Madrid, 28040 Madrid, Spain

¹⁵ Cahill Center for Astrophysics, California Institute of Technology, 1216 East California Boulevard, Pasadena, CA 91125, USA

¹⁶ LERMA, Observatoire de Paris, Collège de France, PSL University, CNRS, Sorbonne University, France, Paris

¹⁷ Institute of Astrophysics, Foundation for Research and Technology-Hellas (FORTH), Heraklion 70013, Greece

¹⁸ School of Sciences, European University Cyprus, Diogenes street, Engomi 1516, Nicosia, Cyprus

¹⁹ School of Physics & Astronomy, University of Southampton, Highfield, Southampton SO17 1BJ, UK

²⁰ Instituto de Radioastronomía y Astrofísica (IRyA), Universidad Nacional Autónoma de México, Antigua Carretera a Pátzcuaro

8701, ExHda. San José de la Huerta, Morelia, Michoacán 58089, Mexico

²¹ National Astronomical Observatory of Japan, National Institutes of Natural Sciences (NINS), 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan

²² Department of Astronomy, School of Science, Graduate University for Advanced Studies (SOKENDAI), Mitaka, Tokyo 181-8588, Japan

²³ Telespazio UK for the European Space Agency (ESA), ESAC, Camino Bajo del Castillo s/n, 28692 Villanueva de la Cañada, Spain

²⁴ Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA

²⁵ Instituto de Estudios Astrofísicos, Facultad de Ingeniería y Ciencias, Universidad Diego Portales, Av. Ejército Libertador 441, Santiago, Chile

²⁶ Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing 100871, China

²⁷ LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Sorbonne Paris Citeé, 5 place Jules Janssen, 92195 Meudon, France

²⁸ Astronomical Observatory, Volgina 7, 11060 Belgrade, Serbia

²⁹ Sterrenkundig Observatorium, Universiteit Gent, Krijgslaan 281-S9, Gent 9000, Belgium

³⁰ Centro de Astrobiología (CAB), CSIC-INTA, Ctra. de Ajalvir, km 4, 28850 Torrejón de Ardoz, Madrid, Spain

³¹ Centre for Extragalactic Astronomy, Durham University, South Road, Durham DH1 3LE, UK
# Appendix A: Additional material of the non-parametric modelling

In this appendix we show the properties of the emission lines from the integrated spectra in different regions estimated using a non-parametric methodology (see Sect. 3.1). An example of the modelling of the integrated nuclear profile of [NeV] line with both methodologies is shown in Fig. A.2. On the one hand, with the non-parametric modelling we derive a $W80 \sim 260\,\text{km s}^{-1}$, $v50 \sim -47\,\text{km s}^{-1}$, and a $v98 \sim 306\,\text{km s}^{-1}$ (maximum in absolute value between v02 and v98). On the other hand, with the parametric modelling we used a single Gaussian with $\sigma \sim 114\,\text{km s}^{-1}$ (i.e. FWHM $\sim 268\,\text{km s}^{-1}$), and mean velocity $\sim -62\,\text{km s}^{-1}$. Thus for this integrated profile of [NeV] the derived parameters are consistent with both methods. We note that we have made a cut of S/N > 3 at the peak of the line for measuring the properties of exclusively the ionised gas emission lines in Table A.1. Additional lines detected in the nuclear spectrum, such as molecular transitions, are reported in Zhang & Packham (2024) using a different methodology.

As mentioned in Sect. 3, for the nuclear region, we additionally applied a Monte Carlo approach to estimate the uncertainty of the W80 measurement (i.e. measuring the v10 and v90 parameters) produced exclusively by the noise in the continuum. We varied the original continuum from the integrated spectra with random Gaussian noise, with $\sigma$ equal to the standard deviation of the original spectra. Then we changed the condition from defining the line as where the S/N was > 3, to instead varying the line integration range from velocities in the range of $\pm 500$ to $\pm 1000\,\text{km s}^{-1}$ from the systemic velocity. We measured the v10 and v90 parameters for 30000 iterations for all the lines except for [Fe VIII] (only 5000 iterations) due to the strong features in the continuum. This approach allows to estimate the impact of the noise into the measurements, constraining the obtained W80 values. For the nuclear region the W80 values derived with both methods are consistent within the errors.

In general the integrated profiles in Fig. A.1 are symmetrical for the nuclear region (middle panels), and show redshifted and blueshifted wings for the two selected regions in the cone (NE and SW, respectively, see Sect. 3). These wings are particularly prominent for the NE region, as can be seen in Table A.1, where the v98 parameter tends to be higher than for the nuclear region. For estimating the significance of the relation of these parameters with respect to the IP and the critical density, we computed the correlation coefficients $\rho$ and p-value (see Table A.2). The Pearson coefficient, $\rho$, quantifies the strength of the relation between the variables, and the p-value indicates the probability of observing a correlation between the variables if there is no actual relation between them. None of the correlations are strong ($\rho > 0.7$), but rather moderate (0.3 < $\rho$ < 0.7) for the nuclear aperture, and weak ($\rho$ < 0.3) for the NE and SW apertures with some exceptions. We see a moderate correlation statistically significant (p-value < 0.05) for the W80 parameter with respect to the critical density in the nuclear region, and for the W80 parameter with respect to the IP in the SW. We did not detect any strong correlation of v98 with respect to the IP or the critical density of the lines in any of the regions (see Table A.2 and Figs. A.3 and 3).

# Appendix B: Additional spatially resolved maps

In this appendix we show the continuum maps of channel 3 (see Fig. B.1), kinematic maps derived using a two component modelling for the emission lines (see Figs. B.2 and B.3, and Sect. 3),

and the channel maps of the neon lines (see Figs. B.4, B.5, and B.6).

When using two Gaussians to fit the [NeV] line (Fig. B.3), the secondary component is distributed in two regions, one at $\sim 1.5''$ NE from the centre, and the second expanding to the SW direction up to $\sim 2.5''$. The velocity dispersion is similar to that of the primary ($\sigma \le 300\,\text{km s}^{-1}$, see Sect. 3.3.2), which is relatively low, as also found for outflows in other AGN (see e.g. Hermosa Muñoz et al. 2024). The region to the NE of the secondary component is completely blueshifted (average velocity $\sim -200\,\text{km s}^{-1}$), whereas the region to the south is both blueshifted and redshifted to the south and west, respectively. We note that these regions are detected with similar properties as the secondary component modelled for the [NeVI], [OIV] 25.89µm (not shown here), and [SIV] emission lines (see García-Bernete & Rigopoulou 2024b). As mentioned in Sect. 3.3.2, we detected approaching and receding velocities everywhere in the bicone, as is seen from the channel maps in Fig. B.6. Thus it is likely that the two components are a consequence of the cone being optically thin in these lines, and we are detecting blue and red wings in the data cubes due to the inclination of the cone with respect to the plane of the sky and its wide angle (see also Roy et al. 2021).
<table><caption>**Table A.1.** Measurements of the W80, v50, and v98 parameters of the ionised gas emission lines in the integrated nuclear spectra (R~0.5").</caption><thead><tr><th>Line</th><th>λ<br>(μm)</th><th>v50<sub>AGN</sub><br>(km s<sup>-1</sup>)</th><th>v50<sub>NE</sub><br>(km s<sup>-1</sup>)</th><th>v50<sub>SW</sub><br>(km s<sup>-1</sup>)</th><th>W80<sub>AGN</sub><br>(km s<sup>-1</sup>)</th><th>W80<sub>NE</sub><br>(km s<sup>-1</sup>)</th><th>W80<sub>SW</sub><br>(km s<sup>-1</sup>)</th><th>v98<sub>AGN</sub><br>(km s<sup>-1</sup>)</th><th>v98<sub>NE</sub><br>(km s<sup>-1</sup>)</th><th>v98<sub>SW</sub><br>(km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>[Fe II]</td><td>5.34</td><td>9.3</td><td>-168.8</td><td>-</td><td>267.1</td><td>89.1</td><td>-</td><td>168.8</td><td>213.3</td><td>-</td></tr><tr><td>[Fe VIII]</td><td>5.45</td><td>27.5</td><td>-</td><td>-</td><td>305.6</td><td>-</td><td>-</td><td>289.4</td><td>-</td><td>-</td></tr><tr><td>[Mg V]</td><td>5.61</td><td>-69.4</td><td>-196.6</td><td>-</td><td>339.1</td><td>296.7</td><td>-</td><td>366.1</td><td>493.2</td><td>-</td></tr><tr><td>[Ar II]</td><td>6.98</td><td>-47.9</td><td>-150.1</td><td>21.0</td><td>272.3</td><td>340.4</td><td>204.2</td><td>252.2</td><td>354.3</td><td>217.3</td></tr><tr><td>[Ne VI]</td><td>7.65</td><td>-58.1</td><td>-159.1</td><td>-6.7</td><td>302.9</td><td>252.5</td><td>252.5</td><td>411.5</td><td>310.5</td><td>259.2</td></tr><tr><td>[Ar V]</td><td>7.90</td><td>-24.0</td><td>-</td><td>-</td><td>342.3</td><td>-</td><td>-</td><td>219.6</td><td>-</td><td>-</td></tr><tr><td>[Ar III]</td><td>8.99</td><td>-33.1</td><td>-250.0</td><td>-32.2</td><td>300.8</td><td>214.9</td><td>257.8</td><td>248.0</td><td>376.9</td><td>248.0</td></tr><tr><td>[S IV]</td><td>10.51</td><td>-62.3</td><td>-172.6</td><td>12.1</td><td>294.1</td><td>257.3</td><td>257.3</td><td>429.9</td><td>319.6</td><td>318.7</td></tr><tr><td>[Ne II]</td><td>12.81</td><td>-24.5</td><td>-140.5</td><td>34.3</td><td>289.9</td><td>289.9</td><td>231.9</td><td>314.5</td><td>488.4</td><td>313.6</td></tr><tr><td>[Ar V]</td><td>13.10</td><td>-50.0</td><td>-</td><td>-</td><td>226.8</td><td>-</td><td>-</td><td>163.5</td><td>-</td><td>-</td></tr><tr><td>[Ne V]</td><td>14.32</td><td>-47.4</td><td>-151.2</td><td>-46.5</td><td>259.4</td><td>363.2</td><td>311.3</td><td>306.8</td><td>358.7</td><td>305.9</td></tr><tr><td>[Ne III]</td><td>15.55</td><td>-26.1</td><td>-169.4</td><td>22.6</td><td>334.4</td><td>334.4</td><td>286.6</td><td>312.7</td><td>408.2</td><td>261.4</td></tr><tr><td>[S III]</td><td>18.71</td><td>-33.1</td><td>-230.5</td><td>56.3</td><td>285.9</td><td>285.9</td><td>190.6</td><td>325.8</td><td>421.1</td><td>246.8</td></tr><tr><td>[Ne V]<sup>†</sup></td><td>24.32</td><td>4.3</td><td>-142.4</td><td>5.2</td><td>366.7</td><td>366.7</td><td>440.0</td><td>362.4</td><td>362.4</td><td>434.8</td></tr><tr><td>[O IV]</td><td>25.89</td><td>44.7</td><td>-93.1</td><td>45.6</td><td>344.4</td><td>344.4</td><td>344.4</td><td>299.7</td><td>344.4</td><td>321.1</td></tr></tbody></table>

**Notes.** † indicates that the continuum near the line for the NE and SW regions contained strong wiggles and thus the measurements are not well constrained (see Fig. A.1). The v98 parameter is the maximum in absolute value between the v02 and v98 measurements. The different integrated regions nuclear-AGN, NE, and SW, are defined in see Sect. 3.

<table><caption>**Table A.2.** Correlation coefficients between the W80 and v98 parameters with the lines IP and critical density at the three selected regions (see Sect. 3).</caption><thead><tr><th rowspan="2">Relation</th><th colspan="2">AGN</th><th colspan="2">NE</th><th colspan="2">SW</th></tr><tr><th>ρ<sub>AGN</sub></th><th>p-value</th><th>ρ<sub>NE</sub></th><th>p-value</th><th>ρ<sub>SW</sub></th><th>p-value</th></tr></thead><tbody><tr><td>W80 vs IP</td><td>0.51</td><td>0.08</td><td>0.39</td><td>0.21</td><td>0.65</td><td>0.03</td></tr><tr><td>v98 vs IP</td><td>0.20</td><td>0.50</td><td>0.24</td><td>0.45</td><td>0.14</td><td>0.68</td></tr><tr><td>W80 vs log n<sub>crit</sub></td><td>0.69</td><td>0.01</td><td>0.13</td><td>0.70</td><td>0.54</td><td>0.09</td></tr><tr><td>v98 vs log n<sub>crit</sub></td><td>0.50</td><td>0.09</td><td>0.64</td><td>0.03</td><td>0.03</td><td>0.94</td></tr></tbody></table>

**Notes.** The Pearson coefficient, *ρ*, quantifies the strength of the relation between the variables, and the p-value indicates the probability of observing a correlation between the variables if there is no actual relation between them.
<figure><img src="image_8.png" /><figcaption>**Fig. A.1.** Same as Fig. 2 but, from top to bottom, we show all the lines with S/N > 3 distributed into the four individual channels of MIRI/MRS data (see Sect. 2). All the velocities were corrected from the systemic value assuming a redshift of 0.00868 (see Sect. 1).</figcaption></figure>
<figure><img src="image_9.png" /></figure>

**Fig. A.2.** Examples on the modelling of the [Ne V]₁₄ line with the non-parametric (left) and the parametric techniques (right). **Left panel:** Normalised flux of the line integrated in the nuclear region (R = 0.5′′) with respect to the velocity. Grey, dashed, vertical lines indicate the limits of the line for doing the non-parametric modelling (see text in Appendix A). Red, dashed, vertical lines indicate the different parameters described in Sect. 2. **Right panel:** Gaussian modelling of the integrated nuclear region. The yellow line is the linear fit to the continuum, the green curve indicates the primary component, and the red curve indicates the global fit. The grey, dashed, horizontal lines in the residuals indicate the 3ε limits (see Sect. 2). The σ is corrected from the instrumental value (see Sect. 2). In both panels the central, grey, dashed, vertical line indicates the expected position of the line in rest frame, corrected from the systemic value with a redshift of 0.00868.

<figure><img src="image_10.png" /></figure>

**Fig. A.3.** Same as Fig. 3 but for the NE (top panels) and SW regions (bottom panels; see Sect. 3.1).
<figure><img src="image_11.png" /><figcaption>**Fig. B.1.** Continuum images for all the bands in channel 3 (upper panels from left to right, short, medium, and long, respectively) and the short band in channel 2 (bottom panel). The flux is in MJy sr⁻¹ and the regions correspond to those in Fig. 1.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Fig. B.2.** Maps of the two components model fit to the emission lines, following the description of Fig. 1. From left to right we show the flux maps in erg s⁻¹ cm⁻², the velocity, and the velocity dispersion in km s⁻¹. The photometric centre measured in the continuum of each channel is marked with a white star, the grey circle in the middle panels is the PSF, and the black line in the left panels is a 1" physical scale.</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>Fig. B.3. Same description as in Fig. B.2 but for [Ne V].</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>Fig. B.4. Channel maps of the [Ne II] line assuming steps of 200 km s⁻¹. The white star is the same as Fig. 1. White circles are the regions in which we obtained the spectra for the NE and SW parts of the ionisation cone, and the red circle is the main clump spectra (see Sect. 3).</figcaption></figure>
<figure><img src="image_15.png" /><figcaption>**Fig. B.5.** Same as Fig. B.4 but for [Ne III] line.</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>**Fig. B.6.** Same as Fig. B.4 but for [Ne V]$_{14}$ line.</figcaption></figure>
<figure><img src="image_17.png" /><figcaption>**Fig. B.7.** Same as Fig. B.4 but for [Ne VI] line. We excluded velocities above/below ± 400 km s⁻¹ to eliminate part of the PAH emission, seen as emission in the direction of the disc, where they are more prominent (see García-Bernete & Rigopoulou 2024b).</figcaption></figure>