BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Massimo Dotti

Massimo Dotti appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-ph.HEgr-qcastro-ph.IMastro-phastro-ph.SRphysics.flu-dyn

Catalog footprint

What is connected

54works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

54 published item(s)

preprint2022arXiv

Black hole spin evolution in warped accretion discs

Massive black holes (BHs) inhabiting galactic nuclei can be described by two parameters only, i.e. mass and spin, that change through cosmic time in response to accretion and merger events. While most numerical simulations accurately track the BH mass, spin evolution is rarely taken into account. In this work, we implement and validate a self-consistent sub-grid model for the evolution of the BH mass and spin via gas accretion in the hydrodynamics code GIZMO. The model assumes that accretion from resolved scales does not occur instantaneously, but is mediated by a sub-grid geometrically thin $α$-disc. After validating our model semi-analytically, we test it in an idealized environment consisting of a circumnuclear disc, where gas accretion onto the accretion disc is consistently determined by GIZMO. In the absence of any accretion-related feedback, the spin evolution closely traces that observed in the semi-analytical models and depends on the free parameters of our implementation, such as the initial BH spin, angular momentum of the accretion disc, and the radius at which the gas inflow circularises. In GIZMO, we also couple our model with the biconical-outflow model presented in a companion paper, wherein the feedback axis is always aligned with the BH spin. In this last case, the evolution of the central BH differs significantly from the previous cases, since the feedback process modifies the gas dynamics and its inflow rates from resolved scales. Such an interaction cannot be modeled by simple semi-analytical models and should be treated using full N-body hydrodynamical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Disc instability and bar formation: view from the IllustrisTNG simulations

We make use of z = 0 samples of strongly barred and unbarred disc galaxies from the TNG100 and TNG50 cosmological hydrodynamical simulations to assess the performance of the simple disc instability criterion proposed by Efstathiou, Lake & Negroponte (1982) (ELN-criterion). We find that strongly barred galaxies generally assemble earlier, are more star-dominated in their central regions, and have more massive and more compact discs than unbarred galaxies. The ELN-criterion successfully identifies ~75% and ~80% of the strongly barred and the unbarred galaxies, respectively. Strongly barred galaxies that the criterion fails to identify tend to have more extended discs, higher spin values and bars that assembled later than is typical for the bulk of the barred population. The bars in many of these cases appear to be produced by an interaction with a close neighbour (i.e. to be externally triggered) rather than to result from secular growth in the disc. On the other hand, we find that unbarred galaxies misclassified as barred by the ELN-criterion typically have stellar discs similar to those of barred galaxies, although more extended in the vertical direction and less star-dominated in their central regions, possibly reflecting later formation times. In addition, the bulge component of these galaxies is significantly more prominent at early times than in the strongly barred sample. Thus, the ELN-criterion robustly identifies secular bar instabilities in most simulated disc galaxies, but additional environmental criteria are needed to account for interaction-induced bar formation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dual AGN in the Horizon-AGN simulation and their link to galaxy and massive black hole mergers, with an excursus on multiple AGN

The occurrence of dual active galactic nuclei (AGN) on scales of a few tens of kpc can be used to study merger-induced accretion on massive black holes (MBHs) and to derive clues on MBH mergers, using dual AGN as a parent population of precursors. We investigate the properties of dual AGN in the cosmological simulation Horizon-AGN. We create catalogs of dual AGN selected with distance and luminosity criteria, plus sub-catalogs where further mass cuts are applied. We divide the sample into dual AGN hosted in different galaxies, on the way to a merger, and into those hosted in one galaxy, after the galaxy merger has happened. We find that the relation between MBH and galaxy mass is similar to that of general AGN population and we compare the properties of dual AGN also with a control sample, discussing differences and similarities in masses and Eddington ratios. The typical mass ratio of galaxy mergers associated to dual AGN is 0.2, with mass loss in the smaller galaxy decreasing the mass ratio as the merger progresses. Between 30 and 80 per cent of dual AGN with separations between 4 and 30~kpc can be matched to an ensuing MBH merger. The dual AGN fraction increases with redshift and with separation threshold, although above 50~kpc the increase of multiple AGN limits that of duals. Multiple AGN are generally associated with massive halos, and mass loss of satellites shapes the galaxy-halo relation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Extreme mass ratio inspirals and tidal disruption events in nuclear clusters. I. Time dependent rates

In this paper we develop a computationally efficient, two-population, time-dependent Fokker-Plank approach in the two dimensions of energy and angular momentum to study the rates of tidal disruption events (TDEs), extreme mass ratio inspirals (EMRIs) and direct plunges occurring around massive black holes (MBHs) in galactic nuclei. We test our code by exploring a wide range of the astrophysically relevant parameter space, including MBH masses, galaxy central densities and inner density slopes. We find that mass segregation and, more in general, the time dependency of the distribution function regulate the event rate: TDEs always decline with time, whereas EMRIs and plunges reach a maximum and undergo a subsequent nearly exponential decay. Once suitably normalized, the rates associated to different choices of MBH mass and galaxy density overlap nearly perfectly. Based on this, we provide a simple scaling that allows to reproduce the time-dependent event rates for any MBH mass and underlying galactic nucleus. Although our peak rates are in general agreement with the literature relying on the steady-state (non-time dependent) assumption, those can be sustained on a timescale that strongly depends on the properties of the system. In particular this can be much shorter than a Gyr for relatively light MBHs residing in dense systems. This warns against using steady state models to compute global TDE, EMRI and plunge rates and calls for a more sophisticated, time dependent treatment of the problem.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Extreme Mass Ratio Inspirals triggered by Massive Black Hole Binaries: from Relativistic Dynamics to Cosmological Rates

Extreme mass ratio inspirals (EMRIs) are compact binary systems characterized by a mass-ratio $q=m/M$ in the range $~10^{-9}-10^{-4}$ and represent primary gravitational wave (GW) sources for the forthcoming Laser Interferometer Space Antenna (LISA). While their standard formation channel involves relaxation processes deflecting compact objects on very low angular momentum orbits around the central massive black hole, a number of alternative formation channels has been proposed, including binary tidal break-up, migration in accretion disks and secular and chaotic dynamics around a massive black hole binary (MBHB). In this work, we take an extensive closer look at this latter scenario, investigating how EMRIs can be triggered by a MBHBs, formed in the aftermath of galaxy mergers. By employing a suite of relativistic three-body simulations, we evaluate the efficiency of EMRI formation for different parameters of the MBHB, assessing the importance of both secular and chaotic dynamics. By modelling the distribution of compact objects in galaxy nuclei, we estimate the resulting EMRI formation rate, finding that EMRI are produced in a sharp burst, with peak rates that are 10-100 times higher than the standard two-body relaxation channel, lasting for 10$^6$--10$^8$ years. By coupling our results with an estimate of the cosmic MBHB merger rate, we finally forecast that LISA could observe ${\cal O}(10)$ EMRIs per year formed by this channel.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Maximally informed Bayesian modelling of disc galaxies

Dissecting the underlying structure of galaxies is of main importance in the framework of galaxy formation and evolution theories. While a classical bulge+disc decomposition of disc galaxies is usually taken as granted, this is only rarely solidly founded upon the full exploitation of the richness of data arising from spectroscopic studies with integral field units. In this work we describe a fully Bayesian estimation method of the global structure of disc galaxies which makes use of the wealth of photometric, kinematic, and mass-to-light ratio data, and that can be seen as a first step towards a machine-learning approach, certainly needed when dealing with larger samples of galaxies. Ours is a novel, hybrid line of action in tackling the problem of galactic parameter estimation, neither purely photometric nor orbit-based. Being rooted on a nested sampler, our code, which is available publicly as an online repository, allows for a statistical assessment of the need for multiple components in the dissecting process. As a first case-study the GPU-optimized code is applied to the S0 galaxy NGC-7683, finding that in this galaxy a pseudo-bulge, possibly the remnant of a bar-like structure, does exist in the center of the system. These results are then tested against the publicly available, orbit-based code DYNAMITE, finding substantial agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Morphological decomposition of TNG50 galaxies: methodology and catalogue

We present MORDOR (MORphological DecOmposeR, a new algorithm for structural decomposition of simulated galaxies based on stellar kinematics. The code measures the properties of up to five structural components (a thin/cold and a thick/warm disc, a classical and a secular bulge, and a spherical stellar halo), and determines the properties of a stellar bar (if present). A comparison with other algorithms presented in the literature yields overall good agreement, with MORDOR displaying a higher flexibility in correctly decomposing systems and identifying bars in crowded environments (e.g. with ongoing fly-bys, often observable in cosmological simulations). We use MORDOR to analyse galaxies in the TNG50 simulation and find the following: ($i$) the thick disc component undergoes the strongest evolution in the binding energy-circularity plane, as expected when disc galaxies decrease their turbulent-rotational support with cosmic time; ($ii$) smaller galaxies (with stellar mass, $10^{9} \lesssim M_{*} / {\rm M_{\odot}} \leq 5 \times 10^{9}$) undergo a major growth in their disc components after $z\sim 1$, whereas ($iii$) the most massive galaxies ($5 \times 10^{10} < M_{*} / {\rm M_{\odot}} \leq 5\times10^{11}$) evolve toward more spheroidal dominated objects down to $z=0$ due to frequent gravitational interactions with satellites; ($iv$) the fraction of barred galaxies grows rapidly at high redshift and stabilizes below $z\sim 2$, except for the most massive galaxies that show a decrease in the bar occupation fraction at low redshift; ($v$) galaxies with $M_{*} \sim 10^{11}~{\rm M_{\odot}}$ exhibit the highest relative occurrence of bars at $z=0$, in agreement with observational studies. We publicly release MORDOR and the morphological catalogue of TNG50 galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Non-isotropic feedback from accreting spinning black holes

Active galactic nuclei (AGN) are massive black holes (BHs) caught in the act of accreting gas at the centre of their host galaxies. Part of the accreting mass is converted to energy and released into the surrounding medium, in a process loosely referred to as AGN feedback. Most numerical simulations include AGN feedback as a sub-grid model, wherein energy or momentum (or both) is coupled to the nearby gas. In this work, we implement a new momentum-driven model in the hydrodynamics code GIZMO, in which accretion from large scales is mediated by a sub-grid accretion disc model, and gas particles are stochastically kicked over a bi-conical region, to mimic observed kinetic winds. The feedback cone's axis can be set either parallel to the angular momentum of the gas surrounding the BH or to the BH spin direction, which is self-consistently evolved within the accretion-disc model. Using a circumnuclear disc (CND) as a test bed, we find that (i) the conical shape of the outflow is always visible and is weakly dependent on the launching orientation and aperture, resulting in comparable mass inflows and outflows; (ii) the cone's orientation is also similar amongst our tests, and it is not always the same as the initial value, due to the interaction with the CND playing a crucial role in shaping the outflow; and (iii) the velocity of the outflow, instead, differs and strongly depends on the interplay with the CND.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the detectability of gravitational waves from primordial black holes orbiting Sgr A*

In this work we characterize the expected gravitational wave signal detectable by the planned space-borne interferometer LISA and the proposed next generation space-borne interferometer $μ$Ares arising from a population of primordial black holes orbiting Sgr A*, the super-massive black hole at the Galactic center. Assuming that such objects indeed form the entire diffuse mass allowed by the observed orbit of S2 in the Galactic center, under the simplified assumption of circular orbits and monochromatic mass function, we assess the expected signal in gravitational waves, either from resolved and non-resolved sources. We estimate a small but non negligible chance of $\simeq$ 10% of detecting one single 1 M$_{\odot}$ primordial black hole with LISA in a 10-year-long data stream, while the background signal due to unresolved sources would essentially elude any reasonable chance of detection. On the contrary, $μ$Ares, with a $\simeq$ 3 orders-of-magnitude better sensitivity at $\simeq$ 10$^{-5}$ Hz, would be able to resolve $\simeq$ 140 solar mass primordial black holes in the same amount of time, while the unresolved background should be observable with an integrated signal-to-noise ratio $\gtrsim$ 100. Allowing the typical PBH mass to be in the range 0.01-10 M$_{\odot}$ would increase LISA chance of detection to $\simeq$ 40% towards the lower limit of the mass spectrum. In the case of $μ$Ares, instead, we find a "sweet spot" just about 1 M$_{\odot}$, a mass for which the number of resolvable events is indeed maximized.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Re-weighting dark matter in disc galaxies: a new general relativistic observational test

A recent analysis of data from the ESA Gaia mission demonstrated that the kinematics of stars in the Milky Way can be modelled without invoking the presence of dark matter whatsoever. Indeed, the higher-than-Keplerian velocities observed in outer stars can be ascribed to the properties of the general relativistic (GR) metric assumed to describe the Galaxy. Here we generalize the concept, and derive the most general exact GR model for a stationary axi-symmetric dynamically cold dust structure. We explicitly show how deviations from the commonly adopted Newtonian dynamics are indeed manifest even at low velocities and low densities. We provide for the first time a detailed description of the frequency shift experienced by photons travelling from any emission site within an external disc galaxy to the observer, relating the outcome of the shift measurement to the gravitational properties of the galaxy. Finally, we devise a novel, groundbreaking observational test potentially able to fully characterize the GR metric under the minimal set of assumptions mentioned above. The proposed experiment exploits the effects non-diagonal GR terms have on the frequency shift of photons, ultimately providing a test to evaluate whether dark matter is actually required by disc galaxy kinematics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The buildup of strongly barred galaxies in the TNG100 simulation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The evolution of the barred galaxy population in the TNG50 simulation

We use the magnetic-hydrodynamical simulation TNG50 to study the evolution of barred massive disc galaxies. Massive spiral galaxies are already present as early as $z=4$, and bar formation takes place already at those early times. The bars grow longer and stronger as the host galaxies evolve, with the bar sizes increasing at a pace similar to that of the disc scale lengths. The bar fraction mildly evolves with redshift for galaxies with $M_{*}\geq10^{10}M\odot$, being greater than $\sim40\%$ at $0.5<z<3$ and $\sim30\%$ at $z=0$. When bars larger than a given physical size ($\geq 2\,\rm kpc$) or the angular resolution limit of twice the I-band angular PSF FWHM of the HST are considered, the bar fraction dramatically decreases with increasing redshift, reconciling the theoretical predictions with observational data. We find that barred galaxies have an older stellar population, lower gas fractions and star formation rates than unbarred galaxies. In most cases, the discs of barred galaxies assembled earlier and faster than the discs of unbarred galaxies. We also find that barred galaxies are typical in haloes with larger concentrations and smaller spin parameters than unbarred galaxies. Furthermore, the inner regions of barred galaxies are more baryon-dominated than those of unbarred galaxies but have comparable global stellar mass fractions. Our findings suggest that the bar population could be used as a potential tracer of the buildup of disc galaxies and their host haloes. With this paper, we release a catalogue of barred galaxies in TNG50 at $6$ redshifts between $z=4$ and $z=0$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The role of bars on the dynamical-friction driven inspiral of massive perturbers

In this paper, we systematically explore the impact of a galactic bar on the inspiral time-scale of a massive object (MO) within a Milky Way-like galaxy. We integrate the orbit of MOs in a multi-component galaxy model via a semi-analytical approach that accounts for dynamical friction generalized to rotationally supported backgrounds. We compare the MO evolution in a galaxy featuring a Milky Way-like rotating bar to the evolution within an analogous axisymmetric galaxy without the bar. In agreement with previous studies, we find that the bar presence may significantly affect the inspiral, sometimes making it shorter by a factor of a few, sometimes hindering it for a Hubble time. The erratic behaviour is mainly impacted by the relative phase at which the MO encounters the stronger bar-induced resonances. In particular, the effect of the bar is more prominent for initially in-plane, prograde MOs, especially those crossing the bar co-rotation radius or outer Lindblad resonance. In the barred galaxy, we find the sinking of the most massive MOs (>~10^7.5 Msun) approaching the galaxy from large separations (>~8 kpc) to be most efficiently hampered. Neglecting the effect of global torques associated with the non-symmetric mass distribution is thus not advisable even within an idealized, smooth galaxy model; we further note that spiral patterns are unlikely to affect the inspiral due to their transient and fluctuating nature. We speculate that the sinking efficiency of massive black holes involved in minor galaxy mergers may be hampered in barred galaxies, making them less likely to host a gravitational wave signal accessible to low-frequency detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Dynamical evolution of massive perturbers in realistic multi-component galaxy models I: implementation and validation

Galaxies are self-gravitating structures composed by several components encompassing spherical, axial and triaxial symmetry. Although real systems feature heterogeneous components whose properties are intimately connected, semi-analytical approaches often exploit the linearity of the Poisson's equation to represent the potential and mass distribution of a multi-component galaxy as the sum of the individual components. In this work, we expand the semi-analytical framework developed in Bonetti et al. (2020) by including both a detailed implementation of the gravitational potential of exponential disc (modelled with a ${\rm sech}^2$ and an exponential vertical profile) and an accurate prescription for the dynamical friction experienced by massive perturbers in composite galaxy models featuring rotating disc structures. Such improvements allow us to evolve arbitrary orbits either within or outside the galactic disc plane. We validate the results obtained by our numerical model against public semi-analytical codes as well as full N-body simulations, finding that our model is in excellent agreement to the codes it is compared with. The ability to reproduce the relevant physical processes responsible for the evolution of massive perturber orbits and its computational efficiency make our framework perfectly suited for large parameter-space exploration studies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Black hole mergers from dwarf to massive galaxies with the NewHorizon and Horizon-AGN simulations

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. We use the same method to analyse the two simulations, post-processing 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dynamical friction-driven orbital circularisation in rotating discs: a semi-analytical description

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

From galactic nuclei to the halo outskirts: tracing supermassive black holes across cosmic history and environments

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Global torques and stochasticity as the drivers of massive black hole pairing in the young Universe

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the eccentricity evolution of massive black hole binaries in stellar backgrounds

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Real galaxy mergers from galaxy pair catalogs

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 $Δz$) gives similar results than the most complex one (based on a neural network analyzing $d$, $Δ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 $Δz_\mathrm{th} \sim 10^{-3}$, in agreement with recent results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Quest for Dual and Binary Supermassive Black Holes: A Multi-Messenger View

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bar-driven evolution and quenching of spiral galaxies in cosmological simulations

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]

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nuclear discs as clocks for the assembly history of early-type galaxies: the case of NGC4458

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Retrograde binaries of massive black holes in circum-binary accretion discs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Robust automatic photometry of local galaxies from SDSS Dissecting the color magnitude relation with color profiles

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Black hole accretion versus star formation rate: theory confronts observations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Growing black holes and galaxies: black hole accretion versus star formation rate

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Linking the fate of massive black hole binaries to the active galactic nuclei luminosity function

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Massive black hole and gas dynamics in mergers of galaxy nuclei - II. Black hole sinking in star-forming nuclear discs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Massive black holes in merging galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Overlapping inflows as catalysts of AGN activity - II: Relative importance of turbulence and inflow-disc interaction

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quasi-periodicities at year-like timescales in Blazars

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 $γ$-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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Growth and activity of black holes in galaxy mergers with varying mass ratios

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Linking the spin evolution of massive black holes to galaxy kinematics

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$α$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Massive black hole and gas dynamics in galaxy nuclei mergers. I. Numerical implementation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

New insights on the recoiling/binary black hole candidate J0927+2943 via molecular gas observations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nuclear coups: dynamics of black holes in galaxy mergers

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Super-Critical Growth of Massive Black Holes from Stellar-Mass Seeds

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The nature of massive black hole binary candidates: II. Spectral energy distribution atlas

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$μ$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~$μ$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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effects Of Circum-nuclear Disk Gas Evolution And The Spin Of Central Black Holes

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 ε= 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Halpha3: an Halpha imaging survey of HI selected galaxies from ALFALFA. III. Nurture shapes up the Hubble sequence in the Great Wall

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Overlapping Inflow Events as Catalysts for Supermassive Black Hole Growth

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Red-channel (6000-8000 Å) nuclear spectra of 376 local galaxies

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 λ [NII] to Hα lines and the equivalent width (EW) of Hα versus [NII]/Hα (WHAN diagnostic introduced by Cid Fernandes and collaborators) after correcting Hα for underlying absorption. The obtained spectra are made available in machine readable format via the Strasbourg Astronomical Data Center (CDS) and NED.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The nature of massive black hole binary candidates: I. Spectral properties and evolution

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Evolution of binary black holes in self gravitating discs: dissecting the torques

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gravitational Recoil From Accretion-Aligned Black-Hole Binaries

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Observability of Dual Active Galactic Nuclei in Merging Galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Limiting eccentricity of sub-parsec massive black hole binaries surrounded by self-gravitating gas discs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Massive black hole binary eccentricity in rotating stellar systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Massive black hole binary plane reorientation in rotating stellar systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Massive black holes in stellar systems: 'quiescent' accretion and luminosity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Geometry and inclination of the broad line region in blazars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gravitational recoil: effects on massive black hole occupation fraction over cosmic time

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).

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Double quasars: probes of black hole scaling relationships and merger scenarios

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.

0 citations0 reviewsNo code linkedOpen access