Shaken, not stirred: inefficient mixing of CM- and CI-like materials
Pith reviewed 2026-05-15 13:03 UTC · model grok-4.3
The pith
Saturn's growth scatters fewer than 2 percent of CM-like planetesimals beyond 15 au because gas drag damps their eccentricities.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Outward scattering during Saturn's growth and migration implants fewer than about 2 percent of CM-like planetesimals beyond 15 au, even under gas-rich conditions, because gas drag damps their eccentricities and drives them back toward their perihelia rather than allowing them to circularize at larger distances. Adding an ice-giant core modestly increases the outward reach (up to about 4 percent in the most gas-rich case), but Type-I migration further lowers perihelia, making long-term retention at large distances difficult. For a CM mass budget M_CM,tot about 1 M_Earth, this implies at most M_CM less than 0.02-0.04 M_Earth reaches 15-25 au, corresponding to a diluted mass fraction less than
What carries the argument
N-body integrations of 100 km planetesimals launched from the outer edge of Jupiter's gap, including gas drag and gravitational perturbations of growing Jupiter and Saturn with optional nearby ice-giant embryo, across varied gas surface-density profiles and growth timescales; gas drag is the mechanism that damps eccentricities and prevents sustained outward implantation.
If this is right
- Saturn's growth scatters CM-like planetesimals efficiently yet gas drag keeps implantation beyond 15 au below 2 percent.
- An ice-giant embryo raises the implanted fraction modestly to 4 percent at most in gas-rich runs, but Type-I migration reduces perihelia and retention.
- At most 0.02-0.04 Earth masses of CM material reaches 15-25 au for a total budget of 1 Earth mass.
- This yields a diluted fraction below 1-2 times 10 to the minus 3 in the outer ring and negligible CI contamination.
- Combined with distinct radial distributions of CM- and CI-like asteroids, mixing of the carbonaceous reservoirs remains limited.
Where Pith is reading between the lines
- Weaker gas drag than modeled or different migration histories could raise implantation fractions and permit more mixing.
- Trace CM signatures in outer solar system small bodies would test the predicted isolation of the CI reservoir.
- The result supports separate formation zones for CM and CI materials with little cross-contamination across the early solar system.
Load-bearing premise
The simulations assume 100 km planetesimals experience the modeled gas-drag regime for the explored range of gas surface-density profiles and growth timescales.
What would settle it
Detection of a CM-like mass fraction exceeding roughly 0.002 in bodies or meteorites from the 15-25 au region would falsify the low implantation efficiency.
Figures
read the original abstract
A recent study suggests that CM chondrite-like planetesimals formed in the vicinity of Saturn, in a pressure bump outside the gap carved by proto-Jupiter. While a fraction of these objects was implanted into the asteroid belt as a consequence of Saturn's growth, it remains unclear whether the scattered remainder could reach the ice-giant region and mix with more distant carbonaceous reservoirs. We test whether outward scattering during Saturn's growth and migration can implant CM-like bodies onto long-lived orbits in the Uranus-Neptune region, where they could contaminate the CI reservoir. We performed N-body integrations of 100 km planetesimals launched from the outer edge of Jupiter's gap, including gas drag and the gravitational perturbations of growing Jupiter and Saturn, with optional inclusion of a nearby ice-giant embryo. We explored a range of gas surface-density profiles and growth timescales. While Saturn's growth efficiently scatters CM-like planetesimals, fewer than about 2 percent are implanted beyond 15 au, even under gas-rich conditions, because gas drag damps their eccentricities and drives them back toward their perihelia rather than allowing them to circularize at larger distances. Adding an ice-giant core modestly increases the outward reach (up to about 4 percent in the most gas-rich case), but Type-I migration further lowers perihelia, making long-term retention at large distances difficult. For a CM mass budget M_CM,tot about 1 M_Earth, this implies at most M_CM < 0.02-0.04 M_Earth reaches 15-25 au, corresponding to a diluted mass fraction < (1-2) x 10^-3 in the outer ring, hence negligible contamination of the CI reservoir. Combined with the distinct radial distributions of CM- and CI-like asteroids in the belt, these results imply limited mixing of carbonaceous reservoirs and isolation of the CI reservoir.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses N-body integrations of 100 km planetesimals launched from the outer edge of Jupiter's gap, incorporating gas drag and perturbations from growing Jupiter and Saturn (with optional ice-giant embryo), to test outward implantation during Saturn's growth and migration. It reports that fewer than ~2% of CM-like bodies reach long-lived orbits beyond 15 au even under gas-rich conditions, because gas drag damps eccentricities and returns bodies to smaller perihelia; adding an ice-giant core raises this to ~4% at most, but Type-I migration hinders retention. For a total CM mass budget of ~1 M_Earth this implies at most 0.02-0.04 M_Earth reaches the 15-25 au region, yielding a diluted mass fraction <(1-2)×10^{-3} and thus negligible contamination of the CI reservoir.
Significance. If the implantation fractions hold, the work supplies quantitative limits on radial mixing between inner (CM) and outer (CI) carbonaceous reservoirs during giant-planet growth, reinforcing the interpretation that CM and CI chondrites sample distinct formation zones with limited cross-contamination. The forward N-body approach with explicitly stated initial conditions and drag prescriptions provides a reproducible basis for the <2% figure and the derived mass constraints.
major comments (2)
- [Methods / simulation setup] Simulation methods (gas-drag implementation): the headline result of <2% implantation beyond 15 au relies on the explored range of gas surface-density profiles and growth timescales, yet it is not shown that the weakest-drag end (lowest Σ or largest effective planetesimal size) was sampled at sufficient resolution to bound the upper implantation fraction; if weaker drag permits >5% retention at 15-25 au the M_CM <0.02-0.04 M_Earth limit and negligible-contamination conclusion would not hold for all plausible disks.
- [Results section] Results on ice-giant embryo runs: the modest increase to ~4% implantation is reported only for the most gas-rich case, but the subsequent effect of Type-I migration on perihelia is described qualitatively; quantitative statistics on final semi-major-axis and eccentricity distributions after migration are needed to confirm that long-term retention at large distances remains difficult.
minor comments (2)
- [Abstract] Abstract: the reported implantation percentages lack accompanying uncertainties, particle counts, or convergence tests, which would help readers assess statistical robustness of the <2% and ~4% figures.
- [Discussion] Notation: the mass budget M_CM,tot is introduced as ~1 M_Earth without an explicit reference to the source of this value or its uncertainty range.
Simulated Author's Rebuttal
We thank the referee for their careful reading and constructive comments, which help strengthen the presentation of our results. We respond to each major comment below and indicate the revisions planned for the manuscript.
read point-by-point responses
-
Referee: [Methods / simulation setup] Simulation methods (gas-drag implementation): the headline result of <2% implantation beyond 15 au relies on the explored range of gas surface-density profiles and growth timescales, yet it is not shown that the weakest-drag end (lowest Σ or largest effective planetesimal size) was sampled at sufficient resolution to bound the upper implantation fraction; if weaker drag permits >5% retention at 15-25 au the M_CM <0.02-0.04 M_Earth limit and negligible-contamination conclusion would not hold for all plausible disks.
Authors: We thank the referee for this observation. Our parameter survey explicitly included the weakest-drag regime through the lowest gas surface densities (0.1×MMSN) and the 100 km planetesimal size, which minimizes the drag force. In these runs the implantation fraction beyond 15 au stayed below 2% (or ~4% with the ice-giant embryo), consistent with the headline result. To make the bounding of the upper limit fully transparent, we will add a supplementary figure showing implantation efficiency as a function of the effective drag parameter across the full range of Σ and growth timescales. This will confirm that even at the weakest drag the fraction remains well below 5%. revision: partial
-
Referee: [Results section] Results on ice-giant embryo runs: the modest increase to ~4% implantation is reported only for the most gas-rich case, but the subsequent effect of Type-I migration on perihelia is described qualitatively; quantitative statistics on final semi-major-axis and eccentricity distributions after migration are needed to confirm that long-term retention at large distances remains difficult.
Authors: We agree that quantitative distributions would strengthen the ice-giant embryo section. In the revised manuscript we will add a new figure (or panels to an existing figure) presenting the cumulative distributions of final semi-major axes and eccentricities for the embryo runs, both prior to and after Type-I migration. These statistics will show that the majority of bodies experience perihelion lowering sufficient to prevent stable retention beyond 15 au, thereby providing the requested quantitative support for the conclusion that long-term outward implantation remains inefficient. revision: yes
Circularity Check
No circularity: results are direct outputs of forward N-body integrations
full rationale
The paper derives its central claims (<2% implantation beyond 15 au, M_CM < 0.02-0.04 M_Earth) exclusively from explicit N-body simulations with stated initial conditions (100 km planetesimals at outer edge of Jupiter's gap), gas-drag prescriptions (Epstein/Stokes regimes), and a scanned parameter space of gas surface-density profiles plus growth timescales. No step defines a quantity in terms of its own output, renames a fitted parameter as a prediction, or relies on a self-citation chain for a uniqueness theorem or ansatz. The implantation statistics are computed quantities, not tautological re-expressions of the inputs.
Axiom & Free-Parameter Ledger
free parameters (2)
- gas surface-density profiles
- growth timescales
axioms (1)
- standard math Newtonian N-body gravity plus standard aerodynamic gas drag on 100 km bodies
Reference graph
Works this paper leans on
- [1]
-
[2]
Anderson, S. E., Vernazza, P., & Brož, M. 2025, Nat. Astro [arXiv:2507.22649]
-
[3]
Baguet, D., Morbidelli, A., & Petit, J. M. 2019, Icarus, 334, 99 Benítez-Llambay, P., Masset, F., Koenigsberger, G., & Szulágyi, J. 2015, Nature, 520, 63
work page 2019
-
[4]
Bierson, C. J. & Nimmo, F. 2019, Icarus, 326, 10
work page 2019
-
[5]
Bitsch, B., Johansen, A., Lambrechts, M., & Morbidelli, A. 2015, A&A, 575, A28
work page 2015
-
[6]
Bottke, W. F., Nesvorný, D., Grimm, R. E., Morbidelli, A., & O’Brien, D. P. 2006, Nature, 439, 821
work page 2006
-
[7]
Bradley, J. P. 1994, Geochim. Cosmochim. Acta, 58, 2123 Brož, M., Vernazza, P., Marsset, M., et al. 2024, A&A, 689, A183
work page 1994
- [8]
- [9]
-
[10]
Carter, P. J. & Stewart, S. T. 2020, Planet. Sci. J., 1, 45
work page 2020
- [11]
-
[12]
Cornejo, S., Masset, F. S., & Sánchez-Salcedo, F. J. 2023, MNRAS, 523, 936
work page 2023
- [13]
- [14]
-
[15]
Davidsson, B. J. R., Sierks, H., Güttler, C., et al. 2016, A&A, 592, A63 de Sousa, R. R., Morbidelli, A., Raymond, S. N., et al. 2020, Icarus, 339, 113605
work page 2016
-
[16]
Deienno, R., Morbidelli, A., Gomes, R. S., & Nesvorný, D. 2017, AJ, 153, 153
work page 2017
- [17]
-
[18]
DeMeo, F. E. & Carry, B. 2013, Icarus, 226, 723
work page 2013
-
[19]
J., Kalyaan, A., & Alexander, C
Desch, S. J., Kalyaan, A., & Alexander, C. M. O. 2018, ApJS, 238, 11
work page 2018
-
[20]
Dodson-Robinson, S. E. & Bodenheimer, P. 2010, Icarus, 207, 491
work page 2010
- [21]
-
[22]
Emery, J. P., Binzel, R. P., Britt, D. T., et al. 2024, Space Sci. Rev., 220, 28
work page 2024
-
[23]
Fendyke, S. M. & Nelson, R. P. 2014, MNRAS, 437, 96
work page 2014
-
[24]
Gomes, R. S. 2003, Icarus, 161, 404
work page 2003
-
[25]
Gonzalez, J. F., Laibe, G., Maddison, S. T., Pinte, C., & Ménard, F. 2015, Planet. Space Sci., 116, 48
work page 2015
-
[26]
M., Fortier, A., Brunini, A., & Benvenuto, O
Guilera, O. M., Fortier, A., Brunini, A., & Benvenuto, O. G. 2011, A&A, 532, A142
work page 2011
-
[27]
Guillot, T., Ida, S., & Ormel, C. W. 2014, A&A, 572, A72
work page 2014
-
[28]
Hartmann, L., Calvet, N., Gullbring, E., & D’Alessio, P. 1998, ApJ, 495, 385
work page 1998
- [29]
-
[30]
Hellmann, J. L., Schneider, J. M., Wölfer, E., et al. 2023, ApJ, 946, L34
work page 2023
-
[31]
J., Brunetto, R., Cruikshank, D
Holler, B. J., Brunetto, R., Cruikshank, D. P., et al. 2025, RNAAS, 9, 241
work page 2025
-
[32]
2022, Science Advances, 8, eadd8141
Hopp, T., Dauphas, N., Abe, Y ., et al. 2022, Science Advances, 8, eadd8141
work page 2022
- [33]
-
[34]
Izidoro, A., Morbidelli, A., Raymond, S. N., Hersant, F., & Pierens, A. 2015, A&A, 582, A99
work page 2015
-
[35]
M., Lacerda, P., & Bizzarro, M
Johansen, A., Low, M.-M. M., Lacerda, P., & Bizzarro, M. 2015, Science Ad- vances, 1, e1500109
work page 2015
-
[36]
Kaib, N. A. & Quinn, T. 2008, Icarus, 197, 221
work page 2008
- [37]
-
[38]
S., Burkhardt, C., Budde, G., & Kleine, T
Kruijer, T. S., Burkhardt, C., Budde, G., & Kleine, T. 2017, Proceedings of the National Academy of Sciences, 114, 6712
work page 2017
-
[39]
Lambrechts, M., Johansen, A., & Morbidelli, A. 2014, A&A, 572, A35
work page 2014
-
[40]
Lau, T. C. H., Birnstiel, T., Dr˛ a˙ zkowska, J., & Stammler, S. M. 2024, A&A, 688, A22
work page 2024
-
[41]
Lee, M. H. & Peale, S. J. 2002, ApJ, 567, 596
work page 2002
-
[42]
Levison, H. F., Bottke, W. F., Gounelle, M., et al. 2009, Nature, 460, 364
work page 2009
-
[43]
Levison, H. F. & Duncan, M. J. 1997, Icarus, 127, 13
work page 1997
-
[44]
F., Morbidelli, A., Vanlaerhoven, C., Gomes, R., & Tsiganis, K
Levison, H. F., Morbidelli, A., Vanlaerhoven, C., Gomes, R., & Tsiganis, K. 2008, Icarus, 196, 258
work page 2008
- [45]
-
[46]
Lin, D. N. C. & Papaloizou, J. 1986, ApJ, 309, 846
work page 1986
-
[47]
Meech, K. & Raymond, S. N. 2020, in Planetary Astrobiology, ed. V . S. Mead- ows, G. N. Arney, B. E. Schmidt, & D. J. Des Marais, 325
work page 2020
- [48]
-
[49]
Morbidelli, A., Levison, H. F., Tsiganis, K., & Gomes, R. 2005, Nature, 435, 462
work page 2005
-
[50]
Nakamura, T., Matsumoto, M., Amano, K., et al. 2023, Science, 379, abn8671
work page 2023
- [51]
-
[52]
Nelson, R. P., Papaloizou, J. C. B., Masset, F., & Kley, W. 2000, MNRAS, 318, 18 Nesvorný, D. 2018a, ARA&A, 56, 137 Nesvorný, D. 2018b, ARA&A, 56, 137 Nesvorný, D., Dauphas, N., V okrouhlický, D., Deienno, R., & Hopp, T. 2024, EPSL, 626, 118521
work page 2000
- [53]
-
[54]
J., Baruteau, C., Crida, A., & Kley, W
Paardekooper, S. J., Baruteau, C., Crida, A., & Kley, W. 2010, MNRAS, 401, 1950
work page 2010
-
[55]
Paardekooper, S. J., Baruteau, C., & Kley, W. 2011, MNRAS, 410, 293
work page 2011
-
[56]
Perets, H. B. & Murray-Clay, R. A. 2011, AJ, 733, 56 Article number, page 12 S.E. Anderson, P. Vernazza, and M. Brož: Shaken, not stirred: inefficient mixing of CM- and CI-like materials
work page 2011
- [57]
-
[58]
Pitjeva, E. V . & Pitjev, N. P. 2018, Astronomy Letters, 44, 554
work page 2018
-
[59]
B., Hubickyj, O., Bodenheimer, P., et al
Pollack, J. B., Hubickyj, O., Bodenheimer, P., et al. 1996, Icarus, 124, 62
work page 1996
-
[60]
Raymond, S. N. & Izidoro, A. 2017, Science Advances, 3, e1701138
work page 2017
-
[61]
Raymond, S. N. & Izidoro, A. 2017, Icarus, 297, 134
work page 2017
-
[62]
Raymond, S. N. & Nesvorný, D. 2022, in Vesta and Ceres. Insights from the Dawn Mission for the Origin of the Solar System, 227
work page 2022
- [63]
-
[64]
Rivkin, A. S. 2012, Icarus, 221, 744
work page 2012
- [65]
-
[66]
Ronnet, T., Mousis, O., Vernazza, P., Lunine, J. I., & Crida, A. 2018, AJ, 155, 224
work page 2018
-
[67]
Shakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337
work page 1973
-
[68]
Sharkey, B. N. L., Rivkin, A. S., Cartwright, R. J., et al. 2025, Planet. Sci. J., 6, 242
work page 2025
-
[69]
2024, Science Advances, 10, eadp2426
Spitzer, F., Kleine, T., Burkhardt, C., et al. 2024, Science Advances, 10, eadp2426
work page 2024
-
[70]
Supulver, K. D. & Lin, D. 2000, Icarus, 146, 525
work page 2000
- [71]
-
[72]
Tsiganis, K., Gomes, R., Morbidelli, A., & Levison, H. F. 2005, Nature, 435, 459
work page 2005
- [73]
-
[74]
Vernazza, P., Simon, P. N., Jorda, L., et al. 2025, A&A, 699, A167
work page 2025
-
[75]
J., Morbidelli, A., Raymond, S
Walsh, K. J., Morbidelli, A., Raymond, S. N., O’Brien, D. P., & Mandell, A. M. 2011, Nature, 475, 206
work page 2011
-
[76]
Warren, P. H. 2011, EEPSL, 311, 93
work page 2011
- [77]
-
[78]
Zega, T. J., McCoy, T. J., Russell, S. S., et al. 2025, Nature Geoscience, 18, 832 Article number, page 13 A&A proofs:manuscript no. aa57394-25 Appendix A: Collisional Efficiencies A.1. Planetesimal-planetesimal collisions. Our simulations neglect planetesimal-planetesimal collisions, which could, in principle, modify implantation outcomes by fragment- in...
work page 2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.