skip to main content

Title: Prediction of H α and [O iii] emission line galaxy number counts for future galaxy redshift surveys

We perform a simulation with Galacticus, a semi-analytical galaxy formation model, to predict the number counts of H α and [O iii] emitting galaxies. With a state-of-the-art N-body simulation, UNIT, we first calibrate Galacticus with the current observation of H α luminosity function. The resulting model coupled with a dust attenuation model, can reproduce the current observations, including the H α luminosity function from HiZELS and number density from WISP. We extrapolate the model prediction to higher redshift and the result is found to be consistent with previous investigations. We then use the same galaxy formation model to predict the number counts for [O iii] emitting galaxies. The result provides further validation of our galaxy formation model and dust model. We present number counts of H α and [O iii] emission line galaxies for three different line flux limits: 5 × 10−17erg s−1 cm−2, 1 × 10−16 erg s−1 cm−2 (6.5σ nominal depth for WFIRST GRS), and 2 × 10−16 erg s−1 cm−2 (3.5σ depth of Euclid GRS). At redshift 2 < z < 3, our model predicts that WFIRST can observe hundreds of [O iii] emission line galaxies per square degree with a line flux limit of 1 × 10−16 erg s−1 cm−2. This will provide accurate measurement of large-scale structure to probe dark energy over a huge cosmic volume to an unprecedented high redshift. more » Finally, we compare the flux ratio of H α/[O iii] within the redshift range of 0 < z < 3. Our results show the known trend of increasing H α/[O iii] flux ratio with H α flux at low redshift, which becomes a weaker trend at higher redshifts.

« less
 ;  ;  ;  ;  
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 3667-3678
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Euclid and the Roman Space Telescope (Roman) will soon use grism spectroscopy to detect millions of galaxies via their Hαand [Oiii]λ5007 emission. To better constrain the expected galaxy counts from these instruments, we use a vetted sample of 4239 emission-line galaxies from the 3D Hubble Space Telescope survey to measure the Hαand [Oiii]λ5007 luminosity functions between 1.16 <z< 1.90; this sample is ∼4 times larger than previous studies at this redshift. We find very good agreement with previous measurements for Hα, but for [Oiii], we predict a higher number of intermediate-luminosity galaxies than from previous works. We find that, for both lines, the characteristic luminosity,*, increases monotonically with redshift, and use the Hαluminosity function to calculate the epoch’s cosmic star formation rate density. We find that Hα-visible galaxies account for ∼81% of the epoch’s total star formation rate, and this value changes very little over the 1.16 <z< 1.56 redshift range. Finally, we derive the surface density of galaxies as a function of limiting flux and find that previous predictions for galaxy counts for the Euclid Wide Survey are unchanged, but there may be more [Oiii] galaxies in the Roman High Latitude Survey than previouslymore »estimated.

    « less
  2. Exploiting the sensitivity of the IRAM NOrthern Extended Millimeter Array (NOEMA) and its ability to process large instantaneous bandwidths, we have studied the morphology and other properties of the molecular gas and dust in the star forming galaxy, H-ATLAS J131611.5+281219 (HerBS-89a), at z = 2.95. High angular resolution (0 . ″3) images reveal a partial 1 . ″0 diameter Einstein ring in the dust continuum emission and the molecular emission lines of 12 CO(9−8) and H 2 O(2 02  − 1 11 ). Together with lower angular resolution (0 . ″6) images, we report the detection of a series of molecular lines including the three fundamental transitions of the molecular ion OH + , namely (1 1  − 0 1 ), (1 2  − 0 1 ), and (1 0  − 0 1 ), seen in absorption; the molecular ion CH + (1 − 0) seen in absorption, and tentatively in emission; two transitions of amidogen (NH 2 ), namely (2 02  − 1 11 ) and (2 20  − 2 11 ) seen in emission; and HCN(11 − 10) and/or NH(1 2  − 0 1 ) seen in absorption. The NOEMA data are complemented with Very Large Array data tracing the 12 CO(1 − 0) emission line, which provides a measurement ofmore »the total mass of molecular gas and an anchor for a CO excitation analysis. In addition, we present Hubble Space Telescope imaging that reveals the foreground lensing galaxy in the near-infrared (1.15  μ m). Together with photometric data from the Gran Telescopio Canarias, we derive a photometric redshift of z phot = 0.9 −0.5 +0.3 for the foreground lensing galaxy. Modeling the lensing of HerBS-89a, we reconstruct the dust continuum (magnified by a factor μ  ≃ 5.0) and molecular emission lines (magnified by μ  ∼ 4 − 5) in the source plane, which probe scales of ∼0 . ″1 (or 800 pc). The 12 CO(9 − 8) and H 2 O(2 02  − 1 11 ) emission lines have comparable spatial and kinematic distributions; the source-plane reconstructions do not clearly distinguish between a one-component and a two-component scenario, but the latter, which reveals two compact rotating components with sizes of ≈1 kpc that are likely merging, more naturally accounts for the broad line widths observed in HerBS-89a. In the core of HerBS-89a, very dense gas with n H 2  ∼ 10 7 − 9 cm −3 is revealed by the NH 2 emission lines and the possible HCN(11 − 10) absorption line. HerBS-89a is a powerful star forming galaxy with a molecular gas mass of M mol  = (2.1 ± 0.4) × 10 11   M ⊙ , an infrared luminosity of L IR  = (4.6 ± 0.4) × 10 12   L ⊙ , and a dust mass of M dust  = (2.6 ± 0.2) × 10 9   M ⊙ , yielding a dust-to-gas ratio δ GDR  ≈ 80. We derive a star formation rate SFR = 614 ± 59  M ⊙ yr −1 and a depletion timescale τ depl  = (3.4 ± 1.0) × 10 8 years. The OH + and CH + absorption lines, which trace low (∼100 cm −3 ) density molecular gas, all have their main velocity component red-shifted by Δ V  ∼ 100 km s −1 relative to the global CO reservoir. We argue that these absorption lines trace a rare example of gas inflow toward the center of a galaxy, indicating that HerBS-89a is accreting gas from its surroundings.« less
  3. ABSTRACT We present predictions for high redshift (z = 2−10) galaxy populations based on the IllustrisTNG simulation suite and a full Monte Carlo dust radiative transfer post-processing. Specifically, we discuss the H α and H β + $[\rm O \,{\small III}]$ luminosity functions up to z = 8. The predicted H β + $[\rm O \,{\small III}]$ luminosity functions are consistent with present observations at z ≲ 3 with ${\lesssim} 0.1\, {\rm dex}$ differences in luminosities. However, the predicted H α luminosity function is ${\sim }0.3\, {\rm dex}$ dimmer than the observed one at z ≃ 2. Furthermore, we explore continuum spectral indices, the Balmer break at 4000 Å; (D4000) and the UV continuum slope β. The median D4000 versus specific star formation rate relation predicted at z = 2 is in agreement with the local calibration despite a different distribution pattern of galaxies in this plane. In addition, we reproduce the observed AUV versus β relation and explore its dependence on galaxy stellar mass, providing an explanation for the observed complexity of this relation. We also find a deficiency in heavily attenuated, UV red galaxies in the simulations. Finally, we provide predictions for the dust attenuation curves of galaxies at z = 2−6 and investigate their dependence on galaxy colours andmore »stellar masses. The attenuation curves are steeper in galaxies at higher redshifts, with bluer colours, or with lower stellar masses. We attribute these predicted trends to dust geometry. Overall, our results are consistent with present observations of high-redshift galaxies. Future James Webb Space Telecope observations will further test these predictions.« less
  4. Abstract

    The Baldwin, Philips, & Terlevich diagram of [Oiii]/Hβversus [Nii]/Hα(hereafter N2-BPT) has long been used as a tool for classifying galaxies based on the dominant source of ionizing radiation. Recent observations have demonstrated that galaxies atz∼ 2 reside offset from local galaxies in the N2-BPT space. In this paper, we conduct a series of controlled numerical experiments to understand the potential physical processes driving this offset. We model nebular line emission in a large sample of galaxies, taken from thesimbacosmological hydrodynamic galaxy formation simulation, using thecloudyphotoionization code to compute the nebular line luminosities from Hiiregions. We find that the observed shift toward higher [Oiii]/Hβand [Nii]/Hαvalues at high redshift arises from sample selection: when we consider only the most massive galaxiesM*∼ 1010–11M, the offset naturally appears, due to their high metallicities. We predict that deeper observations that probe lower-mass galaxies will reveal galaxies that lie on a locus comparable toz∼ 0 observations. Even when accounting for samples-selection effects, we find that there is a subtle mismatch between simulations and observations. To resolve this discrepancy, we investigate the impact of varying ionization parameters, Hiiregion densities, gas-phase abundance patterns, and increasing radiation field hardness on N2-BPT diagrams. We find that either decreasing themore »ionization parameter or increasing the N/O ratio of galaxies at fixed O/H can move galaxies along a self-similar arc in N2-BPT space that is occupied by high-redshift galaxies.

    « less
  5. ABSTRACT We present a systematic investigation of physical conditions and elemental abundances in four optically thick Lyman-limit systems (LLSs) at z = 0.36–0.6 discovered within the cosmic ultraviolet baryon survey (CUBS). Because intervening LLSs at z < 1 suppress far-UV (ultraviolet) light from background QSOs, an unbiased search of these absorbers requires a near-UV-selected QSO sample, as achieved by CUBS. CUBS LLSs exhibit multicomponent kinematic structure and a complex mix of multiphase gas, with associated metal transitions from multiple ionization states such as C ii, C iii, N iii, Mg ii, Si ii, Si iii, O ii, O iii, O vi, and Fe ii absorption that span several hundred km s−1 in line-of-sight velocity. Specifically, higher column density components (log N(H i)/cm−2≳ 16) in all four absorbers comprise dynamically cool gas with $\langle T \rangle =(2\pm 1) \times 10^4\,$K and modest non-thermal broadening of $\langle b_\mathrm{nt} \rangle =5\pm 3\,$km s−1. The high quality of the QSO absorption spectra allows us to infer the physical conditions of the gas, using a detailed ionization modelling that takes into account the resolved component structures of H i and metal transitions. The range of inferred gas densities indicates that these absorbers consist of spatially compact clouds with a median line-of-sight thickness of $160^{+140}_{-50}$ pc. While obtaining robust metallicitymore »constraints for the low density, highly ionized phase remains challenging due to the uncertain $N\mathrm{(H\, {\small I})}$, we demonstrate that the cool-phase gas in LLSs has a median metallicity of $\mathrm{[\alpha /H]_{1/2}}=-0.7^{+0.1}_{-0.2}$, with a 16–84 percentile range of [α/H] = (−1.3, −0.1). Furthermore, the wide range of inferred elemental abundance ratios ([C/α], [N/α], and [Fe/α]) indicate a diversity of chemical enrichment histories. Combining the absorption data with deep galaxy survey data characterizing the galaxy environment of these absorbers, we discuss the physical connection between star-forming regions in galaxies and diffuse gas associated with optically thick absorption systems in the z < 1 circumgalactic medium.« less