{"id":"a6923ba7-bad5-4aa7-9097-c5759c6ce750","arxiv_id":"2501.09963","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A coupled interior-structure and evolution model applied to GJ 486b predicts a core 1.34 times Earth's size, ~0.93 Myr mantle solidification, and a water- and oxygen-rich secondary atmosphere.","lead":"Scientists used a new 1-D computer model, SERPINT, to predict what the interior of the super-Earth GJ 486b looks like and how it cooled over time. The model suggests the planet has a core about 1.34 times larger than Earth's, a central pressure near 1171 GPa, and formed a thick, water-rich atmosphere early in its history.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Headline structural and thermal numbers are single-point outputs for M=3.00 M⊕, R=1.343 R⊕, yet the paper itself cites Meier et al. (2024) with M=2.81 M⊕, R=1.31 R⊕ yielding a 4900 km core; central claims are thus not robust to adopted mass-radius values.","rationale":"The reader's weakest assumption was the Earth-like composition and 20% impurity calibration. That is a valid modeling caveat, but the paper is explicit about it, and a compositional grid is a standard follow-up rather than an internal flaw. The more pressing issue is that even within the stated Earth-like framework, the headline numbers depend on the chosen M and R, and the paper's own text provides an alternative published parameter set that changes the core radius by ~6%. Since hydrostatic equilibrium is solved for given M and R, this is a direct, checkable dependency, not a philosophical degeneracy. The PREM validation in Sec. 3.1 supports the numerical solver, so the concern is not about the code's internal consistency but about the presentation of point estimates as definitive results. The paper's failures to propagate uncertainties and to present the thermal results as a range (Appendix C) compound this. The appropriate fix is to recompute with the adopted modern values, add a sensitivity table, and phrase the conclusions conditionally; hence CONDITIONAL remains the right verdict, with revisions required before acceptance.","tokens_in":27426,"tokens_out":15038,"duration_ms":140686,"concrete_test":"Recompute the SERPINT structure module with M = 2.81 M⊕ and R = 1.31 R⊕, keeping all EOS parameters and the 20% core impurity fixed; record the new core radius and central pressure. If the core radius differs from 4600 km by more than 5% (i.e., outside [4370, 4830] km) or the central pressure differs by more than 10% from 1171 GPa, then the headline numbers must be presented as conditional on the specific adopted mass-radius. Also rerun the thermal module at the alternative parameters with 10 EO initial water to check whether the 0.93 Myr solidification time changes by more than a factor of 2.","verdict_should_be":"UNCHANGED","load_bearing_attack":"GJ 486b's adopted mass and radius are load-bearing inputs: the structural module solves hydrostatic equilibrium for core radius and central pressure given M and R, so any change in these inputs propagates directly into the quoted results. The paper uses M=3.00 M⊕, R=1.343 R⊕ (citing Caballero et al. 2022), but in Sec. 3.2 it reports that Meier et al. (2024) use M=2.81 M⊕, R=1.31 R⊕ (Trifonov et al. 2021; Caballero et al. 2022) and obtain a core radius of 4900 km, a ~6.5% shift from the paper's 4600 km. The abstract and conclusions present 'core 1.34 times larger than Earth' and '1171 GPa' as the results, without noting that these numbers would become ~1.43× and a different central pressure under the alternative published parameters. No uncertainty propagation or sensitivity test is provided, so the reader cannot tell whether the quoted precision reflects the true dependence on input parameters. The thermal evolution result is likewise a single point: Appendix C shows solidification time varying from 0.02 Myr (0.2 EO) to 1.82 Myr (20 EO); the 0.93 Myr figure is simply the 10-EO case, not a constrained prediction. If the input mass-radius values shift, the core radius, central pressure, and solidification time all shift, undermining the specific claims as stated.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces SERPINT, a one-dimensional coupled interior-structure and thermal-atmosphere evolution model, and applies it to the super-Earth GJ 486b. Assuming an Earth-like, fully differentiated interior with a 20 wt% light-element impurity in the core and an initial budget of 10 Earth oceans of water, the model predicts a core radius of about 4600 km (1.34 times Earth's), a central pressure of about 1171 GPa, and mantle solidification after about 0.93 Myr. The model also predicts a water-rich steam atmosphere that later becomes enriched in oxygen through photolysis and hydrogen escape, which the authors connect to the water-rich interpretation of JWST transmission spectra by Moran et al. (2023). The structural model is validated against PREM for Earth. The central quantitative claims, however, depend on adopted planetary parameters and on several unconstrained model inputs, and the manuscript does not propagate uncertainties or systematically test the sensitivity of the headline numbers.","tokens_in":27782,"tokens_out":9923,"duration_ms":96396,"significance":"If the results are robust, the paper would provide a useful end-to-end modeling framework linking the interior structure, thermal evolution, and early atmospheric composition of a JWST-observed super-Earth, and it makes concrete, falsifiable predictions about core size, central pressure, and solidification timescale. The PREM validation of the structural module is a genuine strength, as is the transparent presentation of the governing equations and the explicit sensitivity study of the initial water budget in Appendix C. However, the central predictions are conditional on at least three unconstrained or loosely constrained inputs: the adopted mass and radius of GJ 486b, the 20% core impurity fraction calibrated to Earth, and the assumed initial water inventory. Because the manuscript presents the resulting numbers as the paper's headline results, the significance of the specific predictions is limited until those dependencies are quantified.","major_comments":[{"comment":"The adopted mass and radius of GJ 486b, 3.00 M⊕ and 1.343 R⊕, are load-bearing inputs to the structural module, yet the paper itself quotes the alternative published values of 2.81 M⊕ and 1.31 R⊕ when comparing with Meier et al. (2024) and reports that those values give a core radius of 4900 km rather than 4600 km. The abstract and conclusions present the 4600 km core and 1171 GPa central pressure as the results without stating which published parameter set they use or propagating the published uncertainties. Since the hydrostatic structure equations (Eqs. 2-4) use M and R as boundary conditions, the quoted structural predictions are not robust to this parameter choice; the authors should either justify the adopted values relative to the alternative, or present results for both published parameter sets and for the marginalized mass-radius posteriors.","section":"Sec. 1 and Sec. 3.2"},{"comment":"The thermal-evolution headline, 'the planet's mantle cools and solidifies over approximately 0.93 million years,' is a single-point output for the assumed initial water inventory of 10 Earth oceans. Appendix C shows that this time ranges from 0.02 Myr at 0.2 EO to 1.82 Myr at 20 EO, i.e., an order-of-magnitude sensitivity. The abstract and conclusions do not state that 0.93 Myr is conditional on an assumed, not observationally inferred, water budget. The initial water inventory should be treated as a free parameter with a stated prior or constrained from formation models, and the conclusions should present the solidification time as a function of that parameter rather than as a point prediction.","section":"Sec. 2.5.1 and Appendix C"},{"comment":"The 20 wt% core impurity factor f is calibrated by fitting the model to Earth's PREM profile, and the same f is then applied to GJ 486b with no sensitivity analysis. The core radius and central pressure depend directly on the core EOS via the impurity factor (Appendix A, Eq. A5), so a different—and equally plausible—impurity fraction would change the predicted 4600 km core radius and 1171 GPa central pressure. The authors should vary f over a reasonable range (e.g., 10-30 wt%) and report the resulting spread in the structural predictions before presenting these values as definitive results.","section":"Sec. 3.1 and Sec. 2.2.1"},{"comment":"Two equations in the methods section need correction. First, Eq. (1) as written yields g = 4πGρr for a constant-density sphere, which is a factor of 3 larger than the correct hydrostatic result g = (4πG/3)ρr. If this equation were implemented literally, the reported PREM validation (central pressure 371 GPa for Earth) would not hold, so it is likely a typographical error in the displayed formula; nevertheless the equation must be fixed because it is central to the structure calculation. Second, Eq. (13) states that P is in GPa, but the numerical coefficients are consistent with P in Pa; for example, the quoted 11747 K condition at the core-mantle boundary follows from the lower-mantle solidus only if P is in Pa. This unit inconsistency must be resolved in the text, as it directly affects the solidification-radius evolution and the reported thermal history.","section":"Eq. (1) and Eq. (13)"}],"minor_comments":[{"comment":"The reported densities are missing a factor of 1000: Earth's central density should be about 1.277×10^4 kg m^-3, not 12.77 kg m^-3, and GJ 486b's central density should be about 1.757×10^4 kg m^-3, not 17.57 kg m^-3.","section":"Sec. 3.1 and Sec. 3.2"},{"comment":"The comparison with Meier et al. (2024) reads 'a core size of 0.586 R⊕, corresponding to 4900 km'; since 0.586 × R⊕ is about 3734 km, the intended quantity is presumably 0.586 Rp (planet radius), giving 4900 km for Rp ≈ 1.31 R⊕. This should be clarified.","section":"Sec. 3.2"},{"comment":"The text says saturation timescales in the 7-9 Gyr range were sampled, while Sec. 3.3 states that the XUV luminosity decreases after one billion years, when the saturation time is exceeded. These statements are inconsistent and should be reconciled.","section":"Sec. 2.3.1 and Sec. 3.3"},{"comment":"The caption lists '238U, 238U, 238Th and 40K'; the second species should be 235U.","section":"Fig. 2 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper's internal comparison with Meier et al. (2024) makes the mass-radius sensitivity directly visible but does not resolve it; I would ask the authors to run their structure and evolution model for the alternative published parameters and report both cases, rather than only the adopted set. The unit errors in Eqs. (1) and (13) are easily fixable but should be treated seriously because they affect the method sections that readers will use to reproduce the model."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper builds SERPINT, a 1-D coupled interior structure, thermal evolution, volatile, and escape model, and applies it to GJ 486b. The model is largely assembled from existing parameterizations (Schaefer et al. 2016, Barth et al. 2021, Krissansen-Totton & Fortney 2022, Seager et al. 2007), but putting them together in one code with a PREM validation for Earth is genuinely useful. The code gives concrete numbers: core radius 4600 km, central pressure 1171 GPa, solidification time 0.93 Myr, and a water- and oxygen-rich secondary atmosphere. These are the first SERPINT-style results for this planet, so there is real value here for people working on magma ocean evolution and secondary atmospheres.\n\nNow the soft spots, in proportion. The biggest one: the adopted mass and radius of 3.00 M⊕ and 1.343 R⊕ are internally inconsistent with values cited in the paper itself. Sec. 3.2 notes that Meier et al. (2024) use 2.81 M⊕ and 1.31 R⊕ and get a 4900 km core, while the paper's headline is 4600 km. That is a ~6.5% shift in the central structural result with no sensitivity analysis or uncertainty propagation. The abstract and conclusions present single numbers as predictions. The thermal evolution result is equally fragile: solidification time ranges from 0.02 Myr at 0.2 EO to 1.82 Myr at 20 EO, and the 0.93 Myr is simply the 10 EO case. The initial water inventory is an input assumption, not a constraint. Similarly, the 20% core impurity is calibrated to Earth's PREM profile, so applying it to GJ 486b transfers a free parameter from Earth to the target without any planetary data.\n\nThe water-rich atmosphere conclusion is also stated more strongly than the data warrant. The paper says it is consistent with Moran et al. (2023), but that paper explicitly leaves stellar contamination as an alternative explanation. The model does not rule that out. To be fair, the authors do include a limitations section that acknowledges many simplifications: line-by-line climate, ocean formation, 3D convection, and core-mantle coupling are excluded. That honesty counts in their favor.\n\nWho is this for? Groups studying magma ocean solidification and volatile outgassing on super-Earths, or anyone wanting a compact, well-documented numerical testbed. The framework deserves a serious referee; the specific claims about GJ 486b need major revision: update the mass-radius to a consistent set, run a sensitivity grid, and present the headline numbers as a nominal case with a range, not as firm predictions.\n\nMy recommendation: send it to peer review with a request for major revision. The model is worth the referees' time, and the issues are correctable.","headline":"Useful coupled interior-thermal-volatile framework, but the headline numbers for GJ 486b rest on inconsistent mass-radius inputs and a chosen water inventory, so the specific claims need revision before they can be trusted.","tokens_in":28393,"tokens_out":1437,"would_cite":false,"duration_ms":16223,"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":"This paper argues that GJ 486b, modeled with an Earth-like fully differentiated interior, has a core 1.34 times larger than Earth's, a central pressure of about 1171 GPa, a magma ocean that solidifies in roughly 0.93 million years, and an…","keywords":["super-Earths","GJ 486b","planetary interior structure","thermal evolution","magma ocean","atmospheric escape","secondary atmosphere","JWST observations"],"falsifier":"Recompute the interior with a different core impurity fraction, such as 10% or 30%, or with a solid inner core, and check whether any Earth-like parameter choice reproduces both the planet's bulk density and a central pressure consistent with experimental iron equations of state at terapascal pressures; failure would rule out the predicted 4600 km core and 1171 GPa central pressure. On the observational side, a JWST spectrum that definitively attributes the 2.87–5.14 µm slope to starspots rather than water would falsify the specific water-rich evolutionary branch.","tokens_in":27180,"feed_emoji":"🪐","tokens_out":7802,"duration_ms":74145,"temperature":0.7,"pith_summary":"GJ 486b is a warm, 3-Earth-mass planet orbiting an M dwarf, and its interior cannot be observed directly. This paper asks what the planet would look like inside if its composition and layering resemble Earth's, and what that interior would do as it cooled. The model predicts an iron-rich core about 1.34 times the size of Earth's, a central pressure near 1171 GPa, and a magma ocean that solidifies within about 0.93 million years. As the ocean freezes, water is outgassed and then split by starlight, leaving a thick water- and oxygen-rich atmosphere. The authors present this as a self-consistent picture that connects JWST's ambiguous water detection to the planet's unseen interior.","feed_headline":"Model finds GJ 486b's core is 1.34 times Earth's size","feed_subtitle":"Earth-like interior gives a 1171 GPa center, a 0.93-Myr magma ocean, and a water-and-oxygen sky.","key_machinery":"The load-bearing object is SERPINT (Structure and Evolution model for Rocky Planet INTeriors), a 1-D spherical model that couples four modules: structure, thermal, volatile, and escape. The structure module solves hydrostatic equilibrium with layer-specific equations of state: a Holzapfel EOS for liquid ε-Fe modified by an impurity factor for the core, and a Vinet EOS for a two-layer silicate mantle with phase-transition boundaries at 23.5 GPa and 125 GPa. The thermal module evolves mantle potential temperature, melt fraction, solidification radius, and viscosity using radiogenic heating and a parameterized solidus-liquidus profile. The volatile and escape modules move water between melt, solid, and atmosphere, and couple photolysis, hydrogen and oxygen escape through energy-limited XUV flux, and crossover-mass drag. The identity that carries the argument is self-consistency: structure sets pressures and radii, which set solidus conditions, which set outgassing, which sets atmospheric greenhouse and escape, which feeds back on cooling.","core_discovery":"The paper's central claim is that an Earth-like structural template—a fully liquid iron core carrying 20% by weight of light elements (silicon, oxygen, sulfur) and a two-layer silicate mantle with no crust or water layer—can account for GJ 486b's observed mass and radius and yields a specific interior: a core-radius fraction of 0.537, a core radius near 4600 km, and a central pressure of about 1171 GPa, roughly three times Earth's. Coupled to this structure, the thermal model starts the planet as a 5000 K magma ocean and finds that the mantle cools and becomes rheologically solid in about 0.93 million years, with the lower mantle freezing first. During that freeze, water dissolved in the melt is expelled into the atmosphere, building a steam atmosphere that photolyzes under XUV irradiation; hydrogen escapes while oxygen is first consumed oxidizing FeO to Fe2O3 in the melt and then, once the melt is oxidized, accumulates in the atmosphere. The paper argues this scenario matches the water-rich reading of the JWST transmission spectrum and predicts a secondary atmosphere whose composition could be checked by future observations.","pith_inferences":["The 20% by weight core impurity and the 10-Earth-ocean initial water budget are calibrated to Earth, not to GJ 486b; the predicted core radius, central pressure, and solidification time should shift if stellar abundance data or formation models constrain different values.","Because the model treats the core as fully liquid and imposes a no-flux core-mantle boundary, it omits core cooling and inner-core solidification; including those would change the heat budget and could lengthen or shorten the magma-ocean lifetime.","The same coupled machinery could be run on other hot super-Earths with measured mass and radius, and the ratio of predicted outgassed water to observed atmospheric water would test whether Earth-like interiors are typical or just one solution among many.","A direct observational test is to search for O2 or other photolysis products in the planet's transmission spectrum; because the model predicts oxygen accumulation only after the melt's FeO is oxidized, an O2 detection would trace the interior oxidation state, not just the atmosphere."],"forward_implications":["If GJ 486b is Earth-like inside, its center sits at roughly 1171 GPa, more than three times Earth's central pressure, placing concrete demands on equations of state for iron-silicon-oxygen-sulfur mixtures at terapascal pressures.","The mantle would have solidified within about a million years, so any magma ocean is long gone; today's surface should be in a solid-state, likely stagnant-lid convective regime.","The early outgassing sequence predicts a transient steam atmosphere with water pressures rising to roughly 1275 bars, followed by an atmosphere enriched in oxygen as hydrogen escapes.","Changing the initial water inventory from 0.2 to 20 Earth oceans shifts the solidification time from 0.02 to 1.82 million years, making water content a first-order control on early thermal history.","JWST observations designed to detect trace gases and improve wavelength coverage could distinguish the modeled water-and-oxygen secondary atmosphere from the stellar-contamination alternative."],"supporting_citations":[{"why":"It supplies GJ 486b's mass, radius, orbital parameters, and context as a prime JWST target.","marker":"Caballero et al. 2022"},{"why":"It provides the JWST NIRSpec transmission spectrum and the water-rich versus starspot interpretation that the atmospheric predictions are compared against.","marker":"Moran et al. 2023"},{"why":"It provides the PREM reference profile used to validate the structural model and calibrate the core impurity.","marker":"Dziewonski & Anderson 1981"},{"why":"It supplies the mantle equation-of-state parameters and the layered interior modeling approach.","marker":"Seager et al. 2007"},{"why":"It provides the Holzapfel EOS parameters and the choice of liquid ε-Fe for the core.","marker":"Hakim et al. 2018"},{"why":"It supplies the magma-ocean thermal and volatile reservoir equations for water outgassing, FeO oxidation, and escape.","marker":"Schaefer et al. 2016"},{"why":"It supplies the radiogenic heating, solidus-liquidus, melt fraction, and viscosity parameterizations.","marker":"Lebrun et al. 2013"},{"why":"It supplies the energy-limited XUV escape formulation and crossover mass condition used for hydrogen and oxygen loss.","marker":"Luger & Barnes 2015"},{"why":"It supplies the stellar evolution grid used for the M dwarf's bolometric and XUV luminosity evolution.","marker":"Baraffe et al. 2015"},{"why":"It provides the coupled thermal-atmospheric evolution framework and comparative solidification timescales.","marker":"Barth et al. 2021"}],"fun_headline_variants":["GJ 486b's core: 1.34x Earth's size, 1171 GPa pressure","Magma ocean on GJ 486b froze in 0.93 million years","Water and oxygen build GJ 486b's secondary atmosphere","Interior model predicts GJ 486b's core and sky","JWST can test GJ 486b's oxygen-rich atmosphere"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model assumes GJ 486b's interior is Earth-like in the specific sense of a fully differentiated planet with an entirely liquid iron core containing 20% by weight light elements and a two-layer silicate mantle with no crust or water layer, with that 20% impurity calibrated to Earth's PREM profile rather than to any observation of GJ 486b.","fun_headline_variants_meta":{"raw":{"variants":["GJ 486b's core: 1.34x Earth's size, 1171 GPa pressure","Magma ocean on GJ 486b froze in 0.93 million years","Water and oxygen build GJ 486b's secondary atmosphere","Interior model predicts GJ 486b's core and sky","JWST can test GJ 486b's oxygen-rich atmosphere"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000582,"raw_usage":{"total_tokens":2799,"prompt_tokens":1068,"completion_tokens":1731,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":684,"completion_tokens_details":{"reasoning_tokens":1631}},"tokens_in":684,"tokens_out":1731,"duration_ms":15710,"temperature":1.0,"reasoning_tokens":1631,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T19:30:04.678987+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the interior with a different core impurity fraction, such as 10% or 30%, or with a solid inner core, and check whether any Earth-like parameter choice reproduces both the planet's bulk density and a central pressure consistent with experimental iron equations of state at terapascal pressures; failure would rule out the predicted 4600 km core and 1171 GPa central pressure. On the observational side, a JWST spectrum that definitively attributes the 2.87–5.14 µm slope to starspots rather than water would falsify the specific water-rich evolutionary branch.","supporting_citations":[{"cited_title":"2007, The Astrophysical Journal, 669, 1279","cited_arxiv_id":null,"evidence_quote":"It supplies the mantle equation-of-state parameters and the layered interior modeling approach."},{"cited_title":"2015, Astrobiology, 15, 119","cited_arxiv_id":null,"evidence_quote":"It supplies the energy-limited XUV escape formulation and crossover mass condition used for hydrogen and oxygen loss."}],"review_version":1}