{"id":"574c220b-d3fb-441d-9a57-997e9e557394","arxiv_id":"2507.15359","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Using a new Gauss-Legendre quadrature transit model and 800 injection-retrieval tests, the authors find Saturn-like oblateness is detectable at 3sigma in favorable cases, while planets flatter than Jupiter are not reliably retrieved at any tested noise level.","lead":"This paper tests whether space telescopes can measure the flattening (oblateness) of rotating exoplanets from the shape of their transit light curves, using a new Gauss-Legendre quadrature model and 800 simulated light curves. It finds that Saturn-like oblateness is detectable around bright stars with precise stellar density knowledge, but that noise above 256 ppm makes detection unreliable.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The stellar-density prior that breaks the b-f-obliquity degeneracy is stated as sigma=0.24 in the text but N(73.26,0.04) in Table 7; if the narrower value was used in the 800 retrievals, the 3-sigma detectability claim is inflated by a factor of six.","rationale":"I agree with the reader's conditional verdict. The reader's weakest assumption was the combined fragility of the stellar-density prior, circular-orbit restriction, and limb-darkening priors; my attack isolates the internal inconsistency in the prior width, which is even more actionable because it can be settled from the released files. I do not claim any misrepresentation; the discrepancy is likely a typo, but the paper must disambiguate it because the whole retrieval strategy leans on this prior. The self-injection design is a secondary concern: the 800 injection tests use the same forward model for injection and retrieval, so they validate the statistical machinery but not the model's fidelity to a real oblate planet; the independent squishyplanet comparison in Sect. 3.3 shows a 13 ppm residual, a nontrivial fraction of the roughly 45 ppm oblateness signal. That comparison bounds the model error rather than eliminating it, and the prior-width ambiguity is the sharpest falsifiable point. I recommend keeping the conditional verdict, with the revision requiring the authors to state the exact sigma used, provide a retrieval with sigma=0.24, and clarify the circular-orbit scope in the abstract.","tokens_in":23834,"tokens_out":5810,"duration_ms":67317,"concrete_test":"Download the released light-curve and configuration files (OSF link in Acknowledgements), identify the actual sigma used for the a/R* Gaussian prior in the DE-MCMC setup, and rerun a subset of the grid (e.g., f=0.09, obliquity=9/45/81 deg, sigma_w=32 and 64 ppm, with and without red noise) replacing sigma=0.04 by sigma=0.24. Then recompute the 3-sigma detection count on those subsets. If detection rates and f/Delta-f distributions are unchanged within sampling noise, Table 7 is a harmless typo and the claim stands; if rates drop materially, the abstract should be reworded to state the required stellar-density precision is about 0.2% rather than the 1% claimed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 4.2 motivates a Gaussian prior on a/R* by assuming an asteroseismic stellar density known to 1%, giving Delta(a/R*)/(a/R*)=0.0033 and therefore sigma(a/R*)=0.24 for the fiducial mean 73.26. The text states this value, but Table 7 lists N(73.26,0.04), six times tighter. The width is not cosmetic: Fig. 17 shows Pearson r=0.988 between a/R* and b in the posterior, so the a/R* prior is the main lever that breaks the b-f-obliquity degeneracy identified in Sect. 4.2. If the retrievals actually used sigma=0.04, they implicitly assume a stellar density precision of about 0.16%, far better than the 0.5-2.6% range cited from Silva Aguirre et al. (2017), and the headline 3-sigma detection rates (e.g., 59% for Saturn-like f about 0.09) are optimistic. If they used sigma=0.24, the detection rates and f/Delta-f contours in Figs. 12-15 may degrade substantially at low f; the manuscript does not report which value was used or give a sensitivity run. In either case the central 3-sigma claim is not yet pinned down. The paper is otherwise transparent: it openly notes the circular-orbit restriction of this prior trick, so eccentric orbits are a separate, already-flagged limitation, not the primary blocker.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a new numerical transit model for oblate (biaxial ellipsoid) exoplanets based on Gauss-Legendre quadrature, implemented in the TLCM framework. The model is benchmarked against the Mandel-Agol analytical model in the spherical limit and against the squishyplanet model for one oblate configuration. The authors then perform 800 injection-and-retrieval tests with different oblateness values f, sky-projected obliquities ϑ, white-noise levels, and injected red noise, and they report 3σ oblateness detection rates as a function of f, ϑ, and noise. Their headline claims are that Saturn-like oblateness (f≈0.09) is detectable in about 59% of tested configurations around bright stars when the stellar density is known from asteroseismology, and that noise levels of 256 ppm or higher make oblateness detection unreliable. The paper also introduces the Δ metric for quantifying light-curve model differences and reports a modest speed advantage over the analytical model.","tokens_in":24114,"tokens_out":4154,"duration_ms":48866,"significance":"If the detectability results are robust, the paper provides a practical roadmap for constraining sky-projected oblateness with CHEOPS, JWST, PLATO, or ARIEL photometry, and the publicly archived input/output light curves support reproducibility. The paper is also transparent about a key limitation: the stellar-density prior that breaks the b–f–ϑ degeneracy is stated to be valid only for circular orbits. However, the central 3σ detectability claim is not yet pinned down because of a quantitative inconsistency in the adopted prior width, and the external validation of the oblate model is thin. The work is potentially valuable but currently requires revision before the headline numbers can be taken at face value.","major_comments":[{"comment":"The text in Section 4.2 states that a Gaussian prior on a/R⋆ is applied with mean 73.26 and standard deviation 0.24, justified by a 1% stellar-density precision via Eq. (30), while Table 7 lists the prior as N(73.26, 0.04), which is six times narrower. This is not cosmetic: the narrower width corresponds to Δρ⋆/ρ⋆ ≈ 0.16%, far better than the 0.5–2.6% range cited from Silva Aguirre et al. (2017), and Fig. 17 shows a very strong a/R⋆–b correlation (Pearson r = 0.988). Because the a/R⋆ prior is the main lever that breaks the b–f–ϑ degeneracy described in Section 4.2, the reported 3σ detection rates (e.g., 59% for f≈0.09) and the claim that f≥0.15 is 'guaranteed' at low noise depend directly on which width was used. The manuscript must state which value was actually used in the 800 retrievals, reconcile the text and table, and provide a sensitivity run at the wider prior.","section":"Section 4.2 and Table 7"},{"comment":"The only external benchmark for the oblate model is a single configuration in which TLCM and squishyplanet differ by roughly 13 ppm in amplitude, while the estimated oblateness signal is roughly 45 ppm. The systematic model difference is therefore about 29% of the signal amplitude, and it is comparable to the lowest tested white-noise levels (σw = 2–16 ppm). Since the injection and retrieval steps both use the same TLCM Gauss-Legendre integrator, systematic integration errors partly cancel in the recovery statistics, making the 3σ detection claim model-conditioned. I recommend adding comparisons over at least a small grid of f and ϑ values, or otherwise quantifying the bias that the 13 ppm-level model difference would induce in the retrieved f and in the detection rates.","section":"Section 3.3"},{"comment":"The paper explicitly notes that the stellar-density prior trick is valid strictly for circular orbits, and all 800 injection tests use circular orbits. The abstract's unqualified statement that a 3σ oblateness detection is possible for a planet orbiting a bright enough star therefore overstates the scope: eccentric orbits with unknown eccentricity and argument of periastron are not covered by the reported detection rates. Since the limitation is already flagged in Section 4.2, this is a scope-and-presentation issue rather than an internal inconsistency, but the abstract should carry the circular-orbit qualifier or otherwise clearly state the restricted applicability of the method.","section":"Abstract and Section 4.2"}],"minor_comments":[{"comment":"There are several typographical errors that should be corrected in a revision, including 'oblanteness' in Section 4.1, 'charaterized' in Section 4.3, 'lowe' in Section 3.3, and 'the the model' in Section 4.1.","section":"Throughout"},{"comment":"The reported speed advantage of the numerical model over the analytical model is inconsistent: the abstract says 'about 25% faster,' Section 3.2 says the analytical model takes about 28% more time, and the conclusion says 'about 20% faster.' These numbers should be harmonized.","section":"Abstract, Section 3.2, Section 5"},{"comment":"The quadratic coefficient in the piecewise fit for log Δ versus log(a/R⋆) is given as −0.314(5) in Eq. (26) but as −0.134(5) in the caption of Fig. 5. Please check which value is correct and make the two consistent.","section":"Eq. (26) and Fig. 5"},{"comment":"The table uses the notation N(73.26, 0.04) for the a/R⋆ prior, while the text in Section 4.2 uses a standard deviation of 0.24; beyond resolving this discrepancy, the table would benefit from a footnote stating the corresponding assumed stellar-density precision.","section":"Table 7"},{"comment":"The limb-darkening priors are described as 'strict' with standard deviation 0.01, but the text in Section 4.1 refers to limb-darkening coefficients 'A = B = 1.3' while Table 7 uses uA and uB; please unify the notation.","section":"Section 4.2"}],"recommendation":"major_revision","confidential_remarks":"I do not see grounds for rejection: the model is useful, the test campaign is extensive, and the authors are unusually explicit about some limitations. The decisive issue is the prior-width inconsistency, which directly controls the headline detectability numbers. If the 800 retrievals actually used the narrower prior from Table 7, the 3σ rates are likely optimistic when recast at the text's 0.24 width; if they used the wider prior, the authors should demonstrate that the reported rates survive a correct tabulation. I would also encourage the editor to weigh the limited external benchmark coverage carefully, since the injection-retrieval design cannot expose systematic errors shared by injection and fitting. A sensitivity run rather than a full rerun of all 800 cases may be sufficient to resolve the prior question, but the external validation concern likely needs additional model comparisons."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a useful, honest paper on single-transit oblateness retrieval, but the central 3-sigma detectability claim is not yet pinned down because the a/R* prior that breaks the degeneracy is stated two different ways (sigma=0.24 in Section 4.2, N(73.26,0.04) in Table 7), and the paper doesn't say which was used or give a sensitivity run. The stress-test note is right; this is the main blocker.\n\nWhat's new: the Gauss-Legendre quadrature light curve model, the systematic 8x10x10 injection-retrieval grid over oblateness, obliquity and noise, and the practical noise floor of 256 ppm. The delta metric for light-curve model comparison is a nice idea. The paper is admirably transparent: Section 2.3 admits that the implemented version doesn't strictly follow GL quadrature premises, and the compliant version is two orders of magnitude slower. It also flags the circular-orbit restriction of the density prior trick. The authors ship the light curves and configs at OSF, which is more than most papers do.\n\nThe soft spots, in order. First, the prior inconsistency above. It matters because Fig. 17 shows a/R* and b are nearly perfectly correlated (r=0.988), so the a/R* prior is the main lever that breaks the b-f-obliquity degeneracy. If the retrievals used sigma=0.04, they assume a stellar density precision of about 0.16%, much better than the 0.5-2.6% from Silva Aguirre et al. cited in the text; if they used 0.24, the detection rates in Figs. 12-15 may degrade. The central claim could change either way. Second, the injection-retrieval uses the same TLCM integrator for both, so systematic integration error partially cancels; the external checks give residuals of ~3.3 ppm vs Mandel-Agol in the spherical limit, but ~13 ppm vs squishyplanet in the oblate case, about 29% of the 45 ppm signal in that example. That's not disqualifying, but without a sensitivity run it leaves the 3-sigma statistics model-conditioned. Third, the abstract says '3sigma detection is possible' while the body reports 27.5% for Jupiter-like, 59% for Saturn-like, 10% for f=0.03; that's a wider spread than the abstract implies. The tight limb-darkening priors are a fourth assumption, but the authors state them openly.\n\nFor an observer planning CHEOPS/PLATO/ARIEL targets, the detectability map is genuinely useful as a first look. But I wouldn't cite the detection rates in the current form. The paper deserves peer review; the issues are addressable, but a sensitivity run on the a/R* prior width and a reanalysis of the detection statistics with the correct consistency check are needed before I'd trust the numbers.","headline":"Useful detectability map for single-transit oblateness, but the central 3-sigma claim rests on an unresolved prior-width inconsistency (0.24 vs 0.04) in the a/R* constraint; needs a sensitivity run before the numbers can be trusted.","tokens_in":24740,"tokens_out":2688,"would_cite":false,"duration_ms":27058,"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":"Space-based photometry can detect a Saturn-like exoplanet's rotational flattening in a single transit, provided the host star's density is known to about one percent.","keywords":["oblateness","exoplanets","transit photometry","rotational flattening","Gauss-Legendre quadrature","asteroseismology","light-curve retrieval","space-based photometry"],"falsifier":"Re-run the same 800-case injection grid with eccentric orbits (for example $e=0.05$, all other parameters and priors unchanged): if the $3\\sigma$ detection rate for Saturn-like oblateness falls substantially below 59%, or the recovered $f$ is biased by more than 0.02, the circular-orbit assumption is the part of the claim that fails.","tokens_in":23596,"feed_emoji":"🪐","tokens_out":13216,"duration_ms":123978,"temperature":0.7,"pith_summary":"The paper asks whether the rotational flattening of an exoplanet, the slight difference between its equatorial and polar radii, leaves a measurable imprint in a space-based transit light curve. The authors build a numerical transit model that uses Gauss-Legendre quadrature to integrate the overlap between a limb-darkened star and a sky-projected ellipse, verify it against an established spherical-transit model, and run 800 injection-and-retrieval experiments with realistic white and time-correlated noise. They find that a single transit of a Saturn-like planet (oblateness $f\\approx 0.09$) around a bright, quiet star yields a $3\\sigma$ detection of sky-projected oblateness in about 59% of tested configurations, provided the stellar mean density is known to about 1% from asteroseismology. They also establish a noise floor: at point-to-point scatter of 256 ppm per 60-second exposure, oblateness retrieval is not reliable, and oblateness values below Jupiter's ($f\\approx 0.06$) are not recoverable at any tested noise level.","feed_headline":"A single transit can reveal a Saturn-like planet's oblateness","feed_subtitle":"With the star's density known to 1%, 3-sigma detections succeed in 59% of tested configurations.","key_machinery":"The carrying mechanism is a numerical transit model in which the sky-projected planetary disk is an ellipse and the blocked stellar flux is computed by two-dimensional Gauss-Legendre quadrature (a numerical integration rule sampling the integrand at Legendre-polynomial roots), evaluating the limb-darkened stellar surface brightness at each quadrature point and summing only points inside the stellar disk. The load-bearing identity is Kepler's third law written as $(a/R_\\star)^3 = P^2 G(1+q)/(3\\pi\\rho_\\star)$, which converts an external stellar-density measurement into a Gaussian prior on the scaled semi-major axis and thereby breaks the $b$\\textendash$f$\\textendash$\\vartheta$ degeneracy that otherwise hides the oblateness signal. In the spherical limit the model matches the standard analytical transit light curve to within a few ppm at 96 quadrature points, while running about 25% faster, and its discrepancy with the analytical model shrinks as $n^{-2}$ with only a mild rise in runtime.","core_discovery":"The central claim is that a rotating exoplanet imprints its non-spherical shape on the ingress and egress of the transit light curve, and that this imprint can be recovered once the degeneracy with the impact parameter is broken. The paper shows that a precise stellar mean density, obtained for example from asteroseismology, supplies that break: Kepler's third law ties $a/R_\\star$ to $\\rho_\\star^{1/3}$, so a 1% density measurement becomes a Gaussian prior of $\\sigma = 0.24$ on $a/R_\\star$, which in turn pins down the transit geometry and separates oblateness $f$ from impact parameter $b$ and sky-projected obliquity $\\vartheta$. In 800 synthetic single-transit light curves spanning $f = 0.03$\\textendash$0.30$, $\\vartheta = 0^\\circ$\\textendash$90^\\circ$, and noise levels $\\sigma_w = 1$\\textendash$256$ ppm, a $3\\sigma$ detection of $f$ is achieved in 559 of 700 cases below 256 ppm; Saturn-like oblateness ($f \\approx 0.09$) is detected in about 59% of configurations, with retrieved $f$ and $\\vartheta$ in $1\\sigma$ agreement with the truth in roughly two-thirds of those detections. No tested configuration yields a reliable detection at $\\sigma_w = 256$ ppm, and values below $f \\approx 0.06$ remain below the method's reach.","pith_inferences":["By extension, co-adding several transits of the same planet should reduce the effective white noise by roughly $\\sqrt{N}$, so targets that fail at 256 ppm in a single transit may become detectable in multi-transit campaigns; the paper's single-transit setting is conservative.","An untested extension is to apply the same density-anchor retrieval to eccentric orbits; since the paper's prior is only strictly valid for circular orbits, such systems would need independent eccentricity and periastron constraints from radial velocities.","The measured sky-projected obliquity $\\vartheta$, combined with independent spin-orbit measurements, could begin to constrain the three-dimensional spin geometry of hot planets; the paper does not pursue that combination.","Existing space-photometry archives of bright stars that already meet the sub-256 ppm noise requirement could be re-analyzed for oblateness without new observations."],"forward_implications":["A single transit of a Saturn-like planet ($f\\approx0.09$) around a bright, asteroseismically characterized star can deliver a $3\\sigma$ oblateness measurement in about 59% of tested geometries, making dedicated follow-up observations of individual transits worthwhile.","Oblateness below about $f=0.06$, the value of Jupiter, is not recoverable at any tested noise level, so single-transit non-detections cannot yet constrain the flattening of most known planets.","At $\\sigma_w=256$ ppm per 60 s exposure, only 54 of 100 configurations reach $3\\sigma$ and the recovered values are frequently inaccurate, establishing a noise floor for planning transit-oblateness programs.","The method relies on prior knowledge of the host star's mean density to about 1% and of limb-darkening coefficients to 0.01, so the target list is restricted to bright, quiet, asteroseismically characterized stars.","Because the model is about 25% faster than the analytical spherical model at full precision, large retrieval grids over many planets are computationally feasible."],"supporting_citations":[{"why":"Supplies the analytical spherical-transit light curve used as the accuracy baseline for the numerical model in the non-oblate limit.","marker":"K. Mandel & E. Agol (2002)"},{"why":"Identified the b-f-vartheta degeneracy that the paper confirms and counters with its density prior.","marker":"S. Dholakia et al. (2024)"},{"why":"Provides the asteroseismic stellar-density precision (0.5-2.6%) that justifies the 1% density assumption behind the a/R-star prior.","marker":"V. Silva Aguirre et al. (2017)"},{"why":"Describes the transit-modelling code into which the new quadrature-based oblate model is implemented and which performs the fits.","marker":"Sz. Csizmadia (2020)"},{"why":"Supplies the wavelet-based time-correlated noise model used in the retrieval tests.","marker":"J. A. Carter & J. N. Winn (2009)"},{"why":"Provides an independent analytical oblate-transit code used to cross-check the new model's light curves.","marker":"B. Cassese et al. (2024)"},{"why":"Establishes the transit-duration relation connecting b and a/R-star, the degeneracy the density prior is designed to break.","marker":"S. Seager & G. Mallén-Ornelas (2003)"},{"why":"Supplies the differential-evolution MCMC sampler used for parameter retrieval in the 800 injection tests.","marker":"E. B. Ford (2006)"}],"fun_headline_variants":["Single transit reveals Saturn-like oblateness when star's density is known","Exoplanet flatness from one transit, gated by precise stellar density","Rotational oblateness detectable in single transits via stellar density prior","One transit, one oblateness: stellar density breaks the degeneracy"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The retrieval depends on an externally known stellar mean density (about 1% precision) and on a strictly circular orbit, since the density-based prior on $a/R_\\star$ is what separates oblateness from the impact-parameter degeneracy and the prior is only strictly valid for circular orbits.","fun_headline_variants_meta":{"raw":{"variants":["Single transit reveals Saturn-like oblateness when star's density is known","Exoplanet flatness from one transit, gated by precise stellar density","Rotational oblateness detectable in single transits via stellar density prior","One transit, one oblateness: stellar density breaks the degeneracy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000569,"raw_usage":{"total_tokens":2780,"prompt_tokens":1120,"completion_tokens":1660,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":736,"completion_tokens_details":{"reasoning_tokens":1583}},"tokens_in":736,"tokens_out":1660,"duration_ms":14245,"temperature":1.0,"reasoning_tokens":1583,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T15:34:24.564768+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the same 800-case injection grid with eccentric orbits (for example $e=0.05$, all other parameters and priors unchanged): if the $3\\sigma$ detection rate for Saturn-like oblateness falls substantially below 59%, or the recovered $f$ is biased by more than 0.02, the circular-orbit assumption is the part of the claim that fails.","supporting_citations":[],"review_version":1}