Analytical Solutions for Planet-Scattering Small Bodies

TL;DR

This work derives an analytical framework for the random scattering of planet-crossing, massless bodies by a planet within Öpik's close-encounter theory, recasting the problem as a two-DOF system governed by the relative-velocity vector ${\bm{U}}_{\infty}$. By averaging over encounter geometries, the authors obtain drift and diffusion coefficients for the normalized energy $x$ and solve the resulting Fokker–Planck equation to yield a closed-form, universal description of the time evolution of the particle distribution, including a characteristic scattering timescale $t_S$ and a typical ejection speed $v_{\text{eje}}$. The framework aligns with full numerical integrations in the regimes where crossing orbits dominate and provides rapid, scalable predictions for Solar System and exoplanetary systems, with broad applicability to Kuiper Belt, Oort Cloud, debris disks, interstellar objects, and free-floating planets. Limitations include single-planet assumptions and neglect of resonant sticking and collisions, but extensions to multi-planet networks are outlined for future work.

Abstract

Gravitational scattering of small bodies (planetesimals) by a planet remains a fundamental problem in celestial mechanics. It is traditionally modeled within the circular restricted three-body problem (CR3BP), where individual particle trajectories are obtained via numerical integrations. Here, we use {Ö}pik's close-encounter framework to study the random walk of the orbital energy $x$ for an ensemble of test particles on planet-crossing orbits. We show that the evolution of each particle's orbital elements $(a, e, i)$ is fully encapsulated by the 3D rotation of the relative velocity vector $\bm{U}_\infty$, whose magnitude remains constant. Consequently, the system can be reduced to two degrees of freedom. By averaging over all possible flyby geometries, we derive explicit expressions for the drift and diffusion coefficients of $x$. We then solve the resulting Fokker--Planck equation to obtain a closed-form solution for the time evolution of the particle distribution. A characteristic scattering timescale naturally emerges, scaling as $(P_{p}/M_{p}^{2})/500$, where $P_{p}$ is the planet's orbital period and $M_{p}$ its mass ratio to the central star. The typical ejection speed of small bodies by a planet is estimated to be $3 v_p M_{p}^{1/3}$, where $v_p$ is the planet's orbital speed. Our analytical solution constitutes a universal law applicable to both the Solar System and exoplanetary systems, providing a computationally efficient alternative to costly $N$-body simulations for studying the orbital distributions and ejection of planetesimals and planets (e.g., Kuiper Belt, Oort Cloud, debris disks, interstellar objects, and free-floating planets).

Figures (14)

• Figure 1: Schematic illustration of the 3D scattering geometry in a planet-centered co-rotating frame. The $z$-axis (${ \bm{\hat{\mathbf{k}}} }$) points out of the orbital plane, while $x$-axis (${ \bm{\hat{\mathbf{i}}} }$) and $y$-axis (${ \bm{\hat{\mathbf{j}}} }$) are aligned respectively with the instantaneous radial and transverse directions of the planet at the time of encounter. In this coordinate system, $\bm{r}_p$ lies entirely along the $x$-axis, and $\bm{v}_p$ (the planet's velocity, solid black arrow) lies along the $y$-axis. The small body's heliocentric velocity vector at the intersection $\bm{v_\text{cross}}$ (solid red arrow) is the vector sum of the relative velocity vector $\bm{v}_\infty$ (solid blue arrow) and $\bm{v}_p$. The vector $\bm{U}_\infty = \bm{v}_\infty/v_p$ is the normalized relative velocity. The two angles $\theta$ (the polar angle from $\bm{v}_p$) and $\phi$ (the azimuth angle of $\bm{v}_\infty$ in the $x$-$z$ plane) fully specify the direction of the relative velocity vector.
• Figure 2: Left: Constant $\mathcal{T}$ curves on the dimensionless $\mathcalo{q}$--$\mathcal{A}$ (perihelion distance--semimajor axis) space, assuming $i=0\degree$ (Eq. \ref{['eq:tisserand-a-q']}). In the large-$\mathcal{A}$ limit, $\mathcalo{q}$ asymptotically approaches the value $(\mathcal{T}/\sqrt{8})^2$ (vertical gray dashed lines, Eq. \ref{['eq:tisserand-a-infinity']}). For $\mathcal{T} < \sqrt{8}$, $\lim_{\mathcal{A} \rightarrow \infty}{\mathcalo{q}} < 1$ and the particle can never decouple from the planet in the co-planar case (blue). For $\mathcal{T} > \sqrt{8}$, $\lim_{\mathcal{A} \rightarrow \infty}{\mathcalo{q}} > 1$ and spontaneous decoupling when the particle reaches large $\mathcal{A}$ is guaranteed (orange). The loosely-coupled case and the closely-coupled cases defined in text are separated by the critical value of $\mathcal{T}=\sqrt{8}$ (red), which has $\mathcalo{q} \rightarrow 1$ as $\mathcal{A} \rightarrow \infty$. Eq. \ref{['eq:tisserand-a-q']} is partially undefined near the planet for $\mathcal{T}>3$ (grey). Right: $\mathcal{A}$ as a function of $\theta$ for various $U_\infty$ (Eq. \ref{['eq:def-A']}). The horizontal dashed line marks $\mathcal{A} = 1$. In a loosely-coupled case (golden curves), there is a maximum $\mathcal{A}$ the particle can reach in the planet-crossing state. For the boundary case $U_\infty = \sqrt{2}-1$ (orange), $\theta \rightarrow 0\degree$ corresponds to $\mathcal{A} \rightarrow \infty$. In closely-coupled cases (blue), the object will escape the system once $\theta < \theta_\text{loss}$, which is defined in Eq. \ref{['eq:def-thetaloss']}. A larger $U_\infty$ corresponds to a wider loss cone (blue line segments).
• Figure 3: Evolution of a particle with $U_\infty = 1.0$ ($\mathcal{A}_0$=1.2, $e_0$=0.61, $i_0$=47.8$^\circ$, closely-coupled case) in orbital elements ($\mathcal{A}$, $e$, $i$, left panels) and relative velocity vectors ($\theta$, $\phi$, $U_\infty$, right panels) scattered by a $m_p = 1\times10^{-4} m_\odot$ planet on a $a_p = 1$ au circular orbit. The gray regions mark $\mathcalo{q} > 1$ (red curves), corresponding to where the Öpik framework fails. The blue region in the top-right panel represents the loss cone $\theta_\text{loss}$ (Eq. \ref{['eq:def-thetaloss']}).
• Figure 4: Same as Fig. \ref{['fig:evolution_U0.3']}, but for a particle with $U_\infty = 0.3$ ($\mathcal{A}_0$=2.4, $e_0$=0.59, $i_0$=5.5$^\circ$, loosely-coupled case). The blue dashed line denoted $\mathcal{A}_\text{max}$ (Eq. \ref{['eq:def-amax']}) for this initial condition.
• Figure 5: 3D rotation of $\bm{U}_\infty$ (red) to $\bm{U}_\infty^\prime$ (blue) after a planetary encounter. The strength of the rotation is given by the angle $\gamma$ between the two vectors, whereas the direction of the rotation is oriented by the the counter-clockwise angle $\psi$.