{"id":"1a43b726-2c16-40fb-bf71-1a565ac0af11","arxiv_id":"2501.18700","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"New transit timing observations reduce LP 791-18 d's mass uncertainty by over a factor of two and, under a tidally damped prior, make its eccentricity consistent with zero.","lead":"Three new Palomar/WIRC transits of the sub-Neptune LP 791-18 c cut the mass uncertainty of the Earth-sized LP 791-18 d to 0.91 ± 0.19 Earth masses and show that previously reported orbital eccentricities were partly an artifact of a coordinate parameterization. The revised orbit sharpens predictions for JWST secondary eclipse timing, which can test whether tidal heating drives volcanism on this temperate rocky planet.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The headline factor-of-two mass-precision gain is conditional on the tidally damped eccentricity prior; the data do not prefer that prior, and the undamped fit yields only a ~1.5x gain.","rationale":"The reader's weakest_assumption correctly identifies the tidally damped free-eccentricity prior as the load-bearing assumption, but focuses its impact on the eccentricity and tidal-heating conclusions. I argue the same assumption also undercuts the paper's headline mass-precision claim: the abstract's 'reduce the mass uncertainty by more than a factor of two' is only true for the damped case, and the data do not actually prefer that case (ΔBIC = −1.3). The undamped fit, which the paper itself presents as equally viable, yields M_d = 0.98+0.29−0.25, only a ~1.5x improvement over P23. The analysis itself is sound and transparent, so the paper should be accepted, but the abstract should be revised to avoid overstating the robustness of the mass improvement. This is a minor framing issue rather than a correctness flaw, hence CONDITIONAL rather than REJECT.","tokens_in":19036,"tokens_out":13062,"duration_ms":120093,"concrete_test":"Run a Bayesian model comparison (e.g., nested sampling for evidence or leave-one-out cross-validation) between the damped and undamped eccentricity models on the full TTV data set. If the undamped model is not disfavored by a Bayes factor >3, the mass headline should be changed to the undamped value (0.98+0.29−0.25 M⊕), and the abstract should explicitly state that the >2x mass-precision improvement is conditional on the damped-eccentricity assumption.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The abstract's headline result, 'We reduce the mass uncertainty by more than a factor of two (M_d = 0.91 ± 0.19 M⊕)', is obtained under the tidally damped free-eccentricity prior described in §3.2, which imposes Gaussian priors centered at zero with σ_e,d = 0.001 and σ_e,c = 0.0001. This prior is physically motivated but not validated by the data: §3.3 reports that the undamped (uniform eccentricity prior) fit gives M_d = 0.98+0.29−0.25 M⊕ and e_d = 0.056+0.015−0.014, and the BIC comparison in §3.3 yields ΔBIC = −1.3, i.e., the data actually slightly favor the undamped model. Compared to P23's M_d = 0.9+0.5−0.4 M⊕, the damped case improves the 68% interval width by ~2.4x, but the undamped case improves it by only ~1.5x. Thus the central quantitative claim is prior-dominated, not a robust data-driven constraint. The paper is transparent about this in the body, so the weakness is in the framing rather than the analysis.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents three new Palomar/WIRC diffuser-assisted transit observations of LP 791-18 c, with timing precisions of 6, 4, and 3 seconds, and combines them with the published P23 transit-timing dataset. The authors fit the combined TTV data with TTVFast, first reproducing P23's results, then switching from the ecos(ω)/esin(ω) parameterization to a √e cos(ω)/√e sin(ω) parameterization with and without a tidally damped free-eccentricity prior. Under the damped prior they find M_d = 0.91 ± 0.19 M⊕ and eccentricities consistent with zero (e_d = 0.0011^{+0.0010}_{-0.0008}); under a uniform eccentricity prior they find M_d = 0.98^{+0.29}_{-0.25} M⊕ and e_d = 0.056^{+0.015}_{-0.014}. They use these results to revise tidal heating estimates and to predict JWST secondary-eclipse timing offsets that can discriminate between the damped and undamped states.","tokens_in":19363,"tokens_out":5302,"duration_ms":50027,"significance":"If the results hold, the paper delivers a significantly improved mass constraint for a temperate Earth-sized planet, identifies a bias in the previously reported P23 eccentricities, and produces a concrete, testable JWST prediction that separates two very different tidal-heating scenarios. The analysis is generally careful: the authors validate their TTV model by reproducing P23's solution, use a well-tested n-body code, check that an e=0 model fits within the timing precision, and are transparent about the two eccentricity priors. However, the headline mass-precision improvement is conditional on a prior that the data do not actually prefer, so the significance statement must be framed with that caveat. The extraneous inserted text in the manuscript is a separate but serious issue.","major_comments":[{"comment":"The headline claim 'We reduce the mass uncertainty by more than a factor of two (M_d = 0.91 ± 0.19 M⊕)' is obtained under the tidally damped free-eccentricity prior described in §3.2. The uniform-eccentricity fit in §3.3 yields M_d = 0.98^{+0.29}_{-0.25} M⊕, only a ~1.5× improvement over P23's width, and the BIC comparison there gives ΔBIC = −1.3, i.e., the data slightly favor the undamped model. As stated, the abstract and summary overstate what is robustly measured. Please qualify the claim (e.g., 'under the tidally damped-state prior') and report both mass constraints prominently.","section":"Abstract and §5"},{"comment":"The manuscript contains a large extraneous passage beginning 'Volcanism on Exoplanets 9' that includes equations numbered (13) and (14), a discussion of HST observations of L 98-59, and parameter values (M1 = 2.2 M⊕, P2 = 7.45 d) that have nothing to do with LP 791-18. This appears to be an accidentally included fragment from another paper (likely Seligman et al. 2023). It must be removed; as submitted, the manuscript is not a coherent paper.","section":"Full text, between §5 and Acknowledgments"},{"comment":"The statement that 'the originally reported values from P23 were an artifact of the chosen e and ω parameterization' is only established for fits that also impose the tidally damped free-eccentricity prior. Under the uniform eccentricity prior, the updated fit still returns e_d = 0.056^{+0.015}_{-0.014} and e_c = 0.062^{+0.017}_{-0.014} (Table 2, §3.3), so the non-zero eccentricities are not purely a parameterization artifact. Please restrict the 'artifact' claim to the damped case, as §5 does.","section":"§3.1"}],"minor_comments":[{"comment":"Introductory paragraph contains a duplicated word: 'with amplitudes of ~2.5 min and ~0.5 min for for planets c and d, respectively'.","section":"§1"},{"comment":"The sentence 'As originally noted by Ford (2006), fitting for esin(ω) and esin(ω) results...' should read 'fitting for ecos(ω) and esin(ω)'.","section":"§3.1"},{"comment":"The damped-fit period for LP 791-18 c is printed as '4.9899100+0.0000012−0.0000014'; this should be formatted with a reasonable number of decimal places (e.g., 4.98991 ± 0.0000013).","section":"Table 2"},{"comment":"The caption says 'using 104 draws from the period, e, and ω distributions'; this should be 10^4 draws to be unambiguous.","section":"Figure 6 caption"},{"comment":"The software list credits 'Claude 3.5 Sonnet' without specifying its role; if an AI assistant was used for code or text, please state its contribution explicitly, and if not, remove it.","section":"§6 Acknowledgments"}],"recommendation":"major_revision","confidential_remarks":"The scientific content is largely sound and the TTV analysis appears careful, but the submitted text contains a clearly extraneous inserted passage from another work, which suggests a compilation error in the manuscript source. The authors should be asked to verify the version they submit. The mass-precision headline should also be reframed so that the prior-dependence is transparent; this is fixable in revision and does not undermine the core TTV analysis."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this if you care about small-planet TTVs or the LP 791-18 system. The paper does two genuinely useful things. First, it adds three Palomar/WIRC transits of planet c with 3–6 s timing precision, a real improvement over the Spitzer data, and uses them to cut the mass uncertainty of LP 791-18 d from about 0.9 M⊕ wide to 0.38 M⊕ (damped case) or 0.54 M⊕ (undamped). Second, it shows that the small eccentricities reported in Peterson et al. (2023) were an artifact of the ecosω/esinω parameterization; switching to √e cosω/√e sinω makes those eccentricities vanish under the damped prior. That is a clean, reproducible point and the paper explains it well.\n\nThe analysis is careful: they reproduce P23's solution before extending it, use TTVFast, test an e=0 model, and report both damped and undamped fits. The mass shift between the two priors is within 1σ, so the headline mass is not held hostage by the prior, only the size of the quoted uncertainty is.\n\nThe soft spot is framing. The abstract's 'reduce the mass uncertainty by more than a factor of two' refers to the tidally damped fit. Without that prior the improvement is closer to 1.5–1.7x, and the data do not actually prefer the damped model (ΔBIC = −1.3, slightly favoring undamped). The paper is transparent in the body and explicitly says current data can't distinguish the two states, so this is not a concealed flaw. But a reader skimming the abstract could easily come away thinking the factor-of-two gain is prior-independent. A revised abstract should either quote both numbers or qualify the headline.\n\nThe eccentricity and tidal heating conclusions are genuinely prior-dependent: e_d = 0.0011 vs 0.056 changes the tidal heat flux by orders of magnitude. The paper handles this honestly by laying out both scenarios and predicting JWST eclipse timing offsets that can discriminate. That is the right way to present it.\n\nWho should read this: anyone working on TTV retrieval, small-planet masses, or tidal heating of rocky planets. The parameterization bias lesson is worth knowing generally. It deserves peer review — the data are new, the analysis is competent, and the claims are clearly scoped. I would suggest the referee push for the abstract fix and for a table column that makes the prior dependence of the mass uncertainty visible at a glance.","headline":"Solid TTV follow-up that halves the mass uncertainty of LP 791-18 d and exposes a parameterization bias in the prior mass/eccentricity claims, though the headline precision gain leans on a physically motivated but data-indifferent eccentricity prior.","tokens_in":19904,"tokens_out":3438,"would_cite":true,"duration_ms":30046,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"New transit timings halve the mass uncertainty of the Earth-sized exoplanet LP 791-18 d and show that its reported eccentricity, the basis for volcanic-heating predictions, was a fitting artifact — a question JWST secondary-eclipse timing…","keywords":["LP 791-18 d","transit timing variations","tidal heating","eccentricity","secondary eclipse timing","M dwarf exoplanets","Earth-sized planets","JWST"],"falsifier":"Measure the secondary eclipse of LP 791-18 d with JWST: the damped model predicts a timing offset of $-0.2^{+2.0}_{-2.7}$ minutes relative to a circular orbit, while the undamped model predicts $-117^{+41}_{-47}$ minutes, so an observed offset with magnitude greater than roughly 10–15 minutes would rule out the tidally damped state and an offset consistent with zero would confirm it.","tokens_in":18851,"feed_emoji":"🌋","tokens_out":18109,"duration_ms":134286,"temperature":0.7,"pith_summary":"LP 791-18 d is one of only a handful of temperate, Earth-sized planets known outside the solar system, and earlier work suggested its slightly eccentric orbit would tidally heat its interior enough to drive volcanic outgassing. This paper adds three high-precision transit timings of its outer companion, LP 791-18 c, and re-fits the whole timing data set. The updated mass is $M_d = 0.91 \\pm 0.19\\,M_\\oplus$, reducing the uncertainty by more than a factor of two. The authors also find that the small non-zero eccentricities reported by the earlier study were an artifact of how the fit parameterized the orbit: with an unbiased parameterization and a physically motivated tidally-damped prior, the orbits are consistent with circular. The current data cannot distinguish this damped state from a higher-eccentricity state, but the paper shows that a single JWST secondary eclipse timing measurement can settle the question.","feed_headline":"Transit timings halve the mass error of an Earth-like exoplanet","feed_subtitle":"The catch: its eccentricity and tidal heating hinge on a prior the data can't yet test.","key_machinery":"The argument rests on three linked pieces of machinery. The first is the TTVFast n-body model, which converts the planet-star mass ratios and osculating orbital elements into transit times; the new data extend the baseline by about three years and reach sub-6-second precision, making them the most precise timing measurements of this system. The second is the eccentricity parameterization: fitting $e\\cos\\omega$ and $e\\sin\\omega$ imposes an effective linear prior on $e$ whose correction slightly overweights low eccentricity, producing spurious small non-zero values, whereas the $\\sqrt{e}\\cos\\omega$ and $\\sqrt{e}\\sin\\omega$ parameterization recovers a uniform prior. The third is the tidally damped free-eccentricity prior, constructed by running short REBOUND simulations at each MCMC step to split the osculating eccentricity into forced and free parts and then penalizing free eccentricities away from zero. The paper's decisive output is the predicted secondary-eclipse timing offset $\\Delta t = 2P e\\cos\\omega/\\pi$, which translates the two eccentricity solutions into observationally distinguishable predictions.","core_discovery":"The paper's central claim is that the three new Palomar/WIRC transit observations of LP 791-18 c, with timing precisions of 3–6 seconds, tighten the dynamical constraints on the innermost Earth-sized planet and expose a bias in the earlier analysis. Re-fitting the complete set of transit times with the TTVFast n-body model, the authors obtain $M_d = 0.91 \\pm 0.19\\,M_\\oplus$ under the tidally damped prior and $M_d = 0.98^{+0.29}_{-0.25}\\,M_\\oplus$ without it. When the eccentricity is parameterized as $\\sqrt{e}\\cos\\omega$ and $\\sqrt{e}\\sin\\omega$ rather than $e\\cos\\omega$ and $e\\sin\\omega$, the previously reported values $e_d \\approx 0.0015$ and $e_c \\approx 0.0008$ collapse to values consistent with zero, leading the authors to conclude those earlier detections were artifacts of the chosen parameterization. Absent the damped prior, the data instead prefer $e_d = 0.056^{+0.015}_{-0.014}$ and $e_c = 0.062^{+0.017}_{-0.014}$. The two models fit the current timing data almost equally well ($\\Delta$BIC $= -1.3$), so the paper's main forward prediction is the secondary eclipse timing offset: essentially zero minutes for the damped state versus roughly $-117$ minutes for the undamped state. An observed offset larger than about 10–15 minutes would confirm that tidal damping has not yet circularized the orbit and would make LP 791-18 d a much stronger volcanic candidate.","pith_inferences":["The same parameterization bias identified here may affect other TTV studies that report small non-zero eccentricities near the detection threshold; re-running those fits with $\\sqrt{e}$-type parameters could revise eccentricity estimates without any new data.","If JWST returns an intermediate eclipse offset, it would map $e\\cos\\omega$ continuously and, combined with the TTV mass constraint, could turn LP 791-18 d into a calibrated probe of terrestrial tidal dissipation across $Q'$ values.","The age of the system is a hidden lever: if LP 791-18 is older than roughly 1 Gyr, retaining a free eccentricity would require even larger $Q'$ values (weaker dissipation) or a recent dynamical disturbance such as a stellar flyby, a scenario the paper mentions only in passing.","Should the eclipse timing confirm the damped state, the TTVs' inability to sense $e \\lesssim 0.01$ means the forced eccentricity of about 0.0015 remains the relevant value for tidal heating, leaving the volcanic-outgassing question to atmospheric observations rather than dynamics."],"forward_implications":["If the tidally damped solution is correct, LP 791-18 d's orbit is effectively circular and its tidal heat flux, while small, is still large enough to matter for silicate melting and planetary evolution, consistent with the earlier Io-like interpretation.","If the undamped solution is correct, the tidal heat flux is comparable to the insolation flux, the interior temperature rises by roughly 145 K, and LP 791-18 d becomes a much more compelling target for volcanic outgassing, similar to L 98-59 d.","JWST secondary eclipse timing can decide between the two: a measured eclipse offset of more than about 10–15 minutes rules out the tidally damped state and constrains the tidal quality factor $Q'$.","The improved mass constraint, $M_d = 0.91 \\pm 0.19\\,M_\\oplus$, places LP 791-18 d on an Earth-like rock-iron composition track and confirms it has at most a modest hydrogen envelope of roughly 2–2.5% by mass."],"supporting_citations":[{"why":"Supplies the original discovery of LP 791-18 d, the full TTV data set, and the REBOUND-derived damped-state free eccentricity prior that this paper re-tests.","marker":"P23"},{"why":"Provides the TTVFast symplectic n-body code used to model and fit all transit times.","marker":"Deck et al. (2014)"},{"why":"Showed that fitting $e\\cos\\omega$ and $e\\sin\\omega$ imposes an effective linear prior on $e$, the bias source claimed for P23's small eccentricities.","marker":"Ford (2006)"},{"why":"Recommended the $\\sqrt{e}\\cos\\omega$ and $\\sqrt{e}\\sin\\omega$ parameterization that yields a uniform prior and removes the near-zero bias.","marker":"Eastman et al. (2013)"},{"why":"Documents the mass-eccentricity degeneracy and chopping TTV behavior that explains why the undamped fit raises both planet masses.","marker":"Lithwick et al. (2012)"},{"why":"Gives the eclipse timing offset formula $\\Delta t = 2P e\\cos\\omega/\\pi$ used to predict JWST secondary eclipse times.","marker":"Deming et al. (2005)"},{"why":"Provides the tidal circularization timescale and $Q'$ thresholds that define the damped versus undamped regimes.","marker":"Puranam & Batygin (2018)"},{"why":"Supplies the REBOUND N-body package used to separate forced from free eccentricity for the damped prior.","marker":"Rein & Tamayo (2015)"},{"why":"Describes the diffuser-assisted WIRC photometry technique that yielded the new 3–6 second transit timing measurements.","marker":"Vissapragada et al. (2020)"}],"fun_headline_variants":["Transit timings halve mass error for an Earth-sized planet","Halved mass error, but eccentricity is still a toss-up","New TTVs pin down mass, leave tidal heating ambiguous","Tighter mass for small exoplanet, JWST to test eccentricity","Mass uncertainty cut in half for LP 791-18 d, but heat source unclear"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that tidal damping has already reduced the planets' free eccentricities to near zero, so the observed eccentricity equals only the small forced value; if tidal dissipation is actually weak ($Q' \\gtrsim 10^4$) or the system was recently disturbed, the data equally allow $e_d \\approx 0.056$ and a tidal heat flux orders of magnitude larger.","fun_headline_variants_meta":{"raw":{"variants":["Transit timings halve mass error for an Earth-sized planet","Halved mass error, but eccentricity is still a toss-up","New TTVs pin down mass, leave tidal heating ambiguous","Tighter mass for small exoplanet, JWST to test eccentricity","Mass uncertainty cut in half for LP 791-18 d, but heat source unclear"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000994,"raw_usage":{"total_tokens":4398,"prompt_tokens":1318,"completion_tokens":3080,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":934,"completion_tokens_details":{"reasoning_tokens":2984}},"tokens_in":934,"tokens_out":3080,"duration_ms":22009,"temperature":1.0,"reasoning_tokens":2984,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T22:47:12.252111+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the secondary eclipse of LP 791-18 d with JWST: the damped model predicts a timing offset of $-0.2^{+2.0}_{-2.7}$ minutes relative to a circular orbit, while the undamped model predicts $-117^{+41}_{-47}$ minutes, so an observed offset with magnitude greater than roughly 10–15 minutes would rule out the tidally damped state and an offset consistent with zero would confirm it.","supporting_citations":[],"review_version":1}