Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

High-Resolution Simulations of Catastrophic Disruptions: Resultant Shape Distributions

High-Resolution Simulations of Catastrophic Disruptions: Resultant Shape Distributions The members of asteroid families have various shapes. We investigate the origin of their shapes by high-resolution impact simulations for catastrophic disrup- tions using a Smoothed Particle Hydrodynamics code. Collisional remnants produced through our simulations of the catastrophic disruptions mainly have spherical or bilobed shapes. However, no at remnants with the ratio of minor to major axis lengths c=a . 0:5 are formed. The results of the simulations provide various shapes of asteroids and explain most of the shapes in asteroid families that are supposed to be produced through catastrophic disruptions. However, the present simulations do not explain signi cantly at asteroids. We suggest that these at asteroids may be interlopers or formed through low-velocity col- lisions between member asteroids. Keywords: Impact processes, Asteroidal shapes, Numerical simulations, SPH method, Accretion 1. Introduction Asteroids in the present solar system are considered to be remnant planetesi- mals that existed in the planet formation era (e.g., Petit et al. 2001; Bottke et al. 2005; Kobayashi et al. 2016). Especially, asteroids with diameters larger than 10 km may keep their shapes formed in the primordial environment because the occurrence timescales of destructive collisions (collisional lifetimes) of those asteroids are estimated to be comparable to or much longer than the age of the solar system (O'Brien & Greenberg 2005). On the other hand, the existence of asteroid families implies recent impacts between asteroids (e.g., Hirayama 1918), which alter the shapes of asteroids. Thus the shapes of large asteroids may tell us information about collisions that occurred in the distant or recent past. Owing to material strengths and small sizes of asteroids, they keep a variety of shapes including irregular shapes di erent from spherical shapes (e.g.,Cibulkov a Corresponding author Email address: sugiuraks@elsi.jp (Keisuke Sugiura) Preprint submitted to Planetary and Space Science November 22, 2019 arXiv:1911.08634v2 [astro-ph.EP] 21 Nov 2019 et al. 2016). Shapes of asteroids are estimated or clari ed through light-curve (e.g., Kaasalainen et al. 2003), radar (e.g., Shepard et al. 2017), and in-situ (e.g., Fujiwara et al. 2006) observations. The Database of Asteroid Models from Inversion Techniques (DAMIT: Durech et al. 2010) provides shape mod- els of about 1,600 asteroids obtained from light-curve observations. Thus, the shape distribution is one of important statistical information for asteroids. Shapes of asteroids are mainly formed through asteroidal collisions, and re- sultant shapes produced through impacts depend on impact conditions such as impact velocities (e.g., Leinhardt et al. 2000; Jutzi & Asphaug 2015; Sugiura et al. 2018). The average impact velocity in the present main belt is estimated to be about 5 km=s (e.g., O'Brien & Sykes 2011; Sugiura et al. 2018), so that impacts between asteroids mainly result in catastrophic disruptions. On the other hand, collisional growth of planetesimals may form the size distribution of 100 km sized main-belt asteroids (Kobayashi et al. 2016). The impact veloc- ities in such formation era do not signi cantly exceed their escape velocities, so that impacts mainly result in non-destructive or merging collisions. Therefore, understanding of shapes produced through non-destructive impacts and catas- trophic disruptions may provide valuable information on the formation era of each asteroid. Sugiura et al. (2018) conducted numerical simulations to investigate as- teroidal shapes produced through non-destructive and equal-mass impacts and showed that the impacts produce various shapes including extremely at and elongated shapes. However, they investigate only the largest remnants produced through non-destructive impacts because the numerical resolution is relatively limited, and shapes of remnants produced through catastrophic disruption are not investigated. Schwartz et al. (2018) found that bilobed shapes such as the shape of the comet 67P/Churyumov-Gerasimenko are formed through catas- trophic disruptions. Asteroid families are groups of asteroids that have similar proper orbital elements, and a family is considered to be formed by the destruction of single parent body (Hirayama 1918; Farinella et al. 1996). The ratio of the mass of the largest remnant in a family to the mass of all member asteroids (or the mass of a parent body) tells us how destructive a family forming impact was. Some families such as the Eos family were formed through catastrophic disruptions (Bro z et al. 2013). Thus the comparison between shapes of family asteroids and shapes of collisional outcomes obtained from numerical simulations may provide some insights into the physics of catastrophic disruptions or the evolution of asteroid families. In this study, we investigate the shapes of not only the largest but also smaller remnants produced by catastrophic disruptions through high-resolution simulations. We use a Smoothed Particle Hydrodynamics (SPH) code through- out simulations although previous studies use N -body codes for reaccumulation phases (e.g., Michel & Richardson 2013; Schwartz et al. 2018). We compare shapes of remnants obtained from our simulations with those of asteroids be- longing to families produced through catastrophic disruptions. 2 2. Method, initial conditions, and analysis of results For numerical simulations of catastrophic collisions, we use the SPH method for elastic dynamics (Libersky & Petschek 1991). We utilize the same version of a SPH code described in Sugiura et al. (2018). We use a pressure dependent yielding criterion (Jutzi 2015) and a fracture model (Grady-Kipp fragmentation model: Benz & Asphaug 1995) for undamaged material, and use a friction model (Drucker-Prager yield criterion without cohesion: Jutzi 2015) for completely damaged material. We also solve the self gravity between SPH particles. We use the Tillotson equation of state (Tillotson 1962). The simulation code is parallelized using Framework for Developing Particle Simulator (Iwasawa et al. 2015; Iwasawa et al. 2016). We use the parameter set of basaltic material (Benz & Asphaug 1999) for the Tillotson equation of state and the fracture model. The uncompressed density in this parameter set is 2:7 g=cm . We ignore the increase of the macroporosity of reaccumulated remnants, and thus the density of remnants is almost the same as the uncompressed density. The cohesion of undamaged material is set to 100 MPa. The friction angle for completely damaged rocks is set to 40 . Regarding to other parameters, we use the same values as described in Sugiura et al. (2018). We investigate collisions between uniform, uncompressed, and intact basaltic spheres with zero rotation and radii R and R . We set the radius of the target t i body R = 50 km. We conduct four simulations of catastrophic collisions with the di erent impact conditions: the ratio of the impactor mass to the target mass q = (R =R ) , the impact velocity v , and the impact angle  ( = 0 i t imp imp imp corresponds to a head-on collision). The detailed impact conditions are listed in Table 1. We use about 4 10 SPH particles for each impact simulation. q v (m/s)  ( ) imp imp Impact 1 1/1 350 15 Impact 2 1/4 700 15 Impact 3 1/8 1000 0 Impact 4 1/16 1700 15 Table 1: Impact parameters (colliders' mass ratio q, impact velocity v , and impact angle imp ) for impact simulations. imp We continue the simulations of the impacts and subsequent reaccumulation of fragments until 1:0 10 s after the time of the impacts, which is much longer than the reaccumulation timescale estimated to be 2R =v  1; 600 s, where t esc v is the two-body escape velocity of the two impacting asteroids. Then, we esc nd remnants produced through reaccumulation of fragments using a friend- of-friend algorithm (e.g., Huchra & Geller 1982). We identify swarms of SPH particles with spacing less than 1.5 smoothing lengths as remnants. Finally, we measure axis lengths of remnants composed of more than 5,000 SPH particles using a procedure described in Sugiura et al. (2018). Note that a remnant with 3 5,000 SPH particles is resolved by about 20 SPH particles along each axis so that the shape of such a remnant is resolved to a certain extent. Two axis ratios b=a and c=a characterize shapes of remnants, where a, b, and c represent the major, intermediate, and minor axis lengths, respectively. 3. Results 500 km (c) 1.5x10 s (a) 0.0 s (b) 5.0x10 s 4 4 5 (d) 2.5x10 s (e) 4.0x10 s (f) 1.0x10 s Figure 1: Snapshots of the simulation of Impact 1 (q = 1, v = 350 m=s,  = 15 ). The imp imp snapshots show the face on views of the ejecta curtain produced through the impact. Figures 1 and 2 show the snapshots of the simulation of Impact 1 (q = 1, v = 350 m=s,  = 15 ). The catastrophic disruption produces a wide and imp imp thin ejecta curtain (Figs. 1b and 2b). The gravitational instability and fragmen- tation of the curtain produce lamentary structures (Fig. 1d). The gravitational reaccumulation of fragments in the lamentary structures (Fig. 1e) eventually forms many small remnants (Figs. 1f and 2f ). Note that these processes (the gravitational fragmentation of the sheets, the formation of the laments, and the gravitational reaccumulation of fragments) are independent of the detailed impact conditions or initial particle placement as long as impacts result in catas- trophic disruption. Figure 3 shows the cumulative mass distributions of the remnants produced through the catastrophic collisions. The mass of the largest remnants are 0:1M for Impact 1, 2, and 4, and  0:01M for Impact 3, where M is the mass t t t of the target asteroid. The mass distributions of the remnants show shallow slopes at remnant mass M . 1 10 M , indicating the lack of remnants with r t 4 4 M . 1  10 M . A remnant with 1  10 M consists of about 500 SPH r t t particles, with which remnants are suciently resolved. Therefore, the smallest size of remnants with 1 10 M might be determined by the scale of the self gravity such as the Jeans length. 4 500 km (c) 1.5x10 s (a) 0.0 s (b) 5.0x10 s 4 4 5 (d) 2.5x10 s (e) 4.0x10 s (f) 1.0x10 s Figure 2: The same as Fig. 1, but these snapshots show the edge on views of the ejecta curtain produced through the impact. For the analysis of remnant shapes, the well-resolved remnants that are composed of more than 5,000 SPH particles are required. We obtain 136 well- resolved remnants from Impact 1 - 4. Figure 4 shows the snapshots of major remnants produced through Impact 1. The most remnants have spherical shapes rather than elongated or at shapes. Some remnants have characteristic shapes with two lobes. The bilobed shapes are formed through coalescence of two rounded remnants. Thus, the bilobed shapes of 67P/Churyumov-Gerasimenko visited by the spacecraft Rosetta or (486958) 2014 MU69 (nicknamed Ultima Thule) visited by the spacecraft New Horizons may be formed through catas- trophic disruptions (see also Schwartz et al. 2018). Figure 5 shows the axis ratio bc=a of the 136 well-resolved remnants. For the larger remnants with M  0:01M , the axis ratio bc=a spreads from 0.5 r t to 0.8. The range of the axis ratio is consistent with that of the remnants obtained from the N -body simulations (Schwartz et al. 2018). We observe the trend that the larger remnants have more irregular shapes, and this trend is also observed for the remnants with the major principal axis length > 3 km obtained from Schwartz et al. (2018). In contrast, the smaller remnants in our simulations mainly have bc=a  0:8, but such irregular remnants are rare in the simulations of Schwartz et al. (2018). This di erence may be caused by the di erence of the numerical simulation methods or the numerical resolution. Figure 6 shows the axis ratios of the 136 well-resolved remnants. The rem- nants have the axis ratios of b=a  c=a, i.e., the destructive impacts mainly produce spherical and bilobed remnants. The most elongated bilobed remnant has b=a  0:5 and c=a  0:4. However, there are a few at remnants with c=a  0:6 and there are no at remnants with c=a . 0:5. Thus at shapes are dicult to form through the catastrophic impacts. 5 1000 Impact 1 Impact 2 Impact 3 Impact 4 -6 -5 -4 -3 -2 -1 1.0x10 1.0x10 1.0x10 1.0x10 1.0x10 1.0x10 M /M r t Figure 3: Cumulative mass distributions of the remnants produced through the catastrophic collisions. The horizontal axis shows the mass of the remnants M normalized by the mass of the target asteroid M , the vertical axis shows the cumulative number of the remnants that have the mass larger than M . The red solid curve, the green dotted curve, the blue dashed curve, and the brown chain curve show the cumulative mass distributions of the remnants produced through Impact 1, 2, 3, and 4, respectively. Shapes of resultant remnants depend on the friction angle, especially for complex shapes (see Sugiura et al. 2019). Hence, the fraction of bilobed rem- nants may be a ected by the choice of the friction angle. Note that the friction angle of lunar soil is 30 50 (Heiken et al. 1991) and thus the friction angle of 40 used in our study is the typical value. Figure 7 shows the rotation periods of the 136 well-resolved remnants. The remnants have the wide range of the spin rates longer than 5 h, and there is a peak at the rotation periods of about 10 h. This peak rotation period is consistent with that obtained from the N -body simulation of the collisional destruction (Michel & Richardson 2013). The remnants with the fastest rotation have the rotation periods of about 5 h, which is longer than those of fast rotating asteroids such as the asteroid Bennu (period = 4:3 h: Lauretta et al. 2019) and 1999 KW4 (period = 2:8 h: Ostro et al. 2006). The fast rotation states may need the additional acceleration of the rotation such as the Yarkovsky-O'Keefe- Radzievskii-Paddack e ect. N (M>M ) r r 10 km Figure 4: Remnants produced through Impact 1. This gure shows the shape of the largest remnant (the remnant at the left upper corner) to that of the 35th largest remnant (the remnant at the right bottom corner). The scale nearby the left upper remnant is valid for all the remnants in this gure. 4. Discussion The average impact velocity in the main belt is about 5 km=s, so that similar- mass impacts in the present main belt mainly result in catastrophic disruptions. As we show in Section 3, it is dicult to produce at shapes with c=a . 0:5 through catastrophic disruptions. Impacts with much smaller mass ratios q do not result in catastrophic disruptions even in the present environment. However, such impacts probably do not produce at shapes because deformation is con- sidered to occur on the scale of impactors. The simulations of non-destructive impacts with q = 1=64 and v = 500 3; 000 m=s show that all the largest imp remnants produced through the impacts have c=a > 0:7 (Sugiura 2019). Thus, we expect that at asteroids with c=a . 0:5 are dicult to form in the present environment. Sugiura et al. (2018) show that merging and nearly head-on collisions be- tween equal-mass asteroids form at shapes with c=a . 0:5. Two-body es- cape velocities of 100 km asteroids are about 100 m=s, and impacts with impact velocities comparable to the escape velocities occur in the planet formation era (e.g., Kobayashi et al. 2016). Thus, we suppose that at asteroids with c=a . 0:5 are likely to be formed through low-velocity and non-destructive impacts in the primordial environment. The probability for nearly head-on col- lisions is not so high for impacts with random impact angles, and asteroids that by chance experience such nearly head-on impacts become at shapes. This implies that asteroids that do not have at shapes might be also primordial. Although Bottke et al. (2005) suggest that many asteroids with the diam- 7 Figure 5: Axis ratio of the remnants that are composed of more than 5,000 SPH particles and are produced through the catastrophic collisions. The horizontal axis shows the mass of the remnants M normalized by the mass of the target asteroid M , and the vertical axis shows r t bc=a. The red circles, the green squares, the blue triangles, and the brown diamonds show the axis ratio of the remnants produced through Impact 1, 2, 3, and 4, respectively. eters < 120 km are byproducts of disruptive collisions, O'Brien & Greenberg (2005) estimate that the timescale of collisional destruction for asteroids with the diameters > 10 km is longer than the age of the solar system. The the- oretical estimate for the fraction of primordial asteroids is uncertain. On the other hand, the classi cation of family asteroids using the Yarkovsky V-shapes suggests that some asteroids with the diameters > 10 km are surviving planetes- imals, i.e., they are primordial (Delbo et al. 2019). Thus, we suggest that at asteroids with c=a . 0:5 and the diameters > 10 km are probably primordial. In our future works, we investigate the survivability of at asteroids in more detail because at shapes may be easy to disrupt due to those large cross sections. To verify our expectation, we investigate shapes of asteroids belonging to asteroid families produced through catastrophic disruptions. We consider as- teroid families with V =V < 0:5 as families produced through catastrophic lr tot disruptions, where V is the volume of the largest member and V is the total lr tot volume of all the members. The volume of asteroids is calculated from those diameters. The well-resolved remnants produced through our simulation of the catastrophic disruption have sizes larger than 10 km, so that we investigate as- 8 1 0.8 0.6 0.4 0.2 Impact 1 Impact 2 Impact 3 Impact 4 0 0.2 0.4 0.6 0.8 1 axis ratio c/a Figure 6: Shape distribution of the remnants that are composed of more than 5,000 SPH particles and are produced through the catastrophic collisions. The horizontal axis shows c=a and the vertical axis shows b=a. Spherical shapes are plotted on the right upper side, at shapes are plotted on the left upper side, and elongated shapes are plotted on the left bottom side. The red circles, the green squares, the blue triangles, and the brown diamonds show the axis ratios of the remnants produced through Impact 1, 2, 3, and 4, respectively. teroids with diameters larger than 10 km for the comparison. We obtain asteroid shapes from the shape model given by DAMIT database . We also obtain infor- mation of asteroid families from AstDyS-2 database and diameters of asteroids from JPL small-body database . Figure 8 shows the axis ratios of the asteroids belonging to the families produced through catastrophic disruptions. The axis ratios of the asteroids are mainly b=a  c=a, i.e., the asteroids mainly have spherical and elongated, or bilobed shapes. The most elongated asteroid has b=a  0:5 and c=a  0:3, which is consistent with the axis ratios of the most elongated remnant produced through our simulation (see Fig. 6). However, some asteroids have at shapes with b=a & 0:7 and c=a . 0:5 that are not observed in our simulations. There are some possibilities for the discrepancy regarding to at shapes. A rst possibility is observational diculties. It is dicult to accurately determine http://astro.troja.m .cuni.cz/projects/asteroids3D/web.php https://newton.spacedys.com/astdys2/index.php https://ssd.jpl.nasa.gov/sbdb query.cgi axis ratio b/a 50 Impact 1 Impact 2 Impact 3 Impact 4 -3 -2 -1 0 10 20 1.0x10 1.0x10 1.0x10 M /M Number r t Figure 7: Rotation periods of the remnants that are composed of more than 5,000 SPH particles and are produced through the catastrophic collisions. The horizontal axis shows the mass of the remnants M normalized by the mass of the target asteroid M , and the vertical r t axis shows the rotation periods of the remnants. The red circles, the green squares, the blue triangles, and the brown diamonds show the rotation periods of the remnants produced through Impact 1, 2, 3, and 4, respectively. Right panel shows the histogram of the rotation periods of these remnants. the shapes of asteroids only by light-curve observations because there are obser- vational errors in faint objects and uncertainties in the inversion technique. For example, the shape of (419) Aurelia was rstly estimated to be the at shape with b=a  0:9 and c=a  0:5 through light-curve observations (Hanu s et al. 2016), but later it was revised to be the spherical shape with b=a  0:9 and c=a  0:8 through observations with adaptive optics and occultation (Hanu s et al. 2017). A second possibility is interlopers. Flat asteroids belonging to families pro- duced through catastrophic disruptions may be interlopers. Those at shapes may be formed through low-velocity impacts in the primordial environment, and later they may be mixed up with family asteroids. A third possibility is collisions between member asteroids. Velocity dis- persion of remnants just after disruptions of parent bodies is estimated to be 100 m=s (e.g., Michel et al. 2001), which is comparable to typical two-body escape velocities of large asteroids. Thus relative velocities between disrupted Rotation Period [h] 1 0.8 0.6 0.4 0.2 0 0.2 0.4 0.6 0.8 1 axis ratio c/a Figure 8: Shape distribution of the asteroids that belong to the asteroid families produced through catastrophic disruptions and have diameters larger than 10 km. The horizontal axis shows c=a and the vertical axis shows b=a. remnants are about escape velocities. They may collide with each other with low impact velocities, which may result in the formation of at shapes if collid- ing two bodies have similar masses and nearly head-on collisions occur (Sugiura et al. 2018). A catastrophic collision ejects remnants around a plane (see Fig. 2b). Nearly head-on collisions between the remnants are reasonably prob- able because their orbital planes are similar (Leleu et al. 2018). The existence of at asteroids belonging to asteroid families may suggest the importance of the evolution of shapes through collisions between member asteroids. 5. Summary Collisional lifetimes of asteroids larger than  10 km are longer than the age of the solar system. Although the existence of asteroid families shows that some large asteroids experience catastrophic disruptions, shapes of asteroids larger than  10 km may be formed in the primordial environment and remain the same until today. Shapes of asteroids formed through catastrophic disruptions provide information on the environment of the most recent collisional events since similar-mass impacts in the present environment mainly result in catas- trophic disruptions, while those in the primordial environment mainly result in non-destructive impacts. axis ratio b/a We conduct the high-resolution simulations of the catastrophic disruptions using the SPH method for elastic dynamics with the self gravity, the fracture model, and the friction model. Our simulations show that catastrophic dis- ruptions mainly produce spherical and bilobed shapes with b=a  c=a (see Figs. 4 and 6). Flat shapes with c=a . 0:5 are not formed through our simu- lations. Thus, we suggest that catastrophic disruptions are dicult to produce at asteroids and at asteroids are likely to be formed through low-velocity and non-destructive impacts in the primordial environment. We investigate shapes of asteroids that have diameters larger than 10 km and belong to asteroid families produced through catastrophic disruptions. Many of them have spherical, elongated, or bilobed shapes with b=a  c=a (see Fig. 8). However, some asteroids have at shapes with b=a & 0:7 and c=a . 0:5, which is not consistent with the shapes of the remnants produced in our simulations of the catastrophic disruptions. We discuss some possibilities for the discrepancy regarding to at shapes. Flat asteroids belonging to families produced through catastrophic disruptions may be interlopers, that is, they may be formed in the primordial environment and mixed up with family asteroids. Low-velocity and similar-mass impacts between asteroids belonging to the same family may produce at shapes, and thus the existence of at asteroids in families may suggest the evolution of asteroidal shapes through collisions between member asteroids. Acknowledgements KS is supported by JSPS KAKENHI Grant Number JP 17J01703. HK and SI are supported by Grant-in-Aid for Scienti c Research (18H05436, 18H05437, 18H05438, 17K05632, 17H01105, 17H01103, 16H02160) and by JSPS Core-to- Core Program \International Network of Planetary Science". Numerical simu- lations in this work were carried out on Cray XC50 at Center for Computational Astrophysics, National Astronomical Observatory of Japan. References Benz, W., & Asphaug, E. 1995. Simulations of brittle solids using smooth particle hydrodynamics. Comput. Phys. Comm., 87(May), 253{265. Benz, W., & Asphaug, E. 1999. Catastrophic Disruptions Revisited. Icarus, 142(Nov.), 5{20. Bottke, W. F., Durda, D. D., Nesvorny,  D., Jedicke, R., Morbidelli, A., Vokrouh- licky,  D., & Levison, H. 2005. The fossilized size distribution of the main asteroid belt. Icarus, 175(May), 111{140. Bro z, M., Morbidelli, A., Bottke, W. F., Rozehnal, J., Vokrouhlicky,  D., & Nesvorny,  D. 2013. Constraining the cometary ux through the asteroid belt during the late heavy bombardment. Astron. Astrophys., 551(Mar.), A117. Cibulkov a, H., Durech, J., Vokrouhlicky,  D., Kaasalainen, M., & Oszkiewicz, D. A. 2016. Distribution of spin-axes longitudes and shape elongations of main-belt asteroids. Astron. Astrophys., 596(Nov.), A57. Delbo, Marco, Avdellidou, Chrysa, & Morbidelli, Alessandro. 2019. Ancient and primordial collisional families as the main sources of X-type asteroids of the inner main belt. Astron. Astrophys., 624(Apr), A69. Durech, J., Sidorin, V., & Kaasalainen, M. 2010. DAMIT: a database of asteroid models. Astron. Astrophys., 513(Apr.), A46. Farinella, P., Davis, D. R., & Marzari, F. 1996. Asteroid Families, Old and Young. Pages 45{55 of: Rettig, T., & Hahn, J. M. (eds), Completing the Inventory of the Solar System. Astronomical Society of the Paci c Conference Series, vol. 107. Fujiwara, A., Kawaguchi, J., Yeomans, D. K., Abe, M., Mukai, T., Okada, T., Saito, J., Yano, H., Yoshikawa, M., Scheeres, D. J., Barnouin-Jha, O., Cheng, A. F., Demura, H., Gaskell, R. W., Hirata, N., Ikeda, H., Kominato, T., Miyamoto, H., Nakamura, A. M., Nakamura, R., Sasaki, S., & Uesugi, K. 2006. The Rubble-Pile Asteroid Itokawa as Observed by Hayabusa. Science, 312(June), 1330{1334. Hanu s, J., Durech, J., Oszkiewicz, D. A., Behrend, R., Carry, B., Delbo, M., Adam, O., Afonina, V., Anquetin, R., Antonini, P., Arnold, L., Audejean, M., Aurard, P., Bachschmidt, M., Baduel, B., Barbotin, E., Barroy, P., Bau- douin, P., Berard, L., Berger, N., Bernasconi, L., Bosch, J.-G., Bouley, S., Bozhinova, I., Brins eld, J., Brunetto, L., Canaud, G., Caron, J., Carrier, F., Casalnuovo, G., Casulli, S., Cerda, M., Chalamet, L., Charbonnel, S., Chinaglia, B., Cikota, A., Colas, F., Coliac, J.-F., Collet, A., Coloma, J., Conjat, M., Conseil, E., Costa, R., Crippa, R., Cristofanelli, M., Damerdji, Y., Deback ere, A., Decock, A., D ehais, Q., D el eage, T., Delmelle, S., De- meautis, C., Dr ozd _ z, _ M., Dubos, G., Dulcamara, T., Dumont, M., Durkee, R., Dymock, R., Escalante del Valle, A., Esseiva, N., Esseiva, R., Esteban, M., Fauchez, T., Fauerbach, M., Fauvaud, M., Fauvaud, S., Forn e, E., Four- nel, C., Fradet, D., Garlitz, J., Gerteis, O., Gillier, C., Gillon, M., Giraud, R., Godard, J.-P., Goncalves, R., Hamanowa, H., Hamanowa, H., Hay, K., Hellmich, S., Heterier, S., Higgins, D., Hirsch, R., Hodosan, G., Hren, M., Hygate, A., Innocent, N., Jacquinot, H., Jawahar, S., Jehin, E., Jerosimic, L., Klotz, A., Ko , W., Korlevic, P., Kosturkiewicz, E., Kra t, P., Krugly, Y., Kugel, F., Labrevoir, O., Lecacheux, J., Lehky,  M., Leroy, A., Lesquerbault, B., Lopez-Gonzales, M. J., Lutz, M., Mallecot, B., Manfroid, J., Manzini, F., Marciniak, A., Martin, A., Modave, B., Montaigut, R., Montier, J., Morelle, E., Morton, B., Mottola, S., Naves, R., Nomen, J., Oey, J., Og loza, W., Paiella, M., Pallares, H., Peyrot, A., Pilcher, F., Pirenne, J.-F., Piron, P., Polinsk  a, M., Polotto, M., Poncy, R., Previt, J. P., Reignier, F., Renauld, D., Ricci, D., Richard, F., Rinner, C., Risoldi, V., Robilliard, D., Romeuf, D., Rousseau, G., Roy, R., Ruthro , J., Salom, P. A., Salvador, L., Sanchez, S., 13 Santana-Ros, T., Scholz, A., S en e, G., Ski , B., Sobkowiak, K., Sogorb, P., Sold an, F., Spiridakis, A., Splanska, E., Sposetti, S., Starkey, D., Stephens, R., Stiepen, A., Stoss, R., Strajnic, J., Teng, J.-P., Tumolo, G., Vagnozzi, A., Vanoutryve, B., Vugnon, J. M., Warner, B. D., Waucomont, M., Wertz, O., Winiarski, M., & Wolf, M. 2016. New and updated convex shape models of asteroids based on optical data from a large collaboration network. Astron. Astrophys., 586(Feb.), A108. Hanu s, J., Viikinkoski, M., Marchis, F., Durech, J., Kaasalainen, M., Delbo', M., Herald, D., Frappa, E., Hayamizu, T., Kerr, S., Preston, S., Timerson, B., Dunham, D., & Talbot, J. 2017. Volumes and bulk densities of forty asteroids from ADAM shape modeling. Astron. Astrophys., 601(May), A114. Heiken, G. H., Vaniman, D. T., & French, B. M. 1991. Lunar sourcebook - A user's guide to the moon. Hirayama, K. 1918. Groups of asteroids probably of common origin. Astron. J., 31(Oct.), 185{188. Huchra, J. P., & Geller, M. J. 1982. Groups of galaxies. I - Nearby groups. Astrophys. J., 257(June), 423{437. Iwasawa, M., Tanikawa, A., Hosono, N., Nitadori, K., Muranushi, T., & Makino, J. 2015. FDPS: A Novel Framework for Developing High-performance Par- ticle Simulation Codes for Distributed-memory Systems. Pages 1:1{1:10 of: Proceedings of the 5th International Workshop on Domain-Speci c Languages and High-Level Frameworks for High Performance Computing. WOLFHPC '15. New York, NY, USA: ACM. Iwasawa, M., Tanikawa, A., Hosono, N., Nitadori, K., Muranushi, T., & Makino, J. 2016. Implementation and performance of FDPS: a framework for devel- oping parallel particle simulation codes. Publ. Astron. Soc. Japan, 68(Aug.), Jutzi, M. 2015. SPH calculations of asteroid disruptions: The role of pressure dependent failure models. Planet. Space Sci., 107(Mar.), 3{9. Jutzi, M., & Asphaug, E. 2015. The shape and structure of cometary nuclei as a result of low-velocity accretion. Science, 348(June), 1355{1358. Kaasalainen, M., Kwiatkowski, T., Abe, M., Piironen, J., Nakamura, T., Ohba, Y., Dermawan, B., Farnham, T., Colas, F., Lowry, S., Weissman, P., White- ley, R. J., Tholen, D. J., Larson, S. M., Yoshikawa, M., Toth, I., & Velichko, F. P. 2003. CCD photometry and model of MUSES-C target (25143) 1998 SF36. Astron. Astrophys., 405(July), L29{L32. Kobayashi, H., Tanaka, H., & Okuzumi, S. 2016. From Planetesimals to Planets in Turbulent Protoplanetary Disks. I. Onset of Runaway Growth. Astrophys. J., 817(Feb.), 105. 14 Lauretta, D. S., Dellagiustina, D. N., Bennett, C. A., Golish, D. R., Becker, K. J., Balram-Knutson, S. S., Barnouin, O. S., Becker, T. L., Bottke, W. F., Boynton, W. V., Campins, H., Clark, B. E., Connolly, H. C., Drouet D'Aubigny, C. Y., Dworkin, J. P., Emery, J. P., Enos, H. L., Hamilton, V. E., Hergenrother, C. W., Howell, E. S., Izawa, M. R. M., Kaplan, H. H., Nolan, M. C., Rizk, B., Roper, H. L., Scheeres, D. J., Smith, P. H., Walsh, K. J., Wolner, C. W. V., & Osiris-Rex Team. 2019. The unexpected surface of asteroid (101955) Bennu. Nature, 568(7750), 55{60. Leinhardt, Z. M., Richardson, D. C., & Quinn, T. 2000. Direct N-body Simu- lations of Rubble Pile Collisions. Icarus, 146(July), 133{151. Leleu, A., Jutzi, M., & Rubin, M. 2018. The peculiar shapes of Saturn's small inner moons as evidence of mergers of similar-sized moonlets. Nature Astron., 2(May), 555{561. Libersky, L. D., & Petschek, A. G. 1991. Smooth particle hydrodynamics with strength of materials. Pages 248{257 of: Trease, H. E., Fritts, M. F., & Crow- ley, W. P. (eds), Advances in the Free-Lagrange Method Including Contribu- tions on Adaptive Gridding and the Smooth Particle Hydrodynamics Method. Lecture Notes in Physics, Berlin Springer Verlag, vol. 395. Michel, P., & Richardson, D. C. 2013. Collision and gravitational reaccumula- tion: Possible formation mechanism of the asteroid Itokawa. Astron. Astro- phys., 554(June), L1. Michel, P., Benz, W., Tanga, P., & Richardson, D. C. 2001. Collisions and Grav- itational Reaccumulation: Forming Asteroid Families and Satellites. Science, 294(Nov.), 1696{1700. O'Brien, D. P., & Greenberg, R. 2005. The collisional and dynamical evolution of the main-belt and NEA size distributions. Icarus, 178(Nov.), 179{212. O'Brien, David P., & Sykes, Mark V. 2011. The Origin and Evolution of the Asteroid Belt|Implications for Vesta and Ceres. Space Sci. Rev., 163(1-4), 41{61. Ostro, Steven. J., Margot, Jean-Luc, Benner, Lance A. M., Giorgini, Jon D., Scheeres, Daniel J., Fahnestock, Eugene G., Broschart, Stephen B., Bellerose, Julie, Nolan, Michael C., Magri, Christopher, Pravec, Petr, Scheirich, Petr, Rose, Randy, Jurgens, Raymond F., De Jong, Eric M., & Suzuki, Shigeru. 2006. Radar Imaging of Binary Near-Earth Asteroid (66391) 1999 KW4. Science, 314(5803), 1276{1280. Petit, J.-M., Morbidelli, A., & Chambers, J. 2001. The Primordial Excitation and Clearing of the Asteroid Belt. Icarus, 153(Oct.), 338{347. Schwartz, S. R., Michel, P., Jutzi, M., Marchi, S., Zhang, Y., & Richardson, D. C. 2018. Catastrophic disruptions as the origin of bilobate comets. Nature Astron., 2(May), 379{382. 15 Shepard, M. K., Richardson, J., Taylor, P. A., Rodriguez-Ford, L. A., Conrad, A., de Pater, I., Adamkovics, M., de Kleer, K., Males, J. R., Morzinski, K. M., Close, L. M., Kaasalainen, M., Viikinkoski, M., Timerson, B., Reddy, V., Magri, C., Nolan, M. C., Howell, E. S., Benner, L. A. M., Giorgini, J. D., Warner, B. D., & Harris, A. W. 2017. Radar observations and shape model of asteroid 16 Psyche. Icarus, 281(Jan.), 388{403. Sugiura, K., Kobayashi, H., & Inutsuka, S. 2018. Toward understanding the ori- gin of asteroid geometries. Variety in shapes produced by equal-mass impacts. Astron. Astrophys., 620(Dec.), A167. Sugiura, Keisuke. 2019. Development of the Numerical Simulation Method for Rocky Body Impacts and Theoretical Analysis of Asteroidal Shapes. PhD Thesis, Nagoya University. Sugiura, Keisuke, Kobayashi, Hiroshi, & Inutsuka, Shu-ichiro. 2019. Collisional elongation: Possible origin of extremely elongated shape of 1I/'Oumuamua. Icarus, 328(Aug), 14{22. Tillotson, J. H. 1962. Metallic equations of state for hypervelocity impact. General Atomic Rept., GA-3216, 1{142. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Astrophysics arXiv (Cornell University)

High-Resolution Simulations of Catastrophic Disruptions: Resultant Shape Distributions

Sugiura, Keisuke; Kobayashi, Hiroshi; Inutsuka, Shu-ichiro
Astrophysics , Volume 2019 (1911) – Nov 20, 2019
Download PDF
16 pages

Loading next page...
 
/lp/arxiv-cornell-university/high-resolution-simulations-of-catastrophic-disruptions-resultant-PNNQJqlkdg

References (38)

ISSN
0032-0633
eISSN
ARCH-3330
DOI
10.1016/j.pss.2019.104807
Publisher site
See Article on Publisher Site

Abstract

The members of asteroid families have various shapes. We investigate the origin of their shapes by high-resolution impact simulations for catastrophic disrup- tions using a Smoothed Particle Hydrodynamics code. Collisional remnants produced through our simulations of the catastrophic disruptions mainly have spherical or bilobed shapes. However, no at remnants with the ratio of minor to major axis lengths c=a . 0:5 are formed. The results of the simulations provide various shapes of asteroids and explain most of the shapes in asteroid families that are supposed to be produced through catastrophic disruptions. However, the present simulations do not explain signi cantly at asteroids. We suggest that these at asteroids may be interlopers or formed through low-velocity col- lisions between member asteroids. Keywords: Impact processes, Asteroidal shapes, Numerical simulations, SPH method, Accretion 1. Introduction Asteroids in the present solar system are considered to be remnant planetesi- mals that existed in the planet formation era (e.g., Petit et al. 2001; Bottke et al. 2005; Kobayashi et al. 2016). Especially, asteroids with diameters larger than 10 km may keep their shapes formed in the primordial environment because the occurrence timescales of destructive collisions (collisional lifetimes) of those asteroids are estimated to be comparable to or much longer than the age of the solar system (O'Brien & Greenberg 2005). On the other hand, the existence of asteroid families implies recent impacts between asteroids (e.g., Hirayama 1918), which alter the shapes of asteroids. Thus the shapes of large asteroids may tell us information about collisions that occurred in the distant or recent past. Owing to material strengths and small sizes of asteroids, they keep a variety of shapes including irregular shapes di erent from spherical shapes (e.g.,Cibulkov a Corresponding author Email address: sugiuraks@elsi.jp (Keisuke Sugiura) Preprint submitted to Planetary and Space Science November 22, 2019 arXiv:1911.08634v2 [astro-ph.EP] 21 Nov 2019 et al. 2016). Shapes of asteroids are estimated or clari ed through light-curve (e.g., Kaasalainen et al. 2003), radar (e.g., Shepard et al. 2017), and in-situ (e.g., Fujiwara et al. 2006) observations. The Database of Asteroid Models from Inversion Techniques (DAMIT: Durech et al. 2010) provides shape mod- els of about 1,600 asteroids obtained from light-curve observations. Thus, the shape distribution is one of important statistical information for asteroids. Shapes of asteroids are mainly formed through asteroidal collisions, and re- sultant shapes produced through impacts depend on impact conditions such as impact velocities (e.g., Leinhardt et al. 2000; Jutzi & Asphaug 2015; Sugiura et al. 2018). The average impact velocity in the present main belt is estimated to be about 5 km=s (e.g., O'Brien & Sykes 2011; Sugiura et al. 2018), so that impacts between asteroids mainly result in catastrophic disruptions. On the other hand, collisional growth of planetesimals may form the size distribution of 100 km sized main-belt asteroids (Kobayashi et al. 2016). The impact veloc- ities in such formation era do not signi cantly exceed their escape velocities, so that impacts mainly result in non-destructive or merging collisions. Therefore, understanding of shapes produced through non-destructive impacts and catas- trophic disruptions may provide valuable information on the formation era of each asteroid. Sugiura et al. (2018) conducted numerical simulations to investigate as- teroidal shapes produced through non-destructive and equal-mass impacts and showed that the impacts produce various shapes including extremely at and elongated shapes. However, they investigate only the largest remnants produced through non-destructive impacts because the numerical resolution is relatively limited, and shapes of remnants produced through catastrophic disruption are not investigated. Schwartz et al. (2018) found that bilobed shapes such as the shape of the comet 67P/Churyumov-Gerasimenko are formed through catas- trophic disruptions. Asteroid families are groups of asteroids that have similar proper orbital elements, and a family is considered to be formed by the destruction of single parent body (Hirayama 1918; Farinella et al. 1996). The ratio of the mass of the largest remnant in a family to the mass of all member asteroids (or the mass of a parent body) tells us how destructive a family forming impact was. Some families such as the Eos family were formed through catastrophic disruptions (Bro z et al. 2013). Thus the comparison between shapes of family asteroids and shapes of collisional outcomes obtained from numerical simulations may provide some insights into the physics of catastrophic disruptions or the evolution of asteroid families. In this study, we investigate the shapes of not only the largest but also smaller remnants produced by catastrophic disruptions through high-resolution simulations. We use a Smoothed Particle Hydrodynamics (SPH) code through- out simulations although previous studies use N -body codes for reaccumulation phases (e.g., Michel & Richardson 2013; Schwartz et al. 2018). We compare shapes of remnants obtained from our simulations with those of asteroids be- longing to families produced through catastrophic disruptions. 2 2. Method, initial conditions, and analysis of results For numerical simulations of catastrophic collisions, we use the SPH method for elastic dynamics (Libersky & Petschek 1991). We utilize the same version of a SPH code described in Sugiura et al. (2018). We use a pressure dependent yielding criterion (Jutzi 2015) and a fracture model (Grady-Kipp fragmentation model: Benz & Asphaug 1995) for undamaged material, and use a friction model (Drucker-Prager yield criterion without cohesion: Jutzi 2015) for completely damaged material. We also solve the self gravity between SPH particles. We use the Tillotson equation of state (Tillotson 1962). The simulation code is parallelized using Framework for Developing Particle Simulator (Iwasawa et al. 2015; Iwasawa et al. 2016). We use the parameter set of basaltic material (Benz & Asphaug 1999) for the Tillotson equation of state and the fracture model. The uncompressed density in this parameter set is 2:7 g=cm . We ignore the increase of the macroporosity of reaccumulated remnants, and thus the density of remnants is almost the same as the uncompressed density. The cohesion of undamaged material is set to 100 MPa. The friction angle for completely damaged rocks is set to 40 . Regarding to other parameters, we use the same values as described in Sugiura et al. (2018). We investigate collisions between uniform, uncompressed, and intact basaltic spheres with zero rotation and radii R and R . We set the radius of the target t i body R = 50 km. We conduct four simulations of catastrophic collisions with the di erent impact conditions: the ratio of the impactor mass to the target mass q = (R =R ) , the impact velocity v , and the impact angle  ( = 0 i t imp imp imp corresponds to a head-on collision). The detailed impact conditions are listed in Table 1. We use about 4 10 SPH particles for each impact simulation. q v (m/s)  ( ) imp imp Impact 1 1/1 350 15 Impact 2 1/4 700 15 Impact 3 1/8 1000 0 Impact 4 1/16 1700 15 Table 1: Impact parameters (colliders' mass ratio q, impact velocity v , and impact angle imp ) for impact simulations. imp We continue the simulations of the impacts and subsequent reaccumulation of fragments until 1:0 10 s after the time of the impacts, which is much longer than the reaccumulation timescale estimated to be 2R =v  1; 600 s, where t esc v is the two-body escape velocity of the two impacting asteroids. Then, we esc nd remnants produced through reaccumulation of fragments using a friend- of-friend algorithm (e.g., Huchra & Geller 1982). We identify swarms of SPH particles with spacing less than 1.5 smoothing lengths as remnants. Finally, we measure axis lengths of remnants composed of more than 5,000 SPH particles using a procedure described in Sugiura et al. (2018). Note that a remnant with 3 5,000 SPH particles is resolved by about 20 SPH particles along each axis so that the shape of such a remnant is resolved to a certain extent. Two axis ratios b=a and c=a characterize shapes of remnants, where a, b, and c represent the major, intermediate, and minor axis lengths, respectively. 3. Results 500 km (c) 1.5x10 s (a) 0.0 s (b) 5.0x10 s 4 4 5 (d) 2.5x10 s (e) 4.0x10 s (f) 1.0x10 s Figure 1: Snapshots of the simulation of Impact 1 (q = 1, v = 350 m=s,  = 15 ). The imp imp snapshots show the face on views of the ejecta curtain produced through the impact. Figures 1 and 2 show the snapshots of the simulation of Impact 1 (q = 1, v = 350 m=s,  = 15 ). The catastrophic disruption produces a wide and imp imp thin ejecta curtain (Figs. 1b and 2b). The gravitational instability and fragmen- tation of the curtain produce lamentary structures (Fig. 1d). The gravitational reaccumulation of fragments in the lamentary structures (Fig. 1e) eventually forms many small remnants (Figs. 1f and 2f ). Note that these processes (the gravitational fragmentation of the sheets, the formation of the laments, and the gravitational reaccumulation of fragments) are independent of the detailed impact conditions or initial particle placement as long as impacts result in catas- trophic disruption. Figure 3 shows the cumulative mass distributions of the remnants produced through the catastrophic collisions. The mass of the largest remnants are 0:1M for Impact 1, 2, and 4, and  0:01M for Impact 3, where M is the mass t t t of the target asteroid. The mass distributions of the remnants show shallow slopes at remnant mass M . 1 10 M , indicating the lack of remnants with r t 4 4 M . 1  10 M . A remnant with 1  10 M consists of about 500 SPH r t t particles, with which remnants are suciently resolved. Therefore, the smallest size of remnants with 1 10 M might be determined by the scale of the self gravity such as the Jeans length. 4 500 km (c) 1.5x10 s (a) 0.0 s (b) 5.0x10 s 4 4 5 (d) 2.5x10 s (e) 4.0x10 s (f) 1.0x10 s Figure 2: The same as Fig. 1, but these snapshots show the edge on views of the ejecta curtain produced through the impact. For the analysis of remnant shapes, the well-resolved remnants that are composed of more than 5,000 SPH particles are required. We obtain 136 well- resolved remnants from Impact 1 - 4. Figure 4 shows the snapshots of major remnants produced through Impact 1. The most remnants have spherical shapes rather than elongated or at shapes. Some remnants have characteristic shapes with two lobes. The bilobed shapes are formed through coalescence of two rounded remnants. Thus, the bilobed shapes of 67P/Churyumov-Gerasimenko visited by the spacecraft Rosetta or (486958) 2014 MU69 (nicknamed Ultima Thule) visited by the spacecraft New Horizons may be formed through catas- trophic disruptions (see also Schwartz et al. 2018). Figure 5 shows the axis ratio bc=a of the 136 well-resolved remnants. For the larger remnants with M  0:01M , the axis ratio bc=a spreads from 0.5 r t to 0.8. The range of the axis ratio is consistent with that of the remnants obtained from the N -body simulations (Schwartz et al. 2018). We observe the trend that the larger remnants have more irregular shapes, and this trend is also observed for the remnants with the major principal axis length > 3 km obtained from Schwartz et al. (2018). In contrast, the smaller remnants in our simulations mainly have bc=a  0:8, but such irregular remnants are rare in the simulations of Schwartz et al. (2018). This di erence may be caused by the di erence of the numerical simulation methods or the numerical resolution. Figure 6 shows the axis ratios of the 136 well-resolved remnants. The rem- nants have the axis ratios of b=a  c=a, i.e., the destructive impacts mainly produce spherical and bilobed remnants. The most elongated bilobed remnant has b=a  0:5 and c=a  0:4. However, there are a few at remnants with c=a  0:6 and there are no at remnants with c=a . 0:5. Thus at shapes are dicult to form through the catastrophic impacts. 5 1000 Impact 1 Impact 2 Impact 3 Impact 4 -6 -5 -4 -3 -2 -1 1.0x10 1.0x10 1.0x10 1.0x10 1.0x10 1.0x10 M /M r t Figure 3: Cumulative mass distributions of the remnants produced through the catastrophic collisions. The horizontal axis shows the mass of the remnants M normalized by the mass of the target asteroid M , the vertical axis shows the cumulative number of the remnants that have the mass larger than M . The red solid curve, the green dotted curve, the blue dashed curve, and the brown chain curve show the cumulative mass distributions of the remnants produced through Impact 1, 2, 3, and 4, respectively. Shapes of resultant remnants depend on the friction angle, especially for complex shapes (see Sugiura et al. 2019). Hence, the fraction of bilobed rem- nants may be a ected by the choice of the friction angle. Note that the friction angle of lunar soil is 30 50 (Heiken et al. 1991) and thus the friction angle of 40 used in our study is the typical value. Figure 7 shows the rotation periods of the 136 well-resolved remnants. The remnants have the wide range of the spin rates longer than 5 h, and there is a peak at the rotation periods of about 10 h. This peak rotation period is consistent with that obtained from the N -body simulation of the collisional destruction (Michel & Richardson 2013). The remnants with the fastest rotation have the rotation periods of about 5 h, which is longer than those of fast rotating asteroids such as the asteroid Bennu (period = 4:3 h: Lauretta et al. 2019) and 1999 KW4 (period = 2:8 h: Ostro et al. 2006). The fast rotation states may need the additional acceleration of the rotation such as the Yarkovsky-O'Keefe- Radzievskii-Paddack e ect. N (M>M ) r r 10 km Figure 4: Remnants produced through Impact 1. This gure shows the shape of the largest remnant (the remnant at the left upper corner) to that of the 35th largest remnant (the remnant at the right bottom corner). The scale nearby the left upper remnant is valid for all the remnants in this gure. 4. Discussion The average impact velocity in the main belt is about 5 km=s, so that similar- mass impacts in the present main belt mainly result in catastrophic disruptions. As we show in Section 3, it is dicult to produce at shapes with c=a . 0:5 through catastrophic disruptions. Impacts with much smaller mass ratios q do not result in catastrophic disruptions even in the present environment. However, such impacts probably do not produce at shapes because deformation is con- sidered to occur on the scale of impactors. The simulations of non-destructive impacts with q = 1=64 and v = 500 3; 000 m=s show that all the largest imp remnants produced through the impacts have c=a > 0:7 (Sugiura 2019). Thus, we expect that at asteroids with c=a . 0:5 are dicult to form in the present environment. Sugiura et al. (2018) show that merging and nearly head-on collisions be- tween equal-mass asteroids form at shapes with c=a . 0:5. Two-body es- cape velocities of 100 km asteroids are about 100 m=s, and impacts with impact velocities comparable to the escape velocities occur in the planet formation era (e.g., Kobayashi et al. 2016). Thus, we suppose that at asteroids with c=a . 0:5 are likely to be formed through low-velocity and non-destructive impacts in the primordial environment. The probability for nearly head-on col- lisions is not so high for impacts with random impact angles, and asteroids that by chance experience such nearly head-on impacts become at shapes. This implies that asteroids that do not have at shapes might be also primordial. Although Bottke et al. (2005) suggest that many asteroids with the diam- 7 Figure 5: Axis ratio of the remnants that are composed of more than 5,000 SPH particles and are produced through the catastrophic collisions. The horizontal axis shows the mass of the remnants M normalized by the mass of the target asteroid M , and the vertical axis shows r t bc=a. The red circles, the green squares, the blue triangles, and the brown diamonds show the axis ratio of the remnants produced through Impact 1, 2, 3, and 4, respectively. eters < 120 km are byproducts of disruptive collisions, O'Brien & Greenberg (2005) estimate that the timescale of collisional destruction for asteroids with the diameters > 10 km is longer than the age of the solar system. The the- oretical estimate for the fraction of primordial asteroids is uncertain. On the other hand, the classi cation of family asteroids using the Yarkovsky V-shapes suggests that some asteroids with the diameters > 10 km are surviving planetes- imals, i.e., they are primordial (Delbo et al. 2019). Thus, we suggest that at asteroids with c=a . 0:5 and the diameters > 10 km are probably primordial. In our future works, we investigate the survivability of at asteroids in more detail because at shapes may be easy to disrupt due to those large cross sections. To verify our expectation, we investigate shapes of asteroids belonging to asteroid families produced through catastrophic disruptions. We consider as- teroid families with V =V < 0:5 as families produced through catastrophic lr tot disruptions, where V is the volume of the largest member and V is the total lr tot volume of all the members. The volume of asteroids is calculated from those diameters. The well-resolved remnants produced through our simulation of the catastrophic disruption have sizes larger than 10 km, so that we investigate as- 8 1 0.8 0.6 0.4 0.2 Impact 1 Impact 2 Impact 3 Impact 4 0 0.2 0.4 0.6 0.8 1 axis ratio c/a Figure 6: Shape distribution of the remnants that are composed of more than 5,000 SPH particles and are produced through the catastrophic collisions. The horizontal axis shows c=a and the vertical axis shows b=a. Spherical shapes are plotted on the right upper side, at shapes are plotted on the left upper side, and elongated shapes are plotted on the left bottom side. The red circles, the green squares, the blue triangles, and the brown diamonds show the axis ratios of the remnants produced through Impact 1, 2, 3, and 4, respectively. teroids with diameters larger than 10 km for the comparison. We obtain asteroid shapes from the shape model given by DAMIT database . We also obtain infor- mation of asteroid families from AstDyS-2 database and diameters of asteroids from JPL small-body database . Figure 8 shows the axis ratios of the asteroids belonging to the families produced through catastrophic disruptions. The axis ratios of the asteroids are mainly b=a  c=a, i.e., the asteroids mainly have spherical and elongated, or bilobed shapes. The most elongated asteroid has b=a  0:5 and c=a  0:3, which is consistent with the axis ratios of the most elongated remnant produced through our simulation (see Fig. 6). However, some asteroids have at shapes with b=a & 0:7 and c=a . 0:5 that are not observed in our simulations. There are some possibilities for the discrepancy regarding to at shapes. A rst possibility is observational diculties. It is dicult to accurately determine http://astro.troja.m .cuni.cz/projects/asteroids3D/web.php https://newton.spacedys.com/astdys2/index.php https://ssd.jpl.nasa.gov/sbdb query.cgi axis ratio b/a 50 Impact 1 Impact 2 Impact 3 Impact 4 -3 -2 -1 0 10 20 1.0x10 1.0x10 1.0x10 M /M Number r t Figure 7: Rotation periods of the remnants that are composed of more than 5,000 SPH particles and are produced through the catastrophic collisions. The horizontal axis shows the mass of the remnants M normalized by the mass of the target asteroid M , and the vertical r t axis shows the rotation periods of the remnants. The red circles, the green squares, the blue triangles, and the brown diamonds show the rotation periods of the remnants produced through Impact 1, 2, 3, and 4, respectively. Right panel shows the histogram of the rotation periods of these remnants. the shapes of asteroids only by light-curve observations because there are obser- vational errors in faint objects and uncertainties in the inversion technique. For example, the shape of (419) Aurelia was rstly estimated to be the at shape with b=a  0:9 and c=a  0:5 through light-curve observations (Hanu s et al. 2016), but later it was revised to be the spherical shape with b=a  0:9 and c=a  0:8 through observations with adaptive optics and occultation (Hanu s et al. 2017). A second possibility is interlopers. Flat asteroids belonging to families pro- duced through catastrophic disruptions may be interlopers. Those at shapes may be formed through low-velocity impacts in the primordial environment, and later they may be mixed up with family asteroids. A third possibility is collisions between member asteroids. Velocity dis- persion of remnants just after disruptions of parent bodies is estimated to be 100 m=s (e.g., Michel et al. 2001), which is comparable to typical two-body escape velocities of large asteroids. Thus relative velocities between disrupted Rotation Period [h] 1 0.8 0.6 0.4 0.2 0 0.2 0.4 0.6 0.8 1 axis ratio c/a Figure 8: Shape distribution of the asteroids that belong to the asteroid families produced through catastrophic disruptions and have diameters larger than 10 km. The horizontal axis shows c=a and the vertical axis shows b=a. remnants are about escape velocities. They may collide with each other with low impact velocities, which may result in the formation of at shapes if collid- ing two bodies have similar masses and nearly head-on collisions occur (Sugiura et al. 2018). A catastrophic collision ejects remnants around a plane (see Fig. 2b). Nearly head-on collisions between the remnants are reasonably prob- able because their orbital planes are similar (Leleu et al. 2018). The existence of at asteroids belonging to asteroid families may suggest the importance of the evolution of shapes through collisions between member asteroids. 5. Summary Collisional lifetimes of asteroids larger than  10 km are longer than the age of the solar system. Although the existence of asteroid families shows that some large asteroids experience catastrophic disruptions, shapes of asteroids larger than  10 km may be formed in the primordial environment and remain the same until today. Shapes of asteroids formed through catastrophic disruptions provide information on the environment of the most recent collisional events since similar-mass impacts in the present environment mainly result in catas- trophic disruptions, while those in the primordial environment mainly result in non-destructive impacts. axis ratio b/a We conduct the high-resolution simulations of the catastrophic disruptions using the SPH method for elastic dynamics with the self gravity, the fracture model, and the friction model. Our simulations show that catastrophic dis- ruptions mainly produce spherical and bilobed shapes with b=a  c=a (see Figs. 4 and 6). Flat shapes with c=a . 0:5 are not formed through our simu- lations. Thus, we suggest that catastrophic disruptions are dicult to produce at asteroids and at asteroids are likely to be formed through low-velocity and non-destructive impacts in the primordial environment. We investigate shapes of asteroids that have diameters larger than 10 km and belong to asteroid families produced through catastrophic disruptions. Many of them have spherical, elongated, or bilobed shapes with b=a  c=a (see Fig. 8). However, some asteroids have at shapes with b=a & 0:7 and c=a . 0:5, which is not consistent with the shapes of the remnants produced in our simulations of the catastrophic disruptions. We discuss some possibilities for the discrepancy regarding to at shapes. Flat asteroids belonging to families produced through catastrophic disruptions may be interlopers, that is, they may be formed in the primordial environment and mixed up with family asteroids. Low-velocity and similar-mass impacts between asteroids belonging to the same family may produce at shapes, and thus the existence of at asteroids in families may suggest the evolution of asteroidal shapes through collisions between member asteroids. Acknowledgements KS is supported by JSPS KAKENHI Grant Number JP 17J01703. HK and SI are supported by Grant-in-Aid for Scienti c Research (18H05436, 18H05437, 18H05438, 17K05632, 17H01105, 17H01103, 16H02160) and by JSPS Core-to- Core Program \International Network of Planetary Science". Numerical simu- lations in this work were carried out on Cray XC50 at Center for Computational Astrophysics, National Astronomical Observatory of Japan. References Benz, W., & Asphaug, E. 1995. Simulations of brittle solids using smooth particle hydrodynamics. Comput. Phys. Comm., 87(May), 253{265. Benz, W., & Asphaug, E. 1999. Catastrophic Disruptions Revisited. Icarus, 142(Nov.), 5{20. Bottke, W. F., Durda, D. D., Nesvorny,  D., Jedicke, R., Morbidelli, A., Vokrouh- licky,  D., & Levison, H. 2005. The fossilized size distribution of the main asteroid belt. Icarus, 175(May), 111{140. Bro z, M., Morbidelli, A., Bottke, W. F., Rozehnal, J., Vokrouhlicky,  D., & Nesvorny,  D. 2013. Constraining the cometary ux through the asteroid belt during the late heavy bombardment. Astron. Astrophys., 551(Mar.), A117. Cibulkov a, H., Durech, J., Vokrouhlicky,  D., Kaasalainen, M., & Oszkiewicz, D. A. 2016. Distribution of spin-axes longitudes and shape elongations of main-belt asteroids. Astron. Astrophys., 596(Nov.), A57. Delbo, Marco, Avdellidou, Chrysa, & Morbidelli, Alessandro. 2019. Ancient and primordial collisional families as the main sources of X-type asteroids of the inner main belt. Astron. Astrophys., 624(Apr), A69. Durech, J., Sidorin, V., & Kaasalainen, M. 2010. DAMIT: a database of asteroid models. Astron. Astrophys., 513(Apr.), A46. Farinella, P., Davis, D. R., & Marzari, F. 1996. Asteroid Families, Old and Young. Pages 45{55 of: Rettig, T., & Hahn, J. M. (eds), Completing the Inventory of the Solar System. Astronomical Society of the Paci c Conference Series, vol. 107. Fujiwara, A., Kawaguchi, J., Yeomans, D. K., Abe, M., Mukai, T., Okada, T., Saito, J., Yano, H., Yoshikawa, M., Scheeres, D. J., Barnouin-Jha, O., Cheng, A. F., Demura, H., Gaskell, R. W., Hirata, N., Ikeda, H., Kominato, T., Miyamoto, H., Nakamura, A. M., Nakamura, R., Sasaki, S., & Uesugi, K. 2006. The Rubble-Pile Asteroid Itokawa as Observed by Hayabusa. Science, 312(June), 1330{1334. Hanu s, J., Durech, J., Oszkiewicz, D. A., Behrend, R., Carry, B., Delbo, M., Adam, O., Afonina, V., Anquetin, R., Antonini, P., Arnold, L., Audejean, M., Aurard, P., Bachschmidt, M., Baduel, B., Barbotin, E., Barroy, P., Bau- douin, P., Berard, L., Berger, N., Bernasconi, L., Bosch, J.-G., Bouley, S., Bozhinova, I., Brins eld, J., Brunetto, L., Canaud, G., Caron, J., Carrier, F., Casalnuovo, G., Casulli, S., Cerda, M., Chalamet, L., Charbonnel, S., Chinaglia, B., Cikota, A., Colas, F., Coliac, J.-F., Collet, A., Coloma, J., Conjat, M., Conseil, E., Costa, R., Crippa, R., Cristofanelli, M., Damerdji, Y., Deback ere, A., Decock, A., D ehais, Q., D el eage, T., Delmelle, S., De- meautis, C., Dr ozd _ z, _ M., Dubos, G., Dulcamara, T., Dumont, M., Durkee, R., Dymock, R., Escalante del Valle, A., Esseiva, N., Esseiva, R., Esteban, M., Fauchez, T., Fauerbach, M., Fauvaud, M., Fauvaud, S., Forn e, E., Four- nel, C., Fradet, D., Garlitz, J., Gerteis, O., Gillier, C., Gillon, M., Giraud, R., Godard, J.-P., Goncalves, R., Hamanowa, H., Hamanowa, H., Hay, K., Hellmich, S., Heterier, S., Higgins, D., Hirsch, R., Hodosan, G., Hren, M., Hygate, A., Innocent, N., Jacquinot, H., Jawahar, S., Jehin, E., Jerosimic, L., Klotz, A., Ko , W., Korlevic, P., Kosturkiewicz, E., Kra t, P., Krugly, Y., Kugel, F., Labrevoir, O., Lecacheux, J., Lehky,  M., Leroy, A., Lesquerbault, B., Lopez-Gonzales, M. J., Lutz, M., Mallecot, B., Manfroid, J., Manzini, F., Marciniak, A., Martin, A., Modave, B., Montaigut, R., Montier, J., Morelle, E., Morton, B., Mottola, S., Naves, R., Nomen, J., Oey, J., Og loza, W., Paiella, M., Pallares, H., Peyrot, A., Pilcher, F., Pirenne, J.-F., Piron, P., Polinsk  a, M., Polotto, M., Poncy, R., Previt, J. P., Reignier, F., Renauld, D., Ricci, D., Richard, F., Rinner, C., Risoldi, V., Robilliard, D., Romeuf, D., Rousseau, G., Roy, R., Ruthro , J., Salom, P. A., Salvador, L., Sanchez, S., 13 Santana-Ros, T., Scholz, A., S en e, G., Ski , B., Sobkowiak, K., Sogorb, P., Sold an, F., Spiridakis, A., Splanska, E., Sposetti, S., Starkey, D., Stephens, R., Stiepen, A., Stoss, R., Strajnic, J., Teng, J.-P., Tumolo, G., Vagnozzi, A., Vanoutryve, B., Vugnon, J. M., Warner, B. D., Waucomont, M., Wertz, O., Winiarski, M., & Wolf, M. 2016. New and updated convex shape models of asteroids based on optical data from a large collaboration network. Astron. Astrophys., 586(Feb.), A108. Hanu s, J., Viikinkoski, M., Marchis, F., Durech, J., Kaasalainen, M., Delbo', M., Herald, D., Frappa, E., Hayamizu, T., Kerr, S., Preston, S., Timerson, B., Dunham, D., & Talbot, J. 2017. Volumes and bulk densities of forty asteroids from ADAM shape modeling. Astron. Astrophys., 601(May), A114. Heiken, G. H., Vaniman, D. T., & French, B. M. 1991. Lunar sourcebook - A user's guide to the moon. Hirayama, K. 1918. Groups of asteroids probably of common origin. Astron. J., 31(Oct.), 185{188. Huchra, J. P., & Geller, M. J. 1982. Groups of galaxies. I - Nearby groups. Astrophys. J., 257(June), 423{437. Iwasawa, M., Tanikawa, A., Hosono, N., Nitadori, K., Muranushi, T., & Makino, J. 2015. FDPS: A Novel Framework for Developing High-performance Par- ticle Simulation Codes for Distributed-memory Systems. Pages 1:1{1:10 of: Proceedings of the 5th International Workshop on Domain-Speci c Languages and High-Level Frameworks for High Performance Computing. WOLFHPC '15. New York, NY, USA: ACM. Iwasawa, M., Tanikawa, A., Hosono, N., Nitadori, K., Muranushi, T., & Makino, J. 2016. Implementation and performance of FDPS: a framework for devel- oping parallel particle simulation codes. Publ. Astron. Soc. Japan, 68(Aug.), Jutzi, M. 2015. SPH calculations of asteroid disruptions: The role of pressure dependent failure models. Planet. Space Sci., 107(Mar.), 3{9. Jutzi, M., & Asphaug, E. 2015. The shape and structure of cometary nuclei as a result of low-velocity accretion. Science, 348(June), 1355{1358. Kaasalainen, M., Kwiatkowski, T., Abe, M., Piironen, J., Nakamura, T., Ohba, Y., Dermawan, B., Farnham, T., Colas, F., Lowry, S., Weissman, P., White- ley, R. J., Tholen, D. J., Larson, S. M., Yoshikawa, M., Toth, I., & Velichko, F. P. 2003. CCD photometry and model of MUSES-C target (25143) 1998 SF36. Astron. Astrophys., 405(July), L29{L32. Kobayashi, H., Tanaka, H., & Okuzumi, S. 2016. From Planetesimals to Planets in Turbulent Protoplanetary Disks. I. Onset of Runaway Growth. Astrophys. J., 817(Feb.), 105. 14 Lauretta, D. S., Dellagiustina, D. N., Bennett, C. A., Golish, D. R., Becker, K. J., Balram-Knutson, S. S., Barnouin, O. S., Becker, T. L., Bottke, W. F., Boynton, W. V., Campins, H., Clark, B. E., Connolly, H. C., Drouet D'Aubigny, C. Y., Dworkin, J. P., Emery, J. P., Enos, H. L., Hamilton, V. E., Hergenrother, C. W., Howell, E. S., Izawa, M. R. M., Kaplan, H. H., Nolan, M. C., Rizk, B., Roper, H. L., Scheeres, D. J., Smith, P. H., Walsh, K. J., Wolner, C. W. V., & Osiris-Rex Team. 2019. The unexpected surface of asteroid (101955) Bennu. Nature, 568(7750), 55{60. Leinhardt, Z. M., Richardson, D. C., & Quinn, T. 2000. Direct N-body Simu- lations of Rubble Pile Collisions. Icarus, 146(July), 133{151. Leleu, A., Jutzi, M., & Rubin, M. 2018. The peculiar shapes of Saturn's small inner moons as evidence of mergers of similar-sized moonlets. Nature Astron., 2(May), 555{561. Libersky, L. D., & Petschek, A. G. 1991. Smooth particle hydrodynamics with strength of materials. Pages 248{257 of: Trease, H. E., Fritts, M. F., & Crow- ley, W. P. (eds), Advances in the Free-Lagrange Method Including Contribu- tions on Adaptive Gridding and the Smooth Particle Hydrodynamics Method. Lecture Notes in Physics, Berlin Springer Verlag, vol. 395. Michel, P., & Richardson, D. C. 2013. Collision and gravitational reaccumula- tion: Possible formation mechanism of the asteroid Itokawa. Astron. Astro- phys., 554(June), L1. Michel, P., Benz, W., Tanga, P., & Richardson, D. C. 2001. Collisions and Grav- itational Reaccumulation: Forming Asteroid Families and Satellites. Science, 294(Nov.), 1696{1700. O'Brien, D. P., & Greenberg, R. 2005. The collisional and dynamical evolution of the main-belt and NEA size distributions. Icarus, 178(Nov.), 179{212. O'Brien, David P., & Sykes, Mark V. 2011. The Origin and Evolution of the Asteroid Belt|Implications for Vesta and Ceres. Space Sci. Rev., 163(1-4), 41{61. Ostro, Steven. J., Margot, Jean-Luc, Benner, Lance A. M., Giorgini, Jon D., Scheeres, Daniel J., Fahnestock, Eugene G., Broschart, Stephen B., Bellerose, Julie, Nolan, Michael C., Magri, Christopher, Pravec, Petr, Scheirich, Petr, Rose, Randy, Jurgens, Raymond F., De Jong, Eric M., & Suzuki, Shigeru. 2006. Radar Imaging of Binary Near-Earth Asteroid (66391) 1999 KW4. Science, 314(5803), 1276{1280. Petit, J.-M., Morbidelli, A., & Chambers, J. 2001. The Primordial Excitation and Clearing of the Asteroid Belt. Icarus, 153(Oct.), 338{347. Schwartz, S. R., Michel, P., Jutzi, M., Marchi, S., Zhang, Y., & Richardson, D. C. 2018. Catastrophic disruptions as the origin of bilobate comets. Nature Astron., 2(May), 379{382. 15 Shepard, M. K., Richardson, J., Taylor, P. A., Rodriguez-Ford, L. A., Conrad, A., de Pater, I., Adamkovics, M., de Kleer, K., Males, J. R., Morzinski, K. M., Close, L. M., Kaasalainen, M., Viikinkoski, M., Timerson, B., Reddy, V., Magri, C., Nolan, M. C., Howell, E. S., Benner, L. A. M., Giorgini, J. D., Warner, B. D., & Harris, A. W. 2017. Radar observations and shape model of asteroid 16 Psyche. Icarus, 281(Jan.), 388{403. Sugiura, K., Kobayashi, H., & Inutsuka, S. 2018. Toward understanding the ori- gin of asteroid geometries. Variety in shapes produced by equal-mass impacts. Astron. Astrophys., 620(Dec.), A167. Sugiura, Keisuke. 2019. Development of the Numerical Simulation Method for Rocky Body Impacts and Theoretical Analysis of Asteroidal Shapes. PhD Thesis, Nagoya University. Sugiura, Keisuke, Kobayashi, Hiroshi, & Inutsuka, Shu-ichiro. 2019. Collisional elongation: Possible origin of extremely elongated shape of 1I/'Oumuamua. Icarus, 328(Aug), 14{22. Tillotson, J. H. 1962. Metallic equations of state for hypervelocity impact. General Atomic Rept., GA-3216, 1{142.

Journal

Astrophysics – arXiv (Cornell University)

Published: Nov 20, 2019

There are no references for this article.