Formation Of Supermassive Black Holes Research Articles

Forming supermassive black holes like J1342+0928 (invoking dark matter) in early universe

Recent discovery of J1342+0928 using data from the WISE telescope and ground based surveys indicate presence of a supermassive black hole (SMBH) having a mass of 800 million solar mass at a redshift of about 7.6. This imply that the black hole grew to this mass only 690 million years after the universe started expanding. Here we suggest that formation of such SMBH’s so early in the universe is consistent with our present understanding of the phenomena involved by invoking dark matter (DM).

Intermediate-mass Black Holes and Dark Matter at the Galactic Center

Abstract Could there be a large population of intermediate-mass black holes (IMBHs) formed in the early universe? Whether primordial or formed in Population III, these are likely to be very subdominant compared to the dark matter density, but could seed early dwarf galaxy/globular cluster and supermassive black hole formation. Via survival of dark matter density spikes, we show here that a centrally concentrated relic population of IMBHs, along with ambient dark matter, could account for the Fermi gamma-ray “excess” in the Galactic center because of dark matter particle annihilations.

Dancing to changa: a self-consistent prediction for close SMBH pair formation time-scales following galaxy mergers

We present the first self-consistent prediction for the distribution of formation timescales for close Supermassive Black Hole (SMBH) pairs following galaxy mergers. Using ROMULUS25, the first large-scale cosmological simulation to accurately track the orbital evolution of SMBHs within their host galaxies down to sub-kpc scales, we predict an average formation rate density of close SMBH pairs of 0.013 cMpc^-3 Gyr^-1. We find that it is relatively rare for galaxy mergers to result in the formation of close SMBH pairs with sub-kpc separation and those that do form are often the result of Gyrs of orbital evolution following the galaxy merger. The likelihood and timescale to form a close SMBH pair depends strongly on the mass ratio of the merging galaxies, as well as the presence of dense stellar cores. Low stellar mass ratio mergers with galaxies that lack a dense stellar core are more likely to become tidally disrupted and deposit their SMBH at large radii without any stellar core to aid in their orbital decay, resulting in a population of long-lived 'wandering' SMBHs. Conversely, SMBHs in galaxies that remain embedded within a stellar core form close pairs in much shorter timescales on average. This timescale is a crucial, though often ignored or very simplified, ingredient to models predicting SMBH mergers rates and the connection between SMBH and star formation activity.

Radiation hydrodynamics simulations of the formation of direct-collapse supermassive stellar systems

Formation of supermassive stars (SMSs) with mass ~10^4 Msun is a promising pathway to seed the formation of supermassive black holes in the early universe. The so-called direct-collapse (DC) model postulates that such an SMS forms in a hot gas cloud irradiated by a nearby star-forming galaxy. We study the DC SMS formation in a fully cosmological context using three-dimensional radiation hydrodynamics simulations. We initialize our simulations using the outputs of the cosmological simulation of Chon et al. (2016), where two DC gas clouds are identified. The long-term evolution over a hundred thousand years is followed from the formation of embryo protostars through their growth to SMSs. We show that the strength of the tidal force by a nearby galaxy determines the multiplicity of the formed stars and affects the protostellar growth. In one case, where a collapsing cloud is significantly stretched by strong tidal force, multiple star-disk systems are formed via filament fragmentation. Small-scale fragmentation occurs in each circumstellar disk, and more than 10 stars with masses of a few times 10^3 Msun are finally formed. Interestingly, about a half of them are found as massive binary stars. In the other case, the gas cloud collapses nearly spherically under a relatively weak tidal field, and a single star-disk system is formed. Only a few SMSs with masses ~ 10^4 Msun are found already after evolution of a hundred thousand years, and the SMSs are expected to grow further by gas accretion and to leave massive blackholes at the end of their lives.

Universal relations with fermionic dark matter

We have recently introduced a new model for the distribution of dark matter (DM) in galaxies, the Ruffini-Argüelles-Rueda (RAR) model, based on a self-gravitating system of massive fermions at finite temperatures. The RAR model, for fermion masses above keV, successfully describes the DM halos in galaxies, and predicts the existence of a denser quantum core towards the center of each configuration. We demonstrate here, for the first time, that the introduction of a cutoff in the fermion phase-space distribution, necessary to account for galaxies finite size and mass, defines a new solution with a compact quantum core which represents an alternative to the central black hole (BH) scenario for SgrA*. For a fermion mass in the range 48keV ≤ mc2 ≤ 345keV, the DM halo distribution fulfills the most recent data of the Milky Way rotation curves while harbors a dense quantum core of 4×106M⊙ within the S2 star pericenter. In particular, for a fermion mass of mc2 ∼ 50keV the model is able to explain the DM halos from typical dwarf spheroidal to normal elliptical galaxies, while harboring dark and massive compact objects from ∼ 103M⊙ tp to 108M⊙ at their respective centers. The model is shown to be in good agreement with different observationally inferred universal relations, such as the ones connecting DM halos with supermassive dark central objects. Finally, the model provides a natural mechanism for the formation of supermassive BHs as heavy as few ∼ 108M⊙. We argue that larger BH masses (few ∼ 109−10M⊙) may be achieved by assuming subsequent accretion processes onto the above heavy seeds, depending on accretion efficiency and environment.

Gravitational wave formation from the collapse of dark energy field configurations

Dark energy is the dominant component of the energy density in the Universe. In a previous paper, we have shown that the collapse of dark energy fields leads to the formation of supermassive black holes with masses comparable to the masses of black holes at the centers of galaxies. Thus, it becomes a pressing issue to investigate the other physical consequences of the collapse of dark energy fields. Given that the primary interactions of dark energy fields with the rest of the Universe are gravitational, it is particularly interesting to investigate the gravitational wave signals emitted during the collapse of dark energy fields. This is the focus of the current work described in this paper. We describe and use the 3+1 BSSN formalism to follow the evolution of the dark energy fields coupled with gravity and to extract the gravitational wave signals. Finally, we describe the results of our numerical computations and the gravitational wave signals produced by the collapse of dark energy fields.

Can Superconducting Cosmic Strings Piercing Seed Black Holes Generate Supermassive Black Holes in the Early Universe?

AbstractThe discovery of a large number of supermassive black holes (SMBH) at redshifts , when the Universe was only 900 million years old, raises the question of how such massive compact objects could form in a cosmologically short time interval. Each of the standard scenarios proposed, involving rapid accretion of seed black holes or black hole mergers, faces severe theoretical difficulties in explaining the short‐time formation of supermassive objects. In this work we propose an alternative scenario for the formation of SMBH in the early Universe, in which energy transfer from superconducting cosmic strings piercing small seed black holes is the main physical process leading to rapid mass increase. As a toy model, the accretion rate of a seed black hole pierced by two antipodal strings carrying constant current is considered. Using an effective action approach, which phenomenologically incorporates a large class of superconducting string models, we estimate the minimum current required to form SMBH with masses of order by . This corresponds to the mass of the central black hole powering the quasar ULAS J112001.48+064124.3 and is taken as a test case scenario for early‐epoch SMBH formation. For GUT scale strings, the required fractional increase in the string energy density, due to the presence of the current, is of order 10−7, so that their existence remains consistent with current observational bounds on the string tension. In addition, we consider an “exotic” scenario, in which an SMBH is generated when a small seed black hole is pierced by a higher‐dimensional string, predicted by string theory. We find that both topological defect strings and fundamental strings are able to carry currents large enough to generate early‐epoch SMBH via our proposed mechanism.

Jetted tidal disruptions of stars as a flag of intermediate mass black holes at high redshifts

Tidal disruption events (TDEs) of stars by single or binary supermassive black holes (SMBHs) brighten galactic nuclei and reveal a population of otherwise dormant black holes. Adopting event rates from the literature, we aim to establish general trends in the redshift evolution of the TDE number counts and their observable signals. We pay particular attention to two types of TDEs which are expected to be observable out to high redshifts, namely (i) jetted TDEs whose luminosity is boosted by relativistic beaming, and (ii) TDEs around binary black holes. We show that the brightest (jetted) TDEs are expected to be produced by massive black hole binaries if the occupancy of intermediate mass black holes (IMBHs) in low mass galaxies is high. The same binary population will also provide gravitational wave sources for eLISA. In addition, we find that the shape of the X-ray luminosity function of TDEs strongly depends on the occupancy of IMBHs and could be used to constrain scenarios of SMBH formation. Finally, we make predictions for the expected number of TDEs observed by future X-ray telescopes as a function of their sensitivity limits. We find that an instrument which is 50 times more sensitive than the Burst Alert Telescope (BAT) on board the Swift satellite is expected to trigger ~10 times more events than BAT with 50% of the events coming from z>2. Because of their long decay times, high-redshift TDEs can be mistaken for fixed point sources in deep field surveys such as the 4Ms survey with the Chandra X-ray Observatory and targeted observations of the same deep field with year-long intervals could reveal TDEs. If the occupation fraction of IMBHs is high, 6-20 TDEs are expected in each deep field observed by a telescope 50 times more sensitive than Chandra.

The sustainable growth of the first black holes

Super-Eddington accretion has been suggested as a possible formation pathway of $10^9 \, M_\odot$ supermassive black holes (SMBHs) 800 Myr after the Big Bang. However, stellar feedback from BH seed progenitors and winds from BH accretion disks may decrease BH accretion rates. In this work, we study the impact of these physical processes on the formation of $z \sim 6$ quasar, including new physical prescriptions in the cosmological, data-constrained semi-analytic model GAMETE/QSOdust. We find that the feedback produced by the first stellar progenitors on the surrounding does not play a relevant role in preventing SMBHs formation. In order to grow the $z \gtrsim 6$ SMBHs, the accreted gas must efficiently lose angular momentum. Moreover disk winds, easily originated in super-Eddington accretion regime, can strongly reduce duty cycles. This produces a decrease in the active fraction among the progenitors of $z\sim6$ bright quasars, reducing the probability to observe them.

The Romulus cosmological simulations: a physical approach to the formation, dynamics and accretion models of SMBHs

We present a novel implementation of supermassive black hole (SMBH) formation, dynamics, and accretion in the massively parallel tree+SPH code, ChaNGa. This approach improves the modeling of SMBHs in fully cosmological simulations, allowing for a more de- tailed analysis of SMBH-galaxy co-evolution throughout cosmic time. Our scheme includes novel, physically motivated models for SMBH formation, dynamics and sinking timescales within galaxies, and SMBH accretion of rotationally supported gas. The sub-grid parameters that regulate star formation (SF) and feedback from SMBHs and SNe are optimized against a comprehensive set of z = 0 galaxy scaling relations using a novel, multi-dimensional parameter search. We have incorporated our new SMBH implementation and parameter optimization into a new set of high resolution, large-scale cosmological simulations called Romulus. We present initial results from our flagship simulation, Romulus25, showing that our SMBH model results in SF efficiency, SMBH masses, and global SF and SMBH accretion histories at high redshift that are consistent with observations. We discuss the importance of SMBH physics in shaping the evolution of massive galaxies and show how SMBH feedback is much more effective at regulating star formation compared to SNe feedback in this regime. Further, we show how each aspect of our SMBH model impacts this evolution compared to more common approaches. Finally, we present a science application of this scheme studying the properties and time evolution of an example dual AGN system, highlighting how our approach allows simulations to better study galaxy interactions and SMBH mergers in the context of galaxy-BH co-evolution.

A rumble in the dark: signatures of self-interacting dark matter in supermassive black hole dynamics and galaxy density profiles

We explore for the first time the effect of self-interacting dark matter (SIDM) on the dark matter (DM) and baryonic distribution in massive galaxies formed in hydrodynamical cosmological simulations, including explicit baryonic physics treatment. A novel implementation of Super-Massive Black Hole (SMBH) formation and evolution is used, as in Tremmel et al.(2015, 2016), allowing to explicitly follow SMBH dynamics at the center of galaxies. A high SIDM constant cross-section is chosen, $\sigma$=10 $\rm cm^2/gr$, to amplify differences from CDM models. Milky Way-like galaxies form a shallower DM density profile in SIDM than they do in CDM, with differences already at 20 kpc scales. This demonstrates that even for the most massive spirals the effect of SIDM dominates over the adiabatic contraction due to baryons. Strikingly, the dynamics of SMBHs differs in the SIDM and reference CDM case. SMBHs in massive spirals have sunk to the centre of their host galaxy in both the SIDM and CDM run, while in less massive galaxies about 80$\%$ of the SMBH population is off-centered in the SIDM case, as opposed to the CDM case in which $\sim$90$\%$ of SMBHs have reached their host's centre. SMBHs are found as far as $\sim$9 kpc away from the centre of their host SIDM galaxy. This difference is due to the increased dynamical friction timescale caused by the lower DM density in SIDM galaxies compared to CDM, resulting in 'core stalling'. This pilot work highlights the importance of simulating in a full hydrodynamical context different DM models combined to SMBH physics to study their influence on galaxy formation.

Supermassive black hole seeds: updates on the “quasi-star model”

LISA will allow us to see the infancy of the universe and gather unique constraints on an outstanding mystery: the formation of supermassive black holes. Here, I review a particular scenario, where black hole seeds are formed from direct collapse of gas at the centre of (proto)-galaxies. I will show that very massive (> 105M⊙) black hole seeds are possible but before confidently claim so the impact of physical ingredients, like internal rotation and super Eddington accretion, need to be fully understood.

On the dynamics of supermassive black holes in gas-rich, star-forming galaxies: the case for nuclear star cluster co-evolution

We introduce a new model for the formation and evolution of supermassive black holes (SMBHs) in the ramses code using sink particles, improving over previous work the treatment of gas accretion and dynamical evolution. This new model is tested against a suite of high-resolution simulations of an isolated, gas-rich, cooling halo. We study the effect of various feedback models on the SMBH growth and its dynamics within the galaxy. In runs without any feedback, the SMBH is trapped within a massive bulge and is therefore able to grow quickly, but only if the seed mass is chosen larger than the minimum Jeans mass resolved by the simulation. We demonstrate that, in the absence of supernovae (SN) feedback, the maximum SMBH mass is reached when active galactic nucleus (AGN) heating balances gas cooling in the nuclear region. When our efficient SN feedback is included, it completely prevents bulge formation, so that massive gas clumps can perturb the SMBH orbit, and reduce the accretion rate significantly. To overcome this issue, we propose an observationally motivated model for the joint evolution of the SMBH and a parent nuclear star cluster (NSC), which allows the SMBH to remain in the nuclear region, grow fast and resist external perturbations. In this scenario, however, SN feedback controls the gas supply and the maximum SMBH mass now depends on the balance between AGN heating and gravity. We conclude that SMBH/NSC co-evolution is crucial for the growth of SMBH in high-z galaxies, the progenitors of massive ellipticals today.

An intermediate-mass black hole in the centre of the globular cluster 47 Tucanae.

Intermediate-mass black holes should help us to understand the evolutionary connection between stellar-mass and super-massive black holes. However, the existence of intermediate-mass black holes is still uncertain, and their formation process is therefore unknown. It has long been suspected that black holes with masses 100 to 10,000 times that of the Sun should form and reside in dense stellar systems. Therefore, dedicated observational campaigns have targeted globular clusters for many decades, searching for signatures of these elusive objects. All candidate signatures appear radio-dim and do not have the X-ray to radio flux ratios required for accreting black holes. Based on the lack of an electromagnetic counterpart, upper limits of 2,060 and 470 solar masses have been placed on the mass of a putative black hole in 47 Tucanae (NGC 104) from radio and X-ray observations, respectively. Here we show there is evidence for a central black hole in 47 Tucanae with a mass of solar masses when the dynamical state of the globular cluster is probed with pulsars. The existence of an intermediate-mass black hole in the centre of one of the densest clusters with no detectable electromagnetic counterpart suggests that the black hole is not accreting at a sufficient rate to make it electromagnetically bright and therefore, contrary to expectations, is gas-starved. This intermediate-mass black hole might be a member of an electromagnetically invisible population of black holes that grow into supermassive black holes in galaxies.

Supermassive black holes from collapsing dark matter Bose–Einstein condensates

The discovery of active galactic nuclei at redshifts 6 suggests that supermassive black holes (SMBHs) formed early on. Growth of the remnants of population III stars by accretion of matter, both baryonic as well as collisionless dark matter (DM), leading up to formation of SMBHs is a very slow process. Therefore, such models encounter difficulties in explaining quasars detected at . Furthermore, massive particles making up collisionless DM have not only so far eluded experimental detection but they also do not satisfactorily explain gravitational structures on small scales.In recent years, there has been a surge in research activities concerning cosmological structure formation that involve coherent, ultra-light bosons in a dark fluid-like or fuzzy cold DM state. In this paper, we study collapse of such ultra-light bosonic halo DM that is in a Bose–Einstein condensate (BEC) phase to give rise to SMBHs on dynamical time scales. Time evolution of such self-gravitating BECs is examined using the Gross–Pitaevskii equation in the framework of time-dependent variational method. Comprised of identical dark bosons of mass m, BECs can collapse to form black holes of mass Meff on time scales ∼108 yrs provided . In particular, ultra-light dark bosons of mass can lead to SMBHs with mass at . Recently observed radio-galaxies in the ELAIS-N1 deep field with aligned jets can also possibly be explained if vortices of a rotating cluster size BEC collapse to form spinning SMBHs with angular momentum , where nW and M are the winding number and mass of a vortex, respectively.

Faint progenitors of luminousz ∼ 6 quasars: Why do not we see them?

Observational searches for faint active nuclei at $z > 6$ have been extremely elusive, with a few candidates whose high-$z$ nature is still to be confirmed. Interpreting this lack of detections is crucial to improve our understanding of high-$z$ supermassive black holes (SMBHs) formation and growth. In this work, we present a model for the emission of accreting BHs in the X-ray band, taking into account super-Eddington accretion, which can be very common in gas-rich systems at high-$z$. We compute the spectral energy distribution for a sample of active galaxies simulated in a cosmological context, which represent the progenitors of a $z \sim 6$ SMBH with $M_{\rm BH} \sim 10^9 \, M_\odot$. We find an average Compton thick fraction of $\sim 45\%$ and large typical column densities ($N_H \gtrsim 10^{23} \rm \, cm^2$). However, faint progenitors are still luminous enough to be detected in the X-ray band of current surveys. Even accounting for a maximum obscuration effect, the number of detectable BHs is reduced at most by a factor 2. In our simulated sample, observations of faint quasars are mainly limited by their very low active fraction ($f_{\rm act} \sim 1 \%$), which is the result of short, super-critical growth episodes. We suggest that to detect high-$z$ SMBHs progenitors, large area surveys with shallower sensitivities, such as Cosmos Legacy and XMM-LSS+XXL, are to be preferred with respect to deep surveys probing smaller fields, such as CDF-S.

Gravitational collapse of dark energy field configurations and supermassive black hole formation

Dark energy is the dominant component of the total energy density of our Universe. The primary interaction of dark energy with the rest of the Universe is gravitational. It is therefore important to understand the gravitational dynamics of dark energy. Since dark energy is a low-energy phenomenon from the perspective of particle physics and field theory, a fundamental approach based on fields in curved space should be sufficient to understand the current dynamics of dark energy. Here, we take a field theory approach to dark energy. We discuss the evolution equations for a generic dark energy field in curved space-time and then discuss the gravitational collapse for dark energy field configurations. We describe the 3 + 1 BSSN formalism to study the gravitational collapse of fields for any general potential for the fields and apply this formalism to models of dark energy motivated by particle physics considerations. We solve the resulting equations for the time evolution of field configurations and the dynamics of space-time. Our results show that gravitational collapse of dark energy field configurations occurs and must be considered in any complete picture of our Universe. We also demonstrate the black hole formation as a result of the gravitational collapse of the dark energy field configurations. The black holes produced by the collapse of dark energy fields are in the supermassive black hole category with the masses of these black holes being comparable to the masses of black holes at the centers of galaxies.

NEW CONSTRAINTS ON THE MOLECULAR GAS IN THE PROTOTYPICAL HyLIRGs BRI 1202–0725 AND BRI 1335–0417

ABSTRACT We present Karl G Jansky Very Large Array observations of CO( ) line emission and rest-frame 250 GHz continuum emission of the Hyper-Luminous IR Galaxies (HyLIRGs) BRI 1202–0725 (z = 4.69) and BRI 1335–0417 (z = 4.41), with an angular resolution as high as 0.″15. Our low-order CO observations delineate the cool molecular gas, the fuel for star formation in the systems, in unprecedented detail. For BRI 1202–0725, line emission is seen from both extreme starburst galaxies: the quasar host and the optically obscured submm galaxy (SMG), in addition to one of the Lyα emitting galaxies in the group. Line emission from the SMG shows an east–west extension of about 0.″6. For Lyα-2, the CO emission is detected at the same velocity as [C ii] and [N ii], indicating a total gas mass ∼4.0 × 1010 M ⊙. The CO emission from BRI 1335–0417 peaks at the nominal quasar position, with a prominent northern extension (∼1″, a possible tidal feature). The gas depletion timescales are ∼107 years for the three HyLIRGs, consistent with extreme starbursts, while that of Lyα-2 may be consistent with main sequence galaxies. We interpret these sources as major star formation episodes in the formation of massive galaxies and supermassive black holes via gas-rich mergers in the early universe.

RICH KOZAI–LIDOV DYNAMICS IN AN INITIALLY THIN AND ECCENTRIC STELLAR DISK AROUND A SUPERMASSIVE BLACK HOLE

ABSTRACT There is growing evidence of star formation in the vicinity of supermassive black holes (SMBHs) in galactic nuclei. A viable scenario for this process assumes infall of a massive gas cloud toward the SMBH and subsequent formation of a dense accretion disk, which gives birth to the young stars. Numerical hydrodynamical models indicate that this star formation process is rather fast and precedes full circularization of the accretion flow, i.e., the new stars are born on elliptic orbits. By means of direct numerical N-body modeling, we show in this paper that the nonzero eccentricity of the stellar disks around the SMBH leads to an onset of various types of the Kozai–Lidov oscillations of a non-negligible subset of individual orbits in the disk, showing a remarkable robustness of this classical mechanism. Among others, we demonstrate that under certain circumstances, the presence of an additional spherical cluster (which is generally known to damp Kozai–Lidov oscillations) may trigger such oscillations as a result of affecting the internal flow of the angular momentum through the disk. We conclude that the Kozai–Lidov oscillations are capable of substantially modifying the initial structure of the disk (its thickness and distribution of eccentricities, in particular).

The New Numerical Galaxy Catalog (ν2GC): An updated semi-analytic model of galaxy and active galactic nucleus formation with large cosmological N-body simulations

Abstract We present a new cosmological galaxy formation model, ν2GC, as an updated version of our previous model νGC. We adopt the so-called “semi-analytic” approach, in which the formation history of dark matter halos is computed by N-body simulations, while the baryon physics such as gas cooling, star formation, and supernova feedback are simply modeled by phenomenological equations. Major updates of the model are as follows: (1) the merger trees of dark matter halos are constructed in state-of-the-art N-body simulations, (2) we introduce the formation and evolution process of supermassive black holes and the suppression of gas cooling due to active galactic nucleus (AGN) activity, (3) we include heating of the intergalactic gas by the cosmic UV background, and (4) we tune some free parameters related to the astrophysical processes using a Markov chain Monte Carlo method. Our N-body simulations of dark matter halos have unprecedented box size and mass resolution (the largest simulation contains 550 billion particles in a 1.12 Gpc h−1 box), enabling the study of much smaller and rarer objects. The model was tuned to fit the luminosity functions of local galaxies and mass function of neutral hydrogen. Local observations, such as the Tully–Fisher relation, the size–magnitude relation of spiral galaxies, and the scaling relation between the bulge mass and black hole mass were well reproduced by the model. Moreover, the model also reproduced well the cosmic star formation history and redshift evolution of rest-frame K-band luminosity functions. The numerical catalog of the simulated galaxies and AGNs is publicly available on the web.