REVIEW 4 major objections 3 minor 58 references
Convective mixing in distant and close-in giant planets -- Dependences on the initial composition, luminosity, bloating and semi-convection
T0 review · 4 major / 3 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Evolution models of giant planets indicate that dilute cores—extended interiors enriched in heavy elements—rarely survive when the planet starts out brighter than a few thousand Jupiter luminosities, so most giants should end up fully…
desk verdict Careful and unusually honest parameter study, but the headline dilute-core conclusion rests on an unresolved mesh-resolution dependence that the paper itself concedes. 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 mechanism that carries the argument is compositional convection governed by the Ledoux criterion in a one-dimensional planetary evolution model. Convection begins where the radiative-conductive gradient $\nabla_{\mathrm{rad}}$ exceeds the sum of the adiabatic gradient and the compositional gradient $\nabla_X$; mixing is then modelled as a diffusion process on a separate equal-mass mesh of $5\times10^4$ points using mixing-length theory for the diffusion coefficient. Steep compositional steps—the 'stairs' in the heavy-element profile—are what hold a dilute core in place, because they make $\nabla_X$ large enough to suppress mixing. The paper's key numerical control is the mesh resolution, since it sets the maximum compositional gradient that can be represented, and the key physical control is the initial luminosity, since $\nabla_{\mathrm{rad}}$ grows roughly linearly with it.
What would settle it
Repeat the paper's Jupiter-like simulation at an initial luminosity of $10^4\,L_J$ with about $3\times10^5$ equally spaced mass points, the resolution implied by the overshooting-length estimate quoted in Sect. 3.6; if a dilute core then survives, the luminosity ceiling is an artifact of the coarser mesh rather than a physical limit.
Extended reading notes
Core claim
The central claim is that dilute cores—deep interior regions where heavy elements are present at moderate enrichment rather than in a pure compact core—are difficult to preserve. The paper finds that for a Jupiter-like planet with the heavy-element profiles taken from formation models, the limiting factor is the initial luminosity, which sets the radiative-conductive gradient that drives convection in the deep envelope. Above about $3\times10^3\,L_J$, mixing destroys the compositional staircase that stabilises the dilute core; at $10^4\,L_J$ the envelope mixes completely. Because hot-start and even many cold-start formation scenarios produce luminosities above this threshold, the paper concludes that retaining a dilute core through the whole evolution is an unlikely outcome for most giant planets, at least under the assumptions and resolution of the model.
Load-bearing premise
The central conclusion depends on assuming that 50,000 equal-mass layers capture the real compositional layering, even though the paper shows the survival of the dilute core changes drastically with mesh resolution and offers no physical argument that this count matches the true layering scale.
Editorial extensions
If this is right
- If dilute cores are as fragile as this model says, JWST measurements of atmospheric metallicity in giant planets can be interpreted with relatively simple core-plus-homogeneous-envelope structures for most planets.
- The first billion years matter almost exclusively: nearly all mixing happens early, so the conditions set by formation, not later evolution, decide whether a dilute core survives.
- Close-in hot Jupiters should not simply be assumed convective: bloating lowers the intrinsic luminosity and can suppress mixing, leaving slightly more dilute cores than at wide orbits, provided they formed in place.
- Strong semi-convection can shrink a dilute core and enrich the outer envelope, but it cannot fully erase a large initial core in the paper's models.
- The threshold of roughly $3\times10^3\,L_J$ provides a formation-property test: planets formed by hot accretion are expected to be fully mixed, and any observed dilute core would point to cold accretion or early radiative zones.
Reading between the lines
- If confirmed, the strong mesh dependence raises the possibility that the qualitative conclusion—that dilute cores are rare—is an artifact of unresolved layering; a physically motivated resolution would require resolving the estimated overshoot length, roughly $3\times10^5$ mesh points, which the paper did not run at the decisive luminosity.
- A corollary beyond the paper is that the atmospheric-versus-bulk metallicity mismatch should be systematically absent in young, bright giant planets and possibly present in old, dim ones.
- The bloating results imply a testable orbital-distance trend: if hot Jupiters form in situ, their interior mixing should be minimal near 0.04–0.05 AU and stronger both farther out and closer in; a survey of atmospheric metallicities versus orbital distance could look for that non-monotonic signature.
- The paper uses water for all heavy elements; if real interiors are richer in heavier molecules, the compositional gradient is larger and mixing is weaker, so the luminosity ceiling would move upward and the claim that most planets lose their dilute cores would be too strong.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript models the 4.5 Gyr evolution of giant planets with the 1D code completo21, adding Ledoux-criterion convection with composition gradients and a separate mass-mesh treatment of convective mixing. Starting from four initial heavy-element profiles (compact, extended, metal-rich, and Jupiter-like), the authors vary orbital distance and bloating, mixing length, semi-convection efficiency, opacity, mesh resolution, and initial luminosity. They find that dilute cores can be retained under some conditions, that semi-convection can shrink them, that bloating has a modest effect, and that mesh resolution strongly controls step formation and core extent. The headline claim is that dilute cores cannot persist at initial luminosities much above ~3e3 L_J for a Jupiter-mass planet, and the paper concludes that retention of dilute cores in a large fraction of giant planets is unlikely.
Significance. If the headline result were robust, it would be an important constraint for post-formation giant-planet evolution and for interpreting atmospheric versus bulk metallicities. The paper also contributes a systematic parameter study of hot Jupiters and explicitly tests the effect of an alternative luminosity distribution. The central result is not a fit: the main parameters are selected a priori, not tuned to produce dilute-core survival. However, the load-bearing claim is not numerically converged. The paper itself concedes that the appropriate number of mesh points is unclear and that step locations and sizes are mesh dependent. Because finer meshes consistently make dilute cores more persistent, the luminosity threshold in Sect. 5 could shift materially. The conclusion about the rarity of retained dilute cores is therefore conditional on an unvalidated numerical parameter; the significance is real but provisional.
major comments (4)
- [§3.6, Fig. 10 and §5] The central claim is not shown to be converged with respect to the mass mesh. At L_init = 1e3 L_J, going from 2e4 to 5e4 to 1e5 mesh points changes whether four initial profiles retain a dilute core, how many steps form, and how far the compact core extends (from ~6% to ~17% of the mass). The resolution study is only performed at the standard L_init = 1e3 L_J; the key threshold of ~3e3 L_J (Sect. 3.7) is never tested at higher resolution. Since finer meshes consistently make dilute cores more persistent, the threshold could shift upward at 1e5 or 2e5 points, potentially above the cold-start luminosities of ~2e4 to 6e4 L_J invoked in Sect. 3.7. The paper's own statement in Sect. 3.6 that 'It is unclear what number of mesh points is most realistic' and the Sect. 5 caveat 'Assuming that the number of mesh points of 5e4 ... provides a good approximation' make this a load-bearing unresolved issue, not a presentation issue.
- [§4.4, Eq. (3)] The luminosity distribution dL/dm = L/M is an ad hoc simplifying assumption, and the paper shows that the alternative dL/dm = -T dS/dt changes the early luminosity and radiative-conductive gradient by up to a factor of 3 in the 0.1-0.2 m/M region, which is precisely the usual extent of the dilute core. This alternative also reduces the number of compositional steps from 7 to 2 for the compact structure. Since the mixing criterion and the ability of a step to inhibit convection depend on the radiative-conductive gradient (Eqs. 5, 6, and 8), the simplified luminosity profile is a possible source of the reported retention threshold. At minimum, the paper should quantify how the step-size and luminosity thresholds change under a more physical entropy-based luminosity profile, or present a physical justification for why dL/dm = L/M is adequate for the central claim.
- [§4.2 and §2.1] The artificial 1 Myr delay in applying the radiative-conductive gradient for mixing purposes is a parameter that conditions the outcome. The manuscript states that without this delay, some initial compositions mix completely and lose their dilute core, whereas with the delay they retain it. The delay is justified qualitatively by hot-start accretion scenarios, but no sensitivity study is presented for its duration or the shape of the ramp. Because the central conclusion that dilute cores can persist at L_init = 1e3 L_J depends on this prescription for the compact and Jupiter-like cases, the paper needs to show that the retention does not hinge on the particular choice of 1 Myr or on the functional form of the delayed onset.
- [§3.7 and §5] The luminosity threshold itself is inferred from runs at only three initial luminosities (1e3, 3e3, and 1e4 L_J), and the paper uses the maximum luminosity reached during mixing rather than the literal initial value. The increase in luminosity during mixing is sizable (e.g., from 3e3 to up to 5.7e3 L_J for the metal-rich profile), but it is not shown how this maximum depends on mesh resolution or on the delayed-convection prescription. A combined resolution study at and above the stated threshold is needed before the conclusion 'it is unlikely that a large number of giant planets retain a dilute core' can be considered robust.
minor comments (3)
- [§2.1] There is a typo: 'the size of the the thermodynamic evolution timestep' should read 'the size of the thermodynamic evolution timestep'.
- [§4.1] The phrase 'for which be do not observe a consistent relation' should read 'for which we do not observe a consistent relation'.
- [Abstract and §3.7] The notation '3 x 1e3 LJ' is awkward; in the published version this should be typeset as 3 × 10^3 L_J consistently throughout the abstract and text.
Circularity Check
No circularity: the dilute-core retention threshold is an emergent simulation outcome, and the numerical-resolution dependence is an acknowledged limitation, not a circular step.
full rationale
The paper's central result (dilute cores do not persist at initial luminosities much above about 3 x 10^3 L_J for the studied heavy-element profiles) is an emergent outcome of integrating Eqs. (5)-(12): the radiative-conductive gradient is proportional to luminosity (Eq. 6), the compositional gradient enters through Eq. (8), and the Ledoux criterion (Eq. 5) determines where mixing occurs. The initial luminosity is a scanned input, not a parameter fitted to the retention outcome; the threshold is found, not imposed. The only calibrated value, the 0.5% bloating efficiency in Sect. 2.2, is fitted to the present-day radius of HD 209458 b and affects only the hot-Jupiter subset, not the luminosity threshold. The mesh-size dependence (Sect. 3.6, Fig. 10) is a numerical convergence limitation that the authors explicitly acknowledge with the sentence 'It is unclear what number of mesh points is most realistic'; however serious for the headline claim, a resolution artifact is not a circularity because the claim does not reduce to its inputs by construction. The formation-luminosity comparison in Sect. 3.7 uses Mordasini (2013), a prior work by a co-author, but it is corroborated by Marley et al. (2007) and is an external premise rather than an imported uniqueness claim or a fitted prediction. No load-bearing step in the derivation is equivalent to its own input, so the analysis is self-contained apart from standard, acknowledged modeling assumptions.
Assumptions & free parameters
free parameters (6)
- Mixing length parameter alpha =
1e-3 (standard), varied to 1e-2 and 1e-1
- Semi-convection efficiency alpha_sc =
0 (standard), varied 1e-2, 1e-1, 1
- Bloating efficiency for hot Jupiters =
0.5%
- Initial luminosity L_init =
1e3 L_J (standard), varied to 3e3 L_J
- Number of mesh points =
5e4 (standard), tested 2e4 and 1e5
- Convective initiation delay =
1 Myr ramp
assumptions (7)
- domain assumption Ledoux criterion (nabla_rad < nabla_ad + nabla_X) determines convective stability.
- domain assumption Convective mixing is diffusive with D = 0.1 (l/H_p) v H_p from the Mihalas (1978) efficiency estimate.
- domain assumption Opacity tables of Bell & Lin (1994) for radiative and Cassisi et al. (2007) for conductive opacity are accurate.
- domain assumption Water EOS (Haldemann et al. 2020) represents all heavy elements; H/He EOS from Chabrier & Debras (2021).
- ad hoc to paper The luminosity is distributed as dl/dm = L/M in the envelope.
- domain assumption The double-grey atmosphere model of Guillot (2010) gives the outer boundary condition.
- domain assumption Convective mixing can be numerically separated from thermal evolution and computed with a Crank-Nicolson scheme on a mass grid.
Cite this review
Pith. "Pith review of Convective mixing in distant and close-in giant planets -- Dependences on the initial composition, luminosity, bloating and semi-convection." pith.science (2026). https://pith.science/paper/HR37OUVJ
@misc{pith2026241118686,
author = {Pith},
title = {Pith review of: Convective mixing in distant and close-in giant planets -- Dependences on the initial composition, luminosity, bloating and semi-convection},
year = {2026},
howpublished = {\url{https://pith.science/paper/HR37OUVJ}},
note = {Machine review of arXiv:2411.18686}
}
read the original abstract
Recent structure models of Jupiter suggest the existence of an extended region in the deep interior with a high heavy element abundance, referred to as a dilute core. This finding has led to increased interest in modelling the formation and evolution processes with the goal of understanding how and under what circumstances such a structure is formed and retained, to in turn better understand the relation between atmospheric and bulk metallicity. We modelled the evolution of giant planets, varying various parameters relevant for the convective mixing process, such as the mixing length parameter and the size of the mesh, and parameters related to the general evolution, such as the orbital distance and the initial luminosity. We in particular studied hot Jupiters and find that the effect of bloating on the mixing process is small but can in some cases inhibit convective mixing by lowering the intrinsic luminosity for a given entropy. Semi-convection can significantly lower the extent of a dilute core if it is strong enough. We find that dilute cores are unable to persist for initial luminosities much higher than 3 x 1e3 LJ for a Jupiter-like planet for the initial heavy element profiles we studied. From this we conclude that, based on our model, it is unlikely that a large number of giant planets retain a dilute core throughout their evolution, although this is dependent on the assumptions and limitations of our method. Future work should focus on improving the link between formation and evolution models so that the mixing process is accurately modelled throughout a planet's lifetime and on improving the understanding of how to model convection near radiative-convective boundaries.
Figures
Figures from the paper (8 more)
Reference graph
Works this paper leans on
-
[1]
R., Alam, M
Alderson, L., Wakeford, H. R., Alam, M. K., et al. 2023, Nature, 614, 664
2023
-
[2]
S., Allard, F., & Hauschildt, P
Baraffe, I., Chabrier, G., Barman, T. S., Allard, F., & Hauschildt, P. H. 2003, A&A, 402, 701
2003
-
[3]
Bell, K. R. & Lin, D. N. C. 1994, ApJ, 427, 987
work page 1994
-
[4]
2024, arXiv e-prints, arXiv:2403.03325
Benneke, B., Roy, P.-A., Coulombe, L.-P., et al. 2024, arXiv e-prints, arXiv:2403.03325
arXiv 2024
- [5]
-
[6]
J., Adriani, A., Adumitroaie, V., et al
Bolton, S. J., Adriani, A., Adumitroaie, V., et al. 2017, Science, 356, 821
2017
- [7]
-
[8]
Y., Pietrinferni, A., Catelan, M., & Salaris, M
Cassisi, S., Potekhin, A. Y., Pietrinferni, A., Catelan, M., & Salaris, M. 2007, ApJ, 661, 1094
work page 2007
Show all 58 references
-
[9]
& Debras, F
Chabrier, G. & Debras, F. 2021, ApJ, 917, 4
2021
-
[10]
& Chabrier, G
Debras, F. & Chabrier, G. 2019, ApJ, 872, 100
2019
-
[11]
& Seager, S
Demory, B.-O. & Seager, S. 2011, ApJS, 197, 12
2011
-
[12]
2024, Nature, 625, 51
Dyrek, A., Min, M., Decin, L., et al. 2024, Nature, 625, 51
2024
-
[13]
2019, New Journal of Physics, 21, 023007
French, M. 2019, New Journal of Physics, 21, 023007
2019
-
[14]
& Redmer, R
French, M. & Redmer, R. 2017, Physics of Plasmas, 24, 092306
2017
-
[15]
R., Anders, E
Fuentes, J. R., Anders, E. H., Cumming, A., & Hindman, B. W. 2023, ApJ, 950, L4
2023
-
[16]
R., Hindman, B
Fuentes, J. R., Hindman, B. W., Fraser, A. E., & Anders, E. H. 2024, ApJ, 975, L1
2024
-
[17]
& Sari, R
Ginzburg, S. & Sari, R. 2016, ApJ, 819, 116
2016
-
[18]
2019, ApJ, 887, 228
Gu, P.-G., Peng, D.-K., & Yen, C.-C. 2019, ApJ, 887, 228
2019
-
[19]
1995, Science, 269, 1697
Guillot, T. 1995, Science, 269, 1697
1995
-
[20]
2005, Annual Review of Earth and Planetary Sciences, 33, 493
Guillot, T. 2005, Annual Review of Earth and Planetary Sciences, 33, 493
2005
-
[21]
2010, A&A, 520, A27 Guzmán-Mesa, A., Kitzmann, D., Fisher, C., et al
Guillot, T. 2010, A&A, 520, A27 Guzmán-Mesa, A., Kitzmann, D., Fisher, C., et al. 2020, AJ, 160, 15
2010
-
[22]
2020, A&A, 643, A105 Helled,R.&Stevenson,D.J.2024,AGUAdvances,5,e2024AV001171
Haldemann, J., Alibert, Y., Mordasini, C., & Benz, W. 2020, A&A, 643, A105 Helled,R.&Stevenson,D.J.2024,AGUAdvances,5,e2024AV001171
2020
-
[23]
Hindman, B. W. & Fuentes, J. R. 2023, ApJ, 957, L23
2023
-
[24]
2023, A&A, 672, A33
Howard, S., Guillot, T., Bazot, M., et al. 2023, A&A, 672, A33
2023
-
[25]
M., Woitke, P., & Dominik, C
Khorshid, N., Min, M., Désert, J. M., Woitke, P., & Dominik, C. 2022, A&A, 667, A147
2022
-
[26]
2012, Stellar Structure and Evolution
Kippenhahn, R., Weigert, A., & Weiss, A. 2012, Stellar Structure and Evolution
2012
- [27]
-
[28]
D., Thorngren, D
Komacek, T. D., Thorngren, D. P., Lopez, E. D., & Ginzburg, S. 2020, ApJ, 893, 36
2020
-
[29]
2009, MNRAS, 395, 1857
Kovetz, A., Yaron, O., & Prialnik, D. 2009, MNRAS, 395, 1857
2009
-
[30]
F., & Fricke, K
Langer, N., El Eid, M. F., & Fricke, K. J. 1985, A&A, 145, 179
1985
-
[31]
J., & Sugimoto, D
Langer, N., Fricke, K. J., & Sugimoto, D. 1983, A&A, 126, 207
1983
-
[32]
& Chabrier, G
Leconte, J. & Chabrier, G. 2012, A&A, 540, A20
2012
-
[33]
2017, A&A, 598, A98
Leconte, J., Selsis, F., Hersant, F., & Guillot, T. 2017, A&A, 598, A98
2017
-
[34]
2018, Nature Communications, 9, 3709
Li, L., Jiang, X., West, R., et al. 2018, Nature Communications, 9, 3709
2018
-
[35]
Lin, D. N. C., Bodenheimer, P., & Richardson, D. C. 1996, Nature, 380, 606
1996
-
[36]
2019, Nature, 572, 355
Liu, S.-F., Hori, Y., Müller, S., et al. 2019, Nature, 572, 355
2019
-
[37]
Lodders, K., Palme, H., & Gail, H. P. 2009, Landolt Börnstein, 4B, 712
2009
-
[38]
2019, ARA&A, 57, 617
Madhusudhan, N. 2019, ARA&A, 57, 617
2019
-
[39]
S., Fortney, J
Marley, M. S., Fortney, J. J., Hubickyj, O., Bodenheimer, P., & Lis- sauer, J. J. 2007, ApJ, 655, 541
2007
-
[40]
2022, A&A, 662, A18
Miguel, Y., Bazot, M., Guillot, T., et al. 2022, A&A, 662, A18
2022
-
[41]
1978, Stellar atmospheres
Mihalas, D. 1978, Stellar atmospheres
1978
-
[42]
2013, A&A, 558, A113
Mordasini, C. 2013, A&A, 558, A113
2013
-
[43]
2012, A&A, 547, A111
Mordasini, C., Alibert, Y., Klahr, H., & Henning, T. 2012, A&A, 547, A111
2012
-
[44]
2016, ApJ, 832, 41 Müller, S., Helled, R., & Cumming, A
Mordasini, C., van Boekel, R., Mollière, P., Henning, T., & Benneke, B. 2016, ApJ, 832, 41 Müller, S., Helled, R., & Cumming, A. 2020, A&A, 638, A121 Öberg, K. I., Murray-Clay, R., & Bergin, E. A. 2011, ApJ, 743, L16
2016
-
[45]
2011, ApJS, 192, 3
Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3
2011
-
[46]
2013, ApJS, 208, 4
Paxton, B., Cantiello, M., Arras, P., et al. 2013, ApJS, 208, 4
2013
-
[47]
2019, MNRAS, 487, 2653
Podolak, M., Helled, R., & Schubert, G. 2019, MNRAS, 487, 2653
2019
-
[48]
Potekhin, A. Y. 1999, A&A, 351, 787
1999
-
[49]
D., & Mollière, P
Sarkis, P., Mordasini, C., Henning, T., Marleau, G. D., & Mollière, P. 2021, A&A, 645, A79
2021
-
[50]
Stevenson, D. J. 1982, Planet. Space Sci., 30, 755
1982
-
[51]
J., Bodenheimer, P., Lissauer, J
Stevenson, D. J., Bodenheimer, P., Lissauer, J. J., & D’Angelo, G. 2022, PSJ, 3, 74
2022
-
[52]
P., Fortney, J
Thorngren, D. P., Fortney, J. J., Murray-Clay, R. A., & Lopez, E. D. 2016, ApJ, 831, 64
2016
-
[53]
J., et al
Tremblin, P., Chabrier, G., Mayne, N. J., et al. 2017, ApJ, 841, 30
2017
-
[54]
2021, ApJ, 909, 40
Turrini, D., Schisano, E., Fonte, S., et al. 2021, ApJ, 909, 40
2021
-
[55]
2018, A&A, 610, L14
Vazan, A., Helled, R., & Guillot, T. 2018, A&A, 610, L14
2018
-
[56]
2015, ApJ, 803, 32
Vazan, A., Helled, R., Kovetz, A., & Podolak, M. 2015, ApJ, 803, 32
2015
-
[57]
2016, ApJ, 829, 118
Vazan, A., Helled, R., Podolak, M., & Kovetz, A. 2016, ApJ, 829, 118
2016
-
[58]
M., Hubbard, W
Wahl, S. M., Hubbard, W. B., Militzer, B., et al. 2017, Geo- phys. Res. Lett., 44, 4649 Article number, page 13 of 13
2017
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.