Main

NASA’s DART successfully impacted Dimorphos, the secondary of the near-Earth binary asteroid (65803) Didymos, on 26 September 2022 at 23:14 UTC. The primary objective of DART was to change the orbital period of Dimorphos around Didymos to demonstrate that a kinetic impactor is a viable method of asteroid deflection1,7. The mission targeted the secondary asteroid in an eclipsing binary system because the experiment could use a single impacting spacecraft and measure the change in the orbit of the secondary through ground-based observations. The Didymos system was selected as the target because it is among the most accessible (low ΔV) of the near-Earth binaries, it has been extremely well characterized4,5,6,8,9,10,11,12 and Dimorphos is in the size range identified as relevant for deflection by a kinetic impactor13,14.

The DART spacecraft collided head-on into the leading hemisphere of Dimorphos to maximize the momentum transfer and reduce the semimajor axis of the Dimorphos orbit, resulting in a shorter orbital period7. If the incident momentum from the impacting spacecraft was simply transferred to the asteroid target with no further momentum enhancement, an orbital period change for Dimorphos of roughly 7 min was expected1. Impact simulations conducted in preparation for DART’s kinetic impact test indicated that, depending on the material strength, impact conditions and other properties, the value of the momentum enhancement factor, β, could be considerable, with predicted values as high as 5 (ref. 2) or 6 (ref. 3), with a resulting orbital period change of more than 40 min (ref. 15).

The Didymos system lightcurve is composed of three parts: the rotational lightcurve of Didymos, the rotational lightcurve of Dimorphos and the mutual events that constrain the orbital period. The Didymos rotational lightcurve can be clearly distinguished because the primary contributes approximately 96% of the light from the system. The Dimorphos rotational period has not been resolved because of its comparatively small size, the oblate shape of Dimorphos16 and the accuracy of the photometric observations necessary for such a detection. Mutual events cause a measurable decrease in the total brightness of the system. We define the primary/secondary occultation or eclipse based on which object is being obscured or shadowed, respectively. We use the timings of the observed mutual events in the determination of the orbital period. For the Didymos–Dimorphos system, mutual events occur when the Didymos–Sun or the Didymos–Earth vector forms an angle less than approximately 17° with the mutual orbit plane of the system. Since the inclination of the mutual orbit to the heliocentric orbit of the binary system is lower than this value, eclipses (mutual shadowing of the components; Fig. 1) always occur. Occultations did not occur during the observing period presented in this paper.

Fig. 1: Post-impact Didymos system geometry.
figure 1

We determine the new orbital period of Dimorphos using the times of mutual events, when a measurable decrease in the system brightness occurs because of an eclipse or occultation. Due to the geometry of the Didymos system during this time period, our lightcurve observations include primary eclipses (left), time outside mutual events (centre) and secondary eclipses (right). These diagrams simulate the view of the system from Earth on 10 October 06:09 (primary eclipse), 10 October 08:47 (outside events) and 10 October 12:06 (secondary eclipse) in geocentric UTC. The primary eclipses observed throughout our post-impact dataset are grazing, which resulted in a subtle decrease in system brightness (Fig. 3). During the secondary eclipse, Dimorphos is completely shadowed.

A precise determination of the pre-impact orbital parameters of the Didymos system was a key goal once the system was chosen as the target of DART. The initial orbit of Dimorphos was first defined following the 2003 apparition when the secondary was discovered11,17. Analyses of lightcurve-derived mutual events obtained during 2003–2022 (ref. 4) led to independent and consistent orbital periods5,6. The data used in the published pre-impact orbit solutions were augmented with photometric data obtained in July 2022 to calculate the pre-impact orbit period for Dimorphos (Extended Data Table 1). Both approaches determined a statistically identical pre-impact orbital period of 11.92148 ± 0.00013 (3σ) h.

To determine the post-impact orbital period, we obtained radar and lightcurve observations of the Didymos system. Our radar observations of Didymos and Dimorphos began about 11 h after impact using the Goldstone X-band (3.5 cm, 8,560 MHz) and continued for 14 dates between UTC 27 September and 13 October (all subsequent dates are in UTC). We also used the Green Bank Telescope to receive radar echoes in a bistatic configuration with transmissions from Goldstone on 2, 6 and 9 October. We obtained echo power spectra during each of the observing windows and range-Doppler images (Fig. 2) on 10 days centred on 4 October, when the signal-to-noise ratios (SNRs) were the highest because Didymos was the closest to Earth. The radar observations of the system are not subject to the same shadowing geometry as the lightcurve photometry. Dimorphos can be seen when illuminated by radar and the system was never in a radar eclipse geometry. We measured the separations between Dimorphos and Didymos in the echo power spectra and the range-Doppler images. We used these measurements in the determination of the orbital parameters of Dimorphos relative to Didymos. We only used data in which the SNRs were strong enough to detect both Didymos and Dimorphos. The first observation of Dimorphos (8σ detection), approximately 12 h after impact, yielded the first estimate of the orbital period change of −36 ± 15 min.

Fig. 2: Radar range-Doppler images of the post-impact Didymos system.
figure 2

Radar range-Doppler images obtained on 4 October using Goldstone and 9 October using Goldstone to transmit and the Green Bank Telescope to receive. In each image, the distance from Earth increases from top to bottom and the Doppler frequency increases to the right, so rotation and orbital motion are anticlockwise. Each image was integrated for 20 min, with 10 min of overlap between successive images. Images have resolutions of 75 m × 0.5 Hz. The broader echo is from Didymos and the smaller, fainter echo shown using arrows is from Dimorphos. The open circles show Dimorphos positions predicted by the pre-impact orbit. The yellow ellipses show the trajectory of Dimorphos. Prediction uncertainties are smaller than the image resolution. On 4 October, the ellipse spans −870 m to +870 m along the y axis and −7 Hz to +7 Hz along the x axis, corresponding to line-of-sight velocity of −12 cm s−1 to +12 cm s−1. On 9 October, the ellipse spans −980 m to +980 m along the y axis and −8 Hz to +8 Hz along the x axis, corresponding to line-of-sight velocity of −14 cm s−1 to +14 cm s−1. The physical extents of the ellipse vary because of the viewing geometry.

Following the DART kinetic impact, ejecta was introduced into the system18. The extra flux and the variable brightness from the rapidly evolving ejecta prevented immediate observations of the mutual events. Lightcurve observations began in the hours after impact and our first successful detection of a mutual event was a secondary eclipse approximately 29.5 h after impact (mid-time at geocentric UTC 28 September 04:50). At the time of the first mutual event detection, the flux from the ejecta dominated the signal within the photometric aperture. This contamination resulted in a reduction in the observed amplitude of the Didymos rotational lightcurve by a factor of 3. The apparent depth of the secondary eclipse was also markedly reduced compared with the predictions6. Pre-impact ejecta models19 suggested that it could take up to several days for our ground-based lightcurve observations to detect the first mutual event as a result of the total ejecta brightness and because the rate of change of that brightness could be comparable with the expected changes in the Didymos system brightness during mutual events.

Photometric observations included in this analysis were obtained from 28 September to 10 October 2022 (Extended Data Table 2). This set of observations ends on 10 October because subsequent observations did not have the required precision because of the bright Moon. On average, our data have photometric accuracy of root mean square (RMS) about 0.006 magnitudes. The exceptional quality of the data included in our analysis has enabled the determination of the Dimorphos orbital period change from lightcurves despite the presence of ejecta in all of our observations (Figs. 3 and 4). At the time of these first observations, the primary eclipses were grazing events (Fig. 1), which required exceptionally precise data to measure.

Fig. 3: Post-impact photometric lightcurve of the Didymos system.
figure 3

Measured photometry from UTC 2 October 2022 phase folded to the 2.26-h rotation period of Didymos (top) and the extracted mutual events (= observed data − 9th order Fourier fit to the rotation of Didymos) phase folded to the new orbit period of Dimorphos (bottom). These lightcurves, collected from five different telescopes, show photometric accuracy similar to all the lightcurve datasets in our analysis. The mutual event times are highly consistent across these datasets, although residual systematics in the photometry result in slightly different event depths.

Fig. 4: Observed mutual events of the Didymos system.
figure 4

The data are marked as circles and the solid curve represents the synthetic lightcurve for the best-fit post-impact solution. The dashed curve is the pre-impact orbit prediction from ref. 6. The primary and secondary events are shown on the left and right sides of the plots, respectively. In some cases, the observations of a secondary event precede those of a primary event (that is, their order in the dataset is the inverse of that shown in the plot). We present these events in reverse order and they are separated by a ‘//’ symbol in the plot (0.4728 days are to be subtracted from the x coordinate of data points to the right from this separator). The y axis shows the magnitude on the night of the observation for each dataset and each tick mark has a range of 0.02 magnitudes.

Two independent methods were used to model the available data for determination of the post-impact orbital period: (1) we use the processes described in ref. 6 to model the lightcurve observations alone and (2) we combine the radar and mutual event timings5,11 plus Didymos-relative astrometry of Dimorphos in optical navigation images from the Didymos Reconnaissance and Asteroid Camera for Optical navigation (DRACO) on the DART spacecraft20. Both methods use the same ground-based photometric datasets but have independent processes for accepting individual data points and mutual events. Ellipsoidal approximations of the shapes of Dimorphos and Didymos are incorporated in the calculation of the orbit period of Dimorphos in both approaches and the axial ratios reported in ref. 16 were used for their calculation.

We determine a post-impact period of 11.372 ± 0.017 (3σ) h with a period change of −33.0 ± 1.0 (3σ) min. Both methods provide statistically identical results for the post-impact orbital period. The rotation period of Didymos is measured during the lightcurve analysis process and shows no variation from its pre-impact value of 2.260 h to an uncertainty of approximately 5 s (3σ). The rotational lightcurve of Dimorphos has not been detected. The new orbital period results in Dimorphos completing an extra full orbit roughly every 9.8 days.

The difference between the pre-impact and post-impact mutual orbit period of the Didymos–Dimorphos system greatly exceeds the approximately 7 min period change calculated for the case of a simple momentum transfer with no momentum enhancement1. Estimates of the change in orbital velocity imparted to Dimorphos require modelling beyond the scope of this paper, but it is evident that the ejecta from the DART impact carried a substantial amount of momentum compared with what the DART spacecraft itself was carrying (for example, ref. 21). To serve as a proof-of-concept for the kinetic impactor technique of planetary defence, DART needed to demonstrate that an asteroid could be targeted during a high-speed encounter16 and that the orbit of the target could be changed. DART has successfully done both.

Methods

The models incorporated three types of observation of the Didymos–Dimorphos system: photometric lightcurves, radar, and Didymos-relative astrometry from the DRACO camera aboard DART20. We determined the post-impact orbital period using two separate models (ref. 6, hereafter SP22, and refs. 5,11, hereafter N22+). Both approaches use the same sets of pre-impact and post-impact lightcurves (Extended Data Tables 1 and 2). The SP22 approach models the lightcurve observations to determine the properties of the orbit. The N22+ approach incorporates Didymos-relative astrometry from DRACO optical navigation images to revise the orbital parameters of the pre-impact orbit and includes lightcurve mutual event timings and radar observations for the post-impact solution (Extended Data Tables 37).

Photometric lightcurve data and reductions

Previous observations of the Didymos system4 demonstrated the need for requirements on the photometry used in the analysis. We define our data quality requirement as an RMS < 0.01 magnitudes, in which the RMS value refers to the consistency over the nightly run and results in a minimum SNR on the individual exposures of about 100. For an accurate decomposition of the lightcurve, we require adequate coverage of the primary lightcurve outside mutual events. We prefer two complete rotation periods of the primary (Prot = 2.26 h) outside the events and estimate this requirement as 6 h of continuous observation. The observations can be split between multiple stations. Four observatories contributed data that met the photometric requirements of the lightcurve dataset for the orbital period change (Extended Data Table 2): Las Campanas Observatory 1-m Swope Telescope, the Las Cumbres Observatory Global Telescope (LCOGT) network 1-m telescopes, the Danish 1.54-m telescope at the European Southern Observatory’s La Silla site and the Lowell Observatory 1.1-m Hall telescope.

The Las Campanas Observatory 1-m Swope telescope is located in the Atacama Desert, Chile22 . The Swope 4K CCD is a visible-wavelength, direct-imaging CCD with a 29.7 × 29.8 arcmin field of view. Swope observations were taken in the Sloan-r’ filter and used sidereal tracking with one or two sky pointings each night. Instrumental aperture photometry was performed on every frame using the Python package SEP23. We use the astroquery Python package to query VizieR24 and Horizons25 databases to identify Gaia stars and to obtain the coordinates of the asteroid for the given date of the images, respectively, and the GaiaXPy Python package, to request and download synthetic photometry of Gaia stars26 in Sloan-r band when available. The Swope data show discrepancies in the photometry (as seen in Fig. 2) at the 0.01–0.02-magnitude level. There are no issues on the timing of the events, which are the key drivers for the derivation of the new orbit period. Additional reductions of these data with optimized apertures will be used to address these discrepancies.

The LCOGT network27 consists of telescopes at seven sites around the world, operated robotically using dynamical scheduling software28. We used the 1-m telescopes at the South Africa and Chile nodes with the telescopes tracking at half of the ephemeris rates. These observations were scheduled and reduced using the NEOexchange Target and Observation Manager and data-reduction pipeline29. Images were pre-processed using the Python-based BANZAI pipeline30. Astrometry and photometry was performed using the Python-based NEOexchange pipeline29. The LCOGT data were primarily obtained in PanSTARRS-w band (equivalent to a broad g+r+i band) and was calibrated to the Gaia DR2 (ref. 31) using calviacat32, with the w band treated as an r band. Calibration stars were constrained to have ‘solar-like’ colours.

The Danish 1.54-m telescope is located at the European Southern Observatory’s La Silla site in Chile. Observations were performed by the MiNDSTEp (Microlensing Network for the Detection of Small Terrestrial Exoplanets) consortium. The Danish Faint Object Spectrograph and Camera (DFOSC) instrument, with a field of view 13.7′ × 13.7′, was used in imaging mode. Images were taken with the Bessell R filter using sidereal tracking. Data reduction used a custom Python pipeline, including alignment of frames using Astrometry.net tools33. Relative photometry was calibrated using the procedure outlined in ref. 34 using the calviacat32 package and the Gaia DR3 star catalogue, with conversion to SDSS-r band magnitudes assuming a colour of (g − r) = 0.52 for Didymos4,35.

The Lowell 1.1-m Hall telescope, located on Anderson Mesa south of Flagstaff, Arizona, is equipped with a 4K × 4K CCD that images a 25-arcmin square field. The telescope was tracked at half of the ephemeris rate. Exposures were taken with a broad VR-band filter. Photometric calibration was based on field star magnitudes from the PanSTARRS catalogue. Only stars with high SNR (>100) and solar-like colours were used for calibration. For the 2 October 2022 data, the photometry was measured using the Canopus software package. For the 5 October 2022, the photometry was measured using the PHOTOMETRYPIPELINE36.

We added lightcurve observations from three telescopes (Extended Data Table 1) to augment the pre-impact lightcurve solutions published in ref. 6 and ref. 5: the 6.5-m Magellan Baade telescope, the SOAR (Southern Astrophysical Research) 4.1-m telescope and the 4.3-m Lowell Discovery Telescope. Both of the revised models confirmed the previous solutions.

Lightcurve decomposition

To model the photometric data of the binary asteroid system, we follow the decomposition methods defined in refs. 17,37 and discussed in ref. 4. Outside mutual events, the largest signal in the Didymos system lightcurve is the flux of the primary, which can be represented by the following Fourier series:

$${F}_{1}\left(t\right)={C}_{1}+\,\mathop{\sum }\limits_{k=1}^{{m}_{1}}\left[{C}_{1k}\cos \frac{2\pi k}{{P}_{1}}\left(t-{t}_{0}\right)+{S}_{1k}\sin \frac{2\pi k}{{P}_{1}}(t-{t}_{0})\right]$$

F1(t) is the flux of the primary, Didymos, at time t, C1 is the mean flux of the primary, C1k and S1k are the Fourier coefficients, P1 is the lightcurve rotational period of Didymos, t0 is the zero-point time and m1 is the maximum significant order. By using this mathematical representation for the system, we assume that Didymos is in principal-axis rotation, that mutual illumination between the objects is negligible and that the rotational lightcurve does not change with time. The lightcurve data are corrected to constant geocentric and heliocentric distances and a consistent solar phase angle. We connect data from different telescopes by scaling them in relative magnitude compared with each other, which has no impact on the timing of the mutual events.

We use observations taken outside mutual events to fit the rotational lightcurve of Didymos. The rapidly changing Earth–Didymos–Sun geometry during this period of Didymos’ close approach to Earth causes observable changes in the primary rotational lightcurve. For our previous work4, we were able to combine data on the timescales of days to weeks. For this dataset, separate decompositions are done for each Julian Day (JD). We correct for the overall fading of the ejecta for each dataset by fitting a linear flux trend before performing the lightcurve decomposition.

Radar observations

We observed Didymos and Dimorphos using the Goldstone X-band radar (3.5 cm, 8,560 MHz) on the 70-m DSS-14 telescope on 14 dates between 27 September and 13 October 2022. On 2, 6 and 9 October, we also used the 100-m Green Bank Telescope to receive radar echoes in a bistatic configuration with transmissions from Goldstone. Typical transmitter power was 430 kW. We obtained echo power spectra during each of the observing windows and range-Doppler images on several days centred on 4 October when the SNRs were the highest. Didymos was clearly detected in all of the data (>3σ) and its maximum bandwidth varied from 22 Hz on 27 September, when the subradar latitude was −50°, to 34 Hz on 13 October, when its subradar latitude was −32° (based on the pole direction estimated by ref. 11).

Detecting Dimorphos was challenging and required experimenting with setups having different frequency resolutions, range resolutions and integration times. This process was a trade-off between obtaining longer integrations with sufficiently high SNRs to detect Dimorphos versus reducing the smearing caused by the orbital motion during the integration. We found that the echo from Dimorphos was most consistently visible at resolutions of 1 Hz in the echo power spectra and at 0.5 Hz in the images. Because of the 11.9-h rotation period, a diameter of about 160 m and a subradar latitude of −50° to −30° (ref. 11), the echo from Dimorphos was expected to have a bandwidth of about 1 Hz (ref. 11), so the data do not resolve Dimorphos in frequency but maximize the SNRs by nearly matching the bandwidth. The contribution of self-noise in the echo power spectra is negligible and does not notably affect the SNRs. We attempted imaging with time delay resolutions of 0.5 μs and 1 μs (corresponding to range resolutions of 75 m and 150 m) and found that the 0.5-μs setup yielded more consistent detections. We experimented with summing data spanning a range of time intervals and found that the echo from Dimorphos was not clearly visible in all the data on any given day. It became more difficult to detect Dimorphos after 4 October, as the distance to Didymos increased and the SNRs correspondingly decreased. Figure 2 shows range-Doppler images and Extended Data Fig. 1 shows selected echo power spectra in which the echo from Dimorphos was seen.

We measured the separations between Dimorphos and Didymos in the echo power spectra and range-Doppler images and used these measurements in the estimation of the orbital parameters of Dimorphos relative to Didymos. The separations in Doppler frequency and range between Didymos and Dimorphos relate to the relative velocity and distance along the line of sight of the observer because of their mutual orbit about each other. We used only data in which both Didymos and Dimorphos were clearly visible for making these measurements. The echo power spectra were processed so that hypothetical echoes from the Didymos system barycentre appear at 0 Hz (ref. 38). Because the reflex motion of Didymos about the system barycentre is <10 m (0.08 Hz)11, we assumed that the Didymos centre of mass (COM) is at 0 Hz, so that the Doppler frequency of Dimorphos represents the relative Doppler shift. The echo from Dimorphos is unresolved, so we assumed that its COM was located in the Doppler bin that contained the strongest spike from the echo from Dimorphos. We assigned uncertainties of ±2 Hz to the Doppler separation measurements to take into account the uncertainties resulting from the frequency resolution of the spectra (1 Hz), the ephemeris errors in the location of the system barycentre (0.24 Hz, 3σ) and the reflex motion of Didymos about the system barycentre (<0.1 Hz). Consequently, the principal source of uncertainty in measurements of the range-Doppler separations are the Doppler frequencies of Dimorphos.

Due to the low SNRs, the COM of Didymos is hard to locate in the range-Doppler images, so we assumed that it is located 375 m (5 range pixels at 75 m per pixel) behind the leading edge, which is the brightest part of the echo and easiest to see. This distance equals the equatorial radius reported from the 3D shape model obtained by ref. 11 and is consistent with preliminary estimates from the DART spacecraft images reported by ref. 16. The echo from Dimorphos extended over one to three range rows and we assumed that its COM is in the trailing row. We assigned uncertainties of 150 m (two range rows) to the range separation measurements. Extended Data Tables 5 and 6 show the range and Doppler frequency of Dimorphos relative to Didymos that were used in the orbit determination. We estimated eight range measurements on 9 October (when reception at Green Bank facilitated detecting echoes from Dimorphos), far more than on any other day, so we inflated their uncertainties by a factor of 3 to mitigate the effects of correlated errors.

Didymos-relative optical astrometry from DRACO images

We measured the positions of Dimorphos relative to Didymos in 16 DRACO images taken in the minutes before impact on 26 September 2022 between 23:10:58.235 and 23:12:39.336 UTC to use in the orbit estimation process. At the time these measurements were made, no shape models estimated from spacecraft images were available to fit to the partially illuminated figures of the two bodies, so we measured the intersections of the limbs with the relative position vectors. These measurements were differenced to estimate the limb-to-limb positions of Dimorphos relative to Didymos. These positions were mapped from image coordinates into right ascension (RA) and declination (DEC) using the camera model and the GNC (guidance, navigation and control) spacecraft attitude knowledge. Measurement uncertainties of 1.13 × 10−3 degrees (3σ) were derived by repeating this process and comparing the different observations. We assumed the equatorial extents of Didymos and Dimorphos to be 425 m and 88 m, respectively, and added an angular distance corresponding to 425 − 88 = 337 m (±20 m (1σ) uncertainty) in the direction of the limb-to-limb separations to estimate the distances between the COMs. Because the measurements covered a very short time span, we de-weighted the uncertainties by 4× (√16) to mitigate effects of correlated measurement errors. We de-weighted the DEC measurements by an extra factor of 2 because they are clearly noisier than the RA measurements. Extended Data Table 7 lists the observations and uncertainties.

Orbital period determination by means of lightcurves (SP22 method)

The numerical model of the Didymos system in ref. 6 was developed using the techniques described in refs. 39,40,41. Didymos and Dimorphos are represented by ellipsoids with axial ratios of a1/c1 = b1/c1 = 1.37, a2/c2 = 1.53 and b2/c2 = 1.50 (ref. 16). The motion of the two bodies is assumed to be Keplerian. The post-impact system was analysed with no a priori assumption on the new binary orbital period. The lightcurve data from 28 and 29 September showed that parts of the data were attenuated with respect to the rotational lightcurve of the primary. Those sections of the data were iteratively masked until all of the data points in the mutual events were identified and the lightcurve decomposition was complete. The first mutual event (0.03 magnitudes deep) was determined to be a secondary eclipse, as the system geometry predicted very shallow or absent primary events.

We adapted the method from ref. 6 to estimate the uncertainty of the post-impact period. When stepping the period over a suitable interval, we computed normalized χ2 for each step. We determined its 3σ uncertainty as an interval in which χ2 is below a certain limit. The adopted limiting P-value corresponds to the probability that the χ2 exceeds a particular value only by chance equal to 0.27%. At each step of the period scanning, the mean anomaly of Dimorphos at the epoch of the impact was also scanned within its 3σ uncertainty interval that was determined by ref. 6 and that we have revised using the extra data taken in July 2022. The SP22 pre-impact period was 11.921478 ± 0.000123 (3σ) h.

The SP22 model determines a post-impact period of 11.372 ± 0.017 (3σ) h, corresponding to an orbit period change of −33.0 ± 1.0 (3σ) min.

Orbital period determination by means of radar and lightcurves (N22+ method)

The lightcurve analysis method described in ref. 5 is a less complicated approach compared with the methods presented in ref. 6. However, it has the advantage of combining information from different data types, such as radar, relative optical astrometry from DRACO images and lightcurve mutual events. The pre-impact orbital period using the N22+ approach was 11.92148 ± 0.00013 (3σ) h.

Lightcurve decomposition was done independently from the SP22 process and required identifying mutual events. The first identified post-impact mutual event was on UTC 28 September 2022. We expected that the head-on impact would decrease the orbital period compared with the pre-impact solution and expected an event with a duration of approximately 1 h. To identify the mutual event, we tested a range of orbit periods from 11 to 12 h in time steps of 0.1 h, with a best match of 11.4 h. Subsequent observations helped refine the initial estimate.

For each mutual event, there are four contact times: when the event begins and flux decreases (T1), when flux reaches a minimum (T2), when the flux begins to increase (T3) and when the event ends and the flux returns to the baseline (T4). We use times T1.5 and T3.5 in the orbit determination. These times are when the flux is at half the total reduction in flux during the event (Fig. 1 in ref. 5). We use 1σ uncertainties of (T1.5 − T1)/2 and (T4 − T3.5)/2 for T1.5 and T3.5, respectively.

We used a least-squares approach, as described in ref. 5, for estimating the orbital parameters of Dimorphos relative to Didymos. Before the DART impact, Dimorphos is assumed to be a point mass on a modified Keplerian orbit around Didymos, with an extra term for modelling the drift in mean motion from nongravitational effects, such as the binary YORP effect and tidal dissipation. The post-impact orbit was assumed to be Keplerian, as the data-arc length is too short to detect a drift in mean motion. We used Δn to capture the change in mean motion owing to the DART impact. The mean anomaly, M, and mean motion, n, of Dimorphos at time, t, are given by:

$$\begin{array}{c}M(t)={M}_{0}+{n}_{0}\left(t-{t}_{0}\right)+\frac{1}{2}\dot{n}{\left(t-{t}_{0}\right)}^{2}\,{\rm{for}}\,t < {t}_{{\rm{imp}}}\\ M(t)={M}_{{\rm{imp}}}+\left({n}_{{\rm{imp}}}+\Delta n\right)\left(t-{t}_{{\rm{imp}}}\right)\,{\rm{for}}\,t > {t}_{{\rm{imp}}}\\ n(t)={n}_{0}+\dot{n}\left(t-{t}_{0}\right)\,{\rm{for}}\,t < {t}_{{\rm{imp}}}\\ n(t)={n}_{{\rm{imp}}}+\Delta n\,{\rm{for}}\,t > {t}_{{\rm{imp}}}\end{array}$$

in which timp is the time of the DART impact, M0 and n0 are the mean anomaly and mean motion at t0, respectively, \(\mathop{n}\limits^{.}\) is the linear drift in mean motion from nongravitational effects and Mimp and nimp are the mean anomaly and mean motion at impact, respectively.

We used differential corrections as described in ref. 5 for estimating the orbital parameters M0, n0, \(\mathop{n}\limits^{.}\), Δn, the pre-impact semimajor axis (a) and the orbit pole longitude (λ) and latitude (β). This requires calculating a computed value corresponding to each observation using a model. We used three kinds of observations: lightcurve mutual event times, radar range and Doppler measurements of Dimorphos relative to Didymos and the separation of Dimorphos from Didymos as seen in spatially resolved DRACO images. The modelling of the first two observables is described in ref. 5. To model the separation of Dimorphos from Didymos in DRACO images, we used SPICE42 to subtract the RA and DEC of the COM of Didymos from those of the COM of Dimorphos as seen from the DART spacecraft.

The N22+ approach results in a post-impact period of 11.371 ± 0.016 (3σ) h and an orbit period change of −33.0 ± 1.0 (3σ) min. The best-fit orbit parameters are presented in Extended Data Table 3.