Search

Quasi-periodic eruptions (QPEs) are repeated X-ray flares from galactic nuclei. Despite some diversity in the recurrence and amplitude of eruptions, their striking regularity has motivated theorists to associate QPEs with orbital systems. Among the known QPE sources, eRO-QPE2 has shown the most regular flare timing and luminosity since its discovery. We report here on its long-term evolution over $\sim3.3\,$yr from discovery and find that: i) the average QPE recurrence time per epoch has decreased over time, albeit not at a uniform rate; ii) the distinct alternation between consecutive long and short recurrence times found at discovery has not been significant since; iii) the spectral properties, namely flux and temperature of both eruptions and quiescence components, have remained remarkably consistent within uncertainties. We attempted to interpret these results as orbital period and eccentricity decay coupled with orbital and disk precession. However, since gaps between observations are too long, we are not able to distinguish between an evolution dominated by just a decreasing trend, or by large modulations (e.g. due to the precession frequencies at play). In the former case, the observed period decrease is roughly consistent with that of a star losing orbital energy due to hydrodynamic gas drag from disk collisions, although the related eccentricity decay is too fast and additional modulations have to contribute too. In the latter case, no conclusive remarks are possible on the orbital evolution and the nature of the orbiter due to the many effects at play. However, these two cases come with distinctive predictions for future X-ray data: in the former, we expect all future observations to show a shorter recurrence time than the latest epoch, while in the latter we expect some future observations to be found with a larger recurrence, hence an apparent temporary period increase., Comment: Version updated after re-submission to A&A with minor revisions. Appendix B in v1 was moved to the main text

We present ALMA CO(1-0) observations of the nearby LIRG galaxy pair IRAS05054+1718 with a new analysis of X-ray data collected between 2012 and 2021 using NuSTAR, Swift, and XMM-Newton. The western component of the pair, NED01, hosts a Seyfert 1.9 nucleus launching a powerful X-ray UFO. Our X-ray spectral analysis suggests the UFO could be variable or multi-component in velocity and constrains its momentum flux to $\dot p^{X-ray}_{out} \sim (4\pm2)\times 10^{34}$ gcms$^{-2}$. ALMA CO(1-0) observations include also the eastern component of the pair, a LIRG with no clear evidence for an AGN. We study the CO(1-0) kinematics in the two galaxies using the 3D-BAROLO code. In both sources, we can model the bulk of the CO(1-0) emission with rotating disks and, after subtracting the best-fit models, we detect compact residual emission at S/N=15 within $\sim3$kpc from the centre. A molecular outflow in NED01, if present, cannot be brighter than such residuals, implying an upper limit on its outflow rate of $\dot{M}^{mol}_{out} \lesssim 19\pm14~M_{\odot}~yr^{-1}$ and on its momentum rate of $\dot p^{mol}_{out} \lesssim (2.7\pm2.4) \times 10^{34}$gcms$^{-1}$. Combined with the revised energetics of the X-ray wind, we derive an upper limit on the momentum rate ratio of $\dot{p}^{mol}_{out}/\dot{p}^{X-ray}_{out}<0.67$. We discuss these results in the context of the expectations of AGN feedback models, and we propose the X-ray disk wind in NED01 has not significantly impacted the molecular gas reservoir (yet), and we can constrain its effect to be much smaller than expectations of AGN ''energy-driven'' feedback models. We also consider and discuss the hypothesis of asymmetries of the molecular disk not properly captured by the 3D-BAROLO code. Our results highlight the challenges in testing the predictions of popular AGN disk-wind feedback theories, even with good quality multi-wavelength observations., Comment: Accepted by A&A

3C186 is a powerful radio loud quasar (QSO) at the center of a cool-core cluster at z=1.06. Previous studies reported evidence for a projected spatial offset of ~1'' between the isophotal center of the galaxy and the point-source QSO as well as a spectral shift of ~2000 km/s between the narrow and broad line region of the system. In this work we report high-resolution molecular gas CO(4-3) observations of the system taken with NOEMA interferometer. We clearly detect a large reservoir of molecular gas, $M_{H_2}\sim8\times10^{10}~M_\odot$, that is cospatial with the host galaxy and likely associated with a rotating disk-like structure. We firmly confirm both the spatial offset of the galaxy's gas reservoir with respect to the continuum emission of the QSO and the spectral offset with respect to the redshift of the broad line region. Our morphological and kinematical analysis confirms that the most likely scenario to explain the 3C186 system is that the QSO is a kicked super-massive black hole (SMBH), which we believe may have resulted from a strong gravitational wave recoil as two SMBHs coalesced after the merger of their host galaxies., Comment: 9 pages, Astronomy & Astrophysics Letters in press

Although thousands of galaxy mergers are known, only a handful of sub-kiloparsec-scale supermassive black hole (SMBH) pairs have been confirmed so far, leaving a huge gap between the observed and predicted numbers of such objects. In this work, we present a detailed analysis of the Sloan Digital Sky Survey optical spectrum and of near-infrared (NIR) diffraction limited imaging of SDSS~J1431+4358. This object is a local radio-quiet type 2 active galactic nucleus (AGN) previously selected as a double AGN candidate on the basis of the double-peaked [OIII] emission line. The NIR adaptive optics-assisted observations were obtained at the Large Binocular Telescope with the LUCI+FLAO camera. We found that most of the prominent optical emission lines are characterized by a double-peaked profile, mainly produced by AGN photoionization. Our spectroscopical analysis disfavors the hypothesis that the double-peaked emission lines in the source are the signatures of outflow kinematics, leaving open the possibility that we are detecting either the rotation of a single narrow-line region or the presence of two SMBHs orbiting around a common central potential. The latter scenario is further supported by the high-spatial resolution NIR imaging: after subtracting the dominant contribution of the stellar bulge component in the host galaxy, we detect two faint nuclear sources at r<0.5 kpc projected separation. Interestingly, the two sources have a position angle consistent with that defined by the two regions where the [OIII] double peaks most likely originate. Aside from the discovery of a promising sub-kiloparsec scale dual AGN, our analysis shows the importance of an appropriate host galaxy subtraction in order to achieve a reliable estimate of the incidence of dual AGNs at small projected separations., Comment: 11 pages, 8 figures, 1 table, accepted by A&A (revised to match the A&A version)

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., Comment: MNRAS, in press

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., Comment: 11 pages, 8 figures, accepted for publication in MNRAS

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., Comment: 33 pages, 25 figures, 3 appendices. Replaced with accepted version in MNRAS. Reference website is: http://qsfit.inaf.it/

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., Comment: 9 pages, 6 figures, 3 tables in Monthly Notices of the Royal Astronomical Society 2016

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., Comment: 12 pages, 8 figures, 2 tables. Accepted for publication in MNRAS

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., Comment: 13 pages, 11 figures, 1 table

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., Comment: 10 pages, 6 figures, accepted for publication on MNRAS

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., Comment: accepted for publication in Astronomy & Astrophysics

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., Comment: 20 pages; submitted to the European Space Agency on May 24th, 2013 for the L2/L3 selection of ESA's Cosmic Vision program

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., Comment: 2 pages, in Proceedings of the Workshop "X-ray astronomy: towards the next 50 years!", Milano, October 1-5, 2012, Mem. S.A.It in press

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., Comment: 32 pages, 11 figures, MNRAS submitted

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., Comment: 14 pages, 7 figures, accepted for publication in ApJ

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., Comment: 7 pages, 3 figure, accepted for publication in MNRAS

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., Comment: MNRAS accepted. 14 pages, 6 figures

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., Comment: 4 Figures. Accepted for publication in Advances in Astronomy

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), Comment: 19 pages, 11 figures, submitted to MNRAS, minor revision of the reference list

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

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., Comment: Proceeding of the workshop: "Narrow-Line Seyfert 1 Galaxies and their place in the Universe", April 04-06, 2011, Milan

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., Comment: 7 pages, 4 figures. Accepted for publication in MNRAS

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., Comment: Accepted for publication in MNRAS 2010 May 08; 21 pages, 33 figures and 6 tables; High resolution version http://star.herts.ac.uk/~hledo/paper/ledo10.pdf

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., Comment: To appear in the proceedings of the conference "Accretion and Ejection in AGN: A global view, June 22-26 2009 - Como, Italy"

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., Comment: 4 pages, 2 figures, submitted to ApJ Letters

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., Comment: 11 pages, 3 figures. Accepted for publication in MNRAS

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., Comment: 4 figures, 1 table. Accepted for publication in ApJ Letters

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., Comment: 20 pages, 15 figures (jpg or png format), 3 tables; to be published in MNRAS

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., Comment: Invited Review to appear on Advanced Science Letters (ASL), Special Issue on Computational Astrophysics, edited by Lucio Mayer

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., Comment: 10 pages, 2 figures. Accepted for publication in MNRAS

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., Comment: 8 pages, 5 figures, to appear in Classical and Quantum Gravity, Lisa-7 Special Issue

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., Comment: ApJ Letters, in press

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., Comment: 8 pages, 4 figures, revised version accepted to MNRAS

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., Comment: Accepted for publication in MNRAS Letters

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., Comment: 4 pages, 3 figures, to appear in the proceedings of the Conference "The Central Kiloparsec: Active Galactic Nuclei and Their Hosts", Ierapetra 4-6 June 2008, Greece, Memorie della Societa Astronomica Italiana

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., Comment: 9 pages, 8 figures, new simulations added, accepted for publication in MNRAS

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., Comment: 5 pages, 3 figures. Accepted for publication in MNRAS Letters. Wrong version was uploaded! Check for the correct one in the replacement

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., Comment: To appear in "2007 STScI Spring Symposium: Black Holes", eds. M. Livio & A. M. Koekemoer, Cambridge University Press, 25 pages, 12 figures

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., Comment: 10 pages, 5 figures, Proceedings of the Conference "The Multicoloured Landscape of Compact Objects and their Explosive Origins", Cefalu` 2006