Available online at www.sciencedirect.com ScienceDirect
Advances in Space Research 63 (2019) 2515–2534 www.elsevier.com/locate/asr
Assessing possible mutual orbit period change by shape deformation of Didymos after a kinetic impact in the NASA-led Double Asteroid Redirection Test
Masatoshi Hirabayashi a,⇑, Alex B. Davis b, Eugene G. Fahnestock c, Derek C. Richardson d, Patrick Michel e, Andrew F. Cheng f, Andrew S. Rivkin f, Daniel J. Scheeres b, Steven R. Chesley c, Yang Yu g, Shantanu P. Naidu c, Stephen R. Schwartz h, Lance A.M. Benner c, Petr Pravec i, Angela M. Stickle f, Martin Jutzi j, the DART Dynamical and Physical Properties (WG3) analysis group
a Auburn University, 319 Davis Hall, Auburn, AL 36849, United States b University of Colorado Boulder, 429 UCB, Boulder, CO 80309, United States c Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109, United States d University of Maryland, College Park, MD 20742, United States e Laboratoire Lagrange, Universite´ Coˆte d’Azur, Observatoire de la Coˆte d’Azur, CNRS, CS 34229, 06304 Nice Cedex 4, France f Applied Physics Laboratory/The Johns Hopkins University, Laurel, MD 20723, United States g Beihang University, Beijing 100191, China h University of Arizona, Tucson, AZ 85287, United States i Astronomical Institute, Academy of Sciences of the Czech Republic, Fricˇova 1, CZ-25165 Ondrˇejov, Czech Republic j Physics Institute, Space Research and Planetary Sciences, University of Bern, Sidlerstrasse 5, 3012 Bern, Switzerland
Received 23 June 2018; received in revised form 1 December 2018; accepted 30 December 2018 Available online 14 January 2019
Abstract
The Asteroid Impact & Deflection Assessment (AIDA) targets binary near-Earth asteroid (65803) Didymos. As part of this mission, the NASA-led Double Asteroid Redirection Test (DART) will make a kinetic impactor collide with the smaller secondary of Didymos to test kinetic impact asteroid deflection technology, while the ESA-led Hera mission will evaluate the efficiency of the deflection by con- ducting detailed on-site observations. Research has shown that the larger primary of Didymos is spinning close to its critical spin, and the DART-impact-driven ejecta would give kinetic energy to the primary. It has been hypothesized that such an energy input might cause structural deformation of the primary, affecting the mutual orbit period, a critical parameter for assessing the kinetic impact deflection by the DART impactor. A key issue in the previous work was that the secondary was assumed to be spherical, which may not be realistic. Here, we use a second-order inertia-integral mutual dynamics model to analyze the effects of the shapes of the primary and the secondary on the mutual orbit period change of the system. We first compare the second-order model with three mutual dynamics models, including a high-order inertia-integral model that takes into account the detailed shapes of Didymos. The comparison tests show that the second- order model may have an error of 10% for computing the mutual orbit period change, compared to the high-order model. We next use the second-order model to analyze how the original shape and shape deformation change the mutual orbit period. The results show that when the secondary is elongated, the mutual orbit period becomes short. Also, shape deformation of the secondary further changes the
⇑ Corresponding author. E-mail address: [email protected] (M. Hirabayashi). https://doi.org/10.1016/j.asr.2018.12.041 0273-1177/Ó 2019 COSPAR. Published by Elsevier Ltd. All rights reserved. 2516 M. Hirabayashi et al. / Advances in Space Research 63 (2019) 2515–2534 mutual orbit period. A better understanding of this mechanism allows for detailed assessment of DART’s kinetic impact deflection capa- bility for Didymos. Ó 2019 COSPAR. Published by Elsevier Ltd. All rights reserved.
Keywords: Planetary defense; Asteroid impact hazards; Kinetic impactor; Mutual orbit dynamics; Binary asteroids
1. Introduction et al., 2018). The Hera mission further assesses the b parameter by determining the physical properties including Binary near-Earth asteroid Didymos is the target of the the mass of Didymos (Michel et al., 2018). Asteroid Impact & Deflection Assessment (AIDA) mission, Determination of the b parameter requires detailed eval- which is an international collaboration between NASA and uation of the mutual dynamics and structural condition of ESA (Cheng et al., 2016, 2018; Michel et al., 2016, 2018). In Didymos. A study that conducted impact experiments with this mission, NASA leads the Double Asteroid Redirection aluminum spheres hitting pumice boulders showed that the Test (DART) mission in which a kinetic impactor will be b parameter ranged between 1 and 1.69 (Walker et al., launched to Didymos. As of May 14, 2018, the launch per- 2017). The b parameter is dependent on the surface topog- iod will be opened after June 15, 2021, and the spacecraft raphy (Feldhacker et al., 2017) and surface material condi- will conduct a flyby of asteroid 2001 CB21 and collide with tions such as porosity (Wu¨nnemann et al., 2006; Collins, the smaller secondary of Didymos on October 5, 2022, 2014; Wiggins et al., 2018) and regolith size frequency dis- which is the apparition of Didymos. The DART impactor tributions (Tatsumi and Sugita, 2018; Stickle et al., 2018). will have the Didymos Reconnaissance and Asteroid Cam- Such impact conditions likely change the ejecta conditions era for OpNav (DRACO) visible imager, which will take such as the speed and direction. Therefore, while a fraction high-resolution images of the surface of the secondary right of the DART-impact-driven ejecta may fall onto the sur- before the impact (Cheng et al., 2018). On the other hand, face of the primary (Yu et al., 2017; Yu and Michel, ESA leads the Hera mission, which will send a spacecraft to 2018), its amount is strongly dependent on these impact Didymos to observe the physical properties of this asteroid. conditions. Also, the primary is rotating at a spin period The currently planned instruments onboard the Hera of 2.26 h (Pravec et al., 2006). Given a top-like shape (an spacecraft include the Asteroid Framing Camera (AFC), oblate shape with a raised equatorial ridge), the primary a LIght Detection And Ranging (LIDAR) instrument, might be sensitive to failure of the internal structure, and a six-unit CubeSat carrying two instruments: one being depending on the bulk density and material heterogeneity the Asteroid SPECTral (ASPECT) imaging instrument, (Zhang et al., 2017). If the primary has the reported nom- and the other being chosen from radio science, seismology, inal size and bulk density (Naidu et al., 2016; Michel et al., gravimetry, or volatile detection (Michel et al., 2018). The 2016), it should have the mechanical strength to hold the planned launch date of the Hera spacecraft is October current shape (Hirabayashi et al., 2017). Detailed studies 2023. using Soft-Sphere Discrete Element Methods (SSDEMs) One of the primary goals of the AIDA mission is to have been conducted in multiple research groups, which quantify the asteroid deflection capability using a kinetic consistently confirmed the result by Hirabayashi and impactor. A critical process is the momentum transfer, Scheeres (2014), who pointed out that given a uniform which is augmented by impact-driven ejecta (Holsapple structure, a spheroidal shape might have two distinguish- and Housen, 2012). The efficiency of the momentum trans- able deformation modes: horizontally outward deforma- fer is described using the so-called b parameter. When tion on the equatorial plane and vertically inward b ¼ 1, the impactor collides with the asteroid without net deformation at other locations (e.g. Sa´nchez and generation of momentum-enhancing ejecta so that its Scheeres, 2012, 2016; Zhang et al., 2017, 2018). momentum is transferred only to the target. If b > 1, the Considering these uncertainties, Hirabayashi et al. target asteroid achieves further acceleration compared to (2017) hypothesized that if the DART impact induced the b ¼ 1 case because the momentum carried by the ejecta shape deformation of the primary, the mutual orbit period launched in the opposite direction to the DART impact change, which is a key parameter for determination of the b direction is added to that of the secondary. On the other parameter, would change significantly (the schematics hand, if b < 1 (requiring antipodal spallation, which is given in Fig. 1). Using the results from earlier structural not anticipated for the DART impact), it is accelerated less works (e.g. Hirabayashi and Scheeres, 2014) to assume that than the b ¼ 1 case. In the AIDA mission, the DART mis- the shape deformation process is axisymmetric along the sion demonstrates kinetic impact deflection, determines the spin axis of the primary, they investigated how the mutual orbit period change of Didymos using Earth-based primary’s deformed shape would change the mutual orbit observations and images from the DRACO imager, and period of the system. They found that the primary’s shape attempts to obtain an estimate of the b value (Cheng deformation process might play a significant role in M. Hirabayashi et al. / Advances in Space Research 63 (2019) 2515–2534 2517
Fig. 1. Schematics of a scenario of the mutual orbit period change caused by shape deformation of Didymos after the DART impact. The shape of the primary used is the radar shape model from Michel et al. (2016) and Naidu et al. (2016). In Step 2, the solid line shows the normal case that experiences neither the DART impact nor shape deformation, the dot-dashed line describes a trajectory after the DART impact without shape deformation, and the dashed line gives a perturbed orbit after the DART impact with shape deformation. changing the mutual orbit period. However, an issue with We use the second-order model to analyze how the orig- their work was that the secondary was assumed to be a per- inal shape and shape deformation of Didymos affect the fect sphere although lightcurve observations did not rule mutual orbit period of the system. The present manuscript out that the secondary would be elongated (Pravec et al., contains the following three components: 2006). In fact, many binary asteroids have elongated secon- daries (Pravec et al., 2016); therefore, the assumption made First, we generalize the effects of the primary’s shape by Hirabayashi et al. (2017) might not have accounted for deformation on the mutual orbit period change. the mutual dynamics of Didymos accurately. We hypothe- Second, we investigate how the elongation of the sec- size that the secondary’s elongation also contributes to the ondary changes the mutual orbit period of the system. mutual orbit period change (Fig. 1) and quantify this prob- Third, we explore how the shape deformation process of lem in the present work. the secondary influences the mutual orbit period change The purpose of this work is to use a model that takes of the system. into account the elongation of the secondary to investigate how the mutual orbit period change would occur in the These exercises extend our understanding of the mutual Didymos system after the DART impact. In our work, orbit period change of Didymos due to the DART impact we employ a second-order inertia-integral mutual dynamics from the previous work by Hirabayashi et al. (2017), pro- model in which the shapes of the primary and the sec- viding a better capability of quantifying the kinetic impact ondary are assumed to be an oblate spheroid (i.e., the long asteroid deflection by DART. axis equal to the intermediate axis) and a prolate ellipsoid We organize our discussion as follows. In Section 2,we (i.e., the intermediate axis equal to the short axis), respec- briefly introduce the model used in this work. In Section 3, tively (Scheeres, 2009; McMahon and Scheeres, 2013). This we show a deformation process that we consider in the pre- model is tested by using three approaches capable of com- sent work. Section 4 presents our numerical investigations. puting the mutual orbit period change, which are given In Section 4.1, we show comparison tests of the second- below: order model with other mutual dynamics models. In Sec- tion 4.2, we discuss how the original shape and shape The first is a sphere-sphere dynamics approach that deformation of the primary and secondary affect the assumes the primary and the secondary to be homoge- mutual orbit period change of the system. Finally, Section 5 neous spheres (Cheng et al., 2016). gives interpretations of our numerical results into the The second is a mutual dynamics approach that uses a mutual orbit period change of Didymos. We note that in polyhedron-shape-based gravity model (Werner and the literature, the smaller and larger components of Didy- Scheeres, 1997) for the primary and assumes the sec- mos have been named in different ways (e.g. Michel et al., ondary to be a homogeneous sphere (Hirabayashi 2018; Cheng et al., 2018; Hirabayashi et al., 2017). The lar- et al., 2017). ger component has been called Didymos, Didymmain, The third is an inertia-integral-based mutual dynamics Didymos A, or the primary; on the other hand, the smaller approach that takes into account the detailed shapes component has been called Didymoon, Didymos B, or the using high-order gravity terms (Hou et al., 2017; Davis secondary. In this work, we simply call them the primary and Scheeres, 2017). and the secondary. 2518 M. Hirabayashi et al. / Advances in Space Research 63 (2019) 2515–2534
2. Mutual dynamics modeling 2.2. Polyhedron-shape-based model
This section introduces the models that are used in this This approach uses a polyhedron-shape gravity model work. We use mathematical notations defined in Table 1. (Werner and Scheeres, 1997) to describe the gravity force Recall the approaches used are a sphere-sphere model, a acting on the secondary under the assumption that the sec- polyhedron model, a second-order inertia-integral model, ondary is perfectly spherical (Hirabayashi et al., 2017). The and a high-order inertia-integral model. For the second- secondary is considered to have a non-negligible mass; the order model, the interaction is computed between two rotational motion of the primary evolves due to the orbital biaxial ellipsoids by considering the gravity terms up to sec- motion of the components. Therefore, it is necessary to ond order. This is the main approach used for quantifying take into account the mutual interaction between the pri- the mutual orbit period change due to shape deformation mary and secondary. of the system. The equation of relative translational motion given in the frame rotating with the primary is given as (Scheeres 2.1. Sphere-sphere model et al., 2006), € _ _ Rc þ 2Xp Rc þ Xp Rc þ Xp ðXp RcÞ Cheng et al. (2016) introduced the sphere-sphere model. This model considers that both the primary and the sec- ms @U ¼ 1 þ ; ð2Þ ondary are perfectly spherical. In this case, because these mp @Rc objects can be assumed to be point masses, we can use where Rc is the center of mass (COM) of the secondary rel- the two-body problem theory. Thus, the rotational motion ative to that of the primary, U is the potential energy, and is independent of the orbital motion. Because this theory Xp is the spin vector of the primary. As discussed above, has an analytical solution, the orbit period change in the the rotational state of the primary is coupled with the orbi- system, dT orbit, can be described as tal motion. Therefore, we describe it as sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi _ A3 A3 IpXp þ Xp IpXp ¼ sp; ð3Þ dT p a p b : orbit ¼ 2 G m m 2 G m m ð1Þ ð p þ sÞ ð p þ sÞ where I p is the inertia tensor of the primary, and sp is the torque acting on the primary. In this form, Ab and Aa are the orbital semi-major axes m before and after the DART impact, respectively. p and 2.3. Second-order inertia-integral model ms are the masses of the primary and the secondary, respec- tively. G is the gravitational constant. This model assumes that the primary is an oblate spher- oid, and the secondary is a prolate ellipsoid, instead of con-
Table 1 Notational definitions of key parameters used in this work. Subscript K is replaced with either p or s, where p stands for the primary and s indicates the secondary. The COM stands for the center of mass. Parameters Symbols Units Gravitational constant (¼ 6:674 10 11) G m3 kg 2 s 2 Orbital semi-major axis A m Geometric semi-major axis of a body aK m Mass mK kg Position of small element in a body relative to its COM rK m Force f K N Torque sK N m 2 Inertia tensor IK kg m 1 Spin vector XK rad s Rotational angle hK rad 1 Speed of synchronous orbit vsyn m sec 2 Moment of inertia along the out-of-plane direction IKz kg m Ratio of primary semi-minor axis to semi-major axis v [–] Ratio of secondary semi-minor axis to semi-major axis n [–] Position of secondary COM relative to primary COM Rc m Position of small element in secondary relative to one in primary R m Position of small element in secondary relative to primary COM R m Position of secondary COM relative to primary COM in inertial frame ½Rcx; Rcy ; 0 m Total energy E J Kinetic energy K J Potential energy U J Out-of-plane component of system angular momentum H kg m2 s 1 M. Hirabayashi et al. / Advances in Space Research 63 (2019) 2515–2534 2519 Z Z sidering particular shapes (Fig. 2). By an oblate spheroid, R f s G dmpdms: ¼ 3 ð5Þ we mean that the object’s semi-major axis is the same as p s R its intermediate axis. On the other hand, for the prolate R ellipsoid, the semi-minor axis is the same as the semi- where is the position vector of a body element in the sec- intermediate axis. While this assumption may not be realis- ondary with respect to one in the primary, and R is its mag- tic, the proposed model may be useful to quantify the nitude. In our model, we consider the inertia integrals up to n mutual orbit period change given the fact that there may second order, i.e., ¼ 2. Earlier works considered the remain observational uncertainties of the physical proper- moment of inertia, instead of the inertia integrals, to ties of Didymos after the DART impact. expand Eq. (5) although these models are fundamentally To describe the mutual gravity force, we use an inertia equal because of the order of expansion (Scheeres, 2009; integral, which is given as (Hou et al., 2017; Davis and McMahon and Scheeres, 2013). Scheeres, 2017) We focus on the planar motion as the motion in the out- Z of-plane direction is expected to be small even after the DART impact (Cheng et al., 2016; Hirabayashi et al., T i;j;k xi y j zk dm ; K ¼ K K K K ð4Þ 2017). This consideration is also reasonable for the high- mK order model. In our problem, the high-order model shows that the ratio of the z component to the norm of the x and where subscript K is replaced with either p (primary) or s T y components is less than 0.007 % even if we account for (secondary). Also, ðxK; yK; zKÞ is the position vector of a the z component of the DART impact velocity. Using Rc body element from the COM in its body-fixed frame (see and f s, we describe the equation of relative translational the details in Fig. 3). For the nth-order inertia integral, motion in the inertial frame as i þ j þ k ¼ n should be satisfied. We expand the gravity force vector acting on the secondary, f s, which is given as € 1 1 Rc ¼ þ f s: ð6Þ mp ms
Fig. 2. Illustration of the formulation. In the second-order model, the primary has an oblate shape while the secondary is a prolate ellipsoid.
Fig. 3. Illustration of the axis orientations of the primary and the secondary. For the primary, the z1 axis is defined in the out-of-plane direction, which x y z corresponds to the maximum principal axis of this object, and the 1 and 1 axes are on the orbital plane. Similarly, for the secondary, the 2 axis is equal to z x y its maximum principal axis and is parallel to the 1 axis. The 2 axis is defined along the minimum principal axis, and the 2 axis is along the intermediate principal axis. 2520 M. Hirabayashi et al. / Advances in Space Research 63 (2019) 2515–2534 2 3 On the other hand, the equation of rotational motion in the 4 þ v2 00 6 7 2 frame rotating with the secondary is described as Ap ¼ 4 04þ v 0 5: ð15Þ _ 2 IsXs þ Xs IsXs ¼ ss; ð7Þ 002þ 3v where Is is the inertia tensor, Xs is the spin vector, and ss is The primary’s rotation does not change this matrix. the torque vector acting on the secondary. In the present The next step is to integrate Eq. (14) over the entire mass model, because the primary is an oblate shape, its rota- of the secondary. It is necessary to consider the attitude of tional motion can be decoupled with our consideration the secondary. Because we prefer to describe the sec- (see below). ondary’s inertia integrals in the frame fixed on the sec- We further discuss our formulation by using the ondary, we use a transformation process from the body- schematics in Fig. 2. O and O0 are the COMs of the primary fixed frame to the inertial frame. We obtain the second- and the secondary, respectively. We also define R to order form of the gravity force acting on the secondary as R describe the position vector from O to a small element in f df ; s ¼ s s rK h the secondary. Vector describes the position vector of a2 v2 Gmpms 3 pð1 Þ URc Rc a body element from the COM. We expand Eq. (5) by R3 þ 10 R2 ð16Þ noc c i using the following expressions, a2RT R a2 d 3 s c As cI 3 s ð sIþ2AsÞ R : þ 4 2 c 2Rc 10Rc R ¼R rp; ð8Þ R ¼Rc þ rs; ð9Þ where as is the geometric semi-major axis of the secondary, n to obtain is the secondary’s elongation, i.e., the ratio of the inter- mediate axis to the long axis, I is the identity matrix, ds 2 2 2 R ¼R þ r 2R rp; ð10Þ 2 p is given as ds ¼ 1 þ 2n , and R 2 R2 r2 R r : 2 3 ¼ c þ s þ 2 c s ð11Þ 100 6 7 We then employ the Taylor expansions of these forms with U ¼ 4 0105: ð17Þ x r2 R r =R 2 respect to ¼ð p 2 pÞ for Eq. (10) and 003 2 2 xc ¼ðr þ 2Rc rsÞ=R for Eq. (11). s c This matrix comes from the fact that the primary is an Applying these formulations, we compute f s and ss.We oblate body that is symmetric along the in-plane directions. first compute the gravity force from the entire mass of the Thus, it describes a force contrast between the in-plane primary acting on a small element in the secondary. This components and the out-of-plane component. Also, As is force, denoted as df s, is given as Z given as R 2 3 2 h n2 2 h n2 h h df s ¼ Gdms dmp: ð12Þ cos s þ sin s ð1 Þ cos s sin s 0 R3 6 7 mp 2 2 2 2 As ¼ 4 ð1 n Þ cos hs sin hs n cos hs þ sin hs 0 5; 2 Using Eq. (10), we expand Eq. (12) with respect to x . Con- 00n sidering the terms up to second order, we obtain ð18Þ Z Gdm hs s 3 15 2 where is the rotational angle of the secondary in the iner- df s x x R rp dmp: ¼ 3 1 þ ð Þ ð13Þ R m 2 8 tial frame. Because the motion of this system is assumed to p T be planar, we write Rc as ½Rcx; Rcy; 0 , and thus the z com- We integrate Eq. (13) over the entire mass of the primary. ponent is always zero in Eq. (16). Because the primary is an oblate spheroid, its rotation Similarly, the second-order expression of the torque act- around the short axis does not affect the gravity acting ing on the secondary is computed using Eq. (13), which is on any elements in the secondary. In other words, the grav- given as Z ity is independent of the spin of the primary. We describe 2 3Gmpmsa the motion of the system in the inertial frame. In this s r df s R R : s ¼ s s 5 c As c ð19Þ T s 5Rc frame, R is given as ðRx ; Ry ; Rz Þ , and we obtain "#The equation of rotational motion is then given as Gm dm 3a2 R 3a2 R p s p Ap p 2 2 2 2 2 df s ¼ R þ fðR þ R Þþv R g ; 3Gmpmsa R 3 10 R 2 2 x y z R 4 X_ X X s R R : Is s þ s Is s ¼ 5 c As c ð20Þ 5Rc ð14Þ Because the primary and secondary are axisymmetric, the where ap is the geometric semi-major axis of the primary, v rotational motion along the x and y axes is always zero. is the primary’s oblateness, i.e., the ratio of the short axis to Therefore, only the motion in the z axis is a critical compo- the long axis, and Ap is a matrix defined as nent in this problem and is given as M. Hirabayashi et al. / Advances in Space Research 63 (2019) 2515–2534 2521
2 2 no 3Gmpmsa 1 n 2.4. High-order inertia-integral model I x_ s ð Þ R R h R2 R2 h ; sz s ¼ 5 2 cx cy cos2 s þð cx cyÞsin2 s 10Rc ð21Þ While computational burden increases, including _ higher-order terms of the inertia integrals in Eq. (4) hs xs; 22 ¼ ð Þ increases the accuracy of the gravity calculation because 2 2 where I sz ¼ msas ð1 þ 2n Þ=5, and xs is the spin rate of the this process accounts for the smaller perturbation due to secondary. the gravity variation. The model was developed by Hou If the spin state is synchronized with the orbital motion, et al. (2017) and applied to the mutual dynamics of Didy- xs is constant and identical to the rate of the orbital mos by Davis and Scheeres (2017). Unlike the second-order motion. Considering that the secondary’s long axis points model discussed in Section 2.3, this model takes into towards the COM of the primary, we obtain the orbit account the irregular shapes of Didymos (Davis and speed of the system. Assuming that the mutual gravity Scheeres, 2017). In this section, we briefly summarize this force computed from Eq. (16) is balanced with the centrifu- high-order gravity calculation technique. We note other gal force, we obtain the relative speed when the secondary’s models that describe mutual dynamics of irregularly rotation is synchronized with its orbit, vsyn,as shaped objects (e.g. von Braun, 1991; Maciejewski, 1995; sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiWerner and Scheeres, 2005; Fahnestock and Scheeres, a2 v2 2 2 Gðmp þ msÞ 3 pð1 Þ 3as ð1 n Þ 2006; Richardson et al., 2009; Hirabayashi and Scheeres, vsyn ¼ 1 þ þ : R 2 2 2013; Naidu et al., 2016). c 10Rc 5Rc Given the inertia integrals defined in Eq. (4), we obtain ð23Þ the mutual potential of the system. If we consider an infi- The energy of this system is conserved and described as nite number of the order terms, we can describe the mutual potential as (Hou et al., 2017) E ¼ K þ U; ð24Þ X1 1 where E is the total energy, K is the kinetic energy, and U is U ¼ G U n; ð28Þ Rnþ1 the potential energy of the system. These quantities are as n¼0 c follows: where U n is the order term of the mutual potential, which is m m 1 p s _ 2 _ 2 1 2 2 2 given as K ¼ ðRcx þ Rcy Þþ msas ð1 þ n Þxs ; ð25Þ 2 mp þ ms 10 Xn X n k n k U n ¼ tk a i ;i ;i i ;i ;i b j ;j ;j j ;j ;j ð 1 2 3Þð 4 5 6Þ ð 1 2 3Þð 4 5 6Þ Gmpms 1 Gmpms 2 2 3 Gmpms 2 k n n i ;i ;i i ;i ;i j ;j ;j j ;j ;j U a v a ð2Þ¼ 2½2 ð 1 2 3Þð 41 5 6Þð 1 2 3Þð 4 5 6Þ ¼ Rc þ 10 R3 pð2 þ Þ 10 R3 p c c i i i i i i i j ; i j ; i j i j ; i j ; i j e 1þ 4 e 2þ 5 e 3þ 6 T ð 1þ 1Þ ð 2þ 2Þ ð 3þ 3ÞT 0ð 4þ 4Þ ð 5þ 5Þ ð 6þ 6Þ; 1 Gmpms 2 2 3 Gmpms 2 2 2 2 1 2 3 p s a n a R hs n hs þ 10 R3 s ð1 þ 2 Þ 10 R3 s f cxðcos þ sin Þ c c ð29Þ 2 2 2 2 þ2RcxRcy ð1 n Þ cos hs sin hs þ Rcxðn cos hs þ sin hsÞ: ð26Þ where n is the order, kð2Þ means that k steps up with a size n= n= For Eq. (25), the first term and the second term on the of 2, and ½ 2 defines the integer part of 2(Hou et al., tn; ak bk e ; e right-hand side are the translational effect of the relative 2017). k , and are recursive parameters. 1 2, and motion and the rotational effect of the secondary, respec- e3 are the components of the unit vector of Rc in the rotat- 0 tively. We do not consider the rotational effect of the pri- ing frame. Also, T s means that the inertia integral is given mary because it is independent of the mutual motion in Pin the frame rotatingP with the primary. In addition, 6 6 this problem. Also, Eq. (26) is obtained by integrating f s i k j n k q¼1 q ¼ and q¼1 q ¼ . in Eq. (16) with respect to Rc. Also, the angular momentum Computation of the inertia integral parameters is in the out-of-plane direction is given as accomplished in a parallel manner to spherical harmonics, m m wherein a shape model is decomposed into a set of tetrahe- p s _ _ 1 2 2 H ¼ ðRcxRcy Rcy RcxÞþ msas ð1 þ n Þxs: ð27Þ mp þ ms 5 dra assumed to have uniform density, and the inertia inte- gral of each tetrahedron is summed to compute the inertia Note that this model can only describe the two- integrals for a given shape model (Hou et al., 2017). For v n dimensional motion if either or has a non-unity value. convenience, the resulting inertia integrals are computed The operations regarding inertia integrals from Eqs. (14)– in the rotating frame of each asteroid. Throughout compu- (16) removed the components for the rotational orientation tation of the mutual gravity potential and torques, the iner- in the out-of-plane direction to simplify the present discus- tia integrals of the secondary are rotated into the frame sion. However, because the out-of-plane mode does not fixed on the primary, which requires additional computa- contribute to the mutual orbit period change (Cheng 0 tion of the inertia integral tensor, T s (Hou et al., 2017). et al., 2016; Hirabayashi et al., 2017), we use the derived Also, the mutual gravity torques are computed applying v n forms as the second-order model. Also, when ¼ ¼ 1, a technique that takes partial derivatives of the potential this model becomes identical to the sphere-sphere model. with respect to the components of a direction cosine matrix 2522 M. Hirabayashi et al. / Advances in Space Research 63 (2019) 2515–2534 and then considers cross products of them with the compo- second-order model to analyze the effect of the secondary’s nents of a direction cosine matrix (Maciejewski, 1995; Hou shape deformation on the mutual orbit period change. et al., 2017). To describe the mutual motion between the Thus, we change the elongation parameter, n, i.e., the ratio primary and the secondary, we apply the theory of the full of the secondary’s semi-minor axis to its semi-major axis, two-body problem (e.g. Scheeres, 2002; Fahnestock and to model the secondary’s deformation. If the secondary Scheeres, 2006). becomes less elongated after the DART impact, n becomes larger, i.e., closer to 1. If the secondary becomes more elon- 3. Shape deformation modeling gated, n should be smaller. We also introduce an assumption that after the DART The goal of this research is to analyze how the shape impact occurs, the bodies immediately change their shapes deformation process would change the mutual orbit period but keep their original spin rate.1 However, this assump- after the DART impact as compared to the case without an tion may not be reasonable because it may take some time impact. We add the shape deformation process into the for Didymos to deform completely. Thus, the mutual mutual orbit dynamics models described in Section 2.A dynamics interaction during the deformation process may key issue is that shape deformation highly depends on the cause additional effects on the mutual orbit period change. DART impact condition and the physical conditions of We leave this analysis as our future work. We also note the primary and the secondary. In this work, instead of that the present work does not take into account the tidal proposing the detailed deformation processes, we analyze dissipation effect on the mutual dynamics. Shape irregular- the mutual orbit period change considering a possible ity in a binary system causes the secondary’s libration deformation scenario of the primary and the secondary. (Scheeres et al., 2006). The timescale of tidal dissipation Importantly, in the following deformation scenario, we is on the order of million years for a binary system whose assume that the volume is conserved, considering a con- the primary’s size is about a few hundred kilometers stant bulk density. (Goldreich and Sari, 2009); however, it is highly dependent The primary may experience centrifugal-force-driven on the mass ratio of the system and tidal dissipation and shape deformation that is induced by collisions of may range from a thousand years to a million years DART-impact-driven ejecta with the surface of the pri- (Jacobson and Scheeres, 2011). mary. The primary’s spin period is reported to be 2.26 h (e.g. Michel et al., 2016), implying that the primary may 4. Results be close to structural failure, depending on the density and strength (Zhang et al., 2017; Hirabayashi et al., We use the second-order mutual dynamics model as our 2017). Once the DART impactor hits the surface of the sec- main numerical tool to analyze how the shape and shape ondary, ejecta from the DART impact site orbit the sys- deformation of Didymos affect the mutual orbit period tem. Some may escape from the system, some may fall change. Before discussing the results, we conduct three back, and the rest may fall on the primary (Yu et al., comparison tests: 2017; Yu and Michel, 2018). Under these conditions, the deformation process of the top-like-shaped primary may First, we compare the second-order model with the be driven by surface shedding (e.g., Walsh et al., 2008) sphere-sphere model. and internal deformation (e.g., Hirabayashi and Scheeres, Second, we compare the second-order model with the 2014), a combination of which may make the shape more polyhedron-shape-based mutual dynamics model. oblate (Hirabayashi, 2015). Thus, we consider this mecha- Third, we compare the second-order model with the nism to be the main deformation process of the primary. In high-order model. Note that in these comparison tests, the mutual dynamics models, we change the oblateness of we assume b ¼ 1; in other words, the discussions do the primary under constant volume. The second-order not consider that ejecta delivers additional momentum model has the oblateness parameter, v, to control this con- to the secondary. dition, while the polyhedron-shape-based model changes the aspect ratio of the primary under constant volume. After these comparison tests, we analyze the shape On the other hand, it is reasonable to consider that the effects on the mutual orbit period change, considering the secondary’s spin state is synchronized with its orbit (Pravec following three cases: et al., 2006; Pravec et al., 2016). Because the mutual orbit period of the system is 11.92 h (Michel et al., 2016), the First, we investigate how the mutual orbit period change spin condition may not contribute to the deformation of is affected by the b parameter and shape deformation of the secondary. Instead, the DART impactor may cause the primary. the secondary to deform. When the DART spacecraft hits the surface of the secondary, the shock wave propagates through the interior of the secondary. Here, we consider the secondary’s deformation only by considering its elon- 1 The DART impactor is planned to hit the center of figure of the gation. In the following discussion, we only use the secondary, minimizing the torque acting on it (Cheng et al., 2018). M. Hirabayashi et al. / Advances in Space Research 63 (2019) 2515–2534 2523