Title: A SPectroscopic Survey of Biased Halos in the Reionization Era (ASPIRE): Spectroscopically Complete Census of Obscured Cosmic Star Formation Rate Density at z  = 4–6
Abstract We present a stringent measurement of the dust-obscured star formation rate density (SFRD) atz= 4–6 from the ASPIRE JWST Cycle-1 medium and ALMA Cycle-9 large program. We obtained JWST/NIRCam grism spectroscopy and ALMA 1.2 mm continuum map along 25 independent quasar sightlines, covering a total survey area of  ∼35 arcmin2where we search for dusty star-forming galaxies (DSFGs) atz= 0–7. We identify eight DSFGs in seven fields atz= 4–6 through the detection of Hαor [O iii]λ5008 lines, including fainter lines such as Hβ, [O iii]λ4960, [N ii]λ6585, and [S ii]λλ6718,6733 for six sources. With this spectroscopically complete DSFG sample atz= 4–6 and negligible impact from cosmic variance (shot noise), we measure the infrared luminosity function (IRLF) down toLIR ∼ 2 × 1011L. We find flattening of IRLF atz= 4–6 towards the faint end (power-law slope α = 0.5 9 0.45 + 0.39 ). We determine the dust-obscured cosmic SFRD at this epoch to be log [ ρ SFR , IR / ( M yr 1 Mpc 3 ) ] = 1.5 2 0.13 + 0.14 . This is significantly higher than previous determinations using ALMA data in the Hubble Ultra Deep Field, which is void of DSFGs atz= 4–6 because of strong cosmic variance (shot noise). We conclude that the majority (66% ± 7%) of cosmic star formation atz ∼ 5 is still obscured by dust. We also discuss the uncertainty of SFRD propagated from far-IR spectral energy distribution and IRLF at the bright end, which will need to be resolved with future ALMA and JWST observations.  more » « less
Award ID(s):
2308258 2513040
PAR ID:
10609795
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; « less
Publisher / Repository:
ApJ
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
980
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
12
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We study the luminosity function (LF) and clustering properties of 888 Hαemitters (HAEs) at 3.75 < z < 6 in the GOODS-N field. The sample, built from JWST CONGRESS and FRESCO NIRCam grism surveys using a novel redshift assignment algorithm, spans ∼62 arcmin2and reachesL ∼ 1041.2erg s−1. We identify two prominent filamentary protoclusters atz ≈ 4.41 andz ≈ 5.19, hosting 98 and 144 HAEs, respectively. The observed HαLFs show similar shallow faint-end slopes for both protocluster and field galaxies at 3.75 < z < 5, and for the protocluster at 5 < z < 6 (α ≈ −1.2 to −1.3). In contrast, the field LF at 5 < z < 6 has a much steeper slope ( α = 1.8 7 0.23 + 0.30 ), suggesting that protocluster galaxies atz > 5 are more evolved, resembling the populations at 3.75 < z < 5. The observed star formation rate density from Hαintegrated down to 0.45Myr−1, is 0.05 0 0.003 + 0.002 M yr 1 Mpc 3 at 3.75 < z < 5 and 0.04 6 0.004 + 0.006 M yr 1 Mpc 3 at 5 < z < 6, with protoclusters contributing about 25% and 55%, respectively. This implies a large fraction of star formation atz> 4 occurs in protoclusters. For the first time, we conduct the star formation-rate-limited three-dimensional clustering analysis atz > 4. We find that the filamentary geometry of protoclusters flattens the power-law shape of the HAE autocorrelation functions, with slopes much shallower than the typically assumed value. The autocorrelation function of field HAEs has a correlation length of r 0 = 4.6 1 0.68 + 1.00 h 1 Mpc atz ≈ 4−5 and r 0 = 6.2 3 1.13 + 1.68 h 1 Mpc atz ≈ 5−6. Comparing the observed correlation functions with the UniverseMachinesimulation, we infer the dark matter (sub-)halo masses of HAEs to be log ( M h / M ) = 11.0 11.2 atz ≈ 4−6, with a scatter of 0.4 dex. 
    more » « less
  2. Abstract We present thez≈ 6 type-1 quasar luminosity function (QLF), based on the Pan-STARRS1 (PS1) quasar survey. The PS1 sample includes 125 quasars atz≈ 5.7–6.2, with −28 ≲M1450≲ −25. With the addition of 48 fainter quasars from the SHELLQs survey, we evaluate thez≈ 6 QLF over −28 ≲M1450≲ −22. Adopting a double power law with an exponential evolution of the quasar density (Φ(z) ∝ 10k(z−6);k= −0.7), we use a maximum likelihood method to model our data. We find a break magnitude of M * = 26.38 0.60 + 0.79 mag , a faint-end slope of α = 1.70 0.19 + 0.29 , and a steep bright-end slope of β = 3.84 1.21 + 0.63 . Based on our new QLF model, we determine the quasar comoving spatial density atz≈ 6 to be n ( M 1450 < 26 ) = 1.16 0.12 + 0.13 cGpc 3 . In comparison with the literature, we find the quasar density to evolve with a constant value ofk≈ −0.7, fromz≈ 7 toz≈ 4. Additionally, we derive an ionizing emissivity of ϵ 912 ( z = 6 ) = 7.23 1.02 + 1.65 × 10 22 erg s 1 Hz 1 cMpc 3 , based on the QLF measurement. Given standard assumptions, and the recent measurement of the mean free path by Becker et al. atz≈ 6, we calculate an Hiphotoionizing rate of ΓH I(z= 6) ≈ 6 × 10−16s−1, strongly disfavoring a dominant role of quasars in hydrogen reionization. 
    more » « less
  3. Abstract Dust-obscured galaxies (DOGs) containing central supermassive black holes (SMBHs) that are rapidly accreting (i.e., having high Eddington ratios,λEdd) may represent a key phase closest to the peak of both the black hole and galaxy growth in the coevolution framework for SMBHs and galaxies. In this work, we present a 68 ks XMM-Newton observation of the high-λEddDOG J1324+4501 atz∼ 0.8, which was initially observed by Chandra. We analyze the XMM-Newton spectra jointly with archival Chandra spectra. In performing a detailed X-ray spectral analysis, we find that the source is intrinsically X-ray luminous with log ( L X /erg s 1 ) = 44.71 0.12 + 0.08 and heavily obscured with log ( N H / cm 2 ) = 23.43 0.13 + 0.09 . We further utilize UV-to-IR archival photometry to measure and fit the source’s spectral energy distribution to estimate its host-galaxy properties. We present a supplementary comparison sample of 21 X-ray luminous DOGs from the XMM-SERVS survey with sufficient (>200) 0.5–10 keV counts to perform a similarly detailed X-ray spectral analysis. Of the X-ray luminous DOGs in our sample, we find that J1324+4501 is the most remarkable, possessing one of the highest X-ray luminosities, column densities, and star formation rates. We demonstrate that J1324+4501 is in an extreme evolutionary stage where SMBH accretion and galaxy growth are at their peaks. 
    more » « less
  4. Abstract We present a systematic study of the environments of 25 luminous quasars atz > 6.5 from the ASPIRE program. Using JWST/NIRCam wide-field slitless spectroscopy data, we identified 487 galaxies at 5.3 ≲ z ≲ 7.0 exhibiting [Oiii] emission. Among these, 122 [Oiii] emitters lie within ∣Δvlos∣ < 1000 km s−1of the quasars, corresponding to a ∼9.4-fold enhancement relative to the average galaxy density at other redshifts. Furthermore, we identified 16 [Cii]-emitting galaxies at the quasar redshifts from Atacama Large Millimeter/submillimeter Array (ALMA) mosaic observations. A cross-correlation function analysis between quasars and [Oiii]+[Cii] emitters yields a cross-correlation length of r 0 QG = 8.6 8 0.55 + 0.51 h 1 cMpc and an autocorrelation of r 0 QQ = 15.7 6 2.70 + 2.48 h 1 cMpc , indicating thatz ∼ 7 quasars reside in dark matter halos with M halo = 1 0 12.2 7 0.26 + 0.21 M and have a quasar lifetime of t Q = 1 0 7.05 1.01 + 0.95 yr . Notably, the number of [Oiii]-emitting galaxies at quasar redshifts varies significantly from field to field, ranging from 0 to 20, highlighting a diverse quasar environment. Remarkably, seven quasars trace significant galaxy overdensities (i.e., protoclusters), withδgal > 5 within a volume ofV ∼ 500 cMpc3. We also find that ∣Δvlos∣ increases rapidly toward smaller galaxy–quasar separations in protocluster fields, consistent with galaxy kinematics around extremely massive halos in cosmological simulations. By combining JWST and ALMA data, we reveal the complex and diverse environments of these early quasars, providing robust evidence that the earliest luminous quasars are effective tracers of galaxy overdensities, albeit with substantial field-to-field variation. 
    more » « less
  5. Abstract We present Atacama Large Millimeter/submillimeter Array observations of the [CI] 492 and 806 GHz fine-structure lines in 25 dusty star-forming galaxies (DSFGs) atz= 4.3 in the core of the SPT2349–56 protocluster. The protocluster galaxies exhibit a median L [ C I ] ( 2 1 ) / L [ C I ] ( 1 0 ) ratio of 0.94, with an interquartile range of 0.81–1.24. These ratios are markedly different to those observed in DSFGs in the field (across a comparable redshift and 850μm flux density range), where the median is 0.55, with an interquartile range of 0.50–0.76, and we show that this difference is driven by an excess of [Ci](2–1) in the protocluster galaxies for a given 850μm flux density. Assuming local thermal equilibrium, we estimate gas excitation temperatures of T ex = 59 . 1 6.8 + 8.1 K for our protocluster sample and T ex = 33 . 9 2.2 + 2.4 K for the field sample. Our main interpretation of this result is that the protocluster galaxies have had their cold gas driven to their cores via close-by interactions within the dense environment, leading to an overall increase in the average gas density and excitation temperature, as well as an elevated [Ci](2–1) luminosity-to-far-infrared-luminosity ratio. 
    more » « less