REVIEW 3 major objections 5 minor 46 references
Prospects of detecting rotational flatness of exoplanets from space-based photometry
T0 review · 3 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read 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.
desk verdict 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. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (3)
- [Section 4.2 and Table 7] 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 3.3] 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.
- [Abstract and Section 4.2] 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.
minor comments (5)
- [Throughout] 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.
- [Abstract, Section 3.2, Section 5] 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.
- [Eq. (26) and Fig. 5] 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.
- [Table 7] 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 4.2] 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.
Circularity Check
No circularity: the central claim is an injection-retrieval sensitivity study with external benchmarks (Mandel-Agol and squishyplanet); the text/Table-7 prior discrepancy is a correctness issue, not a circular reduction.
full rationale
The paper's central derivation chain is a Monte Carlo power study: it injects transit light curves generated with a numerical quadrature model at known oblateness values, then retrieves parameters with the same model and reports 3-sigma detection rates. This is a self-consistency test, not a derivation in which the output is equivalent to the input by construction. The forward model is independently benchmarked in the spherical limit against the analytic Mandel-Agol model and in the oblate case against squishyplanet, so the validation is not purely self-referential. The detection claim is explicitly conditional on stated priors: a 1% stellar-density prior on a/R* (Eq. 30) and 1% limb-darkening priors. These priors are adopted from external asteroseismic precision estimates (Silva Aguirre et al. 2017), not fitted from the target data, so they are not a fitted input renamed as a prediction. The manuscript itself flags the circular-orbit restriction of the a/R* prior trick, which is a limitation but not circularity. The internal inconsistency between the text value sigma(a/R*) = 0.24 and Table 7's N(73.26, 0.04) affects the realism or reproducibility of the assumed prior and the resulting detection rates, but it does not make the derivation circular. Self-citations to TLCM and Kálmán et al. are software and noise-model references, not load-bearing uniqueness arguments. No equation defines the predicted quantity in terms of the input, and no parameter is fitted to a subset and then reported as a prediction. Therefore no significant circularity is present.
Assumptions & free parameters
free parameters (3)
- a/R* Gaussian prior width =
0.24 (text) / 0.04 (Table 7)
- Limb darkening Gaussian prior widths =
0.01
- Red noise amplitude scaling =
1/10 of ARIMA clone from Kalman et al. 2024
assumptions (5)
- standard math Gauss-Legendre quadrature with sum of weights = 2 (Eq. 8)
- domain assumption Planet shape is a biaxial ellipsoid (spheroid) with rigid-body rotation; potential V = -GMp/rPol = -GMp/rEq - (1/2) omega^2 rEq^2 (Eq. 3)
- domain assumption Star is spherical with static limb darkening; the 1.1 factor in the delta^2 < (1.1(1+rEq))^2 transit-phase cut is arbitrary (Sect. 2.2)
- domain assumption ARIMA clone of HST photometry, scaled by 1/10, represents space-based telescope red noise (Sect. 4.1)
- domain assumption Stellar density known to 1% via asteroseismology, implying Delta(a/R*)/(a/R*) = 0.0033 (Eq. 30)
Cite this review
Pith. "Pith review of Prospects of detecting rotational flatness of exoplanets from space-based photometry." pith.science (2026). https://pith.science/paper/HBR4MI3K
@misc{pith2026250715359,
author = {Pith},
title = {Pith review of: Prospects of detecting rotational flatness of exoplanets from space-based photometry},
year = {2026},
howpublished = {\url{https://pith.science/paper/HBR4MI3K}},
note = {Machine review of arXiv:2507.15359}
}
abstract
In the era of photometry with space-based telescopes, such as CHEOPS (CHaracterizing ExOPlanets Satellite), JWST (James Webb Space Telescope), PLATO (PLAnetary Transits and Oscillations of stars), and ARIEL (Atmospheric Remote-sensing Infrared Exoplanet Large-survey), the road has opened for detecting subtle distortions in exoplanet transit light curves -- resulting from their non-spherical shape. We investigate the prospects of retrieval of rotational flatness (oblateness) of exoplanets at various noise levels. We present a novel method for calculating the transit light curves based on the Gauss-Legendre quadrature. We compare it in the non-rotating limit to the available analytical models. We conduct injection-and-retrieval tests to assess the precision and accuracy of the retrievable oblateness values. We find that the light curve calculation technique is about $25$\% faster than a well-known analytical counterpart, while still being precise enough. We show that a $3 \sigma$ oblateness detection is possible for a planet orbiting bright enough stars, by exploiting a precise estimate on the stellar density obtained e.g. from asteroseismology. We also show that for noise levels $\geq 256$ ppm (expressed as point-to-point scatter with a $60$~s exposure time) detection of planetary oblateness is not reliable.
Figures
Figures from the paper (14 more)
Reference graph
Works this paper leans on
-
[1]
Akinsanmi, B., Barros, S. C. C., Santos, N. C., et al. 2019, A&A, 621, A117, doi: 10.1051/0004-6361/201834215
-
[2]
Akinsanmi, B., Lendl, M., Bou´ e, G., & Barros, S. C. C. 2024a, A&A, 682, A15, doi: 10.1051/0004-6361/202347739 22
-
[3]
Akinsanmi, B., Barros, S. C. C., Lendl, M., et al. 2024b, A&A, 685, A63, doi: 10.1051/0004-6361/202348502
-
[4]
2009, A&A, 506, 411, doi: 10.1051/0004-6361/200810860
Auvergne, M., Bodin, P., Boisnard, L., et al. 2009, A&A, 506, 411, doi: 10.1051/0004-6361/200810860
-
[5]
Barnes, J. W. 2009, ApJ, 705, 683, doi: 10.1088/0004-637X/705/1/683
-
[6]
Barnes, J. W., & Fortney, J. J. 2003, ApJ, 588, 545, doi: 10.1086/373893
doi:10.1086/373893 2003
-
[7]
2021, Experimental Astronomy, 51, 109, doi: 10.1007/s10686-020-09679-4
Benz, W., Broeg, C., Fortier, A., et al. 2021, Experimental Astronomy, 51, 109, doi: 10.1007/s10686-020-09679-4
-
[8]
2022, ApJ, 935, 178, doi: 10.3847/1538-4357/ac82b2 Bernab` o, L
Berardo, D., & de Wit, J. 2022, ApJ, 935, 178, doi: 10.3847/1538-4357/ac82b2 Bernab` o, L. M., Csizmadia, S., Smith, A. M. S., et al. 2024, A&A, 684, A78, doi: 10.1051/0004-6361/202346852 Bernab` o, L. M., Csizmadia, S., Smith, A. M. S., et al. 2025, A&A, 694, A233, doi: 10.1051/0004-6361/202451994
Show all 46 references
-
[9]
2024, AJ, 168, 243, doi: 10.3847/1538-3881/ad7d8d
Bhowmick, U., & Khaire, V. 2024, AJ, 168, 243, doi: 10.3847/1538-3881/ad7d8d
2024 doi
-
[10]
J., Koch, D., Basri, G., et al
Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977, doi: 10.1126/science.1185402
2010 doi
-
[11]
A., & Winn, J
Carter, J. A., & Winn, J. N. 2009, ApJ, 704, 51, doi: 10.1088/0004-637X/704/1/51
2009 doi
-
[12]
2024, The Journal of Open Source Software, 9, 6972, doi: 10.21105/joss.06972
Cassese, B., Vega, J., Lu, T., et al. 2024, The Journal of Open Source Software, 9, 6972, doi: 10.21105/joss.06972
2024 doi
-
[13]
2004, A&A, 428, 1001, doi: 10.1051/0004-6361:20041673
Claret, A. 2004, A&A, 428, 1001, doi: 10.1051/0004-6361:20041673
2004 doi
-
[14]
2020, MNRAS, 496, 4442, doi: 10.1093/mnras/staa349
Csizmadia, Sz. 2020, MNRAS, 496, 4442, doi: 10.1093/mnras/staa349
2020 doi
-
[15]
Csizmadia, Sz., Hellard, H., & Smith, A. M. S. 2019, A&A, 623, A45, doi: 10.1051/0004-6361/201834376
2019 doi
-
[16]
Csizmadia, Sz., Smith, A. M. S., K´ alm´ an, S., et al. 2023, A&A, 675, A106, doi: 10.1051/0004-6361/202141302
2023 doi
-
[17]
2011, A&A, 531, A41, doi: 10.1051/0004-6361/201117009
Csizmadia, Sz., Moutou, C., Deleuil, M., et al. 2011, A&A, 531, A41, doi: 10.1051/0004-6361/201117009
2011 doi
-
[18]
I., & Johnson, J
Dawson, R. I., & Johnson, J. A. 2012, ApJ, 756, 122, doi: 10.1088/0004-637X/756/2/122
2012 doi
-
[19]
2018, A&A, 619, A97, doi: 10.1051/0004-6361/201731068
Deleuil, M., Aigrain, S., Moutou, C., et al. 2018, A&A, 619, A97, doi: 10.1051/0004-6361/201731068
2018 doi
- [20]
-
[21]
1992, A&A, 259, 227
Diaz-Cordoves, J., & Gimenez, A. 1992, A&A, 259, 227
1992
-
[22]
Ford, E. B. 2006, ApJ, 642, 505, doi: 10.1086/500802
2006 doi
-
[23]
P., Mather, J
Gardner, J. P., Mather, J. C., Clampin, M., et al. 2006, SSRv, 123, 485, doi: 10.1007/s11214-006-8315-7
2006 doi
-
[24]
2019, ApJ, 878, 119, doi: 10.3847/1538-4357/ab2048
Hellard, H., Csizmadia, S., Padovan, S., et al. 2019, ApJ, 878, 119, doi: 10.3847/1538-4357/ab2048
2019 doi
-
[25]
2020, ApJ, 889, 66, doi: 10.3847/1538-4357/ab616e
Hellard, H., Csizmadia, S., Padovan, S., Sohl, F., & Rauer, H. 2020, ApJ, 889, 66, doi: 10.3847/1538-4357/ab616e
2020 doi
-
[26]
1997, A&A, 327, 199 K´ alm´ an, Sz., Csizmadia, S., Simon, A
Hestroffer, D. 1997, A&A, 327, 199 K´ alm´ an, Sz., Csizmadia, S., Simon, A. E., et al. 2024, MNRAS, 528, L66, doi: 10.1093/mnrasl/slad169 K´ alm´ an, Sz., Szab´ o, G. M., & Csizmadia, S. 2023, A&A, 675, A107, doi: 10.1051/0004-6361/202143017
1997 doi
-
[27]
A., & Sobieski, S
Klinglesmith, D. A., & Sobieski, S. 1970, AJ, 75, 175, doi: 10.1086/110960
1970 doi
- [28]
-
[29]
2002, ApJL, 580, L171, doi: 10.1086/345520
Mandel, K., & Agol, E. 2002, ApJL, 580, L171, doi: 10.1086/345520
2002 doi
-
[30]
2023, A&A, 677, A133, doi: 10.1051/0004-6361/202245287
Matuszewski, F., Nettelmann, N., Cabrera, J., B¨ orner, A., & Rauer, H. 2023, A&A, 677, A133, doi: 10.1051/0004-6361/202245287
2023 doi
-
[31]
W., Feinberg, L
McElwain, M. W., Feinberg, L. D., Perrin, M. D., et al. 2023, PASP, 135, 058001, doi: 10.1088/1538-3873/acada0
2023 doi
-
[32]
Milne, E. A. 1921, MNRAS, 81, 361, doi: 10.1093/mnras/81.5.361
1921 doi
-
[33]
2025, ApJL, 981, L7, doi: 10.3847/2041-8213/adb42b
Vanderburg, A. 2025, ApJL, 981, L7, doi: 10.3847/2041-8213/adb42b
2025 doi
-
[34]
2014, Experimental Astronomy, 38, 249, doi: 10.1007/s10686-014-9383-4
Rauer, H., Catala, C., Aerts, C., et al. 2014, Experimental Astronomy, 38, 249, doi: 10.1007/s10686-014-9383-4
2014 doi
- [35]
-
[36]
2019, MNRAS, 490, 1111, doi: 10.1093/mnras/stz2556
Rein, E., & Ofir, A. 2019, MNRAS, 490, 1111, doi: 10.1093/mnras/stz2556
2019 doi
-
[37]
R., Winn, J
Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, Journal of Astronomical Telescopes, Instruments, and Systems, 1, 014003, doi: 10.1117/1.JATIS.1.1.014003
2015 doi
-
[38]
Russell, H. N. 1912, ApJ, 36, 54, doi: 10.1086/141952
1912 doi
-
[39]
2003, ApJ, 585, 1038, doi: 10.1086/346105
Seager, S., & Mall´ en-Ornelas, G. 2003, ApJ, 585, 1038, doi: 10.1086/346105
2003 doi
-
[40]
S., & Wilk, M
Shapiro, S. S., & Wilk, M. B. 1965, Biometrika, 52, 591, doi: 10.1093/biomet/52.3-4.591 Silva Aguirre, V., Lund, M. N., Antia, H. M., et al. 2017, ApJ, 835, 173, doi: 10.3847/1538-4357/835/2/173
1965 doi
-
[41]
Sing, D. K. 2010, A&A, 510, A21, doi: 10.1051/0004-6361/200913675 Szab´ o, Gy. M., K´ alm´ an, S., Pribulla, T., et al. 2022, Experimental Astronomy, 53, 607, doi: 10.1007/s10686-021-09777-x
2010 doi
-
[42]
2022, in European Planetary Science Congress, EPSC2022–1114, doi: 10.5194/epsc2022-1114
Tinetti, G., Eccleston, P., Lueftinger, T., et al. 2022, in European Planetary Science Congress, EPSC2022–1114, doi: 10.5194/epsc2022-1114
2022 doi
-
[43]
2018, Experimental Astronomy, 46, 135, doi: 10.1007/s10686-018-9598-x
Tinetti, G., Drossart, P., Eccleston, P., et al. 2018, Experimental Astronomy, 46, 135, doi: 10.1007/s10686-018-9598-x
2018 doi
-
[44]
2006,, Tech
Trefethen, L. 2006,, Tech. rep
2006
-
[45]
A., & Rucinski, S
Wade, R. A., & Rucinski, S. M. 1985, A&AS, 60, 471 23
1985
-
[46]
E., & Devinney, E
Wilson, R. E., & Devinney, E. J. 1971, ApJ, 166, 605, doi: 10.1086/150986
1971 doi
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.