{"id":"02897605-861d-4c7c-80d2-ae0847b5a401","arxiv_id":"2505.15580","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Seven new eclipsing double white dwarf binaries are characterized, and two show the first measured nonzero eccentricity in a compact double white dwarf binary.","lead":"Astronomers found seven pairs of white dwarf stars that eclipse each other, with orbital periods between 45 minutes and 3 hours, and measured their masses, sizes, and temperatures. Two of the longest-period pairs show a tiny, never-before-measured orbital eccentricity, which challenges the usual assumption that these binaries are born circular.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (1) is missing a factor of two; the reported e cos(omega) values are internally inconsistent with the printed formula.","rationale":"The reader's CONDITIONAL verdict is appropriate, but the main load-bearing concern I identify is not the third-body alternative highlighted by the reader. A third body would induce a common-mode delay on both primary and secondary eclipses and cannot easily produce a constant differential offset unless its orbital period is comparable to the inner binary period, which would require it to be dynamically unstable and detectable; the paper's SED and spectral checks provide reasonable exclusion. The reader's point that only e cos(omega) is measured is valid and acknowledged by the authors, but it is a framing issue rather than a technical error. The concrete weakness is the factor-of-two discrepancy between Eq. (1) and the tabulated e cos(omega) values. This is a correctness risk because the central claim is a number, and the paper's equations do not reproduce it. However, internal consistency (the table values matching the standard (2P/pi) formula) suggests the printed equation is a typo rather than a computational error, so the detection itself appears robust. Independently, several aspects support the claim: the offset is detected in multiple nights and filters, the two systems with nonzero offsets have the longest periods and longest GW circularization times, and the other long-period system (J1110+7445) is consistent with zero, arguing against a systematic period-dependent effect. The apsidal precession prediction provides a falsifiable path. I therefore do not change the reader's CONDITIONAL verdict; the correction of Eq. (1) should be required before acceptance, but it does not overturn the qualitative result.","tokens_in":25029,"tokens_out":22837,"duration_ms":189732,"concrete_test":"Independently re-derive the secondary-eclipse timing offset to first order in eccentricity. Using Kepler's equation, M = E - e sin E with E ≈ nu - e sin nu, gives M ≈ nu - 2e sin nu. For conjunction true anomalies nu_1 = pi/2 - omega and nu_2 = 3pi/2 - omega, the time difference is P/2 + (2P/pi) e cos(omega). Confirm this factor against Winn (2010) or Kaplan (2010). Then recompute e cos(omega) for both systems from the reported delta-t_2 and delta-t_LT using both the printed Eq. (1), e cos(omega) = pi(delta-t_2 - delta-t_LT)/P, and the corrected formula, e cos(omega) = pi(delta-t_2 - delta-t_LT)/(2P). If the corrected formula reproduces the Table 3 values (1.75e-3 and 1.55e-3), the results stand and Eq. (1) should be corrected; if the printed formula reproduces them, the eccentricities are actually twice the reported values.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim, e cos(omega) ≈ 2×10^-3 for ZTF J0238+0933 and ZTF J1758+7642, depends on the conversion from the measured secondary-eclipse delay to eccentricity. Section 3.2.1 prints delta-t_e ≈ (P/pi) e cos(omega) as Eq. (1). Substituting the Table 3 values for ZTF J0238+0933 (delta-t_2 = 12.7 s, delta-t_LT = -0.8 s, P = 11822 s) gives e cos(omega) = pi*(13.5 s)/11822 s ≈ 3.6e-3, a factor of two larger than the tabulated 1.75e-3. The same factor-of-two discrepancy appears for ZTF J1758+7642 (tabulated 1.55e-3 vs. pi*11.8/11347 ≈ 3.3e-3). The tabulated values match instead delta-t_e ≈ (2P/pi) e cos(omega), which is the standard first-order result from the equation of center (M ≈ nu - 2e sin nu, yielding an interval between conjunctions of P/2 + (2P/pi) e cos(omega)). Thus the manuscript as written is not reproducible: a reader who applies Eq. (1) to the reported timings obtains eccentricities twice the quoted values. If the fitting code used the printed formula, the reported e cos(omega) values are too small by a factor of two, strengthening the detection but changing the comparison to common-envelope simulations; if the code used the correct formula, Eq. (1) is a typo and the tabulated values are correct. Either way, the central numeric claim needs clarification.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"Van Roestel et al. present the discovery and characterisation of seven eclipsing double white dwarfs found in a systematic ZTF search, with orbital periods from 45 minutes to 3.3 hours. For each system they combine ZTF and CHIMERA light curves, archival multi-wavelength photometry, Gaia parallaxes, and radial velocities from LRIS/ESI in an MCMC joint model. They identify six double helium-core systems and one CO-helium system, derive merger timescales, and report the first measurement of a non-zero eccentricity in a compact double white dwarf binary: two systems (ZTF J0238+0933 and ZTF J1758+7642) show secondary-eclipse delays of about 10–12 s, attributed to e cos(omega) ≈ 1.5–1.8 × 10^-3. The paper also predicts that gravitational-wave orbital decay and relativistic apsidal precession will produce measurable eclipse-timing changes within a decade.","tokens_in":25319,"tokens_out":20464,"duration_ms":154106,"significance":"If the eccentricity detection holds, it is an important constraint on common-envelope evolution, as it suggests that post-common-envelope binaries can retain small eccentricities instead of being fully circularised. The full sample of seven well-characterised eclipsing double white dwarfs is also valuable for population studies and for planning LISA verification binaries. The analysis is thorough: it uses global fits to multiple datasets, MCMC parameter estimation, and explicit comparison with theoretical mass-radius relations and gravitational-wave predictions. However, the manuscript as written contains several internal quantitative inconsistencies in the formulas that support the headline numbers, and these must be resolved before the reported values can be used as published.","major_comments":[{"comment":"The printed conversion between the secondary-eclipse time offset and eccentricity is inconsistent with the tabulated values by a factor of two. For ZTF J0238+0933, using the table values Δt_2 = 12.7 s and Δt_LT = −0.8 s gives Δt_e = 13.5 s, and with P = 11 822 s, Eq. (1) yields e cos(omega) = π Δt_e / P = 3.6 × 10^-3, whereas Table 3 lists 1.75 × 10^-3. The same discrepancy appears for ZTF J1758+7642 (3.3 × 10^-3 from Eq. (1) versus 1.55 × 10^-3 in Table 3). The tabulated values match the standard first-order result Δt_e ≈ (2P/π) e cos(omega) rather than Eq. (1). Please correct Eq. (1) or the table, and state explicitly which formula was used in the light-curve fits, in the abstract, and in the discussion of §5.2.","section":"Section 3.2.1, Eq. (1) and Table 3"},{"comment":"The gravitational-wave period-derivative formula is missing the standard (2π)^{8/3} prefactor. With the chirp mass defined as M_c = (M1 M2)^{3/5} / (M1 + M2)^{1/5}, the correct expression is P_dot = −(96/5) (2π)^{8/3} (G M_c)^{5/3} / (c^5 P^{5/3}); the printed formula has only 2π, so all P_dot values in Table 3 are too small by a factor (2π)^{5/3} ≈ 21. For example, ZTF J0720+6439 should have P_dot ≈ 7 × 10^-13 s/s rather than the tabulated 3.4 × 10^-14 s/s. Consequently the predicted cumulative eclipse-time shift in Section 6 is understated by the same factor (about 13 s, not 0.6 s, in 10 years). Please recompute the affected values and predictions.","section":"Section 3.2.1, Eq. (3) and Table 3"},{"comment":"The stated apsidal precession rate of approximately 0.7 deg/yr does not follow from Eq. (4). Substituting the Table 3 masses and periods for ZTF J0238+0933 (M_tot = 0.65 M_sun, P = 0.13683 d) into Eq. (4) gives 1.54 × 10^-3 deg/cycle, which, at 2670 cycles/yr, corresponds to 4.1 deg/yr; ZTF J1758+7642 gives about 3.7 deg/yr. The factor-of-six discrepancy should be resolved, since Section 6 uses this rate to predict a measurable change in Δt_2 over a decade.","section":"Section 5.2, Eq. (4) and following sentence"}],"minor_comments":[{"comment":"The total number of known eclipsing double white dwarfs is stated as 23 in Section 6 and as 26 in Section 7; please harmonise these numbers.","section":"Section 6 vs Section 7"},{"comment":"The sentence 'we can only determine a lower limit to the eccentricity (Eq. 2)' should reference Eq. (1), which gives the eccentricity-dependent time offset, not Eq. (2), which is the Rømer delay.","section":"Section 5.2"},{"comment":"The phrase 'small eccentricity of ≈ 2 × 10^-3' would be more precise as 'e cos(omega) ≈ 2 × 10^-3, implying a lower limit on e'; as written, it could be read as a direct measurement of e.","section":"Abstract and Section 5.2"},{"comment":"The two-step procedure of first ignoring the Rømer delay in the fit and then correcting for it using the fitted masses is mildly circular, although the small size of the correction (≲1.3 s) relative to the 10–12 s offsets makes a significant bias unlikely; it would be helpful to state the size of the correction explicitly for each system.","section":"Section 3.2.1"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a timely and interesting question, and the data set and modelling effort are substantial. The eccentricity claim is likely to survive once the factor-of-two issue between Eq. (1) and Table 3 is resolved, because the secondary-eclipse delays are measured in multiple nights and filters. However, the errors in Eq. (3) and the apparent inconsistency in Eq. (4) versus the quoted precession rate affect published numbers that the community will want to use. I recommend major revision rather than rejection, provided the authors correct the formulas, recompute the derived quantities, and clarify whether the tabulated eccentricities correspond to the printed equation. The inconsistency between the 23 and 26 totals should also be fixed during revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the headline: this is a genuinely good paper. Seven new eclipsing double white dwarfs, systematically found in ZTF, with joint light-curve/SED/RV modeling, and a plausible first detection of nonzero eccentricity in two compact dWDs. The eccentricity rests on secondary-eclipse delays of ~10-12 seconds, seen in multiple nights and filters, and the paper correctly notes that only e cos(omega) is constrained. The comparison to common-envelope simulations is careful, and the apsidal-precession prediction is testable.\n\nThe main flaw is in print, not in physics. Equation (1) gives delta-t_e ~ (P/pi) e cos(omega). That is missing a factor of two; the standard first-order result is (2P/pi) e cos(omega). The tabulated e cos(omega) values match the correct version, so the fitting code presumably used the right formula and the printed equation is a typo. But as written, the paper is not reproducible: a reader applying Eq. (1) to the measured delays gets eccentricities twice the quoted values. That needs to be fixed.\n\nThe reader's weaker assumptions are real but moderate. Only e cos(omega) is measured, not e, and the paper says so. A third body (brown dwarf or planet) could in principle shift the secondary eclipse; the paper rules it out only qualitatively. A short paragraph with a simple upper limit from the SED and spectra would strengthen the claim. The Romer correction uses fitted masses, but the correction is small and the circularity is mild. No public data or code is a shortcoming, though the methods are described in enough detail to re-implement.\n\nOverall: the central claim likely holds, but the typo must be corrected and the third-body discussion deserves a bit more quantitative attention. This is a strong observational paper that deserves serious refereeing.","headline":"Seven new eclipsing double white dwarfs, a credible first eccentricity detection, and one embarrassing typo in Eq. (1).","tokens_in":26010,"tokens_out":2032,"would_cite":true,"duration_ms":18609,"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":"The paper discovers seven eclipsing double white dwarfs in ZTF data and reports the first measurement of non-zero eccentricity in a compact double white dwarf, based on 10-12 second delays of the secondary eclipse in two ~3-hour binaries.","keywords":["eclipsing double white dwarfs","Zwicky Transient Facility","eclipse timing","orbital eccentricity","common envelope evolution","gravitational wave orbital decay","apsidal precession","helium-core white dwarfs"],"falsifier":"Repeated high-speed eclipse timing of ZTF J0238+0933 and ZTF J1758+7642 over the next decade: if the eccentricity interpretation is right, relativistic apsidal precession should make $\\Delta t_2$ change by roughly $1.4$ seconds in ten years, oscillating on a period of approximately 460 years. A delay that stays constant at the one-second level, or drifts instead of oscillating, would rule out eccentricity as the cause and point to a third body or another offset.","tokens_in":24757,"feed_emoji":"💫","tokens_out":10051,"duration_ms":85350,"temperature":0.7,"pith_summary":"This paper searches Zwicky Transient Facility light curves for eclipsing white dwarfs with periods over an hour and finds seven double-white-dwarf binaries. It characterises them with high-speed photometry, spectra, archival photometry, and joint binary modelling, showing that six are low-mass double helium-core systems and one is a carbon-oxygen/helium-core system. For the two longest-period systems, the secondary eclipse arrives 10-12 seconds late, and the paper attributes this to a small eccentricity with $e\\cos\\omega \\approx 2\\times10^{-3}$, the first non-zero eccentricity measured in a compact double white dwarf. If right, common-envelope ejection can leave binaries slightly eccentric, and relativistic apsidal precession should make the delay change measurably within a decade. The paper also predicts that gravitational-wave orbital decay will produce a measurable eclipse arrival time change for the shortest-period system within a decade.","feed_headline":"First eccentricity measured in compact double white dwarfs","feed_subtitle":"Seven new eclipsing white-dwarf pairs, two with a 10-12 s secondary-eclipse delay set to shift within a decade.","key_machinery":"The load-bearing observable is the secondary-eclipse timing offset $\\Delta t_2$. In a circular orbit the secondary eclipse sits at phase 0.5; a small offset is related to eccentricity by $\\Delta t_e \\approx (P/\\pi)e\\cos\\omega$, and the relativistic light-travel-time (Rømer) delay $\\Delta t_{LT} = (P K_1/\\pi c)(1-q^{-1})$ must be subtracted first. The paper measures $\\Delta t_2$ from high-speed, GPS-timed light curves and models the full binary, including eclipse shapes, archival photometry, Gaia parallax, and radial velocities, with the ellc code. For the two eccentric candidates the measured offsets are $12.7^{+1.2}_{-1.1}$ seconds and $10.5^{+0.8}_{-0.7}$ seconds. The same machinery yields a predicted precession rate from the formula of Gimenez (1985), turning the offset into a testable decade-scale prediction.","core_discovery":"On the paper's own terms, the discovery is that two of the seven newly characterised eclipsing double white dwarfs, ZTF J0238+0933 and ZTF J1758+7642, have secondary eclipses arriving about 10 to 12 seconds later than phase 0.5. After subtracting the computed Rømer delay, which is small because the mass ratios are near unity, the paper interprets the residual offset as a small orbital eccentricity with $e\\cos\\omega \\approx 2\\times10^{-3}$. Because the longitude of periastron $\\omega$ is not constrained, this is a lower limit on $e$ rather than a full eccentricity measurement. The paper states that this is the first time a non-zero eccentricity has been measured in a compact double white dwarf binary, and proposes that these systems emerged from the common envelope with this small eccentricity and have not yet been circularised by gravitational-wave emission. It further predicts that relativistic apsidal precession, at roughly $0.7$ degrees per year, will change the secondary-eclipse delay by about $1.4$ seconds within ten years.","pith_inferences":["Beyond the paper: if small eccentricities routinely survive common-envelope ejection, population-synthesis models of double white dwarfs should include non-zero initial eccentricities, which would affect predicted gravitational-wave phasing and the distribution of merger times.","Beyond the paper: because only $e\\cos\\omega$ is measured and the true eccentricity could be larger if periastron lies near $90$ degrees, a decade of timing will settle this: apsidal precession should make $\\Delta t_2$ oscillate on a period of roughly 460 years.","Beyond the paper: a wide-orbit third body could mimic or modify the timing offset through light-travel-time variations, and although the spectra and spectral energy distribution show no sign of such a companion, continued timing of both primary and secondary eclipses will distinguish a constant precession signal from a periodic outer-orbit wobble.","Beyond the paper: if this detection holds, searches of longer-period eclipsing double white dwarfs should find more eccentric systems, since timing sensitivity to $e\\cos\\omega$ improves as the orbital period grows."],"forward_implications":["The two roughly 3-hour binaries retain $e\\cos\\omega \\approx 2\\times10^{-3}$, while the shorter-period systems are consistent with zero eccentricity, supporting the paper's conclusion that common-envelope ejection can leave a small eccentricity that gravitational waves have not yet erased in longer-period binaries.","Relativistic apsidal precession at about $0.7$ degrees per year should change the secondary-eclipse delay by roughly $1.4$ seconds within ten years, allowing $e$ and $\\omega$ to be disentangled from future timing measurements.","The shortest-period system, ZTF J0720+6439, should show about $0.6$ seconds of eclipse arrival time change in ten years from gravitational-wave angular momentum loss, and ZTF J1110+7445 may be detectable by the space-based gravitational-wave detector LISA with a signal-to-noise ratio near 4 after ten years.","All seven binaries are spiraling inward and are predicted to make contact in 36 million to 1.2 billion years, with the carbon-oxygen/helium-core system expected to merge into an R Coronae Borealis-like star and eventually a carbon-oxygen white dwarf.","The seven new systems more than double the number of known eclipsing double white dwarfs with orbital periods longer than one hour, providing a cleaner sample for population studies of these binaries."],"supporting_citations":[{"why":"Supplies the relation $\\Delta t_e \\approx (P/\\pi)e\\cos\\omega$ linking secondary-eclipse timing offsets to eccentricity, the basis of the detection.","marker":"Winn 2010"},{"why":"Establish the Rømer-delay correction and the eclipse-timing technique for measuring eccentricity in double white dwarfs, and provide the earlier upper limits this detection surpasses.","marker":"Kaplan 2010; Kaplan et al. 2014"},{"why":"Gives the gravitational-wave circularisation rate used to estimate that the two long-period binaries have not yet had time to circularise.","marker":"Peters 1964"},{"why":"Provides the relativistic apsidal precession rate used to predict a measurable change in the secondary-eclipse delay within a decade.","marker":"Gimenez 1985"},{"why":"Supplies the ellc light-curve model used to fit the eclipses and derive the binary parameters.","marker":"Maxted 2016"},{"why":"Hydrodynamic common-envelope simulations that produce post-common-envelope eccentricities of order 0.03, supporting the paper's interpretation of the measured small eccentricity.","marker":"Sand et al. 2020; Glanz & Perets 2021; Bronner et al. 2024"},{"why":"Box Least Squares period-search method used to find the eclipses in the Zwicky Transient Facility light curves.","marker":"Kovács et al. 2002"}],"fun_headline_variants":["First eccentricity measured in compact double white dwarfs","Seven new eclipsing white-dwarf pairs, two with tiny eccentricity","Double white dwarf binaries reveal unexpected eccentricity","Tiny orbital eccentricity found in compact white dwarf binaries","Eclipsing white dwarfs show first eccentricity signal"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the 10 to 12 second late arrival of the secondary eclipse comes entirely from the orbit being slightly non-circular, once the small light-travel-time effect is removed; a third body or another timing shift would break this, and only the combination $e\\cos\\omega$ is actually measured.","fun_headline_variants_meta":{"raw":{"variants":["First eccentricity measured in compact double white dwarfs","Seven new eclipsing white-dwarf pairs, two with tiny eccentricity","Double white dwarf binaries reveal unexpected eccentricity","Tiny orbital eccentricity found in compact white dwarf binaries","Eclipsing white dwarfs show first eccentricity signal"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000688,"raw_usage":{"total_tokens":3164,"prompt_tokens":1040,"completion_tokens":2124,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":656,"completion_tokens_details":{"reasoning_tokens":2042}},"tokens_in":656,"tokens_out":2124,"duration_ms":14190,"temperature":1.0,"reasoning_tokens":2042,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T15:13:21.763597+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeated high-speed eclipse timing of ZTF J0238+0933 and ZTF J1758+7642 over the next decade: if the eccentricity interpretation is right, relativistic apsidal precession should make $\\Delta t_2$ change by roughly $1.4$ seconds in ten years, oscillating on a period of approximately 460 years. A delay that stays constant at the one-second level, or drifts instead of oscillating, would rule out eccentricity as the cause and point to a third body or another offset.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the relation $\\Delta t_e \\approx (P/\\pi)e\\cos\\omega$ linking secondary-eclipse timing offsets to eccentricity, the basis of the detection."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establish the Rømer-delay correction and the eclipse-timing technique for measuring eccentricity in double white dwarfs, and provide the earlier upper limits this detection surpasses."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the gravitational-wave circularisation rate used to estimate that the two long-period binaries have not yet had time to circularise."},{"cited_title":"1985, The Astrophysical Journal, 297, 405","cited_arxiv_id":null,"evidence_quote":"Provides the relativistic apsidal precession rate used to predict a measurable change in the secondary-eclipse delay within a decade."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the ellc light-curve model used to fit the eclipses and derive the binary parameters."},{"cited_title":"T., Schneider, F","cited_arxiv_id":null,"evidence_quote":"Hydrodynamic common-envelope simulations that produce post-common-envelope eccentricities of order 0.03, supporting the paper's interpretation of the measured small eccentricity."}],"review_version":1}