Hayabusa2’s Superior Solar Conjunction Phase

Hayabusa2 is the ongoing JAXA’s sample and return mission to the asteroid Ryugu. In late 2018, Ryugu was in superior solar conjunction with the Earth. It is the first time that a spacecraft experiences the blackouts in the communication link with the Earth while hovering around a small celestial body. In this article, the design of the nominal conjunction trajectory flown by the Hayabusa2’s spacecraft is presented. The requirements for the conjunction trajectory were (1) to guarantee a low fuel consumption, (2) to ensure the visibility of the asteroid by the spacecraft’s wide angle camera (60∘\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$60^{\circ }$\end{document} FoV), and (3) to increase the spacecraft altitude to a safety location (∼109km\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\sim109~\mbox{km}$\end{document}) from the nominal BOX-A operation of 20 km (Home Position - HP). Finally, (4) to return at BOX-A after the conjunction phase. Given the mission constraints, the designed conjunction trajectory appears to have a fish-shape in the Hill coordinates therefore we renamed it as “ayu” (sweetfish in Japanese) trajectory. The optNEAR tool was developed for the guidance (ΔVs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\Delta V\mbox{s}$\end{document} planning) and navigation design of the Hayabusa2’s conjunction mission phase. A preliminary sensitivity analysis in the Hill reference frame proved that the ayu trajectory is a good candidate for the conjunction operation of hovering satellite. The solution in the Hill coordinates is refined in the full-body planetary dynamics (optNEAR Tool) before flight. The ayu conjunction trajectory requires (a) two deterministic ΔVs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\Delta V\mbox{s}$\end{document} at the Conjunction Orbit Insertion (COI) point and at the Home-position Recovery Maneuver (HRM) point respectively. (b) Two stochastic ΔVs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\Delta V\mbox{s}$\end{document}, known as Trajectory Correction Manoeuvres (TCMs), before and after the deep conjunction phase are also required. The constraint linear covariance analysis in the full-body dynamics is here derived and used for the preliminary guidance and navigation planning. The results of the covariance analysis were validated in a nonlinear sense with a Monte Carlo approach which proved the validity of the semi-analytic method for the stochastic ΔVs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\Delta V\mbox{s}$\end{document} planning derived in this paper.


Introduction
Hayabusa2 is the Japanese sample and return robotic mission to the Ryugu asteroid (Watanabe et al. 2019). On 27th of July 2018, the Hayabusa2 spacecraft arrived at Ryugu where the spacecraft-asteroid proximity operations have started. The spacecraft is designed to hover above Ryugu at a distance of 20 km, known as Home Position (HP). In early December 2018, the Hayabusa2 spacecraft was located behind the Sun as seen from the Earth. This condition is known as superior solar conjunction because of the solar radio interference constraints. During the solar conjunction the signal pass through the solar corona that causes a substantial increase in the Doppler data noise with consequent degradation in the orbit determination accuracy (Morley and Budnik 2007). The Earth, the Sun and Ryugu were almost aligned on the 11 th of December 2018 by causing a temporary blackout in the communication link between the Hayabusa2 spacecraft and the Earth. The solar radio interference constraints for the Hayabusa2 spacecraft is the region within 3 • of the Sun-Earth-Probe (SEP) angle called Solar Exclusion Zone (SEZ) (Canalias and Masdemont 2004). The superior solar conjunction event represents a challenge for hovering satellite as communications with the Earth are prohibited. During the Earth-Ryugu communication blackout, the spacecraft needed to be placed in a safe "conjunction" trajectory to prevent undesired: close approaches, escape dynamics or, in the worst case, collision with the asteroid.
(1) For spacecrafts in fly-by with the asteroid, the spacecraft orbit is designed in a heliocentric way and, usually, a manoeuvre is given to increase the Spacecraft-Asteroid relative distance (i.e. Rosetta mission) with respect to the Sun (Morabito et al. 2003).
(2) For a spacecraft in orbit around an asteroid, it is possible to either place the spacecraft in a stable orbit or design an orbit that naturally avoids the SEZ regions (i.e. Terminator Orbits (Takahashi and Scheeres 2019) and Halo Orbits (Farres et al. 2013).
(3) A spacecraft in hovering with the asteroid has the advantage to ease communication with the Earth through a relatively simple communication link geometry that is a fundamental asset for Touch-Down operations in deep space. However, the major disadvantage is that a spacecraft in hovering is located at a highly unstable condition and it requires to be artificially constrained at the HP altitude (of 20 km altitude) with daily maneuvers called BOX-A (x hp = y hp = ±0.5 km, z hp = 20 ± 2.5 km) operation for the Hayabusa2 mission. Therefore, the design of a low fuel expenditure conjunction manoeuvre for a hovering spacecraft is not straightforward as for the case of spacecrafts in (1) fly-by or (2) orbiting the asteroid.
The Hayabusa spacecraft experienced the solar conjunction during the heliocentric transfer phase before its arrival to the asteroid Itokawa, which is a common procedure for deep space missions . It is the first time that a spacecraft experiences the superior solar conjunction while in hovering. Thus, to the best of our knowledge, the solar conjunction phase for a hovering spacecraft has never been studied in the past.
In this article, a low energy transfer trajectory is proposed as a candidate of the conjunction trajectory for the Hayabusa2 spacecraft. The main idea here is to make use of the natural dynamics in the Asteroid-Sun system, knowing that in the Hill problem, the spacecraft motion towards the Sun direction can be forbidden by fixing the S/C energy. This approach is similar to ballistic captures (Qi et al. 2014), weak stability boundaries (Belbruno et al. 2010) or low energy transfers (Koon et al. 2001). JAXA's Hiten mission used similar ideas for the design of a recovery trajectory in the patched Circular Restricted Three-Body Problem (CR3BP) of the Sun-Earth and Earth-Moon system (Uesugi et al. 1991). This low energy trajectory aims at inserting the spacecraft in a trajectory towards the forbidden region of motion (known as curve of zero velocity) located at higher altitude with respect to the HP position. Once the spacecraft had reached the boundary of the curve of zero velocity on the 11 th of December 2018, its velocity was naturally inverted in the opposite direction and the spacecraft returned at the HP point by decreasing its altitude. The conjunction orbit is designed in the Hill frame of the Sun-Asteroid system and due to its fish-shape was named "ayu" (sweetfish in Japanese). The time of flight of this trajectory is around 30 days and it requires two deterministic V s. The first maneuver is given at the Conjunction Orbit Insertion (COI) point (Home Position before the conjunction, 23 rd November 2018) and the second one at the Home-position Recovery Maneuver (HRM) point (Home Position after the conjunction, 29 th December 2018) with a total expenditure in the V less than 1 m/s which implies low fuel consumption. optNEAR (optimum trajectory in Near Earth Asteroid Regions) tool was developed for the trajectory, guidance navigation planning of the ayu conjunction trajectory as the JAXA's software JATOPS (Jaxa Approach Trajectory Optimizer with Statistical Constraints) developed for the Hayabusa2 approach phase could not determine the optimum conjunction solution following the required mission constraints.
The article is organised as follow: Sect. 2 provides an overview of the conjunction phase's mission planning and the relevant definitions of the reference frames used. The equations of motion are given in the Hill coordinates in Sect. 3. Section 4 presents the shooting method for the design of the ayu trajectory developed in the optNEAR tool. Section 5 shows the results of the uncertainty analysis in the deterministic maneuvers at COI and HRM. Section 6 investigates the effect of the navigation in BOX-A on the COI maneuvers while Sect. 7 shows the shooting method for the design of the TCM maneuvers. Finally, the results of the guidance and navigation design are presented in Sect. 8.

Overview of the Hayabusa2's Superior Solar conjunction Phase
For the V planning of the ayu conjunction trajectory, the required interfaces between the Orbit Determination (OD), Flight Dynamics (FD), and the Attitude Orbit Control System (AOCS) teams were prepared. The OD team provides the inputs for the trajectory design to the Flight Dynamics (FD) team as the estimates of the spacecraft's state vector, and the initial knowledge of the covariance matrix (after the navigation). The AOCS team receives the product of the trajectory planning from the FD team in a "way" file V sequence format. Different teams require the definition of the state vector of the spacecraft and of the V in a specific reference system. At least four reference frames are required for the solar conjunction mission phase and their mutual rotations have to be derived. Figure 1 shows a schematic overview of the entire solar conjunction phase for the Hayabusa2 mission and of the four main reference frames used. Two main phases are identified here: the Mission Design and Mission Operation. Two reference frames are required for Mission Design as: (1) The Hill reference frame is centered at Ryugu. In this system, the Sun occupies a fixed relative distance with Ryugu and it is always located in the x hill negative coordinates. This reference frame is rotating along the z hill -axis with an angular velocity equal to the Sun-Ryugu mean motion 1 with respect to the J2000EQ reference frame.
(2) The J2000 equatorial reference frame is here either centered at the Solar System Barycenter (SSB) (J2000EQ) or at Ryugu (J2000EQ-Ry), as shown in the blue dot square in Fig. 1. Note that the J2000EQ reference frame is an inertial reference frame centered in the SSB, while the axes of the J2000EQ-Ry reference frame are parallel to the J2000 EQuatorial (J2000EQ) frame but centered at Ryugu asteroid.
In this article, the Hill reference frame (1) is used for the design of the nominal trajectory as a 1 st guess of the ephemeries dynamical system (2) (N-Body Problem) where the gravity of all planets, the Sun, the Moon and Ryugu are taken into account. The solution in the Hill system is refined with the developed N-Body propagator called optNEAR (optimum trajectory in Near Earth Asteroid Regions). optNEAR is written in the J2000EQ-Ry reference frame and it makes use of the NASA's SPICE Toolkit for loading the ephemeries of the planets. The effect of the Solar Radiation Pressure (SRP) perturbation is included in both the Hill propagator (1) (only cannon-ball model) and the J2000EQ-Ry propagator (2) (flat surface model).
(3) Interfaces between the FD team with the OD, and the AOCS teams require the definition of the Home Position (HP) reference frame as shown in Fig. 1. The HP reference frame is used mainly for the mission operation (red dashed rectangle in Fig. 1). The HP reference frame is asteroid-centred. The z hp -axis is along the Earth-Asteroid line, the Sun-Asteroid line belongs to the positive x hp -z hp plane (gray area in Fig. 1) and the y hp -axis is given such that the HP frame is a right-handed coordinate system. Thus: • The z hp -axis is defined as:r 3 = r Earth |r Earth | ; • The y hp -axis is defined as:r 2 = r Earth × r Sun |r Earth × r Sun | ; 1 In this case, the mean motion is derived from the Ryugu-Asteroid distance at the epoch of deep conjunction (11 th December 2018). • The x hp -axis is defined as: This article is devoted to present the activities of the FD team (blue rectangle Fig. 1) while the results of the conjunction mission operation are give in Soldini et al. (2020). The solar conjunction phase is defined by the relative geometry of the Sun, the Earth and the Hayabusa2 Probe. The Sun-Earth-Probe (SEP) angle reaches its minimum value 2 of 0.459 • in deep conjunction on the 11 th of December 2018 as shown in Fig. 2b. Since the Hayabusa2 spacecraft is usually located at 20 km from Ryugu along the Earth-Ryugu line, the SEP angle can be approximated as the angle between the Earth-Sun line and the Earth-Ryugu line (or Earth-Hayabusa2 probe line) as shown in Fig. 2a.

Solar Conjunction Mission Operation
The solar conjunction operation is divided in three phases as shown in Table 1. When the SEP angle is ≤ 3 • , the radio communication is disturbed by the solar corona (Morley and Budnik 2007); therefore, it is not safe to send commands to the spacecraft when in deep conjunction (blue area in the right panel of Fig. 2b). In Fig. 2b, the SEP angle is function of the epoch between late November/December 2018. The blue area shows the epochs where the SEP angle is below 3 • , when the noise in the Earth-Spacecraft communication signal is affected by the solar corona (Morley and Budnik 2007). Table 1 highlight the three main phases for the solar conjunction operation: (a) Preparation Phase, (b) Deep Conjunction Phase and (c) Recovery Phase as shown also in Fig. 2b. (a) During the preparation phase, the spacecraft performed a 180 • slew maneuver around the z hp -axis to ensure the correct orientation of the 12 thrusters after the deep conjunction phase (flip of the HP reference frame). On the 23 rd of November 2018, the COI maneuver is performed when the SEP angle is 5 • and TCM1 is performed on the 30 th of November 2018, when the SEP angle is 3 • . Four days of observations are planned: two days before deep conjunction and two days right after deep conjunction. As one can see, the two TCMs are scheduled before and after the deep conjunction. However, the two TCMs are not scheduled at the same SEP angle to allow the two days of observation campaign required for the TCM2 maneuver after the deep conjunction phase, Fig. 2b.

Equations of Motion in the Hill Problem of the Sun-Ryugu System
The ayu conjunction trajectory is designed in the Hill Reference frame introduced in Sect. 2. We recall that the dynamics of the mother spacecraft are therefore written in a rotating reference frame where the system is centred in the asteroid Ryugu's center of mass. As mentioned in Sect. 2, the Sun and the asteroid are placed along the x hill -axis at a fixed position. The Sun is in the negative x hill coordinates. The photo-gravitation Hill problem is a special approximation of the photo-gravitational Circular Restricted Three Body Problem (CR3BP) where the mass of the secondary body (i.e. asteroid Ryugu) is considered smaller than the mass of the primary body (i.e. the Sun). Therefore, the Hill problem is a subset of the CR3BP where the Sun and the asteroid have fixed relative position along the x hill -axis. When the perturbation of Solar Radiation Pressure (SRP) acceleration is taken into account, the Hill/CR3BP is known as the photo-gravitational Hill/CR3BP. For a Sun-pointing spacecraft, the SRP acceleration has a constant direction as seen from the Hill reference frame. The Hill problem has two equilibrium solutions know as Libration points, L 1 and L 2 . When the effect of the SRP acceleration is taken into account, the x hill coordinate of the equilibrium points move towards the Sun and the equilibrium points are here re-named as SL 1 and SL 2 (pseudo Libration points). Depending on how the initial energy of the Spacecraft (state vector) is set, two regions of motion can be distinguished where the spacecraft dynamics are or are not permitted (known as forbidden regions). Those ideas have been used to increase the spacecraft altitude from the nominal operational 20 km (HP position) to a safe spacecraftasteroid distance grater than 20 km during the conjunction operation. Note that from now on the subscript " hill " will be removed when referring to the Hill equations to ease the notation in the paper. Therefore, the Hill equations in dimensional coordinates are given by: where, a x is the SRP acceleration for a Sun-pointing spacecraft (cannon-ball model) and it is defined as: where a x0 is the SRP acceleration at 1 Astronomical Unit (AU): In Eq. (4), A is the Hayabusa2 spacecraft's reflective area approximated to 13.276 m 2 , and its mass, m, of 580 kg. P 0 is the solar flux of 1366 W/m 2 , c is the speed of light of 2.99792458 · 10 8 m/s, and C r is the reflectivity property of the Hayabusa2 spacecraft assumed as 1.321. The mean anomaly, n, is defined as: Note that, in our problem, d is the distance of the Sun from the asteroid in the Hill frame at the spacecraft highest altitude location of 1.38 AU (2018/12/11) during the deep conjunction phase. Therefore, the ayu trajectory was designed in the instantaneous Hill reference frame at the 11th of December 2018 epoch rather than considering the average motion of Ryugu around the Sun. This assumption is more appropriate for asteroids with high eccentricity as Ryugu (e = 0.1902). The equations of motion of the Hill problem are Hamiltonian and time independent; thus, they admit an energy integral of motion. 3 For a conservative system (i.e. Sun-pointing spacecraft), the energy integral is a constant of motion and it is given by: For the conjunction epoch at 1.38 AU (2018/12/11), a x was computed as 7.1442 · 10 −11 km s −2 . The equations of equilibrium can be found by setting the velocities and the accelerations of Eq.
(2) equal to zero, and in particular we get: For a x equal to zero (non reflective case) the coordinate of the libration point is simply x L 2,1 = ± 3 μa 3n 2 which is ±89.62 km (radius of the Hill sphere, in Fig. 3a). For a Sun-pointing spacecraft with reflectivity coefficient equal to the Hayabusa2 case (C r = 1.321), the 3rd order equation in Eq. (7) was solved which gives the coordinates of the pseudo libration points as x SL 2 = 21.03 km and x SL 2 = −1,606.78 km (Fig. 3b).

Fig. 3
The left panel shows the dynamical objects (ZVC and equilibrium points) when the effect of SRP is neglected, while the right panel shows the same dynamical object for the Hayabusa2 spacecraft. In both cases, the dynamical objects are computed for the instantaneous Hill problem in deep conjunction (2018/12/11) The green areas in Fig. 3 are regions where the spacecraft motion is not permitted. Conversely, the white area in Fig. 3 represents the region of possible motion for the spacecraft. Therefore, if the energy of the spacecraft is below the energy of either L 1,2 (for C r = 0) or SL 2 (for C r = 1.321), the spacecraft motion can be naturally bounded and constrained around Ryugu. The boundary line between the white and the green areas is called the Zero Velocity Curve (ZVC). As the name suggests, points that lie on the ZVC line have zero velocities. To obtain the ZVC in Fig. 3, the velocities in Eq. (6) are set to zero. Therefore, the critical energy is defined as: The ZVC can be found by setting the z-coordinate to zero in Eq. (8) while x and y are free to vary. The intersection of the potential surface E * with a fixed energy constant results in the ZVC appearance in the x-y plane. The ZVC and its forbidden regions disappear if the energy is above the energy level of SL 1 (E SL 1 ) and the spacecraft motion is permitted for any x-y coordinates. In Fig. 3, the Home Position (HP) point at COI and HRM is (C r = 0, Fig. 3a) or isn't (C r = 1.321, Fig. 3a) inside the Hill sphere depending on the spacecraft's reflectively property. Therefore, the energy of the Hayabusa2 spacecraft at HP point must be above the SL 2 energy. This means that due to the reflectively of the Hayabusa2 spacecraft, its motion is not bounded around Ryugu and V maintenance for station-keeping at the HP location must be given every 1-2 days. Using this notion and adapting it to the superior solar conjunction phase, we make use of the energy integral, E, to confine the Hayabusa2 spacecraft motion along a low energy transfer trajectory (ayu-shape) with a consequent low fuel expenditure.

Low Energy Ayu Conjunction Trajectory Design for Hovering Spacecraft
The ayu conjunction trajectory has a fixed Time of Flight (ToF) and it depends on the selected SEP (Sun-Earth-Probe) angle. The SEP angle is function of the epoch and it was introduced in Fig. 2 (Sect. 2). The epoch associated to the minimum SEP angle is the reference  Figure 4 shows a schematic idea behind the proposed conjunction trajectory. Figure 4a shows the ZVC before the V is given at COI, while Fig. 4b shows how the effect of the maneuver moves the ZVC towards the Sun at an altitude of −H in the x-coordinates of the Hill reference frame. As shown in Fig. 4-5, the major design constraint of the conjunction trajectory requires the spacecraft to depart from HP before deep conjunction (COI's red star in Fig. 5) and to return to HP after deep conjunction (HRM's red star in Fig. 5). It is well know that a trajectory approaching the ZVC or the forbidden regions inverts naturally its motion. Our main idea was to look for a trajectory that can encounter the forbidden region at an altitude with x-coordinate, −H , greater than -20 km (HP). H is an unknown variable of our problem and it is shown in both Fig. 4 and Fig. 5. Once the SEP angle has been fixed, the position at COI and HRM epochs are also fixed as shown in Fig. 4 for 6 • in the SEP angle. Note that the velocity of the spacecraft is artificially kept equal to zero at the HP point (COI or HRM). The HP position keeping is known as Hayabusa2's BOX-A operation that was already introduced in Sect. 2. Knowing that the COI and HRM positions are fixed and lie almost in the x-y plane of the Hill reference frame, the conjunction trajectory was expected to be planar with insertion velocity at COI (v z close to zero). This holds true if the spacecraft is kept at HP point in nominal condition of zero velocity. As a first approximation, a symmetric conjunction trajectory with respect to the x-z plane was expected. The insertion velocity at COI is related to the unknown parameter H . We made use of the ZVC notion to lead our thinking towards the design of the shooting method proposed here.
The altitude −H is chosen to belong to the ZVC, therefore −H is a point of zero velocity by definition as shown in Fig. 4-5. This relates the velocity of the spacecraft directly with the coordinate −H through Eq. (6). However, Eq. (6) gives information of the modulus of the velocity and not of its components. To solve the two boundary problem between COI and HRM (red stars in Fig. 5), we introduce the angle α and a local reference frame centered at COI (x , y , z ) as shown in Fig. 5. For navigation purposes, the Hayabusa2 spacecraft has to be on the day side of the asteroid, therefore the planar COI injection velocity, V xy has to be in the II or III quadrant of the x , y coordinates frame (90 • < α < 270 • ). However, due to the geometry of the problem, we can limit α further by knowing that at −H the velocity of the spacecraft will be inverted. Therefore to reach HRM V xy has to be in the III quadrant (180 • < α < 270 • ) if the SEP angles are between 4 • -6 • . The geometry of the conjunction trajectory can be intuitively understood by imagining the game of pool where a ball strikes a rail and bounces off it (here the "rail" is the ZVC of the forbidden regions).
The main parameters used for designing the trajectory are H , α and v z as shown in blue in Fig. 5. Note that V xy is function of v z . Once the Time of Flight (ToF) (or SEP) and μ a are selected, the trajectory has to satisfy a two boundary value problem where the position at the insertion velocity (COI) and the position at braking velocity (HRM) are fixed. The optimum trajectory can be found as a function of H (the maximum trajectory distance from the center of the Hill reference frame), α (the in plane direction of the insertion velocity) and v z (the z component of the insertion velocity in the Hill reference frame). v z is kept as an optimal parameter as the trajectory is not perfectly symmetric along the x-z plane.
To find the nominal conjunction trajectory, the main steps of the proposed single shooting method (optNEAR tool) are as follow: 1) From Eq. (6), the energy associated to the altitude −H is given as: where the energy, E H , is associated to a point that belongs to the ZVC therefore the state vector is {−H, 0, 0, 0, 0, 0}. 2) Once the altitude of the forbidden region is set (read the energy E H ), we can find the modulus of the insertion velocity at COI as followed by re-organizing Eq. (6): note that the initial position of COI in the Hill reference frame is given by {x 0 , y 0 , z 0 }. This position can be found by knowing the initial state at COI in the previously introduced HP reference frame {0 km, 0 km, 20 km} in Sect. 2 and by applying the transformation from the HP reference frame to the Hill reference frame. The details are given in Appendix A. 3) The first guess of the COI insertion velocity is given by: with V xy = V 2 − v 2 z 0 and V being the module of v 0 obtained from Eq. (10). 4) The initial guess in point 3) is used to solve Eq. (2) through ODE integration. The time of integration corresponds to the selected ToF by assuming the following initial state vector:

5) After the integration at
Step 4), we can find the final state, X f , and compare the expected final position (r H RM ) with the one computed at the end of the ODE integration (r f ). We can thus find H , α and v z that minimize the following cost function: The boundary conditions are given as 180 • < α < 270 • , −0.001 < v z < 0.001 km/s and 80 < H < 800 km. As a first guess, H 0 = 300 km, α 0 = 188 • and v z 0 = 0 km/s were assumed. Figure 6 shows the designed conjunction trajectory for the nominal case (μ a = 32 m 3 s −2 and the SEP angle equal to 5 • ). The presented shooting algorithm was developed in the optNEAR tool. The 3D view of the nominal conjunction trajectory is shown in the right panel of Fig. 6. As one can see, the conjunction trajectory is not perfectly planar because the HP z-coordinate at COI (−168 m in z) and HRM (361.7 m in z) epochs do not lie on the x-y plane. The left panel in Fig. 6 shows the forbidden regions. The conjunction trajectory has a "fish-like" shape that we named as "ayu" (sweetfish in Japanese) trajectory. The light gray area represents the energy of the Hayabusa2 spacecraft in the Hill reference frame when in hovering position (Home Position), while the dark gray is the case of the energy of the spacecraft once the insertion velocity is given at COI (energy increase). The expected maximum distance from Ryugu is of 109 km (2018/12/11) for the nominal case. Table 2 shows the optimum parameters for the design of the ayu conjunction trajectory as a function of the Ryugu's gravity constant, μ a , and the SEP angles (or ToFs). As one can see, larger SEP angles correspond to longer ToF and consequently higher altitude, H , can be reached. Moreover, for a fixed SEP angle, higher μ a results in higher H reached. Therefore, uncertainties in the Ryugu's gravity parameter mean that the spacecraft might be closer or farther from the center of the asteroid. For the conservation of the kinetic and potential energy, the theoretical expected V cost in the hill reference frame to increase and decrease the altitude of the spacecraft from 20 km to 109 km and finally back to 20 km can be computed as: The total V computed with the optNEAR tool's shooting method is of 0.2359 m/s which is consistent with the theoretical expected value. Table 3 shows the energy integral for: the libration points, the H coordinate and the HP position at COI or HRM, where E SL 2 < E COI/H RM < E H < E SL 1 . Figure 7 shows a cartoon of the conjunction trajectory where the spacecrafts attitude is assumed to be kept constant and pointing in the x negative coordinate. The ayu trajectory is almost symmetrical along the x-z axis by keeping the velocities in x and z equal to zero at −H and this configuration allows to keep the asteroid in the FoV of the wide angle camera, ONC-W1 (60 • ). Therefore, the ayu conjunction trajectory demonstrates to respect all the required constraints at low fuel costs and it has been selected as a reference solution for the Hayabusa2's solar conjunction phase.

Robustness to Uncertainties in the Deterministic Maneuvers
The ayu conjunction trajectory meets the major operation constraints. In this section, the robustness of the deterministic maneuvers and the required guidance maneuvers (stochastic V s) are investigated as function of the maneuver uncertainties.   Figure 8 shows the effect of the uncertainties in the COI maneuver in the planar (x-y) case where an error in the planar velocity, V xy , orientation angle, α presented in the in Sect. 4 of the optNEAR tool's shooting method was considered. The reference ayu trajectory for a SEP angle of 5 • and μ a of 32 m 3 /s 2 presented in Table 2 was selected and it is marked in red in Fig. 8. As shown in Table 2, the optimum α is for 187.18 • . An uncertainty of ±5 • in α has been considered and the resulting trajectories are shown in gray in Fig. 8. The modules of the in-plane velocity V xy and out-of plane velocity v z 0 were kept constant. Figure 8b shows that for α < 187.18 • , the spacecraft will be in the positive y-coordinates at the HRM epoch while it will be in the negative y-coordinates for α > 187.18 • .  Figure 9 shows the effect of the uncertainties in α and on the in-plane velocity at COI on to the ToF (Fig. 9a) of the ayu trajectory and the final error at the HRM point (Fig. 9b). The red square represents the nominal solution. The ToF increases with the increase of V xy and it has little dependency on α as shown in Fig. 9a. This can be understood by looking at the energy equation, Eq. (6) where a higher module in the velocity increases the energy of the spacecraft and it can be expected that the altitude H will be farther away from the center of the asteroid. This holds also true for lower module in the velocity at COI where it can be expected a shorter ToF to reach HRM position. The highest error in the HRM position is found for high V xy (i.e. 0.2 m/s) and lower α (< 187.18 • ) as shown in Fig. 9b. This can be explained by Fig. 8b where solutions with α < 187.18 • bring the spacecraft to the positive y-coordinates while HRM is positioned in the negative y-coordinates resulting in position errors of 250 km (5.5 km in log scale) for the worst case scenario Figure 10 shows the effect of the COI maneuver uncertainties on the V at HRM. The solution is function of the uncertainties in 182.18 • < α < 193.18 • and in the module of the in-plane velocity 0.1 < V xy < 0.2 m/s. The red square in Fig. 10 highlights the solution for the nominal ayu trajectory of Table 2.

Uncertainty in the COI Maneuver
The module of the V at HRM is primarily influenced by the in-plane velocity V xy and its dependency with α shows that for a given V xy there exist a optimum α angle that minimizes the V at HRM. It is possible to notice that the effect of uncertainties in α and in V xy results in higher V at HRM in the x direction as shown in Fig. 11. From Fig. 11, it is possible to conclude that α and V xy do not affect much the y and z components of the V at HRM.

Uncertainty in the HRM Maneuver
This section investigates the effect of uncertainties in the HRM maneuver at the end of the conjunction phase. Figure 12 shows a schematic representation of the ayu conjunction trajectory at its final arrival at the HRM point (20 km altitude from Ryugu on the 29th of November 2018). Two parameters of uncertainties were considered, the in-plane velocity V xy and its orientation angle, α. Figure 13a shows the map of the Close Approach (CA) distance from Ryugu in km after 1 day of integration time from the HRM location as function of −22 • < α < 22 • and the in-plane velocity 0 < V xy < 40 cm/s. This is the contingency scenario in which we considered the robustness of the ayu conjunction trajectory against failures in executing the HRM maneuver at the HRM point allowing a 1 day margin for rescheduling the HRM maneuver. In Fig. 13, three solutions for V xy are highlighted in white in Fig. 13a and in black in Fig. 13b. The nominal ayu trajectory case with V xy equal to 12.11 cm/s in dashed line, the case of 20 cm/s in solid line and the case of 30 cm/s in dot lines.
For uncertainties in V xy less than 20 cm/s, the spacecraft will not impact on Ryugu after 1 day of simulation starting from the HRM position. Conversely, for V xy > 20 cm/s, the spacecraft can potentially impact on Ryugu when 0 • < α < 9 • (white area of Fig. 13a). Moreover, arrival velocities, V xy , at HRM below 18 cm/s are not of much concern if the HRM maneuver is missed because the CA distance has little dependency from α. Therefore, the solution of the nominal ayu trajectory for V xy = 12 cm/s (white dashed line in Fig. 13a) allows recovery maneuvers for any uncertainties in α. The red square marks the nominal ayu trajectory for α = −7.2 • and V xy = 12.11 cm/s at HRM. Figure 13b shows the dependency of CA for a fixed V xy in this case 12.11 cm/s which correspond to the white dashed line solution in the CA map of Fig. 13a). As one can see, the spacecraft maintains a distance from Ryugu between 9-11 km allowing to reschedule a recovery maneuver to reach the desired 20 km altitude in the BOX-A operation.

Contingency Study: Recovery Trajectory
A contingency study is here carried out to demonstrate the robustness of the selected ayu conjunction trajectory. Figure 14 shows a schematic representation of the possible contingency scenario. Following the result of Sect. 5.2, this case scenario study a possible recovery trajectory when the V maneuver at HRM on 29th of December 2018 is cancelled. The proposed strategy assumes that a recovery V will be given a  Recovery operation: the black trajectory is the designed ayu trajectory between COI (green star) and HRM (black start). The recovery maneuver is executed 1 day after the HRM point. The red star is the HP position after 1 month from the recovery position ordinates after the conjunction. Figure 15 shows the optimum recovery trajectory from the 30th of December 2018 until the 29th of January 2019 (20 km altitude from Ryugu). The optNEAR's shooting method presented in Sect. 4 was modified in point 3) where the recovery trajectory starts from the recovery location rather than at COI and the epochs have been changed accordingly. However, the core idea of the presented method holds true for the recovery case. Figure 15a shows the forbidden region (red area) associated to the recovery trajectory (red line). The energy of the recovery trajectory (forbidden region's red area in Fig. 15a) is lower than the energy of the ayu trajectory (forbidden region's dark gray area in Fig. 6). Finally, Fig. 15b shows the ayu and the recovery trajectories as seen in the Hill reference frame. The boundary conditions for the recovery trajectory are: 188 • < α < 195 • , −0.001 < v z < 0.001, and 40 < H < 90 km. The initial guess in the optimum parameters were set as: H 0 = 70 km, α = 190 • and v z 0 = 0 km/s. The optimum parameters after the optimisation method are given by:ᾱ = 188.36 • ,H = 83.17 km andv z = 0.81 mm/s.

Robustness to Uncertainties in the Navigation at Home Position (HP-NAV)
Once the robustness of the ayu trajectory has been verified against uncertainties in the deterministic maneuvers (at COI and HRM), we are now interested to verify the effect of the navigation uncertainties at HP position on the COI maneuver. For any Hayabusa2 proximity operations, the spacecraft is kept within the nominal operation box named A (BOX-A). BOX-A is centered at the HP altitude around 20 km in the HP reference frame introduced in Sect. 2. When the spacecraft is in BOX-A, the navigation requires a V maneuver every 1/2 days to keep the spacecraft inside the BOX-A (know as HP-NAV or HP maintenance). Therefore, the COI operation has to account for the effects of the navigation in BOX-A. For the uncertainty analysis, we considered a uniform distribution in the position uncertainties (Fig. 16a). The spacecraft can be potentially located anywhere within the BOX-A. A Gaussian distribution was selected for the uncertainty in the COI's insertion velocity (Fig. 16b). The nominal COI position is shown in red in Fig. 16. The following parameters have been considered: • Uniform distribution in position within the BOX-A: x H P = y H P = ±0.5 km z H P = 20 ± 2.5 km • Gaussian distribution in velocity in the HP reference frame: 3σ = [0.5, 0.7, 1] cm/s, for example see Fig. 16b We kept the same assumptions for the SEP angle and μ a as shown in Table 2. Note that the coordinates are given in the HP reference frame and they have to be expressed in the Hill reference frame through the transformation given in Appendix A. The uncertainty analysis is carried out in the Hill problem. We aim to verify if correction maneuvers (stochastic V s) between COI and HRM are required.

Uncertainty in HP-NAV Position
For the HP-NAV uncertainty in position, a uniform distribution of 3000 points within the BOX-A (green box in Fig. 16) was considered. The velocity at COI was kept to its nominal value. Figure 17 shows the final position at HRM in the Hill reference frame, starting from initial states at COI as shown in Fig. 16a. In Figs. 18, Ryugu gravity's is set to 32 m 3 /s 2 . In Fig. 18, figures that belong to the same column as Fig. 18.a, Fig. 18.d and Fig. 18.g show the results for increasing values in the SEP angles from 4 • to 6 • . Rows of figures in Fig. 18 show the same SEP angle solution but with different views for example for a SEP angle of  . 18.a), a planar view in the x-y plane (Fig. 18.b) and a planar view in the x-z plane (Fig. 31.c). As one can notice, if the spacecraft is initially inside the BOX-A at COI, it will be inside the final box at HRM. This suggests that the conjunction trajectory is robust in error in the HP-NAV position. The same analysis was carried out for Ryugu's gravity equal to 11 m 3 /s 2 (nominal case shown in Fig. 31, Appendix B) and  (Fig. 32, Appendix B). From this analysis, we concluded that the corners of the BOX-A are the only representative cases to perform the uncertainty analysis in velocity. Therefore, instead of using the full set of 3000 solutions, only 9 points (the 8 corners of BOX-A and the nominal solution, HP position at COI) will be considered. Figure 17 shows a 3D view of the manifold (blue trajectories) associate to BOX-A when Ryugu's gravity is 32 m 3 /s 2 . Further solutions can be found in Appendix B for μ a = 11 and 92 m 3 /s 2 .

Uncertainty in HP-NAV Velocity
In this section, the effect of velocity uncertainty in the HP-NAV at COI are investigated. The nominal case of 3σ = 5 mm/s is presented here (for the other solutions refer to Appendix C). Based on this results and the one in Appendix C we concluded that at least one stochastic maneuver (TCM) right after the deep conjunction phase is required. We also concluded that error in the execution V s during the mission operation below or equal 4 mm/s can be accepted. Therefore, error in the main V s above 4 mm/s will require a contingency correction V (trim V s).

Design of the Trajectory Correction Maneuver (TCM)
The Trajectory Correction Manoeuvre (TCM) is designed by solving a two boundary value problem where the initial position at TCM and the final position at HRM are fixed. From the  previous analysis, we concluded that at least one TCM is required after deep conjunction (SEP ∼ 0 • ). Once the spacecraft reaches its highest altitude along the ayu trajectory, a TCM is scheduled at SEP = 4 • as shown in Fig. 22. The ToF is fixed and is the time frame between  (2) is done and the error between the final desired position and the expected position is computed. Finally, the insertion velocity at TCM is found through the minimisation of the following Note that the V at TCM is the difference between the TCM insertion velocity and the velocity of the spacecraft before the TCM manoeuvre. We adopted a conservative approach so that a second TCM maneuver is scheduled before the spacecraft enters in deep conjunction.
As part of the guidance in the conjunction phase, we have to compute two TCMs (stochastic V s) as mentioned in Sect. 2.

Preliminary Guidance and Navigation Design: Consider Parameters Covariance Analysis
A complete and exact knowledge of the spacecraft state vector is generally not possible. The individual state variables cannot be measured precisely and available measurements are usually function of the state variables. The spacecraft motion can be affected by nonmodelled forces which are not adequately represented in the dynamical model. Its model parameters can be uncertain too. By definition, the linearized system of equations of the nonlinear dynamics is a further approximation. These sources of error make the knowledge of the spacecraft's state uncertain. The computation of the most likely current state of the spacecraft in the presence of measurements and model uncertainty is the focus of orbit determination which will not be studied in this paper. The error analysis instead involves an investigation of the impact of various sources of error on the orbit determination. The output of an error analysis, as done in this paper, provides the magnitudes of the state vector variances and covariances, thus quantifying the relative contribution of the significant error sources for a preliminary guidance and navigation design (Gordon 1991).
In the first part of this article, the design of the ayu conjunction trajectory has been presented. The second part of this article is devoted to investigate the preliminary OD error analysis for this newly proposed trajectory. The technique used here involves the linearisation of the nonlinear equations of motion around a reference solution (nominal ayu path) and then applying linear estimation techniques (Tapley et al. 2004). The OD process is thus changed from estimating the state of the linear, time-varying, deviations from the reference trajectory. The reference solution used here is generated by numerical integration of a highfidelity nonlinear equations of motion in the optNEAR tool. As OD is out of the scope of this analysis, the current work focuses only on the consider parameters covariance analysis that includes the uncertainties in the dynamical model and in the measurements (Montenbruck  Gordon 1991). The scope of the covariance analysis is to compute the covariance matrix after some tracking and to estimate the stochastic V s at TCMs. In this article, the linear covariance analysis is carried out to account for uncertainties in the HP-NAV (initial knowledge error) at the COI maneuver (20 km from the asteroid), the GCP-NAV or OpNav (two-way Doppler and ONC-T or ONC-W1 cameras) and Vmaneuver execution error. The covariance analysis is carried out considering three trajectory legs: from COI to TCM1 (Leg 1), from TCM1 to TCM2 (Leg 2) and from TCM2 to HRM (Leg 3). The uncertainties in the navigation are considered independently for only Leg 1 and Leg 2. The epochs of the maneouvers for the reference conjunction trajectory are shown in Table 1.

Initial Knowledge at COI
Before the insertion V is given at COI, the HP-NAV navigation error was taken into account through the initial knowledge covariance matrix, P 0, hp , given in the HP coordinates as shown in Table 6. Figure 23a shows the three independent V H P along the X H P , Y H P and Z H P coordinate of the HP reference frame. In Fig. 23b, a local reference frame aligned with the V direction is considered (Wagner and Goodson 2008). This V H P is given by transforming the

V -Maneuver Execution Error
nominal V J 2000EQ from the J2000EQ reference frame to the HP reference frame. The V execution error uncertainties are given in the module and in the pointing direction. Three independent covariance matrices for the X H P , Y H P , Z H P components were considered. The assumption made in the pointing accuracy and in the velocity module are shown in Table 7. The process noise covariance matrix is given by: where P (3×3) ExErr, hp is derived from the assumption made in Fig. 23.

Navigation Error: Consider Parameters Covariance Matrix
The consider parameters covariance analysis is a batch weighted least square formulation that also includes parameters uncertainty. For an unbiased system with a sufficient number of observations, the covariance matrix will asymptotically approach zero and the estimation algorithm will be insensitive to any further observations. This condition can lead to filter divergence. A possible approach to prevent this divergence is to recognize that the linearised equations are in error. To compensate for this error, we assumed that the error in the linearised equations is due to uncertain parameters called considered parameters. They are a priori estimates and the associated covariance matrix is known. Therefore, the nonlinear dynamical model in J2000EQ coordinates is given by (Tapley et al. 2004): The first equation in the system of Eq. (17) is a compact form of the N-Body dynamics written in J2000EQ-Ry coordinates, for further details in the propagator developed in optNAER tool refer to Appendix D. C is the vector of the considered parameters. Note that the vector C can be divided in the dynamical model considered parameters, C d , and measurements model parameters, C m : By expanding Eq. (17) in a first order Taylor series around the nominal ayu trajectory, Eq. (17) turns into:

ẋ(t) = A(t)x(t) + B(t)c(t)
where c(t) = C(t) − C * (t) is the residual of the consider parameters. C is the nominal vector of the consider parameters and C * is the reference vector.
The matrix A is the matrix of the linearised equations. Note that ∂F ∂Cm * = 0 as the measurements model parameters do not appear in the dynamical model. In order to compute the observation at a time t i different from the observation epoch, the state transition matrix has to be computed. The state transition matrices associated to the state and the consider parameters are numerically computed through: with (t 0 , t 0 ) = 0 and (t 0 , t 0 ) = I . Note that: with ∂X(t) ∂Cm(t 0 ) = 0. The sensitivity matrix for both the state vector and the consider parameters are defined as: We make use of the state transition matrices to compute the observations at time t i so that: The sensitivity matrix are finally computed as The consider parameters filter is given by solving the following set of equations:

Observations and the Sensitivity Matrix (H )
In our analysis, the measurements in the range and the range rate are assumed to be along the Earth-Spacecraft line. The azimuth and elevation angles are measured from the spacecraft (Gordon 1991) (−r) as shown in Fig. 24. The observations vector, y, is defined as:

(x−xe)(vx−vxe)+(y−ye)(vy−vye)+(z−ze)(vz−vze)
√ (x−xe) 2 +(y−ye) 2 +(z−ze) 2 tan −1 y ONC-W1/T Case 1  with r = x 2 + y 2 + z 2 (spacecraft-asteroid distance). x e , y e and z e are the coordinates of the Earth 4 as seen from Ryugu in the J2000EQ reference frame. We can now compute the sensitivity matrix following the definition given in Eq. (23): (28) H is the measurement matrix and it is a time-varying matrix function of Eq. (27). The expected measurements are given in Table 8 and Table 9 for Leg 1 and Leg 2 respectively where we assumed a 2-days observation campaign before each TCM. The assumed Gaussian distribution (1σ ) is given in the list below:  ONC-W1/T Case 1

Consider Parameter Covariance Matrix (B)
We have chosen the spacecraft reflectively coefficient, Cr, and the camera accuracy, α and δ as consider parameters. Their standard deviations are given as follow: • 1σ Azimuth = Elevation (ONC-W1) [rad]: 0.00102265 (as before) • 1σ Azimuth = Elevation (ONC-T) [rad]: 0.000102265 (as before) • 3σ Cr: 10% The consider parameters matrix B is derived through (Montenbruck and Gill 2005;Tapley et al. 2004): The matrix B can be easily partitioned into two contributions. The matrix B d is associated to the derivatives with respect of the dynamical model parameters, C d , and the matrix B m is associated to the derivatives with respect of the measurements model. Since the measurements do not appear in the dynamical equations, we can conclude that B m = 0. The partial derivatives of F 4 , F 5 and F 6 (components of F in Eq. (17)) with respect to Cr are derived from the solar radiation pressure acceleration and for a flat plate surface the derivative of the acceleration is: withn being the normal vector to the flat plate, pointing towards the Earth: and θ being the angle betweenn and the Sun-spacecraft direction r ls : cos θ = r ls · r Earth r ls r Earth (32) is the reflectively parameter and K is a constant function of the spacecraft's property (mass/size) and Sun's luminosity. Further details in the definition of the SRP acceleration can be found in Appendix D. By knowing that ∂F 1 ∂Cr = ∂F 2 ∂Cr = ∂F 3 ∂Cr = 0, the matrix B can be rewritten as:

Sensitivity Matrix of the Consider Parameter Covariance Matrix (H ci )
The sensitivity matrix of the consider parameters is derived by linearizing the nonlinear measurement equations. From the definition given in Eq. (23), the sensitivity matrix of the consider parameter covariance matrix is given as:

Consider Parameters Covariance Matrix (P c ): Least Square Batch Filter
The consider parameters covariance matrix, P c , is larger than the noise-only covariance, P x . By adopting a large number of measurements, the impact of the data noise is effectively averaged out and the uncertainty of the estimated parameters is decreased. The consider parameters are assumed to be constant throughout a single orbit determination but affected by a given uncertainty. The matrix P c does not decrease with the increasing in the data rate, but it is essentially constant for a given data arc and tracking configuration. The consider covariance batch filter is suited for assessing the impact of systematic errors in the orbit determination process. For the case of the consider parameters, the covariance matrix is computed as (Montenbruck and Gill 2005): The noise-only covariance matrix, P x , is given by: where the weight matrix, W i , is defined as: We can now derive the following matrices as: and The inverse of the weight matrix associated to the constraint parameters is provided by: Note that, the covariance matrix was computed in J2000EQ-Ry but it can be transformed in the HP reference frame.

Covariance Analysis
The covariance matrix along the reference conjunction trajectory can be computed as a function of time as follow: with (t i , t i ) = I . Note that the state transition matrix is used to update the time of the covariance matrix from t i to t i+1 . We also distinguished between the time before and after the maneuver where the effect of the process noise ( V execution error) is given by P n that was previously introduced in Eq. (16) and here transformed in J2000EQ-Ry coordinates. The effect of the initial knowledge covariance matrix is included in the termP (t i ) while the effect of the navigation at time t i+1 is included in the matrix P c (t i+1 ). Note that the effect of the navigation at t i is given in the termP (t i ). Once the covariances before and after the V maneuver at COI, TCM1, TCM2 and HRM are given, it is possible to compute the stochastic V (guidance). Figure 25 shows the standard deviation expressed in the HP frame and projected to the final HRM epoch for both the positions and velocities. The dashed gray vertical lines highlight the TCM1 (-29 days to HRM) and the TCM2 (-4 days to HRM) epochs. The red dots show the effect of the execution error V . As one can see, the effect of the navigation before TCM1 (Leg 1) and TCM2 (Leg 2) aims to reduce the standard deviation error at HRM.

Semi-analytic Stochastic V Method
For the statistical analysis of control maneuvers, we follow the idea presented in Renaultand and Scheeres (2003). A similar approach has been used for the design of the TCMs during the Hayabusa2's approach phase (Tsuda et al. 2013). The stochastic V can be computed analytically following the definition of the state transition matrix. At time t i , a manoeuvre V i is executed to null the error in position at time t H RM , and a manoeuvre V H RM is scheduled at time t H RM to null the error in velocity  Fig. 26. Therefore, we can write the following formulas: so that in a compact form we have: The covariance matrix of the stochastic V is given by: and we recognise that E[ V V T ] is the covariance matrix in Eq. (41): In Eq. (48), t i is the time of the maneuver (e.g. COI, TCM1 and TCM2) and the covariance matrixP i (t i ) is computed as in E. (41). The diagonal of the covariance matrix provides the standard deviation of the stochastic V : and we can now separate the covariance matrix of the stochastic V at time t i and the one at t H RM so that we can write in a compact way: We are interested only in the standard deviation at time t i , P (3×3) V i,J 2000 (t i ). Figure 27 shows the results of the covariance analysis for the deterministic and the stochastic V s in each components of the HP reference frame. Note that due to the mounting direction of the thrusters attached to the body of the spacecraft, we have to account for a thrust loss factor of cos(75 • ) −1 in V s given along the y hp coordinate. Given the linear covariance analysis, Fig. 28 shows an example of estimated trajectory (in red) vs the ayu nominal trajectory (in blue).

Monte Carlo Analysis in the Nonlinear System
In this section, we verified that the covariance matrix propagated from COI to TCM2 (no navigation) returns the same stochastic V in both the previously presented semi-analytic method ( Fig. 29.b) and the Monte Carlo method for the nonlinear system (Gordon 1991) ( Fig. 29.a). For the Monte Carlo analysis, we performed a simple shooting between TCM2 and HRM to determine the standard deviation of the V at TCM2 under the effect of nonlinearities. Table 10 shows the results of the two analysis compared. Both methods give similar solutions showing that the semi-analytic method is reliable under the effect of nonlinearities. We also demonstrated that our shooting method (optNEAR) can control the spacecraft within a 99% confidence criterion as shown in Table 10 and in Fig. 29.a (red trajectories).

Conclusions
In this article, we have shown the design of the ayu conjunction trajectory for hovering satellite around minor bodies. The case of the Hayabusa2 mission trajectory planning was presented. The ayu trajectory is a low energy transfer trajectory designed in the Hill reference frame. It requires two deterministic V s (COI/HRM) with overall fuel cost of less than 1 m/s. Uncertainties in the insertion V at COI affects the overall V x at HRM. Under contingency scenario at HRM, a recovery trajectory can be rescheduled after one day from the missed HRM maneuver. The effect of the BOX-A's initial knowledge navigation error at COI requires at least one TCM maneuver after deep conjunction (SEP angle ∼ 0 • ). To The effect of uncertainties in the observations and dynamical model are investigated and the stochastic V s are computed using a semi-analytic formulation. When compared to the nonlinear Monte Carlo method, we concluded that the semi-analytic formula works well. Moreover, we are confident that the shooting method developed in optNEAR works well for the V planning during the superior solar conjunction mission operation. The rotation matrix C takes coordinates in the Hill reference frame (1) (in Sect. 2) to coordinates in J2000EQ-Ry reference frame (2) (in Sect. 2): where 1)ĉ 1 = ra |ra | ; 2)ĉ 3 = h |h| = ra ×va |ra ×va | ; 3)ĉ 2 =ˆc 3 ×ĉ 1 |ĉ 3 ×ĉ 1 | , where r a and v a are the position and the velocity of the asteroid with respect to the Sun. 5 We assumed thatĊ = 0. Figure 30 shows the ayu trajectory in the HP reference frame. As the ayu trajectory is designed to departure from HP position (20 km altitude in the z-hp coordinate) at COI and return at HP at HRM, the ayu trajectory appears as a periodic orbit in the HP reference frame.
S. Soldini et al.  Figure 31 shows that all the solutions inside BOX-A at COI are bounded inside the box at HRM therefore, it makes sense to limit the analysis to the corner solutions of BOX-A for μ a = 11 m 3 /s 2 . Fig. 32 shows that all the solutions inside BOX-A at COI are bounded inside the box at HRM therefore, it makes sense to limit the analysis to the corner solutions of BOX-A for μ a = 92 m 3 /s 2 . This behaviour seems to hold true for all the values in μ a analysed. Figure 33 shows the Gaussian distribution in position at HRM for 3σ = 7 mm/s of uncertainty at COI while Fig. 34 shows the distribution of the arrival velocity at HRM. Tables 11  and 12 show the 1σ error in both the final position (Fig. 33) and velocity (Fig. 34) at HRM for all the six cases. The maximum errors are in the Hill x-coordinate for solution c in Fig.  33-34 with a 3σ error of 16 km in position and 4.75 cm/s. Initial navigation errors at around 7 mm/s requires at least one TCM after the deep conjunction point to compensate of the large errors in position and velocity at HRM.

Appendix C: Uncertainty Analysis in Velocity
Initial navigation errors at around 1 cm/s requires at least one TCM after the deep conjunction point to compensate of the large errors in position and velocity at HRM.