{"id":"ab43294a-ac95-4dee-9e76-255da3751bd1","arxiv_id":"2505.10304","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A Bayesian retrieval combining Love numbers and atmospheric metallicity shows hot Jupiter interiors can be constrained if Love number precision is better than 40% (homogeneous) or 15% (dilute core), and finds evidence for a core in WASP-19Ab.","lead":"This paper shows that precisely measuring a hot Jupiter's tidal response (Love number) together with its atmospheric metal content can pin down whether the planet has a core and how much metal it holds. It finds that only one of five planets with measured Love numbers, WASP-19Ab, has data good enough to test this, and that planet appears to have a core of about 80 Earth masses.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The WASP-19Ab core detection depends on the linear-tidal k2 interpretation; the ~19% non-linear bias flagged in §5.1 exceeds the 12% measurement error, so the core claim should be tested against a k2 shifted by that systematic before being treated as settled.","rationale":"The paper is careful and the retrieval machinery is validated on a test planet; the homogeneous/inhomogeneous comparison and the Appendix B prior check show the WASP-19Ab core mass estimate is stable to the luminosity/heat prior. The remaining threat to the central claim is not internal inconsistency but the domain of validity of Eq. (5): it solves the linear, non-rotating, hydrostatic tidal response. The authors themselves quantify in Section 5.1 that WASP-19Ab is outside the regime where this is safe, with q0 ≈ 0.06 and an estimated ~19% k2 bias from Wahl et al. (2021). Because the core mass fraction is essentially mapped from k2 in the homogeneous model (Kramm et al. 2011), a 19% upward shift in k2 is a first-order perturbation to the retrieved mcore, and the 12% observational error is smaller than this systematic. The reader's weakest assumption is exactly this, and I agree with it. The appropriate verdict is a conditional acceptance: the method and test-planet validation are sound, but the WASP-19Ab core detection should not be stated as settled without propagating the non-linear tidal systematic (or demonstrating that the detection survives it). Data/code availability is a secondary reproducibility concern, not the main scientific soft spot.","tokens_in":22095,"tokens_out":4377,"duration_ms":41520,"concrete_test":"Re-run the WASP-19Ab retrieval of Section 4.2.2 with the observed k2 shifted upward by the Wahl et al. (2021) non-linear analog, e.g., using k2 = 0.238 instead of 0.20 and the same 12% relative error, with identical priors; then check whether the homogeneous-model mcore posterior still excludes zero and whether the inhomogeneous-model probability of a coreless interior (mcore and mdilute both below 0.1) stays below 1%. If the core signal disappears or weakens materially, the WASP-19Ab core detection is not robust to the systematic the authors themselves identify in Section 5.1. Ideally this is supplemented by computing the actual non-linear k2 for the retrieved density profiles with a concentric Maclaurin spheroid code (Wahl et al. 2021).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The single most load-bearing assumption is that WASP-19Ab's observed Love number k2 = 0.20+0.02/-0.03 can be compared with the linear, hydrostatic first-order k2 computed from Eq. (5) of Section 2.1.3. Section 5.1 states that WASP-19Ab has rotational parameter q0 ≈ 0.06 (assuming synchronous rotation), exceeding the q0 << 0.01 regime where the linear approximation holds, and that for WASP-12b, which has a similar q0, Wahl et al. (2021) estimate that neglecting non-linear effects underestimates k2 by about 19% — larger than the 12% observational uncertainty. WASP-19Ab was not in Wahl et al.'s sample, so this is an analog rather than a computed correction, but it is the authors' own quantitative warning. Since the homogeneous-model core mass fraction mcore = 0.21+0.05/-0.04, and the inhomogeneous-model claim that the probability of a coreless interior is only 0.02%, are driven primarily by k2, a +19% systematic shift in k2 (i.e., a less centrally condensed planet) could plausibly move the mcore posterior toward zero and change the core-detection statement. The abstract and conclusions present the core detection without this caveat, although Section 5.1 does flag it; the claim is therefore conditional on the linear-tidal interpretation holding for WASP-19Ab.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a Bayesian retrieval framework that combines mass, radius, equilibrium temperature, atmospheric metallicity, and Love number measurements to infer the interior structure of hot Jupiters, using CEPAM-based interior models in two configurations: a homogeneous two-layer model and an inhomogeneous dilute-core model. The framework is validated on a synthetic test planet with known interior properties, and the authors quantify the Love-number precision required to recover bulk metallicity and core mass, finding thresholds of approximately 40% for the homogeneous model and 15% for the dilute-core model. The framework is then applied to five hot Jupiters with measured Love numbers, of which WASP-19Ab is the only planet with both a precise Love number and an atmospheric metallicity constraint; for this planet the paper reports a homogeneous-model core mass fraction of 0.21+0.05/-0.04, corresponding to about 79 Earth masses, and a claimed 0.02% probability of a coreless interior in the inhomogeneous model.","tokens_in":22457,"tokens_out":5386,"duration_ms":53843,"significance":"If the WASP-19Ab result holds, the paper would provide the first exoplanet core detection based on a Love number measurement combined with an atmospheric metallicity constraint, which is a notable advance. The test-planet validation with injected values, the explicit comparison of two interior models, and the quantitative precision-threshold analysis are strengths that make the retrieval methodology useful for planning JWST-era observations. However, the flagship core-detection claim relies on a linear hydrostatic Love number calculation that the authors themselves identify as potentially biased for WASP-19Ab; because the potential systematic is larger than the observational uncertainty, the core detection should be tested against this systematic before the claim is treated as settled.","major_comments":[{"comment":"The WASP-19Ab core-detection claim is not robust to the non-linear tidal systematic that the authors themselves quantify. Section 5.1 states that WASP-19Ab has rotational parameter q0 approximately 0.06, which exceeds the q0 << 0.01 regime where the linear approximation holds, and that Wahl et al. (2021) estimate that neglecting non-linear effects underestimates k2 by about 19% for WASP-12b, a planet with similar q0. This is larger than the 12% observational uncertainty adopted in the retrieval, and k2 is the main constraint on core mass. The abstract and Section 4.2.2 present the core mass of 0.21+0.05/-0.04 and the coreless-probability statement without this caveat. I request that the authors add a retrieval with k2 shifted upward by approximately 19% (or with a non-linearly corrected k2) and report whether the core detection and the 0.02% probability survive; if they do not, the abstract and conclusions should be revised to present the core detection as conditional on the linear-tidal interpretation.","section":"Section 5.1, Table 3, Section 4.2.2, Abstract"},{"comment":"The claim that a coreless interior is \"effectively ruled out\" relies on a probability threshold that is not physically defined. The 0.02% figure is computed as the probability that both m_core and m_dilute are below 0.1, but a diluted-core model with these parameters can still contain substantial heavy elements in the envelope (the retrieved Z_atm for WASP-19Ab is about 0.07), so this is not a direct posterior probability for the absence of a compact core. The authors should either justify the threshold as a physically meaningful definition of coreless or replace the statement with a model-comparison metric that directly contrasts interiors with and without a compact core.","section":"Section 4.2.2"},{"comment":"The systematic tendency for the retrieved Love numbers to fall at the lower end of, or below, the observed values for most planets in the sample is discussed by the authors as a possible observational bias or missing physics. Since WASP-19Ab is the only planet with a precise enough measurement to drive a core detection, the absence of a quantitative test of how this systematic discrepancy would affect the WASP-19Ab inference is a gap. The authors should state explicitly whether the WASP-19Ab agreement is expected to be immune to the same bias, or add a sensitivity test in which the observed k2 is shifted within the plausible systematic range.","section":"Section 4.1"}],"minor_comments":[{"comment":"The notation Z_atm appears in Equation (1) while the text and Table 1 use Z_env for the envelope metal mass fraction; please unify the notation.","section":"Equation (1)"},{"comment":"The middle panel label reads \"Love number k22\" but should read \"Love number k2\".","section":"Figure 5"},{"comment":"The prior notation \"LU(0, l_age_min)/U(1, L_age_min)\" for L_grav is confusing, particularly because the text in Section 2.3 says the lower limit is set to 1 L_J; please clarify the intended ranges in the table.","section":"Table 1"},{"comment":"The sentence \"In Figures B2, B3, B4, B5 and B6 we the results obtained\" is missing the verb \"show\"; please correct it.","section":"Appendix B"},{"comment":"The phrase \"interiors structures\" should be \"interior structures\" for grammatical correctness.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid retrieval-methodology study with a useful precision-threshold analysis and a transparent validation on a test planet. The main issue is that the flagship WASP-19Ab core-detection claim depends on a linear tidal calculation whose potential systematic error the authors themselves quantify as larger than the measurement uncertainty; this is addressable within the manuscript's scope by adding a robustness test or by softening the claim. I therefore recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a solid retrieval study with one result worth remembering: combining a Love number measurement with an atmospheric metallicity constraint can pin down hot Jupiter core mass, and the required precision is about 40% for a homogeneous model and 15% for a dilute-core model. The test-planet validation is genuine and the Bayesian framework is competently put together. The application to the five known Love-number planets is useful even where it only yields upper limits.\n\nThe flagship result, however, is the WASP-19Ab core detection, and it is more fragile than the abstract lets on. The authors use a linear, hydrostatic first-order k2 calculation. In Section 5.1 they note that WASP-19Ab has rotational parameter q0 ≈ 0.06, similar to WASP-12b, for which Wahl et al. (2021) estimate that neglecting non-linear effects underestimates k2 by about 19% — larger than the 12% observational uncertainty. Since k2 is the main constraint on the core, a positive shift of that size would push the retrieved core mass toward zero. The abstract and conclusions say both models 'confirm the presence of a core' without this caveat. The authors do flag it in Section 5.1, so they are not hiding it, but the headline claim is stated more strongly than the modeling supports.\n\nThe rest of the paper holds up. The bulk metallicity results, the comparison to mass-metallicity trends, and the heating efficiency constraints are all reasonable. I also appreciate that they keep the inflation treatment flexible rather than committing to one mechanism. The absence of released code and data is a minor disappointment; \"available on reasonable request\" is weaker than it should be in 2025.\n\nMy verdict: worthwhile paper, but the WASP-19Ab core detection is conditional on the linear-tidal interpretation holding. A straightforward robustness test — re-running the retrieval with k2 shifted by +19%, or better, computing a non-linear Love number for WASP-19Ab — would settle it. If that test moves the core posterior to zero, the paper still works as a methodology and thresholds paper; the core claim just goes away. I'd send it to peer review, and I'd ask for that test before accepting the detection.","headline":"Solid retrieval framework and useful precision thresholds, but the WASP-19Ab core detection leans on a linear-tidal calculation that the authors themselves flag as potentially biased by more than the measurement error.","tokens_in":22943,"tokens_out":2714,"would_cite":true,"duration_ms":26734,"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":"WASP-19Ab's measured tidal response, combined with its atmospheric metallicity, points to a core of about 79 Earth masses.","keywords":["hot Jupiters","Love numbers","interior structure","core mass","atmospheric metallicity","Bayesian retrieval","WASP-19Ab","exoplanet interiors"],"falsifier":"A non-linear tidal calculation for WASP-19Ab (for example with the concentric MacLaurin spheroid method used for WASP-12b) that changes the predicted $k_2$ by the ~19% bias would move the observed value relative to the model grid; if the resulting posterior no longer excludes a coreless interior, the detection claim collapses. A direct JWST phase-curve measurement of $k_2$ for this planet that differs from the Bernabò et al. (2024) value by more than the combined uncertainties would likewise settle the question.","tokens_in":21902,"feed_emoji":"🪐","tokens_out":9250,"duration_ms":70982,"temperature":0.7,"pith_summary":"This paper argues that two observable constraints together—the tidal Love number $k_2$ and the atmospheric metallicity—can lift the degeneracies that usually leave hot Jupiter interior models with upper limits instead of detections. The authors build a Bayesian retrieval on a standard planetary structure code, using both a simple two-layer homogeneous model and a more complex dilute-core inhomogeneous model. Applied to the five hot Jupiters with published Love numbers, only WASP-19Ab has data precise enough to pin down the interior, and both models confirm the presence of a core. Under the homogeneous model the core mass fraction is $0.21^{+0.05}_{-0.04}$, about $79^{+21}_{-18}$ Earth masses. This matters because it shows how upcoming JWST-quality phase curves and atmospheric spectra can turn tides into a direct probe of where heavy elements sit inside giant planets.","feed_headline":"Tidal data reveal a core in hot Jupiter WASP-19Ab","feed_subtitle":"Love number plus atmospheric metallicity point to about 79 Earth masses of heavy elements.","key_machinery":"The load-bearing observable is the second-order Love number $k_2$, the surface value of the Love function $K_n(r)$ that measures how strongly a planet's gravitational potential responds to an external tide. It is computed from the interior density profile by integrating a first-order differential equation for $\\eta_2(r)$ under the assumption of hydrostatic equilibrium and a linear tidal response. The retrieval uses the planetary structure code CEPAM to generate density profiles for a homogeneous two-layer model (compact heavy-element core plus uniform H/He envelope) and an inhomogeneous model with a dilute core whose heavy-element fraction falls off through an error-function gradient; a nested-sampling Bayesian engine then fits predicted radius and $k_2$ to observations, with atmospheric metallicity entering as a Gaussian prior on the envelope metal fraction. The $k_2$ measurement carries the core-mass information because a more centrally concentrated mass distribution (larger core) lowers the Love number.","core_discovery":"The central claim is that accurate interior constraints for hot Jupiters require a Love number measurement with precision better than about 40% for a homogeneous model and 15% for a dilute-core model, and that this must be paired with an atmospheric metallicity measurement to break the remaining degeneracy. For the only planet in the current sample that meets both criteria, WASP-19Ab, the retrieved posteriors rule out a coreless interior: the probability that both the core mass fraction and the dilute-core extent are below 0.1 is only 0.02%. The homogeneous model gives a core mass fraction of $0.21^{+0.05}_{-0.04}$ ($79^{+21}_{-18}$ $M_\\oplus$), while the inhomogeneous model gives a compact core of $55^{+25}_{-29}$ $M_\\oplus$ that extends to zero but is almost never zero. The paper does not claim to distinguish a compact from a diluted core, only that some form of dense core is present.","pith_inferences":["If the non-linear tidal bias for WASP-19Ab is as large as the paper's cited estimate for the similar planet WASP-12b (~19% in $k_2$), the retrieved core mass fraction could shift noticeably; the authors themselves recommend modelling non-linear effects before trusting the exact value.","A natural extension is to apply the retrieval to the upcoming JWST Love-number and atmospheric-metallicity measurements for a larger planet sample, which would test whether the mass–metallicity trend seen in the current five planets persists.","The large spread between atmospheric and bulk metallicity found in Figure 10 suggests atmospheric metallicity alone is a weak proxy for a planet's total heavy-element content, so future atmospheric surveys should not be read as direct interior constraints."],"forward_implications":["Love number measurements with uncertainties above roughly 40% provide little interior information, so future observation campaigns should target phase-curve precisions below that threshold.","Only two current planets, HAT-P-13b and WASP-19Ab, meet the precision threshold; for the other three the retrievals return upper limits rather than detections.","WASP-19Ab becomes a benchmark for formation models: any successful giant-planet formation theory must reproduce a core of roughly 80 Earth masses inside a 1.15 Jupiter-mass planet.","Applying the same retrieval to predicted JWST phase-curve detections (~17% precision for WASP-12b-like planets) would move several planets into the regime where a core is detected rather than merely bounded.","The Love number also sharpens the constraint on the heating efficiency parameter $\\gamma$, linking tidal response to the inflation mechanism."],"supporting_citations":[{"why":"Supplies the planetary structure code CEPAM that generates the density profiles and radii used in every retrieval.","marker":"Guillot & Morel 1995"},{"why":"Provides the H/He equation of state and the inhomogeneous (dilute core) model set-up used in one of the two interior models.","marker":"Miguel et al. 2016"},{"why":"Established the direct relation between $k_2$ and core mass in two-layer models that the retrieval exploits.","marker":"Kramm et al. 2011"},{"why":"Demonstrated the retrieval approach with atmospheric metallicity constraints and the dilute-core parameterisation adopted here.","marker":"Bloot et al. 2023"},{"why":"Supplies the 12%-precision Love number for WASP-19Ab that makes the core detection possible.","marker":"Bernabò et al. 2024"},{"why":"Provides the HAT-P-13b Love number and an earlier core-mass upper limit used for comparison.","marker":"Buhler et al. 2016"},{"why":"Quantifies the non-linear tidal bias (≈19% for WASP-12b-like rotation) that defines the caveat on WASP-19Ab's result.","marker":"Wahl et al. 2021"}],"fun_headline_variants":["WASP-19Ab's core mass pinned by Love number and air chemistry","Tidal Love number reveals 79-Earth-mass core in WASP-19Ab","Core confirmed in hot Jupiter WASP-19Ab from tides and metals","Tidal data reveal WASP-19Ab's core: 79 Earth masses"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The core detection for WASP-19Ab rests on interpreting its measured Love number with a linear, hydrostatic first-order tidal model, even though the planet's rotation parameter ($q_0 \\approx 0.06$) is large enough that non-linear effects could shift the Love number by about 19%, larger than the 12% observational uncertainty.","fun_headline_variants_meta":{"raw":{"variants":["WASP-19Ab's core mass pinned by Love number and air chemistry","Tidal Love number reveals 79-Earth-mass core in WASP-19Ab","Core confirmed in hot Jupiter WASP-19Ab from tides and metals","Tidal data reveal WASP-19Ab's core: 79 Earth masses"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001426,"raw_usage":{"total_tokens":5799,"prompt_tokens":1033,"completion_tokens":4766,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":649,"completion_tokens_details":{"reasoning_tokens":4681}},"tokens_in":649,"tokens_out":4766,"duration_ms":35822,"temperature":1.0,"reasoning_tokens":4681,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:11:43.197396+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A non-linear tidal calculation for WASP-19Ab (for example with the concentric MacLaurin spheroid method used for WASP-12b) that changes the predicted $k_2$ by the ~19% bias would move the observed value relative to the model grid; if the resulting posterior no longer excludes a coreless interior, the detection claim collapses. A direct JWST phase-curve measurement of $k_2$ for this planet that differs from the Bernabò et al. (2024) value by more than the combined uncertainties would likewise settle the question.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the planetary structure code CEPAM that generates the density profiles and radii used in every retrieval."}],"review_version":1}