Blazar Variability with the Vera C. Rubin Legacy Survey of Space and Time

With their emission mainly coming from a relativistic jet pointing toward us, blazars are fundamental sources for studying extragalactic jets and their central engines, consisting of supermassive black holes fed by accretion disks. They are also candidate sources of high-energy neutrinos and cosmic rays. Because of the jet orientation, the nonthermal blazar emission is Doppler beamed; its variability is unpredictable, and it occurs on timescales from less than 1 hr to years. Comprehension of the diverse mechanisms producing the flux and spectral changes requires well-sampled multiband light curves over long time periods. In particular, outbursts are the best test bench for shedding light on the underlying physics, especially when studied in a multiwavelength context. The Vera C. Rubin Legacy Survey of Space and Time (Rubin-LSST) will monitor the southern sky for 10 yr in six photometric bands, offering a formidable tool for studying blazar variability features in a statistical way. The alert system will allow us to trigger follow-up observations of outstanding events, especially at high (keV-to-GeV) and very high (TeV) energies. We here examine the simulated Rubin-LSST survey strategies with the aim of understanding which cadences are more suitable for blazar variability science. Our metrics include light curve and color sampling. We also investigate the problem of saturation, which will affect the brightest and many flaring sources, and will have a detrimental impact on follow-up observations.