pith. sign in

arxiv: 2604.19413 · v2 · pith:ZCUKP2QKnew · submitted 2026-04-21 · 🌌 astro-ph.EP · astro-ph.GA· astro-ph.SR

Oort Cloud Ecology -- IV. Exchanging Asteroids

Erwan Hochart , Simon Portegies Zwart This is my paper

Pith reviewed 2026-05-10 01:43 UTC · model grok-4.3

classification 🌌 astro-ph.EP astro-ph.GAastro-ph.SR
keywords asteroid exchangeOort Cloudstar clustersfractal distributionPlummer modelcaptured asteroidsrogue objectstrans-Neptunian objects
0
0 comments X

The pith

Sub-virial fractal star clusters exchange asteroids more readily than virialised Plummer clusters, and both suppress Oort Cloud formation.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper runs numerical simulations of two star cluster models to determine how the birth environment shapes asteroid orbits, captures, and ejections. A sub-virial fractal distribution produces more frequent high-eccentricity and high-inclination orbits, more captured asteroids, and more rogue objects than a virialised Plummer distribution. In the fractal case, objects that resemble extreme trans-Neptunian objects and Sednoids sit in the same phase-space regions occupied by captured material, while the Plummer case yields mostly native asteroids in those same regions. Neither model assembles a substantial Oort Cloud, showing that dense cluster conditions hinder its formation regardless of the precise dynamical temperature.

Core claim

The sub-virial fractal cluster exhibits richer dynamics, with asteroids and planets more frequently acquiring high eccentricities and inclinations, along with a larger fraction of captured and rogue objects. Additionally, this cluster configuration has its extreme trans-Neptunian object and Sednoid analogues occupy regions of phase-space in semi-major axis, eccentricity and inclination commonly frequented by captured asteroids. Although the virialised Plummer model can produce such objects, by being less dynamically active, the vast majority of asteroids occupying these regions are native rather than captured. Lastly, neither model efficiently form an Oort Cloud, indicating that Oort Cloud

What carries the argument

N-body simulations of two 150-star clusters (sub-virial fractal versus virial Plummer), each star initially carrying 500 asteroids and, for stars below 2 solar masses, 1–8 planets, evolved for 30 Myr.

If this is right

  • Asteroids and planets acquire high eccentricities and inclinations more often in the sub-virial fractal cluster.
  • A larger fraction of asteroids become captured or rogue in the fractal model than in the Plummer model.
  • Extreme trans-Neptunian object and Sednoid analogues in the fractal cluster occupy phase-space regions typical of captured asteroids.
  • The Plummer model produces mostly native asteroids in those same extreme regions.
  • Neither cluster model assembles an Oort Cloud efficiently.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The Sun's birth cluster was probably closer to the more quiescent Plummer type, since an efficient Oort Cloud exists today.
  • Asteroid exchange in dense clusters could mix material between different stellar systems and alter planetary compositions.
  • Future surveys that classify Sednoids as captured or native could distinguish between fractal and smooth birth environments.
  • Extending the integrations past 30 Myr would test whether Oort Cloud assembly can occur later once clusters disperse.

Load-bearing premise

The chosen initial conditions, including a 0.5 pc virial radius, exactly 500 asteroids per star, and integration to 30 Myr, represent typical star-forming regions and capture the dominant processes for asteroid exchange and Oort Cloud assembly.

What would settle it

A census of the fraction of captured versus native extreme trans-Neptunian objects in the solar system or in other young systems that matches the proportions produced by one cluster model but not the other.

Figures

Figures reproduced from arXiv: 2604.19413 by Erwan Hochart, Simon Portegies Zwart.

Figure 1
Figure 1. Figure 1: Initial asteroid density profile for a randomly selected system. The black curve shows the theoretical distribution (where Σ ∝ r −3/2 ). The blue curve shows the asteroid population at that annulus. Shaded red regions show locations devoid of asteroids due to the assumption that planets carve out a region ±3 Hill radii. This system has three plan￾ets initially. Nemesis is a hybrid N-body integrator embedde… view at source ↗
Figure 2
Figure 2. Figure 2: Evolution of cluster density in time. Shaded regions represent the interquartile range. clumps (i.e, swarms of asteroids) and isolated massive objects are added directly to the parent code. 2.4.3. Free Parameters In Nemesis, the user is free to choose the integrator. In this study, children are evolved with the symplectic code Huayno (Pelupessy et al. 2012). The optimised-Kepler (‘OK’) scheme of Huayno is … view at source ↗
Figure 3
Figure 3. Figure 3: Heatmap showing the proportion of captured asteroids in a-i (top) and a-e (bottom) space after 30 Myr. Black scatter points denote the locations of extreme trans-Neptunian objects, while white points represent those of Sednoids. Both axes are sliced into ten equal-sized bins. Regions with no data are shown in white. Tiles occupied by fewer than 100 asteroids per run have their average occupancy population … view at source ↗
Figure 4
Figure 4. Figure 4: Proportion of asteroid fates binned in stellar spectral types for both models. Red is NGC 1333p, blue is NGC 1333f. From lighter to darker tones: Native fraction (i.e., bound to the same host at the initial and final time), captured fraction, other (removed due to tight orbits), and ejected. Values exceed unity (black dashed line) due to the cap￾tured population. Spectral type bins follow the classificatio… view at source ↗
Figure 7
Figure 7. Figure 7: The fraction of captured population for various asteroid families. Hills-Oort objects have Porb ≲ 5.2 Myr and rp ≥ 60, long-period comets have e > 0.2, rp < 10(M∗/ M⊙) 1/3 au and Porb > 200 yr, short-period comets have the same e and rp criteria but Porb < 200 yr. Scatter point colours indicate the mean number of asteroids per spectral classification bin. Bars represent the 25% and 75% quartile range. Case… view at source ↗
Figure 6
Figure 6. Figure 6: Fraction of stars hosting fcap captured objects for model NGC 1333f. Data only considers stars that have at least Nast > 10 at the end time. Solid lines consider all asteroids bound to a star. Dashed lines con￾sider only asteroids with periastron rp < 10 au. We ignore the relevant curve for M < 1 M⊙ to reduce clutter, and has results nearly identical to the global curve. > 200 au remain uncommon in young c… view at source ↗
Figure 8
Figure 8. Figure 8: Two-point correlation function of initial distance between an asteroid’s original host star and its eventual captured star. this difference is due to the inner edge of the primordial debris disk considered here being nearer the host star. Consequently, a greater proportion of the native asteroid population is sheltered from the external environment. Systems that never experience close encounters (nearest s… view at source ↗
read the original abstract

Aims. Investigate the influence of cluster environments on asteroids, with special attention towards captured material. Methods. Using numerical methods, a sub-virial fractally distributed star-forming region and a virialised Plummer distributed star-forming region are simulated. Both models are initialised with a virial radius of 0.5pc and 150 stars. Stellar populations and their corresponding planetary systems are identical between cluster models. Stars initially host 500 asteroids and those with mass M_* <= 2.0 MSun are also orbited by 1 - 8 planets. Clusters are integrated until 30 Myr. Results. The sub-virial fractal cluster exhibits richer dynamics, with asteroids and planets more frequently acquiring high eccentricities and inclinations, along with a larger fraction of captured and rogue objects. Additionally, this cluster configuration has its extreme trans-Neptunian object and Sednoid analogues occupy regions of phase-space in semi-major axis, eccentricity and inclination commonly frequented by captured asteroids. Although the virialised Plummer model can produce such objects, by being less dynamically active, the vast majority of asteroids occupying these regions are native rather than captured. Lastly, neither model efficiently form an Oort Cloud, indicating that Oort Cloud assembly is strongly suppressed in both dynamically hot and more quiescent cluster

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript reports N-body simulations comparing asteroid and planet dynamics in two star-forming cluster models with identical stellar masses, planetary systems (1-8 planets for M_* <= 2 Msun), and asteroid populations (500 per star): a sub-virial fractal distribution versus a virialised Plummer sphere, both with 150 stars and 0.5 pc virial radius, integrated to 30 Myr. The central claim is that the fractal model produces richer dynamics, with asteroids and planets more often reaching high eccentricities and inclinations, higher fractions of captured and rogue objects, and extreme TNO/Sednoid analogues occupying phase-space regions (in a, e, i) commonly populated by captured asteroids; the Plummer model produces such objects but mostly as native rather than captured. Neither model efficiently assembles an Oort Cloud.

Significance. If the numerical results hold, the work demonstrates that initial cluster spatial/kinematic structure modulates asteroid exchange, capture rates, and the origin of detached trans-Neptunian objects, with implications for Solar System formation and exoplanetary small-body populations. The controlled comparison that holds stellar and planetary components fixed while varying only the cluster structure is a clear strength, as is the purely forward simulation from Newtonian gravity with no fitted parameters.

major comments (1)
  1. [Methods] Methods section: the description of the N-body setup does not specify the integrator, time-stepping criteria, softening length, close-encounter handling, or any energy-conservation or convergence tests. These details are load-bearing for the central claim that the fractal model exhibits richer dynamics and higher capture fractions, as the reported differences could be sensitive to numerical choices.
minor comments (2)
  1. [Abstract] Abstract: the statement that 'neither model efficiently form an Oort Cloud' would benefit from a brief quantitative definition of 'efficient' (e.g., fraction of asteroids reaching a > 1000 au and q > 40 au by 30 Myr) to allow direct comparison with longer-term studies.
  2. [Results] Results: the phase-space occupation claim for extreme TNO/Sednoid analogues would be strengthened by a quantitative metric (e.g., overlap fraction or Kolmogorov-Smirnov test between captured and native populations) rather than qualitative description.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for recommending minor revision. The single major comment is addressed below; we have revised the manuscript to incorporate the requested information.

read point-by-point responses

  1. Referee: [Methods] Methods section: the description of the N-body setup does not specify the integrator, time-stepping criteria, softening length, close-encounter handling, or any energy-conservation or convergence tests. These details are load-bearing for the central claim that the fractal model exhibits richer dynamics and higher capture fractions, as the reported differences could be sensitive to numerical choices.

    Authors: We agree that these numerical details were insufficiently described in the original submission and are important for assessing the robustness of the reported differences. In the revised manuscript we have expanded the Methods section to specify the integrator, adaptive time-stepping criteria, softening length, close-encounter treatment, and the results of energy-conservation and convergence tests. These tests confirm that the higher fractions of captured and rogue objects, as well as the elevated eccentricities and inclinations, in the sub-virial fractal model relative to the Plummer model are insensitive to the numerical parameters within the ranges explored. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper reports results from direct N-body integrations of two star-cluster models (sub-virial fractal vs. virialised Plummer) that share identical stellar masses, planetary systems, and asteroid populations. All reported differences in eccentricity/inclination distributions, capture/rogue fractions, and phase-space occupation of extreme TNO/Sednoid analogues emerge from the Newtonian evolution of the chosen initial conditions (0.5 pc virial radius, 150 stars, 500 asteroids per star, 1–8 planets for M_* ≤ 2 M_⊙, 30 Myr integration). No parameters are fitted to data, no equations are defined in terms of their own outputs, and no load-bearing claims reduce to self-citation or ansatz. The central comparison is therefore falsifiable by re-running the simulations with the stated initial conditions.

Axiom & Free-Parameter Ledger

5 free parameters · 3 axioms · 0 invented entities

The central results rest on a large set of hand-chosen initial conditions and standard assumptions about gravitational dynamics; no new physical entities are introduced.

free parameters (5)
  • virial radius
    Set to 0.5 pc for both cluster models as the initial scale.
  • number of stars
    Fixed at 150 stars in both models.
  • asteroids per star
    Each star starts with exactly 500 asteroids.
  • integration time
    Clusters evolved to 30 Myr.
  • planet hosting rule
    Stars with M_* <= 2.0 M_Sun host 1-8 planets.
axioms (3)
  • standard math Newtonian gravity governs all interactions between stars, planets, and asteroids.
    Implicit in all N-body integrations described.
  • domain assumption Initial asteroid and planet populations are identical across both cluster models.
    Stated explicitly to enable direct comparison.
  • ad hoc to paper No gas or other non-gravitational forces are included.
    Omitted from the described numerical setup.

pith-pipeline@v0.9.0 · 5537 in / 1701 out tokens · 38044 ms · 2026-05-10T01:43:06.838061+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.