Binary asteroids are increasingly important targets for both science and space mission planning. Researchers estimate that roughly 16 percent of near-Earth asteroids may exist in binary or multiple systems, although the exact fraction remains uncertain because many companions are difficult to detect. Understanding their future trajectories is essential when assessing impact risks, designing spacecraft encounters, or considering the extraction of resources from small bodies. Unlike the classical two-body problem, in which objects are often treated as ideal point masses, a full two-body model must account for the shapes, orientations, and internal distributions of mass within both bodies.

The challenge becomes especially clear when the asteroids are irregular. Their gravitational attraction does not act as though all of their mass were concentrated at a single point. Instead, the force and torque experienced by each body change as its shape rotates relative to the other. Translation and rotation therefore become tightly coupled: a change in orientation can alter the gravitational force, while that force can change both the orbit and the spin. Even small numerical inaccuracies in the attitude of a simulated asteroid can eventually produce incorrect forces and torques, sending the predicted system onto a physically unrealistic path.

To investigate this problem, a team led by Yongxin Guo of Liaoning University modeled each asteroid pair as a rigid dumbbell. In the double-dumbbell representation, each body consists of two asteroid-like masses joined by a massless rod. This model is simple enough to permit detailed mathematical analysis while retaining the essential features of an extended, rotating object. The researchers described the position and orientation of the dumbbells together using the special Euclidean group, SE(3), a mathematical structure that combines three-dimensional translations with three-dimensional rotations.

Using SE(3) avoids the singularities that can occur when rotations are represented with Euler angles. It also provides a natural language for describing rigid-body motion, but directly integrating equations on a Lie group can be computationally demanding. The researchers therefore introduced a body-fixed coordinate system attached to the second dumbbell and expressed the motion in terms of the first dumbbell’s relative position and relative attitude. This reduction removes redundant information associated with the overall motion of the pair and focuses the calculation on the variables that determine the internal orbital evolution.

The central innovation is the construction of a Hamel variational integrator from the discrete form of Hamilton’s principle. Variational integrators are designed to preserve the geometric structure of mechanical systems rather than simply approximating differential equations step by step. That distinction is critical in long-term orbital simulations. Standard explicit Runge–Kutta methods can produce highly accurate results over short intervals, but they generally do not preserve the symplectic structure associated with Hamiltonian mechanics or the orthogonality of rotation matrices. Over many thousands or millions of time steps, energy and angular momentum may drift, while a rotation matrix can gradually cease to represent a physically valid rotation.

The Hamel approach formulates the discrete equations on the Lie algebra associated with SE(3), rather than directly on the group elements themselves. In practical terms, the calculation works with local velocity-like and momentum-like quantities in the body frame before reconstructing the updated position and orientation. The researchers derived discrete Euler–Lagrange and Hamiltonian equations and used a discrete Legendre transform to create an iterative update procedure. At each time step, implicit equations are solved for the relative translational and rotational variables, followed by updates to the momentum variables and the full system state.

This formulation is intended to combine the advantages of geometric integration with lower computational cost. Lie group variational integrators preserve important structures, but their implicit equations are typically written directly in terms of group elements, which can increase the expense of solving each step. In the Hamel variational integrator, the corresponding equations are expressed on the Lie algebra. According to the study, this leads to a more efficient iterative process while retaining the structure-preserving behavior required for reliable long-duration simulations.

The numerical tests compared the Hamel variational integrator with a Lie group variational integrator and a conventional Runge–Kutta method. In regular-shaped double-dumbbell examples, the simulated kinetic and potential energies exchanged periodically as the bodies moved through their orbit. At approximately the ninth unit of simulated time, the dumbbells reached their closest separation, where the potential energy reached a minimum and the kinetic energy reached a maximum. Despite that exchange, the total energy remained nearly constant in the structure-preserving simulations. The Hamel method showed particularly strong performance in limiting energy error and maintaining the orthogonality of the rotation matrices.

The differences became more pronounced when the researchers introduced irregular shapes. The Runge–Kutta calculation gradually lost the orthogonality condition required of a rotation matrix. Once that condition is violated, the simulated orientation no longer corresponds exactly to a rigid-body rotation, and the resulting force and torque calculations can accumulate substantial errors. The Hamel method, by contrast, maintained low orthogonality and energy errors over the simulated interval. It also performed slightly better than the Lie group variational integrator in the reported CPU-time comparison, suggesting that the lower-dimensional algebraic formulation can provide a practical advantage without sacrificing geometric fidelity.

The simulations also revealed that shape irregularity can influence the motion of the binary system in a measurable way. When the bodies were modeled with nonuniform geometries, the trajectory of their barycenter developed a noticeable deviation in a particular direction, especially along the simulated y-axis. This result underscores a point that is easy to overlook in simplified asteroid models: even when the total mass and broad dimensions of an object remain similar, its detailed geometry can modify the gravitational interaction and alter the long-term path of the system. For missions approaching a binary asteroid, those differences could affect navigation, timing, landing-site selection, and the interpretation of observations.

The study does not suggest that every asteroid mission will immediately require a double-dumbbell Hamel model. Real small bodies can have complex shapes, rubble-pile interiors, deformation, surface mass movement, and non-gravitational forces caused by sunlight or outgassing. Nevertheless, the new method addresses a fundamental numerical problem that appears in increasingly realistic simulations: how to preserve the mathematics of rigid-body motion while keeping calculations efficient enough for extended orbital studies. By linking symplectic integration with Lie group geometry and an algebra-based computational scheme, the work offers a new tool for investigating binary asteroid evolution.

As planetary defense programs and asteroid exploration missions move toward more ambitious targets, reliable long-term prediction will become increasingly important. A tiny numerical bias that seems harmless during a short calculation can become a major trajectory error after repeated orbital encounters and rotations. The Hamel variational integrator provides a way to reduce that risk by preserving energy behavior, angular-momentum structure, and the validity of rotational motion more effectively than conventional approaches. The method could ultimately help scientists model the hidden complexity of binary asteroids—and make predictions about these gravitationally coupled worlds more trustworthy.

Subject of Research: Numerical modeling of the full two-body dynamics and long-term orbital evolution of binary asteroid systems.

Article Title: The Orbital Evolution of 2 Binary Asteroid Systems Based on Hamel’s Variational Integrators

News Publication Date: 16-Jul-2026

Web References: https://doi.org/10.34133/space.0450

References: Space: Science & Technology, DOI: 10.34133/space.0450

Image Credits: Space: Science & Technology

Keywords

Binary asteroids, asteroid dynamics, planetary defense, orbital evolution, full two-body problem, Hamel variational integrator, Lie group methods, symplectic integration, rigid-body dynamics, space science