Showing 1–50 of 52 results for author: Schwartz, S

BYORP and Dissipation in Binary Asteroids: Lessons from DART

Authors: Matija Ćuk, Harrison Agrusa, Rachel H. Cueva, Fabio Ferrari, Masatoshi Hirabayashi, Seth A. Jacobson, Jay McMahon, Patrick Michel, Paul Sánchez, Daniel J. Scheeres, Stephen Schwartz, Kevin J. Walsh, Yun Zhang

Abstract: The Near-Earth binary asteroid Didymos was the target of a planetary defense demonstration mission DART in September 2022. The smaller binary component, Dimorphos, was impacted by the spacecraft in order to measure momentum transfer in kinetic impacts into rubble piles. DART and associated Earth-based observation campaigns have provided a wealth of scientific data on the Didymos-Dimorphos binary.… ▽ More The Near-Earth binary asteroid Didymos was the target of a planetary defense demonstration mission DART in September 2022. The smaller binary component, Dimorphos, was impacted by the spacecraft in order to measure momentum transfer in kinetic impacts into rubble piles. DART and associated Earth-based observation campaigns have provided a wealth of scientific data on the Didymos-Dimorphos binary. DART revealed a largely oblate and ellipsoidal shape of Dimorphos before the impact, while the post-impact observations suggest that Dimorphos now has a prolate shape. Here we add those data points to the known properties of small binary asteroids and propose new paradigms of the radiative binary YORP (BYORP) effect as well as tidal dissipation in small binaries. We find that relatively spheroidal bodies like Dimorphos made of small debris may experience a weaker and more size-dependent BYORP effect than previously thought. This could explain the observed values of period drift in several well-characterized binaries. We also propose that energy dissipation in small binaries is dominated by relatively brief episodes of large-scale movement of (likely surface) materials, rather than long-term steady-state tidal dissipation. We propose that one such episode was triggered on Dimorphos by the DART impact. Depending on the longevity of this high-dissipation regime, it is possible that Dimorphos will be more dynamically relaxed in time for the Hera mission than it was in the weeks following the impact. △ Less

Submitted 28 June, 2024; originally announced June 2024.

Comments: Accepted for PSJ

Earth's Alfvén wings driven by the April 2023 Coronal Mass Ejection

Authors: Li-Jen Chen, Daniel Gershman, Brandon Burkholder, Yuxi Chen, Menelaos Sarantos, Lan Jian, James Drake, Chuanfei Dong, Harsha Gurram, Jason Shuster, Daniel Graham, Olivier Le Contel, Steven Schwartz, Stephen Fuselier, Hadi Madanian, Craig Pollock, Haoming Liang, Matthew Argall, Richard Denton, Rachel Rice, Jason Beedle, Kevin Genestreti, Akhtar Ardakani, Adam Stanier, Ari Le , et al. (11 additional authors not shown)

Abstract: We report a rare regime of Earth's magnetosphere interaction with sub-Alfvénic solar wind in which the windsock-like magnetosphere transforms into one with Alfvén wings. In the magnetic cloud of a Coronal Mass Ejection (CME) on April 24, 2023, NASA's Magnetospheric Multiscale mission distinguishes the following features: (1) unshocked and accelerated cold CME plasma coming directly against Earth's… ▽ More We report a rare regime of Earth's magnetosphere interaction with sub-Alfvénic solar wind in which the windsock-like magnetosphere transforms into one with Alfvén wings. In the magnetic cloud of a Coronal Mass Ejection (CME) on April 24, 2023, NASA's Magnetospheric Multiscale mission distinguishes the following features: (1) unshocked and accelerated cold CME plasma coming directly against Earth's dayside magnetosphere; (2) dynamical wing filaments representing new channels of magnetic connection between the magnetosphere and foot points of the Sun's erupted flux rope; (3) cold CME ions observed with energized counter-streaming electrons, evidence of CME plasma captured due to reconnection between magnetic-cloud and Alfvén-wing field lines. The reported measurements advance our knowledge of CME interaction with planetary magnetospheres, and open new opportunities to understand how sub-Alfvénic plasma flows impact astrophysical bodies such as Mercury, moons of Jupiter, and exoplanets close to their host stars. △ Less

Submitted 3 May, 2024; v1 submitted 12 February, 2024; originally announced February 2024.

Comments: 14 pages, including 4 figures, Under review in Geophys. Res. Lett

A New Database of Giant Impacts over a Wide Range of Masses and with Material Strength: A First Analysis of Outcomes

Authors: Alexandre Emsenhuber, Erik Asphaug, Saverio Cambioni, Travis S. J. Gabriel, Stephen R. Schwartz, Robert E. Melikyan, C. Adeene Denton

Abstract: In the late stage of terrestrial planet formation, planets are predicted to undergo pairwise collisions known as giant impacts. Here we present a high-resolution database of giant impacts for differentiated colliding bodies of iron-silicate composition, with target masses ranging from 10^-4 M_Earth up to super-Earths (5 M_Earth). We vary impactor-to-target mass ratio, core-mantle (iron-silicate) f… ▽ More In the late stage of terrestrial planet formation, planets are predicted to undergo pairwise collisions known as giant impacts. Here we present a high-resolution database of giant impacts for differentiated colliding bodies of iron-silicate composition, with target masses ranging from 10^-4 M_Earth up to super-Earths (5 M_Earth). We vary impactor-to-target mass ratio, core-mantle (iron-silicate) fraction, impact velocity, and impact angle. Strength in the form of friction is included in all simulations. We find that due to strength, collisions with bodies smaller than about 2*10^-3 M_Earth can result in irregular shapes, compound core structures, and captured binaries. We observe that the characteristic escaping velocity of smaller remnants (debris) is approximately half of the impact velocity, significantly faster than currently assumed in N-body simulations of planet formation. Incorporating these results in N-body planet formation studies would provide more realistic debris-debris and debris-planet interactions. △ Less

Submitted 30 January, 2024; originally announced January 2024.

Comments: Accepted for publication in PSJ; Table 2 is available in full in an ancillary file

Multi-point Assessment of the Kinematics of Shocks (MAKOS): A Heliophysics Mission Concept Study

Authors: Katherine A. Goodrich, Lynn B. Wilson III, Steven Schwartz, Ian J. Cohen, Drew L. Turner, Phyllis Whittlesey, Amir Caspi, Randall Rose, Keith Smith

Abstract: Collisionless shocks are fundamental processes that are ubiquitous in space plasma physics throughout the Heliosphere and most astrophysical environments. Earth's bow shock and interplanetary shocks at 1 AU offer the most readily accessible opportunities to advance our understanding of the nature of collisionless shocks via fully-instrumented, in situ observations. One major outstanding question p… ▽ More Collisionless shocks are fundamental processes that are ubiquitous in space plasma physics throughout the Heliosphere and most astrophysical environments. Earth's bow shock and interplanetary shocks at 1 AU offer the most readily accessible opportunities to advance our understanding of the nature of collisionless shocks via fully-instrumented, in situ observations. One major outstanding question pertains to the energy budget of collisionless shocks, particularly how exactly collisionless shocks convert incident kinetic bulk flow energy into thermalization (heating), suprathermal particle acceleration, and a variety of plasma waves, including nonlinear structures. Furthermore, it remains unknown how those energy conversion processes change for different shock orientations (e.g., quasi-parallel vs. quasi-perpendicular) and driving conditions (upstream Alfvénic and fast Mach numbers, plasma beta, etc.). Required to address these questions are multipoint observations enabling direct measurement of the necessary plasmas, energetic particles, and electric and magnetic fields and waves, all simultaneously from upstream, downstream, and at the shock transition layer with observatory separations at ion to magnetohydrodynamic (MHD) scales. Such a configuration of spacecraft with specifically-designed instruments has never been available, and this white paper describes a conceptual mission design -- MAKOS -- to address these outstanding questions and advance our knowledge of the nature of collisionless shocks. △ Less

Submitted 8 June, 2023; originally announced June 2023.

Comments: White paper submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033; 9 pages, 3 figures, 5 tables

Journal ref: Bulletin of the AAS, Vol. 55, Issue 3, Whitepaper #135 (9pp); 2023 July 31

The Persistent Mystery of Collisionless Shocks

Authors: Katherine Goodrich, Steven Schwartz, Lynn Wilson III, Ian Cohen, Drew Turner, Amir Caspi, Keith Smith, Randall Rose, Phyllis Whittlesey, Ferdinand Plaschke, Jasper Halekas, George Hospodarsky, James Burch, Imogen Gingell, Li-Jen Chen, Alessandro Retino, Yuri Khotyaintsev

Abstract: Collisionless shock waves are one of the main forms of energy conversion in space plasmas. They can directly or indirectly drive other universal plasma processes such as magnetic reconnection, turbulence, particle acceleration and wave phenomena. Collisionless shocks employ a myriad of kinetic plasma mechanisms to convert the kinetic energy of supersonic flows in space to other forms of energy (e.… ▽ More Collisionless shock waves are one of the main forms of energy conversion in space plasmas. They can directly or indirectly drive other universal plasma processes such as magnetic reconnection, turbulence, particle acceleration and wave phenomena. Collisionless shocks employ a myriad of kinetic plasma mechanisms to convert the kinetic energy of supersonic flows in space to other forms of energy (e.g., thermal plasma, energetic particles, or Poynting flux) in order for the flow to pass an immovable obstacle. The partitioning of energy downstream of collisionless shocks is not well understood, nor are the processes which perform energy conversion. While we, as the heliophysics community, have collected an abundance of observations of the terrestrial bow shock, instrument and mission-level limitations have made it impossible to quantify this partition, to establish the physics within the shock layer responsible for it, and to understand its dependence on upstream conditions. This paper stresses the need for the first ever spacecraft mission specifically designed and dedicated to the observation of both the terrestrial bow shock as well as Interplanetary shocks in the solar wind. △ Less

Submitted 8 June, 2023; originally announced June 2023.

Comments: White paper submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033; 9 pages, 4 figures

Journal ref: Bulletin of the AAS, Vol. 55, Issue 3, Whitepaper #134 (9pp); 2023 July 31

Successful Kinetic Impact into an Asteroid for Planetary Defense

Authors: R. Terik Daly, Carolyn M. Ernst, Olivier S. Barnouin, Nancy L. Chabot, Andrew S. Rivkin, Andrew F. Cheng, Elena Y. Adams, Harrison F. Agrusa, Elisabeth D. Abel, Amy L. Alford, Erik I. Asphaug, Justin A. Atchison, Andrew R. Badger, Paul Baki, Ronald-L. Ballouz, Dmitriy L. Bekker, Julie Bellerose, Shyam Bhaskaran, Bonnie J. Buratti, Saverio Cambioni, Michelle H. Chen, Steven R. Chesley, George Chiu, Gareth S. Collins, Matthew W. Cox , et al. (76 additional authors not shown)

Abstract: While no known asteroid poses a threat to Earth for at least the next century, the catalog of near-Earth asteroids is incomplete for objects whose impacts would produce regional devastation. Several approaches have been proposed to potentially prevent an asteroid impact with Earth by deflecting or disrupting an asteroid. A test of kinetic impact technology was identified as the highest priority sp… ▽ More While no known asteroid poses a threat to Earth for at least the next century, the catalog of near-Earth asteroids is incomplete for objects whose impacts would produce regional devastation. Several approaches have been proposed to potentially prevent an asteroid impact with Earth by deflecting or disrupting an asteroid. A test of kinetic impact technology was identified as the highest priority space mission related to asteroid mitigation. NASA's Double Asteroid Redirection Test (DART) mission is the first full-scale test of kinetic impact technology. The mission's target asteroid was Dimorphos, the secondary member of the S-type binary near-Earth asteroid (65803) Didymos. This binary asteroid system was chosen to enable ground-based telescopes to quantify the asteroid deflection caused by DART's impact. While past missions have utilized impactors to investigate the properties of small bodies those earlier missions were not intended to deflect their targets and did not achieve measurable deflections. Here we report the DART spacecraft's autonomous kinetic impact into Dimorphos and reconstruct the impact event, including the timeline leading to impact, the location and nature of the DART impact site, and the size and shape of Dimorphos. The successful impact of the DART spacecraft with Dimorphos and the resulting change in Dimorphos's orbit demonstrates that kinetic impactor technology is a viable technique to potentially defend Earth if necessary. △ Less

Submitted 3 March, 2023; originally announced March 2023.

Comments: Accepted by Nature

Ejecta from the DART-produced active asteroid Dimorphos

Authors: Jian-Yang Li, Masatoshi Hirabayashi, Tony L. Farnham, Jessica M. Sunshine, Matthew M. Knight, Gonzalo Tancredi, Fernando Moreno, Brian Murphy, Cyrielle Opitom, Steve Chesley, Daniel J. Scheeres, Cristina A. Thomas, Eugene G. Fahnestock, Andrew F. Cheng, Linda Dressel, Carolyn M. Ernst, Fabio Ferrari, Alan Fitzsimmons, Simone Ieva, Stavro L. Ivanovski, Teddy Kareta, Ludmilla Kolokolova, Tim Lister, Sabina D. Raducan, Andrew S. Rivkin , et al. (39 additional authors not shown)

Abstract: Some active asteroids have been proposed to be the result of impact events. Because active asteroids are generally discovered serendipitously only after their tail formation, the process of the impact ejecta evolving into a tail has never been directly observed. NASA's Double Asteroid Redirection Test (DART) mission, apart from having successfully changed the orbital period of Dimorphos, demonstra… ▽ More Some active asteroids have been proposed to be the result of impact events. Because active asteroids are generally discovered serendipitously only after their tail formation, the process of the impact ejecta evolving into a tail has never been directly observed. NASA's Double Asteroid Redirection Test (DART) mission, apart from having successfully changed the orbital period of Dimorphos, demonstrated the activation process of an asteroid from an impact under precisely known impact conditions. Here we report the observations of the DART impact ejecta with the Hubble Space Telescope (HST) from impact time T+15 minutes to T+18.5 days at spatial resolutions of ~2.1 km per pixel. Our observations reveal a complex evolution of ejecta, which is first dominated by the gravitational interaction between the Didymos binary system and the ejected dust and later by solar radiation pressure. The lowest-speed ejecta dispersed via a sustained tail that displayed a consistent morphology with previously observed asteroid tails thought to be produced by impact. The ejecta evolution following DART's controlled impact experiment thus provides a framework for understanding the fundamental mechanisms acting on asteroids disrupted by natural impact. △ Less

Submitted 2 March, 2023; originally announced March 2023.

Comments: accepted by Nature

After DART: Using the first full-scale test of a kinetic impactor to inform a future planetary defense mission

Authors: Thomas S. Statler, Sabina D. Raducan, Olivier S. Barnouin, Mallory E. DeCoster, Steven R. Chesley, Brent Barbee, Harrison F. Agrusa, Saverio Cambioni, Andrew F. Cheng, Elisabetta Dotto, Siegfried Eggl, Eugene G. Fahnestock, Fabio Ferrari, Dawn Graninger, Alain Herique, Isabel Herreros, Masatoshi Hirabayashi, Stavro Ivanovski, Martin Jutzi, Özgür Karatekin, Alice Lucchetti, Robert Luther, Rahil Makadia, Francesco Marzari, Patrick Michel , et al. (16 additional authors not shown)

Abstract: NASA's Double Asteroid Redirection Test (DART) is the first full-scale test of an asteroid deflection technology. Results from the hypervelocity kinetic impact and Earth-based observations, coupled with LICIACube and the later Hera mission, will result in measurement of the momentum transfer efficiency accurate to ~10% and characterization of the Didymos binary system. But DART is a single experim… ▽ More NASA's Double Asteroid Redirection Test (DART) is the first full-scale test of an asteroid deflection technology. Results from the hypervelocity kinetic impact and Earth-based observations, coupled with LICIACube and the later Hera mission, will result in measurement of the momentum transfer efficiency accurate to ~10% and characterization of the Didymos binary system. But DART is a single experiment; how could these results be used in a future planetary defense necessity involving a different asteroid? We examine what aspects of Dimorphos's response to kinetic impact will be constrained by DART results; how these constraints will help refine knowledge of the physical properties of asteroidal materials and predictive power of impact simulations; what information about a potential Earth impactor could be acquired before a deflection effort; and how design of a deflection mission should be informed by this understanding. We generalize the momentum enhancement factor $β$, showing that a particular direction-specific $β$ will be directly determined by the DART results, and that a related direction-specific $β$ is a figure of merit for a kinetic impact mission. The DART $β$ determination constrains the ejecta momentum vector, which, with hydrodynamic simulations, constrains the physical properties of Dimorphos's near-surface. In a hypothetical planetary defense exigency, extrapolating these constraints to a newly discovered asteroid will require Earth-based observations and benefit from in-situ reconnaissance. We show representative predictions for momentum transfer based on different levels of reconnaissance and discuss strategic targeting to optimize the deflection and reduce the risk of a counterproductive deflection in the wrong direction. △ Less

Submitted 23 September, 2022; originally announced September 2022.

Comments: 30 pages, 7 figures. Planetary Science Journal, in press, accepted 2022 September 22

Effects of impact and target parameters on the results of a kinetic impactor: predictions for the Double Asteroid Redirection Test (DART) mission

Authors: Angela M. Stickle, Mallory E. DeCoster, Christoph Burger, Wendy K. Caldwell, Dawn Graninger, Kathryn M. Kumamoto, Robert Luther, Jens Ormö, Sabina Raducan, Emma Rainey, Christoph M. Schäfer, James D. Walker, Yun Zhang, Patrick Michel, J. Michael Owen, Olivier Barnouin, Andy F. Cheng, Sidney Cochron, Gareth S. Collins, Thomas M. Davison, Elisabetta Dotto, Fabio Ferrari, M. Isabel Herreros, Stavro L. Ivanovski, Martin Jutzi , et al. (8 additional authors not shown)

Abstract: The Double Asteroid Redirection Test (DART) spacecraft will impact into the asteroid Dimorphos on September 26, 2022 as a test of the kinetic impactor technique for planetary defense. The efficiency of the deflection following a kinetic impactor can be represented using the momentum enhancement factor, Beta, which is dependent on factors such as impact geometry and the specific target material pro… ▽ More The Double Asteroid Redirection Test (DART) spacecraft will impact into the asteroid Dimorphos on September 26, 2022 as a test of the kinetic impactor technique for planetary defense. The efficiency of the deflection following a kinetic impactor can be represented using the momentum enhancement factor, Beta, which is dependent on factors such as impact geometry and the specific target material properties. Currently, very little is known about Dimorphos and its material properties that introduces uncertainty in the results of the deflection efficiency observables, including crater formation, ejecta distribution, and Beta. The DART Impact Modeling Working Group (IWG) is responsible for using impact simulations to better understand the results of the DART impact. Pre-impact simulation studies also provide considerable insight into how different properties and impact scenarios affect momentum enhancement following a kinetic impact. This insight provides a basis for predicting the effects of the DART impact and the first understanding of how to interpret results following the encounter. Following the DART impact, the knowledge gained from these studies will inform the initial simulations that will recreate the impact conditions, including providing estimates for potential material properties of Dimorphos and Beta resulting from DARTs impact. This paper summarizes, at a high level, what has been learned from the IWG simulations and experiments in preparation for the DART impact. While unknown, estimates for reasonable potential material properties of Dimorphos provide predictions for Beta of 1-5, depending on end-member cases in the strength regime. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: Accepted to PSJ Didymos-DART Focus Issue

Predictions for the Dynamical States of the Didymos System before and after the Planned DART Impact

Authors: Derek C. Richardson, Harrison F. Agrusa, Brent Barbee, William F. Bottke, Andrew F. Cheng, Siegfried Eggl, Fabio Ferrari, Masatoshi Hirabayashi, Özgür Karatekin, Jay McMahon, Stephen R. Schwartz, Ronald-Louis Ballouz, Adriano Campo Bagatin, Elisabetta Dotto, Eugene G. Fahnestock, Oscar Fuentes-Muñoz, Ioannis Gkolias, Douglas P. Hamilton, Seth A. Jacobson, Martin Jutzi, Josh Lyzhoft, Rahil Makadia, Alex J. Meyer, Patrick Michel, Ryota Nakano , et al. (11 additional authors not shown)

Abstract: NASA's Double Asteroid Redirection Test (DART) spacecraft is planned to impact the natural satellite of (65803) Didymos, Dimorphos, around 23:14 UTC on 26 September 2022, causing a reduction in its orbital period that will be measurable with ground-based observations. This test of kinetic impactor technology will provide the first estimate of the momentum transfer enhancement factor $β$ at a reali… ▽ More NASA's Double Asteroid Redirection Test (DART) spacecraft is planned to impact the natural satellite of (65803) Didymos, Dimorphos, around 23:14 UTC on 26 September 2022, causing a reduction in its orbital period that will be measurable with ground-based observations. This test of kinetic impactor technology will provide the first estimate of the momentum transfer enhancement factor $β$ at a realistic scale, wherein ejecta from the impact provides an additional deflection to the target. Earth-based observations, the LICIACube spacecraft (to be detached from DART prior to impact), and ESA's follow-up Hera mission to launch in 2024, will provide additional characterization of the deflection test. Together Hera and DART comprise the Asteroid Impact and Deflection Assessment (AIDA) cooperation between NASA and ESA. Here the predicted dynamical states of the binary system upon arrival and after impact are presented. The assumed dynamically relaxed state of the system will be excited by the impact, leading to an increase in eccentricity and slight tilt of the orbit together with enhanced libration of Dimorphos with amplitude dependent on the currently poorly known target shape. Free rotation around the moon's long axis may also be triggered and the orbital period will experience variations from seconds to minutes over timescales of days to months. Shape change of either body due to cratering or mass wasting triggered by crater formation and ejecta may affect $β$ but can be constrained through additional measurements. Both BYORP and gravity tides may cause measurable orbital changes on the timescale of Hera's rendezvous. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 23 pages, 13 figures, published in PSJ

Journal ref: Planet. Sci. J. 3 157 (2022)

Collision Chains among the Terrestrial Planets. III. Formation of the Moon

Authors: Erik Asphaug, Alexandre Emsenhuber, Saverio Cambioni, Travis S. J. Gabriel, Stephen R. Schwartz

Abstract: In the canonical model of Moon formation, a Mars-sized protoplanet "Theia" collides with proto-Earth at close to their mutual escape velocity $v_{\rm esc}$ and a common impact angle 45°. The "graze-and-merge" collision strands a fraction of Theia's mantle into orbit, while Earth accretes most of Theia and its momentum. Simulations show that this produces a hot, high angular momentum, silicate-domi… ▽ More In the canonical model of Moon formation, a Mars-sized protoplanet "Theia" collides with proto-Earth at close to their mutual escape velocity $v_{\rm esc}$ and a common impact angle 45°. The "graze-and-merge" collision strands a fraction of Theia's mantle into orbit, while Earth accretes most of Theia and its momentum. Simulations show that this produces a hot, high angular momentum, silicate-dominated protolunar system, in substantial agreement with lunar geology, geochemistry, and dynamics. However, a Moon that derives mostly from Theia's mantle, as angular momentum dictates, is challenged by the fact that O, Ti, Cr, radiogenic W, and other elements are indistinguishable in Earth and lunar rocks. Moreover, the model requires an improbably low initial velocity. Here we develop a scenario for Moon formation that begins with a somewhat faster collision, when proto-Theia impacts proto-Earth at ~1.2 $v_{\rm esc}$, also around 45°. Instead of merging, the bodies come into violent contact for a half-hour and their major components escape, a "hit-and-run collision." N-body evolutions show that the "runner" often returns ~0.1-1 Myr later for a second giant impact, closer to $v_{\rm esc}$; this produces a postimpact disk of ~2-3 lunar masses in smoothed particle hydrodynamics simulations, with angular momentum comparable to canonical scenarios. The disk ends up substantially inclined, in most cases, because the terminal collision is randomly oriented to the first. Proto-Earth contributions to the silicate disk are enhanced by the compounded mixing and greater energy of a collision chain. △ Less

Submitted 1 October, 2021; originally announced October 2021.

Comments: 27 pages, 4 figures, 2 tables. Published in PSJ, article available at https://iopscience.iop.org/article/10.3847/PSJ/ac19b2

Journal ref: Planet. Sci. J. 2:200 (2021)

Collision Chains among the Terrestrial Planets. II. An Asymmetry between Earth and Venus

Authors: Alexandre Emsenhuber, Erik Asphaug, Saverio Cambioni, Travis S. J. Gabriel, Stephen R. Schwartz

Abstract: During the late stage of terrestrial planet formation, hit-and-run collisions are about as common as accretionary mergers, for expected velocities and angles of giant impacts. Average hit-and-runs leave two major remnants plus debris: the target and impactor, somewhat modified through erosion, escaping at lower relative velocity. Here we continue our study of the dynamical effects of such collisio… ▽ More During the late stage of terrestrial planet formation, hit-and-run collisions are about as common as accretionary mergers, for expected velocities and angles of giant impacts. Average hit-and-runs leave two major remnants plus debris: the target and impactor, somewhat modified through erosion, escaping at lower relative velocity. Here we continue our study of the dynamical effects of such collisions. We compare the dynamical fates of intact runners that start from hit-and-runs with proto-Venus at 0.7 AU and proto-Earth at 1.0 AU. We follow the orbital evolutions of the runners, including the other terrestrial planets, Jupiter, and Saturn, in an N-body code. We find that the accretion of these runners can take $\gtrsim$10 Myr (depending on the egress velocity of the first collision) and can involve successive collisions with the original target planet or with other planets. We treat successive collisions that the runner experiences using surrogate models from machine learning, as in previous work, and evolve subsequent hit-and-runs in a similar fashion. We identify asymmetries in the capture, loss, and interchange of runners in the growth of Venus and Earth. Hit-and-run is a more probable outcome at proto-Venus, being smaller and faster orbiting than proto-Earth. But Venus acts as a sink, eventually accreting most of its runners, assuming typical events, whereas proto-Earth loses about half, many of those continuing to Venus. This leads to a disparity in the style of late-stage accretion that could have led to significant differences in geology, composition, and satellite formation at Earth and Venus. △ Less

Submitted 1 October, 2021; originally announced October 2021.

Comments: 18 pages, 9 figures, 4 tables. Published in PSJ, article available at https://iopscience.iop.org/article/10.3847/PSJ/ac19b1

Journal ref: Planet. Sci. J. 2:199 (2021)

Ricochets on Asteroids II: Sensitivity of laboratory experiments of low velocity grazing impacts on substrate grain size

Authors: Esteban Wright, Alice C. Quillen, Paul Sanchez, Stephen R. Schwartz, Miki Nakajima, Hesam Askari, Peter Miklavcic

Abstract: We compare low velocity impacts that ricochet with the same impact velocity and impact angle into granular media with similar bulk density, porosity and friction coefficient but different mean grain size. The ratio of projectile diameter to mean grain length ranges from 4 in our coarsest medium to 50 in our finest sand. Using high speed video and fluorescent markers, we measure the ratio of pre- t… ▽ More We compare low velocity impacts that ricochet with the same impact velocity and impact angle into granular media with similar bulk density, porosity and friction coefficient but different mean grain size. The ratio of projectile diameter to mean grain length ranges from 4 in our coarsest medium to 50 in our finest sand. Using high speed video and fluorescent markers, we measure the ratio of pre- to post-impact horizontal and vertical velocity components, which we refer to as coefficients of restitution, and the angle of deflection caused by the impact in the horizontal plane. Coefficients of restitution are sensitive to mean grain size with the ratio associated with the horizontal velocity component about twice as large for our coarsest gravel as that for our finest sand. This implies that coefficients for hydro-static-like, drag-like and lift-like forces, used in empirical force laws, are sensitive to mean grain size. The coefficient that is most strongly sensitive to grain size is the lift coefficient which decreases by a factor of 3 between our coarsest and finest media. The deflection angles are largest in the coarser media and their size approximately depends on grain size to the 3/2 power. This scaling is matched with a model where momentum transfer takes place via collisions with individual grains. The dependence of impact mechanics on substrate size distribution should be considered in future models for populations of objects that impact granular asteroid surfaces. △ Less

Submitted 22 January, 2022; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: Link to the high speed videos for the 4 experiments tracked and discussed: https://youtube.com/playlist?list=PLYnpHw29go3pCWwm554fwR6jNFvB5qVQn Submitted to the journal Icarus on September 12, 2021

Direct Evidence for Magnetic Reflection of Heavy Ions from High Mach Number Collisionless Shocks

Authors: Hadi Madanian, Steven J. Schwartz, Stephen A. Fuselier, David Burgess, Drew L. Turner, Li-Jen Chen, Mihir I. Desai, Michael J. Starkey

Abstract: Strong shocks in collisionless plasmas, such as supernovae shocks and shocks driven by coronal mass ejections, are known to be a primary source of energetic particles. Due to their different mass per charge ratio, the interaction of heavy ions with the shock layer differs from that of protons, and injection of these ions into acceleration processes is a challenge. Here we show the first direct obs… ▽ More Strong shocks in collisionless plasmas, such as supernovae shocks and shocks driven by coronal mass ejections, are known to be a primary source of energetic particles. Due to their different mass per charge ratio, the interaction of heavy ions with the shock layer differs from that of protons, and injection of these ions into acceleration processes is a challenge. Here we show the first direct observational evidence of magnetic reflection of alpha particles from a high Mach number quasi-perpendicular shock using in-situ spacecraft measurements. The intense magnetic amplification at the shock front associated with nonstationarity modulates the trajectory of alpha particles, some of which travel back upstream as they gyrate in the enhanced magnetic field and experience further acceleration in the upstream region. Our results in particular highlight the important role of high magnetic amplification in seeding heavy ions into the energization processes at nonstationary reforming shocks. △ Less

Submitted 20 August, 2021; v1 submitted 18 June, 2021; originally announced June 2021.

Journal ref: ApJL 915 L19 (2021)

The Effect of Inefficient Accretion on Planetary Differentiation

Authors: Saverio Cambioni, Seth A. Jacobson, Alexandre Emsenhuber, Erik Asphaug, David C. Rubie, Travis S. J. Gabriel, Stephen R. Schwartz, Roberto Furfaro

Abstract: Pairwise collisions between terrestrial embryos are the dominant means of accretion during the last stage of planet formation. Hence, their realistic treatment in N-body studies is critical to accurately model the formation of terrestrial planets and to develop interpretations of telescopic and spacecraft observations. In this work, we compare the effects of two collision prescriptions on the core… ▽ More Pairwise collisions between terrestrial embryos are the dominant means of accretion during the last stage of planet formation. Hence, their realistic treatment in N-body studies is critical to accurately model the formation of terrestrial planets and to develop interpretations of telescopic and spacecraft observations. In this work, we compare the effects of two collision prescriptions on the core-mantle differentiation of terrestrial planets: a model in which collisions are always completely accretionary (``perfect merging'') and a more realistic model based on neural networks that has been trained on hydrodynamical simulations of giant impacts. The latter model is able to predict the loss of mass due to imperfect accretion and the evolution of non-accreted projectiles in hit-and-run collisions. We find that the results of the neural-network model feature a wider range of final core mass fractions and metal-silicate equilibration pressures, temperatures, and oxygen fugacities than the assumption of perfect merging. When used to model collisions in N-body studies of terrestrial planet formation, the two models provide similar answers for planets more massive than 0.1 Earth's masses. For less massive final bodies, however, the inefficient-accretion model predicts a higher degree of compositional diversity. This phenomenon is not reflected in planet formation models of the solar system that use perfect merging to determine collisional outcomes. Our findings confirm the role of giant impacts as important drivers of planetary diversity and encourage a realistic implementation of inefficient accretion in future accretion studies. △ Less

Submitted 14 June, 2021; originally announced June 2021.

Comments: 21 pages, 2 tables, 7 figures. Published open access on PSJ: https://iopscience.iop.org/article/10.3847/PSJ/abf0ad

Journal ref: Planet. Sci. J. 2 93 (2021)

Direct Multipoint Observations Capturing the Reformation of a Supercritical Fast Magnetosonic Shock

Authors: D. L. Turner, L. B. Wilson III, K. A. Goodrich, H. Madanian, S. J. Schwartz, T. Z. Liu, A. Johlander, D. Caprioli, I. J. Cohen, D. Gershman, H. Hietala, J. H. Westlake, B. Lavraud, O. Le Contel, J. L. Burch

Abstract: Using multipoint Magnetospheric Multiscale (MMS) observations in an unusual string-of-pearls configuration, we examine in detail observations of the reformation of a fast magnetosonic shock observed on the upstream edge of a foreshock transient structure upstream of Earth's bow shock. The four MMS spacecraft were separated by several hundred km, comparable to suprathermal ion gyro-radius scales or… ▽ More Using multipoint Magnetospheric Multiscale (MMS) observations in an unusual string-of-pearls configuration, we examine in detail observations of the reformation of a fast magnetosonic shock observed on the upstream edge of a foreshock transient structure upstream of Earth's bow shock. The four MMS spacecraft were separated by several hundred km, comparable to suprathermal ion gyro-radius scales or several ion inertial lengths. At least half of the shock reformation cycle was observed, with a new shock ramp rising up out of the "foot" region of the original shock ramp. Using the multipoint observations, we convert the observed time-series data into distance along the shock normal in the shock's rest frame. That conversion allows for a unique study of the relative spatial scales of the shock's various features, including the shock's growth rate, and how they evolve during the reformation cycle. Analysis indicates that: the growth rate increases during reformation, electron-scale physics play an important role in the shock reformation, and energy conversion processes also undergo the same cyclical periodicity as reformation. Strong, thin electron-kinetic-scale current sheets and large-amplitude electrostatic and electromagnetic waves are reported. Results highlight the critical cross-scale coupling between electron-kinetic- and ion-kinetic-scale processes and details of the nature of nonstationarity, shock-front reformation at collisionless, fast magnetosonic shocks. △ Less

Submitted 2 April, 2021; originally announced April 2021.

Comments: In press and to be published in Astrophysical Journal Letters in 2021

Validating N-body code Chrono for granular DEM simulations in reduced-gravity environments

Authors: Cecily Sunday, Naomi Murdoch, Simon Tardivel, Stephen R. Schwartz, Patrick Michel

Abstract: The Discrete Element Method (DEM) is frequently used to model complex granular systems and to augment the knowledge that we obtain through theory, experimentation, and real-world observations. Numerical simulations are a particularly powerful tool for studying the regolith-covered surfaces of asteroids, comets, and small moons, where reduced-gravity environments produce ill-defined flow behaviors.… ▽ More The Discrete Element Method (DEM) is frequently used to model complex granular systems and to augment the knowledge that we obtain through theory, experimentation, and real-world observations. Numerical simulations are a particularly powerful tool for studying the regolith-covered surfaces of asteroids, comets, and small moons, where reduced-gravity environments produce ill-defined flow behaviors. In this work, we present a method for validating soft-sphere DEM codes for both terrestrial and small-body granular environments. The open-source code Chrono is modified and evaluated first with a series of simple two-body-collision tests, and then, with a set of piling and tumbler tests. In the piling tests, we vary the coefficient of rolling friction to calibrate the simulations against experiments with 1 mm glass beads. Then, we use the friction coefficient to model the flow of 1 mm glass beads in a rotating drum, using a drum configuration from a previous experimental study. We measure the dynamic angle of repose, the flowing layer thickness, and the flowing layer velocity for tests with different particle sizes, contact force models, coefficients of rolling friction, cohesion levels, drum rotation speeds and gravity levels. The tests show that the same flow patterns can be observed at Earth and reduced-gravity levels if the drum rotation speed and the gravity-level are set according to the dimensionless parameter known as the Froude number. Chrono is successfully validated against known flow behaviors at different gravity and cohesion levels, and will be used to study small-body regolith dynamics in future works. △ Less

Submitted 22 September, 2020; originally announced September 2020.

Journal ref: Monthly Notices of the Royal Astronomical Society 498.1 (2020): 1062-1079

Proposed Miami Impact Crater Identified as a Solutional Doline of Oolictic Limestone

Authors: Antonio Paris, Ryan Robinson, Skye Schwartz

Abstract: This investigation addresses the discovery of a proposed impact crater located off the coast of Miami, FL under the North Atlantic Ocean. A preliminary analysis of bathymetry data obtained from the National Oceanic and Atmospheric Administration implied a morphology consistent with a complex crater produced by a hypervelocity impact event of extraterrestrial origin. The proposed impact features in… ▽ More This investigation addresses the discovery of a proposed impact crater located off the coast of Miami, FL under the North Atlantic Ocean. A preliminary analysis of bathymetry data obtained from the National Oceanic and Atmospheric Administration implied a morphology consistent with a complex crater produced by a hypervelocity impact event of extraterrestrial origin. The proposed impact features include a central peak, concentric rings, and an ejecta field to the northwest. Analysis of geological data from the US Geological Survey places the strata overlying the proposed impact site as Miami Limestone, Pleistocene, accumulated during Marine Isotope Stage 5e, thereby placing the maximum age of the proposed impact crater at 80 ka to 130 ka. Three other competing hypotheses for the formation of the structure, namely a controlled maritime explosion, radial lava flow from volcano, or a depressed bioherm, doline, or karst (solutional depression) were explored throughout the investigation. To confirm the proposed structure as an impact crater, an in-situ underwater expedition was organized by Planetary Sciences, Inc. specifically to ascertain whether planar formations, shatter cones, and shock metamorphic and or other meteoritic properties were present. After analyzing the geological samples collected at the proposed impact crater, examining the morphology of analogous geologic structures, and evaluating competing hypotheses, we conclude that the structure is a solutional doline formed by the uneven dissolution of the Miami Limestone, and, accordingly, do not recommend that the structure be indexed in the Earth Impact Database △ Less

Submitted 19 July, 2020; v1 submitted 13 July, 2020; originally announced July 2020.

Comments: 10 Pages, 11 Figures

Journal ref: Journal of the Washington Academy of Sciences, Fall 2020, Volume 106, No. 3

Ricochets on Asteroids: Experimental study of low velocity grazing impacts into granular media

Authors: Esteban Wright, Alice Quillen, Juiliana South, Randal C. Nelson, Paul Sanchez, John Siu, Hesam Askari, Miki Nakajima, Stephen R. Schwartz

Abstract: Spin off events and impacts can eject boulders from an asteroid surface and rubble pile asteroids can accumulate from debris following a collision between large asteroids. These processes produce a population of gravitational bound objects in orbit that can impact an asteroid surface at low velocity and with a distribution of impact angles. We present laboratory experiments of low velocity spheric… ▽ More Spin off events and impacts can eject boulders from an asteroid surface and rubble pile asteroids can accumulate from debris following a collision between large asteroids. These processes produce a population of gravitational bound objects in orbit that can impact an asteroid surface at low velocity and with a distribution of impact angles. We present laboratory experiments of low velocity spherical projectiles into a fine granular medium, sand. We delineate velocity and impact angles giving ricochets, those giving projectiles that roll-out from the impact crater and those that stop within their impact crater. With high speed camera images and fluorescent markers on the projectiles we track spin and projectile trajectories during impact. We find that the projectile only reaches a rolling without slipping condition well after the marble has reached peak penetration depth. The required friction coefficient during the penetration phase of impact is 4-5 times lower than that of the sand suggesting that the sand is fluidized near the projectile surface during penetration. We find that the critical grazing impact critical angle dividing ricochets from roll-outs, increases with increasing impact velocity. The critical angles for ricochet and for roll-out as a function of velocity can be matched by an empirical model during the rebound phase that balances a lift force against gravity. We estimate constraints on projectile radius, velocity and impact angle that would allow projectiles on asteroids to ricochet or roll away from impact, finally coming to rest distant from their initial impact sites. △ Less

Submitted 19 June, 2020; v1 submitted 4 February, 2020; originally announced February 2020.

Comments: Link to the high speed videos for the 3 experiments tracked and discussed: https://www.youtube.com/playlist?list=PLYnpHw29go3rKYOvM4WEQkHWYYIGgRC-E Preliminary results presented at American Geophysical Union (AGU) 2019 Fall Meeting (P43D-3497). Submitted to the journal Icarus on February 3, 2020

Electron energy partition across interplanetary shocks: III. Analysis

Authors: L. B. Wilson III, Li-Jen Chen, Shan Wang, Steven J. Schwartz, Drew L. Turner, Michael L. Stevens, Justin C. Kasper, Adnane Osmane, Damiano Caprioli, Stuart D. Bale, Marc P. Pulupa, Chadi S. Salem, Katherine A. Goodrich

Abstract: Analysis of model fit results of 15,210 electron velocity distribution functions (VDFs), observed within $\pm$2 hours of 52 interplanetary (IP) shocks by the Wind spacecraft near 1 AU, is presented as the third and final part on electron VDFs near IP shocks. The core electrons and protons dominate in the magnitude and change in the partial-to-total thermal pressure ratio, with the core electrons o… ▽ More Analysis of model fit results of 15,210 electron velocity distribution functions (VDFs), observed within $\pm$2 hours of 52 interplanetary (IP) shocks by the Wind spacecraft near 1 AU, is presented as the third and final part on electron VDFs near IP shocks. The core electrons and protons dominate in the magnitude and change in the partial-to-total thermal pressure ratio, with the core electrons often gaining as much or more than the protons. Only a moderate positive correlation is observed between the electron temperature and the kinetic energy change across the shock, while weaker, if any, correlations were found with any other macroscopic shock parameter. No VDF parameter correlated with the shock normal angle. The electron VDF evolves from a narrowly peaked core with flaring suprathermal tails in the upstream to either a slightly hotter core with steeper tails or much hotter flattop core with even steeper tails downstream of the weaker and strongest shocks, respectively. Both quasi-static and fluctuating fields are examined as possible mechanisms modifying the VDF but neither is sufficient alone. For instance, flattop VDFs can be generated by nonlinear ion acoustic wave stochastic acceleration (i.e., inelastic collisions) while other work suggested they result from the combination of quasi-static and fluctuating fields. This three-part study shows that not only are these systems not thermodynamic in nature, even kinetic models may require modification to include things like inelastic collision operators to properly model electron VDF evolution across shocks or in the solar wind. △ Less

Submitted 24 January, 2020; originally announced January 2020.

Comments: 30 pages, 10 figures, 3 tables, to be submitted to Astrophys. J

Realistic On-the-fly Outcomes of Planetary Collisions II: Bringing Machine Learning to N-body Simulations

Authors: Alexandre Emsenhuber, Saverio Cambioni, Erik Asphaug, Travis S. J. Gabriel, Stephen R. Schwartz, Roberto Furfaro

Abstract: Terrestrial planet formation theory is at a bottleneck, with the growing realization that pairwise collisions are treated far too simply. Here, and in our companion paper (Cambioni et al. 2019) that introduces the training methodology, we demonstrate the first application of machine learning to more realistically model the late stage of planet formation by giant impacts. We present surrogate model… ▽ More Terrestrial planet formation theory is at a bottleneck, with the growing realization that pairwise collisions are treated far too simply. Here, and in our companion paper (Cambioni et al. 2019) that introduces the training methodology, we demonstrate the first application of machine learning to more realistically model the late stage of planet formation by giant impacts. We present surrogate models that give fast, reliable answers for the masses and velocities of the two largest remnants of a giant impact, as a function of the colliding masses and their impact velocity and angle, with the caveat that our training data do not yet include pre-impact rotation or variable thermal conditions. We compare canonical N-body scenarios of terrestrial planet formation assuming perfect merger (Chambers 2001) with our more realistic treatment that includes inefficient accretions and hit-and-run collisions. The result is a protracted tail of final events lasting ~200 Myr, and the conversion of about half the mass of the initial population to debris. We obtain profoundly different solar system architectures, featuring a much wider range of terrestrial planet masses and enhanced compositional diversity. △ Less

Submitted 3 January, 2020; originally announced January 2020.

Comments: 20 pages, 10 figures, 3 tables; accepted for publication in ApJ; Table 1 is available in full in an ancillary file; the code developed in this work is available at https://github.com/aemsenhuber/collresolve

Journal ref: ApJ 891:6 (2020)

Investigating Asteroid Surface Geophysics with an Ultra-Low-Gravity Centrifuge in Low-Earth Orbit

Authors: Stephen R. Schwartz, Jekan Thangavelautham, Erik Asphaug, Aman Chandra, Ravi teja Nallapu, Leonard Vance

Abstract: Near-Earth small-body mission targets 162173 Ryugu, 101955 Bennu, and 25143 Itokawa produce gravity fields around 4 orders of magnitude below that of Earth and their irregular shapes, combined with rotational effects produce varying surface potentials. Still, we observe familiar geologic textures and landforms that are the result of the granular physical processes that take place on their surfaces… ▽ More Near-Earth small-body mission targets 162173 Ryugu, 101955 Bennu, and 25143 Itokawa produce gravity fields around 4 orders of magnitude below that of Earth and their irregular shapes, combined with rotational effects produce varying surface potentials. Still, we observe familiar geologic textures and landforms that are the result of the granular physical processes that take place on their surfaces. The nature of these landforms, however, their origins, and how these surfaces react to interrogation by probes, landers, rovers, and penetrators, remain largely unknown, and therefore landing on an asteroid and manipulating its surface material remains a daunting challenge. The AOSAT+ design is a 12U CubeSat that will be in Low-Earth Orbit (LEO) and that will operate as a spinning on-orbit centrifuge. Part of this 12U CubeSat will contain a laboratory that will recreate asteroid surface conditions using crushed meteorite as a regolith proxy. The spinning of the laboratory will simulate the surface gravity of asteroids 2 km and smaller. The result is a bed of realistic regolith, the environment that landers and diggers and maybe astronauts will interact with. A crucial component of this mission involves the reproduction of the experimental results in numerical simulation in order to extract the material parameters of the regolith and its behavior in a sustained, very low but nonzero-gravity environment. △ Less

Submitted 10 October, 2019; originally announced October 2019.

Comments: 7 pages, 2 figures, International Astronautical Congress 2019

Electron energy partition across interplanetary shocks: II. Statistics

Authors: Lynn B. Wilson III, Li-Jen Chen, Shan Wang, Steven J. Schwartz, Drew L. Turner, Michael L. Stevens, Justin C. Kasper, Adnane Osmane, Damiano Caprioli, Stuart D. Bale, Marc P. Pulupa, Chadi S. Salem, Katherine A. Goodrich

Abstract: A statistical analysis of 15,210 electron velocity distribution function (VDF) fits, observed within $\pm$2 hours of 52 interplanetary (IP) shocks by the $Wind$ spacecraft near 1 AU, is presented. This is the second in a three-part series on electron VDFs near IP shocks. The electron velocity moment statistics for the dense, low energy core, tenuous, hot halo, and field-aligned beam/strahl are a s… ▽ More A statistical analysis of 15,210 electron velocity distribution function (VDF) fits, observed within $\pm$2 hours of 52 interplanetary (IP) shocks by the $Wind$ spacecraft near 1 AU, is presented. This is the second in a three-part series on electron VDFs near IP shocks. The electron velocity moment statistics for the dense, low energy core, tenuous, hot halo, and field-aligned beam/strahl are a statistically significant list of values illustrated with both histograms and tabular lists for reference and baselines in future work. The beam/strahl fit results in the upstream are currently the closest thing to a proper parameterization of the beam/strahl electron velocity moments in the ambient solar wind. This work will also serve as a 1 AU baseline and reference for missions like $Parker \ Solar \ Probe$ and $Solar \ Orbiter$. The median density, temperature, beta, and temperature anisotropy values for the core(halo)[beam/strahl] components, with subscripts $ec$($eh$)[$eb$], of all fit results respectively are $n{\scriptstyle_{ec(h)[b]}}$ $\sim$ 11.3(0.36)[0.17] $cm^{-3}$, $T{\scriptstyle_{ec(h)[b], tot}}$ $\sim$ 14.6(48.4)[40.2] $eV$, $β{\scriptstyle_{ec(h)[b], tot}}$ $\sim$ 0.93(0.11)[0.05], and $\mathcal{A}{\scriptstyle_{ec(h)[b]}}$ $\sim$ 0.98(1.03)[0.93]. The nuanced details of the fitting method and data product description were published in Paper I and the detailed analysis of the results will be shown in Paper III. △ Less

Submitted 19 September, 2019; originally announced September 2019.

Comments: 40 pages, 11 figures, 18 tables, to be submitted to Astrophys. J. Suppl

Small Solar System Bodies as granular media

Authors: D. Hestroffer, P. Sánchez, L. Staron, A. Campo Bagatin, S. Eggl, W. Losert, N. Murdoch, E. Opsomer, F. Radjai, D. C. Richardson, M. Salazar, D. J. Scheeres, S. Schwartz, N. Taberlet, H. Yano

Abstract: Asteroids and other Small Solar System Bodies (SSSBs) are of high general and scientific interest in many aspects. The origin, formation, and evolution of our Solar System (and other planetary systems) can be better understood by analysing the constitution and physical properties of small bodies in the Solar System. Currently, two space missions (Hayabusa2, OSIRIS-REx) have recently arrived at the… ▽ More Asteroids and other Small Solar System Bodies (SSSBs) are of high general and scientific interest in many aspects. The origin, formation, and evolution of our Solar System (and other planetary systems) can be better understood by analysing the constitution and physical properties of small bodies in the Solar System. Currently, two space missions (Hayabusa2, OSIRIS-REx) have recently arrived at their respective targets and will bring a sample of the asteroids back to Earth. Other small body missions have also been selected by, or proposed to, space agencies. The threat posed to our planet by near-Earth objects (NEOs) is also considered at the international level, and this has prompted dedicated research on possible mitigation techniques. The DART mission, for example, will test the kinetic impact technique. Even ideas for industrial exploitation have risen during the last years. Lastly, the origin of water and life on Earth appears to be connected to asteroids. Hence, future space mission projects will undoubtedly target some asteroids or other SSSBs. In all these cases and research topics, specific knowledge of the structure and mechanical behaviour of the surface as well as the bulk of those celestial bodies is crucial. In contrast to large telluric planets and dwarf planets, a large proportion of such small bodies is believed to consist of gravitational aggregates ('rubble piles') with no -- or low -- internal cohesion, with varying macro-porosity and surface properties (from smooth regolith covered terrain, to very rough collection of boulders), and varying topography (craters, depressions, ridges) [...]. △ Less

Submitted 4 July, 2019; originally announced July 2019.

Comments: This is a pre-print version of an article to be published in AARv, and available online at https://rdcu.be/bHM9E

The Western Bulge of 162173 Ryugu Formed as a Result of a Rotationally Driven Deformation Process

Authors: Masatoshi Hirabayashi, Eri Tatsumi, Hideaki Miyamoto, Goro Komatsu, Seiji Sugita, Sei-ichiro Watanabe, Daniel J. Scheeres, Olivier S. Barnouin, Patrick Michel, Chikatoshi Honda, Tatsuhiro Michikami, Yuichiro Cho, Tomokatsu Morota, Naru Hirata, Naoyuki Hirata, Naoya Sakatani, Stephen R. Schwartz, Rie Honda, Yasuhiro Yokota, Shingo Kameda, Hidehiko Suzuki, Toru Kouyama, Masahiko Hayakawa, Moe Matsuoka, Kazuo Yoshioka , et al. (4 additional authors not shown)

Abstract: 162173 Ryugu, the target of Hayabusa2, has a round shape with an equatorial ridge, which is known as a spinning top-shape. A strong centrifugal force is a likely contributor to Ryugu's top-shaped features. Observations by Optical Navigation Camera onboard Hayabusa2 show a unique longitudinal variation in geomorphology; the western side of this asteroid, later called the western bulge, has a smooth… ▽ More 162173 Ryugu, the target of Hayabusa2, has a round shape with an equatorial ridge, which is known as a spinning top-shape. A strong centrifugal force is a likely contributor to Ryugu's top-shaped features. Observations by Optical Navigation Camera onboard Hayabusa2 show a unique longitudinal variation in geomorphology; the western side of this asteroid, later called the western bulge, has a smooth surface and a sharp equatorial ridge, compared to the other side. Here, we propose a structural deformation process that generated the western bulge. Applying the mission-derived shape model, we employ a finite element model technique to analyze the locations that experience structural failure within the present shape. Assuming that materials are uniformly distributed, our model shows the longitudinal variation in structurally failed regions when the spin period is shorter than ~3.75 h. Ryugu is structurally intact in the subsurface region of the western bulge while other regions are sensitive to structural failure. We infer that this variation is indicative of the deformation process that occurred in the past, and the western bulge is more relaxed structurally than the other region. Our analysis also shows that this deformation process might occur at a spin period between ~3.5 h and ~3.0 h, providing the cohesive strength ranging between ~4 Pa and ~10 Pa. △ Less

Submitted 6 April, 2019; originally announced April 2019.

Comments: 8 pages, 5 figures, published in the Astrophysical Journal Letters, 874, 1, L10

Realistic On-The-Fly Outcomes of Planetary Collisions: Machine Learning Applied to Simulations of Giant Impacts

Authors: Saverio Cambioni, Erik Asphaug, Alexandre Emsenhuber, Travis S. J. Gabriel, Roberto Furfaro, Stephen R. Schwartz

Abstract: Planet formation simulations are capable of directly integrating the evolution of hundreds to thousands of planetary embryos and planetesimals, as they accrete pairwise to become planets. In principle such investigations allow us to better understand the final configuration and geochemistry of the terrestrial planets, as well as to place our solar system in the context of other exosolar systems. T… ▽ More Planet formation simulations are capable of directly integrating the evolution of hundreds to thousands of planetary embryos and planetesimals, as they accrete pairwise to become planets. In principle such investigations allow us to better understand the final configuration and geochemistry of the terrestrial planets, as well as to place our solar system in the context of other exosolar systems. These simulations, however, classically prescribe collisions to result in perfect mergers, but computational advances have begun to allow for more complex outcomes to be implemented. Here we apply machine learning to a large but sparse database of giant impact studies, streamlining simulations into a classifier of collision outcomes and a regressor of accretion efficiency. The classifier maps a 4-Dimensional parameter space (target mass, projectile-to-target mass ratio, impact velocity, impact angle) into the four major collision types: merger, "graze-and-merge", "hit-and-run", and disruption. The definition of the four regimes and their boundary is fully data-driven; the results do not suffer from any model assumption in the fitting. The classifier maps the structure of the parameter space and provides insights about the outcome regimes. The regressor is a neural network which is trained to closely mimic the functional relationship between the 4-D space of collision parameters, and a real-variable outcome, the mass of the largest remnant. This work is a prototype of a more complete surrogate model, based on extended sets of simulations ("big data"), that will quickly and reliably predict specific collision outcomes for use in realistic N-body dynamical studies of planetary formation. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 18 pages, 7 figures, 3 tables. Accepted for publication on ApJ

Journal ref: ApJ 875:40 (2019)

Guidance, Navigation and Control of Asteroid Mobile Imager and Geologic Observer (AMIGO)

Authors: Greg Wilburn, Himangshu Kalita, Aman Chandra, Stephen Schwartz, Erik Asphaug, Jekan Thangavelautham

Abstract: The science and origins of asteroids is deemed high priority in the Planetary Science Decadal Survey. Major scientific goals for the study of planetesimals are to decipher geological processes in SSSBs not determinable from investigation via in-situ experimentation, and to understand how planetesimals contribute to the formation of planets. Ground based observations are not sufficient to examine S… ▽ More The science and origins of asteroids is deemed high priority in the Planetary Science Decadal Survey. Major scientific goals for the study of planetesimals are to decipher geological processes in SSSBs not determinable from investigation via in-situ experimentation, and to understand how planetesimals contribute to the formation of planets. Ground based observations are not sufficient to examine SSSBs, as they are only able to measure what is on the surface of the body; however, in-situ analysis allows for further, close up investigation as to the surface characteristics and the inner composure of the body. To this end, the Asteroid Mobile Imager and Geologic Observer (AMIGO) an autonomous semi-inflatable robot will operate in a swarm to efficiently characterize the surface of an asteroid. The stowed package is 10x10x10 cm (equivalent to a 1U CubeSat) that deploys an inflatable sphere of ~1m in diameter. Three mobility modes are identified and designed: ballistic hopping, rotation during hops, and up-righting maneuvers. Ballistic hops provide the AMIGO robot the ability to explore a larger portion of the asteroid's surface to sample a larger area than a stationary lander. Rotation during the hop entails attitude control of the robot, utilizing propulsion and reaction wheel actuation. In the event of the robot tipping or not landing up-right, a combination of thrusters and reaction wheels will correct the robot's attitude. The AMIGO propulsion system utilizes sublimate-based micro-electromechanical systems (MEMS) technology as a means of lightweight, low-thrust ballistic hopping and coarse attitude control. Each deployed AMIGO will hop across the surface of the asteroid multiple times. △ Less

Submitted 6 February, 2019; originally announced February 2019.

Comments: 17 pages, 13 figures, AAS Guidance, Navigation and Control Conference 2019

Electron energy partition across interplanetary shocks: I. Methodology and Data Product

Authors: Lynn B. Wilson III, Li-Jen Chen, Shan Wang, Steven J. Schwartz, Drew L. Turner, Michael L. Stevens, Justin C. Kasper, Adnane Osmane, Damiano Caprioli, Stuart D. Bale, Marc P. Pulupa, Chadi S. Salem, Katherine A. Goodrich

Abstract: Analysis of 15314 electron velocity distribution functions (VDFs) within $\pm$2 hours of 52 interplanetary (IP) shocks observed by the \emph{Wind} spacecraft near 1 AU are introduced. The electron VDFs are fit to the sum of three model functions for the cold dense core, hot tenuous halo, and field-aligned beam/strahl component. The best results were found by modeling the core as either a bi-kappa… ▽ More Analysis of 15314 electron velocity distribution functions (VDFs) within $\pm$2 hours of 52 interplanetary (IP) shocks observed by the \emph{Wind} spacecraft near 1 AU are introduced. The electron VDFs are fit to the sum of three model functions for the cold dense core, hot tenuous halo, and field-aligned beam/strahl component. The best results were found by modeling the core as either a bi-kappa or a symmetric (or asymmetric) bi-self-similar velocity distribution function, while both the halo and beam/strahl components were best fit to bi-kappa velocity distribution function. This is the first statistical study to show that the core electron distribution is better fit to a self-similar velocity distribution function than a bi-Maxwellian under all conditions. The self-similar distribution deviation from a Maxwellian is a measure of inelasticity in particle scattering from waves and/or turbulence. The range of values defined by the lower and upper quartiles for the kappa exponents are $κ{\scriptstyle_{ec}}$ $\sim$ 5.40--10.2 for the core, $κ{\scriptstyle_{eh}}$ $\sim$ 3.58--5.34 for the halo, and $κ{\scriptstyle_{eb}}$ $\sim$ 3.40--5.16 for the beam/strahl. The lower-to-upper quartile range of symmetric bi-self-similar core exponents are $s{\scriptstyle_{ec}}$ $\sim$ 2.00--2.04, and asymmetric bi-self-similar core exponents are $p{\scriptstyle_{ec}}$ $\sim$ 2.20--4.00 for the parallel exponent, and $q{\scriptstyle_{ec}}$ $\sim$ 2.00--2.46 for the perpendicular exponent. The nuanced details of the fit procedure and description of resulting data product are also presented. The statistics and detailed analysis of the results are presented in Paper II and Paper III of this three-part study. △ Less

Submitted 2 May, 2019; v1 submitted 4 February, 2019; originally announced February 2019.

Comments: 35 pages, 10 figures, 3 tables, under review in Astrophys. J. Suppl [updated in response to 2nd round of referee comments]

An On-Orbit CubeSat Centrifuge for Asteroid Science and Exploration

Authors: Jekan Thangavelautham, Erik Asphaug, Stephen Schwartz

Abstract: There are thousands of asteroids in near-Earth space and millions expected in the Main Belt. They are diverse in their physical properties and compositions. They are also time capsules of the early Solar System making them valuable for planetary science, and are strategic for resource mining, planetary defense/security and as interplanetary depots. But we lack direct knowledge of the geophysical b… ▽ More There are thousands of asteroids in near-Earth space and millions expected in the Main Belt. They are diverse in their physical properties and compositions. They are also time capsules of the early Solar System making them valuable for planetary science, and are strategic for resource mining, planetary defense/security and as interplanetary depots. But we lack direct knowledge of the geophysical behavior of an asteroid surface under milligravity conditions, and therefore landing on an asteroid and manipulating its surface material remains a daunting challenge. Towards this goal we are putting forth plans for a 12U CubeSat that will be in Low Earth Orbit and that will operate as a spinning centrifuge on-orbit. In this paper, we will present an overview of the systems engineering and instrumentation design on the spacecraft. Parts of this 12U CubeSat will contain a laboratory that will recreate asteroid surface conditions by containing crushed meteorite. The laboratory will spin at 1 to 2 RPM during the primary mission to simulate surface conditions of asteroids 2 km and smaller, followed by an extended mission where the spacecraft will spin at even higher RPM. The result is a bed of realistic regolith, the environment that landers and diggers and maybe astronauts will interact with. The CubeSat is configured with cameras, lasers, actuators and small mechanical instruments to both observe and manipulate the regolith at low simulated gravity conditions. A series of experiments will measure the general behavior, internal friction, adhesion, dilatancy, coefficients of restitution and other parameters that can feed into asteroid surface dynamics simulations. Effective gravity can be varied, and external mechanical forces can be applied. △ Less

Submitted 30 December, 2018; originally announced December 2018.

Comments: 10 pages, 8 figures, IEEE Aerospace Conference 2019

Boulder Stranding in Ejecta Launched by an Impact Generated Seismic Pulse

Authors: Esteban Wright, Alice C. Quillen, Juliana South, Randal C. Nelson, Paul Sanchez, Larkin Martini, Stephen Schwartz, Miki Nakajima, Erik Asphaug

Abstract: We consider how an impact generated seismic pulse affects the surface of an asteroid distant from the impact site. With laboratory experiments on dry polydisperse gravel mixtures, we track the trajectories of particles ejected from the surface by a single strong upward propagating pressure pulse. High speed video images show that ejecta trajectories are independent of particle size, and collisions… ▽ More We consider how an impact generated seismic pulse affects the surface of an asteroid distant from the impact site. With laboratory experiments on dry polydisperse gravel mixtures, we track the trajectories of particles ejected from the surface by a single strong upward propagating pressure pulse. High speed video images show that ejecta trajectories are independent of particle size, and collisions primarily take place upon landing. When they land particles are ballistically sorted, as proposed by Shinbrot et al. (2017), leaving larger particles on the surface and smaller particles more widely dispersed. A single strong pulse can leave previously buried boulders stranded on the surface. Boulder stranding due to an impact excited seismic pulse is an additional mechanism that could leave large boulders present on the surface of rubble asteroids such as 162173 Ryugu, 101955 Bennu and 25143 Itokawa. △ Less

Submitted 13 December, 2019; v1 submitted 4 December, 2018; originally announced December 2018.

Comments: Link to experiment videos: https://www.youtube.com/playlist?list=PLYnpHw29go3rqDIHmaC8BjWKojxr2wSu8

Journal ref: Icarus, Volume 337, 2020, 113424. ISSN 0019-1035

Impact Excitation of a Seismic Pulse and Vibrational Normal Modes on Asteroid Bennu and Associated Slumping of Regolith

Authors: Alice C. Quillen, Yuhui Zhao, YuanYuan Chen, Paul Sanchez, Randal C. Nelson, Stephen R. Schwartz

Abstract: We consider an impact on an asteroid that is energetic enough to cause resurfacing by seismic reverberation and just below the catastrophic disruption threshold, assuming that seismic waves are not rapidly attenuated. In asteroids with diameter less than 1 km we identify a regime where rare energetic impactors can excite seismic waves with frequencies near those of the asteroid's slowest normal mo… ▽ More We consider an impact on an asteroid that is energetic enough to cause resurfacing by seismic reverberation and just below the catastrophic disruption threshold, assuming that seismic waves are not rapidly attenuated. In asteroids with diameter less than 1 km we identify a regime where rare energetic impactors can excite seismic waves with frequencies near those of the asteroid's slowest normal modes. In this regime, the distribution of seismic reverberation is not evenly distributed across the body surface. With mass-spring model elastic simulations, we model impact excitation of seismic waves with a force pulse exerted on the surface and using three different asteroid shape models. The simulations exhibit antipodal focusing and normal mode excitation. If the impulse excited vibrational energy is long lasting, vibrations are highest at impact point, its antipode and at high surface elevations such as an equatorial ridge. A near equatorial impact launches a seismic impulse on a non-spherical body that can be focused on two additional points on an the equatorial ridge. We explore simple flow models for the morphology of vibration induced surface slumping. We find that the initial seismic pulse is unlikely to cause large shape changes. Long lasting seismic reverberation on Bennu caused by a near equatorial impact could have raised the height of its equatorial ridge by a few meters and raised two peaks on it, one near impact site and the other near its antipode. △ Less

Submitted 23 September, 2018; v1 submitted 5 April, 2018; originally announced April 2018.

Comments: submitted to Icarus

Mobility and Science operations On An Asteroid Using a Hopping Small Spacecraft on Stilts

Authors: H. Kalita, S. Schwartz, E. Asphaug, J. Thangavelautham

Abstract: There are thousands of asteroids in near-Earth space and millions in the Main Belt. They are diverse in physical properties and composition and are time capsules of the early solar system. This makes them strategic locations for planetary science, resource mining, planetary defense/security and as interplanetary depots and communication relays. Landing on a small asteroid and manipulating its surf… ▽ More There are thousands of asteroids in near-Earth space and millions in the Main Belt. They are diverse in physical properties and composition and are time capsules of the early solar system. This makes them strategic locations for planetary science, resource mining, planetary defense/security and as interplanetary depots and communication relays. Landing on a small asteroid and manipulating its surface materials remains a major unsolved challenge fraught with high risk. The asteroid surface may contain everything from hard boulders to soft regolith loosely held by cohesion and very low-gravity. Upcoming missions Hayabusa II and OSIRIS-REx will perform touch and go operations to mitigate the risks of landing on an asteroid. This limits the contact time and requires fuel expenditure for hovering. An important unknown is the problem of getting stuck or making a hard impact with the surface. The Spacecraft Penetrator for Increasing Knowledge of NEOs (SPIKE) mission concept will utilize a small-satellite bus that is propelled using a xenon-fueled ion engine and will contain an extendable, low-mass, high-strength boom with a tip containing force-moment sensors. SPIKE will enable contact with the asteroid surface, where it will perform detailed regolith analysis and seismology as well as penetrometry, while keeping the main spacecraft bus at a safe distance. Using one or more long stilts frees the spacecraft from having to hover above the asteroid and thus substantially reduces or eliminates fuel use when doing science operations. This enables much longer missions that include a series of hops to multiple locations on the small-body surface. △ Less

Submitted 29 January, 2018; originally announced January 2018.

Comments: 13 pages, 9 figures, to Appear at AAS GNC 2018/Advances in Astronautical Sciences 2018

Constraints on the perturbed mutual motion in Didymos due to impact-induced deformation of its primary after the DART impact

Authors: Masatoshi Hirabayashi, Stephen R. Schwartz, Yang Yu, Alex B. Davis, Steven R. Chesley, Eugene G. Fahnestock, Patrick Michel, Derek C. Richardson, Shantanu P. Naidu, Daniel J. Scheeres, Andrew F. Cheng, Andrew S. Rivkin, Lance A. M. Benner

Abstract: Binary near-Earth asteroid (65803) Didymos is the target of the proposed NASA Double Asteroid Redirection Test (DART), part of the Asteroid Impact & Deflection Assessment (AIDA) mission concept. In this mission, the DART spacecraft is planned to impact the secondary body of Didymos, perturbing mutual dynamics of the system. The primary body is currently rotating at a spin period close to the spin… ▽ More Binary near-Earth asteroid (65803) Didymos is the target of the proposed NASA Double Asteroid Redirection Test (DART), part of the Asteroid Impact & Deflection Assessment (AIDA) mission concept. In this mission, the DART spacecraft is planned to impact the secondary body of Didymos, perturbing mutual dynamics of the system. The primary body is currently rotating at a spin period close to the spin barrier of asteroids, and materials ejected from the secondary due to the DART impact are likely to reach the primary. These conditions may cause the primary to reshape, due to landslides, or internal deformation, changing the permanent gravity field. Here, we propose that if shape deformation of the primary occurs, the mutual orbit of the system would be perturbed due to a change in the gravity field. We use a numerical simulation technique based on the full two-body problem to investigate the shape effect on the mutual dynamics in Didymos after the DART impact. The results show that under constant volume, shape deformation induces strong perturbation in the mutual motion. We find that the deformation process always causes the orbital period of the system to become shorter. If surface layers with a thickness greater than ~0.4 m on the poles of the primary move down to the equatorial region due to the DART impact, a change in the orbital period of the system and in the spin period of the primary will be detected by ground-based measurement. △ Less

Submitted 6 August, 2017; originally announced August 2017.

Comments: 8 pages, 7 figures, 2 tables, accepted for publication in MNRAS

A Cubesat Centrifuge for Long Duration Milligravity Research

Authors: Erik Asphaug, Jekan Thangavelautham, Andrew Klesh, Aman Chandra, Ravi Nallapu, Laksh Raura, Mercedes Herreras-Martinez, Stephen Schwartz

Abstract: We advocate a low-cost strategy for long-duration research into the 'milligravity' environment of asteroids, comets and small moons, where surface gravity is a vector field typically less than 1/1000 the gravity of Earth. Unlike the microgravity environment of space, there is a directionality that gives rise, over time, to strangely familiar geologic textures and landforms. In addition to advancin… ▽ More We advocate a low-cost strategy for long-duration research into the 'milligravity' environment of asteroids, comets and small moons, where surface gravity is a vector field typically less than 1/1000 the gravity of Earth. Unlike the microgravity environment of space, there is a directionality that gives rise, over time, to strangely familiar geologic textures and landforms. In addition to advancing planetary science, and furthering technologies for hazardous asteroid mitigation and in-situ resource utilization, simplified access to long-duration milligravity offers significant potential for advancing human spaceflight, biomedicine and manufacturing. We show that a commodity 3U ($10\times10\times34$ cm$^3$) cubesat containing a laboratory of loose materials can be spun to 1 rpm = $2π/60$ s$^{-1}$ on its long axis, creating a centrifugal force equivalent to the surface gravity of a kilometer-sized asteroid. We describe the first flight demonstration, where small meteorite fragments will pile up to create a patch of real regolith under realistic asteroid conditions, paving the way for subsequent missions where landing and mobility technology can be flight-proven in the operational environment, in Low-Earth Orbit (LEO). The 3U design can be adapted for use onboard the International Space Station (ISS) to allow for variable gravity experiments under ambient temperature and pressure for a broader range of experiments. △ Less

Submitted 23 May, 2017; v1 submitted 22 May, 2017; originally announced May 2017.

Comments: Accepted to NPJ Microgravity, https://www.nature.com/npjmgrav/

Creep stability of the proposed AIDA mission target 65803 Didymos: I. Discrete cohesionless granular physics model

Authors: Yun Zhang, Derek C. Richardson, Olivier S. Barnouin, Clara Maurel, Patrick Michel, Stephen R. Schwartz, Ronald-Louis Ballouz, Lance A. M. Benner, Shantanu P. Naidu, Junfeng Li

Abstract: As the target of the proposed Asteroid Impact & Deflection Assessment (AIDA) mission, the near-Earth binary asteroid 65803 Didymos represents a special class of binary asteroids, those whose primaries are at risk of rotational disruption. To gain a better understanding of these binary systems and to support the AIDA mission, this paper investigates the creep stability of the Didymos primary by rep… ▽ More As the target of the proposed Asteroid Impact & Deflection Assessment (AIDA) mission, the near-Earth binary asteroid 65803 Didymos represents a special class of binary asteroids, those whose primaries are at risk of rotational disruption. To gain a better understanding of these binary systems and to support the AIDA mission, this paper investigates the creep stability of the Didymos primary by representing it as a cohesionless self-gravitating granular aggregate subject to rotational acceleration. To achieve this goal, a soft-sphere discrete element model (SSDEM) capable of simulating granular systems in quasi-static states is implemented and a quasi-static spin-up procedure is carried out. We devise three critical spin limits for the simulated aggregates to indicate their critical states triggered by reshaping and surface shedding, internal structural deformation, and shear failure, respectively. The failure condition and mode, and shear strength of an aggregate can all be inferred from the three critical spin limits. The effects of arrangement and size distribution of constituent particles, bulk density, spin-up path, and interparticle friction are numerically explored. The results show that the shear strength of a spinning self-gravitating aggregate depends strongly on both its internal configuration and material parameters, while its failure mode and mechanism are mainly affected by its internal configuration. Additionally, this study provides some constraints on the possible physical properties of the Didymos primary based on observational data and proposes a plausible formation mechanism for this binary system. With a bulk density consistent with observational uncertainty and close to the maximum density allowed for the asteroid, the Didymos primary in certain configurations can remain geo-statically stable without including cohesion. △ Less

Submitted 5 March, 2017; originally announced March 2017.

Comments: 66 pages, 24 figures, submitted to Icarus on 25/Aug/2016

Optical Navigation for Interplanetary CubeSats

Authors: Stephen R. Schwartz, Shota Ichikawa, Pranay Gankidi, Nalik Kenia, Graham Dektorand Jekan Thangavelautham

Abstract: CubeSats and small satellites are emerging as low-cost tools for performing science and exploration in deep space. These new classes of satellite exploit the latest advancement in miniaturization of electronics, power systems, and communication technologies to promise reduced launch cost and development cadence. JPL's MarCO CubeSats, part of the Mars Insight mission, will head on an Earth escape t… ▽ More CubeSats and small satellites are emerging as low-cost tools for performing science and exploration in deep space. These new classes of satellite exploit the latest advancement in miniaturization of electronics, power systems, and communication technologies to promise reduced launch cost and development cadence. JPL's MarCO CubeSats, part of the Mars Insight mission, will head on an Earth escape trajectory to Mars in 2018 and serve as communication relays for the Mars Insight Lander during Entry, Descent and Landing. Incremental advancements to the MarCO CubeSats, particularly in propulsion and GNC, could enable these spacecraft to get to another planet or to Near Earth Objects. This can have substantial science return with the right science instrument. We have developed an interplanetary CubeSat concept that includes onboard green monopropellant propulsion system and that can get into a capture orbit around a neighboring planet or chase a small-body. One such candidate is the Martian moon Phobos. Because of the limits of current CubeSat hardware and lack of an accurate ephemeris of Phobos, there will be a 2 to 5 km uncertainty in distance between the spacecraft and Phobos. This presents a major GNC challenge when the CubeSat first attempts to get into visual range of the moon. One solution to this challenge is to develop optical navigation technology that enables the CubeSat to take epicyclic orbits around the most probable location of the target, autonomously search and home-in on the target body. In worst-case scenarios, the technology would narrow down the uncertainty of the small-body location and then use optical flow, a computer vision algorithm to track movement of objects in the field of view. A dimly lit small-body can be detected by the occlusion of one or more surrounding stars. Our studies present preliminary simulations that support the concept. △ Less

Submitted 25 January, 2017; originally announced January 2017.

Comments: in Proceedings of the 40th Annual AAS Guidance, Navigation and Control Conference 2017

Ejecta Cloud from a Kinetic Impact on the Secondary of a Binary Asteroid: I. Mechanical Environment and Dynamic Model

Authors: Yang Yu, Patrick Michel, Stephen R. Schwartz, Shantanu P. Naidu, Lance A. M. Benner

Abstract: An understanding of the post-impact dynamics of ejecta clouds are crucial to the planning of a kinetic impact mission to an asteroid, and also has great implications for the history of planetary formation. The purpose of this article to track the evolution of ejecta produced by AIDA mission, which targets for kinetic impact the secondary of near-Earth binary asteroid 65803 Didymos on 2022, and to… ▽ More An understanding of the post-impact dynamics of ejecta clouds are crucial to the planning of a kinetic impact mission to an asteroid, and also has great implications for the history of planetary formation. The purpose of this article to track the evolution of ejecta produced by AIDA mission, which targets for kinetic impact the secondary of near-Earth binary asteroid 65803 Didymos on 2022, and to feedback essential informations to AIDA's ongoing phase-A study. We present a detailed dynamic model for the simulation of an ejecta cloud from a binary asteroid that synthesizes all relevant forces based on a previous analysis of the mechanical environment. We apply our method to gain insight into the expected response of Didymos to the AIDA impact, including the subsequent evolution of debris and dust. The crater scaling relations from laboratory experiments are employed to approximate the distributions of ejecta mass and launching speed. The size composition of fragments is modeled with a power law fitted from observations of real asteroid surface. A full-scale demonstration is simulated using parameters specified by the mission. We report the results of the simulation, which include the computed spread of the ejecta cloud and the recorded history of ejecta accretion and escape. The violent period of the ejecta evolution is found to be short, and is followed by a stage where the remaining ejecta is gradually cleared. Solar radiation pressure proves to be efficient in cleaning dust-size ejecta, and the simulation results after two weeks shows that large debris on polar orbits (perpendicular to the binary orbital plane) has a survival advantage over smaller ejecta and ejecta that keep to low latitudes. △ Less

Submitted 23 March, 2016; originally announced March 2016.

Comments: 50 pages, 9 figures

MSC Class: 85-02

Small-body deflection techniques using spacecraft: techniques in simulating the fate of ejecta

Authors: Stephen R. Schwartz, Yang Yu, Patrick Michel, Martin Jutzi

Abstract: We define a set of procedures to numerically study the fate of ejecta produced by the impact of an artificial projectile with the aim of deflecting an asteroid. Here we develop a simplified, idealized model of impact conditions that can be adapted to fit the details of specific deflection-test scenarios, such as what is being proposed for the AIDA project. Ongoing studies based upon the methodolog… ▽ More We define a set of procedures to numerically study the fate of ejecta produced by the impact of an artificial projectile with the aim of deflecting an asteroid. Here we develop a simplified, idealized model of impact conditions that can be adapted to fit the details of specific deflection-test scenarios, such as what is being proposed for the AIDA project. Ongoing studies based upon the methodology described here can be used to inform observational strategies and safety conditions for an observing spacecraft. To account for ejecta evolution, the numerical strategies we are employing are varied and include a large N-Body component, a smoothed-particle hydrodynamics (SPH) component, and an application of impactor scaling laws. Simulations that use SPH-derived initial conditions show high-speed ejecta escaping at low angles of inclination, and very slowly moving ejecta lofting off the surface at higher inclination angles, some of which re-impacts the small-body surface. We are currently investigating the realism of this and other models' behaviors. Next steps will include the addition of solar perturbations to the model and applying the protocol developed here directly to specific potential mission concepts such as the proposed AIDA scenario. △ Less

Submitted 21 January, 2016; originally announced January 2016.

Comments: 19 pages, 11 figures, accepted for publication in Advances in Space Research, Special Issue: Asteroids & Space Debris

Dealing with Uncertainties in Asteroid Deflection Demonstration Missions: NEOTwIST

Authors: Siegfried Eggl, Daniel Hestroffer, Juan L. Cano, Javier Martin Avila, Line Drube, Alan W. Harris, Albert Falke, Ulrich Johann, Kilian Engel, Stephen R. Schwartz, Patrick Michel

Abstract: Deflection missions to near-Earth asteroids will encounter non-negligible uncertainties in the physical and orbital parameters of the target object. In order to reliably assess future impact threat mitigation operations such uncertainties have to be quantified and incorporated into the mission design. The implementation of deflection demonstration missions offers the great opportunity to test our… ▽ More Deflection missions to near-Earth asteroids will encounter non-negligible uncertainties in the physical and orbital parameters of the target object. In order to reliably assess future impact threat mitigation operations such uncertainties have to be quantified and incorporated into the mission design. The implementation of deflection demonstration missions offers the great opportunity to test our current understanding of deflection relevant uncertainties and their consequences, e.g., regarding kinetic impacts on asteroid surfaces. In this contribution, we discuss the role of uncertainties in the NEOTwIST asteroid deflection demonstration concept, a low-cost kinetic impactor design elaborated in the framework of the NEOShield project. The aim of NEOTwIST is to change the spin state of a known and well characterized near-Earth object, in this case the asteroid (25143) Itokawa. Fast events such as the production of the impact crater and ejecta are studied via cube-sat chasers and a flyby vehicle. Long term changes, for instance, in the asteroid's spin and orbit, can be assessed using ground based observations. We find that such a mission can indeed provide valuable constraints on mitigation relevant parameters. Furthermore, the here proposed kinetic impact scenarios can be implemented within the next two decades without threatening Earth's safety. △ Less

Submitted 4 February, 2016; v1 submitted 9 January, 2016; originally announced January 2016.

Comments: Accepted for publication in the proceedings of the IAUS 318 - Asteroids: New Observations, New Models, held at the IAU General Assembly in Honolulu, Hawaii, USA 2015

Fire Hose instability driven by alpha particle temperature anisotropy

Authors: Lorenzo Matteini, Petr Hellinger, Steven Schwartz, Simone Landi

Abstract: We investigate properties of a solar wind-like plasma including a secondary alpha particle population exhibiting a parallel temperature anisotropy with respect to the background magnetic field, using linear and quasi-linear predictions and by means of one-dimensional hybrid simulations. We show that anisotropic alpha particles can drive a parallel fire hose instability analogous to that generated… ▽ More We investigate properties of a solar wind-like plasma including a secondary alpha particle population exhibiting a parallel temperature anisotropy with respect to the background magnetic field, using linear and quasi-linear predictions and by means of one-dimensional hybrid simulations. We show that anisotropic alpha particles can drive a parallel fire hose instability analogous to that generated by protons, but that, remarkably, the instability can be triggered also when the parallel plasma beta of alpha particles is below unity. The wave activity generated by the alpha anisotropy affects the evolution of the more abundant protons, leading to their anisotropic heating. When both ion species have sufficient parallel anisotropies both of them can drive the instability, and we observe generation of two distinct peaks in the spectra of the fluctuations, with longer wavelengths associated to alphas and shorter ones to protons. If a non-zero relative drift is present, the unstable modes propagate preferentially in the direction of the drift associated with the unstable species. The generated waves scatter particles and reduce their temperature anisotropy to marginally stable state, and, moreover, they significantly reduce the relative drift between the two ion populations. The coexistence of modes excited by both species leads to saturation of the plasma in distinct regions of the beta/anisotropy parameter space for protons and alpha particles, in good agreement with in situ solar wind observations. Our results confirm that fire hose instabilities are likely at work in the solar wind and limit the anisotropy of different ion species in the plasma. △ Less

Submitted 23 August, 2015; originally announced August 2015.

Comments: 10 pages, 9 figures, Accepted for publication in The Astrophysical Journal

Asteroid Surface Geophysics

Authors: Naomi Murdoch, Paul Sanchez, Stephen R. Schwartz, Hideaki Miyamoto

Abstract: The regolith-covered surfaces of asteroids preserve records of geophysical processes that have occurred both at their surfaces and sometimes also in their interiors. As a result of the unique micro-gravity environment that these bodies posses, a complex and varied geophysics has given birth to fascinating features that we are just now beginning to understand. The processes that formed such feature… ▽ More The regolith-covered surfaces of asteroids preserve records of geophysical processes that have occurred both at their surfaces and sometimes also in their interiors. As a result of the unique micro-gravity environment that these bodies posses, a complex and varied geophysics has given birth to fascinating features that we are just now beginning to understand. The processes that formed such features were first hypothesised through detailed spacecraft observations and have been further studied using theoretical, numerical and experimental methods that often combine several scientific disciplines. These multiple approaches are now merging towards a further understanding of the geophysical states of the surfaces of asteroids. In this chapter we provide a concise summary of what the scientific community has learned so far about the surfaces of these small planetary bodies and the processes that have shaped them. We also discuss the state of the art in terms of experimental techniques and numerical simulations that are currently being used to investigate regolith processes occurring on small-body surfaces and that are contributing to the interpretation of observations and the design of future space missions. △ Less

Submitted 6 March, 2015; originally announced March 2015.

Comments: 28 pages, 8 figures. Chapter to appear in the book ASTEROIDS IV, (University of Arizona Press) Space Science Series, edited by P. Michel, F. DeMeo and W. Bottke

Ion kinetic energy conservation and magnetic field strength constancy in multi-fluid solar wind Alfvénic turbulence

Authors: L. Matteini, T. S. Horbury, F. Pantellini, M. Velli, S. J. Schwartz

Abstract: We investigate properties of the plasma fluid motion in the large amplitude low frequency fluctuations of highly Alfvénic fast solar wind. We show that protons locally conserve total kinetic energy when observed from an effective frame of reference comoving with the fluctuations. For typical properties of the fast wind, this frame can be reasonably identified by alpha particles, which, owing to th… ▽ More We investigate properties of the plasma fluid motion in the large amplitude low frequency fluctuations of highly Alfvénic fast solar wind. We show that protons locally conserve total kinetic energy when observed from an effective frame of reference comoving with the fluctuations. For typical properties of the fast wind, this frame can be reasonably identified by alpha particles, which, owing to their drift with respect to protons at about the Alfvén speed along the magnetic field, do not partake in the fluid low frequency fluctuations. Using their velocity to transform proton velocity into the frame of Alfvénic turbulence, we demonstrate that the resulting plasma motion is characterized by a constant absolute value of the velocity, zero electric fields, and aligned velocity and magnetic field vectors as expected for unidirectional Alfvénic fluctuations in equilibrium. We propose that this constraint, via the correlation between velocity and magnetic field in Alfvénic turbulence, is at the origin of the observed constancy of the magnetic field: while the constant velocity corresponding to constant energy can be only observed in the frame of the fluctuations, the correspondingly constant total magnetic field, invariant for Galilean transformations, remains the observational signature, in the spacecraft frame, of the constant total energy in the Alfvén turbulence frame. △ Less

Submitted 4 January, 2015; originally announced January 2015.

Comments: 6 pages, 6 figures, Accepted for publication in The Astrophysical Journal

Numerical Simulations of Collisional Disruption of Rotating Gravitational Aggregates: Dependence on Material Properties

Authors: Ronald-Louis Ballouz, Derek C Richardson, Patrick Michel, Stephen R. Schwartz, Yang Yu

Abstract: Our knowledge of the strengths of small bodies in the Solar System is limited by our poor understanding of their internal structures, and this, in turn, clouds our understanding of the formation and evolution of these bodies. Observations of the rotational states of asteroids whose diameters are larger than a few hundreds of meters have revealed that they are dominated by gravity and that most are… ▽ More Our knowledge of the strengths of small bodies in the Solar System is limited by our poor understanding of their internal structures, and this, in turn, clouds our understanding of the formation and evolution of these bodies. Observations of the rotational states of asteroids whose diameters are larger than a few hundreds of meters have revealed that they are dominated by gravity and that most are unlikely to be monoliths; however, there is a wide range of plausible internal structures. Numerical and analytical studies of shape and spin limits of gravitational aggregates and their collisional evolution show a strong dependence on shear strength. In order to study this effect, we carry out a systematic exploration of the dependence of collision outcomes on dissipation and friction parameters of the material components making up the bodies. We simulate the catastrophic disruption (leading to the largest remnant retaining 50% of the original mass) of km-size asteroids modeled as gravitational aggregates using pkdgrav, a cosmology N-body code adapted to collisional problems and recently enhanced with a new soft-sphere collision algorithm that includes more realistic contact forces. We find that for a range of three different materials, higher friction and dissipation values increase the catastrophic disruption threshold by about half a magnitude. Furthermore, we find that pre-impact rotation systematically increases mass loss on average, regardless of the target's internal configuration. Our results have important implications for the efficiency of planet formation via planetesimal growth, and also more generally to estimate the impact energy threshold for catastrophic disruption, as this generally has only been evaluated for non-spinning bodies without detailed consideration of material properties. △ Less

Submitted 23 September, 2014; originally announced September 2014.

Comments: 27 pages, 6 figures, Accepted for publication in Planetary and Space Science. arXiv admin note: text overlap with arXiv:1406.5228

Low-speed Impact Simulations into Regolith in Support of Asteroid Sampling Mechanism Design I.: Comparison with 1-g Experiments

Authors: Stephen R. Schwartz, Patrick Michel, Derek C. Richardson, Hajime Yano

Abstract: This study is carried out in the framework of sample-return missions to asteroids that use a low-speed projectile as the primary component of its sampling mechanism (e.g., JAXA's Hayabusa and Hayabusa2 missions). We perform numerical simulations of such impacts into granular materials using different projectile shapes under Earth's gravity. We then compare the amounts of ejected mass obtained in o… ▽ More This study is carried out in the framework of sample-return missions to asteroids that use a low-speed projectile as the primary component of its sampling mechanism (e.g., JAXA's Hayabusa and Hayabusa2 missions). We perform numerical simulations of such impacts into granular materials using different projectile shapes under Earth's gravity. We then compare the amounts of ejected mass obtained in our simulations against what was found in experiments that used similar setups, which allows us to validate our numerical approach. For the targets, we consider 2 different monodisperse grain-diameter sizes: 5 mm and 3 mm. The impact speed of the projectile is 11 m s$^{-1}$ directed downward, perpendicular to the surface of the targets. Using an implementation of the soft-sphere discrete element method (SSDEM) in the $N$-Body gravity tree code PKDGRAV, previously validated in the context of low-speed impacts into sintered glass bead agglomerates, we find a noticeable dependence of the amount of ejected mass on the projectile shape. As found in experiments, in the case of the larger target grain size (5 mm), a conically shaped projectile ejects a greater amount of mass than do projectiles of other shapes, including disks and spheres. We then find that numerically the results are sensitive to the normal coefficient of restitution of the grains, especially for impacts into targets comprised of smaller grains (3 mm). We also find that static friction plays a more important role for impacts into targets comprised of the larger grains. As a preliminary demonstration, one of these considered setups is simulated in a microgravity environment. As expected, a reduction in gravity increases both the amount of ejected mass and the timescale of the impact process. Our methodology is also adaptable to the conditions of sampling mechanisms included in specific mission designs. △ Less

Submitted 8 August, 2014; originally announced August 2014.

Numerical predictions of surface effects during the 2029 close approach of asteroid 99942 Apophis

Authors: Yang Yu, Derek C. Richardson, Patrick Michel, Stephen R. Schwartz, Ronald-Louis Ballouz

Abstract: Asteroid (99942) Apophis' close approach in 2029 will be one of the most significant small-body encounter events in the near future and offers a good opportunity for in situ exploration to determine the asteroid's surface properties and measure any tidal effects that might alter its regolith configuration. Resurfacing mechanics has become a new focus for asteroid researchers due to its important i… ▽ More Asteroid (99942) Apophis' close approach in 2029 will be one of the most significant small-body encounter events in the near future and offers a good opportunity for in situ exploration to determine the asteroid's surface properties and measure any tidal effects that might alter its regolith configuration. Resurfacing mechanics has become a new focus for asteroid researchers due to its important implications for interpreting surface observations, including space weathering effects. This paper provides a prediction for the tidal effects during the 2029 encounter, with an emphasis on whether surface refreshing due to regolith movement will occur. The potential shape modification of the object due to the tidal encounter is first confirmed to be negligibly small with systematic simulations, thus only the external perturbations are taken into account for this work (despite this, seismic shaking induced by shifting blocks might still play a weak role and we will look into this mechanism in future work). A two-stage approach is developed to model the responses of asteroid surface particles (the regolith) based on the soft-sphere implementation of the parallel N-body gravity tree code pkdgrav. A full-body model of Apophis is sent past the Earth on the predicted trajectory to generate the data of all forces acting at a target point on the surface. A sandpile constructed in the local frame is then used to approximate the regolith materials; all the forces the sandpile feels during the encounter are imposed as external perturbations to mimic the regolith's behavior in the full scenario. The local mechanical environment on the asteroid surface is represented in detail, leading to an estimation of the change in global surface environment due to the encounter. Typical patterns of perturbation are presented that depend on the asteroid orientation and sense of rotation at perigee. △ Less

Submitted 1 August, 2014; originally announced August 2014.

The Brazil-nut effect and its application to asteroids

Authors: Soko Matsumura, Derek C. Richardson, Patrick Michel, Stephen R. Schwartz, Ronald-Louis Ballouz

Abstract: Out of the handful of asteroids that have been imaged, some have distributions of blocks that are not easily explained. In this paper, we investigate the possibility that seismic shaking leads to the size sorting of particles in asteroids. In particular, we focus on the so-called Brazil Nut Effect (BNE) that separates large particles from small ones under vibrations. We study the BNE over a wide r… ▽ More Out of the handful of asteroids that have been imaged, some have distributions of blocks that are not easily explained. In this paper, we investigate the possibility that seismic shaking leads to the size sorting of particles in asteroids. In particular, we focus on the so-called Brazil Nut Effect (BNE) that separates large particles from small ones under vibrations. We study the BNE over a wide range of parameters by using the N-body code PKDGRAV, and find that the effect is largely insensitive to the coefficients of restitution, but sensitive to friction constants and oscillation speeds. Agreeing with the previous results, we find that convection drives the BNE, where the intruder rises to the top of the particle bed. For the wide-cylinder case, we also observe a "whale" effect, where the intruder follows the convective current and does not stay at the surface. We show that the non-dimensional critical conditions for the BNE agree well with previous studies. We also show that the BNE is scalable for low-gravity environments and that the rise speed of an intruder is proportional to the square root of the gravitational acceleration. Finally, we apply the critical conditions to observed asteroids, and find that the critical oscillation speeds are comparable to the seismic oscillation speeds that are expected from non-destructive impacts. △ Less

Submitted 10 July, 2014; originally announced July 2014.

Comments: 15 pages, 1 table, 14 figures, accepted for publication in MNRAS

Rotation-Dependent Catastrophic Disruption of Gravitational Aggregates

Authors: Ronald-Louis Ballouz, Derek C. Richardson, Patrick Michel, Stephen R. Schwartz

Abstract: We carry out a systematic exploration of the effect of pre-impact rotation on the outcomes of low-speed collisions between planetesimals modeled as gravitational aggregates. We use pkdgrav, a cosmology code adapted to collisional problems and recently enhanced with a new soft-sphere collision algorithm that includes more realistic contact forces. A rotating body has lower effective surface gravity… ▽ More We carry out a systematic exploration of the effect of pre-impact rotation on the outcomes of low-speed collisions between planetesimals modeled as gravitational aggregates. We use pkdgrav, a cosmology code adapted to collisional problems and recently enhanced with a new soft-sphere collision algorithm that includes more realistic contact forces. A rotating body has lower effective surface gravity than a non-rotating one and therefore might suffer more mass loss as the result of a collision. What is less well understood, however, is whether rotation systematically increases mass loss on average regardless of the impact trajectory. This has important implications for the efficiency of planet formation via planetesimal growth, and also more generally for the determination of the impact energy threshold for catastrophic disruption (leading to the largest remnant retaining 50% of the original mass), as this has generally only been evaluated for non-spinning bodies. We find that for most collision scenarios, rotation lowers the threshold energy for catastrophic dispersal. For head-on collisions, we develop a semi-analytic description of the change in the threshold description as a function of the target's pre-impact rotation rate, and find that these results are consistent with the "universal law" of catastrophic disruption developed by Leinhardt & Stewart. Using this approach, we introduce re-scaled catastrophic disruption variables that take into account the interacting mass fraction of the target and the projectile in order to translate oblique impacts into equivalent head-on collisions. △ Less

Submitted 19 June, 2014; originally announced June 2014.

Comments: 13 pages, 11 figures, 2 tables, Accepted for publication in Astrophysical Journal

Classical-Nova Contribution to the Milky Way's $^{26}$Al Abundance: Exit Channel of the Key $^{25}$Al($p,γ$)$^{26}$Si Resonance

Authors: M. B. Bennett, C. Wrede, K. A. Chipps, J. José, S. N. Liddick, M. Santia, A. Bowe, A. A. Chen, N. Cooper, D. Irvine, E. McNeice, F. Montes, F. Naqvi, R. Ortez, S. D. Pain, J. Pereira, C. Prokop, J. Quaglia, S. J. Quinn, S. B. Schwartz, S. Shanab, A. Simon, A. Spyrou, E. Thiagalingam

Abstract: Classical novae are expected to contribute to the 1809-keV Galactic $γ$-ray emission by producing its precursor $^{26}$Al, but the yield depends on the thermonuclear rate of the unmeasured $^{25}$Al($p,γ$)$^{26}$Si reaction. Using the $β$ decay of $^{26}$P to populate the key $J^π=3^+$ resonance in this reaction, we report the first evidence for the observation of its exit channel via a… ▽ More Classical novae are expected to contribute to the 1809-keV Galactic $γ$-ray emission by producing its precursor $^{26}$Al, but the yield depends on the thermonuclear rate of the unmeasured $^{25}$Al($p,γ$)$^{26}$Si reaction. Using the $β$ decay of $^{26}$P to populate the key $J^π=3^+$ resonance in this reaction, we report the first evidence for the observation of its exit channel via a $1741.6 \pm 0.6 (\textrm{stat}) \pm 0.3 (\textrm{syst})$ keV primary $γ$ ray, where the uncertainties are statistical and systematic, respectively. By combining the measured $γ$-ray energy and intensity with other experimental data on $^{26}$Si, we find the center-of-mass energy and strength of the resonance to be $E_r = 414.9 \pm 0.6(\textrm{stat}) \pm 0.3 (\textrm{syst}) \pm 0.6(\textrm{lit.})$ keV and $ωγ= 23 \pm 6 (\textrm{stat})^{+11}_{-10}(\textrm{lit.})$ meV, respectively, where the last uncertainties are from adopted literature data. We use hydrodynamic nova simulations to model $^{26}$Al production showing that these measurements effectively eliminate the dominant experimental nuclear-physics uncertainty and we estimate that novae may contribute up to 30% of the Galactic $^{26}$Al. △ Less

Submitted 12 December, 2013; originally announced December 2013.

Journal ref: Physical Review Letters 111, 232503 (2013)

Electron-Ion Temperature Equilibration in Collisionless Shocks: the Supernova Remnant-Solar Wind Connection

Authors: Parviz Ghavamian, Steven J. Schwartz, Jeremy Mitchell, Adam Masters, J. Martin Laming

Abstract: Collisionless shocks are loosely defined as shocks where the transition between pre-and post-shock states happens on a length scale much shorter than the collisional mean free path. In the absence of collision to enforce thermal equilibrium post-shock, electrons and ions need not have the same temperatures. While the acceleration of electrons for injection into shock acceleration processes to prod… ▽ More Collisionless shocks are loosely defined as shocks where the transition between pre-and post-shock states happens on a length scale much shorter than the collisional mean free path. In the absence of collision to enforce thermal equilibrium post-shock, electrons and ions need not have the same temperatures. While the acceleration of electrons for injection into shock acceleration processes to produce cosmic rays has received considerable attention, the related problem of the shock heating of quasi-thermal electrons has been relatively neglected. In this paper we review that state of our knowledge of electron heating in astrophysical shocks, mainly associated with supernova remnants (SNRs), shocks in the solar wind associated with the terrestrial and Saturnian bowshocks, and galaxy cluster shocks. The solar wind and SNR samples indicate that the ratio of electron temperature to ion temperature declines with increasing shock speed or Alfvenic Mach number. We discuss the extent to which such behavior can be understood via cosmic ray-generated waves in a shock precursor, which then subsequently damp by heating electrons. Finally, we speculate that a similar mechanism may be at work for both solar wind and SNR shocks. △ Less

Submitted 28 May, 2013; originally announced May 2013.

Comments: 36 pages, 7 figures. accepted for publication in Space Science Reviews

Non-adiabatic electron behaviour due to short-scale electric field structures at collisionless shock waves

Authors: V. See, R. F. Cameron, S. J. Schwartz

Abstract: Under sufficiently high electric field gradients, electron behaviour within exactly perpendicular shocks is unstable to the so-called trajectory instability. We extend previous work paying special attention to shortiscale, high amplitude structures as observed within the electric field profile. Via test particle simulations, we show that such structures can cause the electron distribution to heat… ▽ More Under sufficiently high electric field gradients, electron behaviour within exactly perpendicular shocks is unstable to the so-called trajectory instability. We extend previous work paying special attention to shortiscale, high amplitude structures as observed within the electric field profile. Via test particle simulations, we show that such structures can cause the electron distribution to heat in a manner that violates conservation of the first adiabatic invariant. This is the case even if the overall shock width is larger than the upstream electron gyroradius. The spatial distance over which these structures occur therefore constitutes a new scale length relevant to the shock heating problem. Furthermore, we find that the spatial location of the short-scale structure is important in determining the total effect of non-adiabatic behaviour - a result that has not been previously noted. △ Less

Submitted 17 April, 2013; originally announced April 2013.

Comments: 9 pages, 5 figures

Journal ref: Ann. Geophys., 31, 639-646, 2013