Revisiting the black hole mass of M87* using VLT/MUSE Adaptive Optics Integral Field Unit data I: Ionized gas kinematics

The stellar dynamic-based black hole mass measurements of M87 are twice that determined via ionized gas kinematics; the former is closer to the estimation from the diameter of the gravitationally-lensed ring around the black hole. Using deeper and more comprehensive ionized gas kinematic data, we aim to better constrain the morphology and kinematics of the nuclear ionized gas, thus gaining insights into the reasons behind the disagreement of the measurements. We use both Narrow and Wide Field Mode integral field spectroscopic data from the Multi Unit Spectroscopic Explorer instrument, to model the morphology and kinematics of multiple ionized gas emission lines in the nucleus of M87. The new deep dataset reveals complexities in the nuclear ionized gas kinematics. Several ionized gas filaments can be traced down into the projected sphere of influence. We also found evidence of a partially-filled biconical outflow. The velocity isophotes of the ionized gas disk are twisted and the position angle of the innermost gas disk tends toward a value perpendicular to the radio jet axis. The complexity of the nuclear morphology and kinematics precludes the measurement of an accurate black hole mass. The results support a 6.0 $\times 10^{9}\rm M_{\odot}$ black hole in a 25\deg disk, rather than a 3.5 $\times 10^{9}\rm M_{\odot}$ black hole in a 42\deg disk. The specific RIAF model earlier proposed to reconcile the mass measurement discrepancy was also tested. In general, Keplerian disk models perform better than the RIAF model when fitting the sub-arcsec ionized gas disk. A disk inclination close to 25\deg for the nuclear gas disk, and the warp in the sub-arcsec ionized gas disk, help to reconcile the contradictory nature of the mass discrepancy between stellar and ionized gas black hole masses, and the mis-orientation between the axes of the ionized gas disk and the jet.Comment: 21 pages, 22 figures (5 of them in the appendix). Accepted in Astronomy & Astrophysic

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

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

Towards the prediction of molecular parameters from astronomical emission lines using Neural Networks

Molecular astronomy is a field that is blooming in the era of large observatories such as the Atacama Large Millimeter/Submillimeter Array (ALMA). With modern, sensitive, and high spectral resolution radio telescopes like ALMA and the Square Kilometer Array, the size of the data cubes is rapidly escalating, generating a need for powerful automatic analysis tools. This work introduces MolPred, a pilot study to perform predictions of molecular parameters such as excitation temperature (Tex) and column density (log(N)) from input spectra by the use of neural networks. We used as test cases the spectra of CO, HCO+, SiO and CH3CN between 80 and 400 GHz. Training spectra were generated with MADCUBA, a state-of-the-art spectral analysis tool. Our algorithm was designed to allow the generation of predictions for multiple molecules in parallel. Using neural networks, we can predict the column density and excitation temperature of these molecules with a mean absolute error of 8.5% for CO, 4.1% for HCO+, 1.5% for SiO and 1.6% for CH3CN. The prediction accuracy depends on the noise level, line saturation, and number of transitions. We performed predictions upon real ALMA data. The values predicted by our neural network for this real data differ by 13% from the MADCUBA values on average. Current limitations of our tool include not considering linewidth, source size, multiple velocity components, and line blending

Methanol masers in NGC 253 with ALCHEMI

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 ($\lambda$ = 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.Comment: Accepted for publication in A&A. 29 pages, 17 figures (4 in Appendix

The distribution and origin of C₂H in NGC 253 from ALCHEMI

Context. Observations of chemical species can provide insights into the physical conditions of the emitting gas however it is important to understand how their abundances and excitation vary within different heating environments. C2H is a molecule typically found in PDR regions of our own Galaxy but there is evidence to suggest it also traces other regions undergoing energetic processing in extragalactic environments. / Aims. As part of the ALCHEMI ALMA large program, we map the emission of C2H in the central molecular zone of the nearby starburst galaxy NGC 253 at 1.6″ (28 pc) resolution and characterize it to understand its chemical origins. / Methods. We used spectral modeling of the N = 1−0 through N = 4−3 rotational transitions of C2H to derive the C2H column densities towards the dense clouds in NGC 253. We then use chemical modeling, including photodissociation region (PDR), dense cloud, and shock models to investigate the chemical processes and physical conditions that are producing the molecular emission. / Results. We find high C2H column densities of ∼1015 cm−2 detected towards the dense regions of NGC 253. We further find that these column densities cannot be reproduced if it is assumed that the emission arises from the PDR regions at the edge of the clouds. Instead, we find that the C2H abundance remains high even in the high visual extinction interior of these clouds and that this is most likely caused by a high cosmic-ray ionization rate

ALCHEMI Finds a “Shocking” Carbon Footprint in the Starburst Galaxy NGC 253

The 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 analyse the abundances of CO2, a reliable tracer of ice chemistry, from data collected as part of the Atacama Large Millimeter/submillimeter Array large program ALCHEMI, a wide-frequency spectral scan toward the starburst galaxy NGC 253 with an angular resolution of 1.″6. We constrain the CO2 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 CO2 sublimated from ice. The HOCO+ fractional abundances are found to be (1-2) 7 10−9 at the outer part of the central molecular zone (CMZ), while they are lower (∼10−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+/CO2 ratios suggested from chemical models, the gas-phase CO2 fractional abundance is estimated to be (1-20) 7 10−7 at the outer CMZ, and orders of magnitude lower near the center. We estimate the CO2 ice fractional abundances at the outer CMZ to be (2-5) 7 10−6 from the literature. A comparison between the ice and gas CO2 abundances suggests an efficient sublimation mechanism. This sublimation is attributed to large-scale shocks at the orbital intersections of the bar and CMZ

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

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 the Atacama Large Millimeter/submillimeter Array via the ALCHEMI Large Program. To understand the conditions and associated heating mechanisms within NGC 253\u27s 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 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\u27s central giant molecular clouds (GMCs), where observed HCN/HNC abundance ratios are the lowest, n ∼ 105.5 cm−3 and ζ ∼ 10−12 s−1 (greater than 104 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, H ii 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\u27s 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

Starburst Energy Feedback Seen through HCO+/HOC+Emission in NGC 253 from ALCHEMI

Molecular abundances are sensitive to the UV photon flux and cosmic-ray ionization rate. In starburst environments, the effects of high-energy photons and particles are expected to be stronger. We examine these astrochemical signatures through multiple transitions of HCO+ and its metastable isomer HOC+ in the center of the starburst galaxy NGC 253 using data from the Atacama Large Millimeter/submillimeter Array large program ALMA Comprehensive High-resolution Extragalactic Molecular inventory. The distribution of the HOC+(1-0) integrated intensity shows its association with "superbubbles,"cavities created either by supernovae or expanding H ii regions. The observed HCO+/HOC+ abundance ratios are ∼10-150, and the fractional abundance of HOC+ relative to H2 is ∼1.5 × 10-11-6 × 10-10, which implies that the HOC+ abundance in the center of NGC 253 is significantly higher than in quiescent spiral arm dark clouds in the Galaxy and the Galactic center clouds. Comparison with chemical models implies either an interstellar radiation field of G 0 ⪆ 103 if the maximum visual extinction is ⪆5, or a cosmic-ray ionization rate of ζ ⪆ 10-14 s-1 (3-4 orders of magnitude higher than that within clouds in the Galactic spiral arms) to reproduce the observed results. From the difference in formation routes of HOC+, we propose that a low-excitation line of HOC+ traces cosmic-ray dominated regions, while high-excitation lines trace photodissociation regions. Our results suggest that the interstellar medium in the center of NGC 253 is significantly affected by energy input from UV photons and cosmic rays, sources of energy feedback.N.H. acknowledges support from JSPS KAKENHI grant No. JP21K03634. K.S. has been supported by grants MOST 108-2112-M-001-015 and 109- 2112-M-001-020 from the Ministry of Science and Technology, Taiwan. Y.N. is supported by the NAOJ ALMA Scientific Research grant No. 2017-06B. V.M.R. and L.C. are funded by the Comunidad de Madrid through the Atracción de Talento Investigador (Doctores con experiencia) Grant (COOL: Cosmic Origins Of Life; 2019-T1/TIC-15379)

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

Galaxie

ALCHEMI, an ALMA comprehensive high-resolution extragalactic molecular inventory: survey presentation and first results from the ACA array

Galaxie