Quantifying the Impact of Relativistic Precession on Tidal Disruption Event Light Curves

The disruption of a star by the tidal forces of a spinning black hole causes the stellar stream to precess, affecting the conditions for triggering the tidal disruption event (TDE). In this work, we study the effect that precession imprints on TDE light curves due to the interaction of the TDE wind and luminosity with the stream wrapped around the black hole. We perform two-dimensional radiation-hydrodynamic simulations using the moving-mesh hydrodynamic code jet with its radiation treatment module. We study the impact of black hole mass, accretion efficiency, and inclination between the orbital and spin planes. From our results, we identified two behaviours: (i) models with low-mass black holes (Mh ∼ 106 M⊙), low inclination (i ∼ 0), and low accretion efficiency (η ∼ 0.01) show light curves with a short early peak caused by the interaction of the wind with the inner edge of the stream. The line of sight has little effect on the light curve, since the stream covers a small fraction of the solid angle due to the precession occurring in the orbital plane; and (ii) models with high-mass black holes (Mh ≳ 107 M⊙), high inclination (i ∼ 90°), and high accretion efficiency (η ∼ 0.1) produce light curves with luminosity peaks that can be delayed by up to 50–100 d depending on the line of sight due to presence of the precessed stream blocking the radiation in the early phase of the event. Our results show that black hole spin and misalignment do not imprint recognizable features on the light curves but rather can add complications to their analysis.

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.

Context. The strong tidal force in a supermassive black hole’s (SMBH) vicinity, coupled with a higher stellar density at the center of a galaxy, make it an ideal location to study the interaction between stars and black holes. Two stars moving near the SMBH could collide at a very high speed, which can result in a high energy flare. The resulting debris can then accrete onto the SMBH, which could be observed as a separate event. Aims. We simulate the light curves resulting from the fallback accretion in the aftermath of a stellar collision near a SMBH. We investigate how it varies with physical parameters of the system. Methods. Light curves are calculated by simulating post-collision ejecta as N particles moving along individual orbits which are determined by each particle’s angular momentum, and assuming that all particles start from the distance from the black hole at which the two stars collided. We calculate how long it takes for each particle to reach its distance of closest approach to the SMBH, and from there we add to it the viscous accretion timescale as described by the alpha-disk model for accretion disks. Given a timestamp for each particle to accrete, this can be translated into into a luminosity for a given radiative efficiency. Results. With all other physical parameters of the system held constant, the direction of the relative velocity vector at time of impact plays a large role in determining the overall form of the light curve. One distinctive light curve we notice is characterized by a sustained increase in the luminosity some time after accretion has started. We compare this form to the light curves of some candidate tidal disruption events (TDEs). Conclusions. Stellar collision accretion flares can take on unique appearances that would allow them to be easily distinguished, as well as elucidate underlying physical parameters of the system. There exist several ways to distinguish these events from TDEs, including the much wider range of SMBH masses stellar collisions may exist around. The beginning of the Vera Rubin Observatory Legacy Survey of Space and Time will greatly improve survey abilities and facilitate in the identification of more stellar collision events, particularly in higher-mass SMBH systems.

In a supermassive black hole (BH) tidal disruption event (TDE), the tidally disrupted star feeds the BH via an accretion disk. Most often it is assumed that the accretion rate history, hence the emission light curve, tracks the rate at which new debris mass falls back onto the disk, notably the t−5/3 power law. But this is not the case when the disk evolution due to viscous spreading - the driving force for accretion - is carefully considered. We construct a simple analytical model that comprehensively describes the accretion rate history across 4 different phases of the disk evolution, in the presence of mass fallback and disk wind loss. Accretion rate evolves differently in those phases which are governed by how the disk heat energy is carried away, early on by advection and later by radiation. The accretion rate can decline as steeply as t−5/3 only if copious disk wind loss is present during the early advection-cooled phase. Later, the accretion rate history is t−8/7 or shallower. These have great implications on the TDE flare light curve. A TDE accretion disk is most likely misaligned with the equatorial plane of the spinning BH. Moreover, in the TDE the accretion rate is super- or near-Eddington thus the disk is geometrically thick, for which case the BH’s frame dragging effect may cause the disk precess as a solid body, which may manifest itself as quasi-periodic signal in the TDE light curve. Our disk evolution model predicts the disk precession period increases with time, typically as ∝ t. The results are applied to the recently jetted TDE flare Swift transient J1644 + 57 which shows numerous, quasi-periodic dips in its long-term X-ray light curve. As the current TDE sample increases, the identification of the disk precession signature provides a unique way of measuring BH spin and studying BH accretion physics.

Outflows are inevitably driven from the disk if the vertical component of the black hole (BH) gravity cannot resist the radiation force. We derive the mass loss rate in the outflows by solving a dynamical equation for the vertical gas motion in the disk. The structure of a supercritical accretion disk is calculated with the radial energy advection included. We find that most inflowing gas is driven into outflows if the disk is accreting at a moderate Eddington-scaled rate (up to $\sim 100$) at its outer edge, i.e., only a small fraction of gas is accreted by the BH, which is radiating at several Eddington luminosities, while it reaches around ten for extremely high accretion rate cases ($\dot{m}\equiv\dot{M}/\dot{M}_{\rm Edd}\sim 1000$). Compared with a normal slim disk, the disk luminosity is substantially suppressed due to the mass loss in the outflows. We apply the model to the light curves of the tidal disruption events (TDEs), and find that the disk luminosity declines very slowly with time even if a typical accretion rate $\dot{m}\propto t^{-5/3}$ is assumed at the outer edge of the disk, which is qualitatively consistent with the observed light curves in some TDEs, and helps understanding the energy deficient phenomenon observed in the TDEs. Strong outflows from supercritical accretion disks surrounding super massive BHs may play crucial roles on their host galaxies, which can be taken as an ingredient in the mechanical feedback models. The implications of the results on the growth of supper-massive BHs are also discussed.

Tidal disruption events (TDEs) can uncover the quiescent black holes (BHs) at the center of galaxies and also offer a promising method to study them. In a partial TDE (PTDE), the BH’s tidal force cannot fully disrupt the star, so the stellar core survives and only a varied portion of the stellar mass is bound to the BH and feeds it. We calculate the event rate of PTDEs and full TDEs (FTDEs). In general, the event rate of PTDEs is higher than that of FTDEs, especially for the larger BHs, and the detection rate of PTDEs is approximately dozens per year, as observed by the Zwicky Transient Factory. During the circularization process of the debris stream in PTDEs, no outflow can be launched due to the efficient radiative diffusion. The circularized debris ring then experiences viscous evolution and forms an accretion disk. We calculate the light curves of PTDEs contributed by these two processes, along with their radiation temperature evolution. The light curves have double peaks and peak in the UV spectra. Without obscuration or reprocessing of the radiation by an outflow, PTDEs provide a clean environment to study the circularization and transient disk formation in TDEs.

Tidal Disruption Events (TDEs) provide an opportunity to study supermassive black holes that are otherwise quiescent. The Vera C. Rubin Legacy Survey of Space and Time will be capable of discovering thousands of TDEs each year, allowing for a dramatic increase in the number of discovered TDEs. The optical light curves from TDEs can be used to model the physical parameters of the black hole and disrupted star, but the sampling and photometric uncertainty of the real data will couple with model degeneracies to limit our ability to recover these parameters. In this work, we aim to model the impact of the Rubin survey strategy on simulated TDE light curves to quantify the typical errors in the recovered parameters. black hole masses <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mn>5.5</mml:mn> <mml:mo>&lt;</mml:mo> <mml:mi>log</mml:mi> <mml:msub> <mml:mrow> <mml:mi>M</mml:mi> </mml:mrow> <mml:mrow> <mml:mi mathvariant="normal">BH</mml:mi> </mml:mrow> </mml:msub> <mml:mo>/</mml:mo> <mml:msub> <mml:mrow> <mml:mi>M</mml:mi> </mml:mrow> <mml:mrow> <mml:mo>⊙</mml:mo> </mml:mrow> </mml:msub> <mml:mo>&lt;</mml:mo> <mml:mn>8.2</mml:mn> </mml:math> can be recovered with typical errors of 0.26 dex, with early coverage removing large outliers. Recovery of the mass of the disrupted star is difficult, limited by the degeneracy with the accretion efficiency. Only 57% of the cases have accurate recovery of whether the events are full or partial, so we caution the use this method to assess whether TDEs are partially or fully disrupted systems. black hole mass measurements obtained from Rubin observations of TDEs will provide powerful constraints on the black hole mass function, black hole–galaxy co-evolution, and the population of black hole spins, though continued work to understand the origin of TDE observables and how the TDE rate varies among galaxies will be necessarily to fully utilize the upcoming rich dataset from Rubin.

Outflows are inevitably driven from the disk if the vertical component of the black hole (BH) gravity cannot resist the radiation force. We derive the mass-loss rate in the outflows by solving a dynamical equation for the vertical gas motion in the disk. The structure of a supercritical accretion disk is calculated with the radial energy advection included. We find that most inflowing gas is driven into outflows if the disk is accreting at a moderate Eddington-scaled rate (up to ∼100) at its outer edge, i.e., only a small fraction of gas is accreted by the BH, which is radiating at several Eddington luminosities, while it reaches around ten for extremely high accretion rate cases (). Compared with a normal slim disk, the disk luminosity is substantially suppressed due to the mass loss in the outflows. We apply the model to the light curves of the tidal disruption events (TDEs) and find that the disk luminosity declines very slowly with time even if a typical accretion rate is assumed at the outer edge of the disk, which is qualitatively consistent with the observed light curves in some TDEs and helps us to understand the energy deficient phenomenon observed in the TDEs. Strong outflows from supercritical accretion disks surrounding supermassive BHs may play crucial roles in their host galaxies, which can be taken as an ingredient in the mechanical feedback models. The implications of the results on the growth of supermassive BHs are also discussed.

We report an observational estimate of the rate of stellar tidal disruption flares (TDFs) in inactive galaxies, based on a successful search for these events among transients in galaxies using archival SDSS multi-epoch imaging data (Stripe 82). This search yielded 186 nuclear flares in galaxies, of which two are excellent TDF candidates. Because of the systematic nature of the search, the very large number of galaxies, the long time of observation, and the fact that non-TDFs were excluded without resorting to assumptions about TDF characteristics, this study provides an unparalleled opportunity to measure the TDF rate. To compute the rate of optical stellar tidal disruption events, we simulate our entire pipeline to obtain the efficiency of detection. The rate depends on the light curves of TDFs, which are presently still poorly constrained. Using only the observed part of the SDSS light curves gives a model-independent upper limit to the optical TDF rate: < 2 10^-4 per year per galaxy (90% CL). We develop three empirical models of the light curves, based on the two SDSS light curves and two more recent and better-sampled Pan-STARRS TDF light curves, leading to our best-estimate of the rate: (1.5 - 2.0)_{-1.3}^{+2.7} 10^-5 per year per galaxy. We explore the modeling uncertainties by considering two theoretically motivated light curve models, as well as two different relationships between black hole mass and galaxy luminosity, and two different treatments of the cutoff in the visibility of TDFs at large black hole mass. From this we conclude that these sources of uncertainty are not significantly larger than the statistical ones. Our results are applicable for galaxies hosting black holes with mass in the range of few million to 10^8 solar masses, and translates to a volumetric TDF rate of (4 - 8) 10^-8 per year per cubic Mpc.

We present an analysis of the optical and ultraviolet properties of AT 2020wey, a faint and fast tidal disruption event (TDE) at 124.3 Mpc. The light curve of the object peaked at an absolute magnitude of Mg = −17.45 ± 0.08 mag and a maximum bolometric luminosity of Lpeak = (8.74 ± 0.69)×1042 erg s−1, making it comparable to iPTF16fnl, the faintest TDE to date. The time from the last non-detection to the g-band peak is 23 ± 2 days, and the rise is well described by L ∝ t1.80 ± 0.22. The decline of the bolometric light curve is described by a sharp exponential decay steeper than the canonical t−5/3 power law, making AT 2020wey the fastest declining TDE to date. The multi-band light curve analysis shows first a slowly declining blackbody temperature of TBB ∼ 20 000 K around the peak brightness followed by a gradual temperature increase. The blackbody photosphere is found to expand at a constant velocity (∼1300 km s−1) to a value of RBB ∼ 3.5 × 1014 cm before contracting rapidly. Multi-wavelength fits to the light curve indicate a complete disruption of a star of M⋆ = 0.11−0.02+0.05 M⊙ by a black hole of MBH = 106.46−0.09+0.09 M⊙. Our spectroscopic dataset reveals broad (∼104 km s−1) Balmer and He II 4686 Å lines, with Hα reaching its peak with a lag of ∼8.2 days compared to the continuum. In contrast to previous faint and fast TDEs, there are no obvious Bowen fluorescence lines in the spectra of AT 2020wey. There is a strong correlation between the MOSFIT-derived black hole masses of TDEs and their decline rate. However, AT 2020wey is an outlier in this correlation, which could indicate that its fast early decline may be dictated by a different physical mechanism than fallback. After performing a volumetric correction to a sample of 30 TDEs observed between 2018 and 2020, we conclude that faint TDEs are not rare by nature; they should constitute up to ∼50–60% of the entire population and their numbers could alleviate some of the tension between the observed and theoretical TDE rate estimates. We calculate the optical TDE luminosity function and we find a steep power-law relation dN/dLg ∝ Lg−2.36±0.16.

The unusual transient Swift J1644+57 likely resulted from a collimated relativistic jet powered by accretion onto a massive black hole (BH) following the tidal disruption (TD) of a star. Several mysteries cloud the interpretation of this event: (1) extreme flaring and `plateau' shape of the X-ray/gamma-ray light curve during the first 10 days after the gamma-ray trigger; (2) unexpected rebrightening of the forward shock radio emission months after trigger; (3) no obvious evidence for jet precession, despite misalignment typically expected between the angular momentum of the accretion disk and BH; (4) recent abrupt shut-off in jet X-ray emission after 1.5 years. Here we show that all of these seemingly disparate mysteries are naturally resolved by one assumption: the presence of strong magnetic flux Phi threading the BH. Initially, Phi is weak relative to high fall-back mass accretion rate, Mdot, and the disk and jets precess about the BH axis = our line of sight. As Mdot drops, Phi becomes dynamically important and leads to a magnetically-arrested disk (MAD). MAD naturally aligns disk and jet axis along the BH spin axis, but only after a violent rearrangement phase (jet wobbling). This explains the erratic light curve at early times and the lack of precession at later times. We use our model for Swift J1644+57 to constrain BH and disrupted star properties, finding that a solar-mass main sequence star disrupted by a relatively low mass, M~10^5-10^6 Msun, BH is consistent with the data, while a WD disruption (though still possible) is disfavored. The magnetic flux required to power Swift J1644+57 is too large to be supplied by the star itself, but it could be collected from a quiescent `fossil' accretion disk present in the galactic nucleus prior to the TD. The presence (lack of) of such a fossil disk could be a deciding factor in what TD events are accompanied by powerful jets.[abridged]

Dynamical perturbations from supermassive black hole (SMBH) binaries can increase the rates of tidal disruption events (TDEs). However, most previous work focuses on TDEs from the heavier black hole in the SMBH binary (SMBHB) system. In this work, we focus on the lighter black holes in SMBHB systems and show that they can experience a similarly dramatic increase in their TDE rate due to perturbations from a more massive companion. While the increase in TDEs around the more massive black hole is mostly due to chaotic orbital perturbations, we find that, around the smaller black hole, the eccentric Kozai–Lidov mechanism is dominant and capable of producing a comparably large number of TDEs. In this scenario, the mass derived from the light curve and spectra of TDEs caused by the lighter SMBH companion is expected to be significantly smaller than the SMBH mass estimated from galaxy scaling relations, which are dominated by the more massive companion. This apparent inconsistency can help find SMBHB candidates that are not currently accreting as active galactic nuclei and that are at separations too small for them to be resolved as two distinct sources. In the most extreme cases, these TDEs provide us with the exciting opportunity to study SMBHBs in galaxies where the primary SMBH is too massive to disrupt Sun-like stars.

The origin of thermal optical and UV emission from stellar tidal disruption flares (TDFs) remains an open question. We present Hubble Space Telescope far-UV (FUV) observations of eight optical/UV-selected TDFs 5–10 yr post-peak. Six sources are cleanly detected, showing point-like FUV emission ( ) from the centers of their host galaxies. We discover that the light curves of TDFs from low-mass black holes (&lt;106.5 M ⊙) show significant late-time flattening. Conversely, FUV light curves from high-mass black hole TDFs are generally consistent with an extrapolation from the early-time light curve. The observed late-time emission cannot be explained by existing models for early-time TDF light curves (i.e., reprocessing or circularization shocks), but is instead consistent with a viscously spreading, unobscured accretion disk. These disk models can only reproduce the observed FUV luminosities, however, if they are assumed to be thermally and viscously stable, in contrast to the simplest predictions of α-disk theory. For one TDF in our sample, we measure an upper limit to the UV luminosity that is significantly lower than expectations from theoretical modeling and an extrapolation of the early-time light curve. This dearth of late-time emission could be due to a disk instability/state change absent in the rest of the sample. The disk models that explain the late-time UV detections solve the TDF “missing energy problem” by radiating a rest-mass energy of ∼0.1 M ⊙ over a period of decades, primarily in extreme UV wavelengths.

Aims. The modelling of spectroscopic observations of tidal disruption events (TDEs) to date suggests that the newly formed accretion disks are mostly quasi-circular. In this work we study the transient event AT 2020zso, hosted by an active galactic nucleus (AGN; as inferred from narrow emission line diagnostics), with the aim of characterising the properties of its newly formed accretion flow. Methods. We classify AT 2020zso as a TDE based on the blackbody evolution inferred from UV/optical photometric observations and spectral line content and evolution. We identify transient, double-peaked Bowen (N III), He I, He II, and Hα emission lines. We model medium-resolution optical spectroscopy of the He II (after careful de-blending of the N III contribution) and Hα lines during the rise, peak, and early decline of the light curve using relativistic, elliptical accretion disk models. Results. We find that the spectral evolution before the peak can be explained by optical depth effects consistent with an outflowing, optically thick Eddington envelope. Around the peak, the envelope reaches its maximum extent (approximately 1015 cm, or ∼3000–6000 gravitational radii for an inferred black hole mass of 5−10 × 105 M⊙) and becomes optically thin. The Hα and He II emission lines at and after the peak can be reproduced with a highly inclined (i = 85 ± 5 degrees), highly elliptical (e = 0.97 ± 0.01), and relatively compact (Rin = several 100 Rg and Rout = several 1000 Rg) accretion disk. Conclusions. Overall, the line profiles suggest a highly elliptical geometry for the new accretion flow, consistent with theoretical expectations of newly formed TDE disks. We quantitatively confirm, for the first time, the high inclination nature of a Bowen (and X-ray dim) TDE, consistent with the unification picture of TDEs, where the inclination largely determines the observational appearance. Rapid line profile variations rule out the binary supermassive black hole hypothesis as the origin of the eccentricity; these results thus provide a direct link between a TDE in an AGN and the eccentric accretion disk. We illustrate for the first time how optical spectroscopy can be used to constrain the black hole spin, through (the lack of) disk precession signatures (changes in inferred inclination). We constrain the disk alignment timescale to &gt; 15 days in AT2020zso, which rules out high black hole spin values (a &lt; 0.8) for MBH ∼ 106 M⊙ and disk viscosity α ≳ 0.1.

When a star passes too close to a supermassive black hole, it gets disrupted by strong tidal forces. The stellar debris then evolves into an elongated stream of gas that partly falls back toward the black hole. We present an analytical model describing for the first time the full stream evolution during such a tidal disruption event (TDE). Our framework consists of dividing the stream into different sections of elliptical geometry, whose properties are independently evolved in their comoving frame under the tidal, pressure, and self-gravity forces. Through an explicit treatment of the tidal force and the inclusion of the gas angular momentum, we can accurately follow the stream evolution near pericenter. Our model evolves the longitudinal stream stretching and both transverse widths simultaneously. For the latter, we identify two regimes depending on whether the dynamics is entirely dominated by the tidal force (ballistic regime) or additionally influenced by pressure and self-gravity (hydrostatic regime). We find that the stream undergoes transverse collapses both shortly after the stellar disruption and upon its return near the black hole, at specific locations determined by the regime of evolution considered. The stream evolution predicted by our model can be used to determine the subsequent interactions experienced by this gas that are at the origin of most of the electromagnetic emission from TDEs. Our results suggest that the accretion disk may be fed at a rate that differs from the standard fallback rate, which would provide novel observational signatures dependent on black hole spin.