skip to main content

Title: Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). VIII. A less biased view of the early co-evolution of black holes and host galaxies

We present ALMA [C ii] line and far-infrared (FIR) continuum observations of three $z \gt 6$ low-luminosity quasars ($M_{\rm 1450} \gt -25$ mag) discovered by our Subaru Hyper Suprime-Cam (HSC) survey. The [C ii] line was detected in all three targets with luminosities of $(2.4\mbox{--}9.5) \times 10^8\, L_{\odot }$, about one order of magnitude smaller than optically luminous ($M_{\rm 1450} \lesssim -25$ mag) quasars. The FIR continuum luminosities range from $\lt 9 \times 10^{10}\, L_{\odot }$ (3 $\sigma$ limit) to ${\sim } 2 \times 10^{12}\, L_{\odot }$, indicating a wide range in star formation rates in these galaxies. Most of the HSC quasars studied thus far show [C ii]/ FIR luminosity ratios similar to local star-forming galaxies. Using the [C ii]-based dynamical mass ($M_{\rm dyn}$) as a surrogate for bulge stellar mass ($M_{\rm\, bulge}$), we find that a significant fraction of low-luminosity quasars are located on or even below the local $M_{\rm\, BH}$–$M_{\rm\, bulge}$ relation, particularly at the massive end of the galaxy mass distribution. In contrast, previous studies of optically luminous quasars have found that black holes are overmassive relative to the local relation. Given the low luminosities of our targets, we are exploring the nature of the early co-evolution of supermassive black holes and more » their hosts in a less biased way. Almost all of the quasars presented in this work are growing their black hole mass at a much higher pace at $z \sim 6$ than the parallel growth model, in which supermassive black holes and their hosts grow simultaneously to match the local $M_{\rm\, BH}$–$M_{\rm\, bulge}$ relation at all redshifts. As the low-luminosity quasars appear to realize the local co-evolutionary relation even at $z \sim 6$, they should have experienced vigorous starbursts prior to the currently observed quasar phase to catch up with the relation.

« less
 ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  more » ;  ;  ;  ;  ;  ;  ;   « less
Publication Date:
Journal Name:
Publications of the Astronomical Society of Japan
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT We report the serendipitous discovery of a dust-obscured galaxy observed as part of the Atacama Large Millimeter Array (ALMA) Large Program to INvestigate [C ii] at Early times (ALPINE). While this galaxy is detected both in line and continuum emissions in ALMA Band 7, it is completely dark in the observed optical/near-infrared bands and only shows a significant detection in the UltraVISTA Ks band. We discuss the nature of the observed ALMA line, that is [C ii] at $z$ ∼ 4.6 or high-J CO transitions at $z$ ∼ 2.2. In the first case, we find a [C ii]/FIR luminosity ratio of $\mathrm{log}{(L_{[\mathrm{ C}\, \rm {\small {II}}]}/L_{\mathrm{ FIR}})} \sim -2.5$, consistent with the average value for local star-forming galaxies (SFGs). In the second case instead, the source would lie at larger CO luminosities than those expected for local SFGs and high-z submillimetre galaxies. At both redshifts, we derive the star formation rate (SFR) from the ALMA continuum and the physical parameters of the galaxy, such as the stellar mass (M*), by fitting its spectral energy distribution. Exploiting the results of this work, we believe that our source is a ‘main-sequence’, dusty SFG at $z$ = 4.6 (i.e. [C ii] emitter) with $\mathrm{log(SFR/M_{\odot }\, yr^{-1})}\sim 1.4$more »and log(M*/M⊙) ∼ 9.9. As a support to this scenario our galaxy, if at this redshift, lies in a massive protocluster recently discovered at $z$ ∼ 4.57, at only ∼1 proper Mpc from its centre. This work underlines the crucial role of the ALPINE survey in making a census of this class of objects, in order to unveil their contribution to the global SFR density at the end of the Reionization epoch.« less
  2. null (Ed.)
    ABSTRACT The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While these simulations produce galaxies in overall good agreement with observations, they employ different sub-grid models for galaxies and supermassive black holes (BHs). We investigate the impact of the sub-grid models on the BH mass properties of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations, focusing on the MBH − M⋆ relation and the BH mass function. All simulations predict tight MBH − M⋆ relations, and struggle to produce BHs of $M_{\rm BH}\leqslant 10^{7.5}\, \rm M_{\odot }$ in galaxies of $M_{\star }\sim 10^{10.5}\!-\!10^{11.5}\, \rm M_{\odot }$. While the time evolution of the mean MBH − M⋆ relation is mild ($\rm \Delta M_{\rm BH}\leqslant 1\, dex$ for 0 $\leqslant z \leqslant$ 5) for all the simulations, its linearity (shape) and normalization varies from simulation to simulation. The strength of SN feedback has a large impact on the linearity and time evolution for $M_{\star }\leqslant 10^{10.5}\, \rm M_{\odot }$. We find that the low-mass end is a good discriminant of the simulation models, and highlights the need for new observational constraints. At the high-mass end, strong AGN feedback can suppress the time evolutionmore »of the relation normalization. Compared with observations of the local Universe, we find an excess of BHs with $M_{\rm BH}\geqslant 10^{9}\, \rm M_{\odot }$ in most of the simulations. The BH mass function is dominated by efficiently accreting BHs ($\log _{10}\, f_{\rm Edd}\geqslant -2$) at high redshifts, and transitions progressively from the high-mass to the low-mass end to be governed by inactive BHs. The transition time and the contribution of active BHs are different among the simulations, and can be used to evaluate models against observations.« less
  3. Abstract We report the discovery of a bright ( $g = 14.5$ mag (AB), $K = 11.9$ mag (Vega)) quasar at redshift $z=0.83$ — the optically brightest (unbeamed) quasar at $z>0.4$ . SMSS J114447.77-430859.3, at a Galactic latitude of $b=+18.1^{\circ}$ , was identified by its optical colours from the SkyMapper Southern Survey (SMSS) during a search for symbiotic binary stars. Optical and near-infrared spectroscopy reveals broad Mg ii , H $\unicode{x03B2}$ , H $\unicode{x03B1}$ , and Pa $\unicode{x03B2}$ emission lines, from which we measure a black hole mass of $\log_{10}\! (M_{\mathrm{BH}}/\mathrm{M}_{\odot}) = 9.4 \pm 0.5$ . With its high luminosity, $L_{\mathrm{bol}} = (4.7\pm1.0)\times10^{47}\,\mathrm{erg\,s}^{-1}$ or $M_{i}(z=2) = -29.74$ mag (AB), we estimate an Eddington ratio of $\approx1.4$ . As the most luminous quasar known over the last ${\sim}$ 9 Gyr of cosmic history, having a luminosity $8\times$ greater than 3C 273, the source offers a range of potential follow-up opportunities.

    Direct collapse black holes (BHs) are promising candidates for producing massive z ≳ 6 quasars, but their formation requires fine-tuned conditions. In this work, we use cosmological zoom simulations to study systematically the impact of requiring: (1) low gas angular momentum (spin), and (2) a minimum incident Lyman–Werner (LW) flux in order to form BH seeds. We probe the formation of seeds (with initial masses of $M_{\rm seed} \sim 10^4\!-\!10^6\, \mathrm{M}_{\odot }\, h^{-1})$ in haloes with a total mass >3000 × Mseed and a dense, metal-poor gas mass >5 × Mseed. Within this framework, we find that the seed-forming haloes have a prior history of star formation and metal enrichment, but they also contain pockets of dense, metal-poor gas. When seeding is further restricted to haloes with low gas spins, the number of seeds formed is suppressed by factors of ∼6 compared to the baseline model, regardless of the seed mass. Seed formation is much more strongly impacted if the dense, metal-poor gas is required to have a critical LW flux (Jcrit). Even for Jcrit values as low as 50J21, no $8\times 10^{5}~\mathrm{M}_{\odot }\, h^{-1}$ seeds are formed. While lower mass ($1.25\times 10^{4},1\times 10^{5}~\mathrm{M}_{\odot }\, h^{-1}$) seeds do form, they are strongly suppressed (by factors of ∼10–100) comparedmore »to the baseline model at gas mass resolutions of $\sim 10^4~\mathrm{M}_{\odot }\, h^{-1}$ (with even stronger suppression at higher resolutions). As a result, BH merger rates are also similarly suppressed. Since early BH growth is dominated by mergers in our models, none of the seeds are able to grow to the supermassive regime ($\gtrsim 10^6~\mathrm{M}_{\odot }\, h^{-1}$) by z = 7. Our results hint that producing the bulk of the z ≳ 6 supermassive BH population may require alternate seeding scenarios that do not depend on the LW flux, early BH growth dominated by rapid or super-Eddington accretion, or a combination of these possibilities.

    « less

    Previous studies of fueling black holes in galactic nuclei have argued (on scales ${\sim}0.01{-}1000\,$pc) accretion is dynamical with inflow rates $\dot{M}\sim \eta \, M_{\rm gas}/t_{\rm dyn}$ in terms of gas mass Mgas, dynamical time tdyn, and some η. But these models generally neglected expulsion of gas by stellar feedback, or considered extremely high densities where expulsion is inefficient. Studies of star formation, however, have shown on sub-kpc scales the expulsion efficiency fwind = Mejected/Mtotal scales with the gravitational acceleration as $(1-f_{\rm wind})/f_{\rm wind}\sim \bar{a}_{\rm grav}/\langle \dot{p}/m_{\ast }\rangle \sim \Sigma _{\rm eff}/\Sigma _{\rm crit}$ where $\bar{a}_{\rm grav}\equiv G\, M_{\rm tot}(\lt r)/r^{2}$ and $\langle \dot{p}/m_{\ast }\rangle$ is the momentum injection rate from young stars. Adopting this as the simplest correction for stellar feedback, $\eta \rightarrow \eta \, (1-f_{\rm wind})$, we show this provides a more accurate description of simulations with stellar feedback at low densities. This has immediate consequences, predicting the slope and normalization of the MBH − σ and MBH − Mbulge relation, LAGN −SFR relations, and explanations for outliers in compact Es. Most strikingly, because star formation simulations show expulsion is efficient (fwind ∼ 1) below total-mass surface density $M_{\rm tot}/\pi \, r^{2}\lt \Sigma _{\rm crit}\sim 3\times 10^{9}\, \mathrm{M}_{\odotmore »}\, {\rm kpc^{-2}}$ (where $\Sigma _{\rm crit}=\langle \dot{p}/m_{\ast }\rangle /(\pi \, G)$), BH mass is predicted to specifically trace host galaxy properties above a critical surface brightness Σcrit (B-band $\mu _{\rm B}^{\rm crit}\sim 19\, {\rm mag\, arcsec^{-2}}$). This naturally explains why BH masses preferentially reflect bulge properties or central surface densities (e.g. $\Sigma _{1\, {\rm kpc}}$), not ‘total’ galaxy properties.

    « less