The Thesan-Zoom project: bursty star formation is incompatible with prolonged dust survival

Cosmic dust is a key regulator of galaxy evolution, but its build-up and survival in the first billion years remain poorly constrained. We present a systematic analysis of dust in the thesan-zoom suite of radiation-hydrodynamical zoom-in simulations, which self-consistently model dust formation, growth, destruction, and its coupling to radiative transfer in galaxies at z at least 3, a multi-phase ISM and bursty star formation histories. The simulated galaxies reproduce the observed trends of dust-to-gas and dust-to-metal ratios with gas metallicity, while showing a dust deficit at high specific star-formation rates. They also broadly match observed dust temperatures and UV-IR spatial offsets. We find that dust and its properties are strongly time-variable and tightly linked to bursty star formation, with short-lived IR-bright phases (median duration of 20.3+2.3-2.4 Myr) and longer dust-poor phases, naturally producing a correlation between dust temperature and distance from the star-forming main sequence. The predicted attenuation at 1500 is low compared to observations, even when including unresolved dust through post processing, indicating that a mechanism able to shield dust from strong feedback events is necessary to reconcile our galaxy formation model with observations. In our model, bursty star formation prevents the survival of large dust reservoirs (M dust / M star at least 10-3) over a significant fraction of cosmic time. This implies that bursty star formation can produce the observed overabundance of UV-bright galaxies at z at least 10 only if it rapidly settles down by z about 8 (where large dust reservoirs are detected). It is also possible that our models lack physical ingredients or emergent phenomena that aid the survival of dust. Future observations of high-redshift dust will be key to diagnose the physical mechanism at play in the first galaxies.