scholarly journals The Black Hole Mass Function Across Cosmic Times. I. Stellar Black Holes and Light Seed Distribution

2022 ◽  
Vol 924 (2) ◽  
pp. 56
Author(s):  
Alex Sicilia ◽  
Andrea Lapi ◽  
Lumen Boco ◽  
Mario Spera ◽  
Ugo N. Di Carlo ◽  
...  

Abstract This is the first paper in a series aimed at modeling the black hole (BH) mass function, from the stellar to the intermediate to the (super)massive regime. In the present work, we focus on stellar BHs and provide an ab initio computation of their mass function across cosmic times; we mainly consider the standard, and likely dominant, production channel of stellar-mass BHs constituted by isolated single/binary star evolution. Specifically, we exploit the state-of-the-art stellar and binary evolutionary code SEVN, and couple its outputs with redshift-dependent galaxy statistics and empirical scaling relations involving galaxy metallicity, star formation rate and stellar mass. The resulting relic mass function dN / dVd log m • as a function of the BH mass m • features a rather flat shape up to m • ≈ 50 M ⊙ and then a log-normal decline for larger masses, while its overall normalization at a given mass increases with decreasing redshift. We highlight the contribution to the local mass function from isolated stars evolving into BHs and from binary stellar systems ending up in single or binary BHs. We also include the distortion on the mass function induced by binary BH mergers, finding that it has a minor effect at the high-mass end. We estimate a local stellar BH relic mass density of ρ • ≈ 5 × 107 M ⊙ Mpc−3, which exceeds by more than two orders of magnitude that in supermassive BHs; this translates into an energy density parameter Ω• ≈ 4 × 10−4, implying that the total mass in stellar BHs amounts to ≲1% of the local baryonic matter. We show how our mass function for merging BH binaries compares with the recent estimates from gravitational wave observations by LIGO/Virgo, and discuss the possible implications for dynamical formation of BH binaries in dense environments like star clusters. We address the impact of adopting different binary stellar evolution codes (SEVN and COSMIC) on the mass function, and find the main differences to occur at the high-mass end, in connection with the numerical treatment of stellar binary evolution effects. We highlight that our results can provide a firm theoretical basis for a physically motivated light seed distribution at high redshift, to be implemented in semi-analytic and numerical models of BH formation and evolution. Finally, we stress that the present work can constitute a starting point to investigate the origin of heavy seeds and the growth of (super)massive BHs in high-redshift star-forming galaxies, that we will pursue in forthcoming papers.

2019 ◽  
Vol 485 (3) ◽  
pp. 4413-4422 ◽  
Author(s):  
Daniel J D’Orazio ◽  
Abraham Loeb ◽  
James Guillochon

ABSTRACT The rate of tidal disruption flares (TDFs) per mass of the disrupting black hole encodes information on the present-day mass function (PDMF) of stars in the clusters surrounding super massive black holes. We explore how the shape of the TDF rate with black hole mass can constrain the PDMF, with only weak dependence on black hole spin. We show that existing data can marginally constrain the minimum and maximum masses of stars in the cluster, and the high-mass end of the PDMF slope, as well as the overall TDF rate. With $\mathcal {O}(100)$ TDFs expected to be identified with the Zwicky Transient Facility, the overall rate can be highly constrained, but still with only marginal constraints on the PDMF. However, if ${\lesssim } 10 {{\ \rm per\ cent}}$ of the TDFs expected to be found by LSST over a decade ($\mathcal {O}(10^3)$ TDFs) are identified, then precise and accurate estimates can be made for the minimum stellar mass (within a factor of 2) and the average slope of the high-mass PDMF (to within $\mathcal {O}(10{{\ \rm per\ cent}})$) in nuclear star clusters. This technique could be adapted in the future to probe, in addition to the PDMF, the local black hole mass function and possibly the massive black hole binary population.


2019 ◽  
Vol 489 (1) ◽  
pp. 487-496 ◽  
Author(s):  
Boyan K Stoychev ◽  
Keri L Dixon ◽  
Andrea V Macciò ◽  
Marvin Blank ◽  
Aaron A Dutton

ABSTRACT We use 38 high-resolution simulations of galaxy formation between redshift 10 and 5 to study the impact of a 3 keV warm dark matter (WDM) candidate on the high-redshift Universe. We focus our attention on the stellar mass function and the global star formation rate and consider the consequences for reionization, namely the neutral hydrogen fraction evolution and the electron scattering optical depth. We find that three different effects contribute to differentiate warm and cold dark matter (CDM) predictions: WDM suppresses the number of haloes with mass less than few 109 M⊙; at a fixed halo mass, WDM produces fewer stars than CDM, and finally at halo masses below 109 M⊙, WDM has a larger fraction of dark haloes than CDM post-reionization. These three effects combine to produce a lower stellar mass function in WDM for galaxies with stellar masses at and below 107 M⊙. For z > 7, the global star formation density is lower by a factor of two in the WDM scenario, and for a fixed escape fraction, the fraction of neutral hydrogen is higher by 0.3 at z ∼ 6. This latter quantity can be partially reconciled with CDM and observations only by increasing the escape fraction from 23 per cent to 34 per cent. Overall, our study shows that galaxy formation simulations at high redshift are a key tool to differentiate between dark matter candidates given a model for baryonic physics.


2020 ◽  
Vol 495 (3) ◽  
pp. 3065-3078 ◽  
Author(s):  
Pratika Dayal ◽  
Marta Volonteri ◽  
Tirthankar Roy Choudhury ◽  
Raffaella Schneider ◽  
Maxime Trebitsch ◽  
...  

ABSTRACT In this work we investigate the properties of the sources that reionized the intergalactic medium (IGM) in the high-redshift Universe. Using a semi-analytical model aimed at reproducing galaxies and black holes in the first ∼1.5 Gyr of the Universe, we revisit the relative role of star formation and black hole accretion in producing ionizing photons that can escape into the IGM. Both star formation and black hole accretion are regulated by supernova feedback, resulting in black hole accretion being stunted in low-mass haloes. We explore a wide range of combinations for the escape fraction of ionizing photons (redshift-dependent, constant, and scaling with stellar mass) from both star formation ($\langle f_{\rm esc}^{\rm sf} \rangle$) and AGN ($f_{\rm esc}^{\rm bh}$) to find: (i) the ionizing budget is dominated by stellar radiation from low stellar mass ($M_*\lt 10^9 \, {\rm \rm M_\odot }$) galaxies at z > 6 with the AGN contribution (driven by $M_{bh}\gt 10^6 \, {\rm \rm M_\odot }$ black holes in $M_* \gtrsim 10^9\, {\rm \rm M_\odot }$ galaxies) dominating at lower redshifts; (ii) AGN only contribute $10-25{{\ \rm per\ cent}}$ to the cumulative ionizing emissivity by z = 4 for the models that match the observed reionization constraints; (iii) if the stellar mass dependence of $\langle f_{\rm esc}^{\rm sf} \rangle$ is shallower than $f_{\rm esc}^{\rm bh}$, at z < 7 a transition stellar mass exists above which AGN dominate the escaping ionizing photon production rate; (iv) the transition stellar mass decreases with decreasing redshift. While AGN dominate the escaping emissivity above the knee of the stellar mass function at z ∼ 6.8, they take-over at stellar masses that are a tenth of the knee mass by z = 4.


2020 ◽  
Vol 501 (2) ◽  
pp. 1568-1590
Author(s):  
Lukas J Furtak ◽  
Hakim Atek ◽  
Matthew D Lehnert ◽  
Jacopo Chevallard ◽  
Stéphane Charlot

ABSTRACT We present new measurements of the very low mass end of the galaxy stellar mass function (GSMF) at z ∼ 6−7 computed from a rest-frame ultraviolet selected sample of dropout galaxies. These galaxies lie behind the six Hubble Frontier Field clusters and are all gravitationally magnified. Using deep Spitzer/IRAC and Hubble Space Telescope imaging, we derive stellar masses by fitting galaxy spectral energy distributions and explore the impact of different model assumptions and parameter degeneracies on the resulting GSMF. Our sample probes stellar masses down to $M_{\star }\gt 10^{6}\, \text{M}_{\odot}$ and we find the z ∼ 6−7 GSMF to be best parametrized by a modified Schechter function that allows for a turnover at very low masses. Using a Monte Carlo Markov chain analysis of the GSMF, including accurate treatment of lensing uncertainties, we obtain a relatively steep low-mass end slope $\alpha \simeq -1.96_{-0.08}^{+0.09}$ and a turnover at $\log (M_T/\text{M}_{\odot})\simeq 7.10_{-0.56}^{+0.17}$ with a curvature of $\beta \simeq 1.00_{-0.73}^{+0.87}$ for our minimum assumption model with constant star formation history (SFH) and low dust attenuation, AV ≤ 0.2. We find that the z ∼ 6−7 GSMF, in particular its very low mass end, is significantly affected by the assumed functional form of the star formation history and the degeneracy between stellar mass and dust attenuation. For example, the low-mass end slope ranges from $\alpha \simeq -1.82_{-0.07}^{+0.08}$ for an exponentially rising SFH to $\alpha \simeq -2.34_{-0.10}^{+0.11}$ when allowing AV of up to 3.25. Future observations at longer wavelengths and higher angular resolution with the James Webb Space Telescope are required to break these degeneracies and to robustly constrain the stellar mass of galaxies on the extreme low-mass end of the GSMF.


2016 ◽  
Vol 458 (4) ◽  
pp. 3478-3478 ◽  
Author(s):  
Alice Mortlock ◽  
Christopher. J. Conselice ◽  
William G. Hartley ◽  
Ken Duncan ◽  
Caterina Lani ◽  
...  

2021 ◽  
Vol 503 (4) ◽  
pp. 5115-5133
Author(s):  
A A Khostovan ◽  
S Malhotra ◽  
J E Rhoads ◽  
S Harish ◽  
C Jiang ◽  
...  

ABSTRACT The H α equivalent width (EW) is an observational proxy for specific star formation rate (sSFR) and a tracer of episodic, bursty star-formation activity. Previous assessments show that the H α EW strongly anticorrelates with stellar mass as M−0.25 similar to the sSFR – stellar mass relation. However, such a correlation could be driven or even formed by selection effects. In this study, we investigate how H α EW distributions correlate with physical properties of galaxies and how selection biases could alter such correlations using a z = 0.47 narrow-band-selected sample of 1572 H α emitters from the Ly α Galaxies in the Epoch of Reionization (LAGER) survey as our observational case study. The sample covers a 3 deg2 area of COSMOS with a survey comoving volume of 1.1 × 105 Mpc3. We assume an intrinsic EW distribution to form mock samples of H α emitters and propagate the selection criteria to match observations, giving us control on how selection biases can affect the underlying results. We find that H α EW intrinsically correlates with stellar mass as W0∝M−0.16 ± 0.03 and decreases by a factor of ∼3 from 107 M⊙ to 1010 M⊙, while not correcting for selection effects steepens the correlation as M−0.25 ± 0.04. We find low-mass H α emitters to be ∼320 times more likely to have rest-frame EW>200 Å compared to high-mass H α emitters. Combining the intrinsic W0–stellar mass correlation with an observed stellar mass function correctly reproduces the observed H α luminosity function, while not correcting for selection effects underestimates the number of bright emitters. This suggests that the W0–stellar mass correlation when corrected for selection effects is physically significant and reproduces three statistical distributions of galaxy populations (line luminosity function, stellar mass function, EW distribution). At lower stellar masses, we find there are more high-EW outliers compared to high stellar masses, even after we take into account selection effects. Our results suggest that high sSFR outliers indicative of bursty star formation activity are intrinsically more prevalent in low-mass H α emitters and not a byproduct of selection effects.


2019 ◽  
Vol 15 (S352) ◽  
pp. 234-238
Author(s):  
Donatella Romano ◽  
Zhi-Yu Zhang ◽  
Francesca Matteucci ◽  
Rob J. Ivison ◽  
Padelis P. Papadopoulos

AbstractDetermining the shape of the stellar initial mass function (IMF) and whether it is constant or varies in space and time is the Holy Grail of modern astrophysics, with profound implications for all theories of star and galaxy formation. On a theoretical ground, the extreme conditions for star formation (SF) encountered in the most powerful starbursts in the Universe are expected to favour the formation of massive stars. Direct methods of IMF determination, however, cannot probe such systems, because of the severe dust obscuration affecting their starlight. The next best option is to observe CNO bearing molecules in the interstellar medium at millimetre/ submillimetre wavelengths, which, in principle, provides the best indirect evidence for IMF variations. In this contribution, we present our recent findings on this issue. First, we reassess the roles of different types of stars in the production of CNO isotopes. Then, we calibrate a proprietary chemical evolution code using Milky Way data from the literature, and extend it to discuss extragalactic data. We show that, though significant uncertainties still hamper our knowledge of the evolution of CNO isotopes in galaxies, compelling evidence for an IMF skewed towards high-mass stars can be found for galaxy-wide starbursts. In particular, we analyse a sample of submillimetre galaxies observed by us with the Atacama Large Millimetre Array at the peak of the SF activity of the Universe, for which we measure 13C/18O⋍1. This isotope ratio is especially sensitive to IMF variations, and is little affected by observational uncertainties. At the end, ongoing developments of our work are briefly outlined.


2019 ◽  
Vol 492 (2) ◽  
pp. 1706-1712
Author(s):  
Anton Vikaeus ◽  
Erik Zackrisson ◽  
Christian Binggeli

ABSTRACT The upcoming James Webb Space Telescope (JWST) will allow observations of high-redshift galaxies at fainter detection levels than ever before, and JWST surveys targeting gravitationally lensed fields are expected to bring z ≳ 6 objects with very low star formation rate (SFR) within reach of spectroscopic studies. As galaxies at lower and lower star formation activity are brought into view, many of the standard methods used in the analysis of integrated galaxy spectra are at some point bound to break down, due to violation of the assumptions of a well-sampled stellar initial mass function (IMF) and a slowly varying SFR. We argue that galaxies with SFR ∼ 0.1 M⊙ yr−1 are likely to turn up at the spectroscopic detection limit of JWST in lensed fields, and investigate to what extent star formation sampling may affect the spectral analysis of such objects. We use the slug spectral synthesis code to demonstrate that such effects are likely to have significant impacts on spectral diagnostics of, for example, the Balmer emission lines. These effects are found to stem primarily from SFRs varying rapidly on short (∼Myr) time-scales due to star formation in finite units (star clusters), whereas the effects of an undersampled IMF is deemed insignificant in comparison. In contrast, the ratio between the He ii- and H i-ionizing flux is found to be sensitive to IMF-sampling as well as ICMF-sampling (sampling of the initial cluster mass function), which may affect interpretations of galaxies containing Population III stars or other sources of hard ionizing radiation.


Author(s):  
Aldo Rodríguez-Puebla ◽  
A. R. Calette ◽  
Vladimir Avila-Reese ◽  
Vicente Rodriguez-Gomez ◽  
Marc Huertas-Company

Abstract We report the bivariate $\rm HI$ - and $\rm H_{2}$ -stellar mass distributions of local galaxies in addition of an inventory of galaxy mass functions, MFs, for $\rm HI$ , $\rm H_{2}$ , cold gas, and baryonic mass, separately into early- and late-type galaxies. The MFs are determined using the $\rm HI$ and $\rm H_{2}$ conditional distributions and the galaxy stellar mass function (GSMF). For the conditional distributions we use the results from the compilation presented in Calette et al. [(2018) RMxAA, 54, 443.]. For determining the GSMF from $M_{*}\sim3\times10^{7}$ to $3\times10^{12}\ \text{M}_{\odot}$ , we combine two spectroscopic samples from the Sloan Digital Sky Survey at the redshift range $0.0033<z<0.2$ . We find that the low-mass end slope of the GSMF, after correcting from surface brightness incompleteness, is $\alpha\approx-1.4$ , consistent with previous determinations. The obtained $\rm HI\,$ MFs agree with radio blind surveys. Similarly, the $\rm H_{2}\,$ MFs are consistent with CO follow-up optically-selected samples. We estimate the impact of systematics due to mass-to-light ratios and find that our MFs are robust against systematic errors. We deconvolve our MFs from random errors to obtain the intrinsic MFs. Using the MFs, we calculate cosmic density parameters of all the baryonic components. Baryons locked inside galaxies represent 5.4% of the universal baryon content, while $\sim\! 96\%$ of the $\rm HI$ and $\rm H_{2}$ mass inside galaxies reside in late-type morphologies. Our results imply cosmic depletion times of $\rm H_{2}$ and total neutral H in late-type galaxies of $\sim\!1.3$ and 7.2 Gyr, respectively, which shows that late type galaxies are on average inefficient in converting $\rm H_{2}$ into stars and in transforming $\rm HI$ gas into $\rm H_{2}$ . Our results provide a fully self-consistent empirical description of galaxy demographics in terms of the bivariate gas–stellar mass distribution and their projections, the MFs. This description is ideal to compare and/or to constrain galaxy formation models.


2019 ◽  
Vol 627 ◽  
pp. A131 ◽  
Author(s):  
M. Cousin ◽  
P. Guillard ◽  
M. D. Lehnert

Context. Star formation in galaxies is inefficient, and understanding how star formation is regulated in galaxies is one of the most fundamental challenges of contemporary astrophysics. Radiative cooling, feedback from supernovae and active galactic nuclei (AGN), and large-scale dynamics and dissipation of turbulent energy act over various time and spatial scales and all regulate star formation in a complex gas cycle. Aims. This paper presents the physics implemented in a new semi-analytical model of galaxy formation and evolution called the Galaxy Assembler from dark-matter Simulation (G.A.S.). Methods. The fundamental underpinning of our new model is the development of a multiphase interstellar medium (ISM) in which energy produced by supernovae and AGN maintains an equilibrium between a diffuse, hot, and stable gas and a cooler, clumpy, and low-volume filling factor gas. The hot gas is susceptible to thermal and dynamical instabilities. We include a description of how turbulence leads to the formation of giant molecular clouds through an inertial turbulent energy cascade, assuming a constant kinetic energy transfer per unit volume. We explicitly modelled the evolution of the velocity dispersion at different scales of the cascade and accounted for thermal instabilities in the hot halo gas. Thermal instabilities effectively reduce the impact of radiative cooling and moderates accretion rates onto galaxies, and in particular, for those residing in massive haloes. Results. We show that rapid and multiple exchanges between diffuse and unstable gas phases strongly regulates star formation rates in galaxies because only a small fraction of the unstable gas is forming stars. We checked that the characteristic timescales describing the gas cycle, gas depletion timescale, and star-forming laws at different scales are in good agreement with observations. For high-mass haloes and galaxies, cooling is naturally regulated by the growth of thermal instabilities, so we do not need to implement strong AGN feedback in this model. Our results are also in good agreement with the observed stellar mass function from z ≃ 6.0 to z ≃ 0.5. Conclusion. Our model offers the flexibility to test the impact of various physical processes on the regulation of star formation on a representative population of galaxies across cosmic times. Thermal instabilities and the cascade of turbulent energy in the dense gas phase introduce a delay between gas accretion and star formation, which keeps galaxy growth inefficient in the early Universe. The main results presented in this paper, such as stellar mass functions, are available in the GALAKSIENN library.


Sign in / Sign up

Export Citation Format

Share Document