Euclid Quick Data Release (Q1): Hunting for luminous z > 6 galaxies in the Euclid Deep Fields -- forecasts and first bright detections

TL;DR

The Euclid Deep Fields will significantly enhance the detection and study of luminous galaxies at z > 6, enabling better constraints on the ultraviolet luminosity function's bright end through forecasts and initial detections.

Contribution

This paper provides forecasts for high-redshift galaxy counts in Euclid Deep Fields and demonstrates effective selection methods for bright z > 6 galaxies, including first candidate detections.

Findings

Forecasted >100,000 galaxies at z=6-12 in Euclid data.

High-fidelity selection method recovers >76% of z > 6 galaxies.

Detection of ultra-luminous candidates supports a double power-law LF model.

Abstract

The evolution of the rest-frame ultraviolet luminosity function (UV LF) is a powerful probe of early star formation and stellar mass build-up. At z > 6, its bright end (MUV < -21) remains poorly constrained due to the small volumes of existing near-infrared (NIR) space-based surveys. The Euclid Deep Fields (EDFs) will cover 53 deg^2 with NIR imaging down to 26.5 AB, increasing area by a factor of 100 over previous space-based surveys. They thus offer an unprecedented opportunity to select bright z > 6 Lyman break galaxies (LBGs) and constrain the UV LF's bright end. With NIR coverage extending to 2um, Euclid can detect galaxies out to z = 13. We present forecasts for the number densities of z > 6 galaxies expected in the final EDF dataset. Using synthetic photometry from spectral energy distribution (SED) templates of z = 5--15 galaxies, z = 1--4 interlopers, and Milky Way MLT dwarfs, we explore optimal selection methods for high-z LBGs. A combination of S/N cuts with SED fitting (from optical to MIR) yields the highest-fidelity sample, recovering >76% of input z > 6 LBGs while keeping low-z contamination <10%. This excludes instrumental artefacts, which will affect early Euclid releases. Auxiliary data are critical: optical imaging from the Hyper Suprime-Cam and Vera C. Rubin Observatory distinguishes genuine Lyman breaks, while Spitzer/IRAC data help recover z > 10 sources. Based on empirical double power-law LF models, we expect >100,000 LBGs at z = 6-12 and >100 at z > 12 in the final Euclid release. In contrast, steeper Schechter models predict no z > 12 detections. We also present two ultra-luminous (MUV < -23.5) candidates from the EDF-N Q1 dataset. If their redshifts are confirmed, their magnitudes support a DPL LF model at z > 9, highlighting Euclid's power to constrain the UV LF's bright end and identify the most luminous early galaxies for follow-up.