Hilda Population Model

Updated 3 December 2025

The Hilda population model is defined by asteroids locked in a stable 3:2 resonance with Jupiter, characterized by precise orbital elements and bounded libration behavior.
It employs bias-corrected surveys and synthetic proper elements to decompose the population into collisional families and background, revealing distinct magnitude and size-frequency distributions.
The model integrates dynamical simulations, spectral analysis, and collisional evolution studies to elucidate planetary migration impacts and test hypotheses such as free-floating planet flybys.

The Hilda population model characterizes the distribution, dynamics, history, and physical properties of asteroids in the 3:2 mean-motion resonance with Jupiter. This population exhibits distinctive dynamical confinement, collisional evolution, compositional diversity, and resonance-induced structure—serving as a benchmark for understanding Solar System dynamical instability, outer planet migration scenarios, and the interplay between collisional physics and resonance dynamics.

1. Dynamical and Resonant Definition

The Hilda group is defined by stable libration in the 3:2 mean-motion resonance with Jupiter, centered near $a_{\mathrm{J3/2}} \simeq 3.98\,\mathrm{AU}$ , with long-term stability for proper elements $a_p\in[3.95,4.05]$ AU, $e_p\lesssim0.32$ , $I_p\lesssim20^\circ$ , and resonant angle $\sigma=3\lambda_J-2\lambda-\varpi$ undergoing bounded libration. Hamiltonian models using the planar Circular and Elliptic Restricted Three-Body Problems (CRTBP/ERTBP) demonstrate the existence of families of elliptic periodic orbits surrounded by KAM tori, producing robust dynamical islands that confine the Hilda orbits over $>10^8$ yr. Frequency analysis identifies typical dominant spectral lines for Hildas at $\omega_1\in[0.50,0.51]$ , $\omega_2\in[0.45,0.47]$ (in units of Jupiter's synodic frequency), with the combination $(\omega_1,\,\omega_2)$ providing a sharper membership criterion than classical (a,e,i)-boxes or two-body orbital elements. Quasi-periodic approximations to Hilda orbital evolution require combinational frequencies capturing both secular and libration phenomena, with libration frequencies $\nu\approx0.059\pm0.001\,\mathrm{kyr}^{-1}$ and secular frequencies $a_p\in[3.95,4.05]$ 0 for $a_p\in[3.95,4.05]$ 1 and $a_p\in[3.95,4.05]$ 2 for $a_p\in[3.95,4.05]$ 3 (Jorba et al., 2024, Rosaev, 2023).

2. Population Structure: Families, Background, and Magnitude Distribution

Using a decade of bias-corrected observational data, the Hilda population is decomposed into background and collisional family components by computing synthetic proper elements and hierarchical clustering in $a_p\in[3.95,4.05]$ 4-space. Three major families dominate: Hilda (153), Schubart (1911), and Potomac—together exceeding 60% of the population for $a_p\in[3.95,4.05]$ 5. The cumulative magnitude distribution $a_p\in[3.95,4.05]$ 6 of the background is described by piecewise-linear slopes $a_p\in[3.95,4.05]$ 7, with a mean $a_p\in[3.95,4.05]$ 8 for $a_p\in[3.95,4.05]$ 9, significantly shallower than the Jupiter Trojans ( $e_p\lesssim0.32$ 0). Family slopes (mid-range) are steeper: $e_p\lesssim0.32$ 1, $e_p\lesssim0.32$ 2, $e_p\lesssim0.32$ 3. For $e_p\lesssim0.32$ 4, background and family components each contain $e_p\lesssim0.32$ 5 objects, with the total sample $e_p\lesssim0.32$ 6. As the sample extends to fainter magnitudes, the family fraction rises (Vokrouhlický et al., 6 Mar 2025).

3. Size-Frequency and Albedo Distributions

The size distribution of Hildas between $e_p\lesssim0.32$ 7 and $e_p\lesssim0.32$ 8 km is characterized by a single-slope power law. Subaru Hyper Suprime-Cam data yield $e_p\lesssim0.32$ 9 (differential in $I_p\lesssim20^\circ$ 0), giving a cumulative slope $I_p\lesssim20^\circ$ 1 in diameter, with $I_p\lesssim20^\circ$ 2 for $I_p\lesssim20^\circ$ 3 km. This closely matches the Jupiter Trojan population for $I_p\lesssim20^\circ$ 4 km and is distinct from the "wavy" structure of the main-belt asteroids, indicating different formation and early collisional histories. Infrared and Spitzer data show a mean geometric albedo $I_p\lesssim20^\circ$ 5 for $I_p\lesssim20^\circ$ 6 km and $I_p\lesssim20^\circ$ 7 for $I_p\lesssim20^\circ$ 8 km, with a significant anti-correlation between size and albedo. The range $I_p\lesssim20^\circ$ 9 among small Hildas encompasses the C-, D-, and X-type classes and a high-albedo tail attributed to outer solar system contamination (Terai et al., 2018, Ryan et al., 2011).

4. Collisional and Dynamical Evolution

The large-end Hilda SFD ( $\sigma=3\lambda_J-2\lambda-\varpi$ 0 km, slope $\sigma=3\lambda_J-2\lambda-\varpi$ 1) is primordial; collisional models and 4 Gyr Monte Carlo simulations reveal very low disruption rates and only limited evolution for multikilometer bodies. For $\sigma=3\lambda_J-2\lambda-\varpi$ 2 km, SFDs steepen (slope $\sigma=3\lambda_J-2\lambda-\varpi$ 3) due to catastrophic disruption of a small number of larger bodies, but detailed structures depend on initial assumptions about the small-end slope. The impactors responsible for the largest craters on (334) Chicago ( $\sigma=3\lambda_J-2\lambda-\varpi$ 4 km) are themselves only $\sigma=3\lambda_J-2\lambda-\varpi$ 5– $\sigma=3\lambda_J-2\lambda-\varpi$ 6 km in size, producing maximum craters of $\sigma=3\lambda_J-2\lambda-\varpi$ 7– $\sigma=3\lambda_J-2\lambda-\varpi$ 8 km. Subcatastrophic impact timescales for quasi-Hilda objects are $\sigma=3\lambda_J-2\lambda-\varpi$ 9 yr, much longer than their typical dynamical lifetimes, implying collisional activity is not the mechanism for observed cometary activity in these bodies. The current collisional probability ( $>10^8$ 0, $>10^8$ 1) is insufficient to significantly alter the SFD above a few kilometers; therefore, the major population structure is set by primordial implantation and early events (Zain et al., 18 Jan 2025).

5. Color Bimodality and Compositional Interpretations

Multiband photometry and Sloan Digital Sky Survey (SDSS) data reveal a robust bimodality in the visible spectral slope among Hildas, with two Gaussian subpopulations: a less-red (LR) peak at $>10^8$ 2 and a red (R) peak at $>10^8$ 3, in ratio $>10^8$ 4. Collisional families are exclusively LR, explained by volatile loss in parent-body disruption. This bimodality is mirrored in the Jupiter Trojans (almost identical means and ratios), with family fragments also being only LR. The bimodality's invariance under further collisional evolution, and its match with Trojans, supports a common origin in a trans-Neptunian planetesimal reservoir and subsequent migration (Wong et al., 2017).

6. Resonance Amplitude Structure and FFP Flyby Hypothesis

High-precision modeling uncovers an observed "desert" of Hildas with resonant amplitudes $>10^8$ 5 at $>10^8$ 6 and a nearly complete lack of any $>10^8$ 7 orbits across all $>10^8$ 8. Standard migration/capture models reproduce Hilda $>10^8$ 9-distributions but not this unusual amplitude cutoff. Numerical simulations show that a flyby of a free-floating planet (FFP) with $\omega_1\in[0.50,0.51]$ 0, $\omega_1\in[0.50,0.51]$ 1, $\omega_1\in[0.50,0.51]$ 2 can instantaneously shift Jupiter's orbit by $\omega_1\in[0.50,0.51]$ 3 AU, moving the 3:2 resonance by $\omega_1\in[0.50,0.51]$ 4 AU. This projects surviving Hildas across amplitude space, producing the observed $\omega_1\in[0.50,0.51]$ 5– $\omega_1\in[0.50,0.51]$ 6 pattern. The pattern arises independently of the primordial amplitude PDF, persists for a wide range of FFP parameters, and is not replicated by smooth migration alone. The model also accounts for the Trojan L4:L5 asymmetry. The FFP flyby hypothesis predicts minor inclination excitation ( $\omega_1\in[0.50,0.51]$ 7), possible depletions in the Cybele region, and “scars” in the high-inclination main belt. Constraining $\omega_1\in[0.50,0.51]$ 8– $\omega_1\in[0.50,0.51]$ 9 phase-space boundaries in the Hildas and further comparison with Cybele/JT populations can further test this scenario (Li et al., 2024).

7. Integration into Solar System Evolution

The current Hilda model is consistent with the late-stage dynamical reshaping of the Solar System. Major collisional families (e.g., Hilda, Schubart) likely formed during or shortly after the Late Heavy Bombardment (LHB), supported by age estimates $\omega_2\in[0.45,0.47]$ 0 and requirements for rapid Jupiter migration timescales ( $\omega_2\in[0.45,0.47]$ 1– $\omega_2\in[0.45,0.47]$ 2 Myr) to match observed orbital dispersions. The model reproduces the observed "ears" in the $\omega_2\in[0.45,0.47]$ 3 plane of families, the SFD shape, the magnitude and phase-space distributions, and is incompatible with a scenario of high ongoing collisional rates. The synthetic proper-element and family decomposition enables direct quantitative comparison with population synthesis predictions from giant planet migration/instability models (Brož et al., 2011, Vokrouhlický et al., 6 Mar 2025).

This synthesis defines the Hilda population in dynamical, physical, and evolutionary terms, grounded in resonance mapping, collisional modeling, bias-corrected surveys, and comprehensive phase-space treatment, with strong links to planetary migration scenarios and Solar System evolution models.