Showing 1–50 of 74 results for author: Harada, N

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

Internal 1000 AU-scale Structures of the R CrA Cluster-forming Cloud -- I: Filamentary Structures

Authors: Kengo Tachihara, Naofumi Fukaya, Kazuki Tokuda, Yasumasa Yamasaki, Takeru Nishioka, Daisei Abe, Tsuyoshi Inoue, Naoto Harada, Ayumu Shoshi, Shingo Nozaki, Asako Sato, Mitsuki Omura, Kakeru Fujishiro, Misato Fukagawa, Masahiro N. Machida, Takahiro Kanai, Yumiko Oasa, Toshikazu Onishi, Kazuya Saigo, Yasuo Fukui

Abstract: We report on ALMA ACA observations of a high-density region of the Corona Australis cloud forming a young star cluster, and the results of resolving internal structures. In addition to embedded Class 0/I protostars in continuum, a number of complex dense filamentary structures are detected in the C18O and SO lines by the 7m array. These are sub-structures of the molecular clump that are detected b… ▽ More We report on ALMA ACA observations of a high-density region of the Corona Australis cloud forming a young star cluster, and the results of resolving internal structures. In addition to embedded Class 0/I protostars in continuum, a number of complex dense filamentary structures are detected in the C18O and SO lines by the 7m array. These are sub-structures of the molecular clump that are detected by the TP array as the extended emission. We identify 101 and 37 filamentary structures with a few thousand AU widths in C18O and SO, respectively, called as feathers. The typical column density of the feathers in C18O is about 10^{22} cm^{-2}, and the volume density and line mass are ~ 10^5 cm^{-3}, and a few times M_{sun} pc^{-1}, respectively. This line mass is significantly smaller than the critical line mass expected for cold and dense gas. These structures have complex velocity fields, indicating a turbulent internal property. The number of feathers associated with Class 0/I protostars is only ~ 10, indicating that most of them do not form stars but rather being transient structures. The formation of feathers can be interpreted as a result of colliding gas flow as the morphology well reproduced by MHD simulations, supported by the the presence of HI shells in the vicinity. The colliding gas flows may accumulate gas and form filaments and feathers, and trigger the active star formation of the R CrA cluster. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 24 pages, 13 figures; Accepted for publication in ApJ

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

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

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

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

Submitted 25 March, 2024; originally announced March 2024.

Comments: 9 pages, 10 figures

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

Discovery of Asymmetric Spike-like Structures of the 10 au Disk around the Very Low-luminosity Protostar Embedded in the Taurus Dense Core MC 27/L1521F with ALMA

Authors: Kazuki Tokuda, Naoto Harada, Mitsuki Omura, Tomoaki Matsumoto, Toshikazu Onishi, Kazuya Saigo, Ayumu Shoshi, Shingo Nozaki, Kengo Tachihara, Naofumi Fukaya, Yasuo Fukui, Shu-ichiro Inutsuka, Masahiro N. Machida

Abstract: Recent Atacama Large Millimeter/submillimeter Array (ALMA) observations have revealed an increasing number of compact protostellar disks with radii of less than a few tens of astronomical units and that young Class 0/I objects have an intrinsic size diversity. To deepen our understanding of the origin of such tiny disks, we performed the highest-resolution configuration observations with ALMA at a… ▽ More Recent Atacama Large Millimeter/submillimeter Array (ALMA) observations have revealed an increasing number of compact protostellar disks with radii of less than a few tens of astronomical units and that young Class 0/I objects have an intrinsic size diversity. To deepen our understanding of the origin of such tiny disks, we performed the highest-resolution configuration observations with ALMA at a beam size of $\sim$0$''$03 (4 au) on the very low-luminosity Class 0 protostar embedded in the Taurus dense core MC 27/L1521F. The 1.3 mm continuum measurement successfully resolved a tiny, faint ($\sim$1 mJy) disk with a major axis length of $\sim$10 au, one of the smallest examples in the ALMA protostellar studies. In addition, we detected spike-like components in the northeastern direction at the disk edge. Gravitational instability or other fragmentation mechanisms cannot explain the structures, given the central stellar mass of $\sim$0.2 $M_{\odot}$ and the disk mass of $\gtrsim$10$^{-4}$ $M_{\odot}$. Instead, we propose that these small spike structures were formed by a recent dynamic magnetic flux transport event due to interchange instability that would be favorable to occur if the parental core has a strong magnetic field. The presence of complex arc-like structures on a larger ($\sim$2000 au) scale in the same direction as the spike structures suggests that the event was not single. Such episodic, dynamical events may play an important role in maintaining the compact nature of the protostellar disk in the complex gas envelope during the main accretion phase. △ Less

Submitted 3 April, 2024; v1 submitted 1 March, 2024; originally announced March 2024.

Comments: 13 pages, 5 figures, Accepted for publication in ApJ

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

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

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

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

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

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

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)

An Extremely Young Protostellar Core, MMS 1/ OMC-3: Episodic Mass Ejection History Traced by the Micro SiO Jet

Authors: Satoko Takahashi, Masahiro N. Machida, Mitsuki Omura, Doug Johnstone, Kazuya Saigo, Naoto Harada, Kohji Tomisaka, Paul T. P. Ho, Luis A. Zapata, Steve Mairs, Gregory J. Herczeg, Kotomi Taniguchi, Yuhua Liu, Asako Sato

Abstract: We present ${\sim}0.2$ arcsec ($\sim$80 au) resolution observations of the CO (2-1) and SiO (5-4) lines made with the Atacama large millimeter/submillimeter array toward an extremely young intermediate-mass protostellar source (t$_{\rm dyn}<$1000 years), MMS 1 located in the Orion Molecular Cloud-3 region. We have successfully imaged a very compact CO molecular outflow associated with MMS 1, havin… ▽ More We present ${\sim}0.2$ arcsec ($\sim$80 au) resolution observations of the CO (2-1) and SiO (5-4) lines made with the Atacama large millimeter/submillimeter array toward an extremely young intermediate-mass protostellar source (t$_{\rm dyn}<$1000 years), MMS 1 located in the Orion Molecular Cloud-3 region. We have successfully imaged a very compact CO molecular outflow associated with MMS 1, having deprojected lobe sizes of $\sim$18000 au (red-shifted lobe) and $\sim$35000 au (blue-shifted lobe). We have also detected an extremely compact ($\lesssim$1000 au) and collimated SiO protostellar jet within the CO outflow. The maximum deprojected jet speed is measured to be as high as 93 km s$^{-1}$. The SiO jet wiggles and displays a chain of knots. Our detection of the molecular outflow and jet is the first direct evidence that MMS 1 already hosts a protostar. The position-velocity diagram obtained from the SiO emission shows two distinct structures: (i) bow-shocks associated with the tips of the outflow, and (ii) a collimated jet, showing the jet velocities linearly increasing with the distance from the driving source. Comparisons between the observations and numerical simulations quantitatively share similarities such as multiple-mass ejection events within the jet and Hubble-like flow associated with each mass ejection event. Finally, while there is a weak flux decline seen in the 850 $μ$m light curve obtained with JCMT/SCUBA 2 toward MMS 1, no dramatic flux change events are detected. This suggests that there has not been a clear burst event within the last 8 years. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: 19 pages, 9 figures, Accepted for publication in ApJ

The ALMaQUEST Survey XII: Dense Molecular Gas as traced by HCN and HCO$^{+}$ in Green Valley Galaxies

Authors: Lihwai Lin, Hsi-An Pan, Sara L. Ellison, Nanase Harada, Maria J. Jimenez-Donaire, K. Decker French, William M. Baker, Bau-Ching Hsieh, Yusei Koyama, Carlos Lopez-Coba, Tomonari Michiyama, Kate Rowlands, Sebastian F. Sanchez, Mallory Thorp

Abstract: We present ALMA observations of two dense gas tracers, HCN(1-0) and HCO$^{+}$(1-0), for three galaxies in the green valley and two galaxies on the star-forming main sequence with comparable molecular gas fractions as traced by the CO(1-0) emissions, selected from the ALMaQUEST survey. We investigate whether the deficit of molecular gas star formation efficiency (SFE$_{\rm mol}$) that leads to the… ▽ More We present ALMA observations of two dense gas tracers, HCN(1-0) and HCO$^{+}$(1-0), for three galaxies in the green valley and two galaxies on the star-forming main sequence with comparable molecular gas fractions as traced by the CO(1-0) emissions, selected from the ALMaQUEST survey. We investigate whether the deficit of molecular gas star formation efficiency (SFE$_{\rm mol}$) that leads to the low specific star formation rate in these green valley galaxies is due to a lack of dense gas (characterized by the dense gas fraction $f_{\rm dense}$) or the low star formation efficiency of dense gas (SFE$_{\rm dense}$). We find that SFE$_{\rm mol}$ as traced by the CO emissions, when considering both star-forming and retired spaxels together, is tightly correlated with SFE$_{\rm dense}$ and depends only weakly on $f_{\rm dense}$. The specific star formation rate (sSFR) on kpc scales is primarily driven by SFE$_{\rm mol}$ and SFE$_{\rm dense}$, followed by the dependence on $f_{\rm mol}$, and is least correlated with $f_{\rm dense}$ or the dense-to-stellar mass ratio ($R_{\rm dense}$). When compared with other works in the literature, we find that our green valley sample shows lower global SFE$_{\rm mol}$ as well as lower SFE$_{\rm dense}$ while exhibiting similar dense gas fractions when compared to star-forming and starburst galaxies. We conclude that the star formation of the 3 green valley galaxies with a normal abundance of molecular gas is suppressed mainly due to the reduced SFE$_{\rm dense}$ rather than the lack of dense gas. △ Less

Submitted 11 January, 2024; originally announced January 2024.

Comments: 20 pages, 13 figures, ApJ accepted

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

Ring Gap Structure around Class I Protostar WL 17

Authors: Ayumu Shoshi, Naoto Harada, Kazuki Tokuda, Yoshihiro Kawasaki, Hayao Yamasaki, Asako Sato, Mitsuki Omura, Masayuki Yamaguchi, Kengo Tachihara, Masahiro N. Machida

Abstract: WL 17 is a Class I object and was considered to have a ring-hole structure. We analyzed the structure around WL 17 to investigate the detailed properties of WL 17. We used ALMA archival data, which have a higher angular resolution than previous observations. We investigated the WL 17 system with the 1.3 mm dust continuum and 12CO and C18O (J = 2-1) line emissions. The dust continuum emission showe… ▽ More WL 17 is a Class I object and was considered to have a ring-hole structure. We analyzed the structure around WL 17 to investigate the detailed properties of WL 17. We used ALMA archival data, which have a higher angular resolution than previous observations. We investigated the WL 17 system with the 1.3 mm dust continuum and 12CO and C18O (J = 2-1) line emissions. The dust continuum emission showed a clear ring structure with inner and outer edges of ~11 and ~21 au, respectively. In addition, we detected an inner disk of < 5 au radius enclosing the central star within the ring, the first observation of this structure. Thus, WL 17 has a ring-gap structure, not a ring-hole structure. We did not detect any marked emission in either the gap or inner disk, indicating that there is no sign of a planet, circumplanetary disk, or binary companion. We identified the base of both blue-shifted and red-shifted outflows based on the 12CO emission, which is clearly associated with the disk around WL 17. The outflow mass ejection rate is ~3.6x10^-7 Msun yr-1 and the dynamical timescale is as short as ~ 10^4 yr. The C18O emission showed that an inhomogeneous infalling envelope, which can induce episodic mass accretion, is distributed in the region within ~1000 au from the central protostar. With these new findings, we can constrain the planet formation and dust growth scenarios in the accretion phase of star formation. △ Less

Submitted 5 December, 2023; originally announced December 2023.

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

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

Secondary outflow driven by the protostar Ser-emb 15 in Serpens

Authors: Asako Sato, Kazuki Tokuda, Masahiro N. Machida, Kengo Tachihara, Naoto Harada, Hayao Yamasaki, Shingo Hirano, Toshikazu Onishi, Yuko Matsushita

Abstract: We present the detection of a secondary outflow associated with a Class I source, Ser-emb 15, in the Serpens Molecular Cloud. We reveal two pairs of molecular outflows consisting of three lobes, namely primary and secondary outflows, using ALMA 12CO and SiO line observations at a resolution of 318 au. The secondary outflow is elongated approximately perpendicular to the axis of the primary outflow… ▽ More We present the detection of a secondary outflow associated with a Class I source, Ser-emb 15, in the Serpens Molecular Cloud. We reveal two pairs of molecular outflows consisting of three lobes, namely primary and secondary outflows, using ALMA 12CO and SiO line observations at a resolution of 318 au. The secondary outflow is elongated approximately perpendicular to the axis of the primary outflow in the plane of the sky. We also identify two compact structures, Sources A and B, within an extended structure associated with Ser-emb 15 in the 1.3 mm continuum emission at a resolution of 40 au. The projected sizes of Sources A and B are 137 au and 60 au, respectively. Assuming a dust temperature of 20 K, we estimate the dust mass to be 0.0024 Msun for Source A and 0.00033 Msun for Source B. C18O line data imply the existence of rotational motion around the extended structure, however, cannot resolve rotational motion in Source A and/or B, due to insufficient angular and frequency resolutions. Therefore, we cannot conclude whether Ser-emb 15 is a single or binary system. Thus, either Source A or B could drive the secondary outflow. We discuss two scenarios to explain the driving mechanism of the primary and secondary outflows: the Ser-emb 15 system is (1) a binary system composed of Source A and B or (2) a single star system composed of only Source A. In either case, the system could be a suitable target for investigating the disk and/or binary formation processes in complicated environments. Detecting these outflows should contribute to understanding complex star-forming environments, which may be common in the star-formation processes. △ Less

Submitted 9 October, 2023; originally announced October 2023.

Comments: 21 pages, 10 figures, Accepted for publication in the Astrophysical Journal

An ALMA-resolved view of 7000 au Protostellar Gas Ring around the Class I source CrA-IRS 2 as a possible sign of magnetic flux advection

Authors: Kazuki Tokuda, Naofumi Fukaya, Kengo Tachihara, Mitsuki Omura, Naoto Harada, Shingo Nozaki, Ayumu Shoshi, Masahiro N. Machida

Abstract: Transferring a significant fraction of the magnetic flux from a dense cloud core is essential in the star formation process. A ring-like structure produced by magnetic flux loss has been predicted theoretically, but no observational identification has been presented. We have performed ALMA observations of the Class I protostar IRS 2 in the Corona Australis star-forming region and resolved a distin… ▽ More Transferring a significant fraction of the magnetic flux from a dense cloud core is essential in the star formation process. A ring-like structure produced by magnetic flux loss has been predicted theoretically, but no observational identification has been presented. We have performed ALMA observations of the Class I protostar IRS 2 in the Corona Australis star-forming region and resolved a distinctive gas ring in the C$^{18}$O ($J$ = 2-1) line emission. The center of this gas ring is $\sim$5,000 au away from the protostar, with a diameter of $\sim$7,000 au. The radial velocity of the gas is $\lesssim1$ km s$^{-1}$ blueshifted from that of the protostar, with a possible expanding feature judged from the velocity-field (moment 1) map and position-velocity diagram. These features are either observationally new or have been discovered but not discussed in depth because they are difficult to explain by well-studied protostellar phenomena such as molecular outflows and accretion streamers. A plausible interpretation is a magnetic wall created by the advection of magnetic flux which is theoretically expected in the Class 0/I phase during star formation as a removal mechanism of magnetic flux. Similar structures reported in the other young stellar sources could likely be candidates formed by the same mechanism, encouraging us to revisit the issue of magnetic flux transport in the early stages of star formation from an observational perspective. △ Less

Submitted 15 October, 2023; v1 submitted 24 September, 2023; originally announced September 2023.

Comments: 9 pages, 3 figures, Accepted for publication in the Astronomical Journal Letters

Journal ref: 10.3847/2041-8213/acfca9

FOREVER22: Gas and metal outflow from massive galaxies in protocluster regions

Authors: Naoki Harada, Hidenobu Yajima, Makito Abe

Abstract: We study gas and metal outflow from massive galaxies in protocluster regions at $z=3-9$ by using the results of the FOREVER22 simulation project. Our simulations contain massive haloes with $M_{\rm h} \gtrsim 10^{13}~\rm M_{\odot}$, showing high star formation rates of $> 100~\rm M_{\odot}~yr^{-1}$ and hosting supermassive black holes with $M_{\rm BH} \gtrsim 10^{8}~\rm M_{\odot}$. We show that th… ▽ More We study gas and metal outflow from massive galaxies in protocluster regions at $z=3-9$ by using the results of the FOREVER22 simulation project. Our simulations contain massive haloes with $M_{\rm h} \gtrsim 10^{13}~\rm M_{\odot}$, showing high star formation rates of $> 100~\rm M_{\odot}~yr^{-1}$ and hosting supermassive black holes with $M_{\rm BH} \gtrsim 10^{8}~\rm M_{\odot}$. We show that the mass loading factor ($η_{\rm M}$) sensitively depends on the halo mass and it is $η_{\rm M} = 1.2~(9.2)$ for $M_{\rm h} = 10^{13}~(10^{11})~\rm M_{\odot}$. Once the halo mass exceeds $\sim 10^{12.5}~\rm M_{\odot}$, the outflow velocity of the gas rapidly decreases near a virial radius, and the gas returns to a galactic centre finally as a fountain flow. Also, the metal inflow and outflow rates sensitively depend on the halo mass and redshift. At $z=3$, the inflow rate becomes larger than the outflow one if $M_{\rm h} \gtrsim 10^{13.0}~\rm M_{\odot}$. Thus, we suggest that massive haloes cannot be efficient metal enrichment sources beyond virial radii that will be probed in future observations, e.g., studies of metal absorption lines with the Prime Focus Spectrograph on the Subaru telescope. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: 13 pages, 10 figures, accepted for publication in MNRAS

The Detection of Higher-Order Millimeter Hydrogen Recombination Lines in the Large Magellanic Cloud

Authors: Marta Sewiło, Kazuki Tokuda, Stan E. Kurtz, Steven B. Charnley, Thomas Möller, Jennifer Wiseman, C. -H. Rosie Chen, Remy Indebetouw, Álvaro Sánchez-Monge, Kei E. I. Tanaka, Peter Schilke, Toshikazu Onishi, Naoto Harada

Abstract: We report the first extragalactic detection of the higher-order millimeter hydrogen recombination lines ($Δn>2$). The $γ$-, $ε$-, and $η$-transitions have been detected toward the millimeter continuum source N105-1A in the star-forming region N105 in the Large Magellanic Cloud (LMC) with the Atacama Large Millimeter/submillimeter Array (ALMA). We use the H40$α$ line, the brightest of the detected… ▽ More We report the first extragalactic detection of the higher-order millimeter hydrogen recombination lines ($Δn>2$). The $γ$-, $ε$-, and $η$-transitions have been detected toward the millimeter continuum source N105-1A in the star-forming region N105 in the Large Magellanic Cloud (LMC) with the Atacama Large Millimeter/submillimeter Array (ALMA). We use the H40$α$ line, the brightest of the detected recombination lines (H40$α$, H36$β$, H50$β$, H41$γ$, H57$γ$, H49$ε$, H53$η$, and H54$η$), and/or the 3 mm free-free continuum emission to determine the physical parameters of N105-1A (the electron temperature, emission measure, electron density, and size) and study ionized gas kinematics. We compare the physical properties of N105-1A to a large sample of Galactic compact and ultracompact (UC) H II regions and conclude that N105-1A is similar to the most luminous ($L>10^5$ $L_{\odot}$) UC H II regions in the Galaxy. N105-1A is ionized by an O5.5 V star, it is deeply embedded in its natal molecular clump, and likely associated with a (proto)cluster. We incorporate high-resolution molecular line data including CS, SO, SO$_2$, and CH$_3$OH ($\sim$0.12 pc), and HCO$^{+}$ and CO ($\sim$0.087 pc) to explore the molecular environment of N105-1A. Based on the CO data, we find evidence for a cloud-cloud collision that likely triggered star formation in the region. We find no clear outflow signatures, but the presence of filaments and streamers indicates on-going accretion onto the clump hosting the UC H II region. Sulfur chemistry in N105-1A is consistent with the accretion shock model predictions. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: 51 pages, 30 figures, 2 tables (including appendices); accepted for publication in The Astrophysical Journal (ApJ)

An ALMA Glimpse of Dense Molecular Filaments Associated with High-mass Protostellar Systems in the Large Magellanic Cloud

Authors: Kazuki Tokuda, Naoto Harada, Kei E. I. Tanaka, Tsuyoshi Inoue, Takashi Shimonishi, Yichen Zhang, Marta Sewiło, Yuri Kunitoshi, Ayu Konishi, Yasuo Fukui, Akiko Kawamura, Toshikazu Onishi, Masahiro N. Machida

Abstract: Recent millimeter/sub-millimeter facilities have revealed the physical properties of filamentary molecular clouds in relation to high-mass star formation. A uniform survey of the nearest, face-on star-forming galaxy, the Large Magellanic Cloud (LMC), complements the Galactic knowledge. We present ALMA survey data with a spatial resolution of $\sim$0.1 pc in the 0.87 mm continuum and HCO$^{+}$(4-3)… ▽ More Recent millimeter/sub-millimeter facilities have revealed the physical properties of filamentary molecular clouds in relation to high-mass star formation. A uniform survey of the nearest, face-on star-forming galaxy, the Large Magellanic Cloud (LMC), complements the Galactic knowledge. We present ALMA survey data with a spatial resolution of $\sim$0.1 pc in the 0.87 mm continuum and HCO$^{+}$(4-3) emission toward 30 protostellar objects with luminosities of 10$^4$-10$^{5.5}$ $L_{\odot}$ in the LMC. The spatial distributions of the HCO$^{+}$(4-3) line and thermal dust emission are well correlated, indicating that the line effectively traces dense, filamentary gas with an H$_2$ volume density of $\gtrsim$10$^5$ cm$^{-3}$ and a line mass of $\sim$10$^3$-10$^{4}$ $M_{\odot}$ pc$^{-1}$. Furthermore, we obtain an increase in the velocity linewidths of filamentary clouds, which follows a power-law dependence on their H$_2$ column densities with an exponent of $\sim$0.5. This trend is consistent with observations toward filamentary clouds in nearby star-forming regions withiin $ \lesssim$1 kpc from us and suggests enhanced internal turbulence within the filaments owing to surrounding gas accretion. Among the 30 sources, we find that 14 are associated with hub-filamentary structures, and these complex structures predominantly appear in protostellar luminosities exceeding $\sim$5 $\times$10$^4$ $L_{\odot}$. The hub-filament systems tend to appear in the latest stages of their natal cloud evolution, often linked to prominent H$\;${\sc ii} regions and numerous stellar clusters. Our preliminary statistics suggest that the massive filaments accompanied by hub-type complex features may be a necessary intermediate product in forming extremely luminous high-mass stellar systems capable of ultimately dispersing the parent cloud. △ Less

Submitted 10 August, 2023; originally announced August 2023.

Comments: 21 pages, 8 figures, 4 tables, accepted for publication in ApJ

Molecular Abundance of the Circumnuclear Region Surrounding an Active Galactic Nucleus in NGC 1068 based on Imaging Line Survey in the 3-mm Band with ALMA

Authors: Taku Nakajima, Shuro Takano, Tomoka Tosaki, Akio Taniguchi, Nanase Harada, Toshiki Saito, Masatoshi Imanishi, Yuri Nishimura, Takuma Izumi, Yoichi Tamura, Kotaro Kohno, Eric Herbst

Abstract: We present an imaging molecular line survey in the 3-mm band (85-114 GHz) focused on one of the nearest galaxies with an active galactic nucleus (AGN), NGC 1068, based on observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA). Distributions of 23 molecular transitions are obtained in the central ~3 kpc region, including both the circumnuclear disk (CND) and starburst ring… ▽ More We present an imaging molecular line survey in the 3-mm band (85-114 GHz) focused on one of the nearest galaxies with an active galactic nucleus (AGN), NGC 1068, based on observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA). Distributions of 23 molecular transitions are obtained in the central ~3 kpc region, including both the circumnuclear disk (CND) and starburst ring (SBR) with 60 and 350 pc resolution. The column densities and relative abundances of all the detected molecules are estimated under the assumption of local thermodynamic equilibrium in the CND and SBR. Then, we discuss the physical and chemical effects of the AGN on molecular abundance corresponding to the observation scale. We found that H13CN, SiO, HCN, and H13CO+ are abundant in the CND relative to the SBR. In contrast, 13CO is more abundant in the SBR. Based on the calculated column density ratios of N(HCN)/N(HCO+), N(HCN)/N(CN), and other molecular distributions, we conclude that the enhancement of HCN in the CND may be due to high-temperature environments resulting from strong shocks, which are traced by the SiO emission. Moreover, the abundance of CN in the CND is significantly lower than the expected value of the model calculations in the region affected by strong radiation. The expected strong X-ray irradiation from the AGN has a relatively lower impact on the molecular abundance in the CND than mechanical feedback. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: 33 pages, 20 figures, 4 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

Crescent-Shaped Molecular Outflow from the Intermediate-mass Protostar DK Cha Revealed by ALMA

Authors: Naoto Harada, Kazuki Tokuda, Hayao Yamasaki, Asako Sato, Mitsuki Omura, Shingo Hirano, Toshikazu Onishi, Kengo Tachihara, Masahiro N. Machida

Abstract: We report on an Atacama Large Millimeter/submillimeter Array (ALMA) study of the Class I or II intermediate-mass protostar DK Cha in the Chamaeleon II region. The 12CO (J=2-1) images have an angular resolution of ~1'' (~250 au) and show high-velocity blueshifted (>70 km s-1) and redshifted (>50 km s-1) emissions which have 3000 au scale crescent-shaped structures around the protostellar disk trace… ▽ More We report on an Atacama Large Millimeter/submillimeter Array (ALMA) study of the Class I or II intermediate-mass protostar DK Cha in the Chamaeleon II region. The 12CO (J=2-1) images have an angular resolution of ~1'' (~250 au) and show high-velocity blueshifted (>70 km s-1) and redshifted (>50 km s-1) emissions which have 3000 au scale crescent-shaped structures around the protostellar disk traced in the 1.3mm continuum. Because the high-velocity components of the CO emission are associated with the protostar, we concluded that the emission traces the pole-on outflow. The blueshifted outflow lobe has a clear layered velocity gradient with a higher velocity component located on the inner side of the crescent shape, which can be explained by a model of an outflow with a higher velocity in the inner radii. Based on the directly driven outflow scenario, we estimated the driving radii from the observed outflow velocities and found that the driving region extends over two orders of magnitude. The 13CO emission traces a complex envelope structure with arc-like substructures with lengths of ~1000au. We identified the arc-like structures as streamers because they appear to be connected to a rotating infalling envelope. DK Cha is useful for understanding characteristics that are visible by looking at nearly face-on configurations of young protostellar systems, providing an alternative perspective for studying the star-formation process. △ Less

Submitted 3 February, 2023; originally announced February 2023.

Comments: Accepted for publication in ApJ. 12 pages, 5 figures

Twisted magnetic field in star formation processes of L1521 F revealed by submillimeter dual band polarimetry using James Clerk Maxwell Telescope

Authors: Sakiko Fukaya, Hiroko Shinnaga, Ray S. Furuya, Kohji Tomisaka, Masahiro N. Machida, Naoto Harada

Abstract: Understanding the initial conditions of star formation requires both observational studies and theoretical works taking into account the magnetic field, which plays an important role in star formation processes. Herein, we study the young nearby dense cloud core L1521 F ($n$(H$_2$) $\sim 10^{4-6}$ cm$^{-3}$) in the Taurus Molecular Cloud. This dense core hosts a 0.2 $M_\odot$ protostar, categorize… ▽ More Understanding the initial conditions of star formation requires both observational studies and theoretical works taking into account the magnetic field, which plays an important role in star formation processes. Herein, we study the young nearby dense cloud core L1521 F ($n$(H$_2$) $\sim 10^{4-6}$ cm$^{-3}$) in the Taurus Molecular Cloud. This dense core hosts a 0.2 $M_\odot$ protostar, categorized as a Very Low Luminosity Objects with complex velocity structures, particularly in the vicinity of the protostar. To trace the magnetic field within the dense core, we conducted high sensitivity submillimeter polarimetry of the dust continuum at $λ$= 850 $μ$m and 450 $μ$m using the POL-2 polarimeter situated in front of the SCUBA-2 submillimeter bolometer camera on James Clerk Maxwell Tetescope. This was compared with millimeter polarimetry taken at $λ$= 3.3 mm with ALMA. The magnetic field was detected at $λ$= 850 $μ$m in the peripheral region, which is threaded in a north-south direction, while the central region traced at $λ$= 450 $μ$m shows a magnetic field with an east-west direction, i.e., orthogonal to that of the peripheral region. Magnetic field strengths are estimated to be $\sim$70 $μ$G and 200 $μ$G in the peripheral- and central-regions, respectively, using the Davis-Chandrasekhar-Fermi method. The resulting mass-to-flux ratio of 3 times larger than that of magnetically critical state for both regions indicates that L1521 F is magnetically supercritical, i.e., gravitational forces dominate over magnetic turbulence forces. Combining observational data with MHD simulations, detailed parameters of the morphological properties of this puzzling object are derived for the first time. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: 9 pages, 7 figures

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

The First Detection of a Protostellar CO Outflow in the Small Magellanic Cloud with ALMA

Authors: Kazuki Tokuda, Sarolta Zahorecz, Yuri Kunitoshi, Kosuke Higashino, Kei E. I. Tanaka, Ayu Konishi, Taisei Suzuki, Naoya Kitano, Naoto Harada, Takashi Shimonishi, Naslim Neelamkodan, Yasuo Fukui, Akiko Kawamura, Toshikazu Onishi, Masahiro N. Machida

Abstract: Protostellar outflows are one of the most outstanding features of star formation. Observational studies over the last several decades have successfully demonstrated that outflows are ubiquitously associated with low- and high-mass protostars in the solar-metallicity Galactic conditions. However, the environmental dependence of protostellar outflow properties is still poorly understood, particularl… ▽ More Protostellar outflows are one of the most outstanding features of star formation. Observational studies over the last several decades have successfully demonstrated that outflows are ubiquitously associated with low- and high-mass protostars in the solar-metallicity Galactic conditions. However, the environmental dependence of protostellar outflow properties is still poorly understood, particularly in the low-metallicity regime. Here we report the first detection of a molecular outflow in the Small Magellanic Cloud with 0.2 $Z_{\odot}$, using Atacama Large Millimeter/submillimeter Array observations at a spatial resolution of 0.1 pc toward the massive protostar Y246. The bipolar outflow is nicely illustrated by high-velocity wings of CO(3-2) emission at $\gtrsim$15 km s$^{-1}$. The evaluated properties of the outflow (momentum, mechanical force, etc.) are consistent with those of the Galactic counterparts. Our results suggest that the molecular outflows, i.e., the guidepost of the disk accretion at the small scale, might be universally associated with protostars across the metallicity range of $\sim$0.2-1 $Z_{\odot}$. △ Less

Submitted 7 August, 2022; v1 submitted 18 July, 2022; originally announced July 2022.

Comments: 7 pages, 2 figures, Accepted for publication in ApJL

AGN-driven Cold Gas Outflow of NGC 1068 Characterized by Dissociation-Sensitive Molecules

Authors: Toshiki Saito, Shuro Takano, Nanase Harada, Taku Nakajima, Eva Schinnerer, Daizhong Liu, Akio Taniguchi, Takuma Izumi, Yumi Watanabe, Kazuharu Bamba, Kotaro Kohno, Yuri Nishimura, Sophia Stuber, Tomoka Tosaki

Abstract: Recent developments in (sub-)millimeter facilities have drastically changed the amount of information obtained from extragalactic spectral scans. In this paper, we present a feature extraction technique using principal component analysis (PCA) applied to arcsecond-resolution (1.0-2.0 arcsec = 72-144 pc) spectral scan datasets for the nearby type-2 Seyfert galaxy, NGC 1068, using Band 3 of the Atac… ▽ More Recent developments in (sub-)millimeter facilities have drastically changed the amount of information obtained from extragalactic spectral scans. In this paper, we present a feature extraction technique using principal component analysis (PCA) applied to arcsecond-resolution (1.0-2.0 arcsec = 72-144 pc) spectral scan datasets for the nearby type-2 Seyfert galaxy, NGC 1068, using Band 3 of the Atacama Large Millimeter/submillimeter Array. We apply PCA to 16 well-detected molecular line intensity maps convolved to a common 150 pc resolution. In addition, we include the [SIII]/[SII] line ratio and [CI] $^3P_1$-$^3P_0$ maps in the literature, both of whose distributions show remarkable resemblance with that of a kpc-scale biconical outflow from the central AGN. We identify two prominent features: (1) central concentration at the circumnuclear disk (CND) and (2) two peaks across the center that coincide with the biconical outflow peaks. The concentrated molecular lines in the CND are mostly high-dipole molecules (e.g., H$^{13}$CN, HC$_3$N, and HCN). Line emissions from molecules known to be enhanced in irradiated interstellar medium, CN, C$_2$H, and HNC, show similar concentrations and extended components along the bicone, suggesting that molecule dissociation is a dominant chemical effect of the cold molecular outflow of this galaxy. Although further investigation should be made, this scenario is consistent with the faintness or absence of the emission lines from CO isotopologues, CH$_3$OH, and N$_2$H$^+$, in the outflow, which are easily destroyed by dissociating photons and electrons. △ Less

Submitted 13 July, 2022; originally announced July 2022.

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

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

APEX and NOEMA observations of H$_{2}$S in nearby luminous galaxies and the ULIRG Mrk~231 -- Is there a relation between dense gas properties and molecular outflows?

Authors: M. T. Sato, S. Aalto, K. Kohno, S. König, N. Harada, S. Viti, T. Izumi, Y. Nishimura, M. Gorski

Abstract: In order to understand the evolution and feedback of Active Galactic Nuclei (AGN) and star formation it is important to use molecular lines as probes of physical conditions and chemistry. We use H$_{2}$S to investigate the impact of starburst and AGN activity on the chemistry of the molecular interstellar medium in luminous infrared galaxies. Using the APEX single dish telescope, we have observed… ▽ More In order to understand the evolution and feedback of Active Galactic Nuclei (AGN) and star formation it is important to use molecular lines as probes of physical conditions and chemistry. We use H$_{2}$S to investigate the impact of starburst and AGN activity on the chemistry of the molecular interstellar medium in luminous infrared galaxies. Using the APEX single dish telescope, we have observed the $1_{10}$--$1_{01}$ transition of ortho-H$_{2}$S at 168 GHz towards the centres of twelve nearby luminous infrared galaxies. We have also observed the same line towards the ultra luminous infrared galaxy (ULIRG) Mrk~231 with the NOEMA interferometer. We have detected H$_{2}$S towards NGC~253, NGC~1068, NGC~3256, NGC~4418, NGC~4826, NGC~4945, Circinus, M~83 and Mrk~231. Four galaxies show elevated H$_{2}$S emission relative to HCN. We suggest that the high line ratios are caused by elevated H$_{2}$S abundances in the dense gas. However, we do not find any clear connection between the H$_{2}$S/HCN line intensity ratio, and the presence (or speed) of molecular outflows in the sample galaxies. Therefore H$_{2}$S abundances do not seem to be globally affected by the large-scale outflows. We discuss possible mechanisms behind the suggested H$_{2}$S abundance enhancements in NGC~4418, Circinus, NGC~3256 and NGC~4826. These include radiative processes (for example X-rays or cosmic-rays) or smaller scale shocks. We suggest that $L_{\mathrm{H_{2}S}}$ serves as a tracer of the dense gas content, similar to $L_{\mathrm{HCN}}$, and that the correlation between $L_{\mathrm{H_{2}S}}$ and $M_{\rm outflow}$(H$_2$) implies a relation between the dense gas reservoir and the properties and evolution of the molecular feedback. This potential link requires further study since it holds important keys to our understanding of how the properties of molecular outflows relate to that of their host galaxies. △ Less

Submitted 11 March, 2022; originally announced March 2022.

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

The Kiloparsec-scale Neutral Atomic Carbon Outflow in the Nearby Type-2 Seyfert Galaxy NGC 1068: Evidence for Negative AGN Feedback

Authors: Toshiki Saito, Shuro Takano, Nanase Harada, Taku Nakajima, Eva Schinnerer, Daizhong Liu, Akio Taniguchi, Takuma Izumi, Yumi Watanabe, Kazuharu Bamba, Eric Herbst, Kotaro Kohno, Yuri Nishimura, Sophia Stuber, Yoichi Tamura, Tomoka Tosaki

Abstract: Active galactic nucleus (AGN) feedback is postulated as a key mechanism for regulating star formation within galaxies. Studying the physical properties of the outflowing gas from AGN is thus crucial for understanding the co-evolution of galaxies and supermassive black holes. Here we report 55 pc resolution ALMA neutral atomic carbon [CI] $^3P_1\text{-}^3P_0$ observations toward the central 1 kpc o… ▽ More Active galactic nucleus (AGN) feedback is postulated as a key mechanism for regulating star formation within galaxies. Studying the physical properties of the outflowing gas from AGN is thus crucial for understanding the co-evolution of galaxies and supermassive black holes. Here we report 55 pc resolution ALMA neutral atomic carbon [CI] $^3P_1\text{-}^3P_0$ observations toward the central 1 kpc of the nearby type-2 Seyfert galaxy NGC 1068, supplemented by 55 pc resolution CO($J=1\text{-}0$) observations. We find that [CI] emission within the central kpc is strongly enhanced by a factor of $>$5 compared to the typical [CI]/CO intensity ratio of $\sim$0.2 for nearby starburst galaxies (in units of brightness temperature). The most [CI]-enhanced gas (ratio $>$ 1) exhibits a kpc-scale elongated structure centered at the AGN that matches the known biconical ionized gas outflow entraining molecular gas in the disk. A truncated, decelerating bicone model explains well the kinematics of the elongated structure, indicating that the [CI] enhancement is predominantly driven by the interaction between the ISM in the disk and the highly inclined ionized gas outflow (which is likely driven by the radio jet). Our results strongly favor the "CO dissociation scenario" rather than the "in-situ C formation" one which prefers a perfect bicone geometry. We suggest that the high [CI]/CO intensity ratio gas in NGC 1068 directly traces ISM in the disk that is currently dissociated and entrained by the jet and the outflow, i.e., the "negative" effect of the AGN feedback. △ Less

Submitted 2 March, 2022; originally announced March 2022.

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

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)

Impact of Magnetic Braking on High-mass Close Binary Formation

Authors: Naoto Harada, Shingo Hirano, Masahiro N. Machida, Takashi Hosokawa

Abstract: Combining numerical simulations and analytical modeling, we investigate whether close binary systems form by the effect of magnetic braking. Using magnetohydrodynamics simulations, we calculate the cloud evolution with a sink, for which we do not resolve the binary system or binary orbital motion to realize long-term time integration. Then, we analytically estimate the binary separation using the… ▽ More Combining numerical simulations and analytical modeling, we investigate whether close binary systems form by the effect of magnetic braking. Using magnetohydrodynamics simulations, we calculate the cloud evolution with a sink, for which we do not resolve the binary system or binary orbital motion to realize long-term time integration. Then, we analytically estimate the binary separation using the accreted mass and angular momentum obtained from the simulation. In unmagnetized clouds, wide binary systems with separations of >100 au form, in which the binary separation continues to increase during the main accretion phase. In contrast, close binary systems with separations of <100 au can form in magnetized clouds. Since the efficiency of magnetic braking strongly depends on both the strength and configuration of the magnetic field, they also affect the formation conditions of a close binary. In addition, the protostellar outflow has a negative impact on close binary formation, especially when the rotation axis of the prestellar cloud is aligned with the global magnetic field. The outflow interrupts the accretion of gas with small angular momentum, which is expelled from the cloud, while gas with large angular momentum preferentially falls from the side of the outflow onto the binary system and widens the binary separation. This study shows that a cloud with a magnetic field that is not parallel to the rotation axis is a favorable environment for the formation of close binary systems. △ Less

Submitted 24 September, 2021; originally announced September 2021.

Comments: Accepted for publication in MNRAS

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)

The Circumnuclear Disk Revealed by ALMA. I. Dense Clouds and Tides in the Galactic Center

Authors: Pei-Ying Hsieh, Patrick M. Koch, Woong-Tae Kim, Sergio Martín, Hsi-Wei Yen, John Carpenter, Nanase Harada, Jean L. Turner, Paul T. P. Ho, Ya-Wen Tang, Sara C. Beck

Abstract: Utilizing the Atacama Large Millimeter/submillimeter Array (ALMA), we present CS line maps in five rotational lines ($J_{\rm u}=7, 5, 4, 3, 2$) toward the circumnuclear disk (CND) and streamers of the Galactic Center. Our primary goal is to resolve the compact structures within the CND and the streamers, in order to understand the stability conditions of molecular cores in the vicinity of the supe… ▽ More Utilizing the Atacama Large Millimeter/submillimeter Array (ALMA), we present CS line maps in five rotational lines ($J_{\rm u}=7, 5, 4, 3, 2$) toward the circumnuclear disk (CND) and streamers of the Galactic Center. Our primary goal is to resolve the compact structures within the CND and the streamers, in order to understand the stability conditions of molecular cores in the vicinity of the supermassive black hole (SMBH) Sgr A*. Our data provide the first homogeneous high-resolution ($1.3'' = 0.05$ pc) observations aiming at resolving density and temperature structures. The CS clouds have sizes of $0.05-0.2$ pc with a broad range of velocity dispersion ($σ_{\rm FWHM}=5-40$ km s$^{-1}$). The CS clouds are a mixture of warm ($T_{\rm k}\ge 50-500$ K, n$_{\rm H_2}$=$10^{3-5}$ cm$^{-3}$) and cold gas ($T_{\rm k}\le 50$ K, n$_{\rm H_2}$=$10^{6-8}$ cm$^{-3}$). A stability analysis based on the unmagnetized virial theorem including tidal force shows that $84^{+16}_{-37}$ % of the total gas mass is tidally stable, which accounts for the majority of gas mass. Turbulence dominates the internal energy and thereby sets the threshold densities $10-100$ times higher than the tidal limit at distance $\ge 1.5$ pc to Sgr A*, and therefore, inhibits the clouds from collapsing to form stars near the SMBH. However, within the central $1.5$ pc, the tidal force overrides turbulence and the threshold densities for a gravitational collapse quickly grow to $\ge 10^{8}$ cm$^{-3}$. △ Less

Submitted 5 April, 2021; originally announced April 2021.

Comments: accepted for publication in ApJ

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 observations of CS in NGC 1068: chemistry and excitation

Authors: M. Scourfield, S. Viti, S. Garcia-Burillo, A. Saintonge, F. Combes, A. Fuente, C. Henkel, A. Alonso-Herrero, N. Harada, S. Takano, T. Nakajima, S. Martin, M. Krips, P. P. van der Werf, S. Aalto, A. Usero, K. Kohno

Abstract: We present results from Atacama Large Millimeter/submillimeter Array (ALMA) observations of CS from the nearby galaxy NGC 1068 ($\sim14$ Mpc). This Seyfert 2 barred galaxy possesses a circumnuclear disc (CND, $r\sim200$ pc) and a starburst ring (SB ring, $r\sim1.3$ kpc). These high-resolution maps ($\sim0.5$", $\sim35$ pc) allow us to analyse specific sub-regions in the galaxy and investigate diff… ▽ More We present results from Atacama Large Millimeter/submillimeter Array (ALMA) observations of CS from the nearby galaxy NGC 1068 ($\sim14$ Mpc). This Seyfert 2 barred galaxy possesses a circumnuclear disc (CND, $r\sim200$ pc) and a starburst ring (SB ring, $r\sim1.3$ kpc). These high-resolution maps ($\sim0.5$", $\sim35$ pc) allow us to analyse specific sub-regions in the galaxy and investigate differences in line intensity ratios and physical conditions, particularly those between the CND and SB ring. Local thermodynamic equilibrium (LTE) analysis of the gas is used to calculate CS densities in each sub-region, followed by non-LTE analysis conducted using the radiative transfer code RADEX to fit observations and constrain gas temperature, CS column density and hydrogen density. Finally, the chemical code UCLCHEM is used to reconstruct the gas, allowing an insight into its origin and chemical history. The density of hydrogen in the CND is found to be $\geq10^5$ cm$^{-2}$, although exact values vary, reaching $10^6$ cm$^{-2}$ at the AGN. The conditions in the two arms of the SB ring appear similar to one another, though the density found ($\sim10^4$ cm$^{-2}$) is lower than in the CND. The temperature in the CND increases from east to west, and is also overall greater than found in the SB ring. These modelling methods indicate the requirement for multi-phase gas components in order to fit the observed emission over the galaxy. A larger number of high resolution transitions across the SLED may allow for further constraining of the conditions, particularly in the SB ring. △ Less

Submitted 26 June, 2020; originally announced June 2020.

Comments: 32 pages, 19 figures, 6 tables, accepted for publication in MNRAS

ALMA Observations of Multiple-CO and C Lines Toward the Active Galactic Nucleus of NGC 7469: X-Ray-dominated Region Caught in the Act

Authors: Takuma Izumi, Dieu D. Nguyen, Masatoshi Imanishi, Taiki Kawamuro, Shunsuke Baba, Suzuka Nakano, Kotaro Kohno, Satoki Matsushita, David S. Meier, Jean L. Turner, Tomonari Michiyama, Nanase Harada, Sergio Martín, Kouichiro Nakanishi, Shuro Takano, Tommy Wiklind, Naomasa Nakai, Pei-Ying Hsieh

Abstract: We used the Atacama Large Millimeter/submillimeter Array (ALMA) to map $^{12}$CO($J$ = 1-0), $^{12}$CO($J$ = 2-1), $^{12}$CO($J$ = 3-2), $^{13}$CO($J$ = 2-1), and [CI]($^3P_1$-$^3P_0$) emission lines around the type 1 active galactic nucleus (AGN) of NGC 7469 ($z = 0.0164$) at $\sim 100$ pc resolutions. The CO lines are bright both in the circumnuclear disk (central $\sim 300$ pc) and the surround… ▽ More We used the Atacama Large Millimeter/submillimeter Array (ALMA) to map $^{12}$CO($J$ = 1-0), $^{12}$CO($J$ = 2-1), $^{12}$CO($J$ = 3-2), $^{13}$CO($J$ = 2-1), and [CI]($^3P_1$-$^3P_0$) emission lines around the type 1 active galactic nucleus (AGN) of NGC 7469 ($z = 0.0164$) at $\sim 100$ pc resolutions. The CO lines are bright both in the circumnuclear disk (central $\sim 300$ pc) and the surrounding starburst (SB) ring ($\sim 1$ kpc diameter), with two bright peaks on either side of the AGN. By contrast, the [CI]($^3P_1$-$^3P_0$) line is strongly peaked on the AGN. Consequently, the brightness temperature ratio of [CI]($^3P_1$-$^3P_0$) to $^{13}$CO(2-1) is $\sim 20$ at the AGN, as compared to $\sim 2$ in the SB ring. Our local thermodynamic equilibrium (LTE) and non-LTE models indicate that the enhanced line ratios (or CI enhancement) are due to an elevated C$^0$/CO abundance ratio ($\sim 3-10$) and temperature ($\sim 100-500$ K) around the AGN as compared to the SB ring (abundance ratio $\sim 1$, temperature $\lesssim 100$ K), which accords with the picture of the X-ray-dominated Region (XDR). Based on dynamical modelings, we also provide CO(1-0)-to- and [CI]($^3P_1$-$^3P_0$)-to-molecular mass conversion factors at the central $\sim 100$ pc of this AGN as $α_{\rm CO} = 4.1$ and $α_{\rm CI} = 4.4~M_\odot$ (K km s$^{-1}$ pc$^2$)$^{-1}$, respectively. Our results suggest that the CI enhancement is potentially a good marker of AGNs that could be used in a new submillimeter diagnostic method toward dusty environments. △ Less

Submitted 3 July, 2020; v1 submitted 16 June, 2020; originally announced June 2020.

Comments: 25 pages, 12 figures, 8 tables. Accepted for publication in ApJ. Minor updates in a replacement

On the effects of UV photons/X-rays on the chemistry of the Sgr B2 cloud

Authors: J. Armijos-Abendaño, J. Martín-Pintado, E. López, M. Llerena, N. Harada, M. A. Requena-Torres, S. Martín, V. M. Rivilla, D. Riquelme, F. Aldas

Abstract: The lines of HOC$^+$, HCO and CO$^+$ are considered good tracers of photon-dominated regions (PDRs) and X-ray dominated regions. We study these tracers towards regions of the Sgr B2 cloud selected to be affected by different heating mechanisms. We find the lowest values of the column density ratios of HCO$^+$ versus HOC$^+$, HCO and CO$^+$ in dense HII gas, where UV photons dominate the heating an… ▽ More The lines of HOC$^+$, HCO and CO$^+$ are considered good tracers of photon-dominated regions (PDRs) and X-ray dominated regions. We study these tracers towards regions of the Sgr B2 cloud selected to be affected by different heating mechanisms. We find the lowest values of the column density ratios of HCO$^+$ versus HOC$^+$, HCO and CO$^+$ in dense HII gas, where UV photons dominate the heating and chemistry of gas. HOC$^+$, HCO and CO$^+$ abundances and the above ratios are compared with those of chemical modeling, finding that high temperature chemistry, a cosmic-ray ionization rate of 10$^{-16}$ s$^{-1}$ and timescales $>$10$^{5.0}$ years explain well the HOC$^+$ abundances in quiescent Sgr B2 regions, while shocks are also needed to explain the highest HCO abundances derived for these regions. CO$^+$ is mainly formed in PDRs since the highest CO$^+$ abundances of $\sim$(6-10)$\times$10$^{-10}$ are found in HII regions with electron densities $>$540 cm$^{-3}$ and that CO$^+$ emission is undetected in quiescent gas. Between the ratios, the HCO$^+$/HCO ratio is sensitive to the electron density as it shows different values in dense and diffuse HII regions. We compare SiO J=2-1 emission maps of Sgr B2 with X-ray maps from 2004 and 2012. One known spot shown on the 2012 X-ray map is likely associated with molecular gas at velocities of 15-25 km s$^{-1}$. We also derive the X-ray ionization rate of $\sim$10$^{-19}$ s$^{-1}$ for Sgr B2 regions pervaded by X-rays in 2004, which is quite low to affect the chemistry of the molecular gas. △ Less

Submitted 27 April, 2020; originally announced April 2020.

The MALATANG Survey: Dense Gas and Star Formation from High Transition HCN and HCO+ maps of NGC253

Authors: Xue-Jian Jiang, Thomas R. Greve, Yu Gao, Zhi-Yu Zhang, Qinghua Tan, Richard de Grijs, Luis C. Ho, Michal J. Michalowski, Malcolm J. Currie, Christine D. Wilson, Elias Brinks, Yiping Ao, Yinghe Zhao, Jinhua He, Nanase Harada, Chentao Yang, Qian Jiao, Aeree Chung, Bumhyun Lee, Matthew W. L. Smith, Daizhong Liu, Satoki Matsushita, Yong Shi, Masatoshi Imanishi, Mark G. Rawlings , et al. (4 additional authors not shown)

Abstract: To study the high-transition dense-gas tracers and their relationships to the star formation of the inner $\sim$ 2 kpc circumnuclear region of NGC253, we present HCN $J=4-3$ and HCO$^+ J=4-3$ maps obtained with the James Clerk Maxwell Telescope (JCMT). With the spatially resolved data, we compute the concentration indices $r_{90}/r_{50}$ for the different tracers. HCN and HCO$^+$ 4-3 emission feat… ▽ More To study the high-transition dense-gas tracers and their relationships to the star formation of the inner $\sim$ 2 kpc circumnuclear region of NGC253, we present HCN $J=4-3$ and HCO$^+ J=4-3$ maps obtained with the James Clerk Maxwell Telescope (JCMT). With the spatially resolved data, we compute the concentration indices $r_{90}/r_{50}$ for the different tracers. HCN and HCO$^+$ 4-3 emission features tend to be centrally concentrated, which is in contrast to the shallower distribution of CO 1-0 and the stellar component. The dense-gas fraction ($f_\text{dense}$, traced by the velocity-integrated-intensity ratios of HCN/CO and HCO$^+$/CO) and the ratio $R_\text{31}$ (CO 3-2/1-0) decline towards larger galactocentric distances, but increase with higher SFR surface density. The radial variation and the large scatter of $f_\text{dense}$ and $R_\text{31}$ imply distinct physical conditions in different regions of the galactic disc. The relationships of $f_\text{dense}$ versus $Σ_\text{stellar}$, and SFE$_\text{dense}$ versus $Σ_\text{stellar}$ are explored. SFE$_\text{dense}$ increases with higher $Σ_\text{stellar}$ in this galaxy, which is inconsistent with previous work that used HCN 1-0 data. This implies that existing stellar components might have different effects on the high-$J$ HCN and HCO$^+$ than their low-$J$ emission. We also find that SFE$_\text{dense}$ seems to be decreasing with higher $f_\text{dense}$, which is consistent with previous works, and it suggests that the ability of the dense gas to form stars diminishes when the average density of the gas increases. This is expected in a scenario where only the regions with high-density contrast collapse and form stars. △ Less

Submitted 14 March, 2020; originally announced March 2020.

Comments: accepted to MNRAS

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

The Nuclear Filaments inside the Circumnuclear Disk in the Central 0.5 pc of the Galactic center

Authors: Pei-Ying Hsieh, Patrick M. Koch, Woong-Tae Kim, Paul T. P. Ho, Hsi-Wei Yen, Nanase Harada, Ya-Wen Tang

Abstract: We present CS(7-6) line maps toward the central parsec of the Galactic Center (GC), conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The primary goal is to find and characterize the gas structure in the inner cavity of the circumnuclear disk (CND) in high resolution (1.3"=0.05 pc). Our large field-of-view mosaic maps -- combining interferometric and single-dish data that rec… ▽ More We present CS(7-6) line maps toward the central parsec of the Galactic Center (GC), conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The primary goal is to find and characterize the gas structure in the inner cavity of the circumnuclear disk (CND) in high resolution (1.3"=0.05 pc). Our large field-of-view mosaic maps -- combining interferometric and single-dish data that recover extended emission - provide a first homogeneous look to resolve and link the molecular streamers in the CND with the neutral nuclear filaments newly detected within the central cavity of the CND. We find that the nuclear filaments are rotating with Keplerian velocities in a nearly face-on orbit with an inclination angle of ~10-20 degree (radius <= 0.5 pc). This is in contrast to the CND which is highly inclined at ~65-80 degree (radius ~2-5 pc). Our analysis suggests a highly warped structure from the CND to the nuclear filaments. This result may hint that the nuclear filaments and the CND were created by different external clouds passing by Sgr A*. △ Less

Submitted 16 October, 2019; originally announced October 2019.

Comments: 6 pages, 4 figures, accepted in ApJL

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

Molecular Gas Density Measured with H$_2$CO and CS toward a Spiral Arm of M51

Authors: Yuri Nishimura, Yoshimasa Watanabe, Nanase Harada, Kotaro Kohno, Satoshi Yamamoto

Abstract: Observations of various molecular lines toward a disk region of a nearby galaxy are now feasible, and they are being employed as diagnostic tools to study star-formation activities there. However, the spatial resolution attainable for a nearby galaxy with currently available radio telescopes is $10-1000$ pc, which is much larger than the scales of individual star-forming regions and molecular clou… ▽ More Observations of various molecular lines toward a disk region of a nearby galaxy are now feasible, and they are being employed as diagnostic tools to study star-formation activities there. However, the spatial resolution attainable for a nearby galaxy with currently available radio telescopes is $10-1000$ pc, which is much larger than the scales of individual star-forming regions and molecular cloud cores. Hence, it is of fundamental importance to elucidate which part of an interstellar medium such spatially-unresolved observations are tracing. Here we present sensitive measurements of the H$_2$CO ($1_{01}-0_{00}$) line at 72 GHz toward giant molecular clouds (GMCs) in the spiral arm of M51 using the NRO 45 m and IRAM 30 m telescopes. In conjunction with the previously observed H$_2$CO ($2_{02}-1_{01}$) and CS ($2-1$ and $3-2$) lines, we derive the H$_2$ density of the emitting regions to be $(0.6-2.6)\times10^4$ cm$^{-3}$ and $(2.9-12)\times10^4$ cm$^{-3}$ for H$_2$CO and CS, respectively, by the non-LTE analyses, where we assume the source size of $0.8-1$ kpc and the gas kinetic temperature of $10-20$ K. The derived H$_2$ density indicates that the emission of H$_2$CO and CS is not localized to star-forming cores, but is likely distributed over an entire region of GMCs. Such widespread distributions of H$_2$CO and CS are also supported by models assuming lognormal density distributions over the 1 kpc region. Thus, contributions from the widespread less-dense components should be taken into account for interpretation of these molecular emission observed with a GMC-scale resolution. The different H$_2$ densities derived for H$_2$CO and CS imply their different distributions. We discuss this differences in terms of the formation processes of H$_2$CO and CS. △ Less

Submitted 27 May, 2019; originally announced May 2019.

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

Molecular-Cloud-Scale Chemical Composition III: Constraints of Average Physical Properties through Chemical Models

Authors: Nanase Harada, Yuri Nishimura, Yoshimasa Watanabe, Satoshi Yamamoto, Yuri Aikawa, Nami Sakai, Takashi Shimonishi

Abstract: It is important to understand the origin of molecular line intensities and chemical composition in the molecular-cloud scale in the Galactic sources because it serves as a benchmark to compare with the chemical compositions of extragalactic sources. Recent observations of the 3-mm spectra averaged over the 10-pc scale show similar spectral pattern among sources for molecular lines HCN, HCO$^+$, CC… ▽ More It is important to understand the origin of molecular line intensities and chemical composition in the molecular-cloud scale in the Galactic sources because it serves as a benchmark to compare with the chemical compositions of extragalactic sources. Recent observations of the 3-mm spectra averaged over the 10-pc scale show similar spectral pattern among sources for molecular lines HCN, HCO$^+$, CCH, HNC, HNCO, c-C$_3$H$_2$, CS, SO, N$_2$H$^+$, and CN. To constrain the average physical property emitting such spectral pattern, we model molecular spectra using a time-dependent gas-grain chemical model followed by a radiative transfer calculation. We use a grid of physical parameters such as the density $n=3 \times 10^2 - 3\times 10^4$ cm$^{-3}$, the temperature, $T=10-30$ K, the visual extinction $A_{\rm V} = 2,4,10$ mag, the cosmic-ray ionization rate $ζ= 10^{-17} - 10^{-16}$ s$^{-1}$, and the sulfur elemental abundance $S/H = 8\times 10^{-8} - 8\times 10^{-7}$. Comparison with the observed spectra indicates that spectra are well reproduced with the relatively low density of $n=(1-3) \times 10^3\,$cm$^{-3}$, $T=10\,$K, $ζ= 10^{-17}$ s$^{-1}$, and the short chemistry timescale of $10^5$ yrs. This short chemistry timescale may indicate that molecular clouds are constantly affected by the turbulence, and exposed to low-density, low $A_{\rm V}$ regions that "refreshes" the chemical clock by UV radiation. The relatively low density obtained is orders of magnitude lower than the commonly-quoted critical density in the optically thin case. Meanwhile, this range of density is consistent with results from recent observational analysis of molecular-cloud-scale mapping. △ Less

Submitted 8 December, 2018; originally announced December 2018.

Comments: 30 pages, 21 figures, accepted for publication in The Astrophysical Journal

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