Tidal disruption event rates across cosmic time: forecasts for LSST, Roman, and JWST and their constraints on the supermassive black hole mass function

TL;DR

This paper models the evolution of tidal disruption event rates over cosmic time, forecasting observations for upcoming surveys and proposing a method to constrain the supermassive black hole mass function using TDE data.

Contribution

It introduces a semi-empirical model linking TDE rates to SMBH mass function evolution and incorporates galaxy processes affecting TDE rates across redshifts.

Findings

TDE rate peaks near cosmic noon and declines at higher redshift.

Forecasts for TDE detection rates in LSST, Roman, and JWST surveys.

Methodology to constrain SMBH mass function evolution from TDE observations.

Abstract

Measuring the mass distribution of supermassive black holes (SMBHs) over cosmic time remains particularly challenging for the low mass ($M_{\bullet}\lesssim10^8~M_\odot$) population at $z>1$. This population is also the most sensitive to SMBH seeding and early growth models. In this work we construct a semi-empirical model for the redshift evolution of the TDE rate under multiple SMBH mass function prescriptions, and show that the observed redshift-dependent rate of TDEs is very sensitive to the SMBH mass function and its evolution with redshift. We further incorporate galaxy-scale processes that evolve with redshift -- namely, increasing galaxy nuclear stellar densities, enhanced galaxy-galaxy merger rates, dust obscuration, and a possible top-heavy IMF at early cosmic times -- and quantify their combined impact on the TDE rate. We find that including these effects generally results in a volumetric TDE rate that increases with redshift until a maximum near cosmic noon, before declining at higher redshift where SMBHs that can disrupt stars become increasingly scarce. We forecast TDE rates in the Rubin LSST and the Roman High Latitude Time Domain Survey, alongside expectations for serendipitous TDE rates in the JWST COSMOS-Web survey. Finally, we provide a methodology for using a flux-limited survey of TDEs in LSST to directly constrain the redshift evolution of the SMBH mass function.