An 800-million-solar-mass black hole in a significantly neutral Universe at a redshift of 7.5

Detection of stellar light from quasar host galaxies at redshifts above 6

The detection of starlight from the host galaxies of quasars during the reionization epoch (z > 6) has been elusive, even with deep Hubble Space Telescope observations1,2. The current highest redshift quasar host detected3, at z = 4.5, required the magnifying effect of a foreground lensing galaxy. Low-luminosity quasars4-6 from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP)7 mitigate the challenge of detecting their underlying, previously undetected host galaxies. Here we report rest-frame optical images and spectroscopy of two HSC-SSP quasars at z > 6 with the JWST. Using near-infrared camera imaging at 3.6 and 1.5 μm and subtracting the light from the unresolved quasars, we find that the host galaxies are massive (stellar masses of 13 × and 3.4 × 1010 M☉, respectively), compact and disc-like. Near-infrared spectroscopy at medium resolution shows stellar absorption lines in the more massive quasar, confirming the detection of the host. Velocity-broadened gas in the vicinity of these quasars enables measurements of their black hole masses (1.4 × 109 and 2.0 × 108 M☉, respectively). Their location in the black hole mass-stellar mass plane is consistent with the distribution at low redshift, suggesting that the relation between black holes and their host galaxies was already in place less than a billion years after the Big Bang.

A dusty compact object bridging galaxies and quasars at cosmic dawn

Understanding how super-massive black holes form and grow in the early Universe has become a major challenge^1,2 since it was discovered that luminous quasars existed only 700 million years after the Big Bang^3,4. Simulations indicate an evolutionary sequence of dust-reddened quasars emerging from heavily dust-obscured starbursts that then transition to unobscured luminous quasars by expelling gas and dust⁵. Although the last phase has been identified out to a redshift of 7.6 (ref. ⁶), a transitioning quasar has not been found at similar redshifts owing to their faintness at optical and near-infrared wavelengths. Here we report observations of an ultraviolet compact object, GNz7q, associated with a dust-enshrouded starburst at a redshift of 7.1899 ± 0.0005. The host galaxy is more luminous in dust emission than any other known object at this epoch, forming 1,600 solar masses of stars per year within a central radius of 480 parsec. A red point source in the far-ultraviolet is identified in deep, high-resolution imaging and slitless spectroscopy. GNz7q is extremely faint in X-rays, which indicates the emergence of a uniquely ultraviolet compact star-forming region or a Compton-thick super-Eddington black-hole accretion disk at the dusty starburst core. In the latter case, the observed properties are consistent with predictions from cosmological simulations⁷ and suggest that GNz7q is an antecedent to unobscured luminous quasars at later epochs.

Dynamics and Merger Rate of Primordial Black Holes in a Cluster

The PBH clusters can be sources of gravitational waves, and the merger rate depends on the spatial distribution of PBHs in the cluster which changes over time. It is well known that gravitational collisional systems experience the core collapse that leads to significant increase of the central density and shrinking of the core. After core collapse, the cluster expands almost self-similarly (i.e., density profile extends in size without changing its shape). These dynamic processes affect the merger rate of PBHs. In this paper, the dynamics of the PBH cluster is considered using the Fokker–Planck equation. We calculate the merger rate of PBHs on cosmic time scales and show that its time dependence has a unique signature. Namely, it grows by about an order of magnitude at the moment of core collapse which depends on the characteristics of the cluster, and then decreases according to the dependence R∝t−1.48. It was obtained for monochromatic and power-law PBH mass distributions with some fixed parameters. Obtained results can be used to test the model of the PBH clusters via observation of gravitational waves at high redshift.

Peering into the dark (ages) with low-frequency space interferometers: Using the 21-cm signal of neutral hydrogen from the infant universe to probe fundamental (Astro)physics

The Dark Ages and Cosmic Dawn are largely unexplored windows on the infant Universe (z ~ 200-10). Observations of the redshifted 21-cm line of neutral hydrogen can provide valuable new insight into fundamental physics and astrophysics during these eras that no other probe can provide, and drives the design of many future ground-based instruments such as the Square Kilometre Array (SKA) and the Hydrogen Epoch of Reionization Array (HERA). We review progress in the field of high-redshift 21-cm Cosmology, in particular focussing on what questions can be addressed by probing the Dark Ages at z > 30. We conclude that only a space- or lunar-based radio telescope, shielded from the Earth's radio-frequency interference (RFI) signals and its ionosphere, enable the 21-cm signal from the Dark Ages to be detected. We suggest a generic mission design concept, CoDEX, that will enable this in the coming decades.

Probing the nature of black holes: Deep in the mHz gravitational-wave sky

Black holes are unique among astrophysical sources: they are the simplest macroscopic objects in the Universe, and they are extraordinary in terms of their ability to convert energy into electromagnetic and gravitational radiation. Our capacity to probe their nature is limited by the sensitivity of our detectors. The LIGO/Virgo interferometers are the gravitational-wave equivalent of Galileo's telescope. The first few detections represent the beginning of a long journey of exploration. At the current pace of technological progress, it is reasonable to expect that the gravitational-wave detectors available in the 2035-2050s will be formidable tools to explore these fascinating objects in the cosmos, and space-based detectors with peak sensitivities in the mHz band represent one class of such tools. These detectors have a staggering discovery potential, and they will address fundamental open questions in physics and astronomy. Are astrophysical black holes adequately described by general relativity? Do we have empirical evidence for event horizons? Can black holes provide a glimpse into quantum gravity, or reveal a classical breakdown of Einstein's gravity? How and when did black holes form, and how do they grow? Are there new long-range interactions or fields in our Universe, potentially related to dark matter and dark energy or a more fundamental description of gravitation? Precision tests of black hole spacetimes with mHz-band gravitational-wave detectors will probe general relativity and fundamental physics in previously inaccessible regimes, and allow us to address some of these fundamental issues in our current understanding of nature.

Multi-Field versus Single-Field in the Supergravity Models of Inflation and Primordial Black Holes

We review the models unifying inflation and Primordial Black Hole (PBH) formation, which are based on the modified (Starobinsky-type) supergravity. We begin with the basic (Starobinsky) inflationary model of modified gravity and its alpha-attractor-type generalizations for PBH production, and recall how all those single-field models can be embedded into the minimal supergravity. Then, we focus on the effective two-field models arising from the modified (Starobinsky-type) supergravity and compare them to the single-field models under review. Those two-field models describe double inflation whose first stage is driven by Starobinsky’s scalaron and whose second stage is driven by another scalar belonging to the supergravity multiplet. The power spectra are numerically computed, and it is found that the ultra-slow-roll regime gives rise to the enhancement (peak) in the scalar power spectrum leading to an efficient PBH formation. The resulting PBH masses and their density fraction (as part of dark matter) are found to be in agreement with cosmological observations. The PBH-induced gravitational waves, if any, are shown to be detectable by the ground-based and space-based gravitational interferometers under construction.

The limits of cosmology: role of the Moon

The lunar surface allows a unique way forward in cosmology, to go beyond current limits. The far side provides an unexcelled radio-quiet environment for probing the dark ages via 21 cm interferometry to seek elusive clues on the nature of the infinitesimal fluctuations that seeded galaxy formation. Far-infrared telescopes in cold and dark lunar polar craters will probe back to the first months of the Big Bang and study associated spectral distortions in the CMB. Optical and IR megatelescopes will image the first star clusters in the Universe and seek biosignatures in the atmospheres of unprecedented numbers of nearby habitable zone exoplanets. The goals are compelling and a stable lunar platform will enable construction of telescopes that can access trillions of modes in the sky, providing the key to exploration of our cosmic origins. This article is part of a discussion meeting issue 'Astronomy from the Moon: the next decades'.

Cosmological Advection Flows in the Presence of Primordial Black Holes as Dark Matter and Formation of First Sources

In the inflation-based cosmology the dark matter (DM) density component starts moving with respect to the universal expansion at z_{eq}∼3200 while baryons remain frozen until z_{rec}∼1100. It has been suggested that in this case postlinear corrections to the evolution of small fluctuations would result, for the standard Λ-dominated cold DM (CDM) model, in delayed formation of early objects as supersonic advection flows develop after recombination, so baryons are not immediately captured by the DM gravity on small scales. We develop the hydrodynamical description of such two-component advection and show that, in the supersonic regime, the advection within irrotational fluids is governed by the gradient of the difference of the kinetic energies of the two (DM and baryonic here) components. We then apply this formalism to the case where DM is made up of LIGO-type black holes (BHs) and show that there the advection process on scales relevant for early structure collapse will differ significantly from the earlier discussed (CDM) case because of the additional granulation component to the density field produced during inflation. The advection here will lead efficiently to the common motion of the DM and baryon components on scales relevant for collapse and formation of first luminous sources. This leads to early collapse, making it easier to explain the existence of supermassive BHs observed in quasars at high z>7. The resultant net advection rate reaches minimum around ≲10^{9}  M_{⊙} and subsequently rises to a secondary maximum near the typical mass of ∼10^{12}  M_{⊙}, which may be an important consideration for formation of galaxies at z≲(a few).

Leptonic and Hadronic Radiative Processes in Supermassive-Black-Hole Jets

Supermassive black holes lying in the center of galaxies can launch relativistic jets of plasma along their polar axis. The physics of black-hole jets is a very active research topic in astrophysics, owing to the fact that many questions remain open on the physical mechanisms of jet launching, of particle acceleration in the jet, and on the radiative processes. In this work I focus on the last item, and present a review of the current understanding of radiative emission processes in supermassive-black-hole jets.

Fast jet proper motion discovered in a blazar at z=4.72

High-resolution observations of high-redshift (z>4) radio quasars offer a unique insight into jet kinematics at early cosmological epochs, as well as constraints on cosmological model parameters. Due to the general weakness of extremely distant objects and the apparently slow structural changes caused by cosmological time dilation, only a couple of high-redshift quasars (HRQs) have been studied with parsec-scale resolutions, and with limited number of observing epochs. Here we report on very long baseline interferometry (VLBI) observations of a high-redshift blazar J1430 + 4204 (z=4.72) in the 8 GHz frequency band at five different epochs spanning 22 years. The source shows a compact core-jet structure with two jet components being identified within 3 milli-arcsecond (mas) scale. The long time span and multiple-epoch data allow for the kinematic studies of the jet components. That results in a jet proper motion of μ(J1) = 0.017 ± 0.002 mas a^-1 and μ(J2) = 0.156 ± 0.015 mas a^-1, respectively. For the fastest-moving outer jet component J2, the corresponding apparent transverse speed is (19.5±1.9)c. The inferred bulk jet Lorentz factor Γ=14.6±3.8 and viewing angle θ=2.2°±1.6° indicate highly relativistic beaming. The Lorentz factor and apparent proper motion are the highest measured to date among the z>4 jetted radio sources, while the jet kinematics is still consistent with the cosmological interpretation of quasar redshifts.

Beyond Optical Depth: Future Determination of Ionization History from the Cosmic Microwave Background

We explore the fundamental limits to which reionization histories can be constrained using only large-scale cosmic microwave background (CMB) anisotropy measurements. The redshift distribution of the fractional ionization x _e (z) affects the angular distribution of CMB polarization. We project constraints on the reionization history of the universe using low-noise full-sky temperature and E-mode measurements of the CMB. We show that the measured TE power spectrum, C ^ ℓ TE , has roughly one quarter of the constraining power of C ^ ℓ EE on the reionization optical depth τ, and its addition improves the precision on τ by 20% over using C ^ ℓ EE only. We also use a two-step reionization model with an additional high-redshift step, parameterized by an early ionization fraction x e min , and a late reionization step at z _re. We find that future high signal-to-noise measurements of the multipoles 10 ⩽ ℓ < 20 are especially important for breaking the degeneracy between x e min and z _re. In addition, we show that the uncertainties on these parameters determined from a map with sensitivity 10 μK arcmin are less than 5% larger than the uncertainties in the noiseless case, making this noise level a natural target for future large sky area E-mode measurements.

Primordial Black Holes from the QCD Axion

We propose a mechanism to generate primordial black holes (PBHs) that is independent of cosmological inflation and occurs slightly below the QCD phase transition. Our setup relies on the collapse of long-lived string-domain wall networks and is naturally realized in QCD axion models with domain wall number N_{DW}>1 and Peccei-Quinn symmetry broken after inflation. In our framework, dark matter is mostly composed of axions in the meV mass range along with a small fraction, Ω_{PBH}≳10^{-6}Ω_{CDM} of heavy M∼10^{4}-10^{7}  M_{⊙} PBHs. The latter could play a role in alleviating some of the shortcomings of the standard cosmological model on subgalactic scales. The scenario might have distinct signatures in ongoing axion searches as well as gravitational wave observatories.

Formation of massive black holes in rapidly growing pre-galactic gas clouds

The origin of the supermassive black holes that inhabit the centres of massive galaxies remains unclear^1,2. Direct-collapse black holes-remnants of supermassive stars, with masses around 10,000 times that of the Sun-are ideal seed candidates^3-6. However, their very existence and their formation environment in the early Universe are still under debate, and their supposed rarity makes modelling their formation difficult^7,8. Models have shown that rapid collapse of pre-galactic gas (with a mass infall rate above some critical value) in metal-free haloes is a requirement for the formation of a protostellar core that will then form a supermassive star^9,10. Here we report a radiation hydrodynamics simulation of early galaxy formation^11,12 that produces metal-free haloes massive enough and with sufficiently high mass infall rates to form supermassive stars. We find that pre-galactic haloes and their associated gas clouds that are exposed to a Lyman-Werner intensity roughly three times the intensity of the background radiation and that undergo at least one period of rapid mass growth early in their evolution are ideal environments for the formation of supermassive stars. The rapid growth induces substantial dynamical heating^13,14, amplifying the Lyman-Werner suppression that originates from a group of young galaxies 20 kiloparsecs away. Our results strongly indicate that the dynamics of structure formation, rather than a critical Lyman-Werner flux, is the main driver of the formation of massive black holes in the early Universe. We find that the seeds of massive black holes may be much more common than previously considered in overdense regions of the early Universe, with a co-moving number density up to 10^-3 per cubic megaparsec.

Formation of the first generation of stars and blackholes in the Universe

Modern sky surveys using large ground-based telescopes have discovered a variety of celestial objects. Prominent structures such as galaxies and galaxy clusters are found virtually everywhere, and their collective distribution forms the large-scale structure of the Universe. It is thought that all of the rich content in the present-day Universe developed through gravitational amplification of primeval density fluctuations generated in the very early phase of cosmic evolution. The standard theoretical model based on an array of recent observations accurately predicts the physical conditions in the early Universe, and powerful super-computers allow us to simulate in detail the formation and evolution of cosmic structure to the present epoch. We review recent progress in the study on the first generation of stars and blackholes. We focus on the physics of early structure formation, while identifying several key issues and open questions. Finally, we discuss prospects for future observations of the first stars, galaxies and blackholes.

The route to massive black hole formation via merger-driven direct collapse: a review

The direct collapse model for the formation of massive black holes has gained increased support as it provides a natural explanation for the appearance of bright quasars already less than a billion years from the Big Bang. In this paper we review a recent scenario for direct collapse that relies on multi-scale gas inflows initiated by the major merger of massive gas-rich galaxies at z  >  6, where gas has already achieved solar composition. Hydrodynamical simulations undertaken to explore our scenario show that supermassive, gravitationally bound compact gaseous disks weighing a billion solar masses, only a few pc in size, form in the nuclei of merger remnants in less than 10⁵ yr. These could later produce a supermassive protostar or supermassive star at their center via various mechanisms. Moreover, we present a new analytical model, based on angular momentum transport in mass-loaded gravitoturbulent disks. This naturally predicts that a nuclear disk accreting at rates exceeding [Formula: see text] yr^-1, as seen in the simulations, is stable against fragmentation irrespective of its metallicity. This is at variance with conventional direct collapse scenarios, which require the suppression of gas cooling in metal-free protogalaxies for gas collapse to take place. Such high accretion rates reflect the high free-fall velocities in massive halos appearing only at z  <  10, and occur naturally as a result of the efficient angular momentum loss provided by the merger dynamics. We discuss the implications of our scenario on the observed population of high-z quasars and on its evolution to lower redshifts using a semi-analytical galaxy formation model. Finally, we consider the intriguing possibility that the secondary gas inflows in the unstable disks might drive gas to collapse into a supermassive black hole directly via the General Relativistic radial instability. Such dark collapse route could generate gravitational wave emission detectable via the future Laser Interferometer Space Antenna (LISA).

Simulating the magnetorotational collapse of supermassive stars: Incorporating gas pressure perturbations and different rotation profiles

Collapsing supermassive stars (SMSs) with masses M ≳ 104-6 M ⊙ have long been speculated to be the seeds that can grow and become supermassive black holes (SMBHs). We previously performed general relativistic magnetohydrodynamic (GRMHD) simulations of marginally stable Γ = 4/3 polytropes uniformly rotating at the mass-shedding limit and endowed initially with a dynamically unimportant dipole magnetic field to model the direct collapse of SMSs. These configurations are supported entirely by thermal radiation pressure and reliably model SMSs with M ≳ 106 M ⊙. We found that around 90% of the initial stellar mass forms a spinning black hole (BH) remnant surrounded by a massive, hot, magnetized torus, which eventually launches a magnetically-driven jet. SMSs could be therefore sources of ultra-long gamma-ray bursts (ULGRBs). Here we perform GRMHD simulations of Γ ≳ 4/3, polytropes to account for the perturbative role of gas pressure in SMSs with M ≲ 106 M ⊙. We also consider different initial stellar rotation profiles. The stars are initially seeded with a dynamically weak dipole magnetic field that is either confined to the stellar interior or extended from its interior into the stellar exterior. We calculate the gravitational wave burst signal for the different cases. We find that the mass of the black hole remnant is 90%-99% of the initial stellar mass, depending sharply on Γ - 4/3 as well as on the initial stellar rotation profile. After t ~ 250-550M ≈ 1 - 2 × 103(M/106 M ⊙) s following the appearance of the BH horizon, an incipient jet is launched and it lasts for ~104-105(M/106 M ⊙) s, consistent with the duration of long gamma-ray bursts. Our numerical results suggest that the Blandford-Znajek mechanism powers the incipient jet. They are also in rough agreement with our recently proposed universal model that estimates accretion rates and electromagnetic (Poynting) luminosities that characterize magnetized BH-disk remnant systems that launch a jet. This model helps explain why the outgoing electromagnetic luminosities computed for vastly different BH-disk formation scenarios all reside within a narrow range (~1052±1 erg s-1), roughly independent of M.

Maximally rotating supermassive stars at the onset of collapse: the perturbative effects of gas pressure, magnetic fields, dark matter, and dark energy

The discovery of quasars at increasingly large cosmological redshifts may favour 'direct collapse' as the most promising evolutionary route to the formation of supermassive black holes. In this scenario, supermassive black holes form when their progenitors - supermassive stars - become unstable to gravitational collapse. For uniformly rotating stars supported by pure radiation pressure and spinning at the mass-shedding limit, the critical configuration at the onset of collapse is characterized by universal values of the dimensionless spin and radius parameters J/M2 and R/M, independent of mass M. We consider perturbative effects of gas pressure, magnetic fields, dark matter, and dark energy on these parameters, and thereby determine the domain of validity of this universality. We obtain leading-order corrections for the critical parameters and establish their scaling with the relevant physical parameters. We compare two different approaches to approximate the effects of gas pressure, which plays the most important role, find identical results for the above dimensionless parameters, and also find good agreement with recent numerical results.