Pith. Sign in

REVIEW

The independence of the mid-infrared RR Lyrae Period-Luminosity relation from metallicity from a study of three globular clusters in the Large Magellanic Cloud

Solar and Stellar Astrophysics · Submitted 18 Sep 2026 · Published version · 1 reference

Pith's verdict: Careful photometric work with useful new distance moduli, but the PLZ metallicity claim is weaker than it looks — the within-cluster scatter that constrains c is dominated by 0.28 dex measurement noise.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on September 21, 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.21494 (version 1)
Pith IDM44DKKV3, Pith's permanent record of this paper and its review
CitationSteve Ardern, Stijn Wuyts, Vicky Scowcroft. "The independence of the mid-infrared RR Lyrae Period-Luminosity relation from metallicity from a study of three globular clusters in the Large Magellanic Cloud." arXiv:2609.21494 (2026). https://pith.science/paper/M44DKKV3
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

Mid-infrared RR Lyrae period-luminosity relation shows no dependence on metallicity

The reading

This paper aims to establish that the mid-infrared ([3.6 μm]) period-luminosity (PL) relation for fundamental-mode RR Lyrae stars (RRab) is not affected by the star's metal content. By simultaneously fitting the PL relation for three old globular clusters in the Large Magellanic Cloud, the authors find a metallicity coefficient c = −0.033 ± 0.054 mag dex⁻¹, consistent with zero and implying |c| < 0.1 mag dex⁻¹. If this holds true, RR Lyrae distances no longer require accurate metallicity measurements, greatly simplifying the use of these stars as standard candles and strengthening a Cepheid-independent rung of the cosmic distance ladder.

What carries the argument

A Bayesian linear regression model that simultaneously fits the PLZ relation to multiple clusters, incorporating a cluster-level offset parameter and individual stellar metallicities derived from the Muraveva et al. (2025) empirical formula using Gaia DR3 G-band Fourier parameters (φ31) and periods. The fit is weighted by photometric uncertainties and optionally includes an intrinsic width term in the likelihood to capture residual scatter.

Extended reading notes

Core claim

Using Spitzer [3.6] photometry and Gaia DR3 data for RRab stars in three LMC globular clusters (Reticulum, NGC 1841, NGC 1466), the authors perform a simultaneous Bayesian linear regression of the period-luminosity-metallicity (PLZ) relation. They report a metallicity coefficient c = −0.033 ± 0.054 mag dex⁻¹, statistically indistinguishable from zero. They conclude that the mid-IR PL of RRab stars has a negligible dependence on metallicity, with an upper limit |c| < 0.1 mag dex⁻¹. The same analysis yields extinction-corrected distance moduli for the three clusters that agree with literature values and an intrinsic width of about 0.1 mag in the PL and PLZ relations, indicating that the residual scatter is driven by factors other than metallicity.

Load-bearing premise

Individual RRab metallicities derived from the Muraveva et al. (2025) empirical formula using Gaia DR3 Fourier parameters are accurate to within 0.28 dex, and the extinction estimates from Walker (1990, 1992a,b) are correct.

Editorial extensions

If this is right

  • RR Lyrae distances from mid-infrared PL relations will not require metallicity corrections, removing a potential systematic uncertainty.
  • A Cepheid-independent distance ladder using RR Lyrae becomes simpler and cheaper to calibrate because individual metallicity measurements are unnecessary.
  • Distance moduli to the three LMC clusters (first such measurements for NGC 1841 and NGC 1466) are now available from RR Lyrae PL relations, consistent with prior methods.
  • The detected intrinsic width of ≈0.1 mag suggests that non‑metallicity factors (e.g., effective temperature spread, evolutionary state) contribute to the PL scatter, motivating further theoretical work.
  • The methodology can be extended to more clusters and nearby galaxies to test the universality of the metallicity‑independent mid‑IR PL.

Reading between the lines

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

  • If the mid‑IR PL is indeed metallicity‑independent across a wider range of metallicity ([Fe/H] from −2.2 to −1.5 dex), then RR Lyrae in high‑redshift galaxies observed with JWST can be used as distance indicators without needing spectroscopic metallicity follow‑up.
  • This result strengthens the case for an RR Lyrae‑based distance ladder that bypasses Cepheid‑specific systematics, potentially informing the Hubble tension by providing an independent anchor.
  • The residual intrinsic width of 0.1 mag sets a floor on distance precision from single‑epoch mid‑IR observations; future surveys aiming for <1% distances will need to account for or model this scatter.
  • A direct spectroscopic test of the Muraveva et al. (2025) metallcity formula, using medium‑resolution spectra of RRab in these clusters, would either validate or challenge the foundation of the PLZ analysis.
Original abstract

RR Lyrae are pulsating variable stars tracing old ($>10$ Gyr) stellar populations, and exhibit a strong mid-infrared Period-Luminosity (PL) relation that can be calibrated with known parallaxes to infer distances. Here we present Period-Luminosity relations and distance moduli for three isolated globular clusters in the Large Magellanic Cloud: Reticulum, NGC 1841, and NGC 1466. Our analysis uses legacy \textit{Spitzer Space Telescope} images obtained by the Carnegie RR Lyrae Program, and an internally self-consistent sample restricted to RRab stars, with cluster membership confirmed using Gaia DR3 proper motions and photometry. In the Spitzer 3.6~$\mu$m band, we simultaneously fit a PL for these clusters using a slope derived from Galactic Globular Clusters, yielding extinction-corrected distance moduli of $18.47\pm0.09$ mag for Reticulum, $18.29\pm0.09$ mag for NGC 1841, and $18.61\pm0.09$ mag for NGC 1466. The latter two are the first RR Lyrae PL-based distance moduli for these clusters, and all three are consistent with literature values from other techniques. Additionally we fit a PL with slope as a free parameter, and find that this LMC-derived PL is consistent with that derived from Galactic GCs. Simultaneously fitting a PLZ for the three clusters yields a metallicity coefficient $c= -0.03\pm0.05~\mathrm{mag}~\mathrm{dex}^{-1}$ which can be considered a negligible dependence of the PL on metallicity ($|c| <0.1$ mag dex$^{-1}$). Modelling an intrinsic width, $W$, to the PL/ PLZ yields a consistent value $W\approx0.1\pm0.02~\mathrm{mag}$, suggesting intrinsic width is not driven by metallicity.

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,897 input tokens.

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

The abstract, sentence by sentence

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

RR Lyrae are pulsating variable stars tracing old ($>10$ Gyr) stellar populations, and exhibit a strong mid-infrared Period-Luminosity (PL) relation that can be calibrated with known parallaxes to infer distances. Here we present Period-Luminosity relations and distance moduli for three isolated globular clusters in the Large Magellanic Cloud: Reticulum, NGC 1841, and NGC 1466. Our analysis uses legacy \textit{Spitzer Space Telescope} images obtained by the Carnegie RR Lyrae Program, and an internally self-consistent sample restricted to RRab stars, with cluster membership confirmed using Gaia DR3 proper motions and photometry. In the Spitzer 3.6~$\mu$m band, we simultaneously fit a PL for these clusters using a slope derived from Galactic Globular Clusters, yielding extinction-corrected distance moduli of $18.47\pm0.09$ mag for Reticulum, $18.29\pm0.09$ mag for NGC 1841, and $18.61\pm0.09$ mag for NGC 1466. The latter two are the first RR Lyrae PL-based distance moduli for these clusters, and all three are consistent with literature values from other techniques. Additionally we fit a PL with slope as a free parameter, and find that this LMC-derived PL is consistent with that derived from Galactic GCs. Simultaneously fitting a PLZ for the three clusters yields a metallicity coefficient $c= -0.03\pm0.05~\mathrm{mag}~\mathrm{dex}^{-1}$ which can be considered a negligible dependence of the PL on metallicity ($|c| <0.1$ mag dex$^{-1}$). Modelling an intrinsic width, $W$, to the PL/ PLZ yields a consistent value $W\approx0.1\pm0.02~\mathrm{mag}$, suggesting intrinsic width is not driven by metallicity.

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 RR Lyrae are pulsating variable stars tracing old ($>10$ Gyr) stellar populations, and exhibit a strong mid-infrared Period-Luminosity (PL) relation that can be calibrated with known parallaxes to infer distances. background 1.00 0.06 0.04 asserted 1.00
s1 Here we present Period-Luminosity relations and distance moduli for three isolated globular clusters in the Large Magellanic Cloud: Reticulum, NGC 1841, and NGC 1466. contribution 0.99 0.95 0.05 asserted 0.99
s2 Our analysis uses legacy \textit{Spitzer Space Telescope} images obtained by the Carnegie RR Lyrae Program, and an internally self-consistent sample restricted to RRab stars, with cluster membership confirmed using Gaia DR3 proper motions and photometry. method 1.00 0.84 0.02 asserted 0.99
s3 In the Spitzer 3.6~$\mu$m band, we simultaneously fit a PL for these clusters using a slope derived from Galactic Globular Clusters, yielding extinction-corrected distance moduli of $18.47\pm0.09$ mag for Reticulum, $18.29\pm0.09$ mag for NGC 1841, and $18.61\pm0.09$ mag for NGC 1466. result 0.95 0.97 0.07 asserted 1.00
s4 The latter two are the first RR Lyrae PL-based distance moduli for these clusters, and all three are consistent with literature values from other techniques. result 0.67 0.77 0.07 asserted 0.69
s5 Additionally we fit a PL with slope as a free parameter, and find that this LMC-derived PL is consistent with that derived from Galactic GCs. result 0.93 0.97 0.12 asserted 0.62
s6 Simultaneously fitting a PLZ for the three clusters yields a metallicity coefficient $c= -0.03\pm0.05~\mathrm{mag}~\mathrm{dex}^{-1}$ which can be considered a negligible dependence of the PL on metallicity ($|c| <0.1$ mag dex$^{-1}$). result 0.92 0.96 0.80 asserted 0.75
s7 Modelling an intrinsic width, $W$, to the PL/ PLZ yields a consistent value $W\approx0.1\pm0.02~\mathrm{mag}$, suggesting intrinsic width is not driven by metallicity. interpretation 0.52 0.95 0.11 hedged possibility 0.95

What kind of paper the classifier reports this is: empirical study 1.00 · 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. Pulsation Hertzsprung-Russell Diagram (adapted from Jeffery & Saio 2016) with permission. The solid red line shows the Zero Age Main Sequence (ZAMS), while selected evolutionary tracks are shown as dotted lines, with red numbers showing masses (M⊙) for those tracks. The classical instability strip is shown by the dashed black lines. RRL are shown by a yellow-filled region, and are found where the Horizontal Branch crosses the instability strip – the Zero Age Horizontal Branch (ZAHB) is shown as a red dot-dashed line. Cepheids are included as the yellow diagonally-hashed region for…
Figure 2. Sky positions of the 14 old LMC clusters that are identified by Sarajedini (2024) as known hosts of RR Lyrae. Numbers indicate NGC object identifiers. The three clusters used in this work are shown in teal, and the remaining clusters are shown in orange. The position of the central bar of the LMC is shown by the dashed oval (Rathore et al. 2025). The three clusters used in this work are in the LMC’s halo which has a low LMC field star density. The majority of the remaining clusters are located in, or close to, the central bar with its significantly higher LMC field-star density
Figure 3. Sky position, proper motion and colour-magnitude diagrams for the three clusters. Sky positions (left) are plotted in gnomonic projection (Δ𝛼 cos 𝛿 versus Δ𝛿) to give undistorted representation of positions on the sky, grey dots represent all Gaia objects to 6 arcmin, orange triangles rep￾resent Gaia RRab to 3 arcmin of cluster centre. Proper motions (centre) are shown in mas yr−1 , grey dots represent all Gaia objects to 6 arcmin, orange triangles represent Gaia RRab to 3 arcmin of cluster centre. The black dashed circle shows the region of proper-motion space used for cut.…
Figure 4. Pixel-phase correction sampling for [3.6] photometry, showing histograms of 100,000 simulated single images and 100,000 normalised 7- image dithers. The individual images have a mean dimensionless pixel-phase photometric correction of 0.9707 ± 0.0228, while for the 7-image dithers the value is 0.9707 ± 0.0086. errors is ±0.025 mag for single images, for the 7-image dither this was ±0.009 mag. This pixel-phase photometric error for the 7-image dither can be considered negligible in comparison to dominant errors and is not included in our error analysis. 3.4 Light curves and RRab mean…
Figure 5. Example gloess fitted phase-folded light curves for RRab across a range of periods for the three clusters – top row Reticulum, middle row NGC 1841, bottom row NGC 1466. The full set of recovered light curves is provided in Appendix A. The orange data points and errorbars are our photometric results for each epoch, the black light curve is the fit provided by GLOESS to this dataset. A common magnitude range is used on all light curves. Note that some photometric data points do not have errorbars as the EPSF fitter was unable to provide them, however we have included these points for…
Figure 6. Corner plot for the simultaneous three cluster fit PL with fixed slope showing representations of the posterior distributions of the parameters 𝑏 (intercept), Δ𝑏1 (offset for NGC 1841), Δ𝑏2 (offset for NGC 1466) and their cross-correlations. (2026), 𝜇0 = 18.472 ± 0.035 mag, is in excellent agreement with our value for Reticulum of 18.47±0.09 mag, with our slightly larger uncertainty being dominated by the uncertainty in the zero point of Neeley et al. (2017) of ±0.089 mag. Figures 8b and 8c make the same comparison for NGC 1841 and NGC 1466, and while there are fewer historic…
Figure 7. Period-Luminosity relation for RRab in the LMC globular clusters studied in this work. Apparent magnitudes were offset by the clusters’ extinction￾corrected distance moduli relative to Reticulum (Δ𝜇) to allow comparison as a single population: Δ𝜇NGC 1841 = +0.180 mag, and Δ𝜇NGC 1466 = −0.135 mag. It should be noted that there is a difference in the dynamic range of periods between the three clusters which may be reflective of different cluster ages and RGB mass-loss histories
Figure 8. Comparison of distance moduli obtained by a range of methods for Reticulum, NGC 1841, and NGC 1466. Orange data points are from colour￾magnitude diagram (CMD) fitting, black data points are from PL relations, and blue points are from Horizontal Branch/ RRL mags. Open circles denote data points with no uncertainty reported. Reticulum Cluster studies: Bhardwaj et al. (2026), Muraveva et al. (2018) (‘HST’ and ‘Gaia’ labels denote how the zero-point was determined), Wagner-Kaiser et al. (2017), Pieres et al. (2016), Jeon et al. (2014), Kuehn et al. (2013), Sollima et al. (2008),…
Figure 9. Corner plot for the simultaneous three cluster fit PL with slope as a free parameter. 𝑎 is the slope of the PL, 𝑏 is intercept of the PL, while Δ𝑏1 and Δ𝑏2 are the offsets from Reticulum for NGC 1841 and NGC 1466 respectively. percentiles of the distribution. Our finding of negligible dependence of the PL on metallicity (|𝑐| < 0.1 mag dex−1 ) is in line with several empirical studies cited earlier (e.g Sollima et al. 2008; Borissova et al. 2009; Muraveva et al. 2015). 6 THE INTRINSIC WIDTH OF THE RRL PERIOD-LUMINOSITY RELATION The analysis presented earlier in this section follows
Figure 10. Metallicities for the three clusters found using the empirical formula presented in M25. Values are [Fe/H] (dex) on the Crestani et al. (2021) scale. All RRab within 3 arcmin of the cluster centres with 𝜙31 parameter reported in the Gaia DR3 gaiadr3.vari-rrlyrae table are shown. Blue filled circles indicate RRab for which we recovered mean magnitudes and which were used in the PL, other RRab are shown by orange open circles. The dashed line shows mean value for all data points, and grey band indicates standard deviation as reported in
Figure 11. Corner plot for the PLZ showing representations of the posterior distributions of the parameters 𝑎 (slope), 𝑏 (intercept), 𝑐 (metallicity coef￾ficient), Δ𝑏1 (offset for NGC 1841), Δ𝑏2 (offset for NGC 1466), and their cross-correlations. prior to constrain 𝑊 to non-negative values (0 < 𝑊 < 1). Corner plots for the fits with intrinsic width are shown in

Showing 11 of the paper's 14 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:2609.21494 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.