Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Rotational dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

896 records · Page 50

What drives galaxy quenching? A deep connection between galaxy kinematics and quenching in the local Universe

ABSTRACT We develop a 2D inclined rotating disc model, which we apply to the stellar velocity maps of 1862 galaxies taken from the MaNGA survey (SDSS public Data Release 15) . We use a random forest classifier to identify the kinematic parameters that are most connected to galaxy quenching. We find that kinematic parameters that relate predominantly to the disc (such as the mean rotational velocity) and parameters that characterize whether a galaxy is rotation- or dispersion-dominated (such as the ratio of rotational velocity to velocity dispersion) are not fundamentally linked to the quenching of star formation. Instead, we find overwhelmingly that it is the absolute level of velocity dispersion (a property that relates primarily to a galaxy’s bulge/spheroidal component) that is most important for separating star-forming and quenched galaxies. Furthermore, a partial correlation analysis shows that many commonly discussed correlations between galaxy properties and quenching are spurious, and that the fundamental correlation is between quenching and velocity dispersion. In particular, we find that at fixed velocity dispersion, there is only a very weak dependence of quenching on the disc properties, whereby more discy galaxies are slightly more likely to be forming stars. By invoking the tight relationship between black hole mass and velocity dispersion, and noting that black hole mass traces the total energy released by active galactic nuclei (AGNs), we argue that these data support a scenario in which quenching occurs by preventive feedback from AGNs. The kinematic measurements from this work are publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

Use of in Situ Synchrotron Techniques to Probe the Oxidized Surface of Molybdenum Nitride Oxygen Reduction Electrocatalysis

The development of active and stable earth-abundant catalysts for the oxygen reduction reaction (ORR) is needed for widespread, economic development of fuel cell technologies. Designing and optimizing these non-platinum group metals is challenging, however, because they are susceptible to composition and structure changes, including dissolution, oxidation, and corrosion, both in air and under reaction conditions. To identify the active surface, and thus understand the properties that affect activity, the catalyst surface must be characterized in situ. Herein, we utilize a grazing incidence electrochemical cell to investigate in situ composition and morphology changes of a molybdenum nitride (Mo-N) thin film catalyst using grazing incidence x-ray absorption spectroscopy (GI-XAS) and x-ray reflectivity (XRR). In rotating ring disk electrode measurements, we find that the activity, selectivity, stability, and capacitance of the Mo-N catalyst is dependent on the maximum potential to which it has been exposed. Specifically, the overpotential required to reach -2 mA cm-2geo decreases by over 90 mV when the maximum potential is increased from 0.3 to 0.8 V vs RHE (Figure 1). Because the Mo-N oxidizes rapidly in air, ex situ characterization methods including x-ray photoelectron spectroscopy and time-of-flight secondary ion mass spectrometry can provide only limited insight into these in situ catalyst changes. Using in situ GI-XAS measurements at applied potentials between 0.3 and 0.9 V vs RHE, however, we are able to determine that the surface of the film oxidizes and becomes more amorphous when exposed to increasingly higher potentials (Figure 1). Furthermore, the surface remains oxidized on the order of several hours when returned to "ORR relevant potentials" (< 0.6 V vs RHE), indicating that this surface-oxidized nitride is the active surface for ORR. Using in situ XRR measurements, we find that there is no change in surface roughness at potentials below 0.7 V vs RHE, but the film roughens significantly at 0.8 V vs RHE, correlating with ex situ measurements of Mo dissolution at this potential (Figure 1). We therefore conclude that the intrinsic activity of the Mo-N catalyst increases when exposed to potentials up to 0.7 V vs RHE, while above that potential activity enhancements are due to the exposure of more active sites through dissolution. The in situ electrochemical surface-sensitive x-ray characterization as used here is a promising methodology for understanding and leveraging surface dynamics to improve the performance of non-traditional catalysts.

Kreider, Melissa↗

Cluster spin glass correlations and dynamics in Zn 0.5⁢ Mn 0.5⁢ Te

Here, we present a combined magnetometry, muon spin-relaxation (𝜇⁢SR), and neutron-scattering study of the insulating spin glass Zn 0.5 ⁢Mn 0.5 ⁢Te, for which magnetic Mn 2+ and nonmagnetic Zn 2+ ions are randomly distributed on a face-centered cubic lattice. The magnetometry and 𝜇⁢SR results confirm a spin freezing transition around 𝑇 𝑓 ≈ 23 K, with the spin-fluctuation rate decreasing gradually and somewhat inhomogeneously through the sample volume as the temperature decreases toward 𝑇 𝑓 . Characteristic spin-correlation times well above 𝑇 𝑓 are on the order of 10 −10 s, much slower than typically observed in canonical spin glasses but in line with expectations for a cluster spin glass. Using magnetic pair distribution function (mPDF) analysis and reverse Monte Carlo (RMC) modeling of the magnetic diffuse neutron-scattering data, we show that the spin-glass ground state consists of clusters of spins exhibiting short-range-ordered type-III antiferromagnetic correlations with a locally ordered moment of 3.1⁢(1)⁢𝜇 B between nearest-neighbor spins. The type-III correlations decay exponentially as a function of spin separation distance with a correlation length of approximately 5 Å. The diffuse magnetic scattering and corresponding mPDF show no significant changes across 𝑇 𝑓 , indicating that the dynamically fluctuating short-range spin correlations in the paramagnetic state retain the same basic type-III configuration that characterizes the spin-glass state; the only change apparent from the neutron-scattering data is a gradual reduction of the correlation length and locally ordered moment with increasing temperature. Taken together, these results paint a unique and detailed picture of the local magnetic structure and dynamics in Zn 0.5 ⁢Mn 0.5⁢ Te and provide strong evidence that this material is best described as a cluster spin glass. In addition, this work showcases a statistical method for extracting diffuse scattering signals from neutron powder diffraction data, which we developed to facilitate the mPDF and RMC analysis of the neutron data. This method has the potential to be broadly useful for neutron powder diffraction experiments on a variety of materials with short-range atomic or magnetic order.

magnetism↗

Quantum simulations in effective model spaces: Hamiltonian-learning variational quantum eigensolver using digital quantum computers and application to the Lipkin-Meshkov-Glick model

Quantum simulations offer the potential to predict the structure and dynamics of nuclear many-body systems that are beyond the capabilities of classical computing. Generally, preparing the ground state of strongly-interacting many-body systems relevant to nuclear physics is however inefficient, even using ideal quantum computers. In addition, currently available NISQ-era quantum devices possess modest numbers of qubits, limiting the size of quantum many-body systems that can be simulated. In this context, a reformulation of the quantum many-body problems using truncated model spaces and Hamiltonians is desirable to make them more amenable to near-term quantum computers. The importance of symmetries in low-energy theories, including effective field theories (EFTs), lattice quantum chromodynamics (QCD), and effective model spaces for nuclear systems, in particular their interplay with the reduction of active Hilbert spaces, is well known. Lesser known is the fact that the non-commutivity of some symmetries and truncations of the model space can be profitably combined with variational calculations to rearrange the entanglement into localized structures and enable more efficient simulations. Here, the goal of the present study is to explore and utilize the non-commutivity of symmetries and model-space truncations of quantum many-body systems important to nuclear physics, particularly in combination with variational algorithms for quantum simulations and effective Hamiltonian learning. We introduce an iterative hybrid classical-quantum algorithm, Hamiltonian learning variational quantum eigensolver (HL-VQE), that simultaneously optimizes an effective Hamiltonian, thereby rearranging entanglement into the effective model space, and the associated ground-state wavefunction. Quantum simulations, using classical computers and IBM's superconducting-qubit quantum computers, are performed to demonstrate the HL-VQE algorithm, in the context of the Lipkin-Meshkov-Glick (LMG) model of interacting fermions, where the Hamiltonian transformation corresponds to an orbital rotation. We use a mapping where the number of qubits scales with the $\log$ of the size of the effective model space, rather than the particle number. HL-VQE is found to provide an exponential improvement in LMG-model calculations of the ground-state energy and wavefunction, compared to naive truncations without Hamiltonian learning, throughout a significant fraction of the Hilbert space. In the context of EFT, this corresponds to counterterms scaling exponentially with the cut-off as opposed to power law. Implementations on IBM's QExperience quantum computers and simulators for 1- and 2-qubit effective model spaces are shown to provide accurate and precise results, reproducing classical predictions. For a range of parameters defining the LMG model, the HL-VQE algorithm is found to have better scaling of quantum resources requirements than previously explored algorithms. In particular, the HL-VQE scales efficiently over a large fraction of the model space, in contrast to VQE alone. This work constitutes a step in the development of entanglement-driven quantum algorithms for descriptions of nuclear many-body systems. This, in part, leverages the potential of noisy intermediate-scale quantum (NISQ) devices. The exponential scaling of counterterms observed in this study suggests the possibility of more general applicability to other non-perturbative EFTs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Stochastic modeling of x-ray superfluorescence

An approach to modeling the dynamics of x-ray amplified spontaneous emission and superfluorescence, the phenomenon of collective x-ray emission initiated by intense pulses of x-ray free-electron lasers, is developed based on stochastic partial differential equations. The equations are derived from first principles, and the relevant approximations, derivation steps, and extensions specific to stimulated x-ray emission are presented. The resulting equations take the form of three-dimensional generalized Maxwell-Bloch equations augmented with noise terms for both field and atomic variables. The derived noise terms possess specific correlation properties that enable the correct reconstruction of spontaneous emission. Consequently, the developed theoretical formalism is universally suitable for describing all stages of stimulated x-ray emission: spontaneous emission, amplified spontaneous emission, and superfluorescence. Here, we present numerical examples that illustrate various properties of the emitted field, including spatiotemporal coherence and spectral-angular and polarization characteristics. We anticipate that the proposed theoretical framework will establish a robust foundation for interpreting measurements in stimulated x-ray emission spectroscopy, modeling x-ray laser oscillators, and describing other experiments leveraging x-ray superfluorescence.

74 ATOMIC AND MOLECULAR PHYSICS↗

Biomechanical Comparison of Anatomic Restoration of the Ulnar Footprint vs Traditional Ulnar Tunnels in Ulnar Collateral Ligament Reconstruction

Background: Current techniques for ulnar collateral ligament (UCL) reconstruction do not reproduce the anatomic ulnar footprint of the UCL. The purpose of this study was to describe a novel UCL reconstruction technique that utilizes proximal-to-distal ulnar bone tunnels to better re-create the anatomy of the UCL and to compare the biomechanical profile at time zero among this technique, the native UCL, and the traditional docking technique. Hypothesis: The biomechanical profile of the anatomic technique is similar to the native UCL and traditional docking technique. Study Design: Controlled laboratory study. Methods: Ten matched cadaveric elbows were potted with the forearm in neutral rotation. The palmaris longus tendon graft was harvested, and bones were sectioned 14 cm proximal and distal to the elbow joint. Specimen testing included (1) native UCL testing performed at 90° of flexion with 0.5 N·m of valgus moment preload, (2) cyclic loading from 0.5 to 5 N·m of valgus moment for 1000 cycles at 1 Hz, and (3) load to failure at 0.2 mm/s. Elbows then underwent UCL reconstruction with 1 elbow of each pair receiving the classic docking technique using either anatomic (proximal to distal) or traditional (anterior to posterior) tunnel locations. Specimen testing was then repeated as described. Results: There were no differences in maximum load at failure between the anatomic and traditional tunnel location techniques (mean ± SD, 34.90 ± 10.65 vs 37.28 ± 14.26 N·m; P = .644) or when including the native UCL (45.83 ± 17.03 N·m; P = .099). Additionally, there were no differences in valgus angle after 1000 cycles across the anatomic technique (4.58°± 1.47°), traditional technique (4.08°± 1.28°), and native UCL (4.07°± 1.99°). The anatomic group and the native UCL had similar valgus angles at failure (24.13°± 5.86° vs 20.13°± 5.70°; P = .083), while the traditional group had a higher valgus angle at failure when compared with the native UCL (24.88°± 6.18° vs 19.44°± 5.86°; P = .015). Conclusion: In this cadaveric model, UCL reconstruction with the docking technique utilizing proximal-to-distal ulnar tunnels better restored the ulnar footprint while providing valgus stability comparable with reconstruction with the docking technique using traditional anterior-to-posterior ulnar tunnel locations. These results suggest that utilization of the anatomic tunnel location in UCL reconstruction has similar biomechanical properties to the traditional method at the time of initial fixation (ie, not accounting for healing after reconstruction in vivo) while keeping the ulnar tunnels farther from the ulnar nerve. Further studies are warranted to determine if an anatomically based UCL reconstruction results in differing outcomes than traditional reconstruction techniques. Clinical Relevance: Current UCL reconstruction techniques do not accurately re-create the ulnar UCL footprint. The UCL is a dynamic constraint to valgus loads at the elbow, and a more anatomic reconstruction may afford more natural joint kinematics. This more anatomic technique performs similarly to the traditional docking technique at time zero, and the results of this study may offer a starting point for future in vivo studies.

Chang, Edward S.↗

Optically and thermally driven huge lattice orbital and spin angular momenta from spinning fullerenes

Lattice vibration in solids may carry angular momentum. But unlike the intrinsic spin of electrons, the lattice vibration is rarely rotational. To induce angular momentum, one needs to find a material that can accommodate a twisted normal mode, two orthogonal modes, or excitation of magnons. If excitation is too strong, one may exceed the Lindemann limit, so the material melts. Therefore these methods are not ideal. Here, we theoretically propose a route to phonon angular momentum in a molecular crystal C 60 . We find that a single laser pulse is able to inject a significant amount of angular momentum to C 60 , and the momentum transfer is helicity dependent. Changing from right-circularly polarized light to left-circularly polarized light switches the direction of phonon angular momentum. On the ultrafast timescale, the orbital angular momentum change closely resembles the displacive excitation of coherent phonons, with a cosine-function dependence on time, different from the spin counterpart. Atomic displacements, even under strong laser excitation, remain far below the Lindemann criterion. Under thermal excitation, spinning C 60 even at room temperature generates a huge angular momentum close to several hundred h¯. Our finding opens the door to a large group of fullerenes, from C 60 , C 70 to their endohedral derivatives, where angular momentum can be generated through light or temperature. Furthermore, this paves the way to the phononic control electronic spin and harvesting thermal energy through phonon angular momentum.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Theory Guided Fine‐Tune of Strain Effects in Pt Ternary Alloy via Rare Earth Templating: Achieving High Performance PEMFCs Catalysts

The sluggish kinetics and insufficient durability of platinum-based catalysts remain crucial barriers limiting proton-exchange-membrane fuel cells (PEMFCs) deployment. Here, we report a theory-guided synthesis combined with rare-earth templating to realize a previously inaccessible Pt 5 Co-like phase with tailored atomic-scale strain. Guided by density functional theory (DFT) calculations, we identified that a Pt 5 Co-like sublayer can induce a unique mild compressive strain (−1.24%) to the Pt(111) shell and an optimal *OH binding energy shift (ΔE ≈ 0.11 eV). This shift positions the alloy catalyst near the apex of the oxygen reduction reaction activity volcano. This prediction guided the synthesis of ternary alloy Pt 5 (Ce)Co@Pt multilayer nanoparticles, featuring a Ce-stabilized core, a Pt 5 Co-like sublayer, and a Pt-rich shell. This catalyst demonstrates both exceptionally high activity and durability, achieving a mass activity of 2.6 A∙mg Pt −1 in rotating disk electrode testing. In fuel cell membrane electrode assembly tests, Pt 5 (Ce)Co@Pt achieves a current density of 1.9 A∙cm −2 at 0.7 V under heavy-duty vehicle conditions. Remarkably, it maintains 1.2 A∙cm −2 after 1 80 000 AST cycles, doubling the U.S. DOE 2025 target. This work demonstrates a rational design strategy that DFT-guided strain engineering integrates with rare-earth templating to advance Pt-based catalysts for fuel cell applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

The Milky Way’s stellar streams and globular clusters do not align in a Vast Polar Structure

ABSTRACT There is increasing evidence that a substantial fraction of Milky Way satellite galaxies align in a rotationally supported plane of satellites, a rare configuration in cosmological simulations of galaxy formation. It has been suggested that other Milky Way substructures (namely young halo globular clusters and stellar/gaseous streams) similarly tend to align with this plane, accordingly dubbed the Vast Polar Structure (VPOS). Using systemic proper motions inferred from Gaia data, we find that globular cluster orbital poles are not clustered in the VPOS direction, though the population with the highest VPOS membership fraction is the young halo clusters (∼30 per cent). We additionally provide a current census of stellar streams, including new streams discovered using the Dark Energy Survey and Gaia data sets, and find that stellar stream normals are also not clustered in the direction of the VPOS normal. We also find that, based on orbit modelling, there is a likely association between NGC 3201 and the Gjöll stellar stream and that, based on its orbital pole, NGC 4147 is likely not a Sagittarius globular cluster. That the Milky Way’s accreted globular clusters and streams do not align in the same planar configuration as its satellites suggests that the plane of satellites is either a particularly stable orbital configuration or a population of recently accreted satellites. Neither of these explanations is particularly likely in light of other recent studies, leaving the plane of satellites problem as one of the more consequential open problems in galaxy formation and cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

Experimental investigation of turbulent wake flows in a helically wrapped rod bundle in presence of localized blockages

In nuclear sodium fast reactors, bundles of rods are tightly packed into a triangular lattice, enclosed in a hexagonal duct, and each pin is spirally wrapped with a thin wire. Flow blockages can potentially impact the local flow characteristics and heat transfer mechanisms in the bundle due to its small subchannel size. The effects of the blockage on the flow structures and heat transfer mechanisms are important aspects that require an accurate investigation. In this study, the flow-field characteristics in the vicinity of a blockage located in the exterior subchannel of rod bundles with helically wrapped wires were experimentally investigated. The velocity fields in the exterior subchannel were acquired by applying matched-index-of-refraction and time-resolved particle image velocimetry (TR-PIV) techniques for Reynolds numbers of Re 1 = 4000 and Re 2 = 17 000, i.e., equivalent to Rew 1 = 19 600 and Rew 2 = 83 200, respectively, based on the blockage width. The results from the TR-PIV measurements revealed an arch-shaped vortex with a large flow recirculation and a pair of counter-rotating vortices in the wake region downstream of the blockage, which is commonly observed in the wake flow of bluff bodies. The relative lateral distance and angle between the two vortices decreased when the Reynolds numbers increased. Profiles of maximum turbulence intensity along the shear layers illustrated the transition process including the growth, peak, and decay along the flow direction. From the spectral analysis of the turbulent velocities extracted at points along the shear layer, the Strouhal numbers ( St ) representing the vortex shedding frequency were found to be St = 0.25 and St = 0.56 for the left and right shear layers, respectively. Characteristics of shear layers generated by the blockage in the exterior subchannel were investigated via the two-point cross correlation of fluctuating velocities. The patiotemporal cross correlations of turbulent velocities, computed at points in the region where the left shear layer exhibited rolling effects and vortex breakdowns, were considerably wider and longer. The convection velocity U c was estimated to be ~0.82 U m to 0.93 U m . Proper orthogonal decomposition (POD) analysis was applied to the instantaneous velocity fields to extract the statistically dominant flow structures. It was found that POD modes 2–3 and 4–5 formed the pair modes when the corresponding POD temporal coefficients depicted sinusoidal shapes and exhibited nearly circular orbits in the phase space. Spectral analysis of the POD temporal coefficients confirmed the vortex shedding frequencies detected in the analysis of turbulent velocities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Hydrodynamics and Survivability during Post-main-sequence Planetary Engulfment

The engulfment of substellar bodies (SBs), such as brown dwarfs and planets, by giant stars is a possible explanation for rapidly rotating giants, lithium-rich giants, and the presence of SBs in close orbits around subdwarfs and white dwarfs. We perform three-dimensional hydrodynamical simulations of the flow in the vicinity of an engulfed SB. We model the SB as a rigid body with a reflective surface because it cannot accrete. This reflective boundary changes the flow morphology to resemble that of engulfed compact objects with outflows. We measure the drag coefficients for the ram-pressure and gravitational drag forces acting on the SB, and use them to integrate its trajectory inside the star. We find that engulfment can increase the luminosity of a 1 M ⊙ star by up to a few orders of magnitude. The time for the star to return to its original luminosity is up to a few thousand years when the star has evolved to ≈10 R ⊙ and up to a few decades at the tip of the red giant branch (RGB). No SBs can eject the envelope of a 1 M ⊙ star before it evolves to ≈10 R ⊙ if the orbit of the SB is the only energy source contributing to the ejection. In contrast, SBs as small as ≈10 MJup can eject the envelope at the tip of the RGB. The numerical framework we introduce here can be used to study planetary engulfment in a simplified setting that captures the physics of the flow at the scale of the SB.

79 ASTRONOMY AND ASTROPHYSICS↗

Towards modelling AR Sco: calibration – reproducing high-energy pulsar emission and testing convergence to Aristotelian electrodynamics

In recent years, kinetic simulations have been crucial to further our understanding of pulsar electrodynamics. Yet, due to the large-scale separation between the gyro-period and the stellar rotation period, resolving the particle gyration has been computationally unfeasible for realistic pulsar parameters. The main aim of this work is comparing our gyro-phase-resolved model with a gyro-centric pulsar model, where our model solves the general equations of motion with included radiation reaction using a higher order numerical solver with adaptive time-steps. Specifically, we aim to (i) reproduce a pulsar’s high-energy emission maps, namely one with 10 per cent of the surface B-field strength of Vela, and the spectra produced by an independent gyro-centric pulsar emission model; and (ii) test convergence of these results to the radiation-reaction limit of Aristotelian electrodynamics. (iii) Additionally, we identify the effect that a large $E_{\parallel }$-field has on the trajectories and radiation calculations. We find that we can reproduce the curvature radiation emission maps and spectra well, using 10 per cent field strengths of the Vela pulsar and injecting our particles at a higher altitude in the magnetosphere. Using sufficiently large $E_{\parallel }$-fields, our numeric results converge to the analytic radiation-reaction limit trajectories. Additionally, we illustrate the importance of accounting for the $\mathbf {E}\times \mathbf {B}$-drift in the particle trajectories and radiation calculations, validating the Harding and collaborators’ model approach. Lastly, we found that our model deals very well with the high-radiation-reaction and high-field regimes present in pulsars.

79 ASTRONOMY AND ASTROPHYSICS↗

Benchmark for two-dimensional large scale coherent structures in partially magnetized E × B plasmas—community collaboration & lessons learned

Low-temperature plasmas (LTPs) are essential to both fundamental scientific research and critical industrial applications. As in many areas of science, numerical simulations have become a vital tool for uncovering new physical phenomena and guiding technological development. Code benchmarking remains crucial for verifying implementations and evaluating performance. This work continues the Landmark benchmark initiative, a series specifically designed to support the verification of LTP codes. In this study, seventeen simulation codes from a collaborative community of nineteen international institutions modeled a partially magnetized E × B Penning discharge. The emergence of large scale coherent structures, or rotating plasma spokes, endows this configuration with an enormous range of time scales, making it particularly challenging to simulate. The codes showed excellent agreement on the rotation frequency of the spoke as well as key plasma properties, including time-averaged ion density, plasma potential, and electron temperature profiles. Achieving this level of agreement came with challenges, and we share lessons learned on how to conduct future benchmarking campaigns. Comparing code implementations, computational hardware, and simulation runtimes also revealed interesting trends, which are summarized with the aim of guiding future plasma simulation software development.

benchmarking↗

High-speed Boulders and the Debris Field in DART Ejecta

On 2022 September 26 the Double Asteroid Redirection Test (DART) spacecraft collided with Dimorphos, the moon of the near-Earth asteroid 65803 Didymos, in a full-scale demonstration of a kinetic impactor concept. The companion Light Italian Cubesat for Imaging of Asteroids (LICIACube) spacecraft documented the aftermath, capturing images of the expansion and evolution of the ejecta from 29 to 243 s after the impact. We present results from our analyses of these observations, including an improved reduction of the data and new absolute calibration, an updated LICIACube trajectory, and a detailed description of the events and phenomena that were recorded throughout the flyby. One notable aspect of the ejecta was the existence of clusters of boulders, up to 3.6 m in radius, that were ejected at speeds of up to 52 m s −1 . Our analysis of the spatial distribution of 104 of these boulders suggests that they are likely the remnants of larger boulders shattered by the DART spacecraft in the first stages of the impact. The amount of momentum contained in these boulders is more than 3 times that of the DART spacecraft, and it is directed primarily to the south, almost perpendicular to the DART trajectory. Recoil of Dimorphos from the ejection of these boulders has the potential to change its orbital plane by up to a degree and to impart a non-principal-axis component to its rotation state. Damping timescales for these phenomena are such that the Hera spacecraft, arriving at the system in 2026, should be able to measure these effects.

asteroid dynamics↗