Stardust from Neutron Stars: How Kilonovae Unveil the Secrets of Heavy Elements ✛ ✦

Original: "Heavy element dust explains the late-time spectra of kilonovae"
arXiv:2607.00433v1 · 2026-07-01 · CC BY 4.0 · 4 min
Heavy elements born in neutron star mergers condense into dust, producing infrared emission that explains observations by the James Webb Space Telescope.

Context

Mergers of neutron stars are the universe's main factories of elements heavier than iron. In 2017, when gravitational-wave detectors (built with contributions from Rainer Weiss) captured the signal GW170817 and telescopes observed the kilonova AT2017gfo, it was a triumph for the theory of nucleosynthesis, whose foundations were laid by Margaret Burbidge and her colleagues. Neutron stars, predicted after the discovery of pulsars by Jocelyn Bell Burnell, eject neutron-rich material upon collision — an ideal site for the r-process. However, late observations of these events, such as 29 days after the explosion, posed a puzzle: why does the spectrum show a strong infrared continuum with a temperature below 1000 K? Atomic lines could not produce such emission, and the idea of carbon dust was ruled out — the ejecta contains too few light elements. It turned out that dust is born directly from the synthesized heavy nuclei.

Methods

To test the hypothesis of 'metallic' dust formation, researchers built a kinetic model of cluster growth in an expanding gas. Tungsten was chosen as a typical element — it is well studied and condenses readily. The equations included atom association into dimers and larger clusters, thermal dissociation, and destruction by fast electrons from beta-decays of radioactive nuclei. Reaction coefficients were taken from molecular dynamics simulations. Then, using the Sedona radiative transfer code, they simulated the kilonova spectrum, adding dust absorption and emission in the Rayleigh limit (grain size much smaller than the wavelength). The ejecta model had a layered structure: slow inner regions and fast outer ones, with velocities reaching 0.12 the speed of light.

Results

Calculations showed that once the temperature drops below 1800 K (about 10 days after the merger), clusters begin to grow rapidly. By day 20, a significant fraction of the refractory elements has condensed into dust — its mass is on the order of 0.001 solar masses. The key point: dust forms only in the slow layers (velocity less than 0.1 the speed of light), where the density is high enough for atoms to stick together. Comparison with spectroscopy data from James Webb for AT2023vfi showed perfect agreement: the model with dust reproduces a smooth blackbody continuum at 660 K, while the dust-free version predicts a set of narrow atomic lines that bear no resemblance to the actual observation. Interestingly, a small peak around 2.1 μm visible in the data is explained by emission from a thin outer shell that remained dust-free.

Implications

The discovery means that late-time infrared spectroscopy of kilonovae becomes a powerful diagnostic tool. From the intensity and shape of the dust continuum, we can determine the mass and composition of the ejecta, and thus measure the amount of heavy elements produced in the merger. This is especially valuable since many kilonovae are too faint for detailed optical analysis. Moreover, similar mechanisms may operate in supernovae and collapsars, opening a new window into the chemical evolution of galaxies and the origin of the elements.

Future development

The next generation of infrared telescopes, such as the upgraded JWST and future observatories, will be able to detect dust emission from more distant and fainter kilonovae. Improving models to include multi-component alloys, non-equilibrium ionization, and reverse shocks will refine the connection between dust properties and merger parameters. Perhaps we will learn to distinguish 'dusty' and 'dust-free' regimes, revealing the diversity of ejecta types.

Impact

The results will impact high-energy astrophysics, observational cosmology, and laboratory nuclear physics. Infrared dust diagnostics will complement traditional methods for studying neutron stars and gravitational-wave sources.

Next steps

The immediate tasks are to perform non-equilibrium radiative transfer calculations that simultaneously account for condensation, ionization, and heating of the ejecta. Also needed are experiments to measure association cross sections of heavy-element clusters at low pressures.

Key open problems

The formation of dust from r-process elements is directly linked to the unsolved problem of the origin of heavy elements in the Universe. Comparing kilonova observations with elemental abundances in old stars and ultra-diffuse galaxies will help understand whether neutron star mergers are sufficient to explain all the observed diversity. Moreover, the survival of dust during interaction with shock waves is part of the broader question of matter cycling in supernova remnants.

The knowledge machine's view

This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.

What the work does

The study constructs a kinetic model of heavy-element cluster growth in the expanding ejecta of a kilonova, using tungsten as an example. With the Sedona radiative transfer code, a spectrum is modeled accounting for dust in the Rayleigh limit, which reproduces the observed infrared continuum of AT2023vfi at a temperature of about 660 K. It is shown that dust forms in the slow-moving ejecta layers, with a mass of approximately 0.001 solar masses.

Where its strength is

This work provides a concrete physical mechanism explaining the late-time infrared spectra of kilonovae and demonstrates how dust formation modeling can be used to extract ejecta parameters from observations.

Where it could go next

  • It might be interesting to compare the obtained results on dust composition and mass with the conclusions about nucleosynthesis in kilonovae presented in the adjacent review, to clarify which heavy elements condense and how this affects the overall picture of galactic enrichment.based on: Cosmic Kitchen: How Neutron Star Collisions Cook Gold arXiv:2605.26308
  • There is a nearby work that thoroughly examines the role of dust in kilonovae; perhaps comparing modeling approaches and spectral interpretation could extend the applicability of your model to other events or wavelength ranges.based on: Precious Metal Dust: How Kilonovae Conceal and Reveal Secrets arXiv:2607.00433
  • It would be interesting to consider how dust formed in a kilonova interacts with stellar winds and disperses into the interstellar medium—this could link your results to the broader picture of matter cycling discussed in the work on stellar winds.based on: Stellar Wind Blows the Ashes Away arXiv:2606.23818

Worth trying

  • Test how sensitive the predicted spectrum is to the choice of element: replace tungsten with osmium or platinum and compare with observations.
  • Apply the model to other kilonovae with known late-time infrared observations to assess the universality of dust formation parameters.
  • Investigate the effect of grain size distribution on the continuum shape, going beyond the Rayleigh approximation for larger grains.

Why this may matter

Such modeling paves the way for using late-time infrared observations of kilonovae to directly determine the mass and composition of the ejected material, which is crucial for understanding the origin of heavy elements in the universe.

Original abstract

as written by the authors arXiv:2607.00433v1

Neutron star mergers are a leading site of $$r$$-process, producing radioactively powered optical and infrared transients known as kilonovae. Observations of the kilonovae AT2017gfo, associated with the gravitational-wave event GW170817, and AT2023vfi, associated with GRB 230307A, have enabled measurements of the mass of ejected $$r$$-process material and the identification of heavy elements in the ejecta. However, late-time observations reveal strong infrared emission with temperature below 1000 K, which is difficult to explain by atomic absorption and emission processes alone. In this paper, we show that kilonova ejecta provide conditions favorable for the formation of dust grains composed of refractory $$r$$-process elements including Zr, W, and Os. We calculate the kinetic formation of dust grains using reaction rate coefficients of W as a proxy, finding that dust forms efficiently, particularly in slow ejecta. This stands in contrast to a previous study that relied on a classical nucleation framework. By performing radiative transfer simulations that incorporate dust formation, we demonstrate that $$r$$-process dust naturally explains the observed late-time infrared emission. The formation and abundance of $$r$$-process dust are highly sensitive to the ejecta mass, composition, and expansion velocity. Infrared emission from $$r$$-process dust can therefore serve a new probe of heavy-element production in neutron star mergers.

Imagine a dust grain made of a thousand tungsten atoms — it's a million times smaller than a grain of sand, yet thanks to its enormous density (about 20 g/cm³) it can absorb all the infrared radiation inside the ejecta. And if such dust formed in a kilonova in our Galaxy, its thermal glow would be visible even in an amateur telescope!

\tau_d \approx 1 \,\left(\frac{M_d}{10^{-3}\,M_\odot}\right) \left(\frac{v_{\rm ej}}{0.1c}\right)^{-2} \left(\frac{t}{29\,\rm day}\right)^{-2}
Shows that a dust mass of about one thousandth of a solar mass is enough to make the ejecta opaque in the infrared band a month after the merger.

Key numbers

  • dust mass: ~10-3 M☉
  • ejecta velocity: 0.08 c (24,000 km/s)
  • dust emission temperature: ~660 K
  • condensation onset time: ~10 days after merger
  • fraction of refractory elements: 3% of ejecta mass
Scientists
Albert EinsteinHans BetheLise MeitnerMargaret BurbidgeBernhard RiemannJoseph Weber
Concepts
kilonova supernova nucleosynthesis ejecta r process neutrino emission neutron star supernova remnant formation stellar nucleosynthesis infrared excess neutron star binary merger high energy neutrino stellar explosion interstellar dust cold gas dust grain infrared gamma ray burst merger chemical evolution stars red giant helium burning X-ray radiation