Abstract
The dynamics of black hole (BH) seeds in high-redshift galaxies is key to understand their ability to grow via accretion and to pair in close binaries during galactic mergers. To properly ...follow the dynamics of BHs we develop a physically motivated model to capture unresolved dynamical friction from stars, dark matter, and gas. We first validate the model and then we use it to investigate the dynamics of seed BHs born at z ∼ 9 in dwarf proto-galaxies. We perform a suite of zoom cosmological simulations with spatial resolution as high as 10 pc and with a stellar and dark matter mass resolution of $2\times 10^3 \, \, $ and $2\times 10^5 \, \, \mathrm{ M}_{\odot }$, respectively. We first explore the dynamics of a seed BH in the galaxy where it is born and show that it is highly erratic if the seed mass is less than $10^5\, \, \mathrm{ M}_{\odot }$. The dynamics is dominated by the stellar component, whose distribution is irregular and patchy, thus inducing stochasticity in the orbits: the BH may be anywhere in the proto-galaxy. When this dwarf merges into a larger galaxy, it is paramount to simulate the process with very high spatial and mass resolution in order to correctly account for the stripping of the stellar envelope of the satellite BH. The outcome of the encounter could be either a tight binary or, at least temporary, a wandering BH, leading to multiple BHs in a galaxy, each inherited from a different merger.
We consider super-critical accretion with angular momentum onto stellar-mass black holes as a possible mechanism for growing billion-solar-mass black 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 black 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.
ABSTRACT
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 multicomponent 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 sech2 and an exponential vertical profile) and an accurate prescription for the dynamical friction experienced by massive perturbers (MP) 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 MP orbits and its computational efficiency make our framework perfectly suited for large parameter-space exploration studies.
ABSTRACT
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 ...BHs (wBHs) by using the semi-analytical model L-Galaxies 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 pseudo-bulges, classical bulges or ellipticals, being $\overline{a} \sim 0.9$, 0.7 and 0.4, respectively. The outcomes of the model display a good consistency with $z \le 4$ quasar luminosity functions and the $z = 0$ BH mass function, spin values, and BH 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 $\rm {\sim}2\, dex$ smaller than the one of nuclear BHs. At all redshifts, wBHs are typically hosted in $\rm {\it M}_{halo} \gtrsim 10^{13} \, M_{\odot }$ and $\rm {\it M}_{stellar} \gtrsim 10^{10} \, M_{\odot }$, being orphan wBHs the dominant type. Besides, independently of redshift and halo mass, ejected wBHs inhabit the central regions (${\lesssim}\rm 0.3{\it R}_{200}$) of the host DM halo, while orphan wBH linger at larger scales (${\gtrsim}\rm 0.5{\it R}_{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 BH–bulge relation, in particular for BHs in pseudo-bulges, for which the relation is flattened at $\rm {\it M}_{bulge} \gt 10^{10.2}\, M_{\odot }$ and the scatter increase up to ${\sim}\rm 3\, dex$.
The birth of a supermassive black hole binary Pfister, Hugo; Lupi, Alessandro; Capelo, Pedro R ...
Monthly notices of the Royal Astronomical Society,
11/2017, Letnik:
471, Številka:
3
Journal Article
Recenzirano
Odprti dostop
Abstract
We study the dynamical evolution of supermassive black holes, in the late stage of galaxy mergers, from kpc to pc scales. In particular, we capture the formation of the binary, a necessary ...step before the final coalescence, and trace back the main processes causing the decay of the orbit. We use hydrodynamical simulations of galaxy mergers with different resolutions, from 20 pc down to 1 pc, in order to study the effects of the resolution on our results, remove numerical effects, and assess that resolving the influence radius of the orbiting black hole is a minimum condition to fully capture the formation of the binary. Our simulations include the relevant physical processes, namely star formation, supernova feedback, accretion on to the black holes and the ensuing feedback. We find that, in these mergers, dynamical friction from the smooth stellar component of the nucleus is the main process that drives black holes from kpc to pc scales. Gas does not play a crucial role and even clumps do not induce scattering or perturb the orbits. We compare the time needed for the formation of the binary to analytical predictions and suggest how to apply such analytical formalism to obtain estimates of binary formation times in lower resolution simulations.
Abstract
We present results from the ‘Ponos’ simulation suite on the early evolution of a massive, Mvir(z = 0) = 1.2 × 1013 M⊙ galaxy. At z ≳ 6, before feedback from a central supermassive black hole ...becomes dominant, the main galaxy has a stellar mass ∼2 × 109 M⊙ and a star formation rate ∼20 M⊙ yr−1. The galaxy sits near the expected main sequence of star-forming galaxies at those redshifts, and resembles moderately star-forming systems observed at z > 5. The high specific star formation rate results in vigorous heating and stirring of the gas by supernovae feedback, and the galaxy develops a thick and turbulent disc, with gas velocity dispersion ∼40 km s−1, rotation to dispersion ratio ∼2, and with a significant amount of gas at ∼105 K. The Toomre parameter always exceeds the critical value for gravito-turbulence, Q ∼ 1.5–2, mainly due to the contribution of warm/hot gas inside the disc. Without feedback, a nearly gravito-turbulent regime establishes with similar gas velocity dispersion and lower Q. We propose that the ‘hot and turbulent’ disc regime seen in our simulations, unlike the ‘cold and turbulent’ gravito-turbulent regime of massive clumpy disc galaxies at z ∼ 1–2, is a fundamental characterization of the main-sequence galaxies at z ≳ 6, as they can sustain star formation rates comparable to those of low-mass starbursts at z = 0. This results in no sustained coherent gas inflows through the disc, and in fluctuating and anisotropic mass transport, possibly postponing the assembly of the bulge and causing the initial feeding of the central black hole to be highly intermittent.
ABSTRACT
Active galactic nuclei (AGNs) 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 parallel either 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.
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) time-scales, temporal behaviour, and relative magnitude. We find that (i) BHAR and galaxy-wide SFR are typically temporally uncorrelated, and have different variability time-scales, 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, because of the long time spent in this phase. (iv) The time-scale 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.
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 O iii, and strong Fe ii line intensities. The mass MBH 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 on to the MBH larger than the average value appropriate to Seyfert 1 galaxies. In this Letter, we show that NLS1 peculiar MBH and L/LEdd 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 full width at half-maximum 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/LEdd is (on average) comparable to the value of the more common broad-line Seyfert 1 galaxies.