Search

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 obtained with the James Webb Space Telescope (JWST) and 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 H2, 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 H2 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 $\sim$ 160 km/s), in regions devoid of CO(2$-$1). The kinematics of one of these clumps, located at $\sim$ 350 pc westward from the nucleus, are consistent with outflowing gas, possibly driven by localized star formation traced by Polycyclic Aromatic Hydrocarbon (PAH) emission at 11.3 ${\mu}$m. Overall, we observe a stratification of the molecular gas, with the colder gas located in the nuclear spiral, ring, and connecting arms, while most warmer gas with higher velocity-dispersion fills 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 disc., Comment: 19 pages, 15 figures, 2. Accepted for publication in A&A

Bars are remarkable stellar structures that can transport gas toward centers and drive the secular evolution of galaxies. In this context, it is important to locate dynamical resonances associated with bars. For this study, we used ${Spitzer}$ near-infrared images as a proxy for the stellar gravitational potential and the ALMA CO(J=2-1) gas distribution from the PHANGS survey to determine the position of the main dynamical resonances associated with the bars in the PHANGS sample of 74 nearby star-forming galaxies. We used the gravitational torque method to estimate the location of the bar corotation radius ($R_{\rm CR}$), where stars and gas rotate at the same angular velocity as the bar. Of the 46 barred galaxies in PHANGS, we have successfully determined the corotation (CR) for 38 of them. The mean ratio of the $R_{\rm CR}$ to the bar radius ($R_{\rm bar}$) is $\mathcal{R} = R_{\rm CR}/R_{\rm bar} = 1.12$, with a standard deviation of $0.39$. This is consistent with the average value expected from theory and suggests that bars are predominantly fast. We also compared our results with other bar CR measurements from the literature, which employ different methods, and find good agreement ($\rho = 0.64$). Finally, using rotation curves, we have estimated other relevant resonances such as the inner Lindblad resonance (ILR) and the outer Lindblad resonance (OLR), which are often associated with rings. This work provides a useful catalog of resonances for a large sample of nearby galaxies and emphasizes the clear connection between bar dynamics and morphology.

We analyze JWST MIRI/MRS IFU observations of three Seyferts and showcase the intriguing polycyclic aromatic hydrocarbon (PAH) emission characteristics in regions of $\sim 500\,\rm pc$ scales over or around their active galactic nuclei (AGN). Combining the model predictions and the measurements of PAH features and other infrared emission lines, we find that the central regions containing a high fraction of neutral PAHs with small sizes, e.g., those in ESO137-G034, are in highly heated environments, due to collisional shock heating, with hard and moderately intense radiation fields. Such environments are proposed to be associated with inhibited growth or preferential erosion of PAHs, decreasing the average PAH size and the overall abundance of PAHs. We additionally find that the central regions containing a high fraction of ionized PAHs with large sizes, e.g., those in MCG-05-23-016, are likely experiencing severe photo-ionization because of the radiative effects from the radiative shock precursor besides the AGN. The severe photo-ionization can contribute to the ionization of all PAHs and further destruction of small PAHs. Overall, different Seyferts, even different regions in the same galaxy, e.g., those in NGC\,3081, can contain PAH populations of different properties. Specifically, Seyferts that exhibit similar PAH characteristics to ESO137-G034 and MCG-05-23-016 also tend to have similar emission line properties to them, suggesting that the explanations for PAH characteristics of ESO137-G034 and MCG-05-23-016 may also apply generally. These results have promising application in the era of JWST, especially in diagnosing different (i.e., radiative, and kinetic) AGN feedback modes., Comment: ApJL accepted on September 26th, title slightly modified in accordance with ApJL standards, Fig. 2 updated with additional labels

Utilizing JWST MIRI/MRS IFU observations of the kiloparsec scale central regions, we showcase the diversity of ionized gas distributions and kinematics in six nearby Seyfert galaxies included in the GATOS survey. Specifically, we present spatially resolved flux distribution and velocity field maps of six ionized emission lines covering a large range of ionization potentials ($15.8-97.1$ eV). Based on these maps, we showcase the evidence of ionized gas outflows in the six targets, and find some highly disturbed regions in NGC\,5728, NGC\,5506, and ESO137-G034. We propose AGN-driven radio jets plausibly play an important role in triggering these highly disturbed regions. With the outflow rates estimated based on [Ne~{\footnotesize V}] emission, we find the six targets tend to have ionized outflow rates converged to a narrower range than previous finding. These results have important implication for the outflow properties in AGN of comparable luminosity., Comment: 34 pages (11 pages in the appendix), 18 figures in the main text, ApJ in press (accepted on July 26th)

In this paper it is proved that the group $F\left(\frac32\right)$, a Thompson-style group with breaks in $\mathbb{Z}\left[\frac16\right]$ but whose slopes are restricted only to powers of $\frac32$, is finitely generated, with a generating set of two elements.

We analyze JWST MIRI/MRS observations of the infrared PAH bands in the nuclear and circumnuclear regions of local AGN from the GATOS Survey. In this work, we examine the PAH properties in the circumnuclear regions of AGN and AGN-outflows, and compare them to those in star-forming regions and the innermost regions of AGN. This study employs 4.9-28.1 micron sub-arcsecond angular resolution data to investigate the properties of PAH in three nearby sources (DL~30-40 Mpc). Our findings align with previous JWST studies, showing that the central regions of AGN show a larger fraction of neutral PAH molecules (i.e. elevated 11.3/6.2 and 11.3/7.7 PAH ratios) compared to star-forming galaxies. We find that the AGN might affect not only the PAH population in the innermost region but also in the extended regions up to ~kpc scales. By comparing our observations to PAH diagnostic diagrams, we find that, in general, regions located in the projected direction of the AGN-outflow occupy similar positions on the PAH diagnostic diagrams as those of the innermost regions of AGN. Star-forming regions that are not affected by the AGN in these galaxies share the same part of the diagram as Star-forming galaxies. We examine the potential of the PAH-H2 diagram to disentangle AGN versus star-forming activity. Our results suggest that in Sy-like AGN, illumination and feedback from the AGN might affect the PAH population at nuclear and kpc scales, in particular, the ionization state of the PAH grains. However, PAH sizes are rather similar. The carriers of the ionized PAH bands (6.2 and 7.7 micron) are less resilience than those of neutral PAH bands (11.3 micron), which might be particularly important for strongly AGN-host coupled systems. Therefore, caution must be applied when using PAH bands as star-formation rate indicators in these systems even at kpc scales, with the ionized ones being more affected by the AGN., Comment: Accepted for publication in A&A. 21 pages, 13 Figures

We present observations of the type-2 Seyfert NGC7172 obtained with the medium-resolution spectrometer (MRS) of the Mid-Infrared Instrument (MIRI) on board of the JWST. This galaxy hosts one of the lowest ionised gas mass outflow rates (Mout~0.005 M/yr) in a sample of six AGN with similar bolometric luminosities (log Lbol~44erg/s) within the Galactic Activity, Torus and Outflow Survey (GATOS). We aim to understand the properties of the ionised gas outflow, mainly using the emission lines from the neon transitions, that cover a broad range of ionisation potentials (IP) 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<25eV, e.g.[NeII]) trace the rotating disc emission. The high excitation lines (IP>90eV, e.g.[NeV]), 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 of ~350-500 km/s. In particular, [NeV] and [NeVI] lines reveal a biconical ionised gas outflow emerging N-S from the nuclear region, extending at least ~2.5"N and 3.8"S (projected distance of ~450 and 680 pc). 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, NGC7172 may be a case of weak coupling. Nevertheless, we found evidence for positive feedback in two distinct outflowing clumps at projected distances of 3.1" and 4.3" (i.e. ~560 and 780 pc) SW from the AGN. We estimated a star formation rate in these regions using the [NeII] and [NeIII] luminosities of 0.08 M/yr, that 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 ISM of the galaxy., Comment: 20 pages, 6 figures. Accepted for publication in A&A. Abstract adapted for the arxiv version

The ionisation cones of NGC5728 have a deficit of molecular gas based on millimetre observations of CO(2-1) emission. Although photoionisation from the active nucleus may lead to suppression of this transition, warm molecular gas can still be present. We report the detection of eight mid-infrared rotational H$_2$ lines throughout the central kiloparsec, including the ionisation cones, using integral field spectroscopic observations with JWST/MIRI MRS. The H$_2$ line ratios, characteristic of a power-law temperature distribution, indicate that the gas is warmest where it enters the ionisation cone through disk rotation, suggestive of shock excitation. In the nucleus, where the data can be combined with an additional seven ro-vibrational H$_2$ transitions, we find that moderate velocity (30 km s$^{-1}$) shocks in dense ($10^5$ cm$^{-3}$) gas, irradiated by an external UV field ($G_0 = 10^3$), do provide a good match to the full set. The warm molecular gas in the ionisation cone that is traced by the H$_2$ rotational lines has been heated to temperatures $>200$ K. Outside of the ionisation cone the molecular gas kinematics are undisturbed. However, within the ionisation cone, the kinematics are substantially perturbed, indicative of a radial flow, but one that is quantitatively different from the ionised lines. We argue that this outflow is in the plane of the disk, implying a short 50 pc acceleration zone up to speeds of about 400 km s$^{-1}$ followed by an extended deceleration over $\sim$700 pc where it terminates. The deceleration is due to both the radially increasing galaxy mass, and mass-loading as ambient gas in the disk is swept up., Comment: A&A accepted; 16 pages

We use a sample of 64 nearby (D=7-45 Mpc) disk galaxies including 45 AGN and 19 non-AGN, that have high spatial resolution multiline CO observations obtained with the ALMA and/or PdBI arrays to study the distribution of cold molecular gas in their circumunuclear disks (CND). We analyze whether the concentration of cold molecular gas changes as a function of the X-ray luminosity in the 2-10 keV range ($L_{\rm X}$). We also study the concentration of the hot molecular gas using NIR data obtained for the H2 1-0S(1) line. We find a turnover in the distribution of the cold molecular gas concentration as a function of $L_{\rm X}$ with a breakpoint which divides the sample into two branches: the AGN build-up branch ($L_{\rm X}\leq10^{41.5\pm0.3}$erg/s) and the AGN feedback branch ($L_{\rm X}\geq10^{41.5\pm0.3}$erg/s) . Lower luminosity AGN and non-AGN of the AGN build-up branch show high cold molecular gas concentrations and centrally peaked radial profiles on nuclear ($r\leq50$~pc) scales. Higher luminosity AGN of the AGN feedback branch, show a sharp decrease in the concentration of molecular gas and flat or inverted radial profiles. The cold molecular gas concentration index ($CCI$), defined as the ratio of surface densities at $r\leq50$~pc and $r\leq200$~pc , namely $CCI \equiv$~log$_{\rm 10}(\Sigma^{\rm gas}_{\rm 50}/\Sigma^{\rm gas}_{\rm 200}$), spans a factor ~4-5 between the galaxies lying at the high end of the AGN build-up branch and the galaxies of the AGN feedback branch. The concentration and radial distributions of the hot molecular gas in our sample follow less extreme trends as a function of the X-ray luminosity. These observations confirm, on a three times larger sample, previous evidence found by the GATOS survey that the imprint of AGN feedback on the CND-scale distribution of molecular gas is more extreme in higher luminosity Seyfert galaxies of the local universe., Comment: Accepted by Astronomy & Astrophysics (A&A) (14/06/2024), 26 pages, 16 figures

Spiral arms are some of the most spectacular features in disc galaxies, and also present in our own Milky Way. It has been argued that star formation should proceed more efficiently in spiral arms as a result of gas compression. Yet, observational studies have so far yielded contradictory results. Here we examine arm/interarm surface density contrasts at ~100 pc resolution in 28 spiral galaxies from the PHANGS survey. We find that the arm/interarm contrast in stellar mass surface density (Sigma_*) is very modest, typically a few tens of percent. This is much smaller than the contrasts measured for molecular gas (Sigma_mol) or star formation rate (Sigma_SFR) surface density, which typically reach a factor of ~2-3. Yet, Sigma_mol and Sigma_SFR contrasts show a significant correlation with the enhancement in Sigma_*, suggesting that the small stellar contrast largely dictates the stronger accumulation of gas and star formation. All these contrasts increase for grand-design spirals compared to multi-armed and flocculent systems (and for galaxies with high stellar mass). The median star formation efficiency (SFE) of the molecular gas is 16% higher in spiral arms than in interarm regions, with a large scatter, and the contrast increases significantly (median SFE contrast 2.34) for regions of particularly enhanced stellar contrast (Sigma_* contrast >1.97). The molecular-to-atomic gas ratio (Sigma_mol/Sigma_atom) is higher in spiral arms, pointing to a transformation of atomic to molecular gas. In conclusion, the boost in the star formation efficiency of molecular gas in spiral arms is generally modest or absent, except for locations with exceptionally large stellar contrasts. (abridged), Comment: 26 pages, 16 figures. Accepted for publication in A&A

We report the detection of extended (>0.5-1kpc) high-ionization [MgIV] 4.487 $\mu$m (80 eV) emission in four local luminous infrared galaxies observed with JWST/NIRSpec. Excluding the nucleus and outflow of the Type 1 active galactic nucleus (AGN) in the sample, we find that the [MgIV] luminosity is well correlated with that of H recombination lines, which mainly trace star forming clumps in these objects, and that the [ArVI] 4.530 $\mu$m (75 eV), usually seen in AGN, is undetected. On 100-400pc scales, the [MgIV] line profiles are broader (sigma([MgIV])=90 +- 25 km/s) and shifted (Delta_v up to +- 50 km/s) compared to those of the H recombination lines and lower ionization transitions (e.g., sigma(Hu-12)=57 +- 15 km/s). The [MgIV] kinematics follow the large scale rotating velocity field of these galaxies and the broad [MgIV] profiles are compatible with the broad wings detected in the H recombination lines. Based on these observational results, extended highly ionized gas more turbulent than the ambient interstellar medium, possibly as a result of ionizing shocks associated with star-formation, is the most likely origin of the [MgIV] emission. We also computed new grids of photoionization and shock models to investigate where the [MgIV] line originates. Shocks with velocities of 100-130 km/s reproduce the observed line ratios and the [MgIV] luminosity agrees with that expected from the mechanical energy released by supernove (SNe) in these regions. Therefore, these models support shocks induced by SNe as the origin of the [MgIV] line. Future studies on the stellar feedback from SNe will benefit from the [MgIV] line that is little affected by obscuration and, in absence of an AGN, can only be produced by shocks due to its high ionization potential., Comment: Accepted for publication in A&A Letters. 6 pages, 6 figures. Appendix: 8 pages, 10 figures

Polar dust has been discovered in a number of local Active Galactic Nuclei (AGN), with radiation-driven torus models predicting a wind to be its main driver. However, little is known about its characteristics, spatial extent, or connection to the larger scale outflows. We present the first JWST/MIRI study aimed at imaging polar dust by zooming onto the centre of ESO 428-G14, part of the Galaxy Activity, Torus, and Outflow Survey (GATOS) survey of local AGN. We detect extended mid-infrared (MIR) emission within 200 pc from the nucleus. This polar structure is co-linear with a radio jet and lies perpendicular to a molecular gas lane that feeds and obscures the nucleus. Its morphology bears a striking resemblance to that of gas ionised by the AGN in the narrow-line region. We demonstrate that part of this spatial correspondence is due to contamination within the JWST filter bands from strong emission lines. Correcting for the contamination, we find the morphology of the dust continuum to be more compact, though still clearly extended out to ~ 100 pc. We estimate the emitting dust has a temperature of ~ 120 K. Using simple models, we find that the heating of small dust grains by the radiation from the central AGN and/or radiative jet-induced shocks is responsible for the extended MIR emission. Radiation-driven dusty winds from the torus is unlikely to be important. This has important implications for scales to which AGN winds can carry dust and dense gas out into their host galaxies., Comment: 17 pages, 8 figures, 2 tables. Published in MNRAS

Strain is an effective strategy to modulate the optoelectronic properties of 2D materials, but it has been almost unexplored in layered hybrid organic-inorganic metal halide perovskites (HOIPs) due to their complex band structure and mechanical properties. Here, we investigate the temperature-dependent microphotoluminescence (PL) of 2D $(C_6H_5CH_2CH_2NH_3)_2Cs_3Pb_4Br_{13}$ HOIP subject to biaxial strain induced by a $SiO_2$ ring platform on which flakes are placed by viscoelastic stamping. At 80 K, we found that a strain of <1% can change the PL emission from a single peak (unstrained) to three well-resolved peaks. Supported by micro-Raman spectroscopy, we show that the thermomechanically generated strain modulates the bandgap due to changes in the octahedral tilting and lattice expansion. Mechanical simulations demonstrate the coexistence of tensile and compressive strain along the flake. The observed PL peaks add an interesting feature to the rich phenomenology of photoluminescence in 2D HOIPs, which can be exploited in tailored sensing and optoelectronic devices., Comment: 4 figures

The outflowing molecular gas in the circumnuclear disk (CND) of the nearby (D=14 Mpc) AGN-starburst composite galaxy NGC 1068 is considered as a manifestation of ongoing AGN feedback. The large spread of velocities from the outflowing gas is likely driving various kinds of shock chemistry across the CND. We performed a multiline molecular study using CH3OH with the aim of characterizing the gas properties probed by CH3OH in the CND of NGC 1068, and investigating its potential association with molecular shocks. Multi-transition CH3OH were imaged at the resolution of 0.''5-0.''8 with the Atacama Large Millimeter/submillimeter Array (ALMA). We performed non-LTE radiative transfer analysis coupled with a Bayesian inference process in order to determine the gas properties such as the gas volume density and the gas kinetic temperature. The gas densities traced by CH3OH point to $\sim 10^{6}$ cm\textsuperscript{-3} across all the CND regions. The gas kinetic temperature cannot be well constrained in any of the CND regions though the inferred temperature is likely low ($\lesssim 100$ K).The low gas temperature traced by CH3OH suggests shocks and subsequent fast cooling as the origin of the observed gas-phase CH3OH abundance. We also note that the E-/A- isomer column density ratio inferred is fairly close to unity, which is interestingly different from the Galactic measurements in the literature. It remains inconclusive whether CH3OH exclusively traces slow and non-dissociative shocks, or whether the CH3OH abundance can actually be boosted in both fast and slow shocks.

How galaxies regulate nuclear growth through gas accretion by supermassive black holes (SMBHs) is one of the most fundamental questions in galaxy evolution. One potential way to regulate nuclear growth is through a galactic wind that removes gas from the nucleus. It is unclear whether galactic winds are powered by jets, mechanical winds, radiation, or via magnetohydrodynamic (MHD) processes. Compact obscured nuclei (CONs) represent a significant phase of galactic nuclear growth. These galaxies hide growing SMBHs or unusual starbursts in their very opaque, extremely compact (r $<$ 100 pc) centres. They are found in approximately 30 % of the luminous and ultra-luminous infrared galaxy (LIRG and ULIRG) population. Here, we present high-resolution ALMA observations ($\sim$30 mas, $\sim$5 pc) of ground-state and vibrationally excited HCN towards ESO 320-G030 (IRAS 11506-3851). ESO 320-G030 is an isolated luminous infrared galaxy known to host a compact obscured nucleus and a kiloparsec-scale molecular wind. Our analysis of these high-resolution observations excludes the possibility of a starburst-driven wind, a mechanically or energy driven active galactic nucleus (AGN) wind, and exposes a molecular MDH wind. These results imply that the nuclear evolution of galaxies and the growth of SMBHs are similar to the growth of hot cores or protostars where gravitational collapse of the nuclear torus drives a MHD wind. These results mean galaxies are capable, in part, of regulating the evolution of their nuclei without feedback., Comment: 9 pages, 10 figures

We present new optical GTC/MEGARA seeing-limited (0.9") integral-field observations of NGC 5506, together with ALMA observations of the CO(3-2) transition at a 0.2" (25 pc) resolution. NGC 5506 is a luminous (bolometric luminosity of $\sim 10^{44}$ erg/s) nearby (26 Mpc) Seyfert galaxy, part of the Galaxy Activity, Torus, and Outflow Survey (GATOS). We modelled the CO(3-2) kinematics with 3D-Barolo, revealing a rotating and outflowing cold gas ring within the central 1.2 kpc. We derived an integrated cold molecular gas mass outflow rate for the ring of 8 M$_{\odot}$/yr. We fitted the optical emission lines with a maximum of two Gaussian components to separate rotation from non-circular motions. We detected high [OIII]$\lambda$5007 projected velocities (up to 1000 km/s) at the active galactic nucleus (AGN) position, decreasing with radius to an average 330 km/s around 350 pc. We also modelled the [OIII] gas kinematics with a non-parametric method, estimating the ionisation parameter and electron density in every spaxel, from which we derived an ionised mass outflow rate of 0.076 M$_{\odot}$/yr within the central 1.2 kpc. Regions of high CO(3-2) velocity dispersion, extending to projected distances of 350 pc from the AGN, appear to be the result from the interaction of the AGN wind with molecular gas in the galaxy's disc. Additionally, we find the ionised outflow to spatially correlate with radio and soft X-ray emission in the central kiloparsec. We conclude that the effects of AGN feedback in NGC 5506 manifest as a large-scale ionised wind interacting with the molecular disc, resulting in outflows extending to radial distances of 610 pc, Comment: 24 pages, 30 figures, accepted for publication by Astronomy & Astrophysics

We describe standard forms for elements of the higher-dimensional Thompson groups $nV$ arising from gridding subdivision processes. These processes lead to standard normal form descriptions for elements in these groups, and sizes of these standard forms estimate the word length with respect to finite generating sets. These gridded forms lead to standard algebraic descriptions as well, with respect to the both infinite and finite generating sets for these groups., Comment: 18 pages, 6 figures

Nuclear regions of ultraluminous and luminous infrared galaxies (U/LIRGs) are powered by starbursts and/or active galactic nuclei (AGNs). These regions are often obscured by extremely high columns of gas and dust. Molecular lines in the submillimeter windows have the potential to determine the physical conditions of these compact obscured nuclei (CONs). We aim to reveal the distributions of HCN and HCO$^+$ emission in local U/LIRGs and investigate whether and how they are related to galaxy properties. Using ALMA, we have conducted sensitive observations of the HCN J=3--2 and HCO$^+$ J=3--2 lines toward 23 U/LIRGs in the local Universe (z < 0.07) with a spatial resolution of ~0.3" (~50--400 pc). We detected both HCN and HCO$^+$ in 21 galaxies, only HCN in one galaxy, and neither in one galaxy. The global HCN/HCO$^+$ line ratios, averaged over scales of ~0.5--4 kpc, range from 0.4 to 2.3, with an unweighted mean of 1.1. These line ratios appear to have no systematic trend with bolometric AGN luminosity or star formation rate. The line ratio varies with position and velocity within each galaxy, with an average interquartile range of 0.38 on a spaxel-by-spaxel basis. In eight out of ten galaxies known to have outflows and/or inflows, we found spatially and kinematically symmetric structures of high line ratios. These structures appear as a collimated bicone in two galaxies and as a thin spherical shell in six galaxies. Non-LTE analysis suggests that the high HCN/HCO$^+$ line ratio in outflows is predominantly influenced by the abundance ratio. Chemical model calculations indicate that the enhancement of HCN abundance in outflows is likely due to high-temperature chemistry triggered by shock heating. These results imply that the HCN/HCO$^+$ line ratio can aid in identifying the outflow geometry when the shock velocity of the outflows is sufficiently high to heat the gas., Comment: 52 pages, 35 figures, accepted for publication in Astronomy & Astrophysics

The superb image quality, stability and sensitivity of the JWST permit deconvolution techniques to be pursued with a fidelity unavailable to ground-based observations. We present an assessment of several deconvolution approaches to improve image quality and mitigate effects of the complex JWST point spread function (PSF). The optimal deconvolution method is determined by using WebbPSF to simulate JWST's complex PSF and MIRISim to simulate multi-band JWST/Mid-Infrared Imager Module (MIRIM) observations of a toy model of an active galactic nucleus (AGN). Five different deconvolution algorithms are tested: (1) Kraken deconvolution, (2) Richardson-Lucy, (3) Adaptive Imaging Deconvolution Algorithm, (4) Sparse regularization with the Condat-V\~u algorithm, and (5) Iterative Wiener Filtering and Thresholding. We find that Kraken affords the greatest FWHM reduction of the nuclear source of our MIRISim observations for the toy AGN model while retaining good photometric integrity across all simulated wavebands. Applying Kraken to Galactic Activity, Torus, and Outflow Survey (GATOS) multi-band JWST/MIRIM observations of the Seyfert 2 galaxy NGC 5728, we find that the algorithm reduces the FWHM of the nuclear source by a factor of 1.6-2.2 across all five filters. Kraken images facilitate detection of a SE to NW $\thicksim$2".5 ($\thicksim$470 pc, PA $\simeq$115\deg) extended nuclear emission, especially in the longest wavelengths. We demonstrate that Kraken is a powerful tool to enhance faint features otherwise hidden in the complex JWST PSF., Comment: 32 pages, 23 figures, published in AJ 2024 February 7

We present a new physically-motivated model for estimating the molecular line emission in active galaxies. The model takes into account (i) the internal density structure of giant molecular clouds (GMCs), (ii) the heating associated both to stars and to the active galactic nuclei (AGN), respectively producing photodissociation regions (PDRs) and X-ray dominated regions (XDRs) within the GMCs, and (iii) the mass distribution of GMCs within the galaxy volume. The model needs, as input parameters, the radial profiles of molecular mass, far-UV flux and X-ray flux for a given galaxy, and it has two free parameters: the CO-to-H2 conversion factor $\alpha_{CO}$, and the X-ray attenuation column density $N_H$. We test this model on a sample of 24 local ($z \leq 0.06$) AGN-host galaxies, simulating their carbon monoxide spectral line energy distribution (CO SLED). We compare the results with the available observations and calculate, for each galaxy, the best ($\alpha_{CO}$, $N_H$) with a Markov chain Monte Carlo algorithm, finding values consistent with those present in the literature. We find a median $\alpha_{CO} = 4.8$ M$_{\odot}$ (K km s$^{-1}$ pc$^{2}$)$^{-1}$ for our sample. In all the modelled galaxies, we find the XDR component of the CO SLED to dominate the CO luminosity from $J_{\text{upp}} \geq 4$. We conclude that, once a detailed distribution of molecular gas density is taken into account, PDR emission at mid-/high-$J$ becomes negligible with respect to XDR., Comment: 19 pages, 36 figures, resubmitted to MNRAS after addressing referee's comments

We use JWST/MIRI MRS spectroscopy of a sample of six local obscured type 1.9/2 active galactic nuclei (AGN) to compare their nuclear mid-IR absorption bands with the level of nuclear obscuration traced by X-rays. This study is the first to use sub-arcsecond angular resolution data of local obscured AGN to investigate the nuclear mid-IR absorption bands with a wide wavelength coverage (4.9-28.1 $\mu$m). All the nuclei show the 9.7 $\mu$m silicate band in absorption. We compare the strength of the 9.7 and 18 $\mu$m silicate features with torus model predictions. The observed silicate features are generally well explained by clumpy and smooth torus models. We report the detection of the 6 $\mu$m dirty water ice band (i.e., a mix of water and other molecules such as CO and CO$_2$) at sub-arcsecond scales ($\sim$0.26 arcsec at 6 $\mu$m; inner $\sim$50 pc) in a sample of local AGN with different levels of nuclear obscuration in the range log N$_{\rm H}^{\rm X-Ray}$(cm$^{-2}$)$\sim22-25$. We find a good correlation between the 6 $\mu$m water ice optical depths and N$_{\rm H}^{\rm X-Ray}$. This result indicates that the water ice absorption might be a reliable tracer of the nuclear intrinsic obscuration in AGN. The weak water ice absorption in less obscured AGN (log N$_H^{X-ray}$ (cm$^{-2}$)$\lesssim$23.0 cm$^{-2}$) might be related to the hotter dust temperature ($>$T$_{sub}^{H_2O}\sim$110 K) expected to be reached in the outer layers of the torus due to their more inhomogeneous medium. Our results suggest it might be necessary to include the molecular content, such as, H$_2$O, aliphatic hydrocarbons (CH-) and more complex PAH molecules in torus models to better constrain key parameters such as the torus covering factor (i.e. nuclear obscuration)., Comment: Accepted for publication in A&A. 13 pages, 12 Figures

We analyze JWST/NIRSpec observations of the CO rovibrational v=1-0 band at 4.67um around the dust-embedded southern active galactic nucleus (AGN) of NGC3256 (d=40Mpc; L(IR)=10^11.6 Lsun). We classify the CO v=1-0 spectra into three categories based on the behavior of P- and R-branches of the band: (a) both branches in absorption toward the nucleus; (b) P-R asymmetry (P-branch in emission and R-branch in absorption) along the disk of the galaxy; and (c) both branches in emission in the outflow region above and below the disk. In this paper, we focus on the outflow. The CO v=1-0 emission can be explained by the vibrational excitation of CO in the molecular outflow by the bright mid-IR ~4.7um continuum from the AGN up to r~250pc. We model the ratios between the P(J+2) and R(J) transitions of the band to derive the physical properties (column density, kinetic temperature, and CO-to-H2 conversion factor, alpha_CO) of the outflowing gas. We find that the 12CO v=1-0 emission is optically thick for J<4, while the 13CO v=1-0 emission remains optically thin. From the P(2)/R(0) ratio, we identify a temperature gradient in the outflow from >40K in the central 100pc to <15K at 250pc sampling the cooling of the molecular gas in the outflow. We used three methods to derive alpha_CO in eight 100pc (0.5") apertures in the outflow by fitting the P(J+2)/R(J) ratios with non-LTE models. We obtain low median alpha_CO factors (0.40-0.61) x 3.2e-4/[CO/H2] Msun (K km/s/pc2)^-1 in the outflow regions. This implies that outflow rates and energetics might be overestimated if a ULIRG-like alpha_CO, which is 1.3-2 times larger, is assumed. The reduced alpha_CO can be explained if the outflowing molecular clouds are not virialized. We also report the first extragalactic detection of a broad (sigma=0.0091um) spectral feature at 4.645um associated with aliphatic deuterium on polycyclic aromatic hydrocarbons (D_n-PAH)., Comment: Accepted for publication in A&A.17 pages, 13 figures

We introduce a machine-learning approach (denoted Symmetry Seeker Neural Network) capable of automatically discovering discrete symmetry groups in physical systems. This method identifies the finite set of parameter transformations that preserve the system's physical properties. Remarkably, the method accomplishes this without prior knowledge of the system's symmetry or the mathematical relationships between parameters and properties. Demonstrating its versatility, we showcase examples from mathematics, nanophotonics, and quantum chemistry.

We present the first systematic study of the molecular gas and star formation efficiency in a sample of ten narrow-line Seyfert 1 galaxies selected to have X-ray Ultra Fast Outflows and, therefore, to potentially show AGN feedback effects. CO observations were obtained with the IRAM 30m telescope in six galaxies and from the literature for four galaxies. We derived the stellar mass, star formation rate, AGN and FIR dust luminosities by fitting the multi-band spectral energy distributions with the CIGALE code. Most of the galaxies in our sample lie above the main sequence (MS) and the molecular depletion time is one to two orders of magnitude shorter than the one typically measured in local star-forming galaxies. Moreover, we found a promising correlation between the star formation efficiency and the Eddington ratio, as well as a tentative correlation with the AGN luminosity. The role played by the AGN activity in the regulation of star formation within the host galaxies of our sample remains uncertain (little or no effect? positive feedback?). Nevertheless, we can conclude that quenching by the AGN activity is minor and that star formation will likely stop in a short time due to gas exhaustion by the current starburst episode., Comment: Published in MNRAS, Volume 524, Issue 2, Pages 3130-3145

We present hydrodynamic simulations of the interstellar medium (ISM) within the circumnuclear disk (CND) of a typical AGN-dominated galaxy influenced by mechanical feedback from an active galactic nucleus(AGN). The simulations are coupled with the CHIMES non-equilibrium chemistry network to treat the radiative-cooling and AGN-heating. A focus is placed on the central 100 pc scale where AGN outflows are coupled to the ISM and constrained by observational Seyfert-2 galaxies. AGN-feedback models are implemented with different wind-velocity and mass-loading factors. We post-process the simulation snapshots with a radiative-transfer code to obtain the molecular emission lines. We find that the inclusion of an AGN promotes the formation of CO in clumpy and dense regions surrounding supermassive-blackholes (SMBH). The CO(1-0) intensity maps ($<$6 Myr) in the CND seem to match well with observations of NGC 1068 with a best match for a model with 5000 $\rm km/s$ wind-velocity and a high mass-loading factor. We attempt to discern between competing explanations for the apparent counter-rotating gas disk in the NGC 1068 through an analysis of kinematic maps of the CO line emission. We suggest that mechanical AGN-feedback could explain the alignment-stability of position-angle across the different CND radii around the SMBH through momentum and energy loading of the wind. It is the wind-velocity that drives the disk out of alignment on a 100 pc scale for a long period of time. The position-velocity diagrams are in broad agreement with the predicted Keplerian rotation-curve in the model without-AGN, but the AGN models exhibit a larger degree of scatter, in better agreement with NGC 1068 observations., Comment: 16 pages, 13 figures. Accepted for publication in MNRAS

We present new ALMA observations of the CO(3-2) transition and 854micron continuum at 0.06-0.3" resolution, together with new VLT/SINFONI observations of NGC7172. This is a luminous (bolometric luminosity of ~10^44 erg/s) Seyfert galaxy that belongs to the Galaxy Activity, Torus, and Outflow Survey (GATOS). The CO(3-2) observations reveal the presence of a highly inclined cold molecular gas ring with an approximate radius of 3-4"~540-720 pc, which is likely associated with an inner Lindblad resonance of a putative stellar bar. There are noncircular motions in the VLT/SINFONI [SiVI]1.96micron and H2 at 2.12micron, and ALMA CO(3-2) velocity fields. After subtracting the stellar velocity field, we detected [SiVI] blueshifted velocities of a few hundred km/s to the south of the AGN. They trace outflowing ionized gas outside the plane of the galaxy and out to projected distances of ~200 pc. The CO(3-2) position-velocity diagram along the kinematic minor axis displays noncircular motions with observed velocities of up to ~150 km/s. Assuming that these are taking place in the disk of the galaxy, the observed velocity signs imply that the molecular gas ring is not only rotating but also outflowing. We derived an integrated cold molecular gas mass outflow rate of ~40 Msun/yr for the ring. Using the 854micron map, we resolved a 32 pc radius torus with a gas mass of 8x10^5 Msun. These torus properties are similar to other Seyfert galaxies in the GATOS sample. We measured a decreased cold molecular gas concentration in the nuclear-torus region relative to the circumnuclear region when compared to other less luminous Seyfert galaxies. We conclude that the effects of AGN feedback in NGC7172, which are likely caused by the AGN wind and/or the moderate luminosity radio jet, are seen as a large-scale outflowing molecular gas ring and accompanying redistribution of molecular gas in the nuclear regions., Comment: Accepted for publication to A&A

Fits the infrared spectra from the nuclear regions of AGN can place constraints on the dust properties, distribution, and geometry by comparison with models. However, none of the currently available models fully describe the observations of AGN currently available. Among the aspects least explored, here we focus on the role of dust grain size. We offer the community a new spectral energy distribution (SED) library, hereinafter [GoMar23] model, which is based on the two-phase torus model developed before with the inclusion of the grain size as a model parameter, parameterized by the maximum grain size Psize or equivalently the mass-weighted average grain size < P >. We created 691,200 SEDs using the SKIRT code, where the maximum grain size can vary within the range Psize = 0.01 - 10.0um ( < P >= 0.007 - 3.41um). We fit this new and several existing libraries to a sample of 68 nearby and luminous AGNs with Spitzer/IRS spectra dominated by AGN-heated dust. We find that the [GoMar23] model can adequately reproduce up to 85-88% of the spectra. The dust grain size parameter significantly improves the final fit in up to 90% of these spectra. Statistical tests indicate that the grain size is the third most important parameter in the fitting procedure (after the size and half opening angle of the torus). The requirement of a foreground extinction by our model is lower compared to purely clumpy models. We find that 41% of our sample requires that the maximum dust grain size is as large as Psize =10um (< P >= 3.41um). Nonetheless, we also remark that disk+wind and clumpy torus models are still required to reproduce the spectra of a non-negligible fraction of objects, suggesting the need for several dust geometries to explain the infrared continuum of AGN. This work provides tentative evidence for dust grain growth in the proximity of the AGN., Comment: 26 pages, 14 figures, 4 tables, accepted for publication in A&A

We present a spatially-resolved study of the ionised gas in the central 2 kpc of the Seyfert 2 galaxy NGC 2110 and investigate the role of its moderate luminosity radio jet (kinetic radio power of $P_\mathrm{jet} = 2.3 \times 10^{43}\mathrm{erg\ s^{-1}}$). We use new optical integral-field observations taken with the MEGARA spectrograph at GTC. We fit the emission lines with a maximum of two Gaussian components, except at the AGN position where we used three. Aided by existing stellar kinematics, we use the observed velocity and velocity dispersion of the emission lines to classify the different kinematic components. The disc component is characterised by lines with $\sigma \sim 60-200\ \mathrm{km\ s^{-1}}$. The outflow component has typical values of $\sigma \sim 700\ \mathrm{km\ s^{-1}}$ and is confined to the central 400 pc, which is coincident with linear part of the radio jet detected in NGC 2110. At the AGN position, the [O III]$\lambda$5007 line shows high velocity components reaching at least $1000\ \mathrm{km\ s^{-1}}$. This and the high velocity dispersions indicate the presence of outflowing gas outside the galaxy plane. Spatially-resolved diagnostic diagrams reveal mostly LI(N)ER-like excitation in the outflow and some regions in the disc, which could be due to the presence of shocks. However, there is also Seyfert-like excitation beyond the bending of the radio jet, probably tracing the edge of the ionisation cone that intercepts with the disc of the galaxy. NGC 2110 follows well the observational trends between the outflow properties and the jet radio power found for a few nearby Seyfert galaxies. All these pieces of information suggest that part of observed ionised outflow in NGC 2110 might be driven by the radio jet. However, the radio jet was bent at radial distances of 200 pc (in projection) from the AGN, and beyond there, most of the gas in the galaxy disc is rotating., Comment: 25 pages, 22 figures, 3 tables. Accepted for publication in A&A

The factor relating CO emission to molecular hydrogen column density, XCO, is still subject to uncertainty, in particular at low metallicity. Here, to quantify XCO at two different spatial resolutions, we exploit a dust-based method together with ALMA 12-m and ACA data and HI maps of three nearby metal-poor starbursts, NGC625, NGC1705, and NGC5253. Dust opacity at 250pc resolution is derived based on dust temperatures estimated by fitting two-temperature modified blackbodies to Herschel PACS data. By using the HI maps, we are then able to estimate dust-to-gas ratios in the atomic-gas dominated regions, and infer total gas column densities and H2 column densities as the difference with HI. Finally, from the ACA CO(1-0) maps, we derive XCO. We use a similar technique with 40 pc ALMA 12-m data for the three galaxies, but instead derive dust attenuation at 40 pc resolution from reddening maps based on VLT/MUSE data. At 250 pc resolution, XCO $\sim$ 10^22 - 10^23 cm^-2 / K.km/s, 5-1000 times the Milky Way value, with much larger values than would be expected from a simple metallicity dependence. Instead at 40 pc resolution, XCO again shows large variation, but is roughly consistent with a power-law metallicity dependence, given the Z $\sim$ 1/3 Zsun metal abundances of our targets. The large scatter in both estimations could imply additional parameter dependence, that we have investigated by comparing XCO with the observed velocity-integrated brightness temperatures, ICO, as predicted by recent simulations. Indeed, larger XCO is significantly correlated with smaller ICO, but with slightly different slopes and normalizations than predicted by theory. Such behavior can be attributed to the increasing fraction of CO-faint H2 gas with lower spatial resolution. This confirms the idea the XCO is multi-variate, depending not only on metallicity but also on CO brightness temperature and beam size., Comment: A&A, in press (16 figures, 1 appendix)

HNCO and SiO are well known shock tracers and have been observed in nearby galaxies, including the nearby (D=3.5 Mpc) starburst galaxy NGC 253. The simultaneous detection of these two species in regions where the star formation rate is high may be used to study the shock history of the gas. We perform a multi-line molecular study using these two shock tracers (SiO and HNCO) with the aim of characterizing the gas properties. We also explore the possibility of reconstructing the shock history in NGC 253's Central Molecular Zone (CMZ). Six SiO transitions and eleven HNCO transitions were imaged at high resolution $1''.6$ (28 pc) with the Atacama Large Millimeter/submillimeter Array (ALMA) as part of the ALCHEMI Large Programme. Both non-LTE radiative transfer analysis and chemical modelling were performed in order to characterize the gas properties, and to investigate the chemical origin of the emission. The non-LTE radiative transfer analysis coupled with Bayesian inference shows clear evidence that the gas traced by SiO has different densities and temperatures than that traced by HNCO, with an indication that shocks are needed to produce both species. Chemical modelling further confirms such a scenario and suggests that fast and slow shocks are responsible for SiO and HNCO production, respectively, in most GMCs. We are also able to infer the physical characteristics of the shocks traced by SiO and HNCO for each GMC. Radiative transfer and chemical analysis of the SiO and HNCO in the CMZ of NGC 253 reveal a complex picture whereby most of the GMCs are subjected to shocks. We speculate on the possible shock scenarios responsible for the observed emission and provide potential history and timescales for each shock scenario. Higher spatial resolution observations of these two species are required in order to quantitatively differentiate between scenarios.

Strong coupling (SC) between light and matter excitations such as excitons and molecular vibrations bear intriguing potential for controlling chemical reactivity, conductivity or photoluminescence. So far, SC has been typically achieved either between mid-infrared (mid-IR) light and molecular vibrations or between visible light and excitons. Achieving SC simultaneously in both frequency bands may open unexplored pathways for manipulating material properties. Here, we introduce a polaritonic nanoresonator (formed by h-BN layers placed on Al ribbons) hosting surface plasmon polaritons (SPPs) at visible frequencies and phonon polaritons (PhPs) at mid-IR frequencies, which simultaneously couple to excitons and atomic vibration in an adjacent molecular layer (CoPc). Employing near-field optical nanoscopy, we first demonstrate the co-localization of strongly confined near-fields at both visible and mid-IR frequencies. After covering the nanoresonator structure with a layer of CoPc molecules, we observe clear mode splittings in both frequency ranges by far-field transmission spectroscopy, unambiguously revealing simultaneous SPP-exciton and PhP-vibron coupling. Dual-band SC may be exploited for manipulating the coupling between excitons and molecular vibrations in future optoelectronics, nanophotonics, and quantum information applications.

We present the results from the study of the resolved distribution of cold molecular gas around eight young (<10^6 yr), peaked-spectrum radio galaxies. This has allowed us to trace the interplay between the radio jets and the surrounding medium. For three of these sources we present new CO(1-0) observations, obtained with NOEMA. In two targets, we also detected CN lines, both in emission and absorption. Combining the new observations with already published data, we discuss the main results obtained. Although we found that a large fraction of the cold molecular gas is distributed in disc-like rotating structures, in most of the sources high turbulence and deviations from purely quiescent gas (including outflows) were observed in the region co-spatial with the radio continuum emission. This suggests the presence of an interaction between radio plasma and cold molecular gas. We found that newly born and young radio jets, even those with low power i.e., P_jet<10^45 erg/s), can drive massive outflows of cold, molecular gas. The outflows are, however, limited to the sub-kpc regions and likely short lived. On larger scales (a few kpc), we observed cases where the molecular gas appears to avoid the radio lobes and, instead, wraps around them. The results suggest the presence of an evolutionary sequence, consistent with simulations, where the type of impact of the radio plasma changes as the jet expands, going from a direct jet-cloud interaction on sub-kpc scales to a gentler pushing aside of the gas, increasing its turbulence and likely limiting its cooling. This effect can be mediated by the cocoon of shocked gas inflated by the jet-cloud interactions. Building larger samples of young and evolved radio sources for observation at a similar depth and spatial resolution to test this scenario is now needed and may be possible thanks to more data becoming available in the growing public archives., Comment: Published in Galaxies Special Issue "The Symbiosis between Radio Source and Galaxy Evolution" O'Dea & Baum (guest editors), 20 pages, 8 figures; https://doi.org/10.3390/galaxies11010024

In order to investigate the impact of radio jets on the interstellar medium (ISM) of galaxies hosting active galactic nuclei (AGN), we present subarcsecond resolution Atacama Large Millimeter/submillimeter Array (ALMA) CO(2-1) and CO(3-2) observations of the Teacup galaxy. This is a nearby ($D_{\rm L}$=388 Mpc) radio-quiet type-2 quasar (QSO2) with a compact radio jet ($P_{\rm jet}\approx$10$^{43}$ erg s$^{-1}$) that subtends a small angle from the molecular gas disc. Enhanced emission line widths perpendicular to the jet orientation have been reported for several nearby AGN for the ionised gas. For the molecular gas in the Teacup, not only do we find this enhancement in the velocity dispersion but also a higher brightness temperature ratio (T32/T21) perpendicular to the radio jet compared to the ratios found in the galaxy disc. Our results and the comparison with simulations suggest that the radio jet is compressing and accelerating the molecular gas, and driving a lateral outflow that shows enhanced velocity dispersion and higher gas excitation. These results provide further evidence that the coupling between the jet and the ISM is relevant to AGN feedback even in the case of radio-quiet galaxies., Comment: 9 pages, 8 figures. Accepted for publication in A&A Letters

The Narrow Line Seyfert 1 Galaxy IRAS17020+4544 is one of the few AGN where a galaxy-scale energy-conserving outflow was revealed. This paper reports on NOEMA observations addressed to constrain the spatial scale of the CO emission in outflow. The molecular outflowing gas is resolved in five components tracing approaching and receding gas, all located at a distance of 2-3~kpc on the West and East side of the active nucleus. This high velocity gas (up to v_out=~1900 km/s) is not coincident with the rotation pattern of the CO gas in the host galaxy disk. The estimated mass outflow rate shows that with a global mass output of $\dot{M}_{H_2}$=~139$\pm$20$~M_\odot$~yr$^{-1}$, this powerful galaxy-scale outflow is consistent with the wind conserving its energy, and with a momentum rate boost of a factor of ~30 compared to the momentum rate of the nuclear X-ray wind. Preliminary results from ancillary X-ray (Chandra) and radio images (e-MERLIN) are reported. While the nature of the radio source is not conclusive, the Chandra image may tentatively trace extended emission, as expected by an expanding bubble of hot X-ray gas. The outcome of the NOEMA analysis and of past and ongoing publications dedicated to the description of the outflow multi-band phenomenology in IRAS17020+4544 concur to provide compelling reasons to postulate that an outflow shocking with the galaxy interstellar medium is driving the multi-phase wind in this peculiar AGN., Comment: Accepted for publication on MNRAS, 17 pages, 7 figures

We present new HCN and HCO$^+$ ($J$=3-2) images of the nearby star-forming galaxies (SFGs) NGC 3351, NGC 3627, and NGC 4321. The observations, obtained with the Morita ALMA Compact Array, have a spatial resolution of $\sim$290-440 pc and resolve the inner $R_\textrm{gal} \lesssim$ 0.6-1 kpc of the targets, as well as the southern bar end of NGC 3627. We complement this data set with publicly available images of lower excitation lines of HCN, HCO$^+$, and CO and analyse the behaviour of a representative set of line ratios: HCN(3-2)/HCN(1-0), HCN(3-2)/HCO$^+$(3-2), HCN(1-0)/CO(2-1), and HCN(3-2)/CO(2-1). Most of these ratios peak at the galaxy centres and decrease outwards. We compare the HCN and HCO$^+$ observations with a grid of one-phase, non-local thermodynamic equilibrium (non-LTE) radiative transfer models and find them compatible with models that predict subthermally excited and optically thick lines. We study the systematic variations of the line ratios across the targets as a function of the stellar surface density ($\Sigma_\textrm{star}$), the intensity-weighted CO(2-1) ($\langle I_\text{CO}\rangle$), and the star formation rate surface density ($\Sigma_\text{SFR}$). We find no apparent correlation with $\Sigma_\text{SFR}$, but positive correlations with the other two parameters, which are stronger in the case of $\langle I_\text{CO}\rangle$. The HCN/CO-$\langle I_\text{CO}\rangle$ relations show $\lesssim$0.3 dex galaxy-to-galaxy offsets, with HCN(3-2)/CO(2-1)-$\langle I_\text{CO}\rangle$ being $\sim$2 times steeper than HCN(1-0)/CO(2-1). In contrast, the HCN(3-2)/HCN(1-0)-$\langle I_\text{CO}\rangle$ relation exhibits a tighter alignment between galaxies. We conclude that the overall behaviour of the line ratios cannot be ascribed to variations in a single excitation parameter (e.g. density or temperature)., Comment: Accepted for publication in A&A. 14 pages, 8 figures

Molecular line ratios, such as HCN(1-0)/HCO$^+$(1-0) and HCN(4-3)/CS(7-6), are routinely used to identify active galactic nuclei (AGN) activity in galaxies. Such ratios are, however, hard to interpret as they are highly dependent on the physics and energetics of the gas, and hence can seldom be used as a unique, unambiguous diagnostic. We used the composite galaxy NGC 1068 as a `laboratory' to investigate whether molecular line ratios between HCN, HCO$^+$, and CS are useful tracers of AGN-dominated gas and determine the origin of the differences in such ratios across different types of gas. Such a determination will enable a more rigorous use of such ratios. First, we empirically examined the aforementioned ratios at different angular resolutions to quantify correlations. We then used local thermodynamic equilibrium (LTE) and non-LTE analyses coupled with Markov chain Monte Carlo (MCMC) sampling in order to determine the origin of the underlying differences in ratios. We propose that at high spatial resolution (< 50 pc) the HCN(4-3)/CS(2-1) is a reliable tracer of AGN activity. We also find that the variations in ratios are not a consequence of different densities or temperature but of different fractional abundances, yielding to the important result that it is essential to consider the chemical processes at play when drawing conclusions from radiative transfer calculations. From analyses at varying spatial scales, we find that previously proposed molecular line ratios, as well as a new one, have varying levels of consistency. We also determine from an investigation of radiative transfer modelling of our data that it is essential to consider the chemistry of the species when reaching conclusions from radiative transfer calculations., Comment: 15 pages, 13 figures, accepted for publication by Astronomy and Astrophysics on 09/09/2022

We study molecular outflows in a sample of 25 nearby (z< 0.17, d<750 Mpc) ULIRG systems (38 individual nuclei) as part of the "Physics of ULIRGs with MUSE and ALMA" (PUMA) survey, using ~400 pc (0.1-1.0" beam FWHM) resolution ALMA CO(2-1) observations. We used a spectro-astrometry analysis to identify high-velocity (> 300 km/s) molecular gas disconnected from the galaxy rotation, which we attribute to outflows. In 77% of the 26 nuclei with $\log L_{IR}/L_{\odot}>11.8$, we identifid molecular outflows with an average $v_{out}= 490$ km/s, outflow masses $1-35 \times 10^7$ $M_{\odot}$, mass outflow rates $\dot{M}_{out}=6-300$ $M_{\odot}$ yr$^{-1}$, mass-loading factors $\eta = \dot{M}_{out}/SFR = 0.1-1$, and an average outflow mass escape fraction of 45%. The majority of these outflows (18/20) are spatially resolved with radii of 0.2-0.9 kpc and have short dynamical times ($t_{dyn}=R_{out}/v_{out}$) in the range 0.5-2.8 Myr. The outflow detection rate is higher in nuclei dominated by starbursts (SBs, 14/15=93%) than in active galactic nuclei (AGN, 6/11=55%). Outflows perpendicular to the kinematic major axis are mainly found in interacting SBs. We also find that our sample does not follow the $\dot{M}_{out}$ versus AGN luminosity relation reported in previous works. In our analysis, we include a sample of nearby main-sequence galaxies (SFR = 0.3-17 $M_{\odot}$ yr$^{-1}$) with detected molecular outflows from the PHANGS-ALMA survey to increase the $L_{IR}$ dynamic range. Using these two samples, we find a correlation between the outflow velocity and the SFR, as traced by $L_{IR}$ ($v_{out} \propto SFR^{0.25\pm0.01})$, which is consistent with what was found for the atomic ionised and neutral phases. Using this correlation, and the relation between $M_{out}/R_{out}$ and $v_{out}$, we conclude that these outflows are likely momentum-driven., Comment: 27 pages, 23 figures; Appendix: 49 pages, 41 figures. Accepted for publication in Astronomy & Astrophysics. Typo corrected in Table 1 for galaxy 19297-0406