Showing 1–50 of 105 results for author: Gonzalez-Alfonso, E

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

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)

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

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

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

Submitted 25 March, 2024; originally announced March 2024.

Comments: 9 pages, 10 figures

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

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

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

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

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

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

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

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

JWST detection of extremely excited outflowing CO and H2O in VV 114 E SW: a possible rapidly accreting IMBH

Authors: Eduardo González-Alfonso, Ismael García-Bernete, Miguel Pereira-Santaella, David A. Neufeld, Jacqueline Fischer, Fergus R. Donnan

Abstract: Mid-infrared (mid-IR) gas-phase molecular bands are powerful diagnostics of the warm interstellar medium. We report the James Webb Space Telescope detection of the CO v=1-0 (4.4-5.0 um) and H2O nu2=1-0 (5.0-7.8um) ro-vibrational bands, both in absorption, toward the ``s2'' core in the southwest nucleus of the merging galaxy VV 114 E. All ro-vibrational CO lines up to J_low=33 (E_low~3000 K) are de… ▽ More Mid-infrared (mid-IR) gas-phase molecular bands are powerful diagnostics of the warm interstellar medium. We report the James Webb Space Telescope detection of the CO v=1-0 (4.4-5.0 um) and H2O nu2=1-0 (5.0-7.8um) ro-vibrational bands, both in absorption, toward the ``s2'' core in the southwest nucleus of the merging galaxy VV 114 E. All ro-vibrational CO lines up to J_low=33 (E_low~3000 K) are detected, as well as a forest of H2O lines up to 13_{0,13} (E_low~2600 K). The highest-excitation lines are blueshifted by ~180 km s^{-1} relative to the extended molecular cloud, which is traced by the rotational CO J=3-2 346 GHz line observed with the Atacama Large Millimeter/submillimeter Array. The bands also show absorption in a low-velocity component (blueshifted by ~30 km s^{-1}) with lower excitation. The analysis shows that the bands are observed against a continuum with effective temperature of T_bck~550 K extinguished with tau_6um^ext~ 2.5-3 (A_k~6.9-8.3 mag). The high-excitation CO and H2O lines are consistent with v=0 thermalization with T_rot~450 K and column densities of N_CO~(1.7-3.5)x10^{19} cm^{-2} and N_H2O~(1.5-3.0)x10^{19} cm$^{-2}$. Thermalization of the v=0 levels of H2O requires either an extreme density of n_H2>~10^9 cm^{-3}, or radiative excitation by the mid-IR field in a very compact (<1 pc) optically thick source emitting ~10^{10} L_sun. The latter alternative is favored, implying that the observed absorption probes the very early stages of a fully enshrouded active black hole (BH). On the basis of a simple model for BH growth and applying a lifetime constraint to the s2 core, an intermediate-mass BH (IMBH, M_BH~4.5x10^4 M_sun) accreting at super-Eddington rates is suggested, where the observed feedback has not yet been able to break through the natal cocoon. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

The first ground-based detection of the 752 GHz water line in local ultra-luminous infrared galaxies using APEX-SEPIA

Authors: Daysi Quinatoa, Chentao Yang, Edo Ibar, Elizabeth Humphreys, Susanne Aalto, Loreto Barcos-Muñoz, Eduardo González-Alfonso, Violette Impellizzeri, Yara Jaffé, Lijie Liu, Sergio Martín, Axel Weiss, Zhi-Yu Zhang

Abstract: We report the first ground-based detection of the water line p-H2O (211-202) at 752.033 GHz in three z < 0.08 ultra-luminous infrared galaxies (ULIRGs): IRAS 06035-7102, IRAS 17207-0014 and IRAS 09022-3615. Using the Atacama Pathfinder EXperiment (APEX), with its Swedish-ESO PI Instrument for APEX (SEPIA) band-9 receiver, we detect this H2O line with overall signal-to-noise ratios of 8-10 in all t… ▽ More We report the first ground-based detection of the water line p-H2O (211-202) at 752.033 GHz in three z < 0.08 ultra-luminous infrared galaxies (ULIRGs): IRAS 06035-7102, IRAS 17207-0014 and IRAS 09022-3615. Using the Atacama Pathfinder EXperiment (APEX), with its Swedish-ESO PI Instrument for APEX (SEPIA) band-9 receiver, we detect this H2O line with overall signal-to-noise ratios of 8-10 in all three galaxies. Notably, this is the first detection of this line in IRAS 06035-7102. Our new APEX-measured fluxes, between 145 to 705 Jy km s-1, are compared with previous values taken from Herschel SPIRE FTS. We highlight the great capabilities of APEX for resolving the H2O line profiles with high spectral resolutions while also improving by a factor of two the significance of the detection within moderate integration times. While exploring the correlation between the p-H2O(211-202) and the total infrared luminosity, our galaxies are found to follow the trend at the bright end of the local ULIRG's distribution. The p-H2O(211-202) line spectra are compared to the mid-J CO and HCN spectra, and dust continuum previously observed with ALMA. In the complex interacting system IRAS 09022-3615, the profile of the water emission line is offset in velocity with respect to the ALMA CO(J = 4 - 3) emission. For IRAS 17207-0014 and IRAS 06035-7102, the profiles between the water line and the CO lines are spectroscopically aligned. This pilot study demonstrates the feasibility of directly conducting ground-based high-frequency observations of this key water line, opening the possibility of detailed follow-up campaigns to tackle its nature. △ Less

Submitted 3 November, 2023; originally announced November 2023.

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

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

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

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

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

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

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

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

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

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

SUNRISE: The rich molecular inventory of high-redshift dusty galaxies revealed by broadband spectral line surveys

Authors: Chentao Yang, Alain Omont, Sergio Martín, Thomas G. Bisbas, Pierre Cox, Alexandre Beelen, Eduardo González-Alfonso, Raphaël Gavazzi, Susanne Aalto, Paola Andreani, Cecilia Ceccarelli, Yu Gao, Mark Gorski, Michel Guélin, Hai Fu, R. J. Ivison, Kirsten K. Knudsen, Matthew Lehnert, Hugo Messias, Sebastien Muller, Roberto Neri, Dominik Riechers, Paul van der Werf, Zhi-Yu Zhang

Abstract: Understanding the nature of high-$z$ dusty galaxies requires a comprehensive view of their ISM and molecular complexity. However, the molecular ISM at high-$z$ is commonly studied using only a few species beyond CO, limiting our understanding. In this paper, we present the results of deep 3 mm spectral line surveys using the NOEMA targeting two lensed dusty galaxies: APM 08279+5255 (APM), a quasar… ▽ More Understanding the nature of high-$z$ dusty galaxies requires a comprehensive view of their ISM and molecular complexity. However, the molecular ISM at high-$z$ is commonly studied using only a few species beyond CO, limiting our understanding. In this paper, we present the results of deep 3 mm spectral line surveys using the NOEMA targeting two lensed dusty galaxies: APM 08279+5255 (APM), a quasar at redshift $z=3.911$, and NCv1.143 (NC), a $z=3.565$ starburst galaxy. The spectral line surveys cover rest-frame frequencies from about 330-550 GHz. We report the detection of 38 and 25 emission lines in APM and NC, respectively. The spectra reveal the chemical richness and the complexity of the physical properties of the ISM. By comparing the spectra of the two sources and combining the gas excitation analysis, we find that the physical properties and the chemical imprints of the ISM are different between them: the molecular gas is more excited in APM, exhibiting higher molecular-gas temperatures and densities compared to NC; the chemical abundances in APM are akin to the values of local AGN, showing boosted relative abundances of the dense gas tracers that might be related to high-temperature chemistry and/or XDRs, while NC more closely resembles local starburst galaxies. The most significant differences are found in H2O, where the 448GHz H2O line is significantly brighter in APM, likely linked to the intense far-infrared radiation from the dust powered by AGN. Our astrochemical model suggests that at such high column densities, FUV radiation is less important in regulating the ISM, while CRs (X-rays/shocks) are the key players in shaping the abundance of the molecules and the initial conditions of star formation. Such deep spectral line surveys open a new window to study the physical and chemical properties of the ISM and the radiation field of galaxies in the early Universe. (abridged) △ Less

Submitted 22 October, 2023; v1 submitted 14 August, 2023; originally announced August 2023.

Comments: Accepted for publication in A&A. 39 pages (including 9 pages of Appendix). Updated after language editing

Journal ref: A&A 680, A95 (2023)

The importance of radiative pumping on the emission of the H_2O submillimeter lines in galaxies

Authors: Eduardo González-Alfonso, Jacqueline Fischer, Javier R. Goicoechea, Chentao Yang, Miguel Pereira-Santaella, Kenneth P. Stewart

Abstract: H_2O submillimeter emission is a powerful diagnostic of the molecular interstellar medium in a variety of sources, including low- and high-mass star forming regions of the Milky Way, and from local to high redshift galaxies. However, the excitation mechanism of these lines in galaxies has been debated, preventing a basic consensus on the physical information that H_2O provides. Both radiative pump… ▽ More H_2O submillimeter emission is a powerful diagnostic of the molecular interstellar medium in a variety of sources, including low- and high-mass star forming regions of the Milky Way, and from local to high redshift galaxies. However, the excitation mechanism of these lines in galaxies has been debated, preventing a basic consensus on the physical information that H_2O provides. Both radiative pumping due to H_2O absorption of far-infrared photons emitted by dust and collisional excitation in dense shocked gas have been proposed to explain the H_2O emission. Here we propose two basic diagnostics to distinguish between the two mechanisms: 1) in shock excited regions, the ortho-H_2O 3_{21}-2_{12} 75um and the para-H_2O 2_{20}-1_{11} 101um rotational lines are expected to be in emission while, if radiative pumping dominates, both far-infrared lines are expected to be in absorption; 2) based on statistical equilibrium of H_2O level populations, the radiative pumping scenario predicts that the apparent isotropic net rate of far-infrared absorption in the 3_{21}-2_{12} (75um) and 2_{20}-1_{11} (101um) lines should be higher than or equal to the apparent isotropic net rate of submillimeter emission in the 3_{21}-3_{12} (1163 GHz) and 2_20-2_{11} (1229 GHz) lines, respectively. Applying both criteria to all 16 galaxies and several galactic high-mass star-forming regions where the H_2O 75um and submillimeter lines have been observed with Herschel/PACS and SPIRE, we show that in most (extra)galactic sources the H_2O submillimeter line excitation is dominated by far-infrared pumping, with collisional excitation of the low-excitation levels in some of them. Based on this finding, we revisit the interpretation of the correlation between the luminosity of the H_2O 988 GHz line and the source luminosity in the combined galactic and extragalactic sample. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: Accepted for publication in Astronomy & Astrophysics Letters

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

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

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

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

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

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

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

On the thermal structure of the proto-Super Star Cluster 13 in NGC 253

Authors: F. Rico-Villas, E. González-Alfonso, J. Martín-Pintado, V. M. Rivilla, S. Martín

Abstract: Using high angular resolution ALMA observations ($0.02^{\prime\prime}\approx0.34$ pc), we study the thermal structure and kinematics of the proto super star cluster $13$ in the central region of NGC253 through their continuum and vibrationally excited HC$_3$N emission from $J=24-23$ and $J=26-25$ lines arising from vibrational states up to $v_4=1$. We have carried 2D-LTE and non-local radiative tr… ▽ More Using high angular resolution ALMA observations ($0.02^{\prime\prime}\approx0.34$ pc), we study the thermal structure and kinematics of the proto super star cluster $13$ in the central region of NGC253 through their continuum and vibrationally excited HC$_3$N emission from $J=24-23$ and $J=26-25$ lines arising from vibrational states up to $v_4=1$. We have carried 2D-LTE and non-local radiative transfer modelling of the radial profile of the HC$_3$N and continuum emission in concentric rings of $0.1$ pc width. From the 2D-LTE analysis, we found a Super Hot Core (SHC) of $1.5$ pc with very high vibrational temperatures ($>500$ K), and a jump in the radial velocity ($21$ km s$^{-1}$) in the SE-NW direction. From the non-local models, we derive the HC$_3$N column density, H$_2$ density and dust temperature ($T_\text{dust}$) profiles. Our results show that the thermal structure of the SHC is dominated by the greenhouse effect due to the high dust opacity in the IR, leading to an overestimation of the LTE $T_\text{dust}$ and its derived luminosity. The kinematics and $T_\text{dust}$ profile of the SHC suggest that star formation was likely triggered by a cloud-cloud collision. We compare proto-SSC $13$ to other deeply embedded star-forming regions, and discuss the origin of the $L_\text{IR}/M_{\text{H}_2}$ excess above $\sim100$ L$_\odot$ M$_\odot^{-1}$ observed in (U)LIRGs. △ Less

Submitted 3 August, 2022; originally announced August 2022.

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

The ramp-up of interstellar medium enrichment at z>4

Authors: M. Franco, K. E. K. Coppin, J. E. Geach, C. Kobayashi, S. C. Chapman, C. Yang, E. González-Alfonso, J. S. Spilker, A. Cooray, M. J. Michałowski

Abstract: Fluorine is one of the most interesting elements for nuclear and stellar astrophysics. Fluorine abundance was first measured for stars other than the Sun in 1992, then for a handful metal-poor stars, which are likely to have formed in the early Universe. The main production sites of fluorine are under debate and include asymptotic giant branch (AGB) stars, $ν$-process in core-collapse supernovae,… ▽ More Fluorine is one of the most interesting elements for nuclear and stellar astrophysics. Fluorine abundance was first measured for stars other than the Sun in 1992, then for a handful metal-poor stars, which are likely to have formed in the early Universe. The main production sites of fluorine are under debate and include asymptotic giant branch (AGB) stars, $ν$-process in core-collapse supernovae, and Wolf-Rayet (WR) stars. Due to the difference in the mass and lifetime of progenitor stars, high redshift observations of fluorine can help constrain the mechanism of fluorine production in massive galaxies. Here, we report the detection of HF (S/N = 8) in absorption in a gravitationally lensed dusty star-forming galaxy at redshift z=4.4 with $N_{\rm HF}$/$N_{\rm{H_2}}$ as high as $\sim2\times10^{-9}$, indicating a very quick ramp-up of the chemical enrichment in this high-z galaxy. At z=4.4, AGB stars of a few solar masses are very unlikely to dominate the enrichment. Instead, we show that WR stars are required to produce the observed fluorine abundance at this time, with other production mechanisms becoming important at later times. These observations therefore provide an insight into the underlying processes driving the `ramp-up' phase of chemical enrichment alongside rapid stellar mass assembly in a young massive galaxy. △ Less

Submitted 4 November, 2021; originally announced November 2021.

Comments: 21 pages, 4 figures, published in Nature Astronomy

Deeply Buried Nuclei in the Infrared-Luminous Galaxies NGC 4418 and Arp 220: I. ALMA Observations at $λ= $1.4-0.4 mm and Continuum Analysis

Authors: Kazushi Sakamoto, Eduardo Gonzalez-Alfonso, Sergio Martin, David J. Wilner, Susanne Aalto, Aaron S. Evans, Nanase Harada

Abstract: We observed with ALMA three deeply buried nuclei in two galaxies, NGC 4418 and Arp 220, at $\sim$0.2$''$ resolution over a total bandwidth of 67 GHz in $f_{\rm rest}$ = 215 - 697 GHz. Here we (1) introduce our program, (2) describe our data reduction method for wide-band, high-resolution imaging spectroscopy, (3) analyze in visibilities the compact nuclei with line forests, (4) develop a continuum… ▽ More We observed with ALMA three deeply buried nuclei in two galaxies, NGC 4418 and Arp 220, at $\sim$0.2$''$ resolution over a total bandwidth of 67 GHz in $f_{\rm rest}$ = 215 - 697 GHz. Here we (1) introduce our program, (2) describe our data reduction method for wide-band, high-resolution imaging spectroscopy, (3) analyze in visibilities the compact nuclei with line forests, (4) develop a continuum-based estimation method of dust opacity and gas column density in heavily obscured nuclei, which uses the BGN (buried galactic nuclei) model and is sensitive to $\log(N_{\rm H_2}/{\rm cm}^{-2}) \sim $ 25 - 26 at $λ\sim 1$ mm, and (5) present the continuum data and diagnosis of our targets. The three continuum nuclei have major-axis FWHM of $\sim$0.1$''$-0.3$''$ (20-140 pc) aligned to their rotating nuclear disks of molecular gas. However, each nucleus is described better with two or three concentric components than with a single Gaussian. The innermost cores have sizes of 0.05$''$-0.10$''$ (8-40 pc), peak brightness temperatures of ~100-500 K at 350 GHz, and more fractional flux at lower frequencies. The intermediate components correspond to the nuclear disks. They have axial ratios of $\approx$0.5 and hence inclinations $\stackrel{>}{\sim} 60$ deg. The outermost elements include the bipolar outflow from Arp 220W. We estimate 1 mm dust opacity of $τ_{\rm d,1mm} \approx 2.2$, $1.2$, and $\stackrel{<}{\sim} 0.4$ respectively for NGC 4418, Arp 220W, and Arp 220E. The first two correspond to $\log(N_{\rm H}/{\rm cm}^{-2}) \sim 26$ for conventional dust-opacity laws, and hence the nuclei are highly Compton thick. △ Less

Submitted 14 September, 2021; originally announced September 2021.

Comments: 59 pages, 18 figures, 15 tables, accepted for publication in ApJ

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

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

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

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

Comments: Accepted for publication in A&A

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

CON-quest: Searching for the most obscured galaxy nuclei

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

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

Submitted 26 February, 2021; originally announced February 2021.

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

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

A proto-pseudobulge in ESO 320-G030 fed by a massive molecular inflow driven by a nuclear bar

Authors: Eduardo González-Alfonso, Miguel Pereira-Santaella, Jacqueline Fischer, Santiago García-Burillo, Chentao Yang, Almudena Alonso-Herrero, Luis Colina, Matthew L. N. Ashby, Howard A. Smith, Fernando Rico-Villas, Jesús Martín-Pintado, Sara Cazzoli, Kenneth P. Stewart

Abstract: Galaxies with nuclear bars are believed to efficiently drive gas inward, generating a nuclear starburst and possibly an active galactic nucleus (AGN). We confirm this scenario for the isolated, double-barred, luminous infrared galaxy ESO 320-G030 based on an analysis of Herschel and ALMA spectroscopic observations. Herschel/PACS and SPIRE observations of ESO 320-G030 show absorption/emission in 18… ▽ More Galaxies with nuclear bars are believed to efficiently drive gas inward, generating a nuclear starburst and possibly an active galactic nucleus (AGN). We confirm this scenario for the isolated, double-barred, luminous infrared galaxy ESO 320-G030 based on an analysis of Herschel and ALMA spectroscopic observations. Herschel/PACS and SPIRE observations of ESO 320-G030 show absorption/emission in 18 lines of H2O, which we combine with the ALMA H2O 423-330 448 GHz line (Eupper~400 K) and continuum images to study the nuclear region. Radiative transfer models indicate that 3 nuclear components are required to account for the H2O and continuum data. An envelope, with R~130-150 pc, T_dust~50 K, and N_H2~2x10^{23} cm^{-2}, surrounds a nuclear disk with R~40 pc and tau_100um~1.5-3 (N_H2~2x10^{24} cm^{-2}) and an extremely compact (R~12 pc), warm (~100 K), and buried (tau_100um>5, N_H2>~5x10^{24} cm^{-2}) core component. The three nuclear components account for 70% of the galaxy L_IR (SFR~16-18 Msun yr^{-1}). The nucleus is fed by a molecular inflow observed in CO 2-1 with ALMA, which is associated with the nuclear bar. With decreasing radius (r=450-225 pc), the mass inflow rate increases up to ~20 Msun yr^{-1}, which is similar to the nuclear SFR, indicating that the starburst is sustained by the inflow. At lower r, the inflow is best probed by the far-infrared OH ground-state doublets, with an estimated inflow rate of ~30 Msun yr^{-1}. The short timescale of ~20 Myr for nuclear gas replenishment indicates quick secular evolution, and indicates that we are witnessing an intermediate stage (<100 Myr) proto-pseudobulge fed by a massive inflow that is driven by a strong nuclear bar. We also apply the H2O model to the Herschel far-infrared spectroscopic observations of H2^{18}O, OH, $^{18}OH, OH+, H2O^+, H3O^+, NH, NH2, NH3, CH, CH^+, ^{13}CH^+, HF, SH, and C3, and estimate their abundances. △ Less

Submitted 4 November, 2020; v1 submitted 31 October, 2020; originally announced November 2020.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 645, A49 (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)

Vibrationally excited HC3N emission in NGC1068: Tracing the recent star formation in the starburst ring

Authors: Fernando Rico-Villas, Jesús Martín-Pintado, Eduardo González-Alfonso, Victor M. Rivilla, Sergio Martín, Santiago García-Burillo, Izaskun Jiménez-Serra, María Sánchez-García

Abstract: Using ALMA data, we have studied the HC$_3$N and continuum emission in the starburst ring (SB ring) and the circumnuclear disc (CND) of the SB/AGN composite galaxy NGC 1068. We have detected emission from vibrationally excited HC$_3$N (HC$_3$N*) only towards one star-forming region of the SB ring. Remarkably, HC$_3$N* was not detected towards the CND despite its large HC$_3$N $v=0$ column density.… ▽ More Using ALMA data, we have studied the HC$_3$N and continuum emission in the starburst ring (SB ring) and the circumnuclear disc (CND) of the SB/AGN composite galaxy NGC 1068. We have detected emission from vibrationally excited HC$_3$N (HC$_3$N*) only towards one star-forming region of the SB ring. Remarkably, HC$_3$N* was not detected towards the CND despite its large HC$_3$N $v=0$ column density. From LTE and non-LTE modelling of HC$_3$N*, we obtained a dust temperature of $ T_\text{dust} \sim 250$ K and a density of $n_{\text{H}_2}=6\times10^5$ cm$^{-3}$ for this star-forming region. The estimated IR luminosity of $5.8\times10^8$ L$_\odot$ is typical of proto-Super Star Clusters (proto-SSC) observed in the SB galaxy NGC 253. We use the continuum emissions at $147$ GHz and $350$ GHz, along with CO and Pa $α$, to estimate the ages of other $14$ SSCs in the SB ring. We find the youngest SSCs to be associated with the region connecting the nuclear bar with the SB ring, supporting the scenario of sequential star formation. For the CND, our analysis yields $T_\text{dust} \leqslant 100$ K and $n_{\text{H}_2}\sim(3-6)\times10^5$cm$^{-3}$. The very different dust temperatures found for the CND and the proto-SSC indicates that, while the dust in the proto-SSC is being efficiently heated from the inside by the radiation from massive proto-stars, the CND is being heated externally by the AGN, which in the IR optically thin case can only heat the dust to $56$ K. We discuss the implications of the non-detection of HC$_3$N* near the luminous AGN in NGC 1068 on the interpretation of the HC$_3$N* emission observed in the SB/AGN composite galaxies NGC 4418 and Arp 220. △ Less

Submitted 26 January, 2021; v1 submitted 9 August, 2020; originally announced August 2020.

Comments: submitted to MNRAS on 08/08/2020

A Hard X-ray Test of HCN Enhancements as a Tracer of Embedded Black Hole Growth

Authors: G. C. Privon, C. Ricci, S. Aalto, S. Viti, L. Armus, T. Díaz-Santos, E. González-Alfonso, K. Iwasawa, D. L. Jeff, E. Treister, F. Bauer, A. S. Evans, P. Garg, R. Herrero-Illana, J. M. Mazzarella, K. Larson, L. Blecha, L. Barcos-Muñoz, V. Charmandaris, S. Stierwalt, M. A. Pérez-Torres

Abstract: Enhanced emission from the dense gas tracer HCN (relative to HCO$^+ $) has been proposed as a signature of active galactic nuclei (AGN). In a previous single-dish millimeter line survey we identified galaxies with HCN/HCO$ ^+ $ (1-0) intensity ratios consistent with those of many AGN but whose mid-infrared spectral diagnostics are consistent with little to no ( $\lesssim15\% $) contribution of an… ▽ More Enhanced emission from the dense gas tracer HCN (relative to HCO$^+ $) has been proposed as a signature of active galactic nuclei (AGN). In a previous single-dish millimeter line survey we identified galaxies with HCN/HCO$ ^+ $ (1-0) intensity ratios consistent with those of many AGN but whose mid-infrared spectral diagnostics are consistent with little to no ( $\lesssim15\% $) contribution of an AGN to the bolometric luminosity. To search for putative heavily obscured AGN, we present and analyze \nustar hard X-ray (3-79 keV) observations of four such galaxies from the Great Observatories All-sky LIRG Survey. We find no X-ray evidence for AGN in three of the systems and place strong upper limits on the energetic contribution of any heavily obscured ($N_{\rm H}>10^{24}$ cm$^{-2}$) AGN to their bolometric luminosity. The X-ray flux upper limits are presently an order of magnitude below what XDR-driven chemistry model predict are necessary to drive HCN enhancements. In a fourth system we find a hard X-ray excess consistent with the presence of an AGN, but contributing only $\sim3\%$ of the bolometric luminosity. It is also unclear if the AGN is spatially associated with the HCN enhancement. We further explore the relationship between HCN/HCO$^+$ (for several $\mathrm{J}_\mathrm{upper}$ levels) and $L_\mathrm{AGN}/L_\mathrm{IR}$ for a larger sample of systems in the literature. We find no evidence for correlations between the line ratios and the AGN fraction derived from X-rays, indicating that HCN/HCO$^+$ intensity ratios are not driven by the energetic dominance of AGN, nor are they reliable indicators of whether SMBH accretion is ongoing. △ Less

Submitted 24 April, 2020; v1 submitted 6 April, 2020; originally announced April 2020.

Comments: 14 pages, 4 figures. accepted for publication in The Astrophysical Journal

Journal ref: ApJ 893 149 (2020)

Distribution of Water Vapor in Molecular Clouds. II

Authors: Gary J. Melnick, Volker Tolls, Ronald L. Snell, Michael J. Kaufman, Edwin A. Bergin, Javier R. Goicoechea, Paul F. Goldsmith, Eduardo González-Alfonso, David J. Hollenbach, Dariusz C. Lis, David A. Neufeld

Abstract: The depth-dependent abundance of both gas-phase and solid-state water within dense, quiescent, molecular clouds is important to both the cloud chemistry and gas cooling. Where water is in the gas phase, it's free to participate in the network of ion-neutral reactions that lead to a host of oxygen-bearing molecules, and its many ortho and para energy levels make it an effective coolant for gas temp… ▽ More The depth-dependent abundance of both gas-phase and solid-state water within dense, quiescent, molecular clouds is important to both the cloud chemistry and gas cooling. Where water is in the gas phase, it's free to participate in the network of ion-neutral reactions that lead to a host of oxygen-bearing molecules, and its many ortho and para energy levels make it an effective coolant for gas temperatures greater than 20K. Where water is abundant as ice on grain surfaces, and unavailable to cool the gas, significant amounts of oxygen are removed from the gas phase, suppressing the gas-phase chemical reactions that lead to a number of oxygen-bearing species, including O2. Models of FUV-illuminated clouds predict that the gas-phase water abundance peaks in the range Av ~3 and 8mag of the cloud surface, depending on the gas density and FUV field strength. Deeper within such clouds, water is predicted to exist mainly as ice on grain surfaces. More broadly, these models are used to analyze a variety of other regions, including outflow cavities associated with young stellar objects and the surface layers of protoplanetary disks. In this paper, we report the results of observational tests of FUV-illuminated cloud models toward the Orion Molecular Ridge and Cepheus B using data obtained from the Herschel Space Observatory and the Five College Radio Astronomy Observatory. Toward Orion, 2220 spatial positions were observed along the face-on Orion Ridge in the H2O 110-101 557GHz and NH3 J,K=1,0-0,0 572GHz lines. Toward Cepheus B, two strip scans were made in the same lines across the edge-on ionization front. These new observations demonstrate that gas-phase water exists primarily within a few magnitudes of dense cloud surfaces, strengthening the conclusions of an earlier study based on a much smaller data set, and indirectly supports the prediction that water ice is quite abundant in dense clouds. △ Less

Submitted 10 March, 2020; originally announced March 2020.

Comments: 59 pages, including 23 figures, Accepted for publication in ApJ

First detection of the 448 GHz ortho-H2O line at high redshift: probing the structure of a starburst nucleus at z = 3.63

Authors: C. Yang, E. González-Alfonso, A. Omont, M. Pereira-Santaella, J. Fischer, A. Beelen, R. Gavazzi

Abstract: Submillimeter rotational lines of H2O are a powerful probe in warm gas regions of the ISM, tracing scales and structures ranging from kpc disks to the most compact and dust-obscured regions of galactic nuclei. The ortho-H2O(423-330) line at 448 GHz, which was recently detected in a local luminous infrared galaxy (Pereira-Santaella et al. 2017), offers a unique constraint on the excitation conditio… ▽ More Submillimeter rotational lines of H2O are a powerful probe in warm gas regions of the ISM, tracing scales and structures ranging from kpc disks to the most compact and dust-obscured regions of galactic nuclei. The ortho-H2O(423-330) line at 448 GHz, which was recently detected in a local luminous infrared galaxy (Pereira-Santaella et al. 2017), offers a unique constraint on the excitation conditions and ISM properties in deeply buried galaxy nuclei since the line requires high far-IR optical depths to be excited. In this letter, we report the first high-redshift detection of the 448 GHz H2O(423-330) line using ALMA, in a strongly lensed submillimeter galaxy (SMG) at z=3.63. After correcting for magnification, the luminosity of the 448 GHz H2O line is ~10^6 L_sun. In combination with three other previously detected H2O lines, we build a model that "resolves" the dusty ISM structure of the SMG, and find that it is composed of a ~1 kpc optically thin (optical depth at 100μm τ_{100}~0.3) disk component with dust temperature T_{dust} \approx 50 K emitting a total infrared power of 5e12 L_sun with surface density Σ_{IR}=4e11 L_sun kpc^{-2}, and a very compact (0.1 kpc) heavily dust-obscured (τ_{100} \gtrsim 1) nuclear core with very warm dust (100 K) and Σ_{IR}=8e12 L_sun kpc^{-2}. The H2O abundance in the core component, X_{H2O}~(0.3-5)e{-5}, is at least one order of magnitude higher than in the disk component. The optically thick core has the characteristic properties of an Eddington-limited starburst, providing evidence that radiation pressure on dust is capable of supporting the ISM in buried nuclei at high redshifts. The multi-component ISM structure revealed by our models illustrates that dust and molecules such as H2O are present in regions characterized by highly differing conditions and scales, extending from the nucleus to more extended regions of SMGs. △ Less

Submitted 17 January, 2020; originally announced January 2020.

Comments: 1 table and 4 figures, A&A Letter accepted

Journal ref: A&A 634, L3 (2020)

Molecular Gas Inflows and Outflows in Ultraluminous Infrared Galaxies at $z\sim0.2$ and one QSO at $z=6.1$

Authors: R. Herrera-Camus, E. Sturm, J. Graciá-Carpio, S. Veilleux, T. Shimizu, D. Lutz, M. Stone, E. González-Alfonso, R. Davies, J. Fischer, R. Genzel, R. Maiolino, A. Sternberg, L. Tacconi, A. Verma

Abstract: Aims. We aim to search and characterize inflows and outflows of molecular gas in four ultraluminous infrared galaxies (ULIRGs) at $z\sim0.2-0.3$ and one distant QSO at $z=6.13$. Methods. We use Herschel PACS and ALMA Band 7 observations of the hydroxyl molecule (OH) line at rest-frame wavelength 119 $μ$m which in absorption can provide unambiguous evidence for inflows or outflows of molecular ga… ▽ More Aims. We aim to search and characterize inflows and outflows of molecular gas in four ultraluminous infrared galaxies (ULIRGs) at $z\sim0.2-0.3$ and one distant QSO at $z=6.13$. Methods. We use Herschel PACS and ALMA Band 7 observations of the hydroxyl molecule (OH) line at rest-frame wavelength 119 $μ$m which in absorption can provide unambiguous evidence for inflows or outflows of molecular gas in nuclear regions of galaxies. Our study contributes to double the number of OH observations of luminous systems at $z\sim0.2-0.3$, and push the search for molecular outflows based on the OH transition to $z\sim6$. Results. We detect OH high-velocity absorption wings in three of the four ULIRGs. In two cases, IRAS F20036-1547 and IRAS F13352+6402, the blueshifted absorption profiles indicate the presence of powerful and fast molecular gas outflows. Consistent with an inside-out quenching scenario, these outflows are depleting the central reservoir of molecular gas at a similar rate than the intense star formation activity. In the case of the starburst-dominated system IRAS 10091+4704, we detect an inverted P-Cygni profile that is unique among ULIRGs and indicates the presence of a fast ($\sim400$ km s$^{-1}$) inflow of molecular gas at a rate of $\sim100~M_{\odot}~{\rm yr}^{-1}$ towards the central region. Finally, we tentatively detect ($\sim3σ$) the OH doublet in absorption in the $z=6.13$ QSO ULAS J131911+095051. The OH feature is blueshifted with a median velocity that suggests the presence of a molecular outflow, although characterized by a modest molecular mass loss rate of $\sim200~M_{\odot}~{\rm yr}^{-1}$. This value is comparable to the small mass outflow rates found in the stacking of the [CII] spectra of other $z\sim6$ QSOs and suggests that ejective feedback in this phase of the evolution of ULAS J131911+095051 has subsided. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: 9 pages, 6 figures, accepted for publication in A&A as a Letter to the Editor

AGN feedback in a galaxy merger: Multi-phase, galaxy-scale outflows including a fast molecular gas blob ~6 kpc away from IRAS F08572+3915

Authors: R. Herrera-Camus, A. Janssen, E. Sturm, D. Lutz, S. Veilleux, R. Davies, T. Shimizu, E. González-Alfonso, D. S. N. Rupke, L. Tacconi, R. Genzel, C. Cicone, R. Maiolino, A. Contursi, J. Graciá-Carpio

Abstract: To understand the role that AGN feedback plays in galaxy evolution we need in-depth studies of the multi-phase structure and energetics of galaxy-wide outflows. In this work we present new, deep ($\sim$50 hr) NOEMA CO(1-0) line observations of the molecular gas in the powerful outflow driven by the AGN in the ultra-luminous infrared galaxy IRAS F08572+3915. We spatially resolve the outflow, findin… ▽ More To understand the role that AGN feedback plays in galaxy evolution we need in-depth studies of the multi-phase structure and energetics of galaxy-wide outflows. In this work we present new, deep ($\sim$50 hr) NOEMA CO(1-0) line observations of the molecular gas in the powerful outflow driven by the AGN in the ultra-luminous infrared galaxy IRAS F08572+3915. We spatially resolve the outflow, finding that its most likely configuration is a wide-angle bicone aligned with the kinematic major axis of the rotation disk. The molecular gas in the wind reaches velocities up to approximately $\pm$1200 km s$^{-1}$ and transports nearly 20% of the molecular gas mass in the system. We detect a second outflow component located $\sim$6 kpc north-west from the galaxy moving away at $\sim$900 km s$^{-1}$, which could be the result of a previous episode of AGN activity. The total mass and energetics of the outflow, which includes contributions from the ionized, neutral, warm and cold molecular gas phases is strongly dominated by the cold molecular gas. In fact, the molecular mass outflow rate is higher than the star formation rate, even if we only consider the gas in the outflow that is fast enough to escape the galaxy, which accounts for about $\sim$40% of the total mass of the outflow. This results in an outflow depletion time for the molecular gas in the central $\sim$1.5 kpc region of only $\sim3$ Myr, a factor of $\sim2$ shorter than the depletion time by star formation activity. △ Less

Submitted 14 November, 2019; originally announced November 2019.

Comments: 17 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 635, A47 (2020)

Molecular outflows in local galaxies: Method comparison and a role of intermittent AGN driving

Authors: D. Lutz, E. Sturm, A. Janssen, S. Veilleux, S. Aalto, C. Cicone, A. Contursi, R. I. Davies, C. Feruglio, J. Fischer, S. Garcia-Burillo, R. Genzel, E. González-Alfonso, J. Gracía-Carpio R. Herrera-Camus, R. Maiolino, A. Schruba, T. Shimizu, A. Sternberg, L. J. Tacconi, A. Weiß

Abstract: We report new detections and limits from a NOEMA and ALMA CO(1-0) search for molecular outflows in 13 local galaxies with high FIR surface brightness, and combine with results from the literature. CO line ratios and outflow structure provide some constraints on the conversion from observables to quantities such as molecular mass outflow rates. Ratios between outflow emission in higher J CO transit… ▽ More We report new detections and limits from a NOEMA and ALMA CO(1-0) search for molecular outflows in 13 local galaxies with high FIR surface brightness, and combine with results from the literature. CO line ratios and outflow structure provide some constraints on the conversion from observables to quantities such as molecular mass outflow rates. Ratios between outflow emission in higher J CO transitions and in CO(1-0) typically are consistent with excitation Ri1<~1. For IRAS 13120-5453, however, R31=2.10 indicates optically thin CO in the outflow. Like much of the outflow literature, we use alpha(CO) = 0.8, and we present arguments for using C=1 in deriving molecular mass outflow rates Mdot = C*M*v/R. We compare the two main methods for molecular outflow detection: CO mm interferometry and Herschel OH spectroscopy. For 26 sources studied with both methods, we find 80% agreement in detecting vout>~150km/s outflows, and non-matches can be plausibly ascribed to outflow geometry and SNR. For 12 bright ULIRGs with detailed OH-based outflow modeling, CO outflows are detected in all but one. Outflow masses, velocities, and sizes for these 11 sources agree well between the two methods, and modest remaining differences may relate to the different but overlapping regions sampled by CO emission and OH absorption. Outflow properties correlate better with AGN luminosity and with bolometric luminosity than with FIR surface brightness. The most massive outflows are found for systems with current AGN activity, but significant outflows in non-AGN systems must relate to star formation or to AGN activity in the recent past. We report scaling relations for the increase of outflow mass, rate, momentum rate, and kinetic power with bolometric luminosity. Short ~10^6yr flow times and some sources with resolved multiple outflow episodes support a role of intermittent driving, likely by AGN. (abridged) △ Less

Submitted 13 November, 2019; originally announced November 2019.

Comments: 37 pages. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 633, A134 (2020)

Super Hot Cores in NGC 253: Witnessing the formation and early evolution of Super Star Clusters

Authors: Fernando Rico-Villas, Jesus Martin-Pintado, Eduardo Gonzalez-Alfonso, Sergio Martin, Victor M. Rivilla

Abstract: Using $0.2^{\prime \prime}$ ($\sim3$ pc) ALMA images of vibrationally excited HC$_3$N emission (HC$_3$N$^*$) we reveal the presence of $8$ unresolved Super Hot Cores (SHCs) in the inner $160$ pc of NGC\,253. Our LTE and non-LTE modelling of the HC$_3$N$^*$ emission indicate that SHCs have dust temperatures of $200-375$ K, relatively high H$_2$ densities of $1-6\times 10^{6}$ cm$^{-3}$ and high IR… ▽ More Using $0.2^{\prime \prime}$ ($\sim3$ pc) ALMA images of vibrationally excited HC$_3$N emission (HC$_3$N$^*$) we reveal the presence of $8$ unresolved Super Hot Cores (SHCs) in the inner $160$ pc of NGC\,253. Our LTE and non-LTE modelling of the HC$_3$N$^*$ emission indicate that SHCs have dust temperatures of $200-375$ K, relatively high H$_2$ densities of $1-6\times 10^{6}$ cm$^{-3}$ and high IR luminosities of $0.1-1\times 10^8$ L$_\odot$. As expected from their short lived phase ($\sim 10^4$ yr), all SHCs are associated with young Super Star Clusters (SSCs). We use the ratio of luminosities form the SHCs (protostar phase) and from the free-free emission (ZAMS star phase), to establish the evolutionary stage of the SSCs. The youngest SSCs, with the larges ratios, have ages of a few $10^4$ yr (proto-SSCs) and the more evolved SSCs are likely between $10^5$ and $10^6$ yr (ZAMS-SSCs). The different evolutionary stages of the SSCs are also supported by the radiative feedback from the UV radiation as traced by the HNCO/CS ratio, with this ratio being systematically higher in the young proto-SSCs than in the older ZAMS-SSCs. We also estimate the SFR and the SFE of the SSCs. The trend found in the estimated SFE ($\sim40\%$ for proto-SSCs and $>85\%$ for ZAMS-SSCs) and in the gas mass reservoir available for star formation, one order of magnitude higher for proto-SSCs, suggests that star formation is still going on in proto-SSCs. We also find that the most evolved SSCs are located, in projection, closer to the center of the galaxy than the younger proto-SSCs, indicating an inside-out SSC formation scenario. △ Less

Submitted 14 January, 2020; v1 submitted 25 September, 2019; originally announced September 2019.

Comments: typos corrected, references added

The Greenhouse Effect in Buried Galactic Nuclei and the Resonant HCN Vibrational Emission

Authors: Eduardo González-Alfonso, Kazushi Sakamoto

Abstract: Recent interferometric observations have shown bright HCN emission from the nu2=1 vibrational state arising in buried nuclear regions of galaxies, indicating an efficient pumping of the nu2=1 state through absorption of 14 $μ$m continuum photons. We have modeled the continuum and HCN vibrational line emission in these regions, characterized by high column densities of dust and high luminosities, w… ▽ More Recent interferometric observations have shown bright HCN emission from the nu2=1 vibrational state arising in buried nuclear regions of galaxies, indicating an efficient pumping of the nu2=1 state through absorption of 14 $μ$m continuum photons. We have modeled the continuum and HCN vibrational line emission in these regions, characterized by high column densities of dust and high luminosities, with a spherically symmetric approach, simulating both a central heating source (AGN) and a compact nuclear starburst (SB). We find that when the H2 columns become very high, N_{H2}>~10^{25} cm-2, trapping of continuum photons within the nuclear region dramatically enhances the dust temperature (Tdust) in the inner regions, even though the predicted spectral energy distribution as seen from outside becomes relatively cold. The models thus predict bright continuum at millimeter wavelengths for luminosity surface brightness (averaged over the model source) of ~10^{8} Lsun pc^{-2}. This {\it greenhouse} effect significantly enhances the mean mid-infrared intensity within the dusty volume, populating the nu2=1 state to the extent that the HCN vibrational lines become optically thick. AGN models yield higher Tdust in the inner regions and higher peak (sub)millimeter continuum brightness than SB models, but similar HCN vibrational J=3-2 and 4-3 emission owing to both optical depth effects and a moderate impact of high \tdust\ on these low-J lines. The observed HCN vibrational emission in several galaxies can be accounted for with a HCN abundance of ~10^{-6} (relative to H2) and luminosity surface brightness in the range (0.5-2)x10^{8}$ Lsun pc^{-2}, predicting a far-infrared photosphere with Tdust}~80-150 K --in agreement with the values inferred from far-infrared molecular absorption. △ Less

Submitted 12 August, 2019; originally announced August 2019.

Comments: Accepted for publication in The Astrophysical Journal

CO, H2O, H2O+ line and dust emission in a z = 3.63 strongly lensed starburst merger at sub-kiloparsec scales

Authors: C. Yang, R. Gavazzi, A. Beelen, P. Cox, A. Omont, M. D. Lehnert, Y. Gao, R. J. Ivison, A. M. Swinbank, L. Barcos-Muñoz, R. Neri, A. Cooray, S. Dye, S. Eales, H. Fu, E. González-Alfonso, E. Ibar, M. J. Michałowski, H. Nayyeri, M. Negrello, J. Nightingale, I. Pérez-Fournon, D. A. Riechers, I. Smail, P. van der Werf

Abstract: Using ALMA, we report high angular-resolution observations of the redshift z=3.63 galaxy, G09v1.97, one of the most luminous strongly lensed galaxies discovered by the H-ATLAS survey. We present 0"2-0"4 resolution images of the rest-frame 188 and 419$μ$m dust continuum and the CO(6-5), H2O(211-202) and J=2 H2O+ line emission. We also report the detection of H$_2^{18}$O in this source. The dust con… ▽ More Using ALMA, we report high angular-resolution observations of the redshift z=3.63 galaxy, G09v1.97, one of the most luminous strongly lensed galaxies discovered by the H-ATLAS survey. We present 0"2-0"4 resolution images of the rest-frame 188 and 419$μ$m dust continuum and the CO(6-5), H2O(211-202) and J=2 H2O+ line emission. We also report the detection of H$_2^{18}$O in this source. The dust continuum and molecular gas emission are resolved into a nearly complete ~1"5 diameter Einstein ring plus a weaker image in the center, which is caused by a special dual deflector lensing configuration. The observed line profiles of the CO, H2O and H2O+ lines are strikingly similar. In the source plane, we reconstruct the dust continuum images and the spectral cubes of the line emission at sub-kpc scales. The reconstructed dust emission in the source plane is dominated by a compact disk with an effective radius of 0.7kpc plus an overlapping extended disk with a radius twice as large. While the average magnification for the dust continuum is $μ$~10-11, the magnification of the line emission varies 5 to 22 across different velocity components. The emission lines have similar spatial and kinematic distributions. The molecular gas and dust content reveal that G09v1.97 is a gas-rich major merger in its pre-coalescence phase. Both of the merging companions are intrinsically ULIRGs with LIR reaching $\gtrsim 4\times10^{12}L_\odot$, and the total LIR of G09v1.97 is $1.4\times10^{13}L_\odot$. The approaching southern galaxy shows no obvious kinematic structure with a semi-major half-light radius a_s=0.4kpc, while the receding galaxy resembles an a_s=1.2kpc rotating disk. The two galaxies are separated by a projected distance of 1.3kpc, bridged by weak line emission that is co-spatially located with the cold-dust-emission peak, suggesting a large amount of cold ISM in the interacting region. (abridged) △ Less

Submitted 9 March, 2019; v1 submitted 1 March, 2019; originally announced March 2019.

Comments: 19 pages, 12 figures, and 4-page appendix, A&A in press. Minor changes to match language editing

Journal ref: A&A 624, A138 (2019)

Herschel water maps towards the vicinity of the black hole Sgr A*

Authors: J. Armijos-Abendaño, J. Martín-Pintado, M. A. Requena-Torres, E. González-Alfonso, R. Güsten, A. Weiß, A. I. Harris, F. P. Israel, C. Kramer, J. Stutzki, P. van der Werf

Abstract: Aims: We study the spatial distribution and kinematics of water emission in a ~64 pc$^2$ region of the Galactic Center (GC) around Sgr A*. We also analyze the water excitation to derive the physical conditions and water abundances in the CND and the `quiescent clouds'. Methods: We presented the integrated intensity maps of the ortho 1$_{10}-1_{01}$, and para 2$_{02}-1_{11}$ and 1$_{11}-0_{00}$ wat… ▽ More Aims: We study the spatial distribution and kinematics of water emission in a ~64 pc$^2$ region of the Galactic Center (GC) around Sgr A*. We also analyze the water excitation to derive the physical conditions and water abundances in the CND and the `quiescent clouds'. Methods: We presented the integrated intensity maps of the ortho 1$_{10}-1_{01}$, and para 2$_{02}-1_{11}$ and 1$_{11}-0_{00}$ water transitions observed with the HIFI instrument on board Herschel. To study the water excitation we used ground state ortho and para H$_2^{18}$O transitions. In our study, we also used SPIRE continuum measurements of the CND. Using a non-LTE radiative transfer code, the water line profiles and dust continuum were modeled. We also used a rotating ring model to reproduce the CND kinematics represented by the PV diagram. Results: We identify the water emission arising from the CND, the Western Streamer, and the 20 and 50 km s$^{-1}$ clouds. The ortho water maps show absorption structures in the range of [-220,10] km s$^{-1}$. The PV diagram shows that the 2$_{02}-1_{11}$ H$_2$O emission traces the CND. We derive high X$_{H_2O}$ of $\sim$(0.1-1.3)$\times$10$^{-5}$, V$_t$ of 14-23 km s$^{-1}$ and T$_d$ of 15-45 K for the CND, and the lower X$_{\rm H_2O}$ of 4$\times$10$^{-8}$ and V$_t$ of 9 km s$^{-1}$ for the 20 km s$^{-1}$ cloud. Collisional excitation and dust effects are responsible for the water excitation in the southwest lobe of the CND and the 20 km s$^{-1}$ cloud, whereas only collisions can account for the water excitation in the northeast lobe of the CND. We propose that the water vapor in the CND is caused by grain sputtering by shocks of 10-20 km s$^{-1}$, with some contribution of high temperature and cosmic-ray chemistries plus a PDR chemistry. The low X$_{\rm H_2O}$ derived for the 20 km s$^{-1}$ cloud could be partially a consequence of the water freeze-out on grains. △ Less

Submitted 13 February, 2019; originally announced February 2019.

Comments: 15 pages, 11 figures

Journal ref: A&A 624, A112 (2019)

Hidden or missing outflows in highly obscured galaxy nuclei?

Authors: N. Falstad, F. Hallqvist, S. Aalto, S. König, S. Muller, R. Aladro, F. Combes, A. S. Evans, G. A. Fuller, J. S. Gallagher, S. García-Burillo, E. González-Alfonso, T. R. Greve, C. Henkel, M. Imanishi, T. Izumi, J. G. Mangum, S. Martín, G. C. Privon, K. Sakamoto, S. Veilleux, P. P. van der Werf

Abstract: Understanding the nuclear growth and feedback processes in galaxies requires investigating their often obscured central regions. One way to do this is to use (sub)millimeter line emission from vibrationally excited HCN (HCN-vib). It has been suggested that the most intense HCN-vib emission from a galaxy is connected to a phase of nuclear growth that occurs before the nuclear feedback processes hav… ▽ More Understanding the nuclear growth and feedback processes in galaxies requires investigating their often obscured central regions. One way to do this is to use (sub)millimeter line emission from vibrationally excited HCN (HCN-vib). It has been suggested that the most intense HCN-vib emission from a galaxy is connected to a phase of nuclear growth that occurs before the nuclear feedback processes have been fully developed. We aim to investigate if there is a connection between the presence of strong HCN-vib emission and the development of feedback in (U)LIRGs. We collected literature and archival data to compare the luminosities of rotational lines of HCN-vib, normalized to the total infrared luminosity, to the median velocities of 119 μm OH absorption lines, potentially indicating outflows, in a total of 17 (U)LIRGs. The most HCN-vib luminous systems all lack signatures of significant molecular outflows in the far-infrared OH absorption lines. However, at least some of the systems with bright HCN-vib emission do have fast and collimated outflows that can be seen in spectral lines at longer wavelengths. We conclude that the galaxy nuclei with the highest L(HCN-vib)/L(IR) do not drive wide-angle outflows detectable using the median velocities of far-infrared OH absorption lines. It is possible that this is due to an orientation effect where sources which are oriented in such a way that their outflows are not along our line of sight also radiate a smaller proportion of their infrared luminosity in our direction. It could also be that massive wide-angle outflows destroy the deeply embedded regions responsible for bright HCN-vib emission, so that the two phenomena cannot coexist. This would strengthen the idea that vibrationally excited HCN traces a heavily obscured stage of evolution before nuclear feedback mechanisms are fully developed. △ Less

Submitted 20 January, 2019; originally announced January 2019.

Comments: 10 pages, 4 figures, Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 623, A29 (2019)

Probing the High-Redshift Universe with SPICA: Toward the Epoch of Reionization and Beyond

Authors: E. Egami, S. Gallerani, R. Schneider, A. Pallottini, L. Vallini, E. Sobacchi, A. Ferrara, S. Bianchi, M. Bocchio, S. Marassi, L. Armus, L. Spinoglio, A. W. Blain, M. Bradford, D. L. Clements, H. Dannerbauer, J. A. Fernández-Ontiveros, E. González-Alfonso, M. J. Griffin, C. Gruppioni, H. Kaneda, K. Kohno, S. C. Madden, H. Matsuhara, P. Najarro , et al. (14 additional authors not shown)

Abstract: With the recent discovery of a dozen dusty star-forming galaxies and around 30 quasars at z>5 that are hyper-luminous in the infrared ($μ$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $μ$ is a lensing magnification factor), the possibility has opened up for SPICA, the proposed ESA M5 mid-/far-infrared mission, to extend its spectroscopic studies toward the epoch of reionization and beyond. In this paper… ▽ More With the recent discovery of a dozen dusty star-forming galaxies and around 30 quasars at z>5 that are hyper-luminous in the infrared ($μ$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $μ$ is a lensing magnification factor), the possibility has opened up for SPICA, the proposed ESA M5 mid-/far-infrared mission, to extend its spectroscopic studies toward the epoch of reionization and beyond. In this paper, we examine the feasibility and scientific potential of such observations with SPICA's far-infrared spectrometer SAFARI, which will probe a spectral range (35-230 $μ$m) that will be unexplored by ALMA and JWST. Our simulations show that SAFARI is capable of delivering good-quality spectra for hyper-luminous infrared galaxies (HyLIRGs) at z=5-10, allowing us to sample spectral features in the rest-frame mid-infrared and to investigate a host of key scientific issues, such as the relative importance of star formation versus AGN, the hardness of the radiation field, the level of chemical enrichment, and the properties of the molecular gas. From a broader perspective, SAFARI offers the potential to open up a new frontier in the study of the early Universe, providing access to uniquely powerful spectral features for probing first-generation objects, such as the key cooling lines of low-metallicity or metal-free forming galaxies (fine-structure and H2 lines) and emission features of solid compounds freshly synthesized by Population III supernovae. Ultimately, SAFARI's ability to explore the high-redshift Universe will be determined by the availability of sufficiently bright targets (whether intrinsically luminous or gravitationally lensed). With its launch expected around 2030, SPICA is ideally positioned to take full advantage of upcoming wide-field surveys such as LSST, SKA, Euclid, and WFIRST, which are likely to provide extraordinary targets for SAFARI. △ Less

Submitted 25 September, 2018; v1 submitted 20 September, 2018; originally announced September 2018.

Comments: 24 pages, 8 figures, accepted for publication in PASA (SPICA special issue)

Molecular gas in the northern nucleus of Mrk273: Physical and chemical properties of the disk and its outflow

Authors: R. Aladro, S. König, S. Aalto, E. González-Alfonso, N. Falstad, S. Martín, S. Muller, S. García-Burillo, C. Henkel, P. van der Werf, E. Mills, J. Fischer, F. Costagliola, M. Krips, .

Abstract: Aiming to characterise the properties of the molecular gas in the ultraluminous infrared galaxy Mrk273 and its outflow, we used the NOEMA interferometer to image the dense gas molecular tracers HCN, HCO+, HNC, HOC+ and HC3N at 86GHz and 256GHz with angular resolutions of 4.9x4.5 arcsec (3.7x3.4 kpc) and 0.61x0.55 arcsec (460x420 pc). We also modelled the flux of several H2O lines observed with Her… ▽ More Aiming to characterise the properties of the molecular gas in the ultraluminous infrared galaxy Mrk273 and its outflow, we used the NOEMA interferometer to image the dense gas molecular tracers HCN, HCO+, HNC, HOC+ and HC3N at 86GHz and 256GHz with angular resolutions of 4.9x4.5 arcsec (3.7x3.4 kpc) and 0.61x0.55 arcsec (460x420 pc). We also modelled the flux of several H2O lines observed with Herschel using a radiative transfer code that includes excitation by collisions as well as by far-infrared photons. The disk of the Mrk273 north nucleus has two components with decoupled kinematics. The gas in the outer parts (1.5 kpc) rotates with a south-east to north-west direction, while in the inner disk (300 pc) follows a north-east to south-west rotation. The central 300 pc, which hosts a compact starburst region, is filled with dense and warm gas, contains a dynamical mass of (4-5)x10^9M_sun, a luminosity of L'_HCN=(3-4)x10^8 K km/s pc^2, and a dust temperature of 55 K. At the very centre, a compact core with R~50 pc has a luminosity of L_IR=4x10^11L_sun (30% of the total infrared luminosity), and a dust temperature of 95 K. The core is expanding at low velocities ~50-100 km/s, probably affected by the outflowing gas. We detect the blue-shifted component of the outflow, while the red-shifted counterpart remains undetected in our data. Its cold and dense phase reaches fast velocities up to ~1000 km/s, while the warm outflowing gas has more moderate maximum velocities of ~600 km/s. The outflow is detected as far as 460 pc from the centre in the northern direction, and has a mass of dense gas <8x10^8M_sun. The difference between the position angles of the inner disk (~70 degree) and the outflow (~10 degree) indicates that the outflow is likely powered by the AGN, and not by the starburst. Regarding the chemistry, we measure an extremely low HCO+/HOC+ ratio of 10+-5 in the inner disk of Mrk273. △ Less

Submitted 30 May, 2018; v1 submitted 29 May, 2018; originally announced May 2018.

Comments: Accepted for publication in A&A. 21 pages, 17 figures, 7 tables, and a lot of interesting text

Journal ref: A&A 617, A20 (2018)

The dual role of starburst and active galactic nuclei in driving extreme molecular outflows

Authors: Avani Gowardhan, Henrik Spoon, Dominik A. Riechers, Eduardo González-Alfonso, Duncan Farrah, Jacqueline Fischer, Jeremy Darling, Chiara Fergulio, Jose Afonso, Luca Bizzocchi

Abstract: We report molecular gas observations of IRAS 20100-4156 and IRAS 03158+4227, two local ultraluminous infrared galaxies (ULIRGs) hosting some of the fastest and most massive molecular outflows known. Using ALMA and PdBI observations, we spatially resolve the CO(1-0) emission from the outflowing molecular gas in both and find maximum outflow velocities of $ v_{\rm max} \sim 1600$ and $\sim 1700$ km/… ▽ More We report molecular gas observations of IRAS 20100-4156 and IRAS 03158+4227, two local ultraluminous infrared galaxies (ULIRGs) hosting some of the fastest and most massive molecular outflows known. Using ALMA and PdBI observations, we spatially resolve the CO(1-0) emission from the outflowing molecular gas in both and find maximum outflow velocities of $ v_{\rm max} \sim 1600$ and $\sim 1700$ km/s for IRAS 20100-4156 and IRAS 03158+4227, respectively. We find total gas mass outflow rates of $\dot M_{\rm OF} \sim 670$ and $\sim 350$ Msun/yr, respectively, corresponding to molecular gas depletion timescales $τ^{\rm dep}_{\rm OF} \sim 11$ and $\sim 16$ Myr. This is nearly 3 times shorter than the depletion timescales implied by star formation, $τ^{\rm dep}_{\rm SFR} \sim 33$ and $\sim 46$ Myr, respectively. To determine the outflow driving mechanism, we compare the starburst ($L_{*}$) and AGN ($L_{\rm AGN}$) luminosities to the outflowing energy and momentum fluxes, using mid-infrared spectral decomposition to discern $L_{\rm AGN}$. Comparison to other molecular outflows in ULIRGs reveals that outflow properties correlate similarly with $L_{*}$ and $L_{\rm IR}$ as with $L_{\rm AGN}$, indicating that AGN luminosity alone may not be a good tracer of feedback strength and that a combination of AGN and starburst activity may be driving the most powerful molecular outflows. We also detect the OH 1.667 GHz maser line from both sources and demonstrate its utility in detecting molecular outflows. △ Less

Submitted 9 April, 2018; originally announced April 2018.

Comments: 29 pages, 17 figures, 6 tables, accepted to ApJ

Outflowing OH$^+$ in Markarian 231: the ionization rate of the molecular gas

Authors: E. González-Alfonso, J. Fischer, S. Bruderer, M. L. N. Ashby, H. A. Smith, S. Veilleux, H. S. P. Müller, K. P. Stewart, E. Sturm

Abstract: The oxygen-bearing molecular ions OH+, H2O+, and H3O+ are key species that probe the ionization rate of (partially) molecular gas that is ionized by X-rays and cosmic rays permeating the interstellar medium. We report Herschel far-infrared and submillimeter spectroscopic observations of OH+ in Mrk 231, showing both ground-state P-Cygni profiles, and excited line profiles with blueshifted absorptio… ▽ More The oxygen-bearing molecular ions OH+, H2O+, and H3O+ are key species that probe the ionization rate of (partially) molecular gas that is ionized by X-rays and cosmic rays permeating the interstellar medium. We report Herschel far-infrared and submillimeter spectroscopic observations of OH+ in Mrk 231, showing both ground-state P-Cygni profiles, and excited line profiles with blueshifted absorption wings extending up to ~1000 km s^{-1}. In addition, OH+ probes an excited component peaking at central velocities, likely arising from the torus probed by the OH centimeter-wave megamaser. Four lines of H2O+ are also detected at systemic velocities, but H3O+ is undetected. Based on our earlier OH studies, we estimate an abundance ratio of OH/OH+~5-10 for the outflowing components and ~20 for the torus, and an OH+ abundance relative to H nuclei of ~>10^{-7}. We also find high OH+/H2O+ and OH+/H3O+ ratios, both are ~>4 in the torus and ~>10-20 in the outflowing gas components. Chemical models indicate that these high OH+ abundances relative to OH, H2O+, and H3O+ are characteristic of gas with a high ionization rate per unit density, ζ/n_H~(1-5)x10^{-17} cm^3 s^{-1} and ~(1-2)x10^{-16} cm^3 s^{-1} for the above components, respectively, and an ionization rate of ζ~(0.5-2)x10^{-12} s^{-1}. X-rays appear to be unable to explain the inferred ionization rate, and thus we suggest that low-energy (10-400 MeV) cosmic-rays are primarily responsible for the ionization with \dot{M}_{CR}~0.01 M_{sun} yr^{-1} and \dot{E}_{CR}~10^{44} erg s^{-1}, the latter corresponding to 1% of the AGN luminosity and similar to the energetics of the molecular outflow. We suggest that cosmic-rays accelerated in the forward shock associated with the molecular outflow are responsible for the ionization, as they diffuse through the outflowing molecular phase downstream. △ Less

Submitted 13 March, 2018; originally announced March 2018.

Comments: 25 pages, 11 figures

SHINING, A Survey of Far Infrared Lines in Nearby Galaxies. II: Line-Deficit Models, AGN impact, [CII]-SFR Scaling Relations, and Mass-Metallicity Relation in (U)LIRGS

Authors: R. Herrera-Camus, E. Sturm, J. Graciá-Carpio, D. Lutz, A. Contursi, S. Veilleux, J. Fischer, E. González-Alfonso, A. Poglitsch, L. Tacconi, R. Genzel, R. Maiolino, A. Sternberg, R. Davies, A. Verma

Abstract: The SHINING survey (Paper I; Herrera-Camus et al. 2018) offers a great opportunity to study the properties of the ionized and neutral media of galaxies from prototypical starbursts and active galactic nuclei (AGN) to heavily obscured objects. Based on Herschel/PACS observations of the main far-infrared (FIR) fine-structure lines, in this paper we analyze the physical mechanisms behind the observed… ▽ More The SHINING survey (Paper I; Herrera-Camus et al. 2018) offers a great opportunity to study the properties of the ionized and neutral media of galaxies from prototypical starbursts and active galactic nuclei (AGN) to heavily obscured objects. Based on Herschel/PACS observations of the main far-infrared (FIR) fine-structure lines, in this paper we analyze the physical mechanisms behind the observed line deficits in galaxies, the apparent offset of luminous infrared galaxies (LIRGs) from the mass-metallicity relation, and the scaling relations between [CII] 158 $μ$m line emission and star formation rate (SFR). Based on a toy model and the Cloudy code, we conclude that the increase in the ionization parameter with FIR surface brightness can explain the observed decrease in the line-to-FIR continuum ratio of galaxies. In the case of the [CII] line, the increase in the ionization parameter is accompanied by a reduction in the photoelectric heating efficiency and the inability of the line to track the increase in the FUV radiation field as galaxies become more compact and luminous. In the central $\sim$kiloparsec regions of AGN galaxies we observe a significant increase in the [OI] 63 $μ$m/[CII] line ratio; the AGN impact on the line-to-FIR ratios fades on global scales. Based on extinction-insensitive metallicity measurements of LIRGs we confirm that they lie below the mass-metallicity relation, but the offset is smaller than those reported in studies that use optical-based metal abundances. Finally, we present scaling relations between [CII] emission and SFR in the context of the main-sequence of star-forming galaxies. △ Less

Submitted 12 March, 2018; originally announced March 2018.

Comments: 22 pages, 7 Figures, Accepted for publication in the Astrophysical Journal

SHINING, A Survey of Far Infrared Lines in Nearby Galaxies. I: Survey Description, Observational Trends, and Line Diagnostics

Authors: R. Herrera-Camus, E. Sturm, J. Graciá-Carpio, D. Lutz, A. Contursi, S. Veilleux, J. Fischer, E. González-Alfonso, A. Poglitsch, L. Tacconi, R. Genzel, R. Maiolino, A. Sternberg, R. Davies, A. Verma

Abstract: We use the Herschel/PACS spectrometer to study the global and spatially resolved far-infrared (FIR) fine-structure line emission in a sample of 52 galaxies that constitute the SHINING survey. These galaxies include star-forming, active-galactic nuclei (AGN), and luminous infrared galaxies (LIRGs). We find an increasing number of galaxies (and kiloparsec size regions within galaxies) with low line-… ▽ More We use the Herschel/PACS spectrometer to study the global and spatially resolved far-infrared (FIR) fine-structure line emission in a sample of 52 galaxies that constitute the SHINING survey. These galaxies include star-forming, active-galactic nuclei (AGN), and luminous infrared galaxies (LIRGs). We find an increasing number of galaxies (and kiloparsec size regions within galaxies) with low line-to-FIR continuum ratios as a function of increasing FIR luminosity ($L_{\mathrm{FIR}}$), dust infrared color, $L_{\mathrm{FIR}}$ to molecular gas mass ratio ($L_{\mathrm{FIR}}/M_{\mathrm{mol}}$), and FIR surface brightness ($Σ_{\mathrm{FIR}}$). The correlations between the [CII]/FIR or [OI]/FIR ratios with $Σ_{\mathrm{FIR}}$ are remarkably tight ($\sim0.3$ dex scatter over almost four orders of magnitude in $Σ_{\mathrm{FIR}}$). We observe that galaxies with $L_{\mathrm{FIR}}/M_{\mathrm{mol}} \gtrsim 80\,L_{\odot}\,M_{\odot}^{-1}$ and $Σ_{\mathrm{FIR}}\gtrsim10^{11}$ $L_{\odot}$ kpc$^{-2}$ tend to have weak fine-structure line-to-FIR continuum ratios, and that LIRGs with infrared sizes $\gtrsim1$ kpc have line-to-FIR ratios comparable to those observed in typical star-forming galaxies. We analyze the physical mechanisms driving these trends in Paper II (Herrera-Camus et al. 2018). The combined analysis of the [CII], [NII], and [OIII] lines reveals that the fraction of the [CII] line emission that arises from neutral gas increases from 60% to 90% in the most active star-forming regions and that the emission originating in the ionized gas is associated with low-ionization, diffuse gas rather than with dense gas in HII regions. Finally, we report the global and spatially resolved line fluxes of the SHINING galaxies to enable the comparison and planning of future local and high-$z$ studies. △ Less

Submitted 12 March, 2018; originally announced March 2018.

Comments: 30 pages, 13 Figures, Accepted for publication in the Astrophysical Journal

HERUS: The Far-IR/Submm Spectral Energy Distributions of Local ULIRGs & Photometric Atlas

Authors: D. L. Clements, C. Pearson, D. Farrah, J. Greenslade, Jeronimo Bernard-Salas, E. Gonzalez-Alfonso, J. Afonso, A. Efstathiou, D. Rigopoulou, V. Lebouteiller, P. D. Hurley, H. Spoon

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

Submitted 13 December, 2017; originally announced December 2017.

Comments: Accepted for publication in MNRAS

Probing the baryon cycle of galaxies with SPICA mid- and far-infrared observations

Authors: F. F. S. van der Tak, S. C. Madden, P. Roelfsema, L. Armus, M. Baes, J. Bernard-Salas, A. Bolatto, S. Bontemps, C. Bot, C. M. Bradford, J. Braine, L. Ciesla, D. Clements, D. Cormier, J. A. Fernández-Ontiveros, F. Galliano, M. Giard, H. Gomez, E. González-Alfonso, F. Herpin, D. Johnstone, A. Jones, H. Kaneda, F. Kemper, V. Lebouteiller , et al. (7 additional authors not shown)

Abstract: The SPICA mid and far-infrared telescope will address fundamental issues in our understanding of star formation and ISM physics in galaxies. A particular hallmark of SPICA is the outstanding sensitivity enabled by the cold telescope, optimized detectors, and wide instantaneous bandwidth throughout the mid- and far-infrared. The spectroscopic, imaging and polarimetric observations that SPICA will b… ▽ More The SPICA mid and far-infrared telescope will address fundamental issues in our understanding of star formation and ISM physics in galaxies. A particular hallmark of SPICA is the outstanding sensitivity enabled by the cold telescope, optimized detectors, and wide instantaneous bandwidth throughout the mid- and far-infrared. The spectroscopic, imaging and polarimetric observations that SPICA will be able to collect will help in clarifying the complex physical mechanisms which underlie the baryon cycle of galaxies. In particular: (i) The access to a large suite of atomic and ionic fine-structure lines for large samples of galaxies will shed light on the origin of the observed spread in star formation rates within and between galaxies. (ii) Observations of HD rotational lines (out to $\sim$10 Mpc) and fine structure lines such as [CII] 158 $μ$m (out to $\sim$100 Mpc) will clarify the main reservoirs of interstellar matter in galaxies, including phases where CO does not emit. (iii) Far-infrared spectroscopy of dust and ice features will address uncertainties in the mass and composition of dust in galaxies, and the contributions of supernovae to the interstellar dust budget will be quantified by photometry and monitoring of supernova remnants in nearby galaxies. (iv) Observations of far-infrared cooling lines such as [OI] 63 $μ$m from star-forming molecular clouds in our Galaxy will evaluate the importance of shocks to dissipate turbulent energy. The paper concludes with requirements for the telescope and instruments, and recommendations for the observing strategy. △ Less

Submitted 30 November, 2017; originally announced November 2017.

Comments: Accepted for publication in PASA; 19 pages, 8 figures, 1 table

A hidden molecular outflow in the LIRG Zw 049.057

Authors: N. Falstad, S. Aalto, J. G. Mangum, F. Costagliola, J. S. Gallagher, E. González-Alfonso, K. Sakamoto, S. König, S. Muller, A. S. Evans, G. C. Privon

Abstract: Feedback in the form of mass outflows driven by star formation or active galactic nuclei is a key component of galaxy evolution. The luminous infrared galaxy Zw 049.057 harbours a compact obscured nucleus with a possible far-IR signature of outflowing molecular gas. Due to the high optical depths at far-IR wavelengths, the interpretation of the outflow signature is uncertain. At mm and radio wavel… ▽ More Feedback in the form of mass outflows driven by star formation or active galactic nuclei is a key component of galaxy evolution. The luminous infrared galaxy Zw 049.057 harbours a compact obscured nucleus with a possible far-IR signature of outflowing molecular gas. Due to the high optical depths at far-IR wavelengths, the interpretation of the outflow signature is uncertain. At mm and radio wavelengths, the radiation is better able to penetrate the large columns of gas and dust. We used high resolution observations from the SMA, ALMA, and the VLA to image the CO 2-1 and 6-5 emission, the 690 GHz continuum, the radio cm continuum, and absorptions by rotationally excited OH. The CO line profiles exhibit wings extending 300 km/s beyond the systemic velocity. At cm wavelengths, we find a compact (40 pc) continuum component in the nucleus, with weaker emission extending several 100 pc approximately along the major and minor axes of the galaxy. In the OH absorption lines toward the compact continuum, wings extending to a similar velocity as for the CO are seen on the blue side of the profile. The weak cm continuum emission along the minor axis is aligned with a highly collimated, jet-like dust feature previously seen in near-IR images of the galaxy. Comparison of the apparent optical depths in the OH lines indicate that the excitation conditions in Zw 049.057 differ from those in other OH megamaser galaxies. We interpret the wings in the spectral lines as signatures of a molecular outflow. A relation between this outflow and the minor axis radio feature is possible, although further studies are required to investigate this possible association and understand the connection between the outflow and the nuclear activity. Finally, we suggest that the differing OH excitation conditions are further evidence that Zw 049.057 is in a transition phase between megamaser and kilomaser activity. △ Less

Submitted 14 November, 2017; originally announced November 2017.

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

Journal ref: A&A 609, A75 (2018)

Feedback and feeding in the context of galaxy evolution with SPICA: direct characterization of molecular outflows and inflows

Authors: E. González-Alfonso, L. Armus, F. J. Carrera, V. Charmandaris, A. Efstathiou, E. Egami, J. A. Fernández-Ontiveros, J. Fischer, G. L. Granato, C. Gruppioni, E. Hatziminaoglou, M. Imanishi, N. Isobe, H. Kaneda, D. Koziel-Wierzbowska, M. A. Malkan, J. Martin-Pintado, S. Mateos, H. Matsuhara, G. Miniutti, T. Nakagawa, F. Pozzi, F. Rico-Villas, G. Rodighiero, P. Roelfsema , et al. (6 additional authors not shown)

Abstract: A far-infrared observatory such as the {\it SPace Infrared telescope for Cosmology and Astrophysics} ({\it SPICA}), with its unprecedented spectroscopic sensitivity, would unveil the role of feedback in galaxy evolution during the last $\sim10$ Gyr of the Universe ($z=1.5-2$), through the use of far- and mid-infrared molecular and ionic fine structure lines that trace outflowing and infalling gas.… ▽ More A far-infrared observatory such as the {\it SPace Infrared telescope for Cosmology and Astrophysics} ({\it SPICA}), with its unprecedented spectroscopic sensitivity, would unveil the role of feedback in galaxy evolution during the last $\sim10$ Gyr of the Universe ($z=1.5-2$), through the use of far- and mid-infrared molecular and ionic fine structure lines that trace outflowing and infalling gas. Outflowing gas is identified in the far-infrared through P-Cygni line shapes and absorption blueshifted wings in molecular lines with high dipolar moments, and through emission line wings of fine-structure lines of ionized gas. We quantify the detectability of galaxy-scale massive molecular and ionized outflows as a function of redshift in AGN-dominated, starburst-dominated, and main-sequence galaxies, explore the detectability of metal-rich inflows in the local Universe, and describe the most significant synergies with other current and future observatories that will measure feedback in galaxies via complementary tracers at other wavelengths. △ Less

Submitted 6 October, 2017; originally announced October 2017.

Comments: This paper belongs to the SPICA Special Issue on PASA. Accepted for publication in PASA

Tracing the evolution of dust obscured star-formation and accretion back to the reionisation epoch with SPICA

Authors: C. Gruppioni, L. Ciesla, E. Hatziminaoglou, F. Pozzi, G. Rodighiero, P. Santini, L. Armus, M. Baes, J. Braine, V. Charmandaris, D. L. Clements, N. Christopher, H. Dannerbauer, A. Efstathiou, E. Egami, J. A. Fernandez-Ontiveros, F. Fontanot, A. Franceschini, E. Gonzalez-Alfonso, M. Griffin, H. Kaneda, L. Marchetti, P. Monaco, T. Nakagawa, T. Onaka , et al. (13 additional authors not shown)

Abstract: Our current knowledge of star formation and accretion luminosity at high-redshift (z>3-4), as well as the possible connections between them, relies mostly on observations in the rest-frame ultraviolet (UV), which are strongly affected by dust obscuration. Due to the lack of sensitivity of past and current infrared (IR) instrumentation, so far it has not been possible to get a glimpse into the earl… ▽ More Our current knowledge of star formation and accretion luminosity at high-redshift (z>3-4), as well as the possible connections between them, relies mostly on observations in the rest-frame ultraviolet (UV), which are strongly affected by dust obscuration. Due to the lack of sensitivity of past and current infrared (IR) instrumentation, so far it has not been possible to get a glimpse into the early phases of the dust-obscured Universe. Among the next generation of IR observatories, SPICA, observing in the 12-350 micron range, will be the only facility that can enable us to make the required leap forward in understanding the obscured star-formation rate and black-hole accretion rate densities (SFRD and BHARD, respectively) with respect to what Spitzer and Herschel achieved in the mid- and far-IR at z<3. In particular, SPICA will have the unique ability to trace the evolution of the obscured SFRD and BHARD over cosmic time, from the peak of their activity back to the reionisation epoch (i.e., 3<z<6-7), where its predecessors had severe limitations. Here we discuss the potential of both deep and shallow photometric surveys performed with the SPICA mid-IR instrument (SMI), enabled by the very low level of impact of dust obscuration in a band centred at 34 micron. These unique unbiased photometric surveys that SPICA will perform will be followed up by observations both with the SPICA spectrometers and with other facilities at shorter and longer wavelengths, with the aim to fully characterise the evolution of AGNs and star-forming galaxies after re-ionisation. △ Less

Submitted 6 October, 2017; originally announced October 2017.

Comments: This paper belongs to the SPICA Special Issue on PASA

Galaxy evolution studies with the SPace IR telescope for Cosmology and Astrophysics (SPICA): the power of IR spectroscopy

Authors: L. Spinoglio, A. Alonso-Herrero, L. Armus, M. Baes, J. Bernard-Salas, S. Bianchi, M. Bocchio, A. Bolatto, C. M. Bradford, J. Braine, F. J. Carrera, L. Ciesla, D. L. Clements, H. Dannerbauer, Y. Doi, A. Efstathiou, E. Egami, J. A. Fernandez-Ontiveros, A. Ferrara, J. Fischer, A. Franceschini, S. Gallerani, M. Giard, E. Gonzalez-Alfonso, C. Gruppioni , et al. (43 additional authors not shown)

Abstract: IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Larg… ▽ More IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features. △ Less

Submitted 5 October, 2017; originally announced October 2017.

Comments: Paper accepted for publication on PASA on 21st September 2017, as part of the SPICA Special Issue

SPICA and the Chemical Evolution of Galaxies: The Rise of Metals and Dust

Authors: J. A. Fernández-Ontiveros, L. Armus, M. Baes, J. Bernard-Salas, A. D. Bolatto, J. Braine, L. Ciesla, I. De Looze, E. Egami, J. Fischer, M. Giard, E. González-Alfonso, G. L. Granato, C. Gruppioni, M. Imanishi, D. Ishihara, H. Kaneda, S. Madden, M. Malkan, H. Matsuhara, M. Matsuura, T. Nagao, F. Najarro, T. Nakagawa, T. Onaka , et al. (11 additional authors not shown)

Abstract: The physical processes driving the chemical evolution of galaxies in the last $\sim 11\, \rm{Gyr}$ cannot be understood without directly probing the dust-obscured phase of star-forming galaxies and active galactic nuclei. This phase, hidden to optical tracers, represents the bulk of star formation and black hole accretion activity in galaxies at $1 < z < 3$. Spectroscopic observations with a cryog… ▽ More The physical processes driving the chemical evolution of galaxies in the last $\sim 11\, \rm{Gyr}$ cannot be understood without directly probing the dust-obscured phase of star-forming galaxies and active galactic nuclei. This phase, hidden to optical tracers, represents the bulk of star formation and black hole accretion activity in galaxies at $1 < z < 3$. Spectroscopic observations with a cryogenic infrared (IR) observatory like SPICA will be sensitive enough to peer through the dust-obscured regions of galaxies and access the rest-frame mid- to far-IR range in galaxies at high-$z$. This wavelength range contains a unique suite of spectral lines and dust features that serve as proxies for the abundances of heavy elements and the dust composition, providing tracers with a feeble response to both extinction and temperature. In this work, we investigate how SPICA observations could be exploited to understand key aspects in the chemical evolution of galaxies: the assembly of nearby galaxies based on the spatial distribution of heavy element abundances, the global content of metals in galaxies reaching the knee of the luminosity function up to $z \sim 3$, and the dust composition of galaxies at high-$z$. Possible synergies with facilities available in the late 2020s are also discussed. △ Less

Submitted 5 October, 2017; originally announced October 2017.

Comments: In press. This paper belongs to the "SPICA Special Issue" on PASA

Molecular gas in the Herschel-selected strongly lensed submillimeter galaxies at z~2-4 as probed by multi-J CO lines

Authors: C. Yang, A. Omont, A. Beelen, Y. Gao, P. van der Werf, R. Gavazzi, Z. -Y. Zhang, R. Ivison, M. Lehnert, D. Liu, I. Oteo, E. González-Alfonso, H. Dannerbauer, P. Cox, M. Krips, R. Neri, D. Riechers, A. J. Baker, M. J. Michałowski, A. Cooray, I. Smail

Abstract: (abridged) We present the IRAM-30m observations of multiple-J CO and CI line emission in a sample of redshift ~2-4 Herschel-ATLAS SMGs. A non-negligible effect of differential lensing is found for the CO emission lines, which could have caused significant underestimations of the linewidths, hence of the dynamical masses. The CO SLEDs are found to be similar to those of the local starburst-dominate… ▽ More (abridged) We present the IRAM-30m observations of multiple-J CO and CI line emission in a sample of redshift ~2-4 Herschel-ATLAS SMGs. A non-negligible effect of differential lensing is found for the CO emission lines, which could have caused significant underestimations of the linewidths, hence of the dynamical masses. The CO SLEDs are found to be similar to those of the local starburst-dominated ULIRGs and of the previously studied SMGs. After correcting for lensing amplification, we derived the global properties of the bulk of molecular gas in the SMGs using non-LTE radiative transfer modelling. The gas thermal pressure is found to be correlated with star formation efficiency. Further decomposing the CO SLEDs into two excitation components, we find a low-excitation component, which is less correlated with star formation, and a high-excitation one which is tightly related to the on-going star-forming activity. Additionally, tight linear correlations between the FIR and CO line luminosities have been confirmed for the $J \ge 5$ CO lines, implying that these CO lines are good tracers of star formation. The [CI](2-1) lines follow the tight linear correlation between the luminosities of the [CI](2-1) and the CO(1-0) line found in local starbursts, indicating that CI lines could serve as good total molecular gas mass tracers for high-redshift SMGs. The total mass of the molecular gas reservoir, $(1-30) \times 10^{10} M_\odot$, suggests a typical molecular gas depletion time ~20-100 Myr and a gas to dust mass ratio $δ_{\rm GDR}$~30-100. The ratio between CO line luminosity and the dust mass appears to be slowly increasing with redshift for the SMGs, which need to be further confirmed. Finally, through comparing the linewidth of CO and H2O lines, we find that they agree well in almost all our SMGs, confirming that the emitting regions are co-spatially located. △ Less

Submitted 20 September, 2017; v1 submitted 14 September, 2017; originally announced September 2017.

Comments: A&A in press, 22 pages plus 20 pages of appendix (1 long table plus figures), 13 figures and 7 tables in main text; minor edits

Journal ref: A&A 608, A144 (2017)

Quasar Feedback in the Ultraluminous Infrared Galaxy F11119+3257: Connecting the Accretion Disk Wind with the Large-Scale Molecular Outflow

Authors: S. Veilleux, A. Bolatto, F. Tombesi, M. Melendez, E. Sturm, E. Gonzalez-Alfonso, J. Fischer, D. S. N. Rupke

Abstract: In Tombesi et al. (2015), we reported the first direct evidence for a quasar accretion disk wind driving a massive molecular outflow. The target was F11119+3257, an ultraluminous infrared galaxy (ULIRG) with unambiguous type-1 quasar optical broad emission lines. The energetics of the accretion disk wind and molecular outflow were found to be consistent with the predictions of quasar feedback mode… ▽ More In Tombesi et al. (2015), we reported the first direct evidence for a quasar accretion disk wind driving a massive molecular outflow. The target was F11119+3257, an ultraluminous infrared galaxy (ULIRG) with unambiguous type-1 quasar optical broad emission lines. The energetics of the accretion disk wind and molecular outflow were found to be consistent with the predictions of quasar feedback models where the molecular outflow is driven by a hot energy-conserving bubble inflated by the inner quasar accretion disk wind. However, this conclusion was uncertain because the energetics were estimated from the optically thick OH 119 um transition profile observed with Herschel. Here, we independently confirm the presence of the molecular outflow in F11119+3257, based on the detection of broad wings in the CO(1-0) profile derived from ALMA observations. The broad CO(1-0) line emission appears to be spatially extended on a scale of at least ~7 kpc from the center. Mass outflow rate, momentum flux, and mechanical power of (80-200) R_7^{-1} M_sun/yr, (1.5-3.0) R_7^{-1} L_AGN/c, and (0.15-0.40)% R_7^{-1} L_AGN are inferred from these data, assuming a CO-to-H_2 conversion factor appropriate for a ULIRG (R_7 is the radius of the outflow normalized to 7 kpc and L_AGN is the AGN luminosity). These rates are time-averaged over a flow time scale of 7x10^6 yrs. They are similar to the OH-based rates time-averaged over a flow time scale of 4x10^5 yrs, but about a factor 4 smaller than the local ("instantaneous"; <10^5 yrs) OH-based estimates cited in Tombesi et al. The implications of these new results are discussed in the context of time-variable quasar-mode feedback and galaxy evolution. The need for an energy-conserving bubble to explain the molecular outflow is also re-examined. △ Less

Submitted 1 June, 2017; originally announced June 2017.

Comments: 15 pages, 6 figures, 4 tables, accepted for publication in ApJ

First detection of the 448 GHz H2O transition in space

Authors: M. Pereira-Santaella, E. González-Alfonso, A. Usero, S. García-Burillo, J. Martín-Pintado, L. Colina, A. Alonso-Herrero, S. Arribas, S. Cazzoli, F. Rico, D. Rigopoulou, T. Storchi Bergmann

Abstract: We present the first detection of the ortho-H2O 4_23-3_30 transition at 448 GHz in space. We observed this transition in the local (z = 0.010) luminous infrared (IR) galaxy ESO 320-G030 (IRAS F11506-3851) using the Atacama Large Millimeter/submillimeter Array (ALMA). The water 4_23-3_30 emission, which originates in the highly obscured nucleus of this galaxy, is spatially resolved over a region of… ▽ More We present the first detection of the ortho-H2O 4_23-3_30 transition at 448 GHz in space. We observed this transition in the local (z = 0.010) luminous infrared (IR) galaxy ESO 320-G030 (IRAS F11506-3851) using the Atacama Large Millimeter/submillimeter Array (ALMA). The water 4_23-3_30 emission, which originates in the highly obscured nucleus of this galaxy, is spatially resolved over a region of ~65 pc in diameter and shows a regular rotation pattern compatible with the global molecular and ionized gas kinematics. The line profile is symmetric and well fitted by a Gaussian with an integrated flux of 37.0 +- 0.7 Jy km s-1 . Models predict this water transition as a potential collisionally excited maser transition. On the contrary, in this galaxy, we find that the 4_23-3_30 emission is primarily excited by the intense far-IR radiation field present in its nucleus. According to our modeling, this transition is a probe of deeply buried galaxy nuclei thanks to the high dust optical depths (tau_100μm > 1, N_H > 1e24 cm-2) required to efficiently excite it. △ Less

Submitted 18 April, 2017; originally announced April 2017.

Comments: Accepted for publication in A&A Letters; 4 pages, 5 figures

Journal ref: A&A 601, L3 (2017)

Molecular outflows in local ULIRGs: energetics from multi-transition OH analysis

Authors: E. González-Alfonso, J. Fischer, H. W. W. Spoon, K. P. Stewart, M. L. N. Ashby, S. Veilleux, H. A. Smith, E. Sturm, D. Farrah, N. Falstad, M. Meléndez, J. Graciá-Carpio, A. W. Janssen, V. Lebouteiller

Abstract: We report on the energetics of molecular outflows in 14 local Ultraluminous Infrared Galaxies (ULIRGs) that show unambiguous outflow signatures (P-Cygni profiles or high-velocity absorption wings) in the far-infrared lines of OH measured with the Herschel/PACS spectrometer. Detection of both ground-state (at 119 and 79 um) and one or more radiatively-excited (at 65 and 84 um) lines allows us to mo… ▽ More We report on the energetics of molecular outflows in 14 local Ultraluminous Infrared Galaxies (ULIRGs) that show unambiguous outflow signatures (P-Cygni profiles or high-velocity absorption wings) in the far-infrared lines of OH measured with the Herschel/PACS spectrometer. Detection of both ground-state (at 119 and 79 um) and one or more radiatively-excited (at 65 and 84 um) lines allows us to model the nuclear gas (<~300 pc) as well as the more extended components using spherically symmetric radiative transfer models. The highest molecular outflow velocities are found in buried sources, in which slower but massive expansion of the nuclear gas is also observed. With the exception of a few outliers, the outflows have momentum fluxes of (2-5)xL_IR/c and mechanical luminosities of (0.1-0.3)% of L_IR. The moderate momentum boosts in these sources (<~3) suggest that the outflows are mostly momentum-driven by the combined effects of AGN and nuclear starbursts, as a result of radiation pressure, winds, and supernovae remnants. In some sources (~20%), however, powerful (10^{10.5-11} Lsun) AGN feedback and (partially) energy-conserving phases are required, with momentum boosts in the range 3-20. These outflows appear to be stochastic strong-AGN feedback events that occur throughout the merging process. In a few sources, the outflow activity in the innermost regions has subsided in the last ~1 Myr. While OH traces the molecular outflows at sub-kpc scales, comparison of the masses traced by OH with those previously inferred from tracers of more extended outflowing gas suggests that most mass is loaded (with loading factors of Mdot/SFR=1-10) from the central galactic cores (a few x 100 pc). Outflow depletion timescales are <10^8 yr, shorter than the gas consumption timescales by factors of 1.1-15, and are anti-correlated with the AGN luminosity. △ Less

Submitted 11 February, 2017; v1 submitted 24 December, 2016; originally announced December 2016.

Journal ref: The Astrophysical Journal, 2017, 836, 11

The Dense Molecular Gas and Nuclear Activity in the ULIRG IRAS 13120-5453

Authors: G. C. Privon, S. Aalto, N. Falstad, S. Muller, E. González-Alfonso, K. Sliwa, E. Treister, F. Costagliola, L. Armus, A. S. Evans, S. Garcia-Burillo, T. Izumi, K. Sakamoto, P. van der Werf, J. K. Chu

Abstract: We present new ALMA Band 7 ($\sim340$ GHz) observations of the dense gas tracers HCN, HCO$^+$, and CS in the local, single-nucleus, ultraluminous infrared galaxy IRAS 13120-5453. We find centrally enhanced HCN (4-3) emission, relative to HCO$^+$ (4-3), but do not find evidence for radiative pumping of HCN. Considering the size of the starburst (0.5 kpc) and the estimated supernovae rate of… ▽ More We present new ALMA Band 7 ($\sim340$ GHz) observations of the dense gas tracers HCN, HCO$^+$, and CS in the local, single-nucleus, ultraluminous infrared galaxy IRAS 13120-5453. We find centrally enhanced HCN (4-3) emission, relative to HCO$^+$ (4-3), but do not find evidence for radiative pumping of HCN. Considering the size of the starburst (0.5 kpc) and the estimated supernovae rate of $\sim1.2$ yr$^{-1}$, the high HCN/HCO$^+$ ratio can be explained by an enhanced HCN abundance as a result of mechanical heating by the supernovae, though the active galactic nucleus and winds may also contribute additional mechanical heating. The starburst size implies a high $Σ_{IR}$ of $4.7\times10^{12}$ $L_{\odot}$ kpc$^{-2}$, slightly below predictions of radiation-pressure limited starbursts. The HCN line profile has low-level wings, which we tentatively interpret as evidence for outflowing dense molecular gas. However, the dense molecular outflow seen in the HCN line wings is unlikely to escape the galaxy and is destined to return to the nucleus and fuel future star formation. We also present modeling of Herschel observations of the H$_2$O lines and find a nuclear dust temperature of $\sim40$ K. IRAS 13120-5453 has a lower dust temperature and $Σ_{IR}$ than is inferred for the systems termed "compact obscured nuclei" (such as Arp 220 and Mrk 231). If IRAS 13120-5453 has undergone a compact obscured nucleus phase, we are likely witnessing it at a time when the feedback has already inflated the nuclear ISM and diluted star formation in the starburst/AGN core. △ Less

Submitted 20 December, 2016; v1 submitted 13 December, 2016; originally announced December 2016.

Comments: accepted for publication in ApJ, 21 pages, 11 figures