{"id":"bbc00071-1a0f-4516-b638-56656044ffe7","arxiv_id":"2506.14869","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A search of Gaia spectroscopic binaries in TESS data yields 112 heartbeat stars and shows that the fraction of such systems rises sharply with stellar temperature.","lead":"Using Gaia's catalog of binary stars and TESS brightness measurements, the authors found 112 new heartbeat stars, binaries whose light curves resemble an electrocardiogram. The work sharpens the picture of how such binaries form and evolve and measures the masses and radii of two systems.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Fig. 12 claim that the HB fraction rises rapidly with effective temperature is not yet supported: the search's detection efficiency as a function of color is unmodeled, and the steep drop across the Kraft break could be a selection effect.","rationale":"The paper does useful work: a clean Gaia-SB-seeded search, a large catalog, orbital agreement tests, and two PHOEBE mass/radius measurements. My concern is not about internal consistency of the fits but about whether the central population claim can be inferred from the selection. The reader's verdict identifies the same weakest assumption. What would have to be true is that TESS HB detection is nearly complete, or at least color-independent, across the sample. The pipeline itself provides reason to doubt this: the selection threshold is defined in chi2 space, and low-amplitude signals are less likely to cross it; the visual confirmation step adds another sensitivity floor. The authors' own physical argument says amplitudes decrease toward cooler stars, so the selection function is expected to correlate with the very axis on which the fraction is claimed to drop. This is an externally checkable selection effect, not a circularity; a recovery simulation resolves it. In the meantime the headline claim is plausible but conditional. I therefore keep the reader's CONDITIONAL verdict (no change), and agree with the weakest-assumption choice.","tokens_in":17363,"tokens_out":2742,"duration_ms":29718,"concrete_test":"Inject synthetic heartbeat signals into the real TESS QLP light curves of the Gaia SB1 sample. For each BP-RP bin, draw orbital parameters from the observed HB distributions, scale the Kumar-model amplitude using the same tidal/evolutionary scalings invoked in §4 (or from the observed amplitude-color relation), add the signals to the actual noise realizations, and rerun the full SB1 pipeline (LS periodogram, 1-hr binning, Trust-Region Kumar fit, R<0.5 and e>0.15 cut, visual inspection if feasible). Compute recovery fraction f_rec(BP-RP). Recompute Fig. 12 as observed fraction divided by f_rec. If the corrected fraction still rises steeply across the Kraft break, the astrophysical claim stands; if it flattens or reverses, the headline population trend is a selection artifact.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central population claim (abstract; §4, Fig. 12) is that the fraction of Gaia SB1/SB2 systems that are heartbeat stars rises rapidly with effective temperature and that non-giant HBs have evolved off the main sequence. For the fraction claim to hold, the HB detection probability must be roughly flat across BP-RP after the giant and period cuts. The paper does not establish this. The SB1 search pipeline selects systems with R = chi2_HB/chi2_line < 0.5 and e > 0.15 (§2.2, Fig. 3), followed by visual inspection. Cooler, lower-mass stars are expected to have smaller tidal amplitudes and, below the Kraft break, different damping (§4), so the same physical trend that the paper wants to measure also makes their HB signals harder to distinguish from noise in TESS QLP light curves. No injection/recovery or completeness simulation is presented. Consequently the drop in Fig. 12 from roughly 10^-1 near BP-RP < 0.5 to roughly 10^-4 near BP-RP ~ 1 could be dominated by declining sensitivity rather than by the underlying binary population. The luminosity-offset claim is more robust because it is a median comparison, but the fraction-vs-temperature claim is directly at risk. The paper even notes systematic catalog effects for period statistics (§4) but does not model this color-dependent completeness.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper searches for heartbeat stars among Gaia DR3 single-lined (SB1) and double-lined (SB2) spectroscopic binaries using TESS Quick-Look Pipeline light curves. The authors identify 112 new heartbeat systems, fit their phase-folded light curves with the analytic Kumar et al. (1995) model, and compare the fitted periods, eccentricities, and arguments of periastron with the Gaia orbital solutions. For the two SB2 systems with apparently consistent orbits, they use PHOEBE, with Gaia velocity semi-amplitudes as priors, to derive stellar masses and radii. The paper also presents a statistical analysis of the heartbeat-star fraction as a function of color, concluding that non-giant heartbeat stars are evolved off the main sequence and that the fraction of binaries that are heartbeat stars rises rapidly with effective temperature.","tokens_in":17637,"tokens_out":4563,"duration_ms":44762,"significance":"If the population-level conclusions hold, this is one of the largest homogeneous samples of short-period heartbeat stars and provides a new route to discovering them from spectroscopic-binary catalogs. The two PHOEBE mass/radius measurements are a useful addition to the sparse set of heartbeat stars with direct stellar parameters. The main population claim—that the heartbeat fraction rises rapidly with effective temperature and drops across the Kraft break—would be an important constraint on tidal dissipation and binary evolution, but it is not yet securely established because the search's detection efficiency as a function of color is not characterized. The paper's catalog of 112 new HBs and its comparison with Gaia orbits are of independent value.","major_comments":[{"comment":"The central claim that the HB fraction rises rapidly with effective temperature is not yet supported because the search's detection efficiency as a function of color is not modeled. Cooler, lower-mass stars are expected to have smaller tidal amplitudes and, below the Kraft break, different damping timescales, making their heartbeat signals harder to detect in TESS data. The selection region (R<0.5 and e>0.15, §2.2 and Fig. 3) was calibrated using only the 10 SB2 HBs, and the final sample also depends on visual inspection; both steps are likely color-dependent. No injection/recovery or completeness simulation is presented. The drop in Fig. 12 from roughly 10^-1 near BP-RP<0.5 to roughly 10^-4 near BP-RP~1 could therefore be dominated by declining sensitivity rather than by the underlying binary fraction. Please provide a completeness correction or a quantitative argument that detection efficiency is flat across color after the period and giant cuts.","section":"§4, Fig. 12"},{"comment":"The statement that 85% of SB1 HBs and 20% of SB2 HBs have orbital parameters that are 'consistent' with the Gaia solutions is not verifiable because the agreement criterion is not defined. The dashed lines in Fig. 8 are described only as 'the range we consider a reasonable match.' Please specify the tolerance explicitly (e.g., fractional period difference, eccentricity difference) and show how the 85%/20% numbers depend on that choice. This is load-bearing because the comparison underpins the PHOEBE modeling and the discussion of Gaia orbit quality, and because a loose criterion would make the agreement rate trivially high.","section":"§3, Fig. 8"},{"comment":"The definition of the 'fraction' plotted in Fig. 12 is not clearly specified. The axis label 'Fraction (NHB/Nnot)' suggests N_HB / N_non-HB, while the text says 'fraction of HBs' and the figure caption says 'median fraction of HBs.' If the denominator excludes the HBs themselves, the values will differ from N_HB/N_total, especially at the blue end where the fraction is around 10^-1. Please state the exact definition and ensure it is applied consistently in both panels, as this quantity is the basis of the paper's principal statistical conclusion.","section":"§4, Fig. 12"}],"minor_comments":[{"comment":"The formula for the true anomaly is garbled; it should read ν = 2 tan^{-1}( sqrt(1+e)/sqrt(1-e) tan(E/2) ). Please correct.","section":"§2.1, Eq. (2)"},{"comment":"The sentence 'the score statistic seems to be a limited indicator of the Gaia orbit quality, as almost all of the targets with incorrect Gaia P, e, or ω also have S<0.587' is confusing: finding that bad orbits also pass the nominal good-orbit cutoff does show the score is not sufficient, but the phrasing reads as if the low scores themselves are the problem. Please rephrase.","section":"§3"},{"comment":"The legend entry 'Gaia Clean Score Stars' is not defined in the caption; state that this refers to SB1 systems with a Bashi et al. (2022) score S<0.587.","section":"Fig. 5"},{"comment":"The table references the supplementary file 'SBTABLE.full' but gives no description of its columns or how the full table can be accessed; please add a short explanation.","section":"Table 1"},{"comment":"The phrase '85% of the single-line spectroscopic binaries' could be misread as applying to all Gaia SB1s; suggest '85% of the single-line spectroscopic binary heartbeat stars' for clarity.","section":"Abstract"},{"comment":"There is a stray line 'Screenshot from 2025-06-20 14-36-03.png' in the text; please remove it.","section":"Before Fig. 11"}],"recommendation":"major_revision","confidential_remarks":"The completeness issue in Fig. 12 is the main obstacle; if the authors can add injection-recovery tests or a convincing completeness model, the paper would be a strong contribution. The agreement-criterion ambiguity is easy to fix. The paper is within scope for a specialist astronomy journal, and the sample will be useful regardless of the population claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this paper gives the community a genuinely useful new catalog—about 108 truly new heartbeat stars rather than the 112 claimed, since four overlap with Solanki et al. (2025)—plus two new PHOEBE mass/radius measurements for massive stars in HB systems. The catalog and the two stellar parameter measurements are solid. The population trend in Fig. 12 is interesting but not yet proven.\n\nWhat is new: starting from the Gaia SB1/SB2 catalogs instead of searching blindly through TESS light curves is a sensible variant of the Kołaczek-Szymański et al. (2021) approach, and doing it at this scale with a semi-automated preselection is new. The orbital parameter comparison to Gaia, the TEO detections, and the two PHOEBE models all add value. The mass/radius values are consistent with detached eclipsing binary results, which is a good sanity check.\n\nSoft spots, in order of importance. First, the central population claim—that the HB fraction rises rapidly with effective temperature—is undercut by an unmodeled detection efficiency. Cooler, lower-mass stars have smaller tidal amplitudes and longer damping times, so the same physical effects the paper wants to measure also make those systems harder to find in TESS. No injection/recovery or completeness simulation is presented. The drop in Fig. 12 from ~10^-1 to ~10^-4 across the Kraft break could be largely a selection effect. The luminosity-offset claim is more robust because it is a median comparison, but the fraction-vs-temperature statement in the abstract is not yet supported as written. Second, the criterion for 'agreement' between TESS and Gaia orbital parameters is never quantified; the 85% and 20% numbers need a tolerance. Third, the '112 new' is internally inconsistent; four were already in Solanki et al. (2025). Minor: the quoted primary radius for TIC 98552498 in the text (1.41±0.40) does not match the figure posterior (1.21+0.60-0.19); the text and figures should be reconciled.\n\nNone of these are fatal. The catalog will be cited, and the two mass/radius measurements are useful. The paper reads like honest work with a clear understanding of the systematics—it even notes some catalog effects in §4. It just needs a serious revision to tighten the numbers and add (or clearly flag the absence of) a completeness treatment.\n\nRecommendation: send it to peer review. The right referee will ask for a defined agreement metric, a completeness test or an explicit limitation statement, and a cleanup of the catalog and radius inconsistencies. That is a conditional accept, not a desk reject.","headline":"New catalog of ~108 TESS heartbeat stars from Gaia SB1/SB2 samples; the temperature trend is promising but needs completeness modeling.","tokens_in":18205,"tokens_out":4226,"would_cite":true,"duration_ms":38494,"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":"Heartbeat-star pulses are a signature of hotter, slightly evolved binary primaries, not a random slice of close binaries.","keywords":["heartbeat stars","spectroscopic binaries","TESS photometry","Gaia DR3","tidal deformation","stellar evolution","eclipsing binaries","tidally excited oscillations"],"falsifier":"Simulate injecting synthetic heartbeat light curves with realistic period, eccentricity, and amplitude distributions into TESS photometry across the full color range, then run the same search pipeline to measure detection completeness as a function of color; if completeness falls for red stars, the observed drop in heartbeat fraction across the Kraft break could be a selection effect rather than a property of the binary population.","tokens_in":17152,"feed_emoji":"💓","tokens_out":9350,"duration_ms":82661,"temperature":0.7,"pith_summary":"Heartbeat stars are close, eccentric binaries whose mutual tides make their phase-folded light curves look like an electrocardiogram. The paper tries to establish that these systems are not a random slice of short-period binaries but are concentrated in hotter, slightly evolved stars. It does this by starting from known spectroscopic binaries in the Gaia catalog rather than scanning the whole sky, then folding their TESS light curves to find 112 new heartbeat systems. On the color-magnitude diagram the non-giant heartbeat stars lie above the main sequence at fixed color, meaning their primaries have begun to evolve, and the fraction of systems showing heartbeat pulses climbs steeply with stellar temperature, topping 1% among the hottest single-lined binaries. A sympathetic reading of the result is that heartbeat visibility is set by the balance between stellar evolution and tidal dissipation.","feed_headline":"112 new heartbeat stars found, and they favor hot evolved binaries","feed_subtitle":"Gaia binary list plus TESS pulses shows heartbeat fraction jumps with stellar temperature.","key_machinery":"Eq. (1) of the paper: a normalized flux model of the form $Z + S [1 - 3 \\sin^2 i \\sin^2(\\nu + \\omega)] / (1 - e \\cos E)^3$, with the true anomaly $\\nu$ and eccentric anomaly $E$ connected to period and time through Kepler's equation. This analytic tidal model converts the electrocardiogram-like phase-folded pulse into estimates of period, eccentricity, inclination, and argument of periastron without a full binary light-curve fit. The search uses known Gaia spectroscopic binaries as the parent sample, selects candidates by a chi-squared ratio against a straight-line fit plus an eccentricity cut, and then filters visually. The statistical part of the argument compares color-magnitude positions and heartbeat fractions in color bins, interpreted through the competition between stellar evolutionary time scales, which set the tidal amplitude, and dissipation time scales, which set how long the pulse lasts.","core_discovery":"The central claim is that heartbeat stars preferentially appear in hotter, slightly evolved binary systems rather than being a uniform sample of short-period binaries. The authors start from 186,905 spectroscopic binaries in the Gaia DR3 catalog, fold their TESS light curves, and identify 112 heartbeat systems with periods between 1.5 and 12.2 days and eccentricities up to 0.57. For single-lined systems, 85% of the periods and eccentricities agree with the Gaia orbital solutions, while only two of the ten double-lined systems agree, a mismatch the authors attribute to sparse radial-velocity sampling and the hot, early-type nature of these stars. For those two double-lined systems, the light curve together with the Gaia velocity amplitudes yields component masses and radii consistent with massive detached eclipsing binaries. On the color-magnitude diagram the non-giant heartbeat stars are more luminous at fixed color than their parent binary samples, and the heartbeat fraction rises rapidly with effective temperature, exceeding 1% for the hottest single-lined systems.","pith_inferences":["If the temperature trend survives a completeness correction, it predicts that hot binaries of the same period and eccentricity should outnumber cool ones as heartbeat sources by a large factor; that can be tested by injecting synthetic heartbeat signals into the same TESS data.","The same search strategy applied to binaries selected by radial-velocity scatter rather than a full orbital solution could recover longer-period heartbeat systems that short TESS sectors alias.","With better radial velocities, the two successfully modeled double-lined systems show the route to homogeneous masses and radii of massive stars from heartbeat light curves alone.","Because convective-envelope depth sets the dissipation rate, the model implies the heartbeat fraction at fixed temperature should depend on metallicity, a testable prediction once larger samples exist."],"forward_implications":["Among hot short-period Gaia single-lined binaries, more than 1% show heartbeat pulses, so the phenomenon is common rather than rare in that regime.","Because non-giant heartbeat stars sit above the main sequence at fixed color, finding one is a sign that the primary has started evolving off the main sequence.","The rapid drop in heartbeat fraction across the Kraft break implicates convective envelopes and long evolutionary timescales as the suppression mechanism.","Only 2 of 10 double-lined heartbeat systems have Gaia orbital solutions in agreement with the light-curve orbits, so short-period hot SB2 solutions from Gaia DR3 need independent radial velocities or DR4 data.","The 18 eclipsing and 10 tidally oscillating systems among the 112 provide specific targets for measuring tides and stellar interiors."],"supporting_citations":[{"why":"Supplies the analytic tidal flux model used for every phase-folded light-curve fit in the search.","marker":"Kumar et al. (1995)"},{"why":"Defined heartbeat stars and demonstrated that this analytic model reproduces their Kepler light curves.","marker":"Thompson et al. (2012)"},{"why":"Provides the DR3 single-lined and double-lined spectroscopic binary catalogs that define the parent population.","marker":"Gaia Collaboration et al. (2023)"},{"why":"Supplies the score statistic used to gauge which Gaia single-lined orbital solutions are reliable.","marker":"Bashi et al. (2022)"},{"why":"Gives the evolutionary-versus-dissipation scaling used to interpret why heartbeat amplitudes and fractions depend on stellar temperature.","marker":"MacLeod & Loeb (2025)"},{"why":"Defines the effective-temperature break separating convective from radiative envelopes, the boundary across which the heartbeat fraction drops.","marker":"Kraft (1967)"},{"why":"Provides the OGLE heartbeat catalog used as the main comparison sample of longer-period and giant systems.","marker":"Wrona et al. (2022a,b)"},{"why":"Earlier TESS heartbeat search whose catalog overlaps four of the new systems and anchors comparison with prior work.","marker":"Solanki et al. (2025)"},{"why":"Supplies the tidal circularization envelope against which the new periods and eccentricities are checked.","marker":"Mazeh (2008)"}],"fun_headline_variants":["112 new heartbeat stars found in hot evolved binaries","Heartbeat stars favor hot evolved binaries: 112 new ones","Gaia and TESS find 112 heartbeat stars, hottest pulse most","Heartbeat star rate rises with temperature: 112 new finds","Hotter binaries mean more heartbeat stars: 112 discovered"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The temperature trend rests on the assumption that the search detects heartbeat pulses around cool stars about as easily as around hot stars.","fun_headline_variants_meta":{"raw":{"variants":["112 new heartbeat stars found in hot evolved binaries","Heartbeat stars favor hot evolved binaries: 112 new ones","Gaia and TESS find 112 heartbeat stars, hottest pulse most","Heartbeat star rate rises with temperature: 112 new finds","Hotter binaries mean more heartbeat stars: 112 discovered"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00025,"raw_usage":{"total_tokens":1568,"prompt_tokens":975,"completion_tokens":593,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":591,"completion_tokens_details":{"reasoning_tokens":509}},"tokens_in":591,"tokens_out":593,"duration_ms":7095,"temperature":1.0,"reasoning_tokens":509,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:10:04.533876+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate injecting synthetic heartbeat light curves with realistic period, eccentricity, and amplitude distributions into TESS photometry across the full color range, then run the same search pipeline to measure detection completeness as a function of color; if completeness falls for red stars, the observed drop in heartbeat fraction across the Kraft break could be a selection effect rather than a property of the binary population.","supporting_citations":[{"cited_title":"O., Quataert E","cited_arxiv_id":null,"evidence_quote":"Supplies the analytic tidal flux model used for every phase-folded light-curve fit in the search."},{"cited_title":"Tidal Amplitudes in the Magellanic Cloud Population of Heartbeat Stars","cited_arxiv_id":"2503.17133","evidence_quote":"Gives the evolutionary-versus-dissipation scaling used to interpret why heartbeat amplitudes and fractions depend on stellar temperature."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier TESS heartbeat search whose catalog overlaps four of the new systems and anchors comparison with prior work."}],"review_version":1}