Time-Dependent Radiation Transport Simulations of Infrared Echoes from Dust-Shrouded Luminous Transients

A wide range of stellar explosions, including supernovae (SNe), tidal disruption events (TDE), and fast blue optical transients (FBOT), can occur in dusty environments initially opaque to the transient's optical/UV light, becoming visible only once the dust is destroyed by the transient's rising luminosity. We present axisymmetric time-dependent radiation transport simulations of dust-shrouded transients with \texttt{Athena++} and tabulated gray opacities, which predict the light-curves of the dust-reprocessed infrared (IR) radiation. The luminosity and timescale of the IR light-curve depends on whether the transient rises rapidly or slowly compared to the light crossing-time of the photosphere, $t_{\rm lc}$. For slow-rising transients ($t_{\rm rise} \gg t_{\rm lc}$) such as SNe, the reprocessed IR radiation diffuses outwards through the dust shell faster than the sublimation front expands; the IR light-curve therefore begins rising prior to the escape of UV/optical light, but peaks on a timescale $\sim t_{\rm rise}$ shorter than the transient duration. By contrast, for fast-rising transients ($t_{\rm rise} \ll t_{\rm lc}$) such as FBOTs and some TDEs, the finite light-travel time results in the reprocessed radiation arriving as an ``echo'' lasting much longer than the transient itself (despite the dust photosphere having already being destroyed by peak light). We explore the effects of the system geometry by considering a torus-shaped distribution of dust. The IR light-curves seen by observers in the equatorial plane of the torus resemble those for a spherical dust shell, while polar observers see faster-rising, brighter and shorter-lived emission. We successfully model the IR excess seen in AT2018cow as a dust echo, supporting the presence of an opaque dusty medium surrounding FBOTs prior to explosion.