Reionization with galaxies and active galactic nuclei

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_{ m esc}^{ m sf} angle$) and AGN ($f_{ m esc}^{ m bh}$) to find: (i) the ionizing budget is dominated by stellar radiation from low stellar mass ($M_∗lt 10^9 , { m m M_odot }$) galaxies at z > 6 with the AGN contribution (driven by $M_{bh}gt 10^6 , { m m M_odot }$ black holes in $M_∗ gtrsim 10^9, { m m M_odot }$ galaxies) dominating at lower redshifts; (ii) AGN only contribute $10-25{{ m 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_{ m esc}^{ m sf} angle$ is shallower than $f_{ m esc}^{ m 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

Ultra-light Axions: Degeneracies with Massive Neutrinos and Forecasts for Future Cosmological Observations

A generic prediction of string theory is the existence of many axion fields. It has recently been argued that many of these fields should be light and, like the well known QCD axion, lead to observable cosmological consequences. In this paper we study in detail the effect of the so-called string axiverse on large scale structure, focusing on the morphology and evolution of density perturbations, anisotropies in the cosmic microwave background and weak gravitational lensing of distant galaxies. We quantify specific effects that will arise from the presence of the axionic fields and highlight possible degeneracies that may arise in the presence of massive neutrinos. We take particular care understanding the different physical effects and scales that come into play. We then forecast how the string axiverse may be constrained and show that with a combination of different observations, it should be possible to detect a fraction of ultralight axions to dark matter of a few percent.Comment: 24 pages, 16 figures, this version: corrected typos, some comments added, matches published versio

Multimessenger study of merging massive black holes in the Obelisk simulation: gravitational waves, electromagnetic counterparts, and their link to galaxy and black hole populations

Massive black hole (BH) mergers are predicted to be powerful sources of low-frequency gravitational waves (GWs). Coupling the detection of GWs with electromagnetic (EM) detection can provide key information about merging BHs and their environments. We study the high-resolution cosmological radiation-hydrodynamics simulation Obelisk, run to redshift $z=3.5$, to assess the GW and EM detectability of high-redshift BH mergers, modelling spectral energy distribution and obscuration. For EM detectability we further consider sub-grid dynamical delays in postprocessing. We find that most of the merger events can be detected by LISA, except for high-mass mergers with very unequal mass ratios. Intrinsic binary parameters are accurately measured, but the sky localisation is poor generally. Only $\sim 40\%$ of these high-redshift sources have sky localisation better than $10\,\mathrm{deg}^2$. Merging BHs are hard to detect in the restframe UV since they are fainter than the host galaxies, which at high $z$ are star-forming. A significant fraction, $15-35\%$, of BH mergers instead outshines the galaxy in X-rays, and about $5-15\%$ are sufficiently bright to be detected with sensitive X-ray instruments. If mergers induce an Eddington-limited brightening, up to $30\%$ of sources can become observable. The transient flux change originating from such a brightening is often large, allowing $4-20\%$ of mergers to be detected as EM counterparts. A fraction $1-30\%$ of mergers is also detectable at radio frequencies. Observable merging BHs tend to have higher accretion rates and masses and are overmassive at fixed galaxy mass with respect to the full population. Most EM-observable mergers can also be GW-detected with LISA, but their sky localisation is generally poorer. This has to be considered when using EM counterparts to obtain information about the properties of merging BHs and their environment.Comment: 17 pages, 11 figures, submitted to A&

Population statistics of intermediate mass black holes in dwarf galaxies using the NewHorizon simulation

While it is well established that supermassive black holes (SMBHs) co-evolve with their host galaxy, it is currently less clear how lower mass black holes, so-called intermediate mass black holes (IMBHs), evolve within their dwarf galaxy hosts. In this paper, we present results on the evolution of a large sample of IMBHs from the NewHorizon simulation. We show that occupation fractions of IMBHs in dwarf galaxies are at least 50 percent for galaxies with stellar masses down to 1E6 Msun, but BH growth is very limited in dwarf galaxies. In NewHorizon, IMBH growth is somewhat more efficient at high redshift z = 3 but in general IMBH do not grow significantly until their host galaxy leaves the dwarf regime. As a result, NewHorizon under-predicts observed AGN luminosity function and AGN fractions. We show that the difficulties of IMBH to remain attached to the centres of their host galaxies plays an important role in limiting their mass growth, and that this dynamic evolution away from galactic centres becomes stronger at lower redshift.Comment: 15 pages, submitted to MNRA

Gas Accretion and Giant Lyman-alpha Nebulae

Several decades of observations and discoveries have shown that high-redshift AGN and massive galaxies are often surrounded by giant Lyman-alpha nebulae extending in some cases up to 500 kpc in size. In this review, I discuss the properties of the such nebulae discovered at z>2 and their connection with gas flows in and around the galaxies and their halos. In particular, I show how current observations are used to constrain the physical properties and origin of the emitting gas in terms of the Lyman-alpha photon production processes and kinematical signatures. These studies suggest that recombination radiation is the most viable scenario to explain the observed Lyman-alpha luminosities and Surface Brightness for the large majority of the nebulae and imply that a significant amount of dense, ionized and cold clumps should be present within and around the halos of massive galaxies. Spectroscopic studies suggest that, among the giant Lyman-alpha nebulae, the one associated with radio-loud AGN should have kinematics dominated by strong, ionized outflows within at least the inner 30-50 kpc. Radio-quiet nebulae instead present more quiescent kinematics compatible with stationary situation and, in some cases, suggestive of rotating structures. However, definitive evidences for accretion onto galaxies of the gas associated with the giant Lyman-alpha emission are not unambiguously detected yet. Deep surveys currently ongoing using other bright, non-resonant lines such as Hydrogen H-alpha and HeII1640 will be crucial to search for clearer signatures of cosmological gas accretion onto galaxies and AGN.Comment: Invited review to appear in Gas Accretion onto Galaxies, Astrophysics and Space Science Library, eds. A. J. Fox & R. Dave', to be published by Springe

A quasar-galaxy merger at z ~ 6.2: Black hole mass and quasar properties from the NIRSpec spectrum

We present JWST/NIRSpec integral field data of the quasar PJ308-21 at z = 6.2342. As shown by previous ALMA and HST imaging, the quasar has two companion sources, interacting with the quasar host galaxy. The high-resolution G395H/290LP NIRSpec spectrum covers the 2.87 − 5.27 μm wavelength range and shows the rest-frame optical emission of the quasar with exquisite quality (signal-to-noise ratio ∼100 − 400 per spectral element). Based on the Hβ line from the broad line region, we obtain an estimate of the black hole mass MBH, Hβ ∼ 2.7 × 109 M⊙. This value is within a factor ≲1.5 of the Hα-based black hole mass from the same spectrum (MBH, Hα ∼ 1.93 × 109 M⊙) and is consistent with a previous estimate relying on the Mg IIλ2799 line (MBH, MgII ∼ 2.65 × 109 M⊙). All these MBH estimates are within the ∼0.5 dex intrinsic scatter of the adopted mass calibrations. The high Eddington ratio of PJ308-21 λEdd, Hβ ∼ 0.67 (λEdd, Hα ∼ 0.96) is in line with the overall quasar population at z ≳ 6. The relative strengths of the [O III], Fe II, and Hβ lines are consistent with the empirical “Eigenvector 1” correlations as observed for low redshift quasars. We find evidence for blueshifted [O III] λ5007 emission with a velocity offset Δv[O III] = −1922 ± 39 km s−1 from the systemic velocity and a full width at half maximum (FWHM) FWHM([O III]) = 2776−74+75 km s−1. This may be the signature of outflowing gas from the nuclear region, despite the true values of Δv[O III] and FWHM([O III]) likely being more uncertain due to the blending with Hβ and Fe II lines. Our study demonstrates the unique capabilities of NIRSpec in capturing quasar spectra at cosmic dawn and studying their properties in unprecedented detail

The Low-redshift Lyman Continuum Survey: Radio continuum properties of low-zz Lyman continuum emitters

The sources that leak Lyman-continuum (LyC) photons and lead to the reionisation of the universe are intensely studied using multiple observing facilities. Recently, the Low-redshift LyC Survey (LzLCS) has found the first large sample of LyC emitting galaxies at low redshift ($z\sim 0.3$) with the Hubble Space Telescope/Cosmic Origins Spectrograph. The LzLCS sample contains a robust estimate of the LyC escape fraction ($f_\mathrm{esc}^\mathrm{LyC}$) for 66 galaxies spanning a wide range of $f_\mathrm{esc}^\mathrm{LyC}$. Here we, for the first time, aim to study the radio continuum (RC) properties of LzLCS sources and their dependence on $f_\mathrm{esc}^\mathrm{LyC}$. We present Karl G. Jansky Very Large Array RC observations at C (4-8 GHz), S (2-4 GHz) and L (1-2 GHz) bands for a sub-sample of the LzLCS sources. The radio spectral index ($\alpha^{\mathrm{3GHz}}_\mathrm{6GHz}$) spans a wide range from being flat ( $\geq -0.1$) to very steep ($\leq -1.0$). We find that the strongest leakers in our sample show flat $\alpha^{\mathrm{3GHz}}_\mathrm{6GHz}$, weak leakers have $\alpha^{\mathrm{3GHz}}_\mathrm{6GHz}$ close to normal star-forming galaxies, and non-leakers are characterized by steep $\alpha^{\mathrm{3GHz}}_\mathrm{6GHz}$. We argue that a combination of young ages, free-free absorption, and a flat cosmic-ray energy spectrum can altogether lead to a flat $\alpha^{\mathrm{3GHz}}_\mathrm{6GHz}$ for strong leakers. Non-leakers are characterized by steep spectra which can arise due to break/cutoff at high frequencies. Such a cutoff in the spectrum can arise in a single injection model of CRs characteristic of galaxies which have recently stopped star formation. Such a relation between $\alpha^{\mathrm{3GHz}}_\mathrm{6GHz}$ and $f_\mathrm{esc}^\mathrm{LyC}$ hints at the interesting role of supernovae, CRs, and magnetic fields in facilitating the escape (and/or the lack) of LyC photons.Comment: 25 pages, 14 figures, 3 tables, Submitted to Astronomy & Astrophysic

The Low-Redshift Lyman Continuum Survey. Unveiling the ISM properties of low-zz Lyman continuum emitters

Combining 66 ultraviolet (UV) spectra and ancillary data from the Low-Redshift Lyman Continuum Survey (LzLCS) and 23 LyC observations by earlier studies, we form a statistical sample of star-forming galaxies at $z \sim 0.3$ to study the role of the cold interstellar medium (ISM) gas in the leakage of ionizing radiation. We first constrain the massive star content (ages and metallicities) and UV attenuation, by fitting the stellar continuum with a combination of simple stellar population models. The models, together with accurate LyC flux measurements, allow to determine the absolute LyC photon escape fraction for each galaxy ($f_{\rm esc}^{\rm abs}$). We measure the equivalent widths and residual fluxes of multiple HI and low-ionization state (LIS) lines, and the geometrical covering fraction adopting the picket-fence model. The $f_{\rm esc}^{\rm abs}$ spans a wide range, with a median (0.16, 0.84 quantiles) of 0.04 (0.02, 0.20), and 50 out of the 89 galaxies detected in the LyC. The HI and LIS line equivalent widths scale with the UV luminosity and attenuation, and inversely with the residual flux of the lines. The HI and LIS residual fluxes are correlated, indicating that the neutral gas is spatially traced by the LIS transitions. We find the observed trends of the absorption lines and the UV attenuation are primarily driven by the covering fraction. The non-uniform gas coverage demonstrates that LyC photons escape through low-column density channels in the ISM. The equivalent widths and residual fluxes of the UV lines strongly correlate with $f_{\rm esc}^{\rm abs}$: strong LyC leakers show weak absorption lines, low UV attenuation, and large Ly$\alpha$ equivalent widths. We finally show that simultaneous UV absorption line and dust attenuation measurements can predict, on average, the escape fraction of galaxies and the method can be applied to galaxies across a wide redshift range.Comment: 30 pages, 16 figures, 3 tables; accepted for publication in Astronomy and Astrophysics on December 16, 2021. Tables A1 to A4 are part of the LzLCS science products and will be publicly available in a dedicated websit

The COSMOS-Web ring: in-depth characterization of an Einstein ring lensing system at z~2

Aims. We provide an in-depth analysis of the COSMOS-Web ring, an Einstein ring at z=2 that we serendipitously discovered in the COSMOS-Web survey and possibly the most distant lens discovered to date. Methods. We extract the visible and NIR photometry from more than 25 bands and we derive the photometric redshifts and physical properties of both the lens and the source with three different SED fitting codes. Using JWST/NIRCam images, we also produce two lens models to (i) recover the total mass of the lens, (ii) derive the magnification of the system, (iii) reconstruct the morphology of the lensed source, and (iv) measure the slope of the total mass density profile of the lens. Results. The lens is a very massive and quiescent (sSFR < 10^(-13) yr-1) elliptical galaxy at z = 2.02 \pm 0.02 with a total mass Mtot(<thetaE) = (3.66 \pm 0.36) x 10^11 Msun and a stellar mass M* = (1.37 \pm 0.14) x 10^11 Msun. Compared to SHMRs from the literature, we find that the total mass is consistent with the presence of a DM halo of mass Mh = 1.09^(+1.46)_(-0.57) x 10^13 Msun. In addition, the background source is a M* = (1.26 \pm 0.17) x 10^10 Msun star-forming galaxy (SFR=(78 \pm 15) Msun/yr) at z = 5.48 \pm 0.06. Its reconstructed morphology shows two components with different colors. Dust attenuation values from SED fitting and nearby detections in the FIR also suggest it could be partially dust-obscured. Conclusions. We find the lens at z=2. Its total, stellar, and DM halo masses are consistent within the Einstein ring, so we do not need any unexpected changes in our description of the lens (e.g. change its IMF or include a non-negligible gas contribution). The most likely solution for the lensed source is at z = 5.5. Its reconstructed morphology is complex and highly wavelength dependent, possibly because it is a merger or a main sequence galaxy with a heterogeneous dust distribution.Comment: 16 pages, submitted to A&