{"id":"2235c742-e775-44e1-a0c8-b8150cae16f9","arxiv_id":"2602.20015","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Warm Jupiters are frequently eccentric (~59% in a high-eccentricity mode), while warm sub-Neptunes and sub-Saturns are mostly low-eccentricity, with the switch occurring at a break radius of about 9.8 Earth radii.","lead":"Using TESS satellite data, the authors measured how stretched the orbits are for 347 warm, single-transiting planets and found that larger planets are far more likely to have eccentric orbits, with a sharp transition near 9.8 Earth radii. The result separates calm, disk-formed planets from dynamically excited giants and gives follow-up targets for testing formation models.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Continuous model's fixed low-e component mean conflicts with the discrete Jovian fit and may inflate the reported 59% high-e fraction.","rationale":"The reader's flagged grazing-transit cut is a legitimate, untested selection effect and should be checked. However, the paper's own internal results already contain a direct model-robustness test — the discrete size-bin fits — and they disagree with the continuous headline number: the Jovian high-e fraction is 0.37±0.10 in Table 2 versus 0.59±0.13 in Table 3. The mechanism is identifiable: the continuous model fixes µlow across all radii, while the discrete fits show µlow rising from ≈0.03 for sub-Neptunes to ≈0.23 for Jovians. This makes the central quantitative claim vulnerable to an equally plausible model choice, which is more fundamental than the b-threshold. The two conflicting abstracts in the manuscript (219 planets, 8–50 d, Rbr=9.2, 65%, 2.7σ versus 347 planets, 8–200 d, Rbr=9.8, 59%, >4σ) and the internal sample-size mismatch (374 vs 347) are editorial but reinforce that the quantitative summary is not yet stable. The inference framework itself is standard HBM, and the direction of the radius trend is supported by the discrete fits, so this is addressable with a targeted model comparison rather than a rejection.","tokens_in":24457,"tokens_out":12967,"duration_ms":122529,"concrete_test":"Refit the radius-continuous model allowing the low-e component mean to vary with radius, e.g. µlow(log Rp) = µlow,s + (µlow,l − µlow,s) σ(−λ′(log Rp − log Rbr′)) with priors matched to Table 3, or a linear-in-log-Rp form; keep all other components and selection terms unchanged. Compare the posterior Jovian high-e fraction, Rbr, and WAIC/LOO against the fixed-µlow model. If the Jovian high-e fraction moves substantially toward the Table 2 value (≈0.37) or Rbr shifts by more than the reported credible interval, the abstract's quantitative claims must be revised and the varying-µlow model should become the primary result.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing weakness is the radius-continuous sigmoid model's assumption of radius-independent component means. Table 2's Jovian-only two-component Beta fit gives a high-e weight w2 = 0.37+0.10-0.11 with low-e component mean µ1 = 0.225+0.050-0.050. Table 3's continuous model instead forces a single global µlow = 0.070+0.026-0.068 and then reports a Jovian high-e fraction of 0.59+0.13-0.13. Because the continuous low-e component is concentrated at e≈0.07, moderately eccentric Jovians (e≈0.2–0.3) that the discrete fit places in the low-e component are pushed into the high-e component. The headline '59% high-e' is thus not a direct empirical fraction but a consequence of the fixed-µlow parameterization. The same assumption drives the reported break radius Rbr ≈ 9.8 R⊕: the model can only express the rise in mean eccentricity with radius through πlow(Rp), so any radial drift in the low-e component's location is absorbed as an apparent increase in the high-e fraction. The paper does not test a model with radius-dependent µlow, and it explicitly declines to compute model evidence (§5.1). The central quantitative claim is therefore not yet established as robust to an equally plausible parameterization.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper infers the eccentricity distribution of warm, single transiting planets in TESS using the photoeccentric effect within a hierarchical Bayesian framework. It first analyzes discrete radius bins (sub-Neptunes, sub-Saturns, Jovians) and then introduces a radius-continuous three-stage mixture model in which the low-eccentricity fraction varies with planet radius through a logistic sigmoid. The central claims are that the warm-single population is best described by two eccentricity components, that the high-eccentricity fraction increases strongly with radius with a break radius near 9.8 R⊕, and that warm Jovians are frequently eccentric (59% high-e fraction). A non-negligible tail of high-eccentricity sub-Neptunes is also claimed.","tokens_in":24886,"tokens_out":5843,"duration_ms":56019,"significance":"If the results hold, this would be an important homogeneous, all-sky characterization of the radius–eccentricity relation for warm single planets, unifying previously segmented TESS and Kepler studies. The analysis pipeline is state-of-the-art: it uses photoeccentric likelihoods with hierarchical Bayesian modeling, validates a Gaussian-approximation shortcut against full posterior chains, performs a 30% random-dropping robustness test, and checks individual photoeccentric constraints against published radial-velocity eccentricities. These are real strengths. However, the central quantitative claims are conditional on a parameterization that is not formally model-selected and is not tested against equally plausible alternatives; the current manuscript overstates the strength of evidence for bimodality.","major_comments":[{"comment":"The radius-continuous model fixes the component means µ_low and µ_high as global constants. Table 2's independent discrete Beta-mixture fits show that the low-e component mean rises from 0.028 for sub-Neptunes to 0.111 for sub-Saturns to 0.225 for Jovians. Under the continuous model with µ_low = 0.070, moderately eccentric Jovians (e ~ 0.2–0.3) that the discrete fit places in the low-e component are reassigned to the high-e component. This explains the difference between the Jovian high-e fraction of 0.59 in Table 3 and w2 = 0.37 in Table 2. Because the model can express radial change only through π_low(Rp), a radial drift in the location of the low-e component is absorbed as an apparent increase in high-e membership. The paper does not test a model with radius-dependent µ_low. This is load-bearing for the headline 59% high-e fraction and for Rbr. Please rerun with a flexible µ_low(Rp) (","section":"§4.2, Eq. (5), Table 3"},{"comment":"The abstract states that the population is 'best described by two components' and that bimodality is detected at '>4σ', but §5.1 explicitly says 'we do not compute Bayesian evidences for the hierarchical models considered here.' No model comparison between the single-Beta and Beta-mixture models is performed. A credible interval excluding w2 = 0 at 4σ is evidence for a non-zero second component, not for bimodality or for the mixture being the best description. The language in the abstract and conclusions should be softened unless formal model selection (e.g., PSIS-LOO, WAIC, or cross-validated predictive comparison) is added. As it stands, the central 'two-component' claim is not formally supported.","section":"§5.1 and Abstract"},{"comment":"The manuscript contains inconsistent sample definitions. §3.2 says '374 planets are therefore used in our final eccentricity distribution,' while §4.1 and the full-text abstract use N = 347. The initial abstract block reports N = 219, P = 8–50 days, Rbr = 9.2 R⊕, and a 65% Jovian high-e fraction, whereas the full-text abstract and body report N = 347, P = 8–200 days, Rbr = 9.8 R⊕, and a 59% high-e fraction. These are incompatible sets of numbers. The reader cannot tell which sample and which results are the definitive ones. This must be reconciled in any revised version.","section":"§3.2, §4.1, and both abstracts"},{"comment":"The removal of 47 targets with >50% of posterior samples at b ≥ 0.9 is not tested for sensitivity. Grazing/high-impact-parameter transits show shortened durations that can mimic the photoeccentric signature of high eccentricity; if the removed systems are preferentially eccentric, the inferred high-e fraction and break radius would be biased. Unlike the 13 unconstrained-eccentricity systems, which were explicitly re-tested, this larger and arguably more dangerous cut has no robustness test. Please report the sensitivity of the §4.2 results to this cut, or provide a quantitative argument that the bias direction and magnitude are negligible.","section":"§3.2"}],"minor_comments":[{"comment":"The title in the LaTeX source has spacing errors ('W arm', 'inTESS'), and the first abstract block differs from the full-text abstract. This appears to be a version-control issue, but as submitted it is confusing.","section":"Title/Abstracts"},{"comment":"There is a duplicated word: 'find our results are consistent consistent to well within 1σ.'","section":"§4.1"},{"comment":"The phrase 'a 4σ credible interval of w2 > 0' is unclear; credible intervals are two-sided, and a one-sided exclusion of zero is not reported in standard form. Please state the actual posterior probability or interval.","section":"§5.1"},{"comment":"The table caption uses 'Rp < 6 R⊕' while the text mentions 'Rp < 6 R⊕' and 'Rp < 4 R⊕' in different places; please make the selection criterion consistent and explicit.","section":"Table 4"}],"recommendation":"major_revision","confidential_remarks":"The manuscript as submitted appears to mix two versions of the abstract and sample definition (219 vs 347, P=8–50 vs P=8–200, Rbr=9.2 vs 9.8). This is not necessarily a scientific flaw but it substantially complicates review and must be fixed. The fixed-component-mean issue in §4.2 is the main scientific concern: it directly affects the headline Jovian high-e fraction and break radius, so it should be addressed with a flexible-mean model rather than only discussed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nRead this one if you care about warm-planet eccentricities. It's the first homogeneous TESS-wide photoeccentric study spanning sub-Neptunes to Jovians, and the new continuous hierarchical mixture model—with radius uncertainties propagated through a sigmoid-weighted Beta mixture—is a real step beyond the discrete bins used by Dong et al. and the earlier Fairnington work. The paper also does careful light-curve fitting, SED-based stellar densities, a 30% random-drop robustness check, and it produces a concrete list of small high-eccentricity candidates worth follow-up. That is genuine value.\n\nThe soft spots are in the interpretation. First, the abstract numbers don't match the body: the posted arXiv abstract gives N=219, P=8–50 d, different component means and break radius, and a 2.7σ bimodality, while the full-text abstract gives N=347, P=8–200 d, and >4σ. That needs fixing before anything else.\n\nMore substantive: the continuous model forces the low-e component mean μlow to be radius-independent. The discrete Jovian-only fit gets μlow ≈ 0.225, but the continuous model imposes a global μlow ≈ 0.07. So moderate-e Jovians (e ~ 0.2–0.3) that the discrete fit places in the low-e component get reassigned to the high-e component, inflating the reported 59% high-e fraction. The break radius at ~9.8 R⊕ is likewise a property of the sigmoid parameterization; the paper does not test a model where μlow drifts with radius, and it explicitly declines to compute model evidence. The \">4σ bimodality\" in the abstract is not supported by any evidentiary comparison—§5.1 says no Bayesian evidences are computed, and the discussion offers only a 4σ credible interval on w2 > 0, which is a different claim. The qualitative trend—larger planets more eccentric—survives; the quantitative Jovian high-e fraction should not be taken at face value.\n\nThe grazing-transit cut (47 targets with >50% posterior mass at b ≥ 0.9) is also untested. High-impact-parameter transits can mimic the shortened durations associated with eccentricity, and unlike the 13 unconstrained-e systems, this cut receives no sensitivity check. This is a moderate concern, not fatal, but it should be addressed.\n\nVerdict: send to a serious referee. The method is original, the sample is valuable, and the issues are fixable with a radius-dependent μlow test, a proper model comparison, and a cleaner abstract. I would cite it with caveats.","headline":"A genuinely useful new TESS-wide eccentricity–radius analysis, but the headline high-e Jovian fraction depends on a modeling choice the paper doesn't defend.","tokens_in":25515,"tokens_out":2940,"would_cite":true,"duration_ms":27278,"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":"Warm single planets split into two eccentricity families, with the excited share rising steeply above about ten Earth radii.","keywords":["eccentricity distribution","warm Jupiters","photoeccentric effect","hierarchical Bayesian analysis","planet radius","TESS","single-transiting planets","transit duration"],"falsifier":"Re-run the hierarchical mixture including the 47 excluded grazing-transit candidates with an explicit high-impact-parameter model; if their inclusion pulls the Jovian high-e fraction down toward the sub-Saturn level or moves the break radius above 12 Earth radii, the two-component/transition claim fails. Alternatively, measure eccentricities for a sample of 30–50 of these warm Jovians by radial velocity; if the RV eccentricities are predominantly below 0.3, the photoeccentric high-e mode is an artifact.","tokens_in":24428,"feed_emoji":"🪐","tokens_out":6568,"duration_ms":63730,"temperature":0.7,"pith_summary":"The paper argues that the eccentric orbits of warm (8–200 day) single-transiting planets are not drawn from one smooth population. Instead, a hierarchical analysis of 347 TESS planets finds two components: a dominant near-circular mode and a separate dynamically excited mode with mean eccentricity above 0.6. The fraction of planets in the excited mode rises steeply with planet radius, crossing over at a break radius of about 9.8 Earth radii. Below that size, sub-Neptunes and sub-Saturns mostly stay on the low-eccentricity track; above it, warm Jovians are frequently eccentric, with roughly 59% assigned to the high-e mode. This matters because eccentricity encodes whether a planet formed quietly in a disk or was later scattered or perturbed.","feed_headline":"Two eccentricity families emerge among warm giant planets","feed_subtitle":"In 347 TESS planets, the eccentric share jumps from 15% for small worlds to 59% for giants.","key_machinery":"The central object is the photoeccentric effect: the ratio of the stellar density derived from the transit light curve under a circular-orbit assumption (the pseudo-density) to the true stellar density from SED/isochrone fitting yields, via g(e,ω) = (1+e sinω)/sqrt(1−e^2), a joint constraint on eccentricity e and argument of periastron ω. The population-level engine is a three-stage hierarchical Bayesian model: a two-component Beta mixture for eccentricity, with the mixture weight set by a logistic sigmoid function of planet radius, so the transition radius Rbr and the low-e/high-e component means are inferred simultaneously while propagating radius uncertainties.","core_discovery":"On the paper's own terms: using the photoeccentric effect—where a transiting planet's light-curve duration betrays its orbital speed and hence its eccentricity—the authors infer the population-level eccentricity distribution of 347 warm single-planet systems observed by TESS. A two-component Beta mixture, with membership governed by a logistic sigmoid of radius, cleanly separates the population into a low-eccentricity component (mean e ~ 0.07) and a high-eccentricity component (mean e ~ 0.62). The high-e fraction grows from about 15% for 1–4 Earth-radii planets to about 59% for 8–16 Earth-radii Jovians, with a transition at Rbr = 9.8+1.4−1.1 Earth radii. The authors interpret this as evidenc","pith_inferences":["If the 47 removed grazing-transit systems were included, the high-eccentricity fraction could shrink: grazing transits shorten durations and mimic eccentric orbits, so the reported 59% Jovian high-e share and the 9.8 Earth-radii break may be upper bounds rather than intrinsic values.","A natural next test is to apply the same sigmoid mixture to warm multi-planet systems, which should show a much smaller high-e fraction; that would confirm that the high-e mode is tied to single-planet architectures, as the paper expects.","The radius-continuous model predicts a specific, testable conditional distribution: for any newly discovered warm single planet with a measured radius, the probability it belongs to the high-e mode is a smooth function of radius; future TESS samples can check whether the predicted logistic curve reproduces out-of-sample eccentricities."],"forward_implications":["If the claim holds, warm Jupiters are not uniformly circularized; a majority have eccentricities around 0.6, so high-eccentricity migration or scattering must be common for giant planets at 8–200 day periods.","Sub-Saturns (4–8 R⊕) behaving like sub-Neptunes argues against the idea that most sub-Saturns are simply failed gas giants with giant-planet-like dynamics; their dynamical history resembles smaller planets.","The existence of a ~15% eccentric sub-Neptune tail implies a population of small planets excited by unseen companions, predicting detectable outer companions or transit-timing variations around those systems.","The break radius near 10 R⊕ provides a target for formation models: whatever process pumps eccentricity must switch on sharply between 4 and 16 Earth radii.","The observed distribution is transit-selected; accounting for detection completeness will shift the intrinsic distribution, but the qualitative rise of eccentricity with radius should persist."],"fun_headline_variants":["Warm giants mostly eccentric, small planets mostly circular","Eccentricity jumps from 16% to 65% with planet size","Two eccentricity families: calm small, excited large","Warm singles: eccentric share 16% small, 65% giants","Planet size sets eccentricity: small calm, giants eccentric"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The central result relies on removing 47 planets whose transits are almost grazing (impact parameter at least 0.9), because such transits can look like eccentric orbits; if those planets are actually eccentric, the measured high-eccentricity fraction and break radius would be off, and the paper does not test this.","fun_headline_variants_meta":{"raw":{"variants":["Warm giants mostly eccentric, small planets mostly circular","Eccentricity jumps from 16% to 65% with planet size","Two eccentricity families: calm small, excited large","Warm singles: eccentric share 16% small, 65% giants","Planet size sets eccentricity: small calm, giants eccentric"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001938,"raw_usage":{"total_tokens":7497,"prompt_tokens":902,"completion_tokens":6595,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":646,"completion_tokens_details":{"reasoning_tokens":6506}},"tokens_in":646,"tokens_out":6595,"duration_ms":42668,"temperature":1.0,"reasoning_tokens":6506,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T21:25:00.290775+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the hierarchical mixture including the 47 excluded grazing-transit candidates with an explicit high-impact-parameter model; if their inclusion pulls the Jovian high-e fraction down toward the sub-Saturn level or moves the break radius above 12 Earth radii, the two-component/transition claim fails. Alternatively, measure eccentricities for a sample of 30–50 of these warm Jovians by radial velocity; if the RV eccentricities are predominantly below 0.3, the photoeccentric high-e mode is an artifact.","supporting_citations":[],"review_version":1}