Numerical relativity predicts that the coalescence of a black hole-binary causes the newly formed black hole to recoil, and evidence for such recoils has been found in the gravitational waves observed during the merger of stellar-mass black holes. Recoiling (super)massive black holes are expected to reside in hypercompact stellar clusters (HCSCs). Simulations of galaxy assembly predict that hundreds of HCSCs should be present in the halo of a Milky Way-type galaxy, and a fraction of those around the Milky Way should have magnitudes within the sensitivity limit of existing surveys. However, recoiling black holes and their HCSCs are still waiting to be securely identified. With the goal of enabling searches through recent and forthcoming databases, we improve over existing literature to produce realistic renditions of HCSCs bound to black holes with a mass of 10$^{5}$ M$_{\odot}$. Including the effects of a population of blue stragglers, we simulate their appearance in Pan-STARRS and in forthcoming $Euclid$ images. We also derive broad-band spectra and the corresponding multi-wavelength colours, finding that the great majority of the simulated HCSCs fall on the colour-colour loci defined by stars and galaxies, with their spectra resembling those of giant K-type stars. We discuss the clusters properties, search strategies, and possible interlopers.

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.

We report the detection of an intriguing parsec-scale radio source in the "offset AGN" candidate, KISSR 102. The elliptical host galaxy includes two optical nuclei at a projected separation of 1.54 kpc, N1 and N2, to the south-east and north-west, respectively. Phase-referenced VLBA observations at 1.5 and 4.9 GHz of this LINER galaxy, have detected double radio components (A and B) at a projected separation of 4.8 parsec at 1.5 GHz, and another partially-resolved double radio structure at 4.9 GHz coincident with the brighter radio component A. These radio detections are confined to the optical nucleus N1. The brightness temperatures of all the detected radio components are high, $\gtrsim10^8$ K, consistent with them being components of a radio AGN. The 1.5-4.9 GHz spectral index is inverted ($\alpha\sim+0.64\pm0.08$) for component A and steep for component B ($\alpha \lesssim-1.6$). The dramatic change in the spectral indices of A and B is inconsistent with it being a typical "core-jet" structure from a single AGN or the mini-lobes of a compact symmetric object. To be consistent with a "core-jet" structure, the jet in KISSR 102 would need to be undergoing strong jet-medium interaction with dense surrounding media resulting in a drastic spectral steepening of the jet. Alternatively, the results could be consistent with the presence of a parsec-scale binary radio AGN, which is the end result of a three-body interaction involving three supermassive black holes in the centre of KISSR 102.

NGC 3393, a nearby Seyfert 2 galaxy with nuclear radio jets, large-scale and nuclear bars, and a posited secondary super massive black hole, provides an interesting laboratory to test the physics of inflows and outflows. Here we present and analyse the molecular gas (ALMA observations of CO J:2-1 emission over a field of view (FOV) of 45\arcsec \times 45\arcsec, at 0\farcs56 (143 pc) spatial and 5 km/s spectral resolution), ionised gas and stars (GEMINI-GMOS/IFU; over a FOV of 4\arcsec \times 5\arcsec, at 0\farcs62 (159 pc) spatial and 23 km/s spectral resolution) in NGC 3393. The ionised gas emission, detected over the complete GEMINI-GMOS FOV, has three identifiable kinematic components. A narrow (\sigma < 115 km/s) component present in the complete FOV, which is consistent with rotation in the galaxy disk. A broad (\sigma > 115 km/s) redshifted component, detected near the NE and SW radio lobes; which we interpret as a radio jet driven outflow. And a broad (\sigma > 115 km/s) blueshifted component that shows high velocities in a region perpendicular to the radio jet axis; we interpret this as an equatorial outflow. The CO J:2-1 emission is detected in spiral arms on 5\arcsec - 20\arcsec scales, and in two disturbed circumnuclear regions. The molecular kinematics in the spiral arms can be explained by rotation. The highly disturbed kinematics of the inner region can be explained by perturbations induced by the nuclear bar and interactions with the large scale bar. We find no evidence for, but cannot strongly rule out, the presence of the posited secondary black hole.

We present Chandra and optical observations of a candidate dual AGN discovered serendipitously while searching for recoiling black holes via a cross-correlation between the serendipitous XMM source catalog (2XMMi) and SDSS-DR7 galaxies with a separation no larger than ten times the sum of their Petrosian radii. The system has a stellar mass ratio M$_{1}$/M$_{2}\approx 0.7$. One of the galaxies (Source 1) shows clear evidence for AGN activity in the form of hard X-ray emission and optical emission-line diagnostics typical of AGN ionisation. The nucleus of the other galaxy (Source 2) has a soft X-ray spectrum, bluer colours, and optical emission line ratios dominated by stellar photoionisation with a "composite" signature, which might indicate the presence of a weak AGN. When plotted on a diagram with X-ray luminosity vs [OIII] luminosity both nuclei fall within the locus defined by local Seyfert galaxies. From the optical spectrum we estimate the electron densities finding n$_{1} < 27$ e$^{-}$ cm$^{-3}$ and n$_{2} \approx 200$ e$^{-}$ cm$^{-3}$. From a 2D decomposition of the surface brightness distribution we infer that both galaxies host rotationally supported bulges (Sersic index $< 1$). While the active nature of Source 1 can be established with confidence, whether the nucleus of Source 2 is active remains a matter of debate. Evidence that a faint AGN might reside in its nucleus is, however, tantalising.

We aim to map the distribution and kinematics of molecular and ionized gas in a sample of active galaxies, to quantify the nuclear inflows and outflows. Here, we analyze the nuclear kinematics of NGC 1566 via ALMA observations of the CO J:2-1 emission at 24 pc spatial and $\sim$2.6 km s$^{-1}$ spectral resolution, and Gemini-GMOS/IFU observations of ionized gas emission lines and stellar absorption lines at similar spatial resolution, and 123 km s$^{-1}$ of intrinsic spectral resolution. The morphology and kinematics of stellar, molecular (CO) and ionized ([N II]) emission lines are compared to the expectations from rotation, outflows, and streaming inflows. While both ionized and molecular gas show rotation signatures, there are significant non-circular motions in the innermost 200 pc and along spiral arms in the central kpc (CO). The nucleus shows a double-peaked CO profile (Full Width at Zero Intensity of 200 km s$^{-1}$), and prominent ($\sim$80 km s$^{-1}$) blue and redshifted lobes are found along the minor axis in the inner arcseconds. Perturbations by the large-scale bar can qualitatively explain all features in the observed velocity field. We thus favour the presence of a molecular outflow in the disk with true velocities of $\sim$180 km s$^{-1}$ in the nucleus and decelerating to 0 by $\sim$72 pc. The implied molecular outflow rate is $5.6~[M_{o}yr^{-1}]$, with this gas accumulating in the nuclear 2 arcsec arms. The ionized gas kinematics support an interpretation of a similar, but more spherical, outflow in the inner 100 pc, with no signs of deceleration. There is some evidence of streaming inflows of $\sim$50 km s$^{-1}$ along specific spiral arms, and the estimated molecular mass inflow rate, $\sim0.1~[M_{o}yr^{-1}]$, is significantly larger than the SMBH accretion rate ($\dot{m}=4.8\times10^{-5}~[M_{o}yr^{-1}]$).

We present two-dimensional stellar and gaseous kinematics of the inner 0.7 $\times$ 1.2 kpc$^{2}$ of the Seyfert galaxy ESO 362-G18, derived from optical spectra obtained with the GMOS/IFU on the Gemini South telescope at a spatial resolution of $\approx$170 pc and spectral resolution of 36 km s$^{-1}$. ESO 362-G18 is a strongly perturbed galaxy of morphological type Sa or S0/a, with a minor merger approaching along the NE direction. Previous studies have shown that the [OIII] emission shows a fan-shaped extension of $\approx$ 10\arcsec\ to the SE. We detect the [OIII] doublet, [NII] and H${\alpha}$ emission lines throughout our field of view. The stellar kinematics is dominated by circular motions in the galaxy plane, with a kinematic position angle of $\approx$137$^{\circ}$. The gas kinematics is also dominated by rotation, with kinematic position angles ranging from 122$^{\circ}$ to 139$^{\circ}$. A double-Gaussian fit to the [OIII]$\lambda$5007 and H${\alpha}$ lines, which have the highest signal to noise ratios of the emission lines, reveal two kinematic components: (1) a component at lower radial velocities which we interpret as gas rotating in the galactic disk; and (2) a component with line of sight velocities 100-250 km s$^{-1}$ higher than the systemic velocity, interpreted as originating in the outflowing gas within the AGN ionization cone. We estimate a mass outflow rate of 7.4 $\times$ 10$^{-2}$ M$_{\odot}$ yr$^{-1}$ in the SE ionization cone (this rate doubles if we assume a biconical configuration), and a mass accretion rate on the supermassive black hole (SMBH) of 2.2 $\times$ 10$^{-2}$ M$_{\odot}$ yr$^{-1}$. The total ionized gas mass within $\sim$84 pc of the nucleus is 3.3 $\times$ 10$^{5}$ M$_{\odot}$; infall velocities of $\sim$34 km s$^{-1}$ in this gas would be required to feed both the outflow and SMBH accretion.

We present two-dimensional maps of emission-line fluxes and kinematics, as well as of the stellar kinematics of the central few kpc of five bright nearby Seyfert galaxies -- Mrk\,6, Mrk\,79, Mrk\,348, Mrk\,607 and Mrk\,1058 -- obtained from observations with the Gemini Multi-Object Spectrograph (GMOS) Integral Field Unit (IFU) on the Gemini North Telescope. The data cover the inner 3\farcs5$\times$5\farcs0 -- corresponding to physical scales in the range 0.6$\times$0.9 to 1.5$\times$2.2\,kpc$^2$ -- at a spatial resolution ranging from 110 to 280 pc with a spectral coverage of 4300 -- 7100\,\AA\ and velocity resolution of $\approx$ 90\,km\,s$^{-1}$. The gas excitation is Seyfert like everywhere but show excitation, but show excitation gradients that are correlated with the gas kinematics, reddening and/or the gas density. The gas kinematics show in all cases two components: a rotation one similar to that observed in the stellar velocity field, and an outflow component. In the case of Mrk607, the gas is counter-rotating relative to the stars. Enhanced gas velocity dispersion is observed in association to the outflows according to two patterns: at the locations of the highest outflow velocities along the ionization axis or perpendicularly to it in a strip centered at the nucleus that we attribute to an equatorial outflow. Bipolar outflows are observed in Mrk\,348 and Mrk\,79, while in Mrk\,1058 only the blueshifted part is clearly observed, while in the cases of Mrk\,6 and Mrk\,607 the geometry of the outflow needs further constraints from modeling to be presented in a forthcoming study, where the mass flow rate and powers will also be obtained.

We use optical spectra from the inner 1.8 $\times$ 2.5kpc$^2$ of the Seyfert 2 galaxy NGC1358, obtained with the GMOS integral field spectrograph on the Gemini South telescope at a spatial resolution of $\approx$ 165pc, to assess the feeding and feedback processes in this nearby active galaxy. Five gaseous kinematical components are observed in the emission line profiles. One of the components is present in the entire field-of-view and we interpret it as due to gas rotating in the disk of the galaxy. Three of the remaining components we interpret as associated to active galactic nucleus (AGN) feedback: a compact unresolved outflow in the inner 1 arcsec and two gas clouds observed at opposite sides of the nucleus, which we propose have been ejected in a previous AGN burst. The disk component velocity field is strongly disturbed by a large scale bar. The subtraction of a velocity model combining both rotation and bar flows reveals three kinematic nuclear spiral arms: two in inflow and one in outflow. We estimate the mass inflow rate in the inner 180pc obtaining $\dot{M}_{in}$ $\approx$ 1.5 $\times 10^{-2}$M$_{\odot}$yr$^{-1}$, about 160 times larger than the accretion rate necessary to power this AGN.

We used the Gemini Multi-Object Spectrograph Integral Field Unit to map the gas distribution, excitation and kinematics within the inner kiloparsec of four nearby low-luminosity active galaxies: NGC3982, NGC4501, NGC2787 and NGC4450. The observations cover the spectral range 5600-7000{\AA} at a velocity resolution of 120km/s and spatial resolution ranging from 50 to 70pc at the galaxies. Extended emission in H{\alpha}, [NII]{\lambda}{\lambda}6548,6583, [SII]{\lambda}{\lambda}6716,6730 over most of the field-of-view is observed for all galaxies, while only NGC3982 shows [OI]{\lambda}6300 extended emission. The H{\alpha} equivalent widths combined with the [NII]/H{\alpha} line ratios reveal that NGC3982 and NGC4450 harbor Seyfert nuclei surrounded by regions with LINER excitation, while NGC2787 and NGC4501 harbor LINER nuclei. NGC3982 shows a partial ring of recent star-formation at 500pc from the nucleus, while in NGC4501 a region at 500pc west of the nucleus shows LINER excitation but has been interpreted as an aging HII region with the gas excitation dominated by shocks from supernovae. The line-of-sight velocity field of the gas shows a rotation pattern for all galaxies, with deviations from pure disk rotation observed in NGC3982, NGC4501 and NGC4450. For NGC4501 and NGC4450, many of these deviations are spatially coincident with dust structures seen in optical continuum images, leading to the interpretation that the deviations are due to shocks in the gas traced by the dust. A speculation is that these shocks lead to loss of angular momentum, allowing the gas to be transferred inwards to feed the AGN. In the case of NGC2787, instead of deviations in the rotation field, we see a misalignment of 40{^\circ} between the orientation of the line of nodes of the gas rotation and the photometric major axis of the galaxy. Evidence of compact nuclear outflows are seen in NGC4501 and NGC4450.

We observed the nuclear region of the galaxy NGC 1365 with the integral field unit of the Gemini Multi Object Spectrograph mounted on the GEMINI-South telescope. The field of view covers $13^{\prime\prime} \times 6^{\prime\prime}$ ($1173 \times 541$ pc$^{2}$) centered on the nucleus, at a spatial resolution of $52$ pc. The spectral coverage extends from $5600$ \AA\ to $7000$ \AA, at a spectral resolution $R=1918$. NGC 1365 hosts a Seyfert 1.8 nucleus, and exhibits a prominent bar extending out to $100^{\prime\prime}$ (9 kpc) from the nucleus. The field of view lies within the inner Lindblad resonance. Within this region, we found that the kinematics of the ionized gas (as traced by [OI], [NII], H$\alpha$, and [SII]) is consistent with rotation in the large-scale plane of the galaxy. While rotation dominates the kinematics, there is also evidence for a fan-shaped outflow, as found in other studies based on the [OIII] emission lines. Although evidence for gas inflowing along nuclear spirals has been found in a few barred galaxies, we find no obvious signs of such features in the inner kiloparsec of NGC 1365. However, the emission lines exhibit a puzzling asymmetry that could originate from gas which is slower than the gas responsible for the bulk of the narrow-line emission. We speculate that it could be tracing gas which lost angular momentum, and is slowly migrating from the inner Lindblad resonance towards the nucleus of the galaxy.

We present optical integral field spectroscopy of the inner $2.5 \times 3.4$ kpc$^2$ of the broad-line radio galaxy Pictor A, at a spatial resolution of $\approx 400$ pc. Line emission is observed over the whole field-of-view, being strongest at the nucleus and in an elongated linear feature (ELF) crossing the nucleus from the south-west to the north-east along PA $\sim 70^\circ$. Although the broad double-peaked H$\alpha$ line and the [OI]6300/H$\alpha$ and [SII]6717+31/H$\alpha$ ratios are typical of AGNs, the [NII]6584/H$\alpha$ ratio (0.15 - 0.25) is unusually low. We suggest that this is due to the unusually low metallicity of the gas. Centroid velocity maps show mostly blueshifts to the south and redshifts to the north of the nucleus, but the velocity field is not well fitted by a rotation model. Velocity dispersions are low (< 100 km s$^{-1}$) along the ELF, ruling out a jet-cloud interaction as the origin of this structure. The ELF shows both blueshifts and redshifts in channel maps, suggesting that it is close to the plane of the sky. The ELF is evidently photoionized by the AGN, but its kinematics and inferred low metallicity suggest that this structure may have originated in a past merger event with another galaxy. We suggest that the gas acquired in this interaction may be feeding the ELF.

We present two-dimensional gaseous kinematics of the inner 1.2 $\times$ 1.8 kpc$^2$ of the Seyfert 2 galaxy NGC3081, from optical spectra (5600--7000\r{A}) obtained with the GMOS integral field spectrograph on the Gemini North telescope at a spatial resolution of $\approx$ 100pc. We have identified two-components in the line emitting gas. A narrower component (FWHM $\approx$ 60-100km s$^{-1}$), which appears to be gas in the galaxy disk, and which shows a distorted rotation pattern, is observed over the whole field of view. A broader component (FWHM $\approx$150-250 km s$^{-1}$) is present in the inner $\approx$ 2arcsec (200pc) and shows blueshifts and redshifts in the near and far sides of the galaxy, respectively, consistent with a bipolar outflow. Assuming this to be the case, we estimate that the mass outflow rate in ionized gas ($\dot{M}_{out}$) is between 1.9 $\times 10^{-3}$M$_{\odot}$ yr$^{-1}$ and 6.9 $\times 10^{-3}$M$_{\odot}$ yr$^{-1}$. The subtraction of a rotation model from the narrower component velocity field reveals a pattern of excess blueshifts of $\approx$ 50km s$^{-1}$ in the far side of the galaxy and similar excess redshifts in the near side, which are cospatial with a previously known nuclear bar. We interpret these residuals as due to gas following non-circular orbits in the barred potential. Under the assumption that these motions may lead to gas inflows, we estimate an upper limit for the mass inflow rate in ionized gas of $\phi$ $\approx$ 1.3 $\times 10^{-2}$M$_{\odot}$ yr$^{-1}$.

We present optical integral field spectroscopy of the circum-nuclear gas of the Seyfert 2 galaxy NGC 1386. The data cover the central 7$^{\prime\prime} \times 9^{\prime\prime}$ (530 $\times$ 680 pc) at a spatial resolution of 0.9" (68 pc), and the spectral range 5700-7000 \AA\ at a resolution of 66 km s$^{-1}$. The line emission is dominated by a bright central component, with two lobes extending $\approx$ 3$^{\prime\prime}$ north and south of the nucleus. We identify three main kinematic components. The first has low velocity dispersion ($\bar \sigma \approx $ 90 km s$^{-1}$), extends over the whole field-of-view, and has a velocity field consistent with gas rotating in the galaxy disk. We interpret the lobes as resulting from photoionization of disk gas in regions where the AGN radiation cones intercept the disk. The second has higher velocity dispersion ($\bar \sigma \approx$ 200 km s$^{-1}$) and is observed in the inner 150 pc around the continuum peak. This component is double peaked, with redshifted and blueshifted components separated by $\approx$ 500 km s$^{-1}$. Together with previous HST imaging, these features suggest the presence of a bipolar outflow for which we estimate a mass outflow rate of $\mathrm{\dot M} \gtrsim $ 0.1 M$_{\odot}$ yr$^{-1}$. The third component is revealed by velocity residuals associated with enhanced velocity dispersion and suggests that outflow and/or rotation is occurring approximately in the equatorial plane of the torus. A second system of velocity residuals may indicate the presence of streaming motions along dusty spirals in the disk.

The use of integral field spectroscopy is becoming increasingly popular, hovewer data reduction is still a difficult process. Here I present a step-by-step guide to the reduction of integral field data acquired with the Gemini Multi-Object Spectrograph (GMOS) on GEMINI. The reduction process, separately applied to a standard star and to the science data, includes bias and sky subtraction, flat-fielding, trimming, wavelength and flux calibration, creation of the cubes for each exposure and final combination into a single cube. Typical problems encoutered during the reduction process are discussed. The command list has been adapted from IRAF scripts given as tutorials at the South American Gemini Data Workshop (S\~{a}o Jos\'e dos Campos, Brazil, October 27-30, 2011) and scripts kindly provided by collaborators.

The coalescence of a binary black hole can be accompanied by a large gravitational recoil due to anisotropic emission of gravitational waves. A recoiling supermassive black hole (SBH) can subsequently undergo long-lived oscillations in the potential well of its host galaxy, suggesting that offset SBHs may be common in the cores of massive ellipticals. We have analyzed HST archival images of 14 nearby core ellipticals, finding evidence for small ($\lesssim 10$ pc) displacements between the AGN (locating the SBH) and the center of the galaxy (the mean photocenter) in 10 of them. Excluding objects that may be affected by large-scale isophotal asymmetries, we consider six galaxies to have detected displacements, including M87, where a displacement was previously reported by Batcheldor et al. 2010. In individual objects, these displacements can be attributed to residual gravitational recoil oscillations following a major or minor merger within the last few Gyr. For plausible merger rates, however, there is a high probability of larger displacements than those observed, if SBH coalescence took place in these galaxies. Remarkably, the AGN-photocenter displacements are approximately aligned with the radio source axis in four of the six galaxies with displacements, including three of the four having relatively powerful kpc-scale jets. This suggests intrinsic asymmetries in radio jet power as a possible displacement mechanism, although approximate alignments are also expected for gravitational recoil. Orbital motion in SBH binaries and interactions with massive perturbers can produce the observed displacement amplitudes but do not offer a ready explanation for the alignments.

We present two-dimensional gaseous kinematics of the inner 1.1 x 1.6kpc^2 of the Seyfert 2 galaxy NGC2110, from optical spectra obtained with the GMOS integral field spectrograph on the Gemini South telescope at a spatial resolution of 100pc. Gas emission is observed over the whole field-of-view, with complex - and frequently double - emission-line profiles. We have identified four components in the emitting gas, according to their velocity dispersion (sigma), which we refer to as: (1) warm gas disk (sigma = 100-220km/s); (2) cold gas disk (sigma = 60-90km/s); (3) nuclear component (sigma = 220-600km/s); and (4) northern cloud (sigma = 60-80km/s). Both the cold and warm disk components are dominated by rotation and have similar gas densities, but the cold gas disk has lower velocity dispersions and reaches higher rotation velocities. We attribute the warm gas disk to a thick gas layer which encompasses the cold disk as observed in some edge-on spiral galaxies. After subtraction of a rotation model from the cold disk velocity field, we observe excess blueshifts of 50km/s in the far side of the galaxy as well as similar excess redshifts in the near side. These residuals can be interpreted as due to nuclear inflow in the cold gas, with an estimated ionized gas mass inflow rate of 2.2 x 10^(-2)Msun/yr. We have also subtracted a rotating model from the warm disk velocity field and found excess blueshifts of 100km/s to the SW of the nucleus and excess redshifts of 40km/s to the NE, which we attribute to gas disturbed by an interaction with a nuclear spherical outflow. This nuclear outflow is the origin of the nuclear component observed within the inner 300pc and it has a mass outflow rate of 0.9Msun/yr. In a region between 1" and 4" north of the nucleus we find a new low sigma component of ionized gas which we attribute to a high latitude cloud photoionized by the nuclear source.