The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN) -- active SMBHs at projected separations larger than several parsecs -- and binary AGN -- probing the scale where SMBHs are bound in a Keplerian binary -- is an observational challenge. The study of AGN pairs (either dual or binary) also represents an overarching theoretical problem in cosmology and astrophysics. The AGN triggering calls for detailed knowledge of the hydrodynamical conditions of gas in the imminent surroundings of the SMBHs and, at the same time, their duality calls for detailed knowledge on how galaxies assemble through major and minor mergers and grow fed by matter along the filaments of the cosmic web. This review describes the techniques used across the electromagnetic spectrum to detect dual and binary AGN candidates and proposes new avenues for their search. The current observational status is compared with the state-of-the-art numerical simulations and models for formation of dual and binary AGN. Binary SMBHs are among the loudest sources of gravitational waves (GWs) in the Universe. The search for a background of GWs at nHz frequencies from inspiralling SMBHs at low redshifts, and the direct detection of signals from their coalescence by the Laser Interferometer Space Antenna in the next decade, make this a theme of major interest for multi-messenger astrophysics. This review discusses the future facilities and observational strategies that are likely to significantly advance this fascinating field.

The gravitational interaction at parsec to sub-parsec scales between a circumbinary gas disc and a super massive black hole binary (SMBHB) is a promising mechanism to drive the migration of SMBHBs toward coalescence. The typical dynamical evolution can be separated in two main regimes: I) Slow migration ($T_{\rm mig}$ $\sim$ $10^{3-4}\times T_{\rm orb}$), where viscous torques are not efficient enough to redistribute the extracted angular momentum from the binary, leading to the formation of a low density cavity around the binary. II) Fast migration ($T_{\rm mig}$ $\sim$ $10^{1-2} \times T_{\rm orb}$), in which the redistribution of angular momentum is efficient and no low density cavity is formed in the circumbinary disc. Using N-Body/SPH simulations we study the effect of AGN feedback in this phase of a SMBHB evolution. We implement an AGN feedback model in the SPH code Gadget-3 that includes momentum feedback from winds, X-ray heating/radial-momentum and Eddington force. Using this implementation we run a set of simulations of SMBHB+disc in the two main shrinking regimes. From these simulations we conclude that the effect of the AGN mechanical feedback is negligible for SMBHBs in the slowly shrinking regime. However, in the fast shrinking regime the AGN wind excavate a "feedback cavity" leaving the SMBHB almost naked, thus stalling the orbital decay of the binary.

Galaxy observations and N-body cosmological simulations produce conflicting dark matter halo density profiles for galaxy central regions. While simulations suggest a cuspy and universal density profile (UDP) of this region, the majority of observations favor variable profiles with a core in the center. In this paper, we investigate the convergency of standard N-body simulations, especially in the cusp region, following the approach proposed by \protect{\citep{13}}. We simulate the well known Hernquist model using the SPH code Gadget-3 and consider the full array of dynamical parameters of the particles. We find that, although the cuspy profile is stable, all integrals of motion characterizing individual particles suffer strong unphysical variations along the whole halo, revealing an effective interaction between the test bodies. This result casts doubts on the reliability of the velocity distribution function obtained in the simulations. Moreover, we find unphysical Fokker-Planck streams of particles in the cusp region. The same streams should appear in cosmological N-body simulations, being strong enough to change the shape of the cusp or even to create it. Our analysis, based on the Hernquist model and the standard SPH code, strongly suggests that the UDPs generally found by the cosmological N-body simulations may be a consequence of numerical effects. A much better understanding of the N-body simulation convergency is necessary before a 'core-cusp problem' can properly be used to question the validity of the CDM model.

Using N-body/SPH simulations we study the evolution of the separation of a pair of SMBHs embedded in a star forming circumnuclear disk (CND). This type of disk is expected to be formed in the central kilo parsec of the remnant of gas-rich galaxy mergers. Our simulations indicate that orbital decay of the SMBHs occurs more quickly when the mean density of the CND is higher, due to increased dynamical friction. However, in simulations where the CND is fragmented in high density gaseous clumps (clumpy CND), the orbits of the SMBHs are erratically perturbed by the gravitational interaction with these clumps, delaying, in some cases, the orbital decay of the SMBHs. The densities of these gaseous clumps in our simulations and in recent studies of clumpy CNDs are significantly higher than the observed density of molecular clouds in isolated galaxies or ULIRGs, thus, we expect that SMBH orbits are perturbed less in real CNDs than in the simulated CNDs of this study and other recent studies. We also find that the migration timescale has a weak dependence on the star formation rate of the CND. Furthermore, the migration timescale of a SMBH pair in a star-forming clumpy CND is at most a factor three longer than the migration timescale of a pair of SMBHs in a CND modeled with more simple gas physics. Therefore, we estimate that the migration timescale of the SMBHs in a clumpy CND is on the order of $10^7$ yrs.

We study the interaction between an unequal mass binary with an isothermal circumbinary disk, motivated by the theoretical and observational evidence that after a major merger of gas-rich galaxies, a massive gaseous disk with a SMBH binary will be formed in the nuclear region. We focus on the gravitational torques that the binary exerts onto the disk and how these torques can drive the formation of a gap in the disk. This exchange of angular momentum between the binary and the disk is mainly driven by the gravitational interaction between the binary and a strong non-axisymmetric density perturbation that is produced in the disk, as response to the presence of the binary. Using SPH numerical simulations we tested two gap-opening criterion, one that assumes that the geometry of the density perturbation is an ellipsoid/thick-spirals and another that assumes a geometry of flat-spirals for the density perturbation. We find that the flat-spirals gap opening criterion successfully predicts which simulations will have a gap on the disk and which simulations will not have a gap on the disk. We also study the limiting cases predicted by the gap-opening criteria. Since the viscosity in our simulations is considerably smaller than the expected value in the nuclear regions of gas-rich merging galaxies, we conclude that in such environments the formation of a circumbinary gap is unlikely.

We study the interaction of an equal mass binary with an isothermal circumbinary disk motivated by the theoretical and observational evidence of the formation of massive black holes binaries surrounded by gas, after a major merger of gas-rich galaxies. We focus on the torques that the binary produces on the disk and how the exchange of angular momentum can drive the formation of a gap on it. We propose that the angular momentum exchange between the binary and the disk is through the gravitational interaction of the binary and a (tidally formed) global non-axisymmetric perturbation in the disk. Using this gravitational interaction we derive an analytic criterion for the formation of a gap in the disk that can be expressed on the structural parameters h/a and M(< r)/M_{bin}. Using SPH simulations we show that the simulations where the binary opens a gap in the disk and the simulations where the disk does not have a gap are distributed in two well separate regions. Our analytic gap-opening criterion predicts a shape of the threshold between this two regions that is consistent with our simulations and the other ones in the literature. We propose an analogy between the regime without (with) a gap in the disk and the Type I (Type II) migration that is observed in simulations of planet-disk interaction (binaries with extreme mass ratios), emphasizing that the interaction that drives the formation of a gap on the disk is different in the regime that we analyze (comparable mass binary).

We study the gravitational stability of gaseous streams in the complex environment of a galaxy merger, because mergers are known to be places of ongoing massive cluster formation and bursts of star formation. We find an analytic stability parameter for case of gaseous streams orbiting around the merger remnant. We test our stability criteria using hydrodynamical simulations of galaxy mergers, obtaining satisfactory results. We find that our criteria successfully predicts the streams that will be gravitationally unstable to fragment into clumps.

Stochastic impacts by large bodies are, at present, the usually accepted mechanisms able to account for the obliquity of the ice giants. We attempt to set constraints on giant impacts as the cause of Neptune's current obliquity in the framework of modern theories. We also use the present orbital properties of the Neptunian irregular satellites (with the exception of Triton) to set constraints on the scenario of giant impacts at the end of Neptune formation. We model the angular momentum transfer to proto-Neptune and the impulse transfer to its irregular satellites by the last stochastic collision (GC) between the protoplanet and an oligarchic mass at the end of Neptune's formation. We obtain that an impactor mass greather than 4 Earth masses is not possible since it cannot reproduce the present rotational properties of the planet, unless the impact parameter of the collision were very small. On the other hand, if the impactor mass was greather than 1.4 Earth masses, the present Neptunian irregular satellites had to be formed or captured after the end of stochastic impacts. The upper bounds on the oligarchic masses (4 Earth masses from the obliquity of Neptune and 1.4 earth masses from the Neptunian irregular satellites) are independent of unknown parameters, such as the mass and distribution of the planetesimals, the location at which Uranus and Neptune were formed, the Solar Nebula initial surface mass density, and the growth regime. If stochastic impacts had occurred, these results should be understood as upper constraints on the oligarchic masses in the trans-Saturnian region at the end of ice planet formation and may be used to set constraints on planetary formation scenarios.