{"id":"b54a084a-b3d1-4d48-ac82-2c6af91905bd","arxiv_id":"2411.11840","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Eccentric mass transfer modeled in MESA shows that many binaries, especially those with black holes, stay eccentric after mass transfer, unlike the standard instant-circularization assumption.","lead":"Astrophysicists implemented a self-consistent model of mass transfer in eccentric binary orbits into the stellar evolution code MESA, replacing the standard shortcut of instantly circularizing such orbits. The models show many binaries remain eccentric after mass transfer, with black-hole binaries especially likely, which changes predictions for X-ray binaries and gravitational wave sources.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Delta-function-periapse MT is the load-bearing assumption: a finite-width MT profile can flip the sign of de/dt near q≃0.76, shifting the q-e bifurcation and the 33% eccentric-post-MT fraction.","rationale":"The reader's weakest_assumption correctly identifies the delta-function-periapse MT approximation as the most load-bearing concern. My independent reading of the paper reaches the same conclusion: the secular equations (3) and (4) are derived under that assumption, and the sign of de/dt near the critical mass ratio directly controls whether a system remains eccentric or circularizes. The paper's Section 5 explicitly concedes the limitation, but the manuscript provides no quantitative test of its impact on the headline 33% fraction or the fitted bifurcation curve (Eq. 10). I considered other possible concerns — e.g., the normalization of the MT rate in the MESA implementation, the choice of the instant-circularization comparison, and the population priors from POSYDON — but none are as directly tied to the central claim or as under-tested. The MT-rate normalization is likely handled correctly because the secular equations are written in terms of an orbit-averaged rate; the instant-circularization comparison convention is a modeling choice that does not affect the existence of eccentric post-MT systems; and the population dependence is acknowledged and discussed. The lack of released code or data is a reproducibility issue, not a correctness flaw. The paper's use of previously peer-reviewed analytic equations and a widely used stellar evolution code provides some independent support, but the specific implementation and the population-level extrapolations are new. Because the delta-function assumption is acknowledged and could shift the quantitative results, a conditional acceptance is appropriate, but the concern does not invalidate the qualitative result that eccentric mass transfer can prevent circularization. No change to the reader's verdict is needed.","tokens_in":95,"tokens_out":9545,"duration_ms":103510,"concrete_test":"Rerun a representative subset of the f-eMT grid, at least the q–e region bracketing Eq. (10), with a finite-width MT profile centered at periapse, for example Ṁ(υ) ∝ exp[-(υ/σ)^2] with σ chosen to match FWHM≈0.12P as in Lajoie & Sills (2011). Replace the delta-function treatment by numerically integrating the Sepinsky et al. (2007b) secular rates over the true anomaly with this phase-dependent MT rate, while keeping all other MESA physics unchanged. Then recompute the post-MT eccentric fraction and re-fit ecrit(q). If the fraction changes by more than about 10 percentage points or the fitted boundary shifts by more than Δe≈0.05, the headline conclusions are not robust to the delta-function assumption.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claims — that about 33% of CO-hosting binaries remain eccentric post-MT and that the q_i–e_i plane cleanly separates outcomes (Eq. 10) — rest on the Sepinsky et al. (2009) secular equations, which assume all MT occurs as a delta function at periapse. In Eq. (4), the bracket changes sign near q≃0.76 (Eq. 9), and whether a system circularizes or pumps eccentricity hinges on this sign. With a delta function, the sign is evaluated only at periapse, where the Roche lobe is smallest and the instantaneous MT rate is formally a spike. Hydrodynamical simulations (Lajoie & Sills 2011) show the MT rate is spread over a finite arc, FWHM≈0.12P, so contributions from other true anomalies are non-negligible. For systems with q near the critical value, the integral over the actual MT rate profile could yield the opposite sign of de/dt, moving the fitted bifurcation (Eq. 10) and changing which systems remain eccentric. The paper explicitly acknowledges this limitation (Sec. 5) but provides no test of its impact on the headline numbers. This is not an internal inconsistency, but it is the assumption on which the population-level conclusions most directly depend.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper implements the Sepinsky et al. (2009) secular eccentric mass-transfer (eMT) equations into the MESA stellar evolution code, coupling the orbital semi-major axis and eccentricity evolution with a self-consistent mass-transfer rate calculation. The study first isolates eMT effects in a simplified grid of 20 Msun donors with 10 Msun black-hole companions, then runs a full-physics grid whose initial conditions are drawn from a POSYDON binary population synthesis model of CO-hosting binaries that initiate Roche-lobe overflow in eccentric orbits. The central results are: (i) a significant fraction of binaries, especially low-mass-ratio BH-hosting systems, remain eccentric after mass transfer; (ii) the final outcome (eccentric vs circularized) is separated by a clean bifurcation in the initial mass-ratio--eccentricity plane, summarized by the polynomial fit in Eq. (10); and (iii) even binaries that naturally circularize end with systematically different donor masses and orbital periods than predicted by the standard instant-circularization treatment, with roughly 5% of systems switching between stable and unstable mass-transfer outcomes.","tokens_in":26442,"tokens_out":7136,"duration_ms":69343,"significance":"If the results hold, this is the first self-consistent implementation of eccentric RLO mass transfer in a detailed stellar evolution code, and it has direct implications for the interpretation of X-ray binaries, wide BH binaries, and gravitational-wave progenitors. The paper provides a concrete falsifiable prediction (the Eq. (10) bifurcation), a systematic comparison against the standard instant-circularization assumption, and a clear statement of the key assumption (delta-function MT at periapse) on which the results rest. The claimed population-level effect, with about one third of CO-hosting binaries remaining eccentric post-MT, is striking and worth pursuing. However, for the reasons detailed in the major comments, the quantitative headline figures are not yet established to the precision claimed, primarily because the load-bearing delta-function approximation is acknowledged but not tested.","major_comments":[{"comment":"The secular orbital evolution equations (3) and (4) assume that all RLO mass transfer occurs as a delta function at periapse. The sign of de/dt in Eq. (4) changes near q ~ 0.76 (Eq. 9), and the balance between circularization and eccentricity pumping for systems near this mass ratio determines the bifurcation boundary in Eq. (10) and hence the 33% eccentric-post-MT fraction. Hydrodynamical simulations (Lajoie & Sills 2011) show that the mass-transfer rate has a finite width (FWHM ~ 0.12 Porb), so contributions from other true anomalies are non-negligible for systems near the transition. The manuscript acknowledges this in Section 5 but does not quantify the impact on the headline numbers. I request a sensitivity test, for example re-evaluating the integrated de/dt for a Gaussian MT window of the hydrodynamically suggested width on a subset of the f-eMT models, to demonstrate that the bifurcation and the one-third fraction are robust.","section":"Section 2.1 and Section 5"},{"comment":"The procedure for determining the periapse mass-transfer rate Mdot0 used in Eqs. (3) and (4) is not specified. The text says that the 'eccentric orbit-averaged MT rate' from the MESA binary module is used, but it does not state how the orbit average is converted into the delta-function amplitude Mdot0. Since Eqs. (3) and (4) are linear in Mdot0, any mismatch between the orbit-averaged rate and the periapse-normalized rate would rescale da/dt and de/dt, shifting both the qcrit values and the fitted bifurcation of Eq. (10). Please give the exact conversion formula and verify that the total mass lost per orbital period in the secular equations matches the mass lost in the MESA MT calculation.","section":"Section 2.3.1"},{"comment":"The polynomial fit for the critical eccentricity is presented as a predictive tool, but it is an unvalidated empirical fit from a single grid: no uncertainties are quoted, no residuals are shown, and the fit is used at the edge of or beyond its stated range (qi in [0, 5.5]). The authors note that high-qi outliers remain eccentric because of winds, indicating that the clean qi-ei separation is not universal. I recommend showing the scatter about the fit with bootstrap or cross-validation uncertainties and stating the expected error when applying Eq. (10) outside the fitted range.","section":"Section 4.2, Eq. (10), and Fig. 6"}],"minor_comments":[{"comment":"The sentence 'It is through this lenses that we interpret...' should read 'It is through this lens that we interpret...' (grammatical error).","section":"Introduction"},{"comment":"The phrase 'following a flat-in-log distribution in the range [1-0.35] days' is garbled; presumably the range should be [1, 10^3.5] days, matching the previous sentence.","section":"Section 2.4"},{"comment":"'which is only analytical treatment currently available' should read 'which is the only analytical treatment currently available'.","section":"Section 2.1"},{"comment":"The row label 'Mixed TFs' should likely read 'Mixed MT'.","section":"Table A1"},{"comment":"The 33% eccentric-post-MT fraction is derived from a single POSYDON BPS realization; the paper notes in Section 5 that the division depends strongly on assumed BPS physics, but a quantitative statement of the expected sensitivity (e.g., to natal kick dispersion or supernova remnant model) would help readers gauge the robustness of the population-level claim.","section":"Section 4.1 and Section 5"}],"recommendation":"major_revision","confidential_remarks":"This is a solid and timely technical contribution, and I would be willing to accept after a revision that addresses the two central technical gaps: a sensitivity test of the delta-function MT assumption and a clear prescription for the Mdot0 normalization. The heavy citation of Sepinsky et al. is appropriate given that the method directly builds on that work; I see no citation or novelty concern. The main risk is that the headline 'one third remain eccentric' number could shift once the finite width of the MT profile is taken into account, so the requested sensitivity analysis is important for the paper's impact."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a genuinely useful paper. It is the first to put the Sepinsky et al. (2009) secular equations into MESA with a self-consistent mass transfer rate, and it makes a clear case that instant circularization is an unjustified shortcut for CO-hosting binaries.\n\nThe cleanest work is the s-eMT parameter study, which isolates the eMT dynamics without tides or winds, and the full f-eMT grid drawn from POSYDON is representative of the underlying population. The comparison against instant circularization is done honestly, and the finding that even naturally circularizing systems end up with different donor masses and orbital periods is significant. The bifurcation in the qi–ei plane is a useful compact result, and the paper explicitly acknowledges the delta-function limitation in Section 5.\n\nThe stress-test concern is legitimate and lands. The secular equations assume all MT occurs at periapse. Since the sign of de/dt flips near q ≈ 0.76, a finite-width MT rate profile (as in Lajoie & Sills) could shift that boundary and therefore the 33% eccentric-post-MT fraction. The paper acknowledges this but does not test its impact on the headline numbers. That is a genuine quantitative gap, not a fatal one: the qualitative findings—binaries can remain eccentric after MT, and eMT changes stable/unstable outcomes—appear in both simplified and full-physics models and are unlikely to vanish with a more realistic MT profile. The polynomial fit in Eq. 10 has no uncertainties, and the population fraction depends on POSYDON assumptions (SN kicks, common envelope), which the authors note. Lack of released inlists/data also makes reproduction harder.\n\nWho gets value from this: the binary population synthesis community, plus observers hunting for eccentric post-MT systems such as the wide Gaia BH binaries. It deserves a serious referee. I would recommend conditional acceptance: ask for a sensitivity test to the delta-function approximation (even a Gaussian profile in true anomaly) and for release of the MESA inlists and the fit data.","headline":"A solid first self-consistent implementation of eccentric mass transfer in MESA; the delta-function-periapse assumption is a real quantitative caveat that shifts the bifurcation boundary but not the qualitative conclusions.","tokens_in":26966,"tokens_out":2317,"would_cite":true,"duration_ms":24029,"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":"Explicitly modeling eccentric Roche-lobe overflow predicts that about one third of star–compact-object binaries remain eccentric after mass transfer, and that even circularized systems differ from instant-circularization predictions.","keywords":["eccentric mass transfer","Roche-lobe overflow","binary stellar evolution","compact objects","orbital eccentricity","population synthesis","X-ray binaries","gravitational-wave sources"],"falsifier":"Simulate a handful of the paper's representative binaries with three-dimensional hydrodynamics that resolve many orbits and measure the orbit-averaged $da/dt$ and $de/dt$; compare with the delta-function predictions, especially for cases near the mass-ratio thresholds $q\\simeq 1$ and $q\\simeq 0.76$. A systematic difference in eccentricity evolution there would falsify the bifurcation boundary. Observationally, a targeted search for post-mass-transfer binaries with stripped helium donors in eccentric ($e > 0.05$) and wide ($P\\sim 10$ to $10^4$ days) orbits could test whether the predicted eccentric population exists.","tokens_in":26002,"feed_emoji":"💫","tokens_out":7763,"duration_ms":69672,"temperature":0.7,"pith_summary":"Binary stars that exchange mass through Roche-lobe overflow are normally assumed to circularize instantly, but this paper asks what happens if the orbit is eccentric and the transfer occurs mainly at closest approach. Using secular equations for eccentric mass transfer embedded in a full stellar evolution code, it finds that a large fraction of binaries can stay eccentric after mass transfer, with about one third of an astrophysically sampled population of star-plus-compact-object binaries remaining eccentric. Even binaries that do circularize through eccentric mass transfer end up with donor masses and orbital periods that differ from instant-circularization predictions, and a small fraction switch between stable and unstable mass transfer. These results matter because eccentric mass transfer would change the predicted populations of X-ray binaries, gravitational-wave sources, and high-energy transients.","feed_headline":"Eccentric mass transfer leaves one-third of binaries eccentric","feed_subtitle":"Self-consistent modeling shows even circularized binaries end up wider than instant-circularization predicts.","key_machinery":"The load-bearing objects are the secular orbital-evolution equations for $da/dt$ and $de/dt$, derived from a model in which all Roche-lobe overflow occurs as a delta function at periapse. These equations are coupled to a detailed stellar evolution code so that the mass-transfer rate, donor response, tides, winds, magnetic braking, and gravitational-wave losses all feed back into the orbit. The sign of each secular rate depends on the mass ratio $q$: when $q$ drops below about 1 the semi-major axis stops shrinking and starts growing, and when $q$ drops below about 0.76 the eccentricity stops decaying and starts being pumped up. The competition between these thresholds and the time spent transferring mass produces the hook-shaped evolution in the period–eccentricity plane and the empirical bifurcation boundary $e_{\\rm crit}(q)$ that separates binaries that circularize from those that remain eccentric.","core_discovery":"The central claim is that Roche-lobe overflow in eccentric orbits, modeled self-consistently with the star and orbit evolving together, does not universally circularize binaries. The paper implements the analytic secular rates for semi-major axis and eccentricity change under a delta-function mass transfer at periapse into a stellar evolution code and applies it both to a simplified grid and to an astrophysical population of stars with compact-object companions. It finds that about 33% of the population remains eccentric ($e > 0.05$) after mass transfer, and that among black-hole-hosting binaries roughly 64% remain eccentric, while neutron-star-hosting binaries mostly circularize. For binaries that do naturally circularize, the eccentric treatment predicts orbital periods roughly 50 to 100% larger and donor mass differences of order 20% compared to the instant-circularization assumption, with about 5% of these systems showing qualitatively different outcomes such as stable versus unstable mass transfer. The paper concludes that the initial mass ratio and eccentricity separate the two outcomes and provides a fitting function for the boundary.","pith_inferences":["If the delta-function assumption survives hydrodynamical checks, population-synthesis predictions for the Galactic X-ray binary population and for compact-object merger rates would need revision, because a sizable fraction of systems would follow a different orbital path than previously assumed.","The predicted population of wide, eccentric post-mass-transfer binaries with stripped helium donors is observationally accessible: Gaia-type astrometry and wide-orbit spectroscopy could test it directly, though selection effects currently make it difficult to detect.","Extending the same treatment to mass transfer between two non-degenerate stars would connect to observed eccentric post-mass-transfer systems such as barium stars and blue stragglers, suggesting that the qualitative remain-eccentric outcome may be common beyond compact-object binaries.","A hydrodynamically calibrated, finite-width mass-transfer model could be checked against the paper's bifurcation boundary; if the boundary shifts, then even the direction of the current bias (more stable, wider orbits) should be re-examined."],"forward_implications":["About one third of compact-object–star binaries in the sampled population remain eccentric after mass transfer, a population that cannot form under the standard instant-circularization assumption.","Black-hole-hosting binaries, with low initial mass ratios, most often remain eccentric and can end up wider and more eccentric than they started; neutron-star-hosting binaries mostly circularize.","Even binaries that circularize naturally through eccentric mass transfer have orbital periods about 50–100% larger and donor masses differing by roughly 20% relative to instant circularization.","A small but non-negligible fraction of binaries switch between stable and unstable mass-transfer outcomes, with the eccentric treatment usually predicting stability when instant circularization predicts instability.","The initial mass ratio and eccentricity of a binary largely determine which outcome occurs, providing a simple fitting function that population synthesis calculations can use to estimate the impact of eccentric mass transfer."],"supporting_citations":[{"why":"Derives the secular equations for $da/dt$ and $de/dt$ for non-conservative eccentric mass transfer that are the core of the implementation.","marker":"Sepinsky et al. 2009"},{"why":"Provides the delta-function formalism, the critical mass-ratio approximations for the sign of $da/dt$ and $de/dt$, and the geometric definitions used in the equations.","marker":"Sepinsky et al. 2007b"},{"why":"Extends beyond the delta-function approximation but only for conservative transfer, used to justify keeping the non-conservative formalism.","marker":"Hamers & Dosopoulou 2019"},{"why":"Supplies the Roche-lobe overflow mass-transfer rate formula adapted to the eccentric orbit-averaged scheme.","marker":"Ritter 1988"},{"why":"Together with Ritter, gives the mass-transfer rate prescription used to compute the transfer rate in the stellar evolution code.","marker":"Kolb & Ritter 1990"},{"why":"Defines the Roche-lobe radius formula used to determine when overflow begins in eccentric orbits.","marker":"Eggleton 1983"},{"why":"Sets the binary stellar evolution and population-synthesis physics, including tides, winds, magnetic braking, and Eddington-limited accretion, used in the full models.","marker":"Fragos et al. 2023"},{"why":"Provides the population-synthesis infrastructure, initial-condition procedures, and the comparison framework for binary evolution outcomes.","marker":"Andrews et al. 2024"},{"why":"Defines the binary population synthesis setup used to obtain the astrophysical initial conditions and the comparison to observed X-ray binaries.","marker":"Rocha et al. 2024"},{"why":"Supplies the hydrodynamical evidence that the eccentric mass-transfer rate has finite width around periapse, the main caveat against the delta-function assumption.","marker":"Lajoie & Sills 2011"}],"fun_headline_variants":["Self-consistent eMT leaves 1/3 of binaries eccentric","Eccentric mass transfer defies circularization in self-consistent model","Simulating eMT: 33% stay eccentric, rest differ from instant-circularization","Eccentric mass transfer: circularization not inevitable in self-consistent runs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The weakest link is the assumption that all mass transfer happens in an instant at closest approach, while real hydrodynamical flows spread mass transfer over a finite fraction of the orbit; if that spread changes the secular rates, the predicted eccentricity evolution and the boundary between circularizing and eccentric outcomes could shift.","fun_headline_variants_meta":{"raw":{"variants":["Self-consistent eMT leaves 1/3 of binaries eccentric","Eccentric mass transfer defies circularization in self-consistent model","Simulating eMT: 33% stay eccentric, rest differ from instant-circularization","Eccentric mass transfer: circularization not inevitable in self-consistent runs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000201,"raw_usage":{"total_tokens":1409,"prompt_tokens":1006,"completion_tokens":403,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":622,"completion_tokens_details":{"reasoning_tokens":320}},"tokens_in":622,"tokens_out":403,"duration_ms":5220,"temperature":1.0,"reasoning_tokens":320,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T18:04:16.094371+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate a handful of the paper's representative binaries with three-dimensional hydrodynamics that resolve many orbits and measure the orbit-averaged $da/dt$ and $de/dt$; compare with the delta-function predictions, especially for cases near the mass-ratio thresholds $q\\simeq 1$ and $q\\simeq 0.76$. A systematic difference in eccentricity evolution there would falsify the bifurcation boundary. Observationally, a targeted search for post-mass-transfer binaries with stripped helium donors in eccentric ($e > 0.05$) and wide ($P\\sim 10$ to $10^4$ days) orbits could test whether the predicted eccentric population exists.","supporting_citations":[],"review_version":1}