We develop a new model for X-ray emission from tidal disruption events (TDEs), applying stationary general relativistic ``slim disk'' accretion solutions to supermassive black holes (SMBHs) and then ray-tracing the photon trajectories from the image plane to the disk surface, including gravitational redshift, Doppler, and lensing effects self-consistently. We simultaneously and successfully fit the multi-epoch XMM-Newton X-ray spectra for two TDEs: ASASSN-14li and ASASSN-15oi. We test explanations for the observed, unexpectedly slow X-ray brightening of ASASSN-15oi, including delayed disk formation and variable obscuration by a reprocessing layer. We propose a new mechanism that better fits the data: a ``Slimming Disk'' scenario in which accretion onto an edge-on disk slows, reducing the disk height and exposing more X-rays from the inner disk to the sightline over time.For ASASSN-15oi, we constrain the SMBH mass to $4.0^{+2.5}_{-3.1} \times 10^6M_\odot$. For ASASSN-14li, the SMBH mass is $10^{+1}_{-7}\times 10^6M_\odot$ and the spin is $>0.3$. For both TDEs, our fitted masses are consistent with independent estimates; for ASASSN-14li, application of the external mass constraint narrows our spin constraint to $>0.85$. The mass accretion rate of ASASSN-14li decays slowly, as $\propto t^{-1.1}$, perhaps due to inefficient debris circularization. Over $\approx$1100 days, its SMBH has accreted $\Delta M \approx 0.17 M_\odot$, implying a progenitor star mass of $> 0.34 M_\odot$, i.e., no ``missing energy problem.'' For both TDEs, the hydrogen column density declines to the host galaxy plus Milky Way value after a few hundred days, suggesting a characteristic timescale for the depletion or removal of obscuring gas.

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.

Tidal disruption events occur rarely in any individual galaxy. Over the last decade, however, time-domain surveys have begun to accumulate statistical samples of these flares. What dynamical processes are responsible for feeding stars to supermassive black holes? At what rate are stars tidally disrupted in realistic galactic nuclei? What may we learn about supermassive black holes and broader astrophysical questions by estimating tidal disruption event rates from observational samples of flares? These are the questions we aim to address in this Chapter, which summarizes current theoretical knowledge about rates of stellar tidal disruption, and compares theoretical predictions to the current state of observations.

We explore the galaxy structure of four tidal disruption event (TDE) host galaxies on 30 pc to kpc scales using HST WFC3 multi-band imaging. The star formation histories of these hosts are diverse, including one post-starburst galaxy (ASASSN-14li), two hosts with recent weak starbursts (ASASSN-14ae and iPTF15af), and one early type (PTF09ge). Compared to early type galaxies of similar stellar masses, the TDE hosts have higher central surface brightnesses and stellar mass surface densities on 30-100 pc scales. The TDE hosts do not show the large, kpc-scale tidal disruptions seen in some post-starburst galaxies; the hosts have low morphological asymmetries similar to those of early type galaxies. The lack of strong asymmetries are inconsistent with a recent major (~1:1 mass) merger, although minor ($\lesssim$1:3) mergers are possible. Given the time elapsed since the end of the starbursts in the three post-burst TDE hosts and the constraints on the merger mass ratios, it is unlikely that a bound supermassive black hole binary (SMBHB) has had time to coalesce. The TDE hosts have low central (<140 pc) ellipticities compared to early type galaxies. The low central ellipticities disfavor a strong radial anisotropy as the cause for the enhanced TDE rate, although we cannot rule out eccentric disks at the scale of the black hole gravitational radius of influence (~1 pc). These observations suggest that the high central stellar densities are a more important driver than SMBHBs or radial anisotropies in increasing the TDE rate in galaxies with recent starbursts.

We present Chandra X-ray observations of four optically-selected tidal disruption events (TDEs) obtained 4-9 years after discovery. Three sources were detected with luminosities between 9X10^40 and 3X10^42 erg/s. The spectrum of PTF09axc is consistent with a power law of index 2.5+-0.1, whereas the spectrum of PTF09ge is very soft. The power law spectrum of PTF09axc and prior literature findings, provide evidence that TDEs transition from an early-time soft state to a late-time hard state many years after disruption. We propose that the time to peak luminosity for optical and X-ray emission may differ substantially in TDEs, with X-rays being produced or becoming observable later. This delay helps explain the differences in observed properties such as L_opt/L_ X of optically and X-ray selected TDEs. We update TDE rate predictions for the eROSITA instrument: it ranges from 3 per yr to 990 per yr, depending sensitively on the distribution of black hole spins and the time delay between disruption and peak X-ray brightness. We further predict an asymmetry in the number of retrograde and prograde disks in samples of optically and X-ray selected TDEs. The details of the observational biases can contribute to observed differences between optically and X-ray selected TDEs (with optically selected TDEs being fainter in X-rays for retrograde TDE disks).

The three-body problem is arguably the oldest open question in astrophysics, and has resisted a general analytic solution for centuries. Various implementations of perturbation theory provide solutions in portions of parameter space, but only where hierarchies of masses or separations exist. Numerical integrations show that bound, non-hierarchical triples of Newtonian point particles will almost always disintegrate into a single escaping star and a stable, bound binary, but the chaotic nature of the three-body problem prevents the derivation of tractable analytic formulae deterministically mapping initial conditions to final outcomes. However, chaos also motivates the assumption of ergodicity, suggesting that the distribution of outcomes is uniform across the accessible phase volume. Here, we use the ergodic hypothesis to derive a complete statistical solution to the non-hierarchical three-body problem, one which provides closed-form distributions of outcomes (e.g. binary orbital elements) given the conserved integrals of motion. We compare our outcome distributions to large ensembles of numerical three-body integrations, and find good agreement, so long as we restrict ourselves to "resonant" encounters (the ~50% of scatterings that undergo chaotic evolution). In analyzing our scattering experiments, we identify "scrambles" (periods in time where no pairwise binaries exist) as the key dynamical state that ergodicizes a non-hierarchical triple. The generally super-thermal distributions of survivor binary eccentricity that we predict have notable applications to many astrophysical scenarios. For example, non-hierarchical triples produced dynamically in globular clusters are a primary formation channel for black hole mergers, but the rates and properties of the resulting gravitational waves depend on the distribution of post-disintegration eccentricities.

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 years 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 ($<10^{6.5} M_\odot$) 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 alpha-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 solar mass over a period of decades, primarily in extreme UV wavelengths.

We present new medium resolution, optical long-slit spectra of a sample of 6 UV/optical and 15 X-ray selected tidal disruption event candidate host galaxies. We measure emission line ratios from the optical spectra, finding that the large majority of hosts are quiescent galaxies, while those displaying emission lines are generally consistent with star-formation dominated environments; only 3 sources show clear evidence of nuclear activity. We measure bulge velocity dispersions using absorption lines and infer host black hole (BH) masses using the M-$\sigma$ relation. While the optical and X-ray host BH masses are statistically consistent with coming from the same parent distribution, the optical host distribution has a visible peak near $M_{\rm BH} \sim 10^6 M_\odot$, whereas the X-ray host distribution appears flat in $M_{\rm BH}$. We find a subset of X-ray selected candidates that are hosted in galaxies significantly less luminous (M$_{\rm g}$$\sim$-16) and less massive (stellar mass$\sim$10$^{8.5-9}$M$_{\odot}$) than those of optical events. Using statistical tests we find suggestive evidence that, in terms of black hole mass, stellar mass and absolute magnitude, the hard X-ray hosts differ from the UV/optical and soft X-ray samples. Similar to individual studies, we find that the emission region size for the X-ray sample is much smaller than the optical emission region, consistent with a compact accretion disk. We find a typical Eddington ratio of the X-ray emission is $\sim$0.01, as opposed to the optical events which have L$_{\rm BB}$$\sim$L$_{\rm Edd}$.

Gravitational perturbations on an exoplanet from a massive outer body, such as the Kozai-Lidov mechanism, can pump the exoplanet's eccentricity up to values that will destroy it via a collision or strong interaction with its parent star. During the final stages of this process, any exomoons orbiting the exoplanet will be detached by the star's tidal force and placed into orbit around the star. Using ensembles of three and four-body simulations, we demonstrate that while most of these detached bodies either collide with their star or are ejected from the system, a substantial fraction, ~10%, of such "orphaned" exomoons (with initial properties similar to those of the Galilean satellites in our own solar system) will outlive their parent exoplanet. The detached exomoons generally orbit inside the ice line, so that strong radiative heating will evaporate any volatile-rich layers, producing a strong outgassing of gas and dust, analogous to a comet's perihelion passage. Small dust grains ejected from the exomoon may help generate an opaque cloud surrounding the orbiting body but are quickly removed by radiation blow-out. By contrast, larger solid particles inherit the orbital properties of the parent exomoon, feeding an eccentric disk of solids that drains more gradually onto the star via Poynting-Robertson drag, and which could result in longer-timescale dimming of the star. For characteristic exomoon evaporation times of ~ 1e5-1e6 yr, attenuation of the stellar light arising from one or more out-gassing exomoons provides a promising explanation for both the dipping and secular dimming behavior observed from KIC 8462852 (Boyajian's Star).

The growing number of stellar-mass binary black hole mergers discovered by Advanced LIGO and Advanced Virgo are starting to constrain the binaries' origin and environment. However, we still lack sufficiently accurate modeling of binary formation channels to obtain strong constraints, or to identify sub-populations. One promising formation mechanism that could result in different black hole properties is binaries merging within the accretion disks of Active Galactic Nuclei (AGN). Here we show that the black holes' orbital alignment with the AGN disks preferentially selects heavier black holes. We carry out Monte Carlo simulations of orbital alignment with AGN disks, and find that AGNs harden the initial black hole mass function. Assuming an initial power law mass distribution $M_{\rm bh}^{-\beta}$, we find that the power law index changes by $\Delta \beta\sim1.3$, resulting in a more top-heavy population of merging black holes. This change is independent of the mass of, and accretion rate onto, the supermassive black hole in the center of the AGN. Our simulations predict an AGN-assisted merger rate of $\sim4$Gpc$^{-3}$yr$^{-1}$. With its hardened mass spectra, the AGN channel could be responsible for $10-50$% of gravitational-wave detections.

The tidal destruction of a star by a massive black hole, known as a tidal disruption event (TDE), is commonly modeled using the "frozen-in" approximation. Under this approximation, the star maintains exact hydrostatic balance prior to entering the tidal sphere (radius $r_{\rm t}$), after which point its internal pressure and self-gravity become instantaneously negligible and the debris undergoes ballistic free fall. We present a suite of hydrodynamical simulations of TDEs with high penetration factors $\beta \equiv r_{\rm t}/r_{\rm p} = 5-7$, where $r_{\rm p}$ is the pericenter of the stellar center of mass, calculated using a Voronoi-based moving-mesh technique. We show that basic assumptions of the frozen-in model, such as the neglect of self-gravity inside $r_{\rm t}$, are violated. Indeed, roughly equal fractions of the final energy spread accumulate exiting and entering the tidal sphere, though the frozen-in prediction is correct at the order-of-magnitude level. We also show that an $\mathcal{O}(1)$ fraction of the debris mass remains transversely confined by self-gravity even for large $\beta$ which has implications for the radio emission from the unbound debris and, potentially, for the circularization efficiency of the bound streams.

LISA will open the mHz band of gravitational waves (GWs) to the astronomy community. The strong gravity which powers the variety of GW sources in this band is also crucial in a number of important astrophysical processes at the current frontiers of astronomy. These range from the beginning of structure formation in the early universe, through the origin and cosmic evolution of massive black holes in concert with their galactic environments, to the evolution of stellar remnant binaries in the Milky Way and in nearby galaxies. These processes and their associated populations also drive current and future observations across the electromagnetic (EM) spectrum. We review opportunities for science breakthroughs, involving either direct coincident EM+GW observations, or indirect multimessenger studies. We argue that for the US community to fully capitalize on the opportunities from the LISA mission, the US efforts should be accompanied by a coordinated and sustained program of multi-disciplinary science investment, following the GW data through to its impact on broad areas of astrophysics. Support for LISA-related multimessenger observers and theorists should be sized appropriately for a flagship observatory and may be coordinated through a dedicated mHz GW research center.

The question of how supermassive black holes (SMBHs) grow over cosmic time is a major puzzle in high-energy astrophysics. One promising approach to this problem is via the study of tidal disruption flares (TDFs). These are transient events resulting from the disruption of stars by quiescent supermassive black holes at centers of galaxies. A meter-class X-ray observatory with a time resolution $\sim$ a millisecond and a spectral resolution of a few eV at KeV energies would be revolutionary as it will facilitate high signal to noise spectral-timing studies of several cosmological TDFs. It would open a new era of astrophysics where SMBHs in TDFs at cosmic distances can be studied in similar detail as current studies of much nearer, stellar-mass black hole binaries. Using Athena X-ray observatory as an example, we highlight two specific aspects of spectral-timing analysis of TDFs. (1) Detection of X-ray quasi-periodic oscillations (QPOs) over a redshift range and using these signal frequencies to constrain the spin evolution of SMBHs, and (2) Time-resolved spectroscopy of outflows/winds to probe super-Eddington accretion. SMBH spin distributions at various redshifts will directly allow us to constrain their primary mode of growth as higher spins are predicted due to spin-up for prolonged accretion-mode growth, while lower spins are expected for growth via mergers due to angular momentum being deposited from random directions. A meter-class X-ray telescope will also be able to characterize relativistic TDFs, viz., SwJ1644+57-like events, out to a redshift greater than 8, i.e., it would facilitate detailed spectral-timing studies of TDFs by the youngest SMBHs in the Universe.

We present an analytic description of tides raised on a star by a small orbiting body. In particular, we highlight the disproportionate effect of eccentricity and thus the scope for using these tides to detect and characterise the orbits of exoplanets and brown dwarfs. The tidal distortions of the star produced by an eccentric orbit are, in comparison to a circular orbit, much richer in detail, and potentially visible from any viewing angle. The magnitude of these variations is much larger than that in a circular orbit of the same semi-major axis. These variations are visible in both photometric and spectroscopic data, and dominate other regular sources of phase variability (e.g reflection and Doppler beaming) over a particularly interesting portion of parameter space. These tidal signatures will be a useful tool for planet detection on their own, and used in concert with other methods provide powerful constraints on planetary and stellar properties.

The immense tidal forces of massive black holes can rip apart stars that come too close to them. As the resulting stellar debris spirals inwards, it heats up and emits x-rays when near the black hole. Here, we report the discovery of an exceptionally stable 131-second x-ray quasi-periodicity from a black hole after it disrupted a star. Using a black hole mass indicated from host galaxy scaling relations implies that, (1) this periodicity originates from very close to the black hole's event horizon, and (2) the black hole is rapidly spinning. Our findings suggest that other disruption events with similar highly sensitive observations likely also exhibit quasi-periodicities that encode information about the fundamental properties of their black holes.

A large population of X-ray binaries (XRBs) was recently discovered within the central parsec of the Galaxy by Hailey et al. While the presence of compact objects on this scale due to radial mass segregation is, in itself, unsurprising, the fraction of binaries would naively be expected to be small because of how easily primordial binaries are dissociated in the dynamically hot environment of the nuclear star cluster (NSC). We propose that the formation of XRBs in the central parsec is dominated by the tidal capture of stars by black holes (BHs) and neutron stars (NSs). We model the time-dependent radial density profiles of stars and compact objects in the NSC with a Fokker-Planck approach, using the present-day stellar population and rate of in situ massive star (and thus compact object) formation as observational constraints. Of the ~10,000-40,000 BHs that accumulate in the central parsec over the age of the Galaxy, we predict that ~60 - 200 currently exist as BH-XRBs formed from tidal capture, consistent with the population seen by Hailey et al. A somewhat lower number of tidal capture NS-XRBs is also predicted. We also use our observationally calibrated models for the NSC to predict rates of other exotic dynamical processes, such as the tidal disruption of stars by the central supermassive black hole (~0.0001 per year at z=0).

We develop an analytic, steady-state model for the gas environment in quiescent galactic nuclei. We assume that the mass is constantly supplied by a spherically symmetric distribution of wind emitting stars, and that gravity is solely due to a central supermassive black hole. We show that at some finite radius, where the Keplerian velocity is comparable to the wind velocity, the bulk velocity vanishes. Matter generated below that radius will be accreted onto the black hole, while matter outside it will escape the system. Under certain conditions, the flow may become supersonic at both domains. We obtain radial profiles of the hydrodynamic variables and verify them using a time-dependent hydrodynamic simulation. We delineate the conditions under which radiative cooling can be neglected, and predict the luminosity and spectrum of the free-free X-ray emission from such a system. We discuss applications of our solution to our own Galactic Centre and other quiescent galactic nuclei.

A tidal disruption event (TDE) ensues when a star passes too close to the supermassive black hole (SMBH) in a galactic center and is ripped apart by the tidal field of the SMBH. The gaseous debris produced in a TDE can power a bright electromagnetic flare as it is accreted by the SMBH; so far, several dozen TDE candidates have been observed. For SMBHs with masses above $\sim 10^7 M_\odot$, the tidal disruption of solar-type stars occurs within ten gravitational radii of the SMBH, implying that general relativity (GR) is needed to describe gravity. Three promising signatures of GR in TDEs are: (1) a super-exponential cutoff in the volumetric TDE rate for SMBH masses above $\sim 10^8 M_\odot$ due to direct capture of tidal debris by the event horizon, (2) delays in accretion disk formation (and a consequent alteration of the early-time light curve) caused by the effects of relativistic precession on stream circularization, and (3) quasi-periodic modulation of X-ray emission due to global precession of misaligned accretion disks and the jets they launch. We review theoretical models and simulations of TDEs in Newtonian gravity, then describe how relativistic modifications give rise to these proposed observational signatures, as well as more speculative effects of GR. We conclude with a brief summary of TDE observations and the extent to which they show indications of these predicted relativistic signatures.

This paper reports circular spectropolarimetry and X-ray observations of several polluted white dwarfs including WD 1145+017, with the aim to constrain the behavior of disk material and instantaneous accretion rates in these evolved planetary systems. Two stars with previously observed Zeeman splitting, WD 0322-019 and WD 2105-820, are detected above 5 sigma and <Bz> > 1 kG, while WD 1145+017, WD 1929+011, and WD 2326+049 yield (null) detections below this minimum level of confidence. For these latter three stars, high-resolution spectra and atmospheric modeling are used to obtain limits on magnetic field strengths via the absence of Zeeman splitting, finding B* < 20 kG based on data with resolving power R near 40 000. An analytical framework is presented for bulk Earth composition material falling onto the magnetic polar regions of white dwarfs, where X-rays and cyclotron radiation may contribute to accretion luminosity. This analysis is applied to X-ray data for WD 1145+017, WD 1729+371, and WD 2326+049, and the upper bound count rates are modeled with spectra for a range of plasma kT = 1 - 10 keV in both the magnetic and non-magnetic accretion regimes. The results for all three stars are consistent with a typical dusty white dwarf in a steady-state at 1e8 - 1e9 g/s. In particular, the non-magnetic limits for WD 1145+017 are found to be well below previous estimates of up to 1e12 g/s, and likely below 1e10 g/s, thus suggesting the star-disk system may be average in its evolutionary state, and only special in viewing geometry.

The progenitors of astronomical transients are linked to a specific stellar population and galactic environment, and observing their host galaxies hence constrains the physical nature of the transient itself. Here, we use imaging from the Hubble Space Telescope, and spatially-resolved, medium resolution spectroscopy from the Very Large Telescope obtained with X-Shooter and MUSE to study the host of the very luminous transient ASASSN-15lh. The dominant stellar population at the transient site is old (around 1 to 2 Gyr), without signs of recent star-formation. We also detect emission from ionized gas, originating from three different, time-invariable, narrow components of collisionally-excited metal and Balmer lines. The ratios of emission lines in the Baldwin-Phillips-Terlevich diagnostic diagram indicate that the ionization source is a weak Active Galactic Nucleus with a black hole mass of $M_\bullet = 5_{-3}^{+8}\cdot10^{8} M_\odot$, derived through the $M_\bullet$-$\sigma$ relation. The narrow line components show spatial and velocity offsets on scales of 1 kpc and 500 km/s, respectively; these offsets are best explained by gas kinematics in the narrow-line region. The location of the central component, which we argue is also the position of the supermassive black hole, aligns with that of the transient within an uncertainty of 170 pc. Using this positional coincidence as well as other similarities with the hosts of Tidal Disruption Events, we strengthen the argument that the transient emission observed as ASASSN-15lh is related to the disruption of a star around a supermassive black hole, most probably spinning with a Kerr parameter $a_\bullet\gtrsim0.5$.

Using few-body simulations, we investigate the evolution of supermassive black holes (SMBHs) in galaxies ($M_{\star}=10^{10}-10^{12}{\rm M}_{\odot}$ at $z=0$) at $0<z<4$. Following galaxy merger trees from the Millennium simulation, we model BH mergers with two extreme binary decay scenarios for the `hard binary' stage: a full or an empty loss cone. These two models should bracket the true evolution, and allow us to separately explore the role of dynamical friction and that of multi-body BH interactions on BH mergers. Using the computed merger rates, we infer the stochastic gravitational wave background (GWB). Our dynamical approach is a first attempt to study the dynamical evolution of multiple SMBHs in the host galaxies undergoing mergers with various mass ratios ($10^{-4} < q_{\star} < 1$). Our main result demonstrates that SMBH binaries are able to merge in both scenarios. In the empty loss cone case, we find that BHs merge via multi-body interactions, avoiding the `final parsec' problem, and entering the PTA band with substantial orbital eccentricity. Our full loss cone treatment, albeit more approximate, suggests that the eccentricity becomes even higher when GWs become dominant, leading to rapid coalescences (binary lifetime $\lesssim1 {\rm ~Gyr}$). Despite the lower merger rates in the empty loss cone case, due to their higher mass ratios and lower redshifts, the GWB in the full/empty loss cone models are comparable ($0.70\times10^{-15}$ and $0.53\times10^{-15}$ at a frequency of $1~{\rm yr}^{-1}$, respectively). Finally, we compute the effects of high eccentricities on the GWB spectrum.

Recent observations suggest that stellar tidal disruption events (TDE) are strongly overrepresented in rare, post-starburst galaxies. Several dynamical mechanisms have been proposed to elevate their TDE rates, ranging from central stellar overdensities to the presence of supermassive black hole (SMBH) binaries. Another such mechanism, introduced here, is a radial velocity anisotropy in the nuclear star cluster produced during the starburst. These, and other, dynamical hypotheses can be disentangled by comparing observations to theoretical predictions for the TDE delay time distribution (DTD). We show that SMBH binaries are a less plausible solution for the post-starburst preference, as they can only reproduce the observed DTD with extensive fine-tuning. The overdensity hypothesis produces a reasonable match to the observed DTD (based on the limited data currently available), provided that the initial stellar density profile created during the starburst, $\rho(r)$, is exceptional in both steepness and normalization. In particular, explaining the post-starburst preference requires $\rho \propto r^{-\gamma}$ with $\gamma \gtrsim 2.5$, i.e. much steeper than the classic Bahcall-Wolf equilibrium profile of $\gamma = 7/4$. For "ultrasteep" density cusps ($\gamma \ge 9/4$), we show that the TDE rate decays with time measured since the starburst as $\dot{N} \propto t^{-(4\gamma-9)/(2\gamma-3)} / \ln t$. Radial anisotropies also represent a promising explanation, provided that initial anisotropy parameters of $\beta_0 \approx 0.5$ are sustainable against the radial orbit instability. TDE rates in initially anisotropic cusps will decay roughly as $\dot{N} \propto t^{-\beta_0}$. As the sample of TDEs with well-studied host galaxies grows, the DTD will become a powerful tool for constraining the exceptional dynamical properties of post-starburst galactic nuclei.

The mass of the central black hole in a galaxy that hosted a tidal disruption event (TDE) is an important parameter in understanding its energetics and dynamics. We present the first homogeneously measured black hole masses of a complete sample of 12 optically/UV selected TDE host galaxies (down to $g_{host}$$\leq$22 mag and $z$=0.37) in the Northern sky. The mass estimates are based on velocity dispersion measurements, performed on late time optical spectroscopic observations. We find black hole masses in the range 3$\times$10$^5$ M$_{\odot}$$\leq$M$_{\rm BH}$$\leq$2$\times$10$^7$ M$_{\odot}$. The TDE host galaxy sample is dominated by low mass black holes ($\sim$10$^6$ M$_{\odot}$), as expected from theoretical predictions. The blackbody peak luminosity of TDEs with M$_{\rm BH}$$\leq$10$^{7.1}$ M$_{\odot}$ is consistent with the Eddington limit of the SMBH, whereas the two TDEs with M$_{\rm BH}$$\geq$10$^{7.1}$ M$_{\odot}$ have peak luminosities below their SMBH Eddington luminosity, in line with the theoretical expectation that the fallback rate for M$_{\rm BH}$$\geq$10$^{7.1}$ M$_{\odot}$ is sub-Eddington. In addition, our observations suggest that TDEs around lower mass black holes evolve faster. These findings corroborate the standard TDE picture in 10$^6$ M$_{\odot}$ black holes. Our results imply an increased tension between observational and theoretical TDE rates. By comparing the blackbody emission radius with theoretical predictions, we conclude that the optical/UV emission is produced in a region consistent with the stream self-intersection radius of shallow encounters, ruling out a compact accretion disk as the direct origin of the blackbody radiation at peak brightness.

When a main sequence star undergoes Roche lobe overflow onto a supermassive black hole (SMBH) in a circular extreme mass ratio inspiral (EMRI), a phase of steady mass transfer ensues. Over millions of years, the binary evolves to a period minimum before reversing course and migrating outwards. Because the time interval between consecutive EMRIs is comparable to the mass-transfer timescale, the semi-major axes of two consecutive mass-transferring EMRIs will cross on a radial scale < few AU. We show that such EMRI crossing events are inevitably accompanied by a series of mildly relativistic, grazing physical collisions between the stars. Each collision strips a small quantity of mass, primarily from the more massive star, which generally increases their radial separation to set up the next collision after a delay of decades to centuries (or longer) set by further gravitational radiation. Depending on the mass of the SMBH, this interaction can result in N ~ 1-1e4 gas production events of mass Msun/N, thus powering a quasi-periodic sequence of SMBH accretion-powered flares over a total duration of thousands of years or longer. Although the EMRI rate is 2-3 orders of magnitude lower than the rate of tidal disruption events (TDE), the ability of a single interacting EMRI pair to produce a large number of luminous flares - and to make more judicious use of the available stellar fuel - could make their observed rate competitive with the TDE rate, enabling them to masquerade as "TDE Imposters." We predict flares with luminosities that decay both as power laws shallower than t^(-5/3) or as decaying exponentials. Viscous spreading of the gas disks produced by the accumulation of previous mass-stripping events places substantial mass on radial scales > 10-100 AU, providing a reprocessing source required to explain the unexpectedly high optical luminosities of some flares.

The Kepler-field star KIC 8462852, an otherwise apparently ordinary F3 main-sequence star, showed several highly unusual dimming events of variable depth and duration. Adding to the mystery was the discovery that KIC 8462852 faded by 14% from 1890 to 1989, as well as by another 3% over the 4 year Kepler mission. Following an initial suggestion by Wright & Sigurdsson, we propose that the secular dimming behavior is the result of the inspiral of a planetary body or bodies into KIC 8462852, which took place ~10 to 1e4 years ago (depending on the planet mass). Gravitational energy released as the body inspirals into the outer layers of the star caused a temporary and unobserved brightening, from which the stellar flux is now returning to the quiescent state. The transient dimming events could then be due to obscuration by planetary debris from an earlier partial disruption of the same inspiraling bodies, or due to evaporation and out-gassing from a tidally detached moon system. Alternatively, the dimming events could arise from a large number of bodies comet- or planetesimal-mass bodies placed onto high eccentricity orbits by the same mechanism (e.g. Lidov-Kozai oscillations due to the outer M-dwarf companion) responsible for driving the more massive planets into KIC 8462852. The required high occurrence rate of KIC 8462852-like systems which have undergone recent major planet inspiral event(s) is the greatest challenge to the model, placing large lower limits on the mass of planetary systems surrounding F stars and/or requiring an unlikely probability to catch KIC 8462852 in its current state.

We carried out the first multi-wavelength (optical/UV and X-ray) photometric reverberation mapping of a tidal disruption flare (TDF) ASASSN-14li. We find that its X-ray variations are correlated with and lag the optical/UV fluctuations by 32$\pm$4 days. Based on the direction and the magnitude of the X-ray time lag, we rule out X-ray reprocessing and direct emission from a standard circular thin disk as the dominant source of its optical/UV emission. The lag magnitude also rules out an AGN disk-driven instability as the origin of ASASSN-14li and thus strongly supports the tidal disruption picture for this event and similar objects. We suggest that the majority of the optical/UV emission likely originates from debris stream self-interactions. Perturbations at the self-interaction sites produce optical/UV variability and travel down to the black hole where they modulate the X-rays. The time lag between the optical/UV and the X-rays variations thus correspond to the time taken by these fluctuations to travel from the self-interaction site to close to the black hole. We further discuss these time lags within the context of the three variants of the self-interaction model. High-cadence monitoring observations of future TDFs will be sensitive enough to detect these echoes and would allow us to establish the origin of optical/UV emission in TDFs in general.

Stellar-mass binary black hole mergers are poised to represent the majority of gravitational-wave (GW) observations by Advanced LIGO and Virgo. Probing their origin will be difficult due to the expected lack of electromagnetic emission and limited localization accuracy. Associations with rare host galaxy types -- such as active galactic nuclei (AGN) -- can nevertheless be identified statistically through spatial correlation. We show that (i) fractional contributions $f_{\rm agn}=50-100\%$ from AGN hosts to the total BBH merger rate can be statistically established with 70-300 detected events (expected in 0.5-2 years of observation with Advanced LIGO-Virgo at design sensitivity and current rate estimates); (ii) fractional contributions as low as $f_{\rm agn}=25\%$ can be tested with 1000 events ($\sim$ 5\,years of observation); (iii) the $\sim5\%$ best localized GWs drive these constraints. The presented method and results are generally applicable to binary formation channels with rare host populations.

We explore the evolution of stellar mass black hole binaries (BHBs) which are formed in the self-gravitating disks of active galactic nuclei (AGN). Hardening due to three-body scattering and gaseous drag are effective mechanisms that reduce the semi-major axis of a BHB to radii where gravitational waves take over, on timescales shorter than the typical lifetime of the AGN disk. Taking observationally-motivated assumptions for the rate of star formation in AGN disks, we find a rate of disk-induced BHB mergers ($\mathcal{R} \sim 3~{\rm yr}^{-1}~{\rm Gpc}^{-3}$, but with large uncertainties) that is comparable with existing estimates of the field rate of BHB mergers, and the approximate BHB merger rate implied by the recent Advanced LIGO detection of GW150914. BHBs formed thorough this channel will frequently be associated with luminous AGN, which are relatively rare within the sky error regions of future gravitational wave detector arrays. This channel could also possess a (potentially transient) electromagnetic counterpart due to super-Eddington accretion onto the stellar mass black hole following the merger.

Time domain optical surveys have discovered roughly a dozen candidate stellar tidal disruption flares in the last five years, and future surveys like the {\it Large Synoptic Survey Telescope} will likely find hundreds to thousands more. These tidal disruption events (TDEs) present an interesting puzzle: a majority of the current TDE sample is hosted by rare post-starburst galaxies, and tens of percent are hosted in even rarer E+A galaxies, which make up $\sim 0.1\%$ of all galaxies in the local universe. E+As are therefore overrepresented among TDE hosts by 1-2 orders of magnitude, a discrepancy unlikely to be accounted for by selection effects. We analyze {\it Hubble Space Telescope} photometry of one of the nearest E+A galaxies, NGC~3156, to estimate the rate of stellar tidal disruption produced as two-body relaxation diffuses stars onto orbits in the loss cone of the central supermassive black hole. The rate of TDEs produced by two-body relaxation in NGC~3156 is large when compared to other galaxies with similar black hole mass: $\dot{N}_{\rm TDE}\sim 1\times 10^{-3}~{\rm yr}^{-1}$. This suggests that the preference of TDEs for E+A hosts may be due to central stellar overdensities produced in recent starbursts.

When a star passes within the tidal radius of a supermassive black hole, it will be torn apart. For a star with the mass of the Sun ($M_\odot$) and a non-spinning black hole with a mass $<10^8 M_\odot$, the tidal radius lies outside the black hole event horizon and the disruption results in a luminous flare. Here we report observations over a period of 10 months of a transient, hitherto interpreted as a superluminous supernova. Our data show that the transient rebrightened substantially in the ultraviolet and that the spectrum went through three different spectroscopic phases without ever becoming nebular. Our observations are more consistent with a tidal disruption event than a superluminous supernova because of the temperature evolution, the presence of highly ionised CNO gas in the line of sight and our improved localisation of the transient in the nucleus of a passive galaxy, where the presence of massive stars is highly unlikely. While the supermassive black hole has a mass $> 10^8 M_\odot$, a star with the same mass as the Sun could be disrupted outside the event horizon if the black hole were spinning rapidly. The rapid spin and high black hole mass can explain the high luminosity of this event.

Dozens of stellar tidal disruption events (TDEs) have been identified at optical, UV and X-ray wavelengths. A small fraction of these, most notably Swift J1644+57, produce radio synchrotron emission, consistent with a powerful, relativistic jet shocking the surrounding circumnuclear gas. The dearth of similar non-thermal radio emission in the majority of TDEs may imply that powerful jet formation is intrinsically rare, or that the conditions in galactic nuclei are typically unfavorable for producing a detectable signal. Here we explore the latter possibility by constraining the radial profile of the gas density encountered by a TDE jet using a one-dimensional model for the circumnuclear medium which includes mass and energy input from a stellar population. Near the jet Sedov radius of 10$^{18}$ cm, we find gas densities in the range of $n_{18} \sim$ 0.1$-$1000 cm$^{-3}$ across a wide range of plausible star formation histories. Using one- and two-dimensional relativistic hydrodynamical simulations, we calculate the synchrotron radio light curves of TDE jets (as viewed both on and off-axis) across the allowed range of density profiles. We find that bright radio emission would be produced across the plausible range of nuclear gas densities by jets as powerful as Swift J1644+57, and we quantify the relationship between the radio luminosity and jet energy. We use existing radio detections and upper limits to constrain the energy distribution of TDE jets. Radio follow up observations several months to several years after the TDE candidate will strongly constrain the energetics of any relativistic flow.

We study the circularization of tidally disrupted stars on bound orbits around spinning supermassive black holes by performing three-dimensional smoothed particle hydrodynamic simulations with Post-Newtonian corrections. Our simulations reveal that debris circularization depends sensitively on the efficiency of radiative cooling. There are two stages in debris circularization if radiative cooling is inefficient: first, the stellar debris streams self-intersect due to relativistic apsidal precession; shocks at the intersection points thermalize orbital energy and the debris forms a geometrically thick, ring-like structure around the black hole. The ring rapidly spreads via viscous diffusion, leading to the formation of a geometrically thick accretion disk. In contrast, if radiative cooling is efficient, the stellar debris circularizes due to self-intersection shocks and forms a geometrically thin ring-like structure. In this case, the dissipated energy can be emitted during debris circularization as a precursor to the subsequent tidal disruption flare. The possible radiated energy is up to ~2*10^{52} erg for a 1 Msun star orbiting a 10^6 Msun black hole. We also find that a retrograde (prograde) black hole spin causes the shock-induced circularization timescale to be shorter (longer) than that of a non-spinning black hole in both cooling cases. The circularization timescale is remarkably long in the radiatively efficient cooling case, and is also sensitive to black hole spin. Specifically, Lense-Thirring torques cause dynamically important nodal precession, which significantly delays debris circularization. On the other hand, nodal precession is too slow to produce observable signatures in the radiatively inefficient case. We also discuss the relationship between our simulations and the parabolic TDEs that are characteristic of most stellar tidal disruptions.

In this paper, we study the chaotic four-body problem in Newtonian gravity. Assuming point particles and total encounter energies $\le$ 0, the problem has three possible outcomes. We describe each outcome as a series of discrete transformations in energy space, using the diagrams first presented in Leigh \& Geller (2012; see the Appendix). Furthermore, we develop a formalism for calculating probabilities for these outcomes to occur, expressed using the density of escape configurations per unit energy, and based on the Monaghan description originally developed for the three-body problem. We compare this analytic formalism to results from a series of binary-binary encounters with identical point particles, simulated using the \texttt{FEWBODY} code. Each of our three encounter outcomes produces a unique velocity distribution for the escaping star(s). Thus, these distributions can potentially be used to constrain the origins of dynamically-formed populations, via a direct comparison between the predicted and observed velocity distributions. Finally, we show that, for encounters that form stable triples, the simulated single star escape velocity distributions are the same as for the three-body problem. This is also the case for the other two encounter outcomes, but only at low virial ratios. This suggests that single and binary stars processed via single-binary and binary-binary encounters in dense star clusters should have a unique velocity distribution relative to the underlying Maxwellian distribution (provided the relaxation time is sufficiently long), which can be calculated analytically.

We consider the origins of enigmatic stellar populations in four Local Group galactic nuclei, specifically the Milky Way, M31, M32 and M33. These are centrally concentrated blue stars, found in three out of the four nuclear star clusters (NSCs) considered here. Their origins are unknown, but could include blue straggler (BS) stars, extended horizontal branch stars and young recently formed stars. Here, we calculate order-of-magnitude estimates for various collision rates, as a function of the host NSC environment and distance from the cluster centre. These rates are sufficiently high that BSs, formed via collisions between main sequence (MS) stars, could contribute non-negligibly ($\sim$ 1-10% in mass) to every surface brightness profile, with the exception of the Milky Way. Stellar evolution models show that the envelopes of red giant branch (RGB) stars must be nearly completely stripped to significantly affect their photometric appearance, which requires multiple collisions. Hence, the collision rates for individual RGB stars are only sufficiently high in the inner $\lesssim$ 0.1 pc of M31 and M32 for RGB destruction to occur. Collisions between white dwarfs and MS stars, which should ablate the stars, could offer a steady and significant supply of gas in every NSC in our sample. The gas could either fragment to form new stars, or accrete onto old MS stars already present. Thus, collisional processes could contribute significantly to the observed blue excesses in M31 and M33; future studies should be aimed at better constraining theoretical predictions to compliment existing and future observational data.

The magnetic dipole radiation (MDR) model is currently the best approach we have to explain pulsar radiation. However a most characteristic parameter of the observed radiation, the braking index n$_{\rm obs}$ shows deviations for all the eight best studied isolated pulsars, from the simple model prediction n$_{\rm dip}$ = 3. The index depends upon the rotational frequency and its first and second time derivatives, but also on the assumption of that the magnetic dipole moment and inclination angle, and the moment of inertia of the pulsar are constant in time. In a recent paper [Phys. Rev. D 91, 063007 (2015)] we showed conclusively that changes in the moment of inertia with frequency alone, cannot explain the observed braking indices. Possible observational evidence for the magnetic dipole moment migrating away from the rotational axis at a rate $\dot\alpha$ $\sim$ 0.6$^{\circ}$ per 100 years over the life time of the Crab pulsar has been recently suggested by Lyne et al. In this paper, we explore the MDR model with constant moment of inertia and magnetic dipole moment but variable inclination angle $\alpha$. We first discuss the effect of the variation of $\alpha$ on the observed braking indices and show they all can be understood. However, no explanation for the origin of the change in $\alpha$ is provided. After discussion of the possible source(s) of magnetism in pulsars we propose a simple mechanism for the change in $\alpha$ based on a toy model in which the magnetic structure in pulsars consists of two interacting dipoles. We show that such a system can explain the Crab observation and the measured braking indices.

Nuclear star clusters (NSCs) and supermassive black holes (SMBHs) both inhabit galactic nuclei, coexisting in a range of bulge masses, but excluding each other in the largest or smallest galaxies. We propose that the transformation of NSCs into SMBHs occurs via runaway tidal captures, once NSCs exceed a certain critical central density and velocity dispersion. The bottleneck in this process, as with all collisional runaways, is growing the first e-fold in black hole mass. The growth of a stellar mass black hole past this bottleneck occurs as tidally captured stars are consumed in repeated episodes of mass transfer at pericenter. Tidal captures may turn off as a growth channel once the black hole reaches a mass ~100-1000 solar masses, but tidal disruption events will continue and appear capable of growing the seed SMBH to larger sizes. The runaway slows (becomes sub-exponential) once the seed SMBH consumes the core of its host NSC. While the bulk of the cosmic mass density in SMBHs is ultimately produced (via the Soltan-Paczynski argument) by episodic gaseous accretion in very massive galaxies, the smallest SMBHs have probably grown from strong tidal encounters with NSC stars. SMBH seeds that grow for a time $t$ entirely through this channel will follow simple power law relations with the velocity dispersion, $\sigma$, of their host galaxy. In the simplest regime it is $M_\bullet \sim \sigma^{3/2}\sqrt{M_\star t / G} \sim 10^{6}M_\odot (\sigma / 50~{\rm km~s}^{-1})^{3/2}(t/10^{10}~{\rm yr})^{1/2}$, but the exponents and prefactor can differ slightly depending on the details of loss cone refilling. Current tidal disruption event rates predicted from this mechanism are consistent with observations.

Rates of stellar tidal disruption events (TDEs) by supermassive black holes (SMBHs) due to two-body relaxation are calculated using a large galaxy sample (N=146) in order to explore the sensitivity of the TDE rates to observational uncertainties, such as the parametrization of galaxy light profiles and the stellar mass function. The largest uncertainty arises due to the poorly constrained occupation fraction of SMBHs in low-mass galaxies, which otherwise dominate the total TDE rate. The detection rate of TDE flares by optical surveys is calculated as a function of SMBH mass and other observables for several physically-motivated models of TDE emission. We also quantify the fraction of galaxies that produce deeply penetrating disruption events. If the majority of the detected events are characterized by super-Eddington luminosities (such as disk winds, or synchrotron radiation from an off-axis relativistic jet), then the measured SMBH mass distribution will tightly constrain the low-end SMBH occupation fraction. If Eddington-limited emission channels dominate, however, then the occupation fraction sensitivity is much less pronounced in a flux-limited survey (although still present in a volume-complete event sample). The SMBH mass distribution of the current sample of TDEs, though highly inhomogeneous and encumbered by selection effects, already suggests that Eddington-limited emission channels dominate. Even our most conservative rate estimates appear to be in tension with much lower observationally inferred TDE rates, and we discuss several possible resolutions to this discrepancy.

The tidal disruption of a star by a supermassive black hole leads to a short-lived thermal flare. Despite extensive searches, radio follow-up observations of known thermal stellar tidal disruption flares (TDFs) have not yet produced a conclusive detection. We present a detection of variable radio emission from a thermal TDF, which we interpret as originating from a newly-launched jet. The multi-wavelength properties of the source present a natural analogy with accretion state changes of stellar mass black holes, suggesting all TDFs could be accompanied by a jet. In the rest frame of the TDF, our radio observations are an order of magnitude more sensitive than nearly all previous upper limits, explaining how these jets, if common, could thus far have escaped detection.

We calculate steady-state, one-dimensional hydrodynamic profiles of hot gas in slowly accreting ("quiescent") galactic nuclei for a range of central black hole masses $M_{\bullet}$, parametrized gas heating rates, and observationally-motivated stellar density profiles. Mass is supplied to the circumnuclear medium by stellar winds, while energy is injected primarily by stellar winds, supernovae, and black hole feedback. Analytic estimates are derived for the stagnation radius (where the radial velocity of the gas passes through zero) and the large scale gas inflow rate, $\dot{M}$, as a function of $M_{\bullet}$ and the gas heating efficiency, the latter being related to the star-formation history. We assess the conditions under which radiative instabilities develop in the hydrostatic region near the stagnation radius, both in the case of a single burst of star formation and for the average star formation history predicted by cosmological simulations. By combining a sample of measured nuclear X-ray luminosities, $L_x$, of nearby quiescent galactic nuclei with our results for $\dot{M}(M_{\bullet})$ we address whether the nuclei are consistent with accreting in a steady-state, thermally-stable manner for radiative efficiencies predicted for radiatively inefficiency accretion flows. We find thermally-stable accretion cannot explain the short average growth times of low mass black holes in the local Universe, which must instead result from gas being fed in from large radii, due either to gas inflows or thermal instabilities acting on larger, galactic scales. Our results have implications for attempts to constrain the occupation fraction of SMBHs in low mass galaxies using the mean $L_x-M_{\bullet}$ correlation, as well as the predicted diversity of the circumnuclear densities encountered by relativistic outflows from tidal disruption events.

Several lines of observational evidence suggest that white dwarfs receive small birth kicks due to anisotropic mass loss. If other stars possess extrasolar analogues to the Solar Oort cloud, the orbits of comets in such clouds will be scrambled by white dwarf natal kicks. Although most comets will be unbound, some will be placed on low angular momentum orbits vulnerable to sublimation or tidal disruption. The dusty debris from these comets will manifest itself as an IR excess temporarily visible around newborn white dwarfs; examples of such disks may already have been seen in the Helix Nebula, and around several other young white dwarfs. Future observations with the James Webb Space Telescope may distinguish this hypothesis from alternatives such as a dynamically excited Kuiper Belt analogue. Although competing hypotheses exist, the observation that $\gtrsim 15\%$ of young white dwarfs possess such disks, if interpreted as indeed being cometary in origin, provides indirect evidence that low mass gas giants (thought necessary to produce an Oort cloud) are common in the outer regions of extrasolar planetary systems. Hydrogen abundances in the atmospheres of older white dwarfs can, if sufficiently low, also be used to place constraints on the joint parameter space of natal kicks and exo-Oort cloud models.

We present a potential-density pair designed to model nearly isothermal star clusters (and similar self-gravitating systems) with a central core and an outer turnover radius, beyond which density falls off as $r^{-4}$. In the intermediate zone, the profile is similar to that of an isothermal sphere (density $\rho \propto r^{-2}$), somewhat less steep than the King 62 profile, and with the advantage that many dynamical quantities can be written in a simple closed form. We derive new analytic expressions for the cluster binding energy and velocity dispersion, and apply these to create toy models for cluster core collapse and evaporation. We fit our projected surface brightness profiles to observed globular and open clusters, and find that the quality of the fit is generally at least as good as that for the surface brightness profiles of King 62. This model can be used for convenient computation of the dynamics and evolution of globular and nuclear star clusters.

When a star is tidally disrupted by a supermassive black hole (BH), roughly half of its mass falls back to the BH at super-Eddington rates. Being tenuously gravitationally bound and unable to cool radiatively, only a small fraction f_in << 1 of the returning debris will likely be incorporated into the disk and accrete, with the vast majority instead becoming unbound in an outflow of velocity ~1e4 km/s. This slow outflow spreads laterally, encasing the BH. For months or longer, the outflow remains sufficiently neutral to block hard EUV and X-ray radiation from the hot inner disk, which instead becomes trapped in a radiation-dominated nebula. Ionizing nebular radiation heats the inner edge of the ejecta to temperatures of T > few 1e4 K, converting the emission to optical/near-UV wavelengths where photons more readily escape due to the lower opacity. This can explain the unexpectedly low and temporally constant effective temperatures of optically-discovered TDE flares. For BHs with relatively high masses M_BH > 1e7 M_sun the ejecta can become ionized at an earlier stage, or for a wider range of viewing angles, producing a TDE flare which is instead dominated by thermal X-ray emission. We predict total radiated energies consistent with those of observed TDE flares, and ejecta velocities that agree with the measured emission line widths. The peak optical luminosity for M_BH < 1e6 M_sun is suppressed due to adiabatic losses in the inner disk wind, possibly contributing to the unexpected dearth of optical TDEs in galaxies with low mass BHs. In the classical picture, where f_in ~ 1, TDEs de-spin supermassive BHs and cap their maximum spins well below theoretical accretion physics limits. This cap is greatly relaxed in our model, and existing Fe K-alpha spin measurements provide suggestive preliminary evidence that f_in < 1.

The incompressibility (compression modulus) $K_{\rm 0}$ of infinite symmetric nuclear matter at saturation density has become one of the major constraints on mean-field models of nuclear many-body systems as well as of models of high density matter in astrophysical objects and heavy-ion collisions. We present a comprehensive re-analysis of recent data on GMR energies in even-even $^{\rm 112-124}$Sn and $^{\rm 106,100-116}$Cd and earlier data on 58 $\le$ A $\le$ 208 nuclei. The incompressibility of finite nuclei $K_{\rm A}$ is expressed as a leptodermous expansion with volume, surface, isospin and Coulomb coefficients $K_{\rm vol}$, $K_{\rm surf}$, $K_\tau$ and $K_{\rm coul}$. \textit{Assuming} that the volume coefficient $K_{\rm vol}$ is identified with $K_{\rm 0}$, the $K_{\rm coul}$ = -(5.2 $\pm$ 0.7) MeV and the contribution from the curvature term K$_{\rm curv}$A$^{\rm -2/3}$ in the expansion is neglected, compelling evidence is found for $K_{\rm 0}$ to be in the range 250 $ < K_{\rm 0} < $ 315 MeV, the ratio of the surface and volume coefficients $c = K_{\rm surf}/K_{\rm vol}$ to be between -2.4 and -1.6 and $K_{\rm \tau}$ between -840 and -350 MeV. We show that the generally accepted value of $K_{\rm 0}$ = (240 $\pm$ 20) MeV can be obtained from the fits provided $c \sim$ -1, as predicted by the majority of mean-field models. However, the fits are significantly improved if $c$ is allowed to vary, leading to a range of $K_{\rm 0}$, extended to higher values. A self-consistent simple (toy) model has been developed, which shows that the density dependence of the surface diffuseness of a vibrating nucleus plays a major role in determination of the ratio K$_{\rm surf}/K_{\rm vol}$ and yields predictions consistent with our findings.

The precise origin of short gamma ray bursts (SGRBs) remains an important open question in relativistic astrophysics. Increasingly, observational evidence suggests the merger of a binary compact object system as the source for most SGRBs, but it is currently unclear how to distinguish observationally between a binary neutron star progenitor and a black hole-neutron star progenitor. We suggest the quasi-periodic signal of jet precession as an observational signature of SGRBs originating in mixed binary systems, and quantify both the fraction of mixed binaries capable of producing SGRBs, and the distributions of precession amplitudes and periods. The difficulty inherent in disrupting a neutron star outside the horizon of a stellar-mass black hole biases the jet precession signal towards low amplitude and high frequency. Precession periods of ~ 0.01-0.1 s and disk-BH spin misalignments ~10 degrees are generally expected, although sufficiently high viscosity may prevent the accumulation of multiple precession periods during the SGRB. The precessing jet will naturally cover a larger solid angle in the sky than would standard SGRB jets, enhancing observability for both prompt emission and optical afterglows.

The tidal disruption of a star by a supermassive black hole (SMBH) is a highly energetic event with consequences dependent on the degree to which the star plunges inside the SMBH's tidal sphere. We introduce a new analytic model for tidal disruption events (TDEs) to analyze the dependence of these events on beta, the ratio of the tidal radius to the orbital pericenter. We find, contrary to most previous work, that the spread in debris energy for a TDE is largely constant for all beta. This result has important consequences for optical transient searches targeting TDEs, which we discuss. We quantify leading-order general relativistic corrections to this spread in energy and find that they are small. We also examine the role of stellar spin, and find that a combination of spin-orbit misalignment, rapid rotation, and high beta may increase the spread in debris energy. Finally, we quantify for the first time the gravitational wave emission due to the strong compression of a star in a high-beta TDE. Although this signal is unlikely to be detectable for disruptions of main sequence stars, the tidal disruption of a white dwarf by an intermediate mass black hole can produce a strong signal visible to Advanced LIGO at tens of megaparsecs.

We study accretion processes for tidally disrupted stars approaching supermassive black holes on bound orbits, by performing three dimensional Smoothed Particle Hydrodynamics simulations with a pseudo-Newtonian potential. We find that there is a critical value of the orbital eccentricity below which all the stellar debris remains bound to the black hole. For high but sub-critical eccentricities, all the stellar mass is accreted onto the black hole in a finite time, causing a significant deviation from the canonical $t^{-5/3}$ mass fallback rate. When a star is on a moderately eccentric orbit and its pericenter distance is deeply inside the tidal disruption radius, there can be several orbit crossings of the debris streams due to relativistic precession. This dissipates orbital energy in shocks, allowing for rapid circularization of the debris streams and formation of an accretion disk. The resultant accretion rate greatly exceeds the Eddington rate and differs strongly from the canonical rate of $t^{-5/3}$. By contrast, there is little dissipation due to orbital crossings for the equivalent simulation with a purely Newtonian potential. This shows that general relativistic precession is crucial for accretion disk formation via circularization of stellar debris from stars on moderately eccentric orbits.

We study tidal disruption and subsequent mass fallback for stars approaching supermassive black holes on bound orbits, by performing three dimensional Smoothed Particle Hydrodynamics simulations with a pseudo-Newtonian potential. We find that the mass fallback rate decays with the expected -5/3 power of time for parabolic orbits, albeit with a slight deviation due to the self-gravity of the stellar debris. For eccentric orbits, however, there is a critical value of the orbital eccentricity, significantly below which all of the stellar debris is bound to the supermassive black hole. All the mass therefore falls back to the supermassive black hole in a much shorter time than in the standard, parabolic case. The resultant mass fallback rate considerably exceeds the Eddington accretion rate and substantially differs from the -5/3 power of time.

When a star is tidally disrupted by a supermassive black hole (SMBH), the streams of liberated gas form an accretion disk after their return to pericenter. We demonstrate that Lense-Thirring precession in the spacetime around a rotating SMBH can produce significant time evolution of the disk angular momentum vector, due to both the periodic precession of the disk and the nonperiodic, differential precession of the bound debris streams. Jet precession and periodic modulation of disk luminosity are possible consequences. The persistence of the jetted X-ray emission in the Swift J164449.3+573451 flare suggests that the jet axis was aligned with the spin axis of the SMBH during this event.

We analyze stellar tidal disruption events as a possible observational signature of gravitational wave induced recoil of supermassive black holes. As a black hole wanders through its galaxy, it will tidally disrupt bound and unbound stars at rates potentially observable by upcoming optical transient surveys. To quantify these rates, we explore a broad range of host galaxy and black hole kick parameters. We find that emission from a transient accretion disk can produce ~1 event per year which LSST would identify as spatially offset, while super-Eddington tidal flares, if they exist, are likely to produce ~10 spatially offset events per year. A majority of tidal disruption flares, and a large majority of flares with an observable spatial offset, are due to bound rather than unbound stars. The total number of disruption events due to recoiled black holes could be almost 1% of the total stellar tidal disruption rate.

A precise electromagnetic measurement of the sky coordinates and redshift of a coalescing black hole binary holds the key for using its gravitational wave (GW) signal to constrain cosmological parameters and to test general relativity. Here we show that the merger of ~10^{6-7}M_sun black holes is generically followed over a period of years by multiple electromagnetic flares from tidally disrupted stars. The sudden recoil imparted to the merged black hole by GW emission promptly fills its loss cone and results in a tidal disruption rate of stars as high as ~0.1 per year. The prompt disruption of a star within a single galaxy over a short period provides a unique electromagnetic flag of a recent black hole coalescence event, and sequential disruptions could be used on their own to calibrate the expected rate of GW sources for pulsar timing arrays or the proposed Laser Interferometer Space Antenna (LISA).