SN 2018bsz: significant dust formation in a nearby superluminous supernova

We investigate the thermal emission and the extinction from dust associated with the nearby (z = 0.0267) superluminous supernova (SLSN) 2018bsz. Our dataset combines daily-cadence, simultaneous optical and near-infrared coverage up to ~100 days with late-time (+1.7 yr) mid-infrared observations. At 230 days after light-curve peak, the SN is not detected in the optical, but shows a surprisingly strong near-infrared excess, with r - J > 3 mag and r - Ks > 5 mag. The temporal evolution of the infrared light curve allows us to assess whether the mid-infrared emission originates from newly formed dust, either in the SN ejecta or in a post-shock circumstellar shell, or instead from an infrared echo produced by pre-existing circumstellar or interstellar dust heated by the SN radiation. Our radiative-transfer modelling shows that the circumstellar and interstellar echo scenarios cannot simultaneously reproduce the observed infrared luminosity and its temporal evolution, and we therefore disfavour these interpretations. In contrast, a scenario in which new dust forms in the SN ejecta at epochs >200 days provides a self-consistent explanation for the evolution of the SN flux. We can fit the spectral energy distribution well at +230 d with 5x10^-4 solar masses of carbon dust, increasing over the following several hundred days to 10^-2 solar masses by +535 d. SN 2018bsz is the first SLSN showing evidence for dust formation within the SN ejecta, and appears to form ten times more dust than normal core-collapse SNe at similar epochs. Together with their preference for low-mass, low-metallicity host galaxies, we suggest that SLSNe may be a significant contributor to dust formation in the early Universe.