Showing 1–50 of 209 results for author: Rigopoulou, D

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

Authors: L. Hermosa Muñoz, A. Alonso-Herrero, M. Pereira-Santaella, I. García-Bernete, S. García-Burillo, B. García-Lorenzo, R. Davies, T. Shimizu, D. Esparza-Arredondo, E. K. S. Hicks, H. Haidar, M. Leist, E. López-Rodríguez, C. Ramos Almeida, D. Rosario, L. Zhang, A. Audibert, E. Bellocchi, P. Boorman, A. J. Bunker, F. Combes, S. Campbell, T. Díaz-Santos, L. Fuller, P. Gandhi , et al. (13 additional authors not shown)

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

Submitted 22 July, 2024; originally announced July 2024.

Comments: 20 pages, 6 figures. Accepted for publication in A&A. Abstract adapted for the arxiv version

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

Polycyclic Aromatic Hydrocarbon Emission in Galaxies as seen with JWST

Authors: Dimitra Rigopoulou, Fergus R. Donnan, Ismael García-Bernete, Miguel Pereira-Santaella, Almudena Alonso-Herrero, Ric Davies, Leslie K. Hunt, Patrick F. Roche, Taro Shimizu

Abstract: We present a systematic study of mid-infrared spectra of galaxies including star-forming galaxies and Active Galactic Nuclei observed with JWST MIRI-MRS and NIRSpec-IFU. We focus on the relative variations of the 3.3, 6.2, 7.7, 11.3, 12.7 and 17 micron Polycyclic Aromatic Hydrocarbon (PAH) features within spatially resolved regions of galaxies including NGC 3256, NGC 7469, VV 114, II Zw96 and NGC… ▽ More We present a systematic study of mid-infrared spectra of galaxies including star-forming galaxies and Active Galactic Nuclei observed with JWST MIRI-MRS and NIRSpec-IFU. We focus on the relative variations of the 3.3, 6.2, 7.7, 11.3, 12.7 and 17 micron Polycyclic Aromatic Hydrocarbon (PAH) features within spatially resolved regions of galaxies including NGC 3256, NGC 7469, VV 114, II Zw96 and NGC 5728. Using theoretical PAH models and extending our earlier work, we introduce a new PAH diagnostic involving the 17 micron PAH feature. To determine the drivers of PAH band variations in galaxies, we compare observed PAH spectral bands to predictions from theoretical PAH models. We consider extinction, dehydrogenation and PAH size and charge as possible drivers of PAH band variations. We find a surprising uniformity in PAH size distribution among the spatially resolved regions of the galaxies studied here, with no evidence for preferential destruction of the smallest grains, contrary to earlier findings. Neither extinction nor dehydrogenation play a crucial role in setting the observed PAH bands. Instead, we find that PAH charge plays a significant role in PAH inter-band variations. We find a tight relation between PAH charge and the intensity of the radiation field as traced by the [NeIII]/[NeII] maps. In agreement with recent JWST results, we find a predominance of neutral PAH molecules in the nuclei of Active Galaxies and their outflows. Ionised PAHs are the dominant population in star-forming galaxies. We discuss the implications of our findings for the use of PAHs as ISM tracers in high redshift galaxies. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 16 pages, 13 figures, accepted for publication by MNRAS

The fountain of the luminous infrared galaxy Zw049.057 as traced by its OH megamaser

Authors: Boy Lankhaar, Susanne Aalto, Clare Wethers, Javier Moldon, Rob Beswick, Mark Gorski, Sabine König, Chentao Yang, Jeff Mangum, John Gallagher, Francoise Combes, Dimitra Rigopoulou, Eduardo González-Alfonso, Sébastien Muller, Ismael Garcia-Bernete, Christian Henkel, Yuri Nishimura, Claudio Ricci

Abstract: High resolution (0."037-0."13 [10-35 pc]) e-MERLIN ($\lambda6-18$ cm) and (0."024 [6.5 pc]) ALMA ($λ1.1$ mm) observations have been used to image OH (hydroxyl) and H$_2$CO (formaldehyde) megamaser emission, and HCN 3->2 emission towards the nuclear (<100 pc) region of the luminous infrared galaxy Zw049.057. Zw049.057 hosts a compact obscured nucleus (CON), thus representing a class of galaxies tha… ▽ More High resolution (0."037-0."13 [10-35 pc]) e-MERLIN ($\lambda6-18$ cm) and (0."024 [6.5 pc]) ALMA ($λ1.1$ mm) observations have been used to image OH (hydroxyl) and H$_2$CO (formaldehyde) megamaser emission, and HCN 3->2 emission towards the nuclear (<100 pc) region of the luminous infrared galaxy Zw049.057. Zw049.057 hosts a compact obscured nucleus (CON), thus representing a class of galaxies that are often associated with inflow and outflow motions. Formaldehyde megamaser emission is detected towards the nuclear region, <30 pc (<0."1), and traces a structure along the disk major axis. OH megamaser (OHM) emission is detected along the minor axis of the disk, ~30 pc (0."1) from the nucleus, where it exhibits a velocity gradient with extrema of -20 km/s south-east (SE) of the disk and -110 km/s north-west (NW) of the disk. HCN 3->2 emission reveals extended emission, along the disk minor axis out to ~60 pc (0."2). Analysis of the minor axis HCN emission reveals high-velocity features, extending out to 600 km/s, redshifted on the SE side and blueshifted on the NW side. We propose that the high-velocity HCN emission traces a fast >250 km/s and collimated outflow, that is enveloped by a wide-angle and slow ~50 km/s outflow that is traced by the OHM emission. Analysis of the outflow kinematics suggests that the slow wide-angle outflow will not reach escape velocity and instead will fall back to the galaxy disk, evolving as a so-called fountain flow, while the fast collimated outflow traced by HCN emission will likely escape the nuclear region. We suggest that the absence of OHM emission in the nuclear region is due to high densities there. Even though OHMs associated with outflows are an exception to conventional OHM emission, we expect them to be common in CON sources that host both OHM and H$_2$CO megamasers. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: 13 pages, 3 figures. Accepted to A&A

The Kinematics of Polycyclic Aromatic Hydrocarbons (PAHs) in Galaxies revealed by Principal Component Analysis (PCA) tomography with JWST/NIRSpec

Authors: Fergus R. Donnan, Dimitra Rigopoulou, Ismael García-Bernete

Abstract: Polycyclic Aromatic Hydrocarbons (PAHs) are organic molecules which comprise the smallest particles of dust in the interstellar medium (ISM). Due to their broad/complex emission profiles, obtaining kinematics is a challenge with traditional methods, especially before the advent of the JWST. In this work we employ Principal Component Analysis (PCA) tomography to analyse JWST/NIRSpec IFU data of 3 n… ▽ More Polycyclic Aromatic Hydrocarbons (PAHs) are organic molecules which comprise the smallest particles of dust in the interstellar medium (ISM). Due to their broad/complex emission profiles, obtaining kinematics is a challenge with traditional methods, especially before the advent of the JWST. In this work we employ Principal Component Analysis (PCA) tomography to analyse JWST/NIRSpec IFU data of 3 nearby luminous infrared galaxies (LIRGs), namely, NGC 3256 N, NGC 3256 S and NGC 7469. We detect the signature of rotation in the second principal component of the 3.3 $μ$m PAH feature in all three targets. We construct velocity maps from the principal components for the 3.3 $μ$m PAH feature, Br$β$ (2.625$μ$m) and molecular hydrogen, H$_2$ 1-0 S(1) (2.12$μ$m). We find that in each target, the PAHs qualitatively follow the rotation of the galaxy, consistent with the rotational signature derived from Br$β$ and H$_2$. There are however some differences, with the PAH rotation in NGC 3256 N appearing at a different position angle, which suggest differences in the motion of the dust as compared to the gas. This kind of analysis opens a new window into this key component of the ISM. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Accepted for publication in MNRAS Letters. Online supplementary material is included

Single Aperture Large Telescope for Universe Studies (SALTUS): Science Overview

Authors: Gordon Chin, Carrie M. Anderson, Jennifer Bergner, Nicolas Biver, Gordon L. Bjoraker, Thibault Cavalie, Michael DiSanti, Jian-Rong Gao, Paul Hartogh, Leon K. Harding, Qing Hu, Daewook Kim, Craig Kulesa, Gert de Lange, David T. Leisawitz, Rebecca C. Levy, Arthur Lichtenberger, Daniel P. Marronh, Joan Najita, Trent Newswander, George H. Rieke, Dimitra Rigopoulou, Peter Roefsema, Nathan X. Roth, Kamber Schwarz , et al. (11 additional authors not shown)

Abstract: The SALTUS Probe mission will provide a powerful far-infrared (far-IR) pointed space observatory to explore our cosmic origins and the possibility of life elsewhere. The observatory employs an innovative deployable 14-m aperture, with a sunshield that will radiatively cool the off-axis primary to <45K. This cooled primary reflector works in tandem with cryogenic coherent and incoherent instruments… ▽ More The SALTUS Probe mission will provide a powerful far-infrared (far-IR) pointed space observatory to explore our cosmic origins and the possibility of life elsewhere. The observatory employs an innovative deployable 14-m aperture, with a sunshield that will radiatively cool the off-axis primary to <45K. This cooled primary reflector works in tandem with cryogenic coherent and incoherent instruments that span the 34 to 660 micron far-IR range at both high and moderate spectral resolutions. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 49 pages, 10 figures, 3 tables, submitted to SPIE JATIS

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!

Double, double, toil, and trouble: The tails, bubbles, and knots of the local compact obscured nucleus galaxy NGC4418

Authors: C. F. Wethers, S. Aalto, G. C. Privon, F. Stanley, J. Gallagher, M. Gorski, S. König, K. Onishi, M. Sato, C. Yang, R. Beswick, L. Barcos-Munoz F. Combes, T. Diaz-Santos, A. S. Evans, I. Garcia-Bernete, C. Henkel, M. Imanishi, S. Martín, S. Muller, Y. Nishimura, C. Ricci, D. Rigopoulou, S. Viti

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

Submitted 27 February, 2024; originally announced February 2024.

Comments: 17 pages, 15 figures, accepted to A&A

Peeling Back the Layers of Extinction of Dusty Galaxies in the Era of JWST: Modelling Joint NIRSpec + MIRI Spectra at rest-frame 1.5-28 $μ$m

Authors: F. R. Donnan, I. García-Bernete, D. Rigopoulou, M. Pereira-Santaella, P. F. Roche, A. Alonso-Herrero

Abstract: We present an analysis of the combined NIRSpec and MIRI spectra of dusty galaxies between 1.5 - 28 $μ$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 na… ▽ More We present an analysis of the combined NIRSpec and MIRI spectra of dusty galaxies between 1.5 - 28 $μ$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. △ Less

Submitted 27 February, 2024; originally announced February 2024.

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

Structures Of Dust and gAs (SODA): Constraining the innermost dust properties of II Zw96 with JWST observations of H$_2$O and CO

Authors: I. García-Bernete, M. Pereira-Santaella, E. González-Alfonso, D. Rigopoulou, A. Efstathiou, F. R. Donnan, N. Thatte

Abstract: 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 $ν_2$=1-0 ($\sim$5.3-7.2 $μ$m) and $^{12}$CO $ν$=1-0 ($\sim$4.45-4.95 $μ$m) in D1. The CO band shows the R- and P-branches in emission and the spectr… ▽ More 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 $ν_2$=1-0 ($\sim$5.3-7.2 $μ$m) and $^{12}$CO $ν$=1-0 ($\sim$4.45-4.95 $μ$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 $μ$m polycyclic aromatic hydrocarbon band in this source. By including the 4.5-7.0 $μ$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 $μ$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. △ Less

Submitted 5 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A (official date of acceptance: 27/12/2023; submitted 27/11/2023). 11 pages, 10 figures

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

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

Uncovering the MIR emission of quiescent galaxies with $JWST$

Authors: David Blánquez-Sesé, G. E. Magdis, C. Gómez-Guijarro, M. Shuntov, V. Kokorev, G. Brammer, F. Valentino, T. Díaz-Santos, E. -D. Paspaliaris, D. Rigopoulou, J. Hjorth, D. Langeroodi, R. Gobat, S. Jin, N. B. Sillassen, S. Gillman, T. R. Greve, M. Lee

Abstract: 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<z<0.7$) and perform modeling of their rest-frame UV to MIR photometry. We find that QGs exhibit a range of MIR spectra that are composed of a… ▽ More 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<z<0.7$) and perform modeling of their rest-frame UV to MIR photometry. We find that QGs exhibit a range of MIR spectra that are composed of a stellar continuum and a dust component that is 1-2 orders of magnitude fainter to that of star-forming galaxies. The observed scatter in the MIR spectra, especially at $λ_{\rm rest} > 5 μ$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 $μ$m PAHs strengths with respect to those at 6.2- and 7.7$ μ$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$. △ Less

Submitted 2 October, 2023; originally announced October 2023.

Comments: The panchromatic QG SED has been made publicly available at http://www.georgiosmagdis.com/software

Journal ref: A&A 679, L2 (2023)

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

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

Molecular gas content and high excitation of a massive main-sequence galaxy at z = 3

Authors: Han Lei, Francesco Valentino, Georgios E. Magdis, Vasily Kokorev, Daizhong Liu, Dimitra Rigopoulou, Shuowen Jin, Emanuele Daddi

Abstract: We present new CO ($J=5-4$ and $7-6$) and [CI] ($^3P_2\,-\, ^3P_1$ and $^3P_1\,-\, ^3P_0$) emission line observations of the star-forming galaxy D49 at the massive end of the Main Sequence at $z=3$. We incorporate previous CO ($J=3-2$) and optical-to-millimetre continuum observations to fit its spectral energy distribution (SED). Our results hint at high-$J$ CO luminosities exceeding the expected… ▽ More We present new CO ($J=5-4$ and $7-6$) and [CI] ($^3P_2\,-\, ^3P_1$ and $^3P_1\,-\, ^3P_0$) emission line observations of the star-forming galaxy D49 at the massive end of the Main Sequence at $z=3$. We incorporate previous CO ($J=3-2$) and optical-to-millimetre continuum observations to fit its spectral energy distribution (SED). Our results hint at high-$J$ CO luminosities exceeding the expected location on the empirical correlations with the infrared luminosity. [CI] emission fully consistent with the literature trends is found. We do not retrieve any signatures of a bright active galactic nucleus that could boost the $J=5-4,\,7-6$ lines in either the infrared or X-ray bands, but warm photon-dominated regions, shocks or turbulence could in principle do so. We suggest that mechanical heating could be a favourable mechanism able to enhance the gas emission at fixed infrared luminosity in D49 and other main-sequence star-forming galaxies at high redshift, but further investigation is necessary to confirm this explanation. We derive molecular gas masses from dust, CO, and [CI] that all agree within the uncertainties. Given its large star formation rate (SFR) $\sim 500~M_\odot~{\rm yr}^{-1}$ and stellar mass $>10^{11.5}~M_\odot$, the short depletion time scale of $<0.3$ Gyr might indicate that D49 is experiencing its last growth spurt and will soon transit to quiescence. △ Less

Submitted 28 April, 2023; originally announced May 2023.

Comments: 7 pages, 4 figures, accepted for publication in A&A Letter

Journal ref: A&A 673, L13 (2023)

A Diverse Population of z ~ 2 ULIRGs Revealed by JWST Imaging

Authors: J. -S. Huang, Zi-Jian Li, Cheng Cheng, Meicun Hou, Haojing Yan, S. P. Willner, Y. -S. Dai, X. Z. Zheng, J. Pan, D. Rigopoulou, T. Wang, Zhiyuan Li, Piaoran Liang, A. Esamdin, G. G. Fazio

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

Submitted 6 April, 2023; v1 submitted 3 April, 2023; originally announced April 2023.

Comments: 9 pages, 3 figures, 1 table, accepted for publication in ApJ. V2 updates author affiliations and acknowledgments, not scientific content

Absence of nuclear polycyclic aromatic hydrocarbon emission from a compact starburst: The case of the type-2 quasar Mrk 477

Authors: C. Ramos Almeida, D. Esparza-Arredondo, O. Gonzalez-Martin, I. Garcia-Bernete, M. Pereira-Santaella, A. Alonso-Herrero, J. A. Acosta-Pulido, P. S. Bessiere, N. A. Levenson, C. N. Tadhunter, D. Rigopoulou, M. Martinez-Paredes, S. Cazzoli, B. Garcia-Lorenzo

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

Submitted 3 January, 2023; v1 submitted 2 December, 2022; originally announced December 2022.

Comments: 7 pages, 3 figures. A&A Letters accepted version including proof corrections

Journal ref: A&A 669, L5 (2023)

A Detailed Look at the Most Obscured Galactic Nuclei in the Mid-Infrared

Authors: F. R. Donnan, D. Rigopoulou, I. García-Bernete, M. Pereira-Santaella, A. Alonso-Herrero, P. F. Roche, S. Aalto, A. Hernán-Caballero, H. W. W. Spoon

Abstract: 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 Aromat… ▽ More 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, $τ_N$, of the $9.8 μ$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, $τ_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 $τ_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. △ Less

Submitted 17 November, 2022; originally announced November 2022.

Comments: Accepted for publication in A&A. 23 pages, 15 figures, 3 tables

Journal ref: A&A 669, A87 (2023)

A solar metallicity galaxy at $z >$ 7? Possible detection of the [N II] 122 $μ$m and [O III] 52 $μ$m lines

Authors: Meghana Killi, Darach Watson, Seiji Fujimoto, Hollis Akins, Kirsten Knudsen, Johan Richard, Yuichi Harikane, Dimitra Rigopoulou, Francesca Rizzo, Michele Ginolfi, Gergö Popping, Vasily Kokorev

Abstract: We present the first detection of the [N II] 122 $μ$m and [O III] 52 $μ$m lines for a reionisation-epoch galaxy. Based on these lines and previous [C II] 158 $μ$m and [O III] 88 $μ$m measurements, using two different radiative transfer models of the interstellar medium, we estimate an upper limit on electron density of $\lesssim$ 500 cm$^{-3}$ and a gas-phase metallicity… ▽ More We present the first detection of the [N II] 122 $μ$m and [O III] 52 $μ$m lines for a reionisation-epoch galaxy. Based on these lines and previous [C II] 158 $μ$m and [O III] 88 $μ$m measurements, using two different radiative transfer models of the interstellar medium, we estimate an upper limit on electron density of $\lesssim$ 500 cm$^{-3}$ and a gas-phase metallicity $Z/Z_\odot \sim 1.1 \pm 0.2$ for A1689-zD1, a gravitationally-lensed, dusty galaxy at $z$ = 7.133. Other measurements or indicators of metallicity so far in galaxy interstellar media at $z \gtrsim$ 6 are typically an order of magnitude lower than this. The unusually high metallicity makes A1689-zD1 inconsistent with the fundamental metallicity relation, although there is likely significant dust obscuration of the stellar mass, which may partly resolve the inconsistency. Given a solar metallicity, the dust-to-metals ratio is a factor of several lower than expected, hinting that galaxies beyond $z \sim$ 7 may have lower dust formation efficiency. Finally, the inferred nitrogen enrichment compared to oxygen, on which the metallicity measurement depends, indicates that star-formation in the system is older than about 250 Myr, pushing the beginnings of this galaxy to $z >$ 10. △ Less

Submitted 8 March, 2023; v1 submitted 2 November, 2022; originally announced November 2022.

Comments: 10 pages, 7 figures; Published MNRAS: 06 March 2023

The Obscured Nucleus and Shocked Environment of VV 114E Revealed by JWST/MIRI Spectroscopy

Authors: F. R. Donnan, I. García-Bernete, D. Rigopoulou, M. Pereira-Santaella, A. Alonso-Herrero, P. F. Roche, A. Hernán-Caballero, H. W. W. Spoon

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

Submitted 15 December, 2022; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 15 pages, 15 figures. Accepted for publication in MNRAS

Torus and polar dust dependence on AGN properties

Authors: I. García-Bernete, O. González-Martín, C. Ramos Almeida, A. Alonso-Herrero, M. Martínez-Paredes, M. J. Ward, P. F. Roche, J. A. Acosta-Pulido, E. López-Rodríguez, D. Rigopoulou, D. Esparza-Arredondo

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

Submitted 7 October, 2022; originally announced October 2022.

Comments: Accepted for publication in A&A, 43 pages, 43 figures, 16 tables

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

Probing computational methodologies in predicting mid-infrared spectra for large polycyclic aromatic hydrocarbons

Authors: Boutheïna Kerkeni, Ismael García-Bernete, Dimitra Rigopoulou, David P. Tew, Patrick F. Roche, David C. Clary

Abstract: We extend the prediction of vibrational spectra to large sized polycyclic aromatic hydrocarbon (PAH) molecules comprising up to \sim 1500 carbon atoms by evaluating the efficiency of several computational chemistry methodologies. We employ classical mechanics methods (Amber and Gaff) with improved atomic point charges, semi-empirical (PM3, and density functional tight binding), and density functio… ▽ More We extend the prediction of vibrational spectra to large sized polycyclic aromatic hydrocarbon (PAH) molecules comprising up to \sim 1500 carbon atoms by evaluating the efficiency of several computational chemistry methodologies. We employ classical mechanics methods (Amber and Gaff) with improved atomic point charges, semi-empirical (PM3, and density functional tight binding), and density functional theory (B3LYP) and conduct global optimizations and frequency calculations in order to investigate the impact of PAH size on the vibrational band positions. We primarily focus on the following mid-infrared emission bands 3.3, 6.2, 7.7, 8.6, 11.3, 12.7, and 17.0 microns. We developed a general Frequency Scaling Function (FSF) to shift the bands and to provide a systematic comparison versus the three methods for each PAH. We first validate this procedure on IR scaled spectra from the NASA Ames PAH Database, and extend it to new large PAHs. We show that when the FSF is applied to the Amber and Gaff IR spectra, an agreement between the normal mode peak positions with those inferred from the B3LYP/4-31G model chemistry is achieved. As calculations become time intensive for large sized molecules Nc > 450, this proposed methodology has advantages. The FSF has enabled extending the investigations to large PAHs where we clearly see the emergence of the 17.0 microns feature, and the weakening of the 3.3 microns one. We finally investigate the trends in the 3.3 microns 17.0 microns PAH band ratio as a function of PAH size and its response following the exposure to fields of varying radiation intensities. △ Less

Submitted 26 September, 2022; originally announced October 2022.

Comments: 19 pages, 12 figures

Journal ref: 513, 2022, 3663

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)

Stellar and black hole assembly in z<0.3 infrared-luminous mergers: intermittent starbursts vs. super-Eddington accretion

Authors: Duncan Farrah, Andreas Efstathiou, Jose Afonso, Jeronimo Bernard-Salas, Joe Cairns, David L Clements, Kevin Croker, Evanthia Hatziminaoglou, Maya Joyce, Mark Lacy, Vianney Lebouteiller, Alix Lieblich, Carol Lonsdale, Seb Oliver, Chris Pearson, Sara Petty, Lura K Pitchford, Dimitra Rigopoulou, Michael Rowan-Robinson, Jack Runburg, Henrik Spoon, Aprajita Verma, Lingyu Wang

Abstract: We study stellar and black hole mass assembly in a sample of 42 infrared-luminous galaxy mergers at z<0.3 by combining results from radiative transfer modelling with archival measures of molecular gas and black hole mass. The ratios of stellar mass, molecular gas mass, and black hole mass to each other are consistent with those of massive gas-rich galaxies at z<0.3. The advanced mergers may show i… ▽ More We study stellar and black hole mass assembly in a sample of 42 infrared-luminous galaxy mergers at z<0.3 by combining results from radiative transfer modelling with archival measures of molecular gas and black hole mass. The ratios of stellar mass, molecular gas mass, and black hole mass to each other are consistent with those of massive gas-rich galaxies at z<0.3. The advanced mergers may show increased black hole mass to stellar mass ratios, consistent with the transition from AGN to ellipticals and implying substantial black hole mass growth over the course of the merger. Star formation rates are enhanced relative to the local main sequence, by factors of ~100 in the starburst and ~1.8 in the host. The starburst star formation rates appear distinct to star formation in the main sequence at all redshifts up to at least z~5. Starbursts may prefer late-stage mergers, but are observed at any merger stage. We do not find evidence that the starbursts in these low-redshift systems substantially increase the total stellar mass, with a soft upper limit on the stellar mass increase from starburst activity of about a factor of two. In contrast, 12 objects show evidence for super-Eddington accretion, associated with late-stage mergers, suggesting that many AGN in infrared-luminous mergers go through a super-Eddington phase. The super-Eddington phase may increase black hole mass by up to an order of magnitude at an accretion efficiency of 42+/-33% over a period of 44+/-22Myr. Our results imply that super-Eddington accretion is an important black hole growth channel in infrared-luminous galaxies at all redshifts. △ Less

Submitted 29 April, 2022; originally announced May 2022.

Comments: MNRAS Accepted

A new look at local ultraluminous infrared galaxies: the atlas and radiative transfer models of their complex physics

Authors: A. Efstathiou, D. Farrah, J. Afonso, D. L. Clements, E. González-Alfonso, M. Lacy, S. Oliver, V. Papadopoulou Lesta, C. Pearson, D. Rigopoulou, M. Rowan-Robinson, H. W. W. Spoon, A. Verma, L. Wang

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

Submitted 29 April, 2022; originally announced May 2022.

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

Journal ref: 2022MNRAS.512.5183E

On the Viability of Determining Galaxy Properties from Observations I: Star Formation Rates and Kinematics

Authors: Kearn Grisdale, Laurence Hogan, Dimitra Rigopoulou, Niranjan Thatte, Miguel Pereira-Santaella, Julien Devriendt, Adrianne Slyz, Ismael García-Bernete, Yohan Dubois, Sukyoung K. Yi, Katarina Kraljic

Abstract: We explore how observations relate to the physical properties of the emitting galaxies by post-processing a pair of merging $z\sim2$ galaxies from the cosmological, hydrodynamical simulation NewHorizon using LCARS (Light from Cloudy Added to RAMSES) to encode the physical properties of the simulated galaxy into H$α$ emission line. By carrying out mock observations and analysis on these data cubes… ▽ More We explore how observations relate to the physical properties of the emitting galaxies by post-processing a pair of merging $z\sim2$ galaxies from the cosmological, hydrodynamical simulation NewHorizon using LCARS (Light from Cloudy Added to RAMSES) to encode the physical properties of the simulated galaxy into H$α$ emission line. By carrying out mock observations and analysis on these data cubes we ascertain which physical properties of the galaxy will be recoverable with the HARMONI spectrograph on the European Extremely Large Telescope (ELT). We are able to estimate the galaxy's star formation rate and dynamical mass to a reasonable degree of accuracy, with values within a factor of $1.81$ and $1.38$ of the true value. The kinematic structure of the galaxy is also recovered in mock observations. Furthermore, we are able to recover radial profiles of the velocity dispersion and are therefore able to calculate how the dynamical ratio varies as a function of distance from the galaxy centre. Finally, we show that when calculated on galaxy scales the dynamical ratio does not always provide a reliable measure of a galaxy's stability against gravity or act as an indicator of a minor merger. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Comments: 20 pages, 17 figures, Accepted for Publication in MNRAS, Comments and questions welcome

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)

A technique to select the most obscured galaxy nuclei

Authors: I. García-Bernete, D. Rigopoulou, S. Aalto, H. W. W. Spoon, A. Hernán-Caballero, A. Efstathiou, P. F. Roche, S. König

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

Submitted 28 March, 2022; originally announced March 2022.

Comments: 13 pages, 11 figures, 1 table, accepted for publication in A&A

Journal ref: A&A 663, A46 (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

Gas Phase Metallicities of Local Ultra-Luminous Infrared Galaxies Follow Normal Star-Forming Galaxies

Authors: Nima Chartab, Asantha Cooray, Jingzhe Ma, Hooshang Nayyeri, Preston Zilliot, Jonathan Lopez, Dario Fadda, Rodrigo Herrera-Camus, Matthew Malkan, Dimitra Rigopoulou, Kartik Sheth, Julie Wardlow

Abstract: Despite advances in observational data, theoretical models, and computational techniques to simulate key physical processes in the formation and evolution of galaxies, the stellar mass assembly of galaxies still remains an unsolved problem today. Optical spectroscopic measurements appear to show that the gas-phase metallicities of local ultra-luminous infrared galaxies (ULIRGs) are significantly l… ▽ More Despite advances in observational data, theoretical models, and computational techniques to simulate key physical processes in the formation and evolution of galaxies, the stellar mass assembly of galaxies still remains an unsolved problem today. Optical spectroscopic measurements appear to show that the gas-phase metallicities of local ultra-luminous infrared galaxies (ULIRGs) are significantly lower than those of normal star-forming galaxies. This difference has resulted in the claim that ULIRGs are fueled by metal-poor gas accretion from the outskirts\cite{Mannucci10}. Here we report on a new set of gas-phase metallicity measurements making use of the far-infrared spectral lines of [O{\sc iii}]52 $μ$m, [O{\sc iii}]88 $μ$m, and [N{\sc iii}]57 $μ$m instead of the usual optical lines. Photoionization models have resulted in a metallicity diagnostic based on these three lines that break the electron density degeneracy and reduce the scatter of the correlation significantly. Using new data from SOFIA and archival data from Herschel Space Observatory, we find that local ULIRGs lie on the mass-metallicity relation of star-forming galaxies and have metallicities comparable to other galaxies with similar stellar masses and star formation rates. The lack of a departure suggests that ULIRGs follow the same mass assembly mechanism as luminous star-forming galaxies and $\sim 0.3$ dex under-abundance in metallicities derived from optical lines is a result of heavily obscured metal-rich gas which has a negligible effect when using the FIR line diagnostics. △ Less

Submitted 19 January, 2022; originally announced January 2022.

Comments: Preprint submitted for publication 10/29/21

The HASHTAG project: The First Submillimeter Images of the Andromeda Galaxy from the Ground

Authors: Matthew W. L. Smith, Stephen A. Eales, Thomas G. Williams, Bumhyun Lee, Zongnan Li, Pauline Barmby, Martin Bureau, Scott Chapman, Brian S. Cho, Aeree Chung, Eun Jung Chung, Hui-Hsuan Chung, Christopher J. R. Clark, David L. Clements, Timothy A. Davis, Ilse De Looze, David J. Eden, Gayathri Athikkat-Eknath, George P. Ford, Yu Gao, Walter Gear, Haley L. Gomez, Richard de Grijs, Jinhua He, Luis C. Ho , et al. (24 additional authors not shown)

Abstract: Observing nearby galaxies with submillimeter telescopes on the ground has two major challenges. First, the brightness is significantly reduced at long submillimeter wavelengths compared to the brightness at the peak of the dust emission. Second, it is necessary to use a high-pass spatial filter to remove atmospheric noise on large angular scales, which has the unwelcome by-product of also removing… ▽ More Observing nearby galaxies with submillimeter telescopes on the ground has two major challenges. First, the brightness is significantly reduced at long submillimeter wavelengths compared to the brightness at the peak of the dust emission. Second, it is necessary to use a high-pass spatial filter to remove atmospheric noise on large angular scales, which has the unwelcome by-product of also removing the galaxy's large-scale structure. We have developed a technique for producing high-resolution submillimeter images of galaxies of large angular size by using the telescope on the ground to determine the small-scale structure (the large Fourier components) and a space telescope (Herschel or Planck) to determine the large-scale structure (the small Fourier components). Using this technique, we are carrying out the HARP and SCUBA-2 High Resolution Terahertz Andromeda Galaxy Survey (HASHTAG), an international Large Program on the James Clerk Maxwell Telescope, with one aim being to produce the first high-fidelity high-resolution submillimeter images of Andromeda. In this paper, we describe the survey, the method we have developed for combining the space-based and ground-based data, and present the first HASHTAG images of Andromeda at 450 and 850um. We also have created a method to predict the CO(J=3-2) line flux across M31, which contaminates the 850um band. We find that while normally the contamination is below our sensitivity limit, the contamination can be significant (up to 28%) in a few of the brightest regions of the 10 kpc ring. We therefore also provide images with the predicted line emission removed. △ Less

Submitted 30 September, 2021; originally announced October 2021.

Comments: 26 pages, 19 figures. Submitted to ApJS June 2021, Accepted September 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)

Resolved Neutral Outflow from a Lensed Dusty Star Forming Galaxy at z=2.09

Authors: Kirsty M. Butler, Paul P. van der Werf, Matus Rybak, Tiago Costa, Pierre Cox, Axel Weiß, Michał J. Michałowski, Dominik A. Riechers, Dimitra Rigopoulou, Lucia Marchetti, Stephen Eales, Ivan Valtchanov

Abstract: We report the detection of a massive neutral gas outflow in the z=2.09 gravitationally lensed Dusty Star-Forming Galaxy HATLASJ085358.9+015537 (G09v1.40), seen in absorption with the OH+(1_1-1_0) transition using spatially resolved (0.5"x0.4") Atacama Large Millimeter/submillimeter Array (ALMA) observations. The blueshifted OH+ line is observed simultaneously with the CO(9-8) emission line and und… ▽ More We report the detection of a massive neutral gas outflow in the z=2.09 gravitationally lensed Dusty Star-Forming Galaxy HATLASJ085358.9+015537 (G09v1.40), seen in absorption with the OH+(1_1-1_0) transition using spatially resolved (0.5"x0.4") Atacama Large Millimeter/submillimeter Array (ALMA) observations. The blueshifted OH+ line is observed simultaneously with the CO(9-8) emission line and underlying dust continuum. These data are complemented by high angular resolution (0.17"x0.13") ALMA observations of CH+(1-0) and underlying dust continuum, and Keck 2.2 micron imaging tracing the stellar emission. The neutral outflow, dust, dense molecular gas and stars all show spatial offsets from each other. The total atomic gas mass of the observed outflow is 6.7x10^9 M_sun, >25% as massive as the gas mass of the galaxy. We find that a conical outflow geometry best describes the OH+ kinematics and morphology and derive deprojected outflow properties as functions of possible inclination (0.38 deg-64 deg). The neutral gas mass outflow rate is between 83-25400 M_sun/yr, exceeding the star formation rate (788+/-300 M_sun/yr) if the inclination is >3.6 deg (mass-loading factor = 0.3-4.7). Kinetic energy and momentum fluxes span 4.4-290x10^9 L_sun and 0.1-3.7x10^37 dyne, respectively (energy-loading factor = 0.013-16), indicating that the feedback mechanisms required to drive the outflow depend on the inclination assumed. We derive a gas depletion time between 29 and 1 Myr, but find that the neutral outflow is likely to remain bound to the galaxy, unless the inclination is small, and may be re-accreted if additional feedback processes do not occur. △ Less

Submitted 2 July, 2021; v1 submitted 20 April, 2021; originally announced April 2021.

Comments: 33 pages, 20 figures

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)

A Complete 16 micron-Selected Galaxy Sample at $z\sim1$: Mid-infrared Spectral Energy Distributions

Authors: J. -S. Huang, Y. -S. Dai, S. P. Willner, S. M. Faber, C. Cheng, H. Xu, S. Wu, X. Shao, C. Hao, X. Xia, D. Rigopoulou, M. Pereira Santaella, G. Magdis, I. Cortzen, H. Yan, G. Fazio, P. Assmann, N. Araneda, L. Fan, M. Musin, Z. Wang, K. C. Xu, C. He, A. Esamdin

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

Submitted 28 April, 2021; v1 submitted 8 March, 2021; originally announced March 2021.

Comments: 14 pages, 5 tables, 29 figures, ApJ accepted v2 corrects author name formatting

Origins Space Telescope: From First Light to Life -- ESA Voyage 2050 White Paper

Authors: M. C. Wiedner, S. Aalto, J. Birkby, D. Burgarella, P. Caselli, V. Charmandaris, A. Cooray, E. De Beck, J. -M. Desert, M. Gerin, J. Goicoechea, M. Griffin, P. Hartogh, F. Helmich, M. Hogerheijde, L. Hunt, A. Karska, Q. Krall, D. Leisawitz, G. Melnick, M. Meixner, M. Mikako, Ch. Pearson, D. Rigopoulou, T. Roellig , et al. (2 additional authors not shown)

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

Submitted 4 December, 2020; originally announced December 2020.

Comments: 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 in Seyfert and star-forming galaxies

Authors: I. García-Bernete, D. Rigopoulou, A. Alonso-Herrero, M. Pereira-Santaella, P. F. Roche, B. Kerkeni

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

Submitted 26 October, 2021; v1 submitted 21 November, 2020; originally announced November 2020.

Comments: 29 pages, 21 figures, 6 tables, MNRAS, accepted 2021 October 22

The properties of Polycyclic Aromatic Hydrocarbons in galaxies: constraints on PAH sizes, charge and radiation fields

Authors: Dimitra Rigopoulou, Marie Barale, David Clary, Xiao Shan, Almudena Alonso-Herrero, Ismael García-Bernete, Leslie Hunt, Boutheina Kerkeni, Miguel Pereira-Santaella, Patrick Roche

Abstract: Based on theoretical spectra computed using Density Functional Theory we study the properties of Polycyclic Aromatic Hydrocarbons (PAH). In particular using bin-average spectra of PAH molecules with varying number of carbons we investigate how the intensity of the mid-infrared emission bands, 3.3, 6.2, 7.7 and 11.3 microns, respond to changes in the number of carbons, charge of the molecule, and t… ▽ More Based on theoretical spectra computed using Density Functional Theory we study the properties of Polycyclic Aromatic Hydrocarbons (PAH). In particular using bin-average spectra of PAH molecules with varying number of carbons we investigate how the intensity of the mid-infrared emission bands, 3.3, 6.2, 7.7 and 11.3 microns, respond to changes in the number of carbons, charge of the molecule, and the hardness of the radiation field that impinges the molecule. We confirm that the 6.2/7.7 band ratio is a good predictor for the size of the PAH molecule (based on the number of carbons present). We also investigate the efficacy of the 11.3/3.3 ratio to trace the size of PAH molecules and note the dependence of this ratio on the hardness of the radiation field. While the ratio can potentially also be used to trace PAH molecular size, a better understanding of the impact of the underlying radiation field on the 3.3 microns feature and the effect of the extinction on the ratio should be evaluated. The newly developed diagnostics are compared to band ratios measured in a variety of galaxies observed with the Infrared Spectrograph on board the Spitzer Space Telescope. We demonstrate that the band ratios can be used to probe the conditions of the interstellar medium in galaxies and differentiate between environments encountered in normal star forming galaxies and Active Galactic Nuclei. Our work highlights the immense potential that PAH observations with the James Webb Space Telescope will have on our understanding of the PAH emission itself and of the physical conditions in galaxies near and far. △ Less

Submitted 3 May, 2021; v1 submitted 19 November, 2020; originally announced November 2020.

Comments: 14 pages; 12 figures. Final version accepted for publication in MNRAS, accepted 2021 March 31

Rise of the Titans: Gas Excitation and Feedback in a Binary Hyper-Luminous Dusty Starburst Galaxy at z~6

Authors: Dominik A. Riechers, Hooshang Nayyeri, Denis Burgarella, Bjorn H. C. Emonts, David L. Clements, Asantha Cooray, Rob J. Ivison, Seb Oliver, Ismael Perez-Fournon, Dimitra Rigopoulou, Douglas Scott

Abstract: We report new observations toward the hyper-luminous dusty starbursting major merger ADFS-27 (z=5.655), using ATCA and ALMA. We detect CO 2-1, 8-7, 9-8, 10-9 and H2O(321-221) emission, and a P-Cygni-shaped OH+(11-01) absorption/emission feature. We also tentatively detect H2O(321-312) and OH+(12-01) emission and CH+(1-0) absorption. We find a total cold molecular mass of M_gas = (2.1+/-0.2) x 10^1… ▽ More We report new observations toward the hyper-luminous dusty starbursting major merger ADFS-27 (z=5.655), using ATCA and ALMA. We detect CO 2-1, 8-7, 9-8, 10-9 and H2O(321-221) emission, and a P-Cygni-shaped OH+(11-01) absorption/emission feature. We also tentatively detect H2O(321-312) and OH+(12-01) emission and CH+(1-0) absorption. We find a total cold molecular mass of M_gas = (2.1+/-0.2) x 10^11 (alpha_CO/1.0) Msun. We also find that the excitation of the star-forming gas is overall moderate for a z>5 dusty starburst, which is consistent with its moderate dust temperature. A high density, high kinetic temperature gas component embedded in the gas reservoir is required to fully explain the CO line ladder. This component is likely associated with the "maximum starburst" nuclei in the two merging galaxies, which are separated by only (140+/-13) km/s along the line of sight and 9.0 kpc in projection. The kinematic structure of both components is consistent with galaxy disks, but this interpretation remains limited by the spatial resolution of the current data. The OH+ features are only detected towards the northern component, which is also the one that is more enshrouded in dust and thus remains undetected up to 1.6 um even in our sensitive new HST/WFC3 imaging. The absorption component of the OH+ line is blueshifted and peaks near the CO and continuum emission peak while the emission is redshifted and peaks offset by 1.7 kpc from the CO and continuum emission peak, suggesting that the gas is associated with a massive molecular outflow from the intensely star-forming nucleus that supplies 125 Msun/yr of enriched gas to its halo. △ Less

Submitted 30 November, 2020; v1 submitted 28 October, 2020; originally announced October 2020.

Comments: 20 pages, 11 figures, 6 tables; to appear in the Astrophysical Journal (accepted November 29, 2020)

Multi-phase feedback processes in the Sy2 galaxy NGC 5643

Authors: I. García-Bernete, A. Alonso-Herrero, S. García-Burillo, M. Pereira-Santaella, B. García-Lorenzo, F. J. Carrera, D. Rigopoulou, C. Ramos Almeida, M. Villar Martín, O. González-Martín, E. K. S. Hicks, A. Labiano, C. Ricci, S. Mateos

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

Submitted 25 September, 2020; originally announced September 2020.

Comments: 26 pages, 21 figures, 5 tables. Astronomy and Astrophysics, Accepted 2020 September 11, in press

Report number: aa38256-20

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

Excitation and acceleration of molecular outflows in LIRGs: The extended ESO 320-G030 outflow on 200-pc scales

Authors: M. Pereira-Santaella, L. Colina, S. García-Burillo, E. González-Alfonso, A. Alonso-Herrero, S. Arribas, S. Cazzoli, J. Piqueras-López, D. Rigopoulou, A. Usero

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

Submitted 17 September, 2020; originally announced September 2020.

Comments: Accepted for publication in A&A

Journal ref: A&A 643, A89 (2020)

Simulating gas kinematic studies of high-redshift galaxies with the HARMONI Integral Field Spectrograph

Authors: Mark L. A. Richardson, Laurence Routledge, Niranjan Thatte, Matthias. Tecza, Ryan C. W. Houghton, Miguel Pereira-Santaella, Dimitra Rigopoulou

Abstract: We present simulated observations of gas kinematics in galaxies formed in 10 pc resolution cosmological simulations with the hydrodynamical + N-body code RAMSES, using the new RAMSES2HSIM pipeline with the simulated observing pipeline (HSIM) for the ELT HARMONI IFU spectrograph. We post-process the galaxy's gas kinematics and Hα line emission for each simulation cell, and integrate the emission to… ▽ More We present simulated observations of gas kinematics in galaxies formed in 10 pc resolution cosmological simulations with the hydrodynamical + N-body code RAMSES, using the new RAMSES2HSIM pipeline with the simulated observing pipeline (HSIM) for the ELT HARMONI IFU spectrograph. We post-process the galaxy's gas kinematics and Hα line emission for each simulation cell, and integrate the emission to produce an extinction-corrected input cube. We then simulate observations of the input cube with HARMONI, for a range of exposure times, spatial sampling, and spectral resolution. We analyze the mock observations to recover galaxy properties such as its kinematics and compare with the known simulation values. We investigate the cause of biases between the 'real' and 'observed' kinematic values, demonstrating the sensitivity of the inferred rotation curve to knowledge of the instrument's point spread function. △ Less

Submitted 10 September, 2020; originally announced September 2020.

Comments: 29 pages aastex format, 9 figures, 3 tables. Accepted for publication in MNRAS

Extinction in the 11.2 micron PAH band and the low L_11.2/L_IR in ULIRGs

Authors: Antonio Hernan-Caballero, Henrik W. W. Spoon, Almudena Alonso-Herrero, Evanthia Hatziminaoglou, Georgios E. Magdis, Pablo G. Perez-Gonzalez, Miguel Pereira-Santaella, Santiago Arribas, Isabella Cortzen, Alvaro Labiano, Javier Piqueras, Dimitra Rigopoulou

Abstract: We present a method for recovering the intrinsic (extinction-corrected) luminosity of the 11.2 micron PAH band in galaxy spectra. Using 105 high S/N Spitzer/IRS spectra of star-forming galaxies, we show that the equivalent width ratio of the 12.7 and 11.2 micron PAH bands is independent on the optical depth, with small dispersion of ~5% indicative of a nearly constant intrinsic flux ratio R_int =… ▽ More We present a method for recovering the intrinsic (extinction-corrected) luminosity of the 11.2 micron PAH band in galaxy spectra. Using 105 high S/N Spitzer/IRS spectra of star-forming galaxies, we show that the equivalent width ratio of the 12.7 and 11.2 micron PAH bands is independent on the optical depth, with small dispersion of ~5% indicative of a nearly constant intrinsic flux ratio R_int = (f_12.7/f_11.2)_int = 0.377 +/- 0.020. Conversely, the observed flux ratio, R_obs = (f_12.7/f_11.2)_obs strongly correlates with the silicate strength (S_sil) confirming that differences in R_obs reflect variation in the optical depth. The relation between R_obs and S_sil reproduces predictions for the Galactic Centre extinction law but disagrees with other laws. We calibrate the total extinction affecting the 11.2 micron PAH from R_obs, which we apply to another sample of 215 galaxies with accurate measurements of the total infrared luminosity (L_IR) to investigate the impact of extinction on L_11.2/L_IR. Correlation between L_11.2/L_IR and R_obs independently on L_IR suggests that increased extinction explains the well known decrease in the average L_11.2/L_IR at high L_IR. The extinction-corrected L_11.2 is proportional to L_IR in the range L_IR/L_sun = 10^9--10^13. These results consolidate L_11.2 as a robust tracer of star formation in galaxies. △ Less

Submitted 1 September, 2020; originally announced September 2020.

Comments: 14 pages, 9 figures, 2 tables. Accepted for publication in MNRAS

ALMA [N ıı] 205 μm Imaging Spectroscopy of the Lensed Submillimeter galaxy ID 141 at redshift 4.24

Authors: Cheng Cheng, Xiaoyue Cao, Nanyao Lu, Ran Li, Chentao Yang, Dimitra Rigopoulou, Vassilis Charmandaris, Yu Gao, Cong Kevin Xu, Paul van der Werf, Tanio Diaz Santos, George C. Privon, Yinghe Zhao, Tianwen Cao, Y. Sophia Dai, Jia-Sheng Huang, David Sanders, Chunxiang Wang, Zhong Wang, Lei Zhu

Abstract: We present the Atacama Large Millimeter/submillimeter Array (ALMA) observation of the Sub-millimeter galaxy (SMG) ID 141 at z=4.24 in the [N II] 205 $μ$m line (hereafter [N II]) and the underlying continuum at (rest-frame) 197.6 $μ$m. Benefiting from lensing magnification by a galaxy pair at z=0.595, ID 141 is one of the brightest z$>4$ SMGs. At the angular resolutions of $\sim1.2"$ to $1.5"$ (… ▽ More We present the Atacama Large Millimeter/submillimeter Array (ALMA) observation of the Sub-millimeter galaxy (SMG) ID 141 at z=4.24 in the [N II] 205 $μ$m line (hereafter [N II]) and the underlying continuum at (rest-frame) 197.6 $μ$m. Benefiting from lensing magnification by a galaxy pair at z=0.595, ID 141 is one of the brightest z$>4$ SMGs. At the angular resolutions of $\sim1.2"$ to $1.5"$ ($1" \sim6.9$ kpc), our observation clearly separates, and moderately resolves the two lensed images in both continuum and line emission at $\rm S/N>5$ . Our continuum-based lensing model implies an averaged amplification factor of $\sim5.8$ and reveals that the de-lensed continuum image has the Sérsic index $\simeq 0.95$ and the Sérsic radius of $\sim0.18" (\sim 1.24$ kpc). Furthermore, the reconstructed [N II] velocity field in the source plane is dominated by a rotation component with a maximum velocity of $\sim 300$ km/s at large radii, indicating a dark matter halo mass of $\sim 10^{12}M_{\odot}$. This, together with the reconstructed velocity dispersion field being smooth and modest in value ($<100$ km/s) over much of the outer parts of the galaxy, favours the interpretation of ID 141 being a disk galaxy dynamically supported by rotation. The observed [N II]/CO (7-6) and [N II]/[C II] 158 $μ$m line luminosity ratios, which are consistent with the corresponding line ratio vs. far-infrared color correlation from local luminous infrared galaxies, imply a de-lensed star formation rate of ($1.8\pm 0.6)\times10^3M_\odot$/yr and provide an independent estimate on the size of the star-forming region $0.7^{+0.3}_{-0.3}$ kpc in radius. △ Less

Submitted 16 October, 2020; v1 submitted 1 June, 2020; originally announced June 2020.

Comments: 13 pages, 6 figures, 2 tables. Lensing model code can be found here https://gitlab.com/cxylzlx/tiny_lens

Cold molecular gas and PAH emission in Seyfert galaxies

Authors: A. Alonso-Herrero, M. Pereira-Santaella, D. Rigopoulou, I. Garcia-Bernete, S. Garcia-Burillo, A. J. Dominguez-Fernandez, F. Combes, R. I. Davies, T. Diaz-Santos, D. Esparza-Arredondo, O. Gonzalez-Martin, A. Hernan-Caballero, E. K. S. Hicks, S. F. Hoenig, N. A. Levenson, C. Ramos Almeida, P. F. Roche, D. Rosario

Abstract: We investigate the relation between the detection of the $11.3\,μ$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… ▽ More We investigate the relation between the detection of the $11.3\,μ$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\,μ$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\,μ$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. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: Accepted for publication in A&A

Dust and gas content of high-redshift galaxies hosting obscured AGN in the CDF-S

Authors: Q. D'Amato, R. Gilli, C. Vignali, M. Massardi, F. Pozzi, G. Zamorani, C. Circosta, F. Vito, J. Fritz, G. Cresci, V. Casasola, F. Calura, A. Feltre, V. Manieri, D. Rigopoulou, P. Tozzi, C. Norman

Abstract: 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… ▽ More 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 α_CO=0.8 and up to ~2 x 10^11~M_sun for α_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. △ Less

Submitted 19 March, 2020; originally announced March 2020.

Comments: 15 pages, 6 figures, 7 tables. Accepted for publication in A&A

Journal ref: A&A 636, A37 (2020)