The researchers tackled the two-impulse tangential phasing problem for a 2:1 distant retrograde orbit within the planar circular restricted three-body problem, the classical dynamical model that treats Earth and Moon as massive bodies circling their common center while a spacecraft responds to both gravitational fields. Phasing maneuvers allow a spacecraft to change its position along the orbit by a prescribed phase angle, typically using two impulsive burns. Earlier work had produced only discrete, isolated solutions, leaving mission planners without a global picture of how these maneuvers connect, continue, or change character as parameters vary.

The key insight was to impose a tangential constraint on the burns, requiring each impulse to align with the spacecraft’s local velocity direction. This restriction, far from limiting the analysis, collapses the design space into an elegant two-dimensional map that can be traced into continuous families using multiple shooting and pseudo-arclength continuation, numerical techniques that follow solution curves through the full parameter space even where simple iteration would fail.

The mathematical heart of the paper is an analytically established Z2 equivariance property of the phasing equations. Equivariance means that when the solution variables are mirrored through a symmetry transformation—in this case a reflection combined with time reversal—the governing equations transform in a precisely matching way. The authors proved that this property holds not just for the 2:1 distant retrograde orbit but for any planar periodic orbit symmetric about the x-axis in the restricted three-body problem, a result with broad applicability across cislunar trajectory design.

This symmetry does more than look beautiful on paper. It yields a rigorous classification of the phasing families: the team identified two symmetric families and two conjugate pairs of asymmetric families. Crucially, the equivariance predicts symmetry-breaking pitchfork bifurcations—points at which a symmetric solution loses stability or existence and spawns a pair of mirror-image asymmetric solutions. These bifurcations are the connective tissue of the phasing landscape, and the researchers used their symmetry-based insight to build an efficient search pipeline that pinpoints the exact bifurcation locations without exhaustive numerical scanning.

The payoff for mission design is striking. The families interconnect through two pitchfork bifurcations and collectively span the entire phase-shift range from negative pi to pi within a maximum total impulse of just 92.7 meters per second and a transfer time of one sidereal month. In practical terms, a spacecraft can achieve any desired repositioning along the 2:1 distant retrograde orbit within those modest budgets, with every option now documented and reachable through continuous solution curves rather than lucky numerical hits.

Beyond the two-impulse problem, the authors demonstrated a symmetry-based multi-segment construction principle that yields free multi-revisit trajectories—paths that repeatedly return to designated points without additional maneuvering cost in the simplified model. Such trajectories serve as high-quality seeds for higher-fidelity design, particularly for contingency-recovery scenarios in which a spacecraft must reestablish a required phasing condition after an anomaly, a concern that has moved from theoretical to urgent as crewed and robotic lunar operations multiply.

The work lands at a moment when distant retrograde orbits are under intense scrutiny. Their long-duration stability, documented in prior research, makes them attractive for depoting, servicing, and as nodes in cislunar logistics networks. Phasing capability is the operational backbone for any such role: rendezvous, constellation deployment, and recovery from missed burns all depend on knowing exactly which transfer options exist and how much they cost. By transforming a catalog of isolated solutions into a structured, symmetry-organized atlas, the new study gives engineers a tool closer to a map than to a list.

Methodologically, the paper also showcases the power of dynamical systems thinking in astrodynamics. The proof of equivariance, built on the mirror theorem for symmetric periodic orbits and a careful linear-algebraic argument showing the relevant symmetry operators are similar involutions, demonstrates that abstract bifurcation theory can directly accelerate practical trajectory computation. The predicted pitchfork structure guided the numerical continuation, ensuring that no family branch was missed and that transitions between symmetric and asymmetric maneuvering were found to machine precision.

As agencies and companies lay plans for sustained activity between Earth and Moon, studies like this one quietly determine what is possible. With the full family structure of tangential phasing around a 2:1 distant retrograde orbit now characterized, mission designers can select transfers with confidence, explore asymmetric options previously hidden by symmetry-breaking, and construct multi-revisit patterns for resilient operations. The result is a cislunar highway network described not by scattered signposts but by the deep geometric grammar of the three-body problem itself.

Subject of Research: Two-impulse tangential phasing maneuver families for a 2:1 distant retrograde orbit in the Earth–Moon system, characterized through equivariance and pitchfork bifurcations.

Article Title: Characterization of tangential-maneuver phasing families for a 2:1 distant retrograde orbit via equivariance and pitchfork bifurcations

Article References: Wang, Y., Zhang, C., & Zhang, H. (2026). Characterization of tangential-maneuver phasing families for a 2:1 distant retrograde orbit via equivariance and pitchfork bifurcations. Celestial Mechanics and Dynamical Astronomy, 138(5), Article 58. https://doi.org/10.1007/s10569-026-10330-1

Image Credits: AI Generated

DOI: 10.1007/s10569-026-10330-1

Keywords: distant retrograde orbit, cislunar space, orbital phasing, three-body problem, pitchfork bifurcation, equivariance, trajectory design, spacecraft rendezvous, Earth–Moon system, bifurcation analysis, astrodynamics, periodic orbits

The CR3BP reduces the gravitational dance of three bodies to its essential form: an infinitesimally small mass, such as a spacecraft, moving under the influence of two primaries that orbit their shared center of mass in perfect circles. Within this simplified but still chaotic system, trajectories fall into three fundamental categories. Periodic orbits repeat themselves exactly, tracing closed loops through phase space. Quasi-periodic orbits wander across the surface of an invisible torus, never quite repeating but never straying far from home. Chaotic trajectories, by contrast, diverge unpredictably from their initial conditions, sensitive to the smallest perturbation. Distinguishing among these behaviors is the central task of trajectory design, because a mission planner must know whether a candidate path will remain stable or spiral into chaos.

The traditional tool for this task is the Poincaré map, a technique dating back to Henri Poincaré’s foundational work in 1890. Rather than tracking a trajectory continuously, the map records only the points where the path pierces a chosen surface in phase space, converting a continuous curve into a discrete scatter of crossings. In the planar version of the problem, this produces a two-dimensional plot in which periodic orbits appear as fixed points, quasi-periodic orbits form distinctive chains of islands, and chaotic motion scatters into structureless dust. Mission designers have long relied on visual inspection of these maps to identify promising orbits, a process that works well enough in two dimensions but collapses entirely when the problem extends into three spatial dimensions.

Extending the Poincaré map to four dimensions is necessary because real spacecraft do not confine themselves to planes. To capture out-of-plane motion, researchers use a space-plus-color representation: two spatial coordinates occupy the horizontal axes, the vertical coordinate is plotted on a third axis, and the out-of-plane velocity is encoded as color. The result is a rich but visually cluttered dataset in which the familiar island chains of two-dimensional maps vanish, replaced by three-dimensional structures with no established taxonomy. Previous efforts to sift through these maps relied on filtering by the Fast Lyapunov Indicator, a numerical measure of chaos, but this remained a manual, time-intensive operation that yielded only preliminary catalogs of structures. With each map containing thousands to millions of data points, the need for automation became unmistakable.

The new pipeline, developed by Kevin M. Trigg, Daniel J. Broyles, Robert A. Bettinger, and Tyler J. Kapolka, addresses this challenge through a carefully engineered sequence of computational steps. The researchers simulated 10,000 seed trajectories for each of 11 different Jacobi constants, a parameter that acts like an energy level and determines which regions of phase space a spacecraft can reach. Each trajectory was propagated for 1,000 time units using high-precision numerical integration, with early termination for paths that collided with a primary body or escaped the system. The team then faced a fundamental data problem: trajectories produce variable numbers of Poincaré crossings, and clustering algorithms require inputs of fixed length. Their solution was a feature engineering scheme that compresses each trajectory into a 23-dimensional numerical fingerprint.

That fingerprint draws on three complementary strategies. Statistical and dynamical descriptors capture the bounding box of each trajectory’s crossings across position and velocity coordinates, along with the Fast Lyapunov Indicator, which quantifies sensitivity to initial conditions and separates regular from chaotic behavior. Geometric descriptors treat each trajectory’s crossings as a miniature dataset, applying clustering within the trajectory itself to measure how cohesive or fragmented its structure is; a well-formed invariant torus yields dense, orderly crossings, while chaos produces scattered noise. Finally, frequency-domain analysis applies a Fast Fourier Transform to the ordered sequence of crossings, extracting the three strongest oscillatory modes and their amplitudes as signatures of recurrence. Periodic and quasi-periodic trajectories concentrate their spectral energy in sharp peaks, whereas chaotic trajectories spread it broadly.

With feature vectors in hand, the pipeline applies Principal Component Analysis to compress the representation, retaining components that explain at least 95 percent of the variance and typically reducing 23 features to 11. The reduced vectors then feed into HDBSCAN, a hierarchical density-based clustering algorithm chosen for its ability to handle clusters of varying density and shape without requiring the number of clusters to be specified in advance. This flexibility matters because Poincaré maps contain structures with irregular boundaries and non-uniform density, conditions that defeat simpler methods like DBSCAN, which relies on a single global density threshold. HDBSCAN also explicitly labels low-density trajectories as noise, providing a natural mechanism for isolating chaotic orbits that belong to no coherent structure. The researchers compared HDBSCAN against agglomerative clustering, spectral clustering, affinity propagation, and Gaussian mixture models, finding that while competitors achieved slightly higher Silhouette scores, HDBSCAN consistently produced lower Davies-Bouldin scores, better structural similarity, and the crucial advantage of automatic noise detection.

The analysis revealed eight distinct dynamical structures in the four-dimensional maps, which the authors named descriptively in the absence of any formal taxonomy: Figure-Eight, Tube, Scorpion Tail with Dots, Scorpion Tail without Dots, Pillar and Shield, Pillar and Shapes, Two Planes, and Three Pillars. These geometries bear no resemblance to the island chains of planar maps, underscoring how radically the topology changes when out-of-plane motion is included. The Figure-Eight structures, associated with quasi-periodic motion around Lagrange points, dominated the regular regions of the maps. An ablation study confirmed that every feature group contributed essential information: removing the Fast Lyapunov Indicator degraded the separation between regular and chaotic regions, while removing statistical features collapsed intra-cluster cohesion. When the tuned clustering configuration was applied unchanged to maps at other Jacobi constants, it recovered consistent structures with noise proportions holding steady between roughly 20 and 30 percent, demonstrating genuine robustness across the energy landscape.

Perhaps the most ambitious component of the work is its approach to orbit continuation. In classical astrodynamics, tracing how an orbit family evolves requires differential correction and continuation methods that follow a trajectory as parameters change. The researchers instead built a data-driven approximation: they treated each Jacobi constant level as a layer in a directed graph, connected each trajectory to its three nearest neighbors in feature space within the adjacent layer, and used depth-first search to find chains spanning the entire range from C equals 3.18 down to 2.68. Chains were ranked by average feature distance, with the smoothest paths representing candidate continuations of dynamical families. The results confirmed that the learned feature representations preserve meaningful dynamical similarity across energy levels, with Figure-Eight structures persisting coherently across maps. Yet the method also exposed its own limitation: only 24 unique structures appeared among the top 1,000 chains, revealing a strong bias toward the dominant, tightly clustered geometries and a need for additional constraints to encourage exploration of rarer structures.

The implications extend beyond the Earth-Moon system. The authors emphasize that the methodology applies to any multi-body gravitational environment, and that four-dimensional Poincaré maps represent a growing frontier in astrodynamics research. Future work includes cataloging how these structures evolve with mass parameter and energy, comparing the discovered orbits against known periodic orbit families and invariant tori generated by traditional continuation methods, and developing algorithms to locate fixed points in maps where, unlike the planar case, they are not intuitively positioned at the centers of island chains. The pipeline also contributes to a broader scientific goal: automated detection of invariant structures in high-dimensional Hamiltonian systems, a challenge that reaches into plasma physics, celestial mechanics, and accelerator design. For mission planners navigating the increasingly crowded cislunar arena, where spacecraft such as those supporting lunar exploration must exploit subtle gravitational structures to conserve fuel, a tool that converts millions of trajectory crossings into organized, labeled dynamical families could shorten the path from concept to flight-ready trajectory design.

Subject of Research: Unsupervised machine learning for discovering dynamical structures in 4D Poincaré maps of the circular restricted three-body problem

Article Title: Application of machine learning to discover dynamical structures in 4D Poincaré maps in the circular restricted three-body problem

Article References: Trigg, K. M., Broyles, D. J., Bettinger, R. A., & Kapolka, T. J. (2026). Application of machine learning to discover dynamical structures in 4D Poincaré maps in the circular restricted three-body problem. Astrophysics and Space Science, 371(9), Article 107. https://doi.org/10.1007/s10509-026-04638-5

Image Credits: AI Generated

DOI: 10.1007/s10509-026-04638-5

Keywords: Poincaré map, circular restricted three-body problem, machine learning, HDBSCAN, quasi-periodic orbits, chaos, Fast Lyapunov Indicator, principal component analysis, astrodynamics, Earth-Moon system, trajectory design, clustering

The research centers on Milanković cycles, the periodic variations in Earth’s orbital eccentricity, axial tilt, and precession that redistribute the sunlight reaching the planet and thereby pace long-term climate change. In the modern Solar System, these cycles operate on well-known timescales, and they have been used to calibrate the geological time scale for the Cenozoic era with remarkable accuracy. Extending that approach into the Precambrian, however, has been hampered by a fundamental problem: numerical models of the Solar System’s orbital motion become chaotic and lose predictive power over tens of millions of years, and no astronomical solution can currently be trusted beyond roughly the last 100 million years. For intervals more than a billion years in the past, scientists have had to rely on theory alone to estimate how orbital frequencies differed from today’s values.

Sedimentary rocks offer a way around this limitation. When climate cycles driven by orbital variations imprint regular patterns on accumulating sediment—alternations between more resistant and more recessive beds, for example, or rhythmic changes in grain size and composition—the resulting cyclostratigraphy can be read as a recording of the astronomical forcing that produced it. The Hakatai Shale, deposited in shallow-water settings roughly 1.4 to 1.1 billion years ago during the Mesoproterozoic era, preserves such rhythms in striking detail. The team logged and analyzed stratigraphic sections at Red Canyon and Tapeats Creek within Grand Canyon National Park, conducting fieldwork under permit from the National Park Service, and measured the thickness and character of successive sedimentary cycles with centimeter-scale resolution.

The key to interpreting these rhythms lies in the physics of the Earth–Moon system. Tidal friction, the braking effect of lunar tides on Earth’s rotation, has steadily slowed the planet’s spin over geological time while pushing the Moon gradually farther away. As the day lengthens, the frequency of the climatic precession cycle—the wobble in Earth’s axis that changes how seasons align with the planet’s position around the Sun—changes in a predictable way. Because the precession signal modulates the amplitude of the eccentricity cycle, sedimentary records that capture both can be used to solve for the precession constant and, from it, the ancient Earth–Moon distance and length of day. This approach, known as TimeOpt and its Bayesian extension TimeOptBMCMC, was applied to the Hakatai Shale using the Astrochron software package, with 100,000 Monte Carlo samples used to constrain the statistical uncertainty of the reconstruction.

The analysis of the Tapeats Creek composite record, corrected for variations in sediment thickness, revealed a coherent suite of astronomical signals. The team identified cycles corresponding to climatic precession, orbital eccentricity, and obliquity, and used the ratios among them to test which cyclostratigraphic interpretation best fit the data. Among three competing interpretations of the dominant spectral peaks, the preferred option yielded sedimentation rates of a few centimeters per thousand years—values consistent with the quiet, low-energy depositional environments inferred independently from the rock’s lithology, which includes reworked microbial mat fabrics, wind-blown quartz grains, and pseudomorphs after evaporite minerals such as gypsum and anhydrite.

Beyond confirming that Milanković forcing operated in the Mesoproterozoic, the record delivered a surprise. The relative amplitudes of the astronomical forcing frequencies, particularly obliquity—the tilt of Earth’s spin axis—showed anomalous patterns compared with what present-day dynamics would predict. In the spectra of the Hakatai Shale, the strength of individual obliquity-related peaks shifted between different stratigraphic intervals in ways that mirror the behavior of state-of-the-art deep-time astronomical models, specifically the ZB23 solutions developed by Zeebe and colleagues, which extend orbital calculations back 3.5 billion years. In those models, the dominance of particular obliquity cycles changes through time as secular resonances among the planets drift in and out of critical configurations.

One such configuration involves the resonance angle associated with the motions of Mars and the inner planets, which can transiently disrupt the dominant obliquity cycle. Another involves a secular resonance that interferes with the main eccentricity cycle linked to the orbital frequencies of Earth and Jupiter. The Hakatai spectra show amplitude trends—weak expression of one eccentricity peak, enhanced power in a particular obliquity band—that are consistent with the models’ predictions for conditions around 1.2 billion years ago, including the possible influence of a resonance in which combinations of planetary orbital frequencies and Earth’s axial precession frequencies nearly coincide. The authors note that these patterns could also reflect a nonlinear climate response, in which interactions between multiple forcing frequencies generate combination tones that appear in the sedimentary record at sums and differences of the original periods.

Either interpretation carries weighty implications. If the amplitude anomalies record shifts in secular Solar System resonances, then the Grand Canyon strata provide the first empirical evidence from the rock record for how the gravitational architecture of the planetary system has changed over more than a billion years, complementing purely numerical approaches that are limited by chaos. If, instead, the signals arise from nonlinear climate dynamics, they illuminate how the Precambrian climate system responded to astronomical forcing in an atmosphere and ocean very different from today’s, before the rise of complex life and with substantially different greenhouse gas inventories. Distinguishing between these possibilities is a central goal of ongoing work, and the Bayesian inverse modeling framework applied here is designed to weigh such alternatives quantitatively.

The study builds on a growing effort to use geology as a probe of Solar System dynamics, sometimes described as mapping Solar System chaos with the geological record. Previous work by members of the team demonstrated that Milankovitch cycles preserved in 2.46-billion-year-old banded iron formations constrain the Earth–Moon system in the Paleoproterozoic, and theoretical studies have traced how tidal evolution reshaped the lunar orbit through resonant episodes. The Hakatai Shale now extends this empirical reach into the Mesoproterozoic with a record whose internal consistency—matching precession, eccentricity, and obliquity signals across two geographically separated sections—strengthens confidence that the rhythms are genuinely astronomical in origin rather than products of local tectonic or depositional noise.

The practical implications extend well beyond deep-time astronomy. Accurate knowledge of past astronomical frequencies underpins astrochronology, the dating method that uses orbital cycles to refine the geological time scale, and the new results constrain how those frequencies differed in the Precambrian, when shorter days and a closer Moon altered the pacing of climate cycles. The team’s cyclostratigraphic data and analyses have been made openly available through Zenodo, and the ZB23 astronomical solutions are publicly accessible, allowing other researchers to test and extend the reconstruction. As more ancient rhythmically deposited successions are examined with these tools, sedimentary rocks may continue to serve as long-term observatories of the heavens—recording, in ordinary mud, the slow gravitational conversation between Earth, the Moon, and the wandering planets.

Subject of Research: Reconstruction of Mesoproterozoic Earth–Moon dynamics and Solar System orbital evolution from Milanković cycles in Grand Canyon sedimentary strata

Article Title: Earth–Moon and Solar System history recorded in Mesoproterozoic Grand Canyon strata

Article References: Lantink, M. L., Eyster, A., Kocken, I. J., Meyers, S. R., & Zeebe, R. E. (2026). Earth–Moon and Solar System history recorded in Mesoproterozoic Grand Canyon strata. Nature Geoscience. https://doi.org/10.1038/s41561-026-02100-3

Image Credits: AI Generated

DOI: 10.1038/s41561-026-02100-3

Keywords: Milanković cycles, Hakatai Shale, Grand Canyon, Earth–Moon system, Mesoproterozoic, cyclostratigraphy, solar system dynamics, obliquity, lunar distance, astrochronology, secular resonance, tidal evolution

The restricted three-body problem is a classic model in gravitational dynamics. It considers two primary bodies—in this case, Earth and the Moon—that orbit one another, while a third body has negligible mass and does not alter their motion. Despite its apparently simple setup, the equations are nonlinear and can produce a remarkable range of behavior, including escape trajectories, temporary capture, unstable passages, libration-point orbits and long-lived periodic motions. In the circular restricted three-body problem, or CRTBP, Earth and Moon are assumed to travel in circular orbits. The elliptic version, known as the ERTBP, allows their separation and orbital speed to vary as they follow an ellipse. That seemingly modest change removes an important symmetry and makes the search for repeating trajectories substantially more difficult.

Nagai’s work focuses on resonant periodic orbits, in which different components of a spacecraft’s motion return to their original configuration after a precise number of cycles. The key idea is frequency matching. A trajectory can oscillate horizontally around the Earth–Moon system while also moving above and below the orbital plane. When the vertical oscillation frequency bears a rational relationship to the principal orbital frequency, the motion can close after a finite period instead of drifting indefinitely. The study concentrates on 1:n resonances, a class in which the relevant frequencies are related by an integer ratio. Such resonances are central to periodic-orbit theory because they transform what might otherwise be a quasiperiodic, non-repeating path into a closed orbit with a recognizable geometric pattern.

To identify these paths, the study begins with a simplified approximation of the equations of motion near the Moon. In that region, the gravitational influence of the Moon dominates the local motion, while Earth’s field and the rotating reference frame continue to shape the trajectory. The approximation makes it possible to understand the essential structure before confronting the full nonlinear equations. Nagai shows that the in-plane part of the approximate solutions is consistent with the rose curves associated with the Italian mathematician Guido Grandi, who studied these curves in the early eighteenth century. In polar-style form, the radial distance varies sinusoidally while the angular position advances with time. The resulting trajectory repeatedly expands and contracts, creating lobes or “petals” around a central region. The number and arrangement of these petals depend on the ratio of the frequencies and on the initial phases.

The mathematical connection is more than a visual coincidence. A rose curve can be written parametrically so that its radial amplitude follows one sinusoidal function while the direction of motion follows another. In Nagai’s formulation, the coordinates contain a factor of the form (\sin((n/N)t-\phi_1)), multiplied by the rotating directional terms (\sin(t-\phi_2)) and (\cos(t-\phi_2)). Here, (n) and (N) are positive integers, and the phase parameters determine the initial orientation and timing of the pattern. When the frequencies are commensurate—meaning their ratio is rational—the curve repeats. In the restricted three-body setting, however, the physical orbit is not merely a two-dimensional textbook curve. The rose-like form is the projection of a dynamical solution, and the vertical component must satisfy its own resonance condition for the full three-dimensional motion to become periodic.

The approximate trajectories serve as initial guesses for a numerical single-shooting procedure. This is a standard but delicate technique in periodic-orbit computation. A trial state—typically including position and velocity—is integrated forward for a proposed period. At the end of that integration, the numerical state is compared with the starting state. If the position and velocity do not match, the initial conditions and, when necessary, the period are adjusted. An iterative correction process then seeks a solution for which the final and initial states coincide within a specified tolerance. In effect, the method solves a boundary-value problem by repeatedly asking the equations of motion to “shoot” from one point and return precisely to it. Using the rose-like approximation as a guide greatly improves the chances of converging on the desired family rather than landing on an unrelated orbit.

The first accurate solutions are computed in the Earth–Moon CRTBP, where the primaries move on circular paths and the rotating frame provides a comparatively stable environment for numerical analysis. Once the 1:n resonant periodic orbits have been found there, Nagai continues them into the ERTBP by gradually increasing the eccentricity of the Earth–Moon orbit. This continuation strategy avoids trying to discover every elliptic solution from scratch. Instead, a known periodic orbit at zero eccentricity is used as the starting point, and the equations are modified in small increments. At each step, the preceding solution supplies the initial estimate for the next one. The procedure traces how the orbit’s shape, period and stability evolve as the idealized circular model becomes more realistic.

Stability is one of the most important questions surrounding any periodic orbit. A trajectory may close perfectly in a mathematical model but be so sensitive to small disturbances that a spacecraft could not remain near it without frequent correction. Nagai analyzes the linear stability of the rose-like periodic orbits during the continuation in eccentricity. In practical terms, linear stability examines how tiny deviations from the reference orbit grow or shrink over time. This information is commonly extracted from the state-transition or monodromy matrix, which maps a small perturbation through one complete period. Its eigenvalues, often called characteristic multipliers, indicate whether perturbations remain bounded, oscillate or expand. The study therefore does not stop at drawing unusual trajectories; it follows their dynamical response as the Earth–Moon model changes.

The work also places these solutions within a long history of three-dimensional periodic orbits in the restricted three-body problem. Earlier studies identified halo orbits, vertical self-resonant satellite orbits and other families that pass near the Earth–Moon libration points. Such trajectories have influenced both theoretical celestial mechanics and mission design, including concepts for spacecraft operating near gravitational balance regions. Rose-like orbits belong to a different visual and dynamical category, but they emerge from the same fundamental principle: nonlinear gravitational systems can support organized families of repeating motion. Their existence illustrates how planar oscillations and vertical resonances can combine to create geometry that is simultaneously simple to recognize and difficult to derive.

The potential significance of the results lies in the bridge they create between classical geometry, modern numerical dynamics and spaceflight applications. The rose curve was developed centuries ago as a mathematical object; here, a related pattern appears naturally in a gravitational model involving the Earth and Moon. That connection could make complicated resonant behavior easier to classify and communicate, while also supplying useful starting points for searches through the enormous catalogue of possible periodic trajectories. The study does not claim that every rose-like orbit is immediately suitable for a mission, nor does it provide operational designs for a spacecraft. Instead, it establishes a computational pathway: approximate the local dynamics, identify resonant structure, refine the orbit in the circular problem, continue it into the elliptic problem and test its stability. As future missions increasingly explore cislunar space, families of structured periodic orbits may become valuable maps of what gravity can make possible—and of where a spacecraft can repeatedly go without simply following an ordinary Keplerian ellipse.

Subject of Research: Resonant rose-like periodic orbits in the Earth–Moon circular and elliptic restricted three-body problems

Article Title: Rose-like periodic orbits in the restricted three-body problem

Article References: Nagai, Y. “Rose-like periodic orbits in the restricted three-body problem.” Celestial Mechanics and Dynamical Astronomy 138, article 52 (2026).

Image Credits: AI Generated

DOI: 10.1007/s10569-026-10327-w

Keywords: Rose curve; periodic orbit; stability; circular restricted three-body problem (CRTBP); elliptic restricted three-body problem (ERTBP)