pith. sign in

arxiv: 2606.18572 · v1 · pith:MZU2EX6Nnew · submitted 2026-06-17 · 🌌 astro-ph.SR

A Pilot Study to Verify the RR Lyrae Candidates with Vera C. Rubin Observatory Early Alerts

Chow-Choong Ngeow , Anupam Bhardwaj , Sarang Shah , Steven Gough-Kelly , Oleksandra Razim This is my paper

Pith reviewed 2026-06-26 20:12 UTC · model grok-4.3

classification 🌌 astro-ph.SR
keywords RR Lyrae starsVera C. Rubin Observatoryearly alertsvariable star verificationlight curvesPan-STARRS1Dark Energy Surveyalert brokers
0
0 comments X

The pith

Rubin Observatory early alerts verify 32 RR Lyrae candidates as genuine using multiband light curves.

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

This pilot study crossmatches faint RR Lyrae candidates from Pan-STARRS1, the Dark Energy Survey, and the Next-Generation Virgo Cluster Survey catalogs with early alerts from the Vera C. Rubin Observatory in deep drilling fields and the M49 field. From the resulting 40 alerts, the extracted multiband light curves confirm 32 as genuine RR Lyrae stars. Several of these confirmed stars were not flagged as variables by the ALeRCE and Lasair community brokers. The remaining alerts include 20 percent non-RR Lyrae variables plus eight non-typical cases such as two active galactic nuclei and two eclipsing binaries. A small number of cataloged candidates produced no alerts at all.

Core claim

The paper shows that multiband light curves derived from Rubin-LSST early alerts verify 32 RR Lyrae candidates as genuine while identifying contaminants, with ALeRCE and Lasair brokers achieving roughly 70 percent and 40 percent true variable classification rates on the sample and 20 percent of alerts turning out to be other variable types.

What carries the argument

Multiband light curves extracted from early alerts, inspected for characteristic RR Lyrae variability patterns to confirm or reject candidates.

If this is right

  • 32 candidates receive confirmation as RR Lyrae on the basis of their extracted light curves.
  • Community brokers miss some genuine variables that alert light curves can still identify.
  • 20 percent of the alert sample consists of non-RR Lyrae variable stars.
  • Eight candidates display light curves inconsistent with RR Lyrae, including AGN and eclipsing binaries.
  • Some cataloged candidates yield no alerts, implying either misclassification or unavailable template images.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Alert-based verification could scale to screen larger variable-star catalogs once full survey operations begin.
  • The method may reduce the volume of dedicated follow-up needed for initial candidate vetting.
  • Missing alerts for some candidates highlight the role of template-image availability in difference-image pipelines.
  • The same crossmatch-and-inspect approach could apply to verification of other periodic variable classes.

Load-bearing premise

The alert-derived multiband light curves alone provide sufficient information to reliably distinguish genuine RR Lyrae from contaminants without additional photometric or spectroscopic follow-up.

What would settle it

Independent spectroscopic observations or full-survey photometry showing that most of the 32 light-curve-verified objects lack the expected radial-velocity or period-amplitude properties of RR Lyrae stars.

Figures

Figures reproduced from arXiv: 2606.18572 by Anupam Bhardwaj, Chow-Choong Ngeow, Oleksandra Razim, Sarang Shah, Steven Gough-Kelly.

Figure 1
Figure 1. Figure 1: Aitoff projection of the six fields listed in [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: The r-band light curve for an alert associated with an RRL from the DES RRL catalog. The light curve has been folded using its pulsation period of 0.6965 days. The inset figure shows the intranight variability for data points highlighted in the box. AGN or bogus from the other two classifiers. They are likely not RRL variables. 3.2. Constraint from Alert Light Curves Using ALeRCE, we downloaded the multiba… view at source ↗
Figure 4
Figure 4. Figure 4: Multiband light curves for two confirmed RRL in the COSMOS field. The dashed curves are the best-fit template light curves (Braga et al. 2024). Note that the COSMOS field is the only field that is observed in all six ugrizy filters in the early Rubin alerts (even though not all RRL have six-band light curves). The Braga et al. (2024) template light curves did not include the u- and y-band, therefore we fit… view at source ↗
Figure 5
Figure 5. Figure 5: Same as [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Same as [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Multiband light curves for the alert (313853517411385444) associated with a known contact binary. Left and right panels are the light curves folded with periods given in the PS1 RRL catalog and in Drake et al. (2014), respectively. There are five candidates in Cases 2 and 3 that are rejected as RRL based on the multiband alert light curves. One of them is a known contact binary (Drake et al. 2014). Using t… view at source ↗
Figure 8
Figure 8. Figure 8: Multiband light curves for candidates in Cases 2 and 3, located in the COSMOS and the M49 field, that are not RR Lyrae. Note that it is unclear whether the large scatter seen in a given phase is intrinsic or due to the intranight variations. estingly, the Gaia Data Release 3 (DR3) RRL catalog (Clementini et al. 2023) classified this candidate as an ab-type RRL but with a period of 0.33451 days. As demonstr… view at source ↗
Figure 9
Figure 9. Figure 9: The multiband light curves for the three candidates in Case 4, which do not exhibit light curve shapes expected for RRL. Hence, they are rejected as RRL candidates. 17.5 18.0 18.5 u COSMOS OID=313936986529333309 Ndet = 759 AGN(0.63)/VS Nu = 16 16.5 17.0 17.5 g Ng = 164 17.0 17.5 r Nr = 137 17.0 17.5 i Ni = 284 17.0 17.5 z Nz = 135 0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 phase (P = 0.35087 days) 16.75 … view at source ↗
Figure 10
Figure 10. Figure 10: Folded multiband light curves for alert 313936986529333309. The left panels show the resulting light curves when folded with the default period adopted from the PS1 RRL catalog. Light curves in the right panels were folded with a shorter period adopted from Chen et al. (2020) or Clementini et al. (2023). The dashed curves are the best-fit template light curves for a c-type RRL. Note that the Braga et al. … view at source ↗
Figure 11
Figure 11. Figure 11: Comparison of the gri-band light curves ex￾tracted from Rubin alerts (red points) and the zubercal online service (grey points). The light curves have been folded with the period adopted from Sesar et al. (2013), which exhibits the least scatter (as compared to the peri￾ods adopted from the PS1 RRL catalog, or other published periods). Based on the zubercal light curves, this variable is most likely a Bla… view at source ↗
read the original abstract

We present a pilot study using the Vera C. Rubin Observatory early alerts to verify RR Lyrae candidates in the Pan-STARRS1, the Dark Energy Survey, and the Next-Generation Virgo Cluster Survey RR Lyrae catalogs. RR Lyrae candidates fainter than 16 mag in the $g$-band in these catalogs were crossmatched with the alerts observed in several deep drilling fields and the M49 field. After excluding alerts with a low number of detections, there are 40 alerts associated with the RR Lyrae candidates. The multiband Rubin-LSST light curves extracted from the alerts verify 32 variables as genuine RR Lyrae, although several were not classified as variable stars in the ALeRCE and Lasair community alert brokers. While ALeRCE and Lasair provide $\sim70\%$ and $\sim40\%$ true variable classification, respectively, we find that $20\%$ of the alert sample are non-RR Lyrae variables. The remaining eight candidate variables do not show typical RR Lyrae light curves and include two active galactic nuclei and two eclipsing binaries. Additionally, we have also found a small number of known variable candidates with no alerts, which would suggest that they are either not RR Lyrae variables or the template images are not yet available for their difference image analysis.

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.

Referee Report

2 major / 1 minor

Summary. The paper reports a pilot study cross-matching faint RR Lyrae candidates from Pan-STARRS1, DES, and NGVS catalogs against Vera C. Rubin Observatory early alerts in deep-drilling fields and the M49 field. After excluding low-detection alerts, 40 objects remain; multiband alert light curves are inspected to classify 32 as genuine RR Lyrae (some missed by ALeRCE and Lasair brokers), while 8 are deemed non-RR Lyrae (including two AGN and two eclipsing binaries). The work also notes a handful of catalogued candidates lacking alerts, possibly due to missing templates, and compares broker true-positive rates (~70% ALeRCE, ~40% Lasair).

Significance. If the verification procedure can be made reproducible, the study supplies a concrete early demonstration that Rubin alert streams can recover and vet RR Lyrae candidates from existing catalogs, quantifies broker performance on a small but clean sample, and flags the practical impact of difference-imaging template availability. These are useful benchmarks for LSST-era variable-star pipelines even though the sample size precludes broad statistical conclusions.

major comments (2)
  1. [Abstract] Abstract (verification paragraph): the central claim that 32 of 40 alerts are verified as genuine RR Lyrae rests on the statement that they “show typical RR Lyrae light curves,” yet no quantitative criteria (Lomb-Scargle FAP threshold, minimum epochs per band, Fourier coefficients, template-fit χ², or period/amplitude consistency rules) are supplied. Without these, the 80 % verification fraction and the 20 % non-RR Lyrae fraction cannot be independently assessed or reproduced.
  2. [Abstract] Abstract (non-detection paragraph): the observation that some known candidates produce no alerts is attributed to possible missing templates, but the text provides neither the total number of such candidates examined nor an estimate of how template incompleteness biases the recovered sample or the reported verification statistics.
minor comments (1)
  1. [Abstract] The abstract would be clearer if it stated the exact number of deep-drilling fields and the M49 field coordinates or visit counts used for the cross-match.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the constructive comments on our pilot study. We address each major comment below and indicate where revisions will be made.

read point-by-point responses

  1. Referee: [Abstract] Abstract (verification paragraph): the central claim that 32 of 40 alerts are verified as genuine RR Lyrae rests on the statement that they “show typical RR Lyrae light curves,” yet no quantitative criteria (Lomb-Scargle FAP threshold, minimum epochs per band, Fourier coefficients, template-fit χ², or period/amplitude consistency rules) are supplied. Without these, the 80 % verification fraction and the 20 % non-RR Lyrae fraction cannot be independently assessed or reproduced.

    Authors: Classification was performed via visual inspection of multiband light curves for characteristic RR Lyrae morphology (asymmetric sawtooth shapes, periods ~0.2-1 d, amplitudes 0.2-1.5 mag) and consistency across bands. No automated quantitative metrics were applied in this small pilot. We will revise the text to explicitly describe these qualitative criteria and acknowledge the absence of reproducible thresholds as a limitation of the current work. revision: partial

  2. Referee: [Abstract] Abstract (non-detection paragraph): the observation that some known candidates produce no alerts is attributed to possible missing templates, but the text provides neither the total number of such candidates examined nor an estimate of how template incompleteness biases the recovered sample or the reported verification statistics.

    Authors: We will revise the abstract and main text to report the exact number of catalogued candidates examined in the fields. A quantitative bias estimate from template incompleteness cannot be provided because template availability is dynamic and incompletely documented for the early-alert period; we will note this as an unquantified caveat. revision: partial

standing simulated objections not resolved
  • Quantitative estimate of how template incompleteness biases the recovered sample or the reported verification statistics

Circularity Check

0 steps flagged

No circularity: purely observational cross-match and visual classification

full rationale

The paper reports a catalog cross-match between existing RR Lyrae candidate lists and Rubin early alerts, followed by extraction of multiband light curves and qualitative visual inspection to confirm 32 objects show typical RR Lyrae shapes. No equations, fitted parameters, model derivations, or predictions appear in the provided text. The central claim rests on direct comparison to external catalogs and observed light-curve morphology, with no self-referential reduction or load-bearing self-citation. This is standard observational verification and is self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Observational pilot study with no mathematical model; relies on standard domain assumptions about cross-match reliability and light-curve morphology for RR Lyrae classification.

axioms (1)
  • domain assumption Standard astronomical assumptions on positional cross-matching accuracy and on the uniqueness of RR Lyrae light-curve shapes in multiband photometry.
    Invoked when the paper states that extracted light curves verify candidates as genuine RR Lyrae.

pith-pipeline@v0.9.1-grok · 5792 in / 1295 out tokens · 31762 ms · 2026-06-26T20:12:04.464864+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

34 extracted references · 33 canonical work pages · 1 internal anchor

  1. [1]

    2026, A&A, 705, A247, doi: 10.1051/0004-6361/202556258 Astropy Collaboration, Robitaille, T

    Ar´ evalo, P., S´ anchez-S´ aez, P., Sotomayor, B., et al. 2026, A&A, 705, A247, doi: 10.1051/0004-6361/202556258 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 Astro...

    work page doi:10.1051/0004-6361/202556258 2026

  2. [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. [3]

    PASP , author =

    Bellm, E. C., Kulkarni, S. R., Graham, M. J., et al. 2019, PASP, 131, 018002, doi: 10.1088/1538-3873/aaecbe

    work page doi:10.1088/1538-3873/aaecbe 2019

  4. [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 —. 2022, Universe, 8, 122, doi: 10.3390/universe8020122

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

  5. [5]

    , keywords =

    Bianco, F. B., Ivezi´ c,ˇZ., Jones, R. L., et al. 2022, ApJS, 258, 1, doi: 10.3847/1538-4365/ac3e72

    work page doi:10.3847/1538-4365/ac3e72 2022

  6. [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. [7]

    M., Mendes de Oliveira, C., Krabbe, A

    Cardoso, N. M., Mendes de Oliveira, C., Krabbe, A. C., et al. 2026, AJ, 171, 290, doi: 10.3847/1538-3881/ae4fb1

    work page doi:10.3847/1538-3881/ae4fb1 2026

  8. [8]

    2021, AJ, 162, 231, doi: 10.3847/1538-3881/ac0ef1

    Carrasco-Davis, R., Reyes, E., Valenzuela, C., et al. 2021, AJ, 162, 231, doi: 10.3847/1538-3881/ac0ef1

    work page doi:10.3847/1538-3881/ac0ef1 2021

  9. [9]

    A., Cheung, S.-H., et al

    Chan, H.-S., Villar, V. A., Cheung, S.-H., et al. 2022, ApJ, 932, 118, doi: 10.3847/1538-4357/ac69d4

    work page doi:10.3847/1538-4357/ac69d4 2022

  10. [10]

    2023, MNRAS, 525, 3075, doi: 10.1093/mnras/stad2296

    Chen, A., Li, Z., Wang, Y., et al. 2023, MNRAS, 525, 3075, doi: 10.1093/mnras/stad2296

    work page doi:10.1093/mnras/stad2296 2023

  11. [11]

    2020, ApJS, 249, 18, doi: 10.3847/1538-4365/ab9cae

    Chen, X., Wang, S., Deng, L., et al. 2020, ApJS, 249, 18, doi: 10.3847/1538-4365/ab9cae

    work page doi:10.3847/1538-4365/ab9cae 2020

  12. [12]

    and Ripepi, V

    Clementini, G., Ripepi, V., Garofalo, A., et al. 2023, A&A, 674, A18, doi: 10.1051/0004-6361/202243964

    work page doi:10.1051/0004-6361/202243964 2023

  13. [13]

    J., Catelan, M., Djorgovski, S

    Drake, A. J., Catelan, M., Djorgovski, S. G., et al. 2013, ApJ, 763, 32, doi: 10.1088/0004-637X/763/1/32

    work page doi:10.1088/0004-637x/763/1/32 2013

  14. [14]

    J., Graham, M

    Drake, A. J., Graham, M. J., Djorgovski, S. G., et al. 2014, ApJS, 213, 9, doi: 10.1088/0067-0049/213/1/9

    work page doi:10.1088/0067-0049/213/1/9 2014

  15. [15]

    W., et al

    Feng, Y., Guhathakurta, P., Peng, E. W., et al. 2026, ApJ, 998, 157, doi: 10.3847/1538-4357/ae313b —. 2024, ApJ, 966, 159, doi: 10.3847/1538-4357/ad2ae7 F¨ orster, F., Cabrera-Vives, G., Castillo-Navarrete, E., et al. 2021, AJ, 161, 242, doi: 10.3847/1538-3881/abe9bc

    work page doi:10.3847/1538-4357/ae313b 2026

  16. [16]

    and Kulkarni, S

    Graham, M. J., Kulkarni, S. R., Bellm, E. C., et al. 2019, PASP, 131, 078001, doi: 10.1088/1538-3873/ab006c

    work page doi:10.1088/1538-3873/ab006c 2019

  17. [17]

    P., et al., 2025, RTN-011: Rubin Observatory Plans for an Early Science Program, @doi 10.71929/RUBIN/2584021 , https://www.osti.gov/servlets/purl/2584021

    Guy, L. P., Bechtol, K., Bellm, E. C., et al. 2025, RTN-011: Rubin Observatory Plans for an Early Science Program, Tech. rep., NSF-DOE Vera C. Rubin Observatory, doi: 10.71929/RUBIN/2584021

    work page doi:10.71929/rubin/2584021 2025

  18. [18]

    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

    work page doi:10.3847/1538-4357/aabfbb 2018

  19. [19]

    Huang, K.-W., & Koposov, S. E. 2022, MNRAS, 510, 3575, doi: 10.1093/mnras/stab3654

    work page doi:10.1093/mnras/stab3654 2022

  20. [20]

    2021, MNRAS, 502, 5686, doi: 10.1093/mnras/stab005 Ivezi´ c,ˇZ., Kahn, S

    Iorio, G., & Belokurov, V. 2021, MNRAS, 502, 5686, doi: 10.1093/mnras/stab005 Ivezi´ c,ˇZ., Kahn, S. M., Tyson, J. A., et al. 2019, ApJ, 873, 111, doi: 10.3847/1538-4357/ab042c Juri´ c, M., Axelrod, T. S., Becker, A. C., et al. 2023, LSE-163: Data Products Definition Document, Tech. rep., NSF-DOE Vera C. Rubin Observatory, doi: 10.71929/RUBIN/2587118

    work page doi:10.1093/mnras/stab005 2021

  21. [21]

    M., et al

    Lucey, M., Mateu, C., Price-Whelan, A. M., et al. 2026, AJ, 171, 249, doi: 10.3847/1538-3881/ae4aaa

    work page doi:10.3847/1538-3881/ae4aaa 2026

  22. [22]

    E., Mu˜ noz, R

    Medina, G. E., Mu˜ noz, R. R., Carlin, J. L., et al. 2024, MNRAS, 531, 4762, doi: 10.1093/mnras/stae1137

    work page doi:10.1093/mnras/stae1137 2024

  23. [23]

    G., Minniti, D., Capuzzo-Dolcetta, R., et al

    Navarro, M. G., Minniti, D., Capuzzo-Dolcetta, R., et al. 2021, A&A, 646, A45, doi: 10.1051/0004-6361/202038463

    work page doi:10.1051/0004-6361/202038463 2021

  24. [24]

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

    Ngeow, C.-C., & Bhardwaj, A. 2026, arXiv e-prints, arXiv:2605.00344, doi: 10.48550/arXiv.2605.00344

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2605.00344 2026

  25. [25]

    M., Becker, A

    Oluseyi, H. M., Becker, A. C., Culliton, C., et al. 2012, AJ, 144, 9, doi: 10.1088/0004-6256/144/1/9 S´ anchez-S´ aez, P., Reyes, I., Valenzuela, C., et al. 2021, AJ, 161, 141, doi: 10.3847/1538-3881/abd5c1 S´ anchez-S´ aez, P., Arredondo, J., Bayo, A., et al. 2023, A&A, 675, A195, doi: 10.1051/0004-6361/202346077

    work page doi:10.1088/0004-6256/144/1/9 2012

  26. [26]

    S., et al

    Sesar, B., Ivezi´ c,ˇZ., Stuart, J. S., et al. 2013, AJ, 146, 21, doi: 10.1088/0004-6256/146/2/21

    work page doi:10.1088/0004-6256/146/2/21 2013

  27. [27]

    Machine-learned Identification of RR Lyrae Stars from Sparse, Multi-band Data: The PS1 Sample

    Sesar, B., Hernitschek, N., Mitrovi´ c, S., et al. 2017, AJ, 153, 204, doi: 10.3847/1538-3881/aa661b SLAC National Accelerator Laboratory, & NSF-DOE Vera C. Rubin Observatory. 2025, The LSST Camera (LSSTCam), SLAC National Accelerator Laboratory (SLAC), Menlo Park, CA (United States), doi: 10.71929/rubin/2571927. https://www.osti.gov//servlets/purl/2571927 13

    work page doi:10.3847/1538-3881/aa661b 2017

  28. [28]

    M., Long, J

    Stringer, K. M., Long, J. P., Macri, L. M., et al. 2019, AJ, 158, 16, doi: 10.3847/1538-3881/ab1f46

    work page doi:10.3847/1538-3881/ab1f46 2019

  29. [29]

    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

  30. [30]

    2026, arXiv e-prints, arXiv:2603.23786

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

    arXiv 2026

  31. [31]

    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

  32. [32]

    D., Francis, G

    Williams, R. D., Francis, G. P., Lawrence, A., et al. 2024, RAS Techniques and Instruments, 3, 362, doi: 10.1093/rasti/rzae024

    work page doi:10.1093/rasti/rzae024 2024

  33. [33]

    Young, D. R. 2023, Sherlock. Contextual classification of astronomical transient sources, doi: 10.5281/zenodo.8038057. https://zenodo.org/doi/10.5281/zenodo.8038057

    work page doi:10.5281/zenodo.8038057 2023

  34. [34]

    K., et al

    Zinn, R., Horowitz, B., Vivas, A. K., et al. 2014, ApJ, 781, 22, doi: 10.1088/0004-637X/781/1/22

    work page doi:10.1088/0004-637x/781/1/22 2014