# Molecular gas stratification and disturbed kinematics in the Seyfert galaxy MCG-05-23-16 revealed by JWST and ALMA

D. Esparza-Arredondo¹,²,* , C. Ramos Almeida¹,² , A. Audibert¹,², M. Pereira-Santaella³, , I. García-Bernete⁴,
S. García-Burillo⁵ , T. Shimizu⁶, R. Davies⁶ , L. Hermosa Muñoz⁴, A. Alonso-Herrero⁴, F. Combes⁷,
G. Speranza³, L. Zhang⁸, S. Campbell⁹, E. Bellocchi¹⁰,¹¹, A. J. Bunker¹², T. Díaz-Santos¹³,¹⁴,

B. García-Lorenzo¹,², O. González-Martín¹⁵, E. K. S. Hicks⁸,¹⁶,¹⁷, A. Labiano¹⁸, N. A. Levenson¹⁹, C. Ricci²⁰,²¹,
D. Rosario⁹, S. Hoenig²², C. Packham⁸, M. Stalevski²³,²⁴, L. Fuller⁹, T. Izumi²⁵, E. López-Rodríguez²⁶,²⁷,
D. Rigopoulou⁴,¹⁴, D. Rouan²⁷, and M. Ward²⁸,²⁹

*(Affiliations can be found after the references)*

Received 4 October 2024 / Accepted 17 November 2024

## ABSTRACT

Understanding the processes that drive the morphology and kinematics of molecular gas in galaxies is crucial for comprehending star formation and, ultimately, galaxy evolution. Using data from the Galactic Activity, Torus and Outflow Survey (GATOS) obtained with the James Webb Space Telescope (JWST) and the archival data from the Atacama Large Millimeter/submillimeter Array (ALMA), we study the behavior of the warm molecular gas at temperatures of hundreds of Kelvin and the cold molecular gas at tens of Kelvin in the galaxy MCG-05-23-16, which hosts an active galactic nucleus (AGN). Hubble Space Telescope (HST) images of this spheroidal galaxy, classified in the optical as S0, show a dust lane resembling a nuclear spiral and a surrounding ring. These features are also detected in CO(2–1) and H₂, and their morphologies and kinematics are consistent with rotation plus local inward gas motions along the kinematic minor axis in the presence of a nuclear bar. The H₂ transitions 0-0 S(3), 0-0 S(4), and 0-0 S(5), which trace warmer and more excited gas, show more disrupted kinematics than 0-0 S(1) and 0-0 S(2), including clumps of high velocity dispersion (of up to ~160 km s⁻¹), in regions devoid of CO(2–1). The kinematics of one of these clumps, located ~350 pc westward of the nucleus, are consistent with outflowing gas, possibly driven by localized star formation traced by polycyclic aromatic hydrocarbon emission at 11.3 µm. Overall, we observe a stratification of the molecular gas, with the colder gas located in the nuclear spiral, ring, and connecting arms, and most of the warmer gas with a higher velocity dispersion filling the inter-arm space. The compact jet, approximately 200 pc in size, detected with Very Large Array (VLA) observations, does not appear to significantly affect the distribution and kinematics of the molecular gas, possibly due to its limited intersection with the molecular gas disk.

**Key words.** galaxies: active – galaxies: ISM – galaxies: individual: MCG-05-23-16 – galaxies: kinematics and dynamics – galaxies: nuclei – galaxies: Seyfert

## 1. Introduction

A wealth of observational evidence supports the coevolution of active galactic nuclei (AGNs) and their host galaxies. This includes various correlations between the supermassive black hole (SMBH) mass ($M_{BH}$) and different galaxy properties, such as the velocity dispersion and the mass of the bulge (e.g., Kormendy & Richstone 1995; Magorrian et al. 1998; Ferrarese & Merritt 2000; Tremaine et al. 2002). These correlations are thought to be driven by AGN feedback, defined as the impact that growing SMBHs have on the host galaxies in which they reside. This impact can happen through radiation, winds, and jets (see Harrison & Ramos Almeida 2024, for a recent review). In cosmological simulations, AGN feedback plays a crucial role in, for example, quenching star formation, which limits the number of massive galaxies (Dubois et al. 2016). Furthermore, the energy released by the AGN affects the fueling of the SMBH itself, possibly regulating the AGN duty cycle (Ineson et al. 2015; Olivares et al. 2022). Active galactic nuclei are thought to be a phase of limited duration (~0.1–100 Myr; Martini 2004; Novak et al. 2011; Hickox et al. 2014) that might occur in every galaxy hosting a SMBH, provided that gas is supplied to the innermost regions of galaxies. This gas can be provided by galaxy mergers, secular processes, and/or direct gas accretion from filaments. Still, regardless of the mechanism, the gas has to lose most of its angular momentum to reach the sphere

of influence of the SMBH. This can happen through nuclear bars and mini-spirals, as has been shown by high-angular-resolution observations of molecular gas in nearby AGNs (Combes 2003; García-Burillo et al. 2005; García-Burillo & Combes 2012; Davies et al. 2009; Hicks et al. 2009; Audibert et al. 2019; Kolcu et al. 2023; Bianchin et al. 2024).

At low redshifts, molecular gas is the material from which stars form, constituting a significant fraction of the interstellar medium (ISM). In the local Universe, using near-infrared observations, it is possible to observe H₂ emission lines associated with rovibrational transitions that probe molecular gas at temperatures of T > 1000 K. Their ratios permit one to discriminate between thermal and nonthermal excitation, as well as to determine the thermal excitation temperature (Davies et al. 2004, 2005; Rodríguez-Ardila et al. 2005; Ramos Almeida et al. 2009). The problem is that the fraction of gas at these high temperatures is typically very low compared with the total molecular gas mass (hot-to-cold gas mass ratios of 10⁻⁷–10⁻⁵ for nearby star-forming galaxies; Dale et al. 2005; Davies et al. 2006; Pereira-Santaella et al. 2018). For this reason, most studies use the transitions of carbon monoxide (CO) and hydroxyl (OH) present in the (sub)-millimeter and far-infrared as indirect tracers of molecular gas. The CO emission can be accessed, for example, with the Atacama Large Millimeter/submillimeter Array (ALMA) and Northern Extended Millimetre Array (NOEMA) observatories and corresponds to gas at temperatures of approximately 5–20 K for disks and 50–75 K

<footnote>* Corresponding author; dcaesparza@gmail.com</footnote>
for nuclear regions (see Cazaux & Tielens 2002; Dale & Helou 2002; Habart et al. 2005, for reviews). Meanwhile, observations of the absorption lines of OH at 119 µm and 163 µm in nearby galaxies were obtained with the Herschel Space Observatory (Sturm et al. 2011; Stone et al. 2016; Runco et al. 2020). The measured intensities through these indirect molecular gas tracers can then be converted to H₂ column densities using a factor, αCO or XOH. These densities are determined under the assumption of identical physical conditions (temperature and density) in all clouds. As a result, they can vary by a factor of ~10 depending on whether the gas is optically thin or exposed to Galactic excitation conditions (see Herrera-Camus et al. 2020 and references therein). With the advent of space telescopes such as the Infrared Space Observatory (ISO) and Spitzer Space Telescope, the pure rotational H₂ emission lines became observable in the mid-infrared (Roussel et al. 2007; Pereira-Santaella et al. 2013, 2014). These lines trace molecular gas at temperatures between ~100 and 1000 K, commonly known as warm molecular gas. By combining the information they provide with CO, OH, and the hot H₂ gas, it becomes possible to obtain a complete and comprehensive picture of molecular gas emission in galaxies and to study their physical properties, including temperature, density, and column density (e.g., Davies et al. 2006).

The integral field capabilities of the James Webb Space Telescope (JWST; Gardner et al. 2023), combined with ALMA, now give us the opportunity to study the molecular phase of the ISM in galaxies with unprecedented sensitivity, spectral, and spatial resolution, making it possible to characterize the molecular gas kinematics in a spatially resolved manner. In the case of nearby AGNs, this permits us to advance our understanding of the fueling and feedback mechanisms (e.g., Pereira-Santaella et al. 2022; García-Bernete et al. 2022; Bianchin et al. 2024; Davies et al. 2024; García-Burillo et al. 2021, 2024; Zhang & Ho 2023; Hermosa Muñoz et al. 2024).

In this work, we target the nearby active galaxy MCG-05-23-16, which is part of the Galactic Activity, Torus, and Outflow Survey (GATOS). The aim of this collaboration is to understand the gas flow cycle of local AGNs, the properties and evolution of the obscuring material, and the relationship between star formation and nuclear activity in a nearly complete sample of AGNs (García-Burillo et al. 2021, 2024; Alonso-Herrero et al. 2021; García-Bernete et al. 2024a,b; Zhang et al. 2024a).

MCG-05-23-16 is an spheroidal galaxy (S0) hosting a Seyfert 1.9 nucleus (Véron-Cetty & Véron 2006; Davies et al. 2020), with broad recombination lines detected in the near-infrared (Goodrich et al. 1994) and a bolometric luminosity of $log(L_{bol}) [erg s^{-1}] \approx 44.3$ (Davies et al. 2020). This source has been widely studied in X-rays due to its continuum variability (e.g., Reeves et al. 2007; Zoghbi et al. 2014; Marinucci et al. 2022; Liu et al. 2024). The galaxy has also been observed at optical wavelengths with the Hubble Space Telescope (HST; see Fig. 1). Previous studies report that its [O III] 5007 Å and Hα emission are elongated on both sides of the nucleus with a position angle (PA) of 40° (e.g., Ferruit et al. 2000; Prieto et al. 2014). Ruschel-Dutra et al. (2021) explored the [O III] and [N III] emission and kinematics using IFU data from the Gemini Multi-Object Spectrographs (GMOS) and reported a kinematic major axis PA = 66°. They did not find clear signatures of ionized gas outflows for this source. More recently, using the same JWST dataset employed here, Zhang et al. (2024b) reported tentative evidence of ionized gas outflows, especially along the direction of their AGN ionization cones (PA = 172°) with velocities >100 km/s in this target.

Previous studies using mid-infrared images of the continuum emission from ground-based telescopes classified this source as spatially unresolved (e.g., Asmus et al. 2015), and reported good agreement between the ground- and space-based photometric and spectroscopic observations at these wavelengths (see Sasmirala object pages). This indicates that no significant flux or spectral changes happened at these wavelengths over a period of ~25 years, and also that the mid-infrared emission is dominated by the AGN. Studies at radio frequencies found weak diffuse emission at 8.4 GHz with the Very Large Array (VLA), which is marginally resolved with PA ~169° and associated with a compact radio jet (Mundell et al. 2009; Orienti & Prieto 2010). This paper reports and analyzes the distribution and kinematics of the warm and cold molecular gas in MCG-05-23-16. Using JWST and ALMA together for the first time, we study the H₂ pure rotational emission lines and the CO(2-1) emission line at the rest frequency of 230.538 GHz. Throughout this work, we assume a cosmology with $H_0 = 70$ km s⁻¹ Mpc⁻¹, $\Omega_M = 0.27$, and $\Omega_\Lambda = 0.73$. The redshift of the galaxy is z = 0.008486, which corresponds to a spatial scale of 176 pc/arcsec and a luminosity distance of 37 Mpc (Kawamuro et al. 2016).

## 2. Observations and data reduction

### 2.1. JWST/MIRI observations

Observations of MCG-05-23-16 were obtained as part of the Cycle 1 GO program 1670 (PIs: T. Shimizu and R. Davies) using the medium-resolution spectrometer (MRS) mode of the Mid-Infrared Instrument (MIRI, Rieke et al. 2015; Wright et al. 2015, 2023) on board JWST. MCG-05-23-16 was observed on April 3, 2023, with a single pointing and an on-source exposure time of 1121 s in each of the four MRS sub-bands using the FASTR1 read-out mode and four-point dither pattern. Background observations half the exposure time of the science observations were also obtained using a two-point dither pattern.

The MIRI/MRS data covers the spectral range 4.9–28.1 µm, with a spectral resolution of R ~ 3700–1300 (Labiano et al. 2021; Argyriou et al. 2023). The data are split into four channels that cover the following ranges: 4.9–7.65 µm (Channel 1; hereafter Ch1), 7.51–11.71 µm (Channel 2; Ch2), 11.55–18.02 µm (Channel 3; Ch3), and 17.71–28.1 µm (Channel 4; Ch4). Combining these four channels and the three grating settings available (short, medium, and long), there are 12 different wavelength bands. We measured the angular resolutions for the 12 bands from the full-width at half maximum (FWHM) of the star HD 163466 (program ID: 01050; PI: B. Vandenbussche), observed on June 9, 2022. These are reported in Table 1, together with the field of view (FoV) of each channel.

The data were reduced using the JWST Science Calibration Pipeline (version 1.11.4), with the context 1130 for the calibration References Data System (Labiano et al. 2016; Álvarez-Márquez et al. 2023; Bushouse et al. 2023). 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 is described in detail in García-Bernete et al. (2022, 2024a) and Pereira-Santaella et al. (2022). The background subtraction was done spectrum by spectrum for each science cube using background average spectra. These average spectra were generated through cubes created with the dedicated background-only observations of each sub-band for each channel. We also applied a residual fringing correction to each spectrum of the science cubes using the script `RFC1D-UTILS` provided by the JWST help desk (Argyriou et al. 2023). This correction is now available in the latest pipeline version.
<figure><img src="image_1.png" /><figcaption>**Fig. 1.** HST/WFC3 image and color map of MCG-05-23-16. Left: 10" × 10" image obtained using the filter F606W ($\lambda_{\text{ref}} \sim 5889$ Å). The JWST/MIRI/MRS FoV of Channel 4 (Ch4) is shown with a dashed white rectangle. The horizontal white line indicates the physical size of 1". Right: V-H color map using the F606W (left panel) and F160W ($\lambda_{\text{ref}} \sim 15$ 369 Å) images. The innermost region of the map has been masked using a solid black circle of 1.5 × FWHM of the F160W filter (1.5 × 0.18) to avoid artifacts due to the mismatch between the different PSFs of F606W and F160W. The dashed blue circles correspond to the regions used to extract the JWST/MIRI nuclear spectrum considering the correction per aperture (Ch1-s: diameter of 0.8 ≈ 141 pc – Ch4-l: 2.5 ≈ 440 pc). The dashed white circle corresponds to the region that extracts the extended (diameter of 3.6 ≈ 634 pc) JWST/MIRI spectrum. The nuclear and extended spectra are shown in Fig. 2 with solid orange and blue lines, respectively. The *Spitzer*/IRS SL (3.6 × 57", PA = 187°) slit is shown with dashed cyan lines. A dust lane resembling a nuclear spiral can be seen southwest and northeast of the nucleus. See also Fig. A.1.</figcaption></figure>

<table><caption>**Table 1.** FWHM and FoV of MIRI/MRS channels and bands.</caption><thead><tr><th>MRS channel and band</th><th>PSF FWHM (arcsec)</th><th>FoV (arcsec)</th><th>H<sub>2</sub> lines detected</th></tr></thead><tbody><tr><td>Ch1-short (Ch1-s)</td><td>0.278</td><td>3.2 × 3.7</td><td>...</td></tr><tr><td>Ch1-medium (Ch1-m)</td><td>0.290</td><td>3.2 × 3.7</td><td>...</td></tr><tr><td>Ch1-long (Ch1-l)</td><td>0.294</td><td>3.2 × 3.7</td><td>0-0S(5)</td></tr><tr><td>Ch2-short (Ch2-s)</td><td>0.359</td><td>4.0 × 4.8</td><td>0-0S(4)</td></tr><tr><td>Ch2-medium (Ch2-m)</td><td>0.379</td><td>4.0 × 4.8</td><td>0-0S(3)</td></tr><tr><td>Ch2-long (Ch2-l)</td><td>0.423</td><td>4.0 × 4.8</td><td>...</td></tr><tr><td>Ch3-short (Ch3-s)</td><td>0.503</td><td>5.2 × 6.2</td><td>0-0S(2)</td></tr><tr><td>Ch3-medium (Ch3-m)</td><td>0.541</td><td>5.2 × 6.2</td><td>...</td></tr><tr><td>Ch3-long (Ch3-l)</td><td>0.603</td><td>5.2 × 6.2</td><td>0-0S(1)</td></tr><tr><td>Ch4-short (Ch4-s)</td><td>0.760</td><td>6.6 × 7.7</td><td>...</td></tr><tr><td>Ch4-medium (Ch4-m)</td><td>0.823</td><td>6.6 × 7.7</td><td>...</td></tr><tr><td>Ch4-long (Ch4-l)</td><td>0.835</td><td>6.6 × 7.7</td><td>...</td></tr></tbody></table>

**Notes.** FWHM measured from the collapsed spectrum of the star HD 163466 in each of the 12 sub-channels of MIRI/MRS, and FoV of each channel. The pixel sizes employed to convert from pixels to arc-seconds are 0.13 in Ch1, 0.17 in Ch2, 0.20 in Ch3, and 0.35 in Ch4 (see MIRI Observations Modes). The last column indicates the channels where we detect the H<sub>2</sub> lines from S(5) to S(1), which have rest-frame wavelengths of 6.90952, 8.02505, 9.66492, 12.27861, and 17.03484 µm, respectively.

## 2.2. ALMA and VLA data

ALMA observations of MCG-05-23-16 in band 6 were performed in two execution blocks (project ID 2019.1.01742.S, PI: D. Rosario) for a total on-source integration time of 15.6 minutes in two configurations: compact (TC) and extended (TE). The TE configuration (baselines from 15 to 2516 m) was observed with 50 antennas in October 2019 and an on-source integra-

tion time of 11.6 minutes, resulting in a synthesized beam of 0.39 (~70 pc). The TC configuration (baselines from 15 to 314 m) was observed in November 2019 with 47 antennas and an on-source integration time of 4 minutes, providing a synthesized beam of 1.44 (~250 pc). The observations were centered on the nucleus, with a single pointing covering a FoV of 26". The correlator setup was centered at the observed frequency of the CO(2-1) line: 228.597 GHz, with a bandwidth of 1.875 GHz and channel spacing of 7.8 MHz (~10.2 km s<sup>-1</sup>). The phase calibrator J1007-3333 was used both in TE and TC, and J1037-2934 was used as flux calibrator in TE configuration and J0725-0054 in TC, from which we measured the typical 10% flux accuracy.

The data were calibrated, concatenated, and cleaned with the CASA software (version 5.6.1-8; McMullin et al. 2007). Once each configuration was calibrated, we concatenated the two of them using the CONCAT task. The imaging and cleaning were performed with the task TCLEAN. The spectral line map was obtained after subtracting the continuum in the uv-plane using the task UVCONTSUB specifying a zeroth-order polynomial in the channels free from emission lines. The CO(2-1) data cube was produced with a spectral resolution of 15 km s<sup>-1</sup> and using Briggs weighting mode (Briggs 1995) and a robust parameter set to 2 to achieve the best sensitivity, equivalent to a natural weighting. The resulting synthesized beam is 0.65 × 0.53 at PA of 87° and a root mean square (rms) noise of 0.51 mJy beam<sup>-1</sup> per 15 km s<sup>-1</sup> channel. Finally, the datacube was corrected for primary beam attenuation.

We also retrieved the X-band VLA radio image (8.4 GHz) of our target from the NRAO VLA Archive Survey Images Page, observed on August 21, 1999, in A configuration with a beam size of 0.33 (~60 pc; synthesized beam of 0.43 × 0.23) and an rms of 0.1 mJy. The data was first presented in Orienti & Prieto (2010), where the authors reported that the total radio spectrum slightly steepens at 8.4 GHz with a spectral index of 0.8 ± 0.1.
<figure><img src="image_2.png" /></figure>

**Fig. 2.** Mid-infrared integrated spectra of MCG-05-23-16. The solid blue and orange lines correspond to the nuclear and extended JWST/MIRI spectra, extracted in circular apertures of ~0.8–2.5 (146–441 pc) and 3.6 (634 pc) diameter, respectively (black and white circles in Fig. 1). The dotted magenta line represents the spectrum obtained by subtracting the nuclear from the extended spectrum. The dashed black line corresponds to the *Spitzer/IRS* spectrum, extracted in the same aperture as the extended JWST/MIRI spectrum. The middle and bottom panels show zooms of two spectral regions. The black filled circles mark the nuclear 12 µm and 18 µm continuum emission reported by Asmus et al. (2015), using data from VISIR/VLT.

# 3. Analysis and results

## 3.1. Nuclear and extended mid-infrared spectra

The left panel of Fig. 1 shows an HST/WFC3 optical image of the central 10" × 10" (~1.7 × 1.7 kpc²) of MCG-05-23-16 (proposal ID: 15181; PI: D. Rosario). The V-H color map using the F606W and F160W images is shown in the right panel. The blue and white circles show the regions used to extract the nuclear and extended spectra, shown as blue and orange lines in Fig. 2. The nuclear spectrum was extracted as a point source, with increasing form 1.5 × PSF FWHM in each channel (see Table 1). These radii range from 0.42 on Ch1-s to 1.252 on Ch4-l. The extended spectrum was extracted in an aperture of 1.8 radius centered in

the nucleus (diameter of 3.6 ≈ 634 pc) and including the nuclear aperture. The spectra from the four channels were combined using different scaling factors to align them in overlapping areas. These scaling factors were determined by conducting linear fits, and they represent less than 1.2% of the flux measured at any given wavelength in the final spectra shown in Fig. 2.

In Fig. 2, we also show the *Spitzer/IRS* spectra downloaded from CASSIS (Combined Atlas of Sources with *Spitzer* IRS Spectra; Lebouteiller et al. 2011). This spectrum was obtained by putting together the spectra from the IRS Short-Low (SL) and Long-Low (LL) slits, which have resolutions of R ~ 60–127 and 57–126, respectively. The wavelength range between 5.3 and 14 µm is obtained through the SL slits, which have a width of ~3.6 (634 pc), and the 14–38 µm range through the LL slits, of ~10.5 width (1.85 kpc). The size and PA of the SL slit are shown in the right panel of Fig. 1. Furthermore, we include the nuclear 12 µm and 18 µm continuum emission obtained through VISIR/VLT data and reported by Asmus et al. (2015) in Sasmirala catalog.

The continuum shape of the IRS and MIRI extended spectra are almost identical, as was expected, since they were extracted in the same aperture. Thanks to the superior spectral resolution and sensitivity of JWST/MIRI, a large number of previously undetected emission lines are observed in both the nuclear and extended MIRI spectra shown in Fig. 2. In this work, we are particularly interested in the H₂ rotational transitions that trace warm molecular gas at hundreds of Kelvin. In the MIRI spectra shown in Fig. 2, we detect H₂0–0S(1) 17.03 µm, S(2) 12.28 µm, S(3) 9.66 µm, S(4) 8.02 µm, and S(5) 6.90 µm. The more precise rest-frame wavelength of these emission lines and the MRS channel in which they are detected are indicated in Table 1.

The nuclear spectrum, extracted in a circular aperture of 146–634 pc in diameter, shows the same spectral shape as the extended spectrum, consistent with an AGN-dominated mid-infrared emission. The 10 µm silicate feature appears in shallow absorption in both the nuclear and extended-spectrum, indicating that the obscuring material producing it comes from the nuclear region. Considering that the silicate feature lies between Ch2-m and Ch2-l, this corresponds to a diameter of 1.33–1.27 (200–223 pc). The 18 µm silicate feature is practically absent (García-Bernete et al. 2024a).

## 3.2. Morphology and kinematics of the molecular gas

### 3.2.1. Warm versus cold molecular gas from parsec to kiloparsec scales

To understand the behavior of the molecular gas, we produced kinematic maps of the different H₂ lines using the cubes corresponding to Ch3-l, Ch3-s, Ch2-m, Ch2-s, and Ch1-l (see Table 1) and the ALUCINE tool (Peralta de Arriba et al. 2023). This tool creates flux, velocity, and velocity dispersion (σ) maps by fitting the emission lines using single or multiple Gaussian components and a local continuum on a spaxel-by-spaxel basis. The flux, velocity, and velocity dispersion maps of the five rotational transitions of the H₂ molecule considered here are shown in Fig. 3. We used a single Gaussian component to fit the emission line profiles. Since MCG-05-23-16 is a Seyfert 1.9 galaxy, the mid-infrared continuum dominates the central region of the cubes, making the equivalent width (EW) of the emission lines very low. This implies that some of the central spaxels show negative values when we subtract the continuum. To minimize this effect, we applied spatial smoothing to the cubes before running ALUCINE using the smooth tool of QFITSVIEW. We selected a boxcar average of one pixel along both the X and Y axes. After
applying this smoothing, only a few spaxels with negative values remained in the H₂O–OS(1) and S(2) maps; that is, those corresponding to Ch3 (see the first and second rows of Fig. 3).

Unlike the ionized gas traced by the various emission lines detected in the MIRI/MRS spectrum (see Fig. 2), which is dominated by the nuclear emission (see Zhang et al. 2024b for more details on the ionized gas kinematics), the warm molecular gas is extended over the entire FoV of the different channels. For comparison, in Fig. 3, we plot the contours of the CO(2–1) emission obtained from the ALMA data (see moment maps in Fig. B.1) with a solid magenta line and the radio continuum emission at 8.4 GHz from VLA with a solid black line (see Sect. 2.2 for a description of the ALMA and VLA observations). We find that the morphologies of the CO(2–1) and H₂ emission are similar, centrally peaked, and resembling a mini-spiral in the central ~2'' (~350 pc) that coincides with the dust lane seen in the HST optical image of the galaxy (see Fig. 1). This mini-spiral is connected with a ring-like feature of almost the same size as the FoV of Ch3 (~8'' ≃ 1.4 kpc). The northeast side of the ring is the most conspicuous in H₂, as can be seen from the S(1) and S(2) flux maps shown in Fig. 3.

Although Ch1, Ch2, and Ch3 cover different FoVs (see Table 1), the colder H₂ gas traced by the S(1) and S(2) lines (two top rows of Fig. 3) seems to exhibit a more regular rotation pattern than the warmer H₂ gas traced by S(3), S(4), and S(5). If we subtract the S(5) from the S(1) velocity map, we find residual velocities of up to ~100 km/s at certain locations, indicating important differences between the kinematics traced by the two lines. The rotation velocities reach ±250 km s⁻<fnref n="1" />, redshifted to the northeast and blueshifted to the southwest (see the middle column of Fig. 3). These velocities are similar to the low ionization lines (e.g., [Ne II]) shown in Zhang et al. (2024b) for this source, which is consistent with the kinematics of the warm ionized gas traced by the optical [O III] and [N II] emission lines (Ruschel-Dutra et al. 2021). The dispersion maps shown in the last column of the same figure show interesting structures, including knots of high velocity dispersion that reach values of up to 160–170 km s⁻<fnref n="1" />. These knots are mainly located in the regions between the ring and the nuclear spiral. An example is the S(3) knot at ~2'' ~350 pc westward of the nucleus. These areas of high velocity dispersion coincide with regions that appear depleted of CO(2–1). The VLA 8.4 GHz contours show extended emission to the north, perpendicular to the mini-spiral (PA ~170° and size ~200 pc). This slightly extended emission was interpreted as a compact radio jet by Orienti & Prieto (2010), but it does not seem connected with any of the kinematic features detected in the H₂ gas. In Sect. 3.2.3, we compare the kinematics of the warm and cold molecular gas in more detail.

### 3.2.2. H₂ molecular gas in specific regions

This section reports the results from comparing the H₂ emission line profiles extracted in six regions of 0.5' radius (~90 pc) corresponding to areas where high velocity dispersion values are visually identified for one or several H₂ lines. In Fig. 4, we identify these regions with circles of different colors in the S(3) and S(5) velocity dispersion maps. Regions R2 and R4 are the most nuclear ones, with R4 encompassing the nucleus and R2 located westward from it. R3 and R6 are situated northeast and southwest of the nucleus, respectively, and regions R1 and R5 to the east and west of the nucleus, with R5 exhibiting high dispersion in S(3), but the S(5) emission falls outside of the FoV in Ch1-l.

In Fig. 5, we show the corresponding emission line profiles extracted from each region. We modeled each line profile with

a single Gaussian and subtracted a local continuum, defined by interpolating between two regions adjacent to the emission line. In Table 2, we report the flux, velocity shift (ΔV), and intrinsic velocity dispersion (σ) values measured for each line. The integrated flux values of the H₂ lines in different regions range [0.36:62] ×10⁻¹⁶ erg s⁻<fnref n="1" /> cm⁻². The highest flux values are found in the region R4, which includes the nucleus. Notably, in this region, the S(3) and S(5) lines show the highest fluxes, while the lowest are measured for S(2) and S(4). The region R5, centered at the position of the S(3) high-velocity-dispersion knot at ~2.5'' westward of the nucleus, shows the lowest flux values (≤2.5 × 10⁻¹⁶ erg s⁻<fnref n="1" /> cm⁻²). When comparing the flux values of the five H₂ lines in each region, the S(1) emission line is the most luminous in regions R1, R2, and R3, and in R5 and R6, the S(3) flux values are maximum.

We also find that the emission line profiles of different H₂ emission lines show different shapes (see Fig. 5). The lines are blueshifted relative to systemic in regions R2, R4, R5, and R6. The region R6, located southwest of the nucleus, shows the highest velocities, |ΔV| > 120 km s⁻<fnref n="1" />. The maximum blueshift is observed for the S(3) emission in this region (−139 km s⁻<fnref n="1" />). In several radial profiles of these lines, we also observe the presence of wings. The S(1) and S(3) lines in region R5 show red wings (see also S(1) in region R3). We observe slightly blueshifted wings in the region R1 for all transitions, except by the radial profiles of the S(1). In regions R1 and R3, located to the east and northeast of the nucleus, we find velocities closer to zero, some blueshifted and others redshifted, for all transitions (see Table 2). Finally, the velocity dispersion values, corrected for the instrument broadening, range between 40 and 170 km s⁻<fnref n="1" />, increasing from S(1) to S(5) in all regions. The σ values measured for S(1) are ≤64 km s⁻<fnref n="1" />, while for S(5) are ≥140 km s⁻<fnref n="1" /> in all regions (see Table 2). This suggests that the S(1) and S(2) lines, which probe colder molecular gas, may trace gas that is more rotation-dominated, while the S(3) and S(5) lines may be associated with more dispersion-dominated gas.

### 3.2.3. Modeling the molecular gas kinematics

To further characterize and compare the kinematics of the molecular gas from temperatures of tens of Kelvin (CO) to hundreds of Kelvin (H₂), we chose three transitions; namely, CO(2–1), S(1), and S(3). Although S(4) and S(5) would generally be preferable for tracing hotter and more excited molecular gas, we decided not to use them due to the higher number of spaxels with no signal (see Fig. 3). We analyzed the gas kinematics using the “3D-Based Analysis of Rotating Objects from Line Observations” (³DBAROLO) software by Di Teodoro & Fraternali (2015). ³DBAROLO allows the user to model rotation in the observed emission line datacubes by fitting a 3D tilted ring. By then subtracting the fit rotation model from the data, it is possible to investigate the regions that deviate from rotation by inspecting the residuals and position-velocity diagrams (PVDs), as in Domínguez-Fernández et al. (2020) and Ramos Almeida et al. (2022). We used the SEARCH parameter in the 3DFIT task to find spatial and spectral coherent emissions, applying a mask with a threshold of 3σ<sub>rms</sub>. We fixed the central position to the peak of the continuum at 200 GHz for the CO(2–1) emission and to the peak of atomic emission lines lying in the same channels for the H₂ lines<fnref n="1" />: at (RA, Dec) of (9:47:40.1472,−30:56:55.922) for S(1) and (9:47:40.1506,−30:56:56.054) for S(3). We allowed

<footnote n="1">¹ In the case of S(1) we used the position of the peak of [Ne III] and in the case of S(3), of [Ar III].</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 3.** Flux [×10⁻²⁰ erg s⁻¹ cm⁻²] (left), velocity [km/s] (center), and velocity dispersion [km/s] (right) maps of the different H₂ rotational transitions detected in the MIRI/MRS data. The magenta contours represent 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, and 90% of the maximum ALMA CO(2–1) flux. For the S(1), S(2), and S(3) maps we also included the ALMA contours corresponding to 7% and 8% of the maximum flux. The black contours correspond to the 5%, 10%, 15%, and 20% of the average σ (~2.85 × 10⁻⁵) in the VLA 8.4 GHz image. The white star shows the AGN location according to the ALMA continuum. The VLA and ALMA beams, corresponding to 0.′14 × 0.′30 and 0.′68 × 0.′83, are shown in the S(1) and S(5) flux maps.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Fig. 4.** Velocity dispersion maps of the H₂0–0S(3) and S(5) emission lines with the CO(2–1) moment 0 contours superimposed. The units of the color bars are km s⁻¹. The circles of different colors indicate the regions from where we extracted the spectra shown in Fig. 5. These regions have radii of 0.′5 (88 pc). In Sect. 3.2.2, we describe how the regions were selected.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 5.** Emission line profiles of the five H₂ emission lines extracted from the regions with matching colors indicated in Fig. 4. All the velocities refer to the systemic value computed with z = 0.008486, indicated by the dashed vertical black lines.</figcaption></figure>

<table><caption>**Table 2.** Flux, velocity, and velocity dispersion values of the H₂ lines in the integrated regions (r = 0.′5) indicated in Fig. 4.</caption><thead><tr><th rowspan="2">H<sub>2</sub></th><th colspan="3">Region R1</th><th colspan="3">Region R2</th><th colspan="3">Region R3</th><th colspan="3">Region R4</th><th colspan="3">Region R5</th><th colspan="3">Region R6</th></tr><tr><th>Flux</th><th>ΔV</th><th>σ</th><th>Flux</th><th>ΔV</th><th>σ</th><th>Flux</th><th>ΔV</th><th>σ</th><th>Flux</th><th>ΔV</th><th>σ</th><th>Flux</th><th>ΔV</th><th>σ</th><th>Flux</th><th>ΔV</th><th>σ</th></tr></thead><tbody><tr><td>0–0S(1)</td><td>7.67 ± 0.78</td><td>10</td><td>62 ± 10</td><td>14.94 ± 1.50</td><td>−78</td><td>56 ± 10</td><td>11.99 ± 1.21</td><td>10</td><td>64 ± 10</td><td>36.31 ± 3.64</td><td>−34</td><td>64 ± 10</td><td>1.91 ± 0.20</td><td>−121</td><td>57 ± 31</td><td>7.12 ± 0.72</td><td>−121</td><td>41 ± 19</td></tr><tr><td>0–0S(2)</td><td>3.32 ± 0.35</td><td>0</td><td>87 ± 6</td><td>5.31 ± 0.54</td><td>−121</td><td>76 ± 6</td><td>3.21 ± 0.33</td><td>0</td><td>66 ± 6</td><td>20.82 ± 2.09</td><td>−60</td><td>&lt;41</td><td>0.36 ± 0.12</td><td>−121</td><td>&lt;84</td><td>3.27 ± 0.34</td><td>−121</td><td>71 ± 6</td></tr><tr><td>0–0S(3)</td><td>6.64 ± 0.68</td><td>21</td><td>120 ± 4</td><td>12.43 ± 1.26</td><td>−99</td><td>87 ± 25</td><td>7.54 ± 0.77</td><td>−19</td><td>100 ± 4</td><td>61.55 ± 6.17</td><td>−59</td><td>130 ± 4</td><td>2.46 ± 0.28</td><td>−59</td><td>163 ± 15</td><td>8.76 ± 0.89</td><td>−139</td><td>95 ± 4</td></tr><tr><td>0–0S(4)</td><td>2.45 ± 0.28</td><td>22</td><td>140 ± 37</td><td>3.15 ± 0.34</td><td>−75</td><td>112 ± 29</td><td>1.98 ± 0.21</td><td>−27</td><td>&lt;29</td><td>25.18 ± 2.54</td><td>−27</td><td>140 ± 16</td><td>0.59 ± 0.08</td><td>−123</td><td>&lt;131</td><td>1.74 ± 0.20</td><td>−123</td><td>94 ± 27</td></tr><tr><td>0–0S(5)</td><td>3.91 ± 0.41</td><td>−29</td><td>148 ± 37</td><td>6.73 ± 0.69</td><td>−98</td><td>139 ± 30</td><td>3.91 ± 0.41</td><td>−29</td><td>147 ± 3</td><td>61.98 ± 6.22</td><td>−64</td><td>157 ± 3</td><td>...</td><td>...</td><td>...</td><td>3.63 ± 0.38</td><td>−133</td><td>170 ± 3</td></tr></tbody></table>

**Notes.** The flux values are in units of 10⁻¹⁶ erg s⁻¹ cm⁻², and the velocity and velocity dispersion values in km s⁻¹. The fluxes have been corrected from extinction by considering the value of A<sub>V</sub> = 1 mag from Prieto et al. (2014) and the Cardelli et al. (1989) reddening law, and the values of σ have been corrected from instrumental width using the values reported in Argyriou et al. (2023).

the code to fit the observed rotation velocity and velocity dispersion while also varying the PA and the inclination of the disk.

We show the moment maps of the CO(2–1) emission, the corresponding ³D BAROLO model, and the residuals resulting from subtracting the model from the data in Fig. B.1. The

moment 0 map shows the mini-spiral and the ring also detected in H₂, and the rotation pattern can be reproduced to a certain extent with a disk inclination of i = 75° and kinematic major axis PA = 59° (mean values). Significant residuals can be seen in the right panels of Fig. B.1, indicating important deviations
<figure><img src="image_6.png" /></figure>

**Fig. 6.** PVDs of the S(1) line. Left panels: Flux and velocity dispersion maps obtained using bilinear interpolation, with the regions used to extract the PVDs superimposed, which have a size of 0.6 × 7". Middle and right panels: PVDs extracted along the kinematic major and minor axes determined by 3D BAROLO (top) and along PA = 87° and 3° (bottom). The black and magenta contours correspond to the data and rotating disk model above 2σ (see Sect. 3.2.3). The high-velocity emissions observed between −1 and 1 arcsec correspond to artifacts in the nuclear region because of the strong AGN continuum and derived low EW of the molecular emission lines there.

from rotation. We provide an interpretation for these noncircular motions, taking into account the galaxy inclination in Sect. 4.

In Fig. B.2, we show the observed velocity field, rotation model, and corresponding residuals for the three molecular transitions considered here. In the case of the H₂ lines, the disk models fit with 3D BAROLO have inclinations of i = 33° and 37° and major axes PAs of 56° and 61° for S(1) and S(3). The lower inclination of the H₂ disk compared to that derived for the CO is probably related to the smaller FoV of the MIRI/MRS data, which ends in the ring and does not cover the extended and more diffuse emission probed by ALMA on larger scales (see Fig. B.2). The rotation pattern is similar for the three transitions considered here, with blueshifted (redshifted) velocities to the southwest (northeast). If we compare the amplitude of the rotation velocities in the same regions, we find similar values among the five H₂ lines (see Figs. 3 and B.2 and Table 2). The residuals in the right panels of Fig. B.2 show high velocities along the minor axis in the outer part of the ring for CO(2–1) and S(1), indicating the presence of noncircular motions. In the S(3) velocity map, the high-velocity-dispersion knot at ~2.5" westward of the nucleus shows “forbidden” velocities; that is, redshifted velocities in the blueshifted rotating side of the galaxy (see bottom panels in Fig. B.2). The high-velocity residuals in this region, which appear devoid of CO(2–1), highlight how it deviates from regular rotation.

In Fig. B.3, we show the velocity dispersion maps, models, and residuals resulting from 3D BAROLO for the three transitions here considered. The warm molecular gas traced by S(1) and S(3) shows higher velocity dispersion than CO in the nucleus (~130–140 km s⁻¹), and in the case of S(3) also in R5 (~160 km s⁻¹). The CO(2–1) map reveals extended regions of high velocity dispersion throughout the disk, including areas that are outside the MIRI/MRS FoVs. 3D BAROLO is unable to reproduce these high central dispersion regions, and therefore the residuals are large in the center for the warm molecular gas

and across the disk for CO(2–1). In the case of the MIRI data, the high nuclear velocity dispersions might be partly due to beam smearing, and we note that BAROLO only reproduces maximum values of 80–90 km/s. The high-velocity-dispersion knot in R5 stands out in the residual map of S(3).

Finally, to further investigate the kinematics of the high-velocity-dispersion regions, we created PVDs for the S(1) and S(3) emission lines, which are shown in Figs. 6 and 7. We extracted them using a pseudo-slit width of 0.6 (~100 pc). In the first column of Figs. 6 and 7, we show the position and orientation of the regions used to produce the PVDs on top of the flux and velocity dispersion maps of each line. The second and third columns of the top row show the PVDs extracted along the major and minor kinematic axes derived from 3D BAROLO, and those in the bottom correspond to the PVDs extracted along PAs chosen to capture the regions of interest (i.e., of high velocity dispersion). We also included the PVDs generated from the model datacube created by 3D BAROLO: we produced PVDs extracted in slits with the same width and PAs described above and overlaid them as red contours on the S(1) and S(3) PVDs. For S(1) and S(3), we observe a steep velocity gradient that goes from 300 to −250 km s⁻¹ along the kinematic major axis<fnref n="2" />. This velocity gradient is less steep along the minor axis for both lines. The nuclear high-velocity gas that is observed in all the PVDs between −1" and 1" is an artifact of the continuum subtraction. In the left and central bottom panels of Fig. 7, we observe a high-velocity-dispersion blob about ~2.5" westward of the nucleus, which we labeled R5. This blob is located in a “forbidden” quadrant of the rotation pattern and indicates a noncircular gas component that reaches velocities of ~350 km s⁻¹. On the opposite side of the nucleus, at ~1.5 to the east, there is another blob of high velocity dispersion (labeled as R1 in Fig. 7) that can be

<footnote n="2">² We note that PVDs show observed velocities, while the velocity maps shown in Fig. 3 correspond to average velocities.</footnote>
<figure><img src="image_7.png" /><figcaption>Fig. 7. Same as Fig. 6, but for S(3). The black arrows show two high velocity dispersion regions (R1 and R5 in Fig. 4).</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>Fig. 8. Same as Fig. 6, but for CO(2–1). The regions used to extract the PVDs have a size of 0.6 × 7" and PAs = 56° and 146°. The contours correspond to 3, 6, 9, 12, 15, 32, and 64 σ (= 2 × 10⁻<fnref n="4" />).</figcaption></figure>

identified in the same PVD, and that reaches similar velocities. Apart from these two high-velocity and high-velocity-dispersion features, the PVDs of both S(1) and S(3) show extended gas at lower velocities that is outside the model contours, for example in the S(1) and S(3) PVDs extracted along the minor axis and along PA = 3° (see right panels in Figs. 6 and 7). This emission corresponds to gas in the high-velocity-dispersion regions at the north and south edges of the S(1) map.

The PVDs of the CO(2-1) emission are displayed in Fig. 8, extracted along the major and minor kinematic axes derived from ³D BAROLO. The width and orientation are similar to those used to extract the PVDs of S(1). We observe that most of the emission along the major axis can be accounted for by circular motions<fnref n="3" />. However, there is nuclear emission showing high blueshifted velocities (~-200 km⁻¹) to the southwest outside the model contours. Along the minor axis (PA = 146°), we identify noncircular motions in two regions, one to the northwest (~-100 km⁻¹) and the other to the southeast, near the nucleus (~-200 km⁻¹). These regions correspond to noncircular motions, which we discuss in Sect. 4.

### 3.3. Properties of the warm molecular gas

We created rotational diagrams following the methodology described in Rigopoulou et al. (2002) and Pereira-Santaella et al. (2014). For this, we used the line fluxes derived for the six regions defined in Sect. 3.2.2 (see Fig. 4), corrected for extinction by adopting the value of $A_V$ = 1 mag from Prieto et al. (2014) and the Cardelli et al. (1989) reddening law<fnref n="4" />. Rotational diagrams compare the column density in the upper level of each transition ($N_u$), normalized by its statistical weight, $g(u)$, as a function of the temperature, $T(u)$, associated with the upper-level energy, $E_u$, and they are derived from the Boltzmann equation:

$$ \frac{N_u}{N} = \frac{g(u)}{Z(T_{ex})} \times e^{-\frac{T_u}{T_{ex}}}, \quad (1) $$

where $T_{ex}$ is the excitation temperature, $Z(T_{ex})$ is the partition function at $T_{ex}$, and $N_u$ is calculated as

$$ N_u = \frac{\text{flux}(u)}{A(u) \times h \times \nu(u)} \times \frac{4\pi}{\Omega_{\text{beam}}}, \quad (2) $$

<footnote n="3">³ Note that PA = 56° is similar to the PA chosen with BAROLO (PA = 59°).</footnote>

<footnote n="4">⁴ This extinction correction has negligible impact on the line fluxes except for S(3), which is close to the silicate absorption feature.</footnote>
<figure><img src="image_9.png" /><figcaption>**Fig. 9.** H₂ rotational diagrams of the six regions indicated in Fig. 4. The dashed yellow and solid blue lines are the fits using a single PL and a 2T model, respectively. The dash-dotted green and dotted magenta lines correspond to the warm and hot components of the 2T model, respectively.</figcaption></figure>

where flux(*u*) is the flux of a line in the *u*th state, A(*u*) is the Einstein A-coefficient of that transition, *v*(*u*) is the frequency of the transition, *h* is Planck's constant, and Ω<sub>beam</sub> is the area of the region used to extract the flux, calculated as Ω<sub>beam</sub> = πr<sup>2</sup>, where *r* = 0.5". We derived the warm (*T* < 500 K) and hot (*T* > 500 K) molecular gas masses of the different regions from these excitation diagrams, which use the S(1) to S(5) transitions. The critical densities of these transitions are relatively low, with n<sub>H2</sub> (crit) = 10<sup>2</sup>–10<sup>5</sup> cm<sup>-3</sup> at 500 K (Le Bourlot et al. 1999). Therefore, we can assume local thermodynamic equilibrium conditions (Roussel et al. 2007).

We fit the rotational diagram of the six regions using two
models: (1) a single power law (PL) temperature distribution<fnref n="5" />
(solid red line) and (2) a two-temperature (2T) model, which
represents a warm and a hot component (dashed and dotted
green lines in Fig. 9). In the case of region R5, we just fit
one temperature component because there is no S(5) measure-
ment there. At first glance, we notice that the excitation tem-
perature increases as the energy level increases, indicating that
the gas consists of several components with different tempera-
tures (Rigopoulou et al. 2002). In Table 3, we report the χ²/d.o.f.
(degrees of freedom), temperatures, column densities, and gas
masses for each region. These values were obtained by fitting
the S(1) to S(5) transitions using either the single PL or the 2T
model. In the 2T model, the warm component dominates the S(1)
and S(2) emission, and the hotter component dominates the S(3)
to S(5) emission. Generally, the models fit the data reasonably
well (χ²/d.o.f. < 1.9).

Excluding region R5, the average (± the standard devia-
tion of the values reported in Table 3) temperatures for the
warm and hot components of the 2T model are T<sub>w</sub> = 278 ± 38 K
and T<sub>h</sub> = 789 ± 90 K, respectively. The excitation tempera-
tures of the hot component are consistent with those reported
in previous works and are typical temperatures in AGNs
(e.g., Davies et al. 2005). The difference between T<sub>w</sub> and
T<sub>h</sub> is ~500 K, with the warm component being the main

contributor to column density and mass. The average col-
umn density and mass values of all the regions, excluding
R5, are log(N<sub>H</sub><sup>w</sup>/cm<sup>-2</sup>) = 20.3 ± 0.2 and (9.6 ± 5.7) × 10<sup>4</sup> M<sub>☉</sub>
for the warm component, and log(N<sub>H</sub><sup>h</sup>/cm<sup>-2</sup>) = 18.6 ± 0.4 and
(1.5 ± 2.9) × 10<sup>4</sup> M<sub>☉</sub> for the hot component. In regions R1, R2,
R3, and R6, the mass of the hot component is less than 3% of
the mass of the warm component. In region R4, we find the
largest contribution from the hot component, of 1.4 × 10<sup>4</sup> M<sub>☉</sub>
(see Table 3), likely due to its position, closer to the galaxy
nucleus than other regions (see Fig. 4). In the case of region R5,
the knot of high velocity dispersion in S(3), we measure a tem-
perature of 450 K, a column density of log(N<sub>H</sub>/cm<sup>-2</sup>) = 19.0,
and a gas mass of 4 × 10<sup>3</sup> M<sub>☉</sub>, which represents just 4% of the
average warm gas mass obtained for the other regions. This high
temperature comes from the fit using only one temperature com-
ponent, as a S(5) measurement is lacking in this region, resulting
in a high uncertainty for this measurement.

Additionally, in the last two rows of Table 3, labeled as Total,
we report the values obtained using the H₂ fluxes measured in an
aperture of 3′′6 (~630 pc) in diameter. We found that the differ-
ence between the hot and the warm temperatures is ~550 K. The
mass of the warm gas is ~ 6.4 × 10⁵ M☉, which falls within the
range estimated by Roussel et al. (2007) for a sample of low-
ionization nuclear emission-line region (LINER) and Seyfert
galaxies, and the mass of hot gas is ~2.7×10⁴ M☉. We also calcu-
lated the mass of cold molecular gas using the CO(2–1) flux inte-
grated in the same 3′′6 diameter aperture (see dashed white circle
in Fig. 1). We used Eq. (3) from Solomon & Vanden Bout (2005)
to compute the CO(2–1) line luminosity, L'CO, in K km s⁻¹ pc²
units. Then, we assumed that the CO emission is thermal-
ized and optically thick (R<sub>21</sub> = L'CO(2–1)/L'CO(1–0) = 1), and
adopted the Milky-Way CO-to-H₂ conversion factor (αCO =
4.36 ± 1.30 (K km s⁻¹ pc²)⁻¹; Bolatto et al. 2013) to calculate
a mass of M<sub>H2</sub> = (1.36 ± 0.50) × 10⁷ M☉. The warm gas mass
measured from the rotational H₂ lines represents less than
5% of the total molecular gas estimated from CO(2–1). The
major source of uncertainty for this mass gas value is the
assumed conversion factor, which could be up to a factor of

<footnote n="5">⁵ We assumed upper and lower temperatures of 3500 K and 200 K.</footnote>
<table><caption>**Table 3.** Parameters derived from the PL and 2T fits shown in Fig. 9.</caption><thead><tr><th>Region name</th><th>Fit model</th><th>χ<sup>2</sup>/d.o.f.</th><th>β</th><th>T<sub>w</sub> (°K)</th><th>log(N<sub>H</sub><sup>w</sup>) (cm<sup>-2</sup>)</th><th>M<sub>w</sub> 10<sup>5</sup> M<sub>☉</sub></th><th>T<sub>h</sub> (°K)</th><th>log(N<sub>H</sub><sup>h</sup>) (cm<sup>-2</sup>)</th><th>M<sub>h</sub> 10<sup>5</sup> M<sub>☉</sub></th></tr></thead><tbody><tr><td rowspan="2">R1</td><td>PL</td><td>1.28</td><td>5.03</td><td>—</td><td>20.1</td><td>0.55</td><td>—</td><td>—</td><td>—</td></tr><tr><td>2T</td><td>1.84</td><td>—</td><td>271 ± 54</td><td>20.1</td><td>0.50</td><td>770 ± 171</td><td>18.5</td><td>0.01</td></tr><tr><td rowspan="2">R2</td><td>PL</td><td>0.36</td><td>5.23</td><td>—</td><td>20.4</td><td>1.05</td><td>—</td><td>—</td><td>—</td></tr><tr><td>2T</td><td>0.24</td><td>—</td><td>266 ± 16</td><td>20.4</td><td>1.0</td><td>760 ± 46</td><td>18.7</td><td>0.02</td></tr><tr><td rowspan="2">R3</td><td>PL</td><td>1.89</td><td>5.42</td><td>—</td><td>20.3</td><td>0.77</td><td>—</td><td>—</td><td>—</td></tr><tr><td>2T</td><td>0.22</td><td>—</td><td>224 ± 13</td><td>20.5</td><td>1.41</td><td>705 ± 29</td><td>18.7</td><td>0.02</td></tr><tr><td rowspan="2">R4</td><td>PL</td><td>0.41</td><td>4.11</td><td>—</td><td>20.7</td><td>2.16</td><td>—</td><td>—</td><td>—</td></tr><tr><td>2T</td><td>0.42</td><td>—</td><td>308 ± 34</td><td>20.6</td><td>1.62</td><td>943 ± 100</td><td>19.0</td><td>0.67</td></tr><tr><td rowspan="2">R5</td><td>PL</td><td>0.34</td><td>4.78</td><td>—</td><td>19.5</td><td>0.12</td><td>—</td><td>—</td><td>—</td></tr><tr><td>1T</td><td>1.05</td><td>—</td><td>450 ± 55</td><td>19.0</td><td>0.04</td><td>—</td><td>—</td><td>—</td></tr><tr><td rowspan="2">R6</td><td>PL</td><td>0.44</td><td>5.18</td><td>—</td><td>20.2</td><td>0.58</td><td>—</td><td>—</td><td>—</td></tr><tr><td>2T</td><td>0.86</td><td>—</td><td>320 ± 43</td><td>20.0</td><td>0.30</td><td>767 ± 144</td><td>18.0</td><td>0.01</td></tr><tr><td rowspan="2">Total</td><td>PL</td><td>0.64</td><td>4.28</td><td>—</td><td>20.1</td><td>6.98</td><td>—</td><td>—</td><td>—</td></tr><tr><td>2T</td><td>0.48</td><td>—</td><td>277 ± 35</td><td>20.1</td><td>6.38</td><td>834 ± 71</td><td>18.7</td><td>0.27</td></tr></tbody></table>

**Notes.** The results from the fit with the single PL temperature distribution, which represents both the hot and warm gas and assumes that the column density follows a PL ($dN \sim T^{\beta}dT$), are shown in the same columns as those corresponding to the warm component of the two temperatures (2T) fits.

<figure><img src="image_10.png" /><figcaption>**Fig. 10.** H<sub>2</sub> excitation map derived from the S(1)/S(3) line ratio (see Sect. 3.3). The dashed black circle marks the region R5, also indicated in Fig. 4 and with an arrow in Fig. 7. This region has the highest temperature area on this map, with a temperature of ~260 K. The solid black contours correspond to 0.2, 0.4, 0.6, 0.8, and 1σ of the CO(2-1) emission. The white star shows the AGN location according to the ALMA continuum.</figcaption></figure>

3–4 lower<sup>6</sup>. However, our warm-to-cold gas mass ratio is among the highest reported for nearby ULIRGs by Higdon et al. (2006) using data from *Spitzer*, which are typically less than 1%. Furthermore, our fraction of 5% is within the range of values obtained by Pereira-Santaella et al. (2014) for a sample of six local infrared bright Seyfert galaxies also observed with *Spitzer*, of 5–50%. Recently, Pereira-Santaella et al. (2024) calculated the α<sub>CO</sub> value for NGC 3256 using various methods with JWST/NIRSpec observations. They derived α<sub>CO</sub> values for different regions of this galaxy, finding that these values range from 0.26–1.09 (K km s<sup>-1</sup> pc<sup>2</sup>)<sup>-1</sup>. The median α<sub>CO</sub> values were 0.40

and 0.61 (K km s<sup>-1</sup> pc<sup>2</sup>)<sup>-1</sup> using one- and two-temperature models, respectively. These α<sub>CO</sub> values are 1.3–2 times lower than typically assumed for ULIRGs. Using these median α<sub>CO</sub> values, we calculated a warm-to-cold gas mass ratio of 0.51 and 0.33.

Finally, in Fig. 10 we show an excitation map of MCG-05-16-23 using the S(1)/S(3) line ratio in the FoV of Ch2-m, determined by the S(3) line (see Table 1) and following the same procedure as in Davies et al. (2024). The nuclear spiral and the ring show higher temperatures than the inter-arm gas, and the S(3) high-velocity-dispersion knot (R5 in Fig. 4) shows the highest temperature of ~260 K determined from the S(1)/S(3) line ratio. In the next section, we discuss the possible origin of this region.

## 4. Discussion

As can be seen from Fig. 1, MCG-05-23-16 shows a dust lane in the optical HST/WFC3 image that resembles a nuclear-trailing spiral and a larger-scale ring (see also Fig. A.1). The dust lane, of ~2" extent (~360 pc) and oriented along the major axis of the galaxy, was previously reported by Prieto et al. (2014). The CO(2-1) and H<sub>2</sub> moment maps reveal that the nuclear spiral is connected with the ring, of major axis ~8" (1.4 kpc), through a network of arms and filaments (see Figs. B.1 and 3). A nuclear trailing spiral connected with a ring is indicative of gas that is losing angular momentum and feeding the SMBH in the presence of a bar (e.g., Regan & Teuben 2003; Buta 2013; Shimizu et al. 2019; Mondal & Chattopadhyay 2021). According to several studies, circumnuclear rings are formed near the inner Lindblad resonances (ILR, Lindblad 1964) due to the bar potential when the gas loses angular momentum and accumulates in that region (see references in Schmidt et al. 2019). Recently, Sormani et al. (2024) proposed that a nuclear ring is actually the inner edge of an extensive gap that opens around the ILR due to the excitation of waves by a bar potential.

Inside the nuclear ring at the ILR, the gravity torques due to the presence of a nuclear bar are negative and induce the formation of a trailing nuclear spiral (e.g., Kim & Elmegreen 2017), favoring the transport of gas toward the center, as is seen in other observations of nearby Seyfert galaxies (see Combes et al.

<footnote n="6"><sup>6</sup> If instead of using the Milky-Way CO-to-H<sub>2</sub> α<sub>CO</sub> factor we assume the typical value of ULIRGs of 0.8 (K km s<sup>-1</sup> pc<sup>2</sup>)<sup>-1</sup> (Narayanan et al. 2012), we measure M<sub>H<sub>2</sub></sub> = (2.49 ± 0.5) × 10<sup>6</sup> M<sub>☉</sub>. Assuming this mass, we obtain a warm/cold gas mass ratio equal to 0.25, which also agrees with previous works.</footnote>
2014; Audibert et al. 2019, 2021). The winding pattern of the nuclear spiral is predicted due to the presence of a massive black hole in the center (Buta & Combes 1996), and in the case of MCG-05-23-16 it seems to coincide with the trailing spiral arms on a large scale, according to the reported geometry of the galaxy disk (i.e., the northwest side is the far side), and the observed sense of rotation of the gas. The computation of the gravity torques due to the nuclear bar is beyond the scope of this paper. However, the presence of a nuclear bar in MCG-05-23-16 is further supported by the S-shape zero isovelocities observed in the CO(2-1) velocity map shown in Fig. B.1, which are not aligned with the minor axis (see e.g., Sanders & Tubbs 1980; Pence & Blackman 1984; Randriamampandry et al. 2015; Cohen et al. 2020, for references on S-shape velocities as bar tracers).

On the largest scales that can be probed with S(1), which coincide with the ring, we see evidence of strong noncircular motions along the minor axis both in CO and H₂. Considering that the southeast is the near side of the galaxy, based on the higher extinction there (see right panel of Fig. 1 and Prieto et al. 2014) the velocity residuals along the kinematic minor axis (right panels of Fig. B.2) would correspond to the radial projection of local inflowing gas motions in the plane of the galaxy (positive velocities in the southeast and negative velocities in the northwest). This is also seen in the PVD diagram of S(1) along the minor axis (see top right panel of Fig. 6): there is gas within the ring that is not reproduced by the rotating disk model fit with ³D BAROLO that is blueshifted to the northwest and redshifted to the southeast. Something similar is seen in S(3), but to a lesser extent due to the smaller FoV of Ch2. These residuals along the minor axis are relatively far from the nuclear spiral, and they can be better explained by the radial projection of noncircular motions typical of the elliptical gas orbits in the presence of a bar, rather than by being the signature of a truly radial inflow.

Evidence of noncircular motions that could be consistent with outflowing gas, although not aligned with the kinematic minor axis, are found from the S(3) kinematics, and to a lesser extent, from S(1). The velocity maps of these emission lines show red residuals at ~2.5" (440 pc) westward of the nucleus (region R5 in Fig. 4, at PA ~270°; see last column in Fig. B.2) and in the case of S(3), also velocity dispersion residuals (see Fig. B.3). These residuals are consistent, in principle, with outflowing gas, and they can be seen in the PVD along PA = 87° shown in Fig. 7, clearly deviating from the rotation disk model fit with ³D BAROLO.

This area of higher velocity dispersion in S(3) coincides with an area devoid of CO(2-1) emission (see Figs. 3 and 4). This suggests that the more turbulent, hotter, and excited gas detected with S(3) fills regions depleted of cold gas. This is observed when we compare the morphology of CO and those of the different H₂ lines observed by MIRI (see Fig. 3): the CO traces the nuclear spiral, the ring, and the connecting arms, and the regions that are devoid of CO are occupied by warmer molecular gas showing the highest values of the velocity dispersion there. Considering the compact and localized morphology of the R5 knot seen in S(3), it is possible that that gas might have been accelerated or disrupted by star formation there. This R5 knot cannot be explained by noncircular motions in the plane of the galaxy because they are well outside the expected virial range, and they show an opposite sign of the velocities (see middle panel of Fig. 7). A jet-driven outflow could be another possible scenario for explaining the kinematics, but we do not see any spatial correspondence with the jet (black contours in Fig. 4). A localized, star-formation-driven outflow is supported by the

<figure><img src="image_11.png" /><figcaption>**Fig. 11.** Comparison of S(3) velocity dispersion and 11.3 µm PAH feature. The left panel shows the S(3) velocity dispersion map with contours of the 11.3 µm PAH feature overlaid in black. The right panel shows the profiles of the 11.3 µm PAH feature extracted from the regions indicated in the left panel.</figcaption></figure>

presence of polycyclic aromatic hydrocarbon (PAH) emission in the area between the nucleus and the S(3) knot, traced by the 11.3 µm feature. This PAH band has been used as a reliable tracer of recent star formation in AGNs (e.g., Diamond-Stanic & Rieke 2010, 2012; Esparza-Arredondo et al. 2018), which is also found to be the most resilient PAH band in the harsh environment of AGNs and their outflows (García-Bernete et al. 2022, 2024b).

Using the integrated luminosity of the PAH feature extracted from the region indicated with a dashed blue rectangle in Fig. 11 ($L_{11.3\mu m} = 7.34 \times 10^{39} \text{ erg s}^{-1}$)<fnref n="7" /> and Eq. (12) in Shipley et al. (2016), as in Ramos Almeida et al. (2023), we measured an SFR of $0.013_{-0.004}^{+0.005} M_{\odot} yr^{-1}$. Gas compression induced by the inflowing gas might have triggered star formation locally, and the young stars producing the PAH emission there could have driven a localized warm molecular outflow in the region shown with a dashed black rectangle in Fig. 11. The mass of the H₂ knot (R5) calculated in Sect. 3.3 is 4000 $M_{\odot}$ (see Table 3, 1T model), but part of the gas there is low-velocity gas participating in rotation (negative velocities at ~2.5" to the west in the central bottom panel of Fig. 7) and part is outflowing (positive velocities at the same position). Therefore, the outflow mass will be lower than the mass estimated for R5 as a whole, so we consider the latter as an upper limit to the outflow mass. By then using a velocity representative of the outflowing gas from the PVD shown in Fig. 7, of ~200 km s⁻¹, and a radius of ~0.75 (88 pc, which is the R5 radius), we derive a mass outflow rate of $\dot{M}_{\text{out}} = M_{\text{H2}}^{R5} \times v_{\text{out}} \times r_{\text{out}}^{-1} < 0.01 M_{\odot} \text{ yr}^{-1}$ and kinetic energy of <1.2 × 10³⁸ erg s⁻¹ (see Eq. (3) in Speranza et al. 2024). According to the estimated values of the SFR and mass outflow rate, star formation might be enough to drive the outflow (mass loading factor = $\dot{M}_{\text{out}}$/SFR < 0.8), although the uncertainties are significant.

As was mentioned before, despite the presence of extended radio emission that is consistent with a compact jet of ~200 pc in size and PA ~−10° (Orienti & Prieto 2010), we do not find evidence of a strong jet-ISM interaction in MCG-05-23-16, at least in the case of the molecular gas. If the radio emission corresponds to a jet, we would be seeing the approaching side, considering its morphology: elongated and brighter toward the north. In that case, and considering that the northwest is the far side of the galaxy, the jet would be subtending a large angle relative to the CO and H₂ disks, and therefore it would have a small impact on the ISM (Meenakshi et al. 2022; Ramos Almeida et al. 2022; Audibert et al. 2023).

<footnote n="7">⁷ $L_{11.3\mu m}$ was multiplied by a factor of two to make the star formation rate (SFR) comparable with PAHFIT-derived measurements.</footnote>
## 5. Conclusions

In this work, we characterize for the first time the physical properties (temperature, column density, and mass), morphology, and kinematics of the molecular gas in the spheroidal galaxy MCG-05-23-16, using data from the MRS module of JWST/MIRI and ALMA. Thanks to these datasets, we studied the connection between the warm ($T \sim 100-1000$ K) and cold ($T < 100$ K) molecular gas traced by the rotational mid-infrared $H_2$ lines and the sub-millimeter CO(2-1) transition.

Our key findings include the detection of a nuclear trailing spiral ~350 pc in size and a surrounding ring of ~1.8 kpc in both CO(2-1) and $H_2$, which are also identified through optical HST observations. These features indicate the typical gas response to a stellar bar, where the ring may correspond to the ILR of the bar and the trailing nuclear spiral inside the ring is triggering an ongoing fueling episode of the AGN.

The morphology and kinematics of the gas traced with the $H_2$-0S(1), S(2), and CO(2-1) lines are similar, with rotational velocities of up to ±250 km s⁻¹ and noncircular motions on the spatial scales of the ring that are consistent with the projection of noncircular motions associated with the elliptical gas orbits followed by the gas. On the other hand, the warmer and more excited molecular gas traced with the S(3), S(4), and S(5) lines shows more complex kinematics that include clumps of high velocity dispersion (of up to ~160 km s⁻¹) in regions where CO(2-1) is undetected. One of these clumps, located ~350 pc westward of the nucleus, shows kinematics that is consistent with outflowing gas, which could be driven by recent star formation in adjacent regions, traced by the 11.3 µm PAH feature. The compact jet previously detected through VLA observations does not appear to impact the molecular gas, probably because its orientation is not favorable for it to entrain the molecular gas disk.

Overall, we observe a stratification of the molecular gas. The colder gas is located in the nuclear spiral, ring, and connecting arms, while most warmer gas with a higher velocity dispersion fills the space between them. This work highlights the potential synergy between JWST and ALMA observations in advancing the study of molecular gas in galaxies.

**Acknowledgements.** We thank the referee for their useful comments and suggestions. D.E.A., C.R.A., and A.A. acknowledge support from the Agencia Estatal de Investigación of the Ministerio de Ciencia, Innovación y Universidades (MCIU/AEI) under the grant "Tracking active galactic nuclei feedback from parsec to kiloparsec scales", with reference PID2022-141105NB-I00 and the European Regional Development Fund (ERDF). D.E.A. and B.G.L. acknowledge support from the Spanish Ministry of Science and Innovation through the Spanish State Research Agency (AEI-MCINN/10.13039/501100011033) through grants "Participation of the Instituto de Astrofísica de Canarias in the development of HARMONI: D1 and Delta-D1 phases" with references PID2019-107010GB100 and PID2022-140483NB-C21 and the Severo Ochoa Program 2020-2023 (CEX2019-000920-S)". M.P.S. acknowledges support from grants RYC2021-033094-I and CNS2023-145506 funded by MCIN/AEI/10.13039/501100011033 and the European Union NextGenerationEU/PRTR. L.H.M. and A.A.H. 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 ERDF A way of making Europe. A.J.B. acknowledges funding from the "FirstGalaxies" Advanced Grant from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (Grant agreement No. 789056). E.K.S.H., L.Zhang, C.Packhman acknowledges grant support from the Space Telescope Science Institute (ID: JWST-GO-01670.007-A). S.G.B. acknowledges support from the Spanish grant PID2022-138560NB-I00, funded by MCIN/AEI/10.13039/501100011033/FEDER, EU. CR acknowledges support from Fondecyt Regular grant 1230345, ANID BASAL project FB210003, and the China-Chile joint research fund. I.G.B. is supported by the Programa Atención de Talento Investigador "César Nombela" via grant 2023-T1/TEC-29030 funded by the Community of Madrid. E.B. acknowledges support from the

Spanish grants PID2022-138621NB-I00 and PID2021-123417OB-I00, funded by MCIN/AEI/10.13039/501100011033/FEDER, EU. T.D.S. 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: 03382). 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 Telescope Science Institute, operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127 for JWST. The specific observations analyzed here can be accessed via doi:10.17909/vre3-m991. This paper makes use of the following ALMA data: ADS/JAO.ALMA# 2019.1.01742.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC(Canada) and NSC and ASIAA (Taiwan), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under a cooperative agreement by Associated Universities, Inc. We used observations made with the NASA/ESA Hubble Space Telescope, and obtained from the Hubble Legacy Archive, which is a collaboration between the Space Telescope Science Institute (STScI/NASA), the Space Telescope European Coordinating Facility (ST-ECF/ESA), and the Canadian Astronomy Data center (CADC/NRC/CSA).

## References

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

* Á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

* Asmus, D., Gandhi, P., Hönig, S. F., et al. 2015, MNRAS, 454, 766

* Audibert, A., Combes, F., García-Burillo, S., et al. 2019, A&A, 632, A33

* Audibert, A., Combes, F., García-Burillo, S., et al. 2021, A&A, 656, A60

* Audibert, A., Ramos Almeida, C., García-Burillo, S., et al. 2023, A&A, 671, L12

* Bianchin, M., U, V., Song, Y., et al. 2024, ApJ, 965, 103

* Bolatto, A. D., Wolfire, M., & Leroy, A. K. 2013, ARA&A, 51, 207

* Briggs, D. S. 1995, Ph.D. Thesis, New Mexico Institute of Mining and Technology

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

* Buta, R. J. 2013, in *Secular Evolution of Galaxies*, Lecture 2: Barred and spiral galaxies, 155

* Buta, R., & Combes, F. 1996, Fund. Cosmic Phys., 17, 95

* Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245

* Cazaux, S., & Tielens, A. G. G. M. 2002, ApJ, 575, L29

* Cohen, D. P., Turner, J. L., & Consiglio, S. M. 2020, MNRAS, 493, 627

* Combes, F. 2003, *Active Galactic Nuclei: From Central Engine to Host Galaxy*, 290, 411

* Combes, F., García-Burillo, S., Casasola, V., et al. 2014, A&A, 565, A97

* Dale, D. A., & Helou, G. 2002, ApJ, 576, 159

* Dale, D. A., Sheth, K., Helou, G., et al. 2005, AJ, 129, 2197

* Davies, R. I., Tacconi, L. J., & Genzel, R. 2004, ApJ, 602, 148

* Davies, R. I., Sternberg, A., Lehnert, M. D., et al. 2005, ApJ, 633, 105

* Davies, R. I., Thomas, J., Genzel, R., et al. 2006, ApJ, 646, 754

* Davies, R. I., Maciejewski, W., Hicks, E. K. S., et al. 2009, ApJ, 702, 114

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

* Davies, R., Shimizu, T., Pereira-Santaella, M., et al. 2024, A&A, 689, A263

* Diamond-Stanic, A. M., & Rieke, G. H. 2010, ApJ, 724, 140

* Diamond-Stanic, A. M., & Rieke, G. H. 2012, ApJ, 746, 168

* Di Teodoro, E. M., & Fraternali, F. 2015, MNRAS, 451, 3021

* Domínguez-Fernández, A. J., Alonso-Herrero, A., García-Burillo, S., et al. 2020, A&A, 643, A127

* Dubois, Y., Peirani, S., Pichon, C., et al. 2016, MNRAS, 463, 3948

* Esparza-Arredondo, D., González-Martín, O., Dultzin, D., et al. 2018, ApJ, 859, 124

* Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9

* Ferruit, P., Wilson, A. S., & Mulchaey, J. 2000, ApJS, 128, 139

* 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, A&A, 691, A162

* García-Burillo, S., & Combes, F. 2012, J. Phys. Conf. Ser., 372, 012050

* García-Burillo, S., Combes, F., Schinnerer, E., et al. 2005, A&A, 441, 1011

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

* García-Burillo, S., Hicks, E. K. S., Alonso-Herrero, A., et al. 2024, A&A, 689, A347
Gardner, J. P., Mather, J. C., Abbott, R., et al. 2023, PASP, 135, 068001

Goodrich, R. W., Veilleux, S., & Hill, G. J. 1994, ApJ, 422, 521

Habart, E., Walmsley, M., Verstraete, L., et al. 2005, Space Sci. Rev., 119, 71

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

Hermosa Muñoz, L., Alonso-Herrero, A., Pereira-Santaella, M., et al. 2024, A&A, 690, A350

Herrera-Camus, R., Janssen, A., Sturm, E., et al. 2020, A&A, 635, A47

Hickox, R. C., Mullaney, J. R., Alexander, D. M., et al. 2014, ApJ, 782, 9

Hicks, E. K. S., Davies, R. I., Malkan, M. A., et al. 2009, ApJ, 696, 448

Higdon, S. J. U., Armus, L., Higdon, J. L., et al. 2006, ApJ, 648, 323

Ineson, J., Croston, J. H., Hardcastle, M. J., et al. 2015, MNRAS, 453, 2682

Kawamuro, T., Ueda, Y., Tazaki, F., et al. 2016, ApJS, 225, 14

Kim, W.-T., & Elmegreen, B. G. 2017, ApJ, 841, L4

Kolcu, T., Maciejewski, W., Gadotti, D. A., et al. 2023, MNRAS, 524, 207

Kormendy, J., & Richstone, D. 1995, ARA&A, 33, 581

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

Labiano, A., Argyriou, I., Álvarez-Márquez, J., et al. 2021, A&A, 656, A57

Le Bourlot, J., Pineau des Forêts, G., & Flower, D. R. 1999, MNRAS, 305, 802

Lebouteiller, V., Barry, D. J., Spoon, H. W. W., et al. 2011, ApJS, 196, 8

Lindblad, B. 1964, Astrophysica Norvegica, 9, 103

Liu, V., Zoghbi, A., & Miller, J. M. 2024, ApJ, 963, 38

Magorrian, J., Tremaine, S., Richstone, D., et al. 1998, AJ, 115, 2285

Marinucci, A., Muleri, F., Dovciak, M., et al. 2022, MNRAS, 516, 5907

Martini, P. 2004, The Interplay Among Black Holes. Stars and ISM in Galactic Nuclei, 222, 235

McMullin, J. P., Waters, B., Schiebel, D., et al. 2007, Astronomical Data Analysis Software and Systems XVI, 376, 127

Meenakshi, M., Mukherjee, D., Wagner, A. Y., et al. 2022, MNRAS, 516, 766

Mondal, D., & Chattopadhyay, T. 2021, Celest. Mech. Dyn. Astron., 133, 43

Mundell, C. G., Ferruit, P., Nagar, N., et al. 2009, ApJ, 703, 802

Narayanan, D., Krumholz, M. R., Ostriker, E. C., et al. 2012, MNRAS, 421, 3127

Novak, G. S., Ostriker, J. P., & Ciotti, L. 2011, ApJ, 737, 26

Olivares, V., Su, Y., Nulsen, P., et al. 2022, MNRAS, 516, L101

Orienti, M., & Prieto, M. A. 2010, MNRAS, 401, 2599

Pence, W. D., & Blackman, C. P. 1984, MNRAS, 210, 547

Peralta de Arriba, L., Alonso-Herrero, A., García-Burillo, S., et al. 2023, A&A, 675, A58

Pereira-Santaella, M., Spinoglio, L., Busquet, G., et al. 2013, ApJ, 768, 55

Pereira-Santaella, M., Spinoglio, L., van der Werf, P. P., et al. 2014, A&A, 566, A49

Pereira-Santaella, M., Colina, L., García-Burillo, S., et al. 2018, A&A, 616, A171

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

Pereira-Santaella, M., González-Alfonso, E., García-Bernete, I., et al. 2024, A&A, 681, A117

Prieto, M. A., Mezcua, M., Fernández-Ontiveros, J. A., et al. 2014, MNRAS, 442, 2145

Ramos Almeida, C., Pérez García, A. M., & Acosta-Pulido, J. A. 2009, ApJ, 694, 1379

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

Ramos Almeida, C., Esparza-Arredondo, D., González-Martín, O., et al. 2023, A&A, 669, L5

Randriamampandry, T. H., Combes, F., Carignan, C., et al. 2015, MNRAS, 454, 3743

Reeves, J. N., Awaki, H., Dewangan, G. C., et al. 2007, PASJ, 59, 301

Regan, M. W., & Teuben, P. 2003, ApJ, 582, 723

Rieke, G. H., Ressler, M. E., Morrison, J. E., et al. 2015, PASP, 127, 665

Rigopoulou, D., Kunze, D., Lutz, D., et al. 2002, A&A, 389, 374

Rodríguez-Ardila, A., Riffel, R., & Pastoriza, M. G. 2005, MNRAS, 364, 1041

Roussel, H., Helou, G., Hollenbach, D. J., et al. 2007, ApJ, 669, 959

Runco, J. N., Malkan, M. A., Fernández-Ontiveros, J. A., et al. 2020, ApJ, 905, 57

Ruschel-Dutra, D., Storchi-Bergmann, T., Schnorr-Müller, A., et al. 2021, MNRAS, 507, 74

Sanders, R. H., & Tubbs, A. D. 1980, ApJ, 235, 803

Schmidt, E. O., Mast, D., Díaz, R. J., et al. 2019, AJ, 158, 60

Shimizu, T. T., Davies, R. I., Lutz, D., et al. 2019, MNRAS, 490, 5860

Shipley, H. V., Papovich, C., Rieke, G. H., et al. 2016, ApJ, 818, 60

Solomon, P. M., & Vanden Bout, P. A. 2005, ARA&A, 43, 677

Sormani, M. C., Sobacchi, E., & Sanders, J. L. 2024, MNRAS, 528, 5742

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

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

Sturm, E., González-Alfonso, E., Veilleux, S., et al. 2011, ApJ, 733, L16

Tremaine, S., Gebhardt, K., Bender, R., et al. 2002, ApJ, 574, 740

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

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

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

Zhang, L., & Ho, L. C. 2023, ApJ, 953, L9

Zhang, L., García-Bernete, I., Packham, C., et al. 2024a, ApJ, 975, L2

Zhang, L., Packham, C., Hicks, E. K. S., et al. 2024b, ApJ, 974, 195

Zoghbi, A., Cackett, E. M., Reynolds, C., et al. 2014, ApJ, 789, 56

<footnote n="1">¹ Instituto de Astrofísica de Canarias, Calle Vía Láctea, s/n, E-38205 La Laguna, Tenerife, Spain</footnote>

<footnote n="2">² Departamento de Astrofísica, Universidad de La Laguna, E-38206 La Laguna, Tenerife, Spain</footnote>

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

<footnote n="4">⁴ Centro de Astrobiología (CAB), CSIC-INTA, Camino Bajo del Castillo s/n, 28692 Villanueva de la Cañada, Madrid, Spain</footnote>

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

<footnote n="6">⁶ Max Planck Institute for Extraterrestrial Physics (MPE), Giessen-bachstr.1, 85748 Garching, Germany</footnote>

<footnote n="7">⁷ LERMA, Observatoire de Paris, Collège de France, PSL University, 837 CNRS, Sorbonne University, Paris, France</footnote>

<footnote n="8">⁸ Department of Physics and Astronomy, The University of Texas at San Antonio, 1 UTSA Circle, San Antonio, Texas 78249, USA</footnote>

<footnote n="9">⁹ School of Mathematics, Statistics and Physics, Newcastle University, Newcastle upon Tyne NE1 7RU, UK</footnote>

<footnote n="10">¹⁰ Departmento de Física de la Tierra y Astrofísica, Fac. de CC Físicas, Universidad Complutense de Madrid, E-28040 Madrid, Spain</footnote>

<footnote n="11">¹¹ Instituto de Física de Partículas y del Cosmos IPARCOS, Fac. CC Físicas, Universidad Complutense de Madrid, E-28040 Madrid, Spain</footnote>

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

<footnote n="13">¹³ Institute of Astrophysics, Foundation for Research and Technology Hellas (FORTH), Heraklion 70013, Greece</footnote>

<footnote n="14">¹⁴ School of Sciences, European University Cyprus, Diogenes street, Engomi 1516, Nicosia, Cyprus</footnote>

<footnote n="15">¹⁵ Instituto de Radioastronomía y Astrofísica (IRyA-UNAM), 3-72 (Xangari), 8701 Morelia, Mexico</footnote>

<footnote n="16">¹⁶ Department of Physics & Astronomy, University of Alaska Anchorage, Anchorage, Alaska, USA</footnote>

<footnote n="17">¹⁷ Department of Physics, University of Alaska, Fairbanks, Alaska 99775-5920, USA</footnote>

<footnote n="18">¹⁸ Telespazio UK for the European Space Agency (ESA), ESAC, Camino Bajo del Castillo s/n, 28692 Villanueva de la Cañada, Spain</footnote>

<footnote n="19">¹⁹ Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA</footnote>

<footnote n="20">²⁰ Instituto de Estudios Astrofísicos, Facultad de Ingeniería y Ciencias, Universidad Diego Portales, Av. Ejército Libertador 441, Santiago, Chile</footnote>

<footnote n="21">²¹ Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing 100871, China</footnote>

<footnote n="22">²² School of Physics & Astronomy, University of Southampton, High-843 field, Southampton SO17 1BJ, UK</footnote>

<footnote n="23">²³ Astronomical Observatory, Volgina 7, 11060 Belgrade, Serbia</footnote>

<footnote n="24">²⁴ Sterrenkundig Observatorium, Universiteit Gent, Krijgslaan 281-S9, Gent B-9000, Belgium</footnote>

<footnote n="25">²⁵ Department of Astronomy, School of Science, Graduate University for Advanced Studies (SOKENDAI), Mitaka, Tokyo 181-8588, Japan</footnote>

<footnote n="26">²⁶ Kavli Institute for Particle Astrophysics & Cosmology (KIPAC), Stanford University, Stanford, CA 94305, USA</footnote>

<footnote n="27">²⁷ Department of Physics & Astronomy, University of South Carolina, Columbia, SC 29208, USA</footnote>

<footnote n="28">²⁸ LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Sorbonne Paris Cité, 5 place Jules Janssen, F-92195 Meudon, France</footnote>

<footnote n="29">²⁹ Centre for Extragalactic Astronomy, Durham University, South Road, Durham DH1 3LE, UK</footnote>
Appendix A: Comparison between ALMA and HST maps

<figure><img src="image_12.png" /><figcaption>**Fig. A.1.** V-H color map obtained using HST/WFC3 images in the F606W and F160W filters. The solid cyan contours correspond to $2.2 \times 10^{-3}$, $4.4 \times 10^{-3}$, $6.6 \times 10^{-3}$, $8.8 \times 10^{-3}$, $1.1 \times 10^{-2}$, $1.3 \times 10^{-2}$, $1.5 \times 10^{-2}$, $1.7 \times 10^{-2}$, and $2 \times 10^{-2}$ [Jy km/s] ALMA CO(2-1) intensity. The dashed black and blue circles correspond to the regions used to extract the nuclear (diameter of 0.8'' - 2.5'' ≈ 141-440 pc) and extended (diameter of 3.6'' ≈ 634 pc) JWST/MIRI spectra shown in Fig. 2.</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>**Fig. B.2.** From top to bottom, moment 1 maps of the CO(2–1), H₂O–0S(1), and H₂O–0S(3) emission lines. From left to right: data, ³DBAROLO model, and residuals from subtracting the model from the data. The kinematic major axis determined by ³DBAROLO is indicated as a dashed gray line, and the corresponding beam sizes are shown as blue ellipses in the bottom left corner of the left panels. The black cross corresponds to the peak of the ALMA continuum.</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>**Fig. B.3.** Same as Fig. B.2, but for the moment 2 maps.</figcaption></figure>