{"id":"c0350ff4-8784-4eb5-ac56-8c787623a101","arxiv_id":"2608.10909","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Using pulsation period instead of radius as a constraint for Cepheid evolutionary models systematically lowers the predicted radius, exposing a tension that is partially explained by a nonlinear radius increase in full-amplitude pulsation models.","lead":"The paper tests whether pulsation periods or radii of Cepheids in eclipsing binaries give consistent evolutionary model fits, and finds they do not. Period-selected models predict radii several sigma below observations for the most precise stars, and part of the offset is blamed on a nonlinear radius increase in large-amplitude pulsators.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Track-dependent systematic error in RSP periods is not ruled out by the sparse gyre comparison; a varying period offset could rotate the period-matched locus in Fig. 5 and inflate the reported radius tension.","rationale":"The paper is a careful, transparent modeling study, and the central claim is supported by a consistent pattern across several Cepheids, with robustness checks against a constant period shift and enlarged period uncertainty. However, the central claim depends on the computed pulsation periods being accurate across the relevant part of the HR diagram, and the existing validation does not establish that the period error is track-independent. The reader's weakest assumption already identified period accuracy as the load-bearing premise; this stress test sharpens that concern into a specific, testable track-dependence issue. The constant-period and constant-radius loci in Fig. 5 are nearly parallel, so the reported sigma_R offsets are essentially the separation between two lines in the HR diagram. A period error that varies with Teff or with pulsation amplitude along the blue loop would rotate one line relative to the other, moving the separation and changing the quoted significance. The gyre comparison in Sect. 2.4 samples only a few selected points and yields a constant offset, so it cannot detect such rotation. Appendix C applies a constant shift, which is the correct test for a constant systematic but not for a track-dependent one. The proposed dense RSP–gyre comparison along the actual tracks used in Fig. 5 would directly quantify any track dependence. Even if such an error exists, the enormous tension for CEP-0227 (12 sigma_R before and 6 sigma_R after the nonlinear radius correction and a 1.4% shift) suggests the conclusion would likely survive, but the significance could change substantially for the less extreme systems such as CEP-4506. Because the reader's verdict is already CONDITIONAL and hinges on the same premise, this concern does not move the verdict; it provides a concrete way to settle whether the concern actually lands. Until the test is run, the conditional acceptance is the appropriate stance.","tokens_in":52,"tokens_out":10983,"duration_ms":173254,"concrete_test":"Using the same final-grid tracks as in Fig. 5 for CEP-0227 and CEP-4506, compute linear periods with both RSP and gyre at a dense sequence of points spanning the full width of the instability strip (e.g., every 0.005 in log Teff from blue to red edge). Measure the fractional difference (P_RSP − P_gyre)/P_RSP as a function of Teff at fixed M, Z, and f_H. Also run nonlinear RSP models at the blue edge, center, and red edge of each track to measure how δP and the nonlinear radius shift vary along the track. If the RSP–gyre difference varies by more than ~1% across the strip, or if δP changes by more than ~0.3% or the radius shift by more than ~0.1 R_sun between the reference point and the edges, then the constant-correction assumptions in Sects. 2.4 and 3.2 break down and the period-matched locus used to quote the tension must be recomputed with track-dependent corrections.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim in Sect. 3.2 that period- and radius-selected models are mutually inconsistent rests on the accuracy of RSP linear periods plus a per-star constant nonlinear correction (Sect. 2.4). The paper assumes 1% theoretical period uncertainty and validates against gyre at 'a few evolutionary tracks' (Sect. 2.4), finding periods shorter by 0.7–1.4%. Appendix C then applies a constant 1.4% shift and shows the tension persists. However, the constant-period and constant-radius loci in Fig. 5 are nearly parallel, and the headline numbers (e.g., −12.0σ_R for CEP-0227) are the offset between them. A period error that varies with effective temperature, or with position along the blue loop, would rotate the period locus relative to the radius locus, potentially reducing or eliminating the offset. The gyre comparison is too sparse to characterize such track dependence, and the grid-variation tests in Fig. 1 are single-point sensitivity checks, not accuracy tests across the instability strip. The nonlinear period correction is likewise computed only at the reference model and assumed slowly varying, without verification along the track. None of the robustness checks covers a non-constant, Teff-dependent period error, so the quoted tension could be in part an artifact of the period computation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper models four LMC eclipsing binary systems containing classical Cepheids (CEP-0227, CEP-4506, CEP-2532, CEP-1718) and uses CEP-1812 as an additional single-star probe, combining MESA evolutionary tracks with RSP linear pulsation periods and per-star nonlinear period corrections. The central exercise is a chi-square comparison of two families of solutions: those selected using {Teff, log L, R} and those selected using {Teff, log L, P}. The authors find that the two families are mutually inconsistent: models that match the observed pulsation period predict radii that are too small by 3 to 12 sigma_R for the F-mode Cepheids, while models that match the observed radius predict periods that are too long. Part of the offset is attributed to a nonlinear increase in the mean radius of full-amplitude pulsation models, and the paper recommends using the pulsation period rather than the radius as the primary constraint for future modeling.","tokens_in":28355,"tokens_out":10584,"duration_ms":104757,"significance":"If correct, the result is significant for Cepheid modeling: it would imply that static evolutionary radii and RSP-based pulsation periods cannot be simultaneously reconciled for the most precisely measured Cepheids, with consequences for period-radius relations, mass-discrepancy studies, and the use of eclipsing binaries as calibrators. The paper has clear strengths: it uses homogeneous P18 data, presents a detailed sensitivity study of RSP periods (Fig. 1), explicitly tests robustness against a 1.5% period uncertainty and a constant 1.4% period shift (Appendix C), and transparently discusses mass-loss treatment, metallicity degeneracies, and radius definition issues. The period comparison is not circular: periods come from RSP calculations at fixed model parameters, and the nonlinear correction is computed at the observed central parameters rather than tuned to remove the tension. The main weakness is that the headline sigma values assume the RSP periods are accurate at the 1% level and that any remaining period error is approximately constant along the relevant tracks; neither assumption is fully demonstrated.","major_comments":[{"comment":"The headline offsets (e.g., -12.0 sigma_R for CEP-0227) are essentially the separation between the period-matched and radius-matched loci, which the paper states are nearly parallel. A period error that varies with effective temperature along the blue loop would rotate the period locus relative to the radius locus and directly change these offsets. The validation in Sect. 2.4 consists of grid-parameter variations up to 0.6%, a gyre comparison at 'a few evolutionary tracks' yielding periods shorter by 0.7-1.4%, and an Appendix C test that applies a constant 1.4% shift with 1.5% uncertainty. None of these tests constrains a Teff-dependent or crossing-dependent period error. Please either compute RSP and gyre periods along the same tracks at several points spanning the observed error boxes, or provide a conservative model for a Teff-dependent period error and show that the tension survives. Without this, the quantitative tension could be partly an artifact of the period computation.","section":"Sect. 2.4, Fig. 5, Appendix C"},{"comment":"The reported significance is expressed in units of the observational radius uncertainty sigma_R only. The assumed 1% theoretical period uncertainty maps into an uncertainty in the model radius at the period-matched location; through the period-mean-density relation this is roughly (2/3)x1% ~ 0.7% in radius, which is larger than sigma_R for CEP-0227 (0.34%) and comparable to sigma_R for CEP-4506 (0.7%). The quoted deviations (-12.0, -3.2, -6.2 sigma_R for the F-mode stars) therefore overstate the confidence with which the period-selected radius is excluded unless the period-induced radius uncertainty is propagated into the comparison. Appendix C enlarges the period uncertainty to 1.5% and applies a constant 1.4% shift, but it still reports residuals only in units of sigma_R and does not combine the period-induced radius error with the observational radius error. Please report the residuals together with their combined uncertainties.","section":"Sect. 3.2, Fig. 5"},{"comment":"The attribution of part of the tension to a nonlinear radius increase is provisional in a way that matters for the quantitative message. The nonlinear radius correction is computed at one reference model with the observed central parameters and is then applied as a constant along the track, without a direct check that it is slowly varying along the relevant tracks. The correction itself varies by almost a factor of two among convective parameter sets for CEP-4506 (0.8-1.5%), and the paper notes that Farag et al. (2026) suggest sensitivity to the numerical solver. Since the post-correction tension for CEP-4506 is only 1.4 sigma_R, the statement that a nonlinear radius increase explains part of the tension should be presented as an estimate with an unquantified systematic error rather than as a firm accounting. For CEP-0227 the residual remains large (8.8 sigma_R after correction), so the main inconsistency does not depend on this correction, but the two situations should be clearly separated.","section":"Sect. 3.2, Table 1"}],"minor_comments":[{"comment":"The text says 'we modeled five classical Cepheids in DLEBs', while the paper actually analyzes six Cepheids in five systems (with CEP-1812 used only as a single-star probe). Please state the counts more precisely to avoid confusion.","section":"Sect. 4"},{"comment":"The text describes the 1.5% period uncertainty as 'nearly three times' the 1% value estimated in Sect. 2.4; 1.5% is 1.5 times 1%, not three times. Please correct this.","section":"Appendix C"},{"comment":"In the paragraph discussing CEP-0227, the star is referred to as 'CEP-227'; use the consistent identifier CEP-0227.","section":"Sect. 3.3"},{"comment":"The mass ratio q for CEP-1718 is listed only in the row of CEP-1718B; it would be clearer to state explicitly that q refers to the binary mass ratio of the system and to give it once for the system or in both rows.","section":"Table A.1"},{"comment":"The gyre comparison is described only in prose; a small table or figure listing the tracks and the RSP/gyre periods would make the 0.7-1.4% difference auditable and would also help the reader assess whether the difference depends on Teff.","section":"Sect. 2.4"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of A&A and the sample, although small, is the best currently available for this test. The main revision should focus on the period-systematics issue: the qualitative inconsistency is likely robust, but the headline sigma values need proper propagation of the theoretical period uncertainty and a check for Teff-dependent period errors. I do not see concerns about novelty or citation behavior."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is the first systematic demonstration I know of that using the pulsation period versus the radius to select evolutionary models for Cepheids in DLEBs gives mutually inconsistent answers, with period-selected models predicting radii too small by several sigma for the well-measured F-mode stars. The paper is transparent, the pattern across six stars is consistent, and the robustness test with 1.5% period error plus a 1.4% systematic shift reduces but does not remove the tension. That is a genuine result and a useful challenge to the community.\n\nWhat it does well: RSP periods computed along MESA tracks with nonlinear period corrections, honest discussion of radius definition (tau=2/3 vs full-amplitude Fourier mean), explicit flags that the nonlinear radius correction lacks a systematic study, and a clear recommendation to prefer the period as a primary constraint. It also shows the correction is mode-dependent and can be several sigma for F-mode stars. That is new and worth engaging.\n\nThe soft spots are where the reader and the stress test point. The headline sigmas propagate only the observational radius error. The 1% theoretical period uncertainty is an assumption; the gyre comparison covers 'a few' tracks and finds 0.7-1.4% shorter periods, but a constant shift is not the same as a track-dependent error. In Fig. 5 the constant-period and constant-radius loci are nearly parallel; a period error that varies with Teff or position along the blue loop could rotate the period locus and change or even remove the offset. The paper does not test this. The nonlinear radius correction is computed at the central parameters and applied along the track as a constant, with no verification that it is slowly varying. That is a second assumption that feeds the corrected sigmas. These are not fatal flaws, but they mean the result is conditional, not established at the quoted significance.\n\nFor CEP-0227 the tension remains near 6 sigma after the nonlinear correction, so even under the paper's own assumptions something else is needed. The paper says this, which is credit to it.\n\nThis paper deserves a serious referee. The right referee will push for a propagation of theoretical uncertainties into the headline sigmas, and for a denser gyre/alternative-code comparison across the strip. No code or data are shipped, which limits reproducibility and should be noted in review.\n\nI would bring it to a reading group and would cite it if I worked on Cepheid modeling or DLEB constraints.","headline":"A careful, transparent modeling study showing a systematic period-radius tension in Cepheid DLEBs; the pattern is real, but the headline sigmas are softer than advertised because theoretical period and nonlinear radius uncertainties are not propagated into them.","tokens_in":29022,"tokens_out":2028,"would_cite":true,"duration_ms":22566,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["85A15"],"pacs":["97.30.Gj","97.80.Hn"],"model":"deepseek-v4-flash","headline":"Cepheid radii and pulsation periods point to different stars","keywords":["classical Cepheids","eclipsing binaries","pulsation periods","period-mean density relation","stellar radius","nonlinear pulsation models","MESA RSP","Large Magellanic Cloud"],"falsifier":"If a single fundamental-mode Cepheid in an eclipsing binary could be fitted by one MESA+gyre (or MESA+RSP) model that simultaneously matches the observed period and the observed radius within 1σ, with a fully self-consistent nonlinear pulsation-on-evolution calculation including the radius definition used by the eclipse solution, the claimed incompatibility would be refuted; conversely, extending the sample to more F-mode systems with sub-percent radii and confirming the same 3–12σ_R systematic offset would strengthen it.","tokens_in":27901,"feed_emoji":"⭐","tokens_out":3760,"duration_ms":29989,"temperature":0.7,"pith_summary":"Eclipsing binaries give the most precise masses and radii we have for classical Cepheids, and those radii are usually the anchor for testing stellar evolution models. This paper asks what happens if you instead anchor the matching to the pulsation period, the single most precisely measured Cepheid observable. The answer is that the two anchors disagree: models that match the observed period systematically predict radii too small by 3 to 12σ_R for fundamental-mode Cepheids, while models that match the observed radius predict periods too long. The paper finds that part of the offset comes from a neglected nonlinear effect, the small but measureable increase in the mean radius of a star pulsating at full amplitude, and concludes that the radius should not be used as the primary constraint until that effect is systematically mapped.","feed_headline":"Cepheid radii and periods point to different stars","feed_subtitle":"Matching pulsation periods predicts radii too small by up to 12σ; the paper says trust the period, not the radius.","key_machinery":"The central object is the period–mean density relation, $P \\propto \\bar{\\rho}^{-1/2}$, which links the pulsation period to the star's radius and mean density. The paper exploits the fact that the period and the radius encode overlapping but observationally independent information, and compares the constant-period locus and the constant-radius locus along MESA evolutionary tracks in the HR diagram; the gap between these nearly parallel loci is the measured tension. A second load-bearing piece of machinery is the RSP nonlinear pulsation model, whose full-amplitude solution supplies both a nonlinear period correction (≤0.28% in this sample) and a nonlinear radius correction (the zeroth-order Fourier term of the time-dependent radius, up to ~1.1% larger than the static radius), which is applied when comparing model radii to observed radii.","core_discovery":"For classical Cepheids in detached eclipsing binaries, the evolutionary model that reproduces the observed pulsation period does not reproduce the observed radius, and the model that reproduces the observed radius does not reproduce the period. Using χ² matching on grids of MESA evolutionary tracks supplemented with RSP pulsation periods, the paper finds that period-selected solutions place the star at radii 3 to 12σ_R too small for the three fundamental-mode Cepheids (CEP-0227, CEP-4506, CEP-1812), while radius-selected solutions give periods too long by roughly 3 to 8σ under a 1% theoretical period uncertainty. The two constraints trace nearly parallel lines in the HR diagram, so the disagreement cannot be erased by adjusting mass, metallicity, or overshooting. The paper traces part of the effect to a nonlinear radius increase in full-amplitude pulsators: the RSP full-amplitude models have mean radii larger than the static evolutionary radii, up to 1.1% for CEP-0227 (3.2σ_R), and applying this correction reduces but does not eliminate the tension. It recommends treating the pulsation period as the primary constraint and calls for a systematic study of nonlinear radius corrections.","pith_inferences":["A testable corollary not pursued in the paper: Cepheids observed at higher pulsation amplitude should show a proportionally larger nonlinear radius inflation relative to their static models, so a sample spanning a range of amplitudes could directly verify the proposed mechanism.","The paper's nonlinear radius correction suggests that radius determinations of any radially pulsating star from eclipse light curves are subtly biased if static model radii are used as priors; the same correction logic could apply to RR Lyrae stars in eclipsing binaries.","If the tension persists after a systematic nonlinear-radius study, the remaining offset in CEP-0227 would point to an independent physics shortfall, possibly rotationally induced radius inflation, which the paper flags but does not model.","The degeneracy the paper highlights between lower metallicity and stronger core overshooting means that spectroscopic metallicities for these systems would do more to settle the radius question than additional period measurements."],"forward_implications":["If the paper is right, radii of Cepheids in eclipsing binaries should no longer be used as the primary matching constraint in evolutionary-model fits; the pulsation period, with its orders-of-magnitude better observational precision, should be preferred.","Period–radius relations built from static evolutionary models carry a small but real systematic offset relative to observed full-amplitude pulsators, and should be corrected for the nonlinear radius shift before comparisons at sub-percent radius precision.","The known Cepheid mass-discrepancy problem is entangled with a radius discrepancy: part of what looks like a mass problem in eclipsing-binary Cepheids may instead be a mismatch between static model radii and the time-averaged radii of pulsating stars.","First-overtone Cepheids should show the same period–radius tension, but current radius measurements are about ten times less precise; the prediction is that better radii for 1O Cepheids will reveal the same systematic offset.","Systematic uncertainties in computed pulsation periods (up to ~1.4% from RSP–gyre comparison, or a 1.5% period uncertainty) reduce the tension by only a few σ_R, so the discrepancy is robust against period-scale systematics."],"supporting_citations":[{"why":"Supplies the observed physical parameters of the five systems — masses, radii, effective temperatures, luminosities, periods — that are matched by the grids.","marker":"P18"},{"why":"Documents the MESA RSP module used to compute both linear and nonlinear pulsation periods and the nonlinear radius corrections.","marker":"Paxton et al. 2019"},{"why":"Defines the reference MESA evolutionary setup (mixing, opacities, atmosphere, MLT calibration) on which the model grids rest.","marker":"Ziółkowska et al. 2024"},{"why":"Establishes the dependence of blue-loop luminosity and instability-strip crossing on metallicity and overshooting, motivating the grid ranges and the fixed envelope overshooting.","marker":"Smolec et al. 2026b"},{"why":"Supplies the gyre pulsation code used as an independent check on the linear periods and the 0.7–1.4% period offset used in the robustness test.","marker":"Townsend & Teitler 2013"},{"why":"Provides the amplitude-equation basis for the prediction that full-amplitude pulsation modifies the mean radius, referenced as prior theoretical support for the nonlinear radius correction.","marker":"Buchler & Goupil 1984"},{"why":"Introduces a self-consistent TDC pulsation calculation in MESA, cited as evidence that nonlinear radius corrections may depend on the numerical solver and as the path to future study.","marker":"Farag et al. 2026"}],"fun_headline_variants":["Cepheid radii and periods: models can't agree","Period-based Cepheid models fail radius test","Trust the period, not the radius, for Cepheids","Nonlinear pulsation skews Cepheid radius fits"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central claim rests on the assumption that the RSP-computed pulsation periods, after a small nonlinear correction and with a 1% uncertainty, are accurate enough across the relevant part of the HR diagram to be compared with radii at the few-sigma level; if the period calculation carries a systematic error tied to the evolutionary track, the tension could be an artifact.","fun_headline_variants_meta":{"raw":{"variants":["Cepheid radii and periods: models can't agree","Period-based Cepheid models fail radius test","Trust the period, not the radius, for Cepheids","Nonlinear pulsation skews Cepheid radius fits"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000188,"raw_usage":{"total_tokens":1431,"prompt_tokens":1142,"completion_tokens":289,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":758,"completion_tokens_details":{"reasoning_tokens":223}},"tokens_in":758,"tokens_out":289,"duration_ms":3076,"temperature":1.0,"reasoning_tokens":223,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:31:21.834494+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"If a single fundamental-mode Cepheid in an eclipsing binary could be fitted by one MESA+gyre (or MESA+RSP) model that simultaneously matches the observed period and the observed radius within 1σ, with a fully self-consistent nonlinear pulsation-on-evolution calculation including the radius definition used by the eclipse solution, the claimed incompatibility would be refuted; conversely, extending the sample to more F-mode systems with sub-percent radii and confirming the same 3–12σ_R systematic offset would strengthen it.","supporting_citations":[],"review_version":1}