{"id":"8e42f4eb-bf81-4c2a-a4b3-d640250bd78f","arxiv_id":"2608.12474","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A new search of APOGEE binaries with TESS light curves finds 50 heartbeat stars, including 36 new systems, and confirms that the detected fraction rises sharply with stellar temperature.","lead":"Astronomers identified 50 heartbeat stars, close binary pairs that are visibly stretched when they swing past each other, by matching APOGEE binary stars with TESS brightness data. The new sample supports the earlier claim that heartbeat stars are most common among hot, massive stars and are rare among cool, low-mass stars.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The headline population trend is not yet separated from the amplitude-limited selection function: the paper's own Sec. 1 predicts a decreasing detected HB fraction at low T, and Fig. 9 has no completeness correction, so the color trend could be a sensitivity curve.","rationale":"I agree with the reader's weakest-assumption identification. The paper is careful and the 50 detections, RV work, and TEO flags are useful, but the two headline claims in the abstract are population statements that require a known selection function. The manuscript itself supplies the physical reason the selection function must vary with color, and it supplies no completeness correction for Fig. 9 or Fig. 8. This is not a disagreement with consensus; it is an internal consistency issue: the discussion in Sec. 1 predicts exactly the trend that Sec. 4 presents as a confirmation. A flat underlying HB fraction with a steep amplitude sensitivity would produce the same figure. The 'evolved off the main sequence' claim is similarly affected because evolved radius boosts tidal amplitudes, so an amplitude-limited search preferentially finds evolved systems. I would not reject the paper; the data and detections stand, and the claim can likely be repaired with an injection-recovery completeness estimate or by restating the conclusions as describing the detected fraction in an amplitude-limited search. Since the reader already requested this qualification and assigned CONDITIONAL, my assessment does not move the verdict.","tokens_in":18369,"tokens_out":3705,"duration_ms":37917,"concrete_test":"Injection-recovery test on the actual QLP light curves: take a random subset of the 31,548 APOGEE MS binaries with T<13.5 and inject synthetic HB signals drawn from a model with a flat underlying HB fraction across BP-RP. The injected signals should use the fitted period/eccentricity/inclination distribution from Table 1 and physical tidal amplitudes from the Kumar model scaled by stellar parameters, with noise matching each target's QLP errors. Run the exact search pipeline (Lomb-Scargle, PDM, HB fit, R<0.5, e>0.10, and the same visual-inspection criteria) and measure the recovered fraction per 0.1-mag color bin. If the recovered fraction reproduces the Fig. 9 decline, the headline trend is a selection effect; if the recovery fraction is flat or has a different shape and the decline persists after dividing Fig. 9 by it, the physical interpretation is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claims are the elevated CMD position of HB primaries and the rapid rise of the HB fraction with effective temperature (abstract; Sec. 4; Fig. 9). The load-bearing assumption is that the color trend in Fig. 9 reflects the underlying binary population rather than the search's detection efficiency. This assumption is undercut by the paper's own physical discussion: 'HB amplitudes reflect a competition between the rate of stellar evolution and the amount of damping ... Together, these effects lead to a decreasing frequency of HBs at lower temperatures in any amplitude-limited search' (Sec. 1). The APOGEE search is amplitude-limited in practice: candidates must beat a linear fit by R<0.5 and e>0.10 in a single TESS sector, and the detected amplitude is set by tidal deformation, which is larger for evolved, hotter stars. Fig. 9 computes NHB/Nbinaries per color bin with no correction for the fraction of systems whose HB signal is detectable at the typical QLP precision. If the sensitivity is much higher in blue bins than red bins, the observed decline would be entirely instrumental. The same selection also biases the CMD result (Fig. 8): evolved primaries with larger radii produce larger-amplitude HBs and are preferentially included, so the displacement above the equal-mass binary isochrone could be a detection bias rather than evidence about the parent population. The paper acknowledges the amplitude mechanism in Sec. 1 but does not let it qualify the Sec. 4 confirmatory language.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a search for heartbeat stars among 31,548 main-sequence (log g >= 3.25, T < 13.5) APOGEE spectroscopic binaries using TESS QLP light curves from sectors 1-79. The authors fit each phase-folded light curve with the Kumar et al. (1995) analytic model, select candidates with R < 0.5 and e > 0.10, visually vet them, and report 50 heartbeat stars. They measure orbital parameters, identify eclipses and tidally excited oscillations, and use APOGEE radial velocities for eight systems to derive mass functions and secondary masses. The paper's headline claims are that non-giant heartbeat primaries have begun to evolve off the main sequence and that the fraction of main-sequence binaries that are heartbeat stars rises rapidly with effective temperature, reaching about 1% for BP-RP < 0.5 mag.","tokens_in":18622,"tokens_out":9997,"duration_ms":83760,"significance":"If the population trends are robust, this is a valuable independent confirmation of C25 with a different selection channel, adding 36 new HBs, and the TEO detections and RV mass functions provide useful targets for follow-up. The main weakness is that both headline claims are population statements made without correcting for the amplitude-limited nature of the search, which the authors themselves describe in Sec. 1. As a detection paper, the catalog is solid; as a population paper, it requires a completeness correction or a substantially weakened interpretation.","major_comments":[{"comment":"The claim that the HB fraction rises rapidly with effective temperature is not established because the search is amplitude-limited and the detection efficiency is expected to vary strongly with colour. The paper states in Sec. 1 that 'Together, these effects lead to a decreasing frequency of HBs at lower temperatures in any amplitude-limited search, like what was done in C25'; the present search uses the same R < 0.5, e > 0.10 selection on single-sector QLP light curves. Fig. 9 shows NHB/Nbinaries per colour bin with no completeness correction. Because hotter, more evolved stars produce larger tidal deformation amplitudes and pass the threshold more easily, the observed decline toward the Kraft break could be entirely a sensitivity curve. Please add an injection-recovery calculation of the detection efficiency as a function of colour/T_eff, or restrict the conclusion to the detected sample and remove 'confirm' from the abstract.","section":"Sec. 4, Fig. 9"},{"comment":"The CMD result is subject to the same selection bias. If the search preferentially detects primaries with larger radii, the binned median magnitude of the detected HBs will lie above the equal-mass binary isochrone even if the underlying population is unevolved. The argument that unequal masses add further support to the evolved-primary hypothesis assumes the detected sample is representative. Without a selection correction, Fig. 8 cannot distinguish 'primaries have evolved off the MS' from 'the search finds the most easily detectable, hence most evolved, HBs.' The authors should either correct for this bias or soften the claim.","section":"Sec. 4, Fig. 8"},{"comment":"The denominator of the fraction in Fig. 9 is defined by the APOGEE v_scatter >= 3 km/s criterion, but the completeness of this binary parent sample as a function of colour is not demonstrated. APOGEE radial-velocity precision and the ability to detect a given v_scatter depend on spectral type and line broadening, so the parent sample may itself introduce a colour-dependent selection that mimics or masks an HB abundance trend. Please quantify or at least discuss the colour dependence of the binary parent sample.","section":"Sec. 2 / Sec. 4, Fig. 9"}],"minor_comments":[{"comment":"The Fig. 7 caption contains the stray text 'Screenshot from 2026-07-28 15-43-31.png' that should be removed.","section":"Fig. 7 caption"},{"comment":"Table 2 uses 'max(D_RV0,-1)' while the text defines D± and the cut max(D±) > 0.25; define the phase-completeness statistic once and use matching notation.","section":"Table 2"},{"comment":"The statement that 1,167,192 models were tried for 31,548 binaries is not derivable from the described procedure (three period harmonics times two starting points); clarify the source of the factor of about 37, for example the number of sectors.","section":"Sec. 2"},{"comment":"Table 1 lists two solutions for TIC 406749309; marking the adopted solution more prominently (bold row or explicit footnote) would avoid confusion.","section":"Table 1"},{"comment":"The argument about the M2 > M1 systems would benefit from a one-line derivation using Eq. (5) to show why decreasing M1 cannot reverse the mass ratio for these systems.","section":"Sec. 3"}],"recommendation":"major_revision","confidential_remarks":"The catalog and fits are careful and useful. The missing completeness correction is, in my view, the single blocking issue; it is fixable by injection-recovery tests or by reframing the paper as a detection catalog. I would support publication after that work. The paper is within the scope of the journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a good catalog paper that adds 36 new heartbeat stars from APOGEE binaries, with careful light-curve fits, eclipse masks, TEO flags, and eight mass functions. The RV analysis is thoughtful—period grid search, phase-coverage cuts, and the HD 164816 SB2 solution gives a genuinely interesting massive pair. If you work on tidal binaries, this is a useful sample.\n\nThe soft spot is the population claim. The abstract says they confirm that the HB fraction rises rapidly with effective temperature and that non-giant primaries have started to evolve off the MS. But their own Sec 1 explains that an amplitude-limited search will naturally produce a decreasing HB frequency at lower temperatures, because cooler convective envelopes damp tides and slower evolution yields smaller amplitudes. Fig 9 plots the detected fraction versus color with no completeness correction. So the 'confirmation' is partly a restatement of the selection function. The same bias applies to Fig 8: more evolved, larger-radius stars generate larger amplitudes and are preferentially detected, so the displacement above the equal-mass isochrone is not clean evidence about the parent population.\n\nI don't think this sinks the paper. The catalog is solid, and the physical trend could well be real—the agreement between the APOGEE and Gaia samples helps, since the parent samples differ. But the authors should either make a serious attempt at a completeness correction or explicitly restate the conclusions as describing the detected fraction, not the underlying binary population. The two M2>M1 systems are handled with a qualitative argument; it's plausible but not rigorous, and the HD 164816 masses depend on a single light-curve inclination, so treat those numbers as model-dependent.\n\nWho is this for? Observers working on heartbeat stars, tidal physics, and binary evolution. It deserves a serious referee, and the referee report should push on the completeness question. It's a conditional accept, not a reject.","headline":"A solid catalog paper with 36 new heartbeat stars and careful fits, but the headline population trend is not yet separated from the search's own amplitude selection, which the authors themselves describe in Sec 1.","tokens_in":19217,"tokens_out":3516,"would_cite":true,"duration_ms":34692,"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":"From 31,548 binaries, 50 heartbeat stars are found, confirming that heartbeat stars are a hot-star phenomenon that disappears near the Kraft break.","keywords":["heartbeat stars","eccentric binaries","tidal interaction","main-sequence binaries","TESS","APOGEE","tidally excited oscillations","Kraft break"],"falsifier":"Measure the heartbeat fraction in a sample selected without an amplitude threshold—for example by stacking long TESS baselines or using eclipse-detected heartbeat stars, whose detectability does not depend on tidal amplitude—and check whether the roughly 1% fraction for the bluest binaries survives or flattens toward the Kraft break. A flat fraction across colour would falsify the claim that cool binaries are intrinsically less likely to be heartbeat stars.","tokens_in":18093,"feed_emoji":"💓","tokens_out":9456,"duration_ms":75170,"temperature":0.7,"pith_summary":"The paper searches TESS light curves of 31,548 main-sequence binaries flagged as binaries by APOGEE radial-velocity scatter and identifies 50 heartbeat stars—short-period, eccentric binaries whose stars are tidally deformed near periastron. Its central aim is to test two claims from the authors' earlier work: that heartbeat primaries have started to evolve off the main sequence, and that the fraction of main-sequence binaries that are heartbeat stars rises sharply with effective temperature. Both claims are confirmed: roughly 1% of binaries with $B_P-R_P<0.5$ mag are heartbeat stars, and the fraction drops steeply toward the Kraft break, where stellar envelopes become convective. If the result holds, the visibility of heartbeat stars is set by stellar structure and evolution speed, not by binary architecture alone.","feed_headline":"50 new heartbeat stars confirm the hot-star trend","feed_subtitle":"The fraction of binaries with tidal heartbeat signals hits about 1% for blue stars, then vanishes at the Kraft break.","key_machinery":"The load-bearing object is the analytic tidal light-curve model of Kumar, Ao, and Quataert (1995), which writes the flux as $F = Z + S\\,[1 - 3\\sin^2 i\\,\\sin^2(\\nu+\\omega)]/(1-e\\cos E)^3$, where $S$ sets the amplitude, $Z$ the mean flux, $i$ the inclination, $\\omega$ the argument of periastron, $e$ the eccentricity, and $\\nu$ and $E$ the true and eccentric anomalies. The search pipeline uses a Lomb–Scargle periodogram, phase-dispersion minimization, and a $\\chi^2$ ratio against a linear fit, followed by MCMC fitting of the phase-folded light curves with eclipses masked. For the radial velocities, the light-curve values of $e$ and $\\omega$ are held fixed while the velocity semi-amplitude $K$ and centre-of-mass velocity are fitted, and the binary mass function $f(M) = P K^3 (1-e^2)^{3/2}/2\\pi G$ turns $K$ into secondary-mass estimates via StarHorse primary masses. This machinery is what lets a sparse set of APOGEE epochs and TESS photometry yield both orbital elements and enough masses to test population trends.","core_discovery":"Using the APOGEE DR17/19 main-sequence binaries with $T<13.5$, the paper constructs a homogeneous sample of 50 heartbeat stars with orbital periods from 1.48 to 10.92 days and eccentricities from $e=0.085$ to $0.588$. Fourteen of these show eclipses, fifteen show tidally excited oscillations, and eight have enough non-rejected APOGEE radial-velocity epochs to measure binary mass functions and estimate secondary masses. The population-level discovery is that binned heartbeat stars sit at or above the equal-mass binary isochrone on the Gaia colour–magnitude diagram, which the authors read as evidence that the primaries have begun to evolve off the main sequence, and that the fraction of binaries that are heartbeat stars is highest for the bluest stars, reaching about 1% for binaries with $B_P-R_P<0.5$ mag and declining rapidly toward the Kraft break at $T_{\\rm eff}=6550$ K.","pith_inferences":["The paper applies no completeness correction to the colour–fraction plot, so an alternative reading—which the authors do not make—is that the trend tracks the amplitude threshold of the search rather than the intrinsic abundance of heartbeat stars.","If the temperature trend is physical, it predicts that the same binaries observed with longer TESS baselines or in a search with a uniform amplitude limit will show a smoother, less abrupt decline across the Kraft break than Fig. 9 shows.","The two systems with $M_2>M_1$ and the two very massive stars in TIC 406749309 depend on StarHorse primary masses that include some secondary light; future double-lined spectroscopic orbits for these eight systems could confirm or revise the companion masses.","TESS Cycle 8's longer 54-day sectors should reveal heartbeat stars with periods above 13 days, which would test whether the sharp boundary in the period–eccentricity diagram is a detection effect or a tidal-circularization boundary."],"forward_implications":["Heartbeat-star primaries on the upper main sequence are generally near the end of their main-sequence life: the brightness offset above the single-star isochrone persists even though the measured mass functions show a range of mass ratios, not just equal-mass twins.","The drop in heartbeat fraction toward the Kraft break indicates that convective envelopes damp the tidal variability, so cooler stars of the same orbital architecture should produce smaller, harder-to-detect light-curve amplitudes.","Any amplitude-limited heartbeat search will systematically undercount cool main-sequence binaries, so the true occurrence of heartbeat stars below the Kraft break is higher than any current census reports.","Thirty-six newly identified heartbeat stars, including fifteen systems with tidally excited oscillations, extend the sample available for studying dynamical tides, orbital circularization, and companion demographics."],"supporting_citations":[{"why":"The predecessor study that supplies the search methodology, the candidate region, and the two main-sequence trends this paper sets out to confirm.","marker":"Callahan et al. (2025)"},{"why":"The analytic tidal light-curve model used to fit all phase-folded TESS light curves and extract orbital elements.","marker":"Kumar et al. (1995)"},{"why":"Establishes the radial-velocity scatter criterion that selects the binary parent sample from APOGEE.","marker":"Badenes et al. (2018)"},{"why":"Defines the Kraft break separating convective from radiative envelopes, the temperature boundary against which the heartbeat fraction declines.","marker":"Kraft 1967"},{"why":"Provides the evolutionary-rate-versus-damping argument used to explain why amplitude-limited searches find fewer heartbeat stars at low temperature.","marker":"MacLeod & Loeb (2025)"},{"why":"Supplies the theoretical period–eccentricity envelope that all detected heartbeat stars are checked against for tidal circularization.","marker":"Mazeh 2008"},{"why":"The archival spectroscopic orbit for TIC 406749309 used to select the correct solution among the two degenerate light-curve fits.","marker":"Trepl et al. (2012)"}],"fun_headline_variants":["TESS and APOGEE reveal 50 heartbeat stars in hot binaries","Heartbeat stars prefer hot binaries: 50 new examples","Hot binaries more likely to show heartbeat signals: 50 found","50 heartbeat stars: hottest binaries beat the trend"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The colour trend in Fig. 9 is interpreted as a property of the underlying binary population, which requires that the heartbeat detection efficiency does not vary with colour; the paper itself states that an amplitude-limited search should find fewer cool stars even if they occur at the same rate.","fun_headline_variants_meta":{"raw":{"variants":["TESS and APOGEE reveal 50 heartbeat stars in hot binaries","Heartbeat stars prefer hot binaries: 50 new examples","Hot binaries more likely to show heartbeat signals: 50 found","50 heartbeat stars: hottest binaries beat the trend"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00095,"raw_usage":{"total_tokens":4022,"prompt_tokens":880,"completion_tokens":3142,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":496,"completion_tokens_details":{"reasoning_tokens":3073}},"tokens_in":496,"tokens_out":3142,"duration_ms":21810,"temperature":1.0,"reasoning_tokens":3073,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:08:05.476656+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the heartbeat fraction in a sample selected without an amplitude threshold—for example by stacking long TESS baselines or using eclipse-detected heartbeat stars, whose detectability does not depend on tidal amplitude—and check whether the roughly 1% fraction for the bluest binaries survives or flattens toward the Kraft break. A flat fraction across colour would falsify the claim that cool binaries are intrinsically less likely to be heartbeat stars.","supporting_citations":[],"review_version":1}