Chaotic Dynamics and Zero-Velocity Structures in the Pluto-Charon CR3BP

Pluto and Charon are a dwarf binary system with a high mass ratio $μ$, preventing Trojan companions. This instability creates ideal intersections for low-energy pathways that spacecraft can traverse and serves as an important test case for stability at Lagrange points for high $μ$ binaries. This study models the Pluto-Charon system in the planar Circular Restricted Three-Body Problem (CR3BP) and compares the tadpole and horseshoe orbits of a massless particle with known low-$μ$ orbits. Moreover, it compares the trajectories for instances where the L1 neck is opened and closed, and exhibits the corresponding zero-velocity curves; RK4 integration was used to update the position and velocity of the particle. The simulations consistently showed chaotic and unpredictable trajectories, where small changes in initial parameters could completely alter the results; trajectories displayed the influence of higher values of $μ$ on the stability of binary systems. These findings confirm that Pluto-Charon cannot host long-lived Trojan companions and instead behaves as a system of chaotic transport. Future work could extend the study to include the spatial or elliptic three-body problem, further refining the understanding of instability in high-$μ$ binaries.

💡 Research Summary

The paper investigates the dynamical environment of the Pluto–Charon binary using the planar Circular Restricted Three‑Body Problem (CR3BP). Because the mass ratio μ≈0.109 is unusually high for a planetary system, the classical triangular Lagrange points L₄ and L₅ are linearly unstable, and the entire set of equilibrium points (L₁–L₅) behave differently from low‑μ systems such as Sun–Earth or Sun–Jupiter. The authors normalize the system (distance, total mass, and angular velocity set to unity), place Pluto at (‑μ,0) and Charon at (1‑μ,0), and compute the effective potential Ω, the zero‑velocity function Z=2Ω, and the Jacobi constant C=Z−(ẋ²+ẏ²). Zero‑velocity curves (ZVCs) are plotted for several values of C, illustrating how tiny changes (ΔC≈10⁻³–10⁻²) open or close the “necks” around the Lagrange points, thereby altering the connectivity of allowed regions.

The equations of motion in the rotating frame, ẍ = ∂Ω/∂x + 2ẏ and ÿ = ∂Ω/∂y – 2ẋ, are integrated with a fourth‑order Runge‑Kutta (RK4) scheme using a time step h = 5×10⁻⁴. After each step the Jacobi constant is recomputed to verify its conservation, ensuring numerical fidelity.

The study first reproduces the classic tadpole (libration around a single triangular point) and horseshoe (libration encompassing L₃, L₄, and L₅) orbits that are known to exist when μ ≤ 0.03852 (Routh’s critical mass ratio). In those low‑μ cases the trajectories are smooth, periodic, and remain confined for arbitrarily long times. When the same initial conditions are applied to the Pluto–Charon system, however, the trajectories quickly become chaotic.

Key results include:

1. Tadpole attempt – Starting near L₄ with a small negative offset in C (δ = –5×10⁻⁴) leads to rapid divergence; the particle collides with Charon at nondimensional time T ≈ 37, demonstrating the linear instability of the triangular points at high μ.

2. Horseshoe attempt – With C set just below the L₃ value (δ = –1×10⁻³), the particle initially follows a horseshoe‑like path but after only two librations the orbit breaks symmetry, the zero‑velocity curves expand, and the particle either escapes or collides. The larger δ makes the chaotic deviation visually obvious.

3. L₁ neck sensitivity – By varying δ by ±10⁻², the neck at L₁ toggles between open and closed. When the neck is closed the particle lingers near Pluto before reflecting; when open, the same particle transits into Charon’s lobe. The outcome is highly sensitive to the sign of the tiny energy perturbation, a hallmark of chaotic dynamics.

Across all simulations the Jacobi constant drift remains on the order of 10⁻⁹–10⁻¹³, confirming that the observed chaos is not a numerical artifact but an intrinsic property of the high‑μ system. The authors conclude that Pluto–Charon cannot host long‑lived Trojan (L₄/L₅) companions; instead, the binary acts as a chaotic transport network where low‑energy pathways exist but are highly unstable.

The paper suggests extending the analysis to the spatial (3‑D) CR3BP and to the elliptic restricted problem, which would incorporate additional degrees of freedom and time‑varying primaries. Such extensions could refine the identification of transient invariant manifolds and improve low‑fuel trajectory design for missions to high‑μ binaries, while also deepening our theoretical understanding of stability limits in binary dwarf‑planet systems.