Trajectories Derived from Periodic Orbits around the Lagrangian Point L 1 and Lunar Swing-Bys: Application in Transfers to Near-Earth Asteroids

: To present a set of trajectories derived from the retrograde periodic orbits around the Lagrangian equilibrium point L 1 , this paper considers the Circular Restricted Three-body Problem with Earth-Moon masses (CR3BP), the Restricted Bicircular, and Full Four-Body Sun-Earth-Moon-spacecraft Problems (BCR4BP and FR4BP, respectively). These periodic orbits are predicted by the dynamics of the CR3BP. To generate the trajectories of this set, ﬁrst, slightly different increments of velocity ( ∆ V s) from those needed to generate periodic orbits around L 1 are applied to a spacecraft in circular low Earth orbits in the same direction of their motion when the Earth, the spacecraft, and the Moon are aligned in this order. Thus, translunar trajectories derived from the periodic orbits are obtained and they will lead the spacecraft to the vicinity of the Moon. Depending on the values of the | ∆ V s|, which are also functions of the relative positioning between the Sun, the Earth, and the Moon, three types of trajectories of interest are found: Collision with the Moon, escape, and geocentric orbits with large semi-major axes. For a well-deﬁned interval of the | ∆ V s|, the trajectories accomplish swing-bys with the Moon and obtain energy to escape from the Earth–Moon system and reach Near-Earth Asteroids (NEAs) between the orbits of Venus and Mars. This procedure reduces the costs of inserting spacecraft into transfer trajectories to a set of NEAs in terms of the required | ∆ V | by up to 5% when compared to Lambert’s problem, for example. This work also presents analyses of examples of transfers to the NEAs 3361 Orpheus, 99942 Apophis, and 65803 Didymos, from 2025 on.


Introduction
According to the Minor Planet Center (MPC), a division of the International Astronomical Union (IAU), more than 28,600 celestial bodies characterized as Near-Earth Asteroids (NEAs) have been discovered by March 2022 [1]. These bodies are part of a larger class of celestial objects called Near-Earth Objects (NEOs), in which comets with orbits close to the Earth's orbit are also included. NEAs are divided into four groups depending on their semimajor axes (a), aphelion distance (Q), and perihelion distance (q). These groups are Atira (a < 1 au and Q < 0.983 au), Aten (a < 1 au and Q > 0.983 au), Apollo (a > 1 au and q < 1.017 au), and Amor (a > 1 au and 1.017 < q < 1.3 au) [2]. Therefore, the asteroids of the Amor group are called "Earth-approaching NEAs" because their orbits are external to the Earth's orbit; the Atira group covers objects with internal orbits; and the asteroids of the mission [16], and the probe returned to the Earth bringing samples on 6 December 2020. The OSIRIS-Rex mission was launched in September 2016 as part of NASA's New Frontiers Program with the goal of collecting samples from the NEA 101,55 Bennu; it arrived at the asteroid in December 2018 and its return to the Earth is expected in 2023 [17]. Recently, on 24 November 2021, NASA's Double Asteroid Redirection Test (DART) mission was launched toward the double NEA 65803 Didymos as the first demonstration of the kinetic impactor technique to change the motion of an NEA [18]; the goal is to impact the Didymos' moonlet in September 2022, when it will be approximately 11 million kilometers from the Earth. The ESA's Hera mission, scheduled for launch in 2024, will follow the same asteroid starting in 2026 to investigate the effects of the impact of the first probe [19,20].
Regarding swing-by maneuvers, they are performed to change the orbits of spacecraft during close encounters with planets or moons. Research on general approaches to this type of maneuver can be found in [21][22][23]. Research on lunar swing-by maneuvers is also present in many works. For example, Dunham and Davis [24] investigated missions with multiple encounters with the Moon, and Prado [25] considered lunar swing-bys to change the inclination of Earth orbits that leave and return to the same altitudes. As examples of missions that performed lunar swing-bys, the ISEE-3, the Nozomi, STEREO probes, and the HGS1 satellite can be highlighted.
NASA's ISEE3 (3rd International Sun-Earth Explorer) was the first artificial satellite placed into a halo orbit around the Sun-Earth L 1 point, in 1978. After completing its original mission in 1981, it was renamed ICE (International Cometary Explorer), and on 10 June 1982, a maneuver removed it from the halo orbit and placed it in a complex trajectory that involved a series of Earth close encounters, two passages through the deep geomagnetic tail (to study it) and the Sun-Earth L2 point, and five lunar fly-bys between March and December 1983. After the last one, on 22 December, the spacecraft obtained enough energy to escape from the Earth-Moon system and reach the comets Giacobini-Zinner on 11 September 1985, and Halley on 28 March 1986 [26,27].
The Japanese Nozomi spacecraft (or Planet B before launch) was launched on 3 July 1998, and after 2 months and 21 days in a phasing orbit (highly elliptical geocentric orbits with perigees of 340 km and apogees of 400,000 km, therefore, with C3 < 0), on 24 September, it performed the first lunar swing-by. After passage through the boundary of the Earth's gravity field, it returned to the Moon and performed a second lunar swing-by on 18 December and, on 20 December, it conducted a final Earth-propelled swing-by to put it on a flight to Mars (in a Trans-Mars Insertion-TMI). A |∆V| of 420 m/s was planned for this passage, but due to the malfunction of a valve, it was 100 m/s lower than expected, requiring an unscheduled correction speed increment, as well as other increments scheduled for TMI correction [28][29][30]. The probe had many problems during its journey, including in the Earth fly-by, and although various efforts were made to recover it, the spacecraft was not able to fulfil its initial mission.
The twin STEREO (Solar Terrestrial Relations Observatory) probes A and B for solar observation were launched on 10 October 2006 by NASA and placed into highly geocentric orbits with perigees of approximately 500 km and apogees just beyond the Moon's orbit (C 3 < 0). During the passage through the fifth apogee (first periselene), on 15 December 2006, they carried out a lunar swing-by when STEREO A was ejected into a heliocentric orbit with a semimajor axis of approximately 0.962 au and eccentricity of approximately 0.00059. STEREO B carried out a second lunar swing-by on 25 February 2007 and was ejected into a heliocentric orbit with semimajor axis and eccentricity around 1.042 au and 0.042, respectively [31]. After that, they successfully performed their missions.
The HGS-1 communication satellite from Hughes Global Services (firstly called AsiaSat 3) was launched on 15 December 1997. After a failure in the Proton rocket during the insertion into its original orbit, it was put into an operational geosynchronous orbit after a sequence of maneuvers that involved two encounters with the Moon on 13 May and 6 June 1998 [32].
It should be highlighted that this work presents a different approach to those used in the missions aforementioned. While all of them made use of multiple passages through the vicinities of the Moon and/or Earth, in this paper, a single lunar swing-by is explored, and they also had different objectives than those proposed here for the escape trajectories derived from the periodic orbits of family G, that is, their application to decrease the total |∆V| required on missions to NEAs.
The presentation of this work is in the following order: In Section, 2 the periodic orbits of family G are presented as proposed by Broucke [3]; in Section 3, the trajectories derived from the orbits of family G are defined and their main features are shown and discussed; in Section 4, a methodology for planning transfer missions among the Earth and the NEAs from the escape trajectories is developed, and the discussion and results of missions to the NEAs 3361 Orpheus, the 99942 Apophis, and 65803 Didymos are presented. Finally, Section 5 is dedicated to the conclusions of the work.

Circular Restricted Three-Body Problem (CR3BP) and Family G of the Periodic Orbits around the Lagrangian Point L 1
The Circular Restricted Three-body Problem (CR3BP) is well known in the literature [3,33]. It describes the motion of a particle with mass m under the action of the gravitational fields of two other so-called primary masses, m P 1 and m P 2 . The particle's mass is negligible, and it does not affect the motion of the primary masses whose orbits are circular and coplanar around their common center of mass, both with the same angular velocity; therefore, they also maintain a constant distance between them. These properties allow the introduction of an orthogonal reference frame (x, y, z) called a synodic frame whose origin is at the common center of mass of the primary masses and in which these bodies are fixed over the x-axis with coordinates (x 1 , y 1 , z 1 ) = −µ P 2 , 0, 0 and (x 2 , y 2 , z 2 ) = µ P 1 , 0, 0 , i.e., x-and y-axes rotate with the same angular velocity as the primary masses relative to an inertial frame (also fixed at the common center of mass). It is useful to consider the normalization of this frame, such as the sum of the primary gravitational mass parameters being equal to 1, that is, µ = µ P 1 + µ P 2 = 1, with µ P 1 = m P 1 / m P 1 + m P 2 and µ P 2 = m P 2 / m P 1 + m P 2 and the distance between the primary masses is also equal to 1 (unit of measurement). In this way, the mean motion and the orbital period of the primaries are equal to 1 and 2π, respectively.
For the CR3BP with an Earth-Moon mass ratio, m P 1 = m Earth = 5.9724 × 10 24 kg and m P 2 = m Moon = 7.3460 × 10 22 kg. Thus, µ P 1 = 0.987849536 and µ P 2 = 0.012150463. The distance between m P 1 and m P 2 corresponds to the average distance between the Earth and the Moon, 384,400 km, and their orbital period is 27.32166 days.
The equations of motion of the CR3BP have five special solutions that define the so-called Lagrangian equilibrium points. Three of them, L 1 , L 2 , and L 3 , are aligned with the primary masses, while the other two, L 4 and L 5 , form two equilateral triangles with the primaries. These points can be seen in Figure 1a.
Broucke [3] numerically examined several families of periodic orbits considering the CR3BP with the Earth-Moon mass ratio. Family G ("G" in capital letter) of retrograde periodic orbits around L 1 is one of them. Family G corresponds to Class "c" of orbits around L 1 in Stromgren's problem with an equal mass ratio [34]; however, the orbits of Class "c" have different evolution, and their symmetry is with respect to both x-and y-axes, while the orbits of family G have symmetry only with respect to the x-axis.
The periodic orbits of family G can be found by numerical integration considering the following initial condition for a particle, at t 0 = 0, in the synodic frame, For instance, by taking x 0 = 0.049999990 and . y 0 = 5.458017660, the periodic orbit shown in Figure 1 is found. In general, periodic orbits of family G can be obtained for 0.012149617 ≤ x 0 ≤ 0.813471672 and −1.950882153 ≤ . y 0 ≤ 601.045380978 [3]. Periodic orbits such as that shown in Figure 1 for (a) synodic (x, y, z) and (b) geocentric (ξ, η, ζ) frames pass remarkably close to the Earth's surface and, after 14 to 15 days, they pass a few dozens of kilometers away from the Moon's surface, so they form a round-trip link between them. In addition, the periodic orbits of family G are sensitive to minor perturbations in their velocities. These features can be explored to generate interplanetary trajectories, as described in the next section.
According to Equation (1), for any periodic orbit of family G, at t 0 = 0, the Earth, spacecraft, and Moon must be aligned in this order, i.e., in inferior conjunction. Therefore, considering the spacecraft in a circular LEO, the spacecraft and the Moon have parallel velocities. In practical terms, this alignment is a weak constraint, because any spacecraft orbiting the Earth in the same plane as the Moon's orbit will achieve this alignment every orbital period, plus a small period of time to compensate for the Moon's displacement. For instance, by taking 0 = 0.049999990 and ̇0= 5.458017660, the periodic orbit shown in Figure 1 is found. In general, periodic orbits of family G can be obtained for 0.012149617 ≤ 0 ≤ 0.813471672 and−1.950882153 ≤ ̇0 ≤ 601.045380978 [3]. Periodic orbits such as that shown in Figure 1 for (a) synodic (x, y, z) and (b) geocentric (ξ,η,ζ) frames pass remarkably close to the Earth's surface and, after 14 to 15 days, they pass a few dozens of kilometers away from the Moon's surface, so they form a round-trip link between them. In addition, the periodic orbits of family G are sensitive to minor perturbations in their velocities. These features can be explored to generate interplanetary trajectories, as described in the next section.
According to Equation (1), for any periodic orbit of family G, at 0 = 0, the Earth, spacecraft, and Moon must be aligned in this order, i.e., in inferior conjunction. Therefore, considering the spacecraft in a circular LEO, the spacecraft and the Moon have parallel velocities. In practical terms, this alignment is a weak constraint, because any spacecraft orbiting the Earth in the same plane as the Moon's orbit will achieve this alignment every orbital period, plus a small period of time to compensate for the Moon's displacement.

Definition of the Trajectories G
This section describes the trajectories derived from the periodic orbits of family G in more detail, hereinafter referred to as trajectories G, or simply TGs. It covers two approaches to define them, as well as other properties.

The CR3BP and a First Approach for the TGs
Considering a spacecraft in circular parking orbits (LEOs) with altitudes between 200 km and 1000 km, the two-dimensional CR3BP, Earth, spacecraft, and Moon alignment at

Definition of the Trajectories G
This section describes the trajectories derived from the periodic orbits of family G in more detail, hereinafter referred to as trajectories G, or simply TGs. It covers two approaches to define them, as well as other properties.

The CR3BP and a First Approach for the TGs
Considering a spacecraft in circular parking orbits (LEOs) with altitudes between 200 km and 1000 km, the two-dimensional CR3BP, Earth, spacecraft, and Moon alignment at t 0 = 0, and the performance of small variations in the speeds (same direction of motion) of the periodic orbits of the family G [3], a first approach for the insertion speeds (V SC ) of spacecraft into trajectories G is obtained. These TGs are tangential to the LEOs and capable of passing between L 1 and the Moon's surface, colliding with the Moon, or escaping from the Earth-Moon system. The values of V SC are given in m/s by where h 0 is the altitude of the circular LEOs 0 ≤ h 0 ≤ 10 6 m , and ϑ 3B values define the point between L 1 and the Moon's surface on the x-axis of the synodic frame the trajectories pass through, and whether they are collision or escape trajectories. For ϑ 3B = −0.976 m/s, the trajectory crosses the x-axis at point L 1 ; for 0 ≤ ϑ 3B ≤ 0.322 m/s, collisions trajectories  (2) can be used to obtain the ∆V TG for any geocentric orbit, for example, from the apogee of an ellipse. To obtain Equation (2), numerical simulations were conducted. In these, small increments in the spacecraft velocities were added to insert them into different kinds of TGs. The characteristic V SC for each kind of TG was then registered, and linear regression with the points found gave form to Equation (2).

The Three-Dimensional Restricted Full Four-Body Problem Sun-Earth-Moon-Particle (FR4BP)
To better understand the properties of the trajectories G and possible applications, the three-dimensional Restricted Full Four-body Sun-Earth-Moon-spacecraft Problem (FR4BP) is considered. In this dynamical system, indexes 1 through 4 are associated with the Sun, Earth, Moon and a spacecraft with negligible mass (m 4 ), respectively. The same normalization adopted for the CR3BP is also adopted for the FR4BP. Thus, , and µ 4 = m 4 /(m 2 + m 3 ) are the gravitational mass parameters of the Sun, the Earth, the Moon and the spacecraft, respectively; and now µ 2 + µ 3 = 1. The normalization is completed by adopting the average distance between the Earth and the Moon (384,400 km) as a unit of measurement. These four bodies move in three-dimensional space under mutual attraction, and the equation of motion of the i-th body in a coordinate frame (X, Y, Z) with the origin at any point in space is ..
where The introduction of the Sun into the dynamical system makes the definition of V SC more complex due to the relative positioning between the four bodies. To describe this, a particular arrangement of the Four-body Problem, in which the Earth's orbit around the Sun and the Moon's orbit around the Earth is circular (BCR4BP), is adopted as a second approach. The semimajor axes (a), eccentricities (e), inclinations (i), and true anomalies (f ) considered are: a Earth = 149,597,870.8 km (1 au), a Moon = 384,400 km, e Earth = e Moon = 0, The longitude of the ascending node (Ω) is undefined for circular orbits with i = 0 • , and the argument of periapsis (ω) is undefined for circular orbits of any inclination, so Ω Earth = ω Earth = ω Moon = 0 • and 0 ≤ Ω Moon < 360 • . However, these values are taken only to provide a step-by-step method to determine V SC for the trajectories G in this second approach. The numerical method that will be presented in the next section allows for adjusting V SC to the current and future values of the orbital elements of the Earth's and the Moon's orbits as well as their instantaneous variations. Figure 2 presents a general illustration of the relative positioning between the Sun, Earth, a spacecraft (SC), and the Moon at launch time (t 0 ) for the BCR4BP, precisely when the spacecraft is in a circular LEO and acquires the insertion speed (V SC ) into a trajectory derived from a periodic orbit of family G relative to the geocentric frame (ξ, η, ζ). V SC is parallel to the Moon's velocity. Note also that the Earth, spacecraft, and Moon are aligned in this order, as described in Section 2, but the Sun alignment is not necessary. Figure 2 presents a general illustration of the relative positioning between the Sun, Earth, a spacecraft ( ), and the Moon at launch time ( 0 ) for the BCR4BP, precisely when the spacecraft is in a circular LEO and acquires the insertion speed ( ) into a trajectory derived from a periodic orbit of family G relative to the geocentric frame (ξ, η, ζ). is parallel to the Moon's velocity. Note also that the Earth, spacecraft, and Moon are aligned in this order, as described in Section 2, but the Sun alignment is not necessary. Considering the same procedure described in Section 3.1 (for the CR3BP), but now for the BCR4BP, and varying fMoon and , a second approach for the insertion speeds (VSC) into TGs can also be found. For the BCR4BP, VSC depends on the h0, fMoon, and . The dependence of VSC with h0 remains linear as in the CR3BP, while and fMoon are approximately sinusoidal with a period of 14.746107 days (3.392 rad), in agreement with the circular dynamics of the system. As a result of this procedure for the BCR4BP, VSC (m/s) can be obtained by where 4 varies continuously, such that: • For 0 ≤ 4 ≤ 0.312 ± 0.008 m/s, TGs that collide with the Moon are found (they are called "trajectories G of collision", TGCs).

•
For 0.312 ± 0.008 < 4 ≤ 0.667 ± 0.021 m/s , TGs that have encounters with the Moon, gain energy, and escape from the Earth-Moon system are found (they are called "trajectories G of direct escape", TGEs).

•
For 0.667 ± 0.021 < 4 ≤ 1.023 ± 0.030 m/s, TGs that approach and have an encounter with the Moon, but do not gain enough energy to escape from the Earth-Moon system, are found (they are called "trajectories G of inversion", TGIs).   Considering the same procedure described in Section 3.1 (for the CR3BP), but now for the BCR4BP, and varying f Moon and Ω Moon , a second approach for the insertion speeds (V SC ) into TGs can also be found. For the BCR4BP, V SC depends on the h 0 , f Moon , and Ω Moon . The dependence of V SC with h 0 remains linear as in the CR3BP, while Ω Moon and f Moon are approximately sinusoidal with a period of 14.746107 days (3.392 rad), in agreement with the circular dynamics of the system. As a result of this procedure for the BCR4BP, V SC (m/s) can be obtained by

The values of
where ϑ 4B varies continuously, such that:

•
For 0 ≤ ϑ 4B ≤ 0.312 ± 0.008 m/s, TGs that collide with the Moon are found (they are called "trajectories G of collision", TGCs). • For 0.312 ± 0.008 < ϑ 4B ≤ 0.667 ± 0.021 m/s, TGs that have encounters with the Moon, gain energy, and escape from the Earth-Moon system are found (they are called "trajectories G of direct escape", TGEs). • For 0.667 ± 0.021 < ϑ 4B ≤ 1.023 ± 0.030 m/s, TGs that approach and have an encounter with the Moon, but do not gain enough energy to escape from the Earth-Moon system, are found (they are called "trajectories G of inversion", TGIs).
The values of A, B, and C have small variations according to the type of trajectory G obtained, i.e., according to ϑ 4B , and they are determined numerically, being A = 2.227 ± 0.074 m/s, B = 2π/3.392 = 1.852442 ± 0.001036, and C = 121.75514 ± 1.412805. For example, in Figure 3a, for which h 0 = 200 km, Ω Moon = f Earth = 0 • at t 0 = 0, and for 0 ≤ f Moon ≤ 360 • , the curve for the first TGC as a function of f Moon is for the first TGE, for the last TGE, Equations (5)-(8) were obtained by the least-squares method, using the same procedure described in Section 3.1 to obtain the simulation data. The standard deviation for Equations (5)-(8) are 0.028 m/s, 0.028 m/s, 0.033 m/s, and 0.033 m/s, respectively.
The dependence on h0 and fMoon are highlighted in Figure 3a,b, in which h0 is equal 200 km and 700 km, respectively, with the curves obtained for = 0° at t0 = 0. The dependence of VSC on fMoon and for h0 = 200 km is shown in Figure 3c.  (4). From the bottom upwards, the first curve is the first TGC, the second curve is the first TGE, the third is the last TGE, and the fourth is the last TGI.
A TGC, TGE, and TGI obtained for a particular positioning defined by = fMoon = 0, at t0 = 0, can be seen in Figure 4 (remember that ℎ = ωEarth = 0). For these values, the Sun, Earth, spacecraft, and Moon are aligned in this order over the X-axis of the coordinate frame with an origin at the Sun's center of mass. Considering the scale used in Figure 4, the three TGs are overlapping in the initial phase, i.e., during the flight between the Earth and the Moon, with subtle differences, as will be described in the next subsection.  (4). From the bottom upwards, the first curve is the first TGC, the second curve is the first TGE, the third is the last TGE, and the fourth is the last TGI.
Equations (5)-(8) were obtained by the least-squares method, using the same procedure described in Section 3.1 to obtain the simulation data. The standard deviation for Equations (5) Figure 3c.
A TGC, TGE, and TGI obtained for a particular positioning defined by Ω Moon = f Moon = 0, at t 0 = 0, can be seen in Figure 4 (remember that Ω Earth = ω Earth = 0). For these values, the Sun, Earth, spacecraft, and Moon are aligned in this order over the X-axis of the coordinate frame with an origin at the Sun's center of mass. Considering the scale used in Figure 4, the three TGs are overlapping in the initial phase, i.e., during the flight between the Earth and the Moon, with subtle differences, as will be described in the next subsection.
In Section 4, the actual ephemerides of the Earth's and the Moon's orbits will be considered for application in missions to NEAs. In this case, the curves for TGCs, TGEs, and TGIs lose their sinusoidal profile and numerical analysis is required. In Section 4, the actual ephemerides of the Earth's and the Moon's orbits will be considered for application in missions to NEAs. In this case, the curves for TGCs, TGEs, and TGIs lose their sinusoidal profile and numerical analysis is required.

Some Features of the TGs in the BCR4BP
As stated previously, TGs are generated by applying ∆ to spacecraft in circular LEOs whose altitudes vary between 200 km and 1000 km. Therefore, , where VSC is given by Equation (4), ℎ = 6371 km is Earth's average radius, and the square root provides the LEOs speeds as a function of the altitude h0 (Figure 5b). Figure 3, in turn, shows the dependence of VSC with fMoon and , and it is possible to note that for fixed values of h0 and fMoon, the difference between the speeds of the first TGC and last TGI is in the order of 1.050 m/s. On the other hand, the sinusoidal dependence of VSC with fMoon for any h0 and shows that the difference between the minimum VSC to generate the TGC and the maximum VSC to generate the last TGI is in the order of 4.6 m/s (Figure 5a). Consequently, the minimum and maximum |∆ | to generate TG will also vary by 4.6 m/s as shown in Figure 5c. The TGs reach apogees of approximately 540,000 km (first TGC) and 550,000 km (last TGI) prior to their passage through the lunar sphere of influence between 8 and 9 days after the launch. The apogees are positioned approximately on the ξ-axis of the geocentric frame as shown in Figure 1b. At the apogees, the speeds of TGs vary between 155 m/s

Some Features of the TGs in the BCR4BP
As stated previously, TGs are generated by applying ∆V TG to spacecraft in circular LEOs whose altitudes vary between 200 km and 1000 km. Therefore, |∆V TG | = V SC − GM Earth /(R Earth + h 0 ), where V SC is given by Equation (4), R Earth = 6371 km is Earth's average radius, and the square root provides the LEOs speeds as a function of the altitude h 0 (Figure 5b). Figure 3, in turn, shows the dependence of V SC with f Moon and Ω Moon , and it is possible to note that for fixed values of h 0 and f Moon , the difference between the speeds of the first TGC and last TGI is in the order of 1.050 m/s. On the other hand, the sinusoidal dependence of V SC with f Moon for any h 0 and Ω Moon shows that the difference between the minimum V SC to generate the TGC and the maximum V SC to generate the last TGI is in the order of 4.6 m/s (Figure 5a). Consequently, the minimum and maximum |∆V TG | to generate TG will also vary by 4.6 m/s as shown in Figure 5c. In Section 4, the actual ephemerides of the Earth's and the Moon's orbits will be considered for application in missions to NEAs. In this case, the curves for TGCs, TGEs, and TGIs lose their sinusoidal profile and numerical analysis is required.

Some Features of the TGs in the BCR4BP
As stated previously, TGs are generated by applying ∆ to spacecraft in circular LEOs whose altitudes vary between 200 km and 1000 km. Therefore, , where VSC is given by Equation (4), ℎ = 6371 km is Earth's average radius, and the square root provides the LEOs speeds as a function of the altitude h0 (Figure 5b). Figure 3, in turn, shows the dependence of VSC with fMoon and , and it is possible to note that for fixed values of h0 and fMoon, the difference between the speeds of the first TGC and last TGI is in the order of 1.050 m/s. On the other hand, the sinusoidal dependence of VSC with fMoon for any h0 and shows that the difference between the minimum VSC to generate the TGC and the maximum VSC to generate the last TGI is in the order of 4.6 m/s (Figure 5a). Consequently, the minimum and maximum |∆ | to generate TG will also vary by 4.6 m/s as shown in Figure 5c. The TGs reach apogees of approximately 540,000 km (first TGC) and 550,000 km (last TGI) prior to their passage through the lunar sphere of influence between 8 and 9 days after the launch. The apogees are positioned approximately on the ξ-axis of the geocentric frame as shown in Figure 1b. At the apogees, the speeds of TGs vary between 155 m/s The TGs reach apogees of approximately 540,000 km (first TGC) and 550,000 km (last TGI) prior to their passage through the lunar sphere of influence between 8 and 9 days after the launch. The apogees are positioned approximately on the ξ-axis of the geocentric frame as shown in Figure 1b. At the apogees, the speeds of TGs vary between 155 m/s (first TGC) and 149 m/s (last TGI). Therefore, the regions near these apogees could be employed for correction maneuver areas.
The TGs reach their periselenes between 14 and 15 days after the launch, and all of them have their specific orbital energy relative to the Earth (ε 24 -following the indices adopted in Section 3.1) increased during the passage through the lunar sphere of influence. However, the same does not occur for their specific orbital energy relative to the Sun (ε 14 ).
To understand the increase in ε 24 , first we consider the angle between the directions defined by a given position vector relative to the Moon at periselene and the straight line joining the centers of the Earth and the Moon (Ψ). According to Broucke [23], the change in the specific orbital energy of a spacecraft's trajectory relative to the Earth (first body) after a swing-by with the Moon (second body), in three-body dynamics, (∆ε 24 ) is ∆ε 24 = ε 24 (t 2 ) − ε 24 (t 1 ) ∝ −sin(Ψ), where ε 24 (t 1 ) and ε 24 (t 2 ) are the specific orbital energies of the trajectory relative to the Earth when it penetrates and leaves the lunar sphere of influence, respectively. Therefore, the largest ∆ε 24 is obtained for Ψ = 270 • . As an example, Figure 6 presents the relations between ∆ε 24 × Ψ, ∆ε 24 × V SC and periselene radius × Ψ for TGEs and TGIs generated for for h 0 = 200 km, Ω Moon = f Moon = f Earth = 0 at t 0 = 0 and the BCR4BP dynamics. It is possible to verify in Figure 6a that all TGEs and TGIs have ∆ε 24 > 0, as expected, and that the largest value of ∆ε 24 corresponds to Ψ = 266 • (≈ 270 • )-the TGCs have been excluded from this analysis. It is also possible to verify, in Figure 6b, that the largest value of ∆ε 24 occurs for the smallest value of V SC capable of generating a TGE; and in Figure 6c, the smallest periselene corresponds to Ψ = 266 • (≈270 • ) and, in comparison with Figure 6a, consequently, it corresponds to the largest value of ∆ε 24 . These results show that the increase in the specific orbital energy relative to the Earth, ε 24 , in the CR4BP dynamics follows, with a slight difference, that predicted by three-body dynamics.
(first TGC) and 149 m/s (last TGI). Therefore, the regions near these apogees could be employed for correction maneuver areas.
The TGs reach their periselenes between 14 and 15 days after the launch, and all of them have their specific orbital energy relative to the Earth (ε24-following the indices adopted in Section 3.1) increased during the passage through the lunar sphere of influence. However, the same does not occur for their specific orbital energy relative to the Sun (ε14).
To understand the increase in ε24, first we consider the angle between the directions defined by a given position vector relative to the Moon at periselene and the straight line joining the centers of the Earth and the Moon (Ψ). According to Broucke [23], the change in the specific orbital energy of a spacecraft's trajectory relative to the Earth (first body) after a swing-by with the Moon (second body), in three-body dynamics, (Δε24) is ∆ 24 = 24 ( 2 ) − 24 ( 1 ) ∝ − (Ψ), where ε24(t1) and ε24(t2) are the specific orbital energies of the trajectory relative to the Earth when it penetrates and leaves the lunar sphere of influence, respectively. Therefore, the largest Δε24 is obtained for Ψ = 270°. As an example, Figure 6 presents the relations between ∆ 24 × Ψ, ∆ 24 × VSC and periselene radius × Ψ for TGEs and TGIs generated for for h0 = 200 km, = fMoon = fEarth = 0 at t0 = 0 and the BCR4BP dynamics. It is possible to verify in Figure 6a that all TGEs and TGIs have Δε24 > 0, as expected, and that the largest value of Δε24 corresponds to Ψ = 266 o (≈ 270 o )-the TGCs have been excluded from this analysis. It is also possible to verify, in Figure 6b, that the largest value of Δε24 occurs for the smallest value of VSC capable of generating a TGE; and in Figure 6c, the smallest periselene corresponds to Ψ = 266 o (≈ 270 o ) and, in comparison with Figure  6a, consequently, it corresponds to the largest value of Δε24. These results show that the increase in the specific orbital energy relative to the Earth, ε24, in the CR4BP dynamics follows, with a slight difference, that predicted by three-body dynamics. Now, in the BCR4BP, the largest changes in the specific orbital energy of the TGs relative to the Sun (Δε14) can occur for more than one value of VSC and Ψ as shown in Figure 7. Δε14 = ε12 − ε14(300), where ε12 is the Earth's specific orbital energy (= −443.298 MJ/kg) taken as a reference, and ε14(300) is the TG's specific orbital energy at t = 300 days, both relative to the Sun. The interval of 300 days was chosen because, from then on, the osculating orbital elements of trajectories G undergo only small variations in general. Still, in the context of the BCR4BP, Δε14 < 0 corresponds to TGEs with perihelia between Earth's and Venus' orbits and aphelia of approximately 1 au; Δε14 ≈ 0 corresponds to TGIs (therefore, co-orbitals to Earth's orbit relative to the Sun), and Δε14 > 0 are TGEs with aphelia between the Earth's and Mars' orbits and perihelia of approximately 1 au. In the example of Figure 7, the smallest Δε14 (Δε14 = −20.46 MJ/kg) occurs for Ψ ≈ 266 o and the smallest VSC Now, in the BCR4BP, the largest changes in the specific orbital energy of the TGs relative to the Sun (∆ε 14 ) can occur for more than one value of V SC and Ψ as shown in Figure 7. ∆ε 14 = ε 12 − ε 14 (300), where ε 12 is the Earth's specific orbital energy (= −443.298 MJ/kg) taken as a reference, and ε 14 (300) is the TG's specific orbital energy at t = 300 days, both relative to the Sun. The interval of 300 days was chosen because, from then on, the osculating orbital elements of trajectories G undergo only small variations in general. Still, in the context of the BCR4BP, ∆ε 14 < 0 corresponds to TGEs with perihelia between Earth's and Venus' orbits and aphelia of approximately 1 au; ∆ε 14 ≈ 0 corresponds to TGIs (therefore, co-orbitals to Earth's orbit relative to the Sun), and ∆ε 14 > 0 are TGEs with aphelia between the Earth's and Mars' orbits and perihelia of approximately 1 au. In the example of Figure 7 value (VSC = 10,932.280 m/s-the smallest VSC that generates a TGE for h0 = 200 km, = fMoon = 0 at t0 = 0); and the largest Δε14 (Δε14 = 23.49 MJ/kg) is obtained for Ψ ≈ 247 o and VSC = 10,932.467 m/s. The first one has the smallest perihelion (0.84 au), near the orbit of Venus, and the largest internal reach, while the second one has the largest aphelion, near the orbit of Mars (1.156 au), and the largest external reach. Thus, Δε14 can generate internal TGEs (with perihelia between Venus' and Earth's orbits and aphelia at 1 au) and external TGEs (with aphelia between Earth's and Mars' orbits and perihelia at 1 au), and the type of TGE depends on and VSC. The dependence on VSC consequently leads to dependence on fMoon and as seen in Equation (4). Figure  8a shows the distribution of perihelia and aphelia for internal TGEs and Figure 8b shows the distribution for external TGEs, for both = 0° at t0 = 0. For this initial condition, the internal TGEs closest to Venus have perihelia and aphelia in the order of 125 × 10 6 km and 150 × 10 6 km, respectively; while the TGEs farthest from Earth toward Mars have perihelia and aphelia radii in the order of 150 × 10 6 km and 182 × 10 6 km, respectively. These distributions are always similar for any value of and h0. That is, the smallest perihelion and the largest aphelion found for the TGEs are in the order of 125 × 10 6 km and 182 × 10 6 km, respectively. The longitude of the ascending node and the argument of periapsis of the TGs depend on the launch date (because of the relative positioning between the Sun, Earth, and Moon), but their inclinations are always less than 2° with respect to the ecliptic in general.
Returning to Figure 5 and comparing the |Δ | and |Δ |, i.e., the speed changes required to generate TGs and conventional escape trajectories, |Δ | is always smaller than |Δ |. |Δ | is 2.5% smaller than |Δ | for h0 = 200 km, and 4.6% smaller for h0 = 1000 km. On the other hand, the heliocentric trajectory generated with a conventional |Δ | will have perihelia in the order of 145 × 10 6 km and aphelia in the order of 151.5 × 10 6 km, much smaller than those achieved by TGEs. In the next section, some examples Thus, ∆ε 14 can generate internal TGEs (with perihelia between Venus' and Earth's orbits and aphelia at 1 au) and external TGEs (with aphelia between Earth's and Mars' orbits and perihelia at 1 au), and the type of TGE depends on Ψ and V SC . The dependence on V SC consequently leads to dependence on f Moon and Ω Moon as seen in Equation (4). Figure 8a shows the distribution of perihelia and aphelia for internal TGEs and Figure 8b shows the distribution for external TGEs, for both Ω Moon = 0 • at t 0 = 0. For this initial condition, the internal TGEs closest to Venus have perihelia and aphelia in the order of 125 × 10 6 km and 150 × 10 6 km, respectively; while the TGEs farthest from Earth toward Mars have perihelia and aphelia radii in the order of 150 × 10 6 km and 182 × 10 6 km, respectively. These distributions are always similar for any value of Ω Moon and h 0 . That is, the smallest perihelion and the largest aphelion found for the TGEs are in the order of 125 × 10 6 km and 182 × 10 6 km, respectively. The longitude of the ascending node and the argument of periapsis of the TGs depend on the launch date (because of the relative positioning between the Sun, Earth, and Moon), but their inclinations are always less than 2 • with respect to the ecliptic in general.
Returning to Figure 5 and comparing the |∆V TG | and |∆V esc |, i.e., the speed changes required to generate TGs and conventional escape trajectories, |∆V TG | is always smaller than |∆V esc |. |∆V TG | is 2.5% smaller than |∆V esc | for h 0 = 200 km, and 4.6% smaller for h 0 = 1000 km. On the other hand, the heliocentric trajectory generated with a conventional |∆V esc | will have perihelia in the order of 145 × 10 6 km and aphelia in the order of 151.5 × 10 6 km, much smaller than those achieved by TGEs. In the next section, some examples for transferring spacecraft to NEAs will be exposed and the |∆V TG | required for these missions is approximately 5% smaller than the |∆V | required for interplanetary trajectories obtained by Lambert's Method. for transferring spacecraft to NEAs will be exposed and the |Δ | required for these missions is approximately 5% smaller than the |Δ | required for interplanetary trajectories obtained by Lambert's Method.

Applications
Studies and results on transfers to the NEAs 3361 Orpheus (Apollo class), 99,942 Apophis (Aten class), and 65803 Didymos (Apollo class) through TGEs are presented in this section. However, a general introduction to the planning of these transfers is required Figure 9. TGEs with the smallest perihelion (internal TGE) and largest aphelion (external TGE).

Applications
Studies and results on transfers to the NEAs 3361 Orpheus (Apollo class), 99,942 Apophis (Aten class), and 65803 Didymos (Apollo class) through TGEs are presented in this section. However, a general introduction to the planning of these transfers is required first. In this section, all numerical simulations consider the Restricted Full Four-body Problem (FR4BP) describes in Section 3.2 as a dynamical system.

Limits of the TGs
In Figure 8, it is possible to identify 124,814,014.0 km and 181,272,593.0 km as the lower and upper limits achieved by TGEs in the first analysis, respectively. From a geometrical standpoint, these values can be interpreted as the radii of two reference circumferences in the plane of the ecliptic (since the inclination of the TGEs is small), one internal with radius R i = 124,814,014.0 km and one external with radius R e = 181,272,593.0 km. The space between them defines the primary interception region of NEAs by TGEs. Figure 10 presents an illustration of this region, and, in broad terms, the analysis of a transfer starts when an asteroid penetrates it. On the other hand, considering the features of the NEAs orbits, those belonging to the Apollo class can cross the two reference circumferences or only the external one, some NEAs belonging to the Aten class can be located completely inside the primary region, while some orbits of Amor class can only cross the external reference circumference and some orbits of Atira class can only cross the internal reference circumference. . An illustration of the primary intercepting region (between the external and internal circumferences), reference ellipse, and launch window for an NEA whose orbit crosses the external circumference.

Parameters for Planning Transfers to NEAs by TGEs
The step-by-step sequence for planning a transfer is as follows:

•
Step 1: Choose an NEA and an interval of interest for planning a mission. Once this choice is made, data on the orbits of the NEA, Earth, and Moon (orbital elements or vectors state) must be obtained through an orbit propagator or, for example, from JPL's Horizons platform for this interval [35].

•
Step 2: Considering the data obtained for the NEA's orbit in step 1, it is necessary to determine if it passes through the primary intercept region, the dates of entry and exit from this region, and the respective position (relative to the Sun), and if its nodes are also in this region. If the nodes are in this region, there is a great probability of finding a TGE able to intercept the chosen NEA. Otherwise, this probability becomes smaller, so even an interception continues to be possible.

•
Step 3. Firstly, considering the case in which the nodes are in the primary intercept region, from the straight lines segments that join these points to the Earth's orbits through the Sun, it is possible to define the semimajor axes of two references ellipses and their orbital periods. Subtracting half of these orbital periods and adding 20 days (time needed for a TGE to reach the Earth's sphere of influence) from the dates corresponding to the passages of the NEA through its ascending or descending nodes, two new dates are found, and they will be the starting point for definition (or study) of the launch windows. Figure 10 shows a reference ellipse considering that the points in which the NEA's orbit entries and exits the intercept region coincide with Figure 10. An illustration of the primary intercepting region (between the external and internal circumferences), reference ellipse, and launch window for an NEA whose orbit crosses the external circumference.
Due to the small inclinations relative to the ecliptic of the TGEs, the interceptions of the NEAs can only occur in a range near their nodes and if they are within the primary interception region. The lower the inclination of the NEA, the greater the possibility of interception.

Parameters for Planning Transfers to NEAs by TGEs
The step-by-step sequence for planning a transfer is as follows:

•
Step 1: Choose an NEA and an interval of interest for planning a mission. Once this choice is made, data on the orbits of the NEA, Earth, and Moon (orbital elements or vectors state) must be obtained through an orbit propagator or, for example, from JPL's Horizons platform for this interval [35].

•
Step 2: Considering the data obtained for the NEA's orbit in step 1, it is necessary to determine if it passes through the primary intercept region, the dates of entry and exit from this region, and the respective position (relative to the Sun), and if its nodes are also in this region. If the nodes are in this region, there is a great probability of finding a TGE able to intercept the chosen NEA. Otherwise, this probability becomes smaller, so even an interception continues to be possible.

•
Step 3. Firstly, considering the case in which the nodes are in the primary intercept region, from the straight lines segments that join these points to the Earth's orbits through the Sun, it is possible to define the semimajor axes of two references ellipses and their orbital periods. Subtracting half of these orbital periods and adding 20 days (time needed for a TGE to reach the Earth's sphere of influence) from the dates corresponding to the passages of the NEA through its ascending or descending nodes, two new dates are found, and they will be the starting point for definition (or study) of the launch windows. Figure 10 shows a reference ellipse considering that the points in which the NEA's orbit entries and exits the intercept region coincide with its nodes. In this example, the semimajor axis is (149,597,870.7 km + 181,272,593.0 km)/2 = 165,435,231.9 km, and half of the orbital period is 212.4 days, plus 20 days, totaling 232.4 days. Then, the date found can be adopted as the beginning of the launch window, since the Earth is within the launch window interval. An interval of three lunar cycles, approximately 82 days, is generally enough to obtain some TGEs able to carry out an Earth-NEA transfer.

•
Step 4. From the new dates found in step 3, intervals containing these dates are defined, and they will be considered the first choice for the launch window. Following this, a set of simulations are conducted to find TGEs able to intercept the NEA chosen. First, the curves V SC × f Moon for these intervals are found; then, from the V SC that defines the TGEs for these intervals, trajectories are propagated to obtain encounters with the NEA.

•
Step 5. Analyses of the results.

Transfer to the 3361 Orpheus
Considering the 2025-2026 biennium, for example, and applying the steps of the previous subsection, the dates in which the Orpheus's orbit enters and exits the primary intercept region are 4 October 2025 and 16 April 2026, respectively. Three transfers to this NEA via TGEs for two different launch windows can be planned. Subtracting half of the orbital period of the reference ellipse (233.1 days) from these dates, the launch window will be between 15 February and 1 July 2025. Table 2 presents the main information about each transfer between 2025 and 2026. In all three transfers, correction maneuvers are required approximately 25 days after the launch, at the aphelion of the TGEs, to cause impact. For the first launch (on 5 March 2025), a distance of 1.5 × 10 6 km between the spacecraft and the NEA at the closest point approach (CPA) must be overcome. A |∆V| of approximately 48 m/s applied at the aphelion is enough to reduce this distance to zero. For the second launch (on 8 May 2025), the distance to be overcome is approximately 1.3 × 10 6 km, and a |∆V| of approximately 18 m/s applied at the aphelion is also enough. For the third launch (on 1 July 2025), the distance between them is 404,000 km at the CPA, and it is overcome with a |∆V| of approximately 12 m/s applied at the aphelion. The ∆Vs for correction maneuvers are calculated to increase or decrease the energy of the TGE (at this point, a heliocentric ellipse) from the vis-viva equation to ensure the NEA-spacecraft distance is null at the encounter. They are applied to the direction of the spacecraft's motion when it is at the aphelion. Other parameters of these missions, such as |∆V| launch, the day of arrival/impact, the relative speed, C 3 at launch, and the angle of impact, among others, are shown in Table 2. It is interesting to note that C 3 relative to the Earth, at departure time (t 0 = 0), is negative for all TGEs. To find the transfer candidates, a search code propagates all TGEs defined by the V SC × f Moon curve (step 4) in the launch window. Figure 11 shows these curves for two launch windows found. A solution is selected from the analysis of three points: (1) The CPA with the NEA, (2) the lowest relative velocity NEA spacecraft, and (3) the shortest transfer time. For the 3361 Orpheus, the three solutions (transfers) found correspond to the closest point approaches, the lowest relative velocities, and the shortest transfer times. For the simulation of transfers to Apophis, for example, a more detailed analysis is required. Figure 12 shows the three transfers (TGEs) in the heliocentric frame (X, Y, Z) with highlights of their launches from the Earth and encounters with the Orpheus.

Transfer to the 99942 Apophis
Orbits of the Aten class spend most of their orbital period, or all of it, inside the primary interception region due to the values of their semimajor axes and eccentricities. As a result, the launch windows for TGEs are almost continuous, and a larger number of solutions can be found for transfers to NEAs belonging to the Aten class via TGEs with CPA below 10 5 km without applying ∆Vs for correction maneuvers. For example, considering a launch window from 03.03.2028 to 06.27.2028 and the steps described in Section 4.2, an analysis of the NEA-spacecraft distance diagrams only makes sense when this quantity is less than 10 5 km, as shown in Figure 13a. Note that there are solutions with NEA-spacecraft distance of less than 10 5 km for the entire launch window interval.    To find the best transfer, the time of flight and the NEA-spacecraft relative speed need to be analyzed. To drive the search, the NEA-spacecraft distance as a function of time of flight, shown in Figure 13b, and as a function of the NEA-spacecraft relative speed, as shown in Figure 13c, are considered. Comparing these figures, the two TGEs with the shortest NEA-spacecraft distance at CPA have the longest transfer times (over 800 days) and relative speeds of the order of 6000 m/s. Excluding these two solutions, the choice falls to the TGE with the fourth-lowest NEA-spacecraft relative speed (5220 m/s), third-lowest NEA-spacecraft distance (35,228.87 km), and the TGE with the lowest NEA-spacecraft relative speed (4373 m/s); however, the latter has an NEA-spacecraft distance of 94,512.43 km. Both are indicated with arrows in Figure 13. Table 3 presents the parameters of these transfers and Figure 14 shows these two TGEs with highlights of their launches from the Earth and encounters with the Apophis. and relative speeds of the order of 6000 m/s. Excluding these two solutions, the choice falls to the TGE with the fourth-lowest NEA-spacecraft relative speed (5220 m/s), thirdlowest NEA-spacecraft distance (35,228.87 km), and the TGE with the lowest NEA-spacecraft relative speed (4373 m/s); however, the latter has an NEA-spacecraft distance of 94,512.43 km. Both are indicated with arrows in Figure 13. Table 3 presents the parameters of these transfers and Figure 14 shows these two TGEs with highlights of their launches from the Earth and encounters with the Apophis.

Transfer to the 65803 Didymos
Didymos will penetrate the primary interception region on 2 November 2026. Subtracting the half orbital period of the reference ellipse, the beginning of the launch window is found on 6 March 2026, and as the passage for this region is fast, the launch window closes on 29 April 2026. Table 4 presents the parameters of the transfer with CPA of 83,904 km, which is overcome with a |ΔV| of 22 m/s at the perihelion of the selected TGE. Figure  15 shows the TGE with highlights of their launches from the Earth and encounters with Didymos.

Transfer to the 65803 Didymos
Didymos will penetrate the primary interception region on 2 November 2026. Subtracting the half orbital period of the reference ellipse, the beginning of the launch window is found on 6 March 2026, and as the passage for this region is fast, the launch window closes on 29 April 2026. Table 4 presents the parameters of the transfer with CPA of 83,904 km, which is overcome with a |∆V| of 22 m/s at the perihelion of the selected TGE. Figure 15 shows the TGE with highlights of their launches from the Earth and encounters with Didymos.

Conclusions
From a set of periodic orbits around the Lagrangian equilibrium point L1, called family G, predicted by the circular Restricted Three-Body Problem, studies on collisions with the Moon, escape from the Earth-Moon system, and far-reaching geocentric trajectories derived from these periodic orbits were presented. Among these, we highlighted the escape trajectories (named TGEs) whose ranges allow for the planning of missions to Near-Earth Asteroids.
One of the main features of the escape trajectories (TGEs) is the departure from the Earth with C3 (characteristic energy) < 0 and, therefore, a |ΔVTG| at launch for altitudes between 200 km and 1000 km is 2.5% to 5% lower than required to generate direct escape geocentric parabolas (C3 = 0) and up to 5% lower than the value required for transfers to NEAs obtained by conventional methods such as Lambert's method. Since they are derived from the periodic orbits of family G, the escape trajectories make a passage through the Moon, collide, or perform swing-bys. In the latter case, the swing-bys guarantee the necessary energy gain to make C3 relative to the Earth positive and send the spacecraft

Conclusions
From a set of periodic orbits around the Lagrangian equilibrium point L 1 , called family G, predicted by the circular Restricted Three-Body Problem, studies on collisions with the Moon, escape from the Earth-Moon system, and far-reaching geocentric trajectories derived from these periodic orbits were presented. Among these, we highlighted the escape trajectories (named TGEs) whose ranges allow for the planning of missions to Near-Earth Asteroids.
One of the main features of the escape trajectories (TGEs) is the departure from the Earth with C 3 (characteristic energy) < 0 and, therefore, a |∆V TG | at launch for altitudes between 200 km and 1000 km is 2.5% to 5% lower than required to generate direct escape geocentric parabolas (C 3 = 0) and up to 5% lower than the value required for transfers to NEAs obtained by conventional methods such as Lambert's method. Since they are derived from the periodic orbits of family G, the escape trajectories make a passage through the Moon, collide, or perform swing-bys. In the latter case, the swing-bys guarantee the necessary energy gain to make C 3 relative to the Earth positive and send the spacecraft into interplanetary space, called TGEs, or to generate geocentric trajectories with an apogee between 9 × 10 5 km and 2 × 10 6 km. The latter are called TGIs.
Studies on launch speeds (V SC ) required to generate collision (TGCs), escape (TGEs), and inversion (TGIs) trajectories and their dependence on the circular LEO altitudes (h 0 ), Ω Moon , and f Moon , as well as the ranges of their final heliocentric orbits, were also presented. These studies led to a methodology for planning transfers to NEAs.
Finally, missions to the NEAs 3361 Orpheus (Apollo class), 99942 Apophis (Aten class), and Didymos (Apollo class) were studied. They show how launch windows are determined and how missions are investigated. The costs of these missions were analyzed in terms of |∆Vs| needed to launch and correct the TGEs until impact with NEAs, which fall below those required to generate parabolic geocentric escape trajectories and interplanetary trajectories by Lambert's method. Thus, transfers planned via TGEs can be classified as low-cost missions.
In forthcoming works, powered lunar and Earth swing-bys from the TGIs could be planned to generate more far-reaching escape trajectories capable of achieving or crossing the orbits of Venus and Mars, for example.