Showing 1–50 of 381 results for author: Henkel, C

Hyperfine structure of the methanol molecule as traced by Class I methanol masers

Authors: I. I. Agafonova, O. S. Bayandina, Y. Gong, C. Henkel, Kee-Tae Kim, M. G. Kozlov, B. Lankhaar, S. A. Levshakov, K. M. Menten, W. Ubachs, I. E. Val'tts, W. Yang

Abstract: We present results on simultaneous observations of Class~I methanol masers at 25, 36, and 44 GHz towards 22 Galactic targets carried out with the Effelsberg 100-m telescope. The study investigates relations between the hyperfine (HF) structure of the torsion-rotation transitions in CH3OH and maser activity. By analyzing the radial velocity shifts between different maser lines together with the pat… ▽ More We present results on simultaneous observations of Class~I methanol masers at 25, 36, and 44 GHz towards 22 Galactic targets carried out with the Effelsberg 100-m telescope. The study investigates relations between the hyperfine (HF) structure of the torsion-rotation transitions in CH3OH and maser activity. By analyzing the radial velocity shifts between different maser lines together with the patterns of the HF structure based on laboratory measurements and quantum-chemical calculations, we find that in any source only one specific HF transition forms the maser emission and that this transition changes from source to source. The physical conditions leading to this selective behavior are still unclear. Using accurate laboratory rest frequencies for the 25 GHz transitions, we have refined the centre frequencies for the HF multiplets at 36, 44, and 95 GHz: f_36 = (36169.2488 +/- 0.0002_stat +/- 0.0004_sys) MHz. f_44 = (44069.4176 +/- 0.0002_stat +/- 0.0004_sys) MHz, and f_95 = (95169.4414 +/- 0.0003_stat +/- 0.0004_sys) MHz. Comparison with previous observations of 44 GHz masers performed 6-10 years ago with a Korean 21-m KVN telescope towards the same targets confirms the kinematic stability of Class~I maser line profiles during this time interval and reveals a systematic radial velocity shift of 0.013 +/- 0.005 km/s between the two telescopes. △ Less

Submitted 16 July, 2024; originally announced July 2024.

Comments: 24 pages, 8 figures, 9 tables; accepted for publication in MNRAS

First detection of the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ methanol maser transitions at $J=7$ and $10$

Authors: Pedro K. Humire, Gisela N. Ortiz-León, Antonio Hernández-Gómez, Wenjin Yang, Christian Henkel, Sergio Martín

Abstract: Class\,I methanol masers provide sensitive information about the shocked environment around star-forming regions. Among the brightest Class~I methanol masers, we have those in the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series, currently reported for the $J=4-9$ transitions, with the only exception being the $J=7$ one at 181.295~GHz, and never expanded to higher $J$ transitions. We aim to search fo… ▽ More Class\,I methanol masers provide sensitive information about the shocked environment around star-forming regions. Among the brightest Class~I methanol masers, we have those in the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series, currently reported for the $J=4-9$ transitions, with the only exception being the $J=7$ one at 181.295~GHz, and never expanded to higher $J$ transitions. We aim to search for population inversion in the $7_{-1}\rightarrow6_{0}-E$ and $10_{-1}\rightarrow9_{0}-E$ methanol transition lines at 181.295 and 326.961\,GHz, respectively, and also extend the number of known low-mass star-forming sources harboring Class\,I methanol masers. We employed the Atacama Pathfinder Experiment (APEX) 12\,m telescope to survey low-mass Galactic sources, focusing on methanol emission lines. We conducted rotation diagrams for all sources with detected $J=7$ methanol line transitions, while employing radiative transfer modeling (both in and out of local thermodynamic equilibrium) to characterize methanol excitation conditions in detail for one specific source with detected masers. We have detected the $7_{-1}\rightarrow6_{0}-E$ and $10_{-1}\rightarrow9_{0}-E$ methanol transitions in six out of nineteen sources. Among them, we firmly determined the $10_{-1}\rightarrow9_{0}-E$ maser nature in CARMA\,7, L1641N, NGC\,2024, and Serpens FIRS, and we claim for the presence of inverted population emission in the $7_{-1}\rightarrow6_{0}-E$ line toward CARMA\,7 and L1641N. This represents the first report of methanol maser emission in these particular transitions. Our study supports previous works indicating that conditions for Class\,I methanol maser emission are satisfied in low-mass star-forming regions and expands the range of detectable frequencies toward higher values. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Report number: aa51008-24

Journal ref: A&A 688, L1 (2024)

Interstellar Nitrogen Isotope Ratios: Measurements on tracers of C$^{14}$N and C$^{15}$N

Authors: J. L. Chen, J. S. Zhang, C. Henkel, Y. T. Yan, H. Z. Yu, Y. X. Wang, Y. P. Zou, J. Y. Zhao, X. Y. Wang

Abstract: The nitrogen isotope ratio 14N/15N is a powerful tool to trace Galactic stellar nucleosynthesis and constraining Galactic chemical evolution. Previous observations have found lower 14N/15N ratios in the Galactic center and higher values in the Galactic disk. This is consistent with the inside-out formation scenario of our Milky Way. However, previous studies mostly utilized double isotope ratios a… ▽ More The nitrogen isotope ratio 14N/15N is a powerful tool to trace Galactic stellar nucleosynthesis and constraining Galactic chemical evolution. Previous observations have found lower 14N/15N ratios in the Galactic center and higher values in the Galactic disk. This is consistent with the inside-out formation scenario of our Milky Way. However, previous studies mostly utilized double isotope ratios also including 12C/13C, which introduces additional uncertainties. Here we therefore present observations of C14N and its rare isotopologue, C15N, toward a sample of star forming regions, measured by the IRAM 30 m and/or the ARO 12 m telescope at $λ$ ~3 mm wavelength. For those 35 sources detected in both isotopologues, physical parameters are determined. Furthermore we have obtained nitrogen isotope ratios using the strongest hyperfine components of CN and C15N. For those sources showing small deviations from Local Thermodynamical Equilibrium and/or self-absorption, the weakest hyperfine component, likely free of the latter effect, was used to obtain reliable 14N/15N values. Our measured 14N/15N isotope ratios from C14N and C15N measurements are compatible with those from our earlier measurements of NH3 and 15NH3 (Paper I), i.e., increasing ratios to a Galacticentric distance of ~9 kpc. The unweighted second order polynomial fit yields $\frac{{\rm C^{14}N}}{{\rm C^{15}N}} = (-4.85 \pm 1.89)\;{\rm kpc^{-2}} \times R_{\rm GC}^{2} + (82.11 \pm 31.93) \;{\rm kpc^{-1}} \times R_{\rm GC} - (28.12 \pm 126.62)$. Toward the outer galaxy, the isotope ratio tends to decrease, supporting an earlier finding by H13CN/HC15N. Galactic chemical evolution models are consistent with our measurements of the 14N/15N isotope ratio, i.e. a rising trend from the Galactic center region to approximately 9 kpc, followed by a decreasing trend with increasing $R_{\rm GC}$ toward the outer Galaxy. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 34 pages, 9 figures, 6 tables

Journal ref: The Astrophysical Journal (2004)

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

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

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

Submitted 11 June, 2024; originally announced June 2024.

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

Kinetic temperature of massive star-forming molecular clumps measured with formaldehyde V. The massive filament DR21

Authors: X. Zhao, X. D. Tang, C. Henkel, Y. Gong, Y. Lin, D. L. Li, Y. X. He, Y. P. Ao, X. Lu, T. Liu, Y. Sun, K. Wang, X. P. Chen, J. Esimbek, J. J. Zhou, J. W. Wu, J. J. Qiu, X. W. Zheng, J. S. Li, C. S. Luo, Q. Zhao

Abstract: The kinetic temperature structure of the massive filament DR21 has been mapped using the IRAM 30 m telescope. This mapping employed the para-H$_2$CO triplet ($J_{\rm K_aK_c}$ = 3$_{03}$--2$_{02}$, 3$_{22}$--2$_{21}$, and 3$_{21}$--2$_{20}$) on a scale of $\sim$0.1 pc. By modeling the averaged line ratios of para-H$_{2}$CO with RADEX under non-LTE assumptions, the kinetic temperature of the dense g… ▽ More The kinetic temperature structure of the massive filament DR21 has been mapped using the IRAM 30 m telescope. This mapping employed the para-H$_2$CO triplet ($J_{\rm K_aK_c}$ = 3$_{03}$--2$_{02}$, 3$_{22}$--2$_{21}$, and 3$_{21}$--2$_{20}$) on a scale of $\sim$0.1 pc. By modeling the averaged line ratios of para-H$_{2}$CO with RADEX under non-LTE assumptions, the kinetic temperature of the dense gas was derived at a density of $n$(H$_{2}$) = 10$^{5}$ cm$^{-3}$. The para-H$_2$CO lines reveal significantly higher temperatures than NH$_3$ (1,1)/(2,2) and FIR wavelengths. The dense clumps appear to correlate with the notable kinetic temperature. Among the four dense cores (N44, N46, N48, and N54), temperature gradients are observed on a scale of $\sim$0.1-0.3 pc. This suggests that the warm dense gas is influenced by internal star formation activity. With the exception of N54, the temperature profiles of these cores were fitted with power-law indices ranging from $-$0.3 to $-$0.5. This indicates that the warm dense gas is heated by radiation emitted from internally embedded protostar(s) and/or clusters. While there is no direct evidence supporting the idea that the dense gas is heated by shocks resulting from a past explosive event in the DR21 region, our measurements toward the DR21W1 region provide compelling evidence that the dense gas is indeed heated by shocks originating from the western DR21 flow. Higher temperatures appear to be associated with turbulence. The physical parameters of the dense gas in the DR21 filament exhibit a remarkable similarity to the results obtained in OMC-1 and N113. This may imply that the physical mechanisms governing the dynamics and thermodynamics of dense gas traced by H$_{2}$CO in diverse star formation regions may be dominated by common underlying principles despite variations in specific environmental conditions. (abbreviated) △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 16 pages, 8 figures, 3 tabels. Accepted for publication by Astronomy & Astrophysics

Journal ref: A&A 687, A207 (2024)

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

Journal ref: A&A 689, A64 (2024)

First detection of CF$^{+}$ in the Large Magellanic Cloud

Authors: Yan Gong, Karl M. Menten, Arshia M. Jacob, Christian Henkel, C. -H. Rosie Chen

Abstract: CF$^{+}$ has been established as a valuable diagnostic tool for investigating photo-dissociation regions (PDRs) and fluorine abundances in the Milky Way. However, its role in extragalactic environments remains largely uncharted. Our objective is to explore the significance of CF$^{+}$ in the Large Magellanic Cloud (LMC) and assess its utility as a valuable probe for examining C$^{+}$ and fluorine… ▽ More CF$^{+}$ has been established as a valuable diagnostic tool for investigating photo-dissociation regions (PDRs) and fluorine abundances in the Milky Way. However, its role in extragalactic environments remains largely uncharted. Our objective is to explore the significance of CF$^{+}$ in the Large Magellanic Cloud (LMC) and assess its utility as a valuable probe for examining C$^{+}$ and fluorine abundances in external galaxies. We performed pointed CF$^{+}$ observations toward an active star-forming region, N113 in the LMC, using the Atacama Pathfinder EXperiment 12~m sub-millimeter telescope. We report the first discovery of CF$^{+}$ in the LMC through the successful detection of the CF$^{+}$ (2$\to$1) and (3$\to$2) lines. The excitation models indicate that CF$^{+}$ emission originates from dense PDRs characterized by an H$_{2}$ number density of $(0.5-7.9)\times 10^{4}$~cm$^{-3}$ in N113. Our observations provide the first constraint on the fluorine abundance in molecular clouds in the LMC, disclosing a value of $\lesssim 1.7\times 10^{-9}$. This value is about an order of magnitude lower than those previously measured toward red giants in the LMC, indicative of fluorine deficiency in the molecular gas. The estimated column density ratio between C$^{+}$ and CF$^{+}$ appears to be lower than the anticipated equilibrium ratio derived from the fluorine abundance in red giants. Both phenomena can be explained by the deficiency of CF$^{+}$ caused by the freeze-out of its primary chemical precursor, HF, onto dust grains. The deficiency of CF$^{+}$ within molecular clouds suggests that the measurements presented in this work serve exclusively as conservative estimates, establishing lower bounds for both the fluorine abundance and C$^{+}$ column densities in external galaxies. △ Less

Submitted 7 May, 2024; originally announced May 2024.

Comments: 9 pages, 6 figures, 1 table, accepted for publication in A&A

Journal ref: A&A 687, A29 (2024)

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

Discovery of widespread non-metastable ammonia masers in the Milky Way

Authors: Y. T. Yan, C. Henkel, K. M. Menten, T. L. Wilson, A. Wootten, Y. Gong, F. Wyrowski, W. Yang, A. Brunthaler, A. Kraus, B. Winkel

Abstract: We present the results of a search for ammonia maser emission in 119 Galactic high-mass star-forming regions (HMSFRs) known to host 22 GHz H$_2$O maser emission. Our survey has led to the discovery of non-metastable NH$_3$ inversion line masers toward 14 of these sources. This doubles the number of known non-metastable ammonia masers in our Galaxy, including nine new very high excitation ($J,K$)~=… ▽ More We present the results of a search for ammonia maser emission in 119 Galactic high-mass star-forming regions (HMSFRs) known to host 22 GHz H$_2$O maser emission. Our survey has led to the discovery of non-metastable NH$_3$ inversion line masers toward 14 of these sources. This doubles the number of known non-metastable ammonia masers in our Galaxy, including nine new very high excitation ($J,K$)~=~(9,6) maser sources. These maser lines, including NH$_3$ (5,4), (6,4), (6,5), (7,6), (8,6), (9,6), (9,8), (10,8), and (11,9), arise from energy levels of 342 K, 513 K, 465 K, 606 K, 834 K, 1090 K, 942 K, 1226 K, and 1449 K above the ground state. Additionally, we tentatively report a new metastable NH$_3$ (3,3) maser in G048.49 and an NH$_3$ (7,7) maser in G029.95. Our observations reveal that all of the newly detected NH$_3$ maser lines exhibit either blueshifted or redshifted velocities with respect to the source systemic velocities. Among the non-metastable ammonia maser lines, larger velocity distributions, offset from the source systemic velocities, are found in the ortho-NH$_3$ ($K=3n$) than in the para-NH$_3$ ($K\neq3n$) transitions. △ Less

Submitted 12 May, 2024; v1 submitted 26 March, 2024; originally announced March 2024.

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

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

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

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

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

Submitted 25 March, 2024; originally announced March 2024.

Comments: 9 pages, 10 figures

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

Magnetic field morphology and evolution in the Central Molecular Zone and its effect on gas dynamics

Authors: R. G. Tress, M. C. Sormani, P. Girichidis, S. C. O. Glover, R. S. Klessen, R. J. Smith, E. Sobacchi, L. Armillotta, A. T. Barnes, C. Battersby, K. R. J. Bogue, N. Brucy, L. Colzi, C. Federrath, P. García, A. Ginsburg, J. Göller, H P. Hatchfield, C. Henkel, P. Hennebelle, J. D. Henshaw, M. Hirschmann, Y. Hu, J. Kauffmann, J. M. D. Kruijssen , et al. (12 additional authors not shown)

Abstract: The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We f… ▽ More The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We find that: (1) The magnetic field is conveniently decomposed into a regular time-averaged component and an irregular turbulent component. The regular component aligns well with the velocity vectors of the gas everywhere, including within the bar lanes. (2) The field geometry transitions from parallel to the Galactic plane near $z=0$ to poloidal away from the plane. (3) The magneto-rotational instability (MRI) causes an in-plane inflow of matter from the CMZ gas ring towards the central few parsecs of $0.01-0.1$ M$_\odot$ yr$^{-1}$ that is absent in the unmagnetised simulations. However, the magnetic fields have no significant effect on the larger-scale bar-driven inflow that brings the gas from the Galactic disc into the CMZ. (4) A combination of bar inflow and MRI-driven turbulence can sustain a turbulent vertical velocity dispersion of $σ_z \simeq 5$ km s$^{-1}$ on scales of $20$ pc in the CMZ ring. The MRI alone sustains a velocity dispersion of $σ_z \simeq 3$ km s$^{-1}$. Both these numbers are lower than the observed velocity dispersion of gas in the CMZ, suggesting that other processes such as stellar feedback are necessary to explain the observations. (5) Dynamo action driven by differential rotation and the MRI amplifies the magnetic fields in the CMZ ring until they saturate at a value that scales with the average local density as $B \simeq 102 (n/10^3 cm^{-3})^{0.33}$ $μ$G. Finally, we discuss the implications of our results within the observational context in the CMZ. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: submitted to A&A, comments welcome

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

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

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

Submitted 27 February, 2024; originally announced February 2024.

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

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)

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

Ammonia Observations of Planck Cold Cores

Authors: Dilda Berdikhan., Jarken Esimbek., Christian Henkel., Jianjun Zhou., Xindi Tang., Tie Liu., Gang Wu., Dalei Li., Yuxin He., Toktarkhan Komesh., Kadirya Tursun., Dongdong Zhou., Ernar Imanaly., Qaynar Jandaolet

Abstract: Single-pointing observations of NH$_3$ (1,1) and (2,2) were conducted towards 672 Planck Early Release Cold Cores (ECCs) using the Nanshan 26-m radio telescope. Out of these sources, a detection rate of 37% (249 cores) was achieved, with NH$_3$(1,1) hyperfine structure detected in 187 and NH$_3$(2,2) emission lines detected in 76 cores. The detection rate of NH3 is positively correlated with the c… ▽ More Single-pointing observations of NH$_3$ (1,1) and (2,2) were conducted towards 672 Planck Early Release Cold Cores (ECCs) using the Nanshan 26-m radio telescope. Out of these sources, a detection rate of 37% (249 cores) was achieved, with NH$_3$(1,1) hyperfine structure detected in 187 and NH$_3$(2,2) emission lines detected in 76 cores. The detection rate of NH3 is positively correlated with the continuum emission fluxes at a frequency of 857 GHz. Among the observed 672 cores, ~22% have associated stellar and IR objects within the beam size (~2$\arcmin$). This suggests that most of the cores in our sample may be starless. The kinetic temperatures of the cores range from 8.9 to 20.7 K, with an average of 12.3 K, indicating a coupling between gas and dust temperatures. The ammonia column densities range from 0.36 to 6.07$\times10^{15}$ cm$^{-2}$, with a median value of 2.04$\times10^{15}$ cm$^{-2}$. The fractional abundances of ammonia range from 0.3 to 9.7$\times10^{-7}$, with an average of 2.7 $\times10^{-7}$, which is one order of magnitude larger than that of Massive Star-Forming (MSF) regions and Infrared Dark Clouds (IRDCs). The correlation between thermal and non-thermal velocity dispersion of the NH$_3$(1,1) inversion transition indicates the dominance of supersonic non-thermal motions in the dense gas traced by NH$_3$, and the relationship between these two parameters in Planck cold cores is weaker, with lower values observed for both parameters relative to other samples under our examination. The cumulative distribution shapes of line widths in the Planck cold cores closely resemble those of the dense cores found in regions of Cepheus, and Orion L1630 and L1641, with higher values compared to Ophiuchus. A comparison of NH3 line-center velocities with those of $^{13}$CO and C$^{18}$O shows small differences (0.13 and 0.12 km s$^{-1}$ ), suggesting quiescence on small scales. △ Less

Submitted 4 January, 2024; originally announced January 2024.

Comments: 29 pages, 13 figures

MSC Class: 85 ACM Class: A.0

What Determines the Boundaries of H2O Maser Emission in an X-ray Illuminated Gas Disk ?

Authors: C. Y. Kuo, F. Gao, J. A. Braatz, D. W. Pesce, E. M. L. Humphreys, M. J. Reid, C. M. V. Impellizzeri, C. Henkel, J. Wagner, C. E. Wu

Abstract: High precision mapping of H2O megamaser emission from active galaxies has revealed more than a dozen Keplerian H2O maser disks, which enable a ~4% uncertainty estimate of the Hubble constant as well as providing accurate masses for the central black holes. These disks often have well-defined inner and outer boundaries of maser emission on sub-parsec scales. In order to better understand the physic… ▽ More High precision mapping of H2O megamaser emission from active galaxies has revealed more than a dozen Keplerian H2O maser disks, which enable a ~4% uncertainty estimate of the Hubble constant as well as providing accurate masses for the central black holes. These disks often have well-defined inner and outer boundaries of maser emission on sub-parsec scales. In order to better understand the physical conditions that determine the inner and outer radii of a maser disk, we examine the distributions of gas density and X-ray heating rate in a warped molecular disk described by a power-law surface density profile. For a suitable choice of the disk mass, we find that the outer radius R_out of the maser disk predicted from our model can match the observed value, with R_out mainly determined by the maximum heating rate or the minimum density for efficient maser action, depending on the combination of the Eddington ratio, black hole mass, and disk mass. Our analysis also indicates that the inner radius for maser action is comparable to the dust sublimation radius, suggesting that dust may play a role in determining the inner radius of a maser disk. Finally, our model predicts that H2O gigamaser disks could exist at the centers of high-z quasars, with disk sizes of >~ 10-30 pc. △ Less

Submitted 10 July, 2024; v1 submitted 26 December, 2023; originally announced December 2023.

Comments: Accepted by MNRAS, 17 pages, 8 figures, 2 tables

The first Ka-band (26.1-35 GHz) blind line survey towards Orion KL

Authors: Xunchuan Liu, Tie Liu, Zhiqiang Shen, Sheng-Li Qin, Qiuyi Luo, Yan Gong, Yu Cheng, Christian Henkel, Qilao Gu, Fengyao Zhu, Tianwei Zhang, Rongbing Zhao, Yajun Wu, Bin Li, Juan Li, Zhang Zhao, Jinqing Wang, Weiye Zhong, Qinghui Liu, Bo Xia, Li Fu, Zhen Yan, Chao Zhang, Lingling Wang, Qian Ye , et al. (9 additional authors not shown)

Abstract: We conducted a Ka-band (26.1--35 GHz) line survey towards Orion KL using the TianMa 65-m Radio Telescope (TMRT). It is the first blind line survey in the Ka band, and achieves a sensitivity of mK level (1--3 mK at a spectral resolution of $\sim$1 km s$^{-1}$). In total, 592 Gaussian features are extracted. Among them, 257 radio recombination lines (RRLs) are identified. The maximum $Δn$ of RRLs of… ▽ More We conducted a Ka-band (26.1--35 GHz) line survey towards Orion KL using the TianMa 65-m Radio Telescope (TMRT). It is the first blind line survey in the Ka band, and achieves a sensitivity of mK level (1--3 mK at a spectral resolution of $\sim$1 km s$^{-1}$). In total, 592 Gaussian features are extracted. Among them, 257 radio recombination lines (RRLs) are identified. The maximum $Δn$ of RRLs of H, He and C are 20, 15, and 5, respectively. Through stacking, we have detected the $β$ lines of ion RRLs (RRLs of C$^+$ with possible contribution of other ions like O$^+$) for the first time, and tentative signal of the $γ$ lines of ion RRLs can also be seen on the stacked spectrum. Besides, 318 other line features were assigned to 37 molecular species, and ten of these species were not detected in the Q-band survey of TMRT. The vibrationally excited states of nine species were also detected. Emission of most species can be modeled under LTE. A number of transitions of E-CH3OH ($J_2-J_1$) display maser effects, which are confirmed by our modeling, and besides the bumping peak at $J\sim 6$ there is another peak at $J\sim 13$. Methylcyanoacetylene (CH$_3$C$_3$N) is detected in Orion KL for the first time. This work emphasizes that the Ka band, which was long-ignored for spectral line surveys, is very useful for surveying RRLs and molecular lines simultaneously. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Comments: accepted by ApJS

A Systematic Observational Study on Galactic Interstellar Ratio 18O/17O. II. C18O and C17O J=2-1 Data Analysis

Authors: Y. P. Zou, J. S. Zhang, C. Henkel, D. Romano, W. Liu, Y. H. Zheng, Y. T. Yan, J. L. Chen, Y. X. Wang, J. Y. Zhao

Abstract: To investigate the relative amount of ejecta from high-mass versus intermediate-mass stars and to trace the chemical evolution of the Galaxy, we have performed with the IRAM 30m and the SMT 10m telescopes a systematic study of Galactic interstellar 18O/17O ratios toward a sample of 421 molecular clouds, covering a galactocentric distance range of 1-22 kpc. The results presented in this paper are b… ▽ More To investigate the relative amount of ejecta from high-mass versus intermediate-mass stars and to trace the chemical evolution of the Galaxy, we have performed with the IRAM 30m and the SMT 10m telescopes a systematic study of Galactic interstellar 18O/17O ratios toward a sample of 421 molecular clouds, covering a galactocentric distance range of 1-22 kpc. The results presented in this paper are based on the J=2-1 transition and encompass 364 sources showing both C18O and C17O detections. The previously suggested 18O/17O gradient is confirmed. For the 41 sources detected with both facilities, good agreement is obtained. A correlation of 18O/17O ratios with heliocentric distance is not found, indicating that beam dilution and linear beam sizes are not relevant. For the subsample of IRAM 30 m high-mass star-forming regions with accurate parallax distances, an unweighted fit gives 18O/17O = (0.12+-0.02)R_GC+(2.38+-0.13) with a correlation coefficient of R = 0.67. While the slope is consistent with our J=1-0 measurement, ratios are systematically lower. This should be caused by larger optical depths of C18O 2-1 lines, w.r.t the corresponding 1-0 transitions, which is supported by RADEX calculations and the fact that C18O/C17O is positively correlated with 13CO/C18O. After considering optical depth effects with C18O J=2-1 reaching typically an optical depth of 0.5, corrected 18O/17O ratios from the J=1-0 and J=2-1 lines become consistent. A good numerical fit to the data is provided by the MWG-12 model, including both rotating stars and novae. △ Less

Submitted 1 November, 2023; originally announced November 2023.

Comments: 17 pages, 11 figures, published in ApJS

Variations of the HCO$^{+}$, HCN, HNC, N$_2$H$^+$ and NH$_{3}$ deuterium fractionation in high-mass star-forming regions

Authors: A. G. Pazukhin, I. I. Zinchenko, E. A. Trofimova, C. Henkel, D. A. Semenov

Abstract: We use spectra and maps of the $J=1-0$ and $J=2-1$ DCO$^{+}$, DCN, DNC, $\rm N_2D^+$ lines and $1_{11}-1_{01}$ ortho- and para-NH$_{2}$D lines, obtained with the IRAM-30m telescope, as well as observations of their hydrogenated isotopologues to study deuteration processes in five high-mass star-forming regions. The temperature was estimated from CH$_3$CCH lines, also observed with the IRAM-30m tel… ▽ More We use spectra and maps of the $J=1-0$ and $J=2-1$ DCO$^{+}$, DCN, DNC, $\rm N_2D^+$ lines and $1_{11}-1_{01}$ ortho- and para-NH$_{2}$D lines, obtained with the IRAM-30m telescope, as well as observations of their hydrogenated isotopologues to study deuteration processes in five high-mass star-forming regions. The temperature was estimated from CH$_3$CCH lines, also observed with the IRAM-30m telescope, and from NH$_3$ lines, observed with the 100-m radio telescope in Effelsberg, as well as using the integrated intensity ratios of the $J=1-0$ H$^{13}$CN and HN$^{13}$C lines and their main isotopologues. Applying a non-LTE radiative transfer model with RADEX, the gas density and the molecular column densities were estimated. D/H ratios are $0.001-0.05$ for DCO$^{+}$, $0.001-0.02$ for DCN, $0.001-0.05$ for DNC and $0.02-0.4$ for NH$_{2}$D. The D/H ratios decrease with increasing temperature in the range of $\rm 20-40 \,K$ and slightly vary at densities $n(\rm H_2) \sim 10^4-10^6\, cm^{-3}$. The deuterium fraction of $\rm N_2H^{+}$ is $0.008-0.1$ at temperatures in the range of $\rm 20-25\, K$ and at a density of $\sim 10^5\, \rm cm^{-3}$. We also estimate relative abundances and find $ \sim 10^{-11}-10^{-9}$ for DCO$^{+}$ and DNC, $ \sim 10^{-11}-10^{-10}$ for $\rm N_2D^+$ and $ \sim 10^{-10}-10^{-8}$ for NH$_{2}$D. The relative abundances of these species decrease with increasing temperature. However, the DCN/H$_2$ ratio is almost constant ($\sim 10^{-10}$). The observational results agree with the predictions of chemical models (although in some cases there are significant differences). △ Less

Submitted 28 September, 2023; originally announced September 2023.

Comments: 13 pages + appendix (25 pages total), 11 figures. Accepted for publication in MNRAS

Sulfur isotope ratios in the Large Magellanic Cloud

Authors: Y. Gong, C. Henkel, K. M. Menten, C. -H. R. Chen, Z. Y. Zhang, Y. T. Yan, A. Weiss, N. Langer, J. Z. Wang, R. Q. Mao, X. D. Tang, W. Yang, Y. P. Ao, M. Wang

Abstract: Sulfur isotope ratios have emerged as a promising tool for tracing stellar nucleosynthesis, quantifying stellar populations, and investigating the chemical evolution of galaxies. While extensively studied in the Milky Way, in extragalactic environments they remain largely unexplored. We focus on investigating the sulfur isotope ratios in the Large Magellanic Cloud (LMC) to gain insights into sulfu… ▽ More Sulfur isotope ratios have emerged as a promising tool for tracing stellar nucleosynthesis, quantifying stellar populations, and investigating the chemical evolution of galaxies. While extensively studied in the Milky Way, in extragalactic environments they remain largely unexplored. We focus on investigating the sulfur isotope ratios in the Large Magellanic Cloud (LMC) to gain insights into sulfur enrichment in this nearby system and to establish benchmarks for such ratios in metal-poor galaxies. We conducted pointed observations of CS and its isotopologues toward N113, one of the most prominent star-formation regions in the LMC, utilizing the Atacama Pathfinder EXperiment 12~m telescope. We present the first robust detection of C$^{33}$S in the LMC by successfully identifying two C$^{33}$S transitions on a large scale of $\sim$5 pc. Our measurements result in an accurate determination of the $^{34}$S/$^{33}$S isotope ratio, which is 2.0$\pm$0.2. Our comparative analysis indicates that the $^{32}$S/$^{33}$S and $^{34}$S/$^{33}$S isotope ratios are about a factor of 2 lower in the LMC than in the Milky Way. Our findings suggest that the low $^{34}$S/$^{33}$S isotope ratio in the LMC can be attributed to a combination of the age effect, low metallicity, and star formation history. △ Less

Submitted 18 October, 2023; v1 submitted 26 September, 2023; originally announced September 2023.

Comments: 6 pages, 1 figures, 2 tables, accepted for publication in A&A, adjusted to the final version

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

Protonated hydrogen cyanide as a tracer of pristine molecular gas

Authors: Y. Gong, F. J. Du, C. Henkel, A. M. Jacob, A. Belloche, J. Z. Wang, K. M. Menten, W. Yang, D. H. Quan, C. T. Bop, G. N. Ortiz-León, X. D. Tang, M. R. Rugel, S. Liu

Abstract: Protonated hydrogen cyanide, HCNH$^{+}$, plays a fundamental role in astrochemistry because it is an intermediary in gas-phase ion-neutral reactions within cold molecular clouds. However, the impact of the environment on the chemistry of HCNH$^{+}$ remains poorly understood. With the IRAM-30 m and APEX-12 m observations, we report the first robust distribution of HCNH$^{+}$ in the Serpens filament… ▽ More Protonated hydrogen cyanide, HCNH$^{+}$, plays a fundamental role in astrochemistry because it is an intermediary in gas-phase ion-neutral reactions within cold molecular clouds. However, the impact of the environment on the chemistry of HCNH$^{+}$ remains poorly understood. With the IRAM-30 m and APEX-12 m observations, we report the first robust distribution of HCNH$^{+}$ in the Serpens filament and in Serpens South. Our data suggest that HCNH$^{+}$ is abundant in cold and quiescent regions, but is deficit in active star-forming regions. The observed HCNH$^{+}$ fractional abundances relative to H$_{2}$ range from $3.1\times 10^{-11}$ in protostellar cores to $5.9\times 10^{-10}$ in prestellar cores, and the HCNH$^{+}$ abundance generally decreases with increasing H$_{2}$ column density, which suggests that HCNH$^{+}$ coevolves with cloud cores. Our observations and modeling results suggest that the abundance of HCNH$^{+}$ in cold molecular clouds is strongly dependent on the H$_{2}$ number density. The decrease in the abundance of HCNH$^{+}$ is caused by the fact that its main precursors (e.g., HCN and HNC) undergo freeze-out as the number density of H$_{2}$ increases. However, current chemical models cannot explain other observed trends, such as the fact that the abundance of HCNH$^{+}$ shows an anti-correlation with that of HCN and HNC, but a positive correlation with that of N$_{2}$H$^{+}$ in the southern part of the Serpens South northern clump. This indicates that additional chemical pathways have to be invoked for the formation of HCNH$^{+}$ via molecules like N$_{2}$ in regions in which HCN and HNC freeze out. Both the fact that HCNH$^{+}$ is most abundant in molecular cores prior to gravitational collapse and the fact that low-$J$ HCNH$^{+}$ transitions have very low H$_{2}$ critical densities make this molecular ion an excellent probe of pristine molecular gas. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 25 pages, 26 figures, accepted for publication in A&A

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

The 107 GHz methanol transition is a dasar in G0.253+0.016

Authors: Alyssa Bulatek, Adam Ginsburg, Jeremy Darling, Christian Henkel, Karl M. Menten

Abstract: We present observations of population anti-inversion in the $3_1 - 4_0\ A^+$ transition of CH$_3$OH (methanol) at 107.013831 GHz toward the Galactic Center cloud G0.253+0.016 ("The Brick"). Anti-inversion of molecular level populations can result in absorption lines against the cosmic microwave background (CMB) in a phenomenon known as a "dasar." We model the physical conditions under which the 10… ▽ More We present observations of population anti-inversion in the $3_1 - 4_0\ A^+$ transition of CH$_3$OH (methanol) at 107.013831 GHz toward the Galactic Center cloud G0.253+0.016 ("The Brick"). Anti-inversion of molecular level populations can result in absorption lines against the cosmic microwave background (CMB) in a phenomenon known as a "dasar." We model the physical conditions under which the 107 GHz methanol transition dases and determine that dasing occurs at densities below $10^6$ cm$^{-3}$ and column densities between $10^{13}$ and $10^{16}$ cm$^{-2}$. We also find that for this transition, dasing does not strongly depend on the gas kinetic temperature. We evaluate the potential of this tool for future deep galaxy surveys. We note that other works have already reported absorption in this transition (e.g., in NGC 253), but we provide the first definitive evidence that it is absorption against the CMB rather than against undetected continuum sources. △ Less

Submitted 16 August, 2023; originally announced August 2023.

Comments: 19 pages, 9 figures, 1 table. Accepted for publication in ApJ (August 14, 2023)

Star formation in the centre of NGC 1808 as observed by ALMA

Authors: Guangwen Chen, George J. Bendo, Gary A. Fuller, Christian Henkel, Xu Kong

Abstract: We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 85.69 and 99.02 GHz continuum emission and H42$α$ and H40$α$ lines emission from the central 1~kpc of NGC 1808. These forms of emission are tracers of photoionizing stars but unaffected by dust obscuration that we use to test the applicability of other commonly star formation metrics. An analysis of the spectral energy… ▽ More We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 85.69 and 99.02 GHz continuum emission and H42$α$ and H40$α$ lines emission from the central 1~kpc of NGC 1808. These forms of emission are tracers of photoionizing stars but unaffected by dust obscuration that we use to test the applicability of other commonly star formation metrics. An analysis of the spectral energy distributions shows that free-free emission contributes about 60 to 90 per cent of the continuum emission in the 85-100 GHz frequency range, dependent on the region. The star formation rate (SFR) derived from the ALMA free-free emission is $3.1\pm0.3$~M$_\odot$~yr$^{-1}$. This is comparable to the SFRs measured from the infrared emission, mainly because most of the bolometric energy from the heavily obscured region is emitted as infrared emission. The radio 1.5~GHz emission yields a SFR 25 per cent lower than the ALMA value, probably because of the diffusion of the electrons producing the synchrotron emission beyond the star-forming regions. The SFRs measured from the extinction-corrected H$α$ line emission are about 40 to 65 per cent of the SFR derived from the ALMA data, likely because this metric was not calibrated for high extinction regions. Some SFRs based on extinction-corrected ultraviolet emission are similar to those from ALMA and infrared data, but given that the ultraviolet terms in the extinction correction equations are very small, these metrics seem inappropriate to apply to this dusty starburst. △ Less

Submitted 16 August, 2023; originally announced August 2023.

Comments: 18 pages, 7 figures, 10 tables, accepted for publication in MNRAS

A global view on star formation: The GLOSTAR Galactic plane survey VIII. Formaldehyde absorption in Cygnus~X

Authors: Y. Gong, G. N. Ortiz-León, M. R. Rugel, K. M. Menten, A. Brunthaler, F. Wyrowski, C. Henkel, H. Beuther, S. A. Dzib, J. S. Urquhart, A. Y. Yang, J. D. Pandian, R. Dokara, V. S. Veena, H. Nguyen, S. -N. X. Medina, W. D. Cotton, W. Reich, B. Winkel, P. Müller, I. Skretas, T. Csengeri, S. Khan, A. Cheema

Abstract: Cygnus X is one of the closest and most active high-mass star-forming regions in our Galaxy, making it one of the best laboratories for studying massive star formation. As part of the GLOSTAR Galactic plane survey, we performed large scale simultaneous H$_{2}$CO (1$_{1,0}$-1$_{1,1}$) spectral line and radio continuum imaging observations toward Cygnus X at $λ\sim$6 cm with the Karl G. Jansky Very… ▽ More Cygnus X is one of the closest and most active high-mass star-forming regions in our Galaxy, making it one of the best laboratories for studying massive star formation. As part of the GLOSTAR Galactic plane survey, we performed large scale simultaneous H$_{2}$CO (1$_{1,0}$-1$_{1,1}$) spectral line and radio continuum imaging observations toward Cygnus X at $λ\sim$6 cm with the Karl G. Jansky Very Large Array and the Effelsberg-100 m radio telescope. Our Effelsberg observations reveal widespread H$_{2}$CO (1$_{1,0}$-1$_{1,1}$) absorption with a spatial extent of $\gtrsim$50 pc in Cygnus~X for the first time. On large scales of 4.4 pc, the relative orientation between local velocity gradient and magnetic field tends to be more parallel at H$_{2}$ column densities of $\gtrsim$1.8$\times 10^{22}$~cm$^{-2}$. On the smaller scale of 0.17 pc, our VLA+Effelsberg combined data reveal H$_{2}$CO absorption only toward three bright H{\scriptsize II} regions. Our observations demonstrate that H$_{2}$CO (1$_{1,0}$-1$_{1,1}$) is commonly optically thin. Kinematic analysis supports the assertion that molecular clouds generally exhibit supersonic motions on scales of 0.17-4.4 pc. We show a non-negligible contribution of the cosmic microwave background radiation in producing extended absorption features in Cygnus X. Our observations suggest that H$_{2}$CO ($1_{1,0}-1_{1,1}$) can trace molecular gas with H$_{2}$ column densities of $\gtrsim 5 \times 10^{21}$ cm$^{-2}$. The ortho-H$_{2}$CO fractional abundance with respect to H$_{2}$ has a mean value of 7.0$\times 10^{-10}$. A comparison of velocity dispersions on different linear scales suggests that the dominant $-3$ km s$^{-1}$ velocity component in the prominent DR21 region has nearly identical velocity dispersions on scales of 0.17-4.4 pc, which deviates from the expected behavior of classic turbulence. △ Less

Submitted 2 August, 2023; originally announced August 2023.

Comments: 27 pages, 23 figures, accepted for publication in A&A

Journal ref: A&A 678, A130 (2023)

ALMA and VLBA views on the outflow associated with an O-type protostar in G26.50+0.28

Authors: Gang Wu, Christian Henkel, Ye Xu, Andreas Brunthaler, Karl M. Menten, Keping Qiu, Jingjing Li, Bo Zhang, Jarken Esimbek

Abstract: Protostellar jets and outflows are essential ingredients of the star formation process. A better understanding of this phenomenon is important in its own right as well as for many fundamental aspects of star formation. Jets and outflows associated with O-type protostars are rarely studied with observations reaching the close vicinity of the protostars. In this work, we report high-resolution ALMA… ▽ More Protostellar jets and outflows are essential ingredients of the star formation process. A better understanding of this phenomenon is important in its own right as well as for many fundamental aspects of star formation. Jets and outflows associated with O-type protostars are rarely studied with observations reaching the close vicinity of the protostars. In this work, we report high-resolution ALMA and VLBA observations to reveal a clear and consistent picture of an outflow associated with an O-type protostar candidate in the G26.50+0.28 region. These observations reveal, for the first time, a collimated jet located in the middle of the outflow cavity. The jet is found to be perpendicular to an elongated disk/toroid and its velocity gradient. The collimated jet appears to show a small amplitude ($α$$\approx$0$\,.\!\!^{\circ}$06) counterclockwise precession, when looking along the blueshifted jet axis from the strongest continuum source MM1, with a precession length of 0.22 pc. The inclination of the jet is likely to be very low ($\approx$8$^{\circ}$), which makes it a promising target to study its transverse morphologies and kinematics. However, no clear evidence of jet rotation is found in the ALMA and VLBA observations. The three-dimensional velocities of the water maser spots appear to show the same absolute speed with respect to different opening angles, suggesting the jet winds may be launched in a relatively small region. This favors the X-wind model, that is, jets are launched in a small area near the inner disk edge. △ Less

Submitted 18 July, 2023; originally announced July 2023.

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

Journal ref: A&A 677, A80 (2023)

Infall Motions in the Hot Core Associated with Hypercompact HII Region G345.0061+01.794 B

Authors: Toktarkhan Komesh, Guido Garay, Christian Henkel, Aruzhan Omar, Robert Estalella, Zhandos Assembay, Dalei Li, Andrés Guzmán, Jarken Esimbek, Jiasheng Huang, Yuxin He, Nazgul Alimgazinova, Meiramgul Kyzgarina, Shukirgaliyev Bekdaulet, Nurman Zhumabay, Arailym Manapbayeva

Abstract: We report high angular resolution observations, made with the Atacama Large Millimeter Array in band 6, of high excitation molecular lines of $\rm CH_3CN$ and $\rm SO_2$ and of the H29$α$ radio recombination line towards the G345.0061+01.794 B HC H II region, in order to investigate the physical and kinematical characteristics of its surroundings. Emission was detected in all observed components o… ▽ More We report high angular resolution observations, made with the Atacama Large Millimeter Array in band 6, of high excitation molecular lines of $\rm CH_3CN$ and $\rm SO_2$ and of the H29$α$ radio recombination line towards the G345.0061+01.794 B HC H II region, in order to investigate the physical and kinematical characteristics of its surroundings. Emission was detected in all observed components of the J=14$\rightarrow$13 rotational ladder of $\rm CH_3CN$ and in the $30_{4,26}-30_{3,27}$ and $32_{4,28}-32_{3,29}$ lines of $\rm SO_2$. The peak of the velocity integrated molecular emission is located $\sim$0$\,.\!\!^{\prime\prime}$4 northwest of the peak of the continuum emission. The first-order moment images and channel maps show a velocity gradient, of 1.1 km s$^{-1}$ arcsec$^{-1}$, across the source, and a distinctive spot of blueshifted emission towards the peak of the zero-order moment. The rotational temperature is found to decrease from 252$\pm24$ Kelvin at the peak position to 166$\pm16$ Kelvin at its edge, indicating that our molecular observations are probing a hot molecular core that is internally excited. The emission in the H29$α$ line arises from a region of 0$\,.\!\!^{\prime\prime}$65 in size, where its peak coincides with that of the dust continuum. We model the kinematical characteristics of the "central blue spot" feature as due to infalling motions, suggesting a central mass of 172.8$\pm8.8 M_{\odot}$. Our observations indicate that this HC H II region is surrounded by a compact structure of hot molecular gas, which is rotating and infalling toward a central mass, that is most likely confining the ionized region. The observed scenario is reminiscent of a "butterfly pattern" with an approximately edge-on torus and ionized gas roughly parallel to its rotation axis. △ Less

Submitted 4 May, 2024; v1 submitted 14 July, 2023; originally announced July 2023.

Comments: Accepted by ApJ

ATLASGAL: 3-mm class I methanol masers in high-mass star formation regions

Authors: W. Yang, Y. Gong, K. M. Menten, J. S. Urquhart, C. Henkel, F. Wyrowski, T. Csengeri, S. P. Ellingsen, A. R. Bemis, J. Jang

Abstract: We analyzed the 3-mm wavelength spectral line survey of 408 ATLASGAL clumps observed with the IRAM 30m-telescope, focusing on the class I methanol masers with frequencies near 84, 95 and 104.3 GHz. We detect narrow, maser-like features towards 54, 100 and 4 sources in the maser lines near 84, 95 and 104.3 GHz, respectively. Among them, fifty 84 GHz masers, twenty nine 95 GHz masers and four rare 1… ▽ More We analyzed the 3-mm wavelength spectral line survey of 408 ATLASGAL clumps observed with the IRAM 30m-telescope, focusing on the class I methanol masers with frequencies near 84, 95 and 104.3 GHz. We detect narrow, maser-like features towards 54, 100 and 4 sources in the maser lines near 84, 95 and 104.3 GHz, respectively. Among them, fifty 84 GHz masers, twenty nine 95 GHz masers and four rare 104.3 GHz masers are new discoveries. The new detections increase the number of known 104.3 GHz masers from 5 to 9. The 95 GHz class I methanol maser is generally stronger than the 84 GHz maser counterpart. We find 9 sources showing class I methanol masers but no SiO emission, indicating that class I methanol masers might be the only signpost of protostellar outflow activity in extremely embedded objects at the earliest evolutionary stage. Class I methanol masers that are associated with sources that show SiO line wings are more numerous and stronger than those without such wings. The total integrated intensity of class I methanol masers is well correlated with the integrated intensity and velocity coverage of the SiO (2--1) emission. The properties of class I methanol masers are positively correlated with the bolometric luminosity, clump mass, peak H$_2$ column density of their associated clumps but uncorrelated with the luminosity-to-mass ratio, dust temperature, and mean H$_2$ volume density. We suggest that the properties of class I masers are related to shocks traced by SiO. Based on our observations, we conclude that class I methanol masers at 84 and 95 GHz can trace a similar evolutionary stage as H$_2$O maser, and appear prior to 6.7 and 12.2 GHz methanol and OH masers. Despite their small number, the 104.3 GHz class I masers appear to trace a short and more evolved stage compared to the other class I masers. [abridged] △ Less

Submitted 7 May, 2023; originally announced May 2023.

Comments: 23 pages, 27 figures, 8 tables, accepted for publication in A&A

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

Gas, dust, and the CO-to-molecular gas conversion factor in low-metallicity starbursts

Authors: L. K. Hunt, F. Belfiore, F. Lelli, B. T. Draine, A. Marasco, S. Garcia-Burillo, G. Venturi, F. Combes, A. Weiß, C. Henkel, K. M. Menten, F. Annibali, V. Casasola, M. Cignoni, A. McLeod, M. Tosi, M. Beltran, A. Concas, G. Cresci, M. Ginolfi, N. Kumari, F. Mannucci

Abstract: The factor relating CO emission to molecular hydrogen column density, XCO, is still subject to uncertainty, in particular at low metallicity. Here, to quantify XCO at two different spatial resolutions, we exploit a dust-based method together with ALMA 12-m and ACA data and HI maps of three nearby metal-poor starbursts, NGC625, NGC1705, and NGC5253. Dust opacity at 250pc resolution is derived based… ▽ More The factor relating CO emission to molecular hydrogen column density, XCO, is still subject to uncertainty, in particular at low metallicity. Here, to quantify XCO at two different spatial resolutions, we exploit a dust-based method together with ALMA 12-m and ACA data and HI maps of three nearby metal-poor starbursts, NGC625, NGC1705, and NGC5253. Dust opacity at 250pc resolution is derived based on dust temperatures estimated by fitting two-temperature modified blackbodies to Herschel PACS data. By using the HI maps, we are then able to estimate dust-to-gas ratios in the atomic-gas dominated regions, and infer total gas column densities and H2 column densities as the difference with HI. Finally, from the ACA CO(1-0) maps, we derive XCO. We use a similar technique with 40 pc ALMA 12-m data for the three galaxies, but instead derive dust attenuation at 40 pc resolution from reddening maps based on VLT/MUSE data. At 250 pc resolution, XCO $\sim$ 10^22 - 10^23 cm^-2 / K.km/s, 5-1000 times the Milky Way value, with much larger values than would be expected from a simple metallicity dependence. Instead at 40 pc resolution, XCO again shows large variation, but is roughly consistent with a power-law metallicity dependence, given the Z $\sim$ 1/3 Zsun metal abundances of our targets. The large scatter in both estimations could imply additional parameter dependence, that we have investigated by comparing XCO with the observed velocity-integrated brightness temperatures, ICO, as predicted by recent simulations. Indeed, larger XCO is significantly correlated with smaller ICO, but with slightly different slopes and normalizations than predicted by theory. Such behavior can be attributed to the increasing fraction of CO-faint H2 gas with lower spatial resolution. This confirms the idea the XCO is multi-variate, depending not only on metallicity but also on CO brightness temperature and beam size. △ Less

Submitted 11 May, 2023; v1 submitted 5 May, 2023; originally announced May 2023.

Comments: A&A, in press (16 figures, 1 appendix)

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

Kinematics of Galactic Centre clouds shaped by shear-seeded solenoidal turbulence

Authors: Maya A. Petkova, J. M. Diederik Kruijssen, Jonathan D. Henshaw, Steven N. Longmore, Simon C. O. Glover, Mattia C. Sormani, Lucia Armillotta, Ashley T. Barnes, Ralf S. Klessen, Francisco Nogueras-Lara, Robin G. Tress, Jairo Armijos-Abendaño, Laura Colzi, Christoph Federrath, Pablo García, Adam Ginsburg, Christian Henkel, Sergio Martín, Denise Riquelme, Víctor M. Rivilla

Abstract: The Central Molecular Zone (CMZ; the central ~ 500 pc of the Galaxy) is a kinematically unusual environment relative to the Galactic disc, with high velocity dispersions and a steep size-linewidth relation of the molecular clouds. In addition, the CMZ region has a significantly lower star formation rate (SFR) than expected by its large amount of dense gas. An important factor in explaining the low… ▽ More The Central Molecular Zone (CMZ; the central ~ 500 pc of the Galaxy) is a kinematically unusual environment relative to the Galactic disc, with high velocity dispersions and a steep size-linewidth relation of the molecular clouds. In addition, the CMZ region has a significantly lower star formation rate (SFR) than expected by its large amount of dense gas. An important factor in explaining the low SFR is the turbulent state of the star-forming gas, which seems to be dominated by rotational modes. However, the turbulence driving mechanism remains unclear. In this work, we investigate how the Galactic gravitational potential affects the turbulence in CMZ clouds. We focus on the CMZ cloud G0.253+0.016 (`the Brick'), which is very quiescent and unlikely to be kinematically dominated by stellar feedback. We demonstrate that several kinematic properties of the Brick arise naturally in a cloud-scale hydrodynamics simulation that takes into account the Galactic gravitational potential. These properties include the line-of-sight velocity distribution, the steepened size-linewidth relation, and the predominantly solenoidal nature of the turbulence. Within the simulation, these properties result from the Galactic shear in combination with the cloud's gravitational collapse. This is a strong indication that the Galactic gravitational potential plays a crucial role in shaping the CMZ gas kinematics, and is a major contributor to suppressing the SFR by inducing predominantly solenoidal turbulent modes. △ Less

Submitted 1 August, 2023; v1 submitted 21 April, 2023; originally announced April 2023.

Comments: 7 pages, 8 figures; accepted to MNRAS (July 24th 2023)

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)

H$_2$O MegaMaser emission in NGC 4258 indicative of a periodic disc instability

Authors: Wiillem A. Baan, Tao AN, Christian Henkel, Hiroshi Imai, Vladimir Kostenko, Andrej Sobolev

Abstract: H$_2$O MegaMaser emission may arise from thin gas discs surrounding the massive nuclei of galaxies such as NGC\,4258, but the physical conditions responsible for the amplified emission are unclear. A detailed view of these regions is possible using the very high angular resolution afforded by space very long baseline interferometry (SVLBI). Here we report SVLBI experiments conducted using the orbi… ▽ More H$_2$O MegaMaser emission may arise from thin gas discs surrounding the massive nuclei of galaxies such as NGC\,4258, but the physical conditions responsible for the amplified emission are unclear. A detailed view of these regions is possible using the very high angular resolution afforded by space very long baseline interferometry (SVLBI). Here we report SVLBI experiments conducted using the orbiting RadioAstron Observatory that have resulted in detections of the H$_2$O 22 GHz emission in NGC\,4258, with Earth-space baselines of 1.3, 9.5 and 19.5 Earth diameters. Observations at the highest angular resolution of 11 and 23 $μ$as show distinct and regularly spaced regions within the rotating disc, at an orbital radius of about 0.126 pc. These observations at three subsequent epochs also indicate a time evolution of the emission features, with a sudden rise in amplitude followed by a slow decay. The formation of the emission regions, their regular spacing and their time-dependent behaviour appear consistent with the occurrence of a periodic magneto-rotational instability in the disc. This type of shear-driven instability within the differentially rotating disc has been suggested to be the mechanism governing the radial momentum transfer and viscosity within a mass-accreting disc. The connection of the H$_2$O MegaMaser activity with the magneto-rotational instability activity would make it an indicator of the mass-accretion rate in the nuclear disc of the host galaxy. △ Less

Submitted 28 February, 2023; originally announced February 2023.

Comments: 23 pages, 5 figures

183 GHz water megamasers in active galactic nuclei: a new accretion disk tracer

Authors: Dominic W. Pesce, James A. Braatz, Christian Henkel, Elizabeth M. L. Humphreys, C. M. Violette Impellizzeri, Cheng-Yu Kuo

Abstract: We present the results of an ALMA survey to identify 183 GHz H$_2$O maser emission from AGN already known to host 22 GHz megamaser systems. Out of 20 sources observed, we detect significant 183 GHz maser emission from 13; this survey thus increases the number of AGN known to host (sub)millimeter megamasers by a factor of 5. We find that the 183 GHz emission is systematically fainter than the 22 GH… ▽ More We present the results of an ALMA survey to identify 183 GHz H$_2$O maser emission from AGN already known to host 22 GHz megamaser systems. Out of 20 sources observed, we detect significant 183 GHz maser emission from 13; this survey thus increases the number of AGN known to host (sub)millimeter megamasers by a factor of 5. We find that the 183 GHz emission is systematically fainter than the 22 GHz emission from the same targets, with typical flux densities being roughly an order of magnitude lower at 183 GHz than at 22 GHz. However, the isotropic luminosities of the detected 183 GHz sources are comparable to their 22 GHz values. For two of our sources -- ESO 269-G012 and the Circinus galaxy -- we detect rich 183 GHz spectral structure containing multiple line complexes. The 183 GHz spectrum of ESO 269-G012 exhibits the triple-peaked structure characteristic of an edge-on AGN disk system. The Circinus galaxy contains the strongest 183 GHz emission detected in our sample, peaking at a flux density of nearly 5 Jy. The high signal-to-noise ratios achieved by these strong lines enable a coarse mapping of the 183 GHz maser system, in which the masers appear to be distributed similarly to those seen in VLBI maps of the 22 GHz system in the same galaxy and may be tracing the circumnuclear accretion disk at larger orbital radii than are occupied by the 22 GHz masers. This newly identified population of AGN disk megamasers presents a motivation for developing VLBI capabilities at 183 GHz. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Comments: 23 pages, 8 figures, submitted to ApJ

Direct measurements of carbon and sulfur isotope ratios in the Milky Way

Authors: Y. T. Yan, C. Henkel, C. Kobayashi, K. M. Menten, Y. Gong, J. S. Zhang, H. Z. Yu, K. Yang, J. J. Xie, Y. X. Wang

Abstract: With the IRAM 30 meter telescope, we performed observations of the $J$ = 2-1 transitions of CS, C$^{33}$S, C$^{34}$S, C$^{36}$S, $^{13}$CS, $^{13}$C$^{33}$S, and $^{13}$C$^{34}$S as well as the $J$ = 3-2 transitions of C$^{33}$S, C$^{34}$S, C$^{36}$S, and $^{13}$CS toward a large sample of 110 HMSFRs. We measured the $^{12}$C/$^{13}$C, $^{32}$S/$^{34}$S, $^{32}$S/$^{33}$S, $^{32}$S/$^{36}$S,… ▽ More With the IRAM 30 meter telescope, we performed observations of the $J$ = 2-1 transitions of CS, C$^{33}$S, C$^{34}$S, C$^{36}$S, $^{13}$CS, $^{13}$C$^{33}$S, and $^{13}$C$^{34}$S as well as the $J$ = 3-2 transitions of C$^{33}$S, C$^{34}$S, C$^{36}$S, and $^{13}$CS toward a large sample of 110 HMSFRs. We measured the $^{12}$C/$^{13}$C, $^{32}$S/$^{34}$S, $^{32}$S/$^{33}$S, $^{32}$S/$^{36}$S, $^{34}$S/$^{33}$S, $^{34}$S/$^{36}$S, and $^{33}$S/$^{36}$S abundance ratios with rare isotopologs of CS, thus avoiding significant saturation effects. With accurate distances obtained from parallax data, we confirm previously identified $^{12}$C/$^{13}$C and $^{32}$S/$^{34}$S gradients as a function of galactocentric distance (RGC). In the CMZ, $^{12}$C/$^{13}$C ratios are higher than suggested by a linear fit to the disk values as a function of RGC. While $^{32}$S/$^{34}$S ratios near the Galactic center and in the inner disk are similar, this is not the case for $^{12}$C/$^{13}$C, when comparing central values with those near RGC of 5 kpc. As was already known, there is no $^{34}$S/$^{33}$S gradient but the average ratio of 4.35~$\pm$~0.44 derived from the $J$ = 2-1 transition lines of C$^{34}$S and C$^{33}$S is well below previously reported values. A comparison between solar and local interstellar $^{32}$S/$^{34}$S and $^{34}$S/$^{33}$S ratios suggests that the Solar System may have been formed from gas with a particularly high $^{34}$S abundance. For the first time, we report positive gradients of $^{32}$S/$^{33}$S, $^{34}$S/$^{36}$S, $^{33}$S/$^{36}$S, and $^{32}{\rm S}/^{36}{\rm S}$ in our Galaxy. The predicted $^{12}$C/$^{13}$C ratios from the latest GCE models are in good agreement with our results. While $^{32}$S/$^{34}$S and $^{32}$S/$^{36}$S ratios show larger differences at larger RGC, $^{32}$S/$^{33}$S ratios show an offset across the entire inner 12 kpc of the Milky Way. △ Less

Submitted 27 December, 2022; v1 submitted 6 December, 2022; originally announced December 2022.

Comments: HMSFRs: high-mass star-forming regions. CMZ: central molecular zone. GCE: Galactic chemical evolution. 68 pages, 12 tables, 18 figures. Accepted for publication in A&A

Journal ref: A&A 670, A98 (2023)

H$^{13}$CN-HN$^{13}$C intensity ratio as a temperature indicator of interstellar clouds

Authors: A. G. Pazukhin, I. I. Zinchenko, E. A. Trofimova, C. Henkel

Abstract: With the 30-m IRAM radio telescope, we observed several massive star forming regions at wavelengths of 3-4 and 2 mm. The temperature of the gas in the sources was estimated from the lines of CH$_{3}$CCH and from the transitions of the NH$_3$ molecule obtained during observations at the 100-m radio telescope in Effelsberg. As a result, a correlation between the integrated intensity ratios of the… ▽ More With the 30-m IRAM radio telescope, we observed several massive star forming regions at wavelengths of 3-4 and 2 mm. The temperature of the gas in the sources was estimated from the lines of CH$_{3}$CCH and from the transitions of the NH$_3$ molecule obtained during observations at the 100-m radio telescope in Effelsberg. As a result, a correlation between the integrated intensity ratios of the $J=1-0$ transitions of H$^{13}$CN and HN$^{13}$C and the kinetic temperature has been obtained. The obtained results allow us to propose the use of the intensity ratio H$^{13}$CN-HN$^{13}$C as a possible temperature indicator of interstellar clouds. We also compared the obtained estimates of the kinetic temperature with the dust temperature $T_{dust}$. As a result, no significant correlation was found. △ Less

Submitted 25 November, 2022; originally announced November 2022.

Comments: 16 pages, 7 figures, 2 tables, Accepted for publication in Astronomy Reports (2022)

Journal ref: Astron. Rep. 66 (2022) 1302-1310

Star formation timescale in the molecular filament WB 673

Authors: O. L. Ryabukhina, M. S. Kirsanova, C. Henkel, D. S. Wiebe

Abstract: We present observations of ammonia emission lines toward the interstellar filament WB~673 hosting the dense clumps WB~673, WB~668, S233-IR and G173.57+2.43. LTE analysis of the lines allows us to estimate gas kinetic temperature ($\lesssim$ 30~K in all the clumps), number density ($7-17\times10^3$~cm$^{-3}$), and ammonia column density ($\approx 1-1.5\times 10^{15}$~cm$^{-2}$) in the dense clumps.… ▽ More We present observations of ammonia emission lines toward the interstellar filament WB~673 hosting the dense clumps WB~673, WB~668, S233-IR and G173.57+2.43. LTE analysis of the lines allows us to estimate gas kinetic temperature ($\lesssim$ 30~K in all the clumps), number density ($7-17\times10^3$~cm$^{-3}$), and ammonia column density ($\approx 1-1.5\times 10^{15}$~cm$^{-2}$) in the dense clumps. We find signatures of collapse in WB 673 and presence of compact spatially unresolved dense clumps in S233-IR. We reconstruct 1D density and temperature distributions in the clumps and estimate their ages using astrochemical modelling. Considering CO, CS, NH$_3$ and N$_2$H$^+$ molecules (plus HCN and HNC for WB~673), we find a chemical age of $t_{\rm chem}=1-3\times 10^5$~yrs providing the best agreement between the simulated and observed column densities in all the clumps. Therefore, we consider $t_{\rm chem}$ as the chemical age of the entire filament. A long preceding low-density stage of gas accumulation in the astrochemical model would break the agreement between the simulated and observed column densities. We suggest that rapid star formation over a $\sim 10^5$~yrs timescale take place in the filament. △ Less

Submitted 4 October, 2022; originally announced October 2022.

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

Discovery of non-metastable ammonia masers in Sagittarius B2

Authors: Y. T. Yan, C. Henkel, K. M. Menten, Y. Gong, H. Nguyen, J. Ott, A. Ginsburg, T. L. Wilson, A. Brunthaler, A. Belloche, J. S. Zhang, N. Budaiev, D. Jeff

Abstract: We report the discovery of widespread maser emission in non-metastable inversion transitions of NH$_3$ toward various parts of the Sagittarius B2 molecular cloud/star forming region complex: We detect masers in the $J,K = $ (6,3), (7,4), (8,5), (9,6), and (10,7) transitions toward Sgr B2(M) and Sgr B2(N), an NH$_3$ (6,3) maser in Sgr B2(NS), and NH$_3$ (7,4), (9,6), and (10,7) masers in Sgr B2(S).… ▽ More We report the discovery of widespread maser emission in non-metastable inversion transitions of NH$_3$ toward various parts of the Sagittarius B2 molecular cloud/star forming region complex: We detect masers in the $J,K = $ (6,3), (7,4), (8,5), (9,6), and (10,7) transitions toward Sgr B2(M) and Sgr B2(N), an NH$_3$ (6,3) maser in Sgr B2(NS), and NH$_3$ (7,4), (9,6), and (10,7) masers in Sgr B2(S). With the high angular resolution data of the Karl G. Jansky Very Large Array (JVLA) in A-configuration we identify 18 maser spots. Nine maser spots arise from Sgr B2(N), one from Sgr B2(NS), five from Sgr B2(M), and three in Sgr B2(S). Compared to our Effelsberg single dish data, the JVLA data indicate no missing flux. The detected maser spots are not resolved by our JVLA observations. Lower limits to the brightness temperature are $>$3000~K and reach up to several 10$^5$~K, manifesting the lines' maser nature. In view of the masers' velocity differences with respect to adjacent hot molecular cores and/or UCH{\scriptsize II} regions, it is argued that all the measured ammonia maser lines may be associated with shocks caused either by outflows or by the expansion of UCH{\scriptsize II} regions. Overall, Sgr B2 is unique in that it allows us to measure many NH$_3$ masers simultaneously, which may be essential to elucidate their so far poorly understood origin and excitation. △ Less

Submitted 6 October, 2022; v1 submitted 23 September, 2022; originally announced September 2022.

Comments: 17 pages, 14 figures, 7 tables. Accepted for publication in A&A Letter

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

Deuterated molecules in regions of high-mass star formation

Authors: Igor I. Zinchenko, Andrey G. Pazukhin, Elena A. Trofimova, Peter M. Zemlyanukha, Christian Henkel, Magnus Thomasson

Abstract: We present the results of our studies of deuterated molecules (DCN, DNC, DCO$^+$, N$_2$D$^+$ and NH$_2$D) in regions of high-mass star formation, which include a survey of such regions with the 20-m Onsala radio telescope and mapping of several objects in various lines with the 30-m IRAM and 100-m MPIfR radio telescopes. The deuteration degree reaches $\sim$10$^{-2}$ in these objects. We discuss i… ▽ More We present the results of our studies of deuterated molecules (DCN, DNC, DCO$^+$, N$_2$D$^+$ and NH$_2$D) in regions of high-mass star formation, which include a survey of such regions with the 20-m Onsala radio telescope and mapping of several objects in various lines with the 30-m IRAM and 100-m MPIfR radio telescopes. The deuteration degree reaches $\sim$10$^{-2}$ in these objects. We discuss its dependencies on the gas temperature and velocity dispersion, as well as spatial distributions of deuterated molecules. We show that the H$^{13}$CN/HN$^{13}$C intensity ratio may be a good indicator of the gas kinetic temperature and estimate densities of the investigated objects. △ Less

Submitted 16 September, 2022; originally announced September 2022.

Comments: 8 pages, 4 figures, to be published in Proceedings of Science (Proceedings of the conference "The Multifaceted Universe: Theory and Observations - 2022", 23-27 May 2022, SAO RAS, Nizhny Arkhyz, Russia)

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

Interstellar Nitrogen Isotope Ratios: New NH3 Data from the Galactic Center out to the Perseus Arm

Authors: J. L. Chen, J. S. Zhang, C. Henkel, Y. T. Yan, H. Z. Yu, J. J. Qiu, X. D. Tang, J. Wang, W. Liu, Y. X. Wang, Y. H. Zheng, J. Y. Zhao, Y. P. Zou

Abstract: Our aim is to measure the interstellar 14N/15N ratio across the Galaxy, to establish a standard data set on interstellar ammonia isotope ratios, and to provide new constraints on the Galactic chemical evolution. The (J, K ) = (1, 1), (2, 2), and (3, 3) lines of 14NH3 and 15NH3 were observed with the Shanghai Tianma 65 m radio telescope (TMRT) and the Effelsberg 100 m telescope toward a large sampl… ▽ More Our aim is to measure the interstellar 14N/15N ratio across the Galaxy, to establish a standard data set on interstellar ammonia isotope ratios, and to provide new constraints on the Galactic chemical evolution. The (J, K ) = (1, 1), (2, 2), and (3, 3) lines of 14NH3 and 15NH3 were observed with the Shanghai Tianma 65 m radio telescope (TMRT) and the Effelsberg 100 m telescope toward a large sample of 210 sources. One hundred fourty-one of these sources were detected by the TMRT in 14NH3. Eight of them were also detected in 15NH3. For 10 of the 36 sources with strong NH3 emission, the Effelsberg 100 m telescope successfully detected their 15NH3(1, 1) lines, including 3 sources (G081.7522, W51D, and Orion-KL) with detections by the TMRT telescope. Thus, a total of 15 sources are detected in both the 14NH3 and 15NH3 lines. Line and physical parameters for these 15 sources are derived, including optical depths, rotation and kinetic temperatures, and total column densities. 14N/15N isotope ratios were determined from the 14NH3/15NH3 abundance ratios. The isotope ratios obtained from both telescopes agree for a given source within the uncertainties, and no dependence on heliocentric distance and kinetic temperature is seen. 14N/15N ratios tend to increase with galactocentric distance, confirming a radial nitrogen isotope gradient. This is consistent with results from recent Galactic chemical model calculations, including the impact of superasymptotic giant branch stars and novae. △ Less

Submitted 8 August, 2022; originally announced August 2022.

ALMA imaging of the cold molecular and dusty disk in the type 2 active nucleus of the Circinus galaxy

Authors: Konrad R. W. Tristram, C. M. Violette Impellizzeri, Zhi-Yu Zhang, Eric Villard, Christian Henkel, Serena Viti, Leonard Burtscher, Françoise Combes, Santiago García-Burillo, Sergio Martín, Klaus Meisenheimer, Paul P. van der Werf

Abstract: We aim to shed light on the physical properties and kinematics of the molecular material in the nucleus of one of the closest type 2 active galaxies. To this end, we obtained high angular resolution ALMA observations of the nucleus of the Circinus galaxy. The observations map the emission at 350GHz and 690GHz with spatial resolutions of ~3.8pc and ~2.2pc, respectively. The continuum emission t… ▽ More We aim to shed light on the physical properties and kinematics of the molecular material in the nucleus of one of the closest type 2 active galaxies. To this end, we obtained high angular resolution ALMA observations of the nucleus of the Circinus galaxy. The observations map the emission at 350GHz and 690GHz with spatial resolutions of ~3.8pc and ~2.2pc, respectively. The continuum emission traces cold ($T\lesssim100$K) dust in a circumnuclear disk with spiral arms on scales of 25pc, plus a marginally resolved nuclear emission peak. The latter is not extended in polar direction as claimed based on earlier ALMA observations. A significant amount (of the order of 40%) of the 350GHz emission is not related to dust, but most likely free-free emission instead. We detect CO(3-2) and CO(6-5) as well as HCO$^+$(4-3), HCN(4-3), and CS(4-3). The CO emission is extended, showing a spiral pattern, similar to the extended dust emission. Towards the nucleus, CO is excited to higher transitions and its emission is self-absorbed, leading to an apparent hole in the CO(3-2) but not the CO(6-5) emission. On the other hand, the high gas density tracers HCO$^+$, HCN, and CS show a strong, yet unresolved (($\lesssim4$pc) concentration of the emission at the nucleus, pointing at a very small 'torus'. The kinematics are dominated by rotation and point at a geometrically thin disk down to the resolution limit of our observations. In contrast to several other AGNs, no HCN enhancement is found towards the nucleus. The Circinus nucleus is therefore composed of at least two distinct components: (1) an optically thin, warm outflow of ionised gas containing clouds of dust; and (2) a cold molecular and dusty disk. These findings support the most recent radiative transfer calculations of the obscuring structures in AGNs, which find a similar two-component structure. (Abridged) △ Less

Submitted 1 June, 2022; originally announced June 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 664, A142 (2022)

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

The Interstellar Medium of Dwarf Galaxies

Authors: Christian Henkel, Leslie K. Hunt, Yuri I. Izotov

Abstract: Dwarf galaxies are by far the most numerous galaxies in the Universe, showing properties that are quite different from those of their larger and more luminous cousins. This review focuses on the physical and chemical properties of the interstellar medium of those dwarfs that are known to host significant amounts of gas and dust. The neutral and ionized gas components and the impact of the dust wil… ▽ More Dwarf galaxies are by far the most numerous galaxies in the Universe, showing properties that are quite different from those of their larger and more luminous cousins. This review focuses on the physical and chemical properties of the interstellar medium of those dwarfs that are known to host significant amounts of gas and dust. The neutral and ionized gas components and the impact of the dust will be discussed, as well as first indications for the existence of active nuclei in these sources. Cosmological implications are also addressed, considering the primordial helium abundance and the similarity of local Green Pea galaxies with young, sometimes proto-galactic sources in the early Universe. △ Less

Submitted 16 February, 2022; originally announced February 2022.

Comments: 29 pages, 10 figures, review article in Galaxies, volume 10

Journal ref: 2022, Galaxies, 10, 11

Redshifted methanol absorption tracing infall motions of high-mass star formation regions

Authors: W. J. Yang, K. M. Menten, A. Y. Yang, F. Wyrowski, Y. Gong, S. P. Ellingsen, C. Henkel, X. Chen, Y. Xu

Abstract: Gravitational collapse is one of the most important processes in high-mass star formation. Compared with the classic blue-skewed profiles, redshifted absorption against continuum emission is a more reliable method to detect inward motions within high-mass star formation regions. We aim to test if methanol transitions can be used to trace infall motions within high-mass star formation regions. Usin… ▽ More Gravitational collapse is one of the most important processes in high-mass star formation. Compared with the classic blue-skewed profiles, redshifted absorption against continuum emission is a more reliable method to detect inward motions within high-mass star formation regions. We aim to test if methanol transitions can be used to trace infall motions within high-mass star formation regions. Using the Effelsberg-100 m, IRAM-30 m, and APEX-12 m telescopes, we carried out observations of 37 and 16 methanol transitions towards two well-known collapsing dense clumps, W31C (G10.6-0.4) and W3(OH), to search for redshifted absorption features or inverse P-Cygni profiles. Redshifted absorption is observed in 14 and 11 methanol transitions towards W31C and W3(OH), respectively. The infall velocities fitted from a simple two-layer model agree with previously reported values derived from other tracers, suggesting that redshifted methanol absorption is a reliable tracer of infall motions within high-mass star formation regions. Our observations indicate the presence of large-scale inward motions, and the mass infall rates are roughly estimated to be $\gtrsim$10$^{-3}$ $M_{\odot}$yr$^{-1}$, which supports the global hierarchical collapse and clump-fed scenario. With the aid of bright continuum sources and the overcooling of methanol transitions leading to enhanced absorption, redshifted methanol absorption can trace infall motions within high-mass star formation regions hosting bright H{\scriptsize II} regions. △ Less

Submitted 5 January, 2022; originally announced January 2022.

Comments: 21 pages, 13 figures, 5 tables. Accepted for publication in A&A

Discovery of ammonia (9,6) masers in two high-mass star-forming regions

Authors: Y. T. Yan, C. Henkel, K. M. Menten, Y. Gong, J. Ott, T. L. Wilson, A. Wootten, A. Brunthaler, J. S. Zhang, J. L. Chen, K. Yang

Abstract: Molecular maser lines are signposts of high-mass star formation, probing excitation and kinematics of very compact regions in the close environment of young stellar objects and providing useful targets for trigonometric parallax measurements. Only a few NH$_{3}$ (9,6) masers were known so far, and their origin is still poorly understood. Here we aim to find new NH$_{3}$ (9,6) masers to provide a b… ▽ More Molecular maser lines are signposts of high-mass star formation, probing excitation and kinematics of very compact regions in the close environment of young stellar objects and providing useful targets for trigonometric parallax measurements. Only a few NH$_{3}$ (9,6) masers were known so far, and their origin is still poorly understood. Here we aim to find new NH$_{3}$ (9,6) masers to provide a better observational basis to study their role in high-mass star-forming regions. We carried out NH$_{3}$ (9,6) observations toward Cepheus A and G34.26$+$0.15 with the Effelsberg-100 m telescope and the Karl G. Janksy Very Large Array. We discovered new NH$_{3}$ (9,6) masers in Cep A and G34.26$+$0.15, which increases the number of high-mass star-forming regions hosting NH$_{3}$ (9,6) masers from five to seven. Long term monitoring (20 months) at Effelsberg shows that the intensity of the (9,6) maser in G34.26$+$0.15 is decreasing, while the Cep A maser remains stable. Compared to the Effelsberg data and assuming linear variations between the epochs of observation, the JVLA data indicate no missing flux. This suggests that the NH$_3$ (9,6) emission arises from single compact emission regions that are not resolved by the interferometric measurements. As JVLA imaging shows, the NH$_{3}$ (9,6) emission in Cep A originates from a sub-arcsecond sized region, slightly to the west of the peak position of the 1.36\,cm continuum object, HW2. In G34.26$+$0.15, three NH$_{3}$ (9,6) maser spots are observed: one is close to the head of the cometary ultracompact \h2 region C and the other two are emitted from a compact region to the west of the hypercompact \h2 region A. The newly found (9,6) masers appear to be related to outflows. Higher angular resolution of JVLA and VLBI observations are needed to provide more accurate positions and constraints for pumping scenarios. △ Less

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

Comments: 10 pages, 6 figures, 3 tables. Publicated in A&A

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

Probing the electron-to-proton mass ratio gradient in the Milky Way with class I methanol masers

Authors: S. A. Levshakov, I. I. Agafonova, C. Henkel, Kee-Tae Kim, M. G. Kozlov, B. Lankhaar, W. Yang

Abstract: We estimate limits on non-universal coupling of hypothetical hidden fields to standard matter by evaluating the fractional changes in the electron-to-proton mass ratio, mu = m_e/m_p, based on observations of ClassI methanol masers distributed in the Milky Way disk over the range of the galactocentric distances 4 < R < 12 kpc. The velocity offsets DeltaV = V44 - V95 measured between the 44 and 95 G… ▽ More We estimate limits on non-universal coupling of hypothetical hidden fields to standard matter by evaluating the fractional changes in the electron-to-proton mass ratio, mu = m_e/m_p, based on observations of ClassI methanol masers distributed in the Milky Way disk over the range of the galactocentric distances 4 < R < 12 kpc. The velocity offsets DeltaV = V44 - V95 measured between the 44 and 95 GHz methanol lines provide, so far, one of the most stringent constraints on the spatial gradient k_mu = d(Delta mu/mu)/dR < 2x10^-9 kpc-1 and the upper limit on Delta mu/mu < 2x10^-8, where Delta mu/mu = (mu_obs-mu_lab)/mu_lab. We also find that the offsets DeltaV are clustered into two groups which are separated by 0.022 +/- 0.003 km/s (1sigma C.L.). The grouping is most probably due to the dominance of different hyperfine transitions in the 44 and 95 GHz methanol maser emission. Which transition becomes favored is determined by an alignment (polarization) of the nuclear spins of the four hydrogen atoms in the methanol molecule. This result confirms that there are preferred hyperfine transitions involved in the methanol maser action. △ Less

Submitted 29 December, 2021; originally announced December 2021.

Comments: Accepted for publication in MNRAS, 7 pages text, 8 figures, 7 tables

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)

Observations of multiple NH$_3$ transitions in W33

Authors: K. Tursun, C. Henkel, J. Esimbek, X. D. Tang, T. L. Wilson, A. Malawi, E. Alkhuja, F. Wyrowski, R. Mauersberger, K. Immer, H. Asiri, J. J. Zhou, G. Wu

Abstract: At a distance of 2.4kpc, W33 is an outstanding massive and luminous 10pc sized star forming complex containing quiescent infrared dark clouds as well as highly active infrared bright cloud cores heated by young massive stars. We report measurements of ammonia (NH$_3$) inversion lines in the frequency range 18--26GHz, obtained with the 40" resolution of the 100 m Effelsberg telescope. We have detec… ▽ More At a distance of 2.4kpc, W33 is an outstanding massive and luminous 10pc sized star forming complex containing quiescent infrared dark clouds as well as highly active infrared bright cloud cores heated by young massive stars. We report measurements of ammonia (NH$_3$) inversion lines in the frequency range 18--26GHz, obtained with the 40" resolution of the 100 m Effelsberg telescope. We have detected the ($J$, $K$)=(1,1), (2,2), (3,3), (4,4), (5,5), (6,6), (2,1) and (3,2) transitions. There is a maser line in the (3,3) transition towards W33 Main. Brightness temperature and line shape indicate no significant variation during the last $\sim$36yr. We have determined kinetic temperatures, column densities and other physical properties of NH$_3$ and the molecular clouds in W33. For the total-NH$_3$ column density, we find for 40"(0.5pc) sized regions 6.0($\pm$2.1)$\times$10$^{14}$, 3.5($\pm$0.1)$\times$10$^{15}$, 3.4($\pm$0.2)$\times$10$^{15}$, 3.1($\pm$0.2)$\times$10$^{15}$, 2.8($\pm$0.2)$\times$10$^{15}$ and 2.0($\pm$0.2)$\times$10$^{15}$cm$^{-2}$ at the peak positions of W33 Main, W33 A, W33 B, W33 Main1, W33 A1 and W33 B1, respectively. W33 Main has a total-NH$_3$ fractional abundance of 1.3($\pm$0.1)$\times$10$^{-9}$ at the peak position. High values of 1.4($\pm$0.3)$\times$10$^{-8}$, 1.6($\pm$0.3)$\times$10$^{-8}$, 3.4($\pm$0.5)$\times$10$^{-8}$, 1.6($\pm$0.5)$\times$10$^{-8}$ and 4.0($\pm$1.2)$\times$10$^{-8}$ are obtained at the central positions of W33 A, W33 B, W33 Main1, W33 A1, and W33 B1. From this, we confirm the already previously proposed different evolutionary stages of the six W33 clumps and find that there is no hot core in the region approaching the extreme conditions encountered in W51-IRS2 or Sgr B2. The ortho-to-para-NH$_3$ abundance ratios suggest that ammonia should have been formed in the gas phase or on dust grain mantles at kinetic temperatures of $\gtrsim$20K. △ Less

Submitted 11 November, 2021; originally announced November 2021.

Comments: 30 pages, 21 figures, accepted by A&A

Journal ref: A&A 658, A34 (2022)