The retrograde geosynchronous orbit, or RGSO, is prized for a specific operational advantage. Because it runs counter to Earth’s rotation, a satellite in this orbit sweeps across the entire geosynchronous belt roughly every 12 hours, a cadence that makes it exceptionally valuable for space surveillance and space situational awareness missions. Yet the very property that makes the orbit useful also makes it costly to attain. The direction of Earth’s rotation, together with the tracking, telemetry, and control constraints of launch sites, means that a direct launch into RGSO from Earth must overcome an additional velocity increment of more than approximately 4 kilometers per second. That figure translates into enormous fuel consumption and significant engineering difficulty. Conventional workarounds rely on lunar gravity assists to move satellites from low Earth orbit or geosynchronous transfer orbit into RGSO, but even those routes demand relatively high velocity increments and offer limited transfer paths, leaving engineers with few good options.

The Beihang team’s answer begins with a careful choice of dynamical model. Rather than the classical circular restricted three-body problem, which treats only Earth and the Moon, the researchers adopted the bicircular restricted four-body problem, or BCR4BP. In this framework, Earth and the Moon revolve in circular orbits around their common barycenter, while the Sun moves in a coplanar circular orbit around the Earth-Moon barycenter. The addition of solar gravity matters enormously for the transfers in question, because some of the most fuel-efficient routes deliberately exploit the Sun’s perturbing influence. The BCR4BP also provides a closer approximation to the true ephemeris than the CR3BP, which makes solutions found within it more likely to survive the transition to high-fidelity numerical models used in actual mission planning.

Three initial orbits anchor the study. The first is a 2:1 stellar-resonant distant retrograde orbit, the nominal orbit once planned for NASA’s Asteroid Redirect Mission. The other two are near-rectilinear halo orbits: an 11:3 synodic-resonant orbit around the L1 point and a 9:2 synodic-resonant orbit around L2, the latter serving as the nominal orbit for the Lunar Gateway space station. These cislunar orbits are natural staging points because spacecraft stationed there already inhabit the sensitive gravitational regime of the three-body problem, where small maneuvers can produce large changes in trajectory. The destination is fixed and demanding: a retrograde geosynchronous orbit at 36,000 kilometers altitude, inclined at 180 degrees, moving retrograde in the equatorial plane.

To design the transfers, the researchers built a two-stage method. For two-impulse transfers, they integrated trajectories forward from the cislunar orbit and constructed an initial guess database using the departure point position, the solar initial phase, and a velocity ratio as search variables. A nonlinear programming stage then optimized the total fuel consumption of the departure and insertion impulses. For three-impulse transfers, the method adds a perilune impulse: trajectories are integrated forward and backward to build pre-perilune and post-perilune Poincaré section sub-databases, which are then matched to divide the transfer into two segments that are optimized separately, yielding a total of three impulses. This database-plus-optimization architecture allows the team to populate an entire solution plane rather than chasing individual trajectories one at a time.

The resulting solution space reveals four transfer types, each reversing the spacecraft’s orbital direction by a different mechanism. Direct transfers keep the apocenter below three times the Earth-Moon distance and achieve orbital reversal through what the authors call apocenter reversal, tracing a semi-elliptical pattern in the time-versus-velocity-increment plane. Weak stability boundary transfers, by contrast, lift the apocenter beyond three Earth-Moon distances into the weak stability boundary region, where sensitive solar gravitational dynamics lower the perigee altitude. These WSB transfers take more than 60 days but cut fuel consumption to approximately 1.3 kilometers per second. Lunar gravity-assisted transfers use a powered flyby near perilune to reverse the orbit, bringing the total velocity increment below 2.2 kilometers per second for transfer times under 60 days. Finally, combined WSB-plus-lunar-gravity-assisted transfers merge both mechanisms and further optimize the Pareto frontier of the time-fuel tradeoff.

One of the study’s clearest findings is the complementary role of the third impulse. Three-impulse solutions fill gaps in the solution space that two-impulse transfers leave empty, effectively rounding out the Pareto frontier of optimal tradeoffs between transfer time and fuel consumption. Across all three initial orbits, the shortest transfers complete within 8 days, and regardless of the starting orbit, the minimum fuel consumption consistently remains below 1.5 kilometers per second. Set against the more than 4 kilometers per second penalty of direct launch from Earth, those numbers represent a dramatic reduction, and they suggest that a spacecraft already stationed in cislunar space could be repositioned into RGSO with a fuel budget that no ground-launched mission could match.

Designing trajectories in an idealized model is one thing; flying them is another. To test engineering feasibility, the team converted the four transfer types to the ephemeris model, using a multiple shooting method to transform the state vector from the Earth-Moon rotating frame to the Earth-centered J2000 inertial frame while simultaneously accounting for the multi-revolution constraints of the initial orbit and the insertion constraints of the quasi-retrograde geosynchronous orbit. This step is where many elegant three-body designs fail, because the real solar system is not perfectly circular or coplanar. The conversion results therefore serve as a critical reality check on the entire approach.

The epoch sensitivity analysis exposes a fundamental tradeoff among the transfer types. Direct transfers, tested under daily initial epochs in January 2025, showed total velocity increment fluctuations of more than 100 meters per second and transfer time variations within 1 day, but feasible results could be obtained for any initial epoch. Lunar gravity-assisted transfers behave similarly, showing low sensitivity to departure timing. Weak stability boundary transfers tell a different story: examined across daily initial epochs in the first two months of 2025, they offer only two consecutive transfer windows per month, with gray regions in the analysis marking failed transfers or cases where the velocity increment exceeded an upper limit of 2 kilometers per month. The combined WSB-plus-lunar-gravity-assisted transfer is the most constrained of all, offering just one window per month with a narrower duration, because it depends on a specific Sun-Earth-Moon angular configuration that occurs only once per synodic period of the Moon.

Taken together, the results confirm that deploying retrograde geosynchronous orbit satellites from cislunar orbits is feasible and provide a new pathway around the prohibitively high fuel consumption of direct launch from Earth. The practical picture that emerges is one of options rather than a single prescription: direct transfers and lunar gravity-assisted transfers offer high engineering flexibility and can depart at nearly any time, while WSB-type transfers deliver the lowest fuel costs but demand patience and precise scheduling around monthly windows. As cislunar infrastructure grows and spacecraft accumulate in distant retrograde orbits and near-rectilinear halo orbits, the ability to redeploy those assets into surveillance-rich retrograde orbits could transform them from fixed-position platforms into a mobile fleet, with the gravitational interplay of Earth, Moon, and Sun doing much of the work that rockets would otherwise have to do.

Subject of Research: Low-energy trajectory design for spacecraft transfers from cislunar orbits to the retrograde geosynchronous orbit

Article Title: Trajectory design for transfers from cislunar orbits to the retrograde geo-synchronous orbit

Article References: Trajectory design for transfers from cislunar orbits to the retrograde geo-synchronous orbit. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: cislunar space, retrograde geosynchronous orbit, distant retrograde orbit, near-rectilinear halo orbit, bicircular restricted four-body problem, weak stability boundary, lunar gravity assist, trajectory design, velocity increment, space situational awareness, three-body problem, Beihang University

The stakes are considerable. Space Domain Awareness, the practice of observing, tracking, and maintaining custody of objects in orbit, has historically focused on the region near Earth, where ground-based radars and telescopes can readily monitor satellites. Cislunar space is a different beast entirely. The region is enormous compared to the near-Earth regime, and objects within it are governed by the competing gravitational pulls of both Earth and the Moon, with the Sun adding further complications. Small errors in a spacecraft’s estimated state can balloon into vastly different predicted trajectories, making it extraordinarily difficult to know where a satellite was, is, or will be.

At the heart of the problem lies an assumption that underpins much of modern space tracking: that uncertainty in a spacecraft’s state can be described by a Gaussian, or bell-curve, probability distribution. Near Earth, this assumption works reasonably well. In cislunar space, however, the dynamics are so nonlinear that Gaussian distributions quickly distort into curved, elongated, and multi-lobed shapes that no single bell curve can capture. When navigation and risk-assessment systems rely on Gaussian assumptions that no longer hold, errors creep into catalog maintenance, collision avoidance, and maneuver detection, potentially compromising the safety of missions such as NASA’s Gateway lunar outpost.

Gupta and DeMars approached this challenge by changing the language in which spacecraft motion is described. Conventional tracking uses Cartesian coordinates, essentially position and velocity in three-dimensional space. The alternative is to use orbital elements, quantities that describe the size, shape, and orientation of an orbit. The researchers built upon the Generalized Equinoctial Orbital Elements introduced by Baù and colleagues, which are characterized by the total energy of the orbiting system. This energy-based formulation allows perturbing forces, whether conservative forces derived from gravitational potentials, such as lunar and solar gravity, or non-conservative forces, to be embedded directly into the equations that evolve the orbital elements, rather than being bolted on afterward through complicated variational equations.

The modification at the core of the new work overcomes a fundamental limitation of the original generalized elements, which are only valid for trajectories with negative total energy. By replacing the semi-major axis with the generalized semi-latus rectum and the mean longitude with the true longitude, the M-GEqOE set remains robust across the wider variety of trajectories found in cislunar space. The six elements include the generalized semi-latus rectum, two parameters describing a generalized eccentricity vector, two parameters orienting the orbital plane, and the true longitude, which varies rapidly along the orbit. Together, they describe a non-osculating ellipse that carries the imprint of the perturbing environment within its very definition.

To test the approach under realistic conditions, the researchers propagated trajectories under a high-fidelity model incorporating the gravitational influences of Earth, Moon, and Sun, using planetary ephemerides from NASA’s Jet Propulsion Laboratory. One subtle numerical obstacle had to be overcome: the generalized angular momentum, which depends on an effective potential energy, must remain nonnegative, yet the third-body gravitational potential can push this quantity below zero even when the underlying physics remains perfectly valid. The team’s solution was to apply a constant offset to the potential, computed once before propagation, which shifts the reference level of potential energy without altering the equations of motion, since those depend only on the gradient of the potential.

The validation covered four representative orbits spanning the cislunar domain. The first was the 9:2 Near Rectilinear Halo Orbit, the baseline orbit for NASA’s Gateway station, where the M-GEqOE propagation over roughly six and a half days matched a conventional Cartesian N-body ephemeris model closely, with errors remaining bounded apart from a modest increase near the close perilune passage. The second was an Elliptical Lunar Frozen Orbit, a highly inclined, eccentric lunar orbit suited to long-term coverage of the lunar south pole and potential data-relay constellations, propagated over thirty days. The third and fourth were Earth-centered sidereal resonant orbits from the 4:1 and 3:1 families, the latter resembling the orbit used by NASA’s Interstellar Boundary Explorer during its extended mission, each requiring roughly twenty-seven days per revolution under combined Earth-Moon-Sun dynamics.

The most striking results concerned uncertainty. The researchers ran Monte Carlo simulations with ten thousand particles for each scenario, propagating identical dynamics in both M-GEqOE and Cartesian coordinates so that any differences in the resulting uncertainty distributions arose purely from the choice of coordinates. To assess Gaussianity rigorously, they employed the Henze-Zirkler test, a statistical test for multivariate normality that evaluates consistency with a Gaussian distribution directly from sample statistics, without requiring an additional density approximation. This choice matters because alternative measures, such as Kullback-Leibler divergence, depend on assumptions about the approximating distribution that can obscure whether observed deviations are intrinsic or artifacts of the approximation itself.

Across all four orbits, the pattern was consistent: uncertainty propagated in M-GEqOE coordinates preserved Gaussian behavior significantly longer than uncertainty propagated in Cartesian coordinates. Along the Near Rectilinear Halo Orbit, the Cartesian representation showed a sharp departure from Gaussianity at perilune passage, where nonlinearities are strongest, while the M-GEqOE representation maintained Gaussian behavior throughout the entire orbit. For the 4:1 resonant orbit, Gaussianity survived until the third perigee crossing at the fourteen-day mark in generalized coordinates, compared with failure at the very next perigee pass in Cartesian coordinates. Even where both representations eventually departed from Gaussianity, as on the 3:1 resonant orbit, the deviations were markedly smaller and recovered more quickly in the generalized elements. Notably, when the researchers visualized the uncertainty clouds in the eigenspace of the covariance matrix, the curvature and tail-like structures betraying non-Gaussian structure were clearly visible in the Cartesian projections, while the M-GEqOE ensembles retained smooth, Gaussian-like shapes.

The implications extend well beyond a single study. Accurate uncertainty propagation is the foundation upon which estimation, tracking, and data association are built, and the authors emphasize that by reducing departures from Gaussianity, the M-GEqOEs may also simplify non-Gaussian uncertainty representations, for example by requiring fewer components in a Gaussian mixture filter. As traffic in cislunar space grows with planned lunar outposts, commercial missions, and scientific observatories, the ability to predict where spacecraft are and how confident we can be in those predictions will determine whether operations in this new frontier remain safe and coordinated. What Gupta and DeMars have shown is that sometimes the most powerful tool for taming a chaotic environment is not a bigger computer or a more elaborate force model, but a smarter choice of coordinates, one that lets the physics of the three-body problem live inside the mathematics from the very start.

Subject of Research: Cislunar spacecraft state and uncertainty propagation using modified generalized equinoctial orbital elements

Article Title: Cislunar state and uncertainty propagation via the modified generalized equinoctial orbital elements

Article References: Gupta, M., & DeMars, K. J. (2026). Cislunar state and uncertainty propagation via the modified generalized equinoctial orbital elements. Celestial Mechanics and Dynamical Astronomy, 138(5), Article 59. https://doi.org/10.1007/s10569-026-10333-y

Image Credits: AI Generated

DOI: 10.1007/s10569-026-10333-y

Keywords: cislunar space, orbital mechanics, uncertainty propagation, equinoctial orbital elements, Space Domain Awareness, spacecraft navigation, Henze-Zirkler test, Monte Carlo simulation, three-body problem, Gateway lunar outpost, celestial mechanics, Gaussian distributions

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