-
A unified model for the clustering of quasars and galaxies at $z\approx6$
Authors:
Elia Pizzati,
Joseph F. Hennawi,
Joop Schaye,
Matthieu Schaller,
Anna-Christina Eilers,
Feige Wang,
Carlos S. Frenk,
Willem Elbers,
John C. Helly,
Ruari Mackenzie,
Jorryt Matthee,
Rongmon Bordoloi,
Daichi Kashino,
Rohan P. Naidu,
Minghao Yue
Abstract:
Recent observations from the EIGER JWST program have measured for the first time the quasar-galaxy cross-correlation function at $z\approx6$. The auto-correlation function of faint $z\approx6$ quasars was also recently estimated. These measurements provide key insights into the properties of quasars and galaxies at high redshift and their relation with the host dark matter halos. In this work, we…
▽ More
Recent observations from the EIGER JWST program have measured for the first time the quasar-galaxy cross-correlation function at $z\approx6$. The auto-correlation function of faint $z\approx6$ quasars was also recently estimated. These measurements provide key insights into the properties of quasars and galaxies at high redshift and their relation with the host dark matter halos. In this work, we interpret these data building upon an empirical quasar population model that has been applied successfully to quasar clustering and demographic measurements at $z\approx2-4$. We make use of a new, large-volume N-body simulation with more than a trillion particles, FLAMINGO-10k, to model quasars and galaxies simultaneously. We successfully reproduce observations of $z\approx6$ quasars and galaxies (i.e., their clustering properties and luminosity functions), and infer key quantities such as their luminosity-halo mass relation, the mass function of their host halos, and their duty cycle/occupation fraction. Our key findings are: (i) quasars reside on average in $\approx10^{12.4}\,{\rm M}_\odot$ halos (corresponding to $\approx5σ$ fluctuations in the initial conditions of the linear density field), but the distribution of host halo masses is quite broad; (ii) the duty cycle of (UV-bright) quasar activity is relatively low ($\lesssim1\%$); (iii) galaxies (that are bright in [OIII]) live in much smaller halos ($\approx10^{10.9}\,{\rm M}_\odot$) and have a larger duty cycle (occupation fraction) of $\approx20\%$. Finally, we focus on the inferred properties of quasars and present a homogeneous analysis of their evolution with redshift. The picture that emerges reveals a strong evolution of the host halo mass and duty cycle of quasars at $z\approx2-6$, and calls for new investigations of the role of quasar activity across cosmic time.
△ Less
Submitted 18 March, 2024;
originally announced March 2024.
-
EIGER VI. The Correlation Function, Host Halo Mass and Duty Cycle of Luminous Quasars at $z\gtrsim6$
Authors:
Anna-Christina Eilers,
Ruari Mackenzie,
Elia Pizzati,
Jorryt Matthee,
Joseph F. Hennawi,
Haowen Zhang,
Rongmon Bordoloi,
Daichi Kashino,
Simon J. Lilly,
Rohan P. Naidu,
Robert A. Simcoe,
Minghao Yue,
Carlos S. Frenk,
John C. Helly,
Matthieu Schaller,
Joop Schaye
Abstract:
We expect luminous ($M_{1450}\lesssim-26.5$) high-redshift quasars to trace the highest density peaks in our universe, and therefore to reside in proto-clusters encompassing an abundance of galaxies in close vicinity. Here, we present observations of four $z\gtrsim6$ quasar fields using JWST/NIRCam in imaging and widefield slitless spectroscopy mode and report a wide range in the number of detecte…
▽ More
We expect luminous ($M_{1450}\lesssim-26.5$) high-redshift quasars to trace the highest density peaks in our universe, and therefore to reside in proto-clusters encompassing an abundance of galaxies in close vicinity. Here, we present observations of four $z\gtrsim6$ quasar fields using JWST/NIRCam in imaging and widefield slitless spectroscopy mode and report a wide range in the number of detected [OIII]-emitting galaxies in the quasars' environments, ranging between a density enhancement of $δ>100$ within a $2$ cMpc radius - one of the largest proto-clusters during the Epoch of Reionization discovered to date - to a density contrast consistent with zero, indicating the presence of a UV-luminous quasar in a region comparable to the average density of the universe. By measuring the two-point cross-correlation function of quasars and their surrounding galaxies, as well as the galaxy auto-correlation function, we infer a correlation length of quasars at $\langle z\rangle=6.25$ of $r_0^{\rm QQ}=21.3^{+2.7}_{-2.6}~{\rm cMpc}\,h^{-1}$, while we obtain a correlation length of the [OIII]-emitting galaxies of $r_0^{\rm GG}=4.2\pm0.1~{\rm cMpc}\,h^{-1}$. By comparing the correlation functions to dark-matter-only simulations we estimate the minimum mass of the quasars' host dark matter halos to be $\log_{10}(M_{\rm halo, min}/M_\odot)=12.30\pm0.14$ (and $\log_{10}(M_{\rm halo, min}^{\rm [OIII]}/M_\odot) = 10.72\pm0.03$ for the [OIII]-emitters), indicating that (a) luminous quasars do not necessarily reside within the most overdense regions in the early universe, and that (b) the UV-luminous duty cycle of quasar activity at these redshifts is $f_{\rm duty}\ll1$. Such short quasar activity timescales challenge our understanding of early supermassive black hole growth and provide evidence for highly dust-obscured growth phases or episodic, radiatively inefficient accretion rates.
△ Less
Submitted 12 March, 2024;
originally announced March 2024.
-
Revisiting the Extreme Clustering of $z \approx 4$ Quasars with Large Volume Cosmological Simulations
Authors:
Elia Pizzati,
Joseph F. Hennawi,
Joop Schaye,
Matthieu Schaller
Abstract:
Observations from wide-field quasar surveys indicate that the quasar auto-correlation length increases dramatically from $z\approx2.5$ to $z\approx4$. This large clustering amplitude at $z\approx4$ has proven hard to interpret theoretically, as it implies that quasars are hosted by the most massive dark matter halos residing in the most extreme environments at that redshift. In this work, we prese…
▽ More
Observations from wide-field quasar surveys indicate that the quasar auto-correlation length increases dramatically from $z\approx2.5$ to $z\approx4$. This large clustering amplitude at $z\approx4$ has proven hard to interpret theoretically, as it implies that quasars are hosted by the most massive dark matter halos residing in the most extreme environments at that redshift. In this work, we present a model that simultaneously reproduces both the observed quasar auto-correlation and quasar luminosity functions. The spatial distribution of halos and their relative abundance are obtained via a novel method that computes the halo mass and halo cross-correlation functions by combining multiple large-volume dark-matter-only cosmological simulations with different box sizes and resolutions. Armed with these halo properties, our model exploits the conditional luminosity function framework to describe the stochastic relationship between quasar luminosity, $L$, and halo mass, $M$. Assuming a simple power-law relation $L\propto M^γ$ with log-normal scatter, $σ$, we are able to reproduce observations at $z\sim 4$ and find that: (a) the quasar luminosity-halo mass relation is highly non-linear ($γ\gtrsim2$), with very little scatter ($σ\lesssim 0.3$ dex); (b) luminous quasars ($\log_{10} L/{\rm erg s}^{-1} \gtrsim 46.5-47$) are hosted by halos with mass $\log_{10} M/{\rm M}_\odot\gtrsim 13-13.5$; and (c) the implied duty cycle for quasar activity approaches unity ($\varepsilon_{\rm DC}\approx10-60\%$). We also consider observations at $z\approx2.5$ and find that the quasar luminosity-halo mass relation evolves significantly with cosmic time, implying a rapid change in quasar host halo masses and duty cycles, which in turn suggests concurrent evolution in black hole scaling relations and/or accretion efficiency.
△ Less
Submitted 1 February, 2024; v1 submitted 28 November, 2023;
originally announced November 2023.
-
Constraining turbulence in protoplanetary discs using the gap contrast: an application to the DSHARP sample
Authors:
E. Pizzati,
G. P. Rosotti,
B. Tabone
Abstract:
Constraining the strength of gas turbulence in protoplanetary discs is an open problem that has relevant implications for the physics of gas accretion and planet formation. In this work, we gauge the amount of turbulence in 6 of the discs observed in the DSHARP programme by indirectly measuring the vertical distribution of their dust component. We employ the differences in the gap contrasts observ…
▽ More
Constraining the strength of gas turbulence in protoplanetary discs is an open problem that has relevant implications for the physics of gas accretion and planet formation. In this work, we gauge the amount of turbulence in 6 of the discs observed in the DSHARP programme by indirectly measuring the vertical distribution of their dust component. We employ the differences in the gap contrasts observed along the major and the minor axes due to projection effects, and build a radiative transfer model to reproduce these features for different values of the dust scale heights. We find that (a) the scale heights that yield a better agreement with data are generally low ($\lesssim 4$ AU at a radial distance of $100$ AU), and in almost all cases we are only able to place upper limits on their exact values; these conclusions imply (assuming an average Stokes number of $\approx10^{-2}$) low turbulence levels of $α_{\rm SS}\lesssim10^{-3}-10^{-4}$; (b) for the 9 other systems we considered out of the DSHARP sample, our method yields no significant constraints on the disc vertical structure; we conclude that this is because these discs have either a low inclination or gaps that are not deep enough. Based on our analysis we provide an empirical criterion to assess whether a given disc is suitable to measure the vertical scale height.
△ Less
Submitted 20 July, 2023;
originally announced July 2023.
-
[CII] halos in ALPINE galaxies: smoking-gun of galactic outflows?
Authors:
E. Pizzati,
A. Ferrara,
A. Pallottini,
L. Sommovigo,
M. Kohandel,
S. Carniani
Abstract:
ALMA observations have revealed that many high redshift galaxies are surrounded by extended (10-15 kpc) [CII]-emitting halos which are not predicted by even the most advanced zoom-in simulations. Using a semi-analytical model, in a previous work we suggested that such halos are produced by starburst-driven, catastrophically cooling outflows. Here, we further improve the model and compare its predi…
▽ More
ALMA observations have revealed that many high redshift galaxies are surrounded by extended (10-15 kpc) [CII]-emitting halos which are not predicted by even the most advanced zoom-in simulations. Using a semi-analytical model, in a previous work we suggested that such halos are produced by starburst-driven, catastrophically cooling outflows. Here, we further improve the model and compare its predictions with data from 7 star-forming ($10\lesssim \rm SFR/ M_\odot \rm yr^{-1}<100$) galaxies at z=4-6, observed in the ALPINE survey. We find that (a) detected [CII] halos are a natural by-product of starburst-driven outflows; (b) the outflow mass loading factors are in the range $4\lesssimη\lesssim 7$, with higher $η$ values for lower-mass, lower-SFR systems, and scale with stellar mass as $η\propto M_*^{-0.43}$, consistently with the momentum-driven hypothesis. Our model suggests that outflows are widespread phenomena in high-z galaxies. However, in low-mass systems the halo extended [CII] emission is likely too faint to be detected with the current levels of sensitivity.
△ Less
Submitted 15 December, 2022;
originally announced December 2022.
-
A new look at the infrared properties of z $\sim$ 5 galaxies
Authors:
L. Sommovigo,
A. Ferrara,
S. Carniani,
A. Pallottini,
P. Dayal,
E. Pizzati,
M. Ginolfi,
V. Markov,
A. Faisst
Abstract:
Recent ALMA large surveys unveiled the presence of significant dust continuum emission in star-forming galaxies at $z>4$. Unfortunately, such large programs -- i.e. ALPINE ($z\sim 5$) and REBELS ($z \sim 7$) -- only provide us with a single Far-Infrared (FIR) continuum data point for their individual targets. Therefore, high-$z$ galaxies FIR spectral energy densities (SEDs) remain mostly unconstra…
▽ More
Recent ALMA large surveys unveiled the presence of significant dust continuum emission in star-forming galaxies at $z>4$. Unfortunately, such large programs -- i.e. ALPINE ($z\sim 5$) and REBELS ($z \sim 7$) -- only provide us with a single Far-Infrared (FIR) continuum data point for their individual targets. Therefore, high-$z$ galaxies FIR spectral energy densities (SEDs) remain mostly unconstrained, hinging on an assumption for their dust temperature ($T_{\rm d}$) in the SED fitting procedure. This introduces uncertainties in the inferred dust masses ($M_{\rm d }$), infrared luminosities ($L_{\rm IR}$), and obscured Star Formation Rate (SFR) fraction at $z > 4$. In this work we use a method that allows us to constrain $T_{\rm d}$ with a single band measurement by combining the $158\ \mathrm{μm}$ continuum information with the overlying [CII] emission line. We analyse the $21$ [CII] and FIR continuum detected $z\sim 5$ galaxies in ALPINE, finding a range of $T_{\rm d}=25-60\ \mathrm{K}$ and $M_{\rm d} = 0.6-25.1\ \times 10^{7}\ \mathrm{M_{\odot}}$. Given the measured stellar masses of ALPINE galaxies, the inferred dust yields are around $M_{\rm d}/M_{\star} = (0.2-8) \times 10^{-3}$, consistent with theoretical dust-production constraints. We find that $8$ out of $21$ ALPINE galaxies have $L_{\rm IR} \geq 10^{12}\ \mathrm{L_{\odot}}$, comparable to UltraLuminous IR Galaxies (ULIRGs). Relying on ultraviolet-to-optical SED fitting, the SFR was underestimated by up to $2$ orders of magnitude in $4$ of these $8$ ULIRGs-like galaxies. We conclude that these $4$ peculiar sources should be characterised by a two-phase interstellar medium structure with "spatially-segregated" FIR and ultraviolet emitting regions.
△ Less
Submitted 17 October, 2022;
originally announced October 2022.
-
The ALMA REBELS Survey: Average [CII] $158\,\rm{μm}$ sizes of Star-Forming Galaxies from $z\sim 7$ to $z\sim 4$
Authors:
Y. Fudamoto,
R. Smit,
R. A. A. Bowler,
P. A. Oesch,
R. Bouwens,
M. Stefanon,
H. Inami,
R. Endsley,
V. Gonzalez,
S. Schouws,
D. Stark,
H. S. B. Algera,
M. Aravena,
L. Barrufet,
E. da Cunha,
P. Dayal,
A. Ferrara,
L. Graziani,
J. A. Hodge,
A. P. S. Hygate,
A. K. Inoue,
T. Nanayakkara,
A. Pallottini,
E. Pizzati,
R. Schneider
, et al. (14 additional authors not shown)
Abstract:
We present the average [CII] $158\,\rm{μm}$ emission line sizes of UV-bright star-forming galaxies at $z\sim7$. Our results are derived from a stacking analysis of [CII] $158\,\rm{μm}$ emission lines and dust continua observed by ALMA, taking advantage of the large program Reionization Era Bright Emission Line Survey (REBELS). We find that the average [CII] emission at $z\sim7$ has an effective ra…
▽ More
We present the average [CII] $158\,\rm{μm}$ emission line sizes of UV-bright star-forming galaxies at $z\sim7$. Our results are derived from a stacking analysis of [CII] $158\,\rm{μm}$ emission lines and dust continua observed by ALMA, taking advantage of the large program Reionization Era Bright Emission Line Survey (REBELS). We find that the average [CII] emission at $z\sim7$ has an effective radius $r_e$ of $2.2\pm0.2\,\rm{kpc}$. It is $\gtrsim2\times$ larger than the dust continuum and the rest-frame UV emission, in agreement with recently reported measurements for $z\lesssim6$ galaxies. Additionally, we compared the average [CII] size with $4<z<6$ galaxies observed by the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). By analysing [CII] sizes of $4<z<6$ galaxies in two redshift bins, we find an average [CII] size of $r_{\rm e}=2.2\pm0.2\,\rm{kpc}$ and $r_{\rm e}=2.5\pm0.2\,\rm{kpc}$ for $z\sim5.5$ and $z\sim4.5$ galaxies, respectively. These measurements show that star-forming galaxies, on average, show no evolution in the size of the [CII] $158\,{\rm μm}$ emitting regions at redshift between $z\sim7$ and $z\sim4$. This finding suggest that the star-forming galaxies could be morphologically dominated by gas over a wide redshift range.
△ Less
Submitted 15 June, 2022; v1 submitted 3 June, 2022;
originally announced June 2022.
-
A survey of high-$z$ galaxies: SERRA simulations
Authors:
A. Pallottini,
A. Ferrara,
S. Gallerani,
C. Behrens,
M. Kohandel,
S. Carniani,
L. Vallini,
S. Salvadori,
V. Gelli,
L. Sommovigo,
V. D'Odorico,
F. Di Mascia,
E. Pizzati
Abstract:
We introduce SERRA, a suite of zoom-in high-resolution ($\sim 10\,\rm pc$) cosmological simulations including non-equilibrium chemistry and on-the-fly radiative transfer. The outputs are post-processed to derive galaxy UV+FIR continuum and emission line properties. Results are compared with available multi-wavelength data to constrain the physical properties (e.g., star formation rates, stellar/ga…
▽ More
We introduce SERRA, a suite of zoom-in high-resolution ($\sim 10\,\rm pc$) cosmological simulations including non-equilibrium chemistry and on-the-fly radiative transfer. The outputs are post-processed to derive galaxy UV+FIR continuum and emission line properties. Results are compared with available multi-wavelength data to constrain the physical properties (e.g., star formation rates, stellar/gas/dust mass, metallicity) of high-redshift $6 \lesssim z \lesssim 15$ galaxies. This flagship paper focuses on the $z=7.7$ sub-sample, including 202 galaxies with stellar mass $10^7 M_\odot \lesssim M_\star \lesssim 5\times 10^{10}M_\odot$, and specific star formation ranging from ${\rm sSFR} \sim 100\,{\rm Gyr}^{-1}$ in young, low-mass galaxies to $\sim 10\,{\rm Gyr}^{-1}$ for older, massive ones. At this redshift, SERRA galaxies are typically bursty, i.e. they are located above the Schmidt-Kennicutt relation by a factor $κ_s = 3.03^{+4.9}_{-1.8}$, consistent with recent findings for [OIII] and [CII] emitters at high-$z$. They also show relatively large ${\rm IRX} = L_{\rm FIR}/L_{\rm UV}$ values as a result of their compact/clumpy morphology effectively blocking the stellar UV luminosity. Note that this conclusion might be affected by insufficient spatial resolution at the molecular cloud level. We confirm that early galaxies lie on the standard $\rm [CII]-SFR$ relation; their observed $L_{\rm [OIII]}/L_{\rm [CII]} \simeq 1-10$ ratios can be reproduced by a part of the SERRA galaxies without the need of a top-heavy IMF and/or anomalous C/O abundances. [OI] line intensities are similar to local ones, making ALMA high-$z$ detections challenging but feasible ($\sim 6\,\rm hr$ for a SFR of $50\,M_\odot\,{\rm yr}^{-1}$).
△ Less
Submitted 4 May, 2022; v1 submitted 7 January, 2022;
originally announced January 2022.
-
Towards inference of overlapping gravitational wave signals
Authors:
Elia Pizzati,
Surabhi Sachdev,
Anuradha Gupta,
Bangalore Sathyaprakash
Abstract:
Merger rates of binary black holes, binary neutron stars, and neutron star-black hole binaries in the local Universe (i.e., redshift $z=0$), inferred from the Laser Interferometer Gravitational Wave Observatory (LIGO) and Virgo, are $16-130\, \mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$, $13-1900\, \mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$, and $7.4-320\,\mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$, respectively. Thes…
▽ More
Merger rates of binary black holes, binary neutron stars, and neutron star-black hole binaries in the local Universe (i.e., redshift $z=0$), inferred from the Laser Interferometer Gravitational Wave Observatory (LIGO) and Virgo, are $16-130\, \mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$, $13-1900\, \mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$, and $7.4-320\,\mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$, respectively. These rates suggest that there is a significant chance that two or more of these signals will overlap with each other during their lifetime in the sensitivity-band of future gravitational-wave detectors such as the Cosmic Explorer and Einstein Telescope. The detection pipelines provide the coalescence time of each signal with an accuracy $\mathcal{O}(10\,\rm ms)$. We show that using a prior on the coalescence time from a detection pipeline, it is possible to correctly infer the properties of these overlapping signals with the current data-analysis infrastructure. We study different configurations of two overlapping signals created by non-spinning binaries, varying their time and phase at coalescence, as well as their signal-to-noise ratios. We conclude that, for the scenarios considered in this work, parameter inference is robust provided that their coalescence times in the detector frame are more than $\sim 1-2 \,\mathrm{s}$. Signals whose coalescence epochs lie within $\sim 0.5\,\mathrm{s}$ of each other suffer from significant biases in parameter inference, and new strategies and algorithms would be required to overcome such biases.
△ Less
Submitted 27 April, 2022; v1 submitted 15 February, 2021;
originally announced February 2021.
-
Outflows and extended [CII] halos in high redshift galaxies
Authors:
Elia Pizzati,
Andrea Ferrara,
Andrea Pallottini,
Simona Gallerani,
Livia Vallini,
Davide Decataldo,
Seiji Fujimoto
Abstract:
Recent stacked ALMA observations have revealed that normal, star-forming galaxies at $z\approx 6$ are surrounded by extended ($\approx 10\,\mathrm{kpc}$) [CII] emitting halos which are not predicted by the most advanced, zoom-in simulations. We present a model in which these halos are the result of supernova-driven cooling outflows. Our model contains two free parameters, the outflow mass loading…
▽ More
Recent stacked ALMA observations have revealed that normal, star-forming galaxies at $z\approx 6$ are surrounded by extended ($\approx 10\,\mathrm{kpc}$) [CII] emitting halos which are not predicted by the most advanced, zoom-in simulations. We present a model in which these halos are the result of supernova-driven cooling outflows. Our model contains two free parameters, the outflow mass loading factor, $η$, and the parent galaxy dark matter halo circular velocity, $v_c$. The outflow model successfully matches the observed [CII] surface brightness profile if $η= 3.20 \pm 0.10$ and $v_c = 170 \pm 10{\,\rm km\,s^{-1}}$, corresponding to a dynamical mass of $\approx 10^{11}\, \mathrm{M}_\odot$. The predicted outflow rate and velocity range are $128 \pm 5 \,\mathrm{M}_\odot {\rm yr}^{-1}$ and $300-500 {\,\rm km\,s^{-1}}$, respectively. We conclude that: (a) extended halos can be produced by cooling outflows; (b) the large $η$ value is marginally consistent with starburst-driven outflows, but it might indicate additional energy input from AGN; (c) the presence of [CII] halos requires an ionizing photon escape fraction from galaxies $f_{\rm esc} \ll 1$. The model can be readily applied also to individual high-$z$ galaxies, as those observed, e.g., by the ALMA ALPINE survey now becoming available.
△ Less
Submitted 24 April, 2020; v1 submitted 28 January, 2020;
originally announced January 2020.