The rate of tidal disruption events (TDEs) depends sensitively on the stellar properties of the central galactic regions. Simulations show that galaxy mergers cause gas inflows that trigger nuclear starbursts that increase the central stellar density. Motivated by these numerical results, and by the observed over-representation of post-starburst galaxies among TDE hosts, we study the evolution of the TDE rate in high-resolution hydrodynamical simulations of a galaxy merger, in which we capture the evolution of the stellar density around the massive black holes (BHs) in the merging galaxies and in the final remnant. We post-process the simulation to obtain the TDE rate. At the second pericentre, a nuclear starburst enhances the stellar density around the BH in the least massive galaxy, leading to an enhancement of the TDE rate around the secondary BH, although the magnitude and the duration of the increase depend on the stochasticity of star formation on very small scales. The central stellar density around the primary BH remains instead fairly constant, and so is its TDE rate, until the two BHs bind and are surrounded by the same stellar environment. After the formation of the binary, the stellar density decreases, and so does the TDE rate.

Dark matter may be composed of light bosons, ${m_b \sim 10^{-22}\, \mathrm{eV}}$, with a de Broglie wavelength $\lambda \sim 1 \,\mathrm{kpc}$ in typical galactic potentials. Such `fuzzy' dark matter (FDM) behaves like cold dark matter (CDM) on much larger scales than the de Broglie wavelength, but may resolve some of the challenges faced by CDM in explaining the properties of galaxies on small scales ($\lesssim 10\,\mathrm{kpc}$). Because of its wave nature, FDM exhibits stochastic density fluctuations on the scale of the de Broglie wavelength that never damp. The gravitational field from these fluctuations scatters stars and black holes, causing their orbits to diffuse through phase space. We show that this relaxation process can be analyzed quantitatively with the same tools used to analyze classical two-body relaxation in an $N$-body system, and can be described by treating the FDM fluctuations as quasiparticles, with effective mass $\sim 10^7 M_\odot {(1\,\mathrm{kpc}/r)}^2{(10^{-22}\,\mathrm{eV}/m_b)}^3$ in a galaxy with a constant circular speed of $200\,\mathrm{kms}$. This novel relaxation mechanism may stall the inspiral of supermassive black holes or globular clusters due to dynamical friction at radii of a few hundred pc, and can heat and expand the central regions of galaxies. These processes can be used to constrain the mass of the light bosons that might comprise FDM.

In the vicinity of a massive black hole, stars move on precessing Keplerian orbits. The mutual stochastic gravitational torques between the stellar orbits drive a rapid reorientation of their orbital planes, through a process called vector resonant relaxation. We derive, from first principles, the correlation of the potential fluctuations in such a system, and the statistical properties of random walks undergone by the stellar orbital orientations. We compare this new analytical approach with effective $N$-body simulations. We also provide a simple scheme to generate the random walk of a test star's orbital orientation using a stochastic equation of motion. We finally present quantitative estimations of this process for a nuclear stellar cluster such as the one of the Milky Way.

Binaries within the sphere of influence of a massive black hole (MBH) in galactic nuclei are susceptible to the Lidov-Kozai (LK) mechanism, which can drive orbits to high eccentricities and trigger strong interactions within the binary such as the emission of gravitational waves (GWs), and mergers of compact objects. These events are potential sources for GW detectors such as Advanced LIGO and VIRGO. The LK mechanism is only effective if the binary is highly inclined with respect to its orbit around the MBH (within a few degrees of 90 deg), implying low rates. However, close to an MBH, torques from the stellar cluster give rise to the process of vector resonant relaxation (VRR). VRR can bring a low-inclination binary into an `active' LK regime in which high eccentricities and strong interactions are triggered in the binary. Here, we study the coupled LK-VRR dynamics, with implications for LIGO and VIRGO GW sources. We carry out Monte Carlo simulations and find that the merger fraction enhancement due to LK-VRR dynamics is up to a factor of ~10 for the lower end of assumed MBH masses (M_MBH = 1e4 MSun), and decreases sharply with increasing M_MBH. We find that, even in our most optimistic scenario, the baseline BH-BH merger rate is small, and the enhancement by LK-VRR coupling is not large enough to increase the rate to well above the LIGO/VIRGO lower limit, 12 Gpc^{-3} yr^{-1}. For the Galactic Center, the LK-VRR-enhanced rate is ~100 times lower than the LIGO/VIRGO limit, and for M_MBH = 1e4 MSun, the rate barely reaches 12 Gpc^{-3} yr^{-1}.

In nuclear star clusters, the potential is governed by the central massive black hole, so that stars move on nearly Keplerian orbits and the total potential is almost stationary in time. Yet, the deviations of the potential from the Keplerian one, due to the enclosed stellar mass and general relativity, will cause the stellar orbits to precess. Moreover, as a result of the finite number of stars, small deviations of the potential from spherical symmetry induce residual torques that can change the stars' angular momentum faster than the standard two-body relaxation. The combination of these two effects drives a stochastic evolution of orbital angular momentum, a process named "resonant relaxation". Owing to recent developments in the description of the relaxation of self-gravitating systems, we can now fully describe scalar resonant relaxation (relaxation of the magnitude of the angular momentum) as a diffusion process. In this framework, the potential fluctuations due to the complex orbital motion of the stars are described by a random correlated noise, whose statistical properties are fully characterized by the stars' mean field motion. On long timescales, the cluster can be regarded as a diffusive system, whose diffusion coefficients depend explicitly on the mean field stellar distribution through the properties of the noise. We show here, for the first time, how the diffusion coefficients of scalar resonant relaxation, for a spherically symmetric system, can be fully calculated from first principles, without any free parameters. We also provide an open source code that evaluates these diffusion coefficients numerically.

The long timescale evolution of a self-gravitating system is generically driven by two-body encounters. In many cases, the motion of the particles is primarily governed by the mean field potential. When this potential is integrable, particles move on nearly fixed orbits, which can be described in terms of angle-action variables. The mean field potential drives fast orbital motions (angles) whose associated orbits (actions) are adiabatically conserved on short dynamical timescales. The long-term stochastic evolution of the actions is driven by the potential fluctuations around the mean field and in particular by "resonant two-body encounters", for which the angular frequencies of two particles are in resonance. We show that the stochastic gravitational fluctuations acting on the particles can generically be described by a correlated Gaussian noise. Using this approach, the so-called $\eta$-formalism, we derive a diffusion equation for the actions in the test particle limit. We show that in the appropriate limits, this diffusion equation is equivalent to the inhomogeneous Balescu-Lenard and Landau equations. This approach provides a new view of the resonant diffusion processes associated with long-term orbital distortions. Finally, by investigating the example of the Hamiltonian Mean Field Model, we show how the present method generically allows for alternative calculations of the long-term diffusion coefficients in inhomogeneous systems.

Intermediate-mass black holes (IMBHs) of mass $M_{\bullet} \approx 10^{2} - 10^{5}$ solar masses, $M_{\odot}$, are the long-sought missing link between stellar black holes, born of supernovae, and massive black holes, tied to galaxy evolution by the empirical $M_{\bullet}/\sigma_{\star}$ correlation. We show that low-mass black hole seeds that accrete stars from locally dense environments in galaxies following a universal $M_{\bullet}/\sigma_{\star}$ relation grow over the age of the Universe to be above ${\mathcal{M}}_{0}\approx3\times10^{5}M_{\odot}$ ($5\%$ lower limit), independent of the unknown seed masses and formation processes. The mass ${\mathcal{M}}_{0}$ depends weakly on the uncertain formation redshift, and sets a universal minimal mass scale for present-day black holes. This can explain why no IMBHs have yet been found, and it implies that present-day galaxies with ${\sigma_{\star}<{\mathcal{S}}_{0}\approx40\,\mathrm{km\,s}^{-1}}$ lack a central black hole, or formed it only recently. A dearth of IMBHs at low redshifts has observable implications for tidal disruptions and gravitational wave mergers.

A massive black hole (MBH) consumes stars whose orbits evolve into the small phase-space volume of unstable orbits, the "loss-cone", which take them directly into the MBH, or close enough to interact strongly with it. The resulting phenomena: tidal heating and tidal disruption, binary capture and hyper-velocity star ejection, gravitational wave (GW) emission by inspiraling compact remnants, or hydrodynamical interactions with an accretion disk, are of interest as they can produce observable signatures and thereby reveal the existence of the MBH, affect its mass and spin evolution, probe strong gravity, and provide information on stars and gas near the MBH. The continuous loss of stars and the processes that resupply them shape the central stellar distribution. We investigate relativistic stellar dynamics near the loss-cone of a non-spinning MBH in steady-state analytically and by Monte Carlo simulations of the diffusion of the orbital parameters. These take into account Newtonian mass precession due to enclosed stellar mass, in-plane precession due to general relativity, dissipation by GW, uncorrelated two-body relaxation, correlated resonant relaxation (RR) and adiabatic invariance due to secular precession, using a rigorously derived description of correlated post-Newtonian dynamics in the diffusion limit. We argue that general maximal entropy considerations strongly constrain orbital diffusion in steady-state, irrespective of the relaxation mechanism. We identify the exact phase-space separatrix between plunges and inspirals, predict their steady-state rates, and verify they are robust under a wide range of assumptions. We derive the dependence of the rates on the mass of the MBH, show that the contribution of RR is small, and discuss special cases where unquenched RR in restricted volumes of phase-space may affect the steady-state substantially.

Stars around a massive black hole (MBH) move on nearly fixed Keplerian orbits, in a centrally-dominated potential. The random fluctuations of the discrete stellar background cause small potential perturbations, which accelerate the evolution of orbital angular momentum by resonant relaxation. This drives many phenomena near MBHs, such as extreme mass-ratio gravitational wave inspirals, the warping of accretion disks, and the formation of exotic stellar populations. We present here a formal statistical mechanics framework to analyze such systems, where the background potential is described as a correlated Gaussian noise. We derive the leading order, phase-averaged 3D stochastic Hamiltonian equations of motion, for evolving the orbital elements of a test star, and obtain the effective Fokker-Planck equation for a general correlated Gaussian noise, for evolving the stellar distribution function. We show that the evolution of angular momentum depends critically on the temporal smoothness of the background potential fluctuations. Smooth noise has a maximal variability frequency $\nu_{\max}$. We show that in the presence of such noise, the evolution of the normalized angular momentum $j=\sqrt{1-e^{2}}$ of a relativistic test star, undergoing Schwarzschild (in-plane) General Relativistic precession with frequency $\nu_{GR}/j^{2}$, is exponentially suppressed for $j<j_{b}$, where $\nu_{GR}/j_{b}^{2}\sim\nu_{\max}$, due to the adiabatic invariance of the precession against the slowly varying random background torques. This results in an effective Schwarzschild precession-induced barrier in angular momentum. When $j_{b}$ is large enough, this barrier can have significant dynamical implications for processes near the MBH.

[abridged] Energy relaxation around a massive black hole (MBH) is key to establishing the dynamical state of galactic nuclei, and the nature of close stellar interactions with the MBH. The standard description of relaxation as diffusion provides a perturbative 2nd-order solution in the weak two-body interaction limit. We run N-body simulations and find that this solution fails to describe the non-Gaussian relaxation on short timescale, which is strongly influenced by extreme events even in the weak limit, and is thus difficult to characterize and measure. We derive a non-perturbative solution for relaxation as an anomalous diffusion process, and develop a robust estimation technique to measure it in simulations. These enable us to analyze and model our numerical results, and validate in detail, for the first time, this model of energy relaxation around an MBH on all timescales. We derive the relation between the energy diffusion time, t_E, and the time for a small perturbation to return to steady state, t_r, in a relaxed, single mass cusp around a MBH. We constrain the contribution of strong encounters, measure that of the weakest encounters, determine the value of the Coulomb logarithm, and provide a robust analytical estimate for t_E in a finite nuclear stellar cusp. We find that t_r ~ 10t_E ~(5/32)Q^2P_h/N_h log Q, where Q=M_bh/M_* is the MBH to star mass ratio, the orbital period P_h and number of stars N_h are evaluated at the energy scale corresponding to the MBH's sphere of influence, E_h=sigma_inf^2, where sigma_inf is the velocity dispersion far from the MBH. We conclude, using the observed cosmic M_bh/sigma correlation, that cusps around lower-mass MBHs (M_bh<10^7 Mo), which evolved passively over a Hubble time, should be relaxed. We consider the effects of anomalous energy diffusion on orbital perturbations of stars observed near the Galactic MBH.