Pith. Sign in

REVIEW

A reanalysis of the megamaser Hubble constant: from spot catalogues to peculiar velocities

Cosmology and Nongalactic Astrophysics · Submitted 15 Sep 2026 · 6 references

Pith's verdict: A careful, transparent reanalysis of the megamaser H0 that advances the methodology but does not resolve the tension; deserves referee time.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on September 17, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 2 major objections · 6 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2609.17684 (version 1)
Pith ID742577BO, Pith's permanent record of this paper and its review
CitationHarry Desmond, Richard Stiskalek. "A reanalysis of the megamaser Hubble constant: from spot catalogues to peculiar velocities." arXiv:2609.17684 (2026). https://pith.science/paper/742577BO
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

A Bayesian reanalysis of five megamaser galaxies, incorporating warped disc models, sample selection, and peculiar-velocity corrections, yields a geometric Hubble constant of 71.6 ± 3.7 km/s/Mpc (statistical), intermediate between the CMB…

The reading

This paper reanalyzes the five megamaser galaxies of the Megamaser Cosmology Project using a Bayesian forward model that jointly infers each galaxy's angular-diameter distance and warped Keplerian disc parameters directly from VLBI observables. It then combines these distances with host redshifts, including a physical volume prior, a sample selection correction, and a non-linear peculiar-velocity reconstruction (Manticore-Local), to obtain H0 = 71.6 ± 3.7 km/s/Mpc (statistical). This value sits between the Planck CMB measurement (67.4) and the SH0ES distance-ladder measurement (73.2), and its uncertainty does not yet distinguish between the two. The analysis shows that both the peculiar-velocity treatment and the choice of distance prior appreciably shift the result, while the updated maser spot tables from the MCP move it by about 2 km/s/Mpc relative to the original published catalogues.

What carries the argument

The central machinery is a warped Keplerian accretion disc model that maps each maser spot's orbital coordinates (radius and azimuth) to its observed sky position, LOS velocity, and LOS acceleration, with the angular-diameter distance D_A and black hole mass M_BH entering through the Keplerian velocity and centripetal acceleration equations. The authors break the mass–distance degeneracy by sampling the combination η = log[(M_BH/M_⊙)/(D_A/Mpc)] and then infer D_A jointly with disc geometry parameters using a two-step procedure: a differential-evolution global search locates the posterior maximum, after which a Metropolis-within-Gibbs sampler (NUTS for global parameters, adaptive random-walk for per-spot coordinates) explores the full posterior, including a dedicated reflection move to handle the bimodal azimuthal angles.

Extended reading notes

Core claim

From five megamaser-hosting galaxies, the authors obtain a purely geometric, ladder-independent constraint on the Hubble constant of H0 = 71.6 ± 3.7 km/s/Mpc (statistical) using a fiducial configuration that includes a redshift-based selection threshold, the non-linear Manticore-Local peculiar-velocity reconstruction, and a uniform-in-volume distance prior. This result resolves none of the Hubble tension but places the megamaser measurement in the intermediate regime, inconsistent with neither Planck nor SH0ES at the current 5% precision. The paper further demonstrates that the MCP's internal reanalysis of maser spot measurements (affecting mainly NGC 6264 and NGC 6323) raises H0 by 1.8–2.1 km/s/Mpc relative to the published catalogues, and that the non-linear peculiar-velocity correction lowers it by 1.6–2.1 km/s/Mpc relative to ignoring peculiar motions.

Load-bearing premise

The non-linear peculiar-velocity reconstruction (Manticore-Local) and its residual bulk-flow prior, calibrated from Cosmicflows-4, provide unbiased predictions for the line-of-sight velocities of the five maser hosts when evaluated at the sampled comoving distances.

Editorial extensions

If this is right

  • If this result is correct, the megamaser H0 no longer favours the local distance ladder over the CMB, as previous MCP analyses suggested; it is consistent with both at the ~1σ level.
  • The sample selection correction, though crude, matters: omitting it shifts H0 downward by ~2 km/s/Mpc, so future larger maser samples will require well-characterized selection functions.
  • The non-linear peculiar-velocity reconstruction from Manticore-Local gives a larger H0 reduction than the linear Carrick+15 field, and the difference is driven by the residual bulk-flow prior; proper marginalization over the local velocity field is essential.
  • The MCP's updated spot tables change the inferred distances significantly for two galaxies (NGC 6264 and NGC 6323), indicating that re-reduction of the raw VLBI data can shift H0 by ~2 km/s/Mpc, comparable to the current statistical uncertainty.
  • The separate analysis of NGC 4258 with an eccentric quadratic warp yields D_A = 7.558 ± 0.073 Mpc, consistent with the published trigonometric distance to <0.2σ, validating the disc model for the nearest maser galaxy.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If future megamaser samples reach 1% precision, the intermediate H0 value found here could either narrow the Hubble tension if it stays near 71–72 or sharpen it if it moves closer to SH0ES or Planck. The present uncertainty is too large to discriminate.
  • The paper's sensitivity to the choice of distance prior (flat vs. log-flat vs. volume-squared) suggests that any two-stage analysis that first derives distances under a non-physical prior and then reweights may be suboptimal; a full joint inference of all distances and H0 within a single hierarchical model would eliminate this reweighting ambiguity.
  • The reliance on a single, pre-computed peculiar-velocity reconstruction (Manticore-Local) calibrated with a Cosmicflows-4 prior means the H0 error bar does not include uncertainty from the reconstruction method itself. Comparing multiple independent reconstructions would provide a more robust estimate.
  • Because the MCP spot catalogue updates were only revealed through an erratum, this analysis highlights the importance of transparent data releases in high-stakes cosmology: the H0 shift from updated measurements exceeds the formal uncertainty of some earlier analyses.
Original abstract

The Hubble tension motivates careful scrutiny of the redshift-independent distances underpinning the local distance ladder. Water megamasers afford a purely geometric distance from orbital dynamics in an accretion disc. For five of the six Megamaser Cosmology Project (MCP) galaxies, we present a Bayesian forward model that infers the angular-diameter distance and warped Keplerian disc parameters jointly from the very long baseline interferometry quasi-observables: maser spot positions, velocities, and accelerations. We then constrain $H_0$ from these distances and the host redshifts, including a physical distance prior, modelling the galaxy sample selection and testing a linear and a non-linear peculiar-velocity reconstruction. Accounting for peculiar velocities with the non-linear Manticore reconstruction, we infer $H_0 = 71.6 \pm 3.7$ km/s/Mpc, where the uncertainty is statistical. Relative to ignoring peculiar velocities, the linear and non-linear reconstructions lower $H_0$ by $0.8$--$1.4$ and $1.6$--$2.1$ km/s/Mpc, respectively, whereas changing the selection variable alters it by less than $1$ km/s/Mpc. We analyse the sixth galaxy, NGC 4258, separately. Modelling its disc as eccentric and quadratically warped, we infer an angular-diameter distance of $7.558 \pm 0.073$ (stat.) Mpc, within $0.13\sigma$ of the published maser distance.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 6,562 input tokens.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.49 · human participants 0.02 · controlled comparison 0.70 · randomized assignment 0.01 · states quantitative result 0.99 · compares with prior results 0.95 · mentions code or data availability 0.03.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

The Hubble tension motivates careful scrutiny of the redshift-independent distances underpinning the local distance ladder. Water megamasers afford a purely geometric distance from orbital dynamics in an accretion disc. For five of the six Megamaser Cosmology Project (MCP) galaxies, we present a Bayesian forward model that infers the angular-diameter distance and warped Keplerian disc parameters jointly from the very long baseline interferometry quasi-observables: maser spot positions, velocities, and accelerations. We then constrain $H_0$ from these distances and the host redshifts, including a physical distance prior, modelling the galaxy sample selection and testing a linear and a non-linear peculiar-velocity reconstruction. Accounting for peculiar velocities with the non-linear Manticore reconstruction, we infer $H_0 = 71.6 \pm 3.7$ km/s/Mpc, where the uncertainty is statistical. Relative to ignoring peculiar velocities, the linear and non-linear reconstructions lower $H_0$ by $0.8$--$1.4$ and $1.6$--$2.1$ km/s/Mpc, respectively, whereas changing the selection variable alters it by less than $1$ km/s/Mpc. We analyse the sixth galaxy, NGC 4258, separately. Modelling its disc as eccentric and quadratically warped, we infer an angular-diameter distance of $7.558 \pm 0.073$ (stat.) Mpc, within $0.13\sigma$ of the published maser distance.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 The Hubble tension motivates careful scrutiny of the redshift-independent distances underpinning the local distance ladder. background 1.00 0.06 0.03 asserted 0.92
s1 Water megamasers afford a purely geometric distance from orbital dynamics in an accretion disc. background 0.89 0.06 0.02 asserted 0.98
s2 For five of the six Megamaser Cosmology Project (MCP) galaxies, we present a Bayesian forward model that infers the angular-diameter distance and warped Keplerian disc parameters jointly from the very long baseline interferometry quasi-observables: maser spot positions, velocities, and accelerations. method 0.83 0.96 0.04 asserted 1.00
s3 We then constrain $H_0$ from these distances and the host redshifts, including a physical distance prior, modelling the galaxy sample selection and testing a linear and a non-linear peculiar-velocity reconstruction. method 0.96 0.95 0.05 asserted 0.98
s4 Accounting for peculiar velocities with the non-linear Manticore reconstruction, we infer $H_0 = 71.6 \pm 3.7$ km/s/Mpc, where the uncertainty is statistical. result 1.00 0.97 0.79 asserted 1.00
s5 Relative to ignoring peculiar velocities, the linear and non-linear reconstructions lower $H_0$ by $0.8$--$1.4$ and $1.6$--$2.1$ km/s/Mpc, respectively, whereas changing the selection variable alters it by less than $1$ km/s/Mpc. result 0.96 0.95 0.13 asserted 1.00
s6 We analyse the sixth galaxy, NGC 4258, separately. contribution 0.82 0.88 0.02 asserted 0.81
s7 Modelling its disc as eccentric and quadratically warped, we infer an angular-diameter distance of $7.558 \pm 0.073$ (stat.) Mpc, within $0.13\sigma$ of the published maser distance. result 0.99 0.96 0.05 asserted 1.00

What kind of paper the classifier reports this is: empirical study 0.78 · evidence from analysis of existing data 0.99.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. — Position–velocity diagrams of the maser spots in NGC 5765b, CGCG 074-064, and NGC 6264, coloured by spectral class, for the updated-MCP tables (Section 2.1). The vertical axis is the LOS velocity offset from the systemic velocity. For illustration, the horizontal axis is the maser sky position projected onto the disc major axis, which we take empirically as the direction of maximum variance of the spot positions. The systemic spots cluster near the origin, in front of the black hole, while the approaching and receding high-velocity spots populate the disc tangents at opposite ends…
Figure 2. — Illustration of the per-spot azimuthal bimodality in ϕ for NGC 5765b. Each panel shows the data likelihood L(di | r, ϕ, θ) of a single maser spot as a function of its orbital angular radius r and azimuthal angle ϕ, with the disc parameters θ fixed at the baseline linear-warp posterior median and the colour scale normalised to each spot’s peak. The likelihood splits into the two solutions ϕ and π − ϕ, symmetric about the tangent point ϕ = π/2 for receding spots (panels a and b) and ϕ = 3π/2 for approaching spots (panel c). In every panel the two solutions lie on the near and far…
Figure 3. — Observed CMB-frame redshift zobs for each galaxy against its inferred comoving distance DC on the updated-ours spot table (Section 2.1), with the horizontal error bars giving the central 68% interval about the posterior median. The vertical error bars show a 250 km s−1 redshift uncertainty, representing the host peculiar velocity. The dashed and dotted lines are the flat-ΛCDM distance–redshift relations for the SH0ES distance-ladder (H0 = 73.2 km s−1 Mpc−1 ) and CMB (H0 = 67.4 km s−1 Mpc−1 ) values. since omitting them biases the result (Desmond et al. 2026). The per-galaxy…
Figure 4. — Angular-diameter distance posterior of each galaxy under the four spot tables of Section 2.1—original-all (blue), updated-MCP (red), updated-all (brown), and updated-ours (yellow)—together with the P20 posterior (black, dashed). The priors, sampler, and disc model are held fixed across the four inferences, except that the updated-MCP model has no NGC 5765b clump 2 floors (Section 3.1). The updated-MCP and updated-ours tables are based on the current MCP measurements and differ in the outlier treatment. criterion is applied, or, for the NGC 5765b clump 2 spots, on whether they are…
Figure 5. — Inferred H0 (median and central 68% interval) over the reconstruction × selection grid, for the linear warp and the volume prior. Markers denote the baseline input spot table: filled circles the updated-MCP table, filled squares the updated-ours table (Section 2.1). Bands give the CMB (Planck Collaboration and others 2020a) and SH0ES (Breuval et al. 2024) values for reference. values reproduce their published results. Passing the same distance posteriors through our own pop￾ulation model instead (the uniform-in-volume prior, sampled σv, the selection term, and the five-galaxy…
Figure 6. — Posterior on the global disc parameters of NGC 5765b for the baseline linear-warp model on the updated-ours table, typical of the five galaxies. Contours are the 68% and 95% credible regions. The four σ (2) parameters are the clump 2 error floors, which the model contains only on the tables the MCP cut is not applied to. All global disc parameters have R <ˆ 1.01, with ≈6.2 × 103 effective samples on DA
Figure 7. — Posterior on the global disc parameters of NGC 4258 under the eccentric quadratic-warp model, inferred from the full 358-spot table with no outlier rejection. Contours are the 68% and 95% credible regions. The eccentricity enters through the Cartesian components ex and ey, with e = (e 2 x + e 2 y ) 1/2 = 0.0041 ± 0.0009, and dω/dr is the periapsis warp rate
Figure 8. — Joint posterior over the shared population parameters for the redshift-selection runs on the updated-ours baseline spot table (Section 2.1), overlaying the Manticore-Local (blue) and linear (red) reconstructions, both with a linear warp. Contours are the 68% and 95% credible regions. The two reconstructions give consistent H0, while σpec and Vext recover their priors

Showing 8 of the paper's 9 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:2609.17684 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.