{"id":"559f50c7-fe02-41a0-933c-c7f11e5d8437","arxiv_id":"2501.13929","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"F-mode tidal models reproduce late red-giant circularization but under-predict early red-giant circularization by about a factor of two, extending the main-sequence tidal discrepancy to the early red giant branch.","lead":"Using about 30,000 red-giant binaries from Gaia, this paper tests whether tidal theory explains when binary orbits become circular. It finds that standard tidal models work for late red giants but predict too little circularization for younger giants, pointing to missing physics.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The predicted a/R≈4 circularization boundary is set by the frequency-dependent eddy viscosity of Eq. (7); if its fast-tide suppression is too strong in early-RGB envelopes, the claimed factor-of-two discrepancy may vanish.","rationale":"The reader's conditional verdict already identifies the frequency-dependent eddy viscosity, Eq. (7), as the weakest assumption. My stress-test concurs: the central claim is that f-mode tides with this viscosity predict a fixed β≈4 across the RGB, whereas early-RGB Gaia binaries appear circularized out to β≈6 with cool islands to a/R≈10–15. The theoretical boundary derives from a turbulence model whose normalization and high-frequency suppression are explicitly uncertain and contested in the literature. The paper itself flags the possibility in Section 5, so this is a documented limitation rather than an invented one. A direct recomputation with an alternative or unsuppressed viscosity law would settle whether the discrepancy survives. If it does not survive, the paper's conclusion about additional early-RGB circularization loses its force; if it does survive, the conclusion is strengthened. I do not see a reason to move the verdict away from the reader's CONDITIONAL assessment, since the paper is otherwise clearly argued and the needed sensitivity analysis is well defined. The selection-function issue for the Gaia cool island is also worth attention, but it is secondary because the upper-envelope discrepancy at β≈6 does not depend on the cool-island definition alone. For these reasons the appropriate action is to keep the verdict unchanged pending the proposed viscosity test.","tokens_in":20063,"tokens_out":7994,"duration_ms":81424,"concrete_test":"Recompute the population synthesis of Figure 10 with the same MESA/GYRE f-mode machinery but replace Eq. (7) with a frequency-independent ν_c = 5 ℓ_c v_c, and separately with an alternative fast-tide prescription such as that of Terquem (2023). Measure the theoretical β envelope and cool-island edge in the early-RGB radial bins (0.7 < log R/R⊙ < 1.3). If either moves from β≈4 and a/R≈6 toward β≈6 and a/R≈10–15, the claimed factor-of-two discrepancy is an artifact of the chosen viscosity law and the central conclusion is unsupported; if β remains near 4, the viscosity assumption is not the limiting factor.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing assumption is the turbulent viscosity prescription in Eq. (7), specifically the fast-tide suppression branch ν_c ∝ (ω_c/|ω|)^2 adopted from Duguid et al. (2020). The theoretical circularization boundary β ≡ a(1−e)/R ≈ 4 is not a pure prediction of linear f-mode dynamics; it is produced by the dissipation integrals in Eq. (B7), which weight the spatially varying ν_c profile inside the convective envelope. In early-RGB binaries the relevant tidal periods are shorter relative to the local convective turnover frequency, so deep-envelope zones can fall into the suppressed branch, and the authors themselves state in Section 5 that their viscosity prescription may be 'too conservative'. If the true turbulent dissipation is less suppressed, or the eddy-viscosity prefactor is larger, the theoretical boundary should move outward, plausibly reaching the observed early-RGB values of β ≈ 6 and cool islands at a/R ≈ 10–15. That would remove the stated need for an 'additional circularization process' without any new physics. The omission of quantitative uncertainty on the adopted ν_c, combined with the unresolved literature disagreement (Goodman & Oh 1997; Terquem 2021, 2023; Barker & Astoul 2021), makes the size of the claimed discrepancy currently unsecured. A secondary, smaller issue is the period-to-a/R conversion: Section 2.2 assumes total mass 1.5× the primary, while the theoretical model has M + M′ = 2.5 M⊙ for a 1.5 M⊙ primary (1.67×), a ~4% effect on a/R that does not by itself explain the discrepancy.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper confronts state-of-the-art tidal theory with the Gaia DR3 sample of ~31,500 red-giant binaries. The authors build a forward model: a 1.5 M⊙ primary (with a 1.0 M⊙ companion) is evolved through the RGB and AGB in MESA; linear tidal responses are computed in an f-mode/p-mode eigenmode expansion (and in a zero-frequency 'equilibrium' variant) using GYRE; dissipation is supplied by a frequency-dependent eddy-viscosity model (Eq. 7) calibrated to the local simulations of Duguid et al. (2020); and a synthetic population is generated from observationally motivated initial period and eccentricity distributions (Section 4.3). The model predicts that the eccentricity upper envelope in a/R space is pinned at β ≡ a(1−e)/R ≈ 4 at every RGB stage. The Gaia data instead show an envelope that starts at β ≈ 6 in the early RGB and declines to β ≈ 3 near the tip, plus 'cool islands' of circular orbits extending to a/R ≈ 10–15 in early stages, about twice the predicted reach. The paper concludes that f-mode tides with the adopted viscous prescription describe the late RGB well but that 'an additional circularization process' is needed during the early RGB. It also reports theory-driven predictions of strong tidal spin-up of giant primaries (with consequences for mass loss) and of binaries beginning Roche-lobe overflow while still significantly eccentric.","tokens_in":20390,"tokens_out":17501,"duration_ms":146433,"significance":"If substantiated, the early-RGB discrepancy is a significant result: it connects the long-standing main-sequence tidal-circularization problem to evolved stars and points to missing dissipation physics (internal gravity waves, inertial waves, mode locking, or magnetic effects) in early-RGB convective envelopes. The paper's method is a genuine strength: the predicted a/R ≈ 4 boundary is produced from first principles (a linear f-mode response plus an eddy viscosity from independent local simulations), with no parameter fitted to the Gaia circularization data, and the a/R rescaling makes the comparison clean and physically motivated. The paper also makes concrete, falsifiable predictions: spin-up to Ω/ω_dyn ≈ 0.3 at RLOF, eccentric RLOF onset, and synchronization out to ~5000-day periods on the RGB. The main weakness is identified by the authors themselves: Section 5 notes that the adopted viscosity prescription 'may be too conservative,' and the early-RGB discrepancy is concentrated precisely in the regime where the uncertain fast-tide suppression branch of Eq. (7) acts (short-period orbits, with ω/ω_c ≳ 1–5 in parts of the envelope).","major_comments":[{"comment":"The central quantitative claim — that f-mode tides plus the adopted eddy viscosity cap tidal circularization at a/R ≈ 4 on the entire RGB, and hence that 'there needs to be an additional circularization process during the early RGB' (Section 6) — depends on the high-frequency branch of the viscosity prescription in Eq. (7). The early-RGB binaries at issue have the shortest orbital periods (Fig. 4), so their tidal frequencies place a significant fraction of the convective envelope in the frequency-dependent suppressed regimes (1 ≲ ω/ω_c ≲ 5 in the middle branch and ≳ 5 in the fast branch), including deeper zones at the largest ratios. Section 5 concedes both that the Duguid et al. (2020) scaling suppresses dissipation 'stronger in the early RGB' and that the adopted prescription 'may be too conservative.' Since the comparison in Fig. 11 is a factor-of-two mismatch in exactly this regime, the manuscript should include a sensitivity study in which the predicted β and cool-island edges are recomputed for the plausible range of viscosity scalings (e.g., Zahn (1977) without high-frequency suppression, the Goldreich & Nicholson (1977) form with varied prefactor, and the alternatives of Goodman & Oh (1997) and Terquem (2021, 2023)), and the resulting change in the predicted boundary should be reported. Absent this test, the discrepancy that motivates the paper's main conclusion is not established.","section":"§3.3, Eq. (7); §5; §6"},{"comment":"The measured quantities used for the discrepancy claim — the upper-envelope values β ≈ 6.5 → 3.0 and the cool-island edges at a/R ≈ 10–15 (Fig. 4, right panels) — are reported without uncertainties. Section 2.2 states that these uncertainties 'are hard to quantify at this moment,' but the central result is precisely a comparison of these numbers with the theoretical β ≈ 4 boundary (Fig. 11). To support a factor-of-two claim, the paper should provide at least a bootstrap estimate of the sampling uncertainty in β and the island edges (e.g., resampling the binary catalog), propagate the quoted 2–4% eccentricity errors, and test the sensitivity of the derived β values to the adopted Gaia masses in the period-to-a/R conversion. Without error bars on the observational limits, the statistical significance of the early-RGB discrepancy cannot be assessed.","section":"§2.2, Figs. 4 and 11"},{"comment":"The two sides of the comparison use inconsistent mass assumptions. Observed a/R values are computed with a total mass of 1.5× the Gaia primary mass (Section 2.2), while the theoretical population is evolved with M + M′ = 2.5 M⊙ for a 1.5 M⊙ primary, i.e., 1.67× (Section 3.1); this produces a ~4% systematic offset in a/R between data and model. In addition, the model uses a single 1.5 M⊙ track while the observed sample spans roughly 1–10 M⊙ (Fig. 3); the claim in Section 3.1 that radius rather than mass controls the tidal evolution is plausible but untested. Computing even one additional track (e.g., 1.0 M⊙ and 2.5 M⊙ primaries) would show whether the predicted a/R ≈ 4 boundary is robust to the primary mass. These checks are inexpensive and directly bear on the quality of the comparison in Fig. 10.","section":"§2.2 vs. §3.1; Fig. 10"}],"minor_comments":[{"comment":"The boundary-conditions sentence contains a duplicated phrase: 'that the normal and that the normal and tangential stresses vanish'; it should read 'that the normal and tangential stresses vanish.'","section":"Appendix B, text after Eq. (B5)"},{"comment":"The label '1.5M' on the evolutionary track should read '1.5 M⊙' in both the caption and the figure legend.","section":"Fig. 1 caption"},{"comment":"The red-clump excision combines a density-peak criterion with a P < 300 day cut; a brief statement of how the derived β values and cool-island edges change when this cut is varied (e.g., 200 or 400 days) would strengthen confidence in the sample boundary.","section":"§2.2, red-clump excision"},{"comment":"The text lists the software packages used (MESA, GYRE, scipy, numpy) but does not state where the MESA tracks, pre-computed mode properties, or the orbital-integration code can be obtained; making these available would substantially aid reproducibility of the population-level predictions in Figs. 10 and 12.","section":"Software and data availability"}],"recommendation":"major_revision","confidential_remarks":"This is a well-executed and potentially high-impact paper, but the headline conclusion ('an additional circularization process during the early RGB') rests on the a/R ≈ 4 boundary produced by one adopted viscosity model whose uncertainty the authors themselves flag. The authors already hint at the needed test in Section 5 ('It is unclear if a different viscous prescription would help reconcile the early RGB discrepancy'); a quantitative sensitivity study should be a condition of acceptance. Note also the mild tension between Section 3.3, which downplays the fast-tide suppression controversy as 'relatively unimportant,' and Section 5, which invokes it as a possible resolution of the discrepancy; the revision should address this explicitly. If the sensitivity study moves the predicted boundary to β ≈ 6, the paper's contribution shifts from 'evidence for missing physics' to 'a constraint on eddy-viscosity scaling' — still valuable and publishable, but the framing and abstract would then need adjustment."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid, genuinely new forward-model comparison, and the early-RGB discrepancy is worth taking seriously—but its size is not yet secure because the predicted boundary is set by the adopted eddy-viscosity prescription. The paper deserves a serious referee.\n\nWhat is new: rather than fitting a circularization boundary, Dewberry & Wu evolve a 1.5+1.0 Msun binary through the RGB using MESA, computing f-mode tidal dissipation from the linear stellar response with a frequency-dependent eddy viscosity, and compare against ~30,000 Gaia DR3 giants binned by primary radius. That is a step up from the 28 binaries of Verbunt & Phinney and the 234 of Price-Whelan & Goodman. The qualitative successes are real: the model produces roughly the observed late-RGB circularization limit (β ≈ 3–4), the cool island at a/R ≈ 6, and the spin-up near RLOF. The failure at early RGB (observed β ≈ 6 and cool islands out to a/R ~ 10–15) is a localized, testable claim.\n\nThe soft spots are in proportion. The largest is the viscosity. The predicted β ≈ 4 boundary comes from weighting the dissipation integral with the fast-tide suppression in Eq. 7, not from f-modes alone. The authors admit the prescription may be 'too conservative,' and the literature on this suppression is genuinely unsettled. A less suppressed viscosity in early-RGB convective envelopes moves the boundary outward and could erase much of the claimed factor-of-two discrepancy. So the quantitative tension is not yet secured. Second, the observed β and cool-island edges are quoted without uncertainties; the paper says they are hard to quantify, but the central comparison needs those error bars to be definitive. The period-to-a/R conversion assumes Gaia masses and a fixed total-mass multiplier; the 4% inconsistency with the model's 1.67× is minor and not the issue. No code or synthetic catalogs are released, which makes independent verification slower, though the method is described in enough detail to reproduce.\n\nThe paper is honest about its limitations and lays out alternatives (internal gravity waves, mode locking, inertial waves, magnetic fields). It is not overfitting: no parameter is tuned to the giant binaries. This is exactly the kind of paper that should go through peer review, with the referee pushing for quantified uncertainties and a viscosity-sensitivity study. I would bring it to reading group and would cite it.","headline":"A genuinely new forward-model test of tidal circularization against 30,000 Gaia giants; the early-RGB discrepancy is real but its size is not yet secure because it sits on the adopted eddy-viscosity prescription.","tokens_in":20943,"tokens_out":3150,"would_cite":true,"duration_ms":27926,"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":"Gaia red-giant binaries stay circular twice as far out as f-mode tidal theory allows, pointing to an extra circularization process in the early red giant branch.","keywords":["tidal circularization","red giant branch","Gaia DR3","binary stars","eddy viscosity","dynamical tides","f-mode tides","eccentricity distribution"],"falsifier":"Compute the early-RGB population with an eddy-viscosity prescription that omits or weakens the fast-tide suppression in Equation (7); if the predicted circularization boundary shifts from $a/R \\approx 4$ to $\\beta \\approx 6$, the claimed need for an additional circularization process would not be supported.","tokens_in":19856,"feed_emoji":"🔭","tokens_out":11304,"duration_ms":92093,"temperature":0.7,"pith_summary":"This paper uses about 30,000 red-giant binaries from Gaia DR3 to test whether tidal dissipation theory explains where binary orbits become circular as the primary star swells along the red giant branch. The authors compute the tidal evolution of a 1.5-solar-mass primary with f-mode dynamical tides damped by an eddy viscosity in the convective envelope, then compare predicted eccentricity distributions with the observed ones in the dimensionless plane $a/R$ (semi-major axis over stellar radius). The model produces a fixed circularization limit at $a/R \\approx 4$ at every red-giant stage, whereas the data show early red giants circularized out to $\\beta \\equiv a(1-e)/R \\approx 6$, with extended cool islands of circular orbits reaching $a/R \\sim 10$--$15$: roughly twice the predicted reach. The paper concludes that f-mode tides are adequate for the late red giant branch but that an additional circularization process, still unidentified, must operate during the early red giant branch. If correct, the result sharpens both where tidal theory succeeds and where binary evolution models are missing physics.","feed_headline":"Red-giant binaries circularize twice as far as tidal theory allows","feed_subtitle":"A 30,000-binary Gaia study finds early red giants circular out to a/R≈6 where f-mode models stop at ≈4.","key_machinery":"The central machinery is the 'f-mode' tidal model: the tidal response of the evolving 1.5-solar-mass primary is expanded in the star's fundamental and lowest-order pressure modes, and each mode's dissipation is evaluated with a frequency-dependent eddy viscosity from convective turbulence. The resulting dissipative coefficients $\\kappa_{mk}$ enter secular equations for orbital decay, eccentricity damping, and spin evolution that are integrated along the stellar evolutionary track. The organizing diagnostic is the dimensionless ratio $a/R$, with $\\beta \\equiv a(1-e)/R$ marking the eccentricity upper envelope; plotting eccentricity against $a/R$ makes both theory and data nearly stationary across evolutionary stages, so the gap between the predicted boundary at $a/R \\approx 4$ and the observed early-RGB reach at $\\beta \\approx 6$ becomes the measurable quantity that drives the conclusion.","core_discovery":"On its own terms, this paper claims that f-mode tidal theory with a frequency-dependent eddy viscosity predicts a universal circularization boundary of $a/R \\approx 4$ for red giants, independent of the star's position on the red giant branch. Against that prediction, Gaia binaries in the early red giant phase are observed circularized out to $\\beta \\approx 6$ and show an extended cool island of circular orbits out to $a/R \\sim 10$--$15$, about twice the theoretical reach; in the late red giant branch the observed reach shrinks to $\\beta \\approx 3$ and agrees with the model. From this the authors conclude that an additional circularization process must act in the early red giant branch, and they note the same missing physics may underlie the even larger discrepancy long seen for main-sequence binaries. The calculations also produce two by-products: tides can spin giant primaries up to rotation rates that should change their mass loss, and many binaries may enter Roche-lobe overflow while still significantly eccentric.","pith_inferences":["A direct test would be to rerun the population calculation with a weaker fast-tide suppression in the eddy-viscosity law; if the predicted boundary moves from $a/R \\approx 4$ to $\\approx 6$ in the early RGB, the claimed missing process would dissolve.","Because the model predicts a universal $a/R$ limit, observed circularization periods should scale linearly with primary radius; future samples with asteroseismic radii can test this scaling and separate radius errors from genuine tidal effects.","The early-RGB excess could come from internal gravity waves breaking in the still-substantial radiative core, so comparing same-radius giants of different masses (hence different core sizes) could discriminate among the proposed mechanisms.","Tidal spin-up out to periods of roughly 5000 days on the RGB and 10,000 days on the AGB should be measurable in asteroseismic rotation rates of red giants, providing an independent check of the dissipation model."],"forward_implications":["If the f-mode limit is fixed at $a/R \\approx 4$, every red-giant binary with dimensionless pericenter $\\beta \\lesssim 4$ should be effectively circular; the early-RGB binaries with $\\beta \\approx 6$ mark where that rule fails.","Tidal spin-up can bring close giant primaries near orbital synchronization, roughly $\\Omega/\\omega_{\\rm dyn} \\approx 0.3$ at Roche-lobe overflow for a 2/3 mass ratio, so mass-loss prescriptions based on slowly rotating single stars need revision for these binaries.","Binaries with initial periods near $10^3$ to $3\\times 10^3$ days can begin Roche-lobe overflow while still substantially eccentric, contradicting older equilibrium-tide expectations of full circularization beforehand.","If the early-RGB excess is real, it may share a mechanism with the main-sequence cool island, turning two separate tidal-theory failures into one missing-process problem."],"supporting_citations":[{"why":"Supplies the DR3 non-single-star catalogue from which the ~30,000 red-giant binary sample and its mass and radius estimates are drawn.","marker":"Gaia Collaboration et al. 2023"},{"why":"Provides the frequency-dependent eddy-viscosity scaling (Equation 7) whose fast-tide suppression sets the predicted a/R about 4 limit.","marker":"Duguid et al. (2020)"},{"why":"Establishes the equilibrium-tide and eddy-viscosity framework for convective envelopes that the f-mode model generalizes.","marker":"Zahn (1977)"},{"why":"Pioneers the red-giant tidal circularization comparison that this paper extends to the Gaia sample.","marker":"Verbunt & Phinney (1995)"},{"why":"Extends red-giant circularization tests to APOGEE binaries and connects the circularization period to evolutionary stage.","marker":"Price-Whelan & Goodman (2018)"},{"why":"Gives the Rayleigh eccentricity distribution used to generate the initial binary orbits for the population integrations.","marker":"Wu et al. 2024"},{"why":"Supplies the f-mode expansion method described in Appendix B for computing tidal responses at finite frequency.","marker":"Vick & Lai (2020)"},{"why":"Motivates the fast-tide suppression of eddy viscosity incorporated into the adopted prescription.","marker":"Goldreich & Nicholson (1977)"}],"fun_headline_variants":["Red giants circularize twice as far as theory allows","Gaia red-giant binaries beat tidal circularization limit","Early red giants show circular orbits beyond tidal reach","Tidal theory fails to explain red-giant binary orbits","Red-giant binaries expose gap in tidal evolution model"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The predicted circularization limit rests on the assumed eddy-viscosity prescription for convective envelopes, and if the true viscosity is stronger than the fast-tide suppression allows, the theory's boundary would move outward and the early-red-giant discrepancy would shrink or vanish.","fun_headline_variants_meta":{"raw":{"variants":["Red giants circularize twice as far as theory allows","Gaia red-giant binaries beat tidal circularization limit","Early red giants show circular orbits beyond tidal reach","Tidal theory fails to explain red-giant binary orbits","Red-giant binaries expose gap in tidal evolution model"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000185,"raw_usage":{"total_tokens":1323,"prompt_tokens":947,"completion_tokens":376,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":563,"completion_tokens_details":{"reasoning_tokens":299}},"tokens_in":563,"tokens_out":376,"duration_ms":3707,"temperature":1.0,"reasoning_tokens":299,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T15:28:28.494669+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the early-RGB population with an eddy-viscosity prescription that omits or weakens the fast-tide suppression in Equation (7); if the predicted circularization boundary shifts from $a/R \\approx 4$ to $\\beta \\approx 6$, the claimed need for an additional circularization process would not be supported.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Pioneers the red-giant tidal circularization comparison that this paper extends to the Gaia sample."}],"review_version":1}