{"id":"ec287db4-882e-476b-be61-674eed244cdd","arxiv_id":"2602.10187","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A simplified tidal-dissipation prescription implemented in COMPAS predicts equilibrium and dynamical tidal strengths up to orders of magnitude larger than the standard Z77 model for many binary stars.","lead":"This paper fits modern tidal physics into the rapid binary-population code COMPAS, replacing older crude tidal recipes with prescriptions that depend on stellar structure and tidal frequency. The new model predicts tidal dissipation that is often 1–7 orders of magnitude stronger than the standard Zahn-1977 recipe, changing which binaries circularize and synchronize.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Always-on IGW dissipation likely biases the 1–7 order-of-magnitude dynamical tide claim for radiative-core stars; a period-resolved test is needed.","rationale":"The reader's weakest_assumption identifies the always-on IGW dissipation as the load-bearing concern. I agree: the abstract's strongest claim is the 1–7 order-of-magnitude dynamical tide enhancement, which for solar-type stars is driven by IGW (and later IW) dissipation. §2.2.2 explicitly states the wave-breaking condition is omitted. The paper's own text provides the ammunition: the critical period for solar-type stars is ~3 days on MS, yet the illustrative 1 M⊙ binary starts at 10 d and the paper claims IGW is dominant. That is a direct internal tension worth testing. I considered other possible concerns: (1) the equilibrium tide mixing-length overestimate by factor ~5 (acknowledged, and not central to the dynamical tide claim); (2) the fitted spin evolution in Appendix A.2 (a validation weakness but not the central claim); (3) the low-e secular equations used at e=0.5 (acknowledged, and the paper's qualitative interpretation is guarded). The IGW always-on assumption is the most load-bearing because it directly inflates the headline dynamical tide enhancement. The concrete test I propose is a minimal sensitivity run: turn off IGW dissipation when the wave-breaking criterion is not met in the 1 M⊙ example. If the dynamical tide dominance and circularization behavior survive, the claim holds; if not, the abstract's generalization is overstated. The verdict remains CONDITIONAL because the paper is a methods paper with acknowledged approximations; a targeted sensitivity test would either confirm or require qualification, but does not warrant rejection of the framework.","tokens_in":28710,"tokens_out":1856,"duration_ms":18286,"concrete_test":"Re-run the 1 M⊙ + 1 M⊙ P_orb=10 d, e=0.5 COMPAS simulation with a modified IGW dissipation term that multiplies Im[k] by an efficiency factor f(omega_t, L_tide) = 0 when the orbital period exceeds the critical wave-breaking period P_crit(t) from Esseldeurs et al. (2024) Fig. 4 for a 1 M⊙ star, and f=1 below P_crit. If the final circularization time and spin evolution change by less than a factor of ~2, the always-on assumption is safe for this claim; if the dynamical tide enhancement over Z77 drops by an order of magnitude or more, the abstract's 1–7 order-of-magnitude claim needs qualification.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's headline claim that dynamical tides can be 1–7 orders of magnitude stronger than Z77 depends critically on the assumption in §2.2.2 that IGW dissipation is 'always efficient' (wave breaking always achieved), and §2.4 further ignores IGW dissipation in stars with >2 layers. This is not a minor caveat: the 1 M⊙ + 1 M⊙ example (Fig. 5) shows dynamical tides exceeding equilibrium tides by 1–2 orders of magnitude, and the abstract generalizes to 'low mass main sequence and giant type stars.' If the always-on assumption fails, the central comparison against Z77 is systematically overstated, because Z77's dynamical tide is also IGW-based but with a different structure dependence. The paper itself notes the critical orbital period for solar-type stars is ~3 days on MS and decreases with evolution (citing Esseldeurs et al. 2024), yet the simulated 1 M⊙ binary starts at P=10 days and the claim is that IGW dissipation dominates. This appears internally inconsistent: at 10 days, the wave-breaking condition may not be met. Moreover, §2.4 asserts IGW dissipation is 'typically sub-dominant to equilibrium or even IW tides' while Fig. 5 shows it is dominant for most of the evolution. The reader's concern is correct. The concrete test: reproduce Fig. 5 with an IGW efficiency factor that switches off below the critical amplitude/period (e.g., Barker & Ogilvie 2010 criterion, or Esseldeurs et al. 2024 period-dependent efficiency) and measure the resulting change in circularization/synchronization.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a fast, closed-form tidal dissipation prescription for the binary population synthesis code COMPAS. It models equilibrium tides through a frequency-dependent turbulent viscosity in convective envelopes, dynamical tides through internal gravity waves (IGWs) excited at convective-radiative boundaries, and inertial-wave (IW) dissipation in convective envelopes. The authors implement the prescriptions, validate them against a solar-type equilibrium-tide comparison and a Sun+hot-Jupiter dynamical-tide comparison, and evolve four representative binaries (0.3+0.3 Msun, 1+1 Msun, 2.5+2.5 Msun, 3+3 Msun) against the Z77/Hurley reference model. The central claim is that this simple model retains the dominant stellar-structure and frequency dependence of tidal dissipation, agrees with detailed simulations to within an order of magnitude, and predicts equilibrium and dynamical tidal strengths that can exceed commonly used prescriptions by 1-2 and 1-7 orders of magnitude, respectively.","tokens_in":29185,"tokens_out":7668,"duration_ms":85560,"significance":"If the claims hold, this would be a valuable contribution: it would give population synthesis a structure- and frequency-aware tidal module, bridging the gap between expensive MESA/STAREVOL-based tidal calculations and million-system binary population studies. The paper is transparent about its approximations and ships explicit equations and COMPAS integration, which is a strength. However, the headline quantitative claims are currently stronger than the evidence supports. The most load-bearing issue is the \"always efficient\" IGW dissipation assumption, which is internally inconsistent with the quoted critical-period criterion and directly controls the claimed 1-7 order-of-magnitude enhancement for radiative-core stars. The dynamical-tide validation in Appendix A.2 is also partly circular, and the high-eccentricity application of low-eccentricity secular equations needs quantitative justification. These issues are fixable, so the paper warrants a major revision rather than rejection.","major_comments":[{"comment":"The assumption that IGW dissipation is \"always efficient\" contradicts the manuscript's own critical-period statement. §2.2.2 states that the critical orbital period for solar-type stars is about 3 days on the MS and decreases with evolution, and that this critical period \"should always be lower than even the shortest period binaries we consider.\" The 1 Msun + 1 Msun example in §3.2 starts at P_orb=10 days. If P_orb > P_crit, the high-amplitude wave-breaking condition used to justify efficient IGW damping is not satisfied. Yet Fig. 5 shows IGW dissipation dominating the early evolution and the abstract generalizes the 1-7 order-of-magnitude dynamical-tide enhancement to low-mass MS and giant stars. This is a load-bearing internal inconsistency, not just a caveat. I request a period-resolved sensitivity test: for example, switch off IGW dissipation when P_orb exceeds the quoted P_crit, or","section":"§2.2.2, §2.4, §3.2, Fig. 5"},{"comment":"The 6-7 orders-of-magnitude dynamical-tide enhancement is measured against a Z77 model that the authors themselves describe as \"technically irrelevant\" for radiative-core stars with convective envelopes. Comparing a new model to an inapplicable reference does not establish a physical enhancement; it only documents the difference between two different model domains. The abstract's \"1-7 orders of magnitude stronger\" statement therefore overstates the physical significance. I recommend reframing the comparison: either compare against a modern, applicable prescription (e.g., Goodman & Dickson 1998, Terquem et al. 1998, Ahuir et al. 2021, Esseldeurs et al. 2024) or clearly label the Z77 comparison as an inapplicable benchmark and report the actual enhancement relative to an applicable modern model.","section":"§3.2, Figs. 4-5, Abstract"},{"comment":"The dynamical-tide validation is circular in an important way. The spin evolution is fitted piecewise-linearly to Fig. 7 of Ahuir et al. (2021) via Eq. (A3), and the resulting Love numbers are then compared with Fig. 9 of the same paper. Since the tidal frequency entering Im[k] is determined by the fitted spin, this comparison largely confirms that the code can reproduce the input spin history, not that the dissipation model is independently predictive. This is especially problematic because the target spin evolution in Ahuir et al. is itself affected by tides. I request an independent validation: self-consistently evolve the spin with the proposed tidal model and compare the resulting tidal Love numbers and spin evolution, or use a different dataset/observable that does not rely on fitting the target's spin.","section":"Appendix A.2, Eq. (A3)"},{"comment":"The secular tidal equations are formally truncated at O(e^2), yet they are applied at e_ZAMS=0.5 in the headline 1 Msun + 1 Msun example and in the high-eccentricity grid points of Fig. 4. The implementation further drops the O(e^2) terms as a numerical stopgap when spins exceed pseudo-synchronization, which breaks angular-momentum conservation in that regime. The manuscript acknowledges the issue, but the e=0.5 cases are central to the paper's qualitative conclusions, and the top panel of Fig. 4 shows binaries becoming wider and more eccentric, a behavior that may be an artifact of the truncated equations. I request a quantitative assessment of this systematic error, for example by comparing with a calculation that retains higher-order e terms (or with a smaller-e control case) for at least one of the fiducial binaries.","section":"§2.4, Eqs. (11)-(13), §3.2, Fig. 4"}],"minor_comments":[{"comment":"LaTeX encoding artifacts appear in Brunt-Väisälä and in several author names (e.g., \"V¨ais¨al¨a\"); please fix the source to render correctly.","section":"Throughout"},{"comment":"The text says the semi-latus rectum a(1-e^2) \"should remain constant\" under angular momentum conservation, but when angular momentum is exchanged with stellar spins, the orbital angular momentum—and hence the semi-latus rectum—is not constant. The subsequent comparison actually shows a small decrease, consistent with spin angular momentum gain. Please rephrase to avoid this error.","section":"§3.2, p. 14"},{"comment":"The phrase \"the critical period should always be lower than even the shortest period binaries we consider\" is ambiguous and, as written, undermines the immediately following assumption of always-efficient IGW dissipation. Please state explicitly the direction of the inequality and how it justifies (or fails to justify) the modeling choice.","section":"§2.2.2"},{"comment":"The paper acknowledges numerical artifacts after 0.6×10^10 yr in panels (c) and (d). Please consider masking or clearly marking these time ranges in the figure, since they are visually prominent and could be mistaken for physical features.","section":"Fig. 3"},{"comment":"The comparison to observed circularization periods for low-mass binaries is brief. Given that the paper's equilibrium-tide model still underpredicts empirical circularization periods, a more explicit statement about which ingredient (frequency dependence, PMS tides, or missing physics) is most likely responsible would help the reader.","section":"§3.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is a useful methods contribution, and I do not think rejection is warranted if the authors can address the IGW efficiency issue and the validation circularity. The requested sensitivity tests and re-benchmarking are within the scope of a revision. Please note also that the phrase \"Data and configuration files for our simulations are available upon reasonable request\" falls short of current reproducibility expectations for a methods paper; I would encourage the authors to archive the configuration files and scripts in a public repository."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper is worth engaging with. It translates recent tidal dissipation theory into a practical module for COMPAS, and it documents the reference Z77 model with unusual care, which is a service to the field. The new content is the synthesis: frequency-dependent equilibrium tides from Duguid, IGW and inertial-wave parameterizations adapted for rapid synthesis, and the self-consistent coupling in COMPAS with explicit side-by-side comparisons to Z77. The qualitative conclusion that structure-dependent tides can differ from Z77 by orders of magnitude is plausible and supported by the cited literature.\n\nThe soft spots are real but manageable. The biggest is the always-on IGW dissipation assumption. The stress-test is right: for the 1 Msun binary at P=10 days, the paper itself cites a critical wave-breaking period around 3 days, so the dominance of IGW dissipation in Fig. 5 is likely overestimated. The internal tension with Section 2.4's statement that IGW dissipation is \"typically sub-dominant\" should also be cleaned up. This needs a sensitivity test with a period-dependent efficiency factor, not just an acknowledgment.\n\nSecond, the validation in Appendix A.2 fits the spin evolution to Ahuir et al. and then compares Love numbers to the same paper; that is partially circular and should be replaced or supplemented with an independent benchmark. Third, the acknowledged order-of-magnitude approximations (constant-density envelopes, guessed mixing length, low-eccentricity secular equations used at e=0.5) are not quantified with error bars. Fourth, no code or configs are shipped; \"available upon reasonable request\" makes reproduction harder than it should be for a methods paper.\n\nThe central method is not circular: the formulas follow from cited theory, not from fitting to the target claim. The Z77 reference implementation alone is worth having, and the authors are transparent about most of their choices. This is a paper for population synthesis practitioners and anyone modeling tides in binaries who wants a bridge between MESA-type detailed calculations and rapid synthesis.\n\nMy recommendation: send it to peer review. The issues are addressable in revision, and the paper's contribution is important enough to deserve referee time. I would push for the wave-breaking sensitivity test and for public code before acceptance.","headline":"Solid, useful methods paper that brings modern tidal theory into COMPAS; the always-on IGW dissipation is a real caveat but not a deal-breaker for the methods contribution.","tokens_in":29618,"tokens_out":2309,"would_cite":true,"duration_ms":27531,"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":"A structure-aware tidal prescription for rapid binary population synthesis predicts equilibrium tides 1–2 orders of magnitude and dynamical tides up to 7 orders of magnitude stronger than standard prescriptions.","keywords":["binary population synthesis","tidal dissipation","tidal Love number","equilibrium tides","dynamical tides","internal gravity waves","inertial waves","binary circularization"],"falsifier":"A census of circularization periods: the model predicts solar-type binaries circularize out to about 10-day orbits on the main sequence, while standard models predict almost none; a cluster sample showing circularization only below about 2–3 days would rule out the enhancement. Quantitatively, the model gives a tidal quality factor of about 8.8 × 10^9 for a Sun-like star at a 1-day tidal period, so any empirical measurement from heartbeat-star apsidal motion or hot-Jupiter decay that sits well above this would contradict the prescription's strength—the paper itself notes equilibrium tides can","tokens_in":1790,"feed_emoji":"⭐","tokens_out":3887,"duration_ms":166078,"temperature":0.7,"pith_summary":"The paper argues that the full frequency- and structure-dependence of tidal dissipation can be compressed into fast, closed-form expressions suitable for rapid binary population synthesis. Its central claim is that a few stellar-structure quantities—convective-envelope depth and density, convective-core radius, buoyancy frequency at radiative-convective boundaries—together with tidal frequency determine the dissipative tidal response, so the code can evolve millions of binaries with tidal physics faithful to detailed simulations. If true, previous population-synthesis studies have systematically underestimated tidal strength by one to seven orders of magnitude for exactly the stars that matter for binary circularization, synchronization, mass-transfer, and gravitational-wave source spins. The paper demonstrates order-of-magnitude agreement with detailed simulations across mechanisms, and shows that solar-type binaries can circularize at periods up to about 10 days while giant-branch binaries can circularize out to roughly 5000-day orbits—both far beyond what standard models achieve.","feed_headline":"New tidal model: binary star tides up to 10 million times stronger","feed_subtitle":"Standard population synthesis has underestimated tides in solar-type and giant binaries, changing circularization and spin predictions.","key_machinery":"The load-bearing object is the imaginary tidal Love number Im[k_l,n^m(ω_t)]—the dissipative part of a star's gravitational response to the tidal potential—which sets tidal power and torque, and therefore the rates of change of semi-major axis, eccentricity, and spin. The paper feeds this object with three closed-form channels: a constant-density convective envelope with a frequency-dependent turbulent viscosity for equilibrium tides; an E2 excitation factor times the 8/3 power of a dimensionless tidal frequency for internal gravity waves, with separate E2 factors for convective-core and radiative-core stars; and a frequency-averaged inertial-wave formula that turns on when the orbital freque","core_discovery":"The paper's central discovery is that tidal dissipation in binary population synthesis can be encoded in a frequency- and structure-dependent imaginary tidal Love number without sacrificing computational speed. For equilibrium tides, the model replaces a constant tidal-friction efficiency with a convective envelope of estimated density and a frequency-dependent turbulent viscosity. For dynamical tides from internal gravity waves, it uses an excitation factor derived from the stellar structure—scaling with the ninth power of the convective-core radius for convective-core stars, and with the buoyancy frequency at the radiative-convective boundary for radiative-core stars. For inertial waves, i","pith_inferences":["If the paper is right, the observed circularization of solar-type cluster binaries—which equilibrium-tide-only models cannot explain—may be accounted for by the newly strong dynamical tides, without invoking anomalously low tidal quality factors.","A testable extension is to implement the wave-breaking condition that the paper omits: switching off internal-gravity-wave tides above the critical orbital period would produce a sharp drop in circularization efficiency near a few days for solar-type stars, observable in cluster period–eccentricity distributions.","The same structure-dependent Love numbers carry over to exoplanet systems, so hot-Jupiter tidal decay and the ages of circularized planet-hosting stars could be re-derived with structure-aware tides instead of a fixed stellar quality factor.","If the enhancement survives, compact-object binaries inherit it: black-hole and neutron-star progenitors that circularize and synchronize earlier will arrive at compact-binary formation with different spins and separations, changing population predictions for gravitational-wave sources."],"forward_implications":["Solar-type binaries in rapid synthesis will now circularize within their main-sequence lifetimes at initial periods up to about 10 days and eccentricities below about 0.6, where standard recipes predict essentially no tidal evolution.","Equilibrium tidal strength varies by 1–2 orders of magnitude over a star's lifetime as its convective envelope evolves, with a characteristic boost as the binary approaches synchronization.","For convective-core stars, tidal strength decays sharply as the core shrinks on the main sequence (roughly as (R_c/R*)^9), so mass-only excitation coefficients overestimate late-main-sequence dynamical tides.","Giant-branch and AGB binaries circularize at substantially larger orbital periods (about 5000 days) than older models imply (about 3000 days), narrowing but not closing the gap with observed circularization periods.","Because the same structure-dependent Love numbers describe tides in planet-host stars, predictions for hot-Jupiter orbital decay and exoplanet circularization can be re-derived without an arbitrary constant stellar quality factor."],"fun_headline_variants":["Binary tides may be millions of times stronger than thought","New tidal model shows binary tides orders of magnitude stronger","Rapid binary synthesis: tides up to 10 million times stronger","Tidal forces in binary stars underestimated by orders of magnitude","New tidal model shows binary tides underestimated by orders of magnitude"],"cache_read_input_tokens":30720,"weakest_assumption_plain":"The load-bearing premise is that internal gravity waves always dissipate fully at every radiative-convective boundary: the paper omits the wave-breaking condition and drops gravity-wave tides entirely for stars with more than two structural layers, so if real stars quench these waves at low amplitudes or long orbital periods, the claimed dynamical-tide enhancements of 1–7 orders of magnitude are systematically overestimated.","fun_headline_variants_meta":{"raw":{"variants":["Binary tides may be millions of times stronger than thought","New tidal model shows binary tides orders of magnitude stronger","Rapid binary synthesis: tides up to 10 million times stronger","Tidal forces in binary stars underestimated by orders of magnitude","New tidal model shows binary tides underestimated by orders of magnitude"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000936,"raw_usage":{"total_tokens":3795,"prompt_tokens":653,"completion_tokens":3142,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":397,"completion_tokens_details":{"reasoning_tokens":3061}},"tokens_in":397,"tokens_out":3142,"duration_ms":22147,"temperature":1.0,"reasoning_tokens":3061,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T01:13:39.837467+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A census of circularization periods: the model predicts solar-type binaries circularize out to about 10-day orbits on the main sequence, while standard models predict almost none; a cluster sample showing circularization only below about 2–3 days would rule out the enhancement. Quantitatively, the model gives a tidal quality factor of about 8.8 × 10^9 for a Sun-like star at a 1-day tidal period, so any empirical measurement from heartbeat-star apsidal motion or hot-Jupiter decay that sits well above this would contradict the prescription's strength—the paper itself notes equilibrium tides can","supporting_citations":[],"review_version":1}