{"id":"00d2d87b-e927-4892-bbec-f0c59ad85fed","arxiv_id":"1908.05222","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A Kepler-data catalog of surface rotation periods and photometric activity for about 15,900 M and K main-sequence stars, including about 4,431 new rotation periods and cross-checks against the McQuillan et al. catalog.","lead":"This paper measures rotation periods and a magnetic-activity proxy for about 15,900 M and K dwarf stars observed by the Kepler space telescope, using longer light curves than previous catalogs. It adds rotation periods for over 4,400 stars not listed in the earlier McQuillan et al. catalog and flags hundreds of likely non-single or misclassified targets.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The fast-rotator/high-Sph branch that drives the bimodality and the Sph-Prot anti-correlation may still contain unflagged binaries or pulsators; the paper's own cautions in §5.3 and §5.5 make this the load-bearing caveat.","rationale":"The reader's weakest assumption is the same one I would stress: the detected periodicities are assumed to be rotational spot modulation for targets not explicitly flagged as polluters, and the fast-rotator/high-Sph branch is the place where this assumption is least secure. The paper is transparent about this, and the independent McQuillan agreement is strong evidence that the pipeline measures real periodicities for the common sample. The catalog itself can therefore stand. The remaining risk is narrower but real: the two physical conclusions, the bimodal period distribution and the Sph-Prot anti-correlation, are weighted by the fast-rotator/high-Sph tail, and the authors explicitly state that this tail resembles binaries and classical pulsator candidates while failing their own flagging criteria. The proposed excision test uses only the published flags and thresholds, so it is directly feasible. It would not change the verdict from CONDITIONAL: the catalog is accepted with conditions, and the physical interpretation should remain contingent on this check. This is an internally acknowledged correctness risk, not a disagreement with consensus or a criticism of the authors' effort.","tokens_in":24789,"tokens_out":7821,"duration_ms":86182,"concrete_test":"Using the published Tables 3-5 with the Berger, Simonian, and FliPer flags, reconstruct the §5.3 K-dwarf Sph-Prot sample and remove every target with Prot < 10 days and Sph > 7000 ppm that is not already flagged as CP/CB, binary, or multiple-signal, i.e., the exact population §5.3 says may be polluted. Then recompute (a) the Spearman correlation between Sph and Prot and (b) the significance of the fast peak in the Prot histogram. If both survive the excision, the concern is settled; if the correlation weakens or the fast peak disappears, the physical claims are driven by the ambiguous tail. A sanity variant is to repeat with an Sph threshold of 10^4 ppm, the value §5.5 associates with synchronized binaries.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The physical claims that go beyond the catalog are the bimodal rotation-period distribution (§5.1) and the Sph-Prot anti-correlation for K dwarfs (§5.3). Both are most sensitive to the fast-rotator, high-Sph branch. The authors themselves flag this population: §5.3 warns that fast rotators with very large Sph resemble Type 1 CP/CB candidates but show three or fewer harmonics, so they fail the §2.2 flagging criteria and remain in the analyzed sample; §5.5 reports that Simonian et al. (2019) found the fast-rotator population dominated by tidally synchronized binaries, and roughly 70% of the likely synchronized binaries in this sample have Sph above 10^4 ppm. These known cases are flagged but not removed, and the same criteria could miss an unflagged fraction. For the central claim to hold, the unflagged fast-rotator/high-Sph tail must be predominantly single-star spot modulation. The McQuillan agreement in §5.4 validates period recovery for the 11,209 common targets, but it does not validate the physical nature of the extreme fast/high-Sph tail. The catalog entries remain useful periodicity measurements; the load-bearing risk is to the interpretation of the fast branch, not to the period estimates themselves.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a homogeneous rotation-period and photometric-activity catalog for 26,521 Kepler M and K main-sequence stars selected from KSPC DR25. Rotation periods are derived with a pipeline combining wavelet analysis, autocorrelation function, and composite spectrum applied to KEPSEISMIC light curves, with PDC-MAP light curves used as a cross-check. The authors report reliable periods for 15,910 targets, including 4,431 targets not in McQuillan et al. (2013b, 2014), and find 99.4% agreement within 2σ with that catalog for the 11,209 common targets. They also compute the photometric activity proxy Sph and report a mild decrease of rotation period with increasing temperature, an anti-correlation between Sph and rotation period for K dwarfs, and a bimodal rotation-period distribution. Potential polluters are identified and flagged, including classical pulsator/close-in binary candidates, multiple-signal systems, misclassified red giants, eclipsing binaries, and photometrically polluted light curves.","tokens_in":25021,"tokens_out":3038,"duration_ms":32212,"significance":"If the catalog is sound, this is a valuable community resource: it extends rotation measurements to fainter and cooler M and K dwarfs than the McQuillan et al. catalogs, adds roughly 4,400 new period estimates, and provides a homogeneous activity proxy tied to the rotation period. The methodology follows the pipeline previously benchmarked in Aigrain et al. (2015), and the agreement with an independent catalog on common targets is a strong validation of the period measurements. The Sph proxy has also been externally validated against solar and chromospheric activity data in prior work. The main scientific claims about bimodality and the Sph-Prot relation are consistent with earlier studies but are more fragile because they depend on the purity of the fast-rotator, high-Sph branch, as the authors themselves caution.","major_comments":[{"comment":"The physical conclusions of the paper, specifically the bimodal rotation-period distribution and the Sph-Prot anti-correlation for K dwarfs, are most sensitive to the fast-rotator, high-Sph branch. The authors state in §5.3 that fast rotators with very large Sph resemble Type 1 CP/CB candidates but have three or fewer harmonics and therefore are not flagged, and in §5.5 that approximately 70% of likely tidally-synchronized binaries have Sph above 10^4 ppm. Yet these targets remain in the main analysis and in Figs. 5-9 only some are excluded by the Berger/Simonian flags. This means an unknown fraction of the extreme fast-rotator branch may be contamination rather than single-star spot modulation. The paper needs a quantitative robustness test: for example, repeating the Prot and Sph distributions and the Sph-Prot correlation with all Berger et al. (2018) binary candidates, Simonian et al. (2019) synchronized binaries, and Type 1 CP/CB candidates excluded, and comparing with the full-sample results. Without such a test, the physical claims extend beyond what the current analysis can support.","section":"§5.3 and §5.5"},{"comment":"The catalog relies heavily on visual inspection: 6,324 of the 15,910 period estimates come from the visual check described in §3.1.2. The paper does not report any reproducibility or validation statistics for this subset, such as inter-inspector agreement, a blind test on simulated light curves, or a separate comparison with McQuillan et al. (2014) restricted to visually selected targets. Since the automatic selection thresholds are the only objective component and are applied to only ~60% of the catalog, the reliability of the remaining ~40% is not demonstrated. The authors should either provide such validation or flag the visually selected periods as lower confidence in the catalog tables.","section":"§3.1.2"},{"comment":"The red-giant removal, which removes 1,221 targets from the sample, depends on a companion paper described only as 'García et al. in prep' and on neural-network and machine-learning methods that are cited but not described in sufficient detail. Since misclassified red giants are a principal source of spurious periods, the reproducibility and correctness of this step are load-bearing for the catalog's purity. The manuscript should either summarize the red-giant identification criteria, provide a public list of the removed targets, or state that the companion paper will be submitted concurrently so that the criteria can be evaluated. As written, a reader cannot independently reproduce the red-giant mask.","section":"§2.2 and Table 4"},{"comment":"The comparison with McQuillan et al. (2013b, 2014) is reported as 99.4% agreement at 2σ for common targets, but the paper does not separate the agreement for automatically selected versus visually selected targets. Because the visual subset is subjective and forms a large fraction of the catalog, reporting the agreement separately for that subset would directly address whether the visual periods are as reliable as the automatic ones. This is a concrete and feasible analysis that would substantially strengthen the central claim.","section":"§5.4"}],"minor_comments":[{"comment":"The y-axis label in the left panels reads 'Franction' instead of 'Fraction'; please correct the typo.","section":"Figure 5"},{"comment":"The sentence 'we do not perform the rotation analysis for ... and and Type 2 and 3 CP/CB candidates' contains a doubled 'and'; please fix the wording.","section":"Appendix A"},{"comment":"The text states that Type 1 CP/CB candidates and multiple-signal targets are 'neglected in Figs. 5-9', but Fig. 5 shows a 'Full sample' distribution in black whose definition is not precisely specified. Please clarify exactly which targets are included in the 'Full sample' curves of Fig. 5.","section":"§5"},{"comment":"Several references are to works described as 'in prep', including 'Szabó et al. in prep' and 'García et al. in prep'. These should be identified or, if unpublished, the relevant data should be made available in supplementary material.","section":"§2.2 and §5.5"},{"comment":"For the 1% of targets where the Jenkins et al. (2010) correction gives a negative Sph, the paper says the correction is instead computed from the flat component of the power density spectrum but does not specify how this alternative is implemented or whether the resulting Sph values are flagged. Please add a brief description.","section":"§4"},{"comment":"The high-pass filter choices of 20, 55, and 80 days are described, but the criterion for selecting the 'appropriate filter' in §3.1.1 is described in words; a compact table or equation summarizing the period-dependent filter priority would improve clarity.","section":"§2.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a useful catalog paper, but the physical conclusions (bimodality and Sph-Prot relation) outrun what the current validation demonstrates. The authors are appropriately cautious about the fast-rotator branch, but caution is not the same as a quantitative robustness test. I would encourage the editor to request the additional analyses described in the major comments, especially a separation of the McQuillan comparison by selection method and a version of the key figures with all externally and internally flagged binary/CP candidates excluded. The catalog itself is likely publishable even if the physical claims require qualification."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The two things you should know: this is a useful catalog paper, and the part most likely to be over-read is the physical interpretation of the fast-rotator branch. The catalog itself is the real deliverable.\n\nThe genuinely new content is the longer Kepler baseline (through Q17), the use of the KEPSEISMIC light curves with three high-pass filters, and the careful polluter screening. That yields 4,431 rotation periods for M and K dwarfs that McQuillan et al. did not report, mostly for cooler and fainter targets. The 99.4% two-sigma agreement with McQuillan on 11,209 common targets is a strong external consistency check, and the Sph proxy has independent validation from solar data and chromospheric activity comparisons. The multi-method wavelet/ACF/composite-spectrum pipeline was already established, but applying it to this sample and releasing the tables and light curves is a real service to the community.\n\nSoft spots, in proportion. About 40% of the period estimates (6,324 targets) come from visual inspection. That is not disqualifying in this kind of work, but it means the catalog is partly an expert-judgment product, and the paper's Section 3.1.2 does not codify the visual criteria. The red-giant and RR Lyrae removal leans on unpublished companion papers—fine if the lists are released with the catalog, but a reviewer should ask for those lists. The bigger substantive caveat: the paper's own sections 5.3 and 5.5 flag that fast rotators with very large Sph resemble Type 1 CP/CB candidates and that Simonian et al. found the fast-rotator population dominated by tidally synchronized binaries. The bimodality and the Sph-Prot anti-correlation are most sensitive to exactly that branch. The stress-test note is right that the McQuillan agreement validates period recovery, not the physical nature of the extreme fast/high-Sph tail. So the catalog entries stand, but claims like \"fast rotators are spot-dominated\" have a load-bearing caveat.\n\nWho is this for? Anyone doing gyrochronology, stellar activity, or exoplanet host characterization with M and K dwarfs. The catalog will get cited, and should be. It deserves a serious referee. My recommendation: publish after requesting the identification lists, codifying the visual-inspection criteria, and softening the interpretation of the fast-rotator sequence.","headline":"Solid catalog paper whose main asset is 15,910 rotation periods and Sph values, with 4,431 new detections; the interpretation of the fast-rotator branch is shakier than the catalog itself.","tokens_in":25696,"tokens_out":1210,"would_cite":true,"duration_ms":14293,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Wavelet-plus-autocorrelation analysis of Kepler light curves yields reliable rotation periods for 15,910 M and K dwarfs, 4,431 of them new, and shows faster rotators are more active.","keywords":["stars: low-mass","stars: rotation","stars: activity","starspots","techniques: photometric","methods: data analysis","catalogs","gyrochronology"],"falsifier":"Take the fast-rotating K dwarfs with photometric activity above roughly 10,000 ppm — the group the paper itself cautions about — and obtain high-resolution spectra or Gaia astrometric orbits for about 100 of them. If more than a small fraction show binary motion, the fast-rotator branch of the activity–rotation relation, and part of the derived bimodal period distribution, is contaminated by synchronized binaries rather than single-star rotation. A purely photometric version of the same test: require every reported period to be recovered independently from the first and second halves of the light curve; stable spot modulation should persist across both halves, while quarter-dependent pollution would not.","tokens_in":24552,"feed_emoji":"🔄","tokens_out":16229,"duration_ms":128353,"temperature":0.7,"pith_summary":"This paper tries to establish that a rotation-detection pipeline combining wavelet time-frequency analysis, the autocorrelation function, and their product (the composite spectrum) can extract trustworthy surface rotation periods from Kepler light curves of low-mass stars, going beyond what an autocorrelation-only analysis achieves. Applied to 26,521 M and K main-sequence dwarfs, the pipeline reports periods and photometric activity levels for roughly 60 percent of the sample, including 4,431 stars with no previously published period. For the 11,209 targets in common with the earlier McQuillan catalog, the two analyses agree within $2\\sigma$ for 99.4 percent. The resulting sample supports three physical claims: hotter dwarfs rotate faster, faster rotators show higher photometric activity (for K dwarfs in particular), and the rotation-period distribution is bimodal. A careful reader would care because rotation periods feed gyrochronology — the use of spin-down to estimate stellar ages — and a larger, cleaner period catalog sharpens that age clock.","feed_headline":"Spin catalog for 15,910 cool stars adds 4,431 new periods","feed_subtitle":"Wavelet-plus-ACF pipeline recovers rotation for Kepler M and K dwarfs, agreeing with prior values at 99.4 percent.","key_machinery":"The object that carries the argument is the three-diagnostic rotation pipeline. First, a Morlet-wavelet time-frequency decomposition yields the wavelet power spectrum and, summed over time, the global wavelet power spectrum (GWPS); the rotation estimate is the central period of the highest fitted Gaussian peak, with the half-width as its uncertainty. Second, the autocorrelation function (ACF) of the light curve, smoothed with a Gaussian whose width is one tenth of the dominant Lomb-Scargle period, yields the highest significant peak as an independent period. Third, the composite spectrum (CS) is the product of the normalized GWPS and the normalized ACF resampled onto the same period grid, so that periods present in both methods are enhanced. A period is accepted when GWPS, ACF, and CS agree within $2\\sigma$, the estimates agree within 20 percent across the appropriate high-pass filters, and the peak-height thresholds hold ($G_{\\rm ACF} \\geq 0.2$, $H_{\\rm ACF} \\geq 0.3$, $H_{\\rm CS} \\geq 0.15$). The companion object is the activity proxy $S_{\\rm ph}$, defined as the standard deviation of light-curve sub-series of length $5 \\times P_{\\rm rot}$ with photon noise subtracted; it converts a measured period into a magnetic-activity measure and underlies the rotation–activity relation.","core_discovery":"The paper's central claim is that its KEPSEISMIC light curves — its own calibration of Kepler pixel data, high-pass filtered at 20, 55, and 80 days — together with a three-diagnostic rotation pipeline yield reliable rotation-period estimates for 15,290 M and K dwarfs, with a further 620 flagged candidates (possible classical pulsators or close-in binaries, and multi-period sources) reported separately. Reliability is defined internally as agreement among the wavelet global power spectrum, the autocorrelation function, and the composite spectrum within $2\\sigma$, plus agreement within 20 percent between different filters, plus peak-height thresholds; externally it is validated by 99.4 percent agreement within $2\\sigma$ with the McQuillan et al. catalog across 11,209 common targets. The discovery is the catalog itself and the population statements built on it: the newly measured stars are on average cooler and fainter than previously cataloged rotators, rotation periods shorten as effective temperature and mass increase, the photometric activity proxy $S_{\\rm ph}$ spans a wider range for hotter stars, $S_{\\rm ph}$ grows as rotation speeds up, and the rotation-period distribution is bimodal. The paper also claims to clean the sample by removing 1,221 misclassified red giants, eclipsing binaries, RR Lyrae stars, and photometrically polluted light curves.","pith_inferences":["If the pipeline resolves factor-of-two period ambiguities as well as the comparison with the prior catalog suggests, applying the same three-diagnostic scheme to older G and F dwarfs could sharpen the debate over weakened magnetic braking, since spurious half-rotation signals are exactly the kind of error that mimics faster-than-expected spin-down.","A testable extension: the 3,562 targets with visible but unrecoverable spot modulation form a natural completeness sample, so modeling their detection failure as a function of brightness and period would let catalog users correct selection bias in the reported period distribution.","Because the paper flags fast rotators with very large $S_{\\rm ph}$ as possibly tidally synchronized binaries, gyrochronology or activity studies that use the fast-rotator branch should test how their conclusions change when high-activity fast rotators are excluded.","The wider $S_{\\rm ph}$ range at higher temperature, combined with the catalog's detection limits, suggests a re-interpretation test: if intermediate-period stars are under-detected because dark spots and bright faculae cancel, the bimodal gap should widen when the analysis is restricted to low-noise, high-amplitude light curves."],"forward_implications":["The 4,431 newly reported periods extend rotation measurements to cooler and fainter M and K dwarfs, giving gyrochronology relations more calibration points at the low-mass end.","Confirmation of the bimodal rotation-period distribution by an independent pipeline strengthens the case that the bimodality is a real feature of the low-mass field-star population rather than an artifact of one method.","The K-dwarf anti-correlation between $S_{\\rm ph}$ and rotation period ties photometric activity to spin rate and anchors comparisons with solar activity, which the paper quotes at 314.5 ppm at maximum and 67.4 ppm at minimum.","Removal of 1,221 misclassified red giants and flagging of 368 pulsator/close-binary candidates changes the effective M and K dwarf population counts derived from the Kepler catalog.","The catalog tables, with pollution and binarity flags, become a reference data set for stellar spin-down studies, exoplanet-host rotation, and tests of magnetic-braking models."],"supporting_citations":[{"why":"Prior rotation-period catalog for Kepler low-mass stars; provides the 11,209 common targets used for external validation and the baseline this work extends by 4,431 new periods.","marker":"McQuillan et al. (2013b, 2014)"},{"why":"Supplies the wavelet/autocorrelation/composite-spectrum methodology and the automatic-selection thresholds the pipeline applies.","marker":"Ceillier et al. (2016, 2017)"},{"why":"Simulation-based comparison cited as evidence that this combined methodology beats autocorrelation alone in completeness and reliability.","marker":"Aigrain et al. (2015)"},{"why":"Kepler Stellar Properties Catalog (DR25) that defines the M and K dwarf sample and supplies effective temperatures, masses, and surface gravities.","marker":"Mathur et al. (2017)"},{"why":"Defines KADACS, the calibration software that produces the KEPSEISMIC light curves used throughout.","marker":"García et al. (2011)"},{"why":"Defines the photometric activity proxy Sph as the standard deviation of sub-series of length 5 × Prot.","marker":"Mathur et al. (2014)"},{"why":"Gaia-based catalog used to flag misclassified red giants and candidate binaries in the sample.","marker":"Berger et al. (2018)"},{"why":"Identifies fast rotators as likely tidally synchronized binaries, motivating the paper's caution on the fast-rotator branch.","marker":"Simonian et al. (2019)"},{"why":"First published combination of wavelet analysis with the autocorrelation function for rotation periods, the direct antecedent of this pipeline.","marker":"García et al. (2014a)"}],"fun_headline_variants":["4,431 new rotation periods for cool Kepler stars","Bimodal spin periods in 15,910 M/K dwarfs","Cool star spins mapped: 15,910 periods, 4,431 new","Kepler cool stars: 15,910 spin periods, new bimodality","15,910 cool star rotation periods, 4,431 fresh"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that each detected periodic brightness variation is starspot rotation, not pulsation, binarity, eclipses, or light from a nearby star; the paper itself flags 368 classical-pulsator/close-binary candidates, 270 multi-signal light curves, and a set of very fast rotators with unusually large activity as cases where this premise is uncertain.","fun_headline_variants_meta":{"raw":{"variants":["4,431 new rotation periods for cool Kepler stars","Bimodal spin periods in 15,910 M/K dwarfs","Cool star spins mapped: 15,910 periods, 4,431 new","Kepler cool stars: 15,910 spin periods, new bimodality","15,910 cool star rotation periods, 4,431 fresh"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000707,"raw_usage":{"total_tokens":3279,"prompt_tokens":1131,"completion_tokens":2148,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":747,"completion_tokens_details":{"reasoning_tokens":2053}},"tokens_in":747,"tokens_out":2148,"duration_ms":14166,"temperature":1.0,"reasoning_tokens":2053,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:19:58.061719+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the fast-rotating K dwarfs with photometric activity above roughly 10,000 ppm — the group the paper itself cautions about — and obtain high-resolution spectra or Gaia astrometric orbits for about 100 of them. If more than a small fraction show binary motion, the fast-rotator branch of the activity–rotation relation, and part of the derived bimodal period distribution, is contaminated by synchronized binaries rather than single-star rotation. A purely photometric version of the same test: require every reported period to be recovered independently from the first and second halves of the light curve; stable spot modulation should persist across both halves, while quarter-dependent pollution would not.","supporting_citations":[],"review_version":1}