In the local universe, black holes of $10^{5-6}$ M$_{\odot}$ are hosted in galaxies displaying a variety of stellar profiles and morphologies. These black holes are the anticipated targets of LISA, the Laser Interferometer Space Antenna that will detect the low-frequency gravitational-wave signal emitted by binary black holes in this mass interval. In this paper, we infer upper limits on the lifetime of binary black holes of $10^{5-6}$ M$_{\odot}$ and up to $10^8$ M$_{\odot}$, forming in galaxy mergers, exploring two underlying stellar density profiles, by Dehnen and by Prugniel & Simien, and by exploiting local scaling relations between the mass of the black holes and several quantities of their hosts. We focus on the phase of the dynamical evolution when the binary is transitioning from the hardening phase ruled by the interaction with single stars to the phase driven by the emission of gravitational waves. We find that different stellar profiles predict very distinct trends with binary mass, with lifetimes ranging between fractions of a Gyr to more than 10 Gyr, and with a spread of about one order of magnitude, given by the uncertainties in the observed correlations, which are larger in the low-mass tail of the observed black hole population.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

Theoretical models for the expected merger rates of intermediate-mass black holes (IMBHs) are vital for planned gravitational-wave detection experiments such as the Laser Interferometer Space Antenna (LISA). Using collisionless $N$-body simulations of dwarf galaxy (DG) mergers, we examine how the orbital decay of IMBHs and the efficiency of IMBH binary formation depend on the central dark matter (DM) density profile of the merging DGs. Specifically, we explore various asymptotic inner slopes $\gamma$ of the DG's DM density distribution, ranging from steep cusps ($\gamma=1$) to shallower density profiles ($\gamma<1$), motivated by well-known baryonic-feedback effects as well as by DM models that differ from cold DM at the scales of DGs. We find that the inner DM slope is crucial for the formation (or lack thereof) of an IMBH binary; only mergers between DGs with cuspy DM profiles ($\gamma=1$) are favourable to forming a hard IMBH binary, whereas when $\gamma<1$ the IMBHs stall at a separation of 50-100 pc. Consequently, the rate of LISA signals from IMBH coalescence will be determined by the fraction of DGs with a cuspy DM profile. Conversely, the LISA event rates at IMBH mass scales offer in principle a novel way to place constraints on the inner structure of DM halos in DGs and address the core-cusp controversy. We also show that, with spatial resolutions of $\sim$0.1 kpc, as often adopted in cosmological simulations, all IMBHs stall, independent of $\gamma$. This suggests caution in employing cosmological simulations of galaxy formation to study BH dynamics in DGs.

We present the first 2D hydrodynamical finite volume simulations in which dust is fully coupled with the gas, including its back-reaction onto it, and at the same time the dust size is evolving according to coagulation and fragmentation based on a sub-grid model. The aim of this analysis is to present the differences occurring when dust evolution is included relative to simulations with fixed dust size, with and without an embedded Jupiter-mass planet that triggers gap formation. We use the two-fluid polar Godunov-type code RoSSBi developed by Surville et al. 2016 combined with a new local sub-grid method for dust evolution based on the model by Birnstiel et al. 2012. We find striking differences between simulations with variable and fixed dust sizes. The timescales for dust depletion differ significantly and yield a completely different evolution of the dust surface density. In general sharp features such as pile-ups of dust in the inner disk and near gap edges, when a massive planet is present, become much weaker. This has important implications on the interpretation of observed substructure in disks, suggesting that the presence of a massive planet does not necessarily cause sharp gaps and rings in the dust component. Also, particles with different dust sizes show a different distribution, pointing to the importance of multi-wavelength synthetic observations in order to compare with observations by ALMA and other instruments. We also find that simulations adopting fixed intermediate particle sizes, in the range $10^{-2} - 10^{-1}$ cm, best approximate the surface density evolution seen in simulations with dust evolution.