Search for: All records

Despite the increasing number of gravitational wave (GW) detections, the astrophysical origin of binary black hole (BBH) mergers remains elusive. A promising formation channel for BBHs is inside accretion discs around supermassive black holes, that power active galactic nuclei (AGN). In this paper, we test for the first time the spatial correlation between observed GW events and AGN. To this end, we assemble all sky catalogues with 1,412 (242) AGN with a bolometric luminosity greater than 1045.5erg s−1 ($10^{46}\, {\rm erg\, s}^{-1}$) with spectroscopic redshift of z ≤ 0.3 from the Milliquas catalogue, version 7.7b. These AGN are cross-matched with localization volumes of BBH mergers observed in the same redshift range by the LIGO and Virgo interferometers during their first three observing runs. We find that the fraction of the detected mergers originated in AGN brighter than $10^{45.5}\, {\rm erg\, s}^{-1}$ ($10^{46}\, {\rm erg\, s}^{-1}$) cannot be higher than 0.49 (0.17) at a 95 per cent credibility level. Our upper limits imply a limited BBH merger production efficiency of the brightest AGN, while most or all GW events may still come from lower luminosity ones. Alternatively, the AGN formation path for merging stellar-mass BBHs may be actually overall subdominant in the local Universe. To our knowledge, ours are the first observational constraints on the fractional contribution of the AGN channel to the observed BBH mergers.

 

Optical surveys have become increasingly adept at identifying candidate tidal disruption events (TDEs) in large numbers, but classifying these generally requires extensive spectroscopic resources. Here we presenttdescore, a simple binary photometric classifier that is trained using a systematic census of ∼3000 nuclear transients from the Zwicky Transient Facility (ZTF). The sample is highly imbalanced, with TDEs representing ∼2% of the total.tdescoreis nonetheless able to reject non-TDEs with 99.6% accuracy, yielding a sample of probable TDEs with recall of 77.5% for a precision of 80.2%.tdescoreis thus substantially better than any available TDE photometric classifier scheme in the literature, with performance not far from spectroscopy as a method for classifying ZTF nuclear transients, despite relying solely on ZTF data and multiwavelength catalog cross matching. In a novel extension, we use “Shapley additive explanations” to provide a human-readable justification for each individualtdescoreclassification, enabling users to understand and form opinions about the underlying classifier reasoning.tdescorecan serve as a model for photometric identification of TDEs with time-domain surveys, such as the upcoming Rubin observatory.

 

We present fundamental scaling relationships between properties of the optical/UV light curves of tidal disruption events (TDEs) and the mass of the black hole that disrupted the star. We have uncovered these relations from the late-time emission of TDEs. Using a sample of 63 optically selected TDEs, the latest catalogue to date, we observed flattening of the early-time emission into a near-constant late-time plateau for at least two-thirds of our sources. Compared to other properties of the TDE light curves (e.g. peak luminosity or decay rate) the plateau luminosity shows the tightest correlation with the total mass of host galaxy (p-value of 2 × 10−6, with a residual scatter of 0.3 dex). Physically this plateau stems from the presence of an accretion flow. We demonstrate theoretically and numerically that the amplitude of this plateau emission is strongly correlated with black hole mass. By simulating a large population (N = 106) of TDEs, we determine a plateau luminosity-black hole mass scaling relationship well described by $\log _{10} \left({{M_{\bullet }}/M_\odot }\right) = 1.50 \log _{10} \left({ L_{\rm plat}}/10^{43} \, {\rm erg\, s^{-1}}\right) + 9.0$ (here Lplat is measured at 6 × 1014 Hz in the rest frame). The observed plateau luminosities of TDEs and black hole masses in our large sample are in excellent agreement with this simulation. Using the black hole mass predicted from the observed TDE plateau luminosity, we reproduce the well-known scaling relations between black hole mass and galaxy velocity dispersion. The large black hole masses of 10 of the TDEs in our sample allow us to provide constraints on their black hole spins, favouring rapidly rotating black holes. Finally, we also discover two significant correlations between early time properties of optical TDE light curves (the g-band peak luminosity and radiated energy) and the TDEs black hole mass.

 

The origin of cosmic high-energy neutrinos remains largely unexplained. For high-energy neutrino alerts from IceCube, a coincidence with time-variable emission has been seen for three different types of accreting black holes: (1) a gamma-ray flare from a blazar (TXS 0506+056), (2) an optical transient following a stellar tidal disruption event (TDE; AT2019dsg), and (3) an optical outburst from an active galactic nucleus (AGN; AT2019fdr). For the latter two sources, infrared follow-up observations revealed a powerful reverberation signal due to dust heated by the flare. This discovery motivates a systematic study of neutrino emission from all supermassive black hole with similar dust echoes. Because dust reprocessing is agnostic to the origin of the outburst, our work unifies TDEs and high-amplitude flares from AGN into a population that we dub accretion flares. Besides the two known events, we uncover a third flare that is coincident with a PeV-scale neutrino (AT2019aalc). Based solely on the optical and infrared properties, we estimate a significance of 3.6σ for this association of high-energy neutrinos with three accretion flares. Our results imply that at least ∼10 per cent of the IceCube high-energy neutrino alerts could be due to accretion flares. This is surprising because the sum of the fluence of these flares is at least three orders of magnitude lower compared to the total fluence of normal AGN. It thus appears that the efficiency of high-energy neutrino production in accretion flares is increased compared to non-flaring AGN. We speculate that this can be explained by the high Eddington ratio of the flares.

 

We predict late-time optical/UV emission from tidal disruption events (TDEs) from our slim accretion disc model (Wen et al. 2020) and explore the impact of the black hole mass M•, black hole spin a•, and accretion disc size. We use these synthetic spectra to successfully fit the multiband Swift observations of ASASSN-14li at >350 d, setting only the host galaxy extinction and outer disc radius as free parameters and employing the M•, a•, disc inclination, and disc accretion rates derived from fitting 10 epochs of ASASSN-14li’s X-ray spectra with the slim disc. To address the nature of the early-time optical/UV emission, we consider two models: shock dissipation and reprocessing. We find that (1) the predicted late-time optical/UV colour (e.g. u − w2) is insensitive to black hole and disc parameters unless the disc spreads quickly; (2) a starburst galaxy extinction model is required to fit the data, consistent with ASASSN-14li’s post-starburst host; (3) surprisingly, the outer disc radius is ≈2 × the tidal radius and ∼constant at late times, showing that viscous spreading is slow or non-existent; (4) the shock model can be self-consistent if M• ≲ 106.75 M⊙, i.e. on the low end of ASASSN-14li’s M• range (106.5–7.1 M⊙; 1σ CL); larger black hole masses require disruption of an unrealistically massive progenitor star; (5) the gas mass needed for reprocessing, whether by a quasi-static or an outflowing layer, can be <0.5 M⊙, consistent with a (plausible) disruption of a solar-mass star.

 

We present an expansion of FLEET, a machine-learning algorithm optimized to select transients that are most likely tidal disruption events (TDEs). FLEET is based on a random forest algorithm trained on both the light curves and host galaxy information of 4779 spectroscopically classified transients. We find that for transients with a probability of being a TDE,P(TDE) > 0.5, we can successfully recover TDEs with ≈40% completeness and ≈30% purity when using their first 20 days of photometry or a similar completeness and ≈50% purity when including 40 days of photometry, an improvement of almost 2 orders of magnitude compared to random selection. Alternatively, we can recover TDEs with a maximum purity of ≈80% and a completeness of ≈30% when considering only transients withP(TDE) > 0.8. We explore the use of FLEET for future time-domain surveys such as the Legacy Survey of Space and Time on the Vera C. Rubin Observatory (Rubin) and the Nancy Grace Roman Space Telescope (Roman). We estimate that ∼10⁴well-observed TDEs could be discovered every year by Rubin and ∼200 TDEs by Roman. Finally, we run FLEET on the TDEs from our Rubin survey simulation and find that we can recover ∼30% of them at redshiftz< 0.5 withP(TDE) > 0.5, or ∼3000 TDEs yr^–1that FLEET could uncover from the Rubin stream. We have demonstrated that we will be able to run FLEET on Rubin photometry as soon as this survey begins. FLEET is provided as an open source package on GitHub: https://github.com/gmzsebastian/FLEET.

 

We conduct a systematic tidal disruption event (TDE) demographics analysis using the largest sample of optically selected TDEs. A flux-limited, spectroscopically complete sample of 33 TDEs is constructed using the Zwicky Transient Facility over 3 yr (from 2018 October to 2021 September). We infer the black hole (BH) mass (M_BH) with host galaxy scaling relations, showing that the sampleM_BHranges from 10^5.1M_⊙to 10^8.2M_⊙. We developed a survey efficiency corrected maximum volume method to infer the rates. The rest-frameg-band luminosity function can be well described by a broken power law of$\varphi (L_g)\propto {\left[{\left(L_g/L_{\mathrm{bk}}\right)}^{0.3}+{\left(L_g/L_{\mathrm{bk}}\right)}^{2.6}\right]}^{-1}$, withL_bk= 10^43.1erg s⁻¹. In the BH mass regime of 10^5.3≲ (M_BH/M_⊙) ≲ 10^7.3, the TDE mass function follows$\varphi (M_{\mathrm{BH}})\propto M_{\mathrm{BH}}^{-0.25}$, which favors a flat local BH mass function (${\mathit{dn}}_{\mathrm{BH}}/d\log M_{\mathrm{BH}}\approx \mathrm{constant}$). We confirm the significant rate suppression at the high-mass end (M_BH≳ 10^7.5M_⊙), which is consistent with theoretical predictions considering direct capture of hydrogen-burning stars by the event horizon. At a host galaxy mass ofM_gal∼ 10¹⁰M_⊙, the average optical TDE rate is ≈3.2 × 10⁻⁵galaxy⁻¹yr⁻¹. We constrain the optical TDE rate to be [3.7, 7.4, and 1.6] × 10⁻⁵galaxy⁻¹yr⁻¹in galaxies with red, green, and blue colors.

 

Abstract Tidal disruption events (TDEs) offer a unique way to study dormant black holes. While the number of observed TDEs has grown thanks to the emergence of wide-field surveys in the past few decades, questions regarding the nature of the observed optical, UV, and X-ray emission remain. We present a uniformly selected sample of 30 spectroscopically classified TDEs from the Zwicky Transient Facility Phase I survey operations with follow-up Swift UV and X-ray observations. Through our investigation into correlations between light-curve properties, we recover a shallow positive correlation between the peak bolometric luminosity and decay timescales. We introduce a new spectroscopic class of TDE, TDE-featureless, which are characterized by featureless optical spectra. The new TDE-featureless class shows larger peak bolometric luminosities, peak blackbody temperatures, and peak blackbody radii. We examine the differences between the X-ray bright and X-ray faint populations of TDEs in this sample, finding that X-ray bright TDEs show higher peak blackbody luminosities than the X-ray faint subsample. This sample of optically selected TDEs is the largest sample of TDEs from a single survey yet, and the systematic discovery, classification, and follow-up of this sample allows for robust characterization of TDE properties, an important stepping stone looking forward toward the Rubin era.