Showing 1–15 of 15 results for author: Tarsitano, F

Euclid: Early Release Observations -- A preview of the Euclid era through a galaxy cluster magnifying lens

Authors: H. Atek, R. Gavazzi, J. R. Weaver, J. M. Diego, T. Schrabback, N. A. Hatch, N. Aghanim, H. Dole, W. G. Hartley, S. Taamoli, G. Congedo, Y. Jimenez-Teja, J. -C. Cuillandre, E. Bañados, S. Belladitta, R. A. A. Bowler, M. Franco, M. Jauzac, G. Mahler, J. Richard, P. -F. Rocci, S. Serjeant, S. Toft, D. Abriola, P. Bergamini , et al. (178 additional authors not shown)

Abstract: We present the first analysis of the Euclid Early Release Observations (ERO) program that targets fields around two lensing clusters, Abell 2390 and Abell 2764. We use VIS and NISP imaging to produce photometric catalogs for a total of $\sim 500\,000$ objects. The imaging data reach a $5\,σ$ typical depth in the range 25.1-25.4 AB in the NISP bands, and 27.1-27.3 AB in the VIS band. Using the Lyma… ▽ More We present the first analysis of the Euclid Early Release Observations (ERO) program that targets fields around two lensing clusters, Abell 2390 and Abell 2764. We use VIS and NISP imaging to produce photometric catalogs for a total of $\sim 500\,000$ objects. The imaging data reach a $5\,σ$ typical depth in the range 25.1-25.4 AB in the NISP bands, and 27.1-27.3 AB in the VIS band. Using the Lyman-break method in combination with photometric redshifts, we identify $30$ Lyman-break galaxy (LBG) candidates at $z>6$ and 139 extremely red sources (ERSs), most likely at lower redshift. The deeper VIS imaging compared to NISP means we can routinely identify high-redshift Lyman breaks of the order of $3$ magnitudes, which reduces contamination by brown dwarf stars and low-redshift galaxies. Spectroscopic follow-up campaigns of such bright sources will help constrain both the bright end of the ultraviolet galaxy luminosity function and the quasar luminosity function at $z>6$, and constrain the physical nature of these objects. Additionally, we have performed a combined strong lensing and weak lensing analysis of A2390, and demonstrate how Euclid will contribute to better constraining the virial mass of galaxy clusters. From these data, we also identify optical and near-infrared counterparts of known $z>0.6$ clusters, which exhibit strong lensing features, establishing the ability of Euclid to characterize high-redshift clusters. Finally, we provide a glimpse of Euclid's ability to map the intracluster light out to larger radii than current facilities, enabling a better understanding of the cluster assembly history and mapping of the dark matter distribution. This initial dataset illustrates the diverse spectrum of legacy science that will be enabled by the Euclid survey. △ 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. 17 pages, 12 figures

Euclid: Early Release Observations -- Overview of the Perseus cluster and analysis of its luminosity and stellar mass functions

Authors: J. -C. Cuillandre, M. Bolzonella, A. Boselli, F. R. Marleau, M. Mondelin, J. G. Sorce, C. Stone, F. Buitrago, Michele Cantiello, K. George, N. A. Hatch, L. Quilley, F. Mannucci, T. Saifollahi, R. S��nchez-Janssen, F. Tarsitano, C. Tortora, X. Xu, H. Bouy, S. Gwyn, M. Kluge, A. Lançon, R. Laureijs, M. Schirmer, Abdurro'uf , et al. (177 additional authors not shown)

Abstract: The Euclid ERO programme targeted the Perseus cluster of galaxies, gathering deep data in the central region of the cluster over 0.7 square degree, corresponding to approximately 0.25 r_200. The data set reaches a point-source depth of IE=28.0 (YE, JE, HE = 25.3) AB magnitudes at 5 sigma with a 0.16" and 0.48" FWHM, and a surface brightness limit of 30.1 (29.2) mag per square arcsec. The exception… ▽ More The Euclid ERO programme targeted the Perseus cluster of galaxies, gathering deep data in the central region of the cluster over 0.7 square degree, corresponding to approximately 0.25 r_200. The data set reaches a point-source depth of IE=28.0 (YE, JE, HE = 25.3) AB magnitudes at 5 sigma with a 0.16" and 0.48" FWHM, and a surface brightness limit of 30.1 (29.2) mag per square arcsec. The exceptional depth and spatial resolution of this wide-field multi-band data enable the simultaneous detection and characterisation of both bright and low surface brightness galaxies, along with their globular cluster systems, from the optical to the NIR. This study advances beyond previous analyses of the cluster and enables a range of scientific investigations summarised here. We derive the luminosity and stellar mass functions (LF and SMF) of the Perseus cluster in the Euclid IE band, thanks to supplementary u,g,r,i,z and Halpha data from the CFHT. We adopt a catalogue of 1100 dwarf galaxies, detailed in the corresponding ERO paper. We identify all other sources in the Euclid images and obtain accurate photometric measurements using AutoProf or AstroPhot for 138 bright cluster galaxies, and SourceExtractor for half a million compact sources. Cluster membership for the bright sample is determined by calculating photometric redshifts with Phosphoros. Our LF and SMF are the deepest recorded for the Perseus cluster, highlighting the groundbreaking capabilities of the Euclid telescope. Both the LF and SMF fit a Schechter plus Gaussian model. The LF features a dip at M(IE)=-19 and a faint-end slope of alpha_S = -1.2 to -1.3. The SMF displays a low-mass-end slope of alpha_S = -1.2 to -1.35. These observed slopes are flatter than those predicted for dark matter halos in cosmological simulations, offering significant insights for models of galaxy formation and evolution. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Submitted to A&A, 44 pages, 35 figures, 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: Early Release Observations -- Programme overview and pipeline for compact- and diffuse-emission photometry

Authors: J. -C. Cuillandre, E. Bertin, M. Bolzonella, H. Bouy, S. Gwyn, S. Isani, M. Kluge, O. Lai, A. Lançon, D. A. Lang, R. Laureijs, T. Saifollahi, M. Schirmer, C. Stone, Abdurro'uf, N. Aghanim, B. Altieri, F. Annibali, H. Atek, P. Awad, M. Baes, E. Bañados, D. Barrado, S. Belladitta, V. Belokurov , et al. (240 additional authors not shown)

Abstract: The Euclid ERO showcase Euclid's capabilities in advance of its main mission, targeting 17 astronomical objects, from galaxy clusters, nearby galaxies, globular clusters, to star-forming regions. A total of 24 hours observing time was allocated in the early months of operation, engaging the scientific community through an early public data release. We describe the development of the ERO pipeline t… ▽ More The Euclid ERO showcase Euclid's capabilities in advance of its main mission, targeting 17 astronomical objects, from galaxy clusters, nearby galaxies, globular clusters, to star-forming regions. A total of 24 hours observing time was allocated in the early months of operation, engaging the scientific community through an early public data release. We describe the development of the ERO pipeline to create visually compelling images while simultaneously meeting the scientific demands within months of launch, leveraging a pragmatic, data-driven development strategy. The pipeline's key requirements are to preserve the image quality and to provide flux calibration and photometry for compact and extended sources. The pipeline's five pillars are: removal of instrumental signatures; astrometric calibration; photometric calibration; image stacking; and the production of science-ready catalogues for both the VIS and NISP instruments. We report a PSF with a full width at half maximum of 0.16" in the optical and 0.49" in the three NIR bands. Our VIS mean absolute flux calibration is accurate to about 1%, and 10% for NISP due to a limited calibration set; both instruments have considerable colour terms. The median depth is 25.3 and 23.2 AB mag with a SNR of 10 for galaxies, and 27.1 and 24.5 AB mag at an SNR of 5 for point sources for VIS and NISP, respectively. Euclid's ability to observe diffuse emission is exceptional due to its extended PSF nearly matching a pure diffraction halo, the best ever achieved by a wide-field, high-resolution imaging telescope. Euclid offers unparalleled capabilities for exploring the LSB Universe across all scales, also opening a new observational window in the NIR. Median surface-brightness levels of 29.9 and 28.3 AB mag per square arcsec are achieved for VIS and NISP, respectively, for detecting a 10 arcsec x 10 arcsec extended feature at the 1 sigma level. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Submitted to A&A, 44 pages, 36 figures - 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

Weak lensing combined with the kinetic Sunyaev Zel'dovich effect: A study of baryonic feedback

Authors: L. Bigwood, A. Amon, A. Schneider, J. Salcido, I. G. McCarthy, C. Preston, D. Sanchez, D. Sijacki, E. Schaan, S. Ferraro, N. Battaglia, A. Chen, S. Dodelson, A. Roodman, A. Pieres, A. Ferte, A. Alarcon, A. Drlica-Wagner, A. Choi, A. Navarro-Alsina, A. Campos, A. J. Ross, A. Carnero Rosell, B. Yin, B. Yanny , et al. (100 additional authors not shown)

Abstract: Extracting precise cosmology from weak lensing surveys requires modelling the non-linear matter power spectrum, which is suppressed at small scales due to baryonic feedback processes. However, hydrodynamical galaxy formation simulations make widely varying predictions for the amplitude and extent of this effect. We use measurements of Dark Energy Survey Year 3 weak lensing (WL) and Atacama Cosmolo… ▽ More Extracting precise cosmology from weak lensing surveys requires modelling the non-linear matter power spectrum, which is suppressed at small scales due to baryonic feedback processes. However, hydrodynamical galaxy formation simulations make widely varying predictions for the amplitude and extent of this effect. We use measurements of Dark Energy Survey Year 3 weak lensing (WL) and Atacama Cosmology Telescope DR5 kinematic Sunyaev-Zel'dovich (kSZ) to jointly constrain cosmological and astrophysical baryonic feedback parameters using a flexible analytical model, `baryonification'. First, using WL only, we compare the $S_8$ constraints using baryonification to a simulation-calibrated halo model, a simulation-based emulator model and the approach of discarding WL measurements on small angular scales. We find that model flexibility can shift the value of $S_8$ and degrade the uncertainty. The kSZ provides additional constraints on the astrophysical parameters and shifts $S_8$ to $S_8=0.823^{+0.019}_{-0.020}$, a higher value than attained using the WL-only analysis. We measure the suppression of the non-linear matter power spectrum using WL + kSZ and constrain a mean feedback scenario that is more extreme than the predictions from most hydrodynamical simulations. We constrain the baryon fractions and the gas mass fractions and find them to be generally lower than inferred from X-ray observations and simulation predictions. We conclude that the WL + kSZ measurements provide a new and complementary benchmark for building a coherent picture of the impact of gas around galaxies across observations. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Assessing theoretical uncertainties for cosmological constraints from weak lensing surveys

Authors: Ting Tan, Dominik Zuercher, Janis Fluri, Alexandre Refregier, Federica Tarsitano, Tomasz Kacprzak

Abstract: $… ▽ More $ $Weak gravitational lensing is a powerful probe which is used to constrain the standard cosmological model and its extensions. With the enhanced statistical precision of current and upcoming surveys, high accuracy predictions for weak lensing statistics are needed to limit the impact of theoretical uncertainties on cosmological parameter constraints. For this purpose, we present a comparison of the theoretical predictions for the nonlinear matter and weak lensing power spectra, based on the widely used fitting functions ($\texttt{mead}$ and $\texttt{rev-halofit}$), emulators ($\texttt{EuclidEmulator}$, $\texttt{EuclidEmulator2}$, $\texttt{BaccoEmulator}$ and $\texttt{CosmicEmulator}$) and N-body simulations ($\texttt{Pkdgrav3}$). We consider the forecasted constraints on the $Λ\texttt{CDM}$ and $\texttt{wCDM}$ models from weak lensing for stage III and stage IV surveys. We study the relative bias on the constraints and their dependence on the assumed prescriptions. Assuming a $Λ\texttt{CDM}$ cosmology, we find that the relative agreement on the $S_8$ parameter is between $0.2-0.3σ$ for a stage III-like survey between the above predictors. For a stage IV-like survey the agreement becomes $1.4-3.0σ$. In the $\texttt{wCDM}$ scenario, we find broader $S_8$ constraints, and agreements of $0.18-0.26σ$ and $0.7-1.7σ$ for stage III and stage IV surveys, respectively. The accuracies of the above predictors therefore appear adequate for stage III surveys, while the fitting functions would need improvements for future stage IV weak lensing surveys. Furthermore, we find that, of the fitting functions, $\texttt{mead}$ provides the best agreement with the emulators. We discuss the implication of these findings for the preparation of the future weak lensing surveys. △ Less

Submitted 7 July, 2022; originally announced July 2022.

Comments: 18 pages, 10 figures

Symbolic Implementation of Extensions of the $\texttt{PyCosmo}$ Boltzmann Solver

Authors: Beatrice Moser, Christiane S. Lorenz, Uwe Schmitt, Alexandre Refregier, Janis Fluri, Raphael Sgier, Federica Tarsitano, Lavinia Heisenberg

Abstract: $\texttt{PyCosmo}$ is a Python-based framework for the fast computation of cosmological model predictions. One of its core features is the symbolic representation of the Einstein-Boltzmann system of equations. Efficient $\texttt{C/C++}$ code is generated from the $\texttt{SymPy}$ symbolic expressions making use of the $\texttt{sympy2c}… ▽ More $\texttt{PyCosmo}$ is a Python-based framework for the fast computation of cosmological model predictions. One of its core features is the symbolic representation of the Einstein-Boltzmann system of equations. Efficient $\texttt{C/C++}$ code is generated from the $\texttt{SymPy}$ symbolic expressions making use of the $\texttt{sympy2c}$ package. This enables easy extensions of the equation system for the implementation of new cosmological models. We illustrate this with three extensions of the $\texttt{PyCosmo}$ Boltzmann solver to include a dark energy component with a constant equation of state, massive neutrinos and a radiation streaming approximation. We describe the $\texttt{PyCosmo}$ framework, highlighting new features, and the symbolic implementation of the new models. We compare the $\texttt{PyCosmo}$ predictions for the $Λ$CDM model extensions with $\texttt{CLASS}$, both in terms of accuracy and computational speed. We find a good agreement, to better than 0.1% when using high-precision settings and a comparable computational speed. Links to the Python Package Index (PyPI) page of the code release and to the PyCosmo Hub, an online platform where the package is installed, are available at: https://cosmology.ethz.ch/research/software-lab/PyCosmo.html. △ Less

Submitted 17 June, 2022; v1 submitted 15 December, 2021; originally announced December 2021.

Comments: 35 pages including 5 figures and 3 tables. Link to $\texttt{PyCosmo}$ package: https://cosmology.ethz.ch/research/software-lab/PyCosmo.html

Journal ref: Astronomy and Computing 40 (2022) 100603

Combined $13\times2$-point analysis of the Cosmic Microwave Background and Large-Scale Structure: implications for the $S_8$-tension and neutrino mass constraints

Authors: Raphael Sgier, Christiane Lorenz, Alexandre Refregier, Janis Fluri, Dominik Zürcher, Federica Tarsitano

Abstract: We present cosmological constraints for the flat $Λ$CDM model, including the sum of neutrino masses, by performing a multi-probe analysis of a total of 13 tomographic auto- and cross-angular power spectra. This is achieved by combining, at map level, the latest primary CMB and CMB-lensing measurements from the Planck 2018 data release, as well as spectroscopic galaxy samples from BOSS DR12, and th… ▽ More We present cosmological constraints for the flat $Λ$CDM model, including the sum of neutrino masses, by performing a multi-probe analysis of a total of 13 tomographic auto- and cross-angular power spectra. This is achieved by combining, at map level, the latest primary CMB and CMB-lensing measurements from the Planck 2018 data release, as well as spectroscopic galaxy samples from BOSS DR12, and the latest Kilo-Degree Survey (KiDS-1000) tomographic weak lensing shear data release. Our analysis includes auto- and cross-correlations as well as calibration parameters for all cosmological probes, thus providing a self-calibration of the combined data sets. We find a good fit (reduced $χ^2$=1.7) for the combined probes with calibration parameters only moderately different from their nominal value, thus giving a possible interpretation of the tension between the early- and late-Universe probes. The resulting value for the structure growth parameter is $S_8 = 0.754 \pm 0.016$ (68\% CL). We also obtain a $\sim$2.3$σ$ constraint on the neutrino mass sum of $\sum m_ν= 0.51^{+0.21}_{-0.24}$ eV (68\% CL), which is compatible with current particle physics limits. We perform several tests by fixing the neutrino mass sum to a low value, considering narrower priors on the multiplicative bias parameters for cosmic shear, and by fixing all calibration parameters to their expected values. These tests result in worse fits compared to our fiducial run, especially for the case when all calibration parameters are fixed. This latter test also yields a lower upper limit of the neutrino mass sum. We discuss how the interplay between the cosmological and calibration parameters impact the $S_8$-tension and the constraints on the neutrino mass sum. [abridged] △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: 40 pages

Image feature extraction and galaxy classification: a novel and efficient approach with automated machine learning

Authors: F. Tarsitano, C. Bruderer, K. Schawinski, W. G. Hartley

Abstract: In this work we explore the possibility of applying machine learning methods designed for one-dimensional problems to the task of galaxy image classification. The algorithms used for image classification typically rely on multiple costly steps, such as the Point Spread Function (PSF) deconvolution and the training and application of complex Convolutional Neural Networks (CNN) of thousands or even… ▽ More In this work we explore the possibility of applying machine learning methods designed for one-dimensional problems to the task of galaxy image classification. The algorithms used for image classification typically rely on multiple costly steps, such as the Point Spread Function (PSF) deconvolution and the training and application of complex Convolutional Neural Networks (CNN) of thousands or even millions of parameters. In our approach, we extract features from the galaxy images by analysing the elliptical isophotes in their light distribution and collect the information in a sequence. The sequences obtained with this method present definite features allowing a direct distinction between galaxy types, as opposed to smooth Sérsic profiles. Then, we train and classify the sequences with machine learning algorithms, designed through the platform Modulos AutoML, and study how they optimize the classification task. As a demonstration of this method, we use the second public release of the Dark Energy Survey (DES DR2). We show that by applying it to this sample we are able to successfully distinguish between early-type and late-type galaxies, for images with signal-to-noise ratio greater then 300. This yields an accuracy of $86\%$ for the early-type galaxies and $93\%$ for the late-type galaxies, which is on par with most contemporary automated image classification approaches. Our novel method allows for galaxy images to be accurately classified and is faster than other approaches. Data dimensionality reduction also implies a significant lowering in computational cost. In the perspective of future data sets obtained with e.g. Euclid and the Vera Rubin Observatory (VRO), this work represents a path towards using a well-tested and widely used platform from industry in efficiently tackling galaxy classification problems at the peta-byte scale. △ Less

Submitted 3 May, 2021; originally announced May 2021.

Comments: 9 pages, 7 figures

Dark Energy Survey Year 3 Results: Deep Field Optical + Near-Infrared Images and Catalogue

Authors: W. G. Hartley, A. Choi, A. Amon, R. A. Gruendl, E. Sheldon, I. Harrison, G. M. Bernstein, I. Sevilla-Noarbe, B. Yanny, K. Eckert, H. T. Diehl, A. Alarcon, M. Banerji, K. Bechtol, R. Buchs, S. Cantu, C. Conselice, J. Cordero, C. Davis, T. M. Davis, S. Dodelson, A. Drlica-Wagner, S. Everett, A. Ferté, D. Gruen , et al. (93 additional authors not shown)

Abstract: We describe the Dark Energy Survey (DES) Deep Fields, a set of images and associated multi-wavelength catalogue ($ugrizJHKs$) built from Dark Energy Camera (DECam) and Visible and Infrared Survey Telescope for Astronomy (VISTA) data. The DES Deep Fields comprise 11 fields (10 DES supernova fields plus COSMOS), with a total area of $\sim30~$ square degrees in $ugriz$ bands and reaching a maximum… ▽ More We describe the Dark Energy Survey (DES) Deep Fields, a set of images and associated multi-wavelength catalogue ($ugrizJHKs$) built from Dark Energy Camera (DECam) and Visible and Infrared Survey Telescope for Astronomy (VISTA) data. The DES Deep Fields comprise 11 fields (10 DES supernova fields plus COSMOS), with a total area of $\sim30~$ square degrees in $ugriz$ bands and reaching a maximum $i$-band depth of 26.75 (AB, $10σ$, 2 arcsec). We present a catalogue for the DES 3-year cosmology analysis of those four fields with full 8-band coverage, totalling $5.88~$ sq. deg. after masking. Numbering $2.8~$million objects ($1.6~$million post masking), our catalogue is drawn from images coadded to consistent depths of $r=25.7, i=25, z=24.3$ mag. We use a new model-fitting code, built upon established methods, to deblend sources and ensure consistent colours across the $u$-band to $Ks$-band wavelength range. We further detail the tight control we maintain over the point-spread function modelling required for the model fitting, astrometry and consistency of photometry between the four fields. The catalogue allows us to perform a careful star-galaxy separation and produces excellent photometric redshift performance (${\rm NMAD} = 0.023$ at $i<23$). The Deep-Fields catalogue will be made available as part of the cosmology data products release, following the completion of the DES 3-year weak lensing and galaxy clustering cosmology work. △ Less

Submitted 16 February, 2022; v1 submitted 23 December, 2020; originally announced December 2020.

Comments: 32 pages, 23 figures, version accepted by MNRAS. See https://www.darkenergysurvey.org/des-year-3-cosmology-results-papers/ for the full DES Y3 cosmology release

Report number: FERMILAB-PUB-20-670-AE

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 509, Issue 3, pp.3547-3579, 2022

DES Y3 results: Blending shear and redshift biases in image simulations

Authors: N. MacCrann, M. R. Becker, J. McCullough, A. Amon, D. Gruen, M. Jarvis, A. Choi, M. A. Troxel, E. Sheldon, B. Yanny, K. Herner, S. Dodelson, J. Zuntz, K. Eckert, R. P. Rollins, T. N. Varga, G. M. Bernstein, R. A. Gruendl, I. Harrison, W. G. Hartley, I. Sevilla-Noarbe, A. Pieres, S. L. Bridle, J. Myles, A. Alarcon , et al. (71 additional authors not shown)

Abstract: As the statistical power of galaxy weak lensing reaches percent level precision, large, realistic and robust simulations are required to calibrate observational systematics, especially given the increased importance of object blending as survey depths increase. To capture the coupled effects of blending in both shear and photometric redshift calibration, we define the effective redshift distributi… ▽ More As the statistical power of galaxy weak lensing reaches percent level precision, large, realistic and robust simulations are required to calibrate observational systematics, especially given the increased importance of object blending as survey depths increase. To capture the coupled effects of blending in both shear and photometric redshift calibration, we define the effective redshift distribution for lensing, $n_γ(z)$, and describe how to estimate it using image simulations. We use an extensive suite of tailored image simulations to characterize the performance of the shear estimation pipeline applied to the Dark Energy Survey (DES) Year 3 dataset. We describe the multi-band, multi-epoch simulations, and demonstrate their high level of realism through comparisons to the real DES data. We isolate the effects that generate shear calibration biases by running variations on our fiducial simulation, and find that blending-related effects are the dominant contribution to the mean multiplicative bias of approximately $-2\%$. By generating simulations with input shear signals that vary with redshift, we calibrate biases in our estimation of the effective redshfit distribution, and demonstrate the importance of this approach when blending is present. We provide corrected effective redshift distributions that incorporate statistical and systematic uncertainties, ready for use in DES Year 3 weak lensing analyses. △ Less

Submitted 26 January, 2022; v1 submitted 15 December, 2020; originally announced December 2020.

Comments: Submitted to MNRAS. See https://www.darkenergysurvey.org/des-year-3-cosmology-results-papers/ for the full DES Y3 cosmology release. Version accepted by mnras

Pushing automated morphological classifications to their limits with the Dark Energy Survey

Authors: J. Vega-Ferrero, H. Domínguez Sánchez, M. Bernardi, M. Huertas-Company, R. Morgan, B. Margalef, M. Aguena, S. Allam, J. Annis, S. Avila, D. Bacon, E. Bertin, D. Brooks, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, A. Choi, C. Conselice, M. Costanzi, L. N. da Costa, M. E. S. Pereira, J. De Vicente, S. Desai, I. Ferrero, P. Fosalba , et al. (38 additional authors not shown)

Abstract: We present morphological classifications of $\sim$27 million galaxies from the Dark Energy Survey (DES) Data Release 1 (DR1) using a supervised deep learning algorithm. The classification scheme separates: (a) early-type galaxies (ETGs) from late-types (LTGs); and (b) face-on galaxies from edge-on. Our Convolutional Neural Networks (CNNs) are trained on a small subset of DES objects with previousl… ▽ More We present morphological classifications of $\sim$27 million galaxies from the Dark Energy Survey (DES) Data Release 1 (DR1) using a supervised deep learning algorithm. The classification scheme separates: (a) early-type galaxies (ETGs) from late-types (LTGs); and (b) face-on galaxies from edge-on. Our Convolutional Neural Networks (CNNs) are trained on a small subset of DES objects with previously known classifications. These typically have $\mathrm{m}_r \lesssim 17.7~\mathrm{mag}$; we model fainter objects to $\mathrm{m}_r < 21.5$ mag by simulating what the brighter objects with well determined classifications would look like if they were at higher redshifts. The CNNs reach 97\% accuracy to $\mathrm{m}_r<21.5$ on their training sets, suggesting that they are able to recover features more accurately than the human eye. We then used the trained CNNs to classify the vast majority of the other DES images. The final catalog comprises five independent CNN predictions for each classification scheme, helping to determine if the CNN predictions are robust or not. We obtain secure classifications for $\sim$ 87\% and 73\% of the catalog for the ETG vs. LTG and edge-on vs. face-on models, respectively. Combining the two classifications (a) and (b) helps to increase the purity of the ETG sample and to identify edge-on lenticular galaxies (as ETGs with high ellipticity). Where a comparison is possible, our classifications correlate very well with Sérsic index (\textit{n}), ellipticity ($ε$) and spectral type, even for the fainter galaxies. This is the largest multi-band catalog of automated galaxy morphologies to date. △ Less

Submitted 9 March, 2021; v1 submitted 14 December, 2020; originally announced December 2020.

Comments: Accepted for publication in MNRAS (2021 February 22); 17 pages, 16 figures

Predicting Cosmological Observables with PyCosmo

Authors: F. Tarsitano, U. Schmitt, A. Refregier, J. Fluri, R. Sgier, A. Nicola, J. Herbel, A. Amara, T. Kacprzak, L. Heisenberg

Abstract: Current and upcoming cosmological experiments open a new era of precision cosmology, thus demanding accurate theoretical predictions for cosmological observables. Because of the complexity of the codes delivering such predictions, reaching a high level of numerical accuracy is challenging. Among the codes already fulfilling this task, $\textsf{PyCosmo}$ is a Python based framework providing soluti… ▽ More Current and upcoming cosmological experiments open a new era of precision cosmology, thus demanding accurate theoretical predictions for cosmological observables. Because of the complexity of the codes delivering such predictions, reaching a high level of numerical accuracy is challenging. Among the codes already fulfilling this task, $\textsf{PyCosmo}$ is a Python based framework providing solutions to the Einstein-Boltzmann equations and accurate predictions for cosmological observables. In this work, we first describe how the observables are implemented. Then, we check the accuracy of the theoretical predictions for background quantities, power spectra and Limber and beyond-Limber angular power spectra by comparison with other codes: the Core Cosmology Library ($\texttt{CCL}$), $\texttt{CLASS}$, $\texttt{HMCode}$ and $\texttt{iCosmo}$. In our analysis we quantify the agreement of $\textsf{PyCosmo}$ with the other codes, for a range of cosmological models, monitored through a series of $\textit{unit tests}$. $\textsf{PyCosmo}$, conceived as a multi purpose cosmology calculation tool in $\texttt{Python}$, is designed to be interactive and user friendly. A current version of the code (without the Boltzmann Solver) is publicly available and can be used interactively on the platform $\textsf{PyCosmo Hub}$, all accessible from this link: https://cosmology.ethz.ch/research/software-lab/PyCosmo.html . On the hub the users can perform their own computations using $\texttt{Jupyter Notebooks}$ without the need of installing any software, access to the results presented in this work and benefit from tutorial notebooks illustrating the usage of the code. The link above also redirects to the code release and documentation. △ Less

Submitted 1 May, 2020; originally announced May 2020.

Comments: PyCosmo is available online on the PyCosmo Hub: https://pycosmohub.com/hub/login . More information here: https://cosmology.ethz.ch/research/software-lab/PyCosmo.html . 37 pages, 20 figures, 1 table. Submitted to Astronomy and Computing

Monte Carlo Control Loops for cosmic shear cosmology with DES Year 1

Authors: T. Kacprzak, J. Herbel, A. Nicola, R. Sgier, F. Tarsitano, C. Bruderer, A. Amara, A. Refregier, S. L. Bridle, A. Drlica-Wagner, D. Gruen, W. G. Hartley, B. Hoyle, L. F. Secco, J. Zuntz, J. Annis, S. Avila, E. Bertin, D. Brooks, E. Buckley-Geer, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, L. N. da Costa, J. De Vicente , et al. (30 additional authors not shown)

Abstract: Weak lensing by large-scale structure is a powerful probe of cosmology and of the dark universe. This cosmic shear technique relies on the accurate measurement of the shapes and redshifts of background galaxies and requires precise control of systematic errors. The Monte Carlo Control Loops (MCCL) is a forward modelling method designed to tackle this problem. It relies on the Ultra Fast Image Gene… ▽ More Weak lensing by large-scale structure is a powerful probe of cosmology and of the dark universe. This cosmic shear technique relies on the accurate measurement of the shapes and redshifts of background galaxies and requires precise control of systematic errors. The Monte Carlo Control Loops (MCCL) is a forward modelling method designed to tackle this problem. It relies on the Ultra Fast Image Generator (UFig) to produce simulated images tuned to match the target data statistically, followed by calibrations and tolerance loops. We present the first end-to-end application of this method, on the Dark Energy Survey (DES) Year 1 wide field imaging data. We simultaneously measure the shear power spectrum $C_{\ell}$ and the redshift distribution $n(z)$ of the background galaxy sample. The method includes maps of the systematic sources, Point Spread Function (PSF), an Approximate Bayesian Computation (ABC) inference of the simulation model parameters, a shear calibration scheme, and the fast estimation of the covariance matrix. We find a close statistical agreement between the simulations and the DES Y1 data using an array of diagnostics. In a non-tomographic setting, we derive a set of $C_\ell$ and $n(z)$ curves that encode the cosmic shear measurement, as well as the systematic uncertainty. Following a blinding scheme, we measure the combination of $Ω_m$, $σ_8$, and intrinsic alignment amplitude $A_{\rm{IA}}$, defined as $S_8D_{\rm{IA}} = σ_8(Ω_m/0.3)^{0.5}D_{\rm{IA}}$, where $D_{\rm{IA}}=1-0.11(A_{\rm{IA}}-1)$. We find $S_8D_{\rm{IA}}=0.895^{+0.054}_{-0.039}$, where systematics are at the level of roughly 60\% of the statistical errors. We discuss these results in the context of earlier cosmic shear analyses of the DES Y1 data. Our findings indicate that this method and its fast runtime offer good prospects for cosmic shear measurements with future wide-field surveys. △ Less

Submitted 3 June, 2019; originally announced June 2019.

Comments: 31 pages, 15 figures, 2 tables

Report number: DES-2018-0362, FERMILAB-PUB-19-247-AE

Journal ref: Phys. Rev. D 101, 082003 (2020)

A catalogue of structural and morphological measurements for DES Y1

Authors: F. Tarsitano, W. G. Hartley, A. Amara, A. Bluck, C. Bruderer, M. Carollo, C. Conselice, P. Melchior, B. Moraes, A. Refregier, I. Sevilla-Noarbe, J. Woo, T. M. C. Abbott, S. Allam, J. Annis, S. Avila, M. Banerji, E. Bertin, D. Brooks, D. L. Burke, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, C. E. Cunha, C. B. D'Andrea , et al. (38 additional authors not shown)

Abstract: We present a structural and morphological catalogue for 45 million objects selected from the first year of data from the Dark Energy Survey (DES). Single Sersic fits and non-parametric measurements are produced for g, r and i filters. The parameters from the best-fitting Sersic model (total magnitude, half-light radius, Sersic index, axis ratio and position angle) are measured with Galfit; the non… ▽ More We present a structural and morphological catalogue for 45 million objects selected from the first year of data from the Dark Energy Survey (DES). Single Sersic fits and non-parametric measurements are produced for g, r and i filters. The parameters from the best-fitting Sersic model (total magnitude, half-light radius, Sersic index, axis ratio and position angle) are measured with Galfit; the non-parametric coefficients (concentration, asymmetry, clumpiness, Gini, M20) are provided using the Zurich Estimator of Structural Types (ZEST+). To study the statistical uncertainties, we consider a sample of state-of-the-art image simulations with a realistic distribution in the input parameter space and then process and analyse them as we do with real data: this enables us to quantify the observational biases due to PSF blurring and magnitude effects and correct the measurements as a function of magnitude, galaxy size, Sersic index (concentration for the analysis of the non-parametric measurements) and ellipticity. We present the largest structural catalogue to date: we find that accurate and complete measurements for all the structural parameters are typically obtained for galaxies with SExtractor MAG AUTO I < 21. Indeed, the parameters in the filters i and r can be overall well recovered up to MAG AUTO < 21.5, corresponding to a fitting completeness of ~90% below this threshold, for a total of 25 million galaxies. The combination of parametric and non-parametric structural measurements makes this catalogue an important instrument to explore and understand how galaxies form and evolve. The catalogue described in this paper will be publicly released alongside the Dark Energy Survey collaboration Y1 cosmology data products at the following URL: https://des.ncsa.illinois.edu/releases/y1a1/gold/morphology. △ Less

Submitted 2 October, 2018; v1 submitted 27 July, 2018; originally announced July 2018.

Comments: Accepted for publication in MNRAS. 26 pages, 16 figures. Catalogue data are available at https://des.ncsa.illinois.edu/releases/y1a1/gold/morphology

Evidence for Dynamically Driven Formation of the GW170817 Neutron Star Binary in NGC 4993

Authors: A. Palmese, W. Hartley, F. Tarsitano, C. Conselice, O. Lahav, S. Allam, J. Annis, H. Lin, M. Soares-Santos, D. Tucker, D. Brout, M. Banerji, K. Bechtol, H. T. Diehl, A. Fruchter, J. Garcia-Bellido, K. Herner, A. J. Levan, T. S. Li, C. Lidman, K. Misra, M. Sako, D. Scolnic, M. Smith, T. M. C. Abbott , et al. (67 additional authors not shown)

Abstract: We present a study of NGC 4993, the host galaxy of the GW170817 gravitational wave event, the GRB170817A short gamma-ray burst (sGRB) and the AT2017gfo kilonova. We use Dark Energy Camera imaging, AAT spectra and publicly available data, relating our findings to binary neutron star (BNS) formation scenarios and merger delay timescales. NGC4993 is a nearby (40 Mpc) early-type galaxy, with $i$-band… ▽ More We present a study of NGC 4993, the host galaxy of the GW170817 gravitational wave event, the GRB170817A short gamma-ray burst (sGRB) and the AT2017gfo kilonova. We use Dark Energy Camera imaging, AAT spectra and publicly available data, relating our findings to binary neutron star (BNS) formation scenarios and merger delay timescales. NGC4993 is a nearby (40 Mpc) early-type galaxy, with $i$-band Sérsic index $n=4.0$ and low asymmetry ($A=0.04\pm 0.01$). These properties are unusual for sGRB hosts. However, NGC4993 presents shell-like structures and dust lanes indicative of a recent galaxy merger, with the optical transient located close to a shell. We constrain the star formation history (SFH) of the galaxy assuming that the galaxy merger produced a star formation burst, but find little to no on-going star formation in either spatially-resolved broadband SED or spectral fitting. We use the best-fit SFH to estimate the BNS merger rate in this type of galaxy, as $R_{NSM}^{gal}= 5.7^{+0.57}_{-3.3} \times 10^{-6} {\rm yr}^{-1}$. If star formation is the only considered BNS formation scenario, the expected number of BNS mergers from early-type galaxies detectable with LIGO during its first two observing seasons is $0.038^{+0.004}_{-0.022}$, as opposed to $\sim 0.5$ from all galaxy types. Hypothesizing that the binary system formed due to dynamical interactions during the galaxy merger, the subsequent time elapsed can constrain the delay time of the BNS coalescence. By using velocity dispersion estimates and the position of the shells, we find that the galaxy merger occurred $t_{\rm mer}\lesssim 200~{\rm Myr}$ prior to the BNS coalescence. △ Less

Submitted 13 November, 2017; v1 submitted 18 October, 2017; originally announced October 2017.

Comments: 9 pages, 4 figures, 1 table. Edited to match the ApJL version

Report number: FERMILAB-PUB-17-466-AE-CD-PPD

Journal ref: 2017 The Astrophysical Journal Letters, Volume 849, Number 2, L34