No Arabic abstract
Tidal disruption events (TDEs) are among the brightest transients in the optical, ultraviolet, and X-ray sky. These flares are set into motion when a star is torn apart by the tidal field of a massive black hole, triggering a chain of events which is -- so far -- incompletely understood. However, the disruption process has been studied extensively for almost half a century, and unlike the later stages of a TDE, our understanding of the disruption itself is reasonably well converged. In this Chapter, we review both analytical and numerical models for stellar tidal disruption. Starting with relatively simple, order-of-magnitude physics, we review models of increasing sophistication, the semi-analytic ``affine formalism, hydrodynamic simulations of the disruption of polytropic stars, and the most recent hydrodynamic results concerning the disruption of realistic stellar models. Our review surveys the immediate aftermath of disruption in both typical and more unusual TDEs, exploring how the fate of the tidal debris changes if one considers non-main sequence stars, deeply penetrating tidal encounters, binary star systems, and sub-parabolic orbits. The stellar tidal disruption process provides the initial conditions needed to model the formation of accretion flows around quiescent massive black holes, and in some cases may also lead to directly observable emission, for example via shock breakout, gravitational waves or runaway nuclear fusion in deeply plunging TDEs.
When a star approaches a black hole closely, it may be pulled apart by gravitational forces in a tidal disruption event (TDE). The flares produced by TDEs are unique tracers of otherwise quiescent supermassive black holes (SMBHs) located at the centre of most galaxies. In particular, the appearance of such flares and the subsequent decay of the light curve are both sensitive to whether the star is partially or totally destroyed by the tidal field. However, the physics of the disruption and the fall-back of the debris are still poorly understood. We are here modelling the hydrodynamical evolution of realistic stars as they approach a SMBH on parabolic orbits, using for the first time the moving-mesh code AREPO, which is particularly well adapted to the problem through its combination of quasi-Lagrangian behaviour, low advection errors, and high accuracy typical of mesh-based techniques. We examine a suite of simulations with different impact parameters, allowing us to determine the critical distance at which the star is totally disrupted, the energy distribution and the fallback rate of the debris, as well as the hydrodynamical evolution of the stellar remnant in the case of a partial disruption. Interestingly, we find that the internal evolution of the remnants core is strongly influenced by persistent vortices excited in the tidal interaction. These should be sites of strong magnetic field amplification, and the associated mixing may profoundly alter the subsequent evolution of the tidally pruned star.
We analyze the early growth stage of direct-collapse black holes (DCBHs) with $sim 10^{5} rm M_odot$, which are formed by collapse of supermassive stars in atomic-cooling halos at $z gtrsim 10$. A nuclear accretion disk around a newborn DCBH is gravitationally unstable and fragments into clumps with a few $10 rm M_odot$ at $sim 0.01-0.1 rm pc$ from the center. Such clumps evolve into massive population III stars with a few $10-100 rm M_odot$ via successive gas accretion and a nuclear star cluster is formed. Radiative and mechanical feedback from an inner slim disk and the star cluster will significantly reduce the gas accretion rate onto the DCBH within $sim 10^6 rm yr$. Some of the nuclear stars can be scattered onto the loss cone orbits also within $lesssim 10^6 rm yr$ and tidally disrupted by the central DCBH. The jet luminosity powered by such tidal disruption events can be $L_{rm j} gtrsim 10^{50} rm erg s^{-1}$. The prompt emission will be observed in X-ray bands with a peak duration of $delta t_{rm obs} sim 10^{5-6} (1+z) rm s$ followed by a tail $propto t_{rm obs}^{-5/3}$, which can be detectable by Swift BAT and eROSITA even from $z sim 20$. Follow-up observations of the radio afterglows with, e.g., eVLA and the host halos with JWST could probe the earliest AGN feedback from DCBHs.
Optical transient surveys have led to the discovery of dozens of stellar tidal disruption events (TDEs) by massive black hole in the centers of galaxies. Despite extensive searches, X-ray follow-up observations have produced no or only weak X-ray detections in most of them. Here we report the discovery of delayed X-ray brightening around 140 days after the optical outburst in the TDE OGLE16aaa, followed by several flux dips during the decay phase. These properties are unusual for standard TDEs and could be explained by the presence of supermassive black hole binary or patchy obscuration. In either scenario, the X-rays can be produced promptly after the disruption but are blocked in the early phase, possibly by a radiation-dominated ejecta which leads to the bulk of optical and ultraviolet emission. Our findings imply that the reprocessing is important in the TDE early evolution, and X-ray observations are promising in revealing supermassive black hole binaries.
Supermassive black hole binaries (SMBHBs) are products of galaxy mergers, and are important in testing Lambda cold dark matter cosmology and locating gravitational-wave-radiation sources. A unique electromagnetic signature of SMBHBs in galactic nuclei is essential in identifying the binaries in observations from the IR band through optical to X-ray. Recently, the flares in optical, UV, and X-ray caused by supermassive black holes (SMBHs) tidally disrupting nearby stars have been successfully used to observationally probe single SMBHs in normal galaxies. In this Letter, we investigate the accretion of the gaseous debris of a tidally disrupted star by a SMBHB. Using both stability analysis of three-body systems and numerical scattering experiments, we show that the accretion of stellar debris gas, which initially decays with time $propto t^{-5/3}$, would stop at a time $T_{rm tr} simeq eta T_{rm b}$. Here, $eta sim0.25$ and $T_{rm b}$ is the orbital period of the SMBHB. After a period of interruption, the accretion recurs discretely at time $T_{rm r} simeq xi T_b$, where $xi sim 1$. Both $eta$ and $xi$ sensitively depend on the orbital parameters of the tidally disrupted star at the tidal radius and the orbit eccentricity of SMBHB. The interrupted accretion of the stellar debris gas gives rise to an interrupted tidal flare, which could be used to identify SMBHBs in non-active galaxies in the upcoming transient surveys.
We present the first simulations of the tidal disruption of stars with realistic structures and compositions by massive black holes (BHs). We build stars in the stellar evolution code MESA and simulate their disruption in the 3D adaptive-mesh hydrodynamics code FLASH, using an extended Helmholtz equation of state and tracking 49 elements. We study the disruption of a 1$M_odot$ star and 3$M_odot$ star at zero-age main sequence (ZAMS), middle-age, and terminal-age main sequence (TAMS). The maximum BH mass for tidal disruption increases by a factor of $sim$2 from stellar radius changes due to MS evolution; this is equivalent to varying BH spin from 0 to 0.75. The shape of the mass fallback rate curves is different from the results for polytropes of Guillochon & Ramirez-Ruiz (2013). The peak timescale $t_{rm peak}$ increases with stellar age, while the peak fallback rate $dot M_{rm peak}$ decreases with age, and these effects diminish with increasing impact parameter $beta$. For a $beta=1$ disruption of a 1$M_odot$ star by a $10^6 M_odot$ BH, from ZAMS to TAMS, $t_{rm peak}$ increases from 30 to 54 days, while $dot M_{rm peak}$ decreases from 0.66 to 0.14 $M_odot$/yr. Compositional anomalies in nitrogen, helium, and carbon can occur before the peak timescale for disruptions of MS stars, which is in contrast to predictions from the frozen-in model. More massive stars can show stronger anomalies at earlier times, meaning that compositional constraints can be key in determining the mass of the disrupted star. The abundance anomalies predicted by these simulations provide a natural explanation for the spectral features and varying line strengths observed in tidal disruption events.