Showing 1–50 of 107 results for author: Wright, A H

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'

KiDS-1000 and DES-Y1 combined: Cosmology from peak count statistics

Authors: Joachim Harnois-Deraps, Sven Heydenreich, Benjamin Giblin, Nicolas Martinet, Tilman Troester, Marika Asgari, Pierre Burger, Tiago Castro, Klaus Dolag, Catherine Heymans, Hendrik Hildebrandt, Benjamin Joachimi, Angus H. Wright

Abstract: We analyse the fourth data release of the Kilo Degree Survey (KiDS-1000) and extract cosmological parameter constraints based on the cosmic shear peak count statistics. Peaks are identified in aperture mass maps in which the filter is maximally sensitive to angular scales in the range 2-4arcmin, probing deep into the non-linear regime of structure formation. We interpret our results with a simulat… ▽ More We analyse the fourth data release of the Kilo Degree Survey (KiDS-1000) and extract cosmological parameter constraints based on the cosmic shear peak count statistics. Peaks are identified in aperture mass maps in which the filter is maximally sensitive to angular scales in the range 2-4arcmin, probing deep into the non-linear regime of structure formation. We interpret our results with a simulation-based inference pipeline, sampling over a broad $w$CDM prior volume and marginalising over uncertainties on shape calibration, photometric redshift distribution, intrinsic alignment and baryonic feedback. Our measurements constrain the structure growth parameter and the amplitude of the non-linear intrinsic alignment model to $Σ_8 \equiv σ_8\left[Ω_{\rm m}/0.3\right]^{0.60}=0.765^{+0.030}_{-0.030}$ and $A_{\rm IA}= 0.71^{+0.42}_{-0.42}$, respectively, in agreement with previous KiDS-1000 results based on two-point shear statistics. These results are robust against modelling of the non-linear physics, different scale cuts and selections of tomographic bins. The posterior is also consistent with that from the Dark Energy Survey Year-1 peak count analysis presented in Harnois-Déraps et al (2021), and hence we jointly analyse both surveys. We obtain $Σ_8^{\rm joint} \equiv σ_8\left[Ω_{\rm m}/0.3\right]^{0.57}=0.759^{+0.020}_{-0.017}$, in agreement with the Planck $w$CDM results. The shear-CMB tension on this parameter increases to $3.1σ$ when forcing $w=-1.0$, and to $4.1σ$ if comparing instead with $S_{8,Λ{\rm CDM}}^{\rm joint} = 0.736^{+0.016}_{-0.018}$, one of the tightest constraints to date on this quantity. (abridged) △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: 25 pages, 20 figures. comments welcome

KiDS-SBI: Simulation-Based Inference Analysis of KiDS-1000 Cosmic Shear

Authors: Maximilian von Wietersheim-Kramsta, Kiyam Lin, Nicolas Tessore, Benjamin Joachimi, Arthur Loureiro, Robert Reischke, Angus H. Wright

Abstract: We present a simulation-based inference (SBI) cosmological analysis of cosmic shear two-point statistics from the fourth weak gravitational lensing data release of the ESO Kilo-Degree Survey (KiDS-1000). KiDS-SBI efficiently performs non-Limber projection of the matter power spectrum via Levin's method, and constructs log-normal random matter fields on the curved sky for arbitrary cosmologies, inc… ▽ More We present a simulation-based inference (SBI) cosmological analysis of cosmic shear two-point statistics from the fourth weak gravitational lensing data release of the ESO Kilo-Degree Survey (KiDS-1000). KiDS-SBI efficiently performs non-Limber projection of the matter power spectrum via Levin's method, and constructs log-normal random matter fields on the curved sky for arbitrary cosmologies, including effective prescriptions for intrinsic alignments and baryonic feedback. The forward model samples realistic galaxy positions and shapes based on the observational characteristics, incorporating shear measurement and redshift calibration uncertainties, as well as angular anisotropies due to variations in depth and point-spread function. To enable direct comparison with standard inference, we limit our analysis to pseudo-angular power spectra. The SBI is based on sequential neural likelihood estimation to infer the posterior distribution of spatially-flat $Λ$CDM cosmological parameters from 18,000 realisations. We infer a mean marginal of the growth of structure parameter $S_{8} \equiv σ_8 (Ω_\mathrm{m} / 0.3)^{0.5} = 0.731\pm 0.033$ ($68 \%$). We present a measure of goodness-of-fit for SBI and determine that the forward model fits the data well with a probability-to-exceed of $0.42$. For fixed cosmology, the learnt likelihood is approximately Gaussian, while constraints widen compared to a Gaussian likelihood analysis due to cosmology dependence in the covariance. Neglecting variable depth and anisotropies in the point spread function in the model can cause $S_{8}$ to be overestimated by ${\sim}5\%$. Our results are in agreement with previous analysis of KiDS-1000 and reinforce a $2.9 σ$ tension with constraints from cosmic microwave background measurements. This work highlights the importance of forward-modelling systematic effects in upcoming galaxy surveys. △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: 44 pages, 30 figures. Submitted to Astronomy & Astrophysics

Improved Weak Lensing Photometric Redshift Calibration via StratLearn and Hierarchical Modeling

Authors: Maximilian Autenrieth, Angus H. Wright, Roberto Trotta, David A. van Dyk, David C. Stenning, Benjamin Joachimi

Abstract: Discrepancies between cosmological parameter estimates from cosmic shear surveys and from recent Planck cosmic microwave background measurements challenge the ability of the highly successful $Λ$CDM model to describe the nature of the Universe. To rule out systematic biases in cosmic shear survey analyses, accurate redshift calibration within tomographic bins is key. In this paper, we improve phot… ▽ More Discrepancies between cosmological parameter estimates from cosmic shear surveys and from recent Planck cosmic microwave background measurements challenge the ability of the highly successful $Λ$CDM model to describe the nature of the Universe. To rule out systematic biases in cosmic shear survey analyses, accurate redshift calibration within tomographic bins is key. In this paper, we improve photo-$z$ calibration via Bayesian hierarchical modeling of full galaxy photo-$z$ conditional densities, by employing $\textit{StratLearn}$, a recently developed statistical methodology, which accounts for systematic differences in the distribution of the spectroscopic training/source set and the photometric target set. Using realistic simulations that were designed to resemble the KiDS+VIKING-450 dataset, we show that $\textit{StratLearn}$-estimated conditional densities improve the galaxy tomographic bin assignment, and that our $\textit{StratLearn}$-Bayesian framework leads to nearly unbiased estimates of the target population means. This leads to a factor of $\sim 2$ improvement upon the previously best photo-$z$ calibration method. Our approach delivers a maximum bias per tomographic bin of $Δ\langle z \rangle = 0.0095 \pm 0.0089$, with an average absolute bias of $0.0052 \pm 0.0067$ across the five tomographic bins. △ Less

Submitted 12 March, 2024; v1 submitted 9 January, 2024; originally announced January 2024.

Comments: 24 pages, 20 figures, 3 appendices

The PAU Survey: Photometric redshift estimation in deep wide fields

Authors: D. Navarro-Gironés, E. Gaztañaga, M. Crocce, A. Wittje, H. Hildebrandt, A. H. Wright, M. Siudek, M. Eriksen, S. Serrano, P. Renard, E. J. Gonzalez, C. M. Baugh, L. Cabayol, J. Carretero, R. Casas, F. J. Castander, J. De Vicente, E. Fernandez, J. García-Bellido, H. Hoekstra, G. Manzoni, R. Miquel, C. Padilla, E. Sánchez, I. Sevilla-Noarbe , et al. (1 additional authors not shown)

Abstract: We present photometric redshifts (photo-$z$) for the deep wide fields of the Physics of the Accelerating Universe Survey (PAUS), covering an area of $\sim$50 deg$^{2}$, for $\sim$1.8 million objects up to $i_{\textrm{AB}}<23$. The PAUS deep wide fields overlap with the W1 and W3 fields from CFHTLenS and the G09 field from KiDS/GAMA. Photo-$z$ are estimated using the 40 narrow bands (NB) of PAUS an… ▽ More We present photometric redshifts (photo-$z$) for the deep wide fields of the Physics of the Accelerating Universe Survey (PAUS), covering an area of $\sim$50 deg$^{2}$, for $\sim$1.8 million objects up to $i_{\textrm{AB}}<23$. The PAUS deep wide fields overlap with the W1 and W3 fields from CFHTLenS and the G09 field from KiDS/GAMA. Photo-$z$ are estimated using the 40 narrow bands (NB) of PAUS and the broad bands (BB) of CFHTLenS and KiDS. We compute the redshifts with the SED template-fitting code BCNZ, with a modification in the calibration technique of the zero-point between the observed and the modelled fluxes, that removes any dependence on spectroscopic redshift samples. We enhance the redshift accuracy by introducing an additional photo-$z$ estimate ($z_{\textrm{b}}$), obtained through the combination of the BCNZ and the BB-only photo-$z$. Comparing with spectroscopic redshifts estimates ($z_{\textrm{s}}$), we obtain a $σ_{68} \simeq 0.019$ for all galaxies with $i_{\textrm{AB}}<23$ and a typical bias $|z_{\textrm{b}}-z_{\textrm{s}}|$ smaller than 0.01. For $z_{\textrm{b}} \sim (0.10-0.75)$ we find $σ_{68} \simeq (0.003-0.02)$, this is a factor of $10-2$ higher accuracy than the corresponding BB-only results. We obtain similar performance when we split the samples into red (passive) and blue (active) galaxies. We validate the redshift probability $p(z)$ obtained by BCNZ and compare its performance with that of $z_{\textrm{b}}$. These photo-$z$ catalogues will facilitate important science cases, such as the study of galaxy clustering and intrinsic alignment at high redshifts ($z \lesssim 1$) and faint magnitudes. △ Less

Submitted 10 December, 2023; originally announced December 2023.

Comments: 24 pages, 26 figures, submitted to MNRAS

KiDS-1000: Cosmology with improved cosmic shear measurements

Authors: Shun-Sheng Li, Henk Hoekstra, Konrad Kuijken, Marika Asgari, Maciej Bilicki, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Benjamin Joachimi, Lance Miller, Jan Luca van den Busch, Angus H. Wright, Arun Kannawadi, Robert Reischke, HuanYuan Shan

Abstract: We present refined cosmological parameter constraints derived from a cosmic shear analysis of the fourth data release of the Kilo-Degree Survey (KiDS-1000). Our main improvements include enhanced galaxy shape measurements made possible by an updated version of the lensfit code and improved shear calibration achieved with a newly developed suite of multi-band image simulations. Additionally, we inc… ▽ More We present refined cosmological parameter constraints derived from a cosmic shear analysis of the fourth data release of the Kilo-Degree Survey (KiDS-1000). Our main improvements include enhanced galaxy shape measurements made possible by an updated version of the lensfit code and improved shear calibration achieved with a newly developed suite of multi-band image simulations. Additionally, we incorporated recent advancements in cosmological inference from the joint Dark Energy Survey Year 3 and KiDS-1000 cosmic shear analysis. Assuming a spatially flat standard cosmological model, we constrain $S_8\equivσ_8(Ω_{\rm m}/0.3)^{0.5} = 0.776_{-0.027-0.003}^{+0.029+0.002}$, where the second set of uncertainties accounts for the systematic uncertainties within the shear calibration. These systematic uncertainties stem from minor deviations from realism in the image simulations and the sensitivity of the shear measurement algorithm to the morphology of the galaxy sample. Despite these changes, our results align with previous KiDS studies and other weak lensing surveys, and we find a ${\sim}2.3σ$ level of tension with the Planck cosmic microwave background constraints on $S_8$. △ Less

Submitted 2 November, 2023; v1 submitted 19 June, 2023; originally announced June 2023.

Comments: 20 pages, 13 figures, 4 tables, minor revisions to match the final accepted version

Journal ref: A&A 679, A133 (2023)

DES Y3 + KiDS-1000: Consistent cosmology combining cosmic shear surveys

Authors: Dark Energy Survey, Kilo-Degree Survey Collaboration, :, T. M. C. Abbott, M. Aguena, A. Alarcon, O. Alves, A. Amon, F. Andrade-Oliveira, M. Asgari, S. Avila, D. Bacon, K. Bechtol, M. R. Becker, G. M. Bernstein, E. Bertin, M. Bilicki, J. Blazek, S. Bocquet, D. Brooks, P. Burger, D. L. Burke, H. Camacho, A. Campos, A. Carnero Rosell , et al. (138 additional authors not shown)

Abstract: We present a joint cosmic shear analysis of the Dark Energy Survey (DES Y3) and the Kilo-Degree Survey (KiDS-1000) in a collaborative effort between the two survey teams. We find consistent cosmological parameter constraints between DES Y3 and KiDS-1000 which, when combined in a joint-survey analysis, constrain the parameter $S_8 = σ_8 \sqrt{Ω_{\rm m}/0.3}$ with a mean value of… ▽ More We present a joint cosmic shear analysis of the Dark Energy Survey (DES Y3) and the Kilo-Degree Survey (KiDS-1000) in a collaborative effort between the two survey teams. We find consistent cosmological parameter constraints between DES Y3 and KiDS-1000 which, when combined in a joint-survey analysis, constrain the parameter $S_8 = σ_8 \sqrt{Ω_{\rm m}/0.3}$ with a mean value of $0.790^{+0.018}_{-0.014}$. The mean marginal is lower than the maximum a posteriori estimate, $S_8=0.801$, owing to skewness in the marginal distribution and projection effects in the multi-dimensional parameter space. Our results are consistent with $S_8$ constraints from observations of the cosmic microwave background by Planck, with agreement at the $1.7σ$ level. We use a Hybrid analysis pipeline, defined from a mock survey study quantifying the impact of the different analysis choices originally adopted by each survey team. We review intrinsic alignment models, baryon feedback mitigation strategies, priors, samplers and models of the non-linear matter power spectrum. △ Less

Submitted 19 October, 2023; v1 submitted 26 May, 2023; originally announced May 2023.

Comments: 40 pages, 21 figures, 15 tables, accepted Open Journal of Astrophysics. Download the chains from https://des.ncsa.illinois.edu/releases/y3a2/Y3key-joint-des-kids or create your own chains with CosmoSIS using https://github.com/joezuntz/cosmosis-standard-library/blob/main/examples/des-y3_and_kids-1000.ini Watch the core team discuss this analysis at https://cosmologytalks.com/2023/05/26/des-kids

Report number: FERMILAB-PUB-23-267-PPD

ACT-DR5 Sunyaev-Zel'dovich Clusters: weak lensing mass calibration with KiDS

Authors: Naomi Clare Robertson, Cristóbal Sifón, Marika Asgari, Nicholas Battaglia, Maciej Bilicki, J. Richard Bond, Mark J. Devlin, Jo Dunkley, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Matt Hilton, Henk Hoekstra, John P. Hughes, Konrad Kuijken, Thibaut Louis, Maya Mallaby-Kay, Lyman Page, Bruce Partridge, Mario Radovich, Peter Schneider, HuanYuan Shan, David N. Spergel, Tilman Tröster, Edward J. Wollack , et al. (2 additional authors not shown)

Abstract: We present weak gravitational lensing measurements of a sample of 157 clusters within the Kilo Degree Survey (KiDS), detected with a $>5σ$ thermal Sunyaev-Zel'dovich (SZ) signal by the Atacama Cosmology Telescope (ACT). Using a halo-model approach we constrain the average total cluster mass, $M_{\rm WL}$, accounting for the ACT cluster selection function of the full sample. We find that the SZ clu… ▽ More We present weak gravitational lensing measurements of a sample of 157 clusters within the Kilo Degree Survey (KiDS), detected with a $>5σ$ thermal Sunyaev-Zel'dovich (SZ) signal by the Atacama Cosmology Telescope (ACT). Using a halo-model approach we constrain the average total cluster mass, $M_{\rm WL}$, accounting for the ACT cluster selection function of the full sample. We find that the SZ cluster mass estimate $M_{\rm SZ}$, which was calibrated using X-ray observations, is biased with $M_{\rm SZ}/M_{\rm WL} = (1-b_{\rm SZ}) = 0.65\pm 0.05$. Separating the sample into six mass bins, we find no evidence of a strong mass-dependency for the mass bias, $(1-b_{\rm SZ})$. Adopting this ACT-KiDS SZ mass-calibration would bring the Planck SZ cluster count into agreement with the counts expected from the {\it Planck} cosmic microwave background $Λ$CDM cosmological model, although it should be noted that the cluster sample considered in this work has a lower average mass $M_{\rm SZ, uncor} = 3.64 \times 10^{14} M_{\odot}$ compared to the Planck cluster sample which has an average mass in the range $M_{\rm SZ, uncor} = (5.5-8.5) \times 10^{14} M_{\odot}$, depending on the sub-sample used. △ Less

Submitted 20 April, 2023; originally announced April 2023.

Comments: 12 pages, 7 figures

Exploring the mass and redshift dependence of the cluster pressure profile with stacks on thermal SZ maps

Authors: Denis Tramonte, Yin-Zhe Ma, Ziang Yan, Matteo Maturi, Gianluca Castignani, Mauro Sereno, Sandro Bardelli, Carlo Giocoli, Federico Marulli, Lauro Moscardini, Emanuella Puddu, Mario Radovich, Ludovic Van Waerbeke, Angus H. Wright

Abstract: We provide novel constraints on the parameters defining the universal pressure profile (UPP) within clusters of galaxies, and explore their dependence on the cluster mass and redshift, from measurements of Sunyaev-Zel'dovich Compton-$y$ profiles. We employ both the $\textit{Planck}$ 2015 MILCA and the ACT-DR4 $y$ maps over the common $\sim 2,100\,\text{deg}^2$ footprint. We combine existing cluste… ▽ More We provide novel constraints on the parameters defining the universal pressure profile (UPP) within clusters of galaxies, and explore their dependence on the cluster mass and redshift, from measurements of Sunyaev-Zel'dovich Compton-$y$ profiles. We employ both the $\textit{Planck}$ 2015 MILCA and the ACT-DR4 $y$ maps over the common $\sim 2,100\,\text{deg}^2$ footprint. We combine existing cluster catalogs based on KiDS, SDSS and DESI observations, for a total of 23,820 clusters spanning the mass range $10^{14.0}\,\text{M}_{\odot}<M_{500}<10^{15.1}\,\text{M}_{\odot}$ and the redshift range $0.02<z<0.98$. We split the clusters into three independent bins in mass and redshift; for each combination we detect the stacked SZ cluster signal and extract the mean $y$ angular profile. The latter is predicted theoretically adopting a halo model framework, and MCMCs are employed to estimate the UPP parameters, the hydrostatic mass bias $b_{\rm h}$ and possible cluster miscentering effects. We constrain $[P_0,c_{500},α,β]$ to $[5.9,2.0,1.8,4.9]$ with $\textit{Planck}$ and to $[3.8,1.3,1.0,4.4]$ with ACT using the full cluster sample, in agreement with previous findings. We do not find any compelling evidence for a residual mass or redshift dependence, thus expanding the validity of the cluster pressure profile over much larger $M_{500}$ and $z$ ranges; this is the first time the model has been tested on such a large (complete and representative) cluster sample. Finally, we obtain loose constraints on the hydrostatic mass bias in the range 0.2-0.3, again in broad agreement with previous works. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Comments: 39 pages, 22 figures. Accepted for publication in ApjS

Modeling Strong Lenses from Wide-Field Ground-Based Observations in KiDS and GAMA

Authors: Shawn Knabel, B. W. Holwerda, J. Nightingale, T. Treu, M. Bilicki, S. Brough, S. Driver, L. Finnerty, L. Haberzettl, S. Hegde, A. M. Hopkins, K. Kuijken, J. Liske, K. A. Pimbblet, R. C. Steele, A. H. Wright

Abstract: Despite the success of galaxy-scale strong gravitational lens studies with Hubble-quality imaging, the number of well-studied strong lenses remains small. As a result, robust comparisons of the lens models to theoretical predictions are difficult. This motivates our application of automated Bayesian lens modeling methods to observations from public data releases of overlapping large ground-based i… ▽ More Despite the success of galaxy-scale strong gravitational lens studies with Hubble-quality imaging, the number of well-studied strong lenses remains small. As a result, robust comparisons of the lens models to theoretical predictions are difficult. This motivates our application of automated Bayesian lens modeling methods to observations from public data releases of overlapping large ground-based imaging and spectroscopic surveys: Kilo-Degree Survey (KiDS) and Galaxy and Mass Assembly (GAMA), respectively. We use the open-source lens modeling software PyAutoLens to perform our analysis. We demonstrate the feasibility of strong lens modeling with large-survey data at lower resolution as a complementary avenue to studies that utilize more time-consuming and expensive observations of individual lenses at higher resolution. We discuss advantages and challenges, with special consideration given to determining background source redshifts from single-aperture spectra and to disentangling foreground lens and background source light. High uncertainties in the best-fit parameters for the models due to the limits of optical resolution in ground-based observatories and the small sample size can be improved with future study. We give broadly applicable recommendations for future efforts, and with proper application this approach could yield measurements in the quantities needed for robust statistical inference. △ Less

Submitted 17 January, 2023; v1 submitted 12 January, 2023; originally announced January 2023.

Comments: 25 pages, 19 figures, accepted for publication in MNRAS following peer review

Potential scientific synergies in weak lensing studies between the CSST and Euclid space probes

Authors: D. Z. Liu, X. M. Meng, X. Z. Er, Z. H. Fan, M. Kilbinger, G. L. Li, R. Li, T. Schrabback, D. Scognamiglio, H. Y. Shan, C. Tao, Y. S. Ting, J. Zhang, S. H. Cheng, S. Farrens, L. P. Fu, H. Hildebrandt, X. Kang, J. P. Kneib, X. K. Liu, Y. Mellier, R. Nakajima, P. Schneider, J. L. Starck, C. L. Wei , et al. (2 additional authors not shown)

Abstract: Aims. With the next generation of large surveys coming to the stage of observational cosmology soon, it is important to explore their potential synergies and to maximise their scientific outcomes. In this study, we aim to investigate the complementarity of the two upcoming space missions Euclid and the China Space Station Telescope (CSST), focusing on weak lensing (WL) cosmology. In particular, we… ▽ More Aims. With the next generation of large surveys coming to the stage of observational cosmology soon, it is important to explore their potential synergies and to maximise their scientific outcomes. In this study, we aim to investigate the complementarity of the two upcoming space missions Euclid and the China Space Station Telescope (CSST), focusing on weak lensing (WL) cosmology. In particular, we analyse the photometric redshifts (photo-zs) and the galaxy blending effects. For Euclid, WL measurements suffer from chromatic PSF effects. For this, CSST can provide valuable information for Euclid to obtain more accurate PSF, and to calibrate the color and color-gradient biases for WL measurements. Methods. We create image simulations for different surveys, and quantify the photo-z performance. For blending analyses, we employ high-resolution HST/CANDELS data to mock Euclid, CSST, and an LSST-like survey. We analyse the blending fraction for different cases, and the blending effects on galaxy photometry. Furthermore, we demonstrate that CSST can provide a large enough number of high SNR multi-band galaxy images to calibrate the color-gradient biases for Euclid. Results. The sky coverage of Euclid lies entirely within the CSST footprint. The combination of Euclid with CSST data can be done more uniformly than with the various ground-based data. Our studies show that by combining Euclid and CSST, we can reach a photo-z precision of $σ_{\rm NMAD} \approx 0.04$, and an outlier fraction of $η\approx 2.4\%$. Because of the similarly high resolutions, the data combination of Euclid and CSST can be relatively straightforward for photometry. To include ground-based data, however, sophisticated deblending utilizing priors from high-resolution space data is demanded. The color-gradient biases for Euclid can be well calibrated to the level of 0.1% using galaxies from CSST deep survey. △ Less

Submitted 28 October, 2022; originally announced October 2022.

Comments: 18 pages, 19 figures and 2 tables. Accepted for publication in A&A

Journal ref: A&A 669, A128 (2023)

KiDS-Legacy calibration: unifying shear and redshift calibration with the SKiLLS multi-band image simulations

Authors: Shun-Sheng Li, Konrad Kuijken, Henk Hoekstra, Lance Miller, Catherine Heymans, Hendrik Hildebrandt, Jan Luca van den Busch, Angus H. Wright, Mijin Yoon, Maciej Bilicki, Matías Bravo, Claudia del P. Lagos

Abstract: We present SKiLLS, a suite of multi-band image simulations for the weak lensing analysis of the complete Kilo-Degree Survey (KiDS), dubbed KiDS-Legacy analysis. The resulting catalogues enable joint shear and redshift calibration, enhancing the realism and hence accuracy over previous efforts. To create a large volume of simulated galaxies with faithful properties and to a sufficient depth, we int… ▽ More We present SKiLLS, a suite of multi-band image simulations for the weak lensing analysis of the complete Kilo-Degree Survey (KiDS), dubbed KiDS-Legacy analysis. The resulting catalogues enable joint shear and redshift calibration, enhancing the realism and hence accuracy over previous efforts. To create a large volume of simulated galaxies with faithful properties and to a sufficient depth, we integrated cosmological simulations with high-quality imaging observations. We also improved the realism of simulated images by allowing the point spread function (PSF) to differ between CCD images, including stellar density variations and varying noise levels between pointings. Using realistic variable shear fields, we accounted for the impact of blended systems at different redshifts. Although the overall correction is minor, we found a clear redshift-bias correlation in the blending-only variable shear simulations, indicating the non-trivial impact of this higher-order blending effect. We also explored the impact of the PSF modelling errors and found a small yet noticeable effect on the shear bias. Finally, we conducted a series of sensitivity tests, including changing the input galaxy properties. We conclude that our fiducial shape measurement algorithm, lensfit, is robust within the requirements of lensing analyses with KiDS. As for future weak lensing surveys with tighter requirements, we suggest further investments in understanding the impact of blends at different redshifts, improving the PSF modelling algorithm and developing the shape measurement method to be less sensitive to the galaxy properties. △ Less

Submitted 22 January, 2023; v1 submitted 13 October, 2022; originally announced October 2022.

Comments: 28 pages, 31 figures, 2 tables, minor revisions to match the final accepted version

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

KiDS-1000: Combined halo-model cosmology constraints from galaxy abundance, galaxy clustering and galaxy-galaxy lensing

Authors: Andrej Dvornik, Catherine Heymans, Marika Asgari, Constance Mahony, Benjamin Joachimi, Maciej Bilicki, Elisa Chisari, Hendrik Hildebrandt, Henk Hoekstra, Harry Johnston, Konrad Kuijken, Alexander Mead, Hironao Miyatake, Takahiro Nishimichi, Robert Reischke, Sandra Unruh, Angus H. Wright

Abstract: We present constraints on the flat $Λ$CDM cosmological model through a joint analysis of galaxy abundance, galaxy clustering and galaxy-galaxy lensing observables with the Kilo-Degree Survey. Our theoretical model combines a flexible conditional stellar mass function, to describe the galaxy-halo connection, with a cosmological N-body simulation-calibrated halo model to describe the non-linear matt… ▽ More We present constraints on the flat $Λ$CDM cosmological model through a joint analysis of galaxy abundance, galaxy clustering and galaxy-galaxy lensing observables with the Kilo-Degree Survey. Our theoretical model combines a flexible conditional stellar mass function, to describe the galaxy-halo connection, with a cosmological N-body simulation-calibrated halo model to describe the non-linear matter field. Our magnitude-limited bright galaxy sample combines 9-band optical-to-near-infrared photometry with an extensive and complete spectroscopic training sample to provide accurate redshift and stellar mass estimates. Our faint galaxy sample provides a background of accurately calibrated lensing measurements. We constrain the structure growth parameter $S_8=σ_8\sqrt{Ω_{\mathrm{m}}/0.3}=0.773^{+0.028}_{-0.030}$, and the matter density parameter $Ω_{\mathrm{m}}=0.290^{+0.021}_{-0.017}$. The galaxy-halo connection model adopted in the work is shown to be in agreement with previous studies. Our constraints on cosmological parameters are comparable to, and consistent with, joint $3\times2{\mathrm{pt}}$ clustering-lensing analyses that additionally include a cosmic shear observable. This analysis therefore brings attention to the significant constraining power in the often-excluded non-linear scales for galaxy clustering and galaxy-galaxy lensing observables. By adopting a theoretical model that accounts for non-linear halo bias, halo exclusion, scale-dependent galaxy bias and the impact of baryon feedback, this work demonstrates the potential and a way forward to include non-linear scales in cosmological analyses. Varying the width of the satellite galaxy distribution with an additional parameter yields a strong preference for sub-Poissonian variance, improving the goodness of fit by 0.18 in reduced $χ^{2}$ value compared to a fixed Poisson distribution. △ Less

Submitted 13 June, 2023; v1 submitted 6 October, 2022; originally announced October 2022.

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

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

Evolving Complexity is Hard

Authors: Alden H. Wright, Cheyenne L. Laue

Abstract: Understanding the evolution of complexity is an important topic in a wide variety of academic fields. Implications of better understanding complexity include increased knowledge of major evolutionary transitions and the properties of living and technological systems. Genotype-phenotype (G-P) maps are fundamental to evolution, and biologically-oriented G-P maps have been shown to have interesting a… ▽ More Understanding the evolution of complexity is an important topic in a wide variety of academic fields. Implications of better understanding complexity include increased knowledge of major evolutionary transitions and the properties of living and technological systems. Genotype-phenotype (G-P) maps are fundamental to evolution, and biologically-oriented G-P maps have been shown to have interesting and often-universal properties that enable evolution by following phenotype-preserving walks in genotype space. Here we use a digital logic gate circuit G-P map where genotypes are represented by circuits and phenotypes by the functions that the circuits compute. We compare two mathematical definitions of circuit and phenotype complexity and show how these definitions relate to other well-known properties of evolution such as redundancy, robustness, and evolvability. Using both Cartesian and Linear genetic programming implementations, we demonstrate that the logic gate circuit shares many universal properties of biologically derived G-P maps, with the exception of the relationship between one method of computing phenotypic evolvability, robustness, and complexity. Due to the inherent structure of the G-P map, including the predominance of rare phenotypes, large interconnected neutral networks, and the high mutational load of low robustness, complex phenotypes are difficult to discover using evolution. We suggest, based on this evidence, that evolving complexity is hard and we discuss computational strategies for genetic-programming-based evolution to successfully find genotypes that map to complex phenotypes in the search space. △ Less

Submitted 16 September, 2022; originally announced September 2022.

Comments: 15 pages, 13 figures. To be published in the proceedings of the GPTP XIX Workshop

Euclid: Calibrating photometric redshifts with spectroscopic cross-correlations

Authors: K. Naidoo, H. Johnston, B. Joachimi, J. L. van den Busch, H. Hildebrandt, O. Ilbert, O. Lahav, N. Aghanim, B. Altieri, A. Amara, M. Baldi, R. Bender, C. Bodendorf, E. Branchini, M. Brescia, J. Brinchmann, S. Camera, V. Capobianco, C. Carbone, J. Carretero, F. J. Castander, M. Castellano, S. Cavuoti, A. Cimatti, R. Cledassou , et al. (104 additional authors not shown)

Abstract: Cosmological constraints from key probes of the Euclid imaging survey rely critically on the accurate determination of the true redshift distributions, $n(z)$, of tomographic redshift bins. We determine whether the mean redshift, $<z>$, of ten Euclid tomographic redshift bins can be calibrated to the Euclid target uncertainties of $σ(<z>)<0.002\,(1+z)$ via cross-correlation, with spectroscopic sam… ▽ More Cosmological constraints from key probes of the Euclid imaging survey rely critically on the accurate determination of the true redshift distributions, $n(z)$, of tomographic redshift bins. We determine whether the mean redshift, $<z>$, of ten Euclid tomographic redshift bins can be calibrated to the Euclid target uncertainties of $σ(<z>)<0.002\,(1+z)$ via cross-correlation, with spectroscopic samples akin to those from the Baryon Oscillation Spectroscopic Survey (BOSS), Dark Energy Spectroscopic Instrument (DESI), and Euclid's NISP spectroscopic survey. We construct mock Euclid and spectroscopic galaxy samples from the Flagship simulation and measure small-scale clustering redshifts up to redshift $z<1.8$ with an algorithm that performs well on current galaxy survey data. The clustering measurements are then fitted to two $n(z)$ models: one is the true $n(z)$ with a free mean; the other a Gaussian Process modified to be restricted to non-negative values. We show that $<z>$ is measured in each tomographic redshift bin to an accuracy of order 0.01 or better. By measuring the clustering redshifts on subsets of the full Flagship area, we construct scaling relations that allow us to extrapolate the method performance to larger sky areas than are currently available in the mock. For the full expected Euclid, BOSS, and DESI overlap region of approximately 6000 deg$^{2}$, the uncertainties attainable by clustering redshifts exceeds the Euclid requirement by at least a factor of three for both $n(z)$ models considered, although systematic biases limit the accuracy. Clustering redshifts are an extremely effective method for redshift calibration for Euclid if the sources of systematic biases can be determined and removed, or calibrated-out with sufficiently realistic simulations. We outline possible future work, in particular an extension to higher redshifts with quasar reference samples. △ Less

Submitted 24 February, 2023; v1 submitted 22 August, 2022; originally announced August 2022.

Comments: 14 pages, 8 figures, accepted for publication in A&A

Journal ref: A&A 670, A149 (2023)

Extending empirical constraints on the SZ-mass scaling relation to higher redshifts via HST weak lensing measurements of nine clusters from the SPT-SZ survey at $z\gtrsim1$

Authors: Hannah Zohren, Tim Schrabback, Sebastian Bocquet, Martin Sommer, Fatimah Raihan, Beatriz Hernández-Martín, Ole Marggraf, Behzad Ansarinejad, Matthew B. Bayliss, Lindsey E. Bleem, Thomas Erben, Henk Hoekstra, Benjamin Floyd, Michael D. Gladders, Florian Kleinebreil, Michael A. McDonald, Mischa Schirmer, Diana Scognamiglio, Keren Sharon, Angus H. Wright

Abstract: We present a Hubble Space Telescope (HST) weak gravitational lensing study of nine distant and massive galaxy clusters with redshifts $1.0 \lesssim z \lesssim 1.7$ ($z_\mathrm{median} = 1.4$) and Sunyaev Zel'dovich (SZ) detection significance $ξ> 6.0$ from the South Pole Telescope Sunyaev Zel'dovich (SPT-SZ) survey. We measured weak lensing galaxy shapes in HST/ACS F606W and F814W images and used… ▽ More We present a Hubble Space Telescope (HST) weak gravitational lensing study of nine distant and massive galaxy clusters with redshifts $1.0 \lesssim z \lesssim 1.7$ ($z_\mathrm{median} = 1.4$) and Sunyaev Zel'dovich (SZ) detection significance $ξ> 6.0$ from the South Pole Telescope Sunyaev Zel'dovich (SPT-SZ) survey. We measured weak lensing galaxy shapes in HST/ACS F606W and F814W images and used additional observations from HST/WFC3 in F110W and VLT/FORS2 in $U_\mathrm{HIGH}$ to preferentially select background galaxies at $z\gtrsim 1.8$, achieving a high purity. We combined recent redshift estimates from the CANDELS/3D-HST and HUDF fields to infer an improved estimate of the source redshift distribution. We measured weak lensing masses by fitting the tangential reduced shear profiles with spherical Navarro-Frenk-White (NFW) models. We obtained the largest lensing mass in our sample for the cluster SPT-CLJ2040$-$4451, thereby confirming earlier results that suggest a high lensing mass of this cluster compared to X-ray and SZ mass measurements. Combining our weak lensing mass constraints with results obtained by previous studies for lower redshift clusters, we extended the calibration of the scaling relation between the unbiased SZ detection significance $ζ$ and the cluster mass for the SPT-SZ survey out to higher redshifts. We found that the mass scale inferred from our highest redshift bin ($1.2 < z < 1.7$) is consistent with an extrapolation of constraints derived from lower redshifts, albeit with large statistical uncertainties. Thus, our results show a similar tendency as found in previous studies, where the cluster mass scale derived from the weak lensing data is lower than the mass scale expected in a Planck $νΛ$CDM (i.e. $ν$ $Λ$ Cold Dark Matter) cosmology given the SPT-SZ cluster number counts. △ Less

Submitted 22 August, 2022; originally announced August 2022.

Comments: 23 pages, 10 figures + appendix (9 pages, 9 figures); accepted for publication in A&A

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

KiDS-1000 Cosmology: Constraints from density split statistics

Authors: Pierre A. Burger, Oliver Friedrich, Joachim Harnois-Déraps, Peter Schneider, Marika Asgari, Maciej Bilicki, Hendrik Hildebrandt, Angus H. Wright, Tiago Castro, Klaus Dolag, Catherine Heymans, Benjamin Joachimi, Konrad Kuijken, Nicolas Martinet, HuanYuan Shan, Tilman Tröster

Abstract: Context. Weak lensing and clustering statistics beyond two-point functions can capture non-Gaussian information about the matter density field, thereby improving the constraints on cosmological parameters relative to the mainstream methods based on correlation functions and power spectra. Aims. This paper presents a cosmological analysis of the fourth data release of the Kilo Degree Survey based o… ▽ More Context. Weak lensing and clustering statistics beyond two-point functions can capture non-Gaussian information about the matter density field, thereby improving the constraints on cosmological parameters relative to the mainstream methods based on correlation functions and power spectra. Aims. This paper presents a cosmological analysis of the fourth data release of the Kilo Degree Survey based on the density split statistics, which measures the mean shear profiles around regions classified according to foreground densities. The latter is constructed from a bright galaxy sample, which we further split into red and blue samples, allowing us to probe their respective connection to the underlying dark matter density. Methods. We use the state-of-the-art model of the density splitting statistics and validate its robustness against mock data infused with known systematic effects such as intrinsic galaxy alignment and baryonic feedback. Results. After marginalising over the photometric redshift uncertainty and the residual shear calibration bias, we measure for the full KiDS-bright sample a structure growth parameter of $S_8 = σ_8 \sqrt{Ω_\mathrm{m}/0.3} = 0.74^{+0.03}_{-0.02}$ that is competitive to and consistent with two-point cosmic shear results, a matter density of $Ω_\mathrm{m} = 0.28 \pm 0.02$, and a constant galaxy bias of $b = 1.32^{+0.12}_{-0.10}$. △ Less

Submitted 10 January, 2023; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: 18 pages, 23 Figures

Report number: A69

Journal ref: A&A 669, A69 (2023)

The intrinsic reddening of the Magellanic Clouds as traced by background galaxies -- III. The Large Magellanic Cloud

Authors: Cameron P. M. Bell, Maria-Rosa L. Cioni, Angus H. Wright, David L. Nidever, I-Da Chiang, Samyaday Choudhury, Martin A. T. Groenewegen, Clara M. Pennock, Yumi Choi, Richard de Grijs, Valentin D. Ivanov, Pol Massana, Ambra Nanni, Noelia E. D. Noël, Knut Olsen, Jacco Th. van Loon, A. Katherina Vivas, Dennis Zaritsky

Abstract: We present a map of the total intrinsic reddening across ~90 deg$^{2}$ of the Large Magellanic Cloud (LMC) derived using optical (ugriz) and near-infrared (IR; YJKs) spectral energy distributions (SEDs) of background galaxies. The reddening map is created from a sample of 222,752 early-type galaxies based on the LEPHARE $χ^{2}$ minimisation SED-fitting routine. We find excellent agreement between… ▽ More We present a map of the total intrinsic reddening across ~90 deg$^{2}$ of the Large Magellanic Cloud (LMC) derived using optical (ugriz) and near-infrared (IR; YJKs) spectral energy distributions (SEDs) of background galaxies. The reddening map is created from a sample of 222,752 early-type galaxies based on the LEPHARE $χ^{2}$ minimisation SED-fitting routine. We find excellent agreement between the regions of enhanced intrinsic reddening across the central (4x4 deg$^2$) region of the LMC and the morphology of the low-level pervasive dust emission as traced by far-IR emission. In addition, we are able to distinguish smaller, isolated enhancements that are coincident with known star-forming regions and the clustering of young stars observed in morphology maps. The level of reddening associated with the molecular ridge south of 30 Doradus is, however, smaller than in the literature reddening maps. The reduced number of galaxies detected in this region, due to high extinction and crowding, may bias our results towards lower reddening values. Our map is consistent with maps derived from red clump stars and from the analysis of the star formation history across the LMC. This study represents one of the first large-scale categorisations of extragalactic sources behind the LMC and as such we provide the LEPHARE outputs for our full sample of ~2.5 million sources. △ Less

Submitted 31 May, 2022; v1 submitted 9 May, 2022; originally announced May 2022.

Comments: 18 pages, 16 figures, accepted for publication in MNRAS

KiDS+VIKING+GAMA: Halo occupation distributions and correlations of satellite numbers with a new halo model of the galaxy-matter bispectrum for galaxy-galaxy-galaxy lensing

Authors: Laila Linke, Patrick Simon, Peter Schneider, Daniel J. Farrow, Jens Rödiger, Angus H. Wright

Abstract: Halo models and halo occupation distributions (HODs) are important tools to model the galaxy and matter distribution. We present and assess a new method for constraining the parameters of HODs using the gravitational lensing shear around galaxy pairs, galaxy-galaxy-galaxy-lensing (G3L). In contrast to galaxy-galaxy-lensing, G3L is sensitive to correlations between the per-halo numbers of galaxies… ▽ More Halo models and halo occupation distributions (HODs) are important tools to model the galaxy and matter distribution. We present and assess a new method for constraining the parameters of HODs using the gravitational lensing shear around galaxy pairs, galaxy-galaxy-galaxy-lensing (G3L). In contrast to galaxy-galaxy-lensing, G3L is sensitive to correlations between the per-halo numbers of galaxies from different populations. We use G3L to probe these correlations and test the default hypothesis that they are negligible. We derive a halo model for G3L and validate it with realistic mock data from the Millennium Simulation and a semi-analytic galaxy model. Then, we analyse public data from the Kilo-Degree Survey (KiDS), the VISTA Infrared Kilo-Degree Galaxy Survey (VIKING) and data from the Galaxy And Mass Assembly Survey (GAMA) to infer the HODs of galaxies at $z<0.5$ in five different stellar mass bins between $10^{8.5}h^{-2} M_\odot$ and $10^{11.5}h^{-2} M_\odot$ and two colours (red and blue), as well as correlations between satellite numbers. The analysis recovers the true HODs in the simulated data within the $68\%$ credibility range. The inferred HODs vary significantly with colour and stellar mass. There is also strong evidence ($>3σ$) for correlations, increasing with halo mass, between the numbers of red and blue satellites and galaxies with stellar masses below $10^{10} \Msun. Possible causes of these correlations are the selection of similar galaxies in different samples, the survey flux limit, or physical mechanisms like a fixed ratio between the satellite numbers of distinct populations. The decorrelation for halos with smaller masses is probably an effect of shot noise by low-occupancy halos. The inferred HODs can be used to complement galaxy-galaxy-lensing or galaxy clustering HOD studies or as input to cosmological analyses and improved mock galaxy catalogues. △ Less

Submitted 5 April, 2022; originally announced April 2022.

Comments: 20 pages + Appendix, 14 Figures. Submitted to Astronomy & Astrophysics. Abstract is abridged

Journal ref: A&A 665, A38 (2022)

KiDS-1000: cosmic shear with enhanced redshift calibration

Authors: J. L. van den Busch, A. H. Wright, H. Hildebrandt, M. Bilicki, M. Asgari, S. Joudaki, C. Blake, C. Heymans, A. Kannawadi, H. Y. Shan, T. Tröster

Abstract: We present a cosmic shear analysis with an improved redshift calibration for the fourth data release of the Kilo-Degree Survey (KiDS-1000) using self-organising maps (SOMs). Compared to the previous analysis of the KiDS-1000 data, we expand the redshift calibration sample to more than twice its size, now consisting of data of 17 spectroscopic redshift campaigns, and significantly extending the fra… ▽ More We present a cosmic shear analysis with an improved redshift calibration for the fourth data release of the Kilo-Degree Survey (KiDS-1000) using self-organising maps (SOMs). Compared to the previous analysis of the KiDS-1000 data, we expand the redshift calibration sample to more than twice its size, now consisting of data of 17 spectroscopic redshift campaigns, and significantly extending the fraction of KiDS galaxies we are able to calibrate with our SOM redshift methodology. We then enhance the calibration sample with precision photometric redshifts from COSMOS2015 and the Physics of the Accelerated Universe Survey (PAUS), allowing us to fill gaps in the spectroscopic coverage of the KiDS data. Finally we perform a Complete Orthogonal Sets of E/B-Integrals (COSEBIs) cosmic shear analysis of the newly calibrated KiDS sample. We find $S_8 = 0.748_{-0.025}^{+0.021}$, which is in good agreement with previous KiDS studies and increases the tension with measurements of the cosmic microwave background to 3.4σ. We repeat the redshift calibration with different subsets of the full calibration sample and obtain, in all cases, agreement within at most 0.5σ in $S_8$ compared to our fiducial analysis. Including additional photometric redshifts allows us to calibrate an additional 6 % of the source galaxy sample. Even though further systematic testing with simulated data is necessary to quantify the impact of redshift outliers, precision photometric redshifts can be beneficial at high redshifts and to mitigate selection effects commonly found in spectroscopically selected calibration samples. △ Less

Submitted 5 April, 2022; originally announced April 2022.

Comments: 18 pages, 15 figures, 6 tables

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

Cosmic star formation history with tomographic cosmic infrared background-galaxy cross-correlation

Authors: Ziang Yan, Ludovic van Waerbeke, Angus H. Wright, Maciej Bilicki, Shiming Gu, Hendrik Hildebrandt, Abhishek S. Maniyar, Tilman Tröster

Abstract: In this work, we probe the star formation history of the Universe using tomographic cross-correlation between the cosmic infrared background (CIB) and galaxy samples. The galaxy samples are from the Kilo-Degree Survey (KiDS), while the CIB maps are made from \planck\, sky maps. We measure the cross-correlation in harmonic space with a significance of 43$σ$. We model the cross-correlation with a ha… ▽ More In this work, we probe the star formation history of the Universe using tomographic cross-correlation between the cosmic infrared background (CIB) and galaxy samples. The galaxy samples are from the Kilo-Degree Survey (KiDS), while the CIB maps are made from \planck\, sky maps. We measure the cross-correlation in harmonic space with a significance of 43$σ$. We model the cross-correlation with a halo model, which links CIB anisotropies to star formation rates (SFR) and galaxy abundance. We assume that SFR has a lognormal dependence on halo mass, while galaxy abundance follows the halo occupation distribution (HOD) model. The cross-correlations give a best-fit maximum star formation efficiency of $η_{\mathrm{max}}= 0.41^{+0.09}_{-0.14}$ at a halo mass $\log_{10}(M_{\mathrm{peak}}/M_{\odot})= {12.14\pm 0.36}$. The derived star formation rate density (SFRD) is well constrained up to $z\sim 1.5$. The constraining power at high redshift is mainly limited by the KiDS survey depth. A combination with external SFRD measurements from previous studies gives $\log_{10}(M_{\mathrm{peak}}/M_{\odot})=12.42^{+0.35}_{-0.19}$. This tightens the SFRD constraint up to $z=4$, yielding a peak SFRD of $0.09_{-0.004}^{+0.003}\,M_{\odot} \mathrm { year }^{-1} \mathrm{Mpc}^{-3}$ at $z=1.74^{+0.06}_{-0.02}$, corresponding to a lookback time of $10.05^{+0.12}_{-0.03}$ Gyr. Both constraints are consistent, and the derived SFRD agrees with previous studies and simulations. Additionally, we estimate the galaxy bias $b$ of KiDS galaxies from the constrained HOD parameters and yield an increasing bias from $b=1.1_{-0.31}^{+0.17}$ at $z=0$ to $b=1.96_{-0.64}^{+0.18}$ at $z=1.5$. Finally, we provide a forecast for future galaxy surveys and conclude that, due to their considerable depth, future surveys will yield a much tighter constraint on the evolution of the SFRD. △ Less

Submitted 18 November, 2022; v1 submitted 4 April, 2022; originally announced April 2022.

Comments: 22 pages, 14 figures, 3 tables, the abstract is abridged, published on A&A

Journal ref: A&A 665, A52 (2022)

Constraints on galaxy formation from the cosmic-infrared-background$\,-\,$optical-imaging cross-correlation using $\boldsymbol{Herschel}$ and UNIONS

Authors: Seunghwan Lim, Ryley Hill, Douglas Scott, Ludovic van Waerbeke, Jean-Charles Cuillandre, Raymond G. Carlberg, Nora Elisa Chisari, Andrej Dvornik, Thomas Erben, Stephen Gwyn, Alan W. McConnachie, Marc-Antoine Miville-Deschênes, Angus H. Wright, Pierre-Alain Duc

Abstract: Using $\it{Herschel}$-SPIRE imaging and the Canada-France Imaging Survey (CFIS) Low Surface Brightness data products from the Ultraviolet Near-Infrared Optical Northern Survey (UNIONS), we present a cross-correlation between the cosmic infrared background and cosmic optical background fluctuations. With a combined sky area of $91\,{\rm deg}^2$ and a minimum resolved scale of $18\,$arcsec, the cros… ▽ More Using $\it{Herschel}$-SPIRE imaging and the Canada-France Imaging Survey (CFIS) Low Surface Brightness data products from the Ultraviolet Near-Infrared Optical Northern Survey (UNIONS), we present a cross-correlation between the cosmic infrared background and cosmic optical background fluctuations. With a combined sky area of $91\,{\rm deg}^2$ and a minimum resolved scale of $18\,$arcsec, the cross-spectrum is measured for two cases: all galaxies are kept in the images; or all individually-detected galaxies are masked to produce `background' maps. We report the detection of the cross-correlation signal at $\gtrsim 22\,σ$ ($\gtrsim 18\,σ$ for the background map). The part of the optical brightness variations that are correlated with the submm emission translates to an rms brightness of $\simeq 32.5\,{\rm mag}\,{\rm arcsec}^{-2}$ in the $r$ band, a level normally unreachable for individual sources. A critical issue is determining what fraction of the cross-power spectrum might be caused by emission from Galactic cirrus. For one of the fields, the Galactic contamination is approximately a factor of 10 higher than the extragalactic signal, with the contamination being estimated using a linear regression from several external survey maps; however, for the other fields, the contamination is typically around 20 per cent. An additional discriminant is that the cross-power spectrum is of the approximate form $P(k)\propto 1/k$, much shallower than that of Galactic cirrus. We interpret the results in a halo-model framework, which shows good agreement with independent measurements for the scalings of star-formation rates in galaxies. The approach presented in this study holds great promise for future surveys such as FYST/CCAT-prime combined with ${\it Euclid}$ or the Vera Rubin Observatory (LSST), which will enable a detailed exploration of the evolution of star formation in galaxies. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: 35 pages, 27 figures, submitted to MNRAS, comments are welcome

Galaxy And Mass Assembly (GAMA): Self-Organizing Map Application on Nearby Galaxies

Authors: B. W. Holwerda, Dominic Smith, Lori Porter, Chris Henry, Ren Porter-Temple, Kyle Cook, Kevin A. Pimbblet, Andrew M. Hopkins, Maciej Bilicki, Sebastian Turner, Viviana Acquaviva, Lingyu Wang, Angus H. Wright, Lee S. Kelvin, Meiert W. Grootes

Abstract: Galaxy populations show bimodality in a variety of properties: stellar mass, colour, specific star-formation rate, size, and Sérsic index. These parameters are our feature space. We use an existing sample of 7556 galaxies from the Galaxy and Mass Assembly (GAMA) survey, represented using five features and the K-means clustering technique, showed that the bimodalities are the manifestation of a mor… ▽ More Galaxy populations show bimodality in a variety of properties: stellar mass, colour, specific star-formation rate, size, and Sérsic index. These parameters are our feature space. We use an existing sample of 7556 galaxies from the Galaxy and Mass Assembly (GAMA) survey, represented using five features and the K-means clustering technique, showed that the bimodalities are the manifestation of a more complex population structure, represented by between 2 and 6 clusters. Here we use Self Organizing Maps (SOM), an unsupervised learning technique which can be used to visualize similarity in a higher dimensional space using a 2D representation, to map these five-dimensional clusters in the feature space onto two-dimensional projections. To further analyze these clusters, using the SOM information, we agree with previous results that the sub-populations found in the feature space can be reasonably mapped onto three or five clusters. We explore where the "green valley" galaxies are mapped onto the SOM, indicating multiple interstitial populations within the green valley population. Finally, we use the projection of the SOM to verify whether morphological information provided by GalaxyZoo users, for example, if features are visible, can be mapped onto the SOM-generated map. Voting on whether galaxies are smooth, likely ellipticals, or "featured" can reasonably be separated but smaller morphological features (bar, spiral arms) can not. SOMs promise to be a useful tool to map and identify instructive sub-populations in multidimensional galaxy survey feature space, provided they are large enough. △ Less

Submitted 29 March, 2022; originally announced March 2022.

Comments: 14 pages, 14 figures, accepted by MNRAS

Galaxy And Mass Assembly (GAMA): Data Release 4 and the z < 0.1 total and z < 0.08 morphological galaxy stellar mass functions

Authors: Simon P. Driver, Sabine Bellstedt, Aaron S. G. Robotham, Ivan K. Baldry, Luke J. Davies, Jochen Liske, Danail Obreschkow, Edward N. Taylor, Angus H. Wright, Mehmet Alpaslan, Steven P. Bamford, Amanda E. Bauer, Joss Bland-Hawthorn, Maciej Bilicki, Matias Bravo, Sarah Brough, Sarah Casura, Michelle E. Cluver, Matthew Colless, Christopher J. Conselice, Scott M. Croom, Jelte de Jong, Franceso D'Eugenio, Roberto De Propris, Burak Dogruel , et al. (45 additional authors not shown)

Abstract: In Galaxy And Mass Assembly Data Release 4 (GAMA DR4), we make available our full spectroscopic redshift sample. This includes 248682 galaxy spectra, and, in combination with earlier surveys, results in 330542 redshifts across five sky regions covering ~250deg^2. The redshift density, is the highest available over such a sustained area, has exceptionally high completeness (95 per cent to r_KIDS=19… ▽ More In Galaxy And Mass Assembly Data Release 4 (GAMA DR4), we make available our full spectroscopic redshift sample. This includes 248682 galaxy spectra, and, in combination with earlier surveys, results in 330542 redshifts across five sky regions covering ~250deg^2. The redshift density, is the highest available over such a sustained area, has exceptionally high completeness (95 per cent to r_KIDS=19.65mag), and is well suited for the study of galaxy mergers, galaxy groups, and the low redshift (z<0.25) galaxy population. DR4 includes 32 value-added tables or Data Management Units (DMUs) that provide a number of measured and derived data products including GALEX, ESO KiDS, ESO VIKING, WISE and Herschel Space Observatory imaging. Within this release, we provide visual morphologies for 15330 galaxies to z<0.08, photometric redshift estimates for all 18million objects to r_KIDS~25mag, and stellar velocity dispersions for 111830 galaxies. We conclude by deriving the total galaxy stellar mass function (GSMF) and its sub-division by morphological class (elliptical, compact-bulge and disc, diffuse-bulge and disc, and disc only). This extends our previous measurement of the total GSMF down to 10^6.75 M_sol h^-2_70 and we find a total stellar mass density of rho_*=(2.97+/-0.04)x10^8 M_sol h_70 Mpc^-3 or Omega_*=(2.17+/-0.03)x10^-3 h^-1_70. We conclude that at z<0.1, the Universe has converted 4.9+/-0.1 per cent of the baryonic mass implied by Big Bang Nucleosynthesis into stars that are gravitationally bound within the galaxy population. △ Less

Submitted 16 March, 2022; originally announced March 2022.

Comments: Accepted for publication in MNRAS. GAMA Data Release 4 is available at: http://www.gama-survey.org/dr4/

Consistent lensing and clustering in a low-$S_8$ Universe with BOSS, DES Year 3, HSC Year 1 and KiDS-1000

Authors: A. Amon, N. C. Robertson, H. Miyatake, C. Heymans, M. White, J. DeRose, S. Yuan, R. H. Wechsler, T. N. Varga, S. Bocquet, A. Dvornik, S. More, A. J. Ross, H. Hoekstra, A. Alarcon, M. Asgari, J. Blazek, A. Campos, R. Chen, A. Choi, M. Crocce, H. T. Diehl, C. Doux, K. Eckert, J. Elvin-Poole , et al. (83 additional authors not shown)

Abstract: We evaluate the consistency between lensing and clustering probes of large-scale structure based on measurements of projected galaxy clustering from BOSS combined with overlapping galaxy-galaxy lensing from three surveys: DES Y3, HSC Y1, and KiDS-1000. An intra-lensing-survey study finds good agreement between these lensing data. We model the observations using the Dark Emulator and fit the data a… ▽ More We evaluate the consistency between lensing and clustering probes of large-scale structure based on measurements of projected galaxy clustering from BOSS combined with overlapping galaxy-galaxy lensing from three surveys: DES Y3, HSC Y1, and KiDS-1000. An intra-lensing-survey study finds good agreement between these lensing data. We model the observations using the Dark Emulator and fit the data at two fixed cosmologies: Planck, with $S_8=0.83$, and a Lensing cosmology with $S_8=0.76$. For a joint analysis limited to scales with $R>5.25h^{-1}$Mpc, we find that both cosmologies provide an acceptable fit to the data. Full utilisation of the small-scale clustering and lensing measurements is hindered by uncertainty in the impact of baryon feedback and assembly bias, which we account for with a reasoned theoretical error budget. We incorporate a systematic scaling parameter for each redshift bin, $A$, that decouples the lensing and clustering to capture any inconsistency. When a wide range of scales ($0.15<R<60h^{-1}$Mpc) are incorporated, we find different results for the consistency of clustering and lensing between the two cosmologies. Limiting the analysis to the bins for which the impact of the selection of the lens sample is expected to be minimal, for the low-$S_8$ Lensing cosmology, the measurements are consistent with $A$=1; $A=0.91\pm0.04$ using DES+KiDS and $A=0.97\pm0.06$ using HSC. For the Planck cosmology case, we find a discrepancy: $A=0.79\pm0.03$ using DES+KiDS and $A=0.84\pm0.05$ using HSC. We demonstrate that a kSZ-based estimate for baryonic effects alleviates some of the discrepancy in the Planck cosmology. This analysis demonstrates the statistical power of these small-scale measurements, but also indicates that caution is still warranted given current uncertainties in modelling baryonic effects, assembly bias, and selection effects in the foreground sample. △ Less

Submitted 13 October, 2022; v1 submitted 15 February, 2022; originally announced February 2022.

Comments: 28 pages, 11 figures

The PAU Survey: Measurements of the 4000 Å spectral break with narrow-band photometry

Authors: Pablo Renard, Małgorzata Siudek, Martin B. Eriksen, Laura Cabayol, Zheng Cai, Jorge Carretero, Ricard Casas, Francisco J. Castander, Enrique Fernandez, Juan García-Bellido, Enrique Gaztanaga, Henk Hoekstra, Benjamin Joachimi, Ramon Miquel, David Navarro-Girones, Cristóbal Padilla, Eusebio Sanchez, Santiago Serrano, Pau Tallada-Crespí, Juan De Vicente, Anna Wittje, Angus H. Wright

Abstract: The D4000 spectral break index is one of the most important features in the visible spectrum, as it is a proxy for stellar ages and is also used in galaxy classification. However, its direct measurement has always been reserved to spectroscopy. Here, we present a general method to directly measure the D4000 with narrow-band (NB) photometry; it has been validated using realistic simulations, and th… ▽ More The D4000 spectral break index is one of the most important features in the visible spectrum, as it is a proxy for stellar ages and is also used in galaxy classification. However, its direct measurement has always been reserved to spectroscopy. Here, we present a general method to directly measure the D4000 with narrow-band (NB) photometry; it has been validated using realistic simulations, and then evaluated with PAUS NBs, cross-matched with VIPERS spectra ($i_{\rm AB} < 22.5$, $0.562 < z < 0.967$). We also reconstruct the D4000 with the SED-fitting code CIGALE; the use of PAUS NBs instead of broad bands significantly improves the SED fitting results. For D4000$_{\rm n}$, the direct measurement has $\rm \langle SNR \rangle \sim 4$, but we find that for $i_{\rm AB}<21$ all direct D4000 measurements have $\rm SNR>3$. The CIGALE D4000$_{\rm n}$ has $\rm \langle SNR \rangle \sim 20$, but underestimates the error by $>$50\%. Furthermore, the direct method recreates well the D4000-SFR relation, as well as the D4000-mass relation for blue galaxies (for red galaxies, selection effects impact the results). On the other hand, CIGALE accurately classifies galaxies into red and blue populations. We conclude that the direct measurement of D4000 with narrow-band photometry is a promising tool to determine average properties of galaxy samples, with results compatible with spectroscopy. △ Less

Submitted 22 July, 2022; v1 submitted 12 January, 2022; originally announced January 2022.

Comments: 21 pages, 18 figures. Accepted in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 515, Issue 1, September 2022, Pages 146-166

Rubin-Euclid Derived Data Products: Initial Recommendations

Authors: Leanne P. Guy, Jean-Charles Cuillandre, Etienne Bachelet, Manda Banerji, Franz E. Bauer, Thomas Collett, Christopher J. Conselice, Siegfried Eggl, Annette Ferguson, Adriano Fontana, Catherine Heymans, Isobel M. Hook, Éric Aubourg, Hervé Aussel, James Bosch, Benoit Carry, Henk Hoekstra, Konrad Kuijken, Francois Lanusse, Peter Melchior, Joseph Mohr, Michele Moresco, Reiko Nakajima, Stéphane Paltani, Michael Troxel , et al. (95 additional authors not shown)

Abstract: This report is the result of a joint discussion between the Rubin and Euclid scientific communities. The work presented in this report was focused on designing and recommending an initial set of Derived Data products (DDPs) that could realize the science goals enabled by joint processing. All interested Rubin and Euclid data rights holders were invited to contribute via an online discussion forum… ▽ More This report is the result of a joint discussion between the Rubin and Euclid scientific communities. The work presented in this report was focused on designing and recommending an initial set of Derived Data products (DDPs) that could realize the science goals enabled by joint processing. All interested Rubin and Euclid data rights holders were invited to contribute via an online discussion forum and a series of virtual meetings. Strong interest in enhancing science with joint DDPs emerged from across a wide range of astrophysical domains: Solar System, the Galaxy, the Local Volume, from the nearby to the primaeval Universe, and cosmology. △ Less

Submitted 13 October, 2022; v1 submitted 11 January, 2022; originally announced January 2022.

Comments: Report of the Rubin-Euclid Derived Data Products Working Group, 78 pages, 11 figures

The variation of the gas content of galaxy groups and pairs compared to isolated galaxies

Authors: Sambit Roychowdhury, Martin J. Meyer, Jonghwan Rhee, Martin A. Zwaan, Garima Chauhan, Luke J. M. Davies, Sabine Bellstedt, Simon P. Driver, Claudia del P. Lagos, Aaron S. G. Robotham, Joss Bland-Hawthorn, Richard Dodson, Benne W. Holwerda, Andrew M. Hopkins, Maritza A. Lara-Lopez, Angel R. Lopez-Sanchez, Danail Obreschkow, Kristof Rozgonyi, Matthew T. Whiting, Angus H. Wright

Abstract: We measure how the atomic gas (HI) fraction ($f_{HI}={\rm \frac{M_{HI}}{M_{*}}}$) of groups and pairs taken as single units vary with average stellar mass ($\langle {\rm M_*} \rangle$) and average star-formation rate ($\langle {\rm SFR} \rangle$), compared to isolated galaxies. The HI 21 cm emission observation are from (i) archival ALFALFA survey data covering three fields from the GAMA survey (p… ▽ More We measure how the atomic gas (HI) fraction ($f_{HI}={\rm \frac{M_{HI}}{M_{*}}}$) of groups and pairs taken as single units vary with average stellar mass ($\langle {\rm M_*} \rangle$) and average star-formation rate ($\langle {\rm SFR} \rangle$), compared to isolated galaxies. The HI 21 cm emission observation are from (i) archival ALFALFA survey data covering three fields from the GAMA survey (provides environmental and galaxy properties), and (ii) DINGO pilot survey data of one of those fields. The mean $f_{HI}$ for different units (groups/pairs/isolated galaxies) are measured in regions of the log($\langle {\rm M_*} \rangle$) -- log($\langle {\rm SFR} \rangle$) plane, relative to the z $\sim 0$ star-forming main sequence (SFMS) of individual galaxies, by stacking $f_{HI}$ spectra of individual units. For ALFALFA, $f_{HI}$ spectra of units are measured by extracting HI spectra over the full groups/pair areas and dividing by the total stellar mass of member galaxies. For DINGO, $f_{HI}$ spectra of units are measured by co-adding HI spectra of individual member galaxies, followed by division by their total stellar mass. For all units the mean $f_{HI}$ decreases as we move to higher $\langle {\rm M_*} \rangle$ along the SFMS, and as we move from above the SFMS to below it at any $\langle {\rm M_*} \rangle$. From the DINGO-based study, mean $f_{HI}$ in groups appears to be lower compared to isolated galaxies for all $\langle {\rm M_*} \rangle$ along the SFMS. From the ALFALFA-based study we find substantially higher mean $f_{HI}$ in groups compared to isolated galaxies (values for pairs being intermediate) for ${\langle{\rm M_*}\rangle}\lesssim10^{9.5}~{\rm M_{\odot}}$, indicating the presence of substantial amounts of HI not associated with cataloged member galaxies in low mass groups. △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: Accepted for publication in ApJ. Main text: 26 pages, 16 figures, 6 tables

Galaxy and Mass Assembly (GAMA): The Weak Environmental Dependence of Quasar Activity at 0.1<z<0.35

Authors: Clare F. Wethers, Nischal Acharya, Roberto De Propris, Jari Kotilainen, Ivan K. Baldry, Sarah Brough, Simon P. Driver, Alister W. Graham, Benne W. Holwerda, Andrew M. Hopkins Angel R. López-Sánchez, Jonathan Loveday, Steven Phillipps, Kevin A. Pimbblet, Edward Taylor, Lingyu Wang, Angus H. Wright

Abstract: Understanding the connection between nuclear activity and galaxy environment remains critical in constraining models of galaxy evolution. By exploiting extensive catalogued data from the Galaxy and Mass Assembly (GAMA) survey, we identify a representative sample of 205 quasars at 0.1 < z < 0.35 and establish a comparison sample of galaxies, closely matched to the quasar sample in terms of both ste… ▽ More Understanding the connection between nuclear activity and galaxy environment remains critical in constraining models of galaxy evolution. By exploiting extensive catalogued data from the Galaxy and Mass Assembly (GAMA) survey, we identify a representative sample of 205 quasars at 0.1 < z < 0.35 and establish a comparison sample of galaxies, closely matched to the quasar sample in terms of both stellar mass and redshift. On scales <1 Mpc, the galaxy number counts and group membership of quasars appear entirely consistent with those of the matched galaxy sample. Despite this, we find that quasars are ~1.5 times more likely to be classified as the group center, indicating a potential link between quasar activity and cold gas flows or galaxy interactions associated with rich group environments. On scales of ~a few Mpc, the clustering strength of both samples are statistically consistent and beyond 10 Mpc we find no evidence that quasars trace large scale structures any more than the galaxy control sample. Both populations are found to prefer intermediate-density sheets and filaments to either very high- or very low- density environments. This weak dependence of quasar activity on galaxy environment supports a paradigm in which quasars represent a phase in the lifetime of all massive galaxies and in which secular processes and a group-centric location are the dominant trigger of quasars at low redshift. △ Less

Submitted 6 December, 2021; originally announced December 2021.

Comments: 26 pages, 20 figures, accepted to ApJ

Lensing Without Borders. I. A Blind Comparison of the Amplitude of Galaxy-Galaxy Lensing Between Independent Imaging Surveys

Authors: A. Leauthaud, A. Amon, S. Singh, D. Gruen, J. U. Lange, S. Huang, N. C. Robertson, T. N. Varga, Y. Luo, C. Heymans, H. Hildebrandt, C. Blake, M. Aguena, S. Allam, F. Andrade-Oliveira, J. Annis, E. Bertin, S. Bhargava, J. Blazek, S. L. Bridle, D. Brooks, D. L. Burke, A. Carnero Rosell, M. Carrasco Kind, J. Carretero , et al. (82 additional authors not shown)

Abstract: Lensing Without Borders is a cross-survey collaboration created to assess the consistency of galaxy-galaxy lensing signals ($ΔΣ$) across different data-sets and to carry out end-to-end tests of systematic errors. We perform a blind comparison of the amplitude of $ΔΣ$ using lens samples from BOSS and six independent lensing surveys. We find good agreement between empirically estimated and reported… ▽ More Lensing Without Borders is a cross-survey collaboration created to assess the consistency of galaxy-galaxy lensing signals ($ΔΣ$) across different data-sets and to carry out end-to-end tests of systematic errors. We perform a blind comparison of the amplitude of $ΔΣ$ using lens samples from BOSS and six independent lensing surveys. We find good agreement between empirically estimated and reported systematic errors which agree to better than 2.3$σ$ in four lens bins and three radial ranges. For lenses with $z_{\rm L}>0.43$ and considering statistical errors, we detect a 3-4$σ$ correlation between lensing amplitude and survey depth. This correlation could arise from the increasing impact at higher redshift of unrecognised galaxy blends on shear calibration and imperfections in photometric redshift calibration. At $z_{\rm L}>0.54$ amplitudes may additionally correlate with foreground stellar density. The amplitude of these trends is within survey-defined systematic error budgets which are designed to include known shear and redshift calibration uncertainty. Using a fully empirical and conservative method, we do not find evidence for large unknown systematics. Systematic errors greater than 15% (25%) ruled out in three lens bins at 68% (95%) confidence at $z<0.54$. Differences with respect to predictions based on clustering are observed to be at the 20-30% level. Our results therefore suggest that lensing systematics alone are unlikely to fully explain the "lensing is low" effect at $z<0.54$. This analysis demonstrates the power of cross-survey comparisons and provides a promising path for identifying and reducing systematics in future lensing analyses. △ Less

Submitted 26 November, 2021; originally announced November 2021.

Comments: 41 page, 20 figures

The dark matter halo masses of elliptical galaxies as a function of observationally robust quantities

Authors: Alessandro Sonnenfeld, Crescenzo Tortora, Henk Hoekstra, Marika Asgari, Maciej Bilicki, Catherine Heymans, Hendrik Hildebrandt, Koen Kuijken, Nicola R. Napolitano, Nivya Roy, Edwin Valentijn, Angus H. Wright

Abstract: Context. The assembly history of the stellar component of a massive elliptical galaxy is closely related to that of its dark matter halo. Measuring how the properties of galaxies correlate with their halo mass can help understand their evolution. Aims. We investigate how the dark matter halo mass of elliptical galaxies varies as a function of their properties, using weak gravitational lensing ob… ▽ More Context. The assembly history of the stellar component of a massive elliptical galaxy is closely related to that of its dark matter halo. Measuring how the properties of galaxies correlate with their halo mass can help understand their evolution. Aims. We investigate how the dark matter halo mass of elliptical galaxies varies as a function of their properties, using weak gravitational lensing observations. To minimise the chances of biases, we focus on galaxy properties that can be determined robustly: the surface brightness profile and the colour. Methods. We selected 2409 central massive elliptical galaxies from the SDSS spectroscopic sample. We first measured their surface brightness profile and colours by fitting Sersic models to photometric data from the Kilo-Degree Survey (KiDS). We fitted their halo mass distribution as a function of redshift, rest-frame $r-$band luminosity, half-light radius and rest-frame $u-g$ colour, using KiDS weak lensing data and a Bayesian hierarchical approach. For the sake of robustness to assumptions on the large-radii behaviour of the surface brightness, we repeated the analysis replacing total luminosity and half-light radius with the luminosity within a 10~kpc aperture, $L_{r,10}$, and the light-weighted surface brightness slope, $Γ_{10}$. Results. We did not detect any correlation between halo mass and either half-light radius or colour, at fixed redshift and luminosity. Conclusions. Our results indicate that the average star formation efficiency of massive elliptical galaxies has little dependence on their final size or colour. This suggests that the origin of the diversity in the size and colour distribution of these objects lies with properties other than halo mass. △ Less

Submitted 15 April, 2022; v1 submitted 22 October, 2021; originally announced October 2021.

Comments: In press on Astronomy & Astrophysics. A short summary video is available at https://youtu.be/LlvwQNjLT7k

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

Pure-mode correlation functions for cosmic shear and application to KiDS-1000

Authors: Peter Schneider, Marika Asgari, Yasaman Najafi Jozani, Andrej Dvornik, Benjamin Giblin, Joachim Harnois-D/'eraps, Catherine Heymans, Hendrik Hildebrandt, Henk Hoekstra, Konrad Kuijken, HuanYuan Shan, Tilman Tröster, Angus H. Wright

Abstract: One probe for systematic effects in gravitational lensing surveys is the presence of so-called B-modes in the cosmic shear two-point correlation functions ξ_\pm(\vt), since lensing is expected to produce only E-mode shear. Furthermore, there exist ambiguous modes which can not uniquely be assigned to either E- or B-mode shear. We derive explicit equations for the pure-mode shear correlation functi… ▽ More One probe for systematic effects in gravitational lensing surveys is the presence of so-called B-modes in the cosmic shear two-point correlation functions ξ_\pm(\vt), since lensing is expected to produce only E-mode shear. Furthermore, there exist ambiguous modes which can not uniquely be assigned to either E- or B-mode shear. We derive explicit equations for the pure-mode shear correlation functions ξ_\pm^E/B(\vt) and their ambiguous components ξ_\pm^amb(\vt), that can be derived from the measured ξ_\pm(\vt) on a finite angular interval \tmin\le\vt\le\tmax, such that the latter can be decomposed uniquely into pure-mode functions as ξ_+=ξ_+^E+ξ_+^B+ξ_+^amb and ξ_-=ξ_-^E-ξ_-^B+ξ_-^amb. The derivation is obtained by defining a new set of COSEBIs, for which explicit relations are obtained, and which yields a smaller covariance between COSEBI modes. We derive the relation between ξ_\pm^E/B/amb and the underlying E-/B-mode power spectra. The pure-mode correlation functions can provide a diagnostics of systematics in configuration space. We then apply our results to SLICS simulations and the KiDS-1000 cosmic shear data, calculate the new COSEBIs and the pure-mode correlation functions, as well as the corresponding covariances, and show that the new statistics fit equally well to the best-fitting cosmological model as the previous KiDS-1000 analysis and recovers the same level of (insignificant) B-modes. We also consider in some detail the ambiguous modes at first- and second-order level, finding some surprising results; for example, the shear field of a point mass, when cut along a line through the center, can not be ascribed uniquely to an E-mode shear and is thus ambiguous, and the shear correlation functions resulting from a random ensemble of point masses, when measured over a finite angular range, correspond to an ambiguous mode. △ Less

Submitted 17 May, 2022; v1 submitted 19 October, 2021; originally announced October 2021.

Comments: 24 pages, A&A, in press

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

KiDS & Euclid: Cosmological implications of a pseudo angular power spectrum analysis of KiDS-1000 cosmic shear tomography

Authors: A. Loureiro, L. Whittaker, A. Spurio Mancini, B. Joachimi, A. Cuceu, M. Asgari, B. Stölzner, T. Tröster, A. H. Wright, M. Bilicki, A. Dvornik, B. Giblin, C. Heymans, H. Hildebrandt, H. Shan, A. Amara, N. Auricchio, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia, V. Capobianco, C. Carbone, J. Carretero, M. Castellano , et al. (89 additional authors not shown)

Abstract: We present a tomographic weak lensing analysis of the Kilo Degree Survey Data Release 4 (KiDS-1000), using a new pseudo angular power spectrum estimator (pseudo-$C_{\ell}$) under development for the ESA Euclid mission. Over 21 million galaxies with shape information are divided into five tomographic redshift bins, ranging from 0.1 to 1.2 in photometric redshift. We measured pseudo-$C_{\ell}$ using… ▽ More We present a tomographic weak lensing analysis of the Kilo Degree Survey Data Release 4 (KiDS-1000), using a new pseudo angular power spectrum estimator (pseudo-$C_{\ell}$) under development for the ESA Euclid mission. Over 21 million galaxies with shape information are divided into five tomographic redshift bins, ranging from 0.1 to 1.2 in photometric redshift. We measured pseudo-$C_{\ell}$ using eight bands in the multipole range $76<\ell<1500$ for auto- and cross-power spectra between the tomographic bins. A series of tests were carried out to check for systematic contamination from a variety of observational sources including stellar number density, variations in survey depth, and point spread function properties. While some marginal correlations with these systematic tracers were observed, there is no evidence of bias in the cosmological inference. B-mode power spectra are consistent with zero signal, with no significant residual contamination from E/B-mode leakage. We performed a Bayesian analysis of the pseudo-$C_{\ell}$ estimates by forward modelling the effects of the mask. Assuming a spatially flat $Λ$CDM cosmology, we constrained the structure growth parameter $S_8 = σ_8(Ω_{\rm m}/0.3)^{1/2} = 0.754_{-0.029}^{+0.027}$. When combining cosmic shear from KiDS-1000 with baryon acoustic oscillation and redshift space distortion data from recent Sloan Digital Sky Survey (SDSS) measurements of luminous red galaxies, as well as the Lyman-$α$ forest and its cross-correlation with quasars, we tightened these constraints to $S_8 = 0.771^{+0.006}_{-0.032}$. These results are in very good agreement with previous KiDS-1000 and SDSS analyses and confirm a $\sim 3σ$ tension with early-Universe constraints from cosmic microwave background experiments. △ Less

Submitted 4 July, 2022; v1 submitted 13 October, 2021; originally announced October 2021.

Comments: 24 pages, 17 figures, Accepted by A&A

Journal ref: A&A 665, A56 (2022)

High-quality strong lens candidates in the final Kilo Degree survey footprint

Authors: R. Li, N. R. Napolitano, C. Spiniello, C. Tortora, K. Kuijken, L. V. E. Koopmans, P. Schneider, F. Getman, L. Xie, L. Long, W. Shu, G. Vernardos, Z. Huang, G. Covone, A. Dvornik, C. Heymans, H. Hildebrandt, M. Radovich, A. H. Wright

Abstract: We present 97 new high-quality strong lensing candidates found in the final $\sim 350\,\rm deg^2$, that completed the full $\sim 1350\,\rm deg^2$ area of the Kilo-Degree Survey (KiDS). Together with our previous findings, the final list of high-quality candidates from KiDS sums up to 268 systems. The new sample is assembled using a new Convolutional Neural Network (CNN) classifier applied to $r$-b… ▽ More We present 97 new high-quality strong lensing candidates found in the final $\sim 350\,\rm deg^2$, that completed the full $\sim 1350\,\rm deg^2$ area of the Kilo-Degree Survey (KiDS). Together with our previous findings, the final list of high-quality candidates from KiDS sums up to 268 systems. The new sample is assembled using a new Convolutional Neural Network (CNN) classifier applied to $r$-band (best seeing) and $g,~r,~i$ color-composited images separately. This optimizes the complementarity of the morphology and color information on the identification of strong lensing candidates. We apply the new classifiers to a sample of luminous red galaxies (LRGs) and a sample of bright galaxies (BGs) and select candidates that received a high probability to be a lens from the CNN ($P_{\rm CNN}$). In particular, setting $P_{\rm CNN}>0.8$ for the LRGs, the $1$-band CNN predicts 1213 candidates, while the $3$-band classifier yields 1299 candidates, with only $\sim$30\% overlap. For the BGs, in order to minimize the false positives, we adopt a more conservative threshold, $P_{\rm CNN} >0.9$, for both CNN classifiers. This results in 3740 newly selected objects. The candidates from the two samples are visually inspected by 7 co-authors to finally select 97 "high-quality" lens candidates which received mean scores larger than 6 (on a scale from 0 to 10). We finally discuss the effect of the seeing on the accuracy of CNN classification and possible avenues to increase the efficiency of multi-band classifiers, in preparation of next-generation surveys from ground and space. △ Less

Submitted 30 December, 2021; v1 submitted 5 October, 2021; originally announced October 2021.

Comments: Published by APJ

Galaxy And Mass Assembly (GAMA): $\mathbf{z \sim 0}$ Galaxy Luminosity Function down to $\mathbf{L \sim 10^{6}~L_\odot}$ via Clustering Based Redshift Inference

Authors: Geray S. Karademir, Edward N. Taylor, Chris Blake, Ivan K. Baldry, Sabine Bellstedt, Maciej Bilicki, Michael J. I. Brown, Michelle E. Cluver, Simon P. Driver, Hendrik Hildebrandt, Benne W. Holwerda, Andrew M. Hopkins, Jonathan Loveday, Steven Phillipps, Angus H. Wright

Abstract: In this study we present a new experimental design using clustering-based redshift inference to measure the evolving galaxy luminosity function (GLF) spanning 5.5 decades from $L \sim 10^{11.5}$ to $ 10^6 ~ \mathrm{L}_\odot$. We use data from the Galaxy And Mass Assembly (GAMA) survey and the Kilo-Degree Survey (KiDS). We derive redshift distributions in bins of apparent magnitude to the limits of… ▽ More In this study we present a new experimental design using clustering-based redshift inference to measure the evolving galaxy luminosity function (GLF) spanning 5.5 decades from $L \sim 10^{11.5}$ to $ 10^6 ~ \mathrm{L}_\odot$. We use data from the Galaxy And Mass Assembly (GAMA) survey and the Kilo-Degree Survey (KiDS). We derive redshift distributions in bins of apparent magnitude to the limits of the GAMA-KiDS photometric catalogue: $m_r \lesssim 23$; more than a decade in luminosity beyond the limits of the GAMA spectroscopic redshift sample via clustering-based redshift inference. This technique uses spatial cross-correlation statistics for a reference set with known redshifts (in our case, the main GAMA sample) to derive the redshift distribution for the target ensemble. For the calibration of the redshift distribution we use a simple parametrisation with an adaptive normalisation factor over the interval $0.005 < z < 0.48$ to derive the clustering redshift results. We find that the GLF has a relatively constant power-law slope $α\approx -1.2$ for $-17 \lesssim M_r \lesssim -13$, and then appears to steepen sharply for $-13 \lesssim M_r \lesssim -10$. This upturn appears to be where Globular Clusters (GCs) take over to dominate the source counts as a function of luminosity. Thus we have mapped the GLF across the full range of the $z \sim 0$ field galaxy population from the most luminous galaxies down to the GC scale. △ Less

Submitted 3 November, 2021; v1 submitted 8 September, 2021; originally announced September 2021.

Comments: 21 pages, 13 figures, accepted by MNRAS

Joint constraints on cosmology and the impact of baryon feedback: combining KiDS-1000 lensing with the thermal Sunyaev-Zeldovich effect from Planck and ACT

Authors: Tilman Tröster, Alexander J. Mead, Catherine Heymans, Ziang Yan, David Alonso, Marika Asgari, Maciej Bilicki, Andrej Dvornik, Hendrik Hildebrandt, Benjamin Joachimi, Arun Kannawadi, Konrad Kuijken, Peter Schneider, HuanYuan Shan, Ludovic van Waerbeke, Angus H. Wright

Abstract: We conduct a pseudo-$C_\ell$ analysis of the tomographic cross-correlation between 1000 deg$^2$ of weak lensing data from the Kilo-Degree Survey (KiDS-1000) and the thermal Sunyaev-Zeldovich (tSZ) effect measured by Planck and the Atacama Cosmology Telescope (ACT). Using HMx, a halo-model-based approach that consistently models the gas, star, and dark matter components, we are able to derive const… ▽ More We conduct a pseudo-$C_\ell$ analysis of the tomographic cross-correlation between 1000 deg$^2$ of weak lensing data from the Kilo-Degree Survey (KiDS-1000) and the thermal Sunyaev-Zeldovich (tSZ) effect measured by Planck and the Atacama Cosmology Telescope (ACT). Using HMx, a halo-model-based approach that consistently models the gas, star, and dark matter components, we are able to derive constraints on both cosmology and baryon feedback for the first time from these data, marginalising over redshift uncertainties, intrinsic alignment of galaxies, and contamination by the cosmic infrared background (CIB). We find our results to be insensitive to the CIB, while intrinsic alignment provides a small but significant contribution to the lensing--tSZ cross-correlation. The cosmological constraints are consistent with those of other low-redshift probes and prefer strong baryon feedback. The inferred amplitude of the lensing--tSZ cross-correlation signal, which scales as $σ_8(Ω_\mathrm{m}/0.3)^{0.2}$, is low by $\sim 2\,σ$ compared to the primary cosmic microwave background constraints by Planck. The lensing--tSZ measurements are then combined with pseudo-$C_\ell$ measurements of KiDS-1000 cosmic shear into a novel joint analysis, accounting for the full cross-covariance between the probes, providing tight cosmological constraints by breaking parameter degeneracies inherent to both probes. The joint analysis gives an improvement of 40% on the constraint of $S_8=σ_8\sqrt{Ω_\mathrm{m}/0.3}$ over cosmic shear alone, while providing constraints on baryon feedback consistent with hydrodynamical simulations, demonstrating the potential of such joint analyses with baryonic tracers such as the tSZ effect. We discuss remaining modelling challenges that need to be addressed if these baryonic probes are to be included in future precision-cosmology analyses. △ Less

Submitted 8 February, 2022; v1 submitted 9 September, 2021; originally announced September 2021.

Comments: 21 pages, 16 figures. Accepted in A&A. Comments welcome. Code, data vectors, and covariances are available at https://github.com/tilmantroester/KiDS-1000xtSZ

Journal ref: A&A 660, A27 (2022)

KiDS-1000: Constraints on the intrinsic alignment of luminous red galaxies

Authors: Maria Cristina Fortuna, Henk Hoekstra, Harry Johnston, Mohammadjavad Vakili, Arun Kannawadi, Christos Georgiou, Benjamin Joachimi, Angus H. Wright, Marika Asgari, Maciej Bilicki, Catherine Heymans, Hendrik Hildebrandt, Konrad Kuijken, Maximilian Von Wietersheim-Kramsta

Abstract: We constrain the luminosity and redshift dependence of the intrinsic alignment (IA) of a nearly volume-limited sample of luminous red galaxies selected from the fourth public data release of the Kilo-Degree Survey (KiDS-1000). To measure the shapes of the galaxies, we used two complementary algorithms, finding consistent IA measurements for the overlapping galaxy sample. The global significance of… ▽ More We constrain the luminosity and redshift dependence of the intrinsic alignment (IA) of a nearly volume-limited sample of luminous red galaxies selected from the fourth public data release of the Kilo-Degree Survey (KiDS-1000). To measure the shapes of the galaxies, we used two complementary algorithms, finding consistent IA measurements for the overlapping galaxy sample. The global significance of IA detection across our two independent luminous red galaxy samples, with our favoured method of shape estimation, is $\sim10.7σ$. We find no significant dependence with redshift of the IA signal in the range $0.2<z<0.8$, nor a dependence with luminosity below $L_r\lesssim 2.9 \times 10^{10} h^{-2} L_{r,\odot}$. Above this luminosity, however, we find that the IA signal increases as a power law, although our results are also compatible with linear growth within the current uncertainties. This behaviour motivates the use of a broken power law model when accounting for the luminosity dependence of IA contamination in cosmic shear studies. △ Less

Submitted 6 September, 2021; originally announced September 2021.

Comments: 21 pages, 14 figures. Accepted for publication in A&A

The LSST-DESC 3x2pt Tomography Optimization Challenge

Authors: Joe Zuntz, François Lanusse, Alex I. Malz, Angus H. Wright, Anže Slosar, Bela Abolfathi, David Alonso, Abby Bault, Clécio R. Bom, Massimo Brescia, Adam Broussard, Jean-Eric Campagne, Stefano Cavuoti, Eduardo S. Cypriano, Bernardo M. O. Fraga, Eric Gawiser, Elizabeth J. Gonzalez, Dylan Green, Peter Hatfield, Kartheik Iyer, David Kirkby, Andrina Nicola, Erfan Nourbakhsh, Andy Park, Gabriel Teixeira , et al. (3 additional authors not shown)

Abstract: This paper presents the results of the Rubin Observatory Dark Energy Science Collaboration (DESC) 3x2pt tomography challenge, which served as a first step toward optimizing the tomographic binning strategy for the main DESC analysis. The task of choosing an optimal tomographic binning scheme for a photometric survey is made particularly delicate in the context of a metacalibrated lensing catalogue… ▽ More This paper presents the results of the Rubin Observatory Dark Energy Science Collaboration (DESC) 3x2pt tomography challenge, which served as a first step toward optimizing the tomographic binning strategy for the main DESC analysis. The task of choosing an optimal tomographic binning scheme for a photometric survey is made particularly delicate in the context of a metacalibrated lensing catalogue, as only the photometry from the bands included in the metacalibration process (usually riz and potentially g) can be used in sample definition. The goal of the challenge was to collect and compare bin assignment strategies under various metrics of a standard 3x2pt cosmology analysis in a highly idealized setting to establish a baseline for realistically complex follow-up studies; in this preliminary study, we used two sets of cosmological simulations of galaxy redshifts and photometry under a simple noise model neglecting photometric outliers and variation in observing conditions, and contributed algorithms were provided with a representative and complete training set. We review and evaluate the entries to the challenge, finding that even from this limited photometry information, multiple algorithms can separate tomographic bins reasonably well, reaching figures-of-merit scores close to the attainable maximum. We further find that adding the g band to riz photometry improves metric performance by ~15% and that the optimal bin assignment strategy depends strongly on the science case: which figure-of-merit is to be optimized, and which observables (clustering, lensing, or both) are included. △ Less

Submitted 15 October, 2021; v1 submitted 30 August, 2021; originally announced August 2021.

Comments: 30 pages (incl. 12 in appendix), 12 figures. Version accepted for publication in the Open Journal of Astrophysics

Report number: DESC-PUB-00054

Faint polarised sources in the Lockman Hole field at 1.4 GHz

Authors: A. Berger, B. Adebahr, N. Herrera Ruiz, A. H. Wright, I. Prandoni, R. -J. Dettmar

Abstract: We aim to study the nature of the faint, polarised radio source population whose source composition and redshift dependence contain information about the strength, morphology, and evolution of magnetic fields over cosmic timescales. We use a 15 pointing radio continuum L-band mosaic of the Lockman Hole, observed in full polarisation, generated from archival data of the WSRT. The data were analysed… ▽ More We aim to study the nature of the faint, polarised radio source population whose source composition and redshift dependence contain information about the strength, morphology, and evolution of magnetic fields over cosmic timescales. We use a 15 pointing radio continuum L-band mosaic of the Lockman Hole, observed in full polarisation, generated from archival data of the WSRT. The data were analysed using the RM-Synthesis technique. We achieved a noise of 7 μJy/beam in polarised intensity, with a resolution of 15''. Using infrared and optical images and source catalogues, we were able to cross-identify and determine redshifts for one third of our detected polarised sources. We detected 150 polarised sources, most of which are weakly polarised with a mean fractional polarisation of 5.4 %. With a total area of 6.5 deg^2 and a detection threshold of 6.25 σ we find 23 polarised sources per deg^2. Based on our multi wavelength analysis, we find that our sample consists of AGN only. We find a discrepancy between archival number counts and those present in our data, which we attribute to sample variance. Considering the absolute radio luminosty, to distinguish weak and strong sources, we find a general trend of increased probability to detect weak sources at low redshift and strong sources at high redshift. Further, we find an anti-correlation between fractional polarisation and redshift for our strong sources sample at z{\geq}0.6. A decrease in the fractional polarisation of strong sources with increasing redshift cannot be explained by a constant magnetic field and electron density over cosmic scales, however the changing properties of cluster environments over the cosmic timemay play an important role. Disentangling these two effects requires deeper and wider polarisation observations, and better models of the morphology and strength of cosmic magnetic fields. △ Less

Submitted 6 July, 2021; originally announced July 2021.

Comments: 17 pages, 16 figures, to be published in A&A

Journal ref: A&A 653, A155 (2021)

The Weak Lensing Radial Acceleration Relation: Constraining Modified Gravity and Cold Dark Matter theories with KiDS-1000

Authors: Margot M. Brouwer, Kyle A. Oman, Edwin A. Valentijn, Maciej Bilicki, Catherine Heymans, Henk Hoekstra, Nicola R. Napolitano, Nivya Roy, Crescenzo Tortora, Angus H. Wright, Marika Asgari, Jan Luca van den Busch, Andrej Dvornik, Thomas Erben, Benjamin Giblin, Alister W. Graham, Hendrik Hildebrandt, Andrew M. Hopkins, Arun Kannawadi, Konrad Kuijken, Jochen Liske, HuanYuan Shan, Tilman Tröster, Erik Verlinde, Manus Visser

Abstract: We present measurements of the radial gravitational acceleration around isolated galaxies, comparing the expected gravitational acceleration given the baryonic matter with the observed gravitational acceleration, using weak lensing measurements from the fourth data release of the Kilo-Degree Survey. These measurements extend the radial acceleration relation (RAR) by 2 decades into the low-accelera… ▽ More We present measurements of the radial gravitational acceleration around isolated galaxies, comparing the expected gravitational acceleration given the baryonic matter with the observed gravitational acceleration, using weak lensing measurements from the fourth data release of the Kilo-Degree Survey. These measurements extend the radial acceleration relation (RAR) by 2 decades into the low-acceleration regime beyond the outskirts of the observable galaxy. We compare our RAR measurements to the predictions of two modified gravity (MG) theories: MOND and Verlinde's emergent gravity. We find that the measured RAR agrees well with the MG predictions. In addition, we find a difference of at least $6σ$ between the RARs of early- and late-type galaxies (split by Sérsic index and $u-r$ colour) with the same stellar mass. Current MG theories involve a gravity modification that is independent of other galaxy properties, which would be unable to explain this behaviour. The difference might be explained if only the early-type galaxies have significant ($M_{gas} \approx M_*$) circumgalactic gaseous haloes. The observed behaviour is also expected in $Λ$CDM models where the galaxy-to-halo mass relation depends on the galaxy formation history. We find that MICE, a $Λ$CDM simulation with hybrid halo occupation distribution modelling and abundance matching, reproduces the observed RAR but significantly differs from BAHAMAS, a hydrodynamical cosmological galaxy formation simulation. Our results are sensitive to the amount of circumgalactic gas; current observational constraints indicate that the resulting corrections are likely moderate. Measurements of the lensing RAR with future cosmological surveys will be able to further distinguish between MG and $Λ$CDM models if systematic uncertainties in the baryonic mass distribution around galaxies are reduced. △ Less

Submitted 22 June, 2021; originally announced June 2021.

Comments: 31 pages, 15 figures, published in A&A

Journal ref: Astronomy & Astrophysics Volume 650 (2021) ArticleID A113

The PAU Survey: Measurement of Narrow-band galaxy properties with Approximate Bayesian Computation

Authors: Luca Tortorelli, Malgorzata Siudek, Beatrice Moser, Tomasz Kacprzak, Pascale Berner, Alexandre Refregier, Adam Amara, Juan García-Bellido, Laura Cabayol, Jorge Carretero, Francisco J. Castander, Juan De Vicente, Martin Eriksen, Enrique Fernandez, Enrique Gaztanaga, Hendrik Hildebrandt, Benjamin Joachimi, Ramon Miquel, Ignacio Sevilla-Noarbe, Cristóbal Padilla, Pablo Renard, Eusebio Sanchez, Santiago Serrano, Pau Tallada-Crespí, Angus H. Wright

Abstract: Narrow-band imaging surveys allow the study of the spectral characteristics of galaxies without the need of performing their spectroscopic follow-up. In this work, we forward-model the Physics of the Accelerating Universe Survey (PAUS) narrow-band data. The aim is to improve the constraints on the spectral coefficients used to create the galaxy spectral energy distributions (SED) of the galaxy pop… ▽ More Narrow-band imaging surveys allow the study of the spectral characteristics of galaxies without the need of performing their spectroscopic follow-up. In this work, we forward-model the Physics of the Accelerating Universe Survey (PAUS) narrow-band data. The aim is to improve the constraints on the spectral coefficients used to create the galaxy spectral energy distributions (SED) of the galaxy population model in Tortorelli et al. 2020. In that work, the model parameters were inferred from the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) data using Approximate Bayesian Computation (ABC). This led to stringent constraints on the B-band galaxy luminosity function parameters, but left the spectral coefficients only broadly constrained. To address that, we perform an ABC inference using CFHTLS and PAUS data. This is the first time our approach combining forward-modelling and ABC is applied simultaneously to multiple datasets. We test the results of the ABC inference by comparing the narrow-band magnitudes of the observed and simulated galaxies using Principal Component Analysis, finding a very good agreement. Furthermore, we prove the scientific potential of the constrained galaxy population model to provide realistic stellar population properties by measuring them with the SED fitting code \textsc{CIGALE}. We use CFHTLS broad-band and PAUS narrow-band photometry for a flux-limited ($\mathrm{i}<22.5$) sample of galaxies up to redshift $\mathrm{z \sim 0.8}$. We find that properties like stellar masses, star-formation rates, mass-weighted stellar ages and metallicities are in agreement within errors between observations and simulations. Overall, this work shows the ability of our galaxy population model to correctly forward-model a complex dataset such as PAUS and the ability to reproduce the diversity of galaxy properties at the redshift range spanned by CFHTLS and PAUS. △ Less

Submitted 8 November, 2021; v1 submitted 4 June, 2021; originally announced June 2021.

Comments: 53 pages, 17 figures, Accepted for publication in JCAP

Geometry vs growth: Internal consistency of the flat ΛCDM model with KiDS-1000

Authors: J. Ruiz-Zapatero, Benjamin Stölzner, Benjamin Joachimi, Marika Asgari, Maciej Bilicki, Andrej Dvornik, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Arun Kannawadi, Konrad Kuijken, TilmanTröster, Jan Luca van den Busch, Angus H. Wright

Abstract: We carry out a multi-probe self-consistency test of the flat $Λ$CDM model with the aim of exploring potential causes of the reported tensions between high- and low-redshift cosmological observations. We divide the model into two theory regimes determined by the smooth background (geometry) and the evolution of matter density fluctuations (growth), each governed by an independent set of Lambda Cold… ▽ More We carry out a multi-probe self-consistency test of the flat $Λ$CDM model with the aim of exploring potential causes of the reported tensions between high- and low-redshift cosmological observations. We divide the model into two theory regimes determined by the smooth background (geometry) and the evolution of matter density fluctuations (growth), each governed by an independent set of Lambda Cold Dark Matter ($Λ$CDM) cosmological parameters. This extended model is constrained by a combination of weak gravitational lensing measurements from the Kilo-Degree Survey, galaxy clustering signatures extracted from Sloan Digital Sky Survey campaigns and the Six-Degree Field Galaxy Survey, and the angular baryon acoustic scale and the primordial scalar fluctuation power spectrum measured in $\textit{Planck}$ cosmic microwave background (CMB) data. We find strong consistency between the geometry and growth parameters, and with the posterior of standard $Λ$CDM analysis. Tension in the amplitude of matter density fluctuations as measured by the parameter $S_8$ persists at around 3$σ$, with a $1.5\,\%$ constraint of $S_8 = 0.776_{-0.008}^{+0.016}$ for the combined probes. We also observe a less significant preference (at least $2σ$) for higher values of the Hubble constant, $H_0 = 70.5^{+0.7}_{-1.5}\,{\rm km\, s^{-1} Mpc^{-1}}$, as well as for lower values of the total matter density parameter $Ω_{\rm{m}} = 0.289^{+0.007}_{-0.005}$ compared to the full $\textit{Planck}$ analysis. Including the subset of the CMB information in the probe combination enhances these differences rather than alleviate them, which we link to the discrepancy between low and high multipoles in $\textit{Planck}$ data. △ Less

Submitted 8 September, 2021; v1 submitted 20 May, 2021; originally announced May 2021.

Comments: 25 pages, 14 figures, updated after being accepted for publication to Astronomy and Astrophysics

Journal ref: A&A 655, A11 (2021)

Probing galaxy bias and intergalactic gas pressure with KiDS Galaxies-tSZ-CMB lensing cross-correlations

Authors: Ziang Yan, Ludovic van Waerbeke, Tilman Tröster, Angus H. Wright, David Alonso, Marika Asgari, Maciej Bilicki, Thomas Erben, Shiming Gu, Catherine Heymans, Hendrik Hildebrandt, Gary Hinshaw, Nick Koukoufilippas, Arun Kannawadi, Konrad Kuijken, Alexander Mead, HuanYuan Shan

Abstract: We constrain the redshift dependence of gas pressure bias $\left\langle b_{y} P_{\mathrm{e}}\right\rangle$ (bias-weighted average electron pressure), which characterises the thermodynamics of intergalactic gas, through a combination of cross-correlations between galaxy positions and the thermal Sunyaev-Zeldovich (tSZ) effect, as well as galaxy positions and the gravitational lensing of the cosmic… ▽ More We constrain the redshift dependence of gas pressure bias $\left\langle b_{y} P_{\mathrm{e}}\right\rangle$ (bias-weighted average electron pressure), which characterises the thermodynamics of intergalactic gas, through a combination of cross-correlations between galaxy positions and the thermal Sunyaev-Zeldovich (tSZ) effect, as well as galaxy positions and the gravitational lensing of the cosmic microwave background (CMB). The galaxy sample is from the fourth data release of the Kilo-Degree Survey (KiDS). The tSZ $y$ map and the CMB lensing map are from the {\textit{Planck}} 2015 and 2018 data releases, respectively. The measurements are performed in five redshift bins with $z\lesssim1$. With these measurements, combining galaxy-tSZ and galaxy-CMB lensing cross-correlations allows us to break the degeneracy between galaxy bias and gas pressure bias, and hence constrain them simultaneously. In all redshift bins, the best-fit values of $\bpe$ are at a level of $\sim 0.3\, \mathrm{meV/cm^3}$ and increase slightly with redshift. The galaxy bias is consistent with unity in all the redshift bins. Our results are not sensitive to the non-linear details of the cross-correlation, which are smoothed out by the {\textit{Planck}} beam. Our measurements are in agreement with previous measurements as well as with theoretical predictions. We also show that our conclusions are not changed when CMB lensing is replaced by galaxy lensing, which shows the consistency of the two lensing signals despite their radically different redshift ranges. This study demonstrates the feasibility of using CMB lensing to calibrate the galaxy distribution such that the galaxy distribution can be used as a mass proxy without relying on the precise knowledge of the matter distribution. △ Less

Submitted 17 May, 2021; v1 submitted 15 February, 2021; originally announced February 2021.

Comments: 20 pages, 14 figures, 3 tables, accepted for publication on Astronomy & Astrophysics

Journal ref: A&A 651, A76 (2021)

Galaxy and Mass Assembly (GAMA): Tracing galaxy environment using the marked correlation function

Authors: U. Sureshkumar, A. Durkalec, A. Pollo, M. Bilicki, J. Loveday, D. J. Farrow, B. W. Holwerda, A. M. Hopkins, J. Liske, K. A. Pimbblet, E. N. Taylor, A. H. Wright

Abstract: We investigate how different galaxy properties - luminosities in u, g, r, J, K-bands, stellar mass, star formation rate and specific star formation rate trace the environment in the local universe. We also study the effect of survey flux limits on galaxy clustering measurements. We measure the two-point correlation function (2pCF) and marked correlation functions (MCFs) using the aforementioned pr… ▽ More We investigate how different galaxy properties - luminosities in u, g, r, J, K-bands, stellar mass, star formation rate and specific star formation rate trace the environment in the local universe. We also study the effect of survey flux limits on galaxy clustering measurements. We measure the two-point correlation function (2pCF) and marked correlation functions (MCFs) using the aforementioned properties as marks. We use nearly stellar-mass-complete galaxy sample in the redshift range 0.1 < z < 0.16 from the Galaxy And Mass Assembly (GAMA) survey with a flux limit of r < 19.8. Further, we impose a brighter flux limit of r < 17.8 to our sample and repeat the measurements to study how this affects galaxy clustering analysis. We compare our results to measurements from the Sloan Digital Sky Survey (SDSS) with flux limits of r < 17.8 and r < 16.8. We show that the stellar mass is the best tracer of galaxy environment, the K-band luminosity being a good substitute, although such a proxy sample misses close pairs of evolved, red galaxies. We also confirm that the u-band luminosity is a good, but not a perfect proxy of star formation rate in the context of galaxy clustering. We observe an effect of the survey flux limit on clustering studies - samples with a higher flux limit (smaller magnitude) miss some information about close pairs of starburst galaxies. △ Less

Submitted 19 September, 2021; v1 submitted 8 February, 2021; originally announced February 2021.

Comments: 15 pages, 7 figures, matches the version published in A&A

Journal ref: A&A 653, A35 (2021)

Halo shapes constrained from a pure sample of central galaxies in KiDS-1000

Authors: Christos Georgiou, Henk Hoekstra, Konrad Kuijken, Maciej Bilicki, Andrej Dvornik, Thomas Erben, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Jelte T. A. de Jong, Arun Kannawadi, Peter Schneider, Tim Schrabback, HuanYuan Shan, Angus H. Wright

Abstract: We present measurements of $f_h$, the ratio of the aligned components of the projected halo and galaxy ellipticities, for a sample of central galaxies using weak gravitational lensing data from the Kilo-Degree Survey (KiDS). Using a lens galaxy shape estimation that is more sensitive to outer galaxy regions, we find $f_{\rm h}=0.50\pm0.20$ for our full sample and $f_{\rm h}=0.55\pm0.19$ for an int… ▽ More We present measurements of $f_h$, the ratio of the aligned components of the projected halo and galaxy ellipticities, for a sample of central galaxies using weak gravitational lensing data from the Kilo-Degree Survey (KiDS). Using a lens galaxy shape estimation that is more sensitive to outer galaxy regions, we find $f_{\rm h}=0.50\pm0.20$ for our full sample and $f_{\rm h}=0.55\pm0.19$ for an intrinsically red (and therefore higher stellar-mass) sub-sample, rejecting the hypothesis of round halos and/or galaxies being un-aligned with their parent halo at $2.5σ$ and $2.9σ$, respectively. We quantify the 93.4% purity of our central galaxy sample using numerical simulations and overlapping spectroscopy from the Galaxy and Mass Assembly survey. This purity ensures that the interpretation of our measurements is not complicated by the presence of a significant fraction of satellite galaxies. Restricting our central galaxy ellipticity measurement to the inner isophotes, we find $f_{\rm h}=0.34\pm0.17$ for our red sub-sample, suggesting that the outer galaxy regions are more aligned with their dark matter halos compared to the inner regions. Our results are in agreement with previous studies and suggest that lower mass halos are rounder and/or less aligned with their host galaxy than samples of more massive galaxies, studied in galaxy groups and clusters. △ Less

Submitted 17 May, 2021; v1 submitted 6 February, 2021; originally announced February 2021.

Comments: Accepted for publication in A&A

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

Bright galaxy sample in the Kilo-Degree Survey Data Release 4: selection, photometric redshifts, and physical properties

Authors: M. Bilicki, A. Dvornik, H. Hoekstra, A. H. Wright, N. E. Chisari, M. Vakili, M. Asgari, B. Giblin, C. Heymans, H. Hildebrandt, B. W. Holwerda, A. Hopkins, H. Johnston, A. Kannawadi, K. Kuijken, S. J. Nakoneczny, H. Y. Shan, A. Sonnenfeld, E. Valentijn

Abstract: We present a bright galaxy sample with accurate and precise photometric redshifts (photo-zs), selected using $ugriZYJHK_\mathrm{s}$ photometry from the Kilo-Degree Survey (KiDS) Data Release 4 (DR4). The highly pure and complete dataset is flux-limited at $r<20$ mag, covers $\sim1000$ deg$^2$, and contains about 1 million galaxies after artifact masking. We exploit the overlap with Galaxy And Mass… ▽ More We present a bright galaxy sample with accurate and precise photometric redshifts (photo-zs), selected using $ugriZYJHK_\mathrm{s}$ photometry from the Kilo-Degree Survey (KiDS) Data Release 4 (DR4). The highly pure and complete dataset is flux-limited at $r<20$ mag, covers $\sim1000$ deg$^2$, and contains about 1 million galaxies after artifact masking. We exploit the overlap with Galaxy And Mass Assembly (GAMA) spectroscopy as calibration to determine photo-zs with the supervised machine learning neural network algorithm implemented in the ANNz2 software. The photo-zs have mean error of $|\langle δz \rangle| \sim 5 \times 10^{-4}$ and low scatter (scaled mean absolute deviation of $\sim 0.018(1+z)$), both practically independent of the $r$-band magnitude and photo-z at $0.05 < z_\mathrm{phot} < 0.5$. Combined with the 9-band photometry, these allow us to estimate robust absolute magnitudes and stellar masses for the full sample. As a demonstration of the usefulness of these data we split the dataset into red and blue galaxies, use them as lenses and measure the weak gravitational lensing signal around them for five stellar mass bins. We fit a halo model to these high-precision measurements to constrain the stellar-mass--halo-mass relations for blue and red galaxies. We find that for high stellar mass ($M_\star>5\times 10^{11} M_\odot$), the red galaxies occupy dark matter halos that are much more massive than those occupied by blue galaxies with the same stellar mass. The data presented here are publicly released via the KiDS webpage at http://kids.strw.leidenuniv.nl/DR4/brightsample.php. △ Less

Submitted 21 July, 2021; v1 submitted 15 January, 2021; originally announced January 2021.

Comments: Matches the published version. Data available at http://kids.strw.leidenuniv.nl/DR4/brightsample.php

Journal ref: A&A 653, A82 (2021)

Magnification bias in galaxy surveys with complex sample selection functions

Authors: Maximilian von Wietersheim-Kramsta, Benjamin Joachimi, Jan Luca van den Busch, Catherine Heymans, Hendrik Hildebrandt, Marika Asgari, Tilman Tröster, Sandra Unruh, Angus H. Wright

Abstract: Gravitational lensing magnification modifies the observed spatial distribution of galaxies and can severely bias cosmological probes of large-scale structure if not accurately modelled. Standard approaches to modelling this magnification bias may not be applicable in practice as many galaxy samples have complex, often implicit, selection functions. We propose and test a procedure to quantify the m… ▽ More Gravitational lensing magnification modifies the observed spatial distribution of galaxies and can severely bias cosmological probes of large-scale structure if not accurately modelled. Standard approaches to modelling this magnification bias may not be applicable in practice as many galaxy samples have complex, often implicit, selection functions. We propose and test a procedure to quantify the magnification bias induced in clustering and galaxy-galaxy lensing (GGL) signals in galaxy samples subject to a selection function beyond a simple flux limit. The method employs realistic mock data to calibrate an effective luminosity function slope, $α_{\rm{obs}}$, from observed galaxy counts, which can then be used with the standard formalism. We demonstrate this method for two galaxy samples derived from the Baryon Oscillation Spectroscopic Survey (BOSS) in the redshift ranges $0.2 < z \leq 0.5$ and $0.5 < z \leq 0.75$, complemented by mock data built from the MICE2 simulation. We obtain $α_{\rm{obs}} = 1.93 \pm 0.05$ and $α_{\rm{obs}} = 2.62 \pm 0.28$ for the two BOSS samples. For BOSS-like lenses, we forecast a contribution of the magnification bias to the GGL signal between the multipole moments, $\ell$, of 100 and 4600 with a cumulative signal-to-noise ratio between 0.1 and 1.1 for sources from the Kilo-Degree Survey (KiDS), between 0.4 and 2.0 for sources from the Hyper Suprime-Cam survey (HSC), and between 0.3 and 2.8 for ESA Euclid-like source samples. These contributions are significant enough to require explicit modelling in future analyses of these and similar surveys. Our code is publicly available within the \textsc{MagBEt} module (\url{https://github.com/mwiet/MAGBET}). △ Less

Submitted 22 April, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: 15 pages, 13 figures. Accepted for publication in Monthly Notices to the Royal Astronomical Society

Euclid preparation: XI. Mean redshift determination from galaxy redshift probabilities for cosmic shear tomography

Authors: Euclid Collaboration, O. Ilbert, S. de la Torre, N. Martinet, A. H. Wright, S. Paltani, C. Laigle, I. Davidzon, E. Jullo, H. Hildebrandt, D. C. Masters, A. Amara, C. J. Conselice, S. Andreon, N. Auricchio, R. Azzollini, C. Baccigalupi, A. Balaguera-Antolínez, M. Baldi, A. Balestra, S. Bardelli, R. Bender, A. Biviano, C. Bodendorf, D. Bonino , et al. (140 additional authors not shown)

Abstract: The analysis of weak gravitational lensing in wide-field imaging surveys is considered to be a major cosmological probe of dark energy. Our capacity to constrain the dark energy equation of state relies on the accurate knowledge of the galaxy mean redshift $\langle z \rangle$. We investigate the possibility of measuring $\langle z \rangle$ with an accuracy better than $0.002\,(1+z)$, in ten tomogr… ▽ More The analysis of weak gravitational lensing in wide-field imaging surveys is considered to be a major cosmological probe of dark energy. Our capacity to constrain the dark energy equation of state relies on the accurate knowledge of the galaxy mean redshift $\langle z \rangle$. We investigate the possibility of measuring $\langle z \rangle$ with an accuracy better than $0.002\,(1+z)$, in ten tomographic bins spanning the redshift interval $0.2<z<2.2$, the requirements for the cosmic shear analysis of Euclid. We implement a sufficiently realistic simulation to understand the advantages, complementarity, but also shortcoming of two standard approaches: the direct calibration of $\langle z \rangle$ with a dedicated spectroscopic sample and the combination of the photometric redshift probability distribution function (zPDF) of individual galaxies. We base our study on the Horizon-AGN hydrodynamical simulation that we analyse with a standard galaxy spectral energy distribution template-fitting code. Such procedure produces photometric redshifts with realistic biases, precision and failure rate. We find that the Euclid current design for direct calibration is sufficiently robust to reach the requirement on the mean redshift, provided that the purity level of the spectroscopic sample is maintained at an extremely high level of $>99.8\%$. The zPDF approach could also be successful if we debias the zPDF using a spectroscopic training sample. This approach requires deep imaging data, but is weakly sensitive to spectroscopic redshift failures in the training sample. We improve the debiasing method and confirm our finding by applying it to real-world weak-lensing data sets (COSMOS and KiDS+VIKING-450). △ Less

Submitted 6 January, 2021; originally announced January 2021.

Comments: 21 pages, 10 figures, accepted in A&A

Organised Randoms: Learning and correcting for systematic galaxy clustering patterns in KiDS using self-organising maps

Authors: Harry Johnston, Angus H. Wright, Benjamin Joachimi, Maciej Bilicki, Nora Elisa Chisari, Andrej Dvornik, Thomas Erben, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Henk Hoekstra, Shahab Joudaki, Mohammadjavad Vakili

Abstract: We present a new method for the mitigation of observational systematic effects in angular galaxy clustering via corrective random galaxy catalogues. Real and synthetic galaxy data, from the Kilo Degree Survey's (KiDS) 4$^{\rm{th}}$ Data Release (KiDS-$1000$) and the Full-sky Lognormal Astro-fields Simulation Kit (FLASK) package respectively, are used to train self-organising maps (SOMs) to learn t… ▽ More We present a new method for the mitigation of observational systematic effects in angular galaxy clustering via corrective random galaxy catalogues. Real and synthetic galaxy data, from the Kilo Degree Survey's (KiDS) 4$^{\rm{th}}$ Data Release (KiDS-$1000$) and the Full-sky Lognormal Astro-fields Simulation Kit (FLASK) package respectively, are used to train self-organising maps (SOMs) to learn the multivariate relationships between observed galaxy number density and up to six systematic-tracer variables, including seeing, Galactic dust extinction, and Galactic stellar density. We then create `organised' randoms, i.e. random galaxy catalogues with spatially variable number densities, mimicking the learnt systematic density modes in the data. Using realistically biased mock data, we show that these organised randoms consistently subtract spurious density modes from the two-point angular correlation function $w(\vartheta)$, correcting biases of up to $12σ$ in the mean clustering amplitude to as low as $0.1σ$, over a high signal-to-noise angular range of 7-100 arcmin. Their performance is also validated for angular clustering cross-correlations in a bright, flux-limited subset of KiDS-$1000$, comparing against an analogous sample constructed from highly-complete spectroscopic redshift data. Each organised random catalogue object is a `clone' carrying the properties of a real galaxy, and is distributed throughout the survey footprint according to the parent galaxy's position in systematics-space. Thus, sub-sample randoms are readily derived from a single master random catalogue via the same selection as applied to the real galaxies. Our method is expected to improve in performance with increased survey area, galaxy number density, and systematic contamination, making organised randoms extremely promising for current and future clustering analyses of faint samples. △ Less

Submitted 21 April, 2021; v1 submitted 15 December, 2020; originally announced December 2020.

Comments: 19 pages (7 appendix pages), 12 figures (8 appendix figures), accepted to A&A

Journal ref: A&A 648, A98 (2021)

Self-calibration and robust propagation of photometric redshift distribution uncertainties in weak gravitational lensing

Authors: B. Stölzner, B. Joachimi, A. Korn, H. Hildebrandt, A. H. Wright

Abstract: We present a method that accurately propagates residual uncertainties in photometric redshift distributions into the cosmological inference from weak lensing measurements. The redshift distributions of tomographic redshift bins are parameterised using a flexible modified Gaussian mixture model. We fit this model to pre-calibrated redshift distributions and implement an analytic marginalisation ove… ▽ More We present a method that accurately propagates residual uncertainties in photometric redshift distributions into the cosmological inference from weak lensing measurements. The redshift distributions of tomographic redshift bins are parameterised using a flexible modified Gaussian mixture model. We fit this model to pre-calibrated redshift distributions and implement an analytic marginalisation over the potentially several hundred redshift nuisance parameters in the weak lensing likelihood, which is demonstrated to accurately recover the cosmological posterior. By iteratively fitting cosmological and nuisance parameters arising from the redshift distribution model, we perform a self-calibration of the redshift distributions via the tomographic cosmic shear measurements. Our method is applied to the third data release of the Kilo-Degree Survey combined with the VISTA Kilo-Degree Infrared Galaxy Survey (KV450). We find constraints on cosmological parameters that are in very good agreement with the fiducial KV450 cosmic shear analysis and investigate the effects of the more flexible model on the self-calibrated redshift distributions. We observe posterior shifts of the medians of the five tomographic redshift distributions of up to $Δz \approx 0.02$, which are however degenerate with an observed decrease of the amplitude of intrinsic galaxy alignments by about $10\%$. △ Less

Submitted 7 May, 2021; v1 submitted 14 December, 2020; originally announced December 2020.

Comments: 16 pages, 10 figures, 4 appendices, version accepted by A&A

Journal ref: A&A 650, A148 (2021)