Stardust from Neutron Stars: How Kilonovae Unveil the Secrets of Heavy Elements ✦
Neutron star mergers are the main source of r-process elements, producing kilonovae — optical and infrared transients powered by radioactive heating. Observations of kilonovae AT2017gfo (associated with GW170817) and AT2023vfi (associated with GRB 230307A) have allowed measuring the mass of ejected r-process material and identifying heavy elements. However, late observations reveal strong infrared emission with temperatures below 1000 K, which is difficult to explain by atomic processes alone. This work demonstrates that conditions in kilonova ejecta are favorable for the formation of dust grains from refractory r-process elements, including Zr, W, and Os. Using kinetic calculations with reaction rate coefficients for tungsten, it is shown that dust formation is efficient, especially in slow ejecta, contrary to a previous study based on classical nucleation theory. Radiation transfer modeling that includes dust naturally explains the observed late infrared emission. The formation and abundance of r-process dust depend strongly on the mass, composition, and expansion velocity of the ejecta.
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.
🎯 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!
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