{"id":"59f27084-47fc-43b1-a6d3-e6080708270d","arxiv_id":"2608.07792","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"δ Scuti pulsator occurrence drops from 88% in clusters younger than 200 Myr to 62% in older clusters, with surviving pulsators rotating faster, evidence that age and rotation control pulsation.","lead":"Using TESS space telescope data, the authors measured how often A- and F-type stars in 20 nearby open clusters show pulsations, correcting for detection limits. They found that pulsation is more common in clusters younger than 200 million years and that the pulsators that survive in older clusters are spinning faster.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Occurrence correction assumes the Kepler amplitude distribution shape; if cluster amplitude distributions differ with age, the 88% vs 62% gap could be an artifact.","rationale":"The reader's verdict is CONDITIONAL and identifies the same weakest assumption: the amplitude-distribution shape in the §3.4 occurrence correction. I agree with that identification. My stress-test focuses on why this is the most load-bearing point: the central quantitative claim (the 26-point drop in occurrence) is mediated by the correction, while the uncorrected fraction contrast is smaller and noisier. The paper already shows sensitivity to membership catalogs (§3.4, Figure 5) and discusses binary effects (§3.5), but it does not test the amplitude-shape assumption with cluster data. The remedy is straightforward and computational: recompute occurrences under alternative amplitude distributions and check the young-old gap. This does not require rejecting the paper; the conditional acceptance remains appropriate. The rotation claim (Figures 9 and 10) is secondary and I do not object to it as strongly, since vbroad caveats are acknowledged and the trend is independently visible in fractions. No internal inconsistency or fraudulent practice is present; the concern is model dependence in a correction factor. Therefore the verdict stays CONDITIONAL with a specific additional test required.","tokens_in":18559,"tokens_out":1552,"duration_ms":13246,"concrete_test":"Re-run the occurrence calculation while replacing the fixed Kepler amplitude distribution in §3.4 with two alternative distributions: (a) a log-normal fit to each cluster's own detected amplitudes (allowing shape, not just median, to vary), and (b) the Kepler distribution broadened/narrowed by factors of 2 in log-amplitude. If the young-old occurrence gap (88±3% vs 62±3%) changes by more than the quoted ±3% uncertainties under (a) or (b), the age trend is not robust to the amplitude-shape assumption.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim that occurrence drops from 88±3% to 62±3% across the 200 Myr split rests on the §3.4 completeness correction. That correction maps an observed pulsator fraction to a true occurrence by assuming that intrinsic δ Scuti amplitudes in every cluster follow the Kepler field-star distribution of Murphy et al. (2019), shifted uniformly in log-amplitude to match each cluster's observed median amplitude. The recovery-rate grid (§3.3) is amplitude- and magnitude-dependent, so if young clusters have intrinsically larger or broader amplitude distributions—physically plausible given stronger driving, different mode excitation, or rotation—the correction will differentially over- or under-correct old and young clusters. The paper computes occurrence as a function of median amplitude (panel f of Figure 4), but no test checks whether the underlying shape, not just the median, varies with age. Notably, Table 1 shows young clusters like BH 99 and NGC 6405 hit occurrence 100 +0/−14/−15, while the old-cluster average is driven down by Stock 2 (41±8%) and Mamajek 4 (49±15%); a modest shape-change bias could plausibly close part of this gap. The §3.4 method is internally consistent, but the assumption that a single Kepler amplitude-shape applies at all ages is external and untested. This is the load-bearing point because the pulsator fractions alone (74±4% vs 52±4%) show a weaker contrast, and the occurrence correction amplifies it.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses TESS photometry to search for δ Scuti pulsators in 20 open clusters within 500 pc and with ages between roughly 20 and 900 Myr, identifying 487 pulsators in total. The authors define a pulsator occurrence that corrects the observed pulsator fraction for incompleteness using injection-recovery tests and an assumed Kepler amplitude distribution, and they report average occurrences of 88±3% for clusters younger than 200 Myr and 62±3% for older clusters. They also report that pulsators in older clusters rotate faster on average and interpret both trends as evidence that helium settling suppresses pulsations over time unless counteracted by rapid rotation.","tokens_in":18824,"tokens_out":7397,"duration_ms":69878,"significance":"If the age trend survives scrutiny, this is an important result: it turns previously anecdotal cluster-to-cluster differences (Pleiades, NGC 2516, Cep-Her, NGC 3532) into a systematic age sequence and directly constrains models of helium diffusion and rotational mixing in A/F stars. The paper's strengths include a homogeneous TESS analysis of 20 clusters, explicit injection-recovery simulations for each cluster, a membership-catalog robustness check, a treatment of equal-mass binaries, and a machine-readable star table. The central weakness is that the headline occurrence numbers inherit an untested assumption about the universal shape of the pulsation amplitude distribution, together with an inconsistency between the detection threshold used for real stars and the threshold used in the recovery simulations.","major_comments":[{"comment":"The occurrence correction assumes that the intrinsic δ Scuti amplitude distribution of every cluster has the same shape as the Kepler field-star distribution of Murphy et al. (2019), with only its median shifted in log-amplitude. This is an external, untested input. Because the correction is applied per cluster, and the young and old subsamples in Table 1 have very different occurrence values (e.g., BH 99 and NGC 6405 at 100+0−14/−15 versus Stock 2 at 41±8 and Mamajek 4 at 49±15), a modest age dependence in the width or high-amplitude tail of the amplitude distribution would translate into a bias of several percentage points, which is the same order as the reported 88% vs 62% difference. The authors should test this assumption with cluster-internal data, for example by comparing the observed amplitude distributions of detected pulsators in young and old clusters, or by recomputing occurrence under alternative plausible distribution shapes and showing that the age contrast is robust.","section":"Sec. 3.4"},{"comment":"Pulsators are identified in §3.1 using a variable frequency boundary, log skewness ≥0.4, and visual inspection, but the recovery criterion in §3.3 is log skewness ≥0.75. The injection-recovery completeness is therefore measured for a stricter detection rule than the one used to build the real sample. This makes the inferred completeness too low and the occurrence too high. Because the correction depends on apparent magnitude (Figure 7), and the young and old subsamples have different distance and magnitude distributions, the bias need not cancel in the 88% vs 62% comparison. The authors should rerun the recovery tests with the actual classification protocol, including the visual-confirmation step, or quantify how the occurrence-versus-age result changes when the recovery threshold is varied.","section":"Secs. 3.1 and 3.3"},{"comment":"The 200 Myr threshold is introduced as the basis for the central 88±3% vs 62±3% averages, but it appears to be selected after inspecting the data, and NGC 6475, with an age of 200±50 Myr, is placed on the older side. The reported significance therefore depends on a data-informed boundary and on the averaging scheme. The authors should report how the averages and their difference change when the threshold is varied over the cluster-age uncertainties (for example 150-300 Myr), when NGC 6475 is moved to the young group, and when the average is computed with inverse-variance weights rather than as a simple unweighted mean.","section":"Sec. 4.1"}],"minor_comments":[{"comment":"In the first bullet of the conclusions, 'which an average occurrence of 88±3%' should read 'with an average occurrence of 88±3%'.","section":"Sec. 5"},{"comment":"The asymmetric errors in Table 1 (e.g., 100+0−14) should be accompanied by a sentence explaining how they are propagated into the quoted ±3% averages.","section":"Table 1"},{"comment":"The text should quantify the documented non-linearity of Gaia vbroad as a vsini estimator (underestimating slow rotators and overestimating the fastest rotators) and state how this calibration uncertainty affects the mean rotation comparison in Figure 9.","section":"Sec. 4.4"},{"comment":"The paper never explicitly defines that the 'pulsator occurrence' quoted for each cluster is the maximum of the occurrence-versus-T_eff curve rather than an occurrence integrated over the instability strip; this should be stated when the statistic is first used.","section":"Sec. 3.2"}],"recommendation":"major_revision","confidential_remarks":"The dataset and overall approach are publishable, but the central age comparison rests on an untested amplitude-distribution shape assumption and on an inconsistency between the real detection threshold and the injection-recovery threshold. These are addressable with additional analyses rather than being fundamental flaws in the sample or the method. I would be comfortable accepting after the robustness tests requested above are added."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Berry et al. report that δ Scuti pulsator occurrence drops from 88±3% in clusters younger than 200 Myr to 62±3% in older clusters, and that surviving pulsators rotate faster. This is a new result: nobody has measured occurrence across a 20-cluster age sequence. The sample construction is careful—TESS photometry, membership catalogs, injection/recovery corrections—and the authors are appropriately tentative about metallicity and about the disagreement with Wang et al. (2025). I see no fatal flaw in the equations.\n\nThe soft spots are real but not fatal. The completeness correction assumes every cluster has the same amplitude-distribution shape as the Kepler field sample, just shifted in log-amplitude. If young and old clusters have different shapes, the correction could bias the 88 vs 62 difference. The paper does not test this. However, the raw pulsator fractions already show a weaker but same-direction trend (74±4% vs 52±4%), so the age trend is not purely an artifact of the correction; the correction amplifies the contrast. The 200 Myr split is data-informed; NGC 6475 sits right on the boundary. And the rotation proxy (Gaia vbroad) has known biases, so the rotation maintenance story should be framed more cautiously.\n\nIf I were referee, I'd ask for a sensitivity test on the amplitude-shape assumption—e.g., using a broader or skewed distribution, or checking whether the inferred occurrence changes across age bins. The authors should also propagate vbroad uncertainties or rely on independent vsini. These are fixable. The paper is worth a serious referee, and the age trend is worth citing.","headline":"A solid new age trend for δ Scuti occurrence, but the completeness correction's amplitude-shape assumption is the load-bearing soft spot.","tokens_in":19393,"tokens_out":3395,"would_cite":true,"duration_ms":28856,"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":"δ Scuti pulsators become rarer as open clusters age, with occurrence dropping from 88% in clusters younger than 200 Myr to 62% in older clusters.","keywords":["delta Scuti stars","pulsator occurrence","open clusters","stellar pulsation","stellar rotation","helium diffusion","instability strip","TESS photometry"],"falsifier":"Redo the occurrence correction for each cluster using the amplitude distribution measured from that cluster's own detected pulsators instead of the Kepler field distribution; if the average occurrence gap between young and old clusters shrinks below statistical significance, then the stated age trend depends on the assumed amplitude distribution rather than on stellar physics.","tokens_in":18342,"feed_emoji":"🌟","tokens_out":3489,"duration_ms":32549,"temperature":0.7,"pith_summary":"This paper uses TESS photometry to identify 487 δ Scuti pulsators in 20 nearby open clusters spanning ages from about 20 to 900 Myr, then measures pulsator occurrence — the fraction of instability-strip stars that actually pulsate, corrected for detection incompleteness. The central claim is that pulsator occurrence declines with cluster age: clusters younger than 200 Myr average 88±3% occurrence, while older clusters average 62±3%, a statistically significant difference. The paper also finds that pulsators in older clusters rotate faster on average than those in younger clusters, and it interprets this as evidence that rapid rotation helps maintain δ Scuti pulsations by counteracting helium settling out of the ionization zone. If correct, the age trend explains why field-star samples show intermediate pulsator fractions and identifies rotation, not just position in the instability strip, as a key factor controlling whether A/F stars pulsate. The paper is careful to separate the raw observed fraction from the completeness-corrected occurrence, arguing that the age dependence is only clearly visible after this correction.","feed_headline":"δ Scuti pulsation drops from 88% to 62% as clusters age","feed_subtitle":"TESS survey of 20 open clusters links pulsation lifetime to age and rapid rotation.","key_machinery":"The load-bearing tool is the pulsator occurrence calculation: injection-and-recovery tests on non-detection light curves, combined with the field-star δ Scuti amplitude distribution from Kepler, shifted in log-amplitude to match each cluster's observed median amplitude, yield a map from observed pulsator fraction and median amplitude to a completeness-corrected occurrence. Pulsator fraction and median amplitude are computed as continuous functions of $T_{\\rm eff}$ using a Gaussian kernel whose width is optimized per cluster by maximizing a leave-one-out Bernoulli log-likelihood. This machinery converts the raw fraction into an occurrence that can be compared fairly across clusters of different distances and TESS coverage.","core_discovery":"The paper establishes that δ Scuti pulsator occurrence decreases with age across coeval stellar populations. From 20 open clusters within 500 pc, the average occurrence in clusters younger than 200 Myr is 88±3%, while clusters older than 200 Myr average 62±3%; the paper states that this difference is statistically significant and shows that the pulsator occurrence decreases with age. It further reports that pulsators in older clusters rotate more rapidly than their younger counterparts and that hotter pulsators (≳8500 K) may stop pulsating earlier than cooler ones. The physical interpretation is that gravitational settling depletes helium from the near-surface ionization zone that drives the κ-mechanism, while rapid rotation mixes helium back into that zone, so the surviving pulsators in old populations are preferentially fast rotators.","pith_inferences":["An untested implication is that angular momentum history, not age alone, sets the observable pulsator fraction; stellar models coupling rotation and helium diffusion should predict a two-dimensional occurrence surface in age and rotation that this sample could directly test.","A testable extension is to repeat the occurrence calculation using amplitude distributions measured from each cluster's own detected pulsators rather than the Kepler field distribution; if the young-versus-old gap narrows or vanishes, the age trend is sensitive to that assumption.","The slow-rotating, metal-rich Am stars in the instability strip may account for a large share of non-pulsators; estimating the Am fraction per cluster could separate chemical-composition effects from age effects in the occurrence decline.","The same occurrence machinery could be applied to more distant clusters observed with PLATO, extending the age baseline beyond 1 Gyr and testing whether the occurrence continues to fall or plateaus."],"forward_implications":["If the central claim is right, the raw pulsator fraction underestimates how many A/F stars pulsate in young clusters, and completeness corrections are necessary before comparing populations.","The age trend implies that many δ Scuti stars stop pulsating on main-sequence timescales, so the instability strip is not a static boundary for pulsation presence.","Rapid rotation becomes a longevity factor: stars that remain pulsating in older clusters are preferentially fast rotators, linking pulsation to angular momentum evolution.","Hotter δ Scuti stars turning off earlier means effective temperature, age, and rotation must be considered together when modeling pulsator populations.","The observed decline around 200 Myr is consistent with models predicting roughly 50% helium depletion from the ionization zone by 100 Myr, supporting helium diffusion as the driving mechanism."],"supporting_citations":[{"why":"Supplies the Kepler field-star δ Scuti amplitude distribution used for the completeness correction and the empirical instability-strip bounds, as well as the field pulsator fraction baseline.","marker":"Murphy et al. (2019)"},{"why":"Introduced the pulsator occurrence concept with injection-and-recovery tests and the amplitude-distribution shift; this paper extends that method to 20 clusters.","marker":"Berry et al. (2025)"},{"why":"Provides the adopted δ Scuti classifications for the Pleiades and the young-cluster pulsator fraction of about 80% used as a comparison point.","marker":"Bedding et al. (2023)"},{"why":"Provides the adopted δ Scuti classifications for NGC 2516 and another young-cluster pulsator fraction of about 80%.","marker":"Li et al. (2024)"},{"why":"Supplies the Cep-Her complex occurrence near 100% and the Gaia vbroad rotation proxy calibration used for the rotation analysis.","marker":"Murphy et al. (2024)"},{"why":"Demonstrates the increase in pulsator fraction with rotation rate and defines the empirical red edge that the rotation comparison in this paper reproduces.","marker":"Gootkin et al. (2024)"},{"why":"Provides the Gaia DR3 open-cluster membership catalogs and cluster ages used to define the sample and its tidal-radius membership flag.","marker":"Hunt & Reffert (2023, 2024)"}],"fun_headline_variants":["δ Scuti pulsator occurrence drops from 88% to 62% with age","Age kills δ Scuti pulsations unless stars spin fast, TESS finds","Young clusters pulsate more: δ Scuti occurrence falls as they age","Rapid rotation preserves δ Scuti pulsations in older star clusters"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole age comparison rests on the assumption that the shape of the δ Scuti amplitude distribution is the same in every cluster as in the Kepler field-star sample, shifted only in median amplitude; if young and old clusters intrinsically differ in amplitude shape, the completeness corrections are biased differently by age and the 88% versus 62% gap could be an artifact.","fun_headline_variants_meta":{"raw":{"variants":["δ Scuti pulsator occurrence drops from 88% to 62% with age","Age kills δ Scuti pulsations unless stars spin fast, TESS finds","Young clusters pulsate more: δ Scuti occurrence falls as they age","Rapid rotation preserves δ Scuti pulsations in older star clusters"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000736,"raw_usage":{"total_tokens":3283,"prompt_tokens":931,"completion_tokens":2352,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":547,"completion_tokens_details":{"reasoning_tokens":2270}},"tokens_in":547,"tokens_out":2352,"duration_ms":14829,"temperature":1.0,"reasoning_tokens":2270,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T04:13:00.631247+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Redo the occurrence correction for each cluster using the amplitude distribution measured from that cluster's own detected pulsators instead of the Kepler field distribution; if the average occurrence gap between young and old clusters shrinks below statistical significance, then the stated age trend depends on the assumed amplitude distribution rather than on stellar physics.","supporting_citations":[{"cited_title":"2025, ApJ, 995, 128, doi: 10.3847/1538-4357/ae18c6","cited_arxiv_id":null,"evidence_quote":"Introduced the pulsator occurrence concept with injection-and-recovery tests and the amplitude-distribution shift; this paper extends that method to 20 clusters."}],"review_version":1}