{"id":"bae6a7b6-f17f-4075-bd8d-7e8a62c4854c","arxiv_id":"2504.19866","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Simulations predict that detectable two-star pulsators have near-equal initial masses, that double red-clump binaries are absent below about 500 solar radii, and that roughly one percent of red giants have gained or lost significant mass.","lead":"This paper uses computer simulations of the Kepler field to predict when binary stars show pulsation signals from both components, and which binary histories leave detectable traces. It gives observers concrete expectations for what asteroseismology should find in current and future space-based surveys.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"","rationale":"","tokens_in":26179,"tokens_out":9185,"duration_ms":89864,"concrete_test":"","verdict_should_be":"UNCHANGED","load_bearing_attack":"","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses the TRILEGAL population synthesis code with its BinaPSE module to simulate the Kepler field, comparing three prescriptions for initial binary parameters: Eggleton (2006), Moe & Di Stefano (2017), and a non-interacting binary model. It applies published models for the detectability of solar-like oscillations in single stars (Chaplin et al. 2011b) and in unresolved binaries (Miglio et al. 2014), and predicts the occurrence rates, mass ratios, orbital properties, and evolutionary states of asteroseismic binaries and of products of binary interactions (under-massive and over-massive stars). The central claims are that asteroseismic binaries are near-equal-mass systems that have not exchanged mass, that double red-clump binaries should be absent at separations below roughly 500 R_sun, and that about 1% of Kepler red giants with detectable oscillations may have experienced significant mass accretion or loss.","tokens_in":26076,"tokens_out":10693,"duration_ms":92743,"significance":"If the predictions are reliable, this paper provides useful quantitative expectations for the ongoing and future systematic searches for asteroseismic binaries in Kepler, K2, TESS, and PLATO data. Its strengths are the use of a well-established population synthesis code, the explicit comparison of three initial-condition prescriptions, the explicit treatment of detection probability for both components of unresolved binaries, and several clearly stated caveats (target selection, tidal mode suppression, uncertainty in the binary fraction). The main predictions are falsifiable with existing and upcoming data, which makes the paper a valuable benchmark even if the exact rates are model-dependent. However, the quantitative claims need to be internally consistent and properly quantified before they can serve as robust expectations.","major_comments":[{"comment":"The derivation of the 500 R_sun boundary for double red-clump binaries is arithmetically inconsistent. With r1/A = 0.38, a separation of A = 500 R_sun gives r1 = 190 R_sun, not the stated r1 ≲ 150 R_sun. If the relevant RGB-tip radius is about 150 R_sun, the same equation yields a boundary of A ≈ 395 R_sun. The claimed threshold of 500 R_sun in the abstract, Fig. 6, and Conclusions therefore does not follow from Eq. (9) as written; please correct the threshold or the radius estimate and update the associated claims.","section":"Section 3.5, Eq. (9)"},{"comment":"The reported occurrence rates are internally inconsistent with the numbers in Tables 2 and 3. The text states that for every 1000 red giant stars with detectable oscillations, about 140 are in binary systems and 2.7 are in asteroseismic binaries, but the quoted counts (11544 in binaries and 227 in asteroseismic binaries out of 41537) correspond to about 278 and 5.5 per 1000, respectively. Furthermore, the abstract's claim of a minimum fraction of 0.06% for the Eggleton prescription does not match Table 6, where summing the Nseismo entries for RGB, CHeB, and EAGB gives 27 systems (54 stars) against roughly 41,800 detectable red giants, i.e., about 0.13%. Please reconcile the abstract and text with the tables.","section":"Section 3.3 and Abstract"},{"comment":"The occurrence rates are reported without any uncertainties. Many of the key counts are small (e.g., Nseismo = 53 for double CHeB in MDS17, 11 in E06, 7 for double RGB in E06), so Poisson fluctuations are substantial. Without error bars or confidence intervals, the comparisons among prescriptions in Fig. 13 and the assertions that one prescription yields the 'highest' or 'lowest' fraction are not statistically robust. Please add at least Poisson confidence intervals and explicitly discuss the small-number limitations in the affected rates.","section":"Tables 3, 6, and A.1-A.3"}],"minor_comments":[{"comment":"In the discussion of KIC 9246715, 'expected Roche lobe r1≈ 80 M⊙' should read '80 R⊙'; the units are incorrect as printed.","section":"Section 3.5"},{"comment":"The formula for the initial orbital period appears garbled in the typeset version; please check the placement of the exponent and parentheses in Eq. (1).","section":"Section 2.3, Eq. (1)"},{"comment":"The note defining Nseismo_pri and Nseismo_sec refers to 'Nseismo_pri' twice in a way that is confusing; please clarify which column corresponds to primaries and which to secondaries.","section":"Table 2, note"},{"comment":"The comparison of the simulated number of red giants with detectable oscillations (~42,000) to the observed number (~30,000) cites 'Garcia et al., in preparation'; please provide a more specific reference or at least a description of the data source.","section":"Section 3.2"},{"comment":"Please clarify in the text and figure caption that the 500 R_sun boundary is based on the initial semi-major axis, since the discussion elsewhere sometimes refers to 'orbital separation' without specifying whether initial or final values are meant.","section":"Figure 6 and Section 3.5"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a timely question and uses appropriate tools, but the internal numerical inconsistencies and the arithmetic error in the 500 R_sun derivation make the current version unsuitable for publication without substantial revision. The issues are localized and should be fixable; the central simulation approach and the scientific questions are sound. I would be willing to review a revised version that reconciles the tables with the abstract and text and provides uncertainties for the small-number rates."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is the first paper to run interacting binary evolution in TRILEGAL/BinaPSE to set quantitative expectations for asteroseismic binaries, and it will be a useful reference when the companion search paper lands. The central claims are honestly derived and mostly robust; the specific rates are model-dependent and some small-number tables need error bars, but the paper deserves a serious referee.\n\nWhat's new: Miglio et al. (2014) used non-interacting binaries; here they compare Eggleton (2006) and Moe & Di Stefano (2017) initial distributions, include mass transfer, CE, mergers, and give rates per 1000 detectable red giants for double-RGB, RGB+CHeB, double-CHeB, and under/over-massive products. The sharpest output is the prediction that low-mass double red-clump binaries should not exist below ~500 R_sun because the primary overflows its Roche lobe near the RGB tip. That is a clean, falsifiable statement, and it is supported by a simple analytic argument (r_L ~ 0.38A vs RGB-tip radius ~150 R_sun), not just by the simulation.\n\nThey are also honest about limitations: detectability ignores tidal mode suppression (they quantify the affected systems—small numbers), blending, and the Kepler target selection. They report that their fiducial MDS17 simulation underpredicts observed over-massive giant fractions (0.7% vs ~5-18% in the literature) and say so.\n\nSoft spots, in proportion:\n\n- No data/code release. The simulations are 'available upon reasonable request.' For a paper whose main product is a target list, that dampens the value.\n- Several headline counts are tiny (e.g., N=4 double-EAGB systems, N=65 double-CHeB) with no Poisson errors. The rates like 3 per 1000 CHeB are probably a factor of ~2 uncertain, and that should be shown.\n- The 'require q near 1' conclusion is partly a selection effect of the detectability model (flux dilution) and the input q distributions, which include a twin excess in MDS17. Worth discussing what is physical vs baked in.\n- The double-RC boundary depends on BSE's RLOF/CE stability criteria. It is a prediction, not a measurement; the authors frame it correctly, but the 500 R_sun number could move if mass transfer is more conservative.\n\nOn circularity: not a concern. The rates come from external initial-condition prescriptions and published detectability models, not fitted to the target results.\n\nBottom line: solid simulation-expectations paper, no load-bearing flaw. I'd send it to review. My main asks would be Poisson errors on small counts, a clear separation of robust conclusions from prescription-dependent rates, and at least a data table in a public repo. Anyone planning an asteroseismic binary search in Kepler/TESS/PLATO data will want this paper.","headline":"Solid expectations paper with sharp falsifiable predictions; rates are prescription-dependent, but the near-equal-mass and no-short-period-double-RC conclusions are worth taking seriously.","tokens_in":26663,"tokens_out":4985,"would_cite":true,"duration_ms":48847,"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":"Asteroseismic binaries are near-equal-mass systems that never exchanged mass, and double red-clump binaries should be absent at separations below about 500 solar radii.","keywords":["asteroseismic binaries","solar-like oscillations","binary population synthesis","red clump stars","mass transfer","Kepler field","red giant oscillations","age-metallicity relation"],"falsifier":"A confirmed double red-clump binary with total mass below about $4\\,M_\\odot$ and semi-major axis below about $500\\,R_\\odot$ (initial period below roughly 1000 days) would contradict the prediction; look for a low-mass core-helium-burning pair whose orbit is short enough that the primary would have filled its Roche lobe at the red-giant tip, with both components showing solar-like oscillations.","tokens_in":25957,"feed_emoji":"🔭","tokens_out":9057,"duration_ms":85977,"temperature":0.7,"pith_summary":"The paper asks which binary systems in the Milky Way can be seen as asteroseismic binaries — pairs where both stars show detectable solar-like oscillations — and what products of binary interactions look like in Kepler-style observations. By synthesizing the Kepler field with three prescriptions for initial binary properties, including interacting and non-interacting populations, it arrives at a clear expectation: systems detected in both components are born with nearly equal masses and never exchange mass, because even small mass-transfer events push one component's oscillation signal below detectability. The simulation also predicts that double red-clump binaries made of two low-mass giants should be absent at separations below about 500 $R_\\odot$, since the primary overflows its Roche lobe near the tip of the red giant branch. These predictions matter because they set quantitative baselines for interpreting Kepler and future PLATO data, and because the same simulations quantify how many field red giants are over- or under-massive products of binary evolution that can distort apparent age-metallicity relations.","feed_headline":"No double red-clump binaries expected under 500 solar radii","feed_subtitle":"Simulations show low-mass pairs that close must swap mass on the red giant branch and hide both oscillation signals.","key_machinery":"The argument is carried by a population-synthesis pipeline that combines TRILEGAL stellar populations, PARSEC evolutionary tracks, and the BinaPSE binary-evolution module, which applies analytic prescriptions from BSE for mass loss, Roche-lobe overflow, common-envelope evolution, tidal circularization, and mergers. On top of the resolved stars, the paper applies a probability of detecting solar-like oscillations from stellar parameters, and for unresolved binaries multiplies the two components' probabilities, $p_{\\mathrm{seismo,bin}} = p_{\\mathrm{seismo,pri}} \\times p_{\\mathrm{seismo,sec}}$, so that any reduction in one component's signal due to a companion suppresses the chance of seeing both. The key physical mechanism producing the sharpest prediction is Roche-lobe geometry at the red-giant tip: with the volume-equivalent Roche-lobe radius $r_1/A \\approx 0.38$ for equal-mass pairs, a low-mass binary with separation below about $500\\,R_\\odot$ will have its primary overflow its Roche lobe near the RGB tip, initiating mass transfer that destroys the double red-clump configuration.","core_discovery":"The central claim is that asteroseismic binaries are a clean, restricted subset of the binary population: systems with initial mass ratio close to unity that avoided Roche-lobe overflow and common-envelope episodes throughout their evolution. The paper argues this from three simulations of $121\\,\\mathrm{deg}^2$ of the Kepler field using the Moe & Di Stefano (2017, MDS17), Eggleton (2006, E06), and non-interacting prescriptions. In the MDS17 simulation, 95% of detected asteroseismic binaries have initial mass ratio above 0.96, and the close-to-one final mass ratio is preserved apart from single-star wind mass loss or tidal circularization. The consequence is that the occurrence rate of double-giant asteroseismic binaries is low — for red giants with detectable oscillations, the fraction ranges from 0.06% under E06 to 0.46% in the non-interacting case — and that short-period double red-clump binaries with total mass below $4\\,M_\\odot$ should not exist, because at separations below roughly $500\\,R_\\odot$ the more massive component fills its Roche lobe at the red-giant tip and the system evolves away from a double red-clump configuration.","pith_inferences":["A direct search for short-period double red-clump binaries in Kepler, TESS, or spectroscopic surveys would test the adopted Roche-lobe overflow stability criteria: finding one would require more conservative mass transfer than the BSE-style prescriptions assume.","The dilution argument implies that asteroseismic catalogs of binaries will be biased toward unevolved, near-twin systems, so interaction rates inferred from such catalogs alone would be lower limits.","The 500 $R_\\odot$ boundary could be turned into a quantitative probe of the red-giant-tip radius and RGB mass-loss efficiency, since the exact cutoff location encodes the maximum radius reached by low-mass stars before the helium flash.","The paper's distance-mismatch toy model suggests a practical identification pipeline: unresolved equal-luminosity binaries produce a characteristic about 29% underestimate of asteroseismic distances, which could flag candidate asteroseismic binaries in combined Gaia and seismic samples."],"forward_implications":["In Kepler-like data, double-giant asteroseismic binaries should be rare: roughly one double-RGB system per 1000 detectable RGB stars, and three double-CHeB systems per 1000 detectable CHeB stars under the MDS17 prescription.","Any observed double red-clump binary with combined mass below about $4\\,M_\\odot$ and separation below about $500\\,R_\\odot$ would directly challenge the adopted mass-transfer and common-envelope prescriptions.","Over- and under-massive giants make up roughly 1.5% of detectable red giants under MDS17; if not identified, these stars bias age estimates and add scatter to Galactic age-metallicity relations.","The predicted rates differ strongly between the E06 and MDS17 initial-binary prescriptions, so counting asteroseismic binaries is a clean way to constrain the initial binary fraction and period distribution.","Because asteroseismic binaries preserve their birth mass ratio and orbital separation, they give direct observational access to the initial mass-ratio and orbital-parameter distributions of the binary population."],"supporting_citations":[{"why":"Supplies the TRILEGAL population-synthesis machinery that generates the synthetic Kepler-field populations.","marker":"Girardi et al. 2005"},{"why":"Provides the BinaPSE module that evolves binary components and applies binary-interaction prescriptions on top of TRILEGAL tracks.","marker":"Dal Tio et al. 2021"},{"why":"Supplies one of the two initial binary parameter distributions, with correlated periods, mass ratios, and eccentricities.","marker":"Moe & Di Stefano 2017"},{"why":"Supplies the alternative Monte Carlo initial binary parameter distribution used for the comparison simulation.","marker":"Eggleton 2006"},{"why":"Provides the analytic BSE prescriptions for Roche-lobe overflow, common-envelope evolution, and mergers that BinaPSE adopts.","marker":"Hurley et al. 2002"},{"why":"Gives the probability-of-detection formula for solar-like oscillations from stellar parameters and Kepler magnitude.","marker":"Chaplin et al. 2011b"},{"why":"Provides the binary detection-probability model that accounts for dilution of the seismic signal by the companion.","marker":"Miglio et al. 2014"},{"why":"Supplies the Roche-lobe radius formula used to derive the roughly 500 $R_\\odot$ boundary for double red-clump binaries.","marker":"Paczyński (1971)"},{"why":"Supplies the stellar tracks showing red-giant-tip radii near 150 $R_\\odot$ for low-mass stars, establishing the Roche-lobe condition.","marker":"Nguyen et al. 2022"},{"why":"Supplies the PARSEC evolutionary tracks used to set stellar mass, radius, luminosity, and evolutionary state in the simulations.","marker":"Bressan et al. 2012"}],"fun_headline_variants":["Double red-clump binaries ruled out below 500 solar radii","Asteroseismic binaries demand near-equal initial masses","Mass transfer hides both oscillation signals in red giants","Interacting binaries never show double red-clump oscillations","Only near-unity mass ratios yield detectable double oscillators"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The short-period double red-clump prediction rests on the assumption that when a low-mass red giant fills its Roche lobe near the tip of the red giant branch, the resulting mass transfer disrupts the system so that it does not become a double red-clump binary; if the models allow more conservative mass transfer, those binaries could survive below 500 solar radii.","fun_headline_variants_meta":{"raw":{"variants":["Double red-clump binaries ruled out below 500 solar radii","Asteroseismic binaries demand near-equal initial masses","Mass transfer hides both oscillation signals in red giants","Interacting binaries never show double red-clump oscillations","Only near-unity mass ratios yield detectable double oscillators"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000287,"raw_usage":{"total_tokens":1761,"prompt_tokens":1098,"completion_tokens":663,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":714,"completion_tokens_details":{"reasoning_tokens":584}},"tokens_in":714,"tokens_out":663,"duration_ms":6614,"temperature":1.0,"reasoning_tokens":584,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:41:26.976432+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A confirmed double red-clump binary with total mass below about $4\\,M_\\odot$ and semi-major axis below about $500\\,R_\\odot$ (initial period below roughly 1000 days) would contradict the prediction; look for a low-mass core-helium-burning pair whose orbit is short enough that the primary would have filled its Roche lobe at the red-giant tip, with both components showing solar-like oscillations.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the TRILEGAL population-synthesis machinery that generates the synthetic Kepler-field populations."},{"cited_title":"2006, Evolutionary Processes in Binary and Multiple Stars (Cam- bridge University Press)","cited_arxiv_id":null,"evidence_quote":"Supplies the alternative Monte Carlo initial binary parameter distribution used for the comparison simulation."},{"cited_title":"J., Farmer, R., et al","cited_arxiv_id":null,"evidence_quote":"Provides the binary detection-probability model that accounts for dilution of the seismic signal by the companion."},{"cited_title":"T., Costa, G., Girardi, L., et al","cited_arxiv_id":null,"evidence_quote":"Supplies the stellar tracks showing red-giant-tip radii near 150 $R_\\odot$ for low-mass stars, establishing the Roche-lobe condition."}],"review_version":1}