arxiv: 2605.00344 · v1 · submitted 2026-05-01 · 🌌 astro-ph.SR

Recognition: unknown

Exploring RR Lyrae Variable Stars in the Vera C. Rubin Observatory Data Preview 1

Chow-Choong Ngeow , Anupam Bhardwaj

Authors on Pith no claims yet

Pith reviewed 2026-05-09 19:27 UTC · model grok-4.3

classification 🌌 astro-ph.SR

keywords RR LyraeRubin ObservatoryDP1variable starsperiod-luminosity-metallicitydistancesstellar pulsation modelsphotometric metallicities

0 comments

The pith

RR Lyrae distances from Rubin DP1 data agree with literature using PWZ relations but are systematically larger with PLZ relations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper tests theoretical stellar pulsation models for LSST filters by analyzing about 600 known RR Lyrae stars in the Vera C. Rubin Observatory Data Preview 1 fields. It derives photometric metallicities from a gri-based metallicity-color relation and finds good agreement with literature values when light curve sampling allows accurate template fits. Distances are computed via period-luminosity-metallicity and period-Wesenheit-metallicity relations, revealing that the PWZ version matches known values closely while the PLZ version overestimates them, likely from model uncertainties tied to evolved horizontal branch stars. The work also shows that predicted period-amplitude relations from those models disagree with DP1 observations and concludes that sparse sampling makes current estimates sensitive to fitting choices.

Core claim

Distances to RR Lyrae stars in the DP1 fields based on the W_gr period-Wesenheit-metallicity relation match literature values with a mean offset of 0.01 plus or minus 0.36 mag, whereas PLZ-based distance moduli are systematically larger, possibly because of theoretical calibration uncertainties that include evolved horizontal branch models. Predicted period-amplitude relations from the evolved models are also inconsistent with amplitudes measured from the DP1 light curves.

What carries the argument

Theoretical period-luminosity-metallicity (PLZ) and period-Wesenheit-metallicity (PWZ) relations from stellar pulsation models tailored to LSST gri filters, applied after template fitting to derive metallicities and distances.

Load-bearing premise

The theoretical PLZ and PWZ relations calibrated on evolved horizontal branch models accurately represent RR Lyrae stars in the LSST gri filters and that template fitting yields reliable parameters even with sparse DP1 light curve sampling.

What would settle it

Future LSST data releases providing densely sampled light curves for the same stars would show whether PLZ distances continue to overestimate literature values or converge once sampling improves.

Figures

Figures reproduced from arXiv: 2605.00344 by Anupam Bhardwaj, Chow-Choong Ngeow.

**Figure 1.** Figure 1: Examples of two RR Lyrae, one from the 47Tuc and Fornax fields, respectively, with light curves exhibiting a large scatter at a given epoch. Both of these RR Lyrae only have two or four epochs of data available in the DP1. Upper and lower panels present multi-band light curves in time and in pulsation phases (after folded with published period), respectively. DiaObject catalog, two types of nested light cu… view at source ↗

**Figure 2.** Figure 2: The multiband light curves for the three known RR Lyrae in the EDFS field and the Rubin95 field. The best-fit template light curves to the objectForcedSource and the diaObjectForcedSource data are shown as black solid curves and cyan dashed curves, respectively. We did not fit the template light curves to the zy-band data for the RRc star in the Rubin95 field (the middle panels) because Braga et al. (2024)… view at source ↗

**Figure 3.** Figure 3: Comparison of the best-fit gri-band template light curves for CSS J023424.4+072337 in the Rubin38 field. The left and middle panels show the objectForcedSource and the diaObjectForcedSource light curves, respectively, together with the best fit template light curves. In the right panels, both sets of best-fit template light curves were plotted together, along with the improved i-band template light curves … view at source ↗

**Figure 4.** Figure 4: Comparison of the DP1 r-band and the Gaia DR3 G-band light curves for the only known RR Lyrae, SSS J035520.3-484728, in the EDFS field. 0.0 0.1 0.2 0.3 0.4 (g r) [mag] 0.10 0.05 0.00 0.05 0.10 0.15 (r i) [m a g] 3.59 ± 0.35 +1.41 ± 0.34 Z=0.02 Z=0.008 Z=0.004 Z=0.001 Z=0.0006 Z=0.0003 Z=0.0001 EDFS field Rubin95 field Rubin38 field view at source ↗

**Figure 5.** Figure 5: Comparison of the extinction-corrected colors for the four known RR Lyrae listed in view at source ↗

**Figure 6.** Figure 6: Comparison of the light curve amplitudes listed in view at source ↗

**Figure 7.** Figure 7: Left Panel: Comparison of the r-band amplitudes for the known RR Lyrae in the 47Tuc field (star symbols) to the theoretical Bailey diagrams based on the ZAHB models (Marconi et al. 2022). The r-band light curve for the outlier marked in the plot can be found in the upper-left panels of view at source ↗

**Figure 8.** Figure 8: The gri-band diaObjectForcedSource light curves for the nine RR Lyrae in the 47Tuc field. The cyan curves represent the best-fit template light curves. The units for period P and metallicity [Fe/H] are in days and dex, respectively. 20.50 20.75 21.00 21.25 21.50 21.75 22.00 22.25 diaObjectForcedSource mean mags 20.50 20.75 21.00 21.25 21.50 21.75 22.00 22.25 objectForcedSource mean mags r-band i-band 0.2 0… view at source ↗

**Figure 9.** Figure 9: Left Panel: Comparison of the mean magnitudes returned from fitting the objectForcedSource light curves and the diaObjectForcedSource light curves with the template light curves to the known RR Lyrae in the Fornax field. The solid magenta line represents the 1:1 relation. Middle Panel: Similar to the left panel, but for the r-band amplitudes determined from the best-fit template light curves. We did not co… view at source ↗

**Figure 10.** Figure 10: Left Panel: Comparison of the derived µ, using the mean magnitudes obtained from the diaObjectForcedSource light curves and the theoretical PWZ and PLZ relations (Marconi et al. 2022), to the literature values adopted from Li et al. (2023, whenever available). The straight line represents the 1 : 1 relation. Right Panel: Similar to the left panel, but for the RR Lyrae in the DES sample with literature val… view at source ↗

**Figure 11.** Figure 11: Comparison of µW derived in this work with literature values, including the eight RR Lyrae in 47Tuc field (after excluding the extreme outlier as shown in the left panel of view at source ↗

read the original abstract

We investigate the properties of known RR Lyrae in the Vera C. Rubin Observatory Data Preview 1 (DP1) fields and compare those with the predictions based on stellar pulsation models tailored to the Legacy Survey of Space and Time (LSST) filters. The cross-match of the DP1 data with two public variable star catalogs resulted in $\sim 600$ RR Lyrae with adequate light curve sampling in five (out of seven) DP1 fields. The majority of RR Lyrae are in the 47 Tucanae and Fornax fields. We estimated photometric metallicities for these RR Lyrae using the theoretical metallicity-color relation based on $gri$-band data, and find a good agreement with literature values where the light curve sampling is sufficient for fitting template light curves accurately. The distance modulus to all RR Lyrae in DP1 fields were determined using the theoretical period-luminosity-metallicity (PLZ) relations and the $W_{gr}$ period-Wesenheit-metallicity (PWZ) relation which has the smallest metallicity term. The distances based on PWZ relations are in good agreement with the literature values with a mean offset of $0.01\pm0.36$~mag. However, the PLZ-based distance moduli are systematically large which could be due to the theoretical calibration uncertainties that include evolved horizontal branch models. The predicted period-amplitude relations based on evolved models are also inconsistent with the amplitudes based on DP1 light curves. We conclude that the metallicity and distance estimates are sensitive to the template fitting to sparsely sampled light curves in DP1 data and future data release will significantly improve these determinations for RR Lyrae stars.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This is a useful early check on RR Lyrae distances in Rubin DP1 data that confirms PWZ relations work but flags real problems with sparse sampling and PLZ models.

read the letter

The paper cross-matches known RR Lyrae catalogs against the DP1 fields and ends up with about 600 stars that have enough points for template fitting, mostly in 47 Tucanae and Fornax. They derive photometric metallicities from the gri colors using an existing theoretical relation, then compute distances with both PLZ and PWZ relations calibrated on stellar models. The PWZ distances line up with literature values at a mean offset of 0.01 ± 0.36 mag, while the PLZ ones are systematically too large and the observed amplitudes do not match the evolved horizontal-branch predictions. The authors note that both metallicity and distance results are sensitive to how well the templates fit the sparsely sampled light curves, and they stop there without claiming the numbers are definitive for LSST use. That is the main new content: a concrete, quantitative comparison on actual early survey photometry rather than simulations. The work is straightforward and appropriately caveated, which is what makes it worth reading. The soft spots are exactly what you would expect from preview data. The sample is small and field-biased, the light curves are undersampled by design, and everything rests on pre-existing theoretical calibrations instead of fitting new parameters to the DP1 points. No new technique or derivation appears. The reported scatter and systematic offsets are presented clearly, but without the full fitting code, covariance treatment, or per-star diagnostics it is difficult to judge how much of the PLZ discrepancy is real versus fitting artifact. This is the kind of paper that belongs in the LSST preparatory literature. Readers who are building pipelines or distance ladders for the full survey will get practical warnings about template fitting and model choice. It is not a high-impact theoretical result, but the empirical sanity check is honest and timely. I would send it to peer review. The central claims are modest and internally consistent, and a referee can ask for clearer documentation of the fitting steps without needing major new data.

Referee Report

2 major / 2 minor

Summary. The manuscript analyzes ~600 RR Lyrae stars identified in five Vera C. Rubin Observatory DP1 fields via cross-matching with public variable-star catalogs. Photometric metallicities are derived from theoretical gri-band metallicity-color relations, and distance moduli are computed using theoretical PLZ and PWZ relations (with the latter having the smallest metallicity coefficient). PWZ distances show a mean offset of 0.01±0.36 mag relative to literature values, while PLZ distances are systematically larger; period-amplitude predictions from evolved horizontal-branch models also mismatch the DP1 amplitudes. The authors conclude that both metallicity and distance estimates are sensitive to template fitting on the sparsely sampled DP1 light curves.

Significance. If the quantitative comparisons hold, the work supplies early empirical benchmarks for theoretical pulsation models calibrated to LSST gri filters against real preview photometry. The small PWZ offset and the explicit caveats on sampling limitations provide practical guidance for RR Lyrae science with the full LSST survey, where denser light-curve coverage will mitigate the reported sensitivities.

major comments (2)

[Abstract and distance results] The central explanation that PLZ distances are systematically larger because of theoretical calibration uncertainties (including evolved horizontal-branch models) is load-bearing for the interpretation, yet the manuscript provides no quantitative comparison of the PLZ zero-point or slope shift between evolved and non-evolved model sets for the specific LSST gri passbands (see abstract and the distance-modulus results paragraph).
[Methods and conclusions] The conclusion that metallicity and distance estimates are sensitive to template fitting on sparse DP1 light curves is well-motivated by the reported scatter, but the paper does not quantify the fitting uncertainties (e.g., via bootstrap or Monte-Carlo resampling of the template parameters) or show how the derived periods, amplitudes, and mean magnitudes change with different sampling thresholds.

minor comments (2)

The abstract states that the majority of the ~600 stars lie in the 47 Tucanae and Fornax fields; adding the exact field-by-field counts and the two fields excluded for inadequate sampling would improve transparency.
The reported mean offset of 0.01±0.36 mag for PWZ distances is quoted without specifying whether the uncertainty is the standard error of the mean or the rms scatter; clarifying this and showing the distribution of residuals would strengthen the claim of 'good agreement'.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful review and the recommendation for minor revision. We address the major comments below and indicate the changes we will make to the manuscript.

read point-by-point responses

Referee: [Abstract and distance results] The central explanation that PLZ distances are systematically larger because of theoretical calibration uncertainties (including evolved horizontal-branch models) is load-bearing for the interpretation, yet the manuscript provides no quantitative comparison of the PLZ zero-point or slope shift between evolved and non-evolved model sets for the specific LSST gri passbands (see abstract and the distance-modulus results paragraph).

Authors: We agree that a direct quantitative comparison of the PLZ zero-point and slope between evolved and non-evolved model sets for the LSST gri passbands would strengthen the interpretation. The PLZ relations used here are taken from pulsation models that incorporate evolved horizontal-branch stars, and the observed systematic offset is presented only as a possible cause. In the revised manuscript we will qualify the relevant sentences in the abstract and distance-modulus paragraph to make this explicit, citing the model papers for the expected differences, rather than leaving the statement unqualified. revision: partial
Referee: [Methods and conclusions] The conclusion that metallicity and distance estimates are sensitive to template fitting on sparse DP1 light curves is well-motivated by the reported scatter, but the paper does not quantify the fitting uncertainties (e.g., via bootstrap or Monte-Carlo resampling of the template parameters) or show how the derived periods, amplitudes, and mean magnitudes change with different sampling thresholds.

Authors: We acknowledge that explicit quantification of the uncertainties would improve the robustness of the conclusion. The present analysis infers sensitivity from the scatter relative to literature values. In the revised version we will add a Monte-Carlo resampling of the template fits for a representative subset of stars and will include a brief analysis (text and/or table) showing how periods, amplitudes, and mean magnitudes vary with sampling quality thresholds. These additions will be placed in the methods and results sections. revision: yes

Circularity Check

0 steps flagged

No significant circularity; external theoretical relations applied to new data with independent comparisons

full rationale

The paper applies pre-existing theoretical PLZ, PWZ, and metallicity-color relations (calibrated on stellar pulsation models including evolved horizontal branch stars) to DP1 photometry for ~600 RR Lyrae stars. Metallicities and distance moduli are computed from these relations, then compared to independent literature values, yielding a PWZ offset of 0.01±0.36 mag while noting PLZ discrepancies and amplitude mismatches. No parameters of the relations are fitted to the DP1 data; the analysis instead flags sensitivity to sparse light-curve template fitting and attributes differences to model uncertainties. The derivation chain therefore remains self-contained against external benchmarks and does not reduce any claimed result to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Analysis rests on pre-existing theoretical relations for RR Lyrae pulsation and metallicity-color; no new free parameters or entities introduced.

axioms (2)

domain assumption Theoretical metallicity-color relation based on gri-band data for RR Lyrae
Invoked to estimate photometric metallicities from DP1 photometry.
domain assumption Theoretical PLZ and PWZ relations calibrated including evolved horizontal branch models
Used to compute distance moduli; noted as possible source of systematic offset.

pith-pipeline@v0.9.0 · 5613 in / 1314 out tokens · 39841 ms · 2026-05-09T19:27:31.644571+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

29 extracted references · 28 canonical work pages · 1 internal anchor

[1]

1999, XSPEC: An X-ray spectral fitting package,, Astrophysics Source Code Library, record ascl:9910.005 http://ascl.net/9910.005 Astropy Collaboration, Robitaille, T

8 http://www.astropy.org Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33, doi: 10.1051/0004-6361/201322068 Astropy Collaboration, Price-Whelan, A. M., Sip˝ ocz, B. M., et al. 2018, AJ, 156, 123, doi: 10.3847/1538-3881/aabc4f 14 Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022, ApJ, 935, 167, doi: 10...

work page doi:10.1051/0004-6361/201322068 2013
[2]

L., Bono, G., Braga, V

Beaton, R. L., Bono, G., Braga, V. F., et al. 2018, SSRv, 214, 113, doi: 10.1007/s11214-018-0542-1

work page doi:10.1007/s11214-018-0542-1 2018
[3]

Bersier, D., & Wood, P. R. 2002, AJ, 123, 840, doi: 10.1086/338315

work page doi:10.1086/338315 2002
[4]

2020, Journal of Astrophysics and Astronomy, 41, 23, doi: 10.1007/s12036-020-09640-z

Bhardwaj, A. 2020, Journal of Astrophysics and Astronomy, 41, 23, doi: 10.1007/s12036-020-09640-z

work page doi:10.1007/s12036-020-09640-z 2020
[5]

2024, AJ, 167, 247, doi: 10.3847/1538-3881/ad38b6

Bhardwaj, A., Rejkuba, M., Ngeow, C.-C., et al. 2024, AJ, 167, 247, doi: 10.3847/1538-3881/ad38b6

work page doi:10.3847/1538-3881/ad38b6 2024
[6]

F., Monelli, M., Dall’Ora, M., et al

Braga, V. F., Monelli, M., Dall’Ora, M., et al. 2024, A&A, 689, A349, doi: 10.1051/0004-6361/202450971

work page doi:10.1051/0004-6361/202450971 2024
[7]

F., Fiorentino, G., Bono, G., et al

Braga, V. F., Fiorentino, G., Bono, G., et al. 2022, MNRAS, 517, 5368, doi: 10.1093/mnras/stac2813

work page doi:10.1093/mnras/stac2813 2022
[8]

2025, arXiv e-prints, arXiv:2501.02103, doi: 10.48550/arXiv.2501.02103

Caplar, N., Beebe, W., Branton, D., et al. 2025, arXiv e-prints, arXiv:2501.02103, doi: 10.48550/arXiv.2501.02103

work page doi:10.48550/arxiv.2501.02103 2025
[9]

A., Clayton, G

Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245, doi: 10.1086/167900

work page doi:10.1086/167900 1989
[10]

Choi, Y., Olsen, K. A. G., Carlin, J. L., et al. 2025, ApJ, 992, 47, doi: 10.3847/1538-4357/adfb70 Dark Energy Survey Collaboration, Abbott, T., Abdalla, F. B., et al. 2016, MNRAS, 460, 1270, doi: 10.1093/mnras/stw641 D´ ek´ any, I., & Grebel, E. K. 2022, ApJS, 261, 33, doi: 10.3847/1538-4365/ac74ba

work page doi:10.3847/1538-4357/adfb70 2025
[11]

W., et al

Feng, Y., Guhathakurta, P., Peng, E. W., et al. 2024, ApJ, 966, 159, doi: 10.3847/1538-4357/ad2ae7

work page doi:10.3847/1538-4357/ad2ae7 2024
[12]

B., et al

Fiorentino, G., Monelli, M., Stetson, P. B., et al. 2017, A&A, 599, A125, doi: 10.1051/0004-6361/201629501 Gaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016, A&A, 595, A1, doi: 10.1051/0004-6361/201629272 Gaia Collaboration, Vallenari, A., Brown, A. G. A., et al. 2023, A&A, 674, A1, doi: 10.1051/0004-6361/202243940

work page doi:10.1051/0004-6361/201629501 2017
[13]

D., Clayton, G

Gordon, K. D., Clayton, G. C., Decleir, M., et al. 2023, ApJ, 950, 86, doi: 10.3847/1538-4357/accb59

work page doi:10.3847/1538-4357/accb59 2023
[14]

Green, G. M. 2018, The Journal of Open Source Software, 3, 695, doi: 10.21105/joss.00695

work page doi:10.21105/joss.00695 2018
[15]

M., Bianco, F

Hambleton, K. M., Bianco, F. B., Street, R., et al. 2023, PASP, 135, 105002, doi: 10.1088/1538-3873/acdb9a

work page doi:10.1088/1538-3873/acdb9a 2023
[16]

G., Rix, H.-W., et al

Hernitschek, N., Cohen, J. G., Rix, H.-W., et al. 2018, ApJ, 859, 31, doi: 10.3847/1538-4357/aabfbb Ivezi´ c,ˇZ., Kahn, S. M., Tyson, J. A., et al. 2019, ApJ, 873, 111, doi: 10.3847/1538-4357/ab042c

work page doi:10.3847/1538-4357/aabfbb 2018
[17]

Photometric Metallicity and Distance Estimates for 136,000 RR Lyrae Stars from Gaia Data Release 3.ApJ2023,944, 88, [arXiv:astro-ph.SR/2206.07668]

Li, X.-Y., Huang, Y., Liu, G.-C., Beers, T. C., & Zhang, H.-W. 2023, ApJ, 944, 88, doi: 10.3847/1538-4357/acadd5

work page doi:10.3847/1538-4357/acadd5 2023
[18]

2025, arXiv e-prints, arXiv:2506.23955, doi: 10.48550/arXiv.2506.23955

Malanchev, K., DeLucchi, M., Caplar, N., et al. 2025, arXiv e-prints, arXiv:2506.23955, doi: 10.48550/arXiv.2506.23955

work page doi:10.48550/arxiv.2506.23955 2025
[19]

2022, ApJ, 934, 29, doi: 10.3847/1538-4357/ac78ee

Marconi, M., Molinaro, R., Dall’Ora, M., et al. 2022, ApJ, 934, 29, doi: 10.3847/1538-4357/ac78ee

work page doi:10.3847/1538-4357/ac78ee 2022
[20]

2022, Universe, 8, 191, doi: 10.3390/universe8030191

Monelli, M., & Fiorentino, G. 2022, Universe, 8, 191, doi: 10.3390/universe8030191

work page doi:10.3390/universe8030191 2022
[21]

R., & El-Badry, K

Nagarajan, P., Weisz, D. R., & El-Badry, K. 2022, ApJ, 932, 19, doi: 10.3847/1538-4357/ac69e6

work page doi:10.3847/1538-4357/ac69e6 2022
[22]

2026, arXiv e-prints, arXiv:2603.24249, doi: 10.48550/arXiv.2603.24249 NSF-DOE Vera C

Ngeow, C.-C., Bhardwaj, A., Nishad, P., & Susmita, D. 2026, arXiv e-prints, arXiv:2603.24249, doi: 10.48550/arXiv.2603.24249 NSF-DOE Vera C. Rubin Observatory. 2025a, Legacy Survey of Space and Time Data Preview 1 [Data set], NSF-DOE Vera C. Rubin Observatory, doi: 10.71929/RUBIN/2570308. https://www.osti.gov//servlets/purl/2570308

work page doi:10.48550/arxiv.2603.24249 2026
[23]

and Finkbeiner, Douglas P

Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103, doi: 10.1088/0004-637X/737/2/103

work page internal anchor Pith review doi:10.1088/0004-637x/737/2/103 2011
[24]

J., Finkbeiner, D

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525, doi: 10.1086/305772 SLAC National Accelerator Laboratory, & NSF-DOE Vera C. Rubin Observatory. 2024, LSST Commissioning

work page Pith review doi:10.1086/305772 1998
[25]

https://www.osti.gov//servlets/purl/2561361

Camera, SLAC National Accelerator Laboratory (SLAC), Menlo Park, CA (United States), doi: 10.71929/RUBIN/2561361. https://www.osti.gov//servlets/purl/2561361

work page doi:10.71929/rubin/2561361
[26]

M., Drlica-Wagner, A., Macri, L., et al

Stringer, K. M., Drlica-Wagner, A., Macri, L., et al. 2021, ApJ, 911, 109, doi: 10.3847/1538-4357/abe873 Vera C Rubin Observatory Team, Acero Cuellar, T.,

work page doi:10.3847/1538-4357/abe873 2021
[27]

2026, arXiv e-prints, arXiv:2603.23786

Acosta, E., et al. 2026, arXiv e-prints, arXiv:2603.23786. https://arxiv.org/abs/2603.23786

work page arXiv 2026
[28]

W., Xue, X

Wang, F., Zhang, H. W., Xue, X. X., et al. 2022, MNRAS, 513, 1958, doi: 10.1093/mnras/stac874

work page doi:10.1093/mnras/stac874 2022
[29]

L., Henden, A

Watson, C. L., Henden, A. A., & Price, A. 2006, Society for Astronomical Sciences Annual Symposium, 25, 47

2006