Search

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

We present observations of the type-2 Seyfert NGC7172 obtained with the medium-resolution spectrometer (MRS) of the Mid-Infrared Instrument (MIRI) on board of the JWST. This galaxy hosts one of the lowest ionised gas mass outflow rates (Mout~0.005 M/yr) in a sample of six AGN with similar bolometric luminosities (log Lbol~44erg/s) within the Galactic Activity, Torus and Outflow Survey (GATOS). We aim to understand the properties of the ionised gas outflow, mainly using the emission lines from the neon transitions, that cover a broad range of ionisation potentials (IP) from ~20 eV to ~130 eV. We applied parametric and non-parametric methods to characterise the line emission and kinematics. The low excitation lines (IP<25eV, e.g.[NeII]) trace the rotating disc emission. The high excitation lines (IP>90eV, e.g.[NeV]), which are likely photoionised exclusively by the AGN, are expanding in the direction nearly perpendicular to the disc of the galaxy, with maximum projected velocities of ~350-500 km/s. In particular, [NeV] and [NeVI] lines reveal a biconical ionised gas outflow emerging N-S from the nuclear region, extending at least ~2.5"N and 3.8"S (projected distance of ~450 and 680 pc). Most of the emission arising in the northern part of the cone was not previously detected due to obscuration. Given the almost face-on orientation of the outflow and the almost edge-on orientation of the galaxy, NGC7172 may be a case of weak coupling. Nevertheless, we found evidence for positive feedback in two distinct outflowing clumps at projected distances of 3.1" and 4.3" (i.e. ~560 and 780 pc) SW from the AGN. We estimated a star formation rate in these regions using the [NeII] and [NeIII] luminosities of 0.08 M/yr, that is ~10% of that found in the circumnuclear ring. The star formation activity might have been triggered by the interaction between the ionised gas outflow and the ISM of the galaxy., Comment: 20 pages, 6 figures. Accepted for publication in A&A. Abstract adapted for the arxiv version

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

Compact obscured nuclei (CONs) are an extremely obscured (N$_{H2}$ >10$^{25}$ cm$^{-2}$) class of galaxy nuclei thought to exist in 20-40 per cent of nearby (ultra-)luminous infrared galaxies. While they have been proposed to represent a key phase of the active galactic nucleus (AGN) feedback cycle, the nature of these CONs - what powers them, their dynamics, and their impact on the host galaxy - remains unknown. This work analyses the large-scale optical properties of the local CON, NGC4418 (z=0.00727). We present new, targeted integral field unit observations of the galaxy with the Multi-Unit Spectroscopic Explorer (MUSE). For the first time, we map the ionised and neutral gas components of the galaxy, along with their dynamical structure, to reveal several previously unknown features of the galaxy. We confirm the presence of a previously postulated blueshifted outflow along the minor axis of NGC4418. We find this outflow to be decelerating and, for the first time, show it to extend bilaterally from the nucleus. We report the discovery of two further outflow structures: a redshifted southern outflow connected to a tail of ionised gas surrounding the galaxy and a blueshifted bubble to the north. In addition to these features, we find the [OIII] emission reveals the presence of knots across the galaxy, which are consistent with regions of the galaxy that have been photoionised by an AGN. Based on the properties of these features, we conclude that the CON in NGC4418 is most likely powered by AGN activity., Comment: 17 pages, 15 figures, accepted to A&A

We present an analysis of the combined NIRSpec and MIRI spectra of dusty galaxies between 1.5 - 28 $\mu$m restframe by implementing a differential extinction model, where the strength of extinction varies across the spectrum as different layers of the obscuring dust are probed. Our model is able to recover a 2D distribution of dust temperature and extinction which allows inference of the physical nature of the dust in these environments. We show that differential extinction is necessary to reproduce the spectra of 4 highly obscured Luminous Infrared Galaxies observed with NIRSpec IFU and MIRI MRS, where simple screen or uniformly mixed dust distributions fail to fit the data. We additionally compare the extinction of HII regions in these galaxies via hydrogen recombination lines, the extinction of molecular gas via the H$_2$ lines, Polycyclic Aromatic Hydrocarbons via the 12.7/11.3 PAH ratio and the stellar continuum. We find that the molecular gas is deeply buried with the HII regions in star-forming regions, with a similar extinction to the hottest dust components. However we find the cooler dust to be less obscured, at a similar extinction to the stellar continuum and PAHs. The nuclei show a complex dust distribution with VV114 NE, NGC 3256 S, IIZw96 SW showing a deeply buried continuum source relative to the molecular gas/HII regions. Additionally, NGC 3256 S, NGC 7469 and VV114 SW show an isolated hot dust component, indicative of AGN heating, where NGC 3256 S and NGC 7469 are previously known AGN., Comment: 19 pages, 13 figures. Accepted for publication in MNRAS

We analyze JWST NIRSpec$+$MIRI/MRS observations of the infrared (IR) gas-phase molecular bands of the most enshrouded source (D1) within the interacting system and luminous IR galaxy II Zw 096. We report the detection of rovibrational lines of H$_2$O $\nu_2$=1-0 ($\sim$5.3-7.2 $\mu$m) and $^{12}$CO $\nu$=1-0 ($\sim$4.45-4.95 $\mu$m) in D1. The CO band shows the R- and P-branches in emission and the spectrum of the H$_2$O band shows the P-branch in emission and the R-branch in absorption. The H$_2$O R-branch in absorption unveils an IR-bright embedded compact source in D1 and the CO broad component features a highly turbulent environment. From both bands, we also identified extended intense star-forming (SF) activity associated with circumnuclear photodissociation regions (PDRs), consistent with the strong emission of the ionised 7.7 $\mu$m polycyclic aromatic hydrocarbon band in this source. By including the 4.5-7.0 $\mu$m continuum information derived from the H$_2$O and CO analysis, we modelled the IR emission of D1 with a dusty torus and SF component. The torus is very compact (diameter of $\sim$3 pc at 5 $\mu$m) and characterised by warm dust ($\sim$ 370 K), giving an IR surface brightness of $\sim$3.6$\times$10$^{8}$ L$_{\rm sun}$/pc$^2$. This result suggests the presence of a dust-obscured active galactic nucleus (AGN) in D1, which has an exceptionally high covering factor that prevents the direct detection of AGN emission. Our results open a new way to investigate the physical conditions of inner dusty tori via modelling the observed IR molecular bands., Comment: Accepted for publication in A&A (official date of acceptance: 27/12/2023; submitted 27/11/2023). 11 pages, 10 figures

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

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

We present a study of the mid-IR (MIR) emission of quiescent galaxies (QGs) beyond the local universe. Using deep $JWST$ imaging in the SMACS-0723 cluster field we identify a mass limited ($M_{*} >10^{9}$M$_{\odot}$) sample of intermediate redshift QGs ($0.2 5 \mu$m, can be attributed to different dust continuum levels and/or the presence of Polycyclic Aromatic Hydrocarbons (PAHs) features. The latter would indicate enhanced 11.3- and 12.7 $\mu$m PAHs strengths with respect to those at 6.2- and 7.7$ \mu$m, consistent with the observed spectra of local ellipticals and indicative of soft radiation fields. Finally, we augment the average UV-to-MIR spectrum of the population with cold dust and gas emission in the far-IR/mm and construct a panchromatic UV-to-radio SED that can serve as a template for the future exploration of the interstellar medium of $z>0$ QGs with ALMA and $JWST$., Comment: The panchromatic QG SED has been made publicly available at http://www.georgiosmagdis.com/software

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

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

Four ultra-luminous infrared galaxies (ULIRGs) observed with JWST/NIRcam in the Cosmos Evolution Early Release Science program offer an unbiased preview of the $z\approx2$ ULIRG population. The objects were originally selected at 24 $\mu$m and have strong polycyclic aromatic hydrocarbon emission features observed with Spitzer/IRS. The four objects have similar stellar masses of ${\sim}10^{11}$ M$_\odot$ but otherwise are quite diverse. One is an isolated disk galaxy, but it has an active nucleus as shown by X-ray observations and by a bright point-source nucleus. Two others are merging pairs with mass ratios of 6-7:1. One has active nuclei in both components, while the other has only one active nucleus: the one in the less-massive neighbor, not the ULIRG. The fourth object is clumpy and irregular and is probably a merger, but there is no sign of an active nucleus. The intrinsic spectral energy distributions for the four AGNs in these systems are typical of type-2 QSOs. This study is consistent with the idea that even if internal processes can produce large luminosities at $z\sim2$, galaxy merging may still be necessary for the most luminous objects. The diversity of these four initial examples suggests that large samples will be needed to understand the $z\approx2$ ULIRG population., Comment: 9 pages, 3 figures, 1 table, accepted for publication in ApJ. V2 updates author affiliations and acknowledgments, not scientific content

Mrk 477 is the closest type-2 quasar (QSO2), at a distance of 163 Mpc. This makes it an ideal laboratory for studying the interplay between nuclear activity and star formation with a great level of detail and signal-to-noise. In this Letter we present new mid-infrared (mid-IR) imaging and spectroscopic data with an angular resolution of 0.4 arcsec (~300 pc) obtained with the Gran Telescopio Canarias (GTC) instrument CanariCam. The N-band (8-13 micron) spectrum of the central ~400 pc of the galaxy reveals [S IV]10.51 micron emission, but no 8.6 or 11.3 micron PAH features, which are commonly used as tracers of recent star formation. This is in stark contrast with the presence of a nuclear starburst of ~300 pc in size, an age of 6 Myr, and a mass of 1.1x10^8 Msun, as constrained from ultraviolet Hubble Space Telescope observations. Considering this, we argue that even the more resilient, neutral molecules that mainly produce the 11.3 micron PAH band are most likely being destroyed in the vicinity of the active nucleus despite the relatively large X-ray column density, of log N_H=23.5 cm^-2, and modest X-ray luminosity, of 1.5x10^43 erg/s. This highlights the importance of being cautious when using PAH features as star formation tracers in the central region of galaxies to evaluate the impact of feedback from active galactic nuclei., Comment: 7 pages, 3 figures. A&A Letters accepted version including proof corrections

Context. Compact Obscured Nuclei (CONs) are an extreme phase of galaxy evolution where rapid supermassive black hole growth and$/$or compact star-forming activity is completely obscured by gas and dust. Aims. We investigate the properties of CONs in the mid-infrared and explore techniques aimed at identifying these objects such as through the equivalent width (EW) ratios of their Polycyclic Aromatic Hydrocarbon (PAH) features. Methods. We model Spitzer spectra by decomposing the continua into nuclear and star-forming components from which we then measure the nuclear optical depth, $\tau_N$, of the $9.8 \mu$m silicate absorption feature. We also use Spitzer spectral maps to investigate how PAH EW ratios vary with aperture size for objects hosting CONs. Results. We find that the nuclear optical depth, $\tau_N$, strongly correlates with the HCN-vib emission line in the millimetre for CONs with a Pearson correlation coefficient of 0.91. We find the PAH EW ratios technique to be effective at selecting CONs and robust against highly inclined galaxies where strong dust lanes may mimic a CON like spectrum by producing a high $\tau_N$. Our analysis of the Spitzer spectral maps showed that the efficacy of the PAH EW ratios to isolate CONs is reduced when there is a strong star-forming component from the host galaxy. In addition, we find that the use of the inferred nuclear optical depth is a reliable method to identify CONs in $36^{+8}_{-7}\%$ of ULIRGs and $17^{+3}_{-3}\%$ of LIRGs, consistent with previous work. Conclusions. We confirm mid-IR spectra to be a powerful diagnostic of CONs where the increased sensitivity of JWST will allow identification of CONs at cosmic noon revealing this extreme but hidden phase of galaxy evolution., Comment: Accepted for publication in A&A. 23 pages, 15 figures, 3 tables

Compact Obscured Nuclei (CONs) potentially hide extreme supermassive black hole (SMBH) growth behind large column densities of gas/dust. We present a spectroscopic analysis of the heavily obscured nucleus and the surrounding environment of the eastern region of the nearby ($z = 0.02007$) interacting galaxy VV 114 with the JWST Mid-InfraRed Instrument (MIRI). We model the spectrum from 4.9 - 28 $\mu$m to extract Polycyclic Aromatic Hydrocarbon (PAH) emission and the underlying obscured continuum. We find that the NE nucleus (A) is highly obscured where the low PAH equivalent width (EW) ratio, EW(12.7)/EW(11.3), reveals a dust enshrouded continuum source. This is confirmed by decomposing the continuum into nuclear and star-forming where the nuclear component is found to be typical of CONs. The 11.3/6.2 PAH flux ratio is consistent with originating in star-forming regions rather than typical AGN. The second nucleus (B) is much less obscured, with PAH flux ratios also typical of star-forming regions. We do not detect any high ionisation lines such as [Ne V] or [Ne VI] which suggests that if an AGN is present it must be highly obscured. Additionally, we detect a shock front south of the secondary nucleus (B) in the [Fe II] (5.34 $\mu$m) line and in warm molecular hydrogen. The 6.2 PAH emission does not spatially coincide with the low-J transitions of H$_2$ but rather appears strong at the shock front which may suggest destruction of the ionised PAHs in the post-shock gas behind the shock front., Comment: 15 pages, 15 figures. Accepted for publication in MNRAS

We present a statistical analysis of the properties of the obscuring material around active galactic nuclei (AGN). This study represents the first of its kind for an ultra-hard X-ray (14-195keV; Swift/BAT) volume-limited (DL<40 Mpc) sample of 24 Seyfert (Sy) galaxies (BCS40 sample) using high angular resolution infrared data and various torus models: smooth, clumpy and two-phase torus models and clumpy disc+wind models. We find that the smooth, clumpy and two-phase torus models (i.e. without including the polar dusty wind component) and disc+wind models provide best fits for a comparable number of galaxies, 8/24 (33.3%) and 9/24 (37.5%), respectively. We find that the best-fitted models depend on the hydrogen column density (NH), which is related to the X-ray (unobscured/obscured) and/or optical (Sy1/Sy2) classification. In particular, smooth, clumpy and two-phase torus models best reproduce the infrared (IR) emission of AGN with relatively high hydrogen column density (median value of log (NH)=23.5+-0.8; i.e. Sy2). However, clumpy disc+wind models provide the best fits to the nuclear IR spectral energy distributions (SEDs) of Sy1/1.8/1.9 (median value of log (NH)=21.0+-1.0), specifically in the near-IR (NIR) range. The success of the disc+wind models in fitting the NIR emission of Sy1 galaxies is due to the combination of adding large graphite grains to the dust composition and self-obscuration effects caused by the wind at intermediate inclinations. In general, we find that the Seyfert galaxies having unfavourable (favourable) conditions, i.e. nuclear hydrogen column density and Eddington ratio, for launching IR dusty polar outflows are best-fitted with smooth, clumpy and two-phase torus (disk+wind) models confirming the predictions from simulations. Therefore, our results indicate that the nature of the inner dusty structure in AGN depend on the intrinsic AGN properties., Comment: Accepted for publication in A&A, 43 pages, 43 figures, 16 tables

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

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., Comment: 12 pages, 10 figures. Small changes to match the accepted for publication version in A&A (official date of acceptance: 20/09/2022)

We present the ultraviolet to submillimetre spectral energy distributions (SEDs) of the HERschel Ultra Luminous Infrared Galaxy Survey (HERUS) sample of 42 local ultraluminous infrared galaxies (ULIRGs) and fit them with a Markov chain Monte Carlo (MCMC) code using the CYprus models for Galaxies and their NUclear Spectra (CYGNUS) radiative transfer models for starbursts, active galactic nucleus (AGN) tori and host galaxy. The Spitzer IRS spectroscopy data are included in the fitting. Our bayesian SED fitting method takes comparable time to popular energy balance methods but it is more physically motivated and versatile. All HERUS galaxies harbor high rates of star formation but we also find bolometrically significant AGN in all of the galaxies of the sample. We estimate the correction of the luminosities of the AGN in the ULIRGs due to the anisotropic emission of the torus and find that it could be up to a factor of $\sim10$ for nearly edge-on tori. We present a comparison of our results with the smooth torus model of Fritz et al. and the two-phase models of Siebenmorgen et al. and SKIRTOR. We find that the CYGNUS AGN torus models fit significantly better the SEDs of our sample compared to all other models. We find no evidence that strong AGN appear either at the beginning or end of a starburst episode or that starbursts and AGN affect each other. IRAS 01003-2238 and Mrk 1014 show evidence for dual AGN in their SED fits suggesting a minimum dual AGN fraction in the sample of 5%., Comment: Published in Monthly Notices of the Royal Astronomical Society. The CYGNUS models used in this paper are available at https://arc.euc.ac.cy/cygnus/

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., Comment: 36 pages, 20 figures, 8 tables. Accepted for publication in A&A

Compact obscured nuclei (CONs) are mainly found in local U/LIRGs. In the local Universe, these sources are generally selected through the detection of the HCN-vib (3-2) emission line at submillimetre wavelengths. In this work, we present a diagnostic method to select deeply buried nuclei based on mid-infrared (mid-IR) polycyclic aromatic hydrocarbons (PAHs) and continuum ratios. Using Spitzer/IRS spectra of a representative sample of local ULIRGs (z<0.27), we examine their PAH and underlying continuum emission ratios. For deeply embedded sources, we find that the 9.7 micron silicate absorption band has a particularly pronounced effect on the 11.3 micron PAH feature. The low flux level in the nuclear silicate absorption band enhances the 11.3 micron PAH feature contrast (high PAH equivalent width) compared to that of the other PAH features. The technique has been extended to include the use of the continuum ratios. However, the latter are affected both by the extinction coming from the host galaxy as well as the nuclear region, whereas the foreground extinction is cancelled out when using the PAH equivalent width ratios. We apply our method to the HERUS and GOALS samples and classify as CON candidates 14 ULIRGs and 10 LIRGs, corresponding to 30% of ULIRGs and 7% of LIRGs from these samples. We find that the observed continuum ratios of CON-dominated sources can be explained by assuming torus models with a tapered disk geometry and a smooth dust distribution. This suggests that the nuclear dusty structure of CONs has an extremely high dust coverage. We also demonstrate that the use of mid-IR color-color diagrams is an effective way to select CON-dominated sources at different redshifts. In particular, the combination of filters of the JWST/MIRI will enable the selection of CONs out to z~1.5. This will allow extending the selection of CONs to high redshifts where U/LIRGs are more numerous., Comment: 13 pages, 11 figures, 1 table, accepted for publication in A&A

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$\alpha$, 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., Comment: 19 pages, 17 figures, accepted for publication in MNRAS

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), Comment: 27 pages, 27 figures. Accepted for publication in A&A

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., Comment: 41 pages, new version revised and accepted by A&A on the 22nd of June 2021

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., Comment: Accepted for publication in A&A

We describe a complete, flux-density-limited sample of galaxies at redshift $0.8 < z < 1.3$ selected at 16 micron. At the selection wavelength near 8 micron rest, the observed emission comes both from dust heated by intense star formation and from active galactic nuclei (AGNs). Fitting the spectral energy distributions (SEDs) of the sample galaxies to local-galaxy templates reveals that more than half the galaxies have SEDs dominated by star formation. About one sixth of the galaxy SEDs are dominated by an AGN, and nearly all the rest of the SEDs are composite. Comparison with X-ray and far-infrared observations shows that combinations of luminosities at rest-frame 4.5 and 8 micron give good measures of both AGN luminosity and star-formation rate. The sample galaxies mostly follow the established star-forming main sequence for $z=1$ galaxies, but of the galaxies more than 0.5 dex above that main sequence, more than half have AGN-type SEDs. Similarly, the most luminous AGNs tend to have higher star-formation rates than the main sequence value. Galaxies with stellar masses $>$10$^{11}$\,\Msun\ are unlikely to host an AGN. About 1% of the sample galaxies show an SED with dust emission typical of neither star formation nor an AGN., Comment: 14 pages, 5 tables, 29 figures, ApJ accepted v2 corrects author name formatting

The Origins Space Telescope (Origins) is one of four science and technology definition studies selected by National Aeronautics and Space Administration (NASA) in preparation of the 2020 Astronomy and Astrophysics Decadal survey in the US. Origins will trace the history of our origins from the time dust and heavy elements permanently altered the cosmic landscape to present-day life. It is designed to answer three major science questions: How do galaxies form stars, make metals, and grow their central supermassive black holes from reionization? How do the conditions for habitability develop during the process of planet formation? Do planets orbiting M-dwarf stars support life? Origins operates at mid- to far-infrared wavelengths from ~2.8 to 588 {\mu}m, is more than 1000 times more sensitive than prior far-IR missions due to its cold (~4.5 K) aperture and state-of-the-art instruments., Comment: 24 pages, 12 figures, 5 tables. Submitted as White paper to Voyage 2050 Long-term planning of the ESA Science Programme. arXiv admin note: substantial text overlap with arXiv:1912.06213

Polycyclic Aromatic Hydrocarbons (PAHs) are carbon-based molecules resulting from the union of aromatic rings and related species, which are likely responsible for strong infrared emission features (3.3, 6.2, 7.7, 8.6, 11.3 and 12.7 microns). In this work, using a sample of 50 Seyfert galaxies (DL<100 Mpc) we compare the circumnuclear (inner kpc) PAH emission of AGN to that of a control sample of star-forming galaxies (22 luminous infrared galaxies and 30 HII galaxies), and investigate the differences between central and extended PAH emission. Using Spitzer/InfraRed Spectrograph spectral data of Seyfert and star-forming galaxies and newly developed PAH diagnostic model grids, derived from theoretical spectra, we compare the predicted and observed PAH ratios. We find that star-forming galaxies and AGN-dominated systems are located in different regions of the PAH diagnostic diagrams. This suggests that not only are the size and charge of the PAH molecules different, but also the nature and hardness of the radiation field that excite them. We find tentative evidence that PAH ratios in AGN-dominated systems are consistent with emission from larger PAH molecules (Nc>300-400) as well as neutral species. By subtracting the spectrum of the central source from the total, we compare the PAH emission in the central vs extended region of a small sample of AGN. In contrast to the findings for the central regions of AGN-dominated systems, the PAH ratios measured in the extended regions of both type 1 and type 2 Seyfert galaxies can be explained assuming similar PAH molecular size distribution and ionized fractions of molecules to those seen in central regions of star-forming galaxies (100

We study the multi-phase feedback processes in the central ~3 kpc of the barred Sy 2 galaxy NGC 5643. We use observations of the cold molecular gas (ALMA CO(2-1)) and ionized gas (MUSE IFU). We study different regions along the outflow zone which extends out to ~2.3 kpc in the same direction (east-west) as the radio jet, as well as nuclear/circumnuclear regions in the host galaxy disk. The deprojected outflowing velocities of the cold molecular gas (median Vcentral~189 km s^-1) are generally lower than those of the outflowing ionized gas, which reach deprojected velocities of up to 750 km s^-1 close to the AGN, and their spatial profiles follow those of the ionized phase. This suggests that the outflowing molecular gas in the galaxy disk is being entrained by the AGN wind. We derive molecular and ionized outflow masses of ~5.2x10^7 Msun and 8.5x10^4 Msun and molecular and ionized outflow mass rates of ~51 Msun yr^-1 and 0.14 Msun yr^-1. Therefore, the molecular phase dominates the outflow mass and outflow mass rate, while the outflow kinetic power and momentum are similar in both phases. However, the wind momentum load for the molecular and ionized outflow phases are ~27-5 and <1, which suggests that the molecular phase is not momentum conserving while the ionized one most certainly is. The molecular gas content (~1.5x10^7 Msun) of the eastern spiral arm is approximately 50-70% of the content of the western one. We interpret this as destruction/clearing of the molecular gas produced by the AGN wind impacting in the host galaxy. The increase of the molecular phase momentum implies that part of the kinetic energy from the AGN wind is transmitted to the molecular outflow. This suggest that in Sy-like AGN such as NGC 5643, the radiative/quasar and the kinetic/radio AGN feedback modes coexist and may shape the host galaxies even at kpc-scales via both positive and (mild) negative feedback., Comment: 26 pages, 21 figures, 5 tables. Astronomy and Astrophysics, Accepted 2020 September 11, in press

We used high-spatial resolution (70 pc; 0.3") CO multi-transition (1-0, 2-1, 4-3, and 6-5) ALMA data to study the physical conditions and kinematics of the cold molecular outflow in the local LIRG ESO320-G030 (d=48 Mpc, log LIR/Lsun=11.3). ESO320-G030 is a double-barred isolated spiral, but its compact and obscured nuclear starburst (SFR~15 Msun/yr; Av~40 mag) resembles those of more luminous ULIRGs. In the outflow, the 1-0/2-1 ratio is enhanced with respect to the rest of the galaxy and the CO(4-3) transition is undetected. This indicates that the outflowing molecular gas is less excited than the gas in the nuclear starburst (launching site) and the galaxy disk. Non-LTE radiative transfer modeling reveals that the properties of the outflow molecular clouds differ from those of the nuclear and disk clouds: The kinetic temperature is lower (~9 K) in the outflow, and the outflowing clouds have lower column densities. Assuming a 10^-4 CO abundance, the large internal velocity gradients, 60^+250_-45 km/s/pc, imply that the outflowing molecular clouds are not bound by self-gravity. All this suggests that the life-cycle (formation, collapse, dissipation) of the disk clouds might differ from that of the outflowing clouds which might not be able to form stars. The low Tkin of the molecular outflow remains constant up to 1.7 kpc. This indicates that the heating by the hotter ionized outflow phase is not efficient and may favor the survival of the outflow molecular phase. The velocity structure of the outflow shows a 0.8 km/s/pc velocity gradient between 190-560 pc and then a constant maximum velocity (~750 km/s) up to 1.7 kpc. This is compatible with a pure gravitational evolution of the outflow under certain mass outflow rate and launching velocity variations. Alternatively, ram pressure acceleration and cloud evaporation could explain the observed kinematics and size of the molecular phase., Comment: Accepted for publication in A&A

We investigate the relation between the detection of the $11.3\,\mu$m PAH feature in the nuclear ($\sim 24-230\,$pc) regions of 22 nearby Seyfert galaxies and the properties of the cold molecular gas. For the former we use ground-based (0.3-0.6" resolution) mid-infrared (mid-IR) spectroscopy. The cold molecular gas is traced by ALMA and NOEMA high (0.2-1.1") angular resolution observations of the CO(2-1) transition. Galaxies with a nuclear detection of the $11.3\,\mu$m PAH feature contain more cold molecular gas (median $1.6\times 10^7\,M_\odot$) and have higher column densities ($N({\rm H}_2) = 2 \times 10^{23}\,{\rm cm}^{-2}$) over the regions sampled by the mid-IR slits than those without a detection. This suggests that molecular gas plays a role in shielding the PAH molecules in the harsh environments of Seyfert nuclei. Choosing the PAH molecule naphthalene as an illustration, we compute its half-life in the nuclear regions of our sample when exposed to 2.5keV hard X-ray photons. We estimate shorter half-lives for naphthalene in nuclei without a $11.3\,\mu$m PAH detection than in those with a detection. The Spitzer/IRS PAH ratios on circumnuclear scales ($\sim$ 4" $\sim$ 0.25-1.3kpc) are in between model predictions for neutral and partly ionized PAHs. However, Seyfert galaxies in our sample with the highest nuclear H$_2$ column densities are not generally closer to the neutral PAH tracks. This is because in the majority of our sample galaxies, the CO(2-1) emission in the inner $\sim$ 4" is not centrally peaked and in some galaxies traces circumnuclear sites of strong star formation activity. Spatially resolved observations with the MIRI medium-resolution spectrograph (MRS) on the James Webb Space Telescope will be able to distinguish the effects of an active galactic nucleus (AGN) and star formation on the PAH emission in nearby AGN., Comment: Accepted for publication in A&A

Obscured AGN represent a significant fraction of the entire AGN population, especially at high redshift (~70% at z=3--5). They are often characterized by the presence of large gas and dust reservoirs that are thought to sustain and possibly obscure vigorous star formation processes that make these objects shine at far-IR and sub-mm wavelengths. We exploit ALMA Cycle 4 observations of the continuum (~2.1mm) and high-J CO emission of a sample of six X-ray selected SMGs hosting an obscured AGN at z_spec>2.5 in the 7 Ms CDF-S. We measured the masses and sizes of the dust and molecular gas and we derived the gas density and column density on the basis of a uniform sphere geometry. Finally, we compared the measured column densities with those derived from the Chandra X-ray spectra. We detected both the continuum and line emission for three sources for which we measured both the flux density and size. For the undetected sources, we derived an upper limit on the flux density. We found that the detected galaxies are rich in gas and dust (molecular gas mass in the range <0.5 - 2.7 x 10^10 M_sun for {\alpha}_CO=0.8 and up to ~2 x 10^11~M_sun for {\alpha}_CO=6.5, and dust mass <0.9 - 4.9 x 10^8 M_sun) and compact (gas major axis 2.1-3.0 kpc, dust major axis 1.4-2.7 kpc). The column densities associated with the ISM are on the order of 10^(23-24) cm-2, which is comparable with those derived from the X-ray spectra. For the detected sources we also derived dynamical masses in the range 0.8 - 3.7 x 10^10 M_sun. We conclude that the ISM of high redshift galaxies can substantially contribute to nuclear obscuration up to the Compton-thick (>10^24 cm-2) regime. In addition, we found that all the detected sources show a velocity gradient reminding one rotating system, even though two of them show peculiar features in their morphology that can be associated with a chaotic, possibly merging, structure., Comment: 15 pages, 6 figures, 7 tables. Accepted for publication in A&A

AGN-driven outflows are believed to play an important role in regulating the growth of galaxies mostly via negative feedback. However, their effects on their hosts are far from clear, especially for low and moderate luminosity Seyferts. To investigate this issue, we have obtained cold molecular gas observations, traced by the CO(2-1) transition, using the NOEMA interferometer of five nearby (distances between 19 and 58 Mpc) Seyfert galaxies. The resolution of approx. 0.3-0.8 arcsec (approx. 30-100 pc) and field of view of NOEMA allowed us to study the CO(2-1) morphology and kinematics in the nuclear regions (approx. 100 pc) and up to radial distances of approx. 900 pc. We have detected CO(2-1) emission in all five galaxies with disky or circumnuclear ring like morphologies. We derived cold molecular gas masses on nuclear (approx. 100 pc) and circumnuclear (approx. 650 pc) scales in the range from $10^6$ to $10^7$M$_{\odot}$ and from $10^7$ to $10^8$ $M_{\odot}$, respectively. In all of our galaxies the bulk of this gas is rotating in the plane of the galaxy. However, non-circular motions are also present. In NGC 4253, NGC 4388 and NGC 7465, we can ascribe the streaming motions to the presence of a large-scale bar. In Mrk 1066 and NGC 4388, the non-circular motions in the nuclear regions are explained as outflowing material due to the interaction of the AGN wind with molecular gas in the galaxy disk. We conclude that for an unambiguous and precise interpretation of the kinematics of the cold molecular gas we need a detailed knowledge of the host galaxy (i.e., presence of bars, interactions, etc) as well as of the ionized gas kinematics and the ionization cone geometry., Comment: Accepted for publication in A&A, 23 pages, 17 figures

We present the discovery of HLock01-LAB, a luminous and large Lya nebula at z=3.326. Medium-band imaging and long-slit spectroscopic observations with the Gran Telescopio Canarias reveal extended emission in the Lya 1215\AA, CIV1550\AA, and HeII 1640\AA lines over ~100kpc, and a total luminosity L(Lya)=(6.4+/-0.1)x10^44 erg s^-1. HLock01-LAB presents an elongated morphology aligned with two faint radio sources contained within the central ~8kpc of the nebula. The radio structures are consistent to be faint radio jets or lobes of a central galaxy, whose spectrum shows nebular emission characteristic of a type-II active galactic nucleus (AGN). The continuum emission of the AGN at short wavelengths is, however, likely dominated by stellar emission of the host galaxy, for which we derive a stellar mass M* = 2.3x10^11 Msun. The detection of extended emission in CIV and CIII] indicates that the gas within the nebula is not primordial. Feedback may have enriched the halo at at least 50 kpc from the nuclear region. Using rest-frame UV emission-line diagnostics, we find that the gas in the nebula is likely heated by the AGN. Nevertheless, at the center of the nebula we find extreme emission line ratios of Lya/CIV~60 and Lya/HeII~80, one of the highest values measured to date, and well above the standard values of photoionization models (Lya/HeII~30 for case B photoionization). Our data suggest that jet-induced shocks are likely responsible for the increase of the electron temperature and, thus, the observed Lya enhancement in the center of the nebula. This scenario is further supported by the presence of radio structures and perturbed kinematics in this region. The large Lya luminosity in HLock01-LAB is likely due to a combination of AGN photoionization and jet-induced shocks, highlighting the diversity of sources of energy powering Lya nebulae. [abridged], Comment: 13 pages, 8 figures, 3 tables. Accepted for publication in A&A. Abstract abridged for arXiv submission

We present ALMA observations of the CO(2-1) and CO(3-2) molecular gas transitions and associated (sub)-mm continua of the nearby Seyfert 1.5 galaxy NGC3227 with angular resolutions 0.085-0.21" (7-15pc). On large scales the cold molecular gas shows circular motions as well as streaming motions on scales of a few hundred parsecs associated with a large scale bar. We fitted the nuclear ALMA 1.3mm emission with an unresolved component and an extended component. The 850$\mu$m emission shows at least two extended components, one along the major axis of the nuclear disk and the other along the axis of the ionization cone. The molecular gas in the central region (1" ~73pc) shows several CO clumps with complex kinematics which appears to be dominated by non-circular motions. While we cannot demonstrate conclusively the presence of a warped nuclear disk, we also detected non-circular motions along the kinematic minor axis. They reach line-of-sight velocities of v-vsys =150-200km/s. Assuming that the radial motions are in the plane of the galaxy, then we interpret them as a nuclear molecular outflow due to molecular gas in the host galaxy being entrained by the AGN wind. We derive molecular outflow rates of $5\,M_\odot\,{\rm yr}^{-1}$ and $0.6\,M_\odot\,{\rm yr}^{-1}$ at projected distances of up to 30pc to the northeast and southwest of the AGN, respectively. At the AGN location we estimate a mass in molecular gas of $5\times 10^{5}\,M_\odot$ and an average column density $N({\rm H}_2) = 2-3\times 10^{23}\,{\rm cm}^{-2}$ in the inner 15pc. The nuclear molecular gas and sub-mm continuum emission of NGC3227 do not resemble the classical compact torus. Rather, these emissions extend for several tens of parsecs and appear connected with the circumnuclear ring in the host galaxy disk, as found in other local AGN. (Abridged), Comment: Accepted for publication in Astronomy and Astrophysics, 17 pages, 12 figures

The extreme infrared (IR) luminosity of local luminous and ultra-luminous IR galaxies (U/LIRGs; 11 < log LIR /Lsun < 12 and log LIR /Lsun > 12, respectively) is mainly powered by star-formation processes triggered by mergers or interactions. While U/LIRGs are rare locally, at z > 1, they become more common, they dominate the star-formation rate (SFR) density, and a fraction of them are found to be normal disk galaxies. Therefore, there must be an evolution of the mechanism triggering these intense starbursts with redshift. To investigate this evolution, we present new optical SWIFT integral field spectroscopic H{\alpha}+[NII] observations of a sample of 9 intermediate-z (0.2 < z < 0.4) U/LIRG systems selected from Herschel 250{\mu}m observations. The main results are the following: (a) the ratios between the velocity dispersion and the rotation curve amplitude indicate that 10-25% (1-2 out of 8) might be compatible with being isolated disks while the remaining objects are interacting/merging systems; (b) the ratio between un-obscured and obscured SFR traced by H{\alpha} and LIR, respectively, is similar in both local and these intermediate-z U/LIRGs; and (c) the ratio between 250{\mu}m and the total IR luminosities of these intermediate-z U/LIRGs is higher than that of local U/LIRGs with the same LIR . This indicates a reduced dust temperature in these intermediate-z U/LIRGs. This, together with their already measured enhanced molecular gas content, suggests that the interstellar medium conditions are different in our sample of intermediate-z galaxies when compared to local U/LIRGs., Comment: Accepted for publication in MNRAS

[Abridged] We combine new CO(1-0) line observations of 24 intermediate redshift galaxies (0.03 < z < 0.28) along with literature data of galaxies at 0

We present new CO(2-1) observations of 3 low-z (~350 Mpc) ULIRG systems (6 nuclei) observed with ALMA at high-spatial resolution (~500 pc). We detect massive cold molecular gas outflows in 5 out of 6 nuclei (0.3-5)x10^8 Msun. These outflows are spatially resolved with deprojected radii of 0.25-1 kpc although high-velocity molecular gas is detected up to ~0.5-1.8 kpc (1-6 kpc deprojected). The mass outflow rates are 12-400 Msun/yr and the inclination corrected average velocity of the outflowing gas 350-550 km/s (v_max = 500-900 km/s). The origin of these outflows can be explained by the nuclear starbursts although the contribution of an obscured AGN can not be completely ruled out. The position angle (PA) of the outflowing gas along the kinematic minor axis of the nuclear molecular disk suggests that the outflow axis is perpendicular to the disk for three of these outflows. Only in one case, the outflow PA is clearly not along the kinematic minor axis. The outflow depletion times are 15-80 Myr which are slightly shorter than the star-formation (SF) depletion times (30-80 Myr). However, we estimate that only 15-30% of the outflowing gas will escape the gravitational potential of the nucleus. The majority of the outflowing gas will return to the disk after 5-10 Myr and become available to form new stars. Therefore, these outflows will not likely quench the nuclear starbursts. These outflows would be consistent with being driven by radiation pressure (momentum-driven) only if the coupling between radiation and dust increases with increasing SF rates. This can be achieved if the dust optical depth is higher in objects with higher SF. The relatively small sizes (<1 kpc) and dynamical times (<3 Myr) of the cold molecular outflows suggests that molecular gas cannot survive longer in the outflow environment or that it cannot form efficiently beyond these distances or times. (Abridged), Comment: Accepted for publication in A&A; 20 pages, 18 figures

We present ALMA Band 6 $^{12}$CO(2--1) line and rest-frame 232GHz continuum observations of the nearby Compton-thick Seyfert galaxy NGC5643 with angular resolutions 0.11-0.26arcsec (9-21pc). The CO(2--1) integrated line map reveals emission from the nuclear and circumnuclear region with a two-arm nuclear spiral extending 10arcsec on each side. The circumnuclear CO(2--1) kinematics can be fitted with a rotating disk, although there are regions with large residual velocities and/or velocity dispersions. The CO(2--1) line profiles of these regions show two different velocity components. One is ascribed to the circular component and the other to the interaction of the AGN outflow, as traced by the [O III]5007AA emission, with molecular gas in the disk a few hundred parsecs from the AGN. On nuclear scales, we detected an inclined CO(2--1) disk (diameter 26 pc, FWHM) oriented almost in a north-south direction. The CO(2--1) nuclear kinematics can be fitted with a rotating disk which appears to be tilted with respect to the large scale disk. There are strong non-circular motions in the central 0.2-0.3 arcsec with velocities of up to 110km/s. In the absence of a nuclear bar, these motions could be explained as radial outflows in the nuclear disk. We estimate a total molecular gas mass for the nuclear disk of $M({\rm H}_2)=1.1\times 10^7\,M_\odot$ and an H$_2$ column density toward the location of the AGN of $N({\rm H}_2)\sim 5 \times 10^{23}\,{\rm cm}^{-2}$, for a standard CO-to-H$_2$ conversion factor. We interpret this nuclear molecular gas disk as the obscuring torus of NGC5643 as well as the collimating structure of the ionization cone., Comment: 16 papes, 9 figures. Accepted to ApJ. Updated Fig.4

We present the Herschel-SPIRE photometric atlas for a complete flux limited sample of 43 local Ultraluminous Infrared Galaxies (ULIRGs), selected at 60${\mu}$m by IRAS, as part of the HERschel ULIRG Survey (HERUS). Photometry observations were obtained using the SPIRE instrument at 250, 350 and 500${\mu}$m. We describe these observations, present the results, and combine the new observations with data from IRAS to examine the far-IR spectral energy distributions (SEDs) of these sources. We fit the observed SEDs of HERUS objects with a simple parameterised modified black body model where temperature and emissivity $\beta$ are free parameters. We compare the fitted values to those of non-ULIRG local galaxies, and find, in agreement with earlier results, that HERUS ULIRGs have warmer dust (median temperature T = 37.9+/-4.7 K compared to 21.3+/-3.4 K) but a similar $\beta$ distribution (median $\beta$ = 1.7 compared to 1.8) to the Herschel reference sample (HRS, Cortese et al., 2014) galaxies. Dust masses are found to be in the range of 10^7.5 to 10^9 Msun significantly higher than that of Herschel Reference Sample (HRS) sources. We compare our results for local ULIRGs with higher redshift samples selected at 250 and 850${\mu}$m. These latter sources generally have cooler dust and/or redder 100-to-250 ${\mu}$m colours than our 60${\mu}$m-selected ULIRGs. We show that this difference may in part be the result of the sources being selected at different wavelengths rather than being a simple indication of rapid evolution in the properties of the population., Comment: Accepted for publication in MNRAS

In an earlier paper we modeled the far-infrared emission from a star-forming galaxy using the photoionisation code CLOUDY and presented metallicity sensitive diagnostics based on far-infrared fine structure line ratios. Here, we focus on the applicability of the [OIII]88/[NII]122 microns line ratio as a gas phase metallicity indicator in high redshift submillimetre luminous galaxies. The [OIII]88/[NII]122 microns ratio is strongly dependent on the ionization parameter (which is related to the total number of ionizing photons) as well as the gas electron density. We demonstrate how the ratio of 88/$122 continuum flux measurements can provide a reasonable estimate of the ionization parameter while the availability of the [NII]205 microns line can constrain the electron density. Using the [OIII]88/[NII]122 microns line ratios from a sample of nearby normal and star-forming galaxies we measure their gas phase metallicities and find that their mass metallicity relation is consistent with the one derived using optical emission lines. Using new, previously unpublished, Herschel spectroscopic observations of key far-infrared fine structure lines of the z~3 galaxy HLSW-01 and additional published measurements of far-infrared fine structure lines of high-z submillimetre luminous galaxies we derive gas phase metallicities using their [OIII]88/[NII]122 microns line ratio. We find that the metallicities of these z~3 submm luminous galaxies are consistent with solar metallicities and that they appear to follow the mass-metallicity relation expected for z~3 systems., Comment: 10 pages, 7 figures, MNRAS in press

The abundance of metals in galaxies is a key parameter which permits to distinguish between different galaxy formation and evolution models. Most of the metallicity determinations are based on optical line ratios. However, the optical spectral range is subject to dust extinction and, for high-z objects (z > 3), some of the lines used in optical metallicity diagnostics are shifted to wavelengths not accessible to ground based observatories. For this reason, we explore metallicity diagnostics using far-infrared (IR) line ratios which can provide a suitable alternative in such situations. To investigate these far-IR line ratios, we modeled the emission of a starburst with the photoionization code CLOUDY. The most sensitive far-IR ratios to measure metallicities are the [OIII]52$\mu$m and 88$\mu$m to [NIII]57$\mu$m ratios. We show that this ratio produces robust metallicities in the presence of an AGN and is insensitive to changes in the age of the ionizing stellar. Another metallicity sensitive ratio is the [OIII]88$\mu$m/[NII]122$\mu$m ratio, although it depends on the ionization parameter. We propose various mid- and far-IR line ratios to break this dependency. Finally, we apply these far-IR diagnostics to a sample of 19 local ultraluminous IR galaxies (ULIRGs) observed with Herschel and Spitzer. We find that the gas-phase metallicity in these local ULIRGs is in the range 0.7 < Z_gas/Z_sun < 1.5, which corresponds to 8.5 < 12 + log (O/H) < 8.9. The inferred metallicities agree well with previous estimates for local ULIRGs and this confirms that they lie below the local mass-metallicity relation., Comment: Accepted for publication in MNRAS; 25 pages, 19 figures

We present millimetre dust emission measurements of two Lyman Break Galaxies at z~3 and construct for the first time fully sampled infrared spectral energy distributions (SEDs), from mid-IR to the Rayleigh-Jeans tail, of individually detected, unlensed, UV-selected, main sequence (MS) galaxies at $z=3$. The SED modelling of the two sources confirms previous findings, based on stacked ensembles, of an increasing mean radiation field with redshift, consistent with a rapidly decreasing gas metallicity in z > 2 galaxies. Complementing our study with CO[3-2] emission line observations, we measure the molecular gas mass (M_H2) reservoir of the systems using three independent approaches: 1) CO line observations, 2) the dust to gas mass ratio vs metallicity relation and 3) a single band, dust emission flux on the Rayleigh-Jeans side of the SED. All techniques return consistent M_H2 estimates within a factor of ~2 or less, yielding gas depletion time-scales (tau_dep ~ 0.35 Gyrs) and gas-to-stellar mass ratios (M_H2/M* ~ 0.5-1) for our z~3 massive MS galaxies. The overall properties of our galaxies are consistent with trends and relations established at lower redshifts, extending the apparent uniformity of star-forming galaxies over the last 11.5 billion years., Comment: Accepted for publication at A&A