{"id":"89ac1d7a-2df5-40a5-966d-c237d124618f","arxiv_id":"2607.29597","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"The largest homogeneous catalog of RR Lyrae radial velocities (17,563 stars), with template-based systemic velocities, amplitude scaling relations, and metallicity/Blazhko trends.","lead":"This paper assembles radial-velocity measurements for more than 17,000 RR Lyrae stars and uses template fitting to derive their bulk motions and pulsation amplitudes. It is a reference catalog plus empirical relations that will inform Galactic-structure and stellar-pulsation work.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Template V_gamma values are not consistent across Balmer diagnostics; Fig. 5 shows Hα–Hβ mean offsets of ~7 km/s and Hα–Hγ offsets ~14 km/s, exceeding the claimed 3.8 km/s accuracy, so the catalog's homogeneity claim is unsupported.","rationale":"The reader correctly identified the universal-template assumption as the weakest point. My review finds direct evidence in the paper's own Fig. 5 that this assumption fails across diagnostics: for the same stars, Hα and Hβ fixed-amplitude V_gamma means differ by 7.1 km/s, and Hα vs Hγ differ by 14.3 km/s. These offsets are several times the quoted 3.8 km/s accuracy and cannot be dismissed as random error. The paper's explanation that the differences arise from distinct RVC shapes is a statement of the problem, not a resolution, because V_gamma should be invariant. The final catalog's global mean over diagnostics inherits these offsets, making the catalog's zero-point dependent on which lines were measured for each star. Appendix A resampling validates Fe-line templates against Gaia and du Pont, but it does not test Balmer cross-diagnostic consistency, nor does it test stars with fewer than three measurements, so the 11.3 km/s claim for sparse sampling is unvalidated. Since the central claim of homogeneity is load-bearing for the catalog's use as a reference, this concern warrants the CONDITIONAL verdict. The proposed test—quantifying the distribution of star-by-star Hα-Hβ differences and checking metallicity dependence—would settle whether the offsets are a template artifact or a physical line-formation effect. If the offsets are real, the authors can correct them; if not, the homogeneity claim stands. Thus the reader's conditionality is appropriate, and the verdict remains CONDITIONAL (UNCHANGED).","tokens_in":30857,"tokens_out":7606,"duration_ms":75093,"concrete_test":"Compute, for the sample of stars with both Hα and Hβ V_gamma from fixed-amplitude templates (N≈4,490 in Fig. 5), the star-by-star difference Δ = V_gamma(Hα) - V_gamma(Hβ). If the median/mean of Δ is >2 km/s and the distribution does not center at zero, then the diagnostics are inconsistent. Additionally, bin Δ by [Fe/H] and period; if Δ correlates with [Fe/H], the template universality assumption fails and the Section 4.4 metallicity slopes are suspect. A clean check: use the du Pont full-phase RVCs to measure V_gamma independently from each Balmer line without templates; if the same offsets appear, the offsets are physical and should be corrected in the catalog; if they disappear, the template fitting is the source.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of a homogeneous, accurate V_gamma catalog rests on the universal RVC templates. Section 4.1 and Fig. 5 reveal a direct contradiction: for the ~4,490 stars with both Hα and Hβ fixed-amplitude V_gamma, the mean values differ by 7.1 km/s (-48.1 vs -55.2), and Hα versus Hγ differ by 14.3 km/s (-47.3 vs -61.6). The paper explains these as arising from different RVC shapes, but V_gamma is the barycentric velocity and should be line-independent. If the final adopted V_gamma is the global mean over available diagnostics, stars observed only in Hα will be systematically offset from those observed only in Hβ by ~7 km/s, and residual offsets propagate into the catalog as a function of spectral coverage. The quoted accuracy (3.8 km/s for well-sampled RVCs, 6.5 for 3-7 points) is therefore not supported for stars with Balmer-only measurements. Appendix A resampling tests only validate Fe-line templates against Gaia/du Pont, not cross-diagnostic consistency. This undermines both the homogeneity claim and the metallicity-dependent trends (Section 4.4) if templates themselves are metallicity-dependent. The issue is correctable: per-diagnostic offsets can be calibrated using common stars, but until then the catalog's systematic floor is larger than stated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper compiles radial-velocity (RV) measurements for 17,563 field RR Lyrae stars from DESI, LAMOST, SDSS, Gaia, STELLA, du Pont, and the literature. Using the Braga et al. (2021) RVC templates, it derives barycentric velocities V_gamma and RV amplitudes with fixed- and free-amplitude fitting, and combines them with literature values to produce a final catalog of 15,955 stars. The method is validated against Gaia DR3 and du Pont calibration samples, and resampling tests in Appendix A show that template fitting outperforms simple averaging. The paper then uses the derived amplitudes to construct photometric–RV amplitude scaling relations, spectroscopic Bailey diagrams, metallicity trends, and a Blazhko candidate analysis.","tokens_in":31216,"tokens_out":6461,"duration_ms":67445,"significance":"If the homogeneity and accuracy claims hold, this would indeed be a valuable reference catalog: the largest homogeneous set of RRL V_gamma values and RV amplitudes, with immediate applications to Galactic kinematics and as a training set for future multi-object spectroscopic surveys. The paper's strengths include the use of a free-amplitude template method that does not assume a photometric amplitude, external validation on two calibration samples, and quantitative resampling experiments demonstrating the advantage of template fitting over simple means. However, the central homogeneity claim is weakened by large cross-diagnostic V_gamma offsets that are visible in the paper's own Fig. 5 and are not addressed by the calibration samples. The accuracy figures quoted in the abstract are therefore not fully supported for Balmer-only or sparsely sampled stars.","major_comments":[{"comment":"The claim that different diagnostics give 'very similar V_gamma' is not supported by the paper's own numbers. For the same stars, the mean V_gamma from Hα is -48.1 km/s vs -55.2 from Hβ (Δ=7.1 km/s), -47.3 vs -61.6 from Hγ (Δ=14.3 km/s), and -46.6 vs -57.9 from Hδ (Δ=11.3 km/s). Since the final V_gamma is the global mean over available diagnostics, stars observed only in Hα will be systematically offset from stars observed only in Hβ by ~7 km/s. The Appendix A validation is based on Fe-line templates against Gaia/du Pont and does not test Balmer V_gamma zero-points. The quoted 3.8/6.5/11.3 km/s accuracies are therefore not established for Balmer-only RVCs. Please calibrate per-diagnostic offsets using common stars, or restrict the homogeneity/accuracy claims to the validated Fe-line subset.","section":"§4.1, Fig. 5"},{"comment":"The abstract's accuracy for 'fewer than three phase points' (11.3 km/s) is not supported by the validation presented. Section 3 states that templates can be applied even with a single RV measurement, but Appendix A's resampling test draws subsets of 3, 4, 5, 6, 9, and 12 points; no test covers 1 or 2 points. The 11.3 km/s figure therefore appears to be an extrapolation. Please either provide a validation for 1–2 point cases (e.g., Monte Carlo with known V_gamma) or revise the claimed accuracy and the associated catalog uncertainties for such stars.","section":"§3 and Appendix A"},{"comment":"The scaling relations rest on subjective sample cuts. In particular, the Balmer-line relations are derived after 'neglect[ing] 1195 stars from the SDSS catalog ... since they show a larger scatter' with no quantitative criterion; the Fe-line relations additionally exclude all stars with Amp(RV) > 120 km/s (RRab) or > 50 km/s (RRc). These cuts can bias the fitted slopes and the inferred break locations. Please quantify the impact: show fits with and without the SDSS exclusion, report the number of stars removed by each cut, and justify the amplitude caps by an objective outlier criterion rather than by the range that makes the relation appear linear/logarithmic.","section":"§4.2, Figs. 6–8"}],"minor_comments":[{"comment":"The Blazhko analysis uses non-modulated templates and labels the resulting amplitudes as a 'first-order approximation.' This caveat should be repeated when drawing the metallicity-skew conclusion, since the amplitude bias may correlate with period and metallicity.","section":"§4.3"},{"comment":"The SP relation for RRab stars includes a log(P)^4 term but no log(P)^3 term. Please confirm that this is intentional and not a typographical omission.","section":"Eq. (5)"},{"comment":"The DESI resolution entry is listed as '2000–5500' but the notes state that R varies from 2000 at 3600 Å to 5500 at 9800 Å. Please make the wavelength dependence explicit in the table column heading.","section":"Table C.1"}],"recommendation":"major_revision","confidential_remarks":"The cross-diagnostic V_gamma offsets in Fig. 5 are the main scientific risk. They are likely correctable with per-diagnostic zero-point calibrations, but until that is done the catalog's headline homogeneity and accuracy claims are not supported. I do not see grounds for rejection if the authors address this point and the accuracy-validation gap for sparse sampling."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a look if you work with RR Lyrae kinematics, but don't put too much weight on the homogeneity claim yet. The catalog is real and large, but the paper's own Figure 5 undercuts the 'same V_gamma within the errors' statement: for the same stars, Hα and Hβ give mean systemic velocities 7 km/s apart, and Hα and Hγ are 14 km/s apart. Barycentric velocity shouldn't care which line you fit. Since the adopted V_gamma is the mean over whatever diagnostics each star happens to have, you get systematic offsets as a function of spectral coverage. The 3.8/6.5/11.3 km/s accuracy figures come from Fe-line validation only (Appendix A), so they don't cover Balmer-only stars.\n\nWhat's genuinely good: the assembly of 15,955 RRLs from many surveys, the template-fitting framework, and the validation against Gaia and du Pont. For Fe-line V_gamma, template fitting clearly beats simple averages. That's solid work and the catalog will likely be the training set for future surveys.\n\nSoft spots beyond the cross-diagnostic offsets: the Balmer scaling relations drop 1,195 SDSS stars post-hoc because they 'show larger scatter'—that's a selection cut that needs an objective rationale. The metallicity slopes for Balmer lines rest on 33 stars and are only 2–3 sigma; fine as a hint, not as a result. And the Blazhko metallicity skew: the body says no significant offset in the mean, only that a sizeable fraction sits at intermediate metallicity, but the summary calls it 'clear evidence.' That's overreach.\n\nReproducibility is partial—no code or parameters, and the templates come from earlier same-group papers. Not fatal, but for a catalog paper it would help to ship the fitting code.\n\nVerdict: worth serious refereeing. The catalog deserves to exist, and the systematic issues are correctable—calibrate zero-points per diagnostic on common stars, re-derive the accuracy claims, and tone down the Blazhko statement. If the authors do that, this becomes a standard reference. As is, I'd be careful using the V_gamma values for stars with only Balmer measurements.","headline":"Large, useful RV catalog for RR Lyrae, but the paper's own Fig. 5 shows cross-diagnostic V_gamma offsets (Hα–Hβ ~7 km/s, Hα–Hγ ~14 km/s) that contradict the homogeneity claim.","tokens_in":31918,"tokens_out":4144,"would_cite":false,"duration_ms":37514,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"By fitting universal radial-velocity-curve templates to sparse spectra, the paper builds a homogeneous catalog of barycentric velocities and pulsation amplitudes for 15,955 field RR Lyrae stars — and uses it to show that metallic-line RV am","keywords":["RR Lyrae stars","radial velocities","RV templates","barycentric velocity","Bailey diagram","Blazhko variables","metallicity trends","Galactic kinematics"],"falsifier":"Take the paper's calibration samples of fully phase-covered RR Lyrae curves and randomly thin each curve to 1, 3, and 7 phase points hundreds of times; if the template-recovered V_gamma shows a statistically significant offset that grows with iron abundance or with Blazhko phase, or if the free-amplitude RV amplitude systematically drifts lower for stars with periods above 0.6 days, the shape-universality premise fails and the catalog's metallicity and Bailey-diagram trends would require revision.","tokens_in":30721,"feed_emoji":"🌟","tokens_out":7750,"duration_ms":69766,"temperature":0.7,"pith_summary":"RR Lyrae stars are ancient, radially pulsating tracers of the Milky Way's oldest populations, but using them for Galactic kinematics requires knowing each star's barycentric velocity — the velocity of the stellar center of mass — rather than the pulsating surface velocity. This paper aims to establish that a single set of radial-velocity-curve templates, applied to as few as one to three randomly phased spectroscopic measurements, can recover that barycentric velocity and the pulsation velocity amplitude accurately enough to build the largest homogeneous catalog of such measurements for field RR Lyrae. The paper reports catalog-level accuracies of 3.8 km/s for well-sampled curves, 6.5 km/s for 3–7 phase points, and 11.3 km/s for fewer than three points. Using this catalog, the paper derives three physical results: the RV amplitude from metallic lines is essentially independent of iron abundance while higher Balmer lines show a clear decline with metallicity; the scaling between photometric and RV amplitudes is logarithmic for RRab and linear for RRc stars; and the spectroscopic Bailey diagram shows a smooth, low-scatter transition with candidate Blazhko variables skewed toward metal-rich, short-period populations. A sympathetic reader would care because, if the template universality holds, the catalog becomes the reference sample for RR Lyrae kinematics and a calibration set for next-generation spectroscopic surveys.","feed_headline":"RR Lyrae velocities pinned to 3.8 km/s across 15,955 stars","feed_subtitle":"One set of radial-velocity templates turns sparse spectra into a reference catalog for RR Lyrae kinematics.","key_machinery":"The central object is the radial-velocity-curve (RVC) template: a universal, phase-dependent shape for the pulsation velocity of RR Lyrae stars, separately constructed for metallic lines (Fe, Mg, Na) and Balmer lines (Hα, Hβ, Hγ, Hδ). Two fitting modes are used: a free-amplitude mode, which derives the RV amplitude directly from four or more measurements, and a fixed-amplitude mode, which scales the template using the photometric V-band amplitude for three-phase-point data. The templates convert sparse, randomly phased RV measurements into a barycentric velocity (relative to the Sun) and a pulsation amplitude, and are validated against two calibration samples with full phase coverage.","core_discovery":"The central discovery is that template fitting of radial-velocity curves, based on universal shapes for metallic and Balmer lines, yields V_gamma and RV amplitudes that agree with full-phase measurements to within a few km/s even when the target has only three to seven phase points. Applied to 17,563 field RR Lyrae (12,353 RRab, 5,011 RRc, 199 RRd), the method produces a catalog of 15,955 stars after quality filtering with typical V_gamma accuracies of 3.8, 6.5, and 11.3 km/s depending on phase coverage. From this homogeneous sample, the paper reports three physical results: (i) metallic-line RV amplitudes are essentially independent of metallicity (slope below the 2σ level), whereas higher-","pith_inferences":["If the template shapes are truly universal across period, amplitude, and metallicity, the same pipeline could be applied to sparsely sampled RR Lyrae in external galaxies or to massive upcoming surveys with limited epochs, potentially multiplying the kinematic sample size by an order of magnitude without new dense-phase observations.","The contrast between the metallicity-independent metallic-line amplitudes and the metallicity-dependent Balmer-line amplitudes suggests that the atmospheric layers probed by Hγ and Hδ shift in optical depth as a function of iron abundance in a way that directly affects the measured pulsation velocity; a testable prediction is that this slope should correlate with effective temperature across the i","The break in the RRab amplitude-scaling relation implies a physical turnover in the luminosity-to-radius amplitude ratio; an independent check would be to measure RV amplitudes for a sample of long-period, low-amplitude RRab stars and verify they follow the logarithmic extrapolation rather than the linear relation used in earlier work.","The smooth, low-scatter spectroscopic Bailey diagram could be exploited as a purely kinematic classifier of RR Lyrae subtypes, potentially improving mode identification for stars where photometric classification is ambiguous."],"forward_implications":["The catalog becomes the reference sample of barycentric radial velocities and RV amplitudes for field RR Lyrae, allowing future spectroscopic surveys to calibrate and validate their own RV measurements without re-observing stars across a full pulsation cycle.","The nonlinear (logarithmic) scaling relations for RRab stars mean RV amplitudes can be predicted from photometric amplitudes alone, extending kinematic and pulsation studies to the much larger photometric samples of RR Lyrae.","The near-zero metallicity dependence of metallic-line RV amplitudes implies that the spectroscopic Bailey diagram is largely free of metallicity bias, simplifying the interpretation of the short- and long-period sequences and of the spread at fixed period.","The increasing metallicity sensitivity of the higher Balmer lines provides a new spectroscopic diagnostic of atmospheric structure, since those lines form in higher atmospheric layers than the metallic lines.","The finding that candidate Blazhko variables are skewed toward metal-intermediate and metal-rich, short-period populations, if confirmed, indicates that metal-rich environments favor Blazhko modulation and connects the phenomenon to the high-amplitude short-period (HASP) regime.","The reduced scatter in the spectroscopic Bailey diagram suggests that period–RV-amplitude combinations can be used as a distance- and reddening-independent probe of chemical composition, complementing direct iron-abundance measurements."],"fun_headline_variants":["Template RV curves yield 16k-star RR Lyrae catalog at 3.8 km/s","3.8 km/s radial velocities for 15,955 RR Lyrae from sparse spectra","RR Lyrae RV catalog: template fitting tames sparse phase coverage","Convection clues from RV-photometry scaling in 16k RR Lyrae"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the radial-velocity-curve templates have a universal shape across the full range of period, amplitude, metallicity, and Blazhko phase, so that even a single phase point yields an unbiased barycentric velocity; if template shapes vary with metallicity or are distorted for Blazhko or metal-rich stars, the V_gamma and amplitude values inherit systematic offsets that propagate into every scaling relation and Bailey-diagram trend in the paper.","fun_headline_variants_meta":{"raw":{"variants":["Template RV curves yield 16k-star RR Lyrae catalog at 3.8 km/s","3.8 km/s radial velocities for 15,955 RR Lyrae from sparse spectra","RR Lyrae RV catalog: template fitting tames sparse phase coverage","Convection clues from RV-photometry scaling in 16k RR Lyrae"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00023,"raw_usage":{"total_tokens":1411,"prompt_tokens":926,"completion_tokens":485,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":670,"completion_tokens_details":{"reasoning_tokens":408}},"tokens_in":670,"tokens_out":485,"duration_ms":5233,"temperature":1.0,"reasoning_tokens":408,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T03:50:16.768742+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the paper's calibration samples of fully phase-covered RR Lyrae curves and randomly thin each curve to 1, 3, and 7 phase points hundreds of times; if the template-recovered V_gamma shows a statistically significant offset that grows with iron abundance or with Blazhko phase, or if the free-amplitude RV amplitude systematically drifts lower for stars with periods above 0.6 days, the shape-universality premise fails and the catalog's metallicity and Bailey-diagram trends would require revision.","supporting_citations":[],"review_version":1}