The mass fallback rate of the debris in relativistic stellar tidal disruption events

TL;DR

This study uses simulations to analyze the mass fallback rate in relativistic and Newtonian tidal disruption events, revealing how black hole spin, stellar properties, and gravity description influence debris accretion dynamics.

Contribution

It provides a comprehensive simulation-based analysis of TDE debris fallback rates considering relativistic effects, stellar profiles, and black hole spin orientations, extending previous models.

Findings

Relativistic effects lead to higher peak fallback rates and shorter peak times.

Fallback duration scales primarily with stellar mass and compactness.

Black hole spin orientation significantly influences debris fallback curves.

Abstract

Highly energetic stellar tidal disruption events (TDEs) provide a way to study black hole characteristics and their environment. We simulate TDEs with the PHANTOM code in a general relativistic and Newtonian description of a supermassive black hole's gravity. Stars, which are placed on parabolic orbits with different parameters $\beta$, are constructed with the stellar evolution code MESA and therefore have realistic stellar density profiles. We study the mass fallback rate of the debris $\dot{M}$ and its dependence on the $\beta$, stellar mass and age as well as the black hole's spin and the choice of the gravity's description. We calculate peak value $\dot{M}_\mathrm{peak}$, time to the peak $t_\mathrm{peak}$, duration of the super-Eddington phase $t_\mathrm{Edd}$, time $t_{>0.5\dot{M}_\mathrm{peak}}$ during which $\dot{M}>0.5\dot{M}_\mathrm{peak}$, early rise-time $\tau_\mathrm{rise}$ and late-time slope $n_\infty$. We recover the trends of $\dot{M}_\mathrm{peak}$, $t_\mathrm{peak}$, $\tau_\mathrm{rise}$ and $n_\infty$ with $\beta$, stellar mass and age, which were obtained in previous studies. We find that $t_\mathrm{Edd}$, at a fixed $\beta$, scales primarily with the stellar mass, while $t_{>0.5\dot{M}_\mathrm{peak}}$ scales with the compactness of stars. The effect of SMBH's rotation depends on the orientation of its rotational axis relative to the direction of the stellar motion on the initial orbit. Encounters on prograde orbits result in narrower $\dot{M}$ curves with higher $\dot{M}_\mathrm{peak}$, while the opposite occurs for retrograde orbits. We find that disruptions, at the same pericenter distance, are stronger in a relativistic tidal field than in a Newtonian. Therefore, relativistic $\dot{M}$ curves have higher $\dot{M}_\mathrm{peak}$, and shorter $t_\mathrm{peak}$ and $t_\mathrm{Edd}$.