Showing 1–50 of 273 results for author: Teyssier, R

Euclid preparation. Forecasting the recovery of galaxy physical properties and their relations with template-fitting and machine-learning methods

Authors: Euclid Collaboration, A. Enia, M. Bolzonella, L. Pozzetti, A. Humphrey, P. A. C. Cunha, W. G. Hartley, F. Dubath, S. Paltani, X. Lopez Lopez, S. Quai, S. Bardelli, L. Bisigello, S. Cavuoti, G. De Lucia, M. Ginolfi, A. Grazian, M. Siudek, C. Tortora, G. Zamorani, N. Aghanim, B. Altieri, A. Amara, S. Andreon, N. Auricchio , et al. (238 additional authors not shown)

Abstract: Euclid will collect an enormous amount of data during the mission's lifetime, observing billions of galaxies in the extragalactic sky. Along with traditional template-fitting methods, numerous Machine Learning algorithms have been presented for computing their photometric redshifts and physical parameters (PP), requiring significantly less computing effort while producing equivalent performance me… ▽ More Euclid will collect an enormous amount of data during the mission's lifetime, observing billions of galaxies in the extragalactic sky. Along with traditional template-fitting methods, numerous Machine Learning algorithms have been presented for computing their photometric redshifts and physical parameters (PP), requiring significantly less computing effort while producing equivalent performance measures. However, their performance is limited by the quality and amount of input information, to the point where the recovery of some well-established physical relationships between parameters might not be guaranteed. To forecast the reliability of Euclid photo-$z$s and PPs calculations, we produced two mock catalogs simulating Euclid photometry. We simulated the Euclid Wide Survey (EWS) and Euclid Deep Fields (EDF). We tested the performance of a template-fitting algorithm (Phosphoros) and four ML methods in recovering photo-$z$s, stellar masses, star-formation rates, and the SFMS. To mimic the Euclid processing as closely as possible, the models were trained with Phosphoros-recovered labels. For the EWS, we found that the best results are achieved with a Mixed Labels approach, training the models with Wide survey features and labels from the Phosphoros results on deeper photometry, i.e., with the best possible set of labels for a given photometry. This imposes a prior, helping the models to better discern cases in degenerate regions of feature space, i.e., when galaxies have similar magnitudes and colors but different redshifts and PPs, with performance metrics even better than those found with Phosphoros. We found no more than $3\%$ performance degradation using a COSMOS-like reference sample or removing $u$ band data, which will not be available until after data release DR1. The best results are obtained for the EDF, with appropriate recovery of photo-$z$, PPs, and the SFMS. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 26 pages, 13 figures. Submitted to A&A

Euclid preparation. Sensitivity to non-standard particle dark matter model

Authors: Euclid Collaboration, J. Lesgourgues, J. Schwagereit, J. Bucko, G. Parimbelli, S. K. Giri, F. Hervas-Peters, A. Schneider, M. Archidiacono, F. Pace, Z. Sakr, A. Amara, L. Amendola, S. Andreon, N. Auricchio, H. Aussel, C. Baccigalupi, M. Baldi, S. Bardelli, R. Bender, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia, J. Brinchmann , et al. (227 additional authors not shown)

Abstract: 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 these surveys on the parameters describing four int… ▽ More 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 these surveys 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 experiencing feeble interactions with relativistic relics. We model these scenarios at the level of the non-linear matter power spectrum using emulators trained on dedicated N-body simulations. We use a mock Euclid likelihood 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 CMB 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: inclusion of photometric galaxy clustering data, minimum angular scale taken into account, 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. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 31 pages, 21 figures

Report number: TTK-24-26

Nexus: A framework for simulations of idealised galaxies

Authors: Thor Tepper-Garcia, Joss Bland-Hawthorn, Eugene Vasiliev, Oscar Agertz, Romain Teyssier, Christoph Federrath

Abstract: Motivated by the need for realistic, dynamically self-consistent, evolving galaxy models that avoid the inherent complexity of full, and zoom-in, cosmological simulations, we have developed Nexus, an integral, flexible framework to create synthetic galaxies made of both collisionless and gaseous components. Nexus leverages the power of publicly available, tried-and-tested packages: i) the stellar-… ▽ More Motivated by the need for realistic, dynamically self-consistent, evolving galaxy models that avoid the inherent complexity of full, and zoom-in, cosmological simulations, we have developed Nexus, an integral, flexible framework to create synthetic galaxies made of both collisionless and gaseous components. Nexus leverages the power of publicly available, tried-and-tested packages: i) the stellar-dynamics, action-based library AGAMA; and ii) the Adaptive Mesh Refinement, N-body/hydrodynamical code Ramses, modified to meet our needs. In addition, we make use of a proprietary module to account for realistic galaxy formation (sub-grid) physics, including star formation, stellar feedback, and chemical enrichment. As a framework to perform controlled experiments of idealised galaxies, Nexus' basic functionality consists in the generation of bespoke initial conditions (ICs) for any desired galaxy model, which are advanced in time to simulate the system's evolution. The fully self-consistent ICs are generated with a distribution-function based approach, as implemented in the galaxy modelling module of AGAMA -- up to now restricted to collisionless components, extended in this work to treat two types of gaseous configurations: (i) hot halos; and (ii) gas disks. For the first time, we are able to construct equilibrium models with disc gas fractions in the range $0 < f_{\rm gas} < 1$, needed to model high-redshift galaxies. The framework is ideally suited to the study of galactic ecology, specifically how stars and gas work together over billions of years. As a validation of our framework, we reproduce - and improve upon - several isolated galaxy model setups reported in earlier studies. Finally, we showcase Nexus by presenting an interesting type of `nested bar' galaxy class. Future upgrades of Nexus will include magneto-hydrodynamics and highly energetic particle (`cosmic ray') heating. △ Less

Submitted 1 June, 2024; originally announced June 2024.

Comments: To be submitted to MNRAS; we welcome general comments and pointers to relevant, missing references

Dust dynamics in RAMSES -- II. Equilibrium drift velocity distributions of charged dust grains

Authors: Eric R. Moseley, Romain Teyssier

Abstract: We investigate the gas-grain relative drift velocity distributions of charged astrophysical dust grains in MHD turbulence. We do this using a range of MHD-PIC simulations spanning different plasma-$β$, sonic/Alfvén Mach number, and with grains of varying size and charge-to-mass ratio. We find that the root-mean-square drift velocity is a strong function of the grain size, following a power law wit… ▽ More We investigate the gas-grain relative drift velocity distributions of charged astrophysical dust grains in MHD turbulence. We do this using a range of MHD-PIC simulations spanning different plasma-$β$, sonic/Alfvén Mach number, and with grains of varying size and charge-to-mass ratio. We find that the root-mean-square drift velocity is a strong function of the grain size, following a power law with a 1/2 slope. The r.m.s. value has only a very weak dependence on the charge-to-mass ratio. On the other hand, the shape of the distribution is a strong function of the grain charge-to-mass ratio, and in compressible turbulence, also the grain size. We then compare these results to simple analytic models based upon time-domain quasi-linear theory and solutions to the Fokker-Planck equation. These models explain qualitatively the r.m.s. drift velocity's lack of charge-to-mass ratio dependence, as well as why the shape of the distribution changes as the charge-to-mass ratio increases. Finally we scale our results to astrophysical conditions. As an example, at a length scale of one parsec in the cold neutral medium, 0.1 $μ$m grains should be drifting at roughly 40% of the turbulent velocity dispersion. These findings may serve as a basis for a model for grain velocities in the context of grain-grain collisions, non-thermal sputtering, and accretion of metals. These findings also have implications for the transport of grains through the galaxy, suggesting that grains may have non-negligible random motions at length-scales that many modern galaxy simulations approach. △ Less

Submitted 30 May, 2024; v1 submitted 28 May, 2024; originally announced May 2024.

Comments: 21 pages, 12 figures, submitted to MNRAS

Euclid preparation. Observational expectations for redshift z<7 active galactic nuclei in the Euclid Wide and Deep surveys

Authors: Euclid Collaboration, M. Selwood, S. Fotopoulou, M. N. Bremer, L. Bisigello, H. Landt, E. Bañados, G. Zamorani, F. Shankar, D. Stern, E. Lusso, L. Spinoglio, V. Allevato, F. Ricci, A. Feltre, F. Mannucci, M. Salvato, R. A. A. Bowler, M. Mignoli, D. Vergani, F. La Franca, A. Amara, S. Andreon, N. Auricchio, M. Baldi , et al. (238 additional authors not shown)

Abstract: 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 distribu… ▽ More 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. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: 36 pages, 21 figures, submitted to A&A

Euclid preparation. Detecting globular clusters in the Euclid survey

Authors: Euclid Collaboration, K. Voggel, A. Lançon, T. Saifollahi, S. S. Larsen, M. Cantiello, M. Rejkuba, J. -C. Cuillandre, P. Hudelot, A. A. Nucita, M. Urbano, E. Romelli, M. A. Raj, M. Schirmer, C. Tortora, Abdurro'uf, F. Annibali, M. Baes, P. Boldrini, R. Cabanac, D. Carollo, C. J. Conselice, P. -A. Duc, A. M. N. Ferguson, L. K. Hunt , et al. (247 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 29 May, 2024; v1 submitted 22 May, 2024; originally announced May 2024.

Comments: Submitted to A&A

Euclid. V. The Flagship galaxy mock catalogue: a comprehensive simulation for the Euclid mission

Authors: Euclid Collaboration, F. J. Castander, P. Fosalba, J. Stadel, D. Potter, J. Carretero, P. Tallada-Crespí, L. Pozzetti, M. Bolzonella, G. A. Mamon, L. Blot, K. Hoffmann, M. Huertas-Company, P. Monaco, E. J. Gonzalez, G. De Lucia, C. Scarlata, M. -A. Breton, L. Linke, C. Viglione, S. -S. Li, Z. Zhai, Z. Baghkhani, K. Pardede, C. Neissner , et al. (344 additional authors not shown)

Abstract: 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… ▽ More 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) △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid. IV. The NISP Calibration Unit

Authors: Euclid Collaboration, F. Hormuth, K. Jahnke, M. Schirmer, C. G. -Y. Lee, T. Scott, R. Barbier, S. Ferriol, W. Gillard, F. Grupp, R. Holmes, W. Holmes, B. Kubik, J. Macias-Perez, M. Laurent, J. Marpaud, M. Marton, E. Medinaceli, G. Morgante, R. Toledo-Moreo, M. Trifoglio, Hans-Walter Rix, A. Secroun, M. Seiffert, P. Stassi , et al. (310 additional authors not shown)

Abstract: The near-infrared calibration unit (NI-CU) on board Euclid's Near-Infrared Spectrometer and Photometer (NISP) is the first astronomical calibration lamp based on light-emitting diodes (LEDs) to be operated in space. Euclid 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 da… ▽ More The near-infrared calibration unit (NI-CU) on board Euclid's Near-Infrared Spectrometer and Photometer (NISP) is the first astronomical calibration lamp based on light-emitting diodes (LEDs) to be operated in space. Euclid 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 ~14000 deg^2 of extragalactic sky requires a very detailed characterisation of near-infrared (NIR) detector properties, as well their constant monitoring in flight. To cover two of the main contributions - relative pixel-to-pixel sensitivity and non-linearity characteristics - as well as support other calibration activities, NI-CU was designed to provide spatially approximately homogeneous (<12% variations) and temporally stable illumination (0.1%-0.2% over 1200s) over the NISP detector plane, with minimal power consumption and energy dissipation. NI-CU is covers the spectral range ~[900,1900] nm - at cryo-operating temperature - at 5 fixed independent wavelengths to capture wavelength-dependent behaviour of the detectors, with fluence over a dynamic range of >=100 from ~15 ph s^-1 pixel^-1 to >1500 ph s^-1 pixel^-1. For this functionality, NI-CU is based on LEDs. We describe the rationale behind the decision and design process, describe the challenges in sourcing the right LEDs, as well as 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 Euclid satellite, and since Euclid's launch in July 2023 has started supporting survey operations. △ Less

Submitted 10 July, 2024; v1 submitted 22 May, 2024; originally announced May 2024.

Comments: Paper accepted for publication in A&A as part of the special issue 'Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid. III. The NISP Instrument

Authors: Euclid Collaboration, K. Jahnke, W. Gillard, M. Schirmer, A. Ealet, T. Maciaszek, E. Prieto, R. Barbier, C. Bonoli, L. Corcione, S. Dusini, F. Grupp, F. Hormuth, S. Ligori, L. Martin, G. Morgante, C. Padilla, R. Toledo-Moreo, M. Trifoglio, L. Valenziano, R. Bender, F. J. Castander, B. Garilli, P. B. Lilje, H. -W. Rix , et al. (412 additional authors not shown)

Abstract: The Near-Infrared Spectrometer and Photometer (NISP) on board the Euclid satellite provides multiband photometry and R>=450 slitless grism spectroscopy in the 950-2020nm 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 proc… ▽ More The Near-Infrared Spectrometer and Photometer (NISP) on board the Euclid satellite provides multiband photometry and R>=450 slitless grism spectroscopy in the 950-2020nm 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, and its technical potentials and limitations. Links to articles providing more details and technical background are included. NISP's 16 HAWAII-2RG (H2RG) detectors with a plate scale of 0.3" pix^-1 deliver a field-of-view of 0.57deg^2. In photo mode, NISP reaches a limiting magnitude of ~24.5AB mag in three photometric exposures of about 100s exposure time, for point sources and with a signal-to-noise ratio (SNR) of 5. For spectroscopy, NISP's point-source sensitivity is a SNR = 3.5 detection of an emission line with flux ~2x10^-16erg/s/cm^2 integrated over two resolution elements of 13.4A, in 3x560s grism exposures at 1.6 mu (redshifted Ha). Our calibration includes on-ground and in-flight characterisation and monitoring of detector baseline, dark current, non-linearity, and sensitivity, to guarantee a relative photometric accuracy of better than 1.5%, and relative spectrophotometry to better than 0.7%. The wavelength calibration must be better than 5A. NISP is the state-of-the-art instrument in the NIR for all science beyond small areas available from HST and JWST - and an enormous advance due to its combination of field size and high throughput of telescope and instrument. During Euclid's 6-year survey covering 14000 deg^2 of extragalactic sky, NISP will be the backbone for determining distances of more than a billion galaxies. Its NIR data will become a rich reference imaging and spectroscopy data set for the coming decades. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue 'Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid. II. The VIS Instrument

Authors: Euclid Collaboration, M. Cropper, A. Al-Bahlawan, J. Amiaux, S. Awan, R. Azzollini, K. Benson, M. Berthe, J. Boucher, E. Bozzo, C. Brockley-Blatt, G. P. Candini, C. Cara, R. A. Chaudery, R. E. Cole, P. Danto, J. Denniston, A. M. Di Giorgio, B. Dryer, J. Endicott, J. -P. Dubois, M. Farina, E. Galli, L. Genolet, J. P. D. Gow , et al. (403 additional authors not shown)

Abstract: This paper presents the specification, design, and development of the Visible Camera (VIS) on the ESA Euclid mission. VIS is a large optical-band imager with a field of view of 0.54 deg^2 sampled at 0.1" with an array of 609 Megapixels and spatial resolution of 0.18". It will be used to survey approximately 14,000 deg^2 of extragalactic sky to measure the distortion of galaxies in the redshift ran… ▽ More This paper presents the specification, design, and development of the Visible Camera (VIS) on the ESA Euclid mission. VIS is a large optical-band imager with a field of view of 0.54 deg^2 sampled at 0.1" with an array of 609 Megapixels and spatial resolution of 0.18". It will be used to survey approximately 14,000 deg^2 of extragalactic sky to measure the distortion of galaxies in the redshift range z=0.1-1.5 resulting from weak gravitational lensing, one of the two principal cosmology probes of Euclid. With photometric redshifts, the distribution of dark matter can be mapped in three dimensions, and, from how this has changed with look-back time, the nature of dark energy and theories of gravity can be constrained. The entire VIS focal plane will be transmitted to provide the largest images of the Universe from space to date, reaching m_AB>24.5 with S/N >10 in a single broad I_E~(r+i+z) band over a six year survey. The particularly challenging aspects of the instrument are the control and calibration of observational biases, which lead to stringent performance requirements and calibration regimes. With its combination of spatial resolution, calibration knowledge, depth, and area covering most of the extra-Galactic sky, VIS will also provide a legacy data set for many other fields. This paper discusses the rationale behind the VIS concept and describes the instrument design and development before reporting the pre-launch performance derived from ground calibrations and brief results from the in-orbit commissioning. VIS should reach fainter than m_AB=25 with S/N>10 for galaxies of full-width half-maximum of 0.3" in a 1.3" diameter aperture over the Wide Survey, and m_AB>26.4 for a Deep Survey that will cover more than 50 deg^2. The paper also describes how VIS works with the other Euclid components of survey, telescope, and science data processing to extract the cosmological information. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

Spectral Difference method with a posteriori limiting: II- Application to low Mach number flows

Authors: D. A. Velasco-Romero, R. Teyssier

Abstract: Stellar convection poses two main gargantuan challenges for astrophysical fluid solvers: low-Mach number flows and minuscule perturbations over steeply stratified hydrostatic equilibria. Most methods exhibit excessive numerical diffusion and are unable to capture the correct solution due to large truncation errors. In this paper, we analyze the performance of the Spectral Difference (SD) method un… ▽ More Stellar convection poses two main gargantuan challenges for astrophysical fluid solvers: low-Mach number flows and minuscule perturbations over steeply stratified hydrostatic equilibria. Most methods exhibit excessive numerical diffusion and are unable to capture the correct solution due to large truncation errors. In this paper, we analyze the performance of the Spectral Difference (SD) method under these extreme conditions using an arbitrarily high-order shock capturing scheme with a posteriori limiting. We include both a modification to the HLLC Riemann solver adapted to low Mach number flows (L-HLLC) and a well-balanced scheme to properly evolve perturbations over steep equilibrium solutions. We evaluate the performance of our method using a series of test tailored specifically for stellar convection. We observe that our high-order SD method is capable of dealing with very subsonic flows without necessarily using the modified Riemann solver. We find however that the well-balanced framework is unavoidable if one wants to capture accurately small amplitude convective and acoustic modes. Analyzing the temporal and spatial evolution of the turbulent kinetic energy, we show that our fourth-order SD scheme seems to emerge as an optimal variant to solve this difficult numerical problem. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Comments: 15 pages, 11 figures, Submitted to MNRAS

Euclid preparation. Sensitivity to neutrino parameters

Authors: Euclid Collaboration, M. Archidiacono, J. Lesgourgues, S. Casas, S. Pamuk, N. Schöneberg, Z. Sakr, G. Parimbelli, A. Schneider, F. Hervas Peters, F. Pace, V. M. Sabarish, M. Costanzi, S. Camera, C. Carbone, S. Clesse, N. Frusciante, A. Fumagalli, P. Monaco, D. Scott, M. Viel, A. Amara, S. Andreon, N. Auricchio, M. Baldi , et al. (224 additional authors not shown)

Abstract: 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. We present forecasts from the combination of these surveys on the sensitivity to cosmological parameters including the summed neutrino mass $M_ν$ and the effective number of relativistic species… ▽ More 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. We present forecasts from the combination of these surveys on the sensitivity to cosmological parameters including the summed neutrino mass $M_ν$ and the effective number of relativistic species $N_{\rm eff}$ in the standard $Λ$CDM scenario and in a scenario with dynamical dark energy ($w_0 w_a$CDM). We compare the accuracy of different algorithms predicting the nonlinear matter power spectrum for such models. We then validate several pipelines for Fisher matrix and MCMC forecasts, using different theory codes, algorithms for numerical derivatives, and assumptions concerning the non-linear cut-off scale. The Euclid primary probes alone will reach a sensitivity of $σ(M_ν)=$56meV in the $Λ$CDM+$M_ν$ model, whereas the combination with CMB data from Planck is expected to achieve $σ(M_ν)=$23meV and raise the evidence for a non-zero neutrino mass to at least the $2.6σ$ level. This can be pushed to a $4σ$ detection if future CMB data from LiteBIRD and CMB Stage-IV are included. In combination with Planck, Euclid will also deliver tight constraints on $ΔN_{\rm eff}< 0.144$ (95%CL) in the $Λ$CDM+$M_ν$+$N_{\rm eff}$ model, or $ΔN_{\rm eff}< 0.063$ when future CMB data are included. When floating $(w_0, w_a)$, we find that the sensitivity to $N_{\rm eff}$ remains stable, while that to $M_ν$ degrades at most by a factor 2. This work illustrates the complementarity between the Euclid spectroscopic and imaging/photometric surveys and between Euclid and CMB constraints. Euclid will have a great potential for measuring the neutrino mass and excluding well-motivated scenarios with additional relativistic particles. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: 48 pages, 33 figures

Report number: TTK-24-17

Introducing the DREAMS Project: DaRk mattEr and Astrophysics with Machine learning and Simulations

Authors: Jonah C. Rose, Paul Torrey, Francisco Villaescusa-Navarro, Mariangela Lisanti, Tri Nguyen, Sandip Roy, Kassidy E. Kollmann, Mark Vogelsberger, Francis-Yan Cyr-Racine, Mikhail V. Medvedev, Shy Genel, Daniel Anglés-Alcázar, Nitya Kallivayalil, Bonny Y. Wang, Belén Costanza, Stephanie O'Neil, Cian Roche, Soumyodipta Karmakar, Alex M. Garcia, Ryan Low, Shurui Lin, Olivia Mostow, Akaxia Cruz, Andrea Caputo, Arya Farahi , et al. (5 additional authors not shown)

Abstract: We introduce the DREAMS project, an innovative approach to understanding the astrophysical implications of alternative dark matter models and their effects on galaxy formation and evolution. The DREAMS project will ultimately comprise thousands of cosmological hydrodynamic simulations that simultaneously vary over dark matter physics, astrophysics, and cosmology in modeling a range of systems -- f… ▽ More We introduce the DREAMS project, an innovative approach to understanding the astrophysical implications of alternative dark matter models and their effects on galaxy formation and evolution. The DREAMS project will ultimately comprise thousands of cosmological hydrodynamic simulations that simultaneously vary over dark matter physics, astrophysics, and cosmology in modeling a range of systems -- from galaxy clusters to ultra-faint satellites. Such extensive simulation suites can provide adequate training sets for machine-learning-based analyses. This paper introduces two new cosmological hydrodynamical suites of Warm Dark Matter, each comprised of 1024 simulations generated using the Arepo code. One suite consists of uniform-box simulations covering a $(25~h^{-1}~{\rm M}_\odot)^3$ volume, while the other consists of Milky Way zoom-ins with sufficient resolution to capture the properties of classical satellites. For each simulation, the Warm Dark Matter particle mass is varied along with the initial density field and several parameters controlling the strength of baryonic feedback within the IllustrisTNG model. We provide two examples, separately utilizing emulators and Convolutional Neural Networks, to demonstrate how such simulation suites can be used to disentangle the effects of dark matter and baryonic physics on galactic properties. The DREAMS project can be extended further to include different dark matter models, galaxy formation physics, and astrophysical targets. In this way, it will provide an unparalleled opportunity to characterize uncertainties on predictions for small-scale observables, leading to robust predictions for testing the particle physics nature of dark matter on these scales. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 28 pages, 8 figures, DREAMS website: https://www.dreams-project.org

Euclid preparation. LensMC, weak lensing cosmic shear measurement with forward modelling and Markov Chain Monte Carlo sampling

Authors: Euclid Collaboration, G. Congedo, L. Miller, A. N. Taylor, N. Cross, C. A. J. Duncan, T. Kitching, N. Martinet, S. Matthew, T. Schrabback, M. Tewes, N. Welikala, N. Aghanim, A. Amara, S. Andreon, N. Auricchio, M. Baldi, S. Bardelli, R. Bender, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia, J. Brinchmann, S. Camera , et al. (217 additional authors not shown)

Abstract: LensMC is a weak lensing shear measurement method developed for Euclid and Stage-IV surveys. It is based on forward modelling to deal with convolution by a point spread function with comparable size to many galaxies; sampling the posterior distribution of galaxy parameters via Markov Chain Monte Carlo; and marginalisation over nuisance parameters for each of the 1.5 billion galaxies observed by Eu… ▽ More LensMC is a weak lensing shear measurement method developed for Euclid and Stage-IV surveys. It is based on forward modelling to deal with convolution by a point spread function with comparable size to many galaxies; sampling the posterior distribution of galaxy parameters via Markov Chain Monte Carlo; and marginalisation over nuisance parameters for each of the 1.5 billion galaxies observed by Euclid. The scientific performance is quantified through high-fidelity images based on the Euclid Flagship simulations and emulation of the Euclid VIS images; realistic clustering with a mean surface number density of 250 arcmin$^{-2}$ ($I_{\rm E}<29.5$) for galaxies, and 6 arcmin$^{-2}$ ($I_{\rm E}<26$) for stars; and a diffraction-limited chromatic point spread function with a full width at half maximum of $0.^{\!\prime\prime}2$ and spatial variation across the field of view. Objects are measured with a density of 90 arcmin$^{-2}$ ($I_{\rm E}<26.5$) in 4500 deg$^2$. The total shear bias is broken down into measurement (our main focus here) and selection effects (which will be addressed elsewhere). We find: measurement multiplicative and additive biases of $m_1=(-3.6\pm0.2)\times10^{-3}$, $m_2=(-4.3\pm0.2)\times10^{-3}$, $c_1=(-1.78\pm0.03)\times10^{-4}$, $c_2=(0.09\pm0.03)\times10^{-4}$; a large detection bias with a multiplicative component of $1.2\times10^{-2}$ and an additive component of $-3\times10^{-4}$; and a measurement PSF leakage of $α_1=(-9\pm3)\times10^{-4}$ and $α_2=(2\pm3)\times10^{-4}$. When model bias is suppressed, the obtained measurement biases are close to Euclid requirement and largely dominated by undetected faint galaxies ($-5\times10^{-3}$). Although significant, model bias will be straightforward to calibrate given the weak sensitivity. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 28 pages, 18 figures, 2 tables

Comparison between a priori and a posteriori slope limiters for high-order finite volume schemes

Authors: Jonathan Palafoutas, David A. Velasco Romero, Romain Teyssier

Abstract: High-order finite volume and finite element methods offer impressive accuracy and cost efficiency when solving hyperbolic conservation laws with smooth solutions. However, if the solution contains discontinuities, these high-order methods can introduce unphysical oscillations and severe overshoots/undershoots. Slope limiters are an effective remedy, combating these oscillations by preserving monot… ▽ More High-order finite volume and finite element methods offer impressive accuracy and cost efficiency when solving hyperbolic conservation laws with smooth solutions. However, if the solution contains discontinuities, these high-order methods can introduce unphysical oscillations and severe overshoots/undershoots. Slope limiters are an effective remedy, combating these oscillations by preserving monotonicity. Some limiters can even maintain a strict maximum principle in the numerical solution. They can be classified into one of two categories: \textit{a priori} and \textit{a posteriori} limiters. The former revises the high-order solution based only on data at the current time $t^n$, while the latter involves computing a candidate solution at $t^{n+1}$ and iteratively recomputing it until some conditions are satisfied. These two limiting paradigms are available for both finite volume and finite element methods. In this work, we develop a methodology to compare \textit{a priori} and \textit{a posteriori} limiters for finite volume solvers at arbitrarily high order. We select the maximum principle preserving scheme presented in \cite{zhang2011maximum, zhang2010maximum} as our \textit{a priori} limited scheme. For \textit{a posteriori} limiting, we adopt the methodology presented in \cite{clain2011high} and search for so-called \textit{troubled cells} in the candidate solution. We revise them with a robust MUSCL fallback scheme. The linear advection equation is solved in both one and two dimensions and we compare variations of these limited schemes based on their ability to maintain a maximum principle, solution quality over long time integration and computational cost. ... △ Less

Submitted 27 April, 2024; originally announced April 2024.

Euclid preparation. Improving cosmological constraints using a new multi-tracer method with the spectroscopic and photometric samples

Authors: Euclid Collaboration, F. Dournac, A. Blanchard, S. Ilić, B. Lamine, I. Tutusaus, A. Amara, S. Andreon, N. Auricchio, H. Aussel, M. Baldi, S. Bardelli, C. Bodendorf, D. Bonino, E. Branchini, S. Brau-Nogue, M. Brescia, J. Brinchmann, S. Camera, V. Capobianco, J. Carretero, S. Casas, M. Castellano, S. Cavuoti, A. Cimatti , et al. (218 additional authors not shown)

Abstract: Future data provided by the Euclid mission will allow us to better understand the cosmic history of the Universe. A metric of its performance is the figure-of-merit (FoM) of dark energy, usually estimated with Fisher forecasts. The expected FoM has previously been estimated taking into account the two main probes of Euclid, namely the three-dimensional clustering of the spectroscopic galaxy sample… ▽ More Future data provided by the Euclid mission will allow us to better understand the cosmic history of the Universe. A metric of its performance is the figure-of-merit (FoM) of dark energy, usually estimated with Fisher forecasts. The expected FoM has previously been estimated taking into account the two main probes of Euclid, namely the three-dimensional clustering of the spectroscopic galaxy sample, and the so-called 3x2pt signal from the photometric sample (i.e., the weak lensing signal, the galaxy clustering, and their cross-correlation). So far, these two probes have been treated as independent. In this paper, we introduce a new observable given by the ratio of the (angular) two-point correlation function of galaxies from the two surveys. For identical (normalised) selection functions, this observable is unaffected by sampling noise, and its variance is solely controlled by Poisson noise. We present forecasts for Euclid where this multi-tracer method is applied and is particularly relevant because the two surveys will cover the same area of the sky. This method allows for the exploitation of the combination of the spectroscopic and photometric samples. When the correlation between this new observable and the other probes is not taken into account, a significant gain is obtained in the FoM, as well as in the constraints on other cosmological parameters. The benefit is more pronounced for a commonly investigated modified gravity model, namely the $γ$ parametrisation of the growth factor. However, the correlation between the different probes is found to be significant and hence the actual gain is uncertain. We present various strategies for circumventing this issue and still extract useful information from the new observable. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: 16 pages, 7 figures, submitted to A&A

Euclid preparation. XLII. A unified catalogue-level reanalysis of weak lensing by galaxy clusters in five imaging surveys

Authors: Euclid Collaboration, M. Sereno, S. Farrens, L. Ingoglia, G. F. Lesci, L. Baumont, G. Covone, C. Giocoli, F. Marulli, S. Miranda La Hera, M. Vannier, A. Biviano, S. Maurogordato, L. Moscardini, N. Aghanim, S. Andreon, N. Auricchio, M. Baldi, S. Bardelli, F. Bellagamba, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia, J. Brinchmann , et al. (199 additional authors not shown)

Abstract: Precise and accurate mass calibration is required to exploit galaxy clusters as astrophysical and cosmological probes in the Euclid era. Systematic errors in lensing signals by galaxy clusters can be empirically estimated by comparing different surveys with independent and uncorrelated systematics. To assess the robustness of the lensing results to systematic errors, we carried out end-to-end test… ▽ More Precise and accurate mass calibration is required to exploit galaxy clusters as astrophysical and cosmological probes in the Euclid era. Systematic errors in lensing signals by galaxy clusters can be empirically estimated by comparing different surveys with independent and uncorrelated systematics. To assess the robustness of the lensing results to systematic errors, we carried out end-to-end tests across different data sets. We performed a unified analysis at the catalogue level by leveraging the Euclid combined cluster and weak-lensing pipeline (COMB-CL). COMB-CL will measure weak lensing cluster masses for the Euclid Survey. Heterogeneous data sets from five independent, recent, lensing surveys (CHFTLenS, DES~SV1, HSC-SSP~S16a, KiDS~DR4, and RCSLenS), which exploited different shear and photometric redshift estimation algorithms, were analysed with a consistent pipeline under the same model assumptions. We performed a comparison of the amplitude of the reduced excess surface density and of the mass estimates using lenses from the Planck PSZ2 and SDSS redMaPPer cluster samples. Mass estimates agree with literature results collected in the LC2 catalogues. Mass accuracy was further investigated considering the AMICO detected clusters in the HSC-SSP XXL North field. The consistency of the data sets was tested using our unified analysis framework. We found agreement between independent surveys, at the level of systematic noise in Stage-III surveys or precursors. This indicates successful control over systematics. If such control continues in Stage-IV, Euclid will be able to measure the weak lensing masses of around 13000 (considering shot noise only) or 3000 (noise from shape and large-scale-structure) massive clusters with a signal-to-noise ratio greater than 3. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: 25 pages; in press on A&A

Euclid preparation. Measuring detailed galaxy morphologies for Euclid with Machine Learning

Authors: Euclid Collaboration, B. Aussel, S. Kruk, M. Walmsley, M. Huertas-Company, M. Castellano, C. J. Conselice, M. Delli Veneri, H. Domínguez Sánchez, P. -A. Duc, U. Kuchner, A. La Marca, B. Margalef-Bentabol, F. R. Marleau, G. Stevens, Y. Toba, C. Tortora, L. Wang, N. Aghanim, B. Altieri, A. Amara, S. Andreon, N. Auricchio, M. Baldi, S. Bardelli , et al. (233 additional authors not shown)

Abstract: The Euclid mission is expected to image millions of galaxies with high resolution, providing an extensive dataset to study galaxy evolution. We investigate the application of deep learning to predict the detailed morphologies of galaxies in Euclid using Zoobot a convolutional neural network pretrained with 450000 galaxies from the Galaxy Zoo project. We adapted Zoobot for emulated Euclid images, g… ▽ More The Euclid mission is expected to image millions of galaxies with high resolution, providing an extensive dataset to study galaxy evolution. We investigate the application of deep learning to predict the detailed morphologies of galaxies in Euclid using Zoobot a convolutional neural network pretrained with 450000 galaxies from the Galaxy Zoo project. We adapted Zoobot for emulated Euclid images, generated based on Hubble Space Telescope COSMOS images, and with labels provided by volunteers in the Galaxy Zoo: Hubble project. We demonstrate that the trained Zoobot model successfully measures detailed morphology for emulated Euclid images. It effectively predicts whether a galaxy has features and identifies and characterises various features such as spiral arms, clumps, bars, disks, and central bulges. When compared to volunteer classifications Zoobot achieves mean vote fraction deviations of less than 12% and an accuracy above 91% for the confident volunteer classifications across most morphology types. However, the performance varies depending on the specific morphological class. For the global classes such as disk or smooth galaxies, the mean deviations are less than 10%, with only 1000 training galaxies necessary to reach this performance. For more detailed structures and complex tasks like detecting and counting spiral arms or clumps, the deviations are slightly higher, around 12% with 60000 galaxies used for training. In order to enhance the performance on complex morphologies, we anticipate that a larger pool of labelled galaxies is needed, which could be obtained using crowdsourcing. Finally, our findings imply that the model can be effectively adapted to new morphological labels. We demonstrate this adaptability by applying Zoobot to peculiar galaxies. In summary, our trained Zoobot CNN can readily predict morphological catalogues for Euclid images. △ Less

Submitted 15 February, 2024; originally announced February 2024.

Comments: 27 pages, 26 figures, 5 tables, submitted to A&A

The AGORA High-resolution Galaxy Simulations Comparison Project IV: Halo and Galaxy Mass Assembly in a Cosmological Zoom-in Simulation at $z\le2$

Authors: Santi Roca-Fàbrega, Ji-hoon Kim, Joel R. Primack, Minyong Jung, Anna Genina, Loic Hausammann, Hyeonyong Kim, Alessandro Lupi, Kentaro Nagamine, Johnny W. Powell, Yves Revaz, Ikkoh Shimizu, Clayton Strawn, Héctor Velázquez, Tom Abel, Daniel Ceverino, Bili Dong, Thomas R. Quinn, Eun-jin Shin, Alvaro Segovia-Otero, Oscar Agertz, Kirk S. S. Barrow, Corentin Cadiou, Avishai Dekel, Cameron Hummels , et al. (3 additional authors not shown)

Abstract: In this fourth paper from the AGORA Collaboration, we study the evolution down to redshift $z=2$ and below of a set of cosmological zoom-in simulations of a Milky Way mass galaxy by eight of the leading hydrodynamic simulation codes. We also compare this CosmoRun suite of simulations with dark matter-only simulations by the same eight codes. We analyze general properties of the halo and galaxy at… ▽ More In this fourth paper from the AGORA Collaboration, we study the evolution down to redshift $z=2$ and below of a set of cosmological zoom-in simulations of a Milky Way mass galaxy by eight of the leading hydrodynamic simulation codes. We also compare this CosmoRun suite of simulations with dark matter-only simulations by the same eight codes. We analyze general properties of the halo and galaxy at $z=4$ and 3, and before the last major merger, focusing on the formation of well-defined rotationally-supported disks, the mass-metallicity relation, the specific star formation rate, the gas metallicity gradients, and the non-axisymmetric structures in the stellar disks. Codes generally converge well to the stellar-to-halo mass ratios predicted by semi-analytic models at $z\sim$2. We see that almost all the hydro codes develop rotationally-supported structures at low redshifts. Most agree within 0.5 dex with the observed MZR at high and intermediate redshifts, and reproduce the gas metallicity gradients obtained from analytical models and low-redshift observations. We confirm that the inter-code differences in the halo assembly history reported in the first paper of the collaboration also exist in CosmoRun, making the code-to-code comparison more difficult. We show that such differences are mainly due to variations in code-dependent parameters that control the time-stepping strategy of the gravity solver. We find that variations in the early stellar feedback can also result in differences in the timing of the low-redshift mergers. All the simulation data down to $z=2$ and the auxiliary data will be made publicly available. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Comments: Submitted to ApJ. Visit the AGORA Collaboration website (www.agorasimulations.org <http://www.agorasimulations.org/>) for more information

The AGORA High-resolution Galaxy Simulations Comparison Project. V: Satellite Galaxy Populations In A Cosmological Zoom-in Simulation of A Milky Way-mass Halo

Authors: Minyong Jung, Santi Roca-Fàbrega, Ji-hoon Kim, Anna Genina, Loic Hausammann, Hyeonyong Kim, Alessandro Lupi, Kentaro Nagamine, Johnny W. Powell, Yves Revaz, Ikkoh Shimizu, Héctor Velázquez, Daniel Ceverino, Joel R. Primack, Thomas R. Quinn, Clayton Strawn, Tom Abel, Avishai Dekel, Bili Dong, Boon Kiat Oh, Romain Teyssier

Abstract: We analyze and compare the satellite halo populations at $z\sim2$ in the high-resolution cosmological zoom-in simulations of a $10^{12}\,{\rm M}_{\odot}$ target halo ($z=0$ mass) carried out on eight widely-used astrophysical simulation codes ({\sc Art-I}, {\sc Enzo}, {\sc Ramses}, {\sc Changa}, {\sc Gadget-3}, {\sc Gear}, {\sc Arepo-t}, and {\sc Gizmo}) for the {\it AGORA} High-resolution Galaxy… ▽ More We analyze and compare the satellite halo populations at $z\sim2$ in the high-resolution cosmological zoom-in simulations of a $10^{12}\,{\rm M}_{\odot}$ target halo ($z=0$ mass) carried out on eight widely-used astrophysical simulation codes ({\sc Art-I}, {\sc Enzo}, {\sc Ramses}, {\sc Changa}, {\sc Gadget-3}, {\sc Gear}, {\sc Arepo-t}, and {\sc Gizmo}) for the {\it AGORA} High-resolution Galaxy Simulations Comparison Project. We use slightly different redshift epochs near $z=2$ for each code (hereafter ``$z\sim2$') at which the eight simulations are in the same stage in the target halo's merger history. After identifying the matched pairs of halos between the {\it CosmoRun} simulations and the DMO simulations, we discover that each {\it CosmoRun} halo tends to be less massive than its DMO counterpart. When we consider only the halos containing stellar particles at $z\sim2$, the number of satellite {\it galaxies} is significantly fewer than that of dark matter halos in all participating {\it AGORA} simulations, and is comparable to the number of present-day satellites near the Milky Way or M31. The so-called ``missing satellite problem' is fully resolved across all participating codes simply by implementing the common baryonic physics adopted in {\it AGORA} and the stellar feedback prescription commonly used in each code, with sufficient numerical resolution ($\lesssim100$ proper pc at $z=2$). We also compare other properties such as the stellar mass$-$halo mass relation and the mass$-$metallicity relation. Our work highlights the value of comparison studies such as {\it AGORA}, where outstanding problems in galaxy formation theory are studied simultaneously on multiple numerical platforms. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: Accepted for publication to the ApJ, 19 pages, 10 figures. Visit the AGORA Collaboration website (www.agorasimulations.org) for more information. For summary video, please see https://www.youtube.com/watch?v=HShRJYvPe1k

Journal ref: ApJ 964 123 (2024)

The AGORA High-resolution Galaxy Simulations Comparison Project. VI. Similarities and Differences in the Circumgalactic Medium

Authors: Clayton Strawn, Santi Roca-Fàbrega, Joel R. Primack, Ji-hoon Kim, Anna Genina, Loic Hausammann, Hyeonyong Kim, Alessandro Lupi, Kentaro Nagamine, Johnny W. Powell, Yves Revaz, Ikkoh Shimizu, Héctor Velázquez, Tom Abel, Daniel Ceverino, Bili Dong, Minyong Jung, Thomas R. Quinn, Eun-jin Shin, Kirk S. S. Barrow, Avishai Dekel, Boon Kiat Oh, Nir Mandelker, Romain Teyssier, Cameron Hummels , et al. (4 additional authors not shown)

Abstract: We analyze the circumgalactic medium (CGM) for eight commonly-used cosmological codes in the AGORA collaboration. The codes are calibrated to use identical initial conditions, cosmology, heating and cooling, and star formation thresholds, but each evolves with its own unique code architecture and stellar feedback implementation. Here, we analyze the results of these simulations in terms of the str… ▽ More We analyze the circumgalactic medium (CGM) for eight commonly-used cosmological codes in the AGORA collaboration. The codes are calibrated to use identical initial conditions, cosmology, heating and cooling, and star formation thresholds, but each evolves with its own unique code architecture and stellar feedback implementation. Here, we analyze the results of these simulations in terms of the structure, composition, and phase dynamics of the CGM. We show properties such as metal distribution, ionization levels, and kinematics are effective tracers of the effects of the different code feedback and implementation methods, and as such they can be highly divergent between simulations. This is merely a fiducial set of models, against which we will in the future compare multiple feedback recipes for each code. Nevertheless, we find that the large parameter space these simulations establish can help disentangle the different variables that affect observable quantities in the CGM, e.g., showing that abundances for ions with higher ionization energy are more strongly determined by the simulation's metallicity, while abundances for ions with lower ionization energy are more strongly determined by the gas density and temperature. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 32 pages, 17 figures. Published in ApJ

Journal ref: ApJ 962 29 (2024)

Euclid preparation. Optical emission-line predictions of intermediate-z galaxy populations in GAEA for the Euclid Deep and Wide Surveys

Authors: Euclid Collaboration, L. Scharré, M. Hirschmann, G. De Lucia, S. Charlot, F. Fontanot, M. Spinelli, L. Xie, A. Feltre, V. Allevato, A. Plat, M. N. Bremer, S. Fotopoulou, L. Gabarra, B. R. Granett, M. Moresco, C. Scarlata, L. Pozzetti, L. Spinoglio, M. Talia, G. Zamorani, B. Altieri, A. Amara, S. Andreon, N. Auricchio , et al. (217 additional authors not shown)

Abstract: In anticipation of the Euclid Wide and Deep Surveys, we present optical emission-line predictions at intermediate redshifts from 0.4 to 2.5. Our approach combines a mock light cone from the GAEA semi-analytic model to self-consistently model nebular emission from HII regions, narrow-line regions of active galactic nuclei (AGN), and evolved stellar populations. Our analysis focuses on seven optical… ▽ More In anticipation of the Euclid Wide and Deep Surveys, we present optical emission-line predictions at intermediate redshifts from 0.4 to 2.5. Our approach combines a mock light cone from the GAEA semi-analytic model to self-consistently model nebular emission from HII regions, narrow-line regions of active galactic nuclei (AGN), and evolved stellar populations. Our analysis focuses on seven optical emission lines: H$α$, H$β$, [SII]$λλ6717, 6731$, [NII]$λ6584$, [OI]$λ6300$, [OIII]$λ5007$, and [OII]$λλ3727, 3729$. We find that Euclid will predominantly observe massive, star-forming, and metal-rich line-emitters. Interstellar dust, modelled using a Calzetti law with mass-dependent scaling, may decrease observable percentages by a further 20-30% with respect to our underlying emission-line populations from GAEA. We predict Euclid to observe around 30-70% of H$α$-, [NII]-, [SII]-, and [OIII]-emitting galaxies at redshift below 1 and under 10% at higher redshift. Observability of H$β$-, [OII]-, and [OI]- emission is limited to below 5%. For the Euclid-observable sample, we find that BPT diagrams can effectively distinguish between different galaxy types up to around redshift 1.8, attributed to the bias toward metal-rich systems. Moreover, we show that the relationships of H$α$ and [OIII]+H$β$ to the star-formation rate, and the [OIII]-AGN luminosity relation, exhibit minimal changes with increasing redshift. Based on line ratios [NII]/H$α$, [NII]/[OII], and [NII]/[SII], we further propose novel z-invariant tracers for the black hole accretion rate-to-star formation rate ratio. Lastly, we find that commonly used metallicity estimators display gradual shifts in normalisations with increasing redshift, while maintaining the overall shape of local calibrations. This is in tentative agreement with recent JWST data. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: 32 pages, 15 figures, submitted to A&A

Euclid preparation XLVI. The Near-IR Background Dipole Experiment with Euclid

Authors: Euclid Collaboration, A. Kashlinsky, R. G. Arendt, M. L. N. Ashby, F. Atrio-Barandela, R. Scaramella, M. A. Strauss, B. Altieri, A. Amara, S. Andreon, N. Auricchio, M. Baldi, S. Bardelli, R. Bender, C. Bodendorf, E. Branchini, M. Brescia, J. Brinchmann, S. Camera, V. Capobianco, C. Carbone, J. Carretero, S. Casas, M. Castellano, S. Cavuoti , et al. (195 additional authors not shown)

Abstract: Verifying the fully kinematic nature of the cosmic microwave background (CMB) dipole is of fundamental importance in cosmology. In the standard cosmological model with the Friedman-Lemaitre-Robertson-Walker (FLRW) metric from the inflationary expansion the CMB dipole should be entirely kinematic. Any non-kinematic CMB dipole component would thus reflect the preinflationary structure of spacetime p… ▽ More Verifying the fully kinematic nature of the cosmic microwave background (CMB) dipole is of fundamental importance in cosmology. In the standard cosmological model with the Friedman-Lemaitre-Robertson-Walker (FLRW) metric from the inflationary expansion the CMB dipole should be entirely kinematic. Any non-kinematic CMB dipole component would thus reflect the preinflationary structure of spacetime probing the extent of the FLRW applicability. Cosmic backgrounds from galaxies after the matter-radiation decoupling, should have kinematic dipole component identical in velocity with the CMB kinematic dipole. Comparing the two can lead to isolating the CMB non-kinematic dipole. It was recently proposed that such measurement can be done using the near-IR cosmic infrared background (CIB) measured with the currently operating Euclid telescope, and later with Roman. The proposed method reconstructs the resolved CIB, the Integrated Galaxy Light (IGL), from Euclid's Wide Survey and probes its dipole, with a kinematic component amplified over that of the CMB by the Compton-Getting effect. The amplification coupled with the extensive galaxy samples forming the IGL would determine the CIB dipole with an overwhelming signal/noise, isolating its direction to sub-degree accuracy. We develop details of the method for Euclid's Wide Survey in 4 bands spanning 0.6 to 2 mic. We isolate the systematic and other uncertainties and present methodologies to minimize them, after confining the sample to the magnitude range with negligible IGL/CIB dipole from galaxy clustering. These include the required star-galaxy separation, accounting for the extinction correction dipole using the method newly developed here achieving total separation, accounting for the Earth's orbital motion and other systematic effects. (Abridged) △ Less

Submitted 24 June, 2024; v1 submitted 31 January, 2024; originally announced January 2024.

Comments: Euclid Key Project paper, A&A, in press

Euclid preparation: TBD. The pre-launch Science Ground Segment simulation framework

Authors: Euclid Collaboration, S. Serrano, P. Hudelot, G. Seidel, J. E. Pollack, E. Jullo, F. Torradeflot, D. Benielli, R. Fahed, T. Auphan, J. Carretero, H. Aussel, P. Casenove, F. J. Castander, J. E. Davies, N. Fourmanoit, S. Huot, A. Kara, E. Keihanen, S. Kermiche, K. Okumura, J. Zoubian, A. Ealet, A. Boucaud, H. Bretonniere , et al. (251 additional authors not shown)

Abstract: 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 previous… ▽ More 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. △ Less

Submitted 2 January, 2024; originally announced January 2024.

Comments: 38 pages, 25 figures, A&A submitted

Euclid preparation. TBD. Galaxy power spectrum modelling in real space

Authors: Euclid Collaboration, A. Pezzotta, C. Moretti, M. Zennaro, A. Moradinezhad Dizgah, M. Crocce, E. Sefusatti, I. Ferrero, K. Pardede, A. Eggemeier, A. Barreira, R. E. Angulo, M. Marinucci, B. Camacho Quevedo, S. de la Torre, D. Alkhanishvili, M. Biagetti, M. -A. Breton, E. Castorina, G. D'Amico, V. Desjacques, M. Guidi, M. Kärcher, A. Oddo, M. Pellejero Ibanez , et al. (224 additional authors not shown)

Abstract: We investigate the accuracy of the perturbative galaxy bias expansion in view of the forthcoming analysis of the Euclid spectroscopic galaxy samples. We compare the performance of an Eulerian galaxy bias expansion, using state-of-art prescriptions from the effective field theory of large-scale structure (EFTofLSS), against a hybrid approach based on Lagrangian perturbation theory and high-resoluti… ▽ More We investigate the accuracy of the perturbative galaxy bias expansion in view of the forthcoming analysis of the Euclid spectroscopic galaxy samples. We compare the performance of an Eulerian galaxy bias expansion, using state-of-art prescriptions from the effective field theory of large-scale structure (EFTofLSS), against a hybrid approach based on Lagrangian perturbation theory and high-resolution simulations. These models are benchmarked against comoving snapshots of the Flagship I N-body simulation at $z=(0.9,1.2,1.5,1.8)$, which have been populated with H$α$ galaxies leading to catalogues of millions of objects within a volume of about $58\,h^{-3}\,{\rm Gpc}^3$. Our analysis suggests that both models can be used to provide a robust inference of the parameters $(h, ω_{\rm c})$ in the redshift range under consideration, with comparable constraining power. We additionally determine the range of validity of the EFTofLSS model in terms of scale cuts and model degrees of freedom. From these tests, it emerges that the standard third-order Eulerian bias expansion can accurately describe the full shape of the real-space galaxy power spectrum up to the maximum wavenumber $k_{\rm max}=0.45\,h\,{\rm Mpc}^{-1}$, even with a measurement precision well below the percent level. In particular, this is true for a configuration with six free nuisance parameters, including local and non-local bias parameters, a matter counterterm, and a correction to the shot-noise contribution. Fixing either tidal bias parameters to physically-motivated relations still leads to unbiased cosmological constraints. We finally repeat our analysis assuming a volume that matches the expected footprint of Euclid, but without considering observational effects, as purity and completeness, showing that we can get consistent cosmological constraints over this range of scales and redshifts. △ Less

Submitted 1 December, 2023; originally announced December 2023.

Comments: 38 pages, 19 figures

Euclid Preparation. XXXVII. Galaxy colour selections with Euclid and ground photometry for cluster weak-lensing analyses

Authors: Euclid Collaboration, G. F. Lesci, M. Sereno, M. Radovich, G. Castignani, L. Bisigello, F. Marulli, L. Moscardini, L. Baumont, G. Covone, S. Farrens, C. Giocoli, L. Ingoglia, S. Miranda La Hera, M. Vannier, A. Biviano, S. Maurogordato, N. Aghanim, A. Amara, S. Andreon, N. Auricchio, M. Baldi, S. Bardelli, R. Bender, C. Bodendorf , et al. (216 additional authors not shown)

Abstract: We derived galaxy colour selections from Euclid and ground-based photometry, aiming to accurately define background galaxy samples in cluster weak-lensing analyses. Given any set of photometric bands, we developed a method for the calibration of optimal galaxy colour selections that maximises the selection completeness, given a threshold on purity. We calibrated galaxy selections using simulated g… ▽ More We derived galaxy colour selections from Euclid and ground-based photometry, aiming to accurately define background galaxy samples in cluster weak-lensing analyses. Given any set of photometric bands, we developed a method for the calibration of optimal galaxy colour selections that maximises the selection completeness, given a threshold on purity. We calibrated galaxy selections using simulated ground-based $griz$ and Euclid $Y_{\rm E}J_{\rm E}H_{\rm E}$ photometry. Both selections produce a purity higher than 97%. The $griz$ selection completeness ranges from 30% to 84% in the lens redshift range $z_{\rm l}\in[0.2,0.8]$. With the full $grizY_{\rm E}J_{\rm E}H_{\rm E}$ selection, the completeness improves by up to $25$ percentage points, and the $z_{\rm l}$ range extends up to $z_{\rm l}=1.5$. The calibrated colour selections are stable to changes in the sample limiting magnitudes and redshift, and the selection based on $griz$ bands provides excellent results on real external datasets. The $griz$ selection is also purer at high redshift and more complete at low redshift compared to colour selections found in the literature. We find excellent agreement in terms of purity and completeness between the analysis of an independent, simulated Euclid galaxy catalogue and our calibration sample, except for galaxies at high redshifts, for which we obtain up to 50 percent points higher completeness. The combination of colour and photo-$z$ selections applied to simulated Euclid data yields up to 95% completeness, while the purity decreases down to 92% at high $z_{\rm l}$. We show that the calibrated colour selections provide robust results even when observations from a single band are missing from the ground-based data. Finally, we show that colour selections do not disrupt the shear calibration for stage III surveys. △ Less

Submitted 24 January, 2024; v1 submitted 27 November, 2023; originally announced November 2023.

Comments: 20 pages, 13 figures. Published by A&A

Euclid preparation. Modelling spectroscopic clustering on mildly nonlinear scales in beyond-$Λ$CDM models

Authors: Euclid Collaboration, B. Bose, P. Carrilho, M. Marinucci, C. Moretti, M. Pietroni, E. Carella, L. Piga, B. S. Wright, F. Vernizzi, C. Carbone, S. Casas, G. D'Amico, N. Frusciante, K. Koyama, F. Pace, A. Pourtsidou, M. Baldi, L. F. de la Bella, B. Fiorini, C. Giocoli, L. Lombriser, N. Aghanim, A. Amara, S. Andreon , et al. (207 additional authors not shown)

Abstract: We investigate the approximations needed to efficiently predict the large-scale clustering of matter and dark matter halos in beyond-$Λ$CDM scenarios. We examine the normal branch of the Dvali-Gabadadze-Porrati model, the Hu-Sawicki $f(R)$ model, a slowly evolving dark energy, an interacting dark energy model and massive neutrinos. For each, we test approximations for the perturbative kernel calcu… ▽ More We investigate the approximations needed to efficiently predict the large-scale clustering of matter and dark matter halos in beyond-$Λ$CDM scenarios. We examine the normal branch of the Dvali-Gabadadze-Porrati model, the Hu-Sawicki $f(R)$ model, a slowly evolving dark energy, an interacting dark energy model and massive neutrinos. For each, we test approximations for the perturbative kernel calculations, including the omission of screening terms and the use of perturbative kernels based on the Einstein-de Sitter universe; we explore different infrared-resummation schemes, tracer bias models and a linear treatment of massive neutrinos; we employ two models for redshift space distortions, the Taruya-Nishimishi-Saito prescription and the Effective Field Theory of Large-Scale Structure. This work further provides a preliminary validation of the codes being considered by Euclid for the spectroscopic clustering probe in beyond-$Λ$CDM scenarios. We calculate and compare the $χ^2$ statistic to assess the different modelling choices. This is done by fitting the spectroscopic clustering predictions to measurements from numerical simulations and perturbation theory-based mock data. We compare the behaviour of this statistic in the beyond-$Λ$CDM cases, as a function of the maximum scale included in the fit, to the baseline $Λ$CDM case. We find that the Einstein-de Sitter approximation without screening is surprisingly accurate for all cases when comparing to the halo clustering monopole and quadrupole obtained from simulations. Our results suggest that the inclusion of multiple redshift bins, higher-order multipoles, higher-order clustering statistics (such as the bispectrum) and photometric probes such as weak lensing, will be essential to extract information on massive neutrinos, modified gravity and dark energy. △ Less

Submitted 11 July, 2024; v1 submitted 22 November, 2023; originally announced November 2023.

Comments: 29 pages, 10 figures, 4 tables, 3 appendices. Journal accepted version

Euclid preparation. Spectroscopy of active galactic nuclei with NISP

Authors: Euclid Collaboration, E. Lusso, S. Fotopoulou, M. Selwood, V. Allevato, G. Calderone, C. Mancini, M. Mignoli, M. Scodeggio, L. Bisigello, A. Feltre, F. Ricci, F. La Franca, D. Vergani, L. Gabarra, V. Le Brun, E. Maiorano, E. Palazzi, M. Moresco, G. Zamorani, G. Cresci, K. Jahnke, A. Humphrey, H. Landt, F. Mannucci , et al. (224 additional authors not shown)

Abstract: 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… ▽ More 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$α$+[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$α$+[NII] is visible (0.89<z<1.83 at a line flux $>2x10^{-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. △ Less

Submitted 15 January, 2024; v1 submitted 20 November, 2023; originally announced November 2023.

Comments: 29 pages, 23 figures. Submitted to A&A, revised version

A new star formation recipe for magneto-hydrodynamics simulations of galaxy formation

Authors: Eden Girma, Romain Teyssier

Abstract: Star formation has been observed to occur at globally low yet locally varying efficiencies. As such, accurate capture of star formation in numerical simulations requires mechanisms that can replicate both its smaller-scale variations and larger-scale properties. Magnetic fields are thought to play an essential role within the turbulent interstellar medium (ISM) and affect molecular cloud collapse.… ▽ More Star formation has been observed to occur at globally low yet locally varying efficiencies. As such, accurate capture of star formation in numerical simulations requires mechanisms that can replicate both its smaller-scale variations and larger-scale properties. Magnetic fields are thought to play an essential role within the turbulent interstellar medium (ISM) and affect molecular cloud collapse. However, it remains to be fully explored how a magnetised model of star formation might influence galaxy evolution. We present a new model for a sub-grid star formation recipe that depends on the magnetic field. We run isolated disk galaxy simulations to assess its impact on the regulation of star formation using the code RAMSES. Building upon existing numerical methods, our model derives the star formation efficiency from local properties of the sub-grid magnetised ISM turbulence, assuming a constant Alfvén speed at sub-parsec scales. Compared to its non-magnetised counterpart, our star formation model suppresses the initial starburst by a factor of two, while regulating star formation later on to a nearly constant rate of $\sim 1~M_{\odot}~\mathrm{yr}^{-1}$. Differences also arise in the local Schmidt law with a shallower power law index for the magnetised star formation model. Our results encourage further examination into the notion that magnetic fields are likely to play a non-trivial role in our understanding of star and galaxy formation. △ Less

Submitted 22 November, 2023; v1 submitted 17 November, 2023; originally announced November 2023.

Comments: 15 pages, 11 figures; Accepted for publication in MNRAS

Euclid Preparation. TBD. Impact of magnification on spectroscopic galaxy clustering

Authors: Euclid Collaboration, G. Jelic-Cizmek, F. Sorrenti, F. Lepori, C. Bonvin, S. Camera, F. J. Castander, R. Durrer, P. Fosalba, M. Kunz, L. Lombriser, I. Tutusaus, C. Viglione, Z. Sakr, N. Aghanim, A. Amara, S. Andreon, M. Baldi, S. Bardelli, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia, J. Brinchmann, V. Capobianco , et al. (204 additional authors not shown)

Abstract: In this paper we investigate the impact of lensing magnification on the analysis of Euclid's spectroscopic survey, using the multipoles of the 2-point correlation function for galaxy clustering. We determine the impact of lensing magnification on cosmological constraints, and the expected shift in the best-fit parameters if magnification is ignored. We consider two cosmological analyses: i) a full… ▽ More In this paper we investigate the impact of lensing magnification on the analysis of Euclid's spectroscopic survey, using the multipoles of the 2-point correlation function for galaxy clustering. We determine the impact of lensing magnification on cosmological constraints, and the expected shift in the best-fit parameters if magnification is ignored. We consider two cosmological analyses: i) a full-shape analysis based on the $Λ$CDM model and its extension $w_0w_a$CDM and ii) a model-independent analysis that measures the growth rate of structure in each redshift bin. We adopt two complementary approaches in our forecast: the Fisher matrix formalism and the Markov chain Monte Carlo method. The fiducial values of the local count slope (or magnification bias), which regulates the amplitude of the lensing magnification, have been estimated from the Euclid Flagship simulations. We use linear perturbation theory and model the 2-point correlation function with the public code coffe. For a $Λ$CDM model, we find that the estimation of cosmological parameters is biased at the level of 0.4-0.7 standard deviations, while for a $w_0w_a$CDM dynamical dark energy model, lensing magnification has a somewhat smaller impact, with shifts below 0.5 standard deviations. In a model-independent analysis aiming to measure the growth rate of structure, we find that the estimation of the growth rate is biased by up to $1.2$ standard deviations in the highest redshift bin. As a result, lensing magnification cannot be neglected in the spectroscopic survey, especially if we want to determine the growth factor, one of the most promising ways to test general relativity with Euclid. We also find that, by including lensing magnification with a simple template, this shift can be almost entirely eliminated with minimal computational overhead. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 23 pages, 14 figures

Euclid preparation. XXXIX. The effect of baryons on the Halo Mass Function

Authors: Euclid Collaboration, T. Castro, S. Borgani, M. Costanzi, J. Dakin, K. Dolag, A. Fumagalli, A. Ragagnin, A. Saro, A. M. C. Le Brun, N. Aghanim, A. Amara, S. Andreon, N. Auricchio, M. Baldi, S. Bardelli, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia, J. Brinchmann, S. Camera, V. Capobianco, C. Carbone, J. Carretero , et al. (198 additional authors not shown)

Abstract: The Euclid photometric survey of galaxy clusters stands as a powerful cosmological tool, with the capacity to significantly propel our understanding of the Universe. Despite being sub-dominant to dark matter and dark energy, the baryonic component in our Universe holds substantial influence over the structure and mass of galaxy clusters. This paper presents a novel model to precisely quantify the… ▽ More The Euclid photometric survey of galaxy clusters stands as a powerful cosmological tool, with the capacity to significantly propel our understanding of the Universe. Despite being sub-dominant to dark matter and dark energy, the baryonic component in our Universe holds substantial influence over the structure and mass of galaxy clusters. This paper presents a novel model to precisely quantify the impact of baryons on galaxy cluster virial halo masses, using the baryon fraction within a cluster as proxy for their effect. Constructed on the premise of quasi-adiabaticity, the model includes two parameters calibrated using non-radiative cosmological hydrodynamical simulations and a single large-scale simulation from the Magneticum set, which includes the physical processes driving galaxy formation. As a main result of our analysis, we demonstrate that this model delivers a remarkable one percent relative accuracy in determining the virial dark matter-only equivalent mass of galaxy clusters, starting from the corresponding total cluster mass and baryon fraction measured in hydrodynamical simulations. Furthermore, we demonstrate that this result is robust against changes in cosmological parameters and against varying the numerical implementation of the sub-resolution physical processes included in the simulations. Our work substantiates previous claims about the impact of baryons on cluster cosmology studies. In particular, we show how neglecting these effects would lead to biased cosmological constraints for a Euclid-like cluster abundance analysis. Importantly, we demonstrate that uncertainties associated with our model, arising from baryonic corrections to cluster masses, are sub-dominant when compared to the precision with which mass-observable relations will be calibrated using Euclid, as well as our current understanding of the baryon fraction within galaxy clusters. △ Less

Submitted 16 April, 2024; v1 submitted 25 October, 2023; originally announced November 2023.

Comments: 18 pages, 10 figures, 4 tables, 1 appendix, abstract abridged for arXiv submission; v2 matches published version

Euclid preparation. TBD. Forecast impact of super-sample covariance on 3x2pt analysis with Euclid

Authors: Euclid Collaboration, D. Sciotti, S. Gouyou Beauchamps, V. F. Cardone, S. Camera, I. Tutusaus, F. Lacasa, A. Barreira, A. Gorce, M. Aubert, P. Baratta, R. E. Upham, M. Bonici, C. Carbone, S. Casas, S. Ilić, M. Martinelli, Z. Sakr, A. Schneider, R. Maoli, R. Scaramella, S. Escoffier, W. Gillard, N. Aghanim, A. Amara , et al. (199 additional authors not shown)

Abstract: Deviations from Gaussianity in the distribution of the fields probed by large-scale structure surveys generate additional terms in the data covariance matrix, increasing the uncertainties in the measurement of the cosmological parameters. Super-sample covariance (SSC) is among the largest of these non-Gaussian contributions, with the potential to significantly degrade constraints on some of the pa… ▽ More Deviations from Gaussianity in the distribution of the fields probed by large-scale structure surveys generate additional terms in the data covariance matrix, increasing the uncertainties in the measurement of the cosmological parameters. Super-sample covariance (SSC) is among the largest of these non-Gaussian contributions, with the potential to significantly degrade constraints on some of the parameters of the cosmological model under study -- especially for weak lensing cosmic shear. We compute and validate the impact of SSC on the forecast uncertainties on the cosmological parameters for the Euclid photometric survey, obtained with a Fisher matrix analysis, both considering the Gaussian covariance alone and adding the SSC term -- computed through the public code PySSC. The photometric probes are considered in isolation and combined in the `3$\times$2pt' analysis. We find the SSC impact to be non-negligible -- halving the Figure of Merit of the dark energy parameters ($w_0$, $w_a$) in the 3$\times$2pt case and substantially increasing the uncertainties on $Ω_{{\rm m},0}, w_0$, and $σ_8$ for cosmic shear; photometric galaxy clustering, on the other hand, is less affected due to the lower probe response. The relative impact of SSC does not show significant changes under variations of the redshift binning scheme, while it is smaller for weak lensing when marginalising over the multiplicative shear bias nuisance parameters, which also leads to poorer constraints on the cosmological parameters. Finally, we explore how the use of prior information on the shear and galaxy bias changes the SSC impact. Improving shear bias priors does not have a significant impact, while galaxy bias must be calibrated to sub-percent level to increase the Figure of Merit by the large amount needed to achieve the value when SSC is not included. △ Less

Submitted 24 October, 2023; originally announced October 2023.

Comments: 22 pages, 13 figures

Euclid preparation. XXXI. The effect of the variations in photometric passbands on photometric-redshift accuracy

Authors: Euclid Collaboration, Stéphane Paltani, J. Coupon, W. G. Hartley, A. Alvarez-Ayllon, F. Dubath, J. J. Mohr, M. Schirmer, J. -C. Cuillandre, G. Desprez, O. Ilbert, K. Kuijken, N. Aghanim, B. Altieri, A. Amara, N. Auricchio, M. Baldi, R. Bender, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia, J. Brinchmann, S. Camera, V. Capobianco , et al. (192 additional authors not shown)

Abstract: The technique of photometric redshifts has become essential for the exploitation of multi-band extragalactic surveys. While the requirements on photo-zs for the study of galaxy evolution mostly pertain to the precision and to the fraction of outliers, the most stringent requirement in their use in cosmology is on the accuracy, with a level of bias at the sub-percent level for the Euclid cosmology… ▽ More The technique of photometric redshifts has become essential for the exploitation of multi-band extragalactic surveys. While the requirements on photo-zs for the study of galaxy evolution mostly pertain to the precision and to the fraction of outliers, the most stringent requirement in their use in cosmology is on the accuracy, with a level of bias at the sub-percent level for the Euclid cosmology mission. A separate, and challenging, calibration process is needed to control the bias at this level of accuracy. The bias in photo-zs has several distinct origins that may not always be easily overcome. We identify here one source of bias linked to the spatial or time variability of the passbands used to determine the photometric colours of galaxies. We first quantified the effect as observed on several well-known photometric cameras, and found in particular that, due to the properties of optical filters, the redshifts of off-axis sources are usually overestimated. We show using simple simulations that the detailed and complex changes in the shape can be mostly ignored and that it is sufficient to know the mean wavelength of the passbands of each photometric observation to correct almost exactly for this bias; the key point is that this mean wavelength is independent of the spectral energy distribution of the source}. We use this property to propose a correction that can be computationally efficiently implemented in some photo-z algorithms, in particular template-fitting. We verified that our algorithm, implemented in the new photo-z code Phosphoros, can effectively reduce the bias in photo-zs on real data using the CFHTLS T007 survey, with an average measured bias Delta z over the redshift range 0.4<z<0.7 decreasing by about 0.02, specifically from Delta z~0.04 to Delta z~0.02 around z=0.5. Our algorithm is also able to produce corrected photometry for other applications. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 19 pages, 13 figures; Accepted for publication in A&A

Journal ref: A&A 681, A66 (2024)

Cosmology with Galaxy Photometry Alone

Authors: ChangHoon Hahn, Francisco Villaescusa-Navarro, Peter Melchior, Romain Teyssier

Abstract: We present the first cosmological constraints using only the observed photometry of galaxies. Villaescusa-Navarro et al. (2022; arXiv:2201.02202) recently demonstrated that the internal physical properties of a single simulated galaxy contain a significant amount of cosmological information. These physical properties, however, cannot be directly measured from observations. In this work, we present… ▽ More We present the first cosmological constraints using only the observed photometry of galaxies. Villaescusa-Navarro et al. (2022; arXiv:2201.02202) recently demonstrated that the internal physical properties of a single simulated galaxy contain a significant amount of cosmological information. These physical properties, however, cannot be directly measured from observations. In this work, we present how we can go beyond theoretical demonstrations to infer cosmological constraints from actual galaxy observables (e.g. optical photometry) using neural density estimation and the CAMELS suite of hydrodynamical simulations. We find that the cosmological information in the photometry of a single galaxy is limited. However, we combine the constraining power of photometry from many galaxies using hierarchical population inference and place significant cosmological constraints. With the observed photometry of $\sim$20,000 NASA-Sloan Atlas galaxies, we constrain $Ω_m = 0.323^{+0.075}_{-0.095}$ and $σ_8 = 0.799^{+0.088}_{-0.085}$. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 15 pages, 7 figures, submitted to ApJL, comments welcome

Euclid preparation. XXXIV. The effect of linear redshift-space distortions in photometric galaxy clustering and its cross-correlation with cosmic shear

Authors: Euclid Collaboration, K. Tanidis, V. F. Cardone, M. Martinelli, I. Tutusaus, S. Camera, N. Aghanim, A. Amara, S. Andreon, N. Auricchio, M. Baldi, S. Bardelli, E. Branchini, M. Brescia, J. Brinchmann, V. Capobianco, C. Carbone, J. Carretero, S. Casas, M. Castellano, S. Cavuoti, A. Cimatti, R. Cledassou, G. Congedo, L. Conversi , et al. (185 additional authors not shown)

Abstract: The cosmological surveys that are planned for the current decade will provide us with unparalleled observations of the distribution of galaxies on cosmic scales, by means of which we can probe the underlying large-scale structure (LSS) of the Universe. This will allow us to test the concordance cosmological model and its extensions. However, precision pushes us to high levels of accuracy in the th… ▽ More The cosmological surveys that are planned for the current decade will provide us with unparalleled observations of the distribution of galaxies on cosmic scales, by means of which we can probe the underlying large-scale structure (LSS) of the Universe. This will allow us to test the concordance cosmological model and its extensions. However, precision pushes us to high levels of accuracy in the theoretical modelling of the LSS observables, so that no biases are introduced into the estimation of the cosmological parameters. In particular, effects such as redshift-space distortions (RSD) can become relevant in the computation of harmonic-space power spectra even for the clustering of the photometrically selected galaxies, as has previously been shown in literature. In this work, we investigate the contribution of linear RSD, as formulated in the Limber approximation by a previous work, in forecast cosmological analyses with the photometric galaxy sample of the Euclid survey. We aim to assess their impact and to quantify the bias on the measurement of cosmological parameters that would be caused if this effect were neglected. We performed this task by producing mock power spectra for photometric galaxy clustering and weak lensing, as is expected to be obtained from the Euclid survey. We then used a Markov chain Monte Carlo approach to obtain the posterior distributions of cosmological parameters from these simulated observations. When the linear RSD is neglected, significant biases are caused when galaxy correlations are used alone and when they are combined with cosmic shear in the so-called 3$\times$2pt approach. These biases can be equivalent to as much as $5\,σ$ when an underlying $Λ$CDM cosmology is assumed. When the cosmological model is extended to include the equation-of-state parameters of dark energy, the extension parameters can be shifted by more than $1\,σ$. △ Less

Submitted 22 April, 2024; v1 submitted 31 August, 2023; originally announced September 2023.

Comments: 15 pages, 6 figures. Version matching publication at journal

Euclid Preparation XXXIII. Characterization of convolutional neural networks for the identification of galaxy-galaxy strong lensing events

Authors: Euclid Collaboration, L. Leuzzi, M. Meneghetti, G. Angora, R. B. Metcalf, L. Moscardini, P. Rosati, P. Bergamini, F. Calura, B. Clément, R. Gavazzi, F. Gentile, M. Lochner, C. Grillo, G. Vernardos, N. Aghanim, A. Amara, L. Amendola, S. Andreon, N. Auricchio, S. Bardelli, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia , et al. (194 additional authors not shown)

Abstract: Forthcoming imaging surveys will potentially increase the number of known galaxy-scale strong lenses by several orders of magnitude. For this to happen, images of tens of millions of galaxies will have to be inspected to identify potential candidates. In this context, deep learning techniques are particularly suitable for the finding patterns in large data sets, and convolutional neural networks (… ▽ More Forthcoming imaging surveys will potentially increase the number of known galaxy-scale strong lenses by several orders of magnitude. For this to happen, images of tens of millions of galaxies will have to be inspected to identify potential candidates. In this context, deep learning techniques are particularly suitable for the finding patterns in large data sets, and convolutional neural networks (CNNs) in particular can efficiently process large volumes of images. We assess and compare the performance of three network architectures in the classification of strong lensing systems on the basis of their morphological characteristics. We train and test our models on different subsamples of a data set of forty thousand mock images, having characteristics similar to those expected in the wide survey planned with the ESA mission \Euclid, gradually including larger fractions of faint lenses. We also evaluate the importance of adding information about the colour difference between the lens and source galaxies by repeating the same training on single-band and multi-band images. Our models find samples of clear lenses with $\gtrsim 90\%$ precision and completeness, without significant differences in the performance of the three architectures. Nevertheless, when including lenses with fainter arcs in the training set, the three models' performance deteriorates with accuracy values of $\sim 0.87$ to $\sim 0.75$ depending on the model. Our analysis confirms the potential of the application of CNNs to the identification of galaxy-scale strong lenses. We suggest that specific training with separate classes of lenses might be needed for detecting the faint lenses since the addition of the colour information does not yield a significant improvement in the current analysis, with the accuracy ranging from $\sim 0.89$ to $\sim 0.78$ for the different models. △ Less

Submitted 26 January, 2024; v1 submitted 17 July, 2023; originally announced July 2023.

Comments: 23 pages,12 figures

Journal ref: Astronomy & Astrophysics,2024, 681, A68

Euclid preparation. XXIX. Water ice in spacecraft part I: The physics of ice formation and contamination

Authors: Euclid Collaboration, M. Schirmer, K. Thürmer, B. Bras, M. Cropper, J. Martin-Fleitas, Y. Goueffon, R. Kohley, A. Mora, M. Portaluppi, G. D. Racca, A. D. Short, S. Szmolka, L. M. Gaspar Venancio, M. Altmann, Z. Balog, U. Bastian, M. Biermann, D. Busonero, C. Fabricius, F. Grupp, C. Jordi, W. Löffler, A. Sagristà Sellés, N. Aghanim , et al. (196 additional authors not shown)

Abstract: Molecular contamination is a well-known problem in space flight. Water is the most common contaminant and alters numerous properties of a cryogenic optical system. Too much ice means that Euclid's calibration requirements and science goals cannot be met. Euclid must then be thermally decontaminated, a long and risky process. We need to understand how iced optics affect the data and when a decontam… ▽ More Molecular contamination is a well-known problem in space flight. Water is the most common contaminant and alters numerous properties of a cryogenic optical system. Too much ice means that Euclid's calibration requirements and science goals cannot be met. Euclid must then be thermally decontaminated, a long and risky process. We need to understand how iced optics affect the data and when a decontamination is required. This is essential to build adequate calibration and survey plans, yet a comprehensive analysis in the context of an astrophysical space survey has not been done before. In this paper we look at other spacecraft with well-documented outgassing records, and we review the formation of thin ice films. A mix of amorphous and crystalline ices is expected for Euclid. Their surface topography depends on the competing energetic needs of the substrate-water and the water-water interfaces, and is hard to predict with current theories. We illustrate that with scanning-tunnelling and atomic-force microscope images. Industrial tools exist to estimate contamination, and we must understand their uncertainties. We find considerable knowledge errors on the diffusion and sublimation coefficients, limiting the accuracy of these tools. We developed a water transport model to compute contamination rates in Euclid, and find general agreement with industry estimates. Tests of the Euclid flight hardware in space simulators did not pick up contamination signals; our in-flight calibrations observations will be much more sensitive. We must understand the link between the amount of ice on the optics and its effect on Euclid's data. Little research is available about this link, possibly because other spacecraft can decontaminate easily, quenching the need for a deeper understanding. In our second paper we quantify the various effects of iced optics on spectrophotometric data. △ Less

Submitted 23 May, 2023; v1 submitted 17 May, 2023; originally announced May 2023.

Comments: 35 pages, 22 figures, A&A in press. Changes to previous version: language edits, added Z. Bolag as author in the arxiv PDF (was listed in the ASCII author list and in the journal PDF, but not in the arxiv PDF). This version is identical to the journal version

Journal ref: A&A 675, A142 (2023)

Cosmology with one galaxy? -- The ASTRID model and robustness

Authors: Nicolas Echeverri, Francisco Villaescusa-Navarro, Chaitanya Chawak, Yueying Ni, ChangHoon Hahn, Elena Hernandez-Martinez, Romain Teyssier, Daniel Angles-Alcazar, Klaus Dolag, Tiago Castro

Abstract: Recent work has pointed out the potential existence of a tight relation between the cosmological parameter $Ω_{\rm m}$, at fixed $Ω_{\rm b}$, and the properties of individual galaxies in state-of-the-art cosmological hydrodynamic simulations. In this paper, we investigate whether such a relation also holds for galaxies from simulations run with a different code that made use of a distinct subgrid… ▽ More Recent work has pointed out the potential existence of a tight relation between the cosmological parameter $Ω_{\rm m}$, at fixed $Ω_{\rm b}$, and the properties of individual galaxies in state-of-the-art cosmological hydrodynamic simulations. In this paper, we investigate whether such a relation also holds for galaxies from simulations run with a different code that made use of a distinct subgrid physics: Astrid. We find that also in this case, neural networks are able to infer the value of $Ω_{\rm m}$ with a $\sim10\%$ precision from the properties of individual galaxies while accounting for astrophysics uncertainties as modeled in CAMELS. This tight relationship is present at all considered redshifts, $z\leq3$, and the stellar mass, the stellar metallicity, and the maximum circular velocity are among the most important galaxy properties behind the relation. In order to use this method with real galaxies, one needs to quantify its robustness: the accuracy of the model when tested on galaxies generated by codes different from the one used for training. We quantify the robustness of the models by testing them on galaxies from four different codes: IllustrisTNG, SIMBA, Astrid, and Magneticum. We show that the models perform well on a large fraction of the galaxies, but fail dramatically on a small fraction of them. Removing these outliers significantly improves the accuracy of the models across simulation codes. △ Less

Submitted 12 April, 2023; originally announced April 2023.

Comments: 16 pages, 12 figures

A universal equation to predict $Ω_{\rm m}$ from halo and galaxy catalogues

Authors: Helen Shao, Natalí S. M de Santi, Francisco Villaescusa-Navarro, Romain Teyssier, Yueying Ni, Daniel Angles-Alcazar, Shy Genel, Lars Hernquist, Ulrich P. Steinwandel, Tiago Castro, Elena Hernandez-Martınez, Klaus Dolag, Christopher C. Lovell, Eli Visbal, Lehman H. Garrison, Mihir Kulkarni

Abstract: We discover analytic equations that can infer the value of $Ω_{\rm m}$ from the positions and velocity moduli of halo and galaxy catalogues. The equations are derived by combining a tailored graph neural network (GNN) architecture with symbolic regression. We first train the GNN on dark matter halos from Gadget N-body simulations to perform field-level likelihood-free inference, and show that our… ▽ More We discover analytic equations that can infer the value of $Ω_{\rm m}$ from the positions and velocity moduli of halo and galaxy catalogues. The equations are derived by combining a tailored graph neural network (GNN) architecture with symbolic regression. We first train the GNN on dark matter halos from Gadget N-body simulations to perform field-level likelihood-free inference, and show that our model can infer $Ω_{\rm m}$ with $\sim6\%$ accuracy from halo catalogues of thousands of N-body simulations run with six different codes: Abacus, CUBEP$^3$M, Gadget, Enzo, PKDGrav3, and Ramses. By applying symbolic regression to the different parts comprising the GNN, we derive equations that can predict $Ω_{\rm m}$ from halo catalogues of simulations run with all of the above codes with accuracies similar to those of the GNN. We show that by tuning a single free parameter, our equations can also infer the value of $Ω_{\rm m}$ from galaxy catalogues of thousands of state-of-the-art hydrodynamic simulations of the CAMELS project, each with a different astrophysics model, run with five distinct codes that employ different subgrid physics: IllustrisTNG, SIMBA, Astrid, Magneticum, SWIFT-EAGLE. Furthermore, the equations also perform well when tested on galaxy catalogues from simulations covering a vast region in parameter space that samples variations in 5 cosmological and 23 astrophysical parameters. We speculate that the equations may reflect the existence of a fundamental physics relation between the phase-space distribution of generic tracers and $Ω_{\rm m}$, one that is not affected by galaxy formation physics down to scales as small as $10~h^{-1}{\rm kpc}$. △ Less

Submitted 28 February, 2023; originally announced February 2023.

Comments: 32 pages, 13 figures, summary video: https://youtu.be/STZHvDHkVgo

Robust Field-level Likelihood-free Inference with Galaxies

Authors: Natalí S. M. de Santi, Helen Shao, Francisco Villaescusa-Navarro, L. Raul Abramo, Romain Teyssier, Pablo Villanueva-Domingo, Yueying Ni, Daniel Anglés-Alcázar, Shy Genel, Elena Hernandez-Martinez, Ulrich P. Steinwandel, Christopher C. Lovell, Klaus Dolag, Tiago Castro, Mark Vogelsberger

Abstract: We train graph neural networks to perform field-level likelihood-free inference using galaxy catalogs from state-of-the-art hydrodynamic simulations of the CAMELS project. Our models are rotational, translational, and permutation invariant and do not impose any cut on scale. From galaxy catalogs that only contain $3$D positions and radial velocities of $\sim 1, 000$ galaxies in tiny… ▽ More We train graph neural networks to perform field-level likelihood-free inference using galaxy catalogs from state-of-the-art hydrodynamic simulations of the CAMELS project. Our models are rotational, translational, and permutation invariant and do not impose any cut on scale. From galaxy catalogs that only contain $3$D positions and radial velocities of $\sim 1, 000$ galaxies in tiny $(25~h^{-1}{\rm Mpc})^3$ volumes our models can infer the value of $Ω_{\rm m}$ with approximately $12$ % precision. More importantly, by testing the models on galaxy catalogs from thousands of hydrodynamic simulations, each having a different efficiency of supernova and AGN feedback, run with five different codes and subgrid models - IllustrisTNG, SIMBA, Astrid, Magneticum, SWIFT-EAGLE -, we find that our models are robust to changes in astrophysics, subgrid physics, and subhalo/galaxy finder. Furthermore, we test our models on $1,024$ simulations that cover a vast region in parameter space - variations in $5$ cosmological and $23$ astrophysical parameters - finding that the model extrapolates really well. Our results indicate that the key to building a robust model is the use of both galaxy positions and velocities, suggesting that the network have likely learned an underlying physical relation that does not depend on galaxy formation and is valid on scales larger than $\sim10~h^{-1}{\rm kpc}$. △ Less

Submitted 18 July, 2023; v1 submitted 27 February, 2023; originally announced February 2023.

Comments: 34 pages, 12 figures. For a video summarizing the results, see https://youtu.be/b59ep7cyPOs

Journal ref: Volume 952, Number 1, Year 2023, Pages 69

Euclid preparation. XXX. Performance assessment of the NISP Red-Grism through spectroscopic simulations for the Wide and Deep surveys

Authors: Euclid Collaboration, L. Gabarra, C. Mancini, L. Rodriguez Munoz, G. Rodighiero, C. Sirignano, M. Scodeggio, M. Talia, S. Dusini, W. Gillard, B. R. Granett, E. Maiorano, M. Moresco, L. Paganin, E. Palazzi, L. Pozzetti, A. Renzi, E. Rossetti, D. Vergani, V. Allevato, L. Bisigello, G. Castignani, B. De Caro, M. Fumana, K. Ganga , et al. (210 additional authors not shown)

Abstract: This work focuses on the pilot run of a simulation campaign aimed at investigating the spectroscopic capabilities of the Euclid Near-Infrared Spectrometer and Photometer (NISP), in terms of continuum and emission line detection in the context of galaxy evolutionary studies. To this purpose we constructed, emulated, and analysed the spectra of 4992 star-forming galaxies at $0.3 \leq z \leq 2.5$ usi… ▽ More This work focuses on the pilot run of a simulation campaign aimed at investigating the spectroscopic capabilities of the Euclid Near-Infrared Spectrometer and Photometer (NISP), in terms of continuum and emission line detection in the context of galaxy evolutionary studies. To this purpose we constructed, emulated, and analysed the spectra of 4992 star-forming galaxies at $0.3 \leq z \leq 2.5$ using the NISP pixel-level simulator. We built the spectral library starting from public multi-wavelength galaxy catalogues, with value-added information on spectral energy distribution (SED) fitting results, and from Bruzual and Charlot (2003) stellar population templates. Rest-frame optical and near-IR nebular emission lines were included using empirical and theoretical relations. We inferred the 3.5$σ$ NISP red grism spectroscopic detection limit of the continuum measured in the $H$ band for star-forming galaxies with a median disk half-light radius of \ang{;;0.4} at magnitude $H= 19.5\pm0.2\,$AB$\,$mag for the Euclid Wide Survey and at $H = 20.8\pm0.6\,$AB$\,$mag for the Euclid Deep Survey. We found a very good agreement with the red grism emission line detection limit requirement for the Wide and Deep surveys. We characterised the effect of the galaxy shape on the detection capability of the red grism and highlighted the degradation of the quality of the extracted spectra as the disk size increases. In particular, we found that the extracted emission line signal to noise ratio (SNR) drops by $\sim\,$45$\%$ when the disk size ranges from \ang{;;0.25} to \ang{;;1}. These trends lead to a correlation between the emission line SNR and the stellar mass of the galaxy and we demonstrate the effect in a stacking analysis unveiling emission lines otherwise too faint to detect. △ Less

Submitted 25 August, 2023; v1 submitted 18 February, 2023; originally announced February 2023.

Comments: 24 pages, 21 figures

Journal ref: Astronomy and Astrophysics (2023)

Cosmic Variance and the Inhomogeneous UV Luminosity Function of Galaxies During Reionization

Authors: Taha Dawoodbhoy, Paul R. Shapiro, Pierre Ocvirk, Joseph S. W. Lewis, Dominique Aubert, Jenny G. Sorce, Kyungjin Ahn, Ilian T. Iliev, Hyunbae Park, Romain Teyssier, Gustavo Yepes

Abstract: When the first galaxies formed and starlight escaped into the intergalactic medium to reionize it, galaxy formation and reionization were both highly inhomogeneous in time and space, and fully-coupled by mutual feedback. To show how this imprinted the UV luminosity function (UVLF) of reionization-era galaxies, we use our large-scale, radiation-hydrodynamics simulation CoDa II to derive the time- a… ▽ More When the first galaxies formed and starlight escaped into the intergalactic medium to reionize it, galaxy formation and reionization were both highly inhomogeneous in time and space, and fully-coupled by mutual feedback. To show how this imprinted the UV luminosity function (UVLF) of reionization-era galaxies, we use our large-scale, radiation-hydrodynamics simulation CoDa II to derive the time- and space-varying halo mass function and UVLF, from $z\simeq6$-15. That UVLF correlates strongly with local reionization redshift: earlier-reionizing regions have UVLFs that are higher, more extended to brighter magnitudes, and flatter at the faint end than later-reionizing regions observed at the same $z$. In general, as a region reionizes, the faint-end slope of its local UVLF flattens, and, by $z=6$ (when reionization ended), the global UVLF, too, exhibits a flattened faint-end slope, `rolling-over' at $M_\text{UV}\gtrsim-17$. CoDa II's UVLF is broadly consistent with cluster-lensed galaxy observations of the Hubble Frontier Fields at $z=6$-8, including the faint end, except for the faintest data point at $z=6$, based on one galaxy at $M_\text{UV}=-12.5$. According to CoDa II, the probability of observing the latter is $\sim5\%$. However, the effective volume searched at this magnitude is very small, and is thus subject to significant cosmic variance. We find that previous methods adopted to calculate the uncertainty due to cosmic variance underestimated it on such small scales by a factor of 2-4, primarily by underestimating the variance in halo abundance when the sample volume is small. △ Less

Submitted 8 August, 2023; v1 submitted 16 February, 2023; originally announced February 2023.

Comments: 16 pages, 16 figures, accepted by MNRAS 07/20/23, comments welcome

Euclid preparation: XXVIII. Modelling of the weak lensing angular power spectrum

Authors: Euclid Collaboration, A. C. Deshpande, T. Kitching, A. Hall, M. L. Brown, N. Aghanim, L. Amendola, N. Auricchio, M. Baldi, R. Bender, D. Bonino, E. Branchini, M. Brescia, J. Brinchmann, S. Camera, G. P. Candini, V. Capobianco, C. Carbone, V. F. Cardone, J. Carretero, F. J. Castander, M. Castellano, S. Cavuoti, A. Cimatti, R. Cledassou , et al. (178 additional authors not shown)

Abstract: This work considers which higher-order effects in modelling the cosmic shear angular power spectra must be taken into account for Euclid. We identify which terms are of concern, and quantify their individual and cumulative impact on cosmological parameter inference from Euclid. We compute the values of these higher-order effects using analytic expressions, and calculate the impact on cosmological… ▽ More This work considers which higher-order effects in modelling the cosmic shear angular power spectra must be taken into account for Euclid. We identify which terms are of concern, and quantify their individual and cumulative impact on cosmological parameter inference from Euclid. We compute the values of these higher-order effects using analytic expressions, and calculate the impact on cosmological parameter estimation using the Fisher matrix formalism. We review 24 effects and find the following potentially need to be accounted for: the reduced shear approximation, magnification bias, source-lens clustering, source obscuration, local Universe effects, and the flat Universe assumption. Upon computing these explicitly, and calculating their cosmological parameter biases, using a maximum multipole of $\ell=5000$, we find that the magnification bias, source-lens clustering, source obscuration, and local Universe terms individually produce significant ($\,>0.25σ$) cosmological biases in one or more parameters, and accordingly must be accounted for. In total, over all effects, we find biases in $Ω_{\rm m}$, $Ω_{\rm b}$, $h$, and $σ_{8}$ of $0.73σ$, $0.28σ$, $0.25σ$, and $-0.79σ$, respectively, for flat $Λ$CDM. For the $w_0w_a$CDM case, we find biases in $Ω_{\rm m}$, $Ω_{\rm b}$, $h$, $n_{\rm s}$, $σ_{8}$, and $w_a$ of $1.49σ$, $0.35σ$, $-1.36σ$, $1.31σ$, $-0.84σ$, and $-0.35σ$, respectively; which are increased relative to the $Λ$CDM due to additional degeneracies as a function of redshift and scale. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: 20 pages, submitted to A&A

Euclid Preparation. XXVIII. Forecasts for ten different higher-order weak lensing statistics

Authors: Euclid Collaboration, V. Ajani, M. Baldi, A. Barthelemy, A. Boyle, P. Burger, V. F. Cardone, S. Cheng, S. Codis, C. Giocoli, J. Harnois-Déraps, S. Heydenreich, V. Kansal, M. Kilbinger, L. Linke, C. Llinares, N. Martinet, C. Parroni, A. Peel, S. Pires, L. Porth, I. Tereno, C. Uhlemann, M. Vicinanza, S. Vinciguerra , et al. (189 additional authors not shown)

Abstract: Recent cosmic shear studies have shown that higher-order statistics (HOS) developed by independent teams now outperform standard two-point estimators in terms of statistical precision thanks to their sensitivity to the non-Gaussian features of large-scale structure. The aim of the Higher-Order Weak Lensing Statistics (HOWLS) project is to assess, compare, and combine the constraining power of ten… ▽ More Recent cosmic shear studies have shown that higher-order statistics (HOS) developed by independent teams now outperform standard two-point estimators in terms of statistical precision thanks to their sensitivity to the non-Gaussian features of large-scale structure. The aim of the Higher-Order Weak Lensing Statistics (HOWLS) project is to assess, compare, and combine the constraining power of ten different HOS on a common set of $Euclid$-like mocks, derived from N-body simulations. In this first paper of the HOWLS series, we computed the nontomographic ($Ω_{\rm m}$, $σ_8$) Fisher information for the one-point probability distribution function, peak counts, Minkowski functionals, Betti numbers, persistent homology Betti numbers and heatmap, and scattering transform coefficients, and we compare them to the shear and convergence two-point correlation functions in the absence of any systematic bias. We also include forecasts for three implementations of higher-order moments, but these cannot be robustly interpreted as the Gaussian likelihood assumption breaks down for these statistics. Taken individually, we find that each HOS outperforms the two-point statistics by a factor of around two in the precision of the forecasts with some variations across statistics and cosmological parameters. When combining all the HOS, this increases to a $4.5$ times improvement, highlighting the immense potential of HOS for cosmic shear cosmological analyses with $Euclid$. The data used in this analysis are publicly released with the paper. △ Less

Submitted 10 July, 2023; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: 33 pages, 24 figures, main results in Fig. 19 & Table 5, version published in A&A

Journal ref: A&A 675, A120 (2023)

Cosmological simulations of the same spiral galaxy: connecting the dark matter distribution of the host halo with the subgrid baryonic physics

Authors: A. Núñez-Castiñeyra, E. Nezri, P. Mollitor, J. Devriendt, R. Teyssier

Abstract: The role of baryonic physics, star formation, and stellar feedback, in shaping the galaxies and their host halos is an evolving topic. The dark matter aspects are illustrated in this work by showing distribution features in a Milky-Way-sized halo. We focus on the halo morphology, geometry, and profile as well as the phase space distribution using one dark matter only and five hydrodynamical cosmol… ▽ More The role of baryonic physics, star formation, and stellar feedback, in shaping the galaxies and their host halos is an evolving topic. The dark matter aspects are illustrated in this work by showing distribution features in a Milky-Way-sized halo. We focus on the halo morphology, geometry, and profile as well as the phase space distribution using one dark matter only and five hydrodynamical cosmological high-resolution simulations of the same halo with different subgrid prescriptions for the baryonic physics (Kennicut versus multi-freefall star formation and delayed cooling versus mechanical supernovae feedback). If some general properties like the relative halo-galaxy orientation are similar, the modifications of the gravitational potential due to the presence of baryons are found to induce different dark matter distributions (rounder and more concentrated halo). The mass density profile as well as the velocity distribution are modified distinctively according to the specific resulting baryonic distribution highlighting the variability of those properties (e.g inner power index from 1.3 to 1.8, broader speed distribution). The uncertainties on those features are of paramount importance for dark matter phenomenology, particularly when dealing with dark matter dynamics or direct and indirect detection searches. As a consequence, dark matter properties and prospects using cosmological simulations require improvement on baryonic physics description. Modeling such processes is a key issue not only for galaxy formation but also for dark matter investigations. △ Less

Submitted 30 March, 2023; v1 submitted 15 January, 2023; originally announced January 2023.

Euclid preparation: XXII. Selection of Quiescent Galaxies from Mock Photometry using Machine Learning

Authors: Euclid Collaboration, A. Humphrey, L. Bisigello, P. A. C. Cunha, M. Bolzonella, S. Fotopoulou, K. Caputi, C. Tortora, G. Zamorani, P. Papaderos, D. Vergani, J. Brinchmann, M. Moresco, A. Amara, N. Auricchio, M. Baldi, R. Bender, D. Bonino, E. Branchini, M. Brescia, S. Camera, V. Capobianco, C. Carbone, J. Carretero, F. J. Castander , et al. (184 additional authors not shown)

Abstract: The Euclid Space Telescope will provide deep imaging at optical and near-infrared wavelengths, along with slitless near-infrared spectroscopy, across ~15,000 sq deg of the sky. Euclid is expected to detect ~12 billion astronomical sources, facilitating new insights into cosmology, galaxy evolution, and various other topics. To optimally exploit the expected very large data set, there is the need t… ▽ More The Euclid Space Telescope will provide deep imaging at optical and near-infrared wavelengths, along with slitless near-infrared spectroscopy, across ~15,000 sq deg of the sky. Euclid is expected to detect ~12 billion astronomical sources, facilitating new insights into cosmology, galaxy evolution, and various other topics. To optimally exploit the expected very large data set, there is the need to develop appropriate methods and software. Here we present a novel machine-learning based methodology for selection of quiescent galaxies using broad-band Euclid I_E, Y_E, J_E, H_E photometry, in combination with multiwavelength photometry from other surveys. The ARIADNE pipeline uses meta-learning to fuse decision-tree ensembles, nearest-neighbours, and deep-learning methods into a single classifier that yields significantly higher accuracy than any of the individual learning methods separately. The pipeline has `sparsity-awareness', so that missing photometry values are still informative for the classification. Our pipeline derives photometric redshifts for galaxies selected as quiescent, aided by the `pseudo-labelling' semi-supervised method. After application of the outlier filter, our pipeline achieves a normalized mean absolute deviation of ~< 0.03 and a fraction of catastrophic outliers of ~< 0.02 when measured against the COSMOS2015 photometric redshifts. We apply our classification pipeline to mock galaxy photometry catalogues corresponding to three main scenarios: (i) Euclid Deep Survey with ancillary ugriz, WISE, and radio data; (ii) Euclid Wide Survey with ancillary ugriz, WISE, and radio data; (iii) Euclid Wide Survey only. Our classification pipeline outperforms UVJ selection, in addition to the Euclid I_E-Y_E, J_E-H_E and u-I_E,I_E-J_E colour-colour methods, with improvements in completeness and the F1-score of up to a factor of 2. (Abridged) △ Less

Submitted 5 December, 2022; v1 submitted 26 September, 2022; originally announced September 2022.

Comments: 37 pages (including appendices), 26 figures; accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 671, A99 (2023)

Euclid preparation XXVI. The Euclid Morphology Challenge. Towards structural parameters for billions of galaxies

Authors: Euclid Collaboration, H. Bretonnière, U. Kuchner, M. Huertas-Company, E. Merlin, M. Castellano, D. Tuccillo, F. Buitrago, C. J. Conselice, A. Boucaud, B. Häußler, M. Kümmel, W. G. Hartley, A. Alvarez Ayllon, E. Bertin, F. Ferrari, L. Ferreira, R. Gavazzi, D. Hernández-Lang, G. Lucatelli, A. S. G. Robotham, M. Schefer, L. Wang, R. Cabanac, H. Domínguez Sánchez , et al. (193 additional authors not shown)

Abstract: The various Euclid imaging surveys will become a reference for studies of galaxy morphology by delivering imaging over an unprecedented area of 15 000 square degrees with high spatial resolution. In order to understand the capabilities of measuring morphologies from Euclid-detected galaxies and to help implement measurements in the pipeline, we have conducted the Euclid Morphology Challenge, which… ▽ More The various Euclid imaging surveys will become a reference for studies of galaxy morphology by delivering imaging over an unprecedented area of 15 000 square degrees with high spatial resolution. In order to understand the capabilities of measuring morphologies from Euclid-detected galaxies and to help implement measurements in the pipeline, we have conducted the Euclid Morphology Challenge, which we present in two papers. While the companion paper by Merlin et al. focuses on the analysis of photometry, this paper assesses the accuracy of the parametric galaxy morphology measurements in imaging predicted from within the Euclid Wide Survey. We evaluate the performance of five state-of-the-art surface-brightness-fitting codes DeepLeGATo, Galapagos-2, Morfometryka, Profit and SourceXtractor++ on a sample of about 1.5 million simulated galaxies resembling reduced observations with the Euclid VIS and NIR instruments. The simulations include analytic Sérsic profiles with one and two components, as well as more realistic galaxies generated with neural networks. We find that, despite some code-specific differences, all methods tend to achieve reliable structural measurements (10% scatter on ideal Sérsic simulations) down to an apparent magnitude of about 23 in one component and 21 in two components, which correspond to a signal-to-noise ratio of approximately 1 and 5 respectively. We also show that when tested on non-analytic profiles, the results are typically degraded by a factor of 3, driven by systematics. We conclude that the Euclid official Data Releases will deliver robust structural parameters for at least 400 million galaxies in the Euclid Wide Survey by the end of the mission. We find that a key factor for explaining the different behaviour of the codes at the faint end is the set of adopted priors for the various structural parameters. △ Less

Submitted 28 November, 2022; v1 submitted 26 September, 2022; originally announced September 2022.

Comments: Accepted by A&A. 30 pages, 23+6 figures, Euclid pre-launch key paper. Companion paper: Euclid Collaboration XXV: Merlin et al. 2022 Minor corrections after journal review

Journal ref: A&A 671, A102 (2023)

Euclid preparation. XXV. The Euclid Morphology Challenge -- Towards model-fitting photometry for billions of galaxies

Authors: Euclid Collaboration, E. Merlin, M. Castellano, H. Bretonnière, M. Huertas-Company, U. Kuchner, D. Tuccillo, F. Buitrago, J. R. Peterson, C. J. Conselice, F. Caro, P. Dimauro, L. Nemani, A. Fontana, M. Kümmel, B. Häußler, W. G. Hartley, A. Alvarez Ayllon, E. Bertin, P. Dubath, F. Ferrari, L. Ferreira, R. Gavazzi, D. Hernández-Lang, G. Lucatelli , et al. (196 additional authors not shown)

Abstract: The ESA Euclid mission will provide high-quality imaging for about 1.5 billion galaxies. A software pipeline to automatically process and analyse such a huge amount of data in real time is being developed by the Science Ground Segment of the Euclid Consortium; this pipeline will include a model-fitting algorithm, which will provide photometric and morphological estimates of paramount importance fo… ▽ More The ESA Euclid mission will provide high-quality imaging for about 1.5 billion galaxies. A software pipeline to automatically process and analyse such a huge amount of data in real time is being developed by the Science Ground Segment of the Euclid Consortium; this pipeline will include a model-fitting algorithm, which will provide photometric and morphological estimates of paramount importance for the core science goals of the mission and for legacy science. The Euclid Morphology Challenge is a comparative investigation of the performance of five model-fitting software packages on simulated Euclid data, aimed at providing the baseline to identify the best suited algorithm to be implemented in the pipeline. In this paper we describe the simulated data set, and we discuss the photometry results. A companion paper (Euclid Collaboration: Bretonnière et al. 2022) is focused on the structural and morphological estimates. We created mock Euclid images simulating five fields of view of 0.48 deg2 each in the $I_E$ band of the VIS instrument, each with three realisations of galaxy profiles (single and double Sérsic, and 'realistic' profiles obtained with a neural network); for one of the fields in the double Sérsic realisation, we also simulated images for the three near-infrared $Y_E$, $J_E$ and $H_E$ bands of the NISP-P instrument, and five Rubin/LSST optical complementary bands ($u$, $g$, $r$, $i$, and $z$). To analyse the results we created diagnostic plots and defined ad-hoc metrics. Five model-fitting software packages (DeepLeGATo, Galapagos-2, Morfometryka, ProFit, and SourceXtractor++) were compared, all typically providing good results. (cut) △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: 29 pages, 33 figures. Euclid pre-launch key paper. Companion paper: Bretonniere et al. 2022

Journal ref: A&A 671, A101 (2023)

Robust field-level inference with dark matter halos

Authors: Helen Shao, Francisco Villaescusa-Navarro, Pablo Villanueva-Domingo, Romain Teyssier, Lehman H. Garrison, Marco Gatti, Derek Inman, Yueying Ni, Ulrich P. Steinwandel, Mihir Kulkarni, Eli Visbal, Greg L. Bryan, Daniel Angles-Alcazar, Tiago Castro, Elena Hernandez-Martinez, Klaus Dolag

Abstract: We train graph neural networks on halo catalogues from Gadget N-body simulations to perform field-level likelihood-free inference of cosmological parameters. The catalogues contain $\lesssim$5,000 halos with masses $\gtrsim 10^{10}~h^{-1}M_\odot$ in a periodic volume of $(25~h^{-1}{\rm Mpc})^3$; every halo in the catalogue is characterized by several properties such as position, mass, velocity, co… ▽ More We train graph neural networks on halo catalogues from Gadget N-body simulations to perform field-level likelihood-free inference of cosmological parameters. The catalogues contain $\lesssim$5,000 halos with masses $\gtrsim 10^{10}~h^{-1}M_\odot$ in a periodic volume of $(25~h^{-1}{\rm Mpc})^3$; every halo in the catalogue is characterized by several properties such as position, mass, velocity, concentration, and maximum circular velocity. Our models, built to be permutationally, translationally, and rotationally invariant, do not impose a minimum scale on which to extract information and are able to infer the values of $Ω_{\rm m}$ and $σ_8$ with a mean relative error of $\sim6\%$, when using positions plus velocities and positions plus masses, respectively. More importantly, we find that our models are very robust: they can infer the value of $Ω_{\rm m}$ and $σ_8$ when tested using halo catalogues from thousands of N-body simulations run with five different N-body codes: Abacus, CUBEP$^3$M, Enzo, PKDGrav3, and Ramses. Surprisingly, the model trained to infer $Ω_{\rm m}$ also works when tested on thousands of state-of-the-art CAMELS hydrodynamic simulations run with four different codes and subgrid physics implementations. Using halo properties such as concentration and maximum circular velocity allow our models to extract more information, at the expense of breaking the robustness of the models. This may happen because the different N-body codes are not converged on the relevant scales corresponding to these parameters. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: 25 pages, 11 figures, summary video: https://youtu.be/qkw92Z6owJU