Dynamical Evidence for a Billion Solar Mass Black Hole in Galaxy NGC 4061 from ALMA $^{12}$CO(2-1) Kinematics

TL;DR

The paper presents a robust dynamical measurement of the SMBH mass in the massive ETG NGC 4061 using high-resolution ALMA CO(2-1) kinematics, achieving a synthesized beam close to the BH sphere of influence. Through KinMS forward modeling with data-driven SkySampler and analytic Gaussian gas distributions, combined with an HST-based MGE stellar mass model and an ISM mass component, the authors derive MBH = $\left(1.17^{+0.08}_{-0.10}\,\text{stat.} \right)\times 10^{9}$ M$_{\odot}$ and M/L$_{\rm F814W}$ = $3.46^{+0.07}_{-0.06}$ M$_{\odot}$/L$_{\odot}$, with a careful accounting of systematic uncertainties. The results resolve previous discrepancies from σ-based estimates and illustrate that molecular-gas dynamics yield precise MBH measurements at the high-mass end of the local black hole mass function. The study finds MBH to be broadly consistent with the M–σ relation (and with M_BH–L_K bulge) derived from purely dynamical and molecular-gas samples, while noting some tension with relations that include reverberation-mapped masses. Overall, the work highlights the efficacy of high-resolution ALMA observations for calibrating SMBH demographics and their co-evolution with massive galaxies.

Abstract

We present the first robust dynamical measurement of the supermassive black hole (SMBH) mass in the massive early-type galaxy NGC 4061 using high-spatial-resolution ALMA observations of the $^{12}$CO(2-1) emission. By combining archival Cycle 6 data with new Cycle 7 observations, we achieve a synthesized beam of $0''.16 \times 0''.13$, comparable to the expected sphere of influence of the central black hole. The molecular gas forms a regularly rotating circumnuclear disk aligned with the prominent dust lane seen in HST imaging. We model the full three-dimensional ALMA data cube using the KinMS forward-modeling framework, exploring both data-driven and analytic prescriptions for the gas surface brightness distribution. Our Bayesian analysis yields a best-fitting SMBH mass of $M_{\rm BH} = (1.17^{+0.08}_{-0.10}\,[{\rm stat.}] \pm 0.43\,[{\rm syst.}]) \times 10^{9}$ M$_\odot$ and an $I$-band stellar mass-to-light ratio of $M/L_{\rm F814W} = 3.46^{+0.07}_{-0.06}\,[{\rm stat.}] \pm 0.10\,[{\rm syst.}]$ M$_\odot$/L$_\odot$. The inferred black hole mass is fully consistent across different modeling assumptions and remains insensitive to plausible radial variations in the $M/L_{\rm F814W}$ profile. Our results resolve the long-standing discrepancy between previous indirect mass estimates based on conflicting stellar velocity dispersion measurements and demonstrate that the exceptionally large dispersion reported in the literature is likely spurious. This study highlights the power of high-resolution ALMA molecular gas kinematics for precision SMBH mass measurements at the high-mass end of the local black hole mass function.

Figures (16)

• Figure 1: Left: HST/WFPC2 F814W image of NGC 4061 within a 25$\times$ 25 (or 13 kpc $\times$ 13 kpc) field of view, showing the central dust lane. Right: ALMA $^{12}$CO(2$-$1) integrated intensity contours overlaid on the left zoom-in HST image (7$\times$ 7 or 3.64 kpc $\times$ 3.64 kpc), illustrating the alignment of the molecular gas with the dust lane.
• Figure 2: Panel A: 1.3 mm continuum emission in Band 6. Panels B–D: ALMA $^{12}$CO(2$-$1) moment maps of NGC $4061$—integrated intensity, intensity-weighted mean LOS velocity, and intensity-weighted LOS velocity dispersion, respectively. The synthesized beam is shown as a black ellipse in the lower-left corner of each map. Panel E: Integrated spectrum extracted from a $6\arcsec \times 6\arcsec$ ($3.12 \times 3.12$ kpc) box; the horizontal dot-dashed line marks zero flux. Panel F: Position–velocity diagram along the major axis, adopting a systemic velocity of $v_{\rm sys}=7190$ km s$^{-1}$ and position angle $\Gamma=175\degr$.
• Figure 3: Comparison between the molecular gas mass distribution derived from the zeroth-moment map (panel B of Figure \ref{['fig:moment-maps']}) and its MGE model (contours). The strong correspondence between the data and model demonstrates excellent agreement across matching radii and contour levels.
• Figure 4: Comparison between the HST/WFPC2 F814W image and its MGE model, shown in 2D surface brightness density over the field of $70\arcsec \times 70\arcsec$ (left) and a central zoom ($7\arcsec \times 7\arcsec$; right). Black contours indicate the data, and red contours show the model, demonstrating close agreement across radii and contour levels. The yellow region marks masked areas affected by bad pixels and the central dust disk.
• Figure 5: The $^{12}$CO(2$-$1)-CND surface brightness distribution of NGC 4061 is best described by a simple Gaussian model. The ALMA data are plotted in black, and the model is overlaid in blue.