Main

The DART spacecraft was launched on 24 November 2021. The spacecraft carried a narrow-angle imager called the Didymos Reconnaissance and Asteroid Camera for Optical navigation (DRACO), which was used for optical navigation, terminal guidance and asteroid characterization8. DRACO detected Didymos—the primary asteroid in the binary system—in summed optical navigation images 61 days before impact. On 27 August 2022, 30 days before impact, DRACO began taking optical navigation images of Didymos every 5 hours, which were processed by the ground optical navigation team9.

On 26 September 2022 at 19:09:24 utc, 4 h and 5 min before impact, the spacecraft’s autonomous Small-body Maneuvering Autonomous Real Time Navigation (SMART Nav) system10 took control of spacecraft navigation (Fig. 1a). SMART Nav processed DRACO images onboard11 to identify Didymos, and once resolved, Dimorphos. To achieve impact with Dimorphos, the spacecraft needed to distinguish between the two asteroids and hit the smaller, dimmer one. Owing to the dynamics of the binary system and the spatial resolution of DRACO, it was known that the spacecraft’s ultimate target—Dimorphos—would be hidden from DRACO’s view during most of the autonomous phase. By design, SMART Nav manoeuvred the spacecraft towards Didymos until Dimorphos became reliably detected10. SMART Nav first detected Dimorphos 73 min before impact, and at 50 min before impact, SMART Nav began manoeuvring towards Dimorphos (Fig. 1). As planned, SMART Nav manoeuvring ended at 23:11:52 utc, 2.5 min before impact, to give the spacecraft time to settle to minimize jitter and smear in the final images. The spacecraft impacted Dimorphos at 23:14:24.183 ± 0.004 utc (Methods). From the start of autonomous navigation until impact, the spacecraft continuously streamed images to the ground, which were immediately broadcast to the public. The final full image was acquired 1.818 s before impact and has a pixel scale of 5.5 cm. The final image received on the ground was a partial image acquired 0.855 s before impact with a pixel scale of 2.6 cm.

Fig. 1: Milestones leading to the impact with Dimorphos from the time SMART Nav began targeting until the end of SMART Nav manoeuvring.
figure 1

a–d, Each column corresponds to a milestone: Didymos targeted (a), Dimorphos detected (b), Dimorphos targeted (c) and end of manoeuvring (d). Each row shows, from top to bottom, the raw DRACO image at the time of that milestone with circles indicating the two asteroids detected by onboard processing and identified by SMART Nav (yellow dashed circles, Didymos; green solid circles, Dimorphos), a zoom-in of Didymos and a zoom-in of Dimorphos. The SMART Nav system used information in DRACO images to successfully impact Dimorphos. In all images, the north pole of Dimorphos (+Z) is towards the bottom left. Images from left to right: dart_0401915351_36903_01_raw.fits, dart_0401925635_06853_01_raw.fits, dart_0401927052_23729_01_raw.fits and dart_0401929899_33346_01_raw.fits.

Full size image

Little was known about the shape or surface of Dimorphos until DRACO obtained high-resolution images. Ground-based radar observations12 led to a diameter estimate of 150 ± 30 m. Analysis of telescopic photometric observations yielded a comparable diameter for Dimorphos, 171 ± 11 m (refs. 13,14). Although DRACO imaged only a portion of Dimorphos, and illumination was limited at 60° solar phase, the images were used to construct an asteroid shape model (Methods). The shape model revealed Dimorphos to be an oblate spheroid with a volume-equivalent diameter of 151 ± 5 m (Table 1, Fig. 2 and Extended Data Fig. 1). The shape of Dimorphos is unusual relative to other near-Earth asteroids visited by spacecraft15,16,17,18,19 and differs from other binary asteroid secondaries observed so far that have measured elongations20,21,22,23. However, oblate secondaries show little or no measurable light-curve amplitude, which biases the observational sample towards elongated secondaries. A size estimate for Didymos from DRACO images (Table 1 and Methods) combined with previous telescopic observations4 enables calculation of a more accurate visible (0.55 μm) geometric albedo for the system of 0.15 ± 0.02. This value is on the low side, but within 1σ, of the mean geometric albedo for S-type asteroids24.

Table 1 Properties of the DART impact, Dimorphos and Didymos
Full size table
Fig. 2: The asteroid Dimorphos and the DART impact site as seen in calibrated DRACO images.
figure 2

a, Dimorphos with an appropriately scaled and correctly oriented outline of the DART spacecraft centred on the impact site. Note the size of the spacecraft relative to the asteroid. The spacecraft bus was approximately 1.2 × 1.3 × 1.3 m, from which other structures extended, resulting in dimensions of approximately 1.8 × 1.9 × 2.6 m. The spacecraft also had two large solar arrays that were each 8.5 m long. b, A closer view of the DART impact site showing the outline of the spacecraft bus and solar arrays over the DRACO image. Note the positions of the two solar arrays relative to two large boulders, labelled 1 (6.5 m long) and 2 (6.1 m long). This subframe is from an image taken 2.781 s before impact. c, The spacecraft bus hit between boulders 1 and 2, whereas the solar arrays interacted with these boulders. This subframe is from an image taken 1.818 s before impact. The arrow in the bottom right of a indicates the direction of the Dimorphos +Z (north) axis. The solid white box in a shows the location of the image in b. The dashed white box in b shows the location of the image in c. Panels b and c show subimages of the full frame. Image names: dart_0401930039_14119_01_iof.fits (a), dart_0401930048_45552_01_iof.fits (b) and dart_0401930049_43695_01_iof.fits (c).

Full size image

The spacecraft trajectory and pointing were reconstructed to locate the impact site (Methods and Fig. 2). The spacecraft impacted Dimorphos at 8.84 ± 0.45° S, 264.30 ± 0.47° E, within 25 m of the centre of the figure, which is very near the scenario for maximizing momentum transfer with an impact through the centre of the figure25. The 1σ uncertainty in the impact site location is ±68 cm (Methods), which is smaller than the size of the spacecraft bus (Fig. 3). The impact angle was 73 ± 7° from local horizontal (Methods and Extended Data Fig. 2). The impact site was near two large boulders, labelled boulder 1 (6.5-m long and at its highest point standing about 2.2 m above the surrounding terrain; Extended Data Fig. 3) and boulder 2 (6.1-m long and at its highest point standing about 1.6 m above the surrounding terrain; Extended Data Fig. 3) in Fig. 2b,c. The spacecraft approached the asteroid with its solar arrays canted slightly towards the surface (Fig. 3). The leading edge of the +Y solar array contacted the surface of boulder 1, and this solar array directly hit boulder 1 (Fig. 3). Almost immediately thereafter, the −Y solar array grazed boulder 2, with the leading edge of the −Y array contacting the surface near the base of boulder 2 in downrange direction (Fig. 3). Finally, the spacecraft bus hit between boulders 1 and 2 (Fig. 3). Although the solar arrays contacted Dimorphos just before the spacecraft bus, the bulk of the spacecraft’s energy was transferred by the bus, which accounted for about 88% of the spacecraft mass at the time of impact.

Fig. 3: Relationship between the spacecraft and topography at the DART impact site.
figure 3

a–c, The position of the spacecraft immediately before the impact of the spacecraft bus from different perspectives to visualize the three-dimensional interactions between the spacecraft and surface. a, Dimorphos north is towards the top of the panel. b, Dimorphos north is to the right. c, Dimorphos north is roughly into the page. In all panels, the −Y solar array points to Dimorphos north. Length scales vary in these perspective views; the scale bars shown are approximate. Boulders 1 and 2 correspond to boulders 1 and 2 in Fig. 2. The caption to Fig. 2 gives the spacecraft dimensions.

Full size image

DART images of Dimorphos revealed a boulder-strewn surface (Fig. 2a, and Extended Data Figs. 4 and 5) resembling other small near-Earth asteroids, such as the S-type (25143) Itokawa17, and carbonaceous asteroids (101955) Bennu26 and (162173) Ryugu16, suggesting a rubble-pile structure for Dimorphos. The boulder-rich nature of the surface is apparent in images as coarse as 2–3-m pixel scale (Extended Data Fig. 4). No unambiguous impact crater candidates are observed, which indicates a young surface, although craters can be difficult to identify on boulder-covered terrains27,28,29,30. The appearance of Dimorphos contrasts with impressions from lower-resolution images of the Didymos surface, where regional roughness variations are apparent (Extended Data Fig. 4).

The final complete DRACO image captured an approximately 880-m2 area surrounding the impact site at a 5.5-cm pixel scale. The impact region (Fig. 2b) exhibits blocky terrain resembling the rest of the observed hemisphere (Fig. 2a, and Extended Data Figs. 4 and 5). There is evidence for variation within boulders, ‘rocks on rocks’ similar to observations on Bennu31, and partially buried boulders (Extended Data Fig. 5). The longest axes of boulders counted in the final complete image are 0.16 m to 6.5 m in length. The impact region has fewer boulders in the 0.2–0.5 m size range than expected if the cumulative distribution followed a single power law, even though the pixel scale of the image is sufficient for their identification (Extended Data Fig. 6). There is no evidence for expansive smooth deposits (grain size smaller than the image pixel scale) such as those seen on Itokawa17. The blocky nature of the impact site probably influenced crater formation, ejecta and momentum enhancement, as seen in impact experiments28,32,33,34, numerical simulations25,35 and the Small Carry-on Impactor experiment on Hayabusa25.

DART did not measure the mass of Dimorphos. Instead, the mass of Dimorphos is estimated using the orbital properties of the binary system, the total volume of the system and an assumption that Didymos and Dimorphos have equal bulk densities (Methods). This assumption cannot be rigorously tested from DART data, but this approach leads to a bulk density of Dimorphos of 2,400 kg m−3 (Table 1) with difficult-to-quantify uncertainties.

On the basis of analyses of reflectance spectra, the best meteoritic analogues for Didymos are L and LL chondrites36,37. L and LL chondrites have grain densities38 of 3,580 ± 10 kg m−3 and 3,520 ± 10 kg m−3, respectively. If one assumes Dimorphos has the same composition as Didymos and that these meteorite values represent the grain density of Dimorphos, then a bulk density of about 2,400 kg m−3 implies a Dimorphos bulk porosity of the order of 30% (with a difficult-to-quantify uncertainty; Methods). This level of bulk porosity is not inconsistent with a rubble-pile structure for Dimorphos, a structure suggested by the boulder-rich character of Dimorphos’s surface. This bulk porosity probably exists as a combination of macroporosity between pieces of rubble and microporosity within individual pieces of rubble. L and LL chondrite samples have porosities of 8.0 ± 0.3% and 9.5 ± 0.6%, respectively38, which would imply that macroporosity is substantial on Dimorphos. Estimates of the density and porosity of Dimorphos will  improve when the European Space Agency’s Hera mission arrives at the Didymos system in early 202739.

DART’s successful autonomous targeting of a small asteroid with limited prior knowledge is a key accomplishment on the path to advancing kinetic impactor technology to an operational capability. The impact of DART indicates that a precursor reconnaissance mission is not a prerequisite for intercepting a subkilometre asteroid, although the characterization done by a precursor mission would provide valuable information for optimizing, planning and predicting the outcome with greater certainty. Kinetic impactor technology for asteroid deflection requires having sufficient warning time—at least several years but preferably decades—to prevent an asteroid impact with Earth1,2,3. Nevertheless, this successful step to demonstrate the viability of kinetic impactor technology for planetary defence builds optimism about humanity’s capacity to protect Earth from an asteroid threat.

Methods

Determining the time of impact

The time of impact was computed from spacecraft clock timestamps inserted into each downlinked telemetry frame by the spacecraft’s radio. A new telemetry frame was produced every 2.9853 ms with a precision of 20 μs. Assuming that the impact occurred in the middle of the frame, the final received frame placed the time of impact to within ±half a frame period: 401,930,051.59326 ± 0.00149 s spacecraft clock time. This spacecraft clock time of impact was converted to utc. This conversion increases the uncertainty slightly: 2022-26-9 23:14:24.183 ± 0.004 s utc. This time is the utc of the impact at Dimorphos, not the time on Earth when the last telemetry frame was received. As each DRACO image is timestamped with a spacecraft clock value, the time to impact for each image can be computed with similar accuracy.

Shape modelling of Dimorphos

We built the shape model using stereophotoclinometry (SPC), a technique that has been widely used to model the shapes of small bodies40,41,42. Before the impact of DART, the shape modelling team conducted tests to understand the strengths and limitations of the SPC approach given the data expected from DART43. DART images pose a challenge for any image-based shape modelling technique owing to the single viewing geometry and lighting condition. SPC provides robust results despite these challenges43.

Extended Data Fig. 7 illustrates the steps used to build the Dimorphos shape model. First, we used DRACO images to determine the dimensions of an initial triaxial ellipsoid. The ellipsoid was constrained by the location of the sunlit limb, the location of the limb lit by light scattered off Didymos and the position of the terminator. Together, the sunlit and Didymos-lit limbs revealed the complete extents of the X and Z axes of Dimorphos (Extended Data Fig. 8). The position of the terminator constrained the Y-axis extent. The last DRACO image to contain all of Dimorphos in the field of view was registered to the ellipsoid. Other DRACO images were registered to that image or, in the case of the highest-resolution images of the impact site, to the immediately preceding image. To correctly determine the scale of the shape model, we used the known time of impact and spacecraft speed to set the spacecraft range for each image.

Next, we pursued parallel paths (Extended Data Fig. 7). The first path used standard SPC processes to construct small digital terrain models (called ‘maplets’) using an SPC program called lithos40,42,43. Maplets were tiled all across the surface of the asteroid in the areas seen by DRACO. Maplets went through several iterations to compute the maplet topography40,42,43. After several iterations, the maplet ensemble was averaged to create a global digital terrain model (DTM). Areas without maplet coverage remained unchanged. This cycle was repeated with maplets of increasingly finer ground-sample distances (GSDs) in the areas covered by higher-resolution images surrounding the impact site. The finest-scale maplets had a GSD of 3 cm, comparable to the pixel scale of the final, partial DRACO image. The second path focused on matching the sunlit limb. We built maplets along the sunlit limb that were conditioned by limb points identified in the tilt-to-height integration in an SPC iteration. These maplets were made at only two GSDs owing to the coarser resolution of images that contained limb. Once the topography in the standard maplets and limb-only maplets stabilized, we united them via several SPC iterations and proceeded to build the global DTM.

The areas covered by maplets in the Dimorphos global DTM are shown in Extended Data Fig. 1 as shaded grey regions. The blue and magenta points show the locations of the sunlit and Didymos-lit limbs. The unseen side of the asteroid is roughly centred on the 90° E view. The shape model has the volume and extents reported in Table 1. The global DTM has a typical GSD of 26 cm, but that resolution is only meaningful in the areas covered by maplets.

We also used the standard maplets to construct a DTM of the impact site. The impact site is covered by the highest-resolution images, so the impact site can support finer-scale topography than the global model. The impact site DTM has a GSD of 5 cm and is shown in Fig. 3.

Owing to the short period of time for which we have resolved images of Dimorphos, the rotation pole, prime meridian and rotation rate of Dimorphos could not be updated using SPC. Instead, we used values derived by ground-based observers. In the equatorial J2000 frame, the pole values used for Dimorphos are14: BODY120065803_POLE_RA = (60.936309840897856, 0, 0) and BODY120065803_POLE_DEC = (−71.674565992873852, 0, 0).

The coordinate system of the global DTM of Dimorphos has the prime meridian pointed towards Didymos, consistent with International Astronomical Union convention: BODY120065803_PM = (64.914870949788195, 724.723943017441570, 1.0840372309 × 10−6).

A planetary constants kernel (PCK) with these values is available at https://ssd.jpl.nasa.gov/ftp/eph/small_bodies/dart/dimorphos/archive/ called ‘dimorphos_s501-preimp.tpc’. Documentation for PCK files, including units for the terms listed above, is available at https://naif.jpl.nasa.gov/pub/naif/toolkit_docs/FORTRAN/req/pck.html#Text%20PCK%20Kernel%20Variable%20Names. The north pole of Dimorphos is in the +Z direction of the shape model but oriented towards ecliptic south.

Shape model uncertainties

We quantified uncertainties in the global DTM of Dimorphos using techniques employed to estimate the uncertainties associated with the global DTM of Bennu44. These techniques were also used to assess uncertainties in shape models developed during pre-impact tests done for DART43. The analyses involved comparing DRACO images against the global DTM when rendered with the same lighting and viewing geometry as the image. We used three analyses described in the above references. The first method, referred to as the limb and terminator method, used image thresholding to identify the limb and terminator in the DRACO image and rendered shape model. The thresholded images were subtracted to reveal mismatches between the image and rendered shape model and any size bias in the model (Extended Data Fig. 9a–e). The second and third methods relied on analysis of corresponding surface features, or keypoints, in images and the rendered shape model, to understand errors in the overall size of the shape model (Extended Data Fig. 9f–i). In the second method, referred to as keypoint matching, the rendered image was rotated, translated and scaled to minimize differences in keypoint locations on the DRACO image and rendered shape model. The third method, referred to as keypoint distance, was based on the measured distances between all keypoints in the DRACO image and rendered shape model.

The limb and terminator assessments indicated that the model is 75 cm too small in X, Y and Z, with a limb uncertainty of 1.3–2 m. Most of the mismatch between DRACO images and the rendered shape model occurred along the limb and terminator. We expected larger uncertainties in these areas owing to the limited coverage and single lighting condition in DRACO images. The keypoints indicated that the Dimorphos shape model is 6 ± 29 cm too large and that features in the model have point-separation errors of −11 ± 20 cm. The keypoint assessments indicated that the model performed well in areas where DRACO resolved detailed surface features, which is the area of the model that is most relevant to understanding the topography of the impact site. Comparisons between shadow lengths in DRACO images and shadow lengths in shape model renderings indicated that some boulder heights are too small. This mismatch was expected based on shape modelling tests done for DART and is a consequence of constructing the shape model from a single lighting condition and viewing geometry43. Shadow lengths suggest that the height of boulder 1 in the impact site DTM is about 10% too small, but the height of boulder 2 in the impact site DTM is correct. On the basis of all of these analyses, we assumed uncertainties of 2 m in the extents of Dimorphos in the X and Z directions.

Because DRACO saw a complete outline of the asteroid (Extended Data Fig. 8) and the spacecraft approach geometry was such that this outline was primarily in the X–Z plane, the shape model dimension uncertainties in X and Z are less than 1% (Table 1). Hence, the Y-axis extent is the largest source of uncertainty in the volume. Future work will refine the errors in the Y axis, but for the time being we assume that the uncertainty in Y is two times the uncertainties in X and Z (that is, 4 m).

Measuring the volume of Dimorphos and realistically estimating uncertainties on that quantity is of particular interest because it directly affects our understanding of the mass of the asteroid. In the shape modelling tests43, we took scaled shape models of an ellipsoid, Itokawa, and Bennu, rendered a set of simulated DRACO images from the ‘truth’ shape models, and then used SPC to build shape models from simulated images43. We compared the volumes of the models built using SPC against the volumes of truth models43. The volume errors in those tests ranged from −2% to +23% (ref. 43). The test with the lowest volume error was for an ellipsoidal asteroid. The roughly ellipsoidal shape of Dimorphos as seen by DRACO suggests that the volume error in the Dimorphos global DTM is probably on the order of several per cent, rather than 23%, by virtue of the simplicity of the long-wavelength shape of Dimorphos. Moreover, in the tests with the largest volume errors, the terminator of the shape model—which is the primary constraint on the extent of the shape into and out of the page—did not match the terminator of the original images. In the case of Dimorphos, the terminator matched nicely between the global DTM and the rendered shape model (Extended Data Fig. 9). Given the results from the shape modelling tests43, this agreement indicated that the volume errors for Dimorphos are probably on the order of several per cent. On the basis of these assessments, a volume error of 5% seems probable. But, we adopted a volume uncertainty of 10% to be conservative. This volume error is larger than the value implied by the reported uncertainties in the X, Y and Z extents, but the intent is to be conservative.

Shape modelling of Didymos

A preliminary SPC shape model of Didymos was built from DRACO data. This model led to the preliminary volume and extent of Didymos reported in Table 1. DRACO images revealed a smaller Z axis and showed that the visible portion of the X axis probably needs to be extended by a few tens of metres compared with a radar-derived shape model12. We used the techniques employed to estimate shape model uncertainties described in the ‘Shape modelling of Dimorphos’ section to assess the preliminary SPC shape model of Didymos. Because the full X and Z extents of Didymos cannot be measured in DRACO images (in contrast to Dimorphos), we conservatively attributed an uncertainty of 15 m to the lengths of all the axes of Didymos.

Impact site identification

The construction of an SPC shape model produced a set of surface landmarks that were used to determine the location (at the time of each image) of the spacecraft relative to the surface in the Dimorphos body fixed frame. This procedure used the spacecraft’s Didymos relative velocity computed by the ground navigation team using a combination of radiometric tracking data (Doppler, range) and optical images of Didymos. We used this information and the shape model of Dimorphos to compute the location of the impact site. The positions of the DART spacecraft from SPC were converted to the inertial J2000 frame and corrected for light time and aberration. The velocity of the spacecraft was estimated by fitting a second-order polynomial function to these positions as a function of time. This approach is identical to that employed previously45 to estimate the position of the Hayabusa spacecraft relative to Itokawa. We used the locations from the last 14 complete images collected by DART to determine the spacecraft velocity. These images contained large numbers of landmarks due to their fine pixel scales, which help anchor the spacecraft position for each image. The fit residuals (that is, the difference between the fit of the spacecraft location to that determined from SPC) are <1 m (Extended Data Fig. 10).

We combined our estimated velocity and the SPC-derived spacecraft positions of the last five images to determine the impact location. Only the last five images were used to locate the impact because heater cycling on the spacecraft introduced cyclic error rates in the inertial attitude knowledge. These rates were estimated to be small at the time of impact owing to the timing of impact relative to heater cycling but at their highest approximately 30 s before impact. Using only the last few images reduced the influence of this known (and cumulative) error source. At each of these five spacecraft positions, we propagated the velocity vector until it intersected with the surface of the Dimorphos shape model. We took the mean of these positions as the impact location. We also computed the spacecraft state at a few different heights above the intersection point to determine the order in which the solar arrays and bus contacted the surface. The propagation from each of the last five images provided the same impact point to within 1  cm. The uncertainties in the impact point location (reported in Table 1) are dominated by the residuals to the fits of the spacecraft position obtained by SPC in the Y and Z directions (Extended Data Fig. 10).

Impact-angle assessments

The tilt angle of the surface relative to the impact velocity vector defines the impact angle. Because the impact site DTM resolves topography at a 5-cm GSD, the spacecraft bus would have interacted with on the order of a thousand facets of the DTM. Therefore, we calculated a mean tilt with respect to the impact velocity vector for each facet in the impact site DTM. The mean tilt calculation is based on previous work41; however, we computed mean tilt with respect to the impact velocity vector, v, rather than the radial vector to a given facet to determine the impact angle.

The magnitude of the mean tilt for a facet at the impact point is the angle between v and the average normal, nav, which is the weighted average of the normal vectors of all facets in the user-defined region, nav = ∑(niAi)/∑Ai surrounding the impact point. In this study, this region was 1.5 m in radius to exceed the size of the spacecraft bus. The normal vector of each facet in the region of interest, ni, was weighted by Ai, the area of the facet projected onto a best-fit plane to the region selected to determine the surface tilt. This yields a mean tilt of acos (nav × v/|nav||v|). The impact angle, θ, which is typically defined relative to the local horizontal for planetary impacts, is given by θ = 90 − (mean tilt). The impact angle shown in Table 1 is computed from the mean tilt and θ of the facet closest to the impact point.

As discussed in the ‘Impact site identification’ section, the location of the impact site has an uncertainty of ±68 cm. To understand the range of mean tilts that DART may have encountered given the uncertainty in the impact location, we considered the distribution of mean tilt for all facets within a circle with radius 68 cm centred on the impact site (Extended Data Fig. 2). The standard deviation of this distribution was 7° (1σ), which we attributed as the uncertainty of the impact angle DART may have experienced.

Boulder counting

The longest axis of each boulder in the final full DRACO image was identified as a line (as done for asteroid (101955) Bennu46). The length of the longest axis was determined from the line length and image pixel scale, assuming the last full image as a ‘flat’ scene. The total number of boulders and pebbles identified in the final full DRACO image impact site is 953 and range in size (that is, the length of the longest axis) from 0.16 m (limit of image resolution, assuming ≥3-pixel sampling47) to 6.5 m. The resulting size–frequency distribution is shown in Extended Data Fig. 6.

Estimates for the density of Dimorphos

To first order, the combined system mass, Msys, was estimated using Kepler’s third law

$${M}_{{\rm{sys}}}=({4\pi }^{2}{a}^{3})/({GP}^{2})$$

where a and P are the pre-impact semimajor axis and the orbit period, respectively, and G is the gravitational constant. Neglecting the aspherical shapes of Didymos and Dimorphos and their associated gravitational potentials may lead to an overestimate of the system mass by about 1–2% (ref. 48). However, this error is negligible as the uncertainty in Dimorphos’s semimajor axis dominates the uncertainty in the system mass. Next, the bulk density of the combined system was obtained by dividing the system mass with the combined volume of both bodies

$${\rho }_{{\rm{sys}}}=({3M}_{{\rm{sys}}})/(4\pi [{R}_{{\rm{A}}}^{3}+{R}_{{\rm{B}}}^{3}])$$

where RA and RB are the volume-equivalent radii of Didymos and Dimorphos, respectively. On the basis of the pre-impact orbit period and semimajor axis6, and the volume-equivalent diameters provided in this work, we calculated a nominal system bulk density of 2,400 ± 250 kg m−3 using the uncertainties quoted in the previous sentence. However, to represent additional possible systematic uncertainties, we adopted a slightly larger uncertainty, which gives a system bulk density of 2,400 ± 300 kg m−3.

The porosity of Dimorphos was estimated as follows:

$${\phi }=1-({\rho }_{{\rm{bulk}}}/{\rho }_{{\rm{grain}}})$$

In this work, we assumed that Dimorphos’s bulk density matches the entire system’s bulk density. There are three near-Earth asteroid systems for which the satellite’s bulk density has been independently measured: 66391 Moshup, 2000 DP107 and 2001 SN263 (two satellites). Two of those satellites (Squannit and 2001 SN263 gamma) were measured to be denser than the primary and two were less dense (2001 SN263 beta and 2000 DP107 beta)20,21,22. For all three systems, the 1σ uncertainties for the satellite and the primary densities overlap. Furthermore, recent work49 estimated that the size of Squannit is about 30% larger than estimated previously20, meaning that the bulk density of Squannit may be in better agreement with its primary. Given these other examples, assuming Dimorphos’s density matches the system bulk density is a reasonable starting point, although its true density could differ substantially from this value. Hera will determine the masses and densities of Didymos and Dimorphos and test the validity of this assumption39.