We study the dynamical evolution of eccentric massive black hole binaries (MBHBs) interacting with unbound stars by means of an extensive set of three body scattering experiments. Compared to previous studies, we extend the investigation down to a MBHB mass ratio of $q=m_2/m_1=10^{-4}$, where $m_1$ and $m_2$ are the masses of the primary and secondary hole respectively. Contrary to a simple extrapolation from higher mass ratios, we find that for $q\lesssim 10^{-3}$ the eccentricity growth rate becomes negative, i.e., the binary {\it circularises} as it shrinks. This behaviour is due to the subset of interacting stars captured in metastable counter-rotating orbits; those stars tend to extract angular momentum from the binary, promoting eccentricity growth for $q>10^{-3}$, but tend to inject angular momentum into the binary driving it towards circularisation for $q<10^{-3}$. The physical origin of this behaviour requires a detailed study of the orbits of this subset of stars and is currently under investigation. Our findings might have important consequences for intermediate MBHs (IMBHs) inspiralling onto MBHs (e.g. a putative $10^3\rm M_{\odot}$ black hole inspiralling onto SgrA$^*$).

We consider the formation of binary intermediate black holes (BIMBH) in globular clusters (GC), which could happen either in situ or due to the mergers between clusters. We simulate the evolution of the BIMBH orbit (and its subsequent merger) due to stellar ejections. We also take into account the evaporation of GCs due to the tidal field of the host galaxy and two-body relaxation. Our results show that if at least $10^{-3}$ of all GCs become BIMBH hosts and the BIMBH masses are $\sim1\%$ of the GC mass, at least one of the inspiralling (or merging) BIMBHs will be detected by LISA during its 4-year mission lifetime. Most of the detected BIMBHs come 1) from heavy GCs ($\gtrsim 3\times10^5M_\odot$), as lower-mass GCs end up being disrupted before their BIMBHs have time to merge, and 2) from redshifts $1<z<3$, assuming that most of GCs form around $z\sim4$ and given that the merger timescale for most BIMBHs is $\sim1$ Gyr. If the BIMBH to GC mass ratio is lower ($\sim10^{-3}$) but the fraction of BIMBH hosts among GCs is higher ($\gtrsim10^{-2}$), some of their mergers will also be detected by LIGO, VIRGO, and KAGRA and the proposed Einstein Telescope.

We consider a black hole (BH) density cusp in a nuclear star cluster (NSC) hosting a supermassive back hole (SMBH) at its center. Assuming the stars and BHs inside the SMBH sphere of influence are mass-segregated, we calculate the number of BHs that sink into this region under the influence of dynamical friction. We find that the total number of BHs increases significantly in this region due to this process for lower mass SMBHs by up to a factor of 5, but there is no increase in the vicinity of the highest mass SMBHs. Due to the high BH number density in the NSC, BH-BH binaries form during close approaches due to GW emission. We update the previous estimate of O'Leary et al. for the rate of such GW capture events by estimating the $\langle n^2\rangle/\langle n\rangle^2$ parameter where $n$ is the number density. We find a BH merger rate for this channel to be in the range $\sim0.01-0.1 \, \mathrm{Gpc^{-3}yr^{-1}}$. The total merger rate is dominated by the smallest galaxies hosting SMBHs, and the number of heaviest BHs in the NSC. It is also exponentially sensitive to the radial number density profile exponent, reaching $>100 \, \mathrm{Gpc^{-3}yr^{-1}}$ when the BH mass function is $m^{-2.3}$ or shallower and the heaviest BH radial number density is close to $r^{-3}$. Even if the rate is much lower than the range constrained by the current LIGO detections, the GW captures around SMBHs can be distinguished by their high eccentricity in the LIGO band.

In this paper we consider a scenario where the currently observed hypervelocity stars in our Galaxy have been ejected from the Galactic center as a result of dynamical interactions with an intermediate-mass black hole (IMBH) orbiting the central supermassive black hole (SMBH). By performing 3-body scattering experiments, we calculate the distribution of the ejected stars' velocities given various parameters of the IMBH-SMBH binary: IMBH mass, semimajor axis and eccentricity. We also calculate the rates of change of the BH binary orbital elements due to those stellar ejections. One of our new findings is that the ejection rate depends (although mildly) on the rotation of the stellar nucleus (its total angular momentum). We also compare the ejection velocity distribution with that produced by the Hills mechanism (stellar binary disruption) and find that the latter produces faster stars on average. Also, the IMBH mechanism produces an ejection velocity distribution which is flattened towards the BH binary plane while the Hills mechanism produces a spherically symmetric one. The results of this paper will allow us in the future to model the ejection of stars by an evolving BH binary and compare both models with \textit{Gaia} observations, for a wide variety of environments (galactic nuclei, globular clusters, the Large Magellanic Clouds, etc.).

Pulsar timing array (PTA) collaborations in North America, Australia, and Europe, have been exploiting the exquisite timing precision of millisecond pulsars over decades of observations to search for correlated timing deviations induced by gravitational waves (GWs). PTAs are sensitive to the frequency band ranging just below 1 nanohertz to a few tens of microhertz. The discovery space of this band is potentially rich with populations of inspiraling supermassive black-holes binaries, decaying cosmic string networks, relic post-inflation GWs, and even non-GW imprints of axionic dark matter. This article aims to provide an understanding of the exciting open science questions in cosmology, galaxy evolution, and fundamental physics that will be addressed by the detection and study of GWs through PTAs. The focus of the article is on providing an understanding of the mechanisms by which PTAs can address specific questions in these fields, and to outline some of the subtleties and difficulties in each case. The material included is weighted most heavily towards the questions which we expect will be answered in the near-term with PTAs; however, we have made efforts to include most currently anticipated applications of nanohertz GWs.

The recent discovery of super-massive black holes (SMBHs) in the centers of high-mass ultra compact dwarf galaxies (UCDs) suggests that at least some UCDs are the stripped nuclear star clusters of lower mass galaxies. Tracing these former nuclei of stripped galaxies provides a unique way to track the assembly history of a galaxy or galaxy cluster. In this paper we present a new method to estimate how many UCDs host an SMBH in their center and thus are stripped galaxy nuclei. We revisit the dynamical mass measurements that suggest many UCDs have more mass than expected from stellar population estimates, which recent observations have shown is due to the presence of an SMBH. We revise the stellar population mass estimates using a new empirical relation between the mass-to-light ratio (M/L) and metallicity, and use this to predict which UCDs are most likely to host an SMBH. This enables us to calculate the fraction of UCDs that host SMBHs across their entire luminosity range for the first time. We then apply the SMBH occupation fraction to the observed luminosity function of UCDs and estimate that in the Fornax and Virgo cluster alone there should be $69^{+32}_{-25}$ stripped nuclei with SMBHs. This analysis shows that stripped nuclei with SMBHs are almost as common in clusters as present-day galaxy nuclei. We estimate the local SMBH number density in stripped nuclei to $3-8\times10^{-3}Mpc^{-3}$, which represents a significant fraction (10-40\%) of the SMBH density in the local Universe. These SMBHs hidden in stripped nuclei will increase expected event rates for tidal disruption events and SMBH-SMBH and SMBH-BH mergers. The existence of numerous stripped nuclei with SMBHs are a direct consequence of hierarchical galaxy formation, but until now their impact on the SMBH density had not been quantified.

We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released $11$-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no significant evidence for a GWB, we place constraints on a GWB from a population of supermassive black-hole binaries, cosmic strings, and a primordial GWB. For the first time, we find that the GWB upper limits and detection statistics are sensitive to the Solar System ephemeris (SSE) model used, and that SSE errors can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first PTA constraints that are robust against SSE uncertainties. We thus place a $95\%$ upper limit on the GW strain amplitude of $A_\mathrm{GWB}<1.45\times 10^{-15}$ at a frequency of $f=1$ yr$^{-1}$ for a fiducial $f^{-2/3}$ power-law spectrum, and with inter-pulsar correlations modeled. This is a factor of $\sim 2$ improvement over the NANOGrav $9$-year limit, calculated using the same procedure. Previous PTA upper limits on the GWB will need revision in light of SSE systematic uncertainties. We use our constraints to characterize the combined influence on the GWB of the stellar mass-density in galactic cores, the eccentricity of SMBH binaries, and SMBH--galactic-bulge scaling relationships. We constrain cosmic-string tension using recent simulations, yielding an SSE-marginalized $95\%$ upper limit on the cosmic string tension of $G\mu < 5.3\times 10^{-11}$---a factor of $\sim 2$ better than the published NANOGrav $9$-year constraints. Our SSE-marginalized $95\%$ upper limit on the energy density of a primordial GWB (for a radiation-dominated post-inflation Universe) is $\Omega_\mathrm{GWB}(f)h^2<3.4\times10^{-10}$.

We compute the isotropic gravitational wave (GW) background produced by binary supermassive black holes (SBHs) in galactic nuclei. In our model, massive binaries evolve at early times via gravitational-slingshot interaction with nearby stars, and at later times by the emission of GWs. Our expressions for the rate of binary hardening in the "stellar" regime are taken from the recent work of Vasiliev et al., who show that in the non-axisymmetric galaxies expected to form via mergers, stars are supplied to the center at high enough rates to ensure binary coalescence on Gyr timescales. We also include, for the first time, the extra degrees of freedom associated with evolution of the binary's orbital plane; in rotating nuclei, interaction with stars causes the orientation and the eccentricity of a massive binary to change in tandem, leading in some cases to very high eccentricities (e>0.9) before the binary enters the GW-dominated regime. We argue that previous studies have over-estimated the mean ratio of SBH mass to galaxy bulge mass by factors of 2 - 3. In the frequency regime currently accessible to pulsar timing arrays (PTAs), our assumptions imply a factor 2 - 3 reduction in the characteristic strain compared with the values computed in most recent studies, removing the tension that currently exists between model predictions and the non-detection of GWs.

Interaction of a binary supermassive black hole with stars in a galactic nucleus can result in changes to all the elements of the binary's orbit, including the angles that define its orientation. If the nucleus is rotating, the orientation changes can be large, causing large changes in the binary's orbital eccentricity as well. We present a general treatment of this problem based on the Fokker-Planck equation for f, defined as the probability distribution for the binary's orbital elements. First- and second-order diffusion coefficients are derived for the orbital elements of the binary using numerical scattering experiments, and analytic approximations are presented for some of these coefficients. Solutions of the Fokker-Planck equation are then derived under various assumptions about the initial rotational state of the nucleus and the binary hardening rate. We find that the evolution of the orbital elements can become qualitatively different when we introduce nuclear rotation: 1) the orientation of the binary's orbit evolves toward alignment with the plane of rotation of the nucleus; 2) binary orbital eccentricity decreases for aligned binaries and increases for counter-aligned ones. We find that the diffusive (random-walk) component of a binary's evolution is small in nuclei with non-negligible rotation, and we derive the time-evolution equations for the semimajor axis, eccentricity and inclination in that approximation. The aforementioned effects could influence gravitational wave production as well as the relative orientation of host galaxies and radio jets.