Distributions of the Density and Kinetic Temperature of the Molecular Gas in the Central Region of NGC 613 using Hierarchical Bayesian Inference

TL;DR

This study uses hierarchical Bayesian inference on ALMA molecular line data to map the physical conditions of molecular gas in NGC 613's center, revealing variations in star formation efficiency linked to gas distribution.

Contribution

It applies a novel non-LTE hierarchical Bayesian method to external galaxy data, providing detailed gas property maps and insights into star formation efficiency variations.

Findings

Derived gas densities, temperatures, and column densities consistent with previous studies.

Identified two regions with different star formation efficiencies.

Linked gas deficiency to higher star formation efficiency in one region.

Abstract

We present position-position-velocity (PPV) cubes of the physical and chemical properties of the molecular medium in the central 1.2 kpc region of the active galaxy NGC 613 at a PPV resolution of 0.$^{\prime\prime}$8$\times$0.$^{\prime\prime}$8$\times$10 km s$^{-1}$ (0.$^{\prime\prime}$8 = $\sim$68 pc). We used eight molecular lines obtained with ALMA. Non-LTE calculation with hierarchical Bayesian inference was used to construct PPV cubes of the gas kinetic temperature ($T_\mathrm{kin}$), molecular hydrogen volume density ($n_\mathrm{H_2}$), column densities ($N_\mathrm{H_2}$), and fractional abundances of four molecules ($^{12}$C$^{18}$O, HCN, HCO$^+$, and CS). The derived $n_\mathrm{H_2}$, $N_\mathrm{H_2}$, and $T_\mathrm{kin}$ ranged 10$^{3.21-3.85}$ cm$^{-3}$, 10$^{20.8-22.1}$ cm$^{-2}$, and 10$^{2.33-2.64}$ K, respectively. Our first application of the non-LTE method with the hierarchical Bayesian inference to external galaxies yielded compatible results compared with the previous studies of this galaxy, demonstrating the efficacy of this method for application to other galaxies. We examined the correlation between gas surface density $\Sigma_\mathrm{H_2}$ (converted from $N_\mathrm{H_2}$) and the star formation rate $\Sigma_\mathrm{SFR}$ obtained from the 110 GHz continuum flux map and found two distinct sequences in the $\Sigma_\mathrm{H_2}$-$\Sigma_\mathrm{SFR}$ diagram; the southwestern subregion of the star-forming ring exhibited a $\sim$0.5 dex higher star formation efficiency (SFE; $\Sigma_\mathrm{SFR}/\Sigma_\mathrm{H_2}$) than the eastern subregion. However, they exhibited no systematic difference in $n_\mathrm{H_2}$, which is often argued as a driver of SFE variation. We suggest that the deficiency of molecular gas in the southwestern subregion, where no significant gas supply is evident along the offset ridges in the bar, is responsible for the elevated SFE.