Showing 1–50 of 250 results for author: García-Burillo, S

GATOS: missing molecular gas in the outflow of NGC5728 revealed by JWST

Authors: R. Davies, T. Shimizu, M. Pereira-Santaella, A. Alonso-Herrero, A. Audibert, E. Bellocchi, P. Boorman, S. Campbell, Y. Cao, F. Combes, D. Delaney, T. Diaz-Santos, F. Eisenhauer, D. Esparza Arredondo, H. Feuchtgruber, N. M. Forster Schreiber, L. Fuller, P. Gandhi, I. Garcia-Bernete, S. Garcia-Burillo, B. Garcia-Lorenzo, R. Genzel, S. Gillessen, O. Gonzalez Martin, H. Haidar , et al. (27 additional authors not shown)

Abstract: 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 integr… ▽ More 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. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: A&A accepted; 16 pages

Deciphering the imprint of AGN feedback in Seyfert galaxies: Nuclear-scale molecular gas deficits

Authors: S. García-Burillo, E. K. S. Hicks, A. Alonso-Herrero, M. Pereira-Santaella, A. Usero, M. Querejeta, O. González-Martin, D. Delaney, C. Ramos Almeida, F. Combes, D. Anglés-Alcázar, A. Audibert, E. Bellocchi, R. I. Davies, T. A. Davis, J. S. Elford, I. García-Bernete, S. Hönig, A. Labiano, M. T. Leist, N. A. Levenson, E. López-Rodríguez, J. Mercedes-Feliz, C. Packham, C. Ricci , et al. (4 additional authors not shown)

Abstract: 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 ran… ▽ More 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}(Σ^{\rm gas}_{\rm 50}/Σ^{\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. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: Accepted by Astronomy & Astrophysics (A&A) (14/06/2024), 26 pages, 16 figures

Do spiral arms enhance star formation efficiency?

Authors: Miguel Querejeta, Adam K. Leroy, Sharon E. Meidt, Eva Schinnerer, Francesco Belfiore, Eric Emsellem, Ralf S. Klessen, Jiayi Sun, Mattia Sormani, Ivana Bešlic, Yixian Cao, Mélanie Chevance, Dario Colombo, Daniel A. Dale, Santiago García-Burillo, Simon C. O. Glover, Kathryn Grasha, Brent Groves, Eric. W. Koch, Lukas Neumann, Hsi-An Pan, Ismael Pessa, Jérôme Pety, Francesca Pinna, Lise Ramambason , et al. (10 additional authors not shown)

Abstract: 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 fr… ▽ More 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) △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: 26 pages, 16 figures. Accepted for publication in A&A

Extended high-ionization [MgIV] emission tracing widespread shocks in starbursts seen by JWST /NIRSpec

Authors: Miguel Pereira-Santaella, Ismael García-Bernete, Eduardo González-Alfonso, Almudena Alonso-Herrero, Luis Colina, Santiago García-Burillo, Dimitra Rigopoulou, Santiago Arribas, Michele Perna

Abstract: We report the detection of extended (>0.5-1kpc) high-ionization [MgIV] 4.487 $μ$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 o… ▽ More We report the detection of extended (>0.5-1kpc) high-ionization [MgIV] 4.487 $μ$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 $μ$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. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: Accepted for publication in A&A Letters. 6 pages, 6 figures. Appendix: 8 pages, 10 figures

Journal ref: A&A 685, L13 (2024)

Dust beyond the torus: Revealing the mid-infrared heart of local Seyfert ESO 428-G14 with JWST/MIRI

Authors: Houda Haidar, David J. Rosario, Almudena Alonso-Herrero, Miguel Pereira-Santaella, Ismael García-Bernete, Stephanie Campbell, Sebastian F. Hönig, Cristina Ramos Almeida, Erin Hicks, Daniel Delaney, Richard Davies, Claudio Ricci, Chris M. Harrison, Mason Leist, Enrique Lopez-Rodriguez, Santiago Garcia-Burillo, Lulu Zhang, Chris Packham, Poshak Gandhi, Anelise Audibert, Enrica Bellocchi, Peter Boorman, Andrew Bunker, Françoise Combes, Tanio Diaz Santos , et al. (12 additional authors not shown)

Abstract: 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 heart of ESO 428-G14, part of the GATOS surve… ▽ More 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 heart of ESO 428-G14, part of the GATOS survey of local AGN. We clearly 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. The morphology of the MIR structure bears a striking resemblance to that of gas ionised by the AGN in the narrow-line region (NLR). 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 using ancillary spectroscopy, we find the morphology of the dust continuum to be asymmetric around the nucleus and more compact, though still clearly extended out to r ~ 100 pc. We estimate a temperature of the emitting dust of ~ 120 K. Using simple models, we find that the heating of small dust grains (~ 0.01 microns) by the radiation from the central AGN and/or radiative jet-induced shocks is responsible for the extended MIR emission. Large-grained dust, predicted by models of 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. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 17 pages, 8 figures, 1 table. Submitted to MNRAS, comments are welcome!

Multi-transition study of methanol towards NGC 1068 with ALMA

Authors: K. -Y. Huang, D. Abbink, S. Viti, S. García-Burillo

Abstract: 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 pr… ▽ More 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. △ Less

Submitted 15 April, 2024; originally announced April 2024.

A spectacular galactic scale magnetohydrodynamic powered wind in ESO 320-G030

Authors: M. D. Gorski, S. Aalto, S. König, C. F. Wethers, C. Yang, S. Muller, K. Onishi, M. Sato, N. Falstad, Jeffrey G. Mangum, S. T. Linden, F. Combes, S. Martín, M. Imanishi, Keiichi Wada, L. Barcos-Muñoz, F. Stanley, S. García-Burillo, P. P. van der Werf, A. S. Evans, C. Henkel, S. Viti, N. Harada, T. Díaz-Santos, J. S. Gallagher , et al. (1 additional authors not shown)

Abstract: 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. Compa… ▽ More 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. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 9 pages, 10 figures

Journal ref: A&A 684, L11 (2024)

AGN feedback in the Local Universe: multiphase outflow of the Seyfert galaxy NGC 5506

Authors: Federico Esposito, Almudena Alonso-Herrero, Santiago García-Burillo, Viviana Casasola, Françoise Combes, Daniele Dallacasa, Richard Davies, Ismael García-Bernete, Begoña García-Lorenzo, Laura Hermosa Muñoz, Luis Peralta de Arriba, Miguel Pereira-Santaella, Francesca Pozzi, Cristina Ramos Almeida, Thomas Taro Shimizu, Livia Vallini, Enrica Bellocchi, Omaira González-Martín, Erin K. S. Hicks, Sebastian Hönig, Alvaro Labiano, Nancy A. Levenson, Claudio Ricci, David J. Rosario

Abstract: 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-Barol… ▽ More 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]$λ$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 △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: 24 pages, 30 figures, accepted for publication by Astronomy & Astrophysics

CON-quest II. Spatially and spectrally resolved HCN/HCO+ line ratios in local luminous and ultraluminous infrared galaxies

Authors: Y. Nishimura, S. Aalto, M. D. Gorski, S. König, K. Onishi, C. Wethers, C. Yang, L. Barcos-Muñoz, F. Combes, T. Díaz-Santos, J. S. Gallagher, S. García-Burillo, E. González-Alfonso, T. R. Greve, N. Harada, C. Henkel, M. Imanishi, K. Kohno, S. T. Linden, J. G. Mangum, S. Martín, S. Muller, G. C. Privon, C. Ricci, F. Stanley , et al. (2 additional authors not shown)

Abstract: 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 an… ▽ More 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. △ Less

Submitted 25 April, 2024; v1 submitted 23 February, 2024; originally announced February 2024.

Comments: 52 pages, 35 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 686, A48 (2024)

Deconvolution of JWST/MIRI Images: Applications to an AGN Model and GATOS Observations of NGC 5728

Authors: M. T. Leist, C. Packham, D. J. V. Rosario, D. A. Hope, A. Alonso-Herrero, E. K. S. Hicks, S. Hönig, L. Zhang, R. Davies, T. Díaz-Santos, O. Ganzález-Martín, E. Bellocchi, P. G. Boorman, F. Combes, I. García-Bernete, S. García-Burillo, B. García-Lorenzo, H. Haidar, K. Ichikawa, M. Imanishi, S. M. Jefferies, Á. Labiano, N. A. Levenson, R. Nikutta, M. Pereira-Santaella , et al. (10 additional authors not shown)

Abstract: 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… ▽ More 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ũ 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°) 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. △ Less

Submitted 14 February, 2024; v1 submitted 21 December, 2023; originally announced December 2023.

Comments: 32 pages, 23 figures, published in AJ 2024 February 7

Modelling molecular clouds and CO excitation in AGN-host galaxies

Authors: Federico Esposito, Livia Vallini, Francesca Pozzi, Viviana Casasola, Almudena Alonso-Herrero, Santiago García-Burillo, Roberto Decarli, Francesco Calura, Cristian Vignali, Matilde Mingozzi, Carlotta Gruppioni, Dhrubojyoti Sengupta

Abstract: 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… ▽ More 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 $α_{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 ($α_{CO}$, $N_H$) with a Markov chain Monte Carlo algorithm, finding values consistent with those present in the literature. We find a median $α_{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. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: 19 pages, 36 figures, resubmitted to MNRAS after addressing referee's comments

The Galaxy Activity, Torus, and Outflow Survey (GATOS). III: Revealing the inner icy structure in local AGN

Authors: I. García-Bernete, A. Alonso-Herrero, D. Rigopoulou, M. Pereira-Santaella, T. Shimizu, R. Davies, F. R. Donnan, P. F. Roche, O. González-Martín, C. Ramos Almeida, E. Bellocchi, P. Boorman, F. Combes, A. Efstathiou, D. Esparza-Arredondo, S. García-Burillo, E. González-Alfonso, E. K. S. Hicks, S. Hönig, A. Labiano, N. A. Levenson, E. López-Rodríguez, C. Ricci, C. Packham, D. Rouan , et al. (2 additional authors not shown)

Abstract: 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… ▽ More 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 $μ$m). All the nuclei show the 9.7 $μ$m silicate band in absorption. We compare the strength of the 9.7 and 18 $μ$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 $μ$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 $μ$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 $μ$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). △ Less

Submitted 4 December, 2023; v1 submitted 13 October, 2023; originally announced October 2023.

Comments: Accepted for publication in A&A. 13 pages, 12 Figures

The CO-to-H2 conversion factor of molecular outflows. Rovibrational CO emission in NGC 3256-S resolved by JWST/NIRSpec

Authors: M. Pereira-Santaella, E. González-Alfonso, I. García-Bernete, S. García-Burillo, D. Rigopoulou

Abstract: 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-bra… ▽ More 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). △ Less

Submitted 27 October, 2023; v1 submitted 12 September, 2023; originally announced September 2023.

Comments: Accepted for publication in A&A.17 pages, 13 figures

Star formation efficiency and AGN feedback in narrow-line Seyfert 1 galaxies with fast X-ray nuclear winds

Authors: Quentin Salomé, Yair Krongold, Anna Lia Longinotti, Manuela Bischetti, Santiago García-Burillo, Olga Vega, Miguel Sánchez-Portal, Chiara Feruglio, María Jesús Jiménez-Donaire, Maria Vittoria Zanchettin

Abstract: 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… ▽ More 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. △ Less

Submitted 24 July, 2023; v1 submitted 12 July, 2023; originally announced July 2023.

Comments: Published in MNRAS, Volume 524, Issue 2, Pages 3130-3145

Hydrodynamic simulations of the Disk of Gas Around Supermassive black holes (HDGAS) -I; Molecular Gas Dynamics

Authors: Mojtaba Raouf, Serena Viti, S. García-Burillo, Alexander J. Richings, Joop schaye, Ashley Bemis, Folkert S. J. Nobels, Matteo Guainazzi, Ko-Yun Huang, Matthieu Schaller, Violette Impellizzeri, Jon Holdship

Abstract: 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 outflow… ▽ More 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. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: 16 pages, 13 figures. Accepted for publication in MNRAS

AGN feedback in action in the molecular gas ring of the Seyfert galaxy NGC7172

Authors: A. Alonso-Herrero, S. Garcia-Burillo, M. Pereira-Santaella, T. Shimizu, F. Combes, E. K. S. Hicks, R. Davies, C. Ramos Almeida, I. Garcia-Bernete, S. F. Hoenig, N. A. Levenson, C. Packham, E. Bellocchi, L. K. Hunt, M. Imanishi, C. Ricci, P. Roche

Abstract: 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 w… ▽ More 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. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: Accepted for publication to A&A

The role of grain size in AGN torus dust models

Authors: Omaira González-Martín, Cristina Ramos Almeida, Jacopo Fritz, Almudena Alonso-Herrero, Sebastian F. Hönig, Patrick F. Roche, Donaji Esparza-Arredondo, Ismael García-Bernete, Santiago García-Burillo, Natalia Osorio-Clavijo, Ulises Reyes-Amador, Marko Stalevski, César Victoria-Ceballos

Abstract: 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 distribut… ▽ More 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. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Comments: 26 pages, 14 figures, 4 tables, accepted for publication in A&A

Journal ref: A&A 676, A73 (2023)

A radio-jet driven outflow in the Seyfert 2 galaxy NGC 2110?

Authors: L. Peralta de Arriba, A. Alonso-Herrero, S. García-Burillo, I. García-Bernete, M. Villar-Martín, B. García-Lorenzo, R. Davies, D. J. Rosario, S. F. Hönig, N. A. Levenson, C. Packham, C. Ramos Almeida, M. Pereira-Santaella, A. Audibert, E. Bellocchi, E. K. S. Hicks, A. Labiano, C. Ricci, D. Rigopoulou

Abstract: 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 Gauss… ▽ More 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 $σ\sim 60-200\ \mathrm{km\ s^{-1}}$. The outflow component has typical values of $σ\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]$λ$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. △ Less

Submitted 10 May, 2023; originally announced May 2023.

Comments: 25 pages, 22 figures, 3 tables. Accepted for publication in A&A

Gas, dust, and the CO-to-molecular gas conversion factor in low-metallicity starbursts

Authors: L. K. Hunt, F. Belfiore, F. Lelli, B. T. Draine, A. Marasco, S. Garcia-Burillo, G. Venturi, F. Combes, A. Weiß, C. Henkel, K. M. Menten, F. Annibali, V. Casasola, M. Cignoni, A. McLeod, M. Tosi, M. Beltran, A. Concas, G. Cresci, M. Ginolfi, N. Kumari, F. Mannucci

Abstract: 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… ▽ More 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. △ Less

Submitted 11 May, 2023; v1 submitted 5 May, 2023; originally announced May 2023.

Comments: A&A, in press (16 figures, 1 appendix)

Journal ref: A&A 675, A64 (2023)

Reconstructing the shock history in the CMZ of NGC 253 with ALCHEMI

Authors: K. -Y. Huang, S. Viti, J. Holdship, J. G. Mangum, S. Martín, N. Harada, S. Muller, K. Sakamoto, K. Tanaka, Y. Yoshimura, R. Herrero-Illana, D. S. Meier, E. Behrens, P. P. van der Werf, C. Henkel, S. García-Burillo, V. M. Rivilla, K. L. Emig, L. Colzi, P. K. Humire, R. Aladro, M. Bouvier

Abstract: 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 charact… ▽ More 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. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Journal ref: A&A 675, A151 (2023)

Young Radio Sources Expanding in Gas-Rich ISM: Using Cold Molecular Gas to Trace Their Impact

Authors: Raffaella Morganti, Suma Murthy, Pierre Guillard, Tom Oosterloo, Santiago Garcia-Burillo

Abstract: 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 absor… ▽ More 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. △ Less

Submitted 27 February, 2023; originally announced February 2023.

Comments: 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

Jet-induced molecular gas excitation and turbulence in the Teacup

Authors: A. Audibert, C. Ramos Almeida, S. García-Burillo, F. Combes, M. Bischetti, M. Meenakshi, D. Mukherjee, G. Bicknell, A. Y. Wagner

Abstract: 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… ▽ More 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. △ Less

Submitted 27 February, 2023; originally announced February 2023.

Comments: 9 pages, 8 figures. Accepted for publication in A&A Letters

Journal ref: A&A 671, L12 (2023)

NOEMA spatially resolved view of the multi-phase outflow in IRAS17020+4544: a shocked wind in action?

Authors: Anna Lia Longinotti, Quentin Salomé, Chiara Feruglio, Yair Krongold, Santiago García-Burillo, Marcello Giroletti, Francesca Panessa, Carlo Stanghellini, Olga Vega, Victor Manuel Patiño-Álvarez, Vahram Chavushyan, Mauricio Elías-Chavez, Aitor Robleto-Orús

Abstract: 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… ▽ More 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. △ Less

Submitted 16 February, 2023; originally announced February 2023.

Comments: Accepted for publication on MNRAS, 17 pages, 7 figures

Sub-kiloparsec empirical relations and excitation conditions of HCN and HCO+ J=3-2 in nearby star-forming galaxies

Authors: Axel Garcia-Rodriguez, Antonio Usero, Adam K. Leroy, Frank Bigiel, Maria Jesus Jimenez-Donaire, Daizhong Liu, Miguel Querejeta, Toshiki Saito, Eva Schinnerer, Ashley Barnes, Francesco Belfiore, Ivana Beslic, Yixian Cao, Melanie Chevance, Daniel A. Dale, Jakob S. den Brok, Cosima Eibensteiner, Santiago Garcia-Burillo, Simon C. O. Glover, Ralf S. Klessen, Jerome Pety, Johannes Puschnig, Erik Rosolowsky, Karin Sandstrom, Mattia C. Sormani , et al. (2 additional authors not shown)

Abstract: 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 avail… ▽ More 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 ($Σ_\textrm{star}$), the intensity-weighted CO(2-1) ($\langle I_\text{CO}\rangle$), and the star formation rate surface density ($Σ_\text{SFR}$). We find no apparent correlation with $Σ_\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). △ Less

Submitted 1 February, 2023; originally announced February 2023.

Comments: Accepted for publication in A&A. 14 pages, 8 figures

Journal ref: A&A 672, A96 (2023)

Understanding if molecular ratios can be used as diagnostics of AGN and starburst activity: The case of NGC 1068

Authors: J. Butterworth, J. Holdship, S. Viti, S. García-Burillo

Abstract: 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 inv… ▽ More 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. △ Less

Submitted 21 September, 2022; v1 submitted 13 September, 2022; originally announced September 2022.

Comments: 15 pages, 13 figures, accepted for publication by Astronomy and Astrophysics on 09/09/2022

Journal ref: A&A 667, A131 (2022)

Physics of ULIRGs with MUSE and ALMA: PUMA IV. No tight relation between cold molecular outflow rates and AGN luminosities

Authors: I. Lamperti, M. Pereira-Santaella, M. Perna, L. Colina, S. Arribas, S. García-Burillo, E. González-Alfonso, S. Aalto, A. Alonso-Herrero, F. Combes, A. Labiano, J. Piqueras-López, D. Rigopoulou, P. van der Werf

Abstract: 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 attribut… ▽ More 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 $η= \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. △ Less

Submitted 18 March, 2023; v1 submitted 7 September, 2022; originally announced September 2022.

Comments: 27 pages, 23 figures; Appendix: 49 pages, 41 figures. Accepted for publication in Astronomy & Astrophysics. Typo corrected in Table 1 for galaxy 19297-0406

Journal ref: A&A 668, A45 (2022)

A high angular resolution view of the PAH emission in Seyfert galaxies using JWST/MRS data

Authors: I. García-Bernete, D. Rigopoulou, A. Alonso-Herrero, F. R. Donnan, P. F. Roche, M. Pereira-Santaella, A. Labiano, L. Peralta de Arriba, T. Izumi, C. Ramos Almeida, T. Shimizu, S. Hönig, S. García-Burillo, D. J. Rosario, M. J. Ward, E. Bellocchi, E. K. S. Hicks, L. Fuller, C. Packham

Abstract: Polycyclic Aromatic Hydrocarbons (PAHs) are carbon-based molecules that are ubiquitous in a variety of astrophysical objects and environments. In this work, we use JWST/MIRI MRS spectroscopy of three Seyferts to compare their nuclear PAH emission with that of star-forming regions. This study represents the first of its kind using sub-arcsecond angular resolution data of local luminous Seyferts (Lb… ▽ More Polycyclic Aromatic Hydrocarbons (PAHs) are carbon-based molecules that are ubiquitous in a variety of astrophysical objects and environments. In this work, we use JWST/MIRI MRS spectroscopy of three Seyferts to compare their nuclear PAH emission with that of star-forming regions. This study represents the first of its kind using sub-arcsecond angular resolution data of local luminous Seyferts (Lbol>10^44.46 erg/s) on a wide wavelength coverage (4.9-28.1 micron). We present an analysis of their nuclear PAH properties by comparing the observed ratios with PAH diagnostic model grids, derived from theoretical spectra. Our results show that a suite of PAH features is present in the innermost parts (~0.45 arcsec at 12 micron; in the inner ~142-245 pc) of luminous Seyfert galaxies. We find that the nuclear regions of AGN lie at different positions of the PAH diagnostic diagrams, whereas the SF regions are concentrated around the average values of SF galaxies. In particular, we find that the nuclear PAH emission mainly originates in neutral PAHs. In contrast, PAH emission originating in the SF regions favours ionised PAH grains. The observed PAH ratios in the nuclear region of AGN-dominated galaxy NGC 6552 indicate the presence of larger-sized PAH molecules compared with those of the SF regions. Therefore, our results provide evidence that the AGN have a significant impact on the ionization state (and probably the size) of the PAH grains on scales of ~142-245 pc. △ Less

Submitted 28 September, 2022; v1 submitted 24 August, 2022; originally announced August 2022.

Comments: 12 pages, 10 figures. Small changes to match the accepted for publication version in A&A (official date of acceptance: 20/09/2022)

Journal ref: A&A 666, L5 (2022)

Low-power jet-ISM interaction in NGC 7319 revealed by JWST/MIRI MRS

Authors: M. Pereira-Santaella, J. Álvarez-Márquez, I. García-Bernete, A. Labiano, L. Colina, A. Alonso-Herrero, E. Bellocchi, S. García-Burillo, S. F. Hönig, C. Ramos Almeida, D. Rosario

Abstract: We present JWST/MIRI MRS spectroscopy of NGC7319, the largest galaxy in the Stephan's Quintet, observed as part of the Early Release Observations (ERO). NGC7319 hosts a type 2 active galactic nucleus (AGN) and a low-power radio jet (L_1.4GHz=3.3x10^22 W Hz^-1) with two asymmetric radio hotspots at 430 pc (N2) and 1.5 kpc (S2) projected distances from the unresolved radio core. The MRS data suggest… ▽ More We present JWST/MIRI MRS spectroscopy of NGC7319, the largest galaxy in the Stephan's Quintet, observed as part of the Early Release Observations (ERO). NGC7319 hosts a type 2 active galactic nucleus (AGN) and a low-power radio jet (L_1.4GHz=3.3x10^22 W Hz^-1) with two asymmetric radio hotspots at 430 pc (N2) and 1.5 kpc (S2) projected distances from the unresolved radio core. The MRS data suggest that the molecular material in the disk of the galaxy decelerates the jet and causes this length asymmetry. We find enhanced emission from warm and hot H_2 (T_w=330+-40 K, T_h=900+-60 K) and ionized gas at the intersection between the jet axis and dust lanes in the disk. This emission is coincident with the radio hotspot N2, the hotspot closer to the core, suggesting that the jet-interstellar medium (ISM) interaction decelerates the jet. Conversely, the mid-infrared emission at the more distant hotspot is fainter, more highly ionized, and with lower H_2 excitation, suggesting a more diffuse atomic environment where the jet can progress to farther distances. At the N2 radio hotspot, the ionized gas mass (M_ion=(2.4-12)x10^5 Msun) is comparable to that of the warm H_2, but the former is more turbulent (sigma_ion~300 vs. sigma_H2~150 km/s), so the mechanical energy of the ionized gas is ~1.3-10 times higher. From these estimates, we find that only <1% of the jet energy remains as mechanical energy in these two ISM phases at N2. We also find extended (r>0.3-1.5 kpc) high-ionization emission ([MgV], [NeVI], and [NeV]) close to the radio hotspots. This initial analysis of NGC7319 shows the potential of MIRI/MRS to investigate the AGN feedback mechanisms due to radio jets and their radiation field in the, often heavily dust-enshrouded, central regions of galaxies. Understanding these mechanisms is an essential ingredient in the development of cosmological simulations of galaxy evolution. △ Less

Submitted 16 September, 2022; v1 submitted 9 August, 2022; originally announced August 2022.

Comments: 10 pages, 8 figures. Small changes to match the accepted for publication version in A&A

Journal ref: A&A 665, L11 (2022)

From the Circumnuclear Disk in the Galactic Center to thick, obscuring tori of AGNs -- Modelling the molecular emission of a parsec-scale torus as found in NGC1068

Authors: B. Vollmer, R. I. Davies, P. Gratier, Th. Lizee, M. Imanishi, J. F. Gallimore, C. M. V. Impellizzeri, S. Garcia-Burillo, F. Le Petit

Abstract: The accretion rates needed to fuel the central black hole in a galaxy can be achieved via viscous torques in thick disks and rings, which can be resolved by millimetre interferometry within the inner ~20pc of the active galaxy NGC1068 at comparable scales and sensitivity to single dish observations of the Circumnuclear Disk (CND) in the Galactic Center. To interpret observations of these regions a… ▽ More The accretion rates needed to fuel the central black hole in a galaxy can be achieved via viscous torques in thick disks and rings, which can be resolved by millimetre interferometry within the inner ~20pc of the active galaxy NGC1068 at comparable scales and sensitivity to single dish observations of the Circumnuclear Disk (CND) in the Galactic Center. To interpret observations of these regions and determine the physical properties of their gas distribution, we present a modelling effort that includes (i) a simple dynamical simulations involving partially inelastic collisions between disk gas clouds, (ii) an analytical model of a turbulent clumpy gas disk calibrated by the dynamical model and observations, (iii) local turbulent and cosmic ray gas heating and cooling via H2O, H2, and CO emission, and (iv) determination of the molecular abundances. We also consider photodissociation regions (PDR) where gas is directly illuminated by the central engine. We compare the resulting model datacubes of the CO, HCN, HCO+, and CS brightness temperatures to available observations. In both cases the kinematics can be explained by one or two clouds colliding with a pre-existing ring, in a prograde sense for the CND and retrograde for NGC1068. And, with only dense disk clouds, the line fluxes can be reproduced to within a factor of about two. To avoid self-absorption of the intercloud medium, turbulent heating at the largest scales, comparable to the disk height, has to be decreased by a factor of 50-200. Our models indicate that turbulent mechanical energy input is the dominant gas heating mechanism within the thick gas disks. In N1068, while the bulk of the AGN X-ray radiation is absorbed in a layer of Compton-thick gas inside the dust sublimation radius, the optical/UV radiation may enhance the molecular line emission from photodissociation regions by ~50% at the inner edge of the gas ring. △ Less

Submitted 29 June, 2022; originally announced June 2022.

Comments: 36 pages, 32 figures, accepted for publication in A&A

Journal ref: A&A 665, A102 (2022)

ALMA imaging of the cold molecular and dusty disk in the type 2 active nucleus of the Circinus galaxy

Authors: Konrad R. W. Tristram, C. M. Violette Impellizzeri, Zhi-Yu Zhang, Eric Villard, Christian Henkel, Serena Viti, Leonard Burtscher, Françoise Combes, Santiago García-Burillo, Sergio Martín, Klaus Meisenheimer, Paul P. van der Werf

Abstract: We aim to shed light on the physical properties and kinematics of the molecular material in the nucleus of one of the closest type 2 active galaxies. To this end, we obtained high angular resolution ALMA observations of the nucleus of the Circinus galaxy. The observations map the emission at 350GHz and 690GHz with spatial resolutions of ~3.8pc and ~2.2pc, respectively. The continuum emission t… ▽ More We aim to shed light on the physical properties and kinematics of the molecular material in the nucleus of one of the closest type 2 active galaxies. To this end, we obtained high angular resolution ALMA observations of the nucleus of the Circinus galaxy. The observations map the emission at 350GHz and 690GHz with spatial resolutions of ~3.8pc and ~2.2pc, respectively. The continuum emission traces cold ($T\lesssim100$K) dust in a circumnuclear disk with spiral arms on scales of 25pc, plus a marginally resolved nuclear emission peak. The latter is not extended in polar direction as claimed based on earlier ALMA observations. A significant amount (of the order of 40%) of the 350GHz emission is not related to dust, but most likely free-free emission instead. We detect CO(3-2) and CO(6-5) as well as HCO$^+$(4-3), HCN(4-3), and CS(4-3). The CO emission is extended, showing a spiral pattern, similar to the extended dust emission. Towards the nucleus, CO is excited to higher transitions and its emission is self-absorbed, leading to an apparent hole in the CO(3-2) but not the CO(6-5) emission. On the other hand, the high gas density tracers HCO$^+$, HCN, and CS show a strong, yet unresolved (($\lesssim4$pc) concentration of the emission at the nucleus, pointing at a very small 'torus'. The kinematics are dominated by rotation and point at a geometrically thin disk down to the resolution limit of our observations. In contrast to several other AGNs, no HCN enhancement is found towards the nucleus. The Circinus nucleus is therefore composed of at least two distinct components: (1) an optically thin, warm outflow of ionised gas containing clouds of dust; and (2) a cold molecular and dusty disk. These findings support the most recent radiative transfer calculations of the obscuring structures in AGNs, which find a similar two-component structure. (Abridged) △ Less

Submitted 1 June, 2022; originally announced June 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 664, A142 (2022)

Compact molecular gas emission in local LIRGs among low- and high-z galaxies

Authors: E. Bellocchi, M. Pereira-Santaella, L. Colina, A. Labiano, M. Sánchez-García, A. Alonso-Herrero, S. Arribas, S. García-Burillo, M. Villar-Martín, D. Rigopoulou, F. Valentino, A. Puglisi, T. Díaz-Santos, S. Cazzoli, A. Usero

Abstract: We present new CO(2-1) observations of a representative sample of 24 local (z$<$0.02) luminous infrared galaxies (LIRGs) obtained at high spatial resolution ($<$100 pc) from ALMA. We derive the effective radii of the CO(2-1) and the 1.3 mm continuum emissions using the curve-of-growth method. LIRGs show an extremely compact cold molecular gas distribution (median R$_{CO}$ $\sim$0.7 kpc), which is… ▽ More We present new CO(2-1) observations of a representative sample of 24 local (z$<$0.02) luminous infrared galaxies (LIRGs) obtained at high spatial resolution ($<$100 pc) from ALMA. We derive the effective radii of the CO(2-1) and the 1.3 mm continuum emissions using the curve-of-growth method. LIRGs show an extremely compact cold molecular gas distribution (median R$_{CO}$ $\sim$0.7 kpc), which is a factor 2 smaller than the ionized gas, and 3.5 times smaller than the stellar size. The molecular size of LIRGs is similar to that of early-type galaxies (R$_{CO}\sim$1 kpc) and about a factor of 6 more compact than local Spirals of similar stellar mass. Only the CO emission in low-z ULIRGs is more compact than these local LIRGs by a factor of 2. Compared to high-z (1$<$z$<$6) systems, the stellar sizes and masses of local LIRGs are similar to those of high-z MS star-forming galaxies (SFG) and about a factor of 2-3 lower than sub-mm galaxies (SMG). The molecular sizes of high-z MS SFGs and SMGs are larger than those derived for LIRGs by a factor of $\sim$3 and $\sim$8, respectively. These results indicate that while low-z LIRGs and high-z MS-SFGs have similar stellar masses and sizes, the regions of current star formation (traced by the ionized gas) and of potential star-formation (traced by the molecular gas) are substantially smaller in LIRGs, and constrained to the central kpc region. High-z galaxies represent a wider population but their star-forming regions are more extended, even covering the overall size of the host galaxy. High-z galaxies have larger fractions of gas than low-z LIRGs, and therefore the formation of stars could be induced by interactions and mergers in extended disks or filaments with large enough molecular gas surface density involving physical mechanisms similar to those identified in the central kpc of LIRGs. △ Less

Submitted 5 April, 2022; originally announced April 2022.

Comments: 36 pages, 20 figures, 8 tables. Accepted for publication in A&A

Journal ref: A&A 664, A60 (2022)

Unveiling the Main Sequence to Starburst Transition Region with a Sample of Intermediate Redshift Luminous Infrared Galaxies

Authors: L. Hogan, D. Rigopoulou, S. García-Burillo, A. Alonso-Herrero, L. Barrufet, F. Combes, I. García-Bernete, G. E. Magdis, M. Pereira-Santaella, N. Thatte, A. Weiß

Abstract: We present a CO(3-2) study of four systems composed of six (ultra) luminous infrared galaxies (U/LIRGs), located at 0.28 < z < 0.44, that straddle the transition region between regular star forming galaxies and starbursts. These galaxies benefit from previous multi-wavelength analysis allowing in depth exploration of an understudied population of U/LIRGs at a time when the universe is experiencing… ▽ More We present a CO(3-2) study of four systems composed of six (ultra) luminous infrared galaxies (U/LIRGs), located at 0.28 < z < 0.44, that straddle the transition region between regular star forming galaxies and starbursts. These galaxies benefit from previous multi-wavelength analysis allowing in depth exploration of an understudied population of U/LIRGs at a time when the universe is experiencing a rapid decline in star formation rate density. We detect CO(3-2) emission in four targets and these galaxies fall between the loci of regular star forming galaxies and starbursts on the Kennicutt-Schmidtt relation. Compared to low luminosity LIRGs and high luminosity ULIRGs at similar redshifts, we find they all have similar molecular gas budgets with the difference in their star formation rates (SFR) driven by the star formation efficiency (SFE). This suggests that at these redshifts large molecular gas reservoirs must coincide with an increased SFE to transition a galaxy into the starburst regime. We studied the structure and kinematics and found our four detections are either interacting or have disturbed morphology which may be driving the SFE. One of the CO(3-2) non-detections has a strong continuum detection, and has been previously observed in H$α$, suggesting an unusual interstellar medium for a ULIRG. We conclude that our sample of transitioning U/LIRGs fill the gap between regular star forming galaxies and starbursts, suggest a continuous change in SFE between these two populations and the increased SFE may be driven by morphology and differing stages of interaction. △ Less

Submitted 21 February, 2022; originally announced February 2022.

Comments: 19 pages, 17 figures, accepted for publication in MNRAS

The chemical footprint of AGN feedback in the outflowing circumnuclear disk of NGC 1068

Authors: K. -Y. Huang, S. Viti, J. Holdship, S. García-Burillo, K. Kohno, A. Taniguchi, S. Martín, R. Aladro, A. Fuente, M. Sánchez-García

Abstract: In the nearby (D=14 Mpc) AGN-starburst composite galaxy NGC 1068, it has been found that the molecular gas in the Circum-nuclear Disk (CND) is outflowing, which is a manifestation of ongoing AGN feedback. The outflowing gas has a large spread of velocities, which likely drive different shock chemistry signatures at different locations in the CND. We perform a multi-line molecular study using two s… ▽ More In the nearby (D=14 Mpc) AGN-starburst composite galaxy NGC 1068, it has been found that the molecular gas in the Circum-nuclear Disk (CND) is outflowing, which is a manifestation of ongoing AGN feedback. The outflowing gas has a large spread of velocities, which likely drive different shock chemistry signatures at different locations in the CND. We perform a multi-line molecular study using two shock tracers, SiO and HNCO, with the aim to determine the gas properties traced by these two species, and explore the possibility of reconstructing the shock history in the CND. Five SiO transitions and three HNCO transitions were imaged at high resolution $0''.5-0''.8$ with the Atacama Large Millimeter/submillimeter Array (ALMA). We performed both LTE and non-LTE radiative transfer analysis coupled with Bayesian inference process in order to characterize the gas properties, such as molecular gas density and gas temperature. We found clear evidence of chemical differentiation between SiO and HNCO, with the SiO/HNCO ratio ranging from greater than one on the east of CND to lower than one on the west side. The non-LTE radiative transfer analysis coupled with Bayesian inference confirms that the gas traced by SiO has different densities - and possibly temperatures - than that traced by HNCO. We find that SiO traces gas affected by fast shocks while the gas traced by HNCO is either just affected by slow shocks or not shocked at all. A distinct differentiation between SiO and HNCO has been revealed in our observations and the further analysis of the gas properties traced by both species, which confirms the results from previous chemical modelings. △ Less

Submitted 10 February, 2022; originally announced February 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 666, A102 (2022)

The complex infrared dust continuum emission of NGC1068: ground-based N- and Q-band spectroscopy and new radiative transfer models

Authors: César Ivan Victoria-Ceballos, Omaira González-Martín, Jacopo Fritz, Cristina Ramos Almeida, Enrique López-Rodríguez, Santiago García-Burillo, Almudena Alonso-Herrero, Mariela Martínez-Paredes, Donaji Esparza-Arredondo, Natalia Osorio-Clavijo

Abstract: Thanks to ground-based infrared and sub-mm observations the study of the dusty torus of nearby AGN has greatly advanced in the last years. With the aim of further investigating the nuclear mid-infrared emission of the archetypal Seyfert 2 galaxy NGC1068, here we present a fitting to the N- and Q-band Michelle/Gemini spectra. We initially test several available SED libraries, including a smooth, cl… ▽ More Thanks to ground-based infrared and sub-mm observations the study of the dusty torus of nearby AGN has greatly advanced in the last years. With the aim of further investigating the nuclear mid-infrared emission of the archetypal Seyfert 2 galaxy NGC1068, here we present a fitting to the N- and Q-band Michelle/Gemini spectra. We initially test several available SED libraries, including a smooth, clumpy and two phase torus models, and a clumpy disk plus wind model. Although, the smooth torus model describe the spectra of NGC1068 if we allow to vary some model parameters among the two spectral bands. Motivated by this result, we produced new SEDs using the radiative transfer code SKIRT. We use two concentric tori that allow us to test a more complex geometry. We test different values for the inner and outer radii, half opening angle, radial and polar exponent of the power-law density profile, opacity, and viewing angle. Furthermore, we also test the dust grains size and different optical and calorimetric properties of silicate grains. The best fitting model consists of two concentric components with outer radii of 1.8 and 28 pc, respectively. We find that the size and the optical and calorimetric properties of graphite and silicate grains in the dust structure are key to reproduce the spectra of NGC1068. We conclude that the dust in NGC1068 reaches different scales, where the highest contribution to the mid-infrared is given by a central and compact component. A less dense and extended component is present, which can be either part of the same torus (conforming a flared disk) or can represent the emission of a polar dust component, as already suggested from interferometric observations. △ Less

Submitted 27 January, 2022; originally announced January 2022.

Spatially resolved star-formation relations of dense molecular gas in NGC 1068

Authors: M. Sánchez-García, S. García-Burillo, M. Pereira-Santaella, L. Colina, A. Usero, M. Querejeta, A. Alonso-Herrero, A. Fuente

Abstract: We analyse the influence of the dynamical environment on the star formation (SF) relations of the dense molecular gas in the starburst (SB) ring of the Seyfert 2 galaxy NGC 1068. We used ALMA to image the emission of the 1-0 transitions of HCN and HCO+ with a resolution of 56 pc. We also used ancillary data of CO(1-0) at a resolution of ~100 pc, and CO(3-2) and its underlying continuum emission at… ▽ More We analyse the influence of the dynamical environment on the star formation (SF) relations of the dense molecular gas in the starburst (SB) ring of the Seyfert 2 galaxy NGC 1068. We used ALMA to image the emission of the 1-0 transitions of HCN and HCO+ with a resolution of 56 pc. We also used ancillary data of CO(1-0) at a resolution of ~100 pc, and CO(3-2) and its underlying continuum emission at ~40 pc. These observations allow us to probe a range of molecular gas densities (n(H2)~10$^{3-5}cm^{-3}$). The SF rate (SFR) is derived from Pa$α$ line emission imaged by HST/NICMOS. We analysed how SF relations change depending on the choice of aperture sizes and molecular gas tracer. The scatter in the Kennicutt-Schmidt relation is about a factor of two to three lower for the HCN and HCO+ lines compared to CO(1-0) for a common aperture. Correlations lose statistical significance below a critical spatial scale $\approx$300-400 pc. The SF efficiency of the dense molecular gas (SFEdense) shows a scattered distribution as a function of the HCN luminosity (L'(HCN)) around a mean value of $\simeq0.01$Myr$^{-1}$. An alternative prescription for SF relations, linking the SFEdense and the boundedness of the gas measured by the parameter b$\equivΣ$dense/$σ^2$, where $Σ$dense is the dense molecular gas surface density and $σ$ the velocity dispersion, resolves the degeneracy associated with the SFEdense-L'(HCN) plot. We identify two branches in the SFEdense-b plot that correspond to two dynamical environments in the SB ring, which are defined by their proximity to the bar-ring interface region. This region corresponds to the crossing of two density wave resonances, where an increased rate of cloud-cloud collisions would favour an enhanced compression of molecular gas. Our results suggest that galactic dynamics plays a major role in the efficiency of the gas conversion into stars. △ Less

Submitted 17 January, 2022; originally announced January 2022.

Comments: 30 pages, 34 figures, 4 tables. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 660, A83 (2022)

The diverse cold molecular gas contents, morphologies, and kinematics of type-2 quasars as seen by ALMA

Authors: C. Ramos Almeida, M. Bischetti, S. Garcia-Burillo, A. Alonso-Herrero, A. Audibert, C. Cicone, C. Feruglio, C. N. Tadhunter, J. C. S. Pierce, M. Pereira-Santaella, P. S. Bessiere

Abstract: We present CO(2-1) and adjacent continuum observations of 7 nearby radio-quiet type-2 quasars (QSO2s) obtained with ALMA at ~0.2" resolution (370 pc at z~0.1). The CO morphologies are diverse, including disks and interacting systems. Two of the QSO2s are red early-type galaxies with no CO(2-1) detected. In the interacting galaxies, the central kpc contains 18-25% of the total cold molecular gas, w… ▽ More We present CO(2-1) and adjacent continuum observations of 7 nearby radio-quiet type-2 quasars (QSO2s) obtained with ALMA at ~0.2" resolution (370 pc at z~0.1). The CO morphologies are diverse, including disks and interacting systems. Two of the QSO2s are red early-type galaxies with no CO(2-1) detected. In the interacting galaxies, the central kpc contains 18-25% of the total cold molecular gas, whereas in the spirals it is only 5-12%. J1010+0612 and J1430+1339 show double-peaked CO morphologies which do not have optical counterparts. Based on our analysis of the ionized and molecular kinematics and mm continuum emission, these CO morphologies are most likely produced by AGN feedback in the form of outflows, jets, and/or shocks. The CO kinematics of the QSO2s are dominated by rotation but also reveal noncircular motions. According to our analysis of the kinematics, these noncircular motions correspond to molecular outflows mostly coplanar with the CO discs in four of the QSO2s, and either to a coplanar inflow or vertical outflow in the case of J1010+0612. These outflows represent 0.2-0.7% of the QSO2s' total molecular gas mass and have maximum velocities of 200-350 km/s, radii from 0.4 to 1.3 kpc, and outflow rates of 8-16 Msun/yr. These properties are intermediate between those of the mild molecular outflows measured for Seyferts, and the fast and energetic outflows of ULIRGs. This suggests that it is not only AGN luminosity that drives massive molecular outflows. Other factors such as jet power, coupling between winds, jets, and/or ionized outflows and the CO discs, and amount or geometry of dense gas in the nuclear regions might be also relevant. Thus, although we do not find evidence for a significant impact of quasar feedback on the total molecular gas reservoirs and SFRs, it appears to be modifying the distribution of cold molecular gas in the central kpc of the galaxies. △ Less

Submitted 15 March, 2022; v1 submitted 26 November, 2021; originally announced November 2021.

Comments: 29 pages including 18 figures, 7 tables and appendices. A&A accepted version including minor edits

Journal ref: A&A 658, A155 (2022)

Duality in spatially resolved star-formation relations in local LIRGs

Authors: M. Sánchez-García, M. Pereira-Santaella, S. García-Burillo, L. Colina, A. Alonso-Herrero, M. Villar-Martín, T. Saito, T. Díaz-Santos, J. Piqueras López, S. Arribas, E. Bellocchi, S. Cazzoli, A. Labiano

Abstract: We analyse the star formation (SF) relations in a sample of 16 nearby luminous infrared galaxies (LIRGs) with more than 2800 regions defined on scales of 90 to 500 pc. We used ALMA to map the distribution of the cold molecular gas traced by the J = 2-1 line of CO and archival PaαHST/NICMOS imaging to trace the recent SF. In four objects, we find two different branches in the Kennicutt-Schmidt rela… ▽ More We analyse the star formation (SF) relations in a sample of 16 nearby luminous infrared galaxies (LIRGs) with more than 2800 regions defined on scales of 90 to 500 pc. We used ALMA to map the distribution of the cold molecular gas traced by the J = 2-1 line of CO and archival PaαHST/NICMOS imaging to trace the recent SF. In four objects, we find two different branches in the Kennicutt-Schmidt relation at 90 pc scales, suggesting the existence of a duality in this relation. The two branches correspond to two different dynamical environments within each galaxy. One branch, which corresponds to the central region of these galaxies 90% of the regions are located at radii < 0.85 kpc), shows higher gas and star formation rate surface densities with higher velocity dispersion. The other branch, which shows lower molecular gas and SF rate surface densities, corresponds to the more external disk regions (r $\sim$ 1 kpc). Despite the scatter, the SF efficiency of the galaxies with a dual behaviour increases with increasing boundedness as measured by the b parameter (b $\equiv$ $Σ_{H_{2}}$/$σ^{2}$ $\propto$ $α_{vir}^{-1}$). At larger spatial scales (250 and 500 pc), the duality disappears. The rest of the sample does not show evidence of this dual behaviour at any scale. △ Less

Submitted 6 November, 2021; originally announced November 2021.

Comments: Accepted for publication in A&A, 28 pages, 10 figures, 5 tables

Journal ref: A&A 659, A102 (2022)

ALCHEMI: an ALMA Comprehensive High-resolution Extragalactic Molecular Inventory. Survey presentation and first results from the ACA array

Authors: S. Martín, J. G. Mangum, N. Harada, F. Costagliola, K. Sakamoto, S. Muller, R. Aladro, K. Tanaka, Y. Yoshimura, K. Nakanishi, R. Herrero-Illana, S. Mühle, S. Aalto, E. Behrens, L. Colzi, K. L. Emig, G. A. Fuller, S. García-Burillo, T. R. Greve, C. Henkel, J. Holdship, P. Humire, L. Hunt, T. Izumi, K. Kohno , et al. (10 additional authors not shown)

Abstract: We used the Atacama Large Millimeter/submillimeter Array (ALMA), covering a nearly contiguous 289 GHz frequency range between 84.2 and 373.2 GHz, to image the continuum and spectral line emission at 1.6\arcsec ($\sim 28$ pc) resolution down to a sensitivity of $30-50$ mK. This article describes the ALMA Comprehensive High-resolution Extragalactic Molecular Inventory (ALCHEMI) Large Program. We foc… ▽ More We used the Atacama Large Millimeter/submillimeter Array (ALMA), covering a nearly contiguous 289 GHz frequency range between 84.2 and 373.2 GHz, to image the continuum and spectral line emission at 1.6\arcsec ($\sim 28$ pc) resolution down to a sensitivity of $30-50$ mK. This article describes the ALMA Comprehensive High-resolution Extragalactic Molecular Inventory (ALCHEMI) Large Program. We focus on the analysis of the spectra extracted from the $15''$ ($\sim255$ pc) resolution ALMA Compact Array data. We model the molecular emission assuming local thermodynamic equilibrium with 78 species detected. Additionally, multiple hydrogen and helium recombination lines are identified. Spectral lines contribute 5 to 36\% of the total emission in frequency bins of 50 GHz. We report the first extragalactic detections of C$_2$H$_5$OH, HOCN, HC$_3$HO, and several rare isotopologues. Isotopic ratios of carbon, oxygen, sulfur, nitrogen and silicon were measure with multiple species. Infrared pumped vibrationaly excited HCN, HNC, and HC$_3$N emission, originating in massive star formation locations, is clearly detected at low resolution, while we do not detect it for HCO$^+$. We suggest high temperature conditions in these regions driving a seemingly "carbon-rich" chemistry which may also explain the observed high abundance of organic species close to those in Galactic hot cores. The $L_{vib}/L_{IR}$ ratio is used as a proxy to estimate a $3\%$ contribution from proto super star cluster to the global infrared emission. Measured isotopic ratios with high dipole moment species agree with those within the central kiloparsec of the Galaxy, while those derived from $\rm^{13}C^{18}O$ are a factor of 5 larger, confirming the existence of multiple ISM components within NGC 253 with different degrees of nucleosynthesis enrichment. ALCHEMI provides a template for early Universe galaxies. △ Less

Submitted 17 September, 2021; originally announced September 2021.

Comments: 60 pages, 32 figures, 15 tables, Accepted for publication at A&A

Journal ref: A&A 656, A46 (2021)

Kinetic temperature of massive star-forming molecular clumps measured with formaldehyde IV. The ALMA view of N113 and N159W in the LMC

Authors: X. D. Tang, C. Henkel, K. M. Menten, Y. Gong, C. -H. R. Chen, D. L. Li, M. -Y. Lee, J. G. Mangum, Y. P. Ao, S. Mühle, S. Aalto, S. García-Burillo, S. Martín, S. Viti, S. Muller, F. Costagliola, H. Asiri, S. A. Levshakov, M. Spaans, J. Ott, C. M. V. Impellizzeri, Y. Fukui, Y. X. He, J. Esimbek, J. J. Zhou , et al. (3 additional authors not shown)

Abstract: We mapped the kinetic temperature structure of two massive star-forming regions, N113 and N159W, in the Large Magellanic Cloud (LMC). We have used $\sim$1\hbox{$\,.\!\!^{\prime\prime}$}6\,($\sim$0.4\,pc) resolution measurements of the para-H$_2$CO\,$J_{\rm K_ aK_c}$\,=\,3$_{03}$--2$_{02}$, 3$_{22}$--2$_{21}$, and 3$_{21}$--2$_{20}$ transitions near 218.5\,GHz to constrain RADEX non-LTE models of t… ▽ More We mapped the kinetic temperature structure of two massive star-forming regions, N113 and N159W, in the Large Magellanic Cloud (LMC). We have used $\sim$1\hbox{$\,.\!\!^{\prime\prime}$}6\,($\sim$0.4\,pc) resolution measurements of the para-H$_2$CO\,$J_{\rm K_ aK_c}$\,=\,3$_{03}$--2$_{02}$, 3$_{22}$--2$_{21}$, and 3$_{21}$--2$_{20}$ transitions near 218.5\,GHz to constrain RADEX non-LTE models of the physical conditions. The gas kinetic temperatures derived from the para-H$_2$CO line ratios 3$_{22}$--2$_{21}$/3$_{03}$--2$_{02}$ and 3$_{21}$--2$_{20}$/3$_{03}$--2$_{02}$ range from 28 to 105\,K in N113 and 29 to 68\,K in N159W. Distributions of the dense gas traced by para-H$_2$CO agree with those of the 1.3\,mm dust and \emph{Spitzer}\,8.0\,$μ$m emission, but do not significantly correlate with the H$α$ emission. The high kinetic temperatures ($T_{\rm kin}$\,$\gtrsim$\,50\,K) of the dense gas traced by para-H$_2$CO appear to be correlated with the embedded infrared sources inside the clouds and/or YSOs in the N113 and N159W regions. The lower temperatures ($T_{\rm kin}$\,$<$\,50\,K) are measured at the outskirts of the H$_2$CO-bearing distributions of both N113 and N159W. It seems that the kinetic temperatures of the dense gas traced by para-H$_2$CO are weakly affected by the external sources of the H$α$ emission. The non-thermal velocity dispersions of para-H$_2$CO are well correlated with the gas kinetic temperatures in the N113 region, implying that the higher kinetic temperature traced by para-H$_2$CO is related to turbulence on a $\sim$0.4\,pc scale. The dense gas heating appears to be dominated by internal star formation activity, radiation, and/or turbulence. It seems that the mechanism heating the dense gas of the star-forming regions in the LMC is consistent with that in Galactic massive star-forming regions located in the Galactic plane. △ Less

Submitted 24 August, 2021; originally announced August 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 655, A12 (2021)

The Distribution and Origin of C$_2$H in NGC 253 from ALCHEMI

Authors: J. Holdship, S. Viti, S. Martín, N. Harada, J. Mangum, K. Sakamoto, S. Muller, K. Tanaka, Y. Yoshimura, K. Nakanishi, R. Herrero-Illana, S. Mühle, R. Aladro, L. Colzi, K. L. Emig, S. García-Burillo, C. Henkel, P. Humire, D. S. Meier, V. M. Rivilla, P. van der Werf

Abstract: Observations of chemical species can provide an insight into the physical conditions of the emitting gas but it is important to understand how their abundances and excitation vary within different heating environments. C$_2$H is a molecule typically found in PDR regions of our own Galaxy but there is evidence to suggest it also traces other regions undergoing energetic processing in extragalactic… ▽ More Observations of chemical species can provide an insight into the physical conditions of the emitting gas but it is important to understand how their abundances and excitation vary within different heating environments. C$_2$H is a molecule typically found in PDR regions of our own Galaxy but there is evidence to suggest it also traces other regions undergoing energetic processing in extragalactic environments. As part of the ALCHEMI ALMA large program, the emission of C$_2$H in the central molecular zone of the nearby starburst galaxy NGC 253 was mapped at 1.6 " (28 pc) resolution and characterized to understand its chemical origins. Spectral modelling of the N=1-0 through N=4-3 rotational transitions of C$_2$H was used to derive the C$_2$H column densities towards the dense clouds in NGC 253. Chemical modelling, including PDR, dense cloud, and shock models were then used to investigate the chemical processes and physical conditions that are producing the molecular emission. We find high C$_2$H column densities of $\sim 10^{15} cm^{-3}$ detected towards the dense regions of NGC 253. We further find that these column densities cannot be reproduced by assuming that the emission arises from the PDR regions at the edge of the clouds. Instead, we find that the C$_2$H abundance remains high even in the high visual extinction interior of these clouds and that this is most likely caused by a high cosmic-ray ionization rate. △ Less

Submitted 16 July, 2021; v1 submitted 9 July, 2021; originally announced July 2021.

Comments: 18 pages, 13 figures, accepted for publication in A&A

Journal ref: A&A 654, A55 (2021)

The Galaxy Activity, Torus and Outflow Survey (GATOS): II. Torus and polar dust emission in nearby Seyfert galaxies

Authors: A. Alonso-Herrero, S. García-Burillo, S. F. Hoenig, I. García-Bernete, C. Ramos Almeida, O. González-Martín, E. López-Rodríguez, P. G. Boorman, A. J. Bunker, L. Burtscher, F. Combes, R. Davies, T. Díaz-Santos, P. Gandhi, B. García-Lorenzo, E. K. S. Hicks, L. K. Hunt, K. Ichikawa, M. Imanishi, T. Izumi, A. Labiano, N. A. Levenson, C. Packham, M. Pereira-Santaella, C. Ricci , et al. (8 additional authors not shown)

Abstract: We compare mid-IR and ALMA far-IR images of 12 nearby Seyferts selected from GATOS. The mid-IR unresolved emission contributes more than 60% of the nuclear emission in most galaxies. By contrast, the ALMA 870micron continuum emission is mostly resolved and typically along the torus equatorial direction (Paper I, Garcia-Burillo et al. 2021). The Eddington ratios and nuclear hydrogen column densitie… ▽ More We compare mid-IR and ALMA far-IR images of 12 nearby Seyferts selected from GATOS. The mid-IR unresolved emission contributes more than 60% of the nuclear emission in most galaxies. By contrast, the ALMA 870micron continuum emission is mostly resolved and typically along the torus equatorial direction (Paper I, Garcia-Burillo et al. 2021). The Eddington ratios and nuclear hydrogen column densities NH of half the sample are favorable to launching polar and/or equatorial dusty winds, according to simulations. Six show mid-IR extended emission in the polar direction as traced by the NLR and perpendicular to the ALMA emission. In a few, the nuclear NH might be too high to uplift large quantities of dusty material along the polar direction. Five galaxies have low NH and/or Eddington ratios and thus polar dusty winds are not likely. We generate new CAT3D-WIND disk-wind model images. At low wind-to-disk cloud ratios the far-IR model images have disk- and ring-like morphologies. The X-shape associated with dusty winds is seen better in the far-IR at intermediate-high inclinations for the extended-wind configurations. In most models, the mid-IR emission comes from the inner part of the disk/cone. Extended bi-conical and one-sided polar mid-IR emission is seen in extended-wind configurations and high wind-to-disk cloud ratios. When convolved to our resolution, the model images reproduce qualitative aspects of the observed morphologies. Low-intermediate wind-to-disk ratios are required to account for the large fractions of unresolved mid-IR emission. This work and Paper I provide observational support for the torus+wind scenario. The wind component is more relevant at high Eddington ratios and/or AGN luminosities, and polar dust emission is predicted at NH of up to $10^{24}$cm$^{-2}$. The torus/disk component, on the other hand, prevails at low luminosities and/or Eddington ratios. (Abridged) △ Less

Submitted 2 July, 2021; v1 submitted 1 July, 2021; originally announced July 2021.

Comments: 27 pages, 27 figures. Accepted for publication in A&A

The Galaxy Activity, Torus and Outflow Survey (GATOS) I. ALMA images of dusty molecular tori in Seyfert galaxies

Authors: S. García-Burillo, A. Alonso-Herrero, C. Ramos Almeida, O. González-Martín, F. Combes, A. Usero, S. Hönig, M. Querejeta, E. K. S. Hicks, L. K. Hunt, D. Rosario, R. Davies, P. G. Boorman, A. J. Bunker, L. Burtscher, L. Colina, T. Díaz-Santos, P. Gandhi, I. García-Bernete, B. García-Lorenzo, K. Ichikawa, M. Imanishi, T. Izumi, A. Labiano, N. A. Levenson , et al. (10 additional authors not shown)

Abstract: We present the first results of the Galaxy Activity, Torus and Outflow Survey (GATOS), a project aimed at understanding the properties of the dusty molecular tori and their connection to the host galaxy in nearby Seyfert galaxies. Our project expands the range of AGN luminosities and Eddington ratios covered by previous surveys of Seyferts conducted by ALMA and allows us to study the gas feeding a… ▽ More We present the first results of the Galaxy Activity, Torus and Outflow Survey (GATOS), a project aimed at understanding the properties of the dusty molecular tori and their connection to the host galaxy in nearby Seyfert galaxies. Our project expands the range of AGN luminosities and Eddington ratios covered by previous surveys of Seyferts conducted by ALMA and allows us to study the gas feeding and feedback cycle in a combined sample of 19 Seyferts. We used ALMA to obtain new images of the emission of molecular gas and dust using the CO(3-2) and HCO+(4-3) lines as well as their underlying continuum emission at 870 microns with high spatial resolutions (0.1'' ~ 7 - 13 pc) in the CND of 10 nearby (D < 28 Mpc) Seyfert galaxies. Our new ALMA observations detect 870 micron continuum and CO line emission from spatially resolved disks located around the AGN in all the sources. The bulk of the continuum flux can be accounted for by thermal emission from dust in the majority of the targets. For most of the sources the disks show a preponderant orientation perpendicular to the AGN wind axes, as expected for dusty molecular tori. The median diameters and molecular gas masses of the tori are ~ 42 pc, and ~ 6 x 10**5 Msun, respectively. We find a positive correlation between the line-of-sight gas column densities responsible for the absorption of X-rays and the molecular gas column densities derived from CO towards the AGN in our sources. The radial distributions of molecular gas in the CND of our combined sample show signs of nuclear-scale molecular gas deficits. We also detect molecular outflows in the sources that show the most extreme nuclear-scale gas deficits in our sample. These observations find for the first time supporting evidence that the imprint of AGN feedback is more extreme in higher luminosity and/or higher Eddington ratio Seyfert galaxies. △ Less

Submitted 22 June, 2021; v1 submitted 20 April, 2021; originally announced April 2021.

Comments: 41 pages, new version revised and accepted by A&A on the 22nd of June 2021

Journal ref: A&A 652, A98 (2021)

Are local ULIRGs powered by AGN? The sub-kpc view of the 220 GHz continuum. PUMA II

Authors: M. Pereira-Santaella, L. Colina, S. García-Burillo, I. Lamperti, E. González-Alfonso, M. Perna, S. Arribas, A. Alonso-Herrero, S. Aalto, F. Combes, A. Labiano, J. Piqueras-López, D. Rigopoulou, P. van der Werf

Abstract: We analyze high-resolution (400pc) 220GHz continuum and CO(2-1) ALMA observations of a representative sample of 23 local (z<0.165) ULIRG systems (34 individual nuclei) as part of the "Physics of ULIRGs with MUSE and ALMA" (PUMA) project. The deconvolved half-light radii of the 220GHz continuum sources are between <60-350 pc (median 90pc). We associate these regions with the regions emitting the bu… ▽ More We analyze high-resolution (400pc) 220GHz continuum and CO(2-1) ALMA observations of a representative sample of 23 local (z<0.165) ULIRG systems (34 individual nuclei) as part of the "Physics of ULIRGs with MUSE and ALMA" (PUMA) project. The deconvolved half-light radii of the 220GHz continuum sources are between <60-350 pc (median 90pc). We associate these regions with the regions emitting the bulk of the infrared luminosity. The good agreement, within a factor of 2, between the 220GHz fluxes and the extrapolation of the infrared gray-body, and the small synchrotron and free-free contributions support this assumption. The cold molecular gas emission sizes, r_CO, are 60-700 pc and are similar in advanced mergers and early interacting systems. On average, r_CO are 2.5 times larger than the continuum. We derive L_IR and cold molecular gas surface densities: log Sigma(L_IR)=11.5-14.3 Lsun/kpc^2 and log Sigma(H2)=2.9-4.2 Msun/pc^2. Assuming that the L_IR is produced by star-formation, this corresponds to median Sigma(SFR)=2500 Msun/yr/kpc^2 which would imply extremely short depletion times, <1-15 Myr, and unphysical SF efficiencies >1 for 70% of the sample. Therefore, this favors the presence of obscured AGN that could dominate the L_IR. We also classify the ULIRG nuclei in two groups: (a) compact nuclei (r<130 pc) with high mid-IR excess emission found in optically classified AGN; and (b) nuclei following a relation with decreasing mid-IR excess for decreasing r. 60% of the interacting nuclei lie in the low end (<130 pc) of this relation, while only 30% of the advanced mergers do so, suggesting that in the early interaction phases the activity occurs in more compact and obscured regions. About two thirds of the nuclei are above the Eddington limit which is consistent with the detection of massive outflows in local ULIRGs and the potential role of radiation pressure in the launching process. △ Less

Submitted 10 May, 2021; v1 submitted 16 April, 2021; originally announced April 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 651, A42 (2021)

CON-quest: Searching for the most obscured galaxy nuclei

Authors: N. Falstad, S. Aalto, S. König, K. Onishi, S. Muller, M. Gorski, M. Sato, F. Stanley, F. Combes, E. González-Alfonso, J. G. Mangum, A. S. Evans, L. Barcos-Muñoz, G. C. Privon, S. T. Linden, T. Díaz-Santos, S. Martín, K. Sakamoto, N. Harada, G. A. Fuller, J. S. Gallagher, P. P. van der Werf, S. Viti, T. R. Greve, S. García-Burillo , et al. (6 additional authors not shown)

Abstract: Some luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) host extremely compact and dusty nuclei. The intense infrared radiation arising from warm dust in these sources is prone to excite vibrational levels of molecules such as HCN. This results in emission from the rotational transitions of vibrationally excited HCN (HCN-vib), with the brightest emission found in compact obscured nucl… ▽ More Some luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) host extremely compact and dusty nuclei. The intense infrared radiation arising from warm dust in these sources is prone to excite vibrational levels of molecules such as HCN. This results in emission from the rotational transitions of vibrationally excited HCN (HCN-vib), with the brightest emission found in compact obscured nuclei (CONs). We aim to establish how common CONs are in the local Universe, and whether their prevalence depends on the luminosity or other properties of the host galaxy. We have conducted an Atacama Large Millimeter/submillimeter Array (ALMA) survey of the rotational J=3-2 transition of HCN-vib in a sample of 46 far-infrared luminous galaxies. Compact obscured nuclei are identified in 38 percent of ULIRGs, 21 percent of LIRGs, and 0 percent of lower luminosity galaxies. We find no dependence on the inclination of the host galaxy, but strong evidence of lower IRAS 25 to 60 μm flux density ratios (f25/f60) in CONs compared to the rest of the sample. Furthermore, we find that CONs have stronger silicate features (s9.7μm) but similar PAH equivalent widths (EQW6.2μm) compared to other galaxies. In the local Universe, CONs are primarily found in (U)LIRGs. High resolution continuum observations of the individual nuclei are required to determine if the CON phenomenon is related to the inclinations of the nuclear disks. The lower f25/f60 ratios in CONs as well as the results for the mid-infrared diagnostics investigated are consistent with large dust columns shifting the nuclear radiation to longer wavelengths, making the mid- and far-infrared "photospheres" significantly cooler than the interior regions. To assess the importance of CONs in the context of galaxy evolution, it is necessary to extend this study to higher redshifts where (U)LIRGs are more common. △ Less

Submitted 26 February, 2021; originally announced February 2021.

Comments: 16 pages, 11 figures, Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 649, A105 (2021)

Gas phase Elemental abundances in Molecular cloudS (GEMS). IV. Observational results and statistical trends

Authors: M. Rodríguez-Baras, A. Fuente, P. Riviére-Marichalar, D. Navarro-Almaida, P. Caselli, M. Gerin, C. Kramer, E. Roueff, V. Wakelam, G. Esplugues, S. García-Burillo, R. Le Gal, S. Spezzano, T. Alonso-Albi, R. Bachiller, S. Cazaux, B. Commercon, J. R. Goicoechea, J. C. Loison, S. P. Treviño-Morales, O. Roncero, I. Jiménez-Serra, J. Laas, A. Hacar, J. Kirk , et al. (11 additional authors not shown)

Abstract: Gas phase Elemental abundances in Molecular CloudS (GEMS) is an IRAM 30m Large Program designed to estimate the S, C, N, and O depletions and gas ionization degree, X(e-), in a set of star-forming filaments of Taurus, Perseus and Orion. Our immediate goal is to build up a complete database of molecular abundances that can serve as an observational basis for estimating X(e-) and the C, O, N, and S… ▽ More Gas phase Elemental abundances in Molecular CloudS (GEMS) is an IRAM 30m Large Program designed to estimate the S, C, N, and O depletions and gas ionization degree, X(e-), in a set of star-forming filaments of Taurus, Perseus and Orion. Our immediate goal is to build up a complete database of molecular abundances that can serve as an observational basis for estimating X(e-) and the C, O, N, and S depletions through chemical modeling. We observed and derived the abundances of 14 species (13CO, C18O, HCO+, H13CO+, HC18O+, HCN, H13CN, HNC, HCS+, CS, SO, 34SO, H2S, and OCS) in 244 positions, covering the AV 3 to 100 mag, n(H2) a few 10$^{3}$ to 10$^6$ cm$^{-3}$, and Tk 10 to 30 K ranges in these clouds, avoiding protostars, HII regions, and outflows. A statistical analysis is carried out to identify general trends between different species and with physical parameters. Relations between molecules reveal strong linear correlations which define three different families: (1) 13CO and C18O; (2) H13CO+, HC18O+, H13CN, and HNC; and (3) the S-bearing molecules. The abundances of the CO isotopologs increase with the gas kinetic temperature until TK 15 K. For higher temperatures, the abundance remains constant with a scatter of a factor of 3. The abundances of H13CO+, HC18O+, H13CN, and HNC are well correlated with each other, and all of them decrease with molecular hydrogen density, following the law n(H2)$^{-0.8\pm0.2}$. The abundances of S-bearing species also decrease with n(H2) at a rate of (S-bearing/H)gas n(H2)$^{-0.6\pm0.1}$. The abundances of molecules belonging to groups 2 and 3 do not present any clear trend with gas temperature. At scales of molecular clouds, the C18O abundance is the quantity that better correlates with the cloud mass. We discuss the utility of the 13CO/C18O, HCO+/H13CO+, and H13CO+/H13CN abundance ratios as chemical diagnostics of star formation in external galaxies. △ Less

Submitted 25 February, 2021; originally announced February 2021.

Comments: 28 pages, 22 figures. Accepted for publication in A&A

Journal ref: A&A 648, A120 (2021)

The Diverse Molecular Gas Content of Massive Galaxies Undergoing Quenching at z~1

Authors: Sirio Belli, Alessandra Contursi, Reinhard Genzel, Linda J. Tacconi, Natascha M. Förster-Schreiber, Dieter Lutz, Françoise Combes, Roberto Neri, Santiago García-Burillo, Karl F. Schuster, Rodrigo Herrera-Camus, Ken-ichi Tadaki, Rebecca L. Davies, Richard I. Davies, Benjamin D. Johnson, Minju M. Lee, Joel Leja, Erica J. Nelson, Sedona H. Price, Jinyi Shangguan, T. Taro Shimizu, Sandro Tacchella, Hannah Übler

Abstract: We present a detailed study of the molecular gas content and stellar population properties of three massive galaxies at 1 < z < 1.3 that are in different stages of quenching. The galaxies were selected to have a quiescent optical/near-infrared spectral energy distribution and a relatively bright emission at 24 micron, and show remarkably diverse properties. CO emission from each of the three galax… ▽ More We present a detailed study of the molecular gas content and stellar population properties of three massive galaxies at 1 < z < 1.3 that are in different stages of quenching. The galaxies were selected to have a quiescent optical/near-infrared spectral energy distribution and a relatively bright emission at 24 micron, and show remarkably diverse properties. CO emission from each of the three galaxies is detected in deep NOEMA observations, allowing us to derive molecular gas fractions Mgas/Mstar of 13-23%. We also reconstruct the star formation histories by fitting models to the observed photometry and optical spectroscopy, finding evidence for recent rejuvenation in one object, slow quenching in another, and rapid quenching in the third system. To better constrain the quenching mechanism we explore the depletion times for our sample and other similar samples at z~0.7 from the literature. We find that the depletion times are highly dependent on the method adopted to measure the star formation rate: using the UV+IR luminosity we obtain depletion times about 6 times shorter than those derived using dust-corrected [OII] emission. When adopting the star formation rates from spectral fitting, which are arguably more robust, we find that recently quenched galaxies and star-forming galaxies have similar depletion times, while older quiescent systems have longer depletion times. These results offer new, important constraints for physical models of galaxy quenching. △ Less

Submitted 15 February, 2021; originally announced February 2021.

Comments: 9 pages, 4 figures. Accepted for publication in ApJ Letters

Gas-phase Elemental abundances in Molecular cloudS (GEMS) III. Unlocking the CS chemistry: the CS+O reaction

Authors: N. Bulut, O. Roncero, A. Aguado, J. -C. Loison, D. Navarro-Almaida, V. Wakelam, A. Fuente, E. Roueff, R. Le Gal, P. Caselli, M. Gerinm K. M. Hickson, S. Spezzano, P. Riviere-Marichalar, T. Alonso-Albi, R. Bachiller, I. Jimenez-Serra, C. Kramer, B. Tercero, M. Rodriguez-Baras, S. Garcia-Burillo, J. R. Goicoechea, S. P. Treviño-Morales, G. Esplugues, S. Cazaux, B. Commercon , et al. (14 additional authors not shown)

Abstract: CS is among the most abundant gas-phase S-bearing molecules in cold dark molecular clouds. It is easily observable with several transitions in the millimeter wavelength range, and has been widely used as a tracer of the gas density in the interstellar medium in our Galaxy and external galaxies. Chemical models fail to account for the observed CS abundances when assuming the cosmic value for the el… ▽ More CS is among the most abundant gas-phase S-bearing molecules in cold dark molecular clouds. It is easily observable with several transitions in the millimeter wavelength range, and has been widely used as a tracer of the gas density in the interstellar medium in our Galaxy and external galaxies. Chemical models fail to account for the observed CS abundances when assuming the cosmic value for the elemental abundance of sulfur. The CS+O -> CO + S reaction has been proposed as a relevant CS destruction mechanism at low temperatures, and could explain the discrepancy between models and observations. Its reaction rate has been experimentally measured at temperatures of 150-400 K, but the extrapolation to lower temperatures is doubtful. Here we calculate the CS+O reaction rate at temperatures <150 K which are prevailing in the interstellar medium. We performed ab initio calculations to obtain the three lowest PES of the CS+O system. These PESs are used to study the reaction dynamics, using several methods to eventually calculate the CS+O thermal reaction rates. We compare the results of our theoretical calculations for 150-400 K with those obtained in the laboratory. Our detailed theoretical study on the CS+O reaction, which is in agreement with the experimental data obtained at 150-400 K, demonstrates the reliability of our approach. After a careful analysis at lower temperatures, we find that the rate constant at 10 K is negligible, which is consistent with the extrapolation of experimental data using the Arrhenius expression. We use the updated chemical network to model the sulfur chemistry in TMC1 based on molecular abundances determined from GEMS project observations. In our model, we take into account the expected decrease of the cosmic ray ionization rate along the cloud. The abundance of CS is still overestimated when assuming the cosmic value for the sulfur abundance. △ Less

Submitted 18 December, 2020; originally announced December 2020.

Comments: 11 pages, 10 figures, accepted in Astronomy & Astrophysics

Journal ref: Astronomy & Astrophysics 646, A5 (2021)

Black hole feeding and star formation in NGC 1808

Authors: A. Audibert, F. Combes, S. García-Burillo, L. Hunt, A. Eckart, S. Aalto, V. Casasola, F. Boone, M. Krips, S. Viti, S. Muller, K. Dasyra, P. van der Werf, S. Martín

Abstract: We report on Atacama Large Millimeter Array (ALMA) observations of CO(3-2) emission in the Seyfert2/starburst galaxy NGC1808, at a spatial resolution of 4pc. Our aim is to investigate the morphology and dynamics of the gas inside the central 0.5kpc and to probe the nuclear feeding and feedback phenomena. We discovered a nuclear spiral of radius 1"=45pc. Within it, we found a decoupled circumnuclea… ▽ More We report on Atacama Large Millimeter Array (ALMA) observations of CO(3-2) emission in the Seyfert2/starburst galaxy NGC1808, at a spatial resolution of 4pc. Our aim is to investigate the morphology and dynamics of the gas inside the central 0.5kpc and to probe the nuclear feeding and feedback phenomena. We discovered a nuclear spiral of radius 1"=45pc. Within it, we found a decoupled circumnuclear disk or molecular torus of a radius of 0.13"=6pc. The HCN(4-3) and HCO$\rm^+$(4-3) and CS(7-6) dense gas line tracers were simultaneously mapped and detected in the nuclear spiral and they present the same misalignment in the molecular torus. At the nucleus, the HCN/HCO$^+$ and HCN/CS ratios indicate the presence of an active galactic nucleus (AGN). The molecular gas shows regular rotation, within a radius of 400pc, except for the misaligned disk inside the nuclear spiral arms. The computations of the torques exerted on the gas by the barred stellar potential reveal that the gas within a radius of 100pc is feeding the nucleus on a timescale of five rotations or on an average timescale of ~60Myr. Some non-circular motions are observed towards the center, corresponding to the nuclear spiral arms. We cannot rule out that small extra kinematic perturbations could be interpreted as a weak outflow attributed to AGN feedback. The molecular outflow detected at $\geqslant$250pc in the NE direction is likely due to supernovae feedback and it is connected to the kpc-scale superwind. △ Less

Submitted 17 September, 2021; v1 submitted 17 November, 2020; originally announced November 2020.

Comments: 15 pages, 18 figures, accepted to A&A. arXiv admin note: text overlap with arXiv:1905.01979

Evidence of galaxy interaction in the Narrow-line Seyfert 1 galaxy IRAS17020+4544 seen by NOEMA

Authors: Q. Salomé, A. L. Longinotti, Y. Krongold, C. Feruglio, V. Chavushyan, O. Vega, S. García-Burillo, A. Fuente, A. Olguín-Iglesias, V. M. Patiño-Álvarez, I. Puerari, A. Robleto-Orús

Abstract: The narrow-line Seyfert 1 galaxy IRAS17020+4544 is one of the few sources where both an X-ray ultra-fast outflow and a molecular outflow were observed to be consistent with energy conservation. However, IRAS17020+4544 is less massive and has a much more modest active galactic nucleus (AGN) luminosity than the other examples. Using recent CO(1-0) observations with the NOrthern Extended Millimeter A… ▽ More The narrow-line Seyfert 1 galaxy IRAS17020+4544 is one of the few sources where both an X-ray ultra-fast outflow and a molecular outflow were observed to be consistent with energy conservation. However, IRAS17020+4544 is less massive and has a much more modest active galactic nucleus (AGN) luminosity than the other examples. Using recent CO(1-0) observations with the NOrthern Extended Millimeter Array (NOEMA), we characterised the molecular gas content of the host galaxy for the first time. We found that the molecular gas is distributed into an apparent central disc of 1.1x10^9 Msun, and a northern extension located up to 8 kpc from the centre with a molecular gas mass M_H2~10^8 Msun. The molecular gas mass and the CO dynamics in the northern extension reveal that IRAS 17020+4544 is not a standard spiral galaxy, instead it is interacting with a dwarf object corresponding to the northern extension. This interaction possibly triggers the high accretion rate onto the super massive black hole. Within the main galaxy, which hosts the AGN, a simple analytical model predicts that the molecular gas may lie in a ring, with less molecular gas in the nuclear region. Such distribution may be the result of the AGN activity which removes or photodissociates the molecular gas in the nuclear region (AGN feedback). Finally, we have detected a molecular outflow of mass M_H2=(0.7-1.2)x10^7 Msun in projection at the location of the northern galaxy, with a similar velocity to that of the massive outflow reported in previous millimeter data obtained by the Large Millimeter Telescope. △ Less

Submitted 6 January, 2021; v1 submitted 17 November, 2020; originally announced November 2020.

Comments: Published in MNRAS, Volume 501, Issue 1, Pages 219-228

ALMA Resolves Giant Molecular Clouds in a Tidal Dwarf Galaxy

Authors: M. Querejeta, F. Lelli, E. Schinnerer, D. Colombo, U. Lisenfeld, C. G. Mundell, F. Bigiel, S. García-Burillo, C. N. Herrera, A. Hughes, J. M. D. Kruijssen, S. E. Meidt, T. J. T. Moore, J. Pety, A. J. Rigby

Abstract: Tidal dwarf galaxies (TDGs) are gravitationally bound condensations of gas and stars formed during galaxy interactions. Here we present multi-configuration ALMA observations of J1023+1952, a TDG in the interacting system Arp 94, where we resolve CO(2-1) emission down to giant molecular clouds (GMCs) at 0.64" ~ 45pc resolution. We find a remarkably high fraction of extended molecular emission (~80-… ▽ More Tidal dwarf galaxies (TDGs) are gravitationally bound condensations of gas and stars formed during galaxy interactions. Here we present multi-configuration ALMA observations of J1023+1952, a TDG in the interacting system Arp 94, where we resolve CO(2-1) emission down to giant molecular clouds (GMCs) at 0.64" ~ 45pc resolution. We find a remarkably high fraction of extended molecular emission (~80-90%), which is filtered out by the interferometer and likely traces diffuse gas. We detect 111 GMCs that give a similar mass spectrum as those in the Milky Way and other nearby galaxies (a truncated power law with slope of -1.76+/-0.13). We also study Larson's laws over the available dynamic range of GMC properties (~2 dex in mass and ~1 dex in size): GMCs follow the size-mass relation of the Milky Way, but their velocity dispersion is higher such that the size-linewidth and virial relations appear super-linear, deviating from the canonical values. The global molecular-to-atomic gas ratio is very high (~1) while the CO(2-1)/CO(1-0) ratio is quite low (~0.5), and both quantities vary from north to south. Star formation is predominantly taking place in the south of the TDG, where we observe projected offsets between GMCs and young stellar clusters ranging from ~50pc to ~200pc; the largest offsets correspond to the oldest knots, as seen in other galaxies. In the quiescent north, we find more molecular clouds and a higher molecular-to-atomic gas ratio (~1.5); atomic and diffuse molecular gas also have a higher velocity dispersion there. Overall, the organisation of the molecular ISM in this TDG is quite different from other types of galaxies on large scales, but the properties of GMCs seem fairly similar, pointing to near universality of the star-formation process on small scales. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Comments: 21 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 645, A97 (2021)