The Distance to NGC 4258 from Individual Maser Component Tracking

TL;DR

This work reexamines the geometric distance to NGC 4258 by moving beyond the traditional averaged maser datasets and a single-orbit approximation. It tracks individual maser components across epochs, incorporates intrinsic linewidths, and marginalizes nuisance coordinates $r$ and $\phi$ to reduce parameter dimensionality, using a dynesty-based framework. The analysis shows that the current observational cadence is insufficient to reliably track maser components, with distance inferences highly sensitive to maser selection (yielding $D\approx 6.7$–$8.1$ Mpc across configurations). Reproductions of prior results on the averaged dataset reveal biases from the single-orbit assumption and line-width neglect, and the method yields a robust path toward percent-level accuracy only with high-cadence monitoring (e.g., $\sim$40 epochs at 10-day cadence for $\sim$10 tracked masers). Overall, NGC 4258 does not yet provide a percent-level geometric anchor, but the proposed cadence and component-tracking approach delineate a clear route to achieving that precision and strengthening its role in the extragalactic distance scale.

Abstract

We present a reanalysis of the water maser system in NGC 4258 to reassess its geometric distance, commonly reported as approximately 7.6 Mpc with percent-level accuracy, a key anchor in extragalactic distance ladder calibrations and recent determinations of the Hubble constant. We introduce a method that relies exclusively on tracking individual maser components, rather than assuming a single averaged trajectory as in previous works, thereby avoiding arbitrary data averaging that can bias interpretations of the disk's geometry and dynamics. This approach requires spatially resolved individual maser components; consequently, the majority of observational epochs were excluded, as they lack sufficient spatial resolution to localize the maser position. We track individual maser components across multiple epochs and introduce an efficient marginalization method over nuisance parameters (angular radius and azimuth of each maser), reducing the number of free parameters from hundreds to 14. Our analysis reveals that the current observational cadence is insufficient to reliably track the individual masers, which our method relies on, given their intrinsic variability. Across a range of maser selections and model configurations, inferred distances span approximately between 6.7 to 8.1 Mpc, demonstrating significant sensitivity to how our method selects individual masers. Even visually and statistically robust fits can differ by several standard deviations, reflecting ambiguity in component identification across sparsely sampled epochs. We evaluate the impact of observational cadence on tracking fidelity and distance precision, and show that high-cadence monitoring is needed for our method to track individual masers and produce a robust anchor for cosmology.

Figures (22)

• Figure 1: Comparison between the Argon_2007 dataset (grey) and the averaged dataset from Humphreys_2013 (pink), showing the east–west offset, $X$, as a function of LOS velocity, $v_{\text{los}}$. Although the averaged dataset includes uncertainties---typically around $1\,\text{$\mu$as}$---the error bars are barely visible, even after incorporating the $1.6\,\text{$\mu$as}$ error floors adopted by Reid_2019. Green lines show synthetic trajectories of two maser components, computed using initial angular radii of $3.5$ and $4.0\,\text{mas}$ (within the efficient amplification region), based on the disk parameters from Reid_2019. These trajectories demonstrate that multiple $X$ values can correspond to the same $v_{\text{los}}$. The inset in the upper left zooms in on the region near $v_{\text{los}} \approx 425\,\textrm{km}\,\textrm{s}^{-1}$, where two distinct $X$ values are observed in the Argon_2007 data. This example highlights a case in which the averaged dataset fails to capture the true spread of the measurements, and the assumption of a Gaussian distribution is clearly invalid.
• Figure 2: Distribution of normalized residuals $\hat{z}_i$ (defined in Equation \ref{['eq:z_i']}) for the blueshifted (left, blue), systemic (middle, green), and redshifted (right, red) masers, comparing $X$ positions from the averaged and Argon_2007 datasets. The dashed curves represent the standard normal distribution, $\mathcal{N}(0, 1)$, expected under the Gaussian assumptions of the averaging procedure. In all three cases, the observed residuals exhibit significant deviations from this model, indicating that the averaged dataset does not accurately capture the true distribution of positional measurements.
• Figure 3: Integrated spectra from Argon_2007 for epoch F (systemic, top), epoch J (redshifted, middle), and epoch G (blueshifted, bottom). The plotted quantity is normalized flux density (i.e., flux density divided by the maximum value in each spectrum) as a function of relativistic LSR velocity. Spectral peaks and their corresponding FWHMs are identified using the procedure described in Section \ref{['subsec:thermal_width']}. Gray circle markers ('o') denote the initial set of raw peaks, while red cross markers ('×') indicate the final set of clean, well-isolated peaks. FWHMs are shown as continuous gray lines for the raw peaks and dashed red lines for the final set. Prominence values exceeding 0.02 are labeled above the corresponding peaks in black. To standardize prominence calculations across epochs and velocity ranges, synthetic zero-flux values were added at the baseband edges, ensuring all peaks are evaluated relative to a common zero baseline.
• Figure 4: Results of the RANSAC algorithm, as described in Section \ref{['subsec:new_dataset']}, for the dataset corresponding to entry #8 in Table \ref{['tab:results_newdata']}. The top panel shows LOS velocity as a function of time; the middle and bottom panels display east–west ($X$) and north–south ($Y$) sky-plane positions, respectively, also plotted against time. Red crosses indicate observations rejected by the algorithm. Colored lines and marker edges correspond to individual maser components identified by RANSAC, each assigned a group ID consistent with the identifiers in the public catalogfn:url. Marker face color encodes $\log_{10}$(prominence).
• Figure 5: Posterior distributions (blue shaded regions) of the distance-fitting results from Table \ref{['tab:results_newdata']}, excluding statistically less likely fits. The red curve shows the distance posterior from Reid_2019 for comparison.