{"id":"bd862569-b663-42cb-a3d7-d2d091961362","arxiv_id":"2507.07169","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"M dwarf planets larger than about 3.5 Earth radii show elevated orbital eccentricities, while smaller planets remain on nearly circular orbits, a pattern similar to planets around Sun-like stars.","lead":"This paper measures the orbital eccentricities of 236 small planets around M dwarf stars using transit shapes, and finds that planets larger than about 3.5 Earth radii have more stretched-out orbits than smaller planets. The pattern matches planets around Sun-like stars, which helps scientists understand how close-in planets form and lose their atmospheres.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central transition claim lacks a formal significance test; the observed sigmoid may be consistent with a flat line given the small high-radius sample.","rationale":"The reader's weakest assumption is the monotonic Beta distribution (α<1, β>1) for the underlying eccentricity distribution. I agree this is a modeling constraint, but it is partially mitigated by the paper's own empirical histogram model (Appendix B, Figure 12), which is presented as an alternative way to compute ⟨e⟩ and would presumably show the same radius trend if the Beta assumption were driving the result. The more load-bearing gap is the absence of any significance test for the global transition itself. The central claim—that there is 'marked evidence' for a transition near 3.5 R⊕—is supported only by a sigmoid fit whose preference over a flat line is never quantified. The paper does quantify the radius-gap peak significance (A>0 for 69% of posterior samples, ~1σ), but no equivalent statistic is reported for the transition. Given that the large-radius bins contain only 32 planets and are dominated by TESS single-transit detections, the possibility that the apparent transition arises from noise or small-sample fluctuations is real and testable. This does not change the conditional verdict—rather, it sharpens the condition: the paper should either add a formal model comparison or soften the 'marked evidence' language. A minor additional inconsistency is that Section 3.4 writes 'α<1 and β<1' where Section 3.2.1 correctly specifies 'α<1 and β>1'; this typo should be fixed but is not load-bearing. If the proposed significance test shows the sigmoid is strongly preferred, the central claim is validated and the paper would be suitable for acceptance after minor revisions.","tokens_in":25327,"tokens_out":6931,"duration_ms":82530,"concrete_test":"Using the released eccentricity posterior samples, recompute the five coarse-bin ⟨e⟩ values and their 16th/84th percentiles. Fit two models to these five points with a Gaussian likelihood: (i) a constant ⟨e⟩ = B, and (ii) the logistic sigmoid of Eq. 10 with uniform priors on B, L, k, and x_t. Report Δχ² and the Akaike/Bayesian information criterion, or compute a Bayes factor via nested sampling. If Δχ² < 4 or the sigmoid is not preferred at >2σ, the 'marked evidence' claim should be downgraded to 'consistent with a transition but not significant.' Repeat using the Appendix B empirical histogram ⟨e⟩ values (Figure 12) to verify the result is not an artifact of the Beta assumption.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 4.2 reports the logistic sigmoid fit (Eq. 10) to binned ⟨e⟩ and quotes a transition radius of 3.1+1.5/−1.2 R⊕ with e_high/e_low = 4.6+5.8/−1.9, but provides no statistical comparison to a null model with constant ⟨e⟩. The high-radius evidence rests on only 14 planets in the 3.5–7.5 R⊕ bin and 18 in the 7.5–16 R⊕ bin (Table 3), the latter composed entirely of single-transit TESS planets. With five bins and asymmetric error bars, a steep sigmoid can fit noise. The paper applies a flat-line null only to the radius-gap peak (Eq. 12, yielding ~1σ for singles), not to the global transition. Absent a Δχ², Bayes factor, or posterior-predictive p-value, the Conclusions statement 'marked evidence for a transition' is not quantitatively supported.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constrains orbital eccentricities for 236 transiting planets orbiting M dwarfs using the photoeccentric effect on TESS and Kepler light curves, then applies a hierarchical Bayesian model with Beta distributions to infer the underlying eccentricity distribution in radius bins. The authors report a transition from low to high eccentricity at approximately 3.5 R_Earth, with e_high/e_low about 4.6, and interpret this as analogous to the FGK-dwarf eccentricity-radius relation found by G25. They also investigate eccentricity near the radius gap, finding only about 1-sigma evidence for elevated eccentricities among single-transit planets and no evidence among multi-transit planets, and discuss implications for photoevaporation versus giant-impact atmospheric loss.","tokens_in":25576,"tokens_out":8242,"duration_ms":85556,"significance":"If the transition is real, this would be the first clear demonstration that the eccentricity-radius relation for M-dwarf planets mirrors that of FGK dwarfs, implying common formation and evolution channels across spectral types. The paper is methodologically careful: it validates the importance-sampling transit fit against the prior Kepler-based analysis (Appendix A), uses an external stellar density prior, and directly compares to G25 with identical binning. The main limitation is that the central transition claim lacks a formal significance test, and the robustness of the result to the assumed Beta distribution shape is not quantified. With a formal null-model comparison and additional robustness checks, this would be a valuable contribution to the demographics of small exoplanets.","major_comments":[{"comment":"The logistic sigmoid fit to the five binned <e> points is presented without any statistical comparison to a null model with constant <e> across radius. The paper only applies a flat-line test to the Gaussian radius-gap peak (Eq. 12), not to the global transition. Given that the two highest-radius bins contain only 14 and 18 planets, respectively (Table 3), and the 7.5-16 R_Earth bin is entirely composed of single-transit TESS planets, a steep sigmoid can easily fit noise. Please add a formal significance test, such as a Delta-chi^2 or Delta-BIC comparison between the sigmoid and a constant model, or a posterior predictive p-value derived from the hierarchical model, and report the resulting significance. As written, the Conclusions statement 'We show marked evidence for a transition' is not quantitatively supported.","section":"Section 4.2, Eq. (10)"},{"comment":"The assumed Beta distribution with alpha<1 and beta>1 is monotonically decreasing and peaks at e=0; it cannot represent a separate high-eccentricity population or a bump away from zero. The authors acknowledge they 'lack sufficient physical understanding of the true shape' but adopt this form for reproducibility and efficiency. The empirical histogram model in Appendix B is used only for qualitative comparison (Figure 12), and the transition parameters of Eq. (10) are not re-derived from that model. The robustness of the central transition claim to this modeling choice is therefore untested. Please demonstrate that the Beta assumption does not bias the binned <e> values, for example by fitting a two-component mixture or by quantifying the transition with the empirical histogram model and showing that the result is unchanged.","section":"Section 3.2.1"},{"comment":"The combined-sample transition may be driven by the survey and multiplicity composition of the high-radius bins: the 7.5-16 R_Earth bin contains only single-transit TESS planets, and single-transit systems are known to have higher eccentricities (Section 4.1). The paper compares singles and multis in Figure 3, but does not test whether the sigmoid transition is present within the single-transit subsample alone or after controlling for multiplicity and survey. The physical interpretation in Section 5.1 assumes a radius-driven effect, but a multiplicity-driven selection effect is a competing explanation. Please add a test that isolates the radius effect from the multiplicity effect, such as fitting the sigmoid to the single-transit sample only, or including multiplicity as a covariate in the hierarchical model.","section":"Section 4.2, Table 3"}],"minor_comments":[{"comment":"The text states that the Beta distribution is modeled with alpha < 1 and beta < 1, which is inconsistent with Section 3.2.1, where beta > 1 is required for a monotonically decreasing distribution peaking at e=0. This appears to be a typo and should be corrected.","section":"Section 3.4"},{"comment":"The 'Gelman-Ruban statistic' should read 'Gelman-Rubin statistic.'","section":"Section 3.1.1"},{"comment":"The sentence 'We take the mean of these median values as <e>' is confusing: it is unclear whether the authors take the mean or the median of the posterior draws of the Beta distribution's mean, and the wording conflates the two operations.","section":"Section 3.4"},{"comment":"The binned <e> values and the sigmoid fit parameters are only presented in figures; a machine-readable table with the numerical values of <e>, their uncertainties, and the sigmoid parameters (B, L, k, x_t) would allow readers to reproduce the test requested in the first major comment.","section":"Section 4.2, Figure 2"},{"comment":"The abstract and conclusions quote the transition at 3.5 R_Earth, while the fitted transition in Section 4.2 is 3.1(+1.5/-1.2) R_Earth; please clarify that 3.5 R_Earth is the bin edge rather than the fitted transition location, or quote the fitted value consistently throughout.","section":"Abstract and Section 6"}],"recommendation":"major_revision","confidential_remarks":"The paper is a direct extension of the authors' prior work (S23 and G25), and the novelty lies in applying the same framework to M-dwarf planets. The method is carefully described and the validation against the Kepler sample is a strength. The main issue is the unsupported strength of the central claim: the transition is asserted as 'marked evidence' without a formal significance test, and the model-choice limitation is acknowledged but not addressed quantitatively. These are fixable within the scope of a revision. There is no concern about circularity; the external stellar density prior and the G25 comparison are appropriate. I recommend major revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Fairly straightforward extension of the G25 eccentricity-radius analysis to M dwarfs, with 79 new TESS planets added to the earlier S23 Kepler sample. The core measurement—a transition from low to high eccentricity somewhere around 3–4.5 R⊕—is plausible, and the authors are careful to make it directly comparable to G25: same bins, same Beta model, direct comparison of binned curves, plus a non-parametric histogram check in the appendix. The importance sampling pipeline is validated against S23 on one Kepler planet and against a TESS planet with matched SNR. That is honest, reproducible work.\n\nWhere it wobbles: there is no significance test for the sigmoid transition. The fit gives x_t = 3.1+1.5/−1.2 R⊕ and e_high/e_low = 4.6+5.8/−1.9, but five binned points with only 14 and 18 planets in the two large-radius bins can look like a step even if the truth is flat. The paper never compares the sigmoid to a constant-⟨e⟩ null. The Conclusions statement 'marked evidence for a transition' is not supported by the statistics as presented. It may be true, but the evidence is suggestive, not marked. The radius-gap eccentricity peak for singles is 69% positive A, i.e., ~1σ, which the body correctly calls modest, though the abstract leans on it a bit.\n\nThe Beta-distribution assumption (α<1, β>1) is a modeling constraint, but the appendix's empirical histogram model produces the same qualitative relation, so that concern is not fatal. Selection effects from the roughly two dozen excluded TOIs are acknowledged but not quantified; probably minor.\n\nBottom line: a demographic step forward for M dwarfs, squarely in line with G25, with stronger cross-checks than typical for this kind of analysis. The missing null test and the thin high-radius sample are the two things a referee should push on. I would send it to peer review with a request to add a flat-line comparison and soften the conclusions. I'd cite it for the TESS eccentricity posteriors even now, and I'd bring it to reading group to debate whether the step is real.","headline":"Plausible extension of the eccentricity-radius relation to M dwarfs, but the headline transition lacks a null-model significance test and rides on a thin high-radius sample.","tokens_in":26095,"tokens_out":2969,"would_cite":true,"duration_ms":31844,"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":"Using transit photometry for 236 planets orbiting M dwarfs, this paper reports a sharp transition from low to high orbital eccentricity at about 3.5 Earth radii, mirroring the relation previously found for planets around Sun-like stars.","keywords":["exoplanets","orbital eccentricity","M dwarf stars","radius gap","photoeccentric effect","transit photometry","TESS","Kepler"],"falsifier":"Re-fit the same 236 planets with a two-component eccentricity model that allows a secondary peak away from e=0 and compare the inferred transition; if the ~3.5 Earth radii step weakens or moves, the reported relation is an artifact of the Beta prior. The observational counterpart is a radial-velocity survey of M dwarf planets straddling 3.5 Earth radii, which would show whether eccentricities really are higher above that radius than below it.","tokens_in":25135,"feed_emoji":"🪐","tokens_out":11537,"duration_ms":102224,"temperature":0.7,"pith_summary":"Using transit photometry for 236 confirmed planets and candidates orbiting M dwarf stars, this paper asks whether the orbital eccentricity-radius relation seen for planets around Sun-like stars also holds around the most common stars in the galaxy. The authors extract individual eccentricity constraints from transit shapes together with a stellar-density prior, then infer the underlying eccentricity distribution in radius bins with a hierarchical Bayesian model. They report a clear transition from low to high eccentricity at about 3.5 Earth radii, with planets above that size roughly five times more eccentric on average than smaller planets. They also find no evidence that M dwarf planets near the radius gap have elevated eccentricities, which they interpret as supporting photoevaporation rather than giant impacts as the main atmospheric-loss mechanism for these planets.","feed_headline":"236 M dwarf planets reveal an eccentricity jump at 3.5 Earth radii","feed_subtitle":"The same low-to-high eccentricity jump seen for Sun-like stars appears at 3.5 Earth radii around M dwarfs.","key_machinery":"The photoeccentric effect is the load-bearing technique: a planet's transit duration and ingress/egress shape depend on its speed across the stellar disk, which for a known stellar density maps to a joint constraint on eccentricity e and argument of periastron ω. The paper implements it by fitting transit light curves with duration as a free parameter and then importance-sampling the (e, ω) posterior against a stellar-density prior, requiring the density inferred from the transit (assuming a circular orbit) to match the independently known stellar density. Population-level inference is carried out by a hierarchical model that assumes the parent eccentricity distribution in each radius bin is a Beta distribution constrained to peak at e=0 and decrease monotonically (alpha<1, beta>1), reparametrized to avoid prior biases. A logistic sigmoid is then fitted to the binned mean eccentricities to locate and characterize the low-to-high-eccentricity transition.","core_discovery":"The paper's central claim is that the eccentricity-radius relation for M dwarf planets is a rising step function: small planets orbit on nearly circular paths while planets larger than about 3.5 Earth radii have systematically higher eccentricities. Fitting a logistic sigmoid to the binned mean eccentricities locates the transition at 3.1+1.5/-1.2 Earth radii, with the ratio of high- to low-eccentricity levels at 4.6+5.8/-1.9, consistent within 1σ with the transition measured for planets orbiting FGK dwarfs. The authors further claim that, unlike the FGK case, M dwarf planets near the radius gap show no significant elevation in eccentricity: multi-transit planets stay at low eccentricity at all radii, and the modest excess seen for single-transit planets between 1.9 and 3.0 Earth radii is consistent with a flat line at roughly 1σ. This asymmetry, if physical, is taken as evidence that photoevaporation or core-powered mass loss, rather than giant impacts, dominate atmospheric stripping for M dwarf planets.","pith_inferences":["A testable extension would be to run the same hierarchical inference on a larger joint Kepler-TESS sample; if the transition stays pinned near 3.5 Earth radii as the sample grows, the photoevaporation-versus-giant-impact interpretation would gain strength.","The authors leave implicit that their Beta-distribution model cannot represent a separate high-eccentricity population, so a re-analysis with a two-component mixture would test whether the e_high/e_low ratio and transition location are artifacts of the assumed distributional shape.","A prediction of their interpretation, not stated in the paper, is that radial-velocity eccentricities of M dwarf planets above roughly 3.5 Earth radii should be systematically higher than those below, and that the single-transit radius-gap excess, if real, should grow with a larger sample.","Because TESS single-transit systems are likely contaminated by undetected multi-planet systems, the reported single-versus-multi eccentricity contrast could sharpen as longer TESS baselines reveal hidden companions, changing the inferred radius-gap behaviour for singles."],"forward_implications":["If the transition is real, planet formation around M dwarfs and FGK dwarfs produces two distinct radius regimes, with the boundary near 3.5 Earth radii, and the dynamical evolution of larger planets leaves them on eccentric orbits regardless of host-star mass.","The absence of elevated eccentricities at the radius gap for multi-transit M dwarf planets would point to photoevaporation or core-powered mass loss, mechanisms that remove atmospheres without changing orbital dynamics, as the dominant sculptors of the M dwarf radius valley.","The consistency of the transition radius between M dwarf and FGK dwarf planets (3.1+1.5/-1.2 versus 3.3±0.4 and 4.2±0.9 Earth radii) suggests the boundary between rocky and gas-rich planet formation channels is set by planet properties rather than by host star mass.","If single-transit M dwarf planets near the radius gap do have modestly elevated eccentricities, giant impacts or planet-planet scattering may still play a role in atmospheric loss for dynamically hot, single-planet systems."],"supporting_citations":[{"why":"It supplies the FGK dwarf eccentricity-radius relation this work extends and the logistic sigmoid fitting procedure used to locate the transition.","marker":"G25"},{"why":"It provides the Kepler M dwarf planet sample, its eccentricity posteriors, and the hierarchical inference framework reused here.","marker":"S23"},{"why":"It established the photoeccentric effect formalism used to turn transit durations into eccentricity constraints.","marker":"R. I. Dawson & J. A. Johnson (2012)"},{"why":"It supplies the importance-sampling method used to generate individual (e, omega) posteriors from transit fits.","marker":"M. G. MacDougall et al. (2023a)"},{"why":"It gives the eccentricity-dependent transit probability correction applied in the hierarchical likelihood.","marker":"D. M. Kipping (2014)"},{"why":"It defines the period-radius diagonal bins used to place 'radius gap' planets in the adjusted-radius analysis.","marker":"C. S. K. Ho et al. (2024)"},{"why":"It provides simulations showing that giant-impact-driven atmospheric loss produces elevated eccentricities at the radius gap, the interpretive link used in the discussion.","marker":"Q. Chance & S. Ballard (2024)"},{"why":"It supplies the revised TIC stellar densities that act as the density prior for TESS host stars in the photoeccentric measurement.","marker":"K. G. Stassun et al. (2019)"}],"fun_headline_variants":["M dwarf planets get eccentric above 3.5 Earth radii","Eccentricity jumps at 3.5 Earth radii for M dwarfs","Same eccentricity step as Sun-like stars seen for M dwarfs","Larger M dwarf planets orbit on more eccentric paths","M dwarf planets show step in eccentricity at 3.5 Earth radii"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The chain of inference assumes that in every radius bin the true eccentricity distribution has exactly one shape—highest at zero eccentricity and falling monotonically—so a population of moderately eccentric planets sitting away from zero would be invisible to the model and could bias the reported average eccentricities and the location of the transition.","fun_headline_variants_meta":{"raw":{"variants":["M dwarf planets get eccentric above 3.5 Earth radii","Eccentricity jumps at 3.5 Earth radii for M dwarfs","Same eccentricity step as Sun-like stars seen for M dwarfs","Larger M dwarf planets orbit on more eccentric paths","M dwarf planets show step in eccentricity at 3.5 Earth radii"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000596,"raw_usage":{"total_tokens":2836,"prompt_tokens":1038,"completion_tokens":1798,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":654,"completion_tokens_details":{"reasoning_tokens":1706}},"tokens_in":654,"tokens_out":1798,"duration_ms":13488,"temperature":1.0,"reasoning_tokens":1706,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T18:47:25.307282+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-fit the same 236 planets with a two-component eccentricity model that allows a secondary peak away from e=0 and compare the inferred transition; if the ~3.5 Earth radii step weakens or moves, the reported relation is an artifact of the Beta prior. The observational counterpart is a radial-velocity survey of M dwarf planets straddling 3.5 Earth radii, which would show whether eccentricities really are higher above that radius than below it.","supporting_citations":[],"review_version":1}