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.

We report the discovery of the first unambiguous case of resolved double-peaked Balmer emission in a tidal disruption event (TDE), AT~2018hyz, that serves as a clear indication of an accretion disk. We modelled the \Ha\ emission line with a low eccentricity ($e\approx0.1$) accretion disk extending out to $\sim$100 $R_{p}$ and a Gaussian component originating from non-disk clouds. Our analysis indicates that in AT~2018hyz, disk formation took place promptly after the most-bound debris returned to pericenter, which we estimate to be roughly tens of days before the first detection. Redistribution of angular momentum and mass transport, possibly through shocks, must occur on the observed timescale of about a month to create the large \Ha-emitting disk that comprises $\lesssim$10\% of the initial stellar mass. With these new insights from AT~2018hyz, we infer that efficient circularization may be common among optically-detected TDEs despite their lack of double-peaked emission, which can be largely affected by the disk inclination angle and the relative strength of the disk contribution to the non-disk component.

Recent studies of Tidal Disruption Events (TDEs) have revealed unexpected correlations between the TDE rate and the large-scale properties of the host galaxies. In this review, we present the host galaxy properties of all TDE candidates known to date and quantify their distributions. We consider throughout the differences between observationally-identified types of TDEs and differences from spectroscopic control samples of galaxies. We focus here on the black hole and stellar masses of TDE host galaxies, their star formation histories and stellar populations, the concentration and morphology of the optical light, the presence of AGN activity, and the extra-galactic environment of the TDE hosts. We summarize the state of several possible explanations for the links between the TDE rate and host galaxy type. We present estimates of the TDE rate for different host galaxy types and quantify the degree to which rate enhancement in some types results in rate suppression in others. We discuss the possibilities for using TDE host galaxies to assist in identifying TDEs in upcoming large transient surveys and possibilities for TDE observations to be used to study their host galaxies.

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.

We use a simple framework to calculate the time evolution of the composition of the fallback material onto a supermassive black hole arising from the tidal disruption of main sequence stars. We study stars with masses between 0.8 and 3.0 $M_\odot$, at evolutionary stages from zero-age main sequence to terminal-age main sequence, built using the Modules for Experiments in Stellar Astrophysics code. We show that most stars develop enhancements in nitrogen ($^{14}$N) and depletions in carbon ($^{12}$C) and oxygen ($^{16}$O) over their lifetimes, and that these features are more pronounced for higher mass stars. We find that, in an accretion-powered tidal disruption flare, these features become prominent only after the time of peak of the fallback rate and appear at earlier times for stars of increasing mass. We postulate that no severe compositional changes resulting from the fallback material should be expected near peak for a wide range of stellar masses and, as such, are unable to explain the extreme helium-to-hydrogen line ratios observed in some TDEs. On the other hand, the resulting compositional changes could help explain the presence of nitrogen-rich features, which are currently only detected after peak. When combined with the shape of the light curve, the time evolution of the composition of the fallback material provides a clear method to help constrain the nature of the disrupted star. This will enable a better characterization of the event by helping break the degeneracy between the mass of the star and the mass of the black hole when fitting tidal disruption light curves.

We study the properties of tidal disruption event (TDE) host galaxies in the context of a catalog of ~500,000 galaxies from the Sloan Digital Sky Survey. We explore whether selection effects can account for the overrepresentation of TDEs in E+A/post-starburst galaxies by creating matched galaxy samples. Accounting for possible selection effects due to black hole (BH) mass, redshift completeness, strong AGN presence, bulge colors, and surface brightness can reduce the apparent overrepresentation of TDEs in E+A host galaxies by a factor of ~4 (from ~$\times$100-190 to ~$\times$25-48), but cannot fully explain the preference. We find that TDE host galaxies have atypical photometric properties compared to similar, "typical" galaxies. In particular, TDE host galaxies tend to live in or near the "green valley" between star-forming and passive galaxies, and have bluer bulge colors ($\Delta (g-r) \approx 0.3$ mag), lower half-light surface brightnesses (by ~1 mag/arcsec$^2$), higher Sersic indices ($\Delta n_{\rm g} \approx 3$), and higher bulge-to-total-light ratios ($\Delta B/T \approx 0.5$) than galaxies with matched BH masses. We find that TDE host galaxies appear more centrally concentrated and that all have high galaxy Sersic indices and $B/T$ fractions---on average in the top 10% of galaxies of the same BH mass---suggesting a higher nuclear stellar density. We identify a region in Sersic index and BH mass parameter space that contains ~2% of our reference catalog galaxies but $\ge\!60\%$ of TDE host galaxies. The unique photometric properties of TDE host galaxies may be useful for selecting candidate TDEs for spectroscopic follow-up observations in large transient surveys.

We construct a menu of objects that can give rise to bright flares when disrupted by massive black holes (BHs), ranging from planets to evolved stars. Through their tidal disruption, main sequence and evolved stars can effectively probe the existence of otherwise quiescent supermassive BHs and white dwarfs can probe intermediate mass BHs. Many low-mass white dwarfs possess extended hydrogen envelopes, which allow for the production of prompt flares in disruptive encounters with moderately massive BHs of $10^5$ to $10^7~M_\odot$--masses that may constitute the majority of massive BHs by number. These objects are a missing link in two ways: (1) for probing moderately massive BHs and (2) for understanding the hydrodynamics of the disruption of objects with tenuous envelopes. A flare arising from the tidal disruption of a $0.17~M_\odot$ white dwarf by a $10^5~M_\odot$ BH reaches a maximum between 0.6 and 11 days, with a peak fallback rate that is usually super-Eddington and results in a flare that is likely brighter than a typical tidal disruption event. Encounters stripping only the envelope can provide hydrogen-only fallback, while encounters disrupting the core evolve from H- to He-rich fallback. While most tidal disruption candidates observed thus far are consistent with the disruptions of main sequence stars, the rapid timescales of nuclear transients such as Dougie and PTF10iya are naturally explained by the disruption of low-mass white dwarfs. As the number of observed flares continues to increase, the menu presented here will be essential for characterizing nuclear BHs and their environments through tidal disruptions.

We measure the color and stellar mass dependence of clustering in spectroscopic galaxies at $0.6 < z < 0.65$ using data from the Baryon Oscillation Spectroscopic Survey component of the Sloan Digital Sky Survey. We greatly increase the statistical precision of our clustering measurements by using the cross-correlation of 66,657 spectroscopic galaxies to a sample of 6.6 million fainter photometric galaxies. The clustering amplitude $w(R)$ is measured as the ratio of the mean excess number of photometric galaxies found within a specified radius annulus around a spectroscopic galaxy to that from a random photometric galaxy distribution. We recover many of the familiar trends at high signal-to-noise ratio. We find the ratio of the clustering amplitudes of red and blue massive galaxies to be $w_\text{red}/w_\text{blue} = 1.92 \pm 0.11$ in our smallest annulus of 75-125 kpc. At our largest radii (2-4 Mpc), we find $w_\text{red}/w_\text{blue} = 1.24 \pm 0.05$. Red galaxies therefore have denser environments than their blue counterparts at $z \sim 0.625$, and this effect increases with decreasing radius. Irrespective of color, we find that $w(R)$ does not obey a simple power-law relation with radius, showing a dip around 1 Mpc. Holding stellar mass fixed, we find a clear differentiation between clustering in red and blue galaxies, showing that clustering is not solely determined by stellar mass. Holding color fixed, we find that clustering increases with stellar mass, especially for red galaxies at small scales (more than a factor of 2 effect over 0.75 dex in stellar mass).