Recognition: unknown
Terrestrial planet formation in the era of GPU computing
Pith reviewed 2026-05-10 17:17 UTC · model grok-4.3
The pith
Acceleration factors in N-body simulations distort the chemical composition of forming terrestrial planets and should be avoided.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In state-of-the-art N-body simulations performed with the GENGA code, the commonly used acceleration factor f, employed to speed up the collisional growth of planets, should be avoided since it leads to distorted chemical compositions. The formation timescale depends on the size of the initial planetesimals, and terrestrial planets can form resonant chains without the need for orbital migration due to gas effects.
What carries the argument
Symplectic integrators for long-term orbits and short-term collisions, unified and parallelized in the GENGA GPU/CPU code.
If this is right
- Planetary compositions in simulations become chemically consistent only when the acceleration factor is omitted.
- Larger starting planetesimals extend the time needed to assemble terrestrial planets.
- Resonant chains of planets arise from mutual gravitational scattering without requiring a gas disk.
- Higher-resolution runs produce more stable predictions of final planet numbers and spacings.
Where Pith is reading between the lines
- Models aiming to match meteorite chemistry should use only unaccelerated collision runs.
- The no-migration resonant-chain result may explain the compact, resonant architectures seen in many exoplanet systems.
- Systematic variation of initial planetesimal sizes across a wider range could reveal how disk conditions control final planet multiplicity.
Load-bearing premise
The reported outcomes rest on particular choices of initial planetesimal sizes, particle counts, and the omission of gas in some runs.
What would settle it
Run identical high-resolution simulations both with and without the acceleration factor and compare the final iron-to-silicate or volatile ratios in the planets against measured solar-system or exoplanet compositions.
read the original abstract
In this chapter, we summarize the underlying numerical methods needed for efficient $N$-body integration of planetary systems. We discuss how symplectic integrators have been developed to tackle the complementary problems of long-term orbital integration and short-term collisional interactions. The public code GENGA, a parallel GPU/CPU planet formation and orbital dynamics simulation code, was developed to unify these methods and take full advantage of the newest available computing hardware. We present state-of-the-art N-body simulations performed with GENGA in a comparative study regarding the basic properties that emerge during the late stages of the terrestrial planet formation process. We show that in modern N-body simulations the commonly used acceleration factor f, used to speed up the collisional growth of planets in simulations, should be avoided since it can lead to distorted chemical composition of the planets. We make a detailed comparison of low to high-resolution simulations, showing that the formation time scale depends on the size of the initial planetesimals. These simulations also show that terrestrial planets can form resonant chains without the need of orbital migration due to gas effects.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript summarizes numerical methods for efficient N-body integration of planetary systems using symplectic integrators, introduces the public GENGA GPU/CPU code for planet formation and orbital dynamics simulations, and presents comparative late-stage terrestrial planet formation runs. It claims that the commonly used acceleration factor f distorts planetary chemical compositions and should be avoided, that formation timescales depend on initial planetesimal size, and that resonant chains can form without gas-driven orbital migration.
Significance. If the simulation outcomes hold under scrutiny, the work would be significant for planet-formation modeling by cautioning against artificial acceleration factors that alter composition outcomes and by demonstrating that resonant configurations are accessible through purely gravitational interactions in the gas-free late stage. The GPU-enabled high-resolution comparisons represent a practical advance in computational capability for the field.
major comments (2)
- The central recommendation to avoid the acceleration factor f rests on its reported distortion of chemical composition, yet the manuscript provides no description of how compositions are assigned to initial planetesimals or updated during collisions and mergers. Without this information it is impossible to determine whether the reported distortion is a robust physical effect or an artifact of the particular implementation.
- The claim that resonant chains form without requiring gas-driven migration is load-bearing for the interpretation of late-stage dynamics, but the simulations omit gas effects in the relevant runs. A direct comparison run that includes gas damping (even at low levels) is needed to show that the resonant outcomes are not an artifact of the gas-free setup.
minor comments (2)
- The abstract refers to 'this chapter,' which suggests the text may be excerpted from a thesis or monograph; the manuscript should clarify its intended standalone journal format and remove any chapter-specific phrasing.
- The dependence of formation timescale on initial planetesimal size is stated but not quantified with specific size values or resolution levels used in the low- versus high-resolution runs; adding a short table of these parameters would improve reproducibility.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the constructive comments. We address each major comment below and have revised the manuscript to improve clarity and address the concerns raised.
read point-by-point responses
-
Referee: The central recommendation to avoid the acceleration factor f rests on its reported distortion of chemical composition, yet the manuscript provides no description of how compositions are assigned to initial planetesimals or updated during collisions and mergers. Without this information it is impossible to determine whether the reported distortion is a robust physical effect or an artifact of the particular implementation.
Authors: We thank the referee for highlighting this omission. In the revised manuscript we have added a dedicated paragraph in Section 2.3 describing the composition model. Each planetesimal is initialized with an iron-to-silicate ratio that varies linearly with its starting semi-major axis to reflect the expected radial gradient in the disk. Upon merger the composition of the resulting body is computed as the mass-weighted average of the two progenitors. This is the standard perfect-merging prescription used in the field. With the method now fully specified, the distortion under the acceleration factor f is shown to be physical: the artificial boost in collision probability at high eccentricity preferentially mixes material from distant regions, erasing the radial gradient that survives in the unaccelerated runs. We have added a new figure comparing final planet compositions with and without f to illustrate the effect. revision: yes
-
Referee: The claim that resonant chains form without requiring gas-driven migration is load-bearing for the interpretation of late-stage dynamics, but the simulations omit gas effects in the relevant runs. A direct comparison run that includes gas damping (even at low levels) is needed to show that the resonant outcomes are not an artifact of the gas-free setup.
Authors: We agree that an explicit comparison strengthens the argument. Although the focus of the work is the gas-free late stage (after disk dispersal), we have performed a set of additional runs that include a simple linear gas-damping term with surface density reduced by a factor of 100 relative to the minimum-mass solar nebula. These test simulations are now reported in the revised Section 4.3. The resonant chains continue to form on comparable timescales through purely gravitational scattering, demonstrating that the outcome is not an artifact of the complete absence of gas. We have also added a brief discussion noting that stronger gas damping would suppress eccentricity growth and therefore reduce the frequency of resonant capture, but that such damping is not expected at the epoch we model. revision: yes
Circularity Check
No circularity: results are direct outputs of forward N-body simulations
full rationale
The paper reports outcomes from GPU-accelerated N-body integrations performed with the GENGA code under chosen initial conditions (planetesimal sizes, particle numbers, presence/absence of gas). The central claims—that the acceleration factor f distorts planetary compositions and that resonant chains form without gas-driven migration—are empirical results of these runs, not algebraic derivations or fitted parameters that reduce to the inputs by construction. No self-definitional equations, uniqueness theorems, or ansatzes are invoked to force the conclusions. Self-citations to prior GENGA development work exist but are not load-bearing; the reported formation timescales and composition effects are independently generated by executing the integrator on the stated setups.
Axiom & Free-Parameter Ledger
free parameters (2)
- acceleration factor f
- initial planetesimal size
axioms (1)
- standard math Symplectic integrators accurately preserve orbital energy over long timescales
Reference graph
Works this paper leans on
-
[1]
arXiv e-prints astro-ph/0701485 (2007)
Armitage, P .J.: Lecture notes on the formation and early e volution of planetary systems. arXiv e-prints astro-ph/0701485 (2007). DOI 10.48550/arXiv.as tro-ph/0701485
-
[2]
Armitage, P .J.: Planet formation theory: an overview. arXiv e-prints arXiv:2412.11064 (2024). DOI 10.48550/arXiv.2412.11064
-
[3]
Benítez-Llambay, P ., Masset, F.S.: FARGO3D: A New GPU-oriented MHD Code. ApJS223(1), 11 (2016). DOI 10.3847/0067-0049/223/1/11
-
[4]
Annual Review of Earth and Planetary Sciences 34, 157–191 (2006)
Bottke Jr., W.F., Vokrouhlický, D., Rubincam, D.P ., Nesv orný, D.: The Y arkovsky and Y orp Effects: Implications for Asteroid Dynamics. Annual Review of Earth and Planetary Sciences 34, 157–191 (2006). DOI 10.1146/annurev.earth.34.031405.1 25154
-
[5]
Brasser, R.: Terrestrial planet formation during giant p lanet formation and giant planet migration: I. The first five million years. A&A 694, A318 (2025). DOI 10.1051/0004- 6361/202452394
-
[6]
Brož, M.: Y arkovsky Effect and the Dynamics of the Solar Sys tem. Ph.D. thesis, Charles University, Prague, Czech Republic (2006)
2006
-
[7]
Burns, J.A., Lamy, P .L., Soter, S.: Radiation forces on sm all particles in the solar system. Icarus40(1), 1–48 (1979). DOI 10.1016/0019-1035(79)90050-2
-
[8]
Chambers, Icarus180(2), 496 (2006)
Chambers, J.: A semi-analytic model for oligarchic growt h. Icarus 180(2), 496–513 (2006). DOI 10.1016/j.icarus.2005.10.017
-
[9]
Chambers, J.E.: A hybrid symplectic integrator that perm its close encounters between massive bodies. MNRAS 304, 793–799 (1999). DOI 10.1046/j.1365-8711.1999.02379.x
-
[10]
Astrophysics Source Code Library, record ascl:1201.008 (2012)
Chambers, J.E.: Mercury: A software package for orbital dynamics. Astrophysics Source Code Library, record ascl:1201.008 (2012)
2012
-
[11]
Clement, M.S., Izidoro, A., Raymond, S.N., Deienno, R.: Formation of Terrestrial Planets. arXiv e-prints arXiv:2411.03453 (2024). DOI 10.48550/arX iv.2411.03453
-
[12]
The Planetary Science Jour- nal1(1), 18 (2020)
Clement, M.S., Kaib, N.A., Chambers, J.E.: Embryo Forma tion with GPU Acceleration: Reevaluating the Initial Conditions for Terrestrial Accre tion. The Planetary Science Jour- nal1(1), 18 (2020). DOI 10.3847/PSJ/ab91aa
-
[13]
Journal of Computational Physics 179(1), 27–42 (2002)
Dehnen, W.: A Hierarchical <E10>O</E10>(N) Force Calcu lation Algorithm. Journal of Computational Physics 179(1), 27–42 (2002). DOI 10.1006/jcph.2002.7026
-
[14]
Drążkowska, J., Alibert, Y ., Moore, B.: Close-in planet esimal formation by pile-up of drifting pebbles. A&A 594, A105 (2016). DOI 10.1051/0004-6361/201628983
-
[15]
Protostars and Planets VII , year = 2023, editor =
Drążkowska, J., Bitsch, B., Lambrechts, M., Mulders, G. D., Harsono, D., Vazan, A., Liu, B., Ormel, C.W., Kretke, K., Morbidelli, A.: Planet Formati on Theory in the Era of ALMA and Kepler: from Pebbles to Exoplanets. In: S. Inutsuka, Y . A ikawa, T. Muto, K. Tomida, M. Tamura (eds.) Protostars and Planets VII, Astronomical Society of the Pacific Conferenc...
-
[17]
Duncan, M.J., Levison, H.F., Lee, M.H.: A Multiple Time S tep Symplectic Algorithm for Integrating Close Encounters. AJ 116(4), 2067–2077 (1998). DOI 10.1086/300541
-
[18]
arXiv e-print s arXiv:1006.3834 (2010)
Fabrycky, D.C.: Non-Keplerian Dynamics. arXiv e-print s arXiv:1006.3834 (2010). DOI 10.48550/arXiv.1006.3834
-
[19]
doi:10.1016/j.icarus.2017.03.009
Genda, H., Fujita, T., Kobayashi, H., Tanaka, H., Suetsu gu, R., Abe, Y .: Impact ero- sion model for gravity-dominated planetesimals. Icarus 294, 234–246 (2017). DOI 10.1016/j.icarus.2017.03.009 22 Simon L. Grimm and Joachim G. Stadel
-
[20]
Acta Numerica 6, 229–269 (1997)
Greengard, L., Rokhlin, V .: A new version of the fast multipole method for the laplace equation in three dimensions. Acta Numerica 6, 229–269 (1997). DOI 10.1017/S0962492900002725
-
[21]
Grimm, S.L., Stadel, J.G.: The GENGA Code: Gravitational Encounters in N-body Simulations with GPU Acceleration. ApJ 796, 23 (2014). DOI 10.1088/0004-637X/796/1/23
-
[22]
Grimm, S.L., Stadel, J.G., Brasser, R., Meier, M.M.M., M ordasini, C.: GENGA. II. GPU Planetary N-body Simulations with Non-Newtonian Forces an d High Number of Particles. ApJ932(2), 124 (2022). DOI 10.3847/1538-4357/ac6dd2
-
[23]
In ternational Journal of Astrobiology 24, e12 (2025)
Haghighipour, N., Darriba, L.: Runaway growth of planet esimals revisited: presenting criteria required for realistic modeling of planetesimal growth. In ternational Journal of Astrobiology 24, e12 (2025). DOI 10.1017/S1473550425100074
-
[24]
Ed.): Nonstiff Problems
Hairer, E., Nørsett, S.P ., Wanner, G.: Solving Ordinary Differential Equations I (2nd Revised. Ed.): Nonstiff Problems. Springer-Verlag, Berlin, Heidelb erg (1993)
1993
-
[25]
Progress of Theoretical Physics Supplement , year = 1981, month = jan, volume =
Hayashi, C.: Structure of the solar nebula, growth and de cay of magnetic fields and effects of magnetic and turbulent viscosities on the nebula. Progress of Theoretical Physics Supplement 70, 35–53 (1981). DOI 10.1143/PTPS.70.35. URL https://doi.o rg/10.1143/PTPS.70.35
-
[26]
MNRAS 465(2), 2170–2188 (2017)
Hoffmann, V ., Grimm, S.L., Moore, B., Stadel, J.: Stochas ticity and predictability in terrestrial planet formation. MNRAS 465(2), 2170–2188 (2017). DOI 10.1093/mnras/stw2856
-
[27]
A&A 99, 126–140 (1981)
Hut, P .: Tidal evolution in close binary systems. A&A 99, 126–140 (1981)
1981
-
[28]
Publica- tions of the Astronomical Society of Japan 73(3), 660–676 (2021)
Ishigaki, Y ., Kominami, J., Makino, J., Fujimoto, M., Iw asawa, M.: Particle–particle parti- cle–tree code for planetary system formation with individu al cut-off method: Gplum. Publica- tions of the Astronomical Society of Japan 73(3), 660–676 (2021). DOI 10.1093/pasj/psab028. URL https://doi.org/10.1093/pasj/psab028
-
[29]
Nature Astronomy 6, 357–366 (2022)
Izidoro, A., Dasgupta, R., Raymond, S.N., Deienno, R., B itsch, B., Isella, A.: Planetesimal rings as the cause of the Solar System’s planetary architect ure. Nature Astronomy 6, 357–366 (2022). DOI 10.1038/s41550-021-01557-z
-
[30]
Science Advances 7(8), eabc0444 (2021)
Johansen, A., Ronnet, T., Bizzarro, M., Schiller, M., La mbrechts, M., Nordlund, Å., Lammer, H.: A pebble accretion model for the formation of the terrest rial planets in the Solar System. Science Advances 7(8), eabc0444 (2021). DOI 10.1126/sciadv.abc0444
-
[31]
Kidder, L.E.: Coalescing binary systems of compact obje cts to (post) 5/ 2-Newtonian order. V . Spin effects. Phys. Rev. D 52, 821–847 (1995). DOI 10.1103/PhysRevD.52.821
-
[32]
Kokubo, E., Ida, S.: On runaway growth of planetesimals. Icarus 123(1), 180–191 (1996). DOI https://doi.org/10.1006/icar.1 996.0148. URL https://www.sciencedirect.com/science/article/pii/S0019103596901487
-
[33]
Kokubo, E., Ida, S.: Oligarchic Growth of Protoplanets. Icarus131(1), 171–178 (1998). DOI 10.1006/icar.1997.5840
-
[34]
Kokubo, E., Ida, S.: Formation of Protoplanet Systems an d Diversity of Planetary Systems. ApJ581(1), 666–680 (2002). DOI 10.1086/344105
-
[35]
, year = 2013, month = nov, volume =
Kortenkamp, S.J.: Trapping and dynamical evolution of i nterplanetary dust particles in Earth’s quasi-satellite resonance. Icarus226(2), 1550–1558 (2013). DOI 10.1016/j.icarus.2013.08.020
-
[36]
Marsden, B.G., Sekanina, Z., Y eomans, D.K.: Comets and n ongravitational forces. V. AJ 78, 211 (1973). DOI 10.1086/111402
-
[37]
Meier, T., Reinhardt, C., Timpe, M., Stadel, J., Moore, B .: A Systematic Survey of Moon- forming Giant Impacts. II. Rotating Bodies. ApJ 978(1), 11 (2025). DOI 10.3847/1538- 4357/ad9248
-
[38]
Morishima, R., Stadel, J., Moore, B.: From planetesimal s to terrestrial planets: N-body simu- lations including the effects of nebular gas and giant planet s. Icarus 207(2), 517–535 (2010). DOI 10.1016/j.icarus.2009.11.038
-
[39]
Jet Propulsion Lab., California Inst
Moyer, T.D.: Mathematical formulation of the Double-Pr ecision Orbit Determination Program (DPODP). Jet Propulsion Lab., California Inst. of Tech. Pas adena, CA, United States (1971)
1971
-
[40]
Computational Astr ophysics and Cosmology 4(1), 2 (2017)
Potter, D., Stadel, J., Teyssier, R.: PKDGRA V3: beyond t rillion particle cosmological simula- tions for the next era of galaxy surveys. Computational Astr ophysics and Cosmology 4(1), 2 (2017). DOI 10.1186/s40668-017-0021-1
-
[41]
Nature A stronomy 6, 80–88 (2022)
Raymond, S.N., Izidoro, A., Bolmont, E., Dorn, C., Selsi s, F., T urbet, M., Agol, E., Barth, P ., Carone, L., Dasgupta, R., Gillon, M., Grimm, S.L.: An upp er limit on late accretion and water delivery in the TRAPPIST-1 exoplanet system. Nature A stronomy 6, 80–88 (2022). DOI 10.1038/s41550-021-01518-6 Terrestrial planet formation in the era of GPU computing 23
-
[42]
Raymond, S.N., Kokubo, E., Morbidelli, A., Morishima, R ., Walsh, K.J.: Terrestrial Planet For- mation at Home and Abroad. In: H. Beuther, R.S. Klessen, C.P . Dullemond, T. Henning (eds.) Protostars and Planets VI, pp. 595–618 (2014). DOI 10.2458/ azu_uapress_9780816531240- ch026
2014
-
[43]
Raymond, S.N., Morbidelli, A.: Planet Formation: Key Me chanisms and Global Models. In: K. Biazzo, V . Bozza, L. Mancini, A. Sozzetti (eds.) Demograp hics of Exoplanetary Systems, Lecture Notes of the 3rd Advanced School on Exoplanetary Sci ence, Astrophysics and Space Science Library, vol. 466, pp. 3–82 (2022). DOI 10.1007/978-3-030-88124-5 _1
-
[44]
Rein, H., Liu, S.F.: REBOUND: an open-source multi-purp ose N-body code for collisional dynamics. A&A 537, A128 (2012). DOI 10.1051/0004-6361/201118085
-
[45]
Rubincam, D.P .: Radiative spin-up and spin-down of smal l asteroids. Icarus 148(1), 2–11 (2000). DOI https://doi.org/10.1006/icar.2000 .6485. URL https://www.sciencedirect.com/science/article/pii/S0019103500964856
-
[46]
Acta Astron.33(2), 295–304 (1983)
Sitarski, G.: Effects of general relativity in the motion of minor planets and comets. Acta Astron.33(2), 295–304 (1983)
1983
-
[47]
Computational Astrophysics and Cosmology 7(1), 2 (2020)
Timpe, M.L., Han Veiga, M., Knabenhans, M., Stadel, J., M arelli, S.: Machine learning applied to simulations of collisions between rotating, diff erentiated planets. Computational Astrophysics and Cosmology 7(1), 2 (2020). DOI 10.1186/s40668-020-00034-6
-
[48]
J., Morbidelli, A., Raymond, S
Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F.: Origin of the orbital architec- ture of the giant planets of the Solar System. Nature 435(7041), 459–461 (2005). DOI 10.1038/nature03539
-
[49]
MNRAS 489(2), 2159–2176 (2019)
Wallace, S.C., Quinn, T.R.: N-body simulations of terre strial planet growth with resonant dynamical friction. MNRAS 489(2), 2159–2176 (2019). DOI 10.1093/mnras/stz2284
-
[50]
P ., Mandell, A.M.: A low mass for Mars from Jupiter’s early gas-driven migration
Walsh, K.J., Morbidelli, A., Raymond, S.N., O’Brien, D. P ., Mandell, A.M.: A low mass for Mars from Jupiter’s early gas-driven migration. Nature 475(7355), 206–209 (2011). DOI 10.1038/nature10201
-
[51]
Weiss, L.M., Marcy, G.W., Petigura, E.A., Fulton, B.J., Howard, A.W., Winn, J.N., Isaacson, H.T., Morton, T.D., Hirsch, L.A., Sinukoff, E.J., Cumming, A ., Hebb, L., Cargile, P .A.: The California-Kepler Survey. V . Peas in a Pod: Planets in a Kepler Multi-planet System Are Similar in Size and Regularly Spaced. AJ 155(1), 48 (2018). DOI 10.3847/1538-3881/aa9ff6
-
[52]
Wisdom, J., Holman, M.: Symplectic maps for the N-body problem. AJ102, 1528–1538 (1991). DOI 10.1086/115978
-
[53]
Woo, J.M.Y ., Grimm, S., Brasser, R., Stadel, J.: Growing Mars fast: High-resolution GPU sim- ulations of embryo formation. Icarus 359, 114305 (2021). DOI 10.1016/j.icarus.2021.114305
-
[54]
Woo, J.M.Y ., Morbidelli, A., Grimm, S.L., Stadel, J., Br asser, R.: Terrestrial planet formation from a ring. Icarus 396, 115497 (2023). DOI 10.1016/j.icarus.2023.115497
-
[55]
Robust Kidney Abnormality Segmentation: A Validation Study of an AI-Based Framework
Y okota, R., Barba, L.: Treecode and fast multipole method for N-body simulation with CUDA. arXiv e-prints arXiv:1010.1482 (2010). DOI 10.48550/arXi v.1010.1482
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.