Photometric Redshift Selection Research Articles

SIDE-real: Supernova Ia Dust Extinction with truncated marginal neural ratio estimation applied to real data

Abstract We present the first fully simulation-based hierarchical analysis of the light curves of a population of low-redshift type Ia supernov aelig; (SNaelig; Ia). Our hardware-accelerated forward model, released in the Python package slicsim, includes stochastic variations of each SN’s spectral flux distribution (based on the pre-trained BayeSN model), extinction from dust in the host and in the Milky Way, redshift, and realistic instrumental noise. By utilising truncated marginal neural ratio estimation (TMNRE), a neural network-enabled simulation-based inference technique, we implicitly marginalise over 4000 latent variables (for a set of ≈ 100 SNaelig; Ia to efficiently infer SN Ia absolute magnitudes and host-galaxy dust properties at the population level while also constraining the parameters of individual objects. Amortisation of the inference procedure allows us to obtain coverage guarantees for our results through Bayesian validation and frequentist calibration. Furthermore, we show a detailed comparison to full likelihood-based inference, implemented through Hamiltonian Monte Carlo, on simulated data and then apply TMNRE to the light curves of 86 SNaelig; Ia from the Carnegie Supernova Project,, deriving marginal posteriors in excellent agreement with previous work. Given its ability to accommodate arbitrarily complex extensions to the forward model—e.g. different populations based on host properties, redshift evolution, complicated photometric redshift estimates, selection effects, and non-la contamination—without significant modifications to the inference procedure, TMNRE has the potential to become the tool of choice for cosmological parameter inference from future, large SN Ia samples.

NGDEEP Epoch 1: The Faint End of the Luminosity Function at z ∼ 9–12 from Ultradeep JWST Imaging

We present a robust sample of very high redshift galaxy candidates from the first epoch of JWST/NIRCam imaging from the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) survey. The NGDEEP NIRCam imaging, spanning 9.7 arcmin2 in the Hubble Ultra Deep Field Parallel Field 2, reaches m = 30.4 (5σ, point-source, 2″ diameter apertures corrected to total) in F277W, making it the deepest public JWST GO imaging data set to date. We describe our detailed data reduction process of the six-filter broadband JWST/NIRCam imaging, incorporating custom corrections for systematic effects to produce high-quality calibrated images. Using robust photometric redshift selection criteria, we identify a sample of 38 z ≳ 9 galaxy candidates. These objects span a redshift range of z = 8.5–15.8 and apparent magnitudes of m F277W = 27–30.5 AB mag, reaching ∼1.5 mag deeper than previous public JWST imaging surveys. We calculate the rest-frame ultraviolet luminosity function at z ∼ 9 and 11 and present a new measurement of the luminosity function faint-end slope at z ∼ 11. We find a faint-end slope of α = −2.5 ± 0.4 and −2.2 ± 0.2 at z ∼ 9 and 11, respectively. This is consistent with no significant evolution in the faint-end slope and number density from z = 9 to 11. Comparing our results with theoretical predictions, we find that some models produce better agreement at the faint end than the bright end. These results will help to constrain how stellar feedback impacts star formation at these early epochs.

Physical Characterization of Early Galaxies in the Webb’s First Deep Field SMACS J0723.3-7327

This paper highlights initial photometric analyses of JWST NIRCam imaging data in the sightline of SMACS0723, aiming to identify galaxies at redshift z > 7. By applying a conservative Lyman-break selection followed by photometric-redshift analysis and visual inspection, we identify four F090W-dropout and two F150W-dropout sources, three of which were recently confirmed in an independent spectroscopic analysis to z = 7.663, 7.665, and 8.499. We then supplement our sample with a photometric-redshift selection, and identify five additional candidates at 7 < z phot < 13. The NIRCam images clearly resolve all sources and reveal their subgalactic components that were not resolved/detected in the previous imaging by Hubble Space Telescope. Our spectral energy distribution analysis reveals that the selected galaxies are characterized by young stellar populations (median age of ∼50 Myr) of subsolar metallicity (∼0.2 Z ⊙) and little dust attenuation (A V ∼ 0.5). In several cases, we observe extreme Hβ+[O iii] lines being captured in the F444W band and seen as color excess, which is consistent with their observed high star formation rate surface density. Eight of the 11 sources identified in this study appear in at least one of the recent studies (Adams et al.; Atek et al.; Donnan et al.; Harikane et al.; Yan et al.) of the same fields, implying the high fidelity of our selection. We crossmatch all high-z galaxy candidates presented in the five studies with our catalog and discuss the possible causes of discrepancy in the final lists.

Early Results from GLASS-JWST. X. Rest-frame UV-optical Properties of Galaxies at 7 < z < 9

We present the first James Webb Space Telescope/NIRCam-led determination of 7 < z < 9 galaxy properties based on broadband imaging from 0.8 to 5 μm as part of the GLASS-JWST Early Release Science program. This is the deepest data set acquired at these wavelengths to date, with an angular resolution ≲0.″14. We robustly identify 13 galaxies with signal-to-noise ratio ≳ 8 in F444W from 8 arcmin2 of data at m AB ≤ 28 from a combination of dropout and photometric redshift selection. From simulated data modeling, we estimate the dropout sample purity to be ≳90%. We find that the number density of these F444W-selected sources is broadly consistent with expectations from the UV luminosity function determined from Hubble Space Telescope data. We characterize galaxy physical properties using a Bayesian spectral energy distribution fitting method, finding a median stellar mass of 108.5 M ⊙ and age 140 Myr, indicating they started ionizing their surroundings at redshift z > 9.5. Their star formation main sequence is consistent with predictions from simulations. Lastly, we introduce an analytical framework to constrain main-sequence evolution at z > 7 based on galaxy ages and basic assumptions, through which we find results consistent with expectations from cosmological simulations. While this work only gives a glimpse of the properties of typical galaxies that are thought to drive the reionization of the universe, it clearly shows the potential of JWST to unveil unprecedented details of galaxy formation in the first billion years.

Cosmology with dropout selection: straw-man surveys & CMB lensing

We seek to prove the means, motive and opportunity of 'dropout' selected 2 ⩽ z ⩽ 5 galaxies for large-scale structure. Together with acquired low-z tracers, these samples would map practically every linear mode and facilitate a tomographic decomposition of the Cosmic Microwave Background (CMB) lensing kernel over an unprecedented volume. With this, one may infer (the time evolution of) matter density fluctuations and perform compelling tests of horizon-scale General Relativity, neutrino masses and Inflation—viz., curvature, running of the spectral index and a scale-dependent halo bias induced by (local) primordial non-Gaussianity. This is facilitated by the order-of-magnitude increase in sensitivity achieved by planned CMB, optical-to-near-infrared imaging and spectroscopy. Focusing on traditional color-color—rather than photometric redshift—selection, we estimate the expected completeness, contamination, and spectroscopic survey speed of tailored Lyman-break galaxy (LBG) samples. With these in hand, we forecast the potential of CMB lensing cross-correlation, `clustering redshifts' and Redshift-Space Distortions (RSD) analyses. In particular, we estimate: the depth dependence of interlopers based on CFHTLS-Archive-Research Survey (CARS) data and propagate this to biases in cosmology; a simple relation for the dependence of the linear galaxy bias on redshift and depth; new inferences of non-linear halo bias at these redshifts using legacy data; detailed forecasts of LBG spectra as would be observed by the Dark Energy Spectroscopic Instrument, Prime Focus Spectrograph, and their successors. We further assess the relative competitiveness of these spectroscopic facilities based on an intuitive figure-of-merit and define modern equivalents to traditional color selection criteria for the Large Synoptic Survey Telescope, where necessary. We confirm these science cases to be compelling for achievable facilities in the next decade by defining a LBG sample of increasing Lyman-α equivalent width with redshift, which delivers both percent-level RSD constraints on the growth rate at high-z and measurements of CMB lensing cross-correlation at z=3 and 4 with a significance measured in the hundreds, given sufficient area overlap. Finally, we discuss the limitations of this initial exploration and provide avenues for future investigation.

Consistent cosmic shear in the face of systematics: a B-mode analysis of KiDS-450, DES-SV and CFHTLenS

We analyse three public cosmic shear surveys; the Kilo-Degree Survey (KiDS-450), the Dark Energy Survey (DES-SV) and the Canada France Hawaii Telescope Lensing Survey (CFHTLenS). Adopting the “COSEBIs” statistic to cleanly and completely separate the lensing E-modes from the non-lensing B-modes, we detect B-modes in KiDS-450 and CFHTLenS at the level of ∼2.7σ. For DES-SV we detect B-modes at the level of 2.8σ in a non-tomographic analysis, increasing to a 5.5σB-mode detection in a tomographic analysis. In order to understand the origin of these detected B-modes we measure the B-mode signature of a range of different simulated systematics including PSF leakage, random but correlated PSF modelling errors, camera-based additive shear bias and photometric redshift selection bias. We show that any correlation between photometric-noise and the relative orientation of the galaxy to the point-spread-function leads to an ellipticity selection bias in tomographic analyses. This work therefore introduces a new systematic for future lensing surveys to consider. We find that the B-modes in DES-SV appear similar to a superposition of the B-mode signatures from all of the systematics simulated. The KiDS-450 and CFHTLenS B-mode measurements show features that are consistent with a repeating additive shear bias.

The XXL Survey

Context. The luminosity function (LF) is a powerful statistical tool used to describe galaxies and learn about their evolution. In particular, the LFs of galaxies inside clusters allow us to better understand how galaxies evolve in these dense environments. Knowledge of the LFs of galaxies in clusters is also crucial for clusters studies in the optical and near-infrared (NIR) as they encode, along with their density profiles, most of their observational properties. However, no consensus has been reached yet about the evolution of the cluster galaxy LF with halo mass and redshift. Aims. The main goal of this study is to investigate the LF of a sample of 142 X-ray selected clusters, with spectroscopic redshift confirmation and a well defined selection function, spanning a wide redshift and mass range, and to test the LF dependence on cluster global properties, in a homogeneous and unbiased way. Methods. Our study is based on the Canada–France–Hawaii Telescope Legacy Survey (CFHTLS) photometric galaxy catalogue, associated with photometric redshifts. We constructed LFs inside a scaled radius using a selection in photometric redshift around the cluster spectroscopic redshift in order to reduce projection effects. The width of the photometric redshift selection was carefully determined to avoid biasing the LF and depended on both the cluster redshift and the galaxy magnitudes. The purity was then enhanced by applying a precise background subtraction. We constructed composite luminosity functions (CLFs) by stacking the individual LFs and studied their evolution with redshift and richness, analysing separately the brightest cluster galaxy (BCG) and non-BCG members. We fitted the dependences of the CLFs and BCG distributions parameters with redshift and richness conjointly in order to distinguish between these two effects. Results. We find that the usual photometric redshift selection methods can bias the LF estimate if the redshift and magnitude dependence of the photometric redshift quality is not taken into account. Our main findings concerning the evolution of the galaxy luminosity distribution with redshift and richness are that, in the inner region of clusters and in the redshift-mass range we probe (about 0 &lt; z &lt; 1 and 1013 M⊙ &lt; M500 &lt; 5 × 1014 M⊙), the bright part of the LF (BCG excluded) does not depend much on mass or redshift except for its amplitude, whereas the BCG luminosity increases both with redshift and richness.

VIMOS Ultra-Deep Survey (VUDS): IGM transmission towards galaxies with 2.5 &lt; z &lt; 5.5 and the colour selection of high-redshift galaxies

The observed UV rest-frame spectra of distant galaxies are the result of their intrinsic emission combined with absorption along the line of sight produced by the inter-galactic medium (IGM). Here we analyse the evolution of the mean IGM transmission Tr(Lyα) and its dispersion along the line of sight for 2127 galaxies with 2.5 &lt; z &lt; 5.5 in the VIMOS Ultra Deep Survey (VUDS). We fitted model spectra combined with a range of IGM transmission to the galaxy spectra using the spectral fitting algorithm GOSSIP+. We used these fits to derive the mean IGM transmission towards each galaxy for several redshift slices from z = 2.5 to z = 5.5. We found that the mean IGM transmission defined as Tr(Lyα) = e− τ (with τ as the HI optical depth) is 79%, 69%, 59%, 55%, and 46% at redshifts 2.75, 3.22, 3.70, 4.23, and 4.77, respectively. We compared these results to measurements obtained from quasar lines of sight and found that the IGM transmission towards galaxies is in excellent agreement with quasar values up to redshift z ~ 4. We found tentative evidence for a higher IGM transmission at z ≥ 4 compared to results from QSOs, but a degeneracy between dust extinction and IGM prevents us from firmly concluding whether the internal dust extinction for star-forming galaxies at z &gt; 4 takes a mean value significantly in excess of E(B−V) &gt; 0.15. Most importantly, we found a large dispersion of IGM transmission along the lines of sight towards distant galaxies with 68% of the distribution within 10 to 17% of the median value in δz = 0.5 bins, similar to what is found on the lines of sight towards QSOs. We demonstrate that taking this broad range of IGM transmission into account is important when selecting high-redshift galaxies based on their colour properties (e.g. LBG or photometric redshiftselection) because failing to do so causes a significant incompleteness in selecting high-redshift galaxy populations. We finally discuss the observed IGM properties and speculate that the broad range of observed transmissions might be the result of cosmic variance and clustering along lines of sight. This clearly shows that the sources that cause this extinction need to be more completely modelled.

Intrinsic alignments of group and cluster galaxies in photometric surveys

Intrinsic alignments of galaxies have been shown to contaminate weak gravitational lensing observables on linear scales, r > 10 h−1 Mpc, but studies of alignments in the non-linear regime have thus far been inconclusive. We present an estimator for extracting the intrinsic alignment signal of galaxies around stacked clusters of galaxies from multiband imaging data. Our estimator removes the contamination caused by galaxies that are gravitationally lensed by the clusters and scattered in redshift space due to photometric redshift uncertainties. It uses posterior probability distributions for the redshifts of the galaxies in the sample and it is easily extended to obtain the weak gravitational lensing signal while removing the intrinsic alignment contamination. We apply this algorithm to groups and clusters of galaxies identified in the Sloan Digital Sky Survey ‘Stripe 82’ co-added imaging data over ∼150 deg2. We find that the intrinsic alignment signal around stacked clusters in the redshift range 0.1 < z < 0.4 is consistent with zero. In terms of the tidal alignment model of Catelan et al., we set joint constraints on the strength of the alignment and the bias of the lensing groups and clusters on scales between 0.1 and 10 h−1 Mpc, |$b_{\rm L}C_1\rho _{\rm crit} = -2_{-14}^{+14} \times 10^{-4}$|⁠. This constrains the contamination fraction of alignment to lensing signal to the range between [ − 18, 23] per cent below scales of 1 h−1 Mpc at 95 per cent confidence level, and this result depends on our photometric redshift quality and selection criteria used to identify background galaxies. Our results are robust to the choice of photometric band in which the shapes are measured (i and r) and to centring on the Brightest Cluster Galaxy or on the geometrical centre of the clusters.

The large area KX quasar catalogue - I. Analysis of the photometric redshift selection and the complete quasar catalogue

The results of a large area, ~600 deg^2, K-band flux-limited spectroscopic survey for luminous quasars are presented. The survey utilises the UKIRT Infrared Deep Sky Survey (UKIDSS) Large Area Survey (LAS) in regions of sky within the Sloan Digital Sky Survey (SDSS) footprint. The K-band excess (KX) of all quasars with respect to Galactic stars is exploited in combination with a photometric redshift/classification scheme to identify quasar candidates for spectroscopic follow-up observations. The data contained within this investigation will be able to provide new constraints on the fraction of luminous quasars reddened by dust with E(B-V)<=0.5 mag. The spectroscopic sample is defined using the K-band, 14.0<=K<=16.6, and SDSS i-band limits of i=19.5, 19.7 and 22.0 over sky areas of 287, 150 and 196 deg^2, respectively. The survey includes >3200 known quasars from the SDSS and more than 250 additional confirmed quasars from the KX-selection. A well-defined sub-sample of quasars in the redshift interval 1.0<=z<=3.5 includes 1152 objects from the SDSS and 172 additional KX-selected quasars. The quasar selection is >95 per cent complete with respect to known SDSS quasars and >95 per cent efficient, largely independent of redshift and i-band magnitude. The properties of the new KX-selected quasars confirm the known redshift-dependent effectiveness of the SDSS quasar selection and provide a sample of luminous quasars experiencing intermediate levels of extinction by dust. The catalogue represents an important step towards the assembly of a well-defined sample of luminous quasars that may be used to investigate the properties of quasars experiencing intermediate levels of dust extinction within their host galaxies or due intervening absorption line systems.

Comparison of the properties of two fossil groups of galaxies with the normal group NGC 6034 based on multiband imaging and optical spectroscopy

We collected multiband imaging and spectroscopy for two fossil groups (RX J1119.7+2126 and 1RXS J235814.4+150524) and one normal group (NGC 6034). We computed photometric redshifts in the central zones of each group, combining previous data with the SDSS five-band data. For each group we investigated the red sequence (RS) of the color-magnitude relation and computed the luminosity functions, stellar population ages and distributions of the group members. Spectroscopy allowed us to investigate the large-scale surroundings of these groups and the substructure levels in 1RXS J235814.4+150524 and NGC 6034. The large-scale environment of 1RXS J235814.4+150524 is poor, though its galaxy density map shows a clear signature of the surrounding cosmic web. RX J1119.7+2126 appears to be very isolated, while the cosmic environment of NGC 6034 is very rich. At the group scale, 1RXS J235814.4+150524 shows no substructure. Galaxies with recent stellar populations seem preferentially located in the group outskirts. A RS is discernable for all three groups in a color-magnitude diagram. The luminosity functions based on photometric redshift selection and on statistical background subtraction have comparable shapes, and agree with the few points obtained from spectroscopic redshifts. These luminosity functions show the expected dip between first and second brightest galaxies for the fossil groups only. Their shape is also regular and relatively flat at faint magnitudes down to the completeness level for RX J1119.7+2126 and NGC 6034, while there is a clear lack of faint galaxies for 1RXS J235814.4+150524. RX J1119.7+2126 is definitely classified as a fossil group; 1RXS J235814.4+150524 also has properties very close to those of a fossil group, while we confirm that NGC 6034 is a normal group.

USING THE 1.6 μm BUMP TO STUDY REST-FRAME NEAR-INFRARED-SELECTED GALAXIES AT REDSHIFT 2

We explore the feasibility and limitations of using the 1.6um bump as a photometric redshift indicator and selection technique and use it to study the rest-frame H-band galaxy luminosity and stellar mass functions at redshift z~2. We use publicly available Spitzer/IRAC images in the GOODS fields and find that color selection in the IRAC bandpasses alone is comparable in completeness and contamination to BzK selection. We find that the shape of the 1.6um bump is robust, and photometric redshifts are not greatly affected by choice of model parameters. Comparison with spectroscopic redshifts shows photometric redshifts to be reliable. We create a rest-frame NIR selected catalog of galaxies at z~2 and construct a galaxy stellar mass function (SMF). Comparisons with other SMFs at approximately the same redshift but determined using shorter wavelengths show good agreement. This agreement suggests that selection at bluer wavelengths does not miss a significant amount of stellar mass in passive galaxies. Comparison with SMFs at other redshifts shows evidence for the downsizing scenario of galaxy evolution. We conclude by pointing out the potential for using the 1.6um technique to select high-redshift galaxies with the JWST, whose lambda > 0.6 um coverage will not be well suited to selecting galaxies using techniques that require imaging at shorter wavelengths.

Photometric redshifts as a tool for studying the Coma cluster galaxy populations

Aims. We apply photometric redshift techniques to an investigation of the Coma cluster galaxy luminosity function (GLF) at faint magnitudes, in particular in the u* band where basically no studies are presently available at these magnitudes. Methods: Cluster members were selected based on probability distribution function from photometric redshift calculations applied to deep u^*, B, V, R, I images covering a region of almost 1 deg2 (completeness limit R ~ 24). In the area covered only by the u* image, the GLF was also derived after a statistical background subtraction. Results: Global and local GLFs in the B, V, R, and I bands obtained with photometric redshift selection are consistent with our previous results based on a statistical background subtraction. The GLF in the u* band shows an increase in the faint end slope towards the outer regions of the cluster. The analysis of the multicolor type spatial distribution reveals that late type galaxies are distributed in clumps in the cluster outskirts, where X-ray substructures are also detected and where the GLF in the u* band is steeper. Conclusions: We can reproduce the GLFs computed with classical statistical subtraction methods by applying a photometric redshift technique. The u* GLF slope is steeper in the cluster outskirts, varying from α ~ -1 in the cluster center to α ~ -2 in the cluster periphery. The concentrations of faint late type galaxies in the cluster outskirts could explain these very steep slopes, assuming a short burst of star formation in these galaxies when entering the cluster. Based on observations obtained with MegaPrime/MegaCam, a joint project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii Telescope (CFHT) which is operated by the National Research Council (NRC) of Canada, the Institut National des Sciences de l'Univers of the Centre National de la Recherche Scientifique (CNRS) of France, and the University of Hawaii. This work is also partly based on data products produced at TERAPIX and the Canadian Astronomy Data Centre as part of the Canada-France-Hawaii Telescope Legacy Survey, a collaborative project of NRC and CNRS. Also based on data from W. M. Keck Observatory which is operated as a scientific partnership between the California Institute of Technology, the University of California, and NASA. It was made possible by the generous financial support of the W. M. Keck Foundation.