We present the Flagship galaxy mock, a simulated catalogue of billions of galaxies designed to support the scientific exploitation of the Euclid mission. Euclid is a medium-class mission of the European Space Agency optimised to determine the properties of dark matter and dark energy on the largest scales of the Universe. It probes structure formation over more than 10 billion years primarily from the combination of weak gravitational lensing and galaxy clustering data. The breath of Euclid's data will also foster a wide variety of scientific analyses. The Flagship simulation was developed to provide a realistic approximation to the galaxies that will be observed by Euclid and used in its scientific analyses. We ran a state-of-the-art N-body simulation with four trillion particles, producing a lightcone on the fly. From the dark matter particles, we produced a catalogue of 16 billion haloes in one octant of the sky in the lightcone up to redshift z=3. We then populated these haloes with mock galaxies using a halo occupation distribution and abundance matching approach, calibrating the free parameters of the galaxy mock against observed correlations and other basic galaxy properties. Modelled galaxy properties include luminosity and flux in several bands, redshifts, positions and velocities, spectral energy distributions, shapes and sizes, stellar masses, star formation rates, metallicities, emission line fluxes, and lensing properties. We selected a final sample of 3.4 billion galaxies with a magnitude cut of H_E<26, where we are complete. We have performed a comprehensive set of validation tests to check the similarity to observational data and theoretical models. In particular, our catalogue is able to closely reproduce the main characteristics of the weak lensing and galaxy clustering samples to be used in the mission's main cosmological analysis. (abridged)

The 2-point correlation function of the galaxy spatial distribution is a major cosmological observable that enables constraints on the dynamics and geometry of the Universe. The Euclid mission aims at performing an extensive spectroscopic survey of approximately 20–30 million Hα-emitting galaxies up to about redshift two. This ambitious project seeks to elucidate the nature of dark energy by mapping the 3-dimensional clustering of galaxies over a significant portion of the sky. This paper presents the methodology and software developed for estimating the 3-dimensional 2-point correlation function within the Euclid Science Ground Segment. The software is designed to overcome the significant challenges posed by the large and complex Euclid data set, which involves millions of galaxies. Key challenges include efficient pair counting, managing computational resources, and ensuring the accuracy of the correlation function estimation. The software leverages advanced algorithms, including kd-tree, octree, and linked-list data partitioning strategies, to optimise the pair-counting process. The implementation also includes parallel processing capabilities using shared-memory open multi-processing to further enhance performance and reduce computation times. Extensive validation and performance testing of the software are presented. The results indicate that the software is robust and can reliably estimate the 2-point correlation function, which is essential for deriving cosmological parameters with high precision. Furthermore, the paper discusses the expected performance of the software during different stages of the Euclid Wide Survey observations and forecasts how the precision of the correlation function measurements will improve over the mission's timeline, highlighting the software's capability to handle large data sets efficiently.

The Near-Infrared Spectrometer and Photometer (NISP) on board the satellite provides multiband photometry and $R slitless grism spectroscopy in the 950--2020\,nm wavelength range. In this reference article, we illuminate the background of NISP's functional and calibration requirements, describe the instrument's integral components, and provide all its key properties. We also sketch the processes needed to understand how NISP operates and is calibrated as well as its technical potentials and limitations. Links to articles providing more details and the technical background are included. The NISP's 16 H2RG detectors with a plate scale of $ deliver a field of view of 0.57\,deg$^2$. In photometric mode, NISP reaches a limiting magnitude of sim \,24.5\,AB\,mag in three photometric exposures of about 100\,s in exposure time for point sources and with a S/N of five. For spectroscopy, NISP's point-source sensitivity is a SNR = 3.5 detection of an emission line with flux sim \,$2 $ integrated over two resolution elements of 13.4\ in 3times 560\,s grism exposures at 1.6\ (redshifted Halpha ). Our calibration includes on-ground and in-flight characterisation and monitoring of the pixel-based detector baseline, dark current, non-linearity, and sensitivity to guarantee a relative photometric accuracy better than 1.5 and a relative spectrophotometry better than 0.7. The wavelength calibration must be accurate to 5\ or better. The NISP is the state-of-the-art instrument in the near-infrared for all science beyond small areas available from HST and JWST -- and it represents an enormous advance from any existing instrumentation due to its combination of field size and high throughput of telescope and instrument. During six-year survey covering 14\,000\,deg$^2$ of extragalactic sky, NISP will be the backbone in determining distances of more than a billion galaxies. Its near-infrared data will become a rich reference imaging and spectroscopy data set for the coming decades.

The near-infrared calibration unit (NI-CU) on board NISP is the first astronomical calibration lamp based on LED to be operated in space. is a mission in ESA's Cosmic Vision 2015--2025 framework to explore the dark universe and provide a next-level characterisation of the nature of gravitation, dark matter, and dark energy. Calibrating photometric and spectrometric measurements of galaxies to better than 1.5 accuracy in a survey homogeneously mapping sim \,14\,000\,deg$^2$ of extragalactic sky requires a very detailed characterisation of NIR detector properties as well as constant monitoring of them in flight. To cover two of the main contributions -- relative pixel-to-pixel sensitivity and non-linearity characteristics -- and to support other calibration activities, NI-CU was designed to provide spatially approximately homogeneous ($<$\,12 variations) and temporally stable illumination (0.1--0.2 over 1200\,s) over the NISP detector plane with minimal power consumption and energy dissipation. NI-CU covers the spectral range sim \,nm -- at cryo-operating temperature -- at five fixed independent wavelengths to capture wavelength-dependent behaviour of the detectors, with fluence over a dynamic range of gtrsim \,100 from sim $. For this functionality, NI-CU is based on LED . We describe the rationale behind the decision and design process, the challenges in sourcing the right LED and the qualification process and lessons learned. We also provide a description of the completed NI-CU, its capabilities, and performance as well as its limits. NI-CU has been integrated into NISP and the satellite, and since launch in July 2023, it has started supporting survey operations.

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance.

The statistical distribution and evolution of key properties (e.g. accretion rate, mass, or spin) of active galactic nuclei (AGN), remain an open debate in astrophysics. The ESA Euclid space mission, launched on July 1st 2023, promises a breakthrough in this field. We create detailed mock catalogues of AGN spectra, from the rest-frame near-infrared down to the ultraviolet, including emission lines, to simulate what Euclid will observe for both obscured (type 2) and unobscured (type 1) AGN. We concentrate on the red grisms of the NISP instrument, which will be used for the wide-field survey, opening a new window for spectroscopic AGN studies in the near-infrared. We quantify the efficiency in the redshift determination as well as in retrieving the emission line flux of the H$\alpha$+[NII] complex as Euclid is mainly focused on this emission line as it is expected to be the brightest one in the probed redshift range. Spectroscopic redshifts are measured for 83% of the simulated AGN in the interval where the H$\alpha$+[NII] is visible (0.892x10^{-16}$ erg s$^{-1}$ cm$^{-2}$, encompassing the peak of AGN activity at $z\simeq 1-1.5$) within the spectral coverage of the red grism. Outside this redshift range, the measurement efficiency decreases significantly. Overall, a spectroscopic redshift is correctly determined for ~90% of type 2 AGN down to an emission line flux of $3x10^{-16}$ erg s$^{-1}$ cm$^{-2}$, and for type 1 AGN down to $8.5x10^{-16}$ erg s$^{-1}$ cm$^{-2}$. Recovered black hole mass values show a small offset with respect to the input values ~10%, but the agreement is good overall. With such a high spectroscopic coverage at z<2, we will be able to measure AGN demography, scaling relations, and clustering from the epoch of the peak of AGN activity down to the present-day Universe for hundreds of thousand AGN with homogeneous spectroscopic information.

The European Space Agency's Euclid mission is one of the upcoming generation of large-scale cosmology surveys, which will map the large-scale structure in the Universe with unprecedented precision. The development and validation of the SGS pipeline requires state-of-the-art simulations with a high level of complexity and accuracy that include subtle instrumental features not accounted for previously as well as faster algorithms for the large-scale production of the expected Euclid data products. In this paper, we present the Euclid SGS simulation framework as applied in a large-scale end-to-end simulation exercise named Science Challenge 8. Our simulation pipeline enables the swift production of detailed image simulations for the construction and validation of the Euclid mission during its qualification phase and will serve as a reference throughout operations. Our end-to-end simulation framework starts with the production of a large cosmological N-body & mock galaxy catalogue simulation. We perform a selection of galaxies down to I_E=26 and 28 mag, respectively, for a Euclid Wide Survey spanning 165 deg^2 and a 1 deg^2 Euclid Deep Survey. We build realistic stellar density catalogues containing Milky Way-like stars down to H<26. Using the latest instrumental models for both the Euclid instruments and spacecraft as well as Euclid-like observing sequences, we emulate with high fidelity Euclid satellite imaging throughout the mission's lifetime. We present the SC8 data set consisting of overlapping visible and near-infrared Euclid Wide Survey and Euclid Deep Survey imaging and low-resolution spectroscopy along with ground-based. This extensive data set enables end-to-end testing of the entire ground segment data reduction and science analysis pipeline as well as the Euclid mission infrastructure, paving the way to future scientific and technical developments and enhancements.

Core ellipticals, which are massive early-type galaxies with almost constant inner surface brightness profiles, are the result of dry mergers. During these events, a binary black hole (BBH) is formed, destroying the original cuspy central regions of the merging objects and scattering stars that are not on tangential orbits. The size of the emerging core correlates with the mass of the finally merged black hole (BH). Therefore, the determination of the size of the core of massive early-type galaxies provides key insights not only into the mass of the black hole, but also into the origin and evolution of these objects. In this work, we report the first Euclid-based dynamical mass determination of a supermassive black hole (SMBH). To this end, we study the center of NGC 1272, the second most luminous elliptical galaxy in the Perseus cluster, combining the Euclid Visible Camera (VIS) photometry coming from the Early Release Observations (EROs) of the Perseus cluster with the Visible Integral-field Replicable Unit Spectrograph (VIRUS) spectroscopic observations at the Hobby-Eberly Telescope (HET). The core of NGC 1272 is detected on the Euclid VIS image. Its size is 1.​​″29 ± 0.​​″07 or 0.45 kpc, which was determined by fitting PSF-convolved core-Sérsic and Nuker-law functions. We deproject the surface brightness profile of the galaxy, finding that the galaxy is axisymmetric and nearly spherical. The two-dimensional stellar kinematics of the galaxy is measured from the VIRUS spectra by deriving optimally regularized non-parametric line-of-sight velocity distributions. Dynamical models of the galaxy are constructed using our axisymmetric and triaxial Schwarzschild codes. We measure a BH mass of (5 ± 3)×109 M⊙, which is in line with the expectation from the MBH − rb correlation, but is eight times larger than predicted by the MBH − σ correlation (at 1.8σ significance). The core size, rather than the velocity dispersion, allows one to select galaxies harboring the most massive BHs. The spatial resolution, wide area coverage, and depth of the Euclid (Wide and Deep) surveys allow us to find cores of passive galaxies that are larger than 2 kpc at a redshift of up to 1.

Core ellipticals, massive early-type galaxies have an almost constant inner surface brightness profile. The size of the core region correlates with the mass of the finally merged black hole. Here we report the first Euclid-based dynamical mass determination of a supermassive black hole. We study the centre of NGC 1272, the second most luminous elliptical galaxy in the Perseus cluster, combining the Euclid VIS photometry coming from the Early Release Observations of the Perseus cluster with VIRUS spectroscopic observations at the Hobby-Eberly Telescope. The core of NGC 1272 is detected on the Euclid VIS image. Its size is 1.29± 0.07” or 0.45 kpc, determined by fitting PSF-convolved core-Sérsic and Nuker-law functions. The two-dimensional stellar kinematics of the galaxy is measured from the VIRUS spectra by deriving optimally regularized non-parametric line-of-sight velocity distributions. Dynamical models of the galaxy are constructed using our axisymmetric and triaxial Schwarzschild codes. We measure a black hole mass of (5±3) × 10^9 M_⊙, in line with the expectation from the M_ BH-r_ b correlation, but eight times larger than predicted by the M_ BH-σ correlation (at 1.8σ significance). The core size, rather than the velocity dispersion, allows one to select galaxies harboring the most massive black holes. The spatial resolution, wide area coverage, and depth of the (Wide and Deep) surveys allow us to find cores of passive galaxies larger than 2 kpc up to redshift 1.

Core ellipticals, which are massive early-type galaxies with almost constant inner surface brightness profiles, are the result of dry mergers. During these events, a binary black hole (BBH) is formed, destroying the original cuspy central regions of the merging objects and scattering stars that are not on tangential orbits. The size of the emerging core correlates with the mass of the finally merged black hole (BH). Therefore, the determination of the size of the core of massive early-type galaxies provides key insights not only into the mass of the black hole, but also into the origin and evolution of these objects. In this work, we report the first Euclid-based dynamical mass determination of a supermassive black hole (SMBH). To this end, we study the center of NGC 1272, the second most luminous elliptical galaxy in the Perseus cluster, combining the Euclid Visible Camera (VIS) photometry coming from the Early Release Observations (EROs) of the Perseus cluster with the Visible Integral-field Replicable Unit Spectrograph (VIRUS) spectroscopic observations at the Hobby-Eberly Telescope (HET). The core of NGC 1272 is detected on the Euclid VIS image. Its size is 1.'' 29 +/- 0.'' 07 or 0.45 kpc, which was determined by fitting PSF-convolved core-Sersic and Nuker-law functions. We deproject the surface brightness profile of the galaxy, finding that the galaxy is axisymmetric and nearly spherical. The two-dimensional stellar kinematics of the galaxy is measured from the VIRUS spectra by deriving optimally regularized non-parametric line-of-sight velocity distributions. Dynamical models of the galaxy are constructed using our axisymmetric and triaxial Schwarzschild codes. We measure a BH mass of (5 +/- 3)x10(9) M-circle dot, which is in line with the expectation from the M-BH-r(b) correlation, but is eight times larger than predicted by the M-BH-sigma correlation (at 1.8 sigma significance). The core size, rather than the velocity dispersion, allows one to select galaxies harboring the most massive BHs. The spatial resolution, wide area coverage, and depth of the Euclid (Wide and Deep) surveys allow us to find cores of passive galaxies that are larger than 2 kpc at a redshift of up to 1.

Euclid is an ESA mission designed to constrain the properties of dark energy and gravity via weak gravitational lensing and galaxy clustering. It will carry out a wide area imaging and spectroscopy survey (EWS) in visible and near-infrared, covering roughly 15,000 square degrees of extragalactic sky on six years. The wide-field telescope and instruments are optimized for pristine PSF and reduced straylight, producing very crisp images. This paper presents the building of the Euclid reference survey: the sequence of pointings of EWS, Deep fields, Auxiliary fields for calibrations, and spacecraft movements followed by Euclid as it operates in a step-and-stare mode from its orbit around the Lagrange point L2. Each EWS pointing has four dithered frames; we simulate the dither pattern at pixel level to analyse the effective coverage. We use up-to-date models for the sky background to define the Euclid region-of-interest (RoI). The building of the reference survey is highly constrained from calibration cadences, spacecraft constraints and background levels; synergies with ground-based coverage are also considered. Via purposely-built software optimized to prioritize best sky areas, produce a compact coverage, and ensure thermal stability, we generate a schedule for the Auxiliary and Deep fields observations and schedule the RoI with EWS transit observations. The resulting reference survey RSD_2021A fulfills all constraints and is a good proxy for the final solution. Its wide survey covers 14,500 square degrees. The limiting AB magnitudes ($5\sigma$ point-like source) achieved in its footprint are estimated to be 26.2 (visible) and 24.5 (near-infrared); for spectroscopy, the H$_\alpha$ line flux limit is $2\times 10^{-16}$ erg cm$^{-2}$ s$^{-1}$ at 1600 nm; and for diffuse emission the surface brightness limits are 29.8 (visible) and 28.4 (near-infrared) mag arcsec$^{-2}$.

Context. The ESA Euclid mission will produce photometric galaxy samples over 15 000 square degrees of the sky that will be rich for clustering and weak lensing statistics. The accuracy of the cosmological constraints derived from these measurements will depend on the knowledge of the underlying redshift distributions based on photometric redshift calibrations. Aims. A new approach is proposed to use the stacked spectra from Euclid slitless spectroscopy to augment broad-band photometric information to constrain the redshift distribution with spectral energy distribution fitting. The high spectral resolution available in the stacked spectra complements the photometry and helps to break the colour-redshift degeneracy and constrain the redshift distribution of galaxy samples. Methods. We modelled the stacked spectra as a linear mixture of spectral templates. The mixture may be inverted to infer the underlying redshift distribution using constrained regression algorithms. We demonstrate the method on simulated Vera C. Rubin Observatory and Euclid mock survey data sets based on the Euclid Flagship mock galaxy catalogue. We assess the accuracy of the reconstruction by considering the inference of the baryon acoustic scale from angular two-point correlation function measurements. Results. We selected mock photometric galaxy samples at redshift z>1 using the self-organising map algorithm. Considering the idealised case without dust attenuation, we find that the redshift distributions of these samples can be recovered with 0.5% accuracy on the baryon acoustic scale. The estimates are not significantly degraded by the spectroscopic measurement noise due to the large sample size. However, the error degrades to 2% when the dust attenuation model is left free. We find that the colour degeneracies introduced by attenuation limit the accuracy considering the wavelength coverage of Euclid near-infrared spectroscopy.

Euclid is a major ESA mission scheduled for launch in 2023-2024 to map the geometry of the dark Universe using two primary probes, weak gravitational lensing and galaxy clustering. Euclid’s instruments, a visible imager (VIS) and an infrared spectrometer and photometer (NISP) have both been designed and built by Euclid Consortium teams. The NISP instrument will hold a large focal plane array of 16 near-infrared H2RG detectors, which are key elements to the performance of the NISP, and therefore to the science return of the mission.Euclid NISP H2RG flight detectors have been individually and thoroughly characterized at Centre de Physique des Particules de Marseille (CPPM) during a whole year with a view to producing a reference database of performance pixel maps. Analyses have been ongoing and have shown the relevance of taking into account spatial variations in deriving performance parameters. This paper will concentrate on interpixel capacitance (IPC) and conversion gain. First, per pixel IPC coefficient maps will be derived thanks to single pixel reset (SPR) measurements and a new IPC correction method will be defined and validated. Then, the paper will look into correlation effects of IPC and their impact on the derivation of per super-pixel IPC-free conversion gain maps. Eventually, several conversion gain values will be defined over clearly distinguishable regions.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

Context. Future weak lensing surveys, such as the Euclid mission, will attempt to measure the shapes of billions of galaxies in order to derive cosmological information. These surveys will attain very low levels of statistical error, and systematic errors must be extremely well controlled. In particular, the point spread function (PSF) must be estimated using stars in the field, and recovered with high accuracy.Aims. The aims of this paper are twofold. Firstly, we took steps toward a nonparametric method to address the issue of recovering the PSF field, namely that of finding the correct PSF at the position of any galaxy in the field, applicable to Euclid. Our approach relies solely on the data, as opposed to parametric methods that make use of our knowledge of the instrument. Secondly, we studied the impact of imperfect PSF models on the shape measurement of galaxies themselves, and whether common assumptions about this impact hold true in an Euclid scenario.Methods. We extended the recently proposed resolved components analysis approach, which performs super-resolution on a field of under-sampled observations of a spatially varying, image-valued function. We added a spatial interpolation component to the method, making it a true 2-dimensional PSF model. We compared our approach to PSFEx, then quantified the impact of PSF recovery errors on galaxy shape measurements through image simulations.Results. Our approach yields an improvement over PSFEx in terms of the PSF model and on observed galaxy shape errors, though it is at present far from reaching the required Euclid accuracy. We also find that the usual formalism used for the propagation of PSF model errors to weak lensing quantities no longer holds in the case of an Euclid-like PSF. In particular, different shape measurement approaches can react differently to the same PSF modeling errors.

Forthcoming large photometric surveys for cosmology require precise and accurate photometric redshift (photo-z) measurements for the success of their main science objectives. However, to date, no method has been able to produce photo-zs at the required accuracy using only the broad-band photometry that those surveys will provide. An assessment of the strengths and weaknesses of current methods is a crucial step in the eventual development of an approach to meet this challenge. We report on the performance of 13 photometric redshift code single value redshift estimates and redshift probability distributions (PDZs) on a common set of data, focusing particularly on the 0.2-2.6 redshift range that the Euclid mission will probe. We designed a challenge using emulated Euclid data drawn from three photometric surveys of the COSMOS field. The data was divided into two samples: one calibration sample for which photometry and redshifts were provided to the participants; and the validation sample, containing only the photometry to ensure a blinded test of the methods. Participants were invited to provide a redshift single value estimate and a PDZ for each source in the validation sample, along with a rejection flag that indicates the sources they consider unfit for use in cosmological analyses. The performance of each method was assessed through a set of informative metrics, using cross-matched spectroscopic and highly-accurate photometric redshifts as the ground truth. We show that the rejection criteria set by participants are efficient in removing strong outliers, that is to say sources for which the photo-z deviates by more than 0.15(1+z) from the spectroscopic-redshift (spec-z). We also show that, while all methods are able to provide reliable single value estimates, several machine-learning methods do not manage to produce useful PDZs. We find that no machine-learning method provides good results in the regions of galaxy color-space that are sparsely populated by spectroscopic-redshifts, for example z> 1. However they generally perform better than template-fitting methods at low redshift (z< 0.7), indicating that template-fitting methods do not use all of the information contained in the photometry. We introduce metrics that quantify both photo-z precision and completeness of the samples (post-rejection), since both contribute to the final figure of merit of the science goals of the survey (e.g., cosmic shear from Euclid). Template-fitting methods provide the best results in these metrics, but we show that a combination of template-fitting results and machine-learning results with rejection criteria can outperform any individual method. On this basis, we argue that further work in identifying how to best select between machine-learning and template-fitting approaches for each individual galaxy should be pursued as a priority.

We provide predictions of the yield of 7 < z < 9 quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; Phi) with redshift, Phi proportional to 10(k(z-6)), k = 0:72, and a further steeper rate of decline, k = 0:92; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we make use of an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to J(AB) similar to 23. Quasars at z > 8 may be selected from Euclid OYJH photometry alone, but selection over the redshift interval 7 < z < 8 is greatly improved by the addition of z-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond z = 6. If the decline of the QLF accelerates beyond z = 6, with k = 0.92, Euclid should nevertheless find over 100 quasars with 7.0 < z < 7.5, and similar to 25 quasars beyond the current record of z = 7.5, including similar to 8 beyond z = 8.0. The first Euclid quasars at z > 7.5 should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, 7 < z < 8, M-1450 < 25, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at J(AB) similar to 23. The precision with which k can be determined over 7 < z < 8 depends on the value of k, but assuming k = 0.72 it can be measured to a 1 sigma uncertainty of 0.07.

We provide predictions of the yield of $7 8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7 7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7

We present the XXL Survey, the largest XMM programme totaling some 6.9 Ms to date and involving an international consortium of roughly 100 members. The XXL Survey covers two extragalactic areas of 25 deg2 each at a point-source sensitivity of ~ 5E-15 erg/sec/cm2 in the [0.5-2] keV band (completeness limit). The survey's main goals are to provide constraints on the dark energy equation of state from the space-time distribution of clusters of galaxies and to serve as a pathfinder for future, wide-area X-ray missions. We review science objectives, including cluster studies, AGN evolution, and large-scale structure, that are being conducted with the support of approximately 30 follow-up programmes. We describe the 542 XMM observations along with the associated multi-lambda and numerical simulation programmes. We give a detailed account of the X-ray processing steps and describe innovative tools being developed for the cosmological analysis. The paper provides a thorough evaluation of the X-ray data, including quality controls, photon statistics, exposure and background maps, and sky coverage. Source catalogue construction and multi-lambda associations are briefly described. This material will be the basis for the calculation of the cluster and AGN selection functions, critical elements of the cosmological and science analyses. The XXL multi-lambda data set will have a unique lasting legacy value for cosmological and extragalactic studies and will serve as a calibration resource for future dark energy studies with clusters and other X-ray selected sources. With the present article, we release the XMM XXL photon and smoothed images along with the corresponding exposure maps. The XMM XXL observation list (Table B.1) is available in electronic form at the CDS. The present paper is the first in a series reporting results of the XXL-XMM survey.

We present a technique to measure lightcurves of time-variable point sources on a spatially structured background from imaging data. The technique was developed to measure light curves of SNLS supernovae in order to infer their distances. This photometry technique performs simultaneous PSF photometry at the same sky position on an image series. We describe two implementations of the method: one that resamples images before measuring fluxes, and one which does not. In both instances, we sketch the key algorithms involved and present the validation using semi-artificial sources introduced in real images in order to assess the accuracy of the supernova flux measurements relative to that of surrounding stars. We describe the methods required to anchor these PSF fluxes to calibrated aperture catalogs, in order to derive SN magnitudes. We find a marginally significant bias of 2 mmag of the after-resampling method, and no bias at the mmag accuracy for the non-resampling method. Given surrounding star magnitudes, we determine the systematic uncertainty of SN magnitudes to be less than 1.5 mmag, which represents about one third of the current photometric calibration uncertainty affecting SN measurements. The SN photometry delivers several by-products: bright star PSF flux mea- surements which have a repeatability of about 0.6%, as for aperture measurements; we measure relative astrometric positions with a noise floor of 2.4 mas for a single-image bright star measurement; we show that in all bands of the MegaCam instrument, stars exhibit a profile linearly broadening with flux by about 0.5% over the whole brightness range.