Renzo's rule revisited: A statistical study of galaxies' baryon - dark matter coupling

We present a systematic statistical analysis of an informal astrophysical phenomenon known as Renzo’s rule (or Sancisi’s law), which states that “for any feature in a galaxy’s luminosity profile, there is a corresponding feature in the rotation curve, and vice versa.” This is often posed as a challenge for the standard LCDM model while supporting alternative theories such as MOND. Indeed, we identify clear features in the dwarf spiral NGC 1560 – a prime example for Renzo’s rule – and find correlation statistics which support Renzo’s rule with a slight preference for MOND over LCDM halo fits. However, a broader analysis on galaxies in the SPARC database reveals an excess of features in rotation curves that lack clear baryonic counterparts, with correlation statistics deviating up to $3σ$ on average from that predicted by both MOND and LCDM haloes, challenging the validity of Renzo’s rule. Thus we do not find clear evidence for Renzo’s rule in present galaxy data overall. We additionally perform mock tests, which show that a definitive test of Renzo’s rule is primarily limited by the lack of clearly resolved baryonic features in current galaxy data.

💡 Research Summary

This paper presents the first systematic statistical test of the informal astrophysical observation known as Renzo’s rule (or Sancisi’s law), which claims that any feature in a galaxy’s luminosity profile is mirrored by a corresponding feature in its rotation curve, and vice‑versa. The authors examine whether this rule holds across a range of galaxies and whether it can be used to discriminate between the standard ΛCDM cosmology and Modified Newtonian Dynamics (MOND).

The study uses four data sets: (i) two independent rotation‑curve reductions for the dwarf spiral NGC 1560, (ii) the SPARC catalog of 175 late‑type galaxies (restricted to 60 with ≥20 data points per component), (iii) a subset of the LITTLE THINGS dwarf galaxy sample (18 objects), and (iv) 21 simulated galaxies from the MaGICC and CLUES hydrodynamical ΛCDM simulations. For each galaxy the observed total rotation speed V_obs, the contributions from gas, stellar disc and bulge, and the combined baryonic speed V_bar are assembled.

Two theoretical predictions are generated for each system. In the MOND case the authors apply Milgrom’s formula with the “simple” interpolating function ν(y)=1+√(1+4y)−½, using the observed baryonic distribution to compute V_MOND. For ΛCDM they adopt two approaches: (a) an MCMC fit of V_obs that treats the halo parameters (e.g., NFW concentration and scale radius) and stellar mass‑to‑light ratios as free parameters, and (b) a direct Monte‑Carlo propagation of the prior distributions to produce V_ΛCDM without fitting the data.

To isolate small‑scale structure, a Gaussian Process Regression is applied to each rotation curve, removing the smooth large‑scale trend and leaving residuals that highlight local “wiggles”. An automated peak‑and‑valley detection algorithm identifies features in the residuals of both V_obs and V_bar. The authors then compute Pearson correlation coefficients and Dynamic Time Warping (DTW) alignment costs between the two residual series as quantitative measures of the correspondence predicted by Renzo’s rule.

Results for NGC 1560 show a clear “kink” between 4–6 kpc that appears in both the gas component and the total rotation curve, yielding a higher Pearson r and lower DTW cost for the MOND prediction than for ΛCDM, thus providing a single‑galaxy example that supports the rule. However, when the same analysis is applied to the full SPARC sample, the average correlation drops dramatically (r≈0.3) and the DTW costs deviate by up to 3σ from the expectations of both MOND and ΛCDM. In many cases rotation‑curve wiggles have no obvious counterpart in the baryonic profile, indicating a systematic breakdown of Renzo’s rule across the population.

The simulated galaxies show mixed behavior: a few MaGICC and CLUES objects reproduce the expected correspondence, but the majority display a wide scatter similar to the observational sample, suggesting that current ΛCDM simulations do not universally generate the fine‑grained baryon–dark‑matter coupling implied by Renzo’s rule.

Mock experiments reveal that the principal limitation is data quality. Reducing observational uncertainties by half, tightening stellar mass‑to‑light ratio errors to ≲0.1 dex, and achieving sub‑kiloparsec spatial resolution would allow a 5σ detection of the rule, if it truly exists. Presently, H I resolution, distance and inclination uncertainties, and the lack of sharply defined baryonic features dominate the error budget.

In conclusion, Renzo’s rule is demonstrably present in specific cases such as NGC 1560 but does not hold as a universal statistical property of late‑type galaxies. The modest preference for MOND in the single‑galaxy test does not translate into a decisive discrimination between MOND and ΛCDM on a population level. Future progress will require higher‑resolution kinematic maps, more precise stellar population modeling, and improved ΛCDM hydrodynamical simulations that can capture subtle baryon–dark‑matter interactions.