The core idea is deceptively simple. Instead of solving the full multiple-gravity-assist trajectory optimisation problem for every Kuiper Belt Object, the team trained a Gaussian process regression model to predict how expensive a transfer would be, based only on the orbital characteristics of the target. The training labels for this model came from genuine high-fidelity optimisations: for a set of representative targets, the researchers ran a complete optimisation of an Earth–Jupiter–Kuiper Belt Object transfer, including the deep-space manoeuvres that allow a spacecraft to adjust its course between gravity assists. Each optimisation produced two numbers that matter enormously to mission planners: the total deep-space manoeuvre velocity increment, essentially the fuel cost of the transfer, and the time of flight. Those pairs of numbers became the ground truth from which the surrogate learned.

A crucial design decision was how to represent each target to the machine-learning algorithm. The team compared seven different ways of encoding a Kuiper Belt Object’s orbit, ranging from classical orbital elements such as semimajor axis, eccentricity and inclination, to Cartesian, cylindrical and spherical position-and-velocity states sampled at fixed epochs. This choice matters because a Gaussian process model builds its predictions from similarities between inputs, and a representation that hides the physically relevant structure of the problem will degrade the predictions even if the underlying data are identical. After systematic comparison, a multi-epoch Cartesian representation, in which the target’s heliocentric state is recorded at several fixed times, gave the most balanced performance across both prediction of transfer cost and ranking of targets by accessibility.

Two further refinements sharpened the model’s usefulness. First, the researchers applied a logarithmic scaling to the transfer objective, which concentrated the model’s accuracy in the low-cost region of the distribution, precisely where mission designers care most, because targets with transfer costs below roughly one kilometre per second of deep-space manoeuvring are the ones worth pursuing. Second, they employed uncertainty-guided adaptive sampling, a technique rooted in Bayesian optimisation in which the model itself identifies the regions of input space where its predictions are least confident, and new expensive optimisations are run there to improve the training set. This active-learning strategy measurably improved both the precision and the recall of the screening process at the tested screening fractions of 10 and 20 percent, meaning the model both missed fewer genuinely accessible targets and wasted less attention on inaccessible ones.

The quantitative results are striking. On an independent test set of targets the model had never seen, the Gaussian process regression achieved a global mean absolute error of 0.673 kilometres per second on the transfer objective and 537.5 days on the time of flight. Those global figures, however, mask a far more encouraging picture in the region that matters. For cases whose validated transfer objective was below one kilometre per second, the mean absolute error dropped to just 0.133 kilometres per second, and for cases the surrogate itself selected as having predicted costs below that threshold, the error was 0.245 kilometres per second. The larger global mean reflects the sparse high-cost tail of the objective distribution, where a handful of extremely expensive transfers dominate the average error but are of little practical interest for mission selection.

The demonstration that ties the method together was applied to a catalogue of 2,223 representative Kuiper Belt Object targets. Running the surrogate across this population, the team identified the ten targets with the lowest predicted transfer objectives, candidates that would then be passed on for full multiple-gravity-assist and deep-space-manoeuvre validation. This two-stage workflow, cheap surrogate screening followed by expensive high-fidelity confirmation, is the same philosophy that has transformed other branches of engineering design, and the study shows it can be brought to bear on one of the most computationally demanding problems in astrodynamics. The researchers note that the approach is demonstrated for the prescribed Earth–Jupiter–Kuiper Belt Object transfer scenario, with the gravity-assist sequence fixed in advance.

The technical machinery underneath the method draws on a rich literature. Lambert’s problem, the classical boundary-value problem of orbital mechanics that connects two positions in space with a transfer orbit, underpins the generation of candidate trajectory legs, and modern robust solvers have made such computations fast and reliable. The global optimisation of the full trajectory, with its mixture of continuous variables such as departure epochs, manoeuvre magnitudes and flyby geometry, relies on evolutionary algorithms such as differential evolution, which have proven robust for the multimodal landscapes that trajectory optimisation produces. Gaussian process regression, meanwhile, brings a distinctive advantage over other machine-learning approaches: it provides not just a point prediction but an associated uncertainty estimate for every input, which is exactly what makes uncertainty-guided adaptive sampling possible and what allows a mission designer to know how much to trust a given screening result.

Why does this matter beyond the mathematics? The Kuiper Belt is one of the last great unexplored reservoirs of the solar system, a ring of icy bodies beyond Neptune that preserves a record of the primordial material from which the planets formed. NASA’s New Horizons mission gave humanity its first close look at a Kuiper Belt Object when it flew past Arrokoth in 2019, and concepts for an interstellar probe and follow-on Kuiper Belt missions continue to circulate in the planetary science community. But the sheer number of known and expected objects in this region means that target selection is itself a formidable optimisation problem. Every candidate target implies a different launch window, a different gravity-assist choreography and a different fuel budget, and choosing poorly can cost a mission a decade of flight time or render it infeasible altogether. Rapid, reliable accessibility assessment is therefore not a luxury but a prerequisite for ambitious outer-solar-system exploration.

The study also fits into a broader trend in which machine learning serves as an accelerant for, rather than a replacement of, classical astrodynamics. Earlier work demonstrated that Gaussian process regression could assess the accessibility of main-belt asteroids, and neural networks have been applied to orbit uncertainty propagation and estimation. The present research extends that paradigm to the most distant and most expensive class of targets, and its careful comparison of state representations and sampling strategies offers a practical template for other teams. The authors acknowledge support from the National Natural Science Foundation of China, and the datasets generated during the study are available from the corresponding author on reasonable request.

There remain, of course, limits to what any surrogate can do. The model was trained and validated for a specific transfer scenario with a prescribed gravity-assist sequence, and its accuracy in the low-cost region, while excellent, comes at the price of larger errors in the sparse high-cost tail. The ten targets flagged by the surrogate still require full optimisation before any mission commitment, and the screening fractions of 10 and 20 percent represent a trade-off between computational savings and the risk of overlooking a hidden gem. Yet the headline numbers speak for themselves: a model that predicts transfer costs to within a few hundred metres per second for the targets that matter, ranks thousands of destinations in seconds, and hands mission designers a principled, uncertainty-aware shortlist. As the catalogue of Kuiper Belt Objects continues to grow, tools like this one may well decide which frozen worlds humanity visits next.

Subject of Research: Surrogate-based machine learning for rapid assessment of spacecraft transfer accessibility to Kuiper Belt Objects

Article Title: Surrogate-based rapid accessibility assessment of Kuiper Belt object transfers

Article References: Meng, L., Zhou, X., Qiao, D., & Li, X. (2026). Surrogate-based rapid accessibility assessment of Kuiper Belt object transfers. Astrophysics and Space Science, 371(10), Article 116. https://doi.org/10.1007/s10509-026-04650-9

Image Credits: AI Generated

DOI: 10.1007/s10509-026-04650-9

Keywords: Kuiper Belt objects, trajectory design, Gaussian process regression, surrogate modelling, gravity assists, deep-space manoeuvres, adaptive sampling, astrodynamics, mission design, machine learning, interplanetary transfers, Bayesian optimisation

Recent groundbreaking research from Michigan State University has shed light on the processes that craft these two-lobed formations. Utilizing a state-of-the-art high-performance computing system, graduate student Jackson Barnes has developed the first computational simulation that naturally forms contact binaries through gravitational collapse, without relying on improbable or exotic events. Published in the Monthly Notices of the Royal Astronomical Society, this work opens new avenues for understanding the early evolutionary pathways of small bodies in the outer solar system.

Traditional models faced significant limitations, often approximating these small icy objects as fluid blobs that, upon collision, merged into singular spheres. Such simplifications failed to reproduce the characteristic dual-lobed structure observed in about 10% of Kuiper Belt planetesimals. Barnes’ simulations mark a breakthrough by incorporating the mechanical strength and granular nature of these bodies. His approach allows the simulated planetesimals to rest against each other, maintain their distinct shapes, and ultimately fuse gently rather than violently.

The insights from Barnes’ research are critical because they align with the observed abundance of contact binaries. If 10% of planetesimals exhibit this fused shape, the formation mechanism must be a relatively common event in the early solar system, rather than a product of rare or catastrophic phenomena. Earth and Environmental Science Professor Seth Jacobson, a senior author on the paper, emphasizes that gravitational collapse is a compelling and elegant explanation consistent with empirical data acquired through decades of observation.

NASA’s New Horizons mission provided the first close-up images of a contact binary in January 2019 when it flew past the Kuiper Belt object known as 2014 MU69, nicknamed Ultima Thule. These crisp images revealed a distinctly two-lobed shape with smooth lobes fused at a narrow neck, challenging prior assumptions about planetesimal formation. Following this discovery, astronomers revisited other Kuiper Belt objects and identified that approximately one in ten follows this binary configuration, with little evidence of disruptive collisions owing to the sparse population density in that cosmic neighborhood.

The Kuiper Belt, formed remnant from the protoplanetary disk that once encircled the Sun, is an archive of primordial matter dating back over four billion years. Planetesimals are among the first large solid bodies to arise from this disk, developing through the slow agglomeration of pebble-sized fragments pulled together by mutual gravitational attraction. This formative stage is analogous to compaction of snowflakes into a snowball, except occurring over cosmic time scales and within a rotating circumstellar environment.

Barnes’ simulations highlight a fascinating dynamical process: as a rotating cloud of pebbles collapses under gravity, irregularities often lead to the initial formation of binary systems—two planetesimals orbiting each other. Over time, their orbits decay, spiraling closer until they make contact gently. The simulated binaries retain their smooth, rounded shapes without blending into a single sphere, thus reproducing the iconic snowman-like morphology observed in actual Kuiper Belt objects.

A key question that arises is how these delicate binary structures persist over billions of years without disruption. Barnes explains that the Kuiper Belt’s low-density environment minimizes chances of catastrophic collisions that could separate or shatter these contact binaries. This tranquil setting preserves the integrity of their shapes, consistent with the lack of significant cratering seen on many observed binaries.

While the gravitational collapse hypothesis had been proposed before, quantitative and realistic modeling was lacking due to computational constraints and oversimplifications. Barnes’ work pioneers a physics-rich simulation capable of resolving the mechanical and dynamical subtleties necessary to form and sustain contact binaries. This represents a major advancement in small-body astrophysics.

Looking forward, Barnes anticipates that his model will inspire further studies examining more complex multi-lobed systems, where three or more bodies coalesce through related mechanisms. The research team also aims to refine their simulations by incorporating more detailed physics to replicate the collapse and accretion processes with even greater fidelity.

Moreover, ongoing and future space missions venturing into the outer solar system may uncover additional contact binaries, revealing whether these “cosmic snowmen” have distant, untapped cousins. Such discoveries will further deepen our understanding of the delicate balance between gravitational forces and collisional histories that shape the architecture of our solar system’s frontier.

This new insight into the origin of contact binary planetesimals marks a significant milestone in planetary science. It not only clarifies how these peculiar objects form but also enhances our comprehension of the early conditions and evolutionary processes that govern the distant Kuiper Belt. As computational capabilities continue to expand, such interdisciplinary efforts bridging observation, theory, and simulation promise to unravel even more cosmic mysteries.

Subject of Research: Formation of contact binary planetesimals in the Kuiper Belt through gravitational collapse

Article Title: Direct contact binary planetesimal formation from gravitational collapse

News Publication Date: 19-Feb-2026

Web References: http://dx.doi.org/10.1093/mnras/stag002

Image Credits: NASA

Keywords

Kuiper Belt, contact binaries, planetesimals, gravitational collapse, New Horizons, solar system formation, computational simulation, binary planetesimals, outer solar system, planetary science