Showing 1–50 of 79 results for author: Sakamoto, K

A temperature or FUV tracer? The HNC/HCN ratio in M83 on the GMC scale

Authors: Nanase Harada, Toshiki Saito, Yuri Nishimura, Yoshimasa Watanabe, Kazushi Sakamoto

Abstract: The HNC/HCN ratio is observationally known as a thermometer in Galactic interstellar molecular clouds. A recent study has alternatively suggested that the HNC/HCN ratio is affected by the ultraviolet (UV) field, not by the temperature. We aim to study this ratio on the scale of giant molecular clouds in the barred spiral galaxy M83 towards the southwestern bar end and the central region from ALMA… ▽ More The HNC/HCN ratio is observationally known as a thermometer in Galactic interstellar molecular clouds. A recent study has alternatively suggested that the HNC/HCN ratio is affected by the ultraviolet (UV) field, not by the temperature. We aim to study this ratio on the scale of giant molecular clouds in the barred spiral galaxy M83 towards the southwestern bar end and the central region from ALMA observations, and if possible, distinguish the above scenarios. We compare the high (40-50 pc) resolution HNC/HCN ratios with the star formation rate from the 3-mm continuum intensity and the molecular mass inferred from the HCN intensities. Our results show that the HNC/HCN ratios do not vary with the star formation rates, star formation efficiencies, or column densities in the bar-end region. In the central region, the HNC/HCN ratios become higher with higher star formation rates, which tend to cause higher temperatures. This result is not consistent with the previously proposed scenario in which the HNC/HCN ratio decreases with increasing temperature. Spectral shapes suggest that this trend may be due to optically thick HCN and optically thin HNC. In addition, we compare the large-scale ($\sim 200$ pc) correlation between the dust temperature from the FIR ratio and the HNC/HCN ratio for the southwestern bar-end region. The HNC/HCN ratio is lower when the dust temperatures are higher. We suggest from the above results that the HNC/HCN ratio depends on the UV radiation field that affects the interstellar medium on the $\sim100\,$pc scale where the column densities are low. △ Less

Submitted 14 May, 2024; originally announced May 2024.

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

An ALCHEMI inspection of sulphur-bearing species towards the central molecular zone of NGC 253

Authors: M. Bouvier, S. Viti, E. Behrens, J. Butterworth, K. -Y. Huang, J. G. Mangum, N. Harada, S. Martín, V. M. Rivilla, S. Muller, K. Sakamoto, Y. Yoshimura, K. Tanaka, K. Nakanishi, R. Herrero-Illana, L. Colzi, M. D. Gorski, C. Henkel, P. K. Humire, D. S. Meier, P. P. van der Werf, Y. T. Yan

Abstract: Sulphur-bearing species are detected in various environments within Galactic star-forming regions and are particularly abundant in the gas phase of outflows and shocks, and photo-dissociation regions. In this work, we aim to investigate the nature of the emission from the most common sulphur-bearing species observable at millimetre wavelengths towards the nuclear starburst of the galaxy NGC 253. W… ▽ More Sulphur-bearing species are detected in various environments within Galactic star-forming regions and are particularly abundant in the gas phase of outflows and shocks, and photo-dissociation regions. In this work, we aim to investigate the nature of the emission from the most common sulphur-bearing species observable at millimetre wavelengths towards the nuclear starburst of the galaxy NGC 253. We intend to understand which type of regions are probed by sulphur-bearing species and which process(es) dominate(s) the release of sulphur into the gas phase. We used the high-angular resolution (1.6" or 27 pc) observations from the ALCHEMI ALMA Large Program to image several sulphur-bearing species towards the central molecular zone (CMZ) of NGC 253. We performed local thermodynamic equilibrium (LTE) and non-LTE large velocity gradient (LVG) analyses to derive the physical conditions of the gas in which S-bearing species are emitted, and their abundance ratios across the CMZ. Finally, we compared our results with previous ALCHEMI studies and a few selected Galactic environments. We found that not all sulphur-bearing species trace the same type of gas: strong evidence indicates that H2S and part of the emission of OCS, H2CS, and SO, are tracing shocks whilst part of SO and CS emission rather trace the dense molecular gas. For some species, such as CCS and SO2, we could not firmly conclude on their origin of emission. The present analysis indicates that the emission from most sulphur-bearing species throughout the CMZ is likely dominated by shocks associated with ongoing star formation. In the inner part of the CMZ where the presence of super star clusters was previously indicated, we could not distinguish between shocks or thermal evaporation as the main process releasing the S-bearing species. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 44 pages, 20 figures, Accepted for publication in A&A

The rotation of classical bulges in barred galaxies in the presence of gas

Authors: Rubens E. G. Machado, Kenzo R. Sakamoto, Andressa Wille, Gustavo F. Gonçalves

Abstract: Barred galaxies often develop a box/peanut pseudobulge, but they can also host a nearly spherical classical bulge, which is known to gain rotation due to the bar. We aim to explore how the presence of gas impacts the rotation of classical bulges. We carried out a comprehensive set of hydrodynamical N-body simulations with different combinations of bulge masses and gas fractions. In these models, b… ▽ More Barred galaxies often develop a box/peanut pseudobulge, but they can also host a nearly spherical classical bulge, which is known to gain rotation due to the bar. We aim to explore how the presence of gas impacts the rotation of classical bulges. We carried out a comprehensive set of hydrodynamical N-body simulations with different combinations of bulge masses and gas fractions. In these models, both massive bulges and high gas content tend to inhibit the formation of strong bars. For low-mass bulges, the resulting bar is stronger in cases of low gas content. In the stronger bar models, bulges acquire more angular momentum and thus display considerable rotational velocity. Such bulges also develop anisotropic velocity dispersions and become triaxial in shape. We found that the rotation of the bulge becomes less pronounced as the gas fraction is increased from 0 to 30%. These results indicate that the gas content has a significant effect on the dynamics of the classical bulge, because it influences bar strength. Particularly in the case of the low-mass bulges (10% bulge mass fraction), all of the measured rotational and structural properties of the classical bulge depend strongly and systematically on the gas content of the galaxy. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 25 pages, 22 figures, accepted for publication in Universe

Microscale Hydrogen, Carbon, and Nitrogen Isotopic Diversity of Organic Matter in Asteroid Ryugu

Authors: Larry R Nittler, Jens Barosch, Katherine Burgess, Rhonda M Stroud, Jianhua Wang, Hikaru Yabuta, Yuma Enokido, Megumi Matsumoto, Tomoki Nakamura, Yoko Kebukawa, Shohei Yamashita, Yoshio Takahashi, Laure Bejach, Lydie Bonal, George D Cody, Emmanuel Dartois, Alexandre Dazzi, Bradley De Gregorio, Ariane Deniset-Besseau, Jean Duprat, Cécile Engrand, Minako Hashiguchi, A. L. David Kilcoyne, Mutsumi Komatsu, Zita Martins , et al. (35 additional authors not shown)

Abstract: We report the H, C, and N isotopic compositions of microscale (0.2 to 2$μ$m) organic matter in samples of asteroid Ryugu and the Orgueil CI carbonaceous chondrite. Three regolith particles of asteroid Ryugu, returned by the Hayabusa2 spacecraft, and several fragments of Orgueil were analyzed by NanoSIMS isotopic imaging. The isotopic distributions of the Ryugu samples from two different collection… ▽ More We report the H, C, and N isotopic compositions of microscale (0.2 to 2$μ$m) organic matter in samples of asteroid Ryugu and the Orgueil CI carbonaceous chondrite. Three regolith particles of asteroid Ryugu, returned by the Hayabusa2 spacecraft, and several fragments of Orgueil were analyzed by NanoSIMS isotopic imaging. The isotopic distributions of the Ryugu samples from two different collection spots are closely similar to each other and to the Orgueil samples, strengthening the proposed Ryugu-CI chondrite connection. Most individual sub-$μ$m organic grains have isotopic compositions within error of bulk values, but 2-8% of them are outliers exhibiting large isotopic enrichments or depletions in D, $^{15}$N, and/or $^{13}$C. The H, C and N isotopic compositions of the outliers are not correlated with each other: while some C-rich grains are both D- and $^{15}$N-enriched, many are enriched or depleted in one or the other system. This most likely points to a diversity in isotopic fractionation pathways and thus diversity in the local formation environments for the individual outlier grains. The observation of a relatively small population of isotopic outlier grains can be explained either by escape from nebular and/or parent body homogenization of carbonaceous precursor material or addition of later isotopic outlier grains. The strong chemical similarity of isotopically typical and isotopically outlying grains, as reflected by synchrotron x-ray absorption spectra, suggests a genetic connection and thus favors the former, homogenization scenario. However, the fact that even the least altered meteorites show the same pattern of a small population of outliers on top of a larger population of homogenized grains indicates that some or most of the homogenization occurred prior to accretion of the macromolecular organic grains into asteroidal parent bodies. △ Less

Submitted 12 April, 2024; originally announced April 2024.

Comments: Accepted for publication in Earth and Planetary Science Letters. 8 figures (plus 3 supplementary figs), two tables

Physical Properties of the Southwest Outflow Streamer in the Starburst Galaxy NGC 253 with ALCHEMI

Authors: Min Bao, Nanase Harada, Kotaro Kohno, Yuki Yoshimura, Fumi Egusa, Yuri Nishimura, Kunihiko Tanaka, Kouichiro Nakanishi, Sergio Martín, Jeffrey G. Mangum, Kazushi Sakamoto, Sébastien Muller, Mathilde Bouvier, Laura Colzi, Kimberly L. Emig, David S. Meier, Christian Henkel, Pedro Humire, Ko-Yun Huang, Víctor M. Rivilla, Paul van der Werf, Serena Viti

Abstract: The physical properties of galactic molecular outflows are important as they could constrain outflow formation mechanisms. We study the properties of the southwest (SW) outflow streamer including gas kinematics, optical depth, dense gas fraction, and shock strength in the central molecular zone of the starburst galaxy NGC 253. We image the molecular emission at a spatial resolution of $\sim$27 pc… ▽ More The physical properties of galactic molecular outflows are important as they could constrain outflow formation mechanisms. We study the properties of the southwest (SW) outflow streamer including gas kinematics, optical depth, dense gas fraction, and shock strength in the central molecular zone of the starburst galaxy NGC 253. We image the molecular emission at a spatial resolution of $\sim$27 pc based on data from the ALCHEMI program. We trace the kinematics of molecular gas with CO(1-0) line. We constrain the optical depth of CO emission with CO/$^{13}$CO(1-0) ratio, the dense gas fraction with HCN/CO(1-0) ratio, as well as the shock strength with SiO(2-1)/$^{13}$CO(1-0) ratio. The CO/$^{13}$CO(1-0) integrated intensity ratio is $\sim$21 in the SW streamer region, which approximates the C/$^{13}$C isotopic abundance ratio. The higher integrated intensity ratio compared to the disk can be attributed to the optically thinner environment for CO(1-0) emission inside the SW streamer. The HCN/CO(1-0) and SiO(2-1)/$^{13}$CO(1-0) integrated intensity ratios both approach $\sim$0.2 in three giant molecular clouds (GMCs) at the base of the outflow streamers, which implies the higher dense gas fraction and enhanced strength of fast shocks in those GMCs than in the disk. The contours of those two integrated intensity ratios are extended towards the directions of outflow streamers, which connects the enhanced dense gas fraction and shock strength with molecular outflow. Moreover, the molecular gas with enhanced dense gas fraction and shock strength located at the base of the SW streamer shares the same velocity with the outflow. These phenomena suggest that the star formation inside the GMCs can trigger the shocks and further drive the molecular outflow. △ Less

Submitted 6 April, 2024; originally announced April 2024.

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

Molecular isotopologue measurements toward super star clusters and the relation to their ages in NGC253 with ALCHEMI

Authors: J. Butterworth, S. Viti, P. P. Van der Werf, J. G. Mangum, S. Martín, N. Harada, K. L. Emig, S. Muller, K. Sakamoto, Y. Yoshimura, K. Tanaka, R. Herrero-Illana, L. Colzi, V. M. Rivilla, K. Y. Huang, M. Bouvier, E. Behrens, C. Henkel, Y. T. Yan, D. S. Meier, D. Zhou

Abstract: Determining the evolution of the CNO isotopes in the interstellar medium (ISM) of starburst galaxies can yield important constraints on the ages of superstar clusters (SSCs), or on other aspects and contributing factors of their evolution. Due to the time-dependent nature of the abundances of isotopes within the ISM as they are supplied from processes such as nucleosynthesis or chemical fractionat… ▽ More Determining the evolution of the CNO isotopes in the interstellar medium (ISM) of starburst galaxies can yield important constraints on the ages of superstar clusters (SSCs), or on other aspects and contributing factors of their evolution. Due to the time-dependent nature of the abundances of isotopes within the ISM as they are supplied from processes such as nucleosynthesis or chemical fractionation, this provides the possible opportunity to probe the ability of isotopes ratios to trace the ages of high star forming regions, such as SSCs. The goal of this study is to investigate whether the isotopic variations in SSC regions within NGC253 are correlated with their different ages as derived from stellar population modelling. We have measured abundance ratios of CO, HCN and HCO$^+$ isotopologues in six regions containing SSCs within NGC253 using high spatial resolution (1.6",$\sim 28$pc) data from the ALCHEMI (ALma Comprehensive High-resolution Extragalactic Molecular Inventory) ALMA Large program. We have then analysed these ratios using RADEX radiative transfer modelling, with the parameter space sampled using the nested sampling Monte Carlo algorithm MLFriends. These abundance ratios were then compared to ages predicted in each region via the fitting of observed star formation tracers (such as Br$γ$) to starburst stellar population evolution models. We do not find any significant trend with age for the CO and HCN isotopologue ratios on the timescales for the ages of the SSC* regions observed. The driving factors of these ratios within SSCs could be the Initial Mass Function as well as possibly fractionation effects. To further probe these effects in SSCs over time a larger sample of SSCs must be observed spanning a larger age range. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: 44 pages, 43 Figures, Accepted for Publication to A&A

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

The ALCHEMI atlas: principal component analysis reveals starburst evolution in NGC 253

Authors: Nanase Harada, David S. Meier, Sergio Martín, Sebastien Muller, Kazushi Sakamoto, Toshiki Saito, Mark D. Gorski, Christian Henkel, Kunihiko Tanaka, Jeffrey G. Mangum, Susanne Aalto, Rebeca Aladro, Mathilde Bouvier, Laura Colzi, Kimberly L. Emig, Rubén Herrero-Illana, Ko-Yun Huang, Kotaro Kohno, Sabine König, Kouichiro Nakanishi, Yuri Nishimura, Shuro Takano, Víctor M. Rivilla, Serena Viti, Yoshimasa Watanabe , et al. (2 additional authors not shown)

Abstract: Molecular lines are powerful diagnostics of the physical and chemical properties of the interstellar medium (ISM). These ISM properties, which affect future star formation, are expected to differ in starburst galaxies from those of more quiescent galaxies. We investigate the ISM properties in the central molecular zone of the nearby starburst galaxy NGC 253 using the ultra-wide millimeter spectral… ▽ More Molecular lines are powerful diagnostics of the physical and chemical properties of the interstellar medium (ISM). These ISM properties, which affect future star formation, are expected to differ in starburst galaxies from those of more quiescent galaxies. We investigate the ISM properties in the central molecular zone of the nearby starburst galaxy NGC 253 using the ultra-wide millimeter spectral scan survey from the ALMA Large Program ALCHEMI. We present an atlas of velocity-integrated images at a 1".6 resolution of 148 unblended transitions from 44 species, including the first extragalactic detection of HCNH$^+$ and the first interferometric images of C$_3$H$^+$, NO, HCS$^+$. We conduct a principal component analysis (PCA) on these images to extract correlated chemical species and to identify key groups of diagnostic transitions. To the best of our knowledge, our dataset is currently the largest astronomical set of molecular lines to which PCA has been applied. The PCA can categorize transitions coming from different physical components in NGC 253 such as i) young starburst tracers characterized by high-excitation transitions of HC$_3$N and complex organic molecules (COMs) versus tracers of on-going star formation (radio recombination lines) and high-excitation transitions of CCH and CN tracing PDRs, ii) tracers of cloud-collision-induced shocks (low-excitation transitions of CH$_3$OH, HNCO, HOCO$^+$, and OCS) versus shocks from star-formation-induced outflows (high-excitation transitions of SiO), as well as iii) outflows showing emission from HOC$^+$, CCH, H$_3$O$^+$, CO isotopologues, HCN, HCO$^+$, CS, and CN. Our findings show these intensities vary with galactic dynamics, star formation activities, and stellar feedback. △ Less

Submitted 4 January, 2024; originally announced January 2024.

Comments: 65 pages, 39 figures. Accepted for publication in ApJS

Volume density structure of the NGC253 CMZ through ALCHEMI excitation analysis

Authors: Kunihiko Tanaka, Jeffrey G. Mangum, Serena Viti, Sergio Martin, Nanase Harada, Kazushi Sakamoto, Sebastien Muller, Yuki Yoshimura, Kouichiro Nakanishi, Ruben Herrero Illana, Kimberly L. Emig, S. Muhle, Hiroyuki Kaneko, Tomoka Tosaki, Erica Behrens, Victor M. Rivilla, Laura Colzi, Yuri Nishimura, P. K. Humire, Mathilde Bouvier, Ko-Yun Huang, Joshua Butterworth, David S. Meier, Paul P. van der Werf

Abstract: We present a spatially-resolved excitation analysis for the central molecular zone (CMZ) of the starburst galaxy NGC 253 using the data from the ALMA Large program ALCHEMI, whereby we explore parameters distinguishing NGC 253 from the quiescent Milky Way's Galactic Center (GC). Non-LTE analyses employing a hierarchical Bayesian framework are applied to Band 3-7 transitions from nine molecular spec… ▽ More We present a spatially-resolved excitation analysis for the central molecular zone (CMZ) of the starburst galaxy NGC 253 using the data from the ALMA Large program ALCHEMI, whereby we explore parameters distinguishing NGC 253 from the quiescent Milky Way's Galactic Center (GC). Non-LTE analyses employing a hierarchical Bayesian framework are applied to Band 3-7 transitions from nine molecular species to delineate the position-position-velocity distributions of column density ($N_\mathrm{H_2}$), volume density ($n_\mathrm{H_2}$), and temperature ($T_\mathrm{kin}$) at 27 pc resolution. Two distinct components are detected: a low-density component with $(n_\mathrm{H_2},\ T_\mathrm{kin})\sim(10^{3.3}\ \mathrm{cm}^{-3}, 85 K)$ and a high-density component with $(n_\mathrm{H_2},\ T_\mathrm{kin})\sim (10^{4.4}\ \mathrm{cm}^{-3}, 110\ \mathrm{K})$, separated at $n_\mathrm{H_2}\sim10^{3.8}\ \mathrm{cm}^{-3}$. NGC 253 has $\sim10$ times the high-density gas mass and $\sim3$ times the dense-gas mass fraction of the GC. These properties are consistent with their HCN/CO ratio but cannot alone explain the factor of $\sim30$ difference in their star formation efficiencies (SFEs), contradicting the dense-gas mass to star formation rate scaling law. The $n_\mathrm{H_2}$ histogram toward NGC 253 exhibits a shallow declining slope up to $n_\mathrm{H_2}\sim10^6\ \mathrm{cm}^{-3}$, while that of the GC steeply drops in $n_\mathrm{H_2}\gtrsim10^{4.5}\ \mathrm{cm}^{-3}$ and vanishes at $10^5\ \mathrm{cm}^{-3}$. Their dense-gas mass fraction ratio becomes consistent with their SFEs when the threshold $n_\mathrm{H_2}$ for the dense gas is taken at $\sim 10^{4.2\mbox{-}4.6}\ \mathrm{cm}^{-3}$. The rich abundance of gas above this density range in the NGC 253 CMZ, or its scarcity in the GC, is likely to be the critical difference characterizing the contrasting star formation in the centers of the two galaxies. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Comments: 49 pages, 27 figures, 7 tables, accepted for publication in the Astrophysical Journal

Warm Molecular Gas in the Central Parsecs of the Buried Nucleus of NGC 4418 Traced with the Fundamental CO Ro-vibrational Absorptions

Authors: Youichi Ohyama, Shusuke Onishi, Takao Nakagawa, Kosei Matsumoto, Naoki Isobe, Mai Shirahata, Shunsuke Baba, Kazushi Sakamoto

Abstract: We investigated the inner buried nucleus of a nearby luminous infrared galaxy NGC 4418 using high-resolution spectroscopy of fundamental carbon monoxide (CO) ro-vibrational absorptions around $4.67 μ$m for the first time. This method allowed us to examine the physical and kinematical properties in the hot inner region of this nucleus. We detected a series of both very deep (partly saturated)… ▽ More We investigated the inner buried nucleus of a nearby luminous infrared galaxy NGC 4418 using high-resolution spectroscopy of fundamental carbon monoxide (CO) ro-vibrational absorptions around $4.67 μ$m for the first time. This method allowed us to examine the physical and kinematical properties in the hot inner region of this nucleus. We detected a series of both very deep (partly saturated) $^{12}$CO and moderately deep (optically thin) $^{13}$CO absorption lines and inferred a large column density ($N_\mathrm{H2}=(5\pm3)\times10^{23}$ cm$^{-2}$ in front of the $5 μ$m photosphere) of warm ($T_\mathrm{ex}\simeq170$ K) molecular gas by assuming an isothermal plane-parallel slab illuminated by a compact background MIR-emitting source. We modeled that the warm CO absorber almost covers the central heating source and that it is an inner layer around the $5 μ$m photosphere (at $r=$several pc) of a compact shroud of gas and dust ($d\sim100$ pc). The width of the absorption lines ($110$ km s$^{-1}$) and their small deviation from the systemic velocity ($<10$ km s$^{-1}$) are consistent with a warm and turbulent layer with little bulk motion in the radial direction. △ Less

Submitted 31 May, 2023; v1 submitted 17 May, 2023; originally announced May 2023.

Comments: 18 pages, 3 figures, 2 tables, accepted for publication in ApJ

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

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

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

Submitted 22 March, 2023; originally announced March 2023.

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

Diverse Molecular Structures Across The Whole Star-Forming Disk of M83: High fidelity Imaging at 40pc Resolution

Authors: Jin Koda, Akihiko Hirota, Fumi Egusa, Kazushi Sakamoto, Tsuyoshi Sawada, Mark Heyer, Junichi Baba, Samuel Boissier, Daniela Calzetti, Jennifer Donovan Meyer, Bruce G. Elmegreen, Armando Gil de Paz, Nanase Harada, Luis C. Ho, Masato I. N. Kobayashi, Nario Kuno, Amanda M Lee, Barry F. Madore, Fumiya Maeda, Sergio Martin, Kazuyuki Muraoka, Kouichiro Nakanishi, Sachiko Onodera, Jorge L. Pineda, Nick Scoville , et al. (1 additional authors not shown)

Abstract: We present high-fidelity CO(1-0) imaging of molecular gas across the full star-forming disk of M83, using ALMA's 12m, 7m, and TP arrays and the MIRIAD package. The data have a mass sensitivity and resolution of 10^4Msun and 40 pc. The full disk coverage shows that the characteristics of molecular gas change radially from the center to outer disk. The molecular gas distribution shows coherent large… ▽ More We present high-fidelity CO(1-0) imaging of molecular gas across the full star-forming disk of M83, using ALMA's 12m, 7m, and TP arrays and the MIRIAD package. The data have a mass sensitivity and resolution of 10^4Msun and 40 pc. The full disk coverage shows that the characteristics of molecular gas change radially from the center to outer disk. The molecular gas distribution shows coherent large-scale structures in the inner part, including the central concentration, bar offset ridges, and prominent molecular spiral arms. In the outer disk, the spiral arms appear less spatially coherent, and even flocculent. Massive filamentary gas concentrations are abundant even in the interarm regions. Building up these structures in the interarm regions would require a very long time (~>100Myr). Instead, they must have formed within stellar spiral arms and been released into the interarm regions. For such structures to survive through the dynamical processes, the lifetimes of these structures and their constituent molecules and molecular clouds must be long (~>100Myr). These interarm structures host little or no star formation traced by Halpha. The new map also shows extended CO emission, which likely represents an ensemble of unresolved molecular clouds. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

Tracing Interstellar Heating: An ALCHEMI Measurement of the HCN Isomers in NGC 253

Authors: Erica Behrens, Jeffrey G. Mangum, Jonathan Holdship, Serena Viti, Nanase Harada, Sergio Martin, Kazushi Sakamoto, Sebastien Muller, Kunihiko Tanaka, Kouichiro Nakanishi, Ruben Herrero-Illana, Yuki Yoshimura, Rebeca Aladro, Laura Colzi, Kimberly L. Emig, Christian Henkel, Ko-Yun Huang, P. K. Humire, David S. Meier, Victor M. Rivilla

Abstract: We analyze HCN and HNC emission in the nearby starburst galaxy NGC 253 to investigate its effectiveness in tracing heating processes associated with star formation. This study uses multiple HCN and HNC rotational transitions observed using ALMA via the ALCHEMI Large Program. To understand the conditions and associated heating mechanisms within NGC 253's dense gas, we employ Bayesian nested samplin… ▽ More We analyze HCN and HNC emission in the nearby starburst galaxy NGC 253 to investigate its effectiveness in tracing heating processes associated with star formation. This study uses multiple HCN and HNC rotational transitions observed using ALMA via the ALCHEMI Large Program. To understand the conditions and associated heating mechanisms within NGC 253's dense gas, we employ Bayesian nested sampling techniques applied to chemical and radiative transfer models which are constrained using our HCN and HNC measurements. We find that the volume density $n_{\text{H}_{2}}$ and cosmic ray ionization rate (CRIR) $ζ$ are enhanced by about an order of magnitude in the galaxy's central regions as compared to those further from the nucleus. In NGC 253's central GMCs, where observed HCN/HNC abundance ratios are lowest, $n \sim 10^{5.5}$ cm$^{-3}$ and $ζ\sim 10^{-12}$ s$^{-1}$ (greater than $10^4$ times the average Galactic rate). We find a positive correlation in the association of both density and CRIR with the number of star formation-related heating sources (supernova remnants, HII regions, and super hot cores) located in each GMC, as well as a correlation between CRIRs and supernova rates. Additionally, we see an anticorrelation between the HCN/HNC ratio and CRIR, indicating that this ratio will be lower in regions where $ζ$ is higher. Though previous studies suggested HCN and HNC may reveal strong mechanical heating processes in NGC 253's CMZ, we find cosmic ray heating dominates the heating budget, and mechanical heating does not play a significant role in the HCN and HNC chemistry. △ Less

Submitted 8 November, 2022; v1 submitted 13 September, 2022; originally announced September 2022.

Comments: 33 pages, 23 figures, accepted for publication by the Astrophysical Journal

ALCHEMI finds a "shocking" carbon footprint in the starburst galaxy NGC~253

Authors: Nanase Harada, Sergio Martin, Jeff Mangum, Kazushi Sakamoto, Sebastian Muller, Victor Rivilla, Christian Henkel, David Meier, Laura Colzi, Mitsuyoshi Yamagishi, Kunihiko Tanaka, Kouichiro Nakanishi, Ruben Herrero-Illana, Yuki Yoshimura, Pedro Humire, Rebeca Aladro, Paul van der Werf, Kim Emig

Abstract: Centers of starburst galaxies may be characterized by a specific gas and ice chemistry due to their gas dynamics and the presence of various ice desorption mechanisms. This may result in a peculiar observable composition. We analyze abundances of $CO_2$, a reliable tracer of ice chemistry, from data collected as part of the ALMA large program ALCHEMI, a wide-frequency spectral scan toward the star… ▽ More Centers of starburst galaxies may be characterized by a specific gas and ice chemistry due to their gas dynamics and the presence of various ice desorption mechanisms. This may result in a peculiar observable composition. We analyze abundances of $CO_2$, a reliable tracer of ice chemistry, from data collected as part of the ALMA large program ALCHEMI, a wide-frequency spectral scan toward the starburst galaxy NGC~253 with an angular resolution of 1.6$''$. We constrain the $CO_2$ abundances in the gas phase using its protonated form $HOCO^+$. The distribution of $HOCO^+$ is similar to that of methanol, which suggests that $HOCO^+$ is indeed produced from the protonation of $CO_2$ sublimated from ice. The $HOCO^+$ fractional abundances are found to be $(1-2)\times10^{-9}$ at the outer part of the central molecular zone (CMZ), while they are lower ($\sim10^{-10}$) near the kinematic center. This peak fractional abundance at the outer CMZ is comparable to that in the Milky Way CMZ, and orders of magnitude higher than that in Galactic disk star-forming regions. From the range of $HOCO^+/CO_2$ ratios suggested from chemical models, the gas-phase $CO_2$ fractional abundance is estimated to be $(1-20)\times10^{-7}$ at the outer CMZ, and orders of magnitude lower near the center. We estimate the $CO_2$ ice fractional abundances at the outer CMZ to be $(2-5)\times10^{-6}$ from the literature. A comparison between the ice and gas $CO_2$ abundances suggests an efficient sublimation mechanism. This sublimation is attributed to large-scale shocks at the orbital intersections of the bar and CMZ. △ Less

Submitted 30 August, 2022; originally announced August 2022.

Comments: 22 pages, 9 figures. Accepted for publication in the Astrophysical Journal

Presolar stardust in asteroid Ryugu

Authors: Jens Barosch, Larry R. Nittler, Jianhua Wang, Conel M. O'D. Alexander, Bradley T. De Gregorio, Cécile Engrand, Yoko Kebukawa, Kazuhide Nagashima, Rhonda M. Stroud, Hikaru Yabuta, Yoshinari Abe, Jérôme Aléon, Sachiko Amari, Yuri Amelin, Ken-ichi Bajo, Laure Bejach, Martin Bizzarro, Lydie Bonal, Audrey Bouvier, Richard W. Carlson, Marc Chaussidon, Byeon-Gak Choi, George D. Cody, Emmanuel Dartois, Nicolas Dauphas , et al. (99 additional authors not shown)

Abstract: We have conducted a NanoSIMS-based search for presolar material in samples recently returned from C-type asteroid Ryugu as part of JAXA's Hayabusa2 mission. We report the detection of all major presolar grain types with O- and C-anomalous isotopic compositions typically identified in carbonaceous chondrite meteorites: 1 silicate, 1 oxide, 1 O-anomalous supernova grain of ambiguous phase, 38 SiC, a… ▽ More We have conducted a NanoSIMS-based search for presolar material in samples recently returned from C-type asteroid Ryugu as part of JAXA's Hayabusa2 mission. We report the detection of all major presolar grain types with O- and C-anomalous isotopic compositions typically identified in carbonaceous chondrite meteorites: 1 silicate, 1 oxide, 1 O-anomalous supernova grain of ambiguous phase, 38 SiC, and 16 carbonaceous grains. At least two of the carbonaceous grains are presolar graphites, whereas several grains with moderate C isotopic anomalies are probably organics. The presolar silicate was located in a clast with a less altered lithology than the typical extensively aqueously altered Ryugu matrix. The matrix-normalized presolar grain abundances in Ryugu are 4.8$^{+4.7}_{-2.6}$ ppm for O-anomalous grains, 25$^{+6}_{-5}$ ppm for SiC grains and 11$^{+5}_{-3}$ ppm for carbonaceous grains. Ryugu is isotopically and petrologically similar to carbonaceous Ivuna-type (CI) chondrites. To compare the in situ presolar grain abundances of Ryugu with CI chondrites, we also mapped Ivuna and Orgueil samples and found a total of SiC grains and 6 carbonaceous grains. No O-anomalous grains were detected. The matrix-normalized presolar grain abundances in the CI chondrites are similar to those in Ryugu: 23 $^{+7}_{-6}$ ppm SiC and 9.0$^{+5.3}_{-4.6}$ ppm carbonaceous grains. Thus, our results provide further evidence in support of the Ryugu-CI connection. They also reveal intriguing hints of small-scale heterogeneities in the Ryugu samples, such as locally distinct degrees of alteration that allowed the preservation of delicate presolar material. △ Less

Submitted 16 August, 2022; originally announced August 2022.

Comments: 12 pages, 3 figures, 2 tables. Published in ApJL

Journal ref: 2022, The Astrophysical Journal Letters, 935, L3 (12pp)

Methanol masers in NGC 253 with ALCHEMI

Authors: P. K. Humire, C. Henkel, A. Hernández-Gómez, S. Martín, J. Mangum, N. Harada, S. Muller, K. Sakamoto, K. Tanaka, Y. Yoshimura, K. Nakanishi, S. Mühle, R. Herrero-Illana, D. S. Meier, E. Caux, R. Aladro, R. Mauersberger, S. Viti, L. Colzi, V. M. Rivilla, M. Gorski, K. M. Menten, K. -Y. Huang, S. Aalto, P. P. van der Werf , et al. (1 additional authors not shown)

Abstract: Context: Methanol masers of Class I (collisionally-pumped) and Class II (radiatively-pumped) have been studied in great detail in our Galaxy in a variety of astrophysical environments such as shocks and star-forming regions and are helpful to analyze the properties of the dense interstellar medium. However, the study of methanol masers in external galaxies is still in its infancy. Aims: Our main g… ▽ More Context: Methanol masers of Class I (collisionally-pumped) and Class II (radiatively-pumped) have been studied in great detail in our Galaxy in a variety of astrophysical environments such as shocks and star-forming regions and are helpful to analyze the properties of the dense interstellar medium. However, the study of methanol masers in external galaxies is still in its infancy. Aims: Our main goal is to search for methanol masers in the central molecular zone (CMZ; inner 500 pc) of the nearby starburst galaxy NGC 253. Methods: Covering a frequency range between 84 and 373 GHz ($λ$ = 3.6 to 0.8 mm) at high angular (1.6"$\sim$27 pc) and spectral ($\sim$8--9 km s$^{-1}$) resolution with the ALMA large program ALCHEMI, we have probed different regions across the CMZ of NGC 253. In order to look for methanol maser candidates, we employed the rotation diagram method and a set of radiative transfer models. Results: We detect for the first time masers above 84 GHz in NGC 253, covering an ample portion of the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series (at 84, 132, 229, and 278 GHz) and the $J_{0}\rightarrow(J-$ 1)$_{1}-A$ series (at 95, 146, and 198 GHz). This confirms the presence of the Class I maser line at 84 GHz, already reported but now being detected in more than one location. For the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series, we observe a lack of Class I maser candidates in the central star-forming disk. Conclusions: The physical conditions for maser excitation in the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series can be weak shocks and cloud-cloud collisions as suggested by shock tracers (SiO and HNCO) in bi-symmetric shock/active regions located in the outskirts of the CMZ. On the other hand, the presence of photodissociation regions due to a high star-formation rate would be needed to explain the lack of Class I masers in the very central regions. △ Less

Submitted 6 May, 2022; originally announced May 2022.

Comments: Accepted for publication in A&A. 29 pages, 17 figures (4 in Appendix)

Energizing Star Formation: The Cosmic Ray Ionization Rate in NGC 253 Derived From ALCHEMI Measurements of H$_3$O$^+$ and SO

Authors: Jonathan Holdship, Jeffrey G. Mangum, Serena Viti, Erica Behrens, Nanase Harada, Sergio Martín, Kazushi Sakamoto, Sebastien Muller, Kunihiko Tanaka, Kouichiro Nakanishi, Rubén Herrero-Illana, Yuki Yoshimura, Rebeca Aladro, Laura Colzi, Kimberly L. Emig, Christian Henkel, Yuri Nishimura, Víctor M. Rivilla, Paul P. van der Werf

Abstract: The cosmic ray ionization rate (CRIR) is a key parameter in understanding the physical and chemical processes in the interstellar medium. Cosmic rays are a significant source of energy in star formation regions, which impacts the physical and chemical processes which drive the formation of stars. Previous studies of the circum-molecular zone (CMZ) of the starburst galaxy NGC 253 have found evidenc… ▽ More The cosmic ray ionization rate (CRIR) is a key parameter in understanding the physical and chemical processes in the interstellar medium. Cosmic rays are a significant source of energy in star formation regions, which impacts the physical and chemical processes which drive the formation of stars. Previous studies of the circum-molecular zone (CMZ) of the starburst galaxy NGC 253 have found evidence for a high CRIR value; $10^3-10^6$ times the average cosmic ray ionization rate within the Milky Way. This is a broad constraint and one goal of this study is to determine this value with much higher precision. We exploit ALMA observations towards the central molecular zone of NGC 253 to measure the CRIR. We first demonstrate that the abundance ratio of H$_3$O$^+$ and SO is strongly sensitive to the CRIR. We then combine chemical and radiative transfer models with nested sampling to infer the gas properties and CRIR of several star-forming regions in NGC 253 due to emission from their transitions. We find that each of the four regions modelled has a CRIR in the range $(1-80)\times10^{-14}$ s$^{-1}$ and that this result adequately fits the abundances of other species that are believed to be sensitive to cosmic rays including C$_2$H, HCO$^+$, HOC$^+$, and CO. From shock and PDR/XDR models, we further find that neither UV/X-ray driven nor shock dominated chemistry are a viable single alternative as none of these processes can adequately fit the abundances of all of these species. △ Less

Submitted 11 May, 2022; v1 submitted 7 April, 2022; originally announced April 2022.

Comments: 24 pages, 15 figures, accepted for publication in ApJ, updated to include changes made for proofs

First extragalactic detection of a phosphorus-bearing molecule with ALCHEMI: phosphorus nitride (PN)

Authors: D. Haasler, V. M. Rivilla, S. Martín, J. Holdship, S. Viti, N. Harada, J. Mangum, K. Sakamoto, S. Muller, K. Tanaka, Y. Yoshimura, K. Nakanishi, L. Colzi, L. Hunt, K. L. Emig, R. Aladro, P. Humire, C. Henkel, P. van der Werf

Abstract: Phosphorus (P) is a crucial element for life given its central role in several biomolecules. P-bearing molecules have been discovered in different regions of the Milky Way, but not yet towards an extragalactic environment. We have searched for P-bearing molecules towards the nearby starburst Galaxy NGC 253. Using observations from the ALMA Comprehensive High-resolution Extragalactic Molecular Inve… ▽ More Phosphorus (P) is a crucial element for life given its central role in several biomolecules. P-bearing molecules have been discovered in different regions of the Milky Way, but not yet towards an extragalactic environment. We have searched for P-bearing molecules towards the nearby starburst Galaxy NGC 253. Using observations from the ALMA Comprehensive High-resolution Extragalactic Molecular Inventory (ALCHEMI) project, we used the MADCUBA package to model the emission of P-bearing molecules assuming Local Thermodynamic Equilibrium (LTE) conditions. We have also performed a non-LTE analysis using SpectralRadex. We report the detection of a P-bearing molecule, phosphorus nitride (PN), for the first time in an extragalactic environment, towards two giant molecular clouds (GMCs) of NGC 253. The LTE analysis yields total PN beam-averaged column densities $N$=(1.20$\pm$0.09)$\times$10$^{13}$ cm$^{-2}$ and $N$=(6.5$\pm$1.6)$\times$10$^{12}$ cm$^{-2}$, which translate into abundances with respect to H$_2$ of $χ$=(8.0$\pm$1.0)$\times$10$^{-12}$ and $χ$=(4.4$\pm$1.2)$\times$10$^{-12}$. We derived a low excitation temperature of $T_{\rm ex}$=(4.4$\pm$1.3) K towards the GMC with the brightest PN emission, which indicates that PN is sub-thermally excited. The non-LTE analysis results in column densities consistent with the LTE values. We have also searched for other P-bearing molecules (PO, PH$_{3}$, CP and CCP), and upper limits were derived. The derived PO/PN ratios are $<$1.3 and $<$1.7. The abundance ratio between PN and the shock-tracer SiO derived towards NGC 253 follows the same trend previously found towards Galactic sources. Comparison of the observations with chemical models indicates that the derived molecular abundances of PN in NGC 253 can be explained by shock-driven chemistry followed by cosmic-ray-driven photochemistry. △ Less

Submitted 15 January, 2022; v1 submitted 9 December, 2021; originally announced December 2021.

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

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

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

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

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

Submitted 17 September, 2021; originally announced September 2021.

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

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

Deeply Buried Nuclei in the Infrared-Luminous Galaxies NGC 4418 and Arp 220: II. Line Forests at $λ= $1.4--0.4 mm and Circumnuclear Gas Observed with ALMA

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

Abstract: We present the line observations in our ALMA imaging spectral scan toward three deeply buried nuclei in NGC 4418 and Arp 220. We cover 67 GHz in $f_{\rm rest}$=215-697 GHz at about 0.2$"$ (30, 80 pc) resolution. All the nuclei show dense line forests; we report our initial line identification using 55 species. The line velocities generally indicate gas rotation around each nucleus, tracing nuclear… ▽ More We present the line observations in our ALMA imaging spectral scan toward three deeply buried nuclei in NGC 4418 and Arp 220. We cover 67 GHz in $f_{\rm rest}$=215-697 GHz at about 0.2$"$ (30, 80 pc) resolution. All the nuclei show dense line forests; we report our initial line identification using 55 species. The line velocities generally indicate gas rotation around each nucleus, tracing nuclear disks of $\sim$100 pc sizes. We confirmed the counter-rotation of the nuclear disks in Arp 220 and that of the nuclear disk and the galactic disk in NGC 4418. While the brightest lines exceed 100 K, most of the major lines and many $^{13}$C isotopologues show absorption against even brighter continuum cores of the nuclei. The lines with higher upper-level energies, including those from vibrationally-excited molecules, tend to arise from smaller areas, indicating radially varying conditions in these nuclei. The outflows from the two Arp 220 nuclei cause blueshifted line absorption below the continuum level. The absorption mostly has small spatial offsets from the continuum peaks to indicate the outflow orientations. The bipolar outflow from the western nucleus is also imaged in multiple emission lines, showing the extent of $\sim$1$"$ (400 pc). Redshifted line absorption against the nucleus of NGC 4418 indicates either an inward gas motion or a small collimated outflow slanted to the nuclear disk. We also resolved some previous confusions due to line blending and misidentification. △ Less

Submitted 17 September, 2021; originally announced September 2021.

Comments: 56 pages, 27 figures, 5 tables, accepted for publication in ApJ

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

Starburst Energy Feedback Seen Through HCO$^+$/HOC$^+$ Emission in NGC 253 from ALCHEMI

Authors: Nanase Harada, Sergio Martín, Jeffrey G. Mangum, Kazushi Sakamoto, Sebastien Muller, Kunihiko Tanaka, Kouichiro Nakanishi, Rubén Herrero-Illana, Yuki Yoshimura, Stefanie Mühle, Rebeca Aladro, Laura Colzi, Víctor M. Rivilla, Susanne Aalto, Erica Behrens, Christian Henkel, Jonathan Holdship, P. K. Humire, David S. Meier, Yuri Nishimura, Paul P. van der Werf, Serena Viti

Abstract: Molecular abundances are sensitive to UV-photon flux and cosmic-ray ionization rate. In starburst environments, the effects of high-energy photons and particles are expected to be stronger. We examine these astrochemical signatures through multiple transitions of HCO$^+$ and its metastable isomer HOC$^+$ in the center of the starburst galaxy NGC 253 using data from the ALMA large program ALCHEMI.… ▽ More Molecular abundances are sensitive to UV-photon flux and cosmic-ray ionization rate. In starburst environments, the effects of high-energy photons and particles are expected to be stronger. We examine these astrochemical signatures through multiple transitions of HCO$^+$ and its metastable isomer HOC$^+$ in the center of the starburst galaxy NGC 253 using data from the ALMA large program ALCHEMI. The distribution of the HOC$^+$(1-0) integrated intensity shows its association with "superbubbles", cavities created either by supernovae or expanding HII regions. The observed HCO$^+$/HOC$^+$ abundance ratios are $\sim 10-150$, and the fractional abundance of HOC$^+$ relative to H$_2$ is $\sim 1.5\times 10^{-11} - 6\times 10^{-10}$, which implies that the HOC$^+$ abundance in the center of NGC 253 is significantly higher than in quiescent spiral-arm dark clouds in the Galaxy and the Galactic center clouds. Comparison with chemical models implies either an interstellar radiation field of $G_0\gtrsim 10^3$ if the maximum visual extinction is $\gtrsim 5$, or a cosmic-ray ionization rate of $ζ\gtrsim 10^{-14}$ s$^{-1}$ (3-4 orders of magnitude higher than that within clouds in the Galactic spiral-arms) to reproduce the observed results. From the difference in formation routes of HOC$^+$, we propose that a low-excitation line of HOC$^+$ traces cosmic-ray dominated regions, while high-excitation lines trace photodissociation regions. Our results suggest that the interstellar medium in the center of NGC 253 is significantly affected by energy input from UV-photons and cosmic rays, sources of energy feedback. △ Less

Submitted 14 September, 2021; originally announced September 2021.

Comments: 57 pages, 32 figures, accepted for publication in ApJ

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

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

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

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

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

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

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)

ALMA resolves the remarkable molecular jet and rotating wind in the extremely radio-quiet galaxy NGC 1377

Authors: S. Aalto, N. Falstad, S. Muller, K. Wada, J. S. Gallagher, S. König, K. Sakamoto, W. Vlemmings, C. Ceccobello, K. Dasyra, F. Combes, S. García-Burillo, Y. Oya, S. Martín, P. van der Werf, A. S. Evans, J. Kotilainen

Abstract: Submillimetre and millimetre observations are important in probing the properties of the molecular gas and dust around obscured active galactic nuclei (AGNs) and their feedback. With very high-resolution (0."02x0."03 (2x3 pc)) ALMA 345 GHz observations of CO 3-2, HCO$^+$ 4-3, HCN 4-3 $ν_2$=1$f$, and continuum we have studied the molecular outflow and nucleus of the extremely radio-quiet lenticular… ▽ More Submillimetre and millimetre observations are important in probing the properties of the molecular gas and dust around obscured active galactic nuclei (AGNs) and their feedback. With very high-resolution (0."02x0."03 (2x3 pc)) ALMA 345 GHz observations of CO 3-2, HCO$^+$ 4-3, HCN 4-3 $ν_2$=1$f$, and continuum we have studied the molecular outflow and nucleus of the extremely radio-quiet lenticular galaxy NGC1377. The outflow is resolved, revealing a 150 pc long, clumpy, high-velocity, collimated molecular jet. The molecular emission is emerging from the spine of the jet with an average diameter of 3-7 pc. A narrow-angle, rotating molecular wind surrounds the jet and is enveloped by a larger-scale, slower CO-emitting structure. The jet and narrow wind are turbulent ($σ>$40 kms$^{-1}$) and have steep radial gas excitation gradients. The jet shows velocity reversals that we propose are caused by precession, or episodic directional changes. We suggest that an important process powering the outflow is magneto-centrifugal driving. In contrast, the large-scale CO-envelope may be a slow wind, or cocoon that stems from jet-wind interactions. An asymmetric, nuclear r$\sim$2 pc and hot (>180 K) dust structure with a high molecular column density, N(H$_2$)$\sim1.8 \times 10^{24}$ cm$^{-2}$, is detected in continuum and vibrationally excited HCN. Its luminosity is likely powered by a buried AGN. The mass of the supermassive black hole (SMBH) is estimated to $\sim9\times10^6$ M$_\odot$ and the SMBH of NGC1377 appears to be at the end of an intense phase of accretion. The nuclear growth may be fuelled by low-angular momentum gas inflowing from gas ejected in the molecular jet and wind. Such a feedback-loop of cyclic accretion and outflows would be an effective process in growing the nuclear SMBH. This result invites new questions as to SMBH growth processes in obscured, dusty galaxies. △ Less

Submitted 26 July, 2020; v1 submitted 15 July, 2020; originally announced July 2020.

Comments: 16 pages, 12 figures, accepted for publication in Astronomy and Astrophysics. Updated affiliations, added reference

VV 655 and NGC 4418: Implications of an interaction for the evolution of a LIRG

Authors: Erin Boettcher, John S. Gallagher III, Youichi Ohyama, Eskil Varenius, Susanne Aalto, Niklas Falstad, Sabine König, Kazushi Sakamoto, Tova M. Yoast-Hull

Abstract: VV 655, a dwarf irregular galaxy with HI tidal debris, is a companion to the lenticular luminous infrared galaxy (LIRG) NGC 4418. NGC 4418 stands out among nearby LIRGs due to its dense central concentration of molecular gas and the dusty, bi-polar structures along its minor axis suggestive of a wind driven by a central starburst and possible nuclear activity. We seek to understand the consequence… ▽ More VV 655, a dwarf irregular galaxy with HI tidal debris, is a companion to the lenticular luminous infrared galaxy (LIRG) NGC 4418. NGC 4418 stands out among nearby LIRGs due to its dense central concentration of molecular gas and the dusty, bi-polar structures along its minor axis suggestive of a wind driven by a central starburst and possible nuclear activity. We seek to understand the consequences of the ongoing minor interaction between VV 655 and NGC 4418 for the evolution of the LIRG, including the origin of the gas supply responsible for its unusual nuclear properties. We investigate the structural, kinematic, and chemical properties of VV 655 and NGC 4418 by analyzing archival imaging data and optical spectroscopic observations from the SDSS-III and new spectra from SALT-RSS. We characterize their gas-phase metal abundances and spatially resolved, ionized gas kinematics, and demonstrate that the gas-phase metallicity in NGC 4418 significantly exceeds that in VV 655. No kinematic disturbances in the ionized gas are observed along the minor axis of NGC 4418, but we see evidence for ionized gas outflows from VV 655 that may increase the cross-section for gas stripping in grazing collisions. A faint, asymmetric outer arm is detected in NGC 4418 of the type normally associated with galaxy-galaxy interactions. The simplest model suggests that the minor interaction between VV 655 and NGC 4418 produced the unusual nuclear properties of the LIRG via tidal torquing of the interstellar medium of NGC 4418 rather than through a significant gas transfer event. In addition to inducing a central concentration of gas in NGC 4418, this interaction also produced an enhanced star formation rate and an outer tidal arm in the LIRG. The VV 655-NGC 4418 system offers an example of the potential for minor collisions to alter the evolutionary pathways of giant galaxies. △ Less

Submitted 20 April, 2020; originally announced April 2020.

Comments: 9 pages, 4 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 637, A17 (2020)

Systematic Variations of CO J=2-1/1-0 Ratio in The Barred Spiral Galaxy M83

Authors: J. Koda, T. Sawada, K. Sakamoto, A. Hirota, F. Egusa, S. Boissier, D. Calzetti, J. Donovan Meyer, B. G. Elmegreen, A. Gil de Paz, N. Harada, L. C. Ho, M. I. N. Kobayashi, N. Kuno, S. Martin, K. Muraoka, K. Nakanishi, N. Scoville, M. Seibert, C. Vlahakis, Y. Watanabe

Abstract: We present spatial variations of the CO J=2-1/1-0 line ratio in M83 using Total Power array data from ALMA. While the intensities of these two lines correlate tightly, the ratio varies over the disk, with a disk average ratio of 0.69, and shows the galactic center and a two-arm spiral pattern. It is high (>0.7) in regions of high molecular gas surface density, but ranges from low to high ratios in… ▽ More We present spatial variations of the CO J=2-1/1-0 line ratio in M83 using Total Power array data from ALMA. While the intensities of these two lines correlate tightly, the ratio varies over the disk, with a disk average ratio of 0.69, and shows the galactic center and a two-arm spiral pattern. It is high (>0.7) in regions of high molecular gas surface density, but ranges from low to high ratios in regions of low surface density. The ratio correlates well with the spatial distributions and intensities of FUV and IR emissions, with FUV being the best correlated. It also correlates better with the ratio of IR intensities (70/350mic), a proxy for dust temperature, than with the IR intensities. Taken together, these results suggest either a direct or indirect link between the dust heating by the interstellar radiation field and the condition of GMCs, even though no efficient mechanism is known for a thermal coupling of dust and bulk gas in GMCs. △ Less

Submitted 29 January, 2020; originally announced January 2020.

Comments: Accepted for ApJ Letters

Chemical evolution along the circumnuclear ring of M83

Authors: Nanase Harada, Kazushi Sakamoto, Sergio Martin, Yoshimasa Watanabe, Rebeca Aladro, Denise Riquelme, Akihiko Hirota

Abstract: We report an astrochemical study on the evolution of interstellar molecular clouds and consequent star formation in the center of the barred spiral galaxy M83. We used the Atacama Large Millimeter/submillimeter Array (ALMA) to image molecular species indicative of shocks (SiO, CH$_3$OH), dense cores (N$_2$H$^+$), and photodissociation regions (CN and CCH), as well as a radio recombination line (H4… ▽ More We report an astrochemical study on the evolution of interstellar molecular clouds and consequent star formation in the center of the barred spiral galaxy M83. We used the Atacama Large Millimeter/submillimeter Array (ALMA) to image molecular species indicative of shocks (SiO, CH$_3$OH), dense cores (N$_2$H$^+$), and photodissociation regions (CN and CCH), as well as a radio recombination line (H41$α$) tracing active star-forming regions. M83 has a circumnuclear gas ring that is joined at two areas by gas streams from the leading-edge gas lanes on the bar. We found elevated abundances of the shock and dense-core tracers in one of the orbit-intersecting areas, and found peaks of CN and H41$α$ downstream. At the other orbit-intersection area, we found similar enhancement of the shock tracers, but less variation of other tracers, and no sign of active star formation in the stream. We propose that the observed chemical variation or lack of it is due to the presence or absence of collision-induced evolution of molecular clouds and induced star formation. This work presents the most clear case of the chemical evolution in the circumnuclear rings of barred galaxies, thanks to the ALMA resolution and sensitivity. △ Less

Submitted 4 September, 2019; originally announced September 2019.

Comments: 31 pages, 12 figures, accepted for publication in The Astrophysical Journal

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

The hidden heart of the luminous infrared galaxy IC 860 I. A molecular inflow feeding opaque, extreme nuclear activity

Authors: S. Aalto, S. Muller, S. König, N. Falstad, J. Mangum, K. Sakamoto, G. C. Privon, J. Gallagher, F. Combes, S. García-Burillo, S. Martín, S. Viti, P. van der Werf, A. S. Evans, J. H. Black, E. Varenius, R. Beswick, G. Fuller, C. Henkel, K. Kohno, K. Alatalo, S. Mühle

Abstract: High-resolution (0."03 - 0."09) (9 - 26 pc)) ALMA (100 - 350 GHz (3 - 0.8 mm)) and VLA 45 GHz measurements have been used to image continuum and spectral line emission from the inner region of the nearby infrared luminous galaxy IC 860. We detect compact (r=8-10 pc), luminous, mm continuum emission in the core of IC 860, with brightness temperatures $T_{\rm B}>$ 280 K. The 45 GHz continuum is equa… ▽ More High-resolution (0."03 - 0."09) (9 - 26 pc)) ALMA (100 - 350 GHz (3 - 0.8 mm)) and VLA 45 GHz measurements have been used to image continuum and spectral line emission from the inner region of the nearby infrared luminous galaxy IC 860. We detect compact (r=8-10 pc), luminous, mm continuum emission in the core of IC 860, with brightness temperatures $T_{\rm B}>$ 280 K. The 45 GHz continuum is equally compact, but fainter in flux. We suggest that the 3 to 0.8 mm continuum is opaque at mm-wavelengths, implying a very large H$_2$ column density $N$(H$_2$) of 1e26 cm-2, and that it emerges from hot dust with temperature $T_{\rm d} \sim 280$ K. Vibrationally excited lines of HCN $ν_2$=1f J=4-3 and 3-2 (HCN-VIB) are seen in emission and resolved on scales of 40-50 pc. The HCN-VIB emission reveals a north-south nuclear velocity gradient with projected rotation velocities of v=100 kms$^{-1}$ at r=10 pc. The brightest HCN-VIB emission is oriented perpendicular to the velocity gradient. Ground state lines of HCN 3-2, 4-3, HC$^{15}$N 4-3, HCO$^+$ 3-2, 4-3 and CS 7-6, show complex line absorption and emission features. HCN and HCO$^+$ have red-shifted, reversed P-Cygni profiles consistent with gas inflows of $\sim$50 kms$^{-1}$. The absorption features can be traced from the north-east into the nucleus. In contrast, CS 7-6 shows blue-shifted line wings extending to -180 kms$^{-1}$. We suggest that a dense and slow outflow is hidden behind a foreground layer of inflowing gas. It appears that the centre of IC 860 is in a phase of rapid evolution where an inflow is building up the nuclear column density of gas. A slow, dense outflow may be signalling the onset of feedback. The inner, r=10 pc, IR luminosity can be powered by an accreting black hole and/or a compact starburst with a top-heavy initial mass function. △ Less

Submitted 4 June, 2019; v1 submitted 17 May, 2019; originally announced May 2019.

Comments: 13 pages, 9 figures, Accepted for publication in Astronomy and Astrophysics; Language edited

Spatially resolved carbon and oxygen isotopic ratios in NGC 253 using optically thin tracers

Authors: S. Martin, S. Muller, C. Henkel, D. S. Meier, R. Aladro, K. Sakamoto, P. P. van der Werf

Abstract: To trace chemical evolution in a prototypical starburst environment, we spatially resolve the carbon and oxygen isotope ratios across the central molecular zone (full size ~$\sim 600$ pc) in the nearby starburst galaxy NGC 253. We imaged the emission of the optically thin isotopologues $^{13}$CO, C$^{18}$O, C$^{17}$O, $^{13}$C$^{18}$O at a spatial resolution $\sim50$ pc, comparable to the typical… ▽ More To trace chemical evolution in a prototypical starburst environment, we spatially resolve the carbon and oxygen isotope ratios across the central molecular zone (full size ~$\sim 600$ pc) in the nearby starburst galaxy NGC 253. We imaged the emission of the optically thin isotopologues $^{13}$CO, C$^{18}$O, C$^{17}$O, $^{13}$C$^{18}$O at a spatial resolution $\sim50$ pc, comparable to the typical size of giant molecular associations. Optical depth effects and contamination of $^{13}$C$^{18}$O by C$_4$H is discussed and accounted for to derive column densities. This is the first extragalactic detection of the double isotopologue $^{13}$C$^{18}$O. Derived isotopic ratios $^{12}$C/$^{13}$C$\sim21\pm6$, $^{16}$O/$^{18}$O$\sim130\pm40$, and $^{18}$O/$^{17}$O$\sim4.5\pm0.8$ differ from the generally adopted values in the nuclei of galaxies. The molecular clouds in the central region of NGC 253 show similar rare isotope enrichment to those within the central molecular zone of the Milky way. This enrichment is attributed to stellar nucleosynthesis. Measured isotopic ratios suggest an enhancement of $^{18}$O as compared to our Galactic center, which we attribute to an extra $^{18}$O injection from massive stars. Our observations show evidence for mixing of distinct gas components with different degrees of processing. We observe an extra molecular component of highly processed gas on top of the already proposed less processed gas being transported to the central region of NGC 253. Such multicomponent nature and optical depth effects may hinder the use of isotopic ratios based on spatially unresolved line to infer the star formation history and/or initial stellar mass function properties in the nuclei of galaxies. △ Less

Submitted 7 March, 2019; originally announced March 2019.

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

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

Dusty Superwind from a Galaxy with a Compact Obscured Nucleus: Optical Spectroscopic Study of NGC 4418

Authors: Youichi Ohyama, Kazushi Sakamoto, Susanne Aalto, John S. Gallagher III

Abstract: We report our optical spectroscopic study of the nucleus and its surrounding region of a nearby luminous infrared galaxy NGC 4418. This galaxy has been known to host a compact obscured nucleus, showing distinct characteristics such as a very compact ($\sim 20$ pc) sub-mm and mid-infrared core and dusty circumnuclear region with massive molecular gas concentration. We detected dusty superwind outfl… ▽ More We report our optical spectroscopic study of the nucleus and its surrounding region of a nearby luminous infrared galaxy NGC 4418. This galaxy has been known to host a compact obscured nucleus, showing distinct characteristics such as a very compact ($\sim 20$ pc) sub-mm and mid-infrared core and dusty circumnuclear region with massive molecular gas concentration. We detected dusty superwind outflow at $\gtrsim 1$ kpc scale along the disk semiminor axis in both shock-heated emission lines and enhanced interstellar Na D absorption. This superwind shows basic characteristics similar to those of the prototypical superwind in the starburst galaxy M82, such as a kpc-scale extended structure of gas and dust along the disk minor axis, outflowing components (multiphase gas and dust), physical conditions of the ionized gas, and monotonically blueshifting radial velocity field with increasing distance from the nucleus on the front side of the superwind. We also detected a moderately extinct starburst population in the SDSS nuclear spectrum with the burst age of $\simeq 10$ Myr and stellar mass of $\simeq 1\times 10^7\ M_\mathrm{\odot}$. It is powerful enough to drive the superwind within the dynamical age of the superwind ($\simeq 10$ Myr). On the basis of comparison between this starburst--superwind scenario and the observations in terms of the burst age, stellar mass, infrared luminosity, and obscuration in the optical bands, we argue that this superwind-driving starburst is separate from the sub-mm core even if the core is a very young star cluster. Therefore, this galaxy hosts both the enshrouded compact core and the superwind-driving circumnuclear starburst. △ Less

Submitted 21 December, 2018; originally announced December 2018.

Comments: 44 pages, 12 figures, accepted for publication in ApJ

ALMA observations of molecular tori around massive black holes

Authors: F. Combes, S. Garcia-Burillo, A. Audibert, L. Hunt, A. Eckart, S. Aalto, V. Casasola, F. Boone, M. Krips, S. Viti, K. Sakamoto, S. Muller, K. Dasyra, P. van der Werf, S. Martin

Abstract: We report ALMA observations of CO(3-2) emission in a sample of 7 Seyfert/LINER galaxies at the unprecedented spatial resolution of 0.1'' = 4-9 pc. Our aim is to explore the close environment of AGN, and the dynamical structures leading to its fueling. The selected galaxies host low-luminosity AGN, and have a wide range of activity types, and barred or ringed morphologies. The observed maps reveal… ▽ More We report ALMA observations of CO(3-2) emission in a sample of 7 Seyfert/LINER galaxies at the unprecedented spatial resolution of 0.1'' = 4-9 pc. Our aim is to explore the close environment of AGN, and the dynamical structures leading to its fueling. The selected galaxies host low-luminosity AGN, and have a wide range of activity types, and barred or ringed morphologies. The observed maps reveal the existence of circum-nuclear disk structures, defined by their morphology and decoupled kinematics. We call these structures "molecular tori", even though they appear often as disks, without holes in the center. They have varied orientations, unaligned with their host galaxy. The radius of the tori ranges from 6 to 27 pc, and their mass from 0.7 10$^7$ to 3.9 10$^7$ Msun. At larger scale, the gas is always piled up in a few 100~pc scale resonant rings, that play the role of a reservoir to fuel the nucleus. In some cases, a trailing spiral is observed inside the ring, providing evidence for feeding processes. Most frequently, the torus and the AGN are slightly off-centered, with respect to the bar-resonant ring position, implying that the black hole is wandering by a few 10~pc amplitude around the center of mass of the galaxy. Our spatial resolution allows us to measure gas velocities inside the sphere of influence of the central black holes. By fitting the observations with different simulated cubes, varying the torus inclination and the black hole mass, it is possible to estimate the mass of the central black hole, which is in general difficult for such late-type galaxies, with only a pseudo-bulge. △ Less

Submitted 12 February, 2019; v1 submitted 2 November, 2018; originally announced November 2018.

Comments: 20 pages, 21 figures, accepted in Astron. and Astrophys

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

ALMA Astrochemical Observations of the Infrared-Luminous Merger NGC 3256

Authors: Nanase Harada, Kazushi Sakamoto, Sergio Martin, Susanne Aalto, Rebeca Aladro, Kazimierz Sliwa

Abstract: In external galaxies, molecular composition may be influenced by extreme environments such as starbursts and galaxy mergers. To study such molecular chemistry, we observed the luminous-infrared galaxy and merger NGC 3256 using the Atacama Large Millimeter/sub-millimeter Array. We covered most of the 3-mm and 1.3-mm bands for a multi-species, multi-transition analysis. We first analyzed intensity r… ▽ More In external galaxies, molecular composition may be influenced by extreme environments such as starbursts and galaxy mergers. To study such molecular chemistry, we observed the luminous-infrared galaxy and merger NGC 3256 using the Atacama Large Millimeter/sub-millimeter Array. We covered most of the 3-mm and 1.3-mm bands for a multi-species, multi-transition analysis. We first analyzed intensity ratio maps of selected lines such as HCN/HCO$^+$, which shows no enhancement at an AGN. We then compared the chemical compositions within NGC 3256 at the two nuclei, tidal arms, and positions with influence from galactic outflows. We found the largest variation in SiO and CH$_3$OH, species that are likely to be enhanced by shocks. Next, we compared the chemical compositions in the nuclei of NGC 3256, NGC 253, and Arp 220; these galactic nuclei have varying star formation efficiencies. Arp 220 shows higher abundances of SiO and HC$_3$N than NGC 3256 and NGC 253. Abundances of most species do not show strong correlation with the star formation efficiencies, although the CH$_3$CCH abundance seems to have a weak positive correlation with the star formation efficiency. Lastly, the chemistry of spiral arm positions in NGC 3256 is compared with that of W 51, a Galactic molecular cloud complex in a spiral arm. We found higher fractional abundances of shock tracers, and possibly also higher dense gas fraction in NGC 3256 compared with W 51. △ Less

Submitted 18 January, 2018; originally announced January 2018.

Comments: 54 pages, 19 figures. Accepted for publication in ApJ. Some extensive figures will only be shown online-only in the published version

Fast, Collimated Outflow in the Western Nucleus of Arp 220

Authors: Loreto Barcos-Muñoz, Susanne Aalto, Todd A. Thompson, Kazushi Sakamoto, Sergio Martín, Adam K. Leroy, George C. Privon, Aaron S. Evans, Amanda Kepley

Abstract: We present the first spatially and spectrally resolved image of the molecular outflow in the western nucleus of Arp\,220. The outflow, seen in HCN~(1--0) by ALMA, is compact and collimated, with extension $\lesssim$ 120\,pc. Bipolar morphology emerges along the minor axis of the disk, with redshifted and blueshifted components reaching maximum inclination-corrected velocity of $\sim\,\pm$\,840\,km… ▽ More We present the first spatially and spectrally resolved image of the molecular outflow in the western nucleus of Arp\,220. The outflow, seen in HCN~(1--0) by ALMA, is compact and collimated, with extension $\lesssim$ 120\,pc. Bipolar morphology emerges along the minor axis of the disk, with redshifted and blueshifted components reaching maximum inclination-corrected velocity of $\sim\,\pm$\,840\,km\,s$^{-1}$. The outflow is also seen in CO and continuum emission, the latter implying that it carries significant dust. We estimate a total mass in the outflow of $\geqslant$\,10$^{6}$\,M$_{\odot}$, a dynamical time of $\sim$\,10$^{5}$\,yr, and mass outflow rates of $\geqslant55$\,M$_{\odot}$\,yr$^{-1}$ and $\geqslant\,15$\,M$_{\odot}$\,yr$^{-1}$ for the northern and southern lobes, respectively. Possible driving mechanisms include supernovae energy and momentum transfer, radiation pressure feedback, and a central AGN. The latter could explain the collimated morphology of the HCN outflow, however we need more complex theoretical models, including contribution from supernovae and AGN, to pinpoint the driving mechanism of this outflow. △ Less

Submitted 18 December, 2017; originally announced December 2017.

Comments: Accepted for publication in ApJL

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)

Resolved Structure of Arp 220 Nuclei at λ~3 mm

Authors: Kazushi Sakamoto, Susanne Aalto, Loreto Barcos-Muñoz, Francesco Costagliola, Aaron S. Evans, Nanase Harada, Sergio Martín, Martina Wiedner, David Wilner

Abstract: We analyze 3 mm emission of the ultraluminous infrared galaxy Arp 220 for spatially-resolved structure and spectral properties of the merger nuclei. ALMA archival data at ~0.05" resolution are used for extensive visibility fitting and deep imaging of continuum emission. The data are fitted well with two concentric components for each nucleus, such as two Gaussians or one Gaussian plus one exponent… ▽ More We analyze 3 mm emission of the ultraluminous infrared galaxy Arp 220 for spatially-resolved structure and spectral properties of the merger nuclei. ALMA archival data at ~0.05" resolution are used for extensive visibility fitting and deep imaging of continuum emission. The data are fitted well with two concentric components for each nucleus, such as two Gaussians or one Gaussian plus one exponential disk. The larger components in individual nuclei are similar in shape and extent, ~100-150 pc, to the cm-wave emission due to supernovae. They are therefore identified with the known starburst nuclear disks. The smaller components in both nuclei have about a few 10 pc sizes and peak brightness temperatures (Tb) more than twice higher than in previous single-Gaussian fitting. They correspond to the dust emission that we find centrally concentrated in both nuclei by subtracting the plasma emission measured at 33 GHz. The dust emission in the western nucleus is found to have a peak Tb ~ 530 K and a full width at half maximum of about 20 pc. This component is estimated to have a bolometric luminosity on the order of 10^{12.5} Lsun and a 20 pc-scale luminosity surface density 10^{15.5} Lsun/kpc^2. A luminous AGN is a plausible energy source for these high values while other explanations remain to be explored. Our continuum image also reveals a third structural component of the western nucleus --- a pair of faint spurs perpendicular to the disk major axis. We attribute it to a bipolar outflow from the highly inclined (i ~ 60 deg) western nuclear disk. △ Less

Submitted 25 September, 2017; originally announced September 2017.

Comments: ApJ in press. (16 pages, 5 figures)

Submillimeter Array CO(2-1) Imaging of the NGC 6946 Giant Molecular Clouds

Authors: Ya-Lin Wu, Kazushi Sakamoto, Hsi-An Pan

Abstract: We present a CO(2-1) mosaic map of the spiral galaxy NGC 6946 by combining data from the Submillimeter Array and the IRAM 30 m telescope. We identify 390 giant molecular clouds (GMCs) from the nucleus to 4.5 kpc in the disk. GMCs in the inner 1 kpc are generally more luminous and turbulent, some of which have luminosities >10^6 K km/s pc^2 and velocity dispersions >10 km/s. Large-scale bar-driven… ▽ More We present a CO(2-1) mosaic map of the spiral galaxy NGC 6946 by combining data from the Submillimeter Array and the IRAM 30 m telescope. We identify 390 giant molecular clouds (GMCs) from the nucleus to 4.5 kpc in the disk. GMCs in the inner 1 kpc are generally more luminous and turbulent, some of which have luminosities >10^6 K km/s pc^2 and velocity dispersions >10 km/s. Large-scale bar-driven dynamics likely regulate GMC properties in the nuclear region. Similar to the Milky Way and other disk galaxies, GMC mass function of NGC 6946 has a shallower slope (index>-2) in the inner region, and a steeper slope (index<-2) in the outer region. This difference in mass spectra may be indicative of different cloud formation pathways: gravitational instabilities might play a major role in the nuclear region, while cloud coalescence might be dominant in the outer disk. Finally, the NGC 6946 clouds are similar to those in M33 in terms of statistical properties, but they are generally less luminous and turbulent than the M51 clouds. △ Less

Submitted 9 April, 2017; originally announced April 2017.

Comments: Published in ApJ

Journal ref: The Astrophysical Journal, Volume 839, article id. 6, 10 pp. (2017)

Luminous, pc-scale CO 6-5 emission in the obscured nucleus of NGC1377

Authors: S. Aalto, S. Muller, F. Costagliola, K. Sakamoto, J. S. Gallagher, N. Falstad, S. König, K. Dasyra, K. Wada, F. Combes, S. García-Burillo, L. E. Kristensen, S. Martín, P. van der Werf, A. S. Evans, J. Kotilainen

Abstract: High resolution submm observations are important in probing the morphology, column density and dynamics of obscured active galactic nuclei (AGNs). With high resolution (0.06 x 0.05) ALMA 690 GHz observations we have found bright (TB >80 K) and compact (FWHM 10x7 pc) CO 6-5 line emission in the nucleus of the extremely radio-quiet galaxy NGC1377. The integrated CO 6-5 intensity is aligned with the… ▽ More High resolution submm observations are important in probing the morphology, column density and dynamics of obscured active galactic nuclei (AGNs). With high resolution (0.06 x 0.05) ALMA 690 GHz observations we have found bright (TB >80 K) and compact (FWHM 10x7 pc) CO 6-5 line emission in the nucleus of the extremely radio-quiet galaxy NGC1377. The integrated CO 6-5 intensity is aligned with the previously discovered jet/outflow of NGC1377 and is tracing the dense (n>1e4 cm-3), hot gas at the base of the outflow. The velocity structure is complex and shifts across the jet/outflow are discussed in terms of jet-rotation or separate, overlapping kinematical components. High velocity gas (deltaV +-145 km/s) is detected inside r<2-3 pc and we suggest that it is emerging from an inclined rotating disk or torus of position angle PA=140+-20 deg with a dynamical mass of approx 3e6 Msun. This mass is consistent with that of a supermassive black hole (SMBH), as inferred from the M-sigma relation. The gas mass of the proposed disk/torus constitutes <3% of the nuclear dynamical mass. In contrast to the intense CO 6-5 line emission, we do not detect dust continuum with an upper limit of S(690GHz)<2mJy. The corresponding, 5 pc, H2 column density is estimated to N(H2)<3e23 cm-2, which is inconsistent with a Compton Thick (CT) source. We discuss the possibility that CT obscuration may be occuring on small (subparsec) or larger scales. From SED fitting we suggest that half of the IR emission of NGC1377 is nuclear and the rest (mostly the far-infrared (FIR)) is more extended. The extreme radio quietness, and the lack of emission from other star formation tracers, raise questions on the origin of the FIR emission. We discuss the possibility that it is arising from the dissipation of shocks in the molecular jet/outflow or from irradiation by the nuclear source along the poles. △ Less

Submitted 17 February, 2017; originally announced February 2017.

Comments: 7 pages, 5 figures, submitted to Astronomy and Astrophysics

Sub-arcsecond imaging of the water emission in Arp 220

Authors: S. König, S. Martín, S. Muller, J. Cernicharo, K. Sakamoto, L. K. Zschaechner, E. M. L. Humphreys, T. Mroczkowski, M. Krips, M. Galametz, S. Aalto, W. H. T. Vlemmings, J. Ott, D. S. Meier, A. Fuente, S. García-Burillo, R. Neri

Abstract: Extragalactic observations of water emission can provide valuable insights into the excitation of the interstellar medium. In addition, extragalactic megamasers are powerful probes of kinematics close to active nuclei. Therefore, it is paramount to determine the true origin of the water emission, whether it is excited by processes close to an AGN or in star-forming regions. We use ALMA Band 5 scie… ▽ More Extragalactic observations of water emission can provide valuable insights into the excitation of the interstellar medium. In addition, extragalactic megamasers are powerful probes of kinematics close to active nuclei. Therefore, it is paramount to determine the true origin of the water emission, whether it is excited by processes close to an AGN or in star-forming regions. We use ALMA Band 5 science verification observations to analyse the emission of the 183 GHz water line in Arp 220 on sub-arcsecond scales, in conjunction with new ALMA Band 7 data at 325 GHz. Specifically, the nature of the process leading to the excitation of emission at these water lines is studied in this context. Supplementary 22 GHz VLA observations are used to better constrain the parameter space in the excitation modelling of the water lines. We detect 183 GHz H2O and 325 GHz water emission towards the two compact nuclei at the center of Arp 220, being brighter in Arp 220 West. The emission at these two frequencies is compared to previous single-dish data and does not show evidence of variability. The 183 and 325 GHz lines show similar spectra and kinematics, but the 22 GHz profile is significantly different in both nuclei due to a blend with an NH3 absorption line. Our findings suggest that the most likely scenario to cause the observed water emission in Arp 220 is a large number of independent masers originating from numerous star-forming regions. △ Less

Submitted 8 March, 2017; v1 submitted 22 December, 2016; originally announced December 2016.

Comments: 10 pages, 7 figures, 3 tables, accepted for publication in A&A

Journal ref: A&A 602, A42 (2017)

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

The unbearable opaqueness of Arp 220

Authors: S. Martín, S. Aalto, K. Sakamoto, E. González-Alfonso, S. Muller, C. Henkel, S. García-Burillo, R. Aladro, F. Costagliola, N. Harada, M. Krips, J. Martín-Pintado, S. Mühle, P. van der Werf, S. Viti

Abstract: We explore the potential of imaging vibrationally excited molecular emission at high angular resolution to better understand the morphology and physical structure of the dense gas in Arp~220 and to gain insight into the nature of the nuclear powering sources. Vibrationally excited emission of HCN is detected in both nuclei with a very high ratio relative to the total $L_{FIR}$, higher than in any… ▽ More We explore the potential of imaging vibrationally excited molecular emission at high angular resolution to better understand the morphology and physical structure of the dense gas in Arp~220 and to gain insight into the nature of the nuclear powering sources. Vibrationally excited emission of HCN is detected in both nuclei with a very high ratio relative to the total $L_{FIR}$, higher than in any other observed galaxy and well above what is observed in Galactic hot cores. HCN $v_2=1f$ is observed to be marginally resolved in $\sim60\times50$~pc regions inside the dusty $\sim100$~pc sized nuclear cores. Its emission is centered on our derived individual nuclear velocities based on HCO$^+$ emission ($V_{WN}=5342\pm4$ and $V_{EN}=5454\pm8$~\kms, for the western and eastern nucleus, respectively). With virial masses within $r\sim25-30$~pc based on the HCN~$v_2=1f$ line widths, we estimate gas surface densities (gas fraction $f_g=0.1$) of $3\pm0.3\times10^4~M_\odot~\rm pc^{-2}$ (WN) and $1.1\pm0.1\times10^4~M_\odot~\rm pc^{-2}$ (EN). The $4-3/3-2$ flux density ratio could be consistent with optically thick emission, which would further constrain the size of the emitting region to $>15$~pc (EN) and $>22$~pc (WN). The absorption systems that may hide up to $70\%$ of the HCN and HCO$^+$ emission are found at velocities of $-50$~\kms~(EN) and $6$, $-140$, and $-575$~\kms (WN) relative to velocities of the nuclei. Blueshifted absorptions are the evidence of outflowing motions from both nuclei. The bright vibrational emission implies the existence of a hot dust region radiatively pumping these transitions. We find evidence of a strong temperature gradient that would be responsible for both the HCN $v_2$ pumping and the absorbed profiles from the vibrational ground state as a result of both continuum and self-absorption by cooler foreground gas. △ Less

Submitted 3 March, 2016; originally announced March 2016.

Comments: 14 pages accepted for publication in A&A

Journal ref: A&A 590, A25 (2016)

SMA Observations of the Extended CO(6-5) Emission in the Starburst Galaxy NGC253

Authors: Melanie Krips, Sergio Martin, Alison Peck, Kazushi Sakamoto, Roberto Neri, Mark Gurwell, Glen Petitpas, Jun-Hui Zhao

Abstract: We present observations of the $^{12}$CO(6-5) line and 686GHz continuum emission in NGC253 with the Submillimeter Array at an angular resolution of ~4arcsec. The $^{12}$CO(6-5) emission is clearly detected along the disk and follows the distribution of the lower $^{12}$CO line transitions with little variations of the line ratios in it. A large-velocity gradient analysis suggests a two-temperature… ▽ More We present observations of the $^{12}$CO(6-5) line and 686GHz continuum emission in NGC253 with the Submillimeter Array at an angular resolution of ~4arcsec. The $^{12}$CO(6-5) emission is clearly detected along the disk and follows the distribution of the lower $^{12}$CO line transitions with little variations of the line ratios in it. A large-velocity gradient analysis suggests a two-temperature model of the molecular gas in the disk, likely dominated by a combination of low-velocity shocks and the disk wide PDRs. Only marginal $^{12}$CO(6-5) emission is detected in the vicinity of the expanding shells at the eastern and western edges of the disk. While the eastern shell contains gas even warmer (T$_{\rm kin}$>300~K) than the hot gas component (T$_{\rm kin}$=300K) of the disk, the western shell is surrounded by gas much cooler (T$_{\rm kin}$=60K) than the eastern shell but somewhat hotter than the cold gas component of the disk (for similar H$_2$ and CO column densities), indicative of different (or differently efficient) heating mechansisms. The continuum emission at 686GHz in the disk agrees well in shape and size with that at lower (sub-)millimeter frequencies, exhibiting a spectral index consistent with thermal dust emission. We find dust temperatures of ~10-30K and largely optically thin emission. However, our fits suggest a second (more optically thick) dust component at higher temperatures (T$_{\rm d}$>60K), similar to the molecular gas. We estimate a global dust mass of ~10$^6$Msun for the disk translating into a gas-to-dust mass ratio of a few hundred consistent with other nearby active galaxies. △ Less

Submitted 28 January, 2016; v1 submitted 21 January, 2016; originally announced January 2016.

Comments: accepted for publication in ApJ, 21 pages, 12 Figures, 6 Tables; acknowledgements updated in revised version

A precessing molecular jet signaling an obscured, growing supermassive black hole in NGC1377?

Authors: S. Aalto, F. Costagliola, S. Muller, K. Sakamoto, J. S. Gallagher, K. Dasyra, K. Wada, F. Combes, S. García-Burillo, L. Kristensen, S. Martín, P. van der Werf, A. S. Evans, J. Kotilainen

Abstract: With high resolution (0."25 x 0."18) ALMA CO 3-2 observations of the nearby (D=21 Mpc, 1"=102 pc), extremely radio-quiet galaxy NGC1377, we have discovered a high-velocity, very collimated nuclear outflow which we interpret as a molecular jet with a projected length of +-150 pc. Along the jet axis we find strong velocity reversals where the projected velocity swings from -150 km/s to +150 km/s. A… ▽ More With high resolution (0."25 x 0."18) ALMA CO 3-2 observations of the nearby (D=21 Mpc, 1"=102 pc), extremely radio-quiet galaxy NGC1377, we have discovered a high-velocity, very collimated nuclear outflow which we interpret as a molecular jet with a projected length of +-150 pc. Along the jet axis we find strong velocity reversals where the projected velocity swings from -150 km/s to +150 km/s. A simple model of a molecular jet precessing around an axis close to the plane of the sky can reproduce the observations. The velocity of the outflowing gas is difficult to constrain due to the velocity reversals but we estimate it to be between 240 and 850 km/s and the jet to precess with a period P=0.3-1.1 Myr. The CO emission is clumpy along the jet and the total molecular mass in the high-velocity (+-(60 to 150 km/s)) gas lies between 2e6 Msun (light jet) and 2e7 Msun (massive jet). There is also CO emission extending along the minor axis of NGC1377. It holds >40% of the flux in NGC1377 and may be a slower, wide-angle molecular outflow which is partially entrained by the molecular jet. We discuss the driving mechanism of the molecular jet and suggest that it is either powered by a very faint radio jet or by an accretion disk-wind similar to those found towards protostars. The nucleus of NGC1377 harbours intense embedded activity and we detect emission from vibrationally excited HCN J=4-3 v_2=1f which is consistent with hot gas and dust. We find large columns of H2 in the centre of NGC1377 which may be a sign of a high rate of recent gas infall. The dynamical age of the molecular jet is short (<1 Myr), which could imply that it is young and consistent with the notion that NGC1377 is caught in a transient phase of its evolution. However, further studies are required to determine the age of the molecular jet, its mass and the role it is playing in the growth of the nucleus of NGC1377. △ Less

Submitted 28 April, 2016; v1 submitted 29 October, 2015; originally announced October 2015.

Comments: This is a revised and expanded version of a previous submission which now has 13 pages, 6 figures (+ 4 in the Appendix) and is accepted for publication in Astronomy & Astrophysics

Exploring the molecular chemistry and excitation in obscured luminous infrared galaxies: An ALMA mm-wave spectral scan of NGC 4418

Authors: F. Costagliola, K. Sakamoto, S. Muller, S. Martín, S. Aalto, N. Harada, P. van der Werf, S. Viti, S. Garcia-Burillo, M. Spaans

Abstract: We obtained an ALMA Cycle 0 spectral scan of the dusty LIRG NGC 4418, spanning a total of 70.7 GHz in bands 3, 6, and 7. We use a combined local thermal equilibrium (LTE) and non-LTE (NLTE) fit of the spectrum in order to identify the molecular species and derive column densities and excitation temperatures. We derive molecular abundances and compare them with other Galactic and extragalactic sour… ▽ More We obtained an ALMA Cycle 0 spectral scan of the dusty LIRG NGC 4418, spanning a total of 70.7 GHz in bands 3, 6, and 7. We use a combined local thermal equilibrium (LTE) and non-LTE (NLTE) fit of the spectrum in order to identify the molecular species and derive column densities and excitation temperatures. We derive molecular abundances and compare them with other Galactic and extragalactic sources by means of a principal component analysis. We detect 317 emission lines from a total of 45 molecular species, including 15 isotopic substitutions and six vibrationally excited variants. Our LTE/NLTE fit find kinetic temperatures from 20 to 350 K, and densities between 10$^5$ and 10$^7$ cm$^{-3}$. The spectrum is dominated by vibrationally excited HC$_3$N, HCN, and HNC, with vibrational temperatures from 300 to 450 K. We find high abundances of HC$_3$N, SiO, H$_2$S, and c-HCCCH and a low CH$_3$OH abundance. A principal component analysis shows that NGC 4418 and Arp 220 share very similar molecular abundances and excitation, which clearly set them apart from other Galactic and extragalactic environments. The similar molecular abundances observed towards NCG 4418 and Arp 220 are consistent with a hot gas-phase chemistry, with the relative abundances of SiO and CH$_3$OH being regulated by shocks and X-ray driven dissociation. The bright emission from vibrationally excited species confirms the presence of a compact IR source, with an effective diameter $<$5 pc and brightness temperatures $>$350 K. The molecular abundances and the vibrationally excited spectrum are consistent with a young AGN/starburst system. We suggest that NGC 4418 may be a template for a new kind of chemistry and excitation, typical of compact obscured nuclei (CON). Because of the narrow line widths and bright molecular emission, NGC 4418 is the ideal target for further studies of the chemistry in CONs. △ Less

Submitted 30 June, 2015; originally announced June 2015.

Comments: accepted by A&A on 29/06/2015

Journal ref: A&A 582, A91 (2015)

Probing highly obscured, self-absorbed galaxy nuclei with vibrationally excited HCN

Authors: S. Aalto, S. Martin F. Costagliola, E. Gonzalez-Alfonso, S. Muller, K. Sakamoto, G. A. Fuller, S. Garcia-Burillo, P. van der Werf, R. Neri, M. Spaans, F. Combes, S. Viti, S. Muehle, L. Armus, A. Evans, E. Sturm, J. Cernicharo, C. Henkel, T. R. Greve

Abstract: We present high resolution (0."4) IRAM PdBI and ALMA mm and submm observations of the (ultra) luminous infrared galaxies ((U)LIRGs) IRAS17208-0014, Arp220, IC860 and Zw049.057 that reveal intense line emission from vibrationally excited ($ν_2$=1) J=3-2 and 4-3 HCN. The emission is emerging from buried, compact (r<17-70 pc) nuclei that have very high implied mid-infrared surface brightness $>$… ▽ More We present high resolution (0."4) IRAM PdBI and ALMA mm and submm observations of the (ultra) luminous infrared galaxies ((U)LIRGs) IRAS17208-0014, Arp220, IC860 and Zw049.057 that reveal intense line emission from vibrationally excited ($ν_2$=1) J=3-2 and 4-3 HCN. The emission is emerging from buried, compact (r<17-70 pc) nuclei that have very high implied mid-infrared surface brightness $>$$5\times 10^{13}$ L$_{\odot}$ kpc$^{-2}$. These nuclei are likely powered by accreting supermassive black holes (SMBHs) and/or hot (>200 K) extreme starbursts. Vibrational, $ν_2$=1, lines of HCN are excited by intense 14 micron mid-infrared emission and are excellent probes of the dynamics, masses, and physical conditions of (U)LIRG nuclei when H$_2$ column densities exceed $10^{24}$ cm$^{-2}$. It is clear that these lines open up a new interesting avenue to gain access to the most obscured AGNs and starbursts. Vibrationally excited HCN acts as a proxy for the absorbed mid-infrared emission from the embedded nuclei, which allows for reconstruction of the intrinsic, hotter dust SED. In contrast, we show strong evidence that the ground vibrational state ($ν$=0), J=3-2 and 4-3 rotational lines of HCN and HCO$^+$ fail to probe the highly enshrouded, compact nuclear regions owing to strong self- and continuum absorption. The HCN and HCO$^+$ line profiles are double-peaked because of the absorption and show evidence of non-circular motions - possibly in the form of in- or outflows. Detections of vibrationally excited HCN in external galaxies are so far limited to ULIRGs and early-type spiral LIRGs, and we discuss possible causes for this. We tentatively suggest that the peak of vibrationally excited HCN emission is connected to a rapid stage of nuclear growth, before the phase of strong feedback. △ Less

Submitted 9 August, 2015; v1 submitted 26 April, 2015; originally announced April 2015.

Comments: 13 pages, 7 figures, accepted for publication in Astronomy and Astrophysics

Infrared and X-Ray Evidence of an AGN in the NGC 3256 Southern Nucleus

Authors: Youichi Ohyama, Yuichi Terashima, Kazushi Sakamoto

Abstract: We investigate signs of Active Galactic Nucleus (AGN) in the luminous infrared galaxy NGC 3256 at both infrared and X-ray wavelengths. NGC 3256 has double, the Northern and Southern, nuclei (hereafter, N and S nuclei, respectively). We show that the Spitzer IRAC colors extracted at the S nucleus are AGN-like, and the Spitzer IRS spectrum is bluer at <6um than at the N nucleus. We built for the S n… ▽ More We investigate signs of Active Galactic Nucleus (AGN) in the luminous infrared galaxy NGC 3256 at both infrared and X-ray wavelengths. NGC 3256 has double, the Northern and Southern, nuclei (hereafter, N and S nuclei, respectively). We show that the Spitzer IRAC colors extracted at the S nucleus are AGN-like, and the Spitzer IRS spectrum is bluer at <6um than at the N nucleus. We built for the S nucleus an AGN-starburst composite model with a heavily absorbed AGN to successfully reproduce not only the IRAC and IRS specrophotometries at ~3arcsec but also the very deep silicate 9.7um absorption observed at 0.36" scale by Diaz-Santos et al. We found a 2.2um compact source at the S nucleus in a HST NICMOS image and identified its unresolved core (at 0.26" resolution) with the compact core in previous mid-infrared observations at comparable resolution. The flux of the 2.2umm core is consistent with our AGN spectral energy distribution model. We also analyzed a deeper than ever Chandra X-ray spectrum of the unresolved (at 0.5" resolution) source at the S nucleus. We found that a dual-component power-law model (for primary and scattered ones) fits an apparently very hard spectrum with a moderately large absorption on the primary component. Together with a limit on equivalent width of a fluorescent Fe-K emission line at 6.4 keV, the X-ray spectrum is consistent with a typical Compton-thin Seyfert 2. We therefore suggest that the S nucleus hosts a heavily absorbed low-luminosity AGN. △ Less

Submitted 30 March, 2015; originally announced March 2015.

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

High-Resolution Radio Continuum Measurements of the Nuclear Disks of Arp 220

Authors: Loreto Barcos-Muñoz, Adam K. Leroy, Aaron S. Evans, George C. Privon, Lee Armus, Jim Condon, Joseph M. Mazzarella, David S. Meier, Emmanuel Momjian, Eric J. Murphy, Juerguen Ott, Ashely Reichardt, Kazushi Sakamoto, David B. Sanders, Eva Schinnerer, Sabrina Stierwalt, Jason A. Surace, Todd A. Thompson, Fabian Walter

Abstract: We present new Karl G. Jansky Very Large Array radio continuum images of the nuclei of Arp 220, the nearest ultra-luminous infrared galaxy. These images have both the angular resolution to study detailed morphologies of the two nuclei that power the system and sensitivity to a wide range of spatial scales. At 33 GHz, and with a resolution of 0".081 x 0".063 (29.9 x 23.3 pc), we resolve the emissio… ▽ More We present new Karl G. Jansky Very Large Array radio continuum images of the nuclei of Arp 220, the nearest ultra-luminous infrared galaxy. These images have both the angular resolution to study detailed morphologies of the two nuclei that power the system and sensitivity to a wide range of spatial scales. At 33 GHz, and with a resolution of 0".081 x 0".063 (29.9 x 23.3 pc), we resolve the emission surrounding both nuclei and conclude that is mostly synchrotron in nature. The spatial distributions of radio emission in both nuclei are well described by exponential profiles. These have deconvolved half-light radii of 51 and 35 pc for the eastern and western nuclei, and they match the number density profile of radio supernovae observed with very long baseline interferometry. This similarity might be due to the fast cooling of cosmic rays electrons caused by the presence of a strong (~ mG) magnetic field in this system. We estimate high luminosity surface densities of $\mathrm{Σ_{IR} \sim 4.2^{+1.6}_{-0.7} \times 10^{13}}$ (east) and $\mathrm{\sim 9.7^{+3.7}_{-2.4} \times 10^{13}~(west)~L_{\odot}~kpc^{-2}}$, and star formation rate surface densities of $\mathrm{Σ_{SFR} \sim 10^{3.7\pm0.1}}$ (east) and $\mathrm{\sim 10^{4.1\pm0.1}~(west)~M_{\odot}~yr^{-1}~kpc^{-2}}$. These values, especially for the western nucleus are, to our knowledge, the highest luminosity and star formation rate surface densities measured for any star-forming system. Despite these high values, the nuclei lie below the dusty Eddington limit in which radiation pressure is balanced only by self-gravity. The small measured sizes also imply that the nuclei of Arp 220 are only transparent in the frequency range ~ 5 to 350 GHz. Our results offer no clear evidence that an active galactic nucleus dominates the emission from either nucleus at 33 GHz. △ Less

Submitted 4 November, 2014; originally announced November 2014.

Comments: 15 pages, 3 Figures, 3 Tables and 2 appendices. Accepted for publication in ApJ

A Luminous Infrared Merger with Two Bipolar Molecular Outflows: ALMA and SMA Observations of NGC 3256

Authors: Kazushi Sakamoto, Susanne Aalto, Francoise Combes, Aaron Evans, Alison Peck

Abstract: We report ALMA and SMA observations of the luminous infrared merger NGC 3256, the most luminous galaxy within z=0.01. Both of the two merger nuclei separated by 5" (0.8 kpc) on the sky have a compact concentration of molecular gas, i.e., nuclear disks with Sigma_mol > 10^3 Msun pc^-2. The one at the northern nucleus is face-on while the southern nuclear disk is almost edge-on. The northern nucleus… ▽ More We report ALMA and SMA observations of the luminous infrared merger NGC 3256, the most luminous galaxy within z=0.01. Both of the two merger nuclei separated by 5" (0.8 kpc) on the sky have a compact concentration of molecular gas, i.e., nuclear disks with Sigma_mol > 10^3 Msun pc^-2. The one at the northern nucleus is face-on while the southern nuclear disk is almost edge-on. The northern nucleus is more massive and has molecular arcs and spiral arms around. The high-velocity molecular gas previously found in the system is resolved to two molecular outflows associated with each of the two nuclei. The molecular outflow from the northern nuclear disk is part of a starburst-driven superwind seen nearly pole on. Its maximum velocity is >750 km/s and its mass outflow rate is estimated to be > 60 Msun/yr for a conversion factor N_{H_2}/I_{CO(1-0)}=1x10^20 cm^-2/(K km/s). The outflow from the southern nucleus is a highly collimated bipolar molecular jet seen nearly edge-on. Its line-of-sight velocity increases with distance out to 300 pc from the southern nucleus. Its maximum de-projected velocity is ~2000 km/s for the estimated inclination and should exceed 1000 km/s even allowing for its uncertainty. The mass outflow rate is estimated to be >50 Msun/yr for this outflow. There are possible signs that this southern outflow has been driven by a bipolar radio jet from an AGN that became inactive very recently. The sum of these outflow rates, although subject to the uncertainty in the molecular mass estimate, either exceeds or compares to the total star formation rate in NGC 3256. The feedback from nuclear activities in the form of molecular outflows is therefore significant in the gas consumption budget, and hence evolution, of this luminous infrared galaxy. (abridged) △ Less

Submitted 27 March, 2014; originally announced March 2014.

Comments: 36 pages, 21 figures, submitted to ApJ on Mar. 4, 2014

Journal ref: The Astrophysical Journal 797 (2014) 90

A high-resolution mm and cm study of the obscured LIRG NGC 4418 - A compact obscured nucleus fed by in-falling gas?

Authors: F. Costagliola, S. Aalto, K. Sakamoto, S. Martín, R. Beswick, S. Muller, H. -R. Klöckner

Abstract: The aim of this study is to constrain the dynamics, structure and feeding of the compact nucleous of NGC4418, and to reveal the nature of the main hidden power source: starburst or AGN. We obtained high spatial resolution observations of NGC4418 at 1.4 and 5 GHz with MERLIN, and at 230 and 270 GHz with the SMA very extended configuration. We use the continuum morphology and flux density to estimat… ▽ More The aim of this study is to constrain the dynamics, structure and feeding of the compact nucleous of NGC4418, and to reveal the nature of the main hidden power source: starburst or AGN. We obtained high spatial resolution observations of NGC4418 at 1.4 and 5 GHz with MERLIN, and at 230 and 270 GHz with the SMA very extended configuration. We use the continuum morphology and flux density to estimate the size of the emitting region, the star formation rate and the dust temperature. Emission lines are used to study the kinematics through position-velocity diagrams. Molecular emission is studied with population diagrams and by fitting an LTE synthetic spectrum. We detect bright 1mm line emission from CO, HC3N, HNC and C34S, and 1.4 GHz absorption from HI. The CO 2-1 emission and HI absorption can be fit by two velocity components at 2090 and 2180 km s-1. We detect vibrationally excited HC3N and HNC, with Tvib 300K. Molecular excitation is consistent with a layered temperature structure, with three main components at 80, 160 and 300 K. For the hot component we estimate a source size of less than 5 pc. The nuclear molecular gas surface density of 1e4 Msun pc-2 is extremely high, and similar to that found in the ultra-luminous infrared galaxy (ULIRG) Arp220. Our observations confirm the the presence of a molecular and atomic in-flow, previously suggested by Herschel observations, which is feeding the activity in the center of NGC4418. Molecular excitation confirms the presence of a very compact, hot dusty core. If a starburst is responsible for the observed IR flux, this has to be at least as extreme as the one in Arp220, with an age of 3-10 Myr and a star formation rate >10 Msun yr-1. If an AGN is present, it must be extremely Compton-thick. △ Less

Submitted 10 June, 2013; originally announced June 2013.

Comments: 18 pages, 11 figures, Accepted for publication by A&A on 10/6/2013