Characterization of the redshift distribution of ensembles of galaxies is pivotal for large scale structure cosmological studies. In this work, we focus on improving the Self-Organizing Map (SOM) methodology for photometric redshift estimation (SOMPZ), specifically in anticipation of the Dark Energy Survey Year 6 (DES Y6) data. This data set, featuring deeper and fainter galaxies than DES Year 3 (DES Y3), demands adapted techniques to ensure accurate recovery of the underlying redshift distribution. We investigate three strategies for enhancing the existing SOM-based approach used in DES Y3: 1) Replacing the Y3 SOM algorithm with one tailored for redshift estimation challenges; 2) Incorporating $\textit{g}$-band flux information to refine redshift estimates (i.e. using $\textit{griz}$ fluxes as opposed to only $\textit{riz}$); 3) Augmenting redshift data for galaxies where available. These methods are applied to DES Y3 data, and results are compared to the Y3 fiducial ones. Our analysis indicates significant improvements with the first two strategies, notably reducing the overlap between redshift bins. By combining strategies 1 and 2, we have successfully managed to reduce redshift bin overlap in DES Y3 by up to 66$\%$. Conversely, the third strategy, involving the addition of redshift data for selected galaxies as an additional feature in the method, yields inferior results and is abandoned. Our findings contribute to the advancement of weak lensing redshift characterization and lay the groundwork for better redshift characterization in DES Year 6 and future stage IV surveys, like the Rubin Observatory.

Modeling the intrinsic alignment (IA) of galaxies poses a challenge to weak lensing analyses. Using the Dark Energy Survey Year 3 shape catalog, we expect less impact from IA when we limit the sample to blue, star-forming galaxies. The cosmological parameter constraints from this BLUE cosmic shear sample are stable to IA model choice, unlike passive galaxies in the full DES Y3 sample, the goodness-of-fit is improved and the Omega(m) and S-8 better agree with the observations from Planck on the cosmic microwave background. Mitigating IA with sample selection in DES, rather than flexible model choices, can reduce uncertainty in S-8 by a factor of 1.5.

Using 1533 type Ia supernovae (SNe Ia) from the five-year sample of the Dark Energy Survey (DES), we investigate the relationship between the projected galactocentric separation of the SNe and their host galaxies and their light curves and standardization. We show, for the first time, that the difference in SN Ia post-standardization brightnesses between high and low-mass hosts reduces from 0.078±0.011 mag in the full sample to 0.036 ± 0.018 mag for SNe Ia located in the outer regions of their host galaxies, while increasing to 0.100 ± 0.014 mag for SNe in the inner regions. The difference in the size of the mass step between inner and outer regions is 0.064±0.023 mag. In these inner regions, the step can be reduced (but not removed) using a model where the R_V of dust along the line-of-sight to the SN changes as a function of galaxy properties. We investigate the remaining difference using the distributions of the SN Ia stretch parameter to test the inferred age of SN progenitors. Comparing red (older) environments only, outer regions have a higher proportion of high-stretch SNe and a more homogeneous stretch distribution. However, this effect cannot explain the reduction in significance of any Hubble residual step in outer regions. We conclude that the standardized distances of SNe Ia located in the outer regions of galaxies are less affected by their global host galaxy properties than those in the inner regions.

The present state of cosmology is facing a crisis where there is a fundamental disagreement in measurements of the Hubble constant (H-0), with significant tension between the early and late Universe methods. Type Ia supernovae (SNe Ia) are important to measuring H-0 through the astronomical distance ladder. However, there remains potential to better standardize SN Ia light curves by using known dependencies on host galaxy properties after the standard light curve width and colour corrections have been applied to the peak SN Ia luminosities. To explore this, we use the 5-yr photometrically identified SNe Ia sample obtained by the Dark Energy Survey, along with host galaxy spectra obtained by the Australian Dark Energy Survey. Using host galaxy spectroscopy, we find a significant trend with the equivalent width (EW) of the [O II] lambda lambda 3727, 29 doublet, a proxy for specific star formation rate, and Hubble residuals. We find that the correlation with [O II] EW is a powerful alternative to the commonly used mass step after initial light-curve corrections. Applying this [O II] EW correction to 20 SNe Ia in calibrator galaxies observed with WiFeS, we examined the impact on SN Ia absolute magnitudes and H-0. Our [O II] EW corrections result in H-0 values ranging between 73.04 and 73.51 kms(-1)Mpc(-1), with a combined statistical and systematic uncertainty of similar to 1.31kms(-1)Mpc(-1). However, even with this additional correction, the impact of host galaxy properties in standardizing SNe Ia appears limited in reducing the current tension (similar to 5 sigma) with the cosmic microwave background result for H-0.

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)

We measure the current expansion rate of the Universe, Hubble's constant H-0, by calibrating the absolute magnitudes of supernovae to distances measured by baryon acoustic oscillations (BAO). This 'inverse distance ladder' technique provides an alternative to calibrating supernovae using nearby absolute distance measurements, replacing the calibration with a high-redshift anchor. We use the recent release of 1829 supernovae from the Dark Energy Survey spanning 0.01 < z < 1.13 anchored to the recent baryon acoustic oscillation measurements from Dark Energy Spectroscopic Instrument (DESI) spanning 0.30 < z(eff) < 2.33. To trace cosmology to z = 0, we use the third-, fourth-, and fifth-order cosmographic models, which, by design, are agnostic about the energy content and expansion history of the universe. With the inclusion of the higher redshift DESI-BAO data, the third-order model is a poor fit to both data sets, with the fourth-order model being preferred by the Akaike Information Criterion. Using the fourth-order cosmographic model, we find H-0 = 67. 19(-0.64)(+0.66) km s(-1) Mpc(-1), in agreement with the value found by Planck without the need to assume Flat-Lambda CDM. However, the best-fitting expansion history differs from that of Planck, providing continued motivation to investigate these tensions.

Context. The Euclid mission of the European Space Agency will deliver weak gravitational lensing and galaxy clustering surveys that can be used to constrain the standard cosmological model and extensions thereof. Aims. We present forecasts from the combination of the Euclid photometric galaxy surveys (weak lensing, galaxy clustering, and their crosscorrelations) and its spectroscopic redshift survey with respect to their sensitivity to cosmological parameters. We include the summed neutrino mass, Sigma m (v), and the e ffective number of relativistic species, N-e ff, in the standard Lambda alpha CDM scenario and in the dynamical dark energy (w (0) w(alpha)CDM) scenario. Methods. We compared the accuracy of di fferent algorithms predicting the non-linear matter power spectrum for such models. We then validated several pipelines for Fisher matrix and Markov chain Monte Carlo (MCMC) forecasts, using di fferent theory codes, algorithms for numerical derivatives, and assumptions on the non-linear cut-o ff scale. Results. The Euclid primary probes alone will reach a sensitivity of sigma(Sigma m (v) = 60 meV) = 56 meV in the Lambda CDM +Sigma m (v) model, whereas the combination with cosmic microwave background (CMB) data from Planck is expected to achieve sigma(Sigma m (v)) = 23 meV, o ffering evidence of a non-zero neutrino mass to at least the 2:6 sigma level. This could be pushed to a 4 sigma detection if future CMB data from LiteBIRD and CMB Stage-IV were included. In combination with Planck, Euclid will also deliver tight constraints on Delta N-e ff < 0:144 (95%CL) in the Lambda CDM +Sigma m (v)+N-e ff model or even Delta N-e ff < 0:063 when future CMB data are included. When floating the dark energy parameters, we find that the sensitivity to Ne ff remains stable, but for Sigma m (v), it gets degraded by up to a factor of 2, at most. Conclusions. This work illustrates the complementarity among the Euclid spectroscopic and photometric surveys and among Euclid and CMB constraints. Euclid will o ffer great potential in measuring the neutrino mass and excluding well-motivated scenarios with additional relativistic particles.

We present a study of the weak lensing inferred matter profiles Delta Sigma(R) of 698 South Pole Telescope (SPT) thermal Sunyaev-Zel'dovich effect (tSZE) selected and MCMF optically confirmed galaxy clusters in the redshift range 0.25 < z < 0.94 that have associated weak gravitational lensing shear profiles from the Dark Energy Survey (DES). Rescaling these profiles to account for the mass dependent size and the redshift dependent density produces average rescaled matter profiles Delta Sigma(R/R-200c)/(rho R-crit(200c)) with a lower dispersion than the unscaled Delta Sigma(R) versions, indicating a significant degree of self-similarity. Galaxy clusters from hydrodynamical simulations also exhibit matter profiles that suggest a high degree of self-similarity, with RMS variation among the average rescaled matter profiles with redshift and mass falling by a factor of approximately six and 23, respectively, compared to the unscaled average matter profiles. We employed this regularity in a new Bayesian method for weak lensing mass calibration that employs the so-called cluster mass posterior P(M-200|zeta, lambda, z), which describes the individual cluster masses given their tSZE (zeta) and optical (lambda, z) observables. This method enables simultaneous constraints on richness lambda-mass and tSZE detection significance zeta-mass relations using average rescaled cluster matter profiles. We validated the method using realistic mock datasets and present observable-mass relation constraints for the SPTxDES sample, where we constrained the amplitude, mass trend, redshift trend, and intrinsic scatter. Our observable-mass relation results are in agreement with the mass calibration derived from the recent cosmological analysis of the SPTxDES data based on a cluster-by-cluster lensing calibration. Our new mass calibration technique offers a higher efficiency when compared to the single cluster calibration technique. We present new validation tests of the observable-mass relation that indicate the underlying power-law form and scatter are adequate to describe the real cluster sample but that also suggest a redshift variation in the intrinsic scatter of the lambda-mass relation may offer a better description. In addition, the average rescaled matter profiles offer high signal-to-noise ratio (S/N) constraints on the shape of real cluster matter profiles, which are in good agreement with available hydrodynamical Lambda CDM simulations. This high S/N profile contains information about baryon feedback, the collisional nature of dark matter, and potential deviations from general relativity.

We study the constraint on f(R) gravity that can be obtained by photometric primary probes of the Euclid mission. Our focus is the dependence of the constraint on the theoretical modelling of the nonlinear matter power spectrum. In the Hu-Sawicki f(R) gravity model, we consider four different predictions for the ratio between the power spectrum in f(R) and that in Lambda cold dark matter (Lambda CDM): a fitting formula, the halo model reaction approach, ReACT, and two emulators based on dark matter only N-body simulations, FORGE and e-Mantis. These predictions are added to the MontePython implementation to predict the angular power spectra for weak lensing (WL), photometric galaxy clustering, and their cross-correlation. By running Markov chain Monte Carlo, we compare constraints on parameters and investigate the bias of the recovered f(R) parameter if the data are created by a different model. For the pessimistic setting of WL, one-dimensional bias for the f(R) parameter, log10|fR0|, is found to be 0.5 sigma when FORGE is used to create the synthetic data with log10|fR0| = -5.301 and fitted by e-Mantis. The impact of baryonic physics on WL is studied by using a baryonification emulator, BCemu. For the optimistic setting, the f(R) parameter and two main baryonic parameters are well constrained despite the degeneracies among these parameters. However, the difference in the nonlinear dark matter prediction can be compensated for the adjustment of baryonic parameters, and the one-dimensional marginalised constraint on log10|fR0| is biased. This bias can be avoided in the pessimistic setting at the expense of weaker constraints. For the pessimistic setting, using the Lambda CDM synthetic data for WL, we obtain the prior-independent upper limit of log10|fR0| < -5.6. Finally, we implement a method to include theoretical errors to avoid the bias due to inaccuracies in the nonlinear matter power spectrum prediction.

We forecast the expected population of active galactic nuclei (AGN) observable in the Euclid Wide Survey (EWS) and Euclid Deep Survey (EDS). Starting from an X-ray luminosity function (XLF) we generate volume-limited samples of the AGN expected in the survey footprints. Each AGN is assigned an SED appropriate for its X-ray luminosity and redshift, with perturbations sampled from empirical distributions. The photometric detectability of each AGN is assessed via mock observation of the assigned SED. We estimate 40 million AGN will be detectable in at least one band in the EWS and 0.24 million in the EDS, corresponding to surface densities of 2.8$\times$10$^{3}$ deg$^{-2}$ and 4.7$\times$10$^{3}$ deg$^{-2}$. Employing colour selection criteria on our simulated data we select a sample of 4.8$\times$10$^{6}$ (331 deg$^{-2}$) AGN in the EWS and 1.7$\times$10$^{4}$ (346 deg$^{-2}$) in the EDS, amounting to 10% and 8% of the AGN detectable in the EWS and EDS. Including ancillary Rubin/LSST bands improves the completeness and purity of AGN selection. These data roughly double the total number of selected AGN to comprise 21% and 15% of the detectable AGN in the EWS and EDS. The total expected sample of colour-selected AGN contains 6.0$\times$10$^{6}$ (74%) unobscured AGN and 2.1$\times$10$^{6}$ (26%) obscured AGN, covering $0.02 \leq z \lesssim 5.2$ and $43 \leq \log_{10} (L_{bol} / erg s^{-1}) \leq 47$. With this simple colour selection, expected surface densities are already comparable to the yield of modern X-ray and mid-infrared surveys of similar area. The relative uncertainty on our expectation for detectable AGN is 6.7% for the EWS and 12.5% for the EDS, driven by the uncertainty of the XLF.

Extragalactic globular clusters (EGCs) are an abundant and powerful tracer of galaxy dynamics and formation, and their own formation and evolution is also a matter of extensive debate. The compact nature of globular clusters means that they are hard to spatially resolve and thus study outside the Local Group. In this work we have examined how well EGCs will be detectable in images from the Euclid telescope, using both simulated pre-launch images and the first early-release observations of the Fornax galaxy cluster. The Euclid Wide Survey will provide high-spatial resolution VIS imaging in the broad IE band as well as near-infrared photometry (YE, JE, and HE). We estimate that the galaxies within 100 Mpc in the footprint of the Euclid survey host around 830 000 EGCs of which about 350 000 are within the survey's detection limits. For about half of these EGCs, three infrared colours will be available as well. For any galaxy within 50Mpc the brighter half of its GC luminosity function will be detectable by the Euclid Wide Survey. The detectability of EGCs is mainly driven by the residual surface brightness of their host galaxy. We find that an automated machine-learning EGC-classification method based on real Euclid data of the Fornax galaxy cluster provides an efficient method to generate high purity and high completeness GC candidate catalogues. We confirm that EGCs are spatially resolved compared to pure point sources in VIS images of Fornax. Our analysis of both simulated and first on-sky data show that Euclid will increase the number of GCs accessible with high-resolution imaging substantially compared to previous surveys, and will permit the study of GCs in the outskirts of their hosts. Euclid is unique in enabling systematic studies of EGCs in a spatially unbiased and homogeneous manner and is primed to improve our understanding of many understudied aspects of GC astrophysics.

The number density , redshift evolution of optically selected galaxy clusters offer an independent measurement of the amplitude of matter fluctuations, S 8 . However, recent results have shown that clusters chosen by the redMaPPer algorithm show richness-dependent biases that affect the weak lensing signals and number densities of clusters, increasing uncertainty in the cluster mass calibration and reducing their constraining power. In this work, we evaluate an alternative cluster proxy, outskirt stellar mass, M out , defined as the total stellar mass within a [50, 100] kpc envelope centered on a massive galaxy. This proxy exhibits scatter comparable to redMaPPer richness, A, but is less likely to be subject to projection effects. We compare the Dark Energy Survey Year 3 redMaPPer cluster catalog with a M out selected cluster sample from the Hyper-Suprime Camera survey. We use weak lensing measurements to quantify and compare the scatter of M out and A with halo mass. Our results show M out has a scatter consistent with A, with a similar halo mass dependence , that both proxies contain unique information about the underlying halo mass. We find A-selected samples introduce features into the measured Delta Sigma signal that are not well fit by a log- normal scatter only model, absent in M out selected samples. Our findings suggest that M out offers an alternative for cluster selection with more easily calibrated selection biases, at least at the generally lower richnesses probed here. Combining both proxies may yield a mass proxy with a lower scatter and more tractable selection biases, enabling the use of lower mass clusters in cosmology. Finally, we find the scatter and slope in the A - M out scaling relation to be 0.49 + 0.02 and 0.38 + 0.09.

Context. Cluster cosmology can benefit from combining multi-wavelength studies. In turn, these studies benefit from a characterisation of the correlation coefficients among different mass-observable relations. Aims. In this work, we aim to provide information on the scatter, skewness, and covariance of various mass-observable relations in galaxy clusters in cosmological hydrodynamic simulations. This information will help future analyses improve the general approach to accretion histories and projection effects, as well as to model mass-observable relations for cosmology studies. Methods. We identified galaxy clusters in Magneticum Box2b simulations with masses of M-200c > 10(14) M-circle dot at redshifts of z = 0.24 and z = 0.90. Our analysis included Euclid-derived properties such as richness, stellar mass, lensing mass, and concentration. Additionally, we investigated complementary multi-wavelength data, including X-ray luminosity, integrated Compton-y parameter, gas mass, and temperature. We then examined the impact of projection effects on mass-observable residuals and correlations. Results. We find that at intermediate redshift (z = 0.24), projection effects have the greatest impact of lensing concentration, richness, and gas mass in terms of the scatter and skewness of the log-residuals of scaling relations. The contribution of projection effects can be significant enough to boost a spurious hot- versus cold-baryon correlations and consequently hide underlying correlations due to halo accretion histories. At high redshift (z = 0.9), the richness has a much lower scatter (of log-residuals), while the quantity that is most impacted by projection effects is the lensing mass. The lensing concentration reconstruction, in particular, is affected by deviations of the reduced-shear profile shape from that derived using a Navarro-Frenk-White (NFW) profile; the amount of interlopers in the line of sight, on the other hand, is not as important.

We present the results of a search for gravitational-wave transients associated with core-collapse supernova SN 2023ixf, which was observed in the galaxy Messier 101 via optical emission on 2023 May 19, during the LIGO–Virgo–KAGRA 15th Engineering Run. We define a five-day on-source window during which an accompanying gravitational-wave signal may have occurred. No gravitational waves have been identified in data when at least two gravitational-wave observatories were operating, which covered ∼14% of this five-day window. We report the search detection efficiency for various possible gravitational-wave emission models. Considering the distance to M101 (6.7 Mpc), we derive constraints on the gravitational-wave emission mechanism of core-collapse supernovae across a broad frequency spectrum, ranging from 50 Hz to 2 kHz, where we assume the gravitational-wave emission occurred when coincident data are available in the on-source window. Considering an ellipsoid model for a rotating proto-neutron star, our search is sensitive to gravitational-wave energy 1 × 10 −4 M ⊙ c 2 and luminosity 2.6 × 10 −4 M ⊙ c 2 s −1 for a source emitting at 82 Hz. These constraints are around an order of magnitude more stringent than those obtained so far with gravitational-wave data. The constraint on the ellipticity of the proto-neutron star that is formed is as low as 1.08, at frequencies above 1200 Hz, surpassing past results.

The Euclid mission of the European Space Agency will provide weak gravitational lensing and galaxy clustering surveys that can be used to constrain the standard cosmological model and its extensions, with an opportunity to test the properties of dark matter beyond the minimal cold dark matter paradigm. We present forecasts from the combination of the Euclid weak lensing and photometric galaxy clustering data on the parameters describing four interesting and representative non-minimal dark matter models: a mixture of cold and warm dark matter relics; unstable dark matter decaying either into massless or massive relics; and dark matter undergoing feeble interactions with relativistic relics. We modelled these scenarios at the level of the non-linear matter power spectrum using emulators trained on dedicated N-body simulations. We used a mock Euclid likelihood and Monte Carlo Markov chains to fit mock data and infer error bars on dark matter parameters marginalised over other parameters. We find that the Euclid photometric probe (alone or in combination with cosmic microwave background data from the Planck satellite) will be sensitive to the effect of each of the four dark matter models considered here. The improvement will be particularly spectacular for decaying and interacting dark matter models. With Euclid, the bounds on some dark matter parameters can improve by up to two orders of magnitude compared to current limits. We discuss the dependence of predicted uncertainties on different assumptions: the inclusion of photometric galaxy clustering data, the minimum angular scale taken into account, and modelling of baryonic feedback effects. We conclude that the Euclid mission will be able to measure quantities related to the dark sector of particle physics with unprecedented sensitivity. This will provide important information for model building in high-energy physics. Any hint of a deviation from the minimal cold dark matter paradigm would have profound implications for cosmology and particle physics.

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.

Surveys of cosmological large-scale structure (LSS) are sensitive to the presence of local primordial non-Gaussianity (PNG), and may be used to constrain models of inflation. Local PNG, characterized by fNL, the amplitude of the quadratic correction to the potential of a Gaussian random field, is traditionally measured from LSS two-point and three-point clustering via the power spectrum and bi-spectrum. We propose a framework to measure fNL using the configuration space two-point correlation function (2pcf) monopole and three-point correlation function (3pcf) monopole of survey tracers. Our model estimates the effect of the scale-dependent bias induced by the presence of PNG on the 2pcf and 3pcf from the clustering of simulated dark matter halos. We describe how this effect may be scaled to an arbitrary tracer of the cosmological matter density. The 2pcf and 3pcf of this tracer are measured to constrain the value of fNL. Using simulations of luminous red galaxies observed by the Dark Energy Spectroscopic Instrument (DESI), we demonstrate the accuracy and constraining power of our model, and forecast the ability to constrainfNL to a precision of sigma(fNL) = 22 with one year of DESI survey data.

To date, galaxy image simulations for weak lensing surveys usually approximate the light profiles of all galaxies as a single or double Sérsic profile, neglecting the influence of galaxy substructures and morphologies deviating from such a simplified parametric characterisation. While this approximation may be sufficient for previous data sets, the stringent cosmic shear calibration requirements and the high quality of the data in the upcoming survey demand a consideration of the effects that realistic galaxy substructures and irregular shapes have on shear measurement biases. Here we present a novel deep learning-based method to create such simulated galaxies directly from Hubble Space Telescope (HST) data. We first build and validate a convolutional neural network based on the wavelet scattering transform to learn noise-free representations independent of the point-spread function (PSF) of HST galaxy images. These can be injected into simulations of images from optical instrument VIS without introducing noise correlations during PSF convolution or shearing. Then, we demonstrate the generation of new galaxy images by sampling from the model randomly as well as conditionally. In the latter case, we fine-tune the interpolation between latent space vectors of sample galaxies to directly obtain new realistic objects following a specific Sérsic index and half-light radius distribution. Furthermore, we show that the distribution of galaxy structural and morphological parameters of our generative model matches the distribution of the input HST training data, proving the capability of the model to produce realistic shapes. Next, we quantify the cosmic shear bias from complex galaxy shapes in simulations by comparing the shear measurement biases between a sample of model objects and their best-fit double-Sérsic counterparts, thereby creating two separate branches that only differ in the complexity of their shapes. Using the Kaiser, Squires, and Broadhurst shape measurement algorithm, we find a multiplicative bias difference between these branches with realistic morphologies and parametric profiles on the order of $(6.9± 0.6) $ for a realistic magnitude-Sérsic index distribution. Moreover, we find clear detection bias differences between full image scenes simulated with parametric and realistic galaxies, leading to a bias difference of $(4.0± 0.9) $ independent of the shape measurement method. This makes complex morphology relevant for stage IV weak lensing surveys, exceeding the full error budget of the Euclid Wide Survey (Δμ_