Observationally, kilonovae are astrophysical transients powered by the radioactive decay of nuclei heavier than iron, thought to be synthesized in the merger of two compact objects1–4. Over the first few days, the kilonova evolution is dominated by a large number of radioactive isotopes contributing to the heating rate2,5. On timescales of weeks to months, its behaviour is predicted to differ depending on the ejecta composition and the merger remnant6–8. Previous work has shown that the kilonova associated with gamma-ray burst 230307A is similar to kilonova AT2017gfo (ref. 9), and mid-infrared spectra revealed an emission line at 2.15 micrometres that was attributed to tellurium. Here we report a multi-wavelength analysis, including publicly available James Webb Space Telescope data9 and our own Hubble Space Telescope data, for the same gamma-ray burst. We model its evolution up to two months after the burst and show that, at these late times, the recession of the photospheric radius and the rapidly decaying bolometric luminosity (Lbol ∝ t−2.7±0.4, where t is time) support the recombination of lanthanide-rich ejecta as they cool.
