We study the mass assembly and spin evolution of supermassive black holes (BHs) across cosmic time as well as the impact of gravitational recoil on the population of nuclear and wandering black holes (wBHs) by using the semi-analytical model LGalaxies run on top of Millennium merger trees. We track spin changes that BHs experience during both coalescence events and gas accretion phases. For the latter, we assume that spin changes are coupled with the bulge assembly. This assumption leads to predictions for the median spin values of z=0 BHs that depend on whether they are hosted by pseudobulges, classical bulges or ellipticals, being a ~ 0.9, 0.7 and 0.4, respectively. The outcomes of the model display a good consistency with z<4 quasar luminosity functions and the z=0 BH mass function, spin values and black hole-bulge mass correlation. Regarding the wBHs, we assume that they can originate from both the disruption of satellite galaxies (orphan wBH) and ejections due to gravitational recoils (ejected wBH). The model points to a number density of wBHs that increases with decreasing redshift, although this population is always ~2 dex smaller than the one of nuclear black holes. At all redshifts, wBHs are typically hosted in M_{halo}>10^{13} M_{\odot} and M_{stellar}>10^{10} M_{\odot}, being orphan wBHs the dominant type. Besides, independently of redshift and halo mass, \textit{ejected} wBHs inhabit the central regions (<0.3R_{200}) of the host DM halo, while orphan wBHs linger at larger scales (>0.5R_{200}). Finally, we find that gravitational recoils cause a progressive depletion of nuclear BHs with decreasing redshift and stellar mass. Moreover, ejection events lead to changes in the predicted local black hole-bulge relation, in particular for BHs in pseudobulges, for which the relation is flattened at M_{bulge}>10^{10.2} M{\odot} and the scatter increase up to ~3dex.

Massive black hole (MBH) coalescences are powerful sources of low-frequency gravitational waves. To study these events in the cosmological context we need to trace the large-scale structure and cosmic evolution of a statistical population of galaxies, from dim dwarfs to bright galaxies. To cover such a large range of galaxy masses, we analyse two complementary simulations: Horizon-AGN with a large volume and low resolution which tracks the high-mass (> 1e7 Msun) MBH population, and NewHorizon with a smaller volume but higher resolution that traces the low-mass (< 1e7 Msun) MBH population. While Horizon-AGN can be used to estimate the rate of inspirals for Pulsar Timing Arrays, NewHorizon can investigate MBH mergers in a statistical sample of dwarf galaxies for LISA, which is sensitive to low-mass MBHs. After the numerical MBH merger at the resolution limit, we post-process MBH dynamics to account for time delays mostly determined by dynamical friction and stellar hardening. In both simulations, MBHs typically merge long after the galaxies do, so that the galaxy morphology at the time of the MBH merger is no longer determined by the galaxy merger from which the MBH merger originated. These time delays cause a loss of high-z MBH coalescences, shifting the peak of the MBH merger rate to z~1-2. This study shows how tracking MBH mergers in low-mass galaxies is crucial to probing the MBH merger rate for LISA and investigate the properties of the host galaxies.

The forthcoming Laser Interferometer Space Antenna (LISA) will probe the population of coalescing massive black hole (MBH) binaries up to the onset of structure formation. Here we simulate the galactic-scale pairing of $\sim10^6 M_\odot$ MBHs in a typical, non-clumpy main-sequence galaxy embedded in a cosmological environment at $z = 7-6$. In order to increase our statistical sample, we adopt a strategy that allows us to follow the evolution of six secondary MBHs concomitantly. We find that the magnitude of the dynamical-friction induced torques is significantly smaller than that of the large-scale, stochastic gravitational torques arising from the perturbed and morphologically evolving galactic disc, suggesting that the standard dynamical friction treatment is inadequate for realistic galaxies at high redshift. The dynamical evolution of MBHs is very stochastic, and a variation in the initial orbital phase can lead to a drastically different time-scale for the inspiral. Most remarkably, the development of a galactic bar in the host system either significantly accelerates the inspiral by dragging a secondary MBH into the centre, or ultimately hinders the orbital decay by scattering the MBH in the galaxy outskirts. The latter occurs more rarely, suggesting that galactic bars overall promote MBH inspiral and binary coalescence. The orbital decay time can be an order of magnitude shorter than what would be predicted relying on dynamical friction alone. The stochasticity, and the important role of global torques, have crucial implications for the rates of MBH coalescences in the early Universe: both have to be accounted for when making predictions for the upcoming LISA observatory.

We present and validate a novel semi-analytical approach to study the effect of dynamical friction on the orbits of massive perturbers in rotating stellar discs. We find that dynamical friction efficiently circularises the orbit of co-rotating perturbers, while it constantly increases the eccentricity of counter-rotating ones until their angular momenta reverse, then once again promoting circularisation. Such "drag toward circular corotation" could shape the distribution of orientations of kinematically decoupled cores in disc galaxies, naturally leading to the observed larger fraction of co-rotating cores.

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$^*$).

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.

Mergers of galaxies are extremely violent events shaping their evolution. Such events are thought to trigger starbursts and, possibly, black hole accretion. Nonetheless, it is still not clear how to know the fate of a galaxy pair from the data available at a given time, limiting our ability to constrain the exact role of mergers. In this paper we use the lightcone of the Horizon-AGN simulation, for which we know the fate of each pair, to test three selection processes aiming at identifying true merging pairs. We find that the simplest one (selecting objects within two thresholds on projected distance $d$ and redshift difference $\Delta z$) gives similar results than the most complex one (based on a neural network analyzing $d$, $\Delta z$, redshift of the primary, masses/star formation rates/aspect ratio of both galaxies). Our best thresholds are $d_\mathrm{th}\sim100\mathrm{\, kpc}$ and $\Delta z_\mathrm{th} \sim 10^{-3}$, in agreement with recent results.

We analyse the properties of strongly barred disc galaxies using the TNG100 simulation, a cosmological hydrodynamical realisation of the IllustrisTNG suite. We identify 270 disc galaxies at $z=0$ in the stellar mass range $M_{*}=10^{10.4-11}M_{\odot}$, of which 40 per cent are barred. Of the detected bars, more than half are strong. We find that the fraction of barred galaxies increases with stellar mass, in agreement with observational results. Strongly barred galaxies exhibit, overall, lower gas-to-stellar mass ratios compared to unbarred galaxies. The majority of barred galaxies are quenched (sSFR $\sim10^{-11.7} $yr$^{-1}$), whereas unbarred galaxies continue to be active (sSFR $\sim10^{-10.3}$yr$^{-1}$) on the main sequence of star-forming galaxies. We explore the evolution of strongly barred and unbarred galaxies to investigate their formation and quenching histories. We find that strong bars form between $0.5< z< 1.5$, with more massive galaxies hosting older bars. Strong bars form in galaxies with an early-established prominent disc component, undergoing periods of enhanced star formation and black hole accretion, possibly assisted by cosmological inflows. Unbarred galaxies, on the other hand, assemble most of their mass and disc component at late times. The nuclear region of strongly barred galaxies quenches shortly after bar formation, while unbarred galaxies remain active across time. Our findings are indicative of bar quenching, possibly assisted by nuclear feedback processes. We conclude that the cosmological environment, together with small scale feedback processes, determine the chances of a galaxy to form a bar and to rapidly quench its central region.

Coalescing neutron star binary (NSB) systems are primary candidates for $r$-process enrichment of galaxies. The recent detection of $r$-process elements in ultra-faint dwarf (UFD) galaxies and the abundances measured in classical dwarfs challenges the NSB merger scenario both in terms of coalescence time scales and merger locations. In this paper, we focus on the dynamics of NSBs in the gravitational potentials of different types of host galaxies and on its impact on the subsequent galactic enrichment. We find that, for a $\sim t^{-1}$ delay time distribution, even when receiving a low kick ($\sim 10~{\rm km~s^{-1}}$) from the second supernova explosion, in shallow dwarf galaxy potentials NSBs tend to merge with a large off-set from the host galaxy. This results in a significant geometrical dilution of the amount of produced $r-$process elements that fall back and pollute the host galaxy gas reservoir. The combination of dilution and small number statistics produces a large scatter in the expected $r$-process enrichment within a single UFD or classical dwarf galaxy. Comparison between our results and observed europium abundances reveals a tension that even a systematic choice of optimistic parameters in our models cannot release. Such a discrepancy could point to the need of additional $r$-process production sites that suffer less severe dilution or to a population of extremely fast merging binaries.

We report the first Atacama large millimeter/submillimeter array observations of MCG-03-58-007, a local ($z=0.03236\pm0.00002$, this work) AGN ($L_{AGN}\sim10^{45}~\rm erg~s^{-1}$), hosting a powerful X-ray ultra-fast ($v=0.1c$) outflow (UFO). The CO(1-0) line emission is observed across $\sim18\,$kpc scales with a resolution of $\sim 1\,\rm kpc$. About 78\% of the CO(1-0) luminosity traces a galaxy-size rotating disk. However, after subtracting the emission due to such rotating disk, we detect with a S/N=20 a residual emission in the central $\sim 4\,$kpc. Such residuals may trace a low velocity ($v_{LOS}=170\,\rm km\,s^{-1}$) outflow. We compare the momentum rate and kinetic power of such putative molecular outflow with that of the X-ray UFO and find $\dot{P}_{out}/\dot{P}_{UFO}=0.3\pm0.2$ and $\dot{E}_{mol}/\dot{E}_{UFO}\sim4\cdot10^{-3}$. This result is at odds with the energy-conserving scenario suggested by the large momentum boosts measured in some other molecular outflows. An alternative interpretation of the residual CO emission would be a compact rotating structure, distinct from the main disk, which would be a factor of $\sim10-100$ more extended and massive than typical circumnuclear disks revealed in Seyferts. However, in both scenarios, our results rule out the hypothesis of a momentum-boosted molecular outflow in this AGN, despite the presence of a powerful X-ray UFO. [Abridged]

Bars are a key factor in the long-term evolution of spiral galaxies, in their unique role in redistributing angular momentum and transporting gas and stars on large scales. The Eris-suite simulations are cosmological zoom-in, N-body, smoothed-particle hydrodynamic simulations built to follow the formation and evolution of a Milky Way-sized galaxy across the build-up of the large scale structure. Here we analyse and describe the outcome of two particular simulations taken from the Eris suite - ErisBH and Eris2k - which mainly differ in the prescriptions employed for gas cooling, star formation, and feedback from supernovae and black holes. Our study shows that the enhanced effective feedback in Eris2k, due to the collective effect of the different micro-physics implementations, results in a galaxy which is less massive than its ErisBH counterpart till z~1. However, when the stellar content is large enough so that global dynamical instabilities can be triggered, the galaxy in Eris2k develops a stronger and more extended bar with respect to ErisBH. We demonstrate that he structural properties and time evolution of the two bars are very different. Our results highlight the importance of accurate sub-grid prescriptions in cosmological zoom-in simulations of the process of galaxy formation and evolution, and the possible use of a statistical sample of barred galaxies to assess the strength of the stellar feedback.

We present kpc-scale ALMA and HST imaging of the quasar PJ308-21 at $z$=$6.2342$, tracing dust, gas (via the [CII] 158 $\mu$m line) and young stars. At a resolution of $\sim0.3''$ ($\approx1.7$ kpc), the system is resolved over $>4''$ ($>$20 kpc). In particular, it features a main component, identified to be the quasar host galaxy, centered on the accreting supermassive black hole; and two other extended components on the West and East side, one redshifted and the other blueshifted relative to the quasar. The [CII] emission of the entire system stretches over $>$1500 km/s along the line of sight. All the components of the system are observed in dust, [CII], and rest-frame UV emission. The inferred [CII] luminosities [(0.9-4.6)$\times 10^9$ L$_\odot$], dust luminosities [(0.15-2.6)$\times10^{12}$ L$_\odot$], and rest-frame UV luminosities [(6.6-15)$\times10^{10}$ L$_\odot$], their ratios, and the implied gas/dust masses and star formation rates [11-290 M$_\odot$ yr$^{-1}$] are typical of high-redshift star-forming galaxies. A toy model of a single satellite galaxy that is tidally stripped by the interaction with the quasar host galaxy can account for the observed velocity and spatial extent of the two extended components. An outflow interpretation of the unique features in PJ308-21 is not supported by the data. PJ308-21 is thus one of the earliest galaxy mergers imaged at cosmic dawn.

The dynamics of black hole seeds in high redshift galaxies is key to understand their ability to grow via accretion and to pair in close binaries during galactic mergers. To properly follow the dynamics of black holes we develop a physically motivated model to capture unresolved dynamical friction from stars, dark matter and gas. We first validate the model and then we use it to investigate the dynamics of seed black holes born at $z\sim9$ in dwarf proto-galaxies. We perform a suite of zoom cosmological simulations with spatial resolution as high as 10 pc and with a stellar and dark matter mass resolution of $2\times10^3 \, \textrm{M}_\odot$ and $2\times10^5 \, \textrm{M}_\odot$ respectively. We first explore the dynamics of a seed black hole in the galaxy where it is born and show that it is highly erratic if the seed mass is less than $10^5\, \textrm{M}_\odot$. The dynamics is dominated by the stellar component, whose distribution is irregular and patchy, thus inducing stochasticity in the orbits: the black hole may be anywhere in the proto-galaxy. When this dwarf merges into a larger galaxy, it is paramount to simulate the process with very high spatial and mass resolution in order to correctly account for the stripping of the stellar envelope of the satellite black hole. The outcome of the encounter could be either a tight binary or, at least temporary, a wandering black hole, leading to multiple black holes in a galaxy, each inherited from a different merger.

Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts.

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.

The emergence of a large-scale stellar bar is one of the most striking features in disc galaxies. By means of state-of-the-art cosmological zoom-in simulations, we study the formation and evolution of bars in Milky Way-like galaxies in a fully cosmological context, including the physics of gas dissipation, star formation, and supernova feedback. Our goal is to characterise the actual trigger of the non-axisymmetric perturbation that leads to the strong bar observable in the simulations at z=0, discriminating between an internal/secular versus an external/tidal origin. To this aim, we run a suite of cosmological zoom-in simulations altering the original history of galaxy-satellite interactions at a time when the main galaxy, though already bar-unstable, does not feature any non-axisymmetric structure yet. We find that the main effect of a late minor merger and of a close fly-by is to delay the time of bar formation and those two dynamical events are not directly responsible for the development of the bar and do not alter significantly its global properties (e.g. its final extension). We conclude that, once the disc has grown to a mass large enough to sustain global non-axisymmetric modes, then bar formation is inevitable.

It has been proposed that close interactions with satellite galaxies can significantly perturb the morphology of the main galaxy. However, the dynamics of an already formed bar following the interaction with the external environment has not been studied in detail in a fully cosmological context. In this work, analysing the cosmological zoom-in simulation Eris2k, we study the effects that a very unequal-mass flyby crossing the stellar disc has on the stability of the pre-existing bar. We characterize the evolution of the bar strength and length showing that the perturbation exerted by the flyby shuffles the orbits of stars for less than one Gyr. After this time, the bar shows a remarkable resilience, reforming with properties comparable to those it had before the interaction. Our work shows that close unequal-mass encounters, the most frequent interactions occurring during the evolution of cosmic structures, have (i) an overall minor impact on the global evolution of the bar in the long term, still (ii) the effect is destructive and (iii) a very weak interaction is sufficient to dismantle a strong bar leading to its "apparent death". As a consequence, due to the non-negligible duration of the bar-less period, a fraction of observed spiral galaxies classified as non-barred could be prone to bar formation.

We study a sample of Narrow-line Seyfert 1 galaxies (NLS1) in their optical and radio features, in order to understand the differences between their radio silent, radio-loud and radio-quiet subclasses. We first show that the different redshift and mass distributions of radio-loud and -quiet NLS1s could be ascribed to observational biases. We then present a geometrical model according to which most of the different observational features of radio-loud and radio-quiet NLS1s are ascribed to the orientation of an intrinsically structured jet. We estimate the fraction of intrinsically jetted sources among NLS1s that justifies the observed radio-detected population. Noticeably, under the assumptions of the geometrical model, we derive a fraction of jetted sources significantly larger than in standard AGN.

Surface abundance observations of halo stars hint at the occurrence of $r$-process nucleosynthesis at low metallicity ($\rm{[Fe/H]< -3}$), possibly within the first $10^8$ yr after the formation of the first stars. Possible loci of early-Universe $r$-process nucleosynthesis are the ejecta of either black hole--neutron star or neutron star--neutron star binary mergers. Here we study the effect of the inclination--eccentricity oscillations raised by a tertiary (e.g. a star) on the coalescence time scale of the inner compact object binaries. Our results are highly sensitive to the assumed initial distribution of the inner binary semi-major axes. Distributions with mostly wide compact object binaries are most affected by the third object, resulting in a strong increase (by more than a factor of 2) in the fraction of fast coalescences. If instead the distribution preferentially populates very close compact binaries, general relativistic precession prevents the third body from increasing the inner binary eccentricity to very high values. In this last case, the fraction of coalescing binaries is increased much less by tertiaries, but the fraction of binaries that would coalesce within $10^8$ yr even without a third object is already high. Our results provide additional support to the compact object merger scenario for $r$-process nucleosynthesis.

Massive black hole (MBH) binaries, formed as a result of galaxy mergers, are expected to harden by dynamical friction and three-body stellar scatterings, until emission of gravitational waves (GWs) leads to their final coalescence. According to recent simulations, MBH binaries can efficiently harden via stellar encounters only when the host geometry is triaxial, even if only modestly, as angular momentum diffusion allows an efficient repopulation of the binary loss cone. In this paper, we carry out a suite of N-body simulations of equal-mass galaxy collisions, varying the initial orbits and density profiles for the merging galaxies and running simulations both with and without central MBHs. We find that the presence of an MBH binary in the remnant makes the system nearly oblate, aligned with the galaxy merger plane, within a radius enclosing 100 MBH masses. We never find binary hosts to be prolate on any scale. The decaying MBHs slightly enhance the tangential anisotropy in the centre of the remnant due to angular momentum injection and the slingshot ejection of stars on nearly radial orbits. This latter effect results in about 1% of the remnant stars being expelled from the galactic nucleus. Finally, we do not find any strong connection between the remnant morphology and the binary hardening rate, which depends only on the inner density slope of the remnant galaxy. Our results suggest that MBH binaries are able to coalesce within a few Gyr, even if the binary is found to partially erase the merger-induced triaxiality from the remnant.

Models and observations suggest that luminous quasar activity is triggered by mergers, so it should preferentially occur in the most massive primordial dark matter haloes, where the frequency of mergers is expected to be the highest. Since the importance of galaxy mergers increases with redshift, we identify the high-redshift Universe as the ideal laboratory for studying dual AGN. Here we present the X-ray properties of two systems of dual quasars at z=3.0-3.3 selected from the SDSS-DR6 at separations of 6-8 arcsec (43-65kpc) and observed by Chandra for 65ks each. Both members of each pair are detected with good photon statistics to allow us to constrain the column density, spectral slope and intrinsic X-ray luminosity. We also include a recently discovered dual quasar at z=5 (separation of 21 arcsec, 136kpc) for which XMM-Newton archival data allow us to detect the two components separately. Using optical spectra we derived bolometric luminosities, BH masses and Eddington ratios that were compared to those of luminous SDSS quasars in the same redshift ranges. We find that the brighter component of both pairs at z=3.0-3.3 has high luminosities compared to the distribution of SDSS quasars at similar redshift, with J1622A having an order magnitude higher luminosity than the median. This source lies at the luminous end of the z~3.3 quasar luminosity function. While we cannot conclusively state that the unusually high luminosities of our sources are related to their having a close companion, for J1622A there is only a 3% probability that it is by chance.

Virial-based methods for estimating active supermassive black hole masses are now commonly used on extremely large spectroscopic quasar catalogues. Most spectral analyses, though, do not pay enough attention to a detailed continuum decomposition. To understand how this affects virial mass estimate results, we test the influence of host galaxy light on them, along with Balmer continuum component. A detailed fit with the new spectroscopic analysis software QSFit demonstrated that the presence or absence of continuum components do not affect significantly the virial-based results for our sample. Taking or not in consideration a host galaxy component, instead, affects the emission line fitting in a more pronounced way at lower redshifts, where in fact we observe dimmer quasars and more visible host galaxies.

The origin of the Nuclear Star Cluster in the centre of our Galaxy is still unknown. One possibility is that it formed after the disruption of stellar clusters that spiralled into the Galactic Centre due to dynamical friction. We trace the formation of the Nuclear Star Cluster around the central black hole, using state-of-the-art N-body simulations, and follow the dynamics of the neutron stars born in the clusters. We then estimate the number of Millisecond Pulsars (MSPs) that are released in the Nuclear Star Cluster, during its formation. The assembly and tidal dismemberment of globular clusters lead to a population of MSPs distributed over a radius of about 20 pc, with a peak near 3 pc. No clustering is found on the sub-parsec scale. We simulate the detectability of this population with future radio telescopes like the MeerKAT radio telescope and SKA1, and find that about of order ten MSPs can be observed over this large volume, with a paucity of MSPs within the central parsec. This helps discriminating this scenario from the in-situ formation model for the Nuclear Star Cluster that would predict an over abundance of MSPs closer to the black hole. We then discuss the potential contribution of our MSP population to the gamma-ray excess at the Galactic Centre.

We present QSFit, a new software package to automatically perform the analysis of Active Galactic Nuclei (AGN) optical spectra. The software provides luminosity estimates for the AGN continuum, the Balmer continuum, both optical and UV iron blended complex, host galaxy and emission lines, as well as width, velocity offset and equivalent width of 20 emission lines. Improving on a number of previous studies on AGN spectral analysis, QSFit fits all the components simultaneously, using an AGN continuum model which extends over the entire available spectrum, and is thus a probe of the actual AGN continuum whose estimates are scarcely influenced by localized features (e.g. emission lines) in the spectrum. We used QSFit to analyze 71,251 optical spectra of Type 1 AGN at z < 2 (obtained by the Sloan Digital Sky Survey, SDSS) and to produce a publicly available catalog of AGN spectral properties. Such catalog allowed us (for the first time) to estimate the AGN continuum slope and the Balmer continuum luminosity on a very large sample, and to show that there is no evident correlation between these quantities the redshift. All data in the catalog, the plots with best fitting model and residuals, and the IDL code we used to perform the analysis, are available on a dedicated website. The whole fitting process is customizable for specific needs, and can be extended to analyze spectra from other data sources. The ultimate purpose of QSFit is to allow astronomers to run standardized recipes to analyze the AGN data, in a simple, replicable and shareable way.

We study the dynamical evolution of supermassive black holes, in the late stage of galaxy mergers, from kpc to pc scales. In particular, we capture the formation of the binary, a necessary step before the final coalescence, and trace back the main processes causing the decay of the orbit. We use hydrodynamical simulations of galaxy mergers with different resolutions, from $20\,\rm pc$ down to $1\,\rm pc$, in order to study the effects of the resolution on our results, remove numerical effects, and assess that resolving the influence radius of the orbiting black hole is a minimum condition to fully capture the formation of the binary. Our simulations include the relevant physical processes, namely star formation, supernova feedback, accretion onto the black holes and the ensuing feedback. We find that, in these mergers, dynamical friction from the smooth stellar component of the nucleus is the main process that drives black holes from kpc to pc scales. Gas does not play a crucial role and even clumps do not induce scattering or perturb the orbits. We compare the time needed for the formation of the binary to analytical predictions and suggest how to apply such analytical formalism to obtain estimates of binary formation times in lower resolution simulations.

We investigate the simultaneous triggering of active galactic nuclei (AGN) in merging galaxies, using a large suite of high-resolution hydrodynamical simulations. We compute dual-AGN observability time-scales using bolometric, X-ray, and Eddington-ratio thresholds, confirming that dual activity from supermassive black holes (BHs) is generally higher at late pericentric passages, before a merger remnant has formed, especially at high luminosities. For typical minor and major mergers, dual activity lasts ~20-70 and ~100-160 Myr, respectively. We also explore the effects of X-ray obscuration from gas, finding that the dual-AGN time decreases at most by a factor of ~2, and of contamination from star formation. Using projected separations and velocity differences rather than three-dimensional quantities can decrease the dual-AGN time-scales by up to ~4, and we apply filters which mimic current observational-resolution limitations. In agreement with observations, we find that, for a sample of major and minor mergers hosting at least one AGN, the fraction harbouring dual AGN is ~20-30 and ~1-10 per cent, respectively. We quantify the effects of merger mass ratio (0.1 to 1), geometry (coplanar, prograde, retrograde, and inclined), disc gas fraction, and BH properties, finding that the mass ratio is the most important factor, with the difference between minor and major mergers varying between factors of a few to orders of magnitude, depending on the luminosity and filter used. We also find that a shallow imaging survey will require very high angular resolution, whereas a deep imaging survey will be less resolution-dependent.

A conclusive model for the formation of dwarf spheroidal (dSph) galaxies still remains elusive. Owing to their proximity to the massive spirals Milky Way (MW) and M31, various environmental processes have been invoked to explain their origin. In this context, the tidal stirring model postulates that interactions with MW-sized hosts can transform rotationally supported dwarfs, resembling present-day dwarf irregular (dIrr) galaxies, into systems with the kinematic and structural properties of dSphs. Using N-body+SPH simulations, we investigate the dependence of this transformation mechanism on the gas fraction, fgas, in the disk of the progenitor dwarf. Our numerical experiments incorporate for the first time the combined effects of radiative cooling, ram-pressure stripping, star formation, supernova (SN) winds, and a cosmic UV background. For a given orbit inside the primary galaxy, rotationally supported dwarfs with gas fractions akin to those of observed dIrrs (fgas >= 0.5), demonstrate a substantially enhanced likelihood and efficiency of transformation into dSphs relative to their collisionless (fgas = 0) counterparts. We argue that the combination of ram-pressure stripping and SN winds causes the gas-rich dwarfs to respond more impulsively to tides, augmenting their transformation. When fgas >= 0.5, disky dwarfs on previously unfavorable low-eccentricity or large-pericenter orbits are still able to transform. On the widest orbits, the transformation is incomplete; the dwarfs retain significant rotational support, a relatively flat shape, and some gas, naturally resembling transition-type systems. We conclude that tidal stirring constitutes a prevalent evolutionary mechanism for shaping the structure of dwarf galaxies within the currently favored CDM cosmological paradigm.

We present results from the "Ponos" simulation suite on the early evolution of a massive, $M_{\rm vir}(z=0)=1.2\times 10^{13}$ M$_{\odot}$ galaxy. At $z\gtrsim6$, before feedback from a central supermassive black hole becomes dominant, the main galaxy has a stellar mass $\sim 2\times 10^{9}$ M$_{\odot}$ and a star formation rate $\sim 20$ M$_{\odot}$ yr$^{-1}$. The galaxy sits near the expected main sequence of star-forming galaxies at those redshifts, and resembles moderately star-forming systems observed at $z>5$. The high specific star formation rate results in vigorous heating and stirring of the gas by supernovae feedback, and the galaxy develops a thick and turbulent disc, with gas velocity dispersion $\sim 40$ km s$^{-1}$, rotation to dispersion ratio $\sim 2$, and with a significant amount of gas at $\sim 10^5$ K. The Toomre parameter always exceeds the critical value for gravito-turbulence, $Q\sim 1.5-2$, mainly due to the contribution of warm/hot gas inside the disc. Without feedback, a nearly gravito-turbulent regime establishes with similar gas velocity dispersion and lower $Q$. We propose that the "hot and turbulent" disc regime seen in our simulations, unlike the "cold and turbulent" gravito-turbulent regime of massive clumpy disc galaxies at $z\sim 1-2$, is a fundamental characterisation of main sequence galaxies at $z\gtrsim 6$, as they can sustain star formation rates comparable to those of low-mass starbursts at $z=0$. This results in no sustained coherent gas inflows through the disc, and in fluctuating and anisotropic mass transport, possibly postponing the assembly of the bulge and causing the initial feeding of the central black hole to be highly intermittent.

We present an observational study of the effect of bars on the gas component and on the star formation properties of their host galaxies in a statistically significant sample of resolved objects, the $Herschel$ Reference Sample. The analysis of optical and far--infrared images allows us to identify a clear spatial correlation between stellar bars and the cold-gas distribution mapped by the warm dust emission. We find that the infrared counterparts of optically identified bars are either bar--like structures or dead central regions in which star formation is strongly suppressed. Similar morphologies are found in the distribution of star formation directly traced by H$\alpha$ maps. The sizes of such optical and infrared structures correlate remarkably well, hinting at a causal connection. In the light of previous observations and of theoretical investigations in the literature, we interpret our findings as further evidence of the scenario in which bars drive strong inflows toward their host nuclei: young bars are still in the process of perturbing the gas and star formation clearly delineates the shape of the bars; old bars on the contrary already removed any gas within their extents, carving a dead region of negligible star formation.

We study the dynamics of galaxy mergers, with emphasis on the gas feeding of nuclear regions, using a suite of hydrodynamical simulations of galaxy encounters. The high spatial and temporal resolution of the simulations allows us to not only recover the standard picture of tidal-torque induced inflows, but also to detail another, important feeding path produced by ram pressure. The induced shocks effectively decouple the dynamics of the gas from that of the stars, greatly enhancing the loss of gas angular momentum and leading to increased central inflows. The ram-pressure shocks also cause, in many cases, the entire galactic gas disc of the smaller galaxy to abruptly change its direction of rotation, causing a complete "flip" and, several $10^8$ yr later, a subsequent "counter-flip". This phenomenon results in the existence of long-lived decoupled gas-stellar and stellar-stellar discs, which could hint at a new explanation for the origin of some of the observed kinematically decoupled cores/counter-rotating discs. Lastly, we speculate, in the case of non-coplanar mergers, on the possible existence of a new class of remnant systems similar to some of the observed X-shaped radio galaxies.

We analyse the output of the hi-res cosmological zoom-in simulation ErisBH to study self-consistently the formation of a strong stellar bar in a Milky Way-type galaxy and its effect on the galactic structure, on the central gas distribution and on star formation. The simulation includes radiative cooling, star formation, SN feedback and a central massive black hole which is undergoing gas accretion and is heating the surroundings via thermal AGN feedback. A large central region in the ErisBH disk becomes bar-unstable after z~1.4, but a clear bar-like structure starts to grow significantly only after z~0.4, possibly triggered by the interaction with a massive satellite. At z~0.1 the bar reaches its maximum radial extent of l~2.2 kpc. As the bar grows, it becomes prone to buckling instability, which we quantify based on the anisotropy of the stellar velocity dispersion. The actual buckling event is observable at z~0.1, resulting in the formation of a boxy-peanut bulge clearly discernible in the edge-on view of the galaxy at z=0. The bar in ErisBH does not dissolve during the formation of the bulge but remains strongly non-axisymmetric down to the resolution limit of ~100 pc at z=0. During its early growth, the bar exerts a strong torque on the gas within its extent and drives gas inflows that enhance the nuclear star formation on sub-kpc scales. Later on the infalling gas is nearly all consumed into stars and, to a lesser extent, accreted onto the central black hole, leaving behind a gas-depleted region within the central ~2 kpc. Observations would more likely identify a prominent, large-scale bar at the stage when the galactic central region has already been quenched. Bar-driven quenching may play an important role in disk-dominated galaxies at all redshift. [Abridged]

Massive black hole binaries (BHBs) are expected to be one of the most powerful sources of gravitational waves (GWs) in the frequency range of the pulsar timing array and of forthcoming space-borne detectors. They are believed to form in the final stages of galaxy mergers, and then harden by slingshot ejections of passing stars. However, evolution via the slingshot mechanism may be ineffective if the reservoir of interacting stars is not readily replenished, and the binary shrinking may come to a halt at roughly a parsec separation. Recent simulations suggest that the departure from spherical symmetry, naturally produced in merger remnants, leads to efficient loss cone refilling, preventing the binary from stalling. However, current N-body simulations able to accurately follow the evolution of BHBs are limited to very modest particle numbers. Brownian motion may artificially enhance the loss cone refilling rate in low-N simulations, where the binary encounters a larger population of stars due its random motion. Here we study the significance of Brownian motion of BHBs in merger remnants in the context of the final parsec problem. We simulate mergers with various particle numbers (from 8k to 1M) and with several density profiles. Moreover, we compare simulations where the BHB is fixed at the centre of the merger remnant with simulations where the BHB is free to random walk. We find that Brownian motion does not significantly affect the evolution of BHBs in simulations with particle numbers in excess of one million, and that the hardening measured in merger simulations is due to collisionless loss cone refilling.

We present an automatic procedure to perform reliable photometry of galaxies on SDSS images. We selected a sample of 5853 galaxies in the Coma and Virgo superclusters. For each galaxy, we derive Petrosian g and i magnitudes, surface brightness profiles and color profiles. Unlike the SDSS pipeline, our procedure is not affected by the well known shredding problem and efficiently extracts Petrosian magnitudes for all galaxies. Hence we derived magnitudes even from the population of galaxies missed by the SDSS that represents 25% of all Local supercluster galaxies and ~95% of galaxies with g < 11 mag. After correcting the g and i magnitudes for Galactic and internal extinction, the blue and red sequences in the color magnitude diagram are well separated, with similar slopes. In addition, we study (i) the color-magnitude diagrams in different galaxy regions, the inner (r <= 1 kpc), intermediate (0.2R{Pet} <=r<= 0.3R{Pet}) and outer, disk-dominated (r>=0.35R{Pet})) zone, and (ii), we compute template color profiles, discussing the dependences of the templates on the galaxy masses and on their morphological type. The two analyses consistently lead to a picture where elliptical galaxies show no color gradients, irrespective of their masses. Spirals, instead, display a steeper gradient in their color profiles with increasing mass, consistently with the growing relevance of a bulge and/or a bar component above 10^{10} Modot.

We explore the hardening of a massive black hole binary embedded in a circum-binary gas disc when the binary and the gas are coplanar and the gas is counter-rotating. The secondary black hole, revolving in the direction opposite to the gas, experiences a drag from gas-dynamical friction and from direct accretion of part of it. Using two-dimensional (2D) hydrodynamical grid simulations we investigate the effect of changing the accretion prescriptions on the dynamics of the secondary black hole which in turn affect the binary hardening and eccentricity evolution. We find that realistic accretion prescriptions lead to results that differ from those inferred assuming accretion of all the gas within the Roche Lobe of the secondary black hole. Different accretion prescriptions result in different disc's surface densities which alter the black hole's dynamics back. Full 3D SPH realizations of a number of representative cases, run over a shorter interval of time, validate the general trends observed in the less computationally demanding 2D simulations. Initially circular black hole binaries increase only slightly their eccentricity which then oscillates around small values (<0.1) while they harden. By contrast, initially eccentric binaries become more and more eccentric. A semi-analytical model describing the black hole's dynamics under accretion only explores the late evolution stages of the binary in an otherwise unperturbed retrograde disc to illustrate how eccentricity evolves with time in relation to the shape of the underlying surface density distribution.

Approximately 20% of early-type galaxies host small nuclear stellar discs that are tens to a few hundred parsecs in size. Such discs are expected to be easily disrupted during major galactic encounters, hence their age serve to constrain their assembly history. We use VIMOS integral-field spectroscopic observations for the intermediate-mass E0 galaxy NGC4458 and age-date its nuclear disc via high-resolution fitting of various model spectra. We find that the nuclear disc is at least 6 Gyr old. A clue to gain narrow limits to the stellar age is our knowledge of the nuclear disc contribution to the central surface brightness. The presence of an old nuclear disk, or the absence of disruptive encounters since z~0.6, for a small galaxy such as NGC4458 which belongs to the Virgo cluster, may be consistent with a hierarchical picture for galaxy formation where the smallest galaxies assembles earlier and the crowded galactic environments reduce the incidence of galaxy mergers. On the other hand, NGC4458 displays little or no bulk rotation except for a central kpc-scale kinematically-decoupled core. Slow rotation and decoupled core are usually explained in terms of mergers. The presence and age of the nuclear disc constraint these mergers to have happened at high redshift.

We searched for quasi-periodicities on year-like timescales in the light curves of 6 blazars in the optical - near infrared bands and we made a comparison with the high energy emission. We obtained optical/NIR light curves from REM photometry plus archival SMARTS data and we accessed the Fermi light curves for the $\gamma$-ray data. The periodograms often show strong peaks in the optical and gamma-ray bands, which in some cases may be inter-related. The significance of the revealed peaks is then discussed, taking into account that the noise is frequency dependent. Quasi-periodicities on a year-like timescale appear to occur often in blazars. No straightforward model describing these possible periodicities is yet available, but some plausible interpretations for the physical mechanisms causing periodic variabilities of these sources are examined.

The rapid assembly of the massive black holes that power the luminous quasars observed at $z \sim 6-7$ remains a puzzle. Various direct collapse models have been proposed to head-start black hole growth from initial seeds with masses $\sim 10^5\,\rm M_\odot$, which can then reach a billion solar mass while accreting at the Eddington limit. Here we propose an alternative scenario based on radiatively inefficient super-critical accretion of stellar-mass holes embedded in the gaseous circum-nuclear discs (CNDs) expected to exist in the cores of high redshift galaxies. Our sub-pc resolution hydrodynamical simulations show that stellar-mass holes orbiting within the central 100 pc of the CND bind to very high density gas clumps that arise from the fragmentation of the surrounding gas. Owing to the large reservoir of dense cold gas available, a stellar-mass black hole allowed to grow at super-Eddington rates according to the "slim disc" solution can increase its mass by 3 orders of magnitudes within a few million years. These findings are supported by simulations run with two different hydro codes, RAMSES based on the Adaptive Mesh Refinement technique and GIZMO based on a new Lagrangian Godunov-type method, and with similar, but not identical, sub-grid recipes for star formation, supernova feedback, black hole accretion and feedback. The low radiative efficiency of super-critical accretion flows are instrumental to the rapid mass growth of our black holes, as they imply modest radiative heating of the surrounding nuclear environment.

The dynamics of massive black holes (BHs) in galaxy mergers is a rich field of research that has seen much progress in recent years. In this contribution we briefly review the processes describing the journey of BHs during mergers, from the cosmic context all the way to when BHs coalesce. If two galaxies each hosting a central BH merge, the BHs would be dragged towards the center of the newly formed galaxy. If/when the holes get sufficiently close, they coalesce via the emission of gravitational waves. How often two BHs are involved in galaxy mergers depends crucially on how many galaxies host BHs and on the galaxy merger history. It is therefore necessary to start with full cosmological models including BH physics and a careful dynamical treatment. After galaxies have merged, however, the BHs still have a long journey until they touch and coalesce. Their dynamical evolution is radically different in gas-rich and gas-poor galaxies, leading to a sort of "dichotomy" between high-redshift and low-redshift galaxies, and late-type and early-type, typically more massive galaxies.

Stellar bars are a common feature in massive disc galaxies. On a theoretical ground, the response of gas to a bar is generally thought to cause nuclear starbursts and, possibly, AGN activity once the perturbed gas reaches the central super-massive black hole. By means of high resolution numerical simulations we detail the purely dynamical effects that a forming bar exerts on the gas of an isolated disc galaxy. The galaxy is initially unstable to the formation of non-axisymmetric structures, and within 1 Gyr it develops spiral arms that eventually evolve into a central stellar bar on kpc scale. A first major episode of gas inflow occurs during the formation of the spiral arms while at later times, when the stellar bar is establishing, a low density region is carved between the bar co-rotational and inner Lindblad resonance radii. The development of such "dead zone" inhibits further massive gas inflows. Indeed, the gas inflow reaches its maximum during the relatively fast bar formation phase and not, as often assumed, when the bar is fully formed. We conclude that the low efficiency of long-lived, evolved bars in driving gas toward galactic nuclei is the reason why observational studies have failed to establish an indisputable link between bars and AGNs. On the other hand, the high efficiency in driving strong gas inflows of the intrinsically transient process of bar formation suggests that the importance of bars as drivers of AGN activity in disc galaxies has been overlooked so far. We finally prove that our conclusions are robust against different numerical implementations of the hydrodynamics routinely used in galaxy evolution studies.

Mergers of gas-rich galaxies are key events in the hierarchical built-up of cosmic structures, and can lead to the formation of massive black hole binaries. By means of high-resolution hydrodynamical simulations we consider the late stages of a gas-rich major merger, detailing the dynamics of two circumnuclear discs, and of the hosted massive black holes during their pairing phase. During the merger gas clumps with masses of a fraction of the black hole mass form because of fragmentation. Such high-density gas is very effective in forming stars, and the most massive clumps can substantially perturb the black hole orbits. After $\sim 10$ Myr from the start of the merger a gravitationally bound black hole binary forms at a separation of a few parsecs, and soon after, the separation falls below our resolution limit of $0.39$ pc. At the time of binary formation the original discs are almost completely disrupted because of SNa feedback, while on pc scales the residual gas settles in a circumbinary disc with mass $\sim 10^5 M_\odot$. We also test that binary dynamics is robust against the details of the SNa feedback employed in the simulations, while gas dynamics is not. We finally highlight the importance of the SNa time-scale on our results.

The main challenge for understanding the fuelling of supermassive black holes in active galactic nuclei is not to account for the source of fuel, but rather to explain its delivery from the boundaries of the black hole sphere of influence (10-100 pc) down to sub-parsec scales. In this work, we report on a series of numerical experiments aimed at exploring in further depth our model of "overlapping inflow events" as catalysts for rapid accretion, seeding a turbulent field in the infalling gas. We initially set a gaseous shell in non-equilibrium rotation around a supermassive black hole. After infall, the shell stalls in a disc-like structure. A second shell is then set in either co-rotation or counter-rotation with respect to the first and is let to impinge on the previously-formed disc. We find that combined turbulence and overlap significantly enhance accretion in counter-rotating inflows, while turbulence dominates for co-rotating inflows. The leftovers of overlapping inflows are warped nuclear discs, whose morphology depends on the relative orientation and angular momentum of the disc and the shell. Overlapping inflows leave observational signatures in the gas rotation curves.

Nuclear stellar discs (NSDs) can help to constrain the assembly history of their host galaxies, as long as we can assume them to be fragile structures that are disrupted during merger events. In this work we investigate the fragility of NSDs by means of N-body simulations reproducing the last phases of a galaxy encounter, when the nuclear regions of the two galaxies merge. For this, we exposed a NSD set in the gravitational potential of the bulge and supermassive black hole of a primary galaxy to the impact of the supermassive black hole from a secondary galaxy. We explored merger events of different mass ratios, from major mergers with a 1:1 mass ratio to intermediate and minor interactions with 1:5 and 1:10 ratios, while considering various impact geometries. We analyse the end results of such mergers from different viewing angles and looked for possible photometric and kinematic signatures of the presence of a disc in the remnant surface density and velocity maps, while adopting detection limits from real observations. Our simulations show that indeed NSDs are fragile against major mergers, which leave little trace of NSDs both in images and velocity maps, while signatures of a disc can be found in the majority of the intermediate to minor-merger remnants and in particular when looking at their kinematics. These results show that NSDs could allow to distinguish between these two modes of galaxy assembly, which may indeed pertain to different kinds of galaxies or galactic environments.

A growing body of evidence indicates that the star formation rate per unit stellar mass (sSFR) decreases with increasing mass in normal "main-sequence" star forming galaxies. Many processes have been advocated as responsible for such a trend (also known as mass quenching), e.g., feedback from active galactic nuclei (AGNs), and the formation of classical bulges. We determine a refined star formation versus stellar mass relation in the local Universe. To this aim we use the Halpha narrow-band imaging follow-up survey (Halpha3) of field galaxies selected from the HI Arecibo Legacy Fast ALFA Survey (ALFALFA) in the Coma and Local superclusters. By complementing this local determination with high-redshift measurements from the literature, we reconstruct the star formation history of main-sequence galaxies as a function of stellar mass from the present epoch up to z=3. In agreement with previous studies, our analysis shows that quenching mechanisms occur above a threshold stellar mass M_knee that evolves with redshift as propto (1+z)^{2}. Moreover, visual morphological classification of individual objects in our local sample reveals a sharp increase in the fraction of visually-classified strong bars with mass, hinting that strong bars may contribute to the observed downturn in the sSFR above M_knee. We test this hypothesis using a simple but physically-motivated numerical model for bar formation, finding that strong bars can rapidly quench star formation in the central few kpc of field galaxies. We conclude that strong bars contribute significantly to the red colors observed in the inner parts of massive galaxies, although additional mechanisms are likely required to quench the star formation in the outer regions of massive spiral galaxies. Intriguingly, when we extrapolate our model to higher redshifts, we successfully recover the observed redshift evolution for M_knee.

We use a suite of hydrodynamical simulations of galaxy mergers to compare star formation rate (SFR) and black hole accretion rate (BHAR) for galaxies before the interaction ('stochastic' phase), during the `merger' proper, lasting ~0.2-0.3 Gyr, and in the `remnant' phase. We calculate the bi-variate distribution of SFR and BHAR and define the regions in the SFR-BHAR plane that the three phases occupy. No strong correlation between BHAR and galaxy-wide SFR is found. A possible exception are galaxies with the highest SFR and the highest BHAR. We also bin the data in the same way used in several observational studies, by either measuring the mean SFR for AGN in different luminosity bins, or the mean BHAR for galaxies in bins of SFR. We find that the apparent contradiction or SFR versus BHAR for observed samples of AGN and star forming galaxies is actually caused by binning effects. The two types of samples use different projections of the full bi-variate distribution, and the full information would lead to unambiguous interpretation. We also find that a galaxy can be classified as AGN-dominated up to 1.5 Gyr after the merger-driven starburst took place. Our study is consistent with the suggestion that most low-luminosity AGN hosts do not show morphological disturbances.

We present a new suite of hydrodynamical simulations and use it to study, in detail, black hole and galaxy properties. The high time, spatial and mass resolution, and realistic orbits and mass ratios, down to 1:6 and 1:10, enable us to meaningfully compare star formation rate (SFR) and BH accretion rate (BHAR) timescales, temporal behaviour and relative magnitude. We find that (i) BHAR and galaxy-wide SFR are typically temporally uncorrelated, and have different variability timescales, except during the merger proper, lasting ~0.2-0.3 Gyr. BHAR and nuclear (<100 pc) SFR are better correlated, and their variability are similar. Averaging over time, the merger phase leads typically to an increase by a factor of a few in the BHAR/SFR ratio. (ii) BHAR and nuclear SFR are intrinsically proportional, but the correlation lessens if the long-term SFR is measured. (iii) Galaxies in the remnant phase are the ones most likely to be selected as systems dominated by an active galactic nucleus (AGN), because of the long time spent in this phase. (iv) The timescale over which a given diagnostic probes the SFR has a profound impact on the recovered correlations with BHAR, and on the interpretation of observational data.

Massive black hole binaries are naturally predicted in the context of the hierarchical model of structure formation. The binaries that manage to lose most of their angular momentum can coalesce to form a single remnant. In the last stages of this process, the holes undergo an extremely loud phase of gravitational wave emission, possibly detectable by current and future probes. The theoretical effort towards obtaining a coherent physical picture of the binary path down to coalescence is still underway. In this paper, for the first time, we take advantage of observational studies of active galactic nuclei evolution to constrain the efficiency of gas-driven binary decay. Under conservative assumptions we find that gas accretion toward the nuclear black holes can efficiently lead binaries of any mass forming at high redshift (> 2) to coalescence within the current time. The observed "downsizing" trend of the accreting black hole luminosity function further implies that the gas inflow is sufficient to drive light black holes down to coalescence, even if they bind in binaries at lower redshifts, down to z~0.5 for binaries of ~10 million solar masses, and z~0.2 for binaries of ~1 million solar masses. This has strong implications for the detection rates of coalescing black hole binaries of future space-based gravitational wave experiments.

Numerical effects are known to plague adaptive mesh refinement (AMR) codes when treating massive particles, e.g. representing massive black holes (MBHs). In an evolving background, they can experience strong, spurious perturbations and then follow unphysical orbits. We study by means of numerical simulations the dynamical evolution of a pair MBHs in the rapidly and violently evolving gaseous and stellar background that follows a galaxy major merger. We confirm that spurious numerical effects alter the MBH orbits in AMR simulations, and show that numerical issues are ultimately due to a drop in the spatial resolution during the simulation, drastically reducing the accuracy in the gravitational force computation. We therefore propose a new refinement criterion suited for massive particles, able to solve in a fast and precise way for their orbits in highly dynamical backgrounds. The new refinement criterion we designed enforces the region around each massive particle to remain at the maximum resolution allowed, independently upon the local gas density. Such maximally-resolved regions then follow the MBHs along their orbits, and effectively avoids all spurious effects caused by resolution changes. Our suite of high resolution, adaptive mesh-refinement hydrodynamic simulations, including different prescriptions for the sub-grid gas physics, shows that the new refinement implementation has the advantage of not altering the physical evolution of the MBHs, accounting for all the non trivial physical processes taking place in violent dynamical scenarios, such as the final stages of a galaxy major merger.

We study supermassive black holes (BHs) in merging galaxies, using a suite of hydrodynamical simulations with very high spatial (~10 pc) and temporal (~1 Myr) resolution, where we vary the initial mass ratio, the orbital configuration, and the gas fraction. (i) We address the question of when and why, during a merger, increased BH accretion occurs, quantifying gas inflows and BH accretion rates. (ii) We also quantify the relative effectiveness in inducing AGN activity of merger-related versus secular-related causes, by studying different stages of the encounter: the stochastic (or early) stage, the (proper) merger stage, and the remnant (or late) stage. (iii) We assess which galaxy mergers preferentially enhance BH accretion, finding that the initial mass ratio is the most important factor. (iv) We study the evolution of the BH masses, finding that the BH mass contrast tends to decrease in minor mergers and to increase in major mergers. This effect hints at the existence of a preferential range of mass ratios for BHs in the final pairing stages. (v) In both merging and dynamically quiescent galaxies, the gas accreted by the BH is not necessarily the gas with $low$ angular momentum, but the gas that $loses$ angular momentum.

The peculiar QSO J0927+2943 shows multiple sets of emission lines in its optical spectrum. This signature has been interpreted as the relative motion between a black hole, either recoiling or bound in a binary system, and its host galaxy, or as a superposition of two galaxies along the line of sight. In order to test these scenarios, we have collected 2mm CO(2-1) observations using the IRAM Plateau de Bure Interferometer, and optical images and spectroscopy at the Calar Alto observatory. Together with archival HST images, these data provide unique insights on the nature of this system. The recoiling/binary black hole scenarios are ruled out by the clear detection of a galactic-scale molecular gas reservoir at the same redshift of the QSO broad lines. The observations presented here also disfavour the superposition model, although with less constraints. Thus, the origin of the second, bright set of narrow emission lines in J0927+2943 is still unknown.

We present the results of a semianalytical model that evolves the masses and spins of massive black holes together with the properties of their host galaxies along the cosmic history. As a consistency check, our model broadly reproduces a number of observations, e.g. the cosmic star formation history, the black hole mass and luminosity function and the galaxy mass function at low redshift, the black hole to bulge mass relation, and the morphological distribution at low redshift. For the first time in a semianalytical investigation, we relax the simplifying assumptions of perfect coherency or perfect isotropy of the gas fueling the black holes. The dynamics of gas is instead linked to the morphological properties of the host galaxies, resulting in different spin distributions for black holes hosted in different galaxy types. We compare our results with the observed sample of spin measurements obtained through broad K$\alpha$ iron line fitting. The observational data disfavor both accretion along a fixed direction and isotropic fueling. Conversely, when the properties of the accretion flow are anchored to the kinematics of the host galaxy, we obtain a good match between theoretical expectations and observations. A mixture of coherent accretion and phases of activity in which the gas dynamics is similar to that of the stars in bulges (i.e., with a significant velocity dispersion superimposed to a net rotation) best describes the data, adding further evidence in support to the coevolution of massive black holes and their hosts.

Recoiling supermassive black holes (SMBHs) are considered one plausible physical mechanism to explain high velocity shifts between narrow and broad emission lines sometimes observed in quasar spectra. If the sphere of influence of the recoiling SMBH is such that only the accretion disc is bound, the dusty torus would be left behind, hence the SED should then present distinctive features (i.e. a mid-infrared deficit). Here we present results from fitting the Spectral Energy Distributions (SEDs) of 32 Type-1 AGN with high velocity shifts between broad and narrow lines. The aim is to find peculiar properties in the multi-wavelength SEDs of such objects by comparing their physical parameters (torus and disc luminosity, intrinsic reddening, and size of the 12$\mu$m emitter) with those estimated from a control sample of $\sim1000$ \emph{typical} quasars selected from the Sloan Digital Sky Survey in the same redshift range. We find that all sources, with the possible exception of J1154+0134, analysed here present a significant amount of 12~$\mu$m emission. This is in contrast with a scenario of a SMBH displaced from the center of the galaxy, as expected for an undergoing recoil event.

We consider super-critical accretion with angular momentum onto stellar-mass black holes as a possible mechanism for growing billion-solar-mass holes from light seeds at early times. We use the radiatively-inefficient "slim disk" solution -- advective, optically thick flows that generalize the standard geometrically thin disk model -- to show how mildly super-Eddington intermittent accretion may significantly ease the problem of assembling the first massive black holes when the Universe was less than 0.8 Gyr old. Because of the low radiative efficiencies of slim disks around non-rotating as well as rapidly rotating holes, the mass e-folding timescale in this regime is nearly independent of the spin parameter. The conditions that may lead to super-critical growth in the early Universe are briefly discussed.

We study the dynamical evolution of supermassive black holes (BHs) in merging galaxies on scales of hundreds of kpc to 10 pc, to identify the physical processes that aid or hinder the orbital decay of BHs. We present hydrodynamical simulations of galaxy mergers with a resolution of $\leq$20 pc, chosen to accurately track the motion of the nuclei and provide a realistic environment for the evolution of the BHs. We find that, during the late stages of the merger, tidal shocks inject energy in the nuclei, causing one or both nuclei to be disrupted and leaving their BH `naked', without any bound gas or stars. In many cases, the nucleus that is ultimately disrupted is that of the larger galaxy (`nuclear coup'), as star formation grows a denser nuclear cusp in the smaller galaxy. We supplement our simulations with an analytical estimate of the orbital-decay time required for the BHs to form a binary at unresolved scales, due to dynamical friction. We find that, when a nuclear coup occurs, the time-scale is much shorter than when the secondary's nucleus is disrupted, as the infalling BH is more massive, and it also finds itself in a denser stellar environment.

One of the greatest issues in modelling black hole fuelling is our lack of understanding of the processes by which gas loses angular momentum and falls from galactic scales down to the nuclear region where an accretion disc forms, subsequently guiding the inflow of gas down to the black hole horizon. It is feared that gas at larger scales might still retain enough angular momentum and settle into a larger scale disc with very low or no inflow to form or replenish the inner accretion disc (on ~0.01 pc scales). In this paper we report on hydrodynamical simulations of rotating infalling gas shells impacting at different angles onto a pre-existing, primitive large scale (~10 pc) disc around a super-massive black hole. The aim is to explore how the interaction between the shell and the disc redistributes the angular momentum on scales close to the black hole's sphere of influence. Angular momentum redistribution via hydrodynamical shocks leads to inflows of gas across the inner boundary, enhancing the inflow rate by more than 2-3 orders of magnitude. In all cases, the gas inflow rate across the inner parsec is higher than in the absence of the interaction, and the orientation of the angular momentum of the flow in the region changes with time due to gas mixing. Warped discs or nested misaligned rings form depending on the angular momentum content of the infalling shell relative to the disc. In the cases in which the shell falls in near counter-rotation, part of the resulting flows settle into an inner dense disc which becomes more susceptible to mass transfer.

We obtained long-slit optical spectra of the nuclear regions of 376 galaxies in the local Universe using the 1.5m Cassini telescope of Bologna Observatory. Of these spectra, 164 were either never taken before by the Sloan Digital Sky Survey (SDSS), or given by the Nasa Extragalactic Database (NED). With these new spectra, we contribute investigating the occurrence of active galactic nuclei (AGNs). Nevertheless, we stress that the present sample is by no means complete, thus, it cannot be used to perform any demographic study. Following the method presented in Gavazzi et al. (2011), we classify the nuclear spectra using a six bin scheme: SEY (Seyfert), sAGN (strong AGN), and wAGN (weak AGN) represent active galactic nuclei of different levels of activity; HII accounts for star-forming nuclei; RET (retired) and PAS (passive) refer to nuclei with poor or no star-formation activity. The spectral classification is performed using the ratio of 6584 {\lambda} [NII] to H{\alpha} lines and the equivalent width (EW) of H{\alpha} versus [NII]/H{\alpha} (WHAN diagnostic introduced by Cid Fernandes and collaborators) after correcting H{\alpha} for underlying absorption. The obtained spectra are made available in machine readable format via the Strasbourg Astronomical Data Center (CDS) and NED.

Theoretically, bound binaries of massive black holes are expected as the natural outcome of mergers of massive galaxies. From the observational side, however, massive black hole binaries remain elusive. Velocity shifts between narrow and broad emission lines in quasar spectra are considered a promising observational tool to search for spatially unresolved, dynamically bound binaries. In this series of papers we investigate the nature of such candidates through analyses of their spectra, images and multi-wavelength spectral energy distributions. Here we investigate the properties of the optical spectra, including the evolution of the broad line profiles, of all the sources identified in our previous study. We find a diverse phenomenology of broad and narrow line luminosities, widths, shapes, ionization conditions and time variability, which we can broadly ascribe to 4 classes based on the shape of the broad line profiles: 1) Objects with bell-shaped broad lines with big velocity shifts (>1000 km/s) compared to their narrow lines show a variety of broad line widths and luminosities, modest flux variations over a few years, and no significant change in the broad line peak wavelength. 2) Objects with double-peaked broad emission lines tend to show very luminous and broadened lines, and little time variability. 3) Objects with asymmetric broad emission lines show a broad range of broad line luminosities and significant variability of the line profiles. 4) The remaining sources tend to show moderate to low broad line luminosities, and can be ascribed to diverse phenomena. We discuss the implications of our findings in the context of massive black hole binary searches.

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

We present the analysis of Halpha3, an Halpha imaging survey of galaxies selected from the HI ALFALFA Survey in the Coma Supercluster. By using the Halpha line as a tracer of the "instantaneous" star formation, complemented with optical colors from SDSS we explore the hypothesis that a morphological sequence of galaxies of progressively earlier type, lower gas-content exists in the neighborhood of the Coma cluster, with specific star formation activity decreasing with increasing local galaxy density and velocity dispersion. In the dwarf regime (8.5<\log(M*)< 9.5) we identify a 4-step sequence of galaxies with progressively redder colors, i.e. of decreasing specific star formation, from (1) HI-rich Late-Type Galaxies belonging to the blue-cloud exhibiting extended plus nuclear star formation, (2) HI-poor LTGs with nuclear star formation only, (3) HI-poor galaxies with no star formation either extended or nuclear, but with nuclear Post-Star-Burst signature,(4) Early-type Galaxies in the red-sequence, with no gas or star formation on all scales. Along this sequence the quenching of the star formation proceeds radially outside-in. The progression toward redder colors found along this "morphological" (gas content) sequence is comparable to the one obtained increasing the local galaxy density, from the cosmic filaments (1,2), to the rich clusters (2,3,4). In the dwarf regime we find evidence for evolution of HI-rich LTGs into ETGs, passing through HI-poor and PSB galaxies, driven by the environment. We identify ram-pressure as the mechanism most likely responsible for this transformation. We conclude that infall of low-mass galaxies has proceeded for the last 7.5 Gyr building up the Coma cluster at a rate of approximately 100 galaxies per Gyr.

Mass and spin are the only two parameters needed to completely characterize black holes in General Relativity. However, the interaction between black holes and their environment is where complexity lies, as the relevant physical processes occur over a large range of scales. That is particularly relevant in the case of super-massive black holes (SMBHs), hosted in galaxy centers, and surrounded by swirling gas and various generations of stars. These compete with the SMBH for gas consumption and affect both dynamics and thermodynamics of the gas itself. How the behavior of such fiery environment influence the angular momentum of the gas accreted onto SMBHs, and, hence, black-hole spins is uncertain. We explore the interaction between SMBHs and their environment via first 3D sub-parsec resolution simulations (ranging from 0.1 pc to 1 kpc scales) that study the evolution of the SMBH spin by including the effects of star formation, stellar feedback, radiative transfer, and metal pollution according to the proper stellar yields and lifetimes. This approach is crucial to investigate the impact of star formation processes and feedback effects on the angular momentum of the material that could accrete on the central hole. We find that star formation and feedback mechanisms can locally inject significant amounts of entropy in the surrounding medium, and impact on the inflow inclination angles and Eddington fractions. As a consequence, the resulting trends show upper-intermediate equilibrium values for the spin parameter, a, of about 0.6 - 0.9, corresponding to radiative efficiencies \epsilon = 9% - 15%. These results suggest that star formation feedback taking place in the circum-nuclear disk during the in-fall cannot induce alone very strong chaotic trends in the gas flow, quite independently from the different numerical parameters.

We study the link between the X-ray emission in radio-quiet AGNs and the accretion rate on the central Supermassive Black-Hole (SMBH) using a well-defined and statistically complete sample of 70 type1 AGNs extracted from the XMM-Newton Bright Serendipitous survey (XBS). To this end, we search and quantify the statistical correlations between the main parameters that characterize the X-ray emission (i.e. the X-ray spectral slope and the X-ray loudness), and the accretion rate, both absolute and relative to the Eddington limit (Eddington ratio). Here, we summarize and discuss the main statistical correlations found and their possible implications on current disk-corona models.

We discuss the relationship between a standard Shakura & Sunyaev (1973) accretion disk model and the Big Blue Bump (BBB) observed in Type 1 AGN, and propose a new method to estimate black hole masses. We apply this method to a sample of 23 radio-loud narrow-line Seyfert 1 (RL-NLS1) galaxies, using data from WISE (Wide-field Infrared Survey Explorer), SDSS (Sloan Digital Sky Survey) and GALEX. Our black hole mass estimates are at least a factor $\sim$6 above previous results based on single epoch virial methods, while the Eddington ratios are correspondingly lower. Hence, the black hole masses of RL-NLS1 galaxies are typically above $10^8 M_{\sun}$, in agreement with the typical black hole mass of blazars.

Massive black holes in galactic nuclei vary their mass M and spin vector J due to accretion. In this study we relax, for the first time, the assumption that accretion can be either chaotic, i.e. when the accretion episodes are randomly and isotropically oriented, or coherent, i.e. when they occur all in a preferred plane. Instead, we consider different degrees of anisotropy in the fueling, never confining to accretion events on a fixed direction. We follow the black hole growth evolving contemporarily mass, spin modulus a and spin direction. We discover the occurrence of two regimes. An early phase (M <~ 10 million solar masses) in which rapid alignment of the black hole spin direction to the disk angular momentum in each single episode leads to erratic changes in the black hole spin orientation and at the same time to large spins (a ~ 0.8). A second phase starts when the black hole mass increases above >~ 10 million solar masses and the accretion disks carry less mass and angular momentum relatively to the hole. In the absence of a preferential direction the black holes tend to spin-down in this phase. However, when a modest degree of anisotropy in the fueling process (still far from being coherent) is present, the black hole spin can increase up to a ~ 1 for very massive black holes (M >~ 100 million solar masses), and its direction is stable over the many accretion cycles. We discuss the implications that our results have in the realm of the observations of black hole spin and jet orientations.

We study the interplay between gas accretion and gravity torques in changing a binary elements and its total angular momentum (L) budget. Especially, we analyse the physical origin of the gravity torques (T_g) and their location within the disc. We analyse 3D SPH simulations of the evolution of initially quasi-circular massive black hole binaries (BHBs) residing in the central hollow of massive self-gravitating circumbinary discs. We use different thermodynamics within the cavity and for the numerical size of the black holes to show that (i) the BHB eccentricity growth found previously is a general result, independent of the accretion and the adopted thermodynamics; (ii) the semi-major axis decay depends both on the T_g and on the interplay with the disc-binary L-transfer due to accretion; (iii) the spectral structure of the T_g is predominately caused by disc edge overdensities and spiral arms developing in the body of the disc and, in general, does not reflect directly the period of the binary; (iv) the net T_g changes sign across the BHB corotation radius. We quantify the relative importance of the two, which appear to depend on the thermodynamical properties of the instreaming gas, and which is crucial in assessing the disc-binary L-transfer; (v) the net torque manifests as a purely kinematic (non-resonant) effect as it stems from the cavity, where the material flows in and out in highly eccentric orbits. Both accretion onto the black holes and the interaction with gas streams inside the cavity must be taken into account to assess the fate of the BHB. Moreover, the total torque exerted by the disc affects L(BHB) by changing all the elements (mass, mass ratio, eccentricity, semimajor axis) of the BHB. Common prescriptions equating tidal torque to semi-major axis shrinking might therefore be poor approximations for real astrophysical systems.

Massive black hole binaries at sub-parsec separations may display in their spectra anomalously small flux ratios between the MgII and CIV broad emission lines, i.e. F_MgII/F_CIV <~ 0.1, due to the erosion of the broad line region around the active, secondary black hole, by the tidal field of the primary. In Paper I by Montuori et al. (2011), we focussed on broad lines emitted by gas bound to the lighter accreting member of a binary when the binary is at the center of a hollow density region (the gap) inside a circum-binary disc. The main aim of this new study is at exploring the potential contribution to the broad line emission by the circum-binary disc and by gaseous streams flowing toward the black hole through the gap. We carry out a post-process analysis of data extracted from a SPH simulation of a circum-binary disc around a black hole binary. Our main result is that the MgII to CIV flux ratio can be reduced to ~ 0.1 within an interval of sub-pc binary separations of the order of a ~ (0.01-0.2)(f_Edd/0.1)^(1/2) pc corresponding to orbital periods of ~ (20-200) (f_Edd/0.1)^(3/4) years for a secondary BH mass in the range M_2 ~ 10^7-10^9 M_sun and a binary mass ratio of 0.3. At even closer separations this ratio returns to increase to values that are indistinguishable from the case of a single AGN (typically F_MgII/F_CIV ~ 0.3-0.4) because of the contribution to the MgII line from gas in the circum-binary disc.

B2 0954+25A, detected by the {\it Fermi} satellite, is a blazar with interesting observational properties: it has been observed to transit from a jet dominated to a disk dominated state; its radio spectrum appears flat at all observing frequencies (down to 74 MHz); optically, the H$\beta$ line profile is asymmetric. The flatness of radio spectrum suggests that the isotropic emission from radio lobes is very weak, despite the large size of its jet ($\gtrsim$ 500 kpc). Its broad--band spectral energy distribution is surprisingly similar to that of the prototypical $\gamma$--ray, radio loud, Narrow Line Seyfert 1 ($\gamma$--NLS1) galaxy PMN J0948+0022. In this work we revisit the mass estimates of B2 0954+25A considering only the symmetric component of the H$\beta$ line and find (1--3) $\times 10^8$ M$_{\sun}$. In light of our composite analysis, we propose to classify the source as a transition object between the class of Flat Spectrum Radio Quasar and $\gamma$--ray, radio loud NLS1. A comparison with two members of each class (3C 273 and PMN J0948+0022) is discussed.

We explore the newly discovered "hangup-kick" effect, which greatly amplifies the recoil for configuration with partial spin- orbital-angular momentum alignment, by studying a set of 48 new simulations of equal-mass, spinning black-hole binaries. We propose a phenomenological model for the recoil that takes this new effect into account and then use this model, in conjunction with statistical distributions for the spin magnitude and orientations, based on accretion simulations, to find the probabilities for observing recoils of several thousand km/s. In addition, we provide initial parameters, eccentricities, radiated linear and angular momentum, precession rates and remnant mass, spin, and recoils for all 48 configurations. Our results indicate that surveys exploring peculiar (redshifted or blueshifted) differential line-of-sight velocities should observe at least one case above 2000 km/s out of four thousand merged galaxies. The probability that a remnant BH receives a total recoil exceeding the ~2000 km/s escape velocity of large elliptical galaxies is ten times larger. Probabilities of recoils exceeding the escape velocity quickly rise to 5% for galaxies with escape velocities of 1000 km/s and nearly 20% for galaxies with escape velocities of 500 km/s. In addition the direction of these large recoils is strongly peaked toward the angular momentum axis, with very low probabilities of recoils exceeding 350 km/s for angles larger than 45 deg. with respect to the orbital angular momentum axis.

Supermassive black holes (SMBHs) have been detected in the centers of most nearby massive galaxies. Galaxies today are the products of billions of years of galaxy mergers, but also billions of years of SMBH activity as active galactic nuclei (AGNs) that is connected to galaxy mergers. In this context, detection of AGN pairs should be relatively common. Observationally, however, dual AGN are scant, being just a few percent of all AGN. In this Letter we investigate the triggering of AGN activity in merging galaxies via a suite of high resolution hydrodynamical simulations. We follow the dynamics and accretion onto the SMBHs as they move from separations of tens of kiloparsecs to tens of parsecs. Our resolution, cooling and star formation implementation produce an inhomogeneous, multi-phase interstellar medium, allowing us to accurately trace star formation and accretion onto the SMBHs. We study the impact of gas content, morphology, and mass ratio, allowing us to study AGN activity and dynamics across a wide range of relevant conditions. We test when the two AGN are simultaneously detectable, for how long and at which separations. We find that strong dual AGN activity occurs during the late phases of the mergers, at small separations (<1-10 kpc) below the resolution limit of most surveys. Much of the SMBH accretion is not simultaneous, limiting the dual AGN fraction detectable through imaging and spectroscopy to a few percent, in agreement with observational samples.

We study the evolution of the orientation of the orbital plane of massive black hole binaries (BHBs) in rotating stellar systems in which the total angular momentum of the stellar cusp is misaligned with respect to that of the binary. We compare results from direct summation N-body simulations with predictions from a simple theoretical model. We find that the same encounters between cusp stars and the BHB that are responsible for the hardening and eccentricity evolution of the binary, lead to a reorientation of the binary orbital plane. In particular, binaries whose angular momentum is initially misaligned with respect to that of the stellar cusp tend to realign their orbital planes with the angular momentum of the cusp on a timescale of a few hardening times. This is due to angular momentum exchange between stars and the BHB during close encounters, and may have important implications for the relative orientation of host galaxies and radio jets.

The study of the dynamical evolution of massive black hole pairs in mergers is crucial in the context of a hierarchical galaxy formation scenario. The timescales for the formation and the coalescence of black hole binaries are still poorly constrained, resulting in large uncertainties in the expected rate of massive black hole binaries detectable in the electromagnetic and gravitational wave spectra. Here we review the current theoretical understanding of the black hole pairing in galaxy mergers, with a particular attention to recent developments and open issues. We conclude with a review of the expected observational signatures of massive binaries, and of the candidates discussed in literature to date.

We demonstrate that very massive (>10^8\msun), cosmologically nearby (z<1) black hole binaries (MBHBs), which are primary targets for ongoing and upcoming pulsar timing arrays (PTAs), are particularly appealing multimessenger carriers. According to current models for massive black hole formation and evolution, the planned Square Kilometer Array (SKA) will collect gravitational wave signals from thousands of such massive systems, being able to individually resolve and locate in the sky several of them (maybe up to a hundred). By employing a standard model for the evolution of MBHBs in circumbinary discs, with the aid of dedicated numerical simulations, we characterize the gas-binary interplay, identifying possible electromagnetic signatures of the PTA sources. We concentrate our investigation on two particularly promising scenarios in the high energy domain, namely, the detection of X-ray periodic variability and of double broad K\alpha iron lines. Up to several hundreds of periodic X-ray sources with a flux >10^-13 erg s^-1 cm^-2 will be in the reach of upcoming X-ray observatories. Double relativistic K\alpha lines may be observable in a handful of low redshift (z<0.3) sources by proposed deep X-ray probes, such as Athena. (Abridged)

Narrow Line Seyfert 1 galaxies (NLS1) are generally believed to host "under-massive" black holes with respect to the predictions from the host galaxy -- black hole mass scale relations. Black hole masses in NLS1 are estimated from the continuum luminosity and the width of broad emission lines. Here we show that the "mass deficit" can be canceled out if we assume that the broad line region (BLR) in type-1 AGN has a flat geometry, which is seen face-on in NLS1. The detection of relativistic jets aligned along the line of sight in a number of NLS1 supports this picture. Moreover, a flat geometry of the BLR is also suggested by a general trend of the mass deficit as a function of the line width observed in other type-1 AGN, from quasars to BL Lac objects, and is consistent with a simple extension of the Unified Model of AGN to the BLR geometry.

We present the results of a systematic search for massive black hole binaries in the Sloan Digital Sky Survey spectroscopic database. We focus on bound binaries, under the assumption that one of the black holes is active. In this framework, the broad lines associated to the accreting black hole are expected to show systematic velocity shifts with respect to the narrow lines, which trace the rest-frame of the galaxy. For a sample of 54586 quasars and 3929 galaxies at redshifts 0.1<z<1.5 we brute-force model each spectrum as a mixture of two quasars at two different redshifts. The spectral model is a data-driven dimensionality reduction of the SDSS quasar spectra based on a matrix factorization. We identified 32 objects with peculiar spectra. Nine of them can be interpreted as black hole binaries. This doubles the number of known black hole binary candidates. We also report on the discovery of a new class of extreme double-peaked emitters with exceptionally broad and faint Balmer lines. For all the interesting sources, we present detailed analysis of the spectra, and discuss possible interpretations.

In this letter we study the eccentricity evolution of a massive black hole (MBH) binary (MBHB) embedded in a rotating stellar cusp. Following the observation that stars on counter-rotating (with respect to the MBHB) orbits extract angular momentum from the binary more efficiently then their co-rotating counterparts, the eccentricity evolution of the MBHB must depend on the degree of co-rotation (counter-rotation) of the surrounding stellar distribution. Using an hybrid scheme that couples numerical three-body scatterings to an analytical formalism for the cusp-binary interaction, we verify this hypothesis by evolving the MBHB in spherically symmetric cusps with different fractions F of co-rotating stars. Consistently with previous works, binaries in isotropic cusps (F=0.5) tend to increase their eccentricity, and when F approaches zero (counter-rotating cusps) the eccentricity rapidly increases to almost unity. Conversely, binaries in cusps with a significant degree of co-rotation (F>0.7) tend to become less and less eccentric, circularising quite quickly for F approaching unity. Direct N-body integrations performed to test the theory, corroborate the results of the hybrid scheme, at least at a qualitative level. We discuss quantitative differences, ascribing their origin to the oversimplified nature of the hybrid approach.

We study the dynamics of supermassive black hole binaries embedded in circumbinary gaseous discs, with the SPH code Gadget-2. The sub-parsec binary (of total mass M and mass ratio q=1/3) has excavated a gap and transfers its angular momentum to the self--gravitating disc (M_disc=0.2 M). We explore the changes of the binary eccentricity e, by simulating a sequence of binary models that differ in the initial eccentricity e_0, only. In initially low-eccentric binaries, the eccentricity increases with time, while in high-eccentric binaries e declines, indicating the existence of a limiting eccentricity e_crit that is found to fall in the interval [0.6,0.8]. We also present an analytical interpretation for this saturation limit. An important consequence of the existence of e_crit is the detectability of a significant residual eccentricity e_LISA} by the proposed gravitational wave detector LISA. It is found that at the moment of entering the LISA frequency domain e_LISA ~ 10^{-3}-10^{-2}; a signature of its earlier coupling with the massive circumbinary disc. We also observe large periodic inflows across the gap, occurring on the binary and disc dynamical time scales rather than on the viscous time. These periodic changes in the accretion rate (with amplitudes up to ~100%, depending on the binary eccentricity) can be considered a fingerprint of eccentric sub-parsec binaries migrating inside a circumbinary disc.

Only a small fraction of local galaxies harbor an accreting black hole, classified as an active galactic nucleus (AGN). However, many stellar systems are plausibly expected to host black holes, from globular clusters to nuclear star clusters, to massive galaxies. The mere presence of stars in the vicinity of a black hole provides a source of fuel via mass loss of evolved stars. In this paper we assess the expected luminosities of black holes embedded in stellar systems of different sizes and properties, spanning a large range of masses. We model the distribution of stars and derive the amount of gas available to a central black hole through a geometrical model. We estimate the luminosity of the black holes under simple, but physically grounded, assumptions on the accretion flow. Finally we discuss the detectability of 'quiescent' black holes in the local Universe.

We study the properties of the broad line region in blazars by comparing the virial estimate of black hole masses with that derived from the mass of the host galaxies. The former is sensitive to the width of broad lines, i.e., to the projection of the velocity of line-emitting clouds along the line of sight; the latter is not. This comparison allows us to constrain the deprojection factor f, thus revealing general properties of the geometry of the broad line region. We show that blazars tend to have 1) higher f values than the quasars of our reference sample: <f_{BLLacs}>=6.9 +/- 2.3, <f_{blazars}>=5.6 +/- 1.3 and <f_{quasars}>=2.0 +/- 0.3; 2) relatively narrow broad emission lines; 3) modest equivalent widths, as expected because of the occurrence of jet emission at very low inclination angles. In a disc-like sketch of the broad line region, these results indicate a pole-on view of a flat geometry in blazars. This consistently extends the orientation--dependent unified model of active nuclei to the geometry of the broad line region.

We investigate on the spectral properties of an active black hole, member of a massive (10^7 - 10^9 Msun) sub-parsec black hole binary. We work under the hypothesis that the binary, surrounded by a circum-binary disc, has cleared a gap, and that accretion occurs onto the secondary black hole fed by material closer to the inner edge of the disc. Broad line emission clouds orbit around the active black hole and suffer erosion due to tidal truncation at the Roche Lobe surface, following gap opening and orbital decay. We consider three of the most prominent broad emission lines observed in the spectra of AGNs, i.e. CIV, MgII and H{\beta}, and compute the flux ratios between the lines of MgII and CIV (FMgII/FCIV) and those of MgII and H{\beta} (FMgII/FH{\beta}). We find that close black hole binaries have FMgII/FCIV up to one order of magnitude smaller than single black holes. By contrast FMgII/FH{\beta} may be significantly reduced only at the shortest separations. Peculiarly low values of line flux ratios together with large velocity offsets between the broad and narrow emission lines and/or periodic variability in the continuum (on timescales >= years) would identify genuine sub-pc binary candidates.

In this proceeding we explore a pathway to radio-loudness under the hypothesis that retrograde accretion onto giant spinning black holes leads to the launch of powerful jets, as seen in radio loud QSOs and recently in LAT/Fermi and BAT/Swift Blazars. Counter-rotation of the accretion disc relative to the BH spin is here associated to gas-poor galaxy mergers progenitors of giant (missing-light) ellipticals. The occurrence of retrograde accretion enters as unifying element that may account for the radio-loudness/galaxy morphology dichotomy observed in AGN.

We report the discovery of peculiar features in the optical spectrum of 4C+22.25, a flat spectrum radio quasar at z=0.4183 observed in the SDSS and in a dedicated spectroscopic follow-up from the Nordic Optical Telescope. The Hbeta and Halpha lines show broad profiles (FWHM~12,000 km/s), faint fluxes and extreme offsets (Delta v=8,700+/-1,300 km/s) with respect to the narrow emission lines. These features show no significant variation in a time lag of ~3.1 yr (rest frame). We rule out possible interpretations based on the superposition of two sources or on recoiling black holes, and we discuss the virtues and limitations of a massive black hole binary scenario.

Nuclear Stellar Disks (NSDs), of a few tens to hundreds of parsec across, are a common and yet poorly studied feature of early-type galaxies. Still, such small disks represent a powerful tool to constrain the assembling history of galaxies, since they can be used to trace to the epoch when galaxies experienced their last major merger event. By studying the fraction and stellar age of NSDs it is thus possible to test the predictions for the assembly history of early-type galaxies according the current hierarchical paradigm for galaxy formation. In this paper we have produced the most comprehensive census of NSDs in nearby early-type galaxies by searching for such disks in objects within 100 Mpc and by using archival images from the Hubble Space Telescope. We found that NSDs are present in approximately 20% of early-type galaxies, and that the fraction of galaxies with NSDs does not depend on their Hubble type nor on their galactic environment, whereas the incidence of NSDs appears to decline in the most massive systems. Furthermore, we have separated the light contribution of twelve such disks from that of their surrounding stellar bulge in order to extract their physical properties. This doubles the number of decomposed NSDs and although the derived values for their central surface brightness and scale-length are consistent with previous studies they also give a hint of possible different characteristics due to different formation scenario between nuclear disks and other kinds of large galactic disks.

We assess the influence of massive black hole (MBH) ejections from galaxy centres, due to the gravitational radiation recoil, along the cosmic merger history of the MBH population. We discuss the 'danger' of the recoil for MBHs as a function of different MBH spin/orbit configurations and of the host halo cosmic bias, and on how that reflects on the 'occupation fraction' of MBHs. We assess ejection probabilities for mergers occurring in a gas-poor environment, where the MBH binary coalescence is driven by stellar dynamical processes, and the spin/orbit configuration is expected to be isotropically distributed. We contrast this case with the 'aligned' case. The latter is the most realistic situation for 'wet', gas-rich mergers, which are the expectation for high-redshift galaxies. We find that if all halos at z>5-7 host a MBH, the probability of the Milky Way (or similar size galaxy) to host a MBH today is less than 50%, unless MBHs form continuously in galaxies. The 'occupation fraction' of MBHs, intimately related to halo bias and MBH formation efficiency, plays a crucial role in increasing the retention fraction. Small halos, with shallow potential wells and low escape velocities, have a high ejection probability, but the MBH merger rate is very low along their galaxy formation merger hierarchy: MBH formation processes are likely inefficient in such shallow potential wells. Recoils can decrease the overall frequency of MBHs in small galaxies to ~60%, while they have little effect on the frequency of MBHs in large galaxies (at most a 20% effect).

Several active galactic nuclei (AGN) with multiple sets of emission lines separated by over 2000 km/s have been observed recently. These have been interpreted as being due to massive black hole (MBH) recoil following a black hole merger, MBH binaries, or chance superpositions of AGN in galaxy clusters. Moreover, a number of double-peaked AGN with velocity offsets of ~ a few 100 km/s have also been detected and interpreted as being due to the internal kinematics of the narrow line regions or MBH binary systems. Here we reexamine the superposition model. Using the Millennium Run we estimate the total number of detectable AGN pairs as a function of the emission line offset. We show that AGN pairs with high velocity line separations up to ~2000 km/s are very likely to be chance superpositions of two AGN in clusters of galaxies for reasonable assumptions about the relative fraction of AGN. No superimposed AGN pairs are predicted for velocity offsets in excess of ~3000 km/s as the required AGN fractions would violate observational constraints. The high velocity AGN pair numbers predicted here are competitive with those predicted from the models relying on MBH recoil or MBH binaries. However, the model fails to account for the largest emission line velocity offsets that require the presence of MBH binaries.

Using high resolution hydrodynamical simulations, we explore the spin evolution of massive dual black holes orbiting inside a circumnuclear disc, relic of a gas-rich galaxy merger. The black holes spiral inwards from initially eccentric co or counter-rotating coplanar orbits relative to the disc's rotation, and accrete gas that is carrying a net angular momentum. As the black hole mass grows, its spin changes in strength and direction due to its gravito-magnetic coupling with the small-scale accretion disc. We find that the black hole spins loose memory of their initial orientation, as accretion torques suffice to align the spins with the angular momentum of their orbit on a short timescale (<1-2 Myr). A residual off-set in the spin direction relative to the orbital angular momentum remains, at the level of <10 degrees for the case of a cold disc, and <30 degrees for a warmer disc. Alignment in a cooler disc is more effective due to the higher coherence of the accretion flow near each black hole that reflects the large-scale coherence of the disc's rotation. If the massive black holes coalesce preserving the spin directions set after formation of a Keplerian binary, the relic black hole resulting from their coalescence receives a relatively small gravitational recoil. The distribution of recoil velocities inferred from a simulated sample of massive black hole binaries has median <70 km/s much smaller than the median resulting from an isotropic distribution of spins.

We investigate the theoretical expectations for detections of supermassive binary black holes that can be identified as sub-parsec luminous quasars. To-date, only two candidates have been selected in a sample comprising 17,500 sources selected from the Sloan Digital Sky Survey (SDSS) Quasar Catalog at z<0.70 (Boroson & Lauer 2009) In this Letter, we use models of assembly and growth of supermassive black holes (SMBHs) in hierarchical cosmologies to study the statistics and observability of binary quasars at sub-parsec separations. Our goal is twofold: (1) test if such a scarce number of binaries is consistent with theoretical prediction of SMBH merger rates, and (2) provide additional predictions at higher redshifts, and at lower flux levels. We determine the cumulative number of expected binaries in a complete, volume limited sample. Motivated by Boroson & Lauer (2009), we apply the SDSS Quasar luminosity cut (M_i<-22) to our theoretical sample, deriving an upper limit to the observable binary fraction. We find that sub-parsec quasar binaries are intrinsically rare. Our best models predict ~0.01 deg^-2 sub-parsec binary quasars with separations below ~10^4 Schwarzschild radii (v_orb>2000 km/s) at z<0.7, which represent a fraction ~6x10^-4 of unabsorbed quasars in our theoretical sample. In a complete sample of ~10,000 sources, we therefore predict an upper limit of ~10 sub-parsec binary quasars. The number of binaries increases rapidly with increasing redshift. The decreasing lifetime with SMBH binary mass suggests that lowering the luminosity threshold does not lead to a significant increase in the number of detectable sub-parsec binary quasars.

We present an imaging study of the black hole binary candidate SDSS J1536+0441 (z=0.3893), based on deep, high resolution VzK images collected at the ESO/VLT. The images clearly show an asymmetric elongation, indicating the presence of a companion source at ~1" (~5 kpc projected distance) East from the quasar. The host galaxy of the quasar is marginally resolved. We find that the companion source is a luminous galaxy, the light profile of which suggests the presence of an unresolved, faint nucleus (either an obscured AGN or a compact stellar bulge). The study of the environment around the quasar indicates the occurrence of a significant over-density of galaxies with a redshift compatible with z~0.4. This suggests that it resides in a moderately rich cluster of galaxies.

In this paper, we explore the gravitomagnetic interaction of a black hole (BH) with a misaligned accretion disc to study BH spin precession and alignment jointly with BH mass and spin parameter evolution, under the assumption that the disc is continually fed, in its outer region, by matter with angular momentum fixed on a given direction. We develop an iterative scheme based on the adiabatic approximation to study the BH-disc coevolution: in this approach, the accretion disc transits through a sequence of quasi-steady warped states (Bardeen-Petterson effect) and interacts with the BH until the BH spin aligns with the outer angular momentum direction. For a BH aligning with a co-rotating disc, the fractional increase in mass is typically less than a few percent, while the spin modulus can increase up to a few tens of percent. The alignment timescale is between ~ 100 thousands and ~ 1 millions years for a maximally rotating BH accreting at the Eddington rate. BH-disc alignment from an initially counter-rotating disc tends to be more efficient compared to the specular co-rotating case due to the asymmetry seeded in the Kerr metric: counter-rotating matter carries a larger and opposite angular momentum when crossing the innermost stable orbit, so that the spin modulus decreases faster and so the relative inclination angle.

In this Letter we explore the hypothesis that the quasar SDSSJ092712.65+294344.0 is hosting a massive black hole binary embedded in a circumbinary disc. The lightest, secondary black hole is active, and gas orbiting around it is responsible for the blue-shifted broad emission lines with velocity off-set of 2650 km/s, relative to the galaxy rest frame. As the tidal interaction of the binary with the outer disc is expected to excavate a gap, the blue-shifted narrow emission lines are consistent with being emitted from the low-density inhomogeneous gas of the hollow region. From the observations we infer a binary mass ratio q ~ 0.3, a mass for the primary of M1 ~ 2 billion Msun and a semi-major axis of 0.34 pc, corresponding to an orbital period of 370 years. We use the results of cosmological merger trees to estimate the likely-hood of observing SDSSJ092712.65+294344.0 as recoiling black hole or as a binary. We find that the binary hypothesis is preferred being one hundred times more probable than the ejection hypothesis. If SDSSJ092712.65+294344.0 hosts a binary, it would be the one closest massive black hole binary system ever discovered.

Binary black holes occupy a special place in our quest for understanding the evolution of galaxies along cosmic history. If massive black holes grow at the center of (pre-)galactic structures that experience a sequence of merger episodes, then dual black holes form as inescapable outcome of galaxy assembly. But, if the black holes reach coalescence, then they become the loudest sources of gravitational waves ever in the universe. Nature seems to provide a pathway for the formation of these exotic binaries, and a number of key questions need to be addressed: How do massive black holes pair in a merger? Depending on the properties of the underlying galaxies, do black holes always form a close Keplerian binary? If a binary forms, does hardening proceed down to the domain controlled by gravitational wave back reaction? What is the role played by gas and/or stars in braking the black holes, and on which timescale does coalescence occur? Can the black holes accrete on flight and shine during their pathway to coalescence? N-Body/hydrodynamical codes have proven to be vital tools for studying their evolution, and progress in this field is expected to grow rapidly in the effort to describe, in full realism, the physics of stars and gas around the black holes, starting from the cosmological large scale of a merger. If detected in the new window provided by the upcoming gravitational wave experiments, binary black holes will provide a deep view into the process of hierarchical clustering which is at the heart of the current paradigm of galaxy formation. They will also be exquisite probes for testing General Relativity, as the theory of gravity. The waveforms emitted during the inspiral, coalescence and ring-down phase carry in their shape the sign of a dynamically evolving space-time and the proof of the existence of an horizon.

We present a photometric FUV to Ks-band study of the field around quasar SDSS J092712.65+294344.0. The SDSS spectrum of this object shows various emission lines with two distinct redshifts, at z=0.699 and z=0.712. Because of this peculiar spectroscopic feature this source has been proposed as a candidate recoiling or binary black hole. A third alternative model involves two galaxies moving in the centre of a rich galaxy cluster. Here we present a study addressing the possible presence of such a rich cluster of galaxies in the SDSS J092712.65+294344.0 field. We observed the 3.6x2.6 square arcmin field in the Ks-band and matched the NIR data with the FUV and NUV images in the GALEX archive and the ugriz observations in the SDSS. From various colour-colour diagrams we were able to classify the nature of 32 sources, only 6-11 of which have colours consistent with galaxies at z~0.7. We compare these numbers with the surface density of galaxies, stars & quasars, and the expectations for typical galaxy clusters both at low and high redshift. Our study shows that the galaxy cluster scenario is in clear disagreement with the new observations.

We study the orbital evolution and accretion history of massive black hole (MBH) pairs in rotationally supported circumnuclear discs up to the point where MBHs form binary systems. Our simulations have high resolution in mass and space which, for the first time, makes it feasible to follow the orbital decay of a MBH either counter- or co-rotating with respect to the circumnuclear disc. We show that a moving MBH on an initially counter-rotating orbit experiences an "orbital angular momentum flip" due to the gas-dynamical friction, i.e., it starts to corotate with the disc before a MBH binary forms. We stress that this effect can only be captured in very high resolution simulations. Given the extremely large number of gas particles used, the dynamical range is sufficiently large to resolve the Bondi-Hoyle-Lyttleton radii of individual MBHs. As a consequence, we are able to link the accretion processes to the orbital evolution of the MBH pairs. We predict that the accretion rate is significantly suppressed and extremely variable when the MBH is moving on a retrograde orbit. It is only after the orbital angular momentum flip has taken place that the secondary rapidly "lights up" at which point both MBHs can accrete near the Eddington rate for a few Myr. The separation of the double nucleus is expected to be around ~10 pc at this stage. We show that the accretion rate can be highly variable also when the MBH is co-rotating with the disc (albeit to a lesser extent) provided that its orbit is eccentric. Our results have significant consequences for the expected number of observable double AGNs at separations of <100 pc.

We report on key studies on the dynamics of black holes (BHs) in gas-rich galaxy mergers that underscore the vital role played by gas dissipation in promoting BH inspiral down to the smallest scales ever probed with use of high-resolution numerical simulations. In major mergers, the BHs sink rapidly under the action of gas-dynamical friction while orbiting inside the massive nuclear disc resulting from the merger. The BHs then bind and form a Keplerian binary on a scale of 5 pc. In minor mergers, BH pairing proceeds down to the minimum scale explored of 10-100 pc only when the gas fraction in the less massive galaxy is comparatively large to avoid its tidal and/or ram pressure disruption and the wandering of the light BH in the periphery of the main halo. Binary BHs enter the gravitational wave dominated inspiral only when their relative distance is typically of 0.001 pc. If the gas preserves the degree of dissipation expected in a star-burst environment, binary decay continues down to 0.1 pc, the smallest length-scale ever attained. Stalling versus hardening below 0.1 pc is still matter of deep investigations.

Anisotropic gravitational radiation from a coalescing black hole binary is known to impart recoil velocities of up to ~1000 km/s to the remnant black hole. In this context, we study the motion of a recoiling black hole inside a galaxy modelled as an Hernquist sphere, and the signature that the hole imprints on the hot gas, using N-body/SPH simulations. Ejection of the black hole results in a sudden expansion of the gas ending with the formation of a gaseous core, similarly to what is seen for the stars. A cometary tail of particles bound to the black hole is initially released along its trail. As the black hole moves on a return orbit, a nearly spherical swarm of hot gaseous particles forms at every apocentre: this feature can live up to ~ 100 Myr. If the recoil velocity exceeds the sound speed initially, the black hole shocks the gas in the form of a Mach cone in density near each super-sonic pericentric passage. We find that the X-ray fingerprint of a recoiling black hole can be detected in Chandra X-ray maps out to a distance of Virgo. For exceptionally massive black holes the Mach cone and the wakes could be observed out to a few hundred of Mpc. Detection of the Mach cone is found to become of twofold importance: i) as a probe of high-velocity recoils and ii) as an assessment of the scatter of the mass-sigma relation at large black hole masses.

We analyze the available sample of double quasars, and investigate their physical properties. Our sample comprises 85 pairs, selected from the Sloan Digital Sky Survey (SDSS). We derive physical parameters for the engine and the host, and model the dynamical evolution of the pair. First, we compare different scaling relationships between massive black holes and their hosts (bulge mass, velocity dispersion, and their possible redshift dependences), and discuss their consistency. We then compute dynamical friction timescales for the double quasar systems to investigate their frequency and their agreement with scenarios for quasar triggering. Comparing typical merging timescales to expected quasar lifetimes, the fraction of double quasars should be roughly a factor of 10 larger than observed. Additionally, we find that, depending on the correlations between black holes and their hosts, the occurrence of double quasars could be redshift-dependent. Comparison of our models to the SDSS quasar catalog suggests that double quasars should be more common at high redshift. We compare the typical separations at which double quasars are observed to the predictions of merger simulations. We find that the distribution of physical separations peaks at ~30 kpc, with a tail at larger separations (~100-200 kpc). These large separation pairs are instead consistent with unequal mass mergers where gas is dynamically perturbed during the first pericentric passage, but the gas reaches the black hole only at the next apocenter, where the pair is observed.

We study the evolution of a massive black hole pair in a rotationally supported nuclear disc. The distributions of stars and gas mimic the nuclear region of a gas-rich galaxy merger remnant. Using high-resolution SPH simulations, we follow the black hole dynamics and trace the evolution of the underlying background, until the black holes form a binary. We find that the gravitational perturbation of the pair creates a core in the disc density profile, hence decreasing the gas-dynamical drag. This leads the newly formed binary to stall at a separation of ~5 pc. In the early phases of the sinking, black holes lose memory of their initial orbital eccentricity if they co-rotate with the disc, as rotation of the gaseous background promotes circularization of the black hole orbits. Circularization is efficient until the black holes bind in a binary, though in the latest stages of the simulations a residual eccentricity > 0.1 is still present. Black holes are treated as sink particles, allowing for gas accretion. We find that accretion strongly depends on the dynamical properties of the black holes, and occurs preferentially after circularization.

Narrow Line Seyfert 1 galaxies (NLS1s) are generally considered peculiar objects among the broad class of Type 1 active galactic nuclei, due to the relatively small width of the broad lines, strong X-ray variability, soft X-ray continua, weak [OIII], and strong FeII line intensities. The mass M_BH of the central massive black hole (MBH) is claimed to be lighter than expected from known MBH-host galaxy scaling relations, while the accretion rate onto the MBH larger than the average value appropriate to Seyfert 1 galaxies. In this Letter, we show that NLS1 peculiar M_BH and L/L_Edd turn out to be fairly standard, provided that the broad line region is allowed to have a disc-like, rather than isotropic, geometry. Assuming that NLS1s are rather ``normal'' Seyfert 1 objects seen along the disc axis, we could estimate the typical inclination angles from the fraction of Seyfert 1 classified as NLS1s, and compute the geometrical factor relating the observed FWHM of broad lines to the virial mass of the MBH. We show that the geometrical factor can fully account for the "black hole mass deficit" observed in NLS1s, and that L/L_Edd is (on average) comparable to the value of the more common broad line Seyfert 1 galaxies.

If massive black holes (BHs) are ubiquitous in galaxies and galaxies experience multiple mergers during their cosmic assembly, then BH binaries should be common albeit temporary features of most galactic bulges. Observationally, the paucity of active BH pairs points toward binary lifetimes far shorter than the Hubble time, indicating rapid inspiral of the BHs down to the domain where gravitational waves lead to their coalescence. Here, we review a series of studies on the dynamics of massive BHs in gas-rich galaxy mergers that underscore the vital role played by a cool, gaseous component in promoting the rapid formation of the BH binary. The BH binary is found to reside at the center of a massive self-gravitating nuclear disc resulting from the collision of the two gaseous discs present in the mother galaxies. Hardening by gravitational torques against gas in this grand disc is found to continue down to sub-parsec scales. The eccentricity decreases with time to zero and when the binary is circular, accretion sets in around the two BHs. When this occurs, each BH is endowed with it own small-size (< 0.01 pc) accretion disc comprising a few percent of the BH mass. Double AGN activity is expected to occur on an estimated timescale of < 1 Myr. The double nuclear point-like sources that may appear have typical separation of < 10 pc, and are likely to be embedded in the still ongoing starburst. We note that a potential threat of binary stalling, in a gaseous environment, may come from radiation and/or mechanical energy injections by the BHs. Only short-lived or sub-Eddington accretion episodes can guarantee the persistence of a dense cool gas structure around the binary necessary for continuing BH inspiral.

We present a study on the dynamics of massive BHs in galaxy mergers, obtained from a series of high-resolution N-Body/SPH simulations. The presence of a gaseous component is essential for the rapid formation of an eccentric (Keplerian) BH binary, that resides at the center of a massive (~10^9 Msun) turbulent nuclear disc. Using physically and/or numerically motivated recipes, we follow the accretion history of the BHs during the merger. The mass of the BHs increases as large central inflows of gas occur inside each galaxy, and their mass ratio varies with time. Given the encountered strong degeneracy between numerical resolution and physical assumptions, we suggest here three possible paths followed by the galaxies and the BHs during a merger in order to fulfill the M-sigma relation : Adjustment, Symbiosis, and BH Dominance. In an extremely high resolution run, we resolved the turbulent gas pattern down to parsec scales, and found that BH feedback is expected to be effective near the end of the merger. We then trace the BH binary orbit down to a scale of 0.1 pc modeling the nuclear disc as an equilibrium Mestel disc composed either of gas, gas and stars, or just stars. Under the action of dynamical friction against the rotating gaseous and/or stellar background the orbit circularizes. When this occurs, each BH is endowed with its own small-size (~0.01 pc) accretion disc comprising a few percent of the BH mass. Double AGN activity is expected to occur on an estimated timescale of ~10 Myrs, comparable to the inspiral time. The double nuclear point--like sources that may appear have typical separations of ~10 pc, and are likely to be embedded in the still ongoing starburst.

The dynamics of two massive black holes in a rotationally supported nuclear disc of 10^8 solar masses is explored using N-Body/SPH simulations. Gas and star particles are co-present in the disc. Described by a Mestel profile, the disc has a vertical support provided by turbulence of the gas, and by stellar velocity dispersion. A primary black hole of 4 million solar masses is placed at the centre of the disc, while a secondary black hole is set initially on an eccentric co-rotating orbit in the disc plane. Its mass is in a 1 to 1, 1 to 4, and 1 to 10 ratio, relative to the primary. With this choice, we mimic the dynamics of black hole pairs released in the nuclear region at the end of a gas-rich galaxy merger. It is found that, under the action of dynamical friction, the two black holes form a close binary in ~10 Myrs. The inspiral process is insensitive to the mass fraction in stars and gas present in the disc and is accompanied by the circularization of the orbit. We detail the gaseous mass profile bound to each black hole that can lead to the formation of two small Keplerian discs, weighing ~2 % of the black hole mass, and of size ~0.01 pc. The mass of the tightly (loosely) bound particles increases (decreases) with time as the black holes spiral into closer and closer orbits. Double AGN activity is expected to occur on an estimated timescale of ~10 Myrs, comparable to the inspiral timescale. The double nuclear point-like sources that may appear during dynamical evolution have typical separations of ~10 pc.

We explore the nature of possible electromagnetic counterparts of coalescences of massive black hole binaries at cosmological distances detectable by the Laser Interferometer Space Antenna (LISA). An electromagnetic precursor, during the last year of gravitational wave (GW)-driven inspiral, or an afterglow within few years after coalescence, may highlight the position in the sky of galaxies hosting LISA sources. We show that observations of precursors and afterglows are mutually exclusive, depending on the mass of the primary black hole. Precursors are expected to occur in binaries where the primary (more massive) black hole is heavier than ~10^7 \Msun. They may correspond to on-off states of accretion, i.e., to a bright X-ray source decaying into quiescence before black hole coalescence, and are likely associated to disturbed galaxies showing signs of ongoing starbursts. Coalescences of lighter binaries, with masses <5x10^6 \Msun, lack of any precursor, as gas is expected to be consumed long before the GW-driven orbital decay. Such events would not be hosted by (massive) galaxies with an associated starburst, given the slow binary inspiral time compared to the typical time scale of starbursts. By contrast, coalescence, for such light binaries, is followed by an electromagnetic afterglow, i.e., an off-on accretion state rising in <20 yrs. Using a cosmological merger tree algorithm, we show that future X-ray missions such as XEUS will be able to identify, in 20 yrs operation, almost all the massive BH binary detectable by LISA, and, in only 5 yrs, all the LISA sources at z>6.

We study the inspiral of double black holes, with masses in the LISA window of detectability, orbiting inside a massive circum-nuclear disc. Using high-resolution SPH simulations, we follow the black hole dynamics in the early phase when gas-dynamical friction acts on the black holes individually, and continue our simulation until they form a close binary. We find that in the early sinking the black holes lose memory of their initial orbital eccentricity if they co-rotate with the gaseous disc, forming a binary with a low eccentricity, consistent with zero within our numerical resolution limit. The cause of circularization resides in the rotation present in the gaseous background where dynamical friction operates. Circularization may hinder gravitational waves from taking over and leading the binary to coalescence. In the case of counter-rotating orbits the initial eccentricity does not decrease, and the black holes may bind forming an eccentric binary. When dynamical friction has subsided, for equal mass black holes and regardless their initial eccentricity, angular momentum loss, driven by the gravitational torque exerted on the binary by surrounding gas, is nevertheless observable down to the smallest scale probed. In the case of unequal masses, dynamical friction remains efficient down to our resolution limit, and there is no sign of formation of any ellipsoidal gas distribution that may further harden the binary. During inspiral, gravitational capture of gas by the black holes occurs mainly along circular orbits: eccentric orbits imply high relative velocities and weak gravitational focusing. Thus, AGN activity may be excited during the black hole pairing process and double active nuclei may form when circularization is completed, on distance-scales of tens of pcs.