Measuring the Central Dark Mass in NGC 4258 with JWST/NIRSpec Stellar Kinematics

TL;DR

This study measures the central supermassive black hole mass in NGC 4258 using high-resolution 2D stellar kinematics from JWST/NIRSpec. By applying Jeans Anisotropic Models to a grid of 12 runs that vary PSF, M/L_K, and orbital anisotropy, the authors robustly quantify systematic uncertainties and validate the stellar-dynamical mass against the gold-standard maser result. They report an ensemble median M_BH = (4.08^{+0.19}_{-0.33})×10^7 M⊙, with a 5% offset from the maser mass, and demonstrate that JWST can resolve the SMBH sphere of influence in the presence of significant AGN continuum. The work highlights the importance of ensemble modeling to account for systematics and confirms the reliability of JWST NIRSpec stellar kinematics for precise SMBH mass measurements.

Abstract

We present a new stellar dynamical measurement of the supermassive black hole (SMBH) mass in the nearby spiral galaxy NGC 4258, a critical benchmark for extragalactic mass measurements. We use archival JWST/NIRSpec IFU data (G235H/F170LP grating) to extract high-resolution two-dimensional stellar kinematics from the CO bandhead absorption features within the central $3'' \times 3''$. We extract the stellar kinematics after correcting for instrumental artifacts and separating the stellar light from the non-thermal AGN continuum. We employ Jeans Anisotropic Models (JAM) to fit the observed kinematics, exploring a grid of 12 models to systematically test the impact of different assumptions for the point-spread function, stellar mass-to-light ratio ($M/L$) profile, and orbital anisotropy. All 12 models provide broadly acceptable fits, albeit with minor differences. The ensemble median and 68% (1$σ$) bootstrap confidence intervals of our 12 models yield a black hole mass of $M_{\rm BH} = (4.08^{+0.19}_{-0.33}) \times 10^7$ M$_\odot$. This paper showcases the utility of using the full model ensemble to robustly account for systematic uncertainties, rather than relying on formal errors from a single preferred model, as has been common practice. Our result is just 5% larger than, and consistent with, the benchmark SMBH mass derived from water maser dynamics, validating the use of NIRSpec stellar kinematics for robust SMBH mass determination. Our analysis demonstrates JWST's capability to resolve the SMBH's sphere of influence and deliver precise dynamical masses, even in the presence of significant AGN continuum emission.

Figures (10)

• Figure 1: Modeling of the wiggles in single-spaxel spectra for NGC 4258. Top panel: The integrated spectrum (orange), single-spaxel spectrum (blue), and residual wiggles (gray). The red curve shows the best-fit model to the wiggles. Bottom panel: The corrected single-spaxel spectrum (dark blue), compared to the integrated spectrum (orange); the gray curve shows the residuals after correction. In all panels, red shaded regions indicate emission lines excluded from the fit.
• Figure 2: Left: Logarithmic integrated intensity map of the NIRSpec G235H/F170LP data cube, collapsed along the spectral axis, excluding the detector gap at $\sim$2.41--2.49 . The annular region outlined by the two red rings ($0\farcs4 < r < 1\farcs4$) indicates the area from which the global spectrum was extracted. Right: Optimal stellar template for NGC 4258. The observed global spectrum (black line) is compared to the broadened best-fit template obtained with pPXF (red line). The fit residuals (green points) are vertically offset by +0.62 to compress the $y$-axis range and better illustrate the stellar CO absorption bandheads. This same vertical offset is applied to all subsequent figures of this type.
• Figure 3: Radial variation in the NIRSpec G235H/F170LP spectrum of NGC 4258. Left column: Each panel shows the observed spectrum (black line), obtained by coadding spaxels within circular annuli of radius $r$. The best-fitting pPXF model (red line) includes the global stellar template (left panel of \ref{['fig:Intensity_maps']}, constrained from \ref{['sec:global_template']}) convolved with a Gaussian LOSVD, plus fourth-degree multiplicative and additive polynomials, which represent the nuclear non-thermal spectrum. Residuals are shown as green dots below each fit. Spectral regions affected by emission lines or artifacts were excluded from the fit (gray areas). Right column: The global stellar template convolved with the LOSVD (red line) is compared to the observed spectrum after subtraction of the nuclear non-thermal component (black line).
• Figure 4: The radial surface brightness profile $I(r)$, measured from individual Voronoi bins in the NIRSpec G235H/F170LP datacube of NGC 4258 (filled black dots), is compared to the estimated stellar light profile $\Gamma(r) = I(r)\gamma(r)$ (red open circles). The underlying stellar distribution is smooth and well approximated by a single-power-law (red line with pink region shows its $1\sigma$ uncertainty). The PSF–convolved central interpolation of this single power-law also reproduces the two innermost $\gamma$ measurements, which lie below the intrinsic (unconvolved) profile, indicating that their apparent decline is fully consistent with PSF effects. The gray points represent the pixels where $\gamma < 0.5$.
• Figure 5: Stellar kinematics extracted from the JWST/NIRSpec G235H F170LP IFU of NGC 4258 are presented. Panel A: Same as the left panel of \ref{['fig:Intensity_maps']}. Panel B displays a representative pPXF fit for a central-offset Voronoi bin that is unaffected by wiggle artifacts; its location is marked in Panel A. The observed spectrum, which includes the stellar CO band heads absorption features at $\sim$2.3 , is shown in black, with the best-fit empirical XSL template overlaid in red. The fit residuals, (data - model), are shown in green. The vertical gray lines indicate the spectral range used in fitting the templates to the spectra in all bins. Panels C–F: present maps of the stellar rotation velocity ($V$), velocity dispersion ($\sigma$), root-mean-square velocity ($V_{\rm rms} = \sqrt{V^2 + \sigma^2}$), and stellar light contribution fraction ($\gamma$). White contours trace the intensity distribution, decreasing by 1 magnitude per arcseconds$^{2}$ from the center outward.