Akio Taniguchi

We present a new framework for multiline analysis that combines kernel principal component analysis (Kernel PCA), an unsupervised machine learning method, and Kernel Shapley Additive Explanations (SHAP), an explainable artificial intelligence technique. To enable a comparison with PCA-based studies, which have been widely used in multiline analyses, we apply our framework to integrated intensity maps of 13 molecular lines from Atacama Large Millimeter/submillimeter Array Band 3 archival data of the nearby galaxy NGC 1068. Previous PCA-based studies of NGC 1068 reported that physically meaningful structures are mainly captured up to the second component. In contrast, our framework suggests that the features up to the third component are physically meaningful, whereas the physical meaning of the fourth component still requires further discussion. Furthermore, by comparing the results obtained from our framework with molecular column densities derived from local thermodynamical equilibrium analysis, we suggest that the abundance of HCO ⁺ may be relatively enhanced in the molecular outflow region extending to a radius of about 400 pc from the galactic center, likely due to the effects of ultraviolet radiation and highly dense gas. These results show that our framework can provide data-driven insights into physical and chemical features that have not been clearly identified in previous studies. It also provides an efficient tool for interpreting the rapidly increasing amount of multiline observational data.

We present a new framework for multi-line analysis that combines kernel principal component analysis (Kernel PCA), an unsupervised machine-learning method, and Kernel SHapley Additive exPlanations (Kernel SHAP), an explainable artificial intelligence (XAI) technique. To enable a comparison with PCA-based studies, which have been widely used in multi-line analyses, we apply our framework to integrated intensity maps of 13 molecular lines from Atacama Large Millimeter/submillimeter Array (ALMA) Band 3 archival data of the nearby galaxy NGC 1068. Previous PCA-based studies of NGC 1068 reported that physically meaningful structures are mainly captured up to the second component. In contrast, our framework can interpret physically meaningful features up to the fourth component. Furthermore, by comparing the results obtained from our framework with molecular column densities derived from local thermodynamical equilibrium (LTE) analysis, we suggest that the abundance of HCO+ is relatively enhanced in the molecular outflow region extending to a radius of about 400 pc from the galactic center, likely due to the effects of ultraviolet radiation and highly dense gas. These results show that our framework can provide data-driven insights into physical and chemical features that have not been clearly identified in previous studies. It also provides an efficient tool for interpreting the rapidly increasing amount of multi-line observational data.

We report on the results of the on-sky test and science demonstration conducted with the 2 mm receiver system, B4R, on the 50 m Large Millimeter Telescope (LMT), located at an altitude of 4600 m in Mexico. The B4R receiver was developed based on the dual-polarization sideband-separating mixer technology of the Atacama Large Millimeter/submillimeter Array and is equipped with a fast Fourier transform digital spectrometer, XFFTS. The primary science objective is the spectroscopic redshift identification of high-redshift dusty star-forming galaxies, complementing the existing 3 mm Redshift Search Receiver by enabling the detection of multiple carbon monoxide lines. Additionally, the B4R receiver broadens the range of science cases possible with the LMT, including astrochemistry, as the 2 mm band encompasses unique molecular lines such as deuterated molecules and shock tracers. During on-site commissioning in 2018 and 2019, we successfully demonstrated on-the-fly mapping and position-switching observations toward the Orion Molecular Cloud 1 and bright high-redshift dusty star-forming galaxies, respectively. We confirmed that the installed B4R system largely met its basic performance specifications. Furthermore, we measured the LMT’s aperture efficiencies across the entire B4R frequency range (130–160 GHz), finding them to be roughly consistent with expectations based on a surface accuracy of 100 μ m and the receiver optics design. These results with the B4R will enable the most sensitive single-dish spectroscopic observations at 2 mm using the LMT.

Porous interstellar medium (ISM) structure in galaxies at the epoch of reionization gives us a hint to understand what types of galaxies contribute to reionization. Although recent studies have pointed out the positive correlation between high ionizing photon escape fractions and high [O iii] 88 μm/[C ii] 158 μm ratios found in UV-luminous star-forming galaxies at z > 6 with the Atacama Large Millimeter/submillimeter Array, previous studies have paid little attention to the neutral gas porosity that allows ionizing photons to escape. Here we present a detailed analysis of a z = 8.312 Lyman break galaxy, MACS0416_Y1 with a high L _{[O III] 88 μm}/L _{[C II] 158 μm} ratio (≈9) and dust continuum detection. We construct a multiphase ISM model incorporating the neutral gas covering fraction (c ov _PDR). The best-fit model reveals a cov _PDR ≈ 25%, indicating that ≈75% of the ionized gas region is exposed to intercloud space. We confirm that our conclusions hold even when varying star formation history, stellar age, gas/stellar metallicity, and carbon-to-oxygen abundance ratio. This finding meets one of the necessary conditions for galaxies to have a nonzero escape fraction of ionizing photons and supports recent studies that galaxies with a high [O iii] 88 μm/[C ii] 158 μm ratio, such as MACS0416_Y1, could contribute to cosmic reionization. Furthermore, the modeled H ii region with the best-fitting parameters has a typical size (D = 0.90 pc) and gas density ( log n H , c / cm − 3 = 2.7 ) that are comparable to local compact H ii regions. This suggests that the H ii regions in MACS0416_Y1 are in an early evolutionary stage.

We report the Australian Telescope Compact Array and Nobeyama 45 m telescope detection of a remarkably bright (S _1.1mm = 44 mJy) submillimeter galaxy MM J154506.4−344318 in emission lines at 48.5 and 97.0 GHz, respectively. We also identify part of an emission line at ≈218.3 GHz using the Atacama Large Millimeter/submillimeter Array (ALMA). Together with photometric redshift estimates and the ratio between the line and infrared luminosities, we conclude that the emission lines are most likely to be the J = 2–1, 4–3, and 9–8 transitions of ¹²CO at redshift z = 3.753 ± 0.001. ALMA 1.3 mm continuum imaging reveals an arc and a spot separated by an angular distance of 1 . ″ 6 , indicative of a strongly lensed dusty star-forming galaxy with respective molecular and dust masses of log M mol / M ⊙ ≈ 11.5 and log M dust / M ⊙ ≈ 9.4 after being corrected for ≈6.6× gravitational magnification. The inferred dust-to-gas mass ratio is found to be high (≈0.0083) among coeval dusty star-forming galaxies, implying the presence of a massive, chemically enriched reservoir of cool interstellar medium at z ≈ 4, or 1.6 Gyr after the Big Bang.

The star formation rate (SFR) is a fundamental parameter for describing galaxies and inferring their evolutionary course. H ii regions yield the best measure of instantaneous SFR in galaxies, although the derived SFR can have large uncertainties depending on tracers and assumptions. We present an SFR calibration for the entire molecular gas disk of the nearby Seyfert galaxy NGC 1068, based on our new high-sensitivity Atacama Large Millimeter/submillimeter Array 100 GHz continuum data at 55 pc (= 0.″8) resolution in combination with the Hubble Space Telescope Paα line data. In this calibration, we account for the spatial variations of dust extinction, electron temperature of H ii regions, AGN contamination, and diffuse ionized gas (DIG) based on publicly available multiwavelength data. Especially, given the extended nature and the possible nonnegligible contribution to the total SFR, a careful consideration of DIG is essential. With a cross-calibration between two corrected ionized gas tracers (free–free continuum and Paα), the total SFR of the NGC 1068 disk is estimated to be 3.2 ± 0.5 M _⊙ yr⁻¹, one-third of the SFR without accounting for DIG (9.1 ± 1.4 M _⊙ yr⁻¹). We confirmed a high SFR around the southern bar end and the corotation radius, which is consistent with the previous SFR measurements. In addition, our total SFR exceeds the total SFR based on 8 μm dust emission by a factor of 1.5. We attribute this discrepancy to the differences in the young stars at different stages of evolution traced by each tracer and their respective timescales. This study provides an example to address the various uncertainties in conventional SFR measurements and their potential to lead to significant SFR miscalculations.

Star formation rate (SFR) is a fundamental parameter for describing galaxies and inferring their evolutionary course. HII regions yield the best measure of instantaneous SFR in galaxies, although the derived SFR can have large uncertainties depending on tracers and assumptions. We present an SFR calibration for the entire molecular gas disk of the nearby Seyfert galaxy NGC 1068, based on our new high-sensitivity ALMA 100GHz continuum data at 55pc (=0."8) resolution in combination with the HST Pa{\alpha} line data. In this calibration, we account for the spatial variations of dust extinction, electron temperature of HII regions, AGN contamination, and diffuse ionized gas (DIG) based on publicly available multi-wavelength data. Especially, given the extended nature and the possible non-negligible contribution to the total SFR, a careful consideration of DIG is essential. With a cross-calibration between two corrected ionized gas tracers (free-free continuum&Pa{\alpha}), the total SFR of the NGC 1068 disk is estimated to be 3.2\pm0.5 Msol/yr, one-third of the SFR without accounting for DIG (9.1\pm1.4 Msol/yr). We confirmed high SFR around the southern bar-end and the corotation radius, which is consistent with the previous SFR measurements. In addition, our total SFR exceeds the total SFR based on 8{\mu}m dust emission by a factor of 1.5. We attribute this discrepancy to the differences in the young stars at different stages of evolution traced by each tracer and their respective timescales. This study provides an example to address the various uncertainties in conventional SFR measurements and their potential to lead to significant SFR miscalculations.

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