Euclid Quick Data Release (Q1): Identification of massive galaxy candidates at the end of the Epoch of Reionisation

Abstract: Probing the presence and properties of massive galaxies at high redshift is one of the most critical tests for galaxy formation models. In this work, we search for galaxies with stellar masses M* > 10^10.25 Msun at z in [5,7], i.e., towards the end of the Epoch of Reionisation, over a total of ~23 deg^2 in two of the Euclid Quick Data Release (Q1) fields: the Euclid Deep Field North and Fornax (EDF-N and EDF-F). In addition to the Euclid photometry, we incorporate Spitzer Infrared Camera (IRAC) and ground-based optical data to perform spectral energy distribution (SED) fitting, obtaining photometric redshifts and derived physical parameters. After applying rigorous selection criteria, we identify a conservative sample of 145 candidate massive galaxies with M* > 10^10.25 Msun at z in [5,7], including 5 objects with M* > 10^11 Msun. This makes for a surface density of about 6.3 deg^-2 at z in [5,7], which should be considered a lower limit because of the current depth of the Euclid data (H_E < 24, 5 sigma in Q1). We find that the inferred stellar masses are consistent with galaxy formation models with standard star-formation efficiencies. These massive galaxies have colour excess E(B-V) values up to 0.75, indicating significant dust attenuation in some of them. In addition, half of the massive galaxies have best-fit ages comparable to the age of the Universe at those redshifts, which suggests that their progenitors were formed very early in cosmic time. About 78% of the massive galaxies lie on the star-forming main sequence (MS) in the SFR-M* plane, ~12% are found in the starburst region, and 10% in the transition zone between the MS and starbursts. We find no significant evidence for outshining or AGN contamination that could account for the elevated specific star-formation rates (sSFR) observed in the ~12% of galaxies classified as starbursts.