Stars Crushed by Black Holes. I. On the Energy Distribution of Stellar Debris in Tidal Disruption Events
Abstract The distribution of orbital energies imparted into stellar debris following the close encounter of a star with a supermassive black hole is the principal factor in determining the rate of return of debris to the black hole, and thus in determining the properties of the resulting lightcurves from such events. We present simulations of tidal disruption events for a range of β ≡ r t / r p where r p is the pericenter distance and r t the tidal radius. We perform these simulations at different spatial resolutions to determine the numerical convergence of our models. We compare simulations in which the heating due to shocks is included or excluded from the dynamics. For β ≲ 8, the simulation results are well-converged at sufficiently moderate-to-high spatial resolution, while for β ≳ 8, the breadth of the energy distribution can be grossly exaggerated by insufficient spatial resolution. We find that shock heating plays a non-negligible role only for β ≳ 4, and that typically the effect of shock heating is mild. We show that self-gravity can modify the energy distribution over time after the debris has receded to large distances for all β . Primarily, our results show that more »
Authors:
; ;
Award ID(s):
Publication Date:
NSF-PAR ID:
10350415
Journal Name:
The Astrophysical Journal
Volume:
923
Issue:
2
Page Range or eLocation-ID:
184
ISSN:
0004-637X
National Science Foundation
##### More Like this
1. Abstract Upon entering the tidal sphere of a supermassive black hole, a star is ripped apart by tides and transformed into a stream of debris. The ultimate fate of that debris, and the properties of the bright flare that is produced and observed, depends on a number of parameters, including the energy of the center of mass of the original star. Here we present the results of a set of smoothed particle hydrodynamics simulations in which a 1 M ⊙ , γ = 5/3 polytrope is disrupted by a 10 6 M ⊙ supermassive black hole. Each simulation has a pericenter distance of r p = r t (i.e., β ≡ r t / r p = 1 with r t the tidal radius), and we vary the eccentricity e of the stellar orbit from e = 0.8 up to e = 1.20 and study the nature of the fallback of debris onto the black hole and the long-term fate of the unbound material. For simulations with eccentricities e ≲ 0.98, the fallback curve has a distinct, three-peak structure that is induced by self-gravity. For simulations with eccentricities e ≳ 1.06, the core of the disrupted star reforms following itsmore »
2. Abstract

We present a toy model for the thermal optical/UV/X-ray emission from tidal disruption events (TDEs). Motivated by recent hydrodynamical simulations, we assume that the debris streams promptly and rapidly circularize (on the orbital period of the most tightly bound debris), generating a hot quasi-spherical pressure-supported envelope of radiusRv∼ 1014cm (photosphere radius ∼1015cm) surrounding the supermassive black hole (SMBH). As the envelope cools radiatively, it undergoes Kelvin–Helmholtz contractionRvt−1, its temperature risingTefft1/2while its total luminosity remains roughly constant; the optical luminosity decays as$νLν∝Rv2Teff∝t−3/2$. Despite this similarity to the mass fallback rate$Ṁfb∝t−5/3$, envelope heating from fallback accretion is subdominant compared to the envelope cooling luminosity except near optical peak (where they are comparable). Envelope contraction can be delayed by energy injection from accretion from the inner envelope onto the SMBH in a regulated manner, leading to a late-time flattening of the optical/X-ray light curves, similar to those observed in some TDEs. Eventually, as the envelope contracts to near the circularization radius, the SMBH accretion rate rises to its maximum, in tandem with the decreasing optical luminosity. This cooling-induced (rather than circularization-induced) delay of up to several hundred days may account for themore »

3. Abstract

We develop a Newtonian model of a deep tidal disruption event (TDE), for which the pericenter distance of the star,rp, is well within the tidal radius of the black hole,rt, i.e., whenβrt/rp≫ 1. We find that shocks form forβ≳ 3, but they are weak (with Mach numbers ∼1) for allβ, and that they reach the center of the star prior to the time of maximum adiabatic compression forβ≳ 10. The maximum density and temperature reached during the TDE follow much shallower relations withβthan the previously predicted$ρmax∝β3$and$Tmax∝β2$scalings. Belowβ≃ 10, this shallower dependence occurs because the pressure gradient is dynamically significant before the pressure is comparable to the ram pressure of the free-falling gas, while aboveβ≃ 10, we find that shocks prematurely halt the compression and yield the scalings$ρmax∝β1.62$and$Tmax∝β1.12$. We find excellent agreement between our results and high-resolution simulations. Our results demonstrate that, in the Newtonian limit, the compression experienced by the star is completely independent of the mass of the black hole. We discuss our results in the context of existing (affine) models, polytropic versus non-polytropic stars, and general relativistic effects, which become important when the pericenter ofmore »

4. ABSTRACT

Tidal disruption events (TDEs) occur when a star gets torn apart by the strong tidal forces of a supermassive black hole, which results in the formation of a debris stream that partly falls back towards the compact object. This gas moves along inclined orbital planes that intersect near pericentre, resulting in a so-called ‘nozzle shock’. We perform the first dedicated study of this interaction, making use of a two-dimensional simulation that follows the transverse gas evolution inside a given section of stream. This numerical approach circumvents the lack of resolution encountered near pericentre passage in global three-dimensional simulations using particle-based methods. As it moves inward, we find that the gas motion is purely ballistic, which near pericentre causes strong vertical compression that squeezes the stream into a thin sheet. Dissipation takes place at the resulting nozzle shock, inducing a rise in pressure that causes the collapsing gas to bounce back, although without imparting significant net expansion. As it recedes to larger distances, this matter continues to expand while remaining thin despite the influence of pressure forces. This gas evolution specifies the strength of the subsequent self-crossing shock, which we find to be more affected by black hole spin than previouslymore »

5. Abstract Tidal disruption events with tidal radius r t and pericenter distance r p are characterized by the quantity β = r t / r p , and “deep encounters” have β ≫ 1. It has been assumed that there is a critical β ≡ β c ∼ 1 that differentiates between partial and full disruption: for β < β c a fraction of the star survives the tidal interaction with the black hole, while for β > β c the star is completely destroyed, and hence all deep encounters should be full. Here we show that this assumption is incorrect by providing an example of a β = 16 encounter between a γ = 5/3, solar-like polytrope and a 10 6 M ⊙ black hole—for which previous investigations have found β c ≃ 0.9—that results in the reformation of a stellar core post-disruption that comprises approximately 25% of the original stellar mass. We propose that the core reforms under self-gravity, which remains important because of the compression of the gas both near pericenter, where the compression occurs out of the orbital plane, and substantially after pericenter, where compression is within the plane. We find that the core forms onmore »