Showing 1–48 of 48 results for author: Peel, A

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'

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

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

Microlensing of strongly lensed quasars

Authors: G. Vernardos, D. Sluse, D. Pooley, R. W. Schmidt, M. Millon, L. Weisenbach, V. Motta, T. Anguita, P. Saha, M. O'Dowd, A. Peel, P. L. Schechter

Abstract: Strong gravitational lensing of quasars has the potential to unlock the poorly understood physics of these fascinating objects, as well as serve as a probe of the lensing mass distribution and of cosmological parameters. In particular, gravitational microlensing by compact bodies in the lensing galaxy can enable mapping of quasar structure to $\lt 10^{-6}$ arcsec scales. Some of this potential has… ▽ More Strong gravitational lensing of quasars has the potential to unlock the poorly understood physics of these fascinating objects, as well as serve as a probe of the lensing mass distribution and of cosmological parameters. In particular, gravitational microlensing by compact bodies in the lensing galaxy can enable mapping of quasar structure to $\lt 10^{-6}$ arcsec scales. Some of this potential has been realized over the past few decades, however the upcoming era of large sky surveys promises to bring this to full fruition. Here we review the theoretical framework of this field, describe the prominent current methods for parameter inference from quasar microlensing data across different observing modalities, and discuss the constraints so far derived on the geometry and physics of quasar inner structure. We also review the application of strong lensing and microlensing to constraining the granularity of the lens potential, i.e. the contribution of the baryonic and dark matter components, and the local mass distribution in the lens, i.e. the stellar mass function. Finally, we discuss the future of the field, including the new possibilities that will be opened by the next generation of large surveys and by new analysis methods now being developed. △ Less

Submitted 1 December, 2023; originally announced December 2023.

Comments: To be submitted to Space Science Reviews, Topical Collection "Strong Gravitational Lensing", eds. J. Wambsganss et al

SKA Science Data Challenge 2: analysis and results

Authors: P. Hartley, A. Bonaldi, R. Braun, J. N. H. S. Aditya, S. Aicardi, L. Alegre, A. Chakraborty, X. Chen, S. Choudhuri, A. O. Clarke, J. Coles, J. S. Collinson, D. Cornu, L. Darriba, M. Delli Veneri, J. Forbrich, B. Fraga, A. Galan, J. Garrido, F. Gubanov, H. Håkansson, M. J. Hardcastle, C. Heneka, D. Herranz, K. M. Hess , et al. (83 additional authors not shown)

Abstract: The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed t… ▽ More The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed to familiarise the scientific community with SKAO data and to drive the development of new analysis techniques. We present the results from Science Data Challenge 2 (SDC2), which invited participants to find and characterise 233245 neutral hydrogen (Hi) sources in a simulated data product representing a 2000~h SKA MID spectral line observation from redshifts 0.25 to 0.5. Through the generous support of eight international supercomputing facilities, participants were able to undertake the Challenge using dedicated computational resources. Alongside the main challenge, `reproducibility awards' were made in recognition of those pipelines which demonstrated Open Science best practice. The Challenge saw over 100 participants develop a range of new and existing techniques, with results that highlight the strengths of multidisciplinary and collaborative effort. The winning strategy -- which combined predictions from two independent machine learning techniques to yield a 20 percent improvement in overall performance -- underscores one of the main Challenge outcomes: that of method complementarity. It is likely that the combination of methods in a so-called ensemble approach will be key to exploiting very large astronomical datasets. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: Under review by MNRAS; 28 pages, 16 figures

JAX-COSMO: An End-to-End Differentiable and GPU Accelerated Cosmology Library

Authors: Jean-Eric Campagne, François Lanusse, Joe Zuntz, Alexandre Boucaud, Santiago Casas, Minas Karamanis, David Kirkby, Denise Lanzieri, Yin Li, Austin Peel

Abstract: We present jax-cosmo, a library for automatically differentiable cosmological theory calculations. It uses the JAX library, which has created a new coding ecosystem, especially in probabilistic programming. As well as batch acceleration, just-in-time compilation, and automatic optimization of code for different hardware modalities (CPU, GPU, TPU), JAX exposes an automatic differentiation (autodiff… ▽ More We present jax-cosmo, a library for automatically differentiable cosmological theory calculations. It uses the JAX library, which has created a new coding ecosystem, especially in probabilistic programming. As well as batch acceleration, just-in-time compilation, and automatic optimization of code for different hardware modalities (CPU, GPU, TPU), JAX exposes an automatic differentiation (autodiff) mechanism. Thanks to autodiff, jax-cosmo gives access to the derivatives of cosmological likelihoods with respect to any of their parameters, and thus enables a range of powerful Bayesian inference algorithms, otherwise impractical in cosmology, such as Hamiltonian Monte Carlo and Variational Inference. In its initial release, jax-cosmo implements background evolution, linear and non-linear power spectra (using halofit or the Eisenstein and Hu transfer function), as well as angular power spectra with the Limber approximation for galaxy and weak lensing probes, all differentiable with respect to the cosmological parameters and their other inputs. We illustrate how autodiff can be a game-changer for common tasks involving Fisher matrix computations, or full posterior inference with gradient-based techniques. In particular, we show how Fisher matrices are now fast, exact, no longer require any fine tuning, and are themselves differentiable. Finally, using a Dark Energy Survey Year 1 3x2pt analysis as a benchmark, we demonstrate how jax-cosmo can be combined with Probabilistic Programming Languages to perform posterior inference with state-of-the-art algorithms including a No U-Turn Sampler, Automatic Differentiation Variational Inference,and Neural Transport HMC. We further demonstrate that Normalizing Flows using Neural Transport are a promising methodology for model validation in the early stages of analysis. △ Less

Submitted 27 April, 2023; v1 submitted 10 February, 2023; originally announced February 2023.

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. XXXII. Evaluating the weak lensing cluster mass biases using the Three Hundred Project hydrodynamical simulations

Authors: Euclid Collaboration, C. Giocoli, M. Meneghetti, E. Rasia, S. Borgani, G. Despali, G. F. Lesci, F. Marulli, L. Moscardini, M. Sereno, W. Cui, A. Knebe, G. Yepes, T. Castro, P. -S. Corasaniti, S. Pires, G. Castignani, L. Ingoglia, T. Schrabback, G. W. Pratt, A. M. C. Le Brun, N. Aghanim, L. Amendola, N. Auricchio, M. Baldi , et al. (191 additional authors not shown)

Abstract: The photometric catalogue of galaxy clusters extracted from ESA Euclid data is expected to be very competitive for cosmological studies. Using state-of-the-art hydrodynamical simulations, we present systematic analyses simulating the expected weak lensing profiles from clusters in a variety of dynamic states and at wide range of redshifts. In order to derive cluster masses, we use a model consiste… ▽ More The photometric catalogue of galaxy clusters extracted from ESA Euclid data is expected to be very competitive for cosmological studies. Using state-of-the-art hydrodynamical simulations, we present systematic analyses simulating the expected weak lensing profiles from clusters in a variety of dynamic states and at wide range of redshifts. In order to derive cluster masses, we use a model consistent with the implementation within the Euclid Consortium of the dedicated processing function and find that, when jointly modelling mass and the concentration parameter of the Navarro-Frenk-White halo profile, the weak lensing masses tend to be, on average, biased low by 5-10% with respect to the true mass, up to z=0.5. Using a fixed value for the concentration $c_{200} = 3$, the mass bias is diminished below 5%, up to z=0.7, along with its relative uncertainty. Simulating the weak lensing signal by projecting along the directions of the axes of the moment of inertia tensor ellipsoid, we find that orientation matters: when clusters are oriented along the major axis, the lensing signal is boosted, and the recovered weak lensing mass is correspondingly overestimated. Typically, the weak lensing mass bias of individual clusters is modulated by the weak lensing signal-to-noise ratio, related to the redshift evolution of the number of galaxies used for weak lensing measurements: the negative mass bias tends to be larger toward higher redshifts. However, when we use a fixed value of the concentration parameter, the redshift evolution trend is reduced. These results provide a solid basis for the weak-lensing mass calibration required by the cosmological application of future cluster surveys from Euclid and Rubin. △ Less

Submitted 18 October, 2023; v1 submitted 1 February, 2023; originally announced February 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

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

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)

Modeling lens potentials with continuous neural fields in galaxy-scale strong lenses

Authors: Luca Biggio, Georgios Vernardos, Aymeric Galan, Austin Peel

Abstract: Strong gravitational lensing is a unique observational tool for studying the dark and luminous mass distribution both within and between galaxies. Given the presence of substructures, current strong lensing observations demand more complex mass models than smooth analytical profiles, such as power-law ellipsoids. In this work, we introduce a continuous neural field to predict the lensing potential… ▽ More Strong gravitational lensing is a unique observational tool for studying the dark and luminous mass distribution both within and between galaxies. Given the presence of substructures, current strong lensing observations demand more complex mass models than smooth analytical profiles, such as power-law ellipsoids. In this work, we introduce a continuous neural field to predict the lensing potential at any position throughout the image plane, allowing for a nearly model-independent description of the lensing mass. We apply our method on simulated Hubble Space Telescope imaging data containing different types of perturbations to a smooth mass distribution: a localized dark subhalo, a population of subhalos, and an external shear perturbation. Assuming knowledge of the source surface brightness, we use the continuous neural field to model either the perturbations alone or the full lensing potential. In both cases, the resulting model is able to fit the imaging data, and we are able to accurately recover the properties of both the smooth potential and of the perturbations. Unlike many other deep learning methods, ours explicitly retains lensing physics (i.e., the lens equation) and introduces high flexibility in the model only where required, namely, in the lens potential. Moreover, the neural network does not require pre-training on large sets of labelled data and predicts the potential from the single observed lensing image. Our model is implemented in the fully differentiable lens modeling code Herculens. △ Less

Submitted 17 October, 2022; originally announced October 2022.

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. XXIV. Calibration of the halo mass function in $Λらむだ(νにゅー)$CDM cosmologies

Authors: Euclid Collaboration, T. Castro, A. Fumagalli, R. E. Angulo, S. Bocquet, S. Borgani, C. Carbone, J. Dakin, K. Dolag, C. Giocoli, P. Monaco, A. Ragagnin, A. Saro, E. Sefusatti, M. Costanzi, A. M. C. Le Brun, P. -S. Corasaniti, A. Amara, L. Amendola, M. Baldi, R. Bender, C. Bodendorf, E. Branchini, M. Brescia, S. Camera , et al. (157 additional authors not shown)

Abstract: Euclid's photometric galaxy cluster survey has the potential to be a very competitive cosmological probe. The main cosmological probe with observations of clusters is their number count, within which the halo mass function (HMF) is a key theoretical quantity. We present a new calibration of the analytic HMF, at the level of accuracy and precision required for the uncertainty in this quantity to be… ▽ More Euclid's photometric galaxy cluster survey has the potential to be a very competitive cosmological probe. The main cosmological probe with observations of clusters is their number count, within which the halo mass function (HMF) is a key theoretical quantity. We present a new calibration of the analytic HMF, at the level of accuracy and precision required for the uncertainty in this quantity to be subdominant with respect to other sources of uncertainty in recovering cosmological parameters from Euclid cluster counts. Our model is calibrated against a suite of N-body simulations using a Bayesian approach taking into account systematic errors arising from numerical effects in the simulation. First, we test the convergence of HMF predictions from different N-body codes, by using initial conditions generated with different orders of Lagrangian Perturbation theory, and adopting different simulation box sizes and mass resolution. Then, we quantify the effect of using different halo-finder algorithms, and how the resulting differences propagate to the cosmological constraints. In order to trace the violation of universality in the HMF, we also analyse simulations based on initial conditions characterised by scale-free power spectra with different spectral indexes, assuming both Einstein--de Sitter and standard $Λらむだ$CDM expansion histories. Based on these results, we construct a fitting function for the HMF that we demonstrate to be sub-percent accurate in reproducing results from 9 different variants of the $Λらむだ$CDM model including massive neutrinos cosmologies. The calibration systematic uncertainty is largely sub-dominant with respect to the expected precision of future mass-observation relations; with the only notable exception of the effect due to the halo finder, that could lead to biased cosmological inference. △ Less

Submitted 16 March, 2023; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: 25 pages, 21 figures, 5 tables, 3 appendixes; v2 matches published version

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

Using wavelets to capture deviations from smoothness in galaxy-scale strong lenses

Authors: Aymeric Galan, Georgios Vernardos, Austin Peel, Frédéric Courbin, Jean-Luc Starck

Abstract: Modeling the mass distribution of galaxy-scale strong gravitational lenses is a task of increasing difficulty. The high-resolution and depth of imaging data now available render simple analytical forms ineffective at capturing lens structures spanning a large range in spatial scale, mass scale, and morphology. In this work, we address the problem with a novel multiscale method based on wavelets. W… ▽ More Modeling the mass distribution of galaxy-scale strong gravitational lenses is a task of increasing difficulty. The high-resolution and depth of imaging data now available render simple analytical forms ineffective at capturing lens structures spanning a large range in spatial scale, mass scale, and morphology. In this work, we address the problem with a novel multiscale method based on wavelets. We tested our method on simulated Hubble Space Telescope (HST) imaging data of strong lenses containing the following different types of mass substructures making them deviate from smooth models: (1) a localized small dark matter subhalo, (2) a Gaussian random field (GRF) that mimics a nonlocalized population of subhalos along the line of sight, and (3) galaxy-scale multipoles that break elliptical symmetry. We show that wavelets are able to recover all of these structures accurately. This is made technically possible by using gradient-informed optimization based on automatic differentiation over thousands of parameters, which also allow us to sample the posterior distributions of all model parameters simultaneously. By construction, our method merges the two main modeling paradigms - analytical and pixelated - with machine-learning optimization techniques into a single modular framework. It is also well-suited for the fast modeling of large samples of lenses. All methods presented here are publicly available in our new Herculens package. △ Less

Submitted 21 November, 2022; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: 24 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 668, A155 (2022)

Euclid preparation: XXIII. Derivation of galaxy physical properties with deep machine learning using mock fluxes and H-band images

Authors: Euclid Collaboration, L. Bisigello, C. J. Conselice, M. Baes, M. Bolzonella, M. Brescia, S. Cavuoti, O. Cucciati, A. Humphrey, L. K. Hunt, C. Maraston, L. Pozzetti, C. Tortora, S. E. van Mierlo, N. Aghanim, N. Auricchio, M. Baldi, R. Bender, C. Bodendorf, D. Bonino, E. Branchini, J. Brinchmann, S. Camera, V. Capobianco, C. Carbone , et al. (174 additional authors not shown)

Abstract: Next generation telescopes, like Euclid, Rubin/LSST, and Roman, will open new windows on the Universe, allowing us to infer physical properties for tens of millions of galaxies. Machine learning methods are increasingly becoming the most efficient tools to handle this enormous amount of data, because they are often faster and more accurate than traditional methods. We investigate how well redshift… ▽ More Next generation telescopes, like Euclid, Rubin/LSST, and Roman, will open new windows on the Universe, allowing us to infer physical properties for tens of millions of galaxies. Machine learning methods are increasingly becoming the most efficient tools to handle this enormous amount of data, because they are often faster and more accurate than traditional methods. We investigate how well redshifts, stellar masses, and star-formation rates (SFR) can be measured with deep learning algorithms for observed galaxies within data mimicking the Euclid and Rubin/LSST surveys. We find that Deep Learning Neural Networks and Convolutional Neutral Networks (CNN), which are dependent on the parameter space of the training sample, perform well in measuring the properties of these galaxies and have a better accuracy than methods based on spectral energy distribution fitting. CNNs allow the processing of multi-band magnitudes together with $H_{\scriptscriptstyle\rm E}$-band images. We find that the estimates of stellar masses improve with the use of an image, but those of redshift and SFR do not. Our best results are deriving i) the redshift within a normalised error of less than 0.15 for 99.9$\%$ of the galaxies with S/N>3 in the $H_{\scriptscriptstyle\rm E}$-band; ii) the stellar mass within a factor of two ($\sim0.3 \rm dex$) for 99.5$\%$ of the considered galaxies; iii) the SFR within a factor of two ($\sim0.3 \rm dex$) for $\sim$70$\%$ of the sample. We discuss the implications of our work for application to surveys as well as how measurements of these galaxy parameters can be improved with deep learning. △ Less

Submitted 4 January, 2023; v1 submitted 29 June, 2022; originally announced June 2022.

Comments: accepted for publication in MNRAS, 21 pages, 22 figures, 6 tables

ShapePipe: A modular weak-lensing processing and analysis pipeline

Authors: S. Farrens, A. Guinot, M. Kilbinger, T. Liaudat, L. Baumont, X. Jimenez, A. Peel, A. Pujol, M. Schmitz, J. -L. Starck, A. Z. Vitorelli

Abstract: We present the first public release of ShapePipe, an open-source and modular weak-lensing measurement, analysis, and validation pipeline written in Python. We describe the design of the software and justify the choices made. We provide a brief description of all the modules currently available and summarise how the pipeline has been applied to real Ultraviolet Near-Infrared Optical Northern Survey… ▽ More We present the first public release of ShapePipe, an open-source and modular weak-lensing measurement, analysis, and validation pipeline written in Python. We describe the design of the software and justify the choices made. We provide a brief description of all the modules currently available and summarise how the pipeline has been applied to real Ultraviolet Near-Infrared Optical Northern Survey data. Finally, we mention plans for future applications and development. The code and accompanying documentation are publicly available on GitHub. △ Less

Submitted 29 June, 2022; originally announced June 2022.

Comments: 7 pages, 2 figures, accepted for publication in A&A

Journal ref: A&A 664, A141 (2022)

Euclid preparation. XXI. Intermediate-redshift contaminants in the search for $z>6$ galaxies within the Euclid Deep Survey

Authors: Euclid Collaboration, S. E. van Mierlo, K. I. Caputi, M. Ashby, H. Atek, M. Bolzonella, R. A. A. Bowler, G. Brammer, C. J. Conselice, J. Cuby, P. Dayal, A. Díaz-Sánchez, S. L. Finkelstein, H. Hoekstra, A. Humphrey, O. Ilbert, H. J. McCracken, B. Milvang-Jensen, P. A. Oesch, R. Pello, G. Rodighiero, M. Schirmer, S. Toft, J. R. Weaver, S. M. Wilkins , et al. (181 additional authors not shown)

Abstract: (Abridged) The Euclid mission is expected to discover thousands of z>6 galaxies in three Deep Fields, which together will cover a ~40 deg2 area. However, the limited number of Euclid bands and availability of ancillary data could make the identification of z>6 galaxies challenging. In this work, we assess the degree of contamination by intermediate-redshift galaxies (z=1-5.8) expected for z>6 gala… ▽ More (Abridged) The Euclid mission is expected to discover thousands of z>6 galaxies in three Deep Fields, which together will cover a ~40 deg2 area. However, the limited number of Euclid bands and availability of ancillary data could make the identification of z>6 galaxies challenging. In this work, we assess the degree of contamination by intermediate-redshift galaxies (z=1-5.8) expected for z>6 galaxies within the Euclid Deep Survey. This study is based on ~176,000 real galaxies at z=1-8 in a ~0.7 deg2 area selected from the UltraVISTA ultra-deep survey, and ~96,000 mock galaxies with 25.3$\leq$H<27.0, which altogether cover the range of magnitudes to be probed in the Euclid Deep Survey. We simulate Euclid and ancillary photometry from the fiducial, 28-band photometry, and fit spectral energy distributions (SEDs) to various combinations of these simulated data. Our study demonstrates that identifying z>6 with Euclid data alone will be very effective, with a z>6 recovery of 91(88)% for bright (faint) galaxies. For the UltraVISTA-like bright sample, the percentage of z=1-5.8 contaminants amongst apparent z>6 galaxies as observed with Euclid alone is 18%, which is reduced to 4(13)% by including ultra-deep Rubin (Spitzer) photometry. Conversely, for the faint mock sample, the contamination fraction with Euclid alone is considerably higher at 39%, and minimized to 7% when including ultra-deep Rubin data. For UltraVISTA-like bright galaxies, we find that Euclid (I-Y)>2.8 and (Y-J)<1.4 colour criteria can separate contaminants from true z>6 galaxies, although these are applicable to only 54% of the contaminants, as many have unconstrained (I-Y) colours. In the most optimistic scenario, these cuts reduce the contamination fraction to 1% whilst preserving 81% of the fiducial z>6 sample. For the faint mock sample, colour cuts are infeasible. △ Less

Submitted 31 October, 2022; v1 submitted 5 May, 2022; originally announced May 2022.

Comments: 27 pages, 14 figures, 9 tables; version accepted for publication by A&A

Journal ref: A&A 666, A200 (2022)

Lightweight HI source finding for next generation radio surveys

Authors: Emma Tolley, Damien Korber, Aymeric Galan, Austin Peel, Mark T. Sargent, Jean-Paul Kneib, Frederic Courbin, Jean-Luc Starck

Abstract: Future deep HI surveys will be essential for understanding the nature of galaxies and the content of the Universe. However, the large volume of these data will require distributed and automated processing techniques. We introduce LiSA, a set of python modules for the denoising, detection and characterization of HI sources in 3D spectral data. LiSA was developed and tested on the Square Kilometer A… ▽ More Future deep HI surveys will be essential for understanding the nature of galaxies and the content of the Universe. However, the large volume of these data will require distributed and automated processing techniques. We introduce LiSA, a set of python modules for the denoising, detection and characterization of HI sources in 3D spectral data. LiSA was developed and tested on the Square Kilometer Array Science Data Challenge 2 dataset, and contains modules and pipelines for easy domain decomposition and parallel execution. LiSA contains algorithms for 2D-1D wavelet denoising using the starlet transform and flexible source finding using null-hypothesis testing. These algorithms are lightweight and portable, needing only a few user-defined parameters reflecting the resolution of the data. LiSA also includes two convolutional neural networks developed to analyse data cubes which separate HI sources from artifacts and predict the HI source properties. All of these components are designed to be as modular as possible, allowing users to mix and match different components to create their ideal pipeline. We demonstrate the performance of the different components of LiSA on the SDC2 dataset, which is able to find 95% of HI sources with SNR > 3 and accurately predict their properties. △ Less

Submitted 20 April, 2022; originally announced April 2022.

Comments: 18 pages, 16 figures

ShapePipe: a new shape measurement pipeline and weak-lensing application to UNIONS/CFIS data

Authors: Axel Guinot, Martin Kilbinger, Samuel Farrens, Austin Peel, Arnau Pujol, Morgan Schmitz, Jean-Luc Starck, Thomas Erben, Raphael Gavazzi, Stephen Gwyn, Michael J. Hudson, Hendrik Hiledebrandt, Tobias Liaudat, Lance Miller, Isaac Spitzer, Ludovic Van Waerbeke, Jean-Charles Cuillandre, Sébastien Fabbro, Alan McConnachie

Abstract: UNIONS is an ongoing collaboration that will provide the largest deep photometric survey of the Northern sky in four optical bands to date. As part of this collaboration, CFIS is taking $r$-band data with an average seeing of 0.65 arcsec, which is complete to magnitude 24.5 and thus ideal for weak-lensing studies. We perform the first weak-lensing analysis of CFIS $r$-band data over an area spanni… ▽ More UNIONS is an ongoing collaboration that will provide the largest deep photometric survey of the Northern sky in four optical bands to date. As part of this collaboration, CFIS is taking $r$-band data with an average seeing of 0.65 arcsec, which is complete to magnitude 24.5 and thus ideal for weak-lensing studies. We perform the first weak-lensing analysis of CFIS $r$-band data over an area spanning 1700 deg$^2$ of the sky. We create a catalogue with measured shapes for 40 million galaxies, corresponding to an effective density of 6.8 galaxies per square arcminute, and demonstrate a low level of systematic biases. This work serves as the basis for further cosmological studies using the full UNIONS survey of 4800 deg$^2$ when completed. Here we present ShapePipe, a newly developed weak-lensing pipeline. This pipeline makes use of state-of-the-art methods such as Ngmix for accurate galaxy shape measurement. Shear calibration is performed with metacalibration. We carry out extensive validation tests on the Point Spread Function (PSF), and on the galaxy shapes. In addition, we create realistic image simulations to validate the estimated shear. We quantify the PSF model accuracy and show that the level of systematics is low as measured by the PSF residuals. Their effect on the shear two-point correlation function is sub-dominant compared to the cosmological contribution on angular scales <100 arcmin. The additive shear bias is below 5x$10^{-4}$, and the residual multiplicative shear bias is at most $10^{-3}$ as measured on image simulations. Using COSEBIs we show that there are no significant B-modes present in second-order shear statistics. We present convergence maps and see clear correlations of the E-mode with known cluster positions. We measure the stacked tangential shear profile around Planck clusters at a significance higher than $4σしぐま$. △ Less

Submitted 10 April, 2022; originally announced April 2022.

Journal ref: A&A 666, A162 (2022)

Euclid preparation. XVIII. The NISP photometric system

Authors: Euclid Collaboration, M. Schirmer, K. Jahnke, G. Seidel, H. Aussel, C. Bodendorf, F. Grupp, F. Hormuth, S. Wachter, P. N. Appleton, R. Barbier, J. Brinchmann, J. M. Carrasco, F. J. Castander, J. Coupon, F. De Paolis, A. Franco, K. Ganga, P. Hudelot, E. Jullo, A. Lancon, A. A. Nucita, S. Paltani, G. Smadja, L. M. G. Venancio , et al. (198 additional authors not shown)

Abstract: Euclid will be the first space mission to survey most of the extragalactic sky in the 0.95-2.02 $μみゅー$m range, to a 5$σしぐま$ point-source median depth of 24.4 AB mag. This unique photometric data set will find wide use beyond Euclid's core science. In this paper, we present accurate computations of the Euclid Y_E, J_E and H_E passbands used by the Near-Infrared Spectrometer and Photometer (NISP), and the… ▽ More Euclid will be the first space mission to survey most of the extragalactic sky in the 0.95-2.02 $μみゅー$m range, to a 5$σしぐま$ point-source median depth of 24.4 AB mag. This unique photometric data set will find wide use beyond Euclid's core science. In this paper, we present accurate computations of the Euclid Y_E, J_E and H_E passbands used by the Near-Infrared Spectrometer and Photometer (NISP), and the associated photometric system. We pay particular attention to passband variations in the field of view, accounting among others for spatially variable filter transmission, and variations of the angle of incidence on the filter substrate using optical ray tracing. The response curves' cut-on and cut-off wavelengths - and their variation in the field of view - are determined with 0.8 nm accuracy, essential for the photometric redshift accuracy required by Euclid. After computing the photometric zeropoints in the AB mag system, we present linear transformations from and to common ground-based near-infrared photometric systems, for normal stars, red and brown dwarfs, and galaxies separately. A Python tool to compute accurate magnitudes for arbitrary passbands and spectral energy distributions is provided. We discuss various factors from space weathering to material outgassing that may slowly alter Euclid's spectral response. At the absolute flux scale, the Euclid in-flight calibration program connects the NISP photometric system to Hubble Space Telescope spectrophotometric white dwarf standards; at the relative flux scale, the chromatic evolution of the response is tracked at the milli-mag level. In this way, we establish an accurate photometric system that is fully controlled throughout Euclid's lifetime. △ Less

Submitted 31 March, 2022; v1 submitted 3 March, 2022; originally announced March 2022.

Comments: 33 pages, 25 figures, accepted for publication in A&A

Journal ref: A&A 662, A92 (2022)

Galaxy flows within 8,000 km/s from Numerical Action methods

Authors: Edward Shaya, R. Brent Tully, Daniel Pomarède, Alan Peel

Abstract: The trajectories since z=4 of systems of galaxies (`halos') with cz < 8,000 km/s are found through Numerical Action reconstructions. A set of 9,719 halos from a 2MASS group catalog and Cosmicflows-3 catalogs are given attention. Present distances are adjusted to minimize departures from observed redshifts. For those with the most precisely determined distances, compromises are made between distanc… ▽ More The trajectories since z=4 of systems of galaxies (`halos') with cz < 8,000 km/s are found through Numerical Action reconstructions. A set of 9,719 halos from a 2MASS group catalog and Cosmicflows-3 catalogs are given attention. Present distances are adjusted to minimize departures from observed redshifts. For those with the most precisely determined distances, compromises are made between distance and redshift agreement. $H_0$ is varied from 69 to 77 km s$^{-1}$ Mpc$^{-1}$ with $Ωおめが_m$ set by the baryon acoustic oscillation constraint from the Planck Satellite. A best fitting amplitude of the mass-to-light relation is found. A uniform density associated with the interhalo medium accounts for the matter not in halos. The solution paths provide the histories of the formation of the nearby large structures and depict how the voids emptied. Assuming no local over/underdensity, the best model has $H_0=73$ km s$^{-1}$ Mpc$^{-1}$ with nearly the same density arising from interhalo matter (IHM) as from halos. We examine local over/underdensities by varying the IHM density and find a valley of best fit models along $H_0 = 73.0 (1 + 0.165δでるた)$ km s$^{-1}$ Mpc$^{-1}$. Friedmann models with distinct densities internal and external to the study region give a similar relationship. The fraction of matter in the IHM seen in n-body simulations roughly matches that in our $H_0=72$ scenario. Videos have been created to visualize the complexities of formation of large-scale structures. Standard n-body calculations starting from the first time-steps as tests of the NAM solutions, and continue until cosmic scale factor $a=2$ provide glimpses into the future. △ Less

Submitted 28 January, 2022; originally announced January 2022.

Comments: 29 pages, 28 figures, 4 links to sketchfab.com, and 4 links to vimeo.com videos

Pool samples to efficiently estimate pathogen prevalence dynamics

Authors: Braden Scherting, Alison Peel, Raina Plowright, Andrew Hoegh

Abstract: Estimating the prevalence of a disease is necessary for evaluating and mitigating risks of its transmission within or between populations. Estimates that consider how prevalence changes with time provide more information about these risks but are difficult to obtain due to the necessary sampling intensity and commensurate testing costs. We propose pooling and jointly testing multiple samples to re… ▽ More Estimating the prevalence of a disease is necessary for evaluating and mitigating risks of its transmission within or between populations. Estimates that consider how prevalence changes with time provide more information about these risks but are difficult to obtain due to the necessary sampling intensity and commensurate testing costs. We propose pooling and jointly testing multiple samples to reduce testing costs and use a novel nonparametric, hierarchical Bayesian model to infer population prevalence from the pooled test results. This approach is shown to reduce uncertainty compared to individual testing at the same budget and to produce similar estimates compared to individual testing at a much higher budget through two synthetic studies and two case studies of natural infection data. △ Less

Submitted 11 November, 2021; originally announced November 2021.

Weak lensing mass reconstruction using sparsity and a Gaussian random field

Authors: J. -L. Starck, K. E. Themelis, N. Jeffrey, A. Peel, F. Lanusse

Abstract: We introduce a novel approach to reconstruct dark matter mass maps from weak gravitational lensing measurements. The cornerstone of the proposed method lies in a new modelling of the matter density field in the Universe as a mixture of two components:(1) a sparsity-based component that captures the non-Gaussian structure of the field, such as peaks or halos at different spatial scales; and (2) a G… ▽ More We introduce a novel approach to reconstruct dark matter mass maps from weak gravitational lensing measurements. The cornerstone of the proposed method lies in a new modelling of the matter density field in the Universe as a mixture of two components:(1) a sparsity-based component that captures the non-Gaussian structure of the field, such as peaks or halos at different spatial scales; and (2) a Gaussian random field, which is known to well represent the linear characteristics of the field.Methods. We propose an algorithm called MCALens which jointly estimates these two components. MCAlens is based on an alternating minimization incorporating both sparse recovery and a proximal iterative Wiener filtering. Experimental results on simulated data show that the proposed method exhibits improved estimation accuracy compared to state-of-the-art mass map reconstruction methods. △ Less

Submitted 8 February, 2021; originally announced February 2021.

Comments: A&A, accepted

Journal ref: A&A 649, A99 (2021)

SLITronomy: towards a fully wavelet-based strong lensing inversion technique

Authors: A. Galan, A. Peel, R. Joseph, F. Courbin, J. -L. Starck

Abstract: Strong gravitational lensing provides a wealth of astrophysical information on the baryonic and dark matter content of galaxies. It also serves as a valuable cosmological probe by allowing us to measure the Hubble constant independently of other methods. These applications all require the difficult task of inverting the lens equation and simultaneously reconstructing the mass profile of the lens a… ▽ More Strong gravitational lensing provides a wealth of astrophysical information on the baryonic and dark matter content of galaxies. It also serves as a valuable cosmological probe by allowing us to measure the Hubble constant independently of other methods. These applications all require the difficult task of inverting the lens equation and simultaneously reconstructing the mass profile of the lens along with the original light profile of the unlensed source. As there is no reason for either the lens or the source to be simple, we need methods that both invert the lens equation with a large number of degrees of freedom and also enforce a well-controlled regularisation that avoids the appearance of spurious structures. This can be beautifully accomplished by representing signals in wavelet space. Building on the Sparse Lens Inversion Technique (SLIT), in this work we present an improved sparsity-based method that describes lensed sources using wavelets and optimises over the parameters given an analytical lens mass profile. We apply our technique on simulated HST and E-ELT data, as well as on real HST images of lenses from the Sloan Lens ACS (SLACS) sample, assuming a lens model. We show that wavelets allow us to reconstruct lensed sources containing detailed substructures when using both present-day data and high-resolution images from future thirty-meter-class telescopes. Wavelets moreover provide a much more tractable solution in terms of quality and computation time compared to using a source model that combines smooth analytical profiles and shapelets. Requiring very little human interaction, our pixel-based technique fits into the effort to devise automated modelling schemes. It can be incorporated in the standard workflow of sampling analytical lens model parameters. The method, which we call SLITronomy, is freely available as a new plug-in to the modelling software Lenstronomy. △ Less

Submitted 4 December, 2020; originally announced December 2020.

Journal ref: A&A 647, A176 (2021)

H0LiCOW XI. A weak lensing measurement of the external convergence in the field of the lensed quasar B1608+656 using HST and Subaru deep imaging

Authors: O. Tihhonova, F. Courbin, D. Harvey, S. Hilbert, A. Peel, C. E. Rusu, C. D. Fassnacht, V. Bonvin, P. J. Marshall, G. Meylan, D. Sluse, S. H. Suyu, T. Treu, K. C. Wong

Abstract: We investigate the environment and line of sight of the H0LiCOW lens B1608+656 using Subaru Suprime-Cam and the Hubble Space Telescope (HST) to perform a weak lensing analysis. We compare three different methods to reconstruct the mass map of the field, i.e. the standard Kaiser-Squires inversion coupled with inpainting and Gaussian or wavelet filtering, and $\tt{Glimpse}$ a method based on sparse… ▽ More We investigate the environment and line of sight of the H0LiCOW lens B1608+656 using Subaru Suprime-Cam and the Hubble Space Telescope (HST) to perform a weak lensing analysis. We compare three different methods to reconstruct the mass map of the field, i.e. the standard Kaiser-Squires inversion coupled with inpainting and Gaussian or wavelet filtering, and $\tt{Glimpse}$ a method based on sparse regularization of the shear field. We find no substantial difference between the 2D mass reconstructions, but we find that the ground-based data is less sensitive to small-scale structures than the space-based observations. Marginalising over the results obtained with all the reconstruction techniques applied to the two available HST filters F606W and F814W, we estimate the external convergence, $κかっぱ_{\rm ext}$ at the position of B1608+656 is $κかっぱ_{\rm ext} = 0.11^{+0.06}_{-0.04}$, where the error bars corresponds respectively to the 16th and 84th quartiles. This result is compatible with previous estimates using the number-counts technique, suggesting that B1608+656 resides in an over-dense line of sight, but with a completely different technique. Using our mass reconstructions, we also compare the convergence at the position of several groups of galaxies in the field of B1608+656 with the mass measurements using various analytical mass profiles, and find that the weak lensing results favor truncated halo models. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: Accepted MNRAS

Constraining neutrino masses with weak-lensing multiscale peak counts

Authors: Virginia Ajani, Austin Peel, Valeria Pettorino, Jean-Luc Starck, Zack Li, Jia Liu

Abstract: Massive neutrinos influence the background evolution of the Universe as well as the growth of structure. Being able to model this effect and constrain the sum of their masses is one of the key challenges in modern cosmology. Weak-lensing cosmological constraints will also soon reach higher levels of precision with next-generation surveys like LSST, WFIRST and Euclid. We use the MassiveNus simulati… ▽ More Massive neutrinos influence the background evolution of the Universe as well as the growth of structure. Being able to model this effect and constrain the sum of their masses is one of the key challenges in modern cosmology. Weak-lensing cosmological constraints will also soon reach higher levels of precision with next-generation surveys like LSST, WFIRST and Euclid. We use the MassiveNus simulations to derive constraints on the sum of neutrino masses $M_νにゅー$, the present-day total matter density $Ωおめが_{\rm m}$, and the primordial power spectrum normalization $A_{\rm s}$ in a tomographic setting. We measure the lensing power spectrum as second-order statistics along with peak counts as higher-order statistics on lensing convergence maps generated from the simulations. We investigate the impact of multiscale filtering approaches on cosmological parameters by employing a starlet (wavelet) filter and a concatenation of Gaussian filters. In both cases peak counts perform better than the power spectrum on the set of parameters [$M_νにゅー$, $Ωおめが_{\rm m}$, $A_{\rm s}$] respectively by 63$\%$, 40$\%$ and 72$\%$ when using a starlet filter and by 70$\%$, 40$\%$ and 77$\%$ when using a multiscale Gaussian. More importantly, we show that when using a multiscale approach, joining power spectrum and peaks does not add any relevant information over considering just the peaks alone. While both multiscale filters behave similarly, we find that with the starlet filter the majority of the information in the data covariance matrix is encoded in the diagonal elements; this can be an advantage when inverting the matrix, speeding up the numerical implementation. △ Less

Submitted 5 January, 2021; v1 submitted 29 January, 2020; originally announced January 2020.

Comments: 17 pages, 10 figures, 6 tables

Journal ref: Phys. Rev. D 102, 103531, (2020)

The impact of baryonic physics and massive neutrinos on weak lensing peak statistics

Authors: Matthew Fong, Miyoung Choi, Victoria Catlett, Brandyn Lee, Austin Peel, Rachel Bowyer, Lindsay J. King, Ian G. McCarthy

Abstract: We study the impact of baryonic processes and massive neutrinos on weak lensing peak statistics that can be used to constrain cosmological parameters. We use the BAHAMAS suite of cosmological simulations, which self-consistently include baryonic processes and the effect of massive neutrino free-streaming on the evolution of structure formation. We construct synthetic weak lensing catalogues by ray… ▽ More We study the impact of baryonic processes and massive neutrinos on weak lensing peak statistics that can be used to constrain cosmological parameters. We use the BAHAMAS suite of cosmological simulations, which self-consistently include baryonic processes and the effect of massive neutrino free-streaming on the evolution of structure formation. We construct synthetic weak lensing catalogues by ray-tracing through light-cones, and use the aperture mass statistic for the analysis. The peaks detected on the maps reflect the cumulative signal from massive bound objects and general large-scale structure. We present the first study of weak lensing peaks in simulations that include both baryonic physics and massive neutrinos (summed neutrino mass $M_νにゅー =$ 0.06, 0.12, 0.24, and 0.48 eV assuming normal hierarchy), so that the uncertainty due to physics beyond the gravity of dark matter can be factored into constraints on cosmological models. Assuming a fiducial model of baryonic physics, we also investigate the correlation between peaks and massive haloes, over a range of summed neutrino mass values. As higher neutrino mass tends to suppress the formation of massive structures in the Universe, the halo mass function and lensing peak counts are therefore modified as a function of $M_νにゅー$. Over most of the S/N range, the impact of fiducial baryonic physics is greater (less) than neutrinos for 0.06 and 0.12 (0.24 and 0.48) eV models. Both baryonic physics and massive neutrinos should be accounted for when deriving cosmological parameters from weak lensing observations. △ Less

Submitted 6 October, 2019; v1 submitted 6 July, 2019; originally announced July 2019.

Journal ref: MNRAS 488 (2019) 3340

The Role of Machine Learning in the Next Decade of Cosmology

Authors: Michelle Ntampaka, Camille Avestruz, Steven Boada, Joao Caldeira, Jessi Cisewski-Kehe, Rosanne Di Stefano, Cora Dvorkin, August E. Evrard, Arya Farahi, Doug Finkbeiner, Shy Genel, Alyssa Goodman, Andy Goulding, Shirley Ho, Arthur Kosowsky, Paul La Plante, Francois Lanusse, Michelle Lochner, Rachel Mandelbaum, Daisuke Nagai, Jeffrey A. Newman, Brian Nord, J. E. G. Peek, Austin Peel, Barnabas Poczos , et al. (5 additional authors not shown)

Abstract: In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster an… ▽ More In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors. △ Less

Submitted 14 January, 2021; v1 submitted 26 February, 2019; originally announced February 2019.

Comments: Submitted to the Astro2020 call for science white papers

Distinguishing standard and modified gravity cosmologies with machine learning

Authors: Austin Peel, Florian Lalande, Jean-Luc Starck, Valeria Pettorino, Julian Merten, Carlo Giocoli, Massimo Meneghetti, Marco Baldi

Abstract: We present a convolutional neural network to classify distinct cosmological scenarios based on the statistically similar weak-lensing maps they generate. Modified gravity (MG) models that include massive neutrinos can mimic the standard concordance model ($Λらむだ$CDM) in terms of Gaussian weak-lensing observables. An inability to distinguish viable models that are based on different physics potentially… ▽ More We present a convolutional neural network to classify distinct cosmological scenarios based on the statistically similar weak-lensing maps they generate. Modified gravity (MG) models that include massive neutrinos can mimic the standard concordance model ($Λらむだ$CDM) in terms of Gaussian weak-lensing observables. An inability to distinguish viable models that are based on different physics potentially limits a deeper understanding of the fundamental nature of cosmic acceleration. For a fixed redshift of sources, we demonstrate that a machine learning network trained on simulated convergence maps can discriminate between such models better than conventional higher-order statistics. Results improve further when multiple source redshifts are combined. To accelerate training, we implement a novel data compression strategy that incorporates our prior knowledge of the morphology of typical convergence map features. Our method fully distinguishes $Λらむだ$CDM from its most similar MG model on noise-free data, and it correctly identifies among the MG models with at least 80% accuracy when using the full redshift information. Adding noise lowers the correct classification rate of all models, but the neural network still significantly outperforms the peak statistics used in a previous analysis. △ Less

Submitted 23 May, 2019; v1 submitted 25 October, 2018; originally announced October 2018.

Comments: 12 pages, 10 figures, 1 table. Text updated and expanded throughout

Journal ref: Phys. Rev. D 100, 023508 (2019)

On the dissection of degenerate cosmologies with machine learning

Authors: Julian Merten, Carlo Giocoli, Marco Baldi, Massimo Meneghetti, Austin Peel, Florian Lalande, Jean-Luc Starck, Valeria Pettorino

Abstract: Based on the DUSTGRAIN-pathfinder suite of simulations, we investigate observational degeneracies between nine models of modified gravity and massive neutrinos. Three types of machine learning techniques are tested for their ability to discriminate lensing convergence maps by extracting dimensional reduced representations of the data. Classical map descriptors such as the power spectrum, peak coun… ▽ More Based on the DUSTGRAIN-pathfinder suite of simulations, we investigate observational degeneracies between nine models of modified gravity and massive neutrinos. Three types of machine learning techniques are tested for their ability to discriminate lensing convergence maps by extracting dimensional reduced representations of the data. Classical map descriptors such as the power spectrum, peak counts and Minkowski functionals are combined into a joint feature vector and compared to the descriptors and statistics that are common to the field of digital image processing. To learn new features directly from the data we use a Convolutional Neural Network (CNN). For the mapping between feature vectors and the predictions of their underlying model, we implement two different classifiers; one based on a nearest-neighbour search and one that is based on a fully connected neural network. We find that the neural network provides a much more robust classification than the nearest-neighbour approach and that the CNN provides the most discriminating representation of the data. It achieves the cleanest separation between the different models and the highest classification success rate of 59% for a single source redshift. Once we perform a tomographic CNN analysis, the total classification accuracy increases significantly to 76% with no observational degeneracies remaining. Visualising the filter responses of the CNN at different network depths provides us with the unique opportunity to learn from very complex models and to understand better why they perform so well. △ Less

Submitted 27 March, 2019; v1 submitted 25 October, 2018; originally announced October 2018.

Comments: 21 pages, 14 figures, 10 tables. Associated code and data respository at https://www.bitbucket.org/jmerten82/mydnn . Accepted for publication by the MNRAS

Breaking degeneracies in modified gravity with higher (than 2nd) order weak-lensing statistics

Authors: Austin Peel, Valeria Pettorino, Carlo Giocoli, Jean-Luc Starck, Marco Baldi

Abstract: General relativity (GR) has been well tested up to solar system scales, but it is much less certain that standard gravity remains an accurate description on the largest, that is, cosmological, scales. Many extensions to GR have been studied that are not yet ruled out by the data, including by that of the recent direct gravitational wave detections. Degeneracies among the standard model ($Λらむだ$CDM) an… ▽ More General relativity (GR) has been well tested up to solar system scales, but it is much less certain that standard gravity remains an accurate description on the largest, that is, cosmological, scales. Many extensions to GR have been studied that are not yet ruled out by the data, including by that of the recent direct gravitational wave detections. Degeneracies among the standard model ($Λらむだ$CDM) and modified gravity (MG) models, as well as among different MG parameters, must be addressed in order to best exploit information from current and future surveys and to unveil the nature of dark energy. We propose various higher-order statistics in the weak-lensing signal as a new set of observables able to break degeneracies between massive neutrinos and MG parameters. We have tested our methodology on so-called $f(R)$ models, which constitute a class of viable models that can explain the accelerated universal expansion by a modification of the fundamental gravitational interaction. We have explored a range of these models that still fit current observations at the background and linear level, and we show using numerical simulations that certain models which include massive neutrinos are able to mimic $Λらむだ$CDM in terms of the 3D power spectrum of matter density fluctuations. We find that depending on the redshift and angular scale of observation, non-Gaussian information accessed by higher-order weak-lensing statistics can be used to break the degeneracy between $f(R)$ models and $Λらむだ$CDM. In particular, peak counts computed in aperture mass maps outperform third- and fourth-order moments. △ Less

Submitted 10 October, 2018; v1 submitted 14 May, 2018; originally announced May 2018.

Comments: 19 pages, 13 figures, 2 tables; appendix and references added; matches published version

Journal ref: A&A 619, A38 (2018)

Sparse reconstruction of the merging A520 cluster system

Authors: Austin Peel, François Lanusse, Jean-Luc Starck

Abstract: Merging galaxy clusters present a unique opportunity to study the properties of dark matter in an astrophysical context. These are rare and extreme cosmic events in which the bulk of the baryonic matter becomes displaced from the dark matter halos of the colliding subclusters. Since all mass bends light, weak gravitational lensing is a primary tool to study the total mass distribution in such syst… ▽ More Merging galaxy clusters present a unique opportunity to study the properties of dark matter in an astrophysical context. These are rare and extreme cosmic events in which the bulk of the baryonic matter becomes displaced from the dark matter halos of the colliding subclusters. Since all mass bends light, weak gravitational lensing is a primary tool to study the total mass distribution in such systems. Combined with X-ray and optical analyses, mass maps of cluster mergers reconstructed from weak-lensing observations have been used to constrain the self-interaction cross-section of dark matter. The dynamically complex Abell 520 (A520) cluster is an exceptional case, even among merging systems: multi-wavelength observations have revealed a surprising high mass-to-light concentration of dark mass, the interpretation of which is difficult under the standard assumption of effectively collisionless dark matter. We revisit A520 using a new sparsity-based mass-mapping algorithm to independently assess the presence of the puzzling dark core. We obtain high-resolution mass reconstructions from two separate galaxy shape catalogs derived from Hubble Space Telescope observations of the system. Our mass maps agree well overall with the results of previous studies, but we find important differences. In particular, although we are able to identify the dark core at a certain level in both data sets, it is at much lower significance than has been reported before using the same data. As we cannot confirm the detection in our analysis, we do not consider A520 as posing a significant challenge to the collisionless dark matter scenario. △ Less

Submitted 1 August, 2017; originally announced August 2017.

Comments: 13 pages, 8 figures, submitted to ApJ

Cosmological constraints with weak lensing peak counts and second-order statistics in a large-field survey

Authors: Austin Peel, Chieh-An Lin, Francois Lanusse, Adrienne Leonard, Jean-Luc Starck, Martin Kilbinger

Abstract: Peak statistics in weak lensing maps access the non-Gaussian information contained in the large-scale distribution of matter in the Universe. They are therefore a promising complement to two-point and higher-order statistics to constrain our cosmological models. To prepare for the high-precision data of next-generation surveys, we assess the constraining power of peak counts in a simulated Euclid-… ▽ More Peak statistics in weak lensing maps access the non-Gaussian information contained in the large-scale distribution of matter in the Universe. They are therefore a promising complement to two-point and higher-order statistics to constrain our cosmological models. To prepare for the high-precision data of next-generation surveys, we assess the constraining power of peak counts in a simulated Euclid-like survey on the cosmological parameters $Ωおめが_\mathrm{m}$, $σしぐま_8$, and $w_0^\mathrm{de}$. In particular, we study how the Camelus model--a fast stochastic algorithm for predicting peaks--can be applied to such large surveys. We measure the peak count abundance in a mock shear catalogue of ~5,000 sq. deg. using a multiscale mass map filtering technique. We then constrain the parameters of the mock survey using Camelus combined with approximate Bayesian computation (ABC). We find that peak statistics yield a tight but significantly biased constraint in the $σしぐま_8$-$Ωおめが_\mathrm{m}$ plane, indicating the need to better understand and control the model's systematics. We calibrate the model to remove the bias and compare results to those from the two-point correlation functions (2PCF) measured on the same field. In this case, we find the derived parameter $Σしぐま_8=σしぐま_8(Ωおめが_\mathrm{m}/0.27)^αあるふぁ=0.76_{-0.03}^{+0.02}$ with $αあるふぁ=0.65$ for peaks, while for 2PCF the value is $Σしぐま_8=0.76_{-0.01}^{+0.02}$ with $αあるふぁ=0.70$. We therefore see comparable constraining power between the two probes, and the offset of their $σしぐま_8$-$Ωおめが_\mathrm{m}$ degeneracy directions suggests that a combined analysis would yield tighter constraints than either measure alone. As expected, $w_0^\mathrm{de}$ cannot be well constrained without a tomographic analysis, but its degeneracy directions with the other two varied parameters are still clear for both peaks and 2PCF. (abridged) △ Less

Submitted 7 December, 2016; originally announced December 2016.

Comments: 12 pages, 10 figures, submitted to A&A

Journal ref: A&A 599, A79 (2017)

Effect of inhomogeneities on high precision measurements of cosmological distances

Authors: Austin Peel, M. A. Troxel, Mustapha Ishak

Abstract: We study effects of inhomogeneities on distance measures in an exact relativistic Swiss-cheese model of the universe, focusing on the distance modulus. The model has LCDM background dynamics, and the `holes' are non-symmetric structures described by the Szekeres metric. The Szekeres exact solution of Einstein's equations, which is inhomogeneous and anisotropic, allows us to capture potentially rel… ▽ More We study effects of inhomogeneities on distance measures in an exact relativistic Swiss-cheese model of the universe, focusing on the distance modulus. The model has LCDM background dynamics, and the `holes' are non-symmetric structures described by the Szekeres metric. The Szekeres exact solution of Einstein's equations, which is inhomogeneous and anisotropic, allows us to capture potentially relevant effects on light propagation due to nontrivial evolution of structures in an exact framework. Light beams traversing a single Szekeres structure in different ways can experience either magnification or demagnification, depending on the particular path, and we explore a small additional effect due time evolution of the structures. We study the probability distributions of $Δでるたμみゅー=μみゅー_{LCDM}-μみゅー_{SC}$ for sources at different redshifts in various Swiss-cheese constructions, where the light beams travel through a large number of randomly oriented holes with random impact parameters. We find for $Δでるたμみゅー$ the dispersions $0.004\le σしぐま_{Δでるたμみゅー} \le 0.008$ mag for sources with redshifts $1.0\le z \le 1.5$, which are smaller than the intrinsic dispersion of, for example, magnitudes of type Ia supernovae. The shapes of the distributions we obtain for our Swiss-cheese constructions are peculiar in the sense that they are not consistently skewed toward the demagnification side, as they are in analyses of lensing in cosmological simulations, reflecting a limitation of these constructions. This could be the result of requiring the continuity of Einstein's equations throughout the overall spacetime patchwork, which imposes the condition that compensating overdense shells must accompany the underdense void regions in the holes. The possibility to explore other uses of these constructions that could circumvent this limitation and lead to different statistics remains open. (Abridged) △ Less

Submitted 5 January, 2015; v1 submitted 19 August, 2014; originally announced August 2014.

Comments: 18 pages, 8 figures, revised conclusion, matches version published in PRD

Journal ref: Phys. Rev. D 90, 123536 (2014)

The effects of structure anisotropy on lensing observables in an exact general relativistic setting for precision cosmology

Authors: M. A. Troxel, Mustapha Ishak, Austin Peel

Abstract: The study of relativistic, higher order and nonlinear effects has become necessary in recent years in the pursuit of precision cosmology. We develop and apply here a framework to study gravitational lensing in exact models in general relativity that are not restricted to homogeneity and isotropy, and where full nonlinearity and relativistic effects are included. We apply the framework to a specifi… ▽ More The study of relativistic, higher order and nonlinear effects has become necessary in recent years in the pursuit of precision cosmology. We develop and apply here a framework to study gravitational lensing in exact models in general relativity that are not restricted to homogeneity and isotropy, and where full nonlinearity and relativistic effects are included. We apply the framework to a specific, anisotropic galaxy cluster model which is based on a modified NFW halo density profile and described by the Szekeres metric. We examine the effects of increasing levels of anisotropy in the galaxy cluster on lensing observables like the convergence and shear for various lensing geometries, finding a strong nonlinear response in both the convergence and shear for rays passing through anisotropic regions of the cluster. Deviation from the expected values in a spherically symmetric structure are asymmetric with respect to path direction and thus will persist as a statistical effect when averaged over some ensemble of such clusters. The resulting relative difference in various geometries can be as large as approximately 2%, 8%, and 24% in the measure of convergence for levels of anisotropy of 5%, 10%, and 15%, respectively, as a fraction of total cluster mass. For the total magnitude of shear, the relative difference can grow near the center of the structure to be as large as 15%, 32%, and 44% for the same levels of anisotropy, averaged over the two extreme geometries. The convergence is impacted most strongly for rays which pass in directions along the axis of maximum dipole anisotropy in the structure, while the shear is most strongly impacted for rays which pass in directions orthogonal to this axis, as expected. These effects due to anisotropic structures will affect lensing measurements and must be fully examined in an era of precision cosmology. △ Less

Submitted 9 April, 2014; v1 submitted 22 November, 2013; originally announced November 2013.

Comments: 26 pages, 8 figures, matches version published in JCAP

Journal ref: JCAP 1403:040, 2014

Effects of anisotropy on gravitational infall in galaxy clusters using an exact general relativistic model

Authors: M. A. Troxel, Austin Peel, Mustapha Ishak

Abstract: We study the effects and implications of anisotropies at the scale of galaxy clusters by building an exact general relativistic model of a cluster using the inhomogeneous and anisotropic Szekeres metric. The model is built from a modified Navarro-Frenk-White (NFW) density profile. We compare this to a corresponding spherically symmetric structure in the Lemaitre-Tolman (LT) model and quantify the… ▽ More We study the effects and implications of anisotropies at the scale of galaxy clusters by building an exact general relativistic model of a cluster using the inhomogeneous and anisotropic Szekeres metric. The model is built from a modified Navarro-Frenk-White (NFW) density profile. We compare this to a corresponding spherically symmetric structure in the Lemaitre-Tolman (LT) model and quantify the impact of introducing varying levels of anisotropy. We examine two physical measures of gravitational infall -- the growth rate of density and the velocity of the source dust in the model. We introduce a generalization of the LT dust velocity profile for the Szekeres metric and demonstrate its consistency with the growth rate of density. We find that the growth rate of density in one substructure increases by 0.5%, 1.5%, and 3.75% for 5%, 10%, and 15% levels of introduced anisotropy, which is measured as the fractional displaced mass relative to the spherically symmetric case. The infall velocity of the dust is found to increase by 2.5, 10, and 20 km/s (0.5%, 2%, and 4.5%), respectively, for the same three levels of anisotropy. This response to the anisotropy in a structure is found to be strongly nonlinear with respect to the strength of anisotropy. These relative velocities correspond to an equivalent increase in the total mass of the spherically symmetric structure of 1%, 3.8%, and 8.4%, indicating that not accounting for the presence of anisotropic mass distributions in cluster models can strongly bias the determination of physical properties like the total mass. △ Less

Submitted 27 January, 2014; v1 submitted 21 November, 2013; originally announced November 2013.

Comments: 20 pages, 6 figures, matches version published in JCAP

Journal ref: JCAP 1312:048, 2013

Stringent Restriction from the Growth of Large-Scale Structure on Apparent Acceleration in Inhomogeneous Cosmological Models

Authors: Mustapha Ishak, Austin Peel, M. A. Troxel

Abstract: Probes of cosmic expansion constitute the main basis for arguments to support or refute a possible apparent acceleration due to different expansion rates in the universe as described by inhomogeneous cosmological models. We present in this Letter a separate argument based on results from an analysis of the growth rate of large-scale structure in the universe as modeled by the inhomogeneous cosmolo… ▽ More Probes of cosmic expansion constitute the main basis for arguments to support or refute a possible apparent acceleration due to different expansion rates in the universe as described by inhomogeneous cosmological models. We present in this Letter a separate argument based on results from an analysis of the growth rate of large-scale structure in the universe as modeled by the inhomogeneous cosmological models of Szekeres. We use the models with no assumptions of spherical or axial symmetries. We find that while the Szekeres models can fit very well the observed expansion history without a $Λらむだ$, they fail to produce the observed late-time suppression in the growth unless $Λらむだ$ is added to the dynamics. A simultaneous fit to the supernova and growth factor data shows that the cold dark matter model with a cosmological constant ($Λらむだ$CDM) provides consistency with the data at a confidence level of 99.65% while the Szekeres model without $Λらむだ$ achieves only a 60.46% level. When the data sets are considered separately, the Szekeres with no $Λらむだ$ fits the supernova data as well as the $Λらむだ$CDM does, but provides a very poor fit to the growth data with only 31.31% consistency level compared to 99.99% for the $Λらむだ$CDM. This absence of late-time growth suppression in inhomogeneous models without a $Λらむだ$ is consolidated by a physical explanation. △ Less

Submitted 27 January, 2014; v1 submitted 2 July, 2013; originally announced July 2013.

Comments: 5 pages, 1 figure, matches version published in PRL

Journal ref: Phys. Rev. Lett. 111, 251302, 2013

Large-scale growth evolution in the Szekeres inhomogeneous cosmological models with comparison to growth data

Authors: Austin Peel, Mustapha Ishak, M. A. Troxel

Abstract: We use the Szekeres inhomogeneous cosmological models to study the growth of large-scale structure in the universe including nonzero spatial curvature and a cosmological constant. In particular, we use the Goode and Wainwright formulation, as in this form the models can be considered to represent exact nonlinear perturbations of an averaged background. We identify a density contrast in both classe… ▽ More We use the Szekeres inhomogeneous cosmological models to study the growth of large-scale structure in the universe including nonzero spatial curvature and a cosmological constant. In particular, we use the Goode and Wainwright formulation, as in this form the models can be considered to represent exact nonlinear perturbations of an averaged background. We identify a density contrast in both classes I and II of the models, for which we derive growth evolution equations. By including Lambda, the time evolution of the density contrast as well as kinematic quantities can be tracked through the matter- and Lambda-dominated cosmic eras up to the present and into the future. In various models of class I and class II, the growth rate is found to be stronger than that of the LCDM cosmology, and it is suppressed at later times due to the presence of Lambda. We find that there are Szekeres models able to provide a growth history similar to that of LCDM while requiring less matter content and nonzero curvature, which speaks to the importance of including the effects of large-scale inhomogeneities in analyzing the growth of large-scale structure. Using data for the growth factor f from redshift space distortions and the Lyman-alpha forest, we obtain best fit parameters for class II models and compare their ability to match observations with LCDM. We find that there is negligible difference between best fit Szekeres models with no priors and those for LCDM, both including and excluding Lyman-alpha data. We also find that the growth index gamma parametrization cannot be applied in a simple way to the growth in Szekeres models, so a direct comparison of the function f to the data is performed. We conclude that the Szekeres models can provide an exact framework for the analysis of large-scale growth data that includes inhomogeneities and allows for different interpretations of observations. (abridged) △ Less

Submitted 10 December, 2012; originally announced December 2012.

Comments: 25 pages, 13 figures

Journal ref: Phys.Rev.D86,123508,2012

The growth of structure in the Szekeres inhomogeneous cosmological models and the matter-dominated era

Authors: Mustapha Ishak, Austin Peel

Abstract: This study belongs to a series devoted to using Szekeres inhomogeneous models to develop a theoretical framework where observations can be investigated with a wider range of possible interpretations. We look here into the growth of large-scale structure in the models. The Szekeres models are exact solutions to Einstein's equations that were originally derived with no symmetries. We use a formulati… ▽ More This study belongs to a series devoted to using Szekeres inhomogeneous models to develop a theoretical framework where observations can be investigated with a wider range of possible interpretations. We look here into the growth of large-scale structure in the models. The Szekeres models are exact solutions to Einstein's equations that were originally derived with no symmetries. We use a formulation of the models that is due to Goode and Wainwright, who considered the models as exact perturbations of an FLRW background. Using the Raychaudhuri equation, we write for the two classes of the models, exact growth equations in terms of the under/overdensity and measurable cosmological parameters. The new equations in the overdensity split into two informative parts. The first part, while exact, is identical to the growth equation in the usual linearly perturbed FLRW models, while the second part constitutes exact non-linear perturbations. We integrate numerically the full exact growth rate equations for the flat and curved cases. We find that for the matter-dominated era, the Szekeres growth rate is up to a factor of three to five stronger than the usual linearly perturbed FLRW cases, reflecting the effect of exact Szekeres non-linear perturbations. The growth is also stronger than that of the non-linear spherical collapse model, and the difference between the two increases with time. This highlights the distinction when we use general inhomogeneous models where shear and a tidal gravitational field are present and contribute to the gravitational clustering. Additionally, it is worth observing that the enhancement of the growth found in the Szekeres models during the matter-dominated era could suggest a substitute to the argument that dark matter is needed when using FLRW models to explain the enhanced growth and resulting large-scale structures that we observe today (abridged) △ Less

Submitted 24 February, 2012; v1 submitted 13 April, 2011; originally announced April 2011.

Comments: 18 pages, 4 figures, matches PRD accepted version

Journal ref: Phys. Rev. D85, 083502, 2012

Delay before synchronization and its role in latency of sensory awareness

Authors: Adele Peel, Henrik Jeldtoft Jensen

Abstract: Here we show that for coupled-map systems, the length of the transient prior to synchronization is both dependant on the coupling strength and dynamics of connections: systems with fixed connections and with no self-coupling display quasi-instantaneous synchronization. Too strong tendency for synchronization would in terms of brain dynamics be expected to be a pathological case. We relate how th… ▽ More Here we show that for coupled-map systems, the length of the transient prior to synchronization is both dependant on the coupling strength and dynamics of connections: systems with fixed connections and with no self-coupling display quasi-instantaneous synchronization. Too strong tendency for synchronization would in terms of brain dynamics be expected to be a pathological case. We relate how the time to synchrony depends on coupling strength and connection dynamics to the latency between neuronal stimulation and conscious awareness. We suggest that this latency can be identified with the delay before a threshold level of synchrony is achieved between distinct regions within the brain, as suggested by recent empirical evidence, in which case the latency can easily be understood as the inevitable delay before such synchrony builds-up. This is demonstrated here through the study of simplistic coupled-map models. △ Less

Submitted 9 November, 2009; originally announced November 2009.

Comments: 5 pages and 4 figures

Our Peculiar Motion Away from the Local Void

Authors: R. Brent Tully, Edward J. Shaya, Igor D. Karachentsev, Helene M. Courtois, Dale D. Kocevski, Luca Rizzi, Alan Peel

Abstract: The peculiar velocity of the Local Group of galaxies manifested in the Cosmic Microwave Background dipole is found to decompose into three dominant components. The three components are clearly separated because they arise on distinct spatial scales and are fortuitously almost orthogonal in their influences. The nearest, which is distinguished by a velocity discontinuity at ~7 Mpc, arises from th… ▽ More The peculiar velocity of the Local Group of galaxies manifested in the Cosmic Microwave Background dipole is found to decompose into three dominant components. The three components are clearly separated because they arise on distinct spatial scales and are fortuitously almost orthogonal in their influences. The nearest, which is distinguished by a velocity discontinuity at ~7 Mpc, arises from the evacuation of the Local Void. We lie in the Local Sheet that bounds the void. Random motions within the Local Sheet are small. Our Galaxy participates in the bulk motion of the Local Sheet away from the Local Void. The component of our motion on an intermediate scale is attributed to the Virgo Cluster and its surroundings, 17 Mpc away. The third and largest component is an attraction on scales larger than 3000 km/s and centered near the direction of the Centaurus Cluster. The amplitudes of the three components are 259, 185, and 455 km/s, respectively, adding collectively to 631 km/s in the reference frame of the Local Sheet. Taking the nearby influences into account causes the residual attributed to large scales to align with observed concentrations of distant galaxies and reduces somewhat the amplitude of motion attributed to their pull. On small scales, in addition to the motion of our Local Sheet away from the Local Void, the nearest adjacent filament, the Leo Spur, is seen to be moving in a direction that will lead to convergence with our filament. Finally, a good distance to an isolated galaxy within the Local Void reveals that this dwarf system has a motion of at least 230 km/s away from the void center. Given the velocities expected from gravitational instability theory in the standard cosmological paradigm, the distance to the center of the Local Void must be at least 23 Mpc from our position. The Local Void is large! △ Less

Submitted 13 December, 2007; v1 submitted 28 May, 2007; originally announced May 2007.

Comments: Tentatively scheduled for Astrophysical Journal, 676 (March 20), 2008. 18 figures, 3 tables including web link for 2 tables, web links to 2 videos

Journal ref: Astrophys.J.676:184-205,2008

Sensitivity of Global Dynamics on the microscopic details of a network of dynamically coupled maps

Authors: Adele Peel, Henrik Jeldtoft Jensen

Abstract: Here we analyze the behavior of dynamically coupled maps, based on those introduced in a series of papers by Ito and Kaneko (Phys. Rev. Lett. 88:028701, 2002 and Phys. Rev. E 67:046226, 2003). We show how the microscopic coupling mechanism changes the behavior of the system both by affecting the stability of fixed points and through a more subtle effect in the crossover behavior between differen… ▽ More Here we analyze the behavior of dynamically coupled maps, based on those introduced in a series of papers by Ito and Kaneko (Phys. Rev. Lett. 88:028701, 2002 and Phys. Rev. E 67:046226, 2003). We show how the microscopic coupling mechanism changes the behavior of the system both by affecting the stability of fixed points and through a more subtle effect in the crossover behavior between different regions of the parameter space. This makes it necessary to choose very carefully the exact manner in which one couples maps if they are to be used as a general model of composite systems. △ Less

Submitted 30 March, 2007; originally announced March 2007.

Comments: 5 pages, 9 figures

Mass Selection Bias in Galaxy Cluster Peculiar Velocities from the Kinetic Sunyaev-Zel'dovich Effect

Authors: Alan C. Peel

Abstract: Upcoming surveys for galaxy clusters using the Sunyaev-Zel'dovich effect are potentially sensitive enough to create a peculiar velocity catalog. The statistics of these peculiar velocities are sensitive to cosmological parameters. We develop a method to explore parameter space using N-body simulations in order to quantify dark matter halo velocity statistics which will be useful for cluster pecu… ▽ More Upcoming surveys for galaxy clusters using the Sunyaev-Zel'dovich effect are potentially sensitive enough to create a peculiar velocity catalog. The statistics of these peculiar velocities are sensitive to cosmological parameters. We develop a method to explore parameter space using N-body simulations in order to quantify dark matter halo velocity statistics which will be useful for cluster peculiar velocity observations. We show that mass selection bias from a kinetic Sunyaev-Zel'dovich velocity catalog forecasts rms peculiar velocities with a much more complicated Omega_m dependency than suggested by linear theory. In addition, we show that both two-point functions for velocities disagree with linear theory predictions out to separations of approx. 40 Mpc/h. A pedagogical appendix is included developing linear theory notation with respect to the two-point peculiar velocity functions. △ Less

Submitted 28 July, 2006; v1 submitted 6 January, 2005; originally announced January 2005.

Comments: LaTeX, 11 pages, 10 figures. Changed to match journal paper. Useful pedagogic appendix added for calculations of biased velocity two point functions

Report number: DAMTP-2005-1

Journal ref: Mon.Not.Roy.Astron.Soc. 365 (2006) 1191-1202

COMPASS: An Upper Limit on CMB Polarization at an Angular Scale of 20 arc minutes

Authors: Philip C. Farese, Giorgio Dall'Oglio, Joshua O. Gundersen, Brian G. Keating, Slade Klawikowski, Lloyd Knox, Alan Levy, Philip M. Lubin, Chris W. O'Dell, Alan Peel, Lucio Piccirillo, John Ruhl, Peter T. Timbie

Abstract: COMPASS is an on-axis 2.6 meter telescope coupled to a correlation polarimeter operating at a wavelength of 1 cm. The entire instrument was built specifically for CMB polarization studies. We report here on observations of February 2001 - April 2001 using this system. We set an upper limit on E-mode polarized anisotropies of 33.5 uK (95% confidence limit) in the l-range 200-600. COMPASS is an on-axis 2.6 meter telescope coupled to a correlation polarimeter operating at a wavelength of 1 cm. The entire instrument was built specifically for CMB polarization studies. We report here on observations of February 2001 - April 2001 using this system. We set an upper limit on E-mode polarized anisotropies of 33.5 uK (95% confidence limit) in the l-range 200-600. △ Less

Submitted 19 August, 2003; originally announced August 2003.

Comments: Submitted to ApJ 8/19/2003

Journal ref: Astrophys.J.610:625-634,2004

COMPASS: An Instrument for Measuring the Polarization of the CMB on Intermediate Angular Scales

Authors: Philip C. Farese, Giorgio Dall'Oglio, Josh Gundersen, Brian Keating, Slade Klawikowski, Lloyd Knox, Alan Levy, Chris O'Dell, Alan Peel, Lucio Piccirillo, John Ruhl, Peter Timbie

Abstract: COMPASS is an on-axis 2.6 meter telescope coupled to a correlation polarimeter. The entire instrument was built specifically for CMB polarization studies. Careful attention was given to receiver and optics design, stability of the pointing platform, avoidance of systematic offsets, and development of data analysis techniques. Here we describe the experiment, its strengths and weaknesses, and the… ▽ More COMPASS is an on-axis 2.6 meter telescope coupled to a correlation polarimeter. The entire instrument was built specifically for CMB polarization studies. Careful attention was given to receiver and optics design, stability of the pointing platform, avoidance of systematic offsets, and development of data analysis techniques. Here we describe the experiment, its strengths and weaknesses, and the various things we have learned that may benefit future efforts to measure the polarization of the CMB. △ Less

Submitted 23 June, 2003; v1 submitted 30 May, 2003; originally announced May 2003.

Comments: To be published in the proceedings of "The Cosmic Microwave Background and its Polarization", New Astronomy Reviews, (eds. S. Hanany and K.A. Olive)

Gravitational Potential Reconstruction from Peculiar Velocity and Weak Lensing Measurements

Authors: Olivier Doré, Lloyd Knox, Alan Peel

Abstract: We present an analytic method for rapidly forecasting the accuracy of gravitational potential reconstruction possible from measurement of radial peculiar velocities of every galaxy cluster with M > M_th in solid angle θしーた^2 and over redshift range z_min < z < z_max. These radial velocities can be determined from measurement of the kinetic and thermal Sunyaev-Zeldovich effects. For a shallow survey… ▽ More We present an analytic method for rapidly forecasting the accuracy of gravitational potential reconstruction possible from measurement of radial peculiar velocities of every galaxy cluster with M > M_th in solid angle θしーた^2 and over redshift range z_min < z < z_max. These radial velocities can be determined from measurement of the kinetic and thermal Sunyaev-Zeldovich effects. For a shallow survey with 0.2 < z < 0.4, coincident with the SDSS photometric survey, one mode of the gravitational potential (on length scales > 60 Mpc) can be reconstructed for every ~8 cluster velocity determinations. Deeper surveys require measurement of more clusters per S/N > 1 mode. Accuracy is limited by the ``undersampling noise'' due to our non-observation of the large fraction of mass that is not in galaxy clusters. Determining the gravitational potential will allow for detailed study of the relationship between galaxies and their surrounding large-scale density fields over a wide range of redshifts, and test the gravitational instability paradigm on very large scales. Observation of weak lensing by large-scale structure provides complementary information since lensing is sensitive to the tangential modes that do not affect the velocity. △ Less

Submitted 18 July, 2002; originally announced July 2002.

Comments: 4 pages, 3 figures, submitted to ApJ Letters

Journal ref: Astrophys.J. 585 (2003) L81-L84

Using Galaxy Cluster Peculiar Velocities to Constrain Cosmological Parameters

Authors: A. Peel, L. Knox

Abstract: Galaxy cluster peculiar velocities can be inferred from high-sensitivity, high-resolution multiple-frequency observations in the 30 to 400 GHz range. While galaxy cluster counts and power spectra are sensitive to the growth factor, peculiar velocities are sensitive to the time-derivative of the growth factor and are hence complementary. Using linear perturbation theory, we forecast constraints o… ▽ More Galaxy cluster peculiar velocities can be inferred from high-sensitivity, high-resolution multiple-frequency observations in the 30 to 400 GHz range. While galaxy cluster counts and power spectra are sensitive to the growth factor, peculiar velocities are sensitive to the time-derivative of the growth factor and are hence complementary. Using linear perturbation theory, we forecast constraints on Omega(matter), the Hubble constant today, and the dark-energy equation of state parameter, given 820 densely sampled cluster locations (near redshift 1) from a Lambda-CDM N-body simulation and 820 sparsely sampled cluster locations in a broader redshift range. △ Less

Submitted 24 May, 2002; originally announced May 2002.

Comments: 4 pages, 2 figures, to appear in Proceedings of Dark Matter 2002, requires espcrc2.sty

CMB Power Spectrum Estimation via Hierarchical Decomposition

Authors: O. Doré, L. Knox, A. Peel

Abstract: We have developed a fast, accurate and generally applicable method for inferring the power spectrum and its uncertainties from maps of the cosmic microwave background (CMB) in the presence of inhomogeneous and correlated noise. For maps with 10 to 100 thousand pixels, we apply an exact power spectrum estimation algorithm to submaps of the data at various resolutions, and then combine the results… ▽ More We have developed a fast, accurate and generally applicable method for inferring the power spectrum and its uncertainties from maps of the cosmic microwave background (CMB) in the presence of inhomogeneous and correlated noise. For maps with 10 to 100 thousand pixels, we apply an exact power spectrum estimation algorithm to submaps of the data at various resolutions, and then combine the results in an optimal manner. To analyze larger maps efficiently one must resort to sub-optimal combinations in which cross-map power spectrum error correlations are only calculated approximately. We expect such approximations to work well in general, and in particular for the megapixel maps to come from the next generation of satellite missions. △ Less

Submitted 27 April, 2001; originally announced April 2001.

Comments: 10 pages, 5 figures, to be submitted to Phys. Rev. D

Journal ref: Phys.Rev.D64:083001,2001