Showing 1–50 of 223 results for author: Heymans, C

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-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

The halo model for cosmology: a pedagogical review

Authors: Marika Asgari, Alexander J. Mead, Catherine Heymans

Abstract: We present a pedagogical review of the halo model, a flexible framework that can describe the distribution of matter and its tracers on non-linear scales for both conventional and exotic cosmological models. We start with the premise that the complex structure of the cosmic web can be described by the sum of its individual components: dark matter, gas, and galaxies, all distributed within spherica… ▽ More We present a pedagogical review of the halo model, a flexible framework that can describe the distribution of matter and its tracers on non-linear scales for both conventional and exotic cosmological models. We start with the premise that the complex structure of the cosmic web can be described by the sum of its individual components: dark matter, gas, and galaxies, all distributed within spherical haloes with a range of masses. The halo properties are specified through a series of simulation-calibrated ingredients including the halo mass function, non-linear halo bias and a dark matter density profile that can additionally account for the impact of baryon feedback. By incorporating a model of the galaxy halo occupation distribution, the properties of central and satellite galaxies, their non-linear bias and intrinsic alignment can be predicted. Through analytical calculations of spherical collapse in exotic cosmologies, the halo model also provides predictions for non-linear clustering in beyond-$Λ$CDM models. The halo model has been widely used to model observations of a variety of large-scale structure probes, most notably as the primary technique to model the underlying non-linear matter power spectrum. By documenting these varied and often distinct use cases, we seek to further coherent halo model analyses of future multi-tracer observables. This review is accompanied by the release of pyhalomodel: https://github.com/alexander-mead/pyhalomodel , flexible software to conduct a wide range of halo-model calculations. △ Less

Submitted 6 November, 2023; v1 submitted 15 March, 2023; originally announced March 2023.

Comments: Comments from the community are gratefully received for version 1

KiDS-1000: cross-correlation with Planck cosmic microwave background lensing and intrinsic alignment removal with self-calibration

Authors: Ji Yao, Huanyuan Shan, Pengjie Zhang, Xiangkun Liu, Catherine Heymans, Benjamin Joachimi, Marika Asgari, Maciej Bilicki, Hendrik Hildebrandt, Konrad Kuijken, Tilman Tröster, Jan Luca van den Busch, Angus Wright, Ziang Yan

Abstract: Galaxy shear - cosmic microwave background (CMB) lensing convergence cross-correlations contain additional information on cosmology to auto-correlations. While being immune to certain systematic effects, they are affected by the galaxy intrinsic alignments (IA). This may be responsible for the reported low lensing amplitude of the galaxy shear $\times$ CMB convergence cross-correlations, compared… ▽ More Galaxy shear - cosmic microwave background (CMB) lensing convergence cross-correlations contain additional information on cosmology to auto-correlations. While being immune to certain systematic effects, they are affected by the galaxy intrinsic alignments (IA). This may be responsible for the reported low lensing amplitude of the galaxy shear $\times$ CMB convergence cross-correlations, compared to the standard Planck $Λ$CDM (cosmological constant and cold dark matter) cosmology prediction. In this work, we investigate how IA affects the Kilo-Degree Survey (KiDS) galaxy lensing shear - Planck CMB lensing convergence cross-correlation and compare it to previous treatments with or without IA taken into consideration. More specifically, we compare marginalization over IA parameters and the IA self-calibration (SC) method (with additional observables defined only from the source galaxies) and prove that SC can efficiently break the degeneracy between the CMB lensing amplitude $A_{\rm lens}$ and the IA amplitude $A_{\rm IA}$. We further investigate how different systematics affect the resulting $A_{\rm IA}$ and $A_{\rm lens}$, and validate our results with the MICE2 simulation. We find that by including the SC method to constrain IA, the information loss due to the degeneracy between CMB lensing and IA is strongly reduced. The best-fit values are $A_{\rm lens}=0.84^{+0.22}_{-0.22}$ and $A_{\rm IA}=0.60^{+1.03}_{-1.03}$, while different angular scale cuts can affect $A_{\rm lens}$ by $\sim10\%$. We show that appropriate treatment of the boost factor, cosmic magnification, and photometric redshift modeling is important for obtaining the correct IA and cosmological results. △ Less

Submitted 10 April, 2023; v1 submitted 31 January, 2023; originally announced January 2023.

Comments: match version accepted by A&A

Journal ref: A&A 673, A111 (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)

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)

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)

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 halo model with beyond-linear halo bias: unbiasing cosmological constraints from galaxy-galaxy lensing and clustering

Authors: Constance Mahony, Andrej Dvornik, Alexander Mead, Catherine Heymans, Marika Asgari, Hendrik Hildebrandt, Hironao Miyatake, Takahiro Nishimichi, Robert Reischke

Abstract: We determine the error introduced in a joint halo model analysis of galaxy-galaxy lensing and galaxy clustering observables when adopting the standard approximation of linear halo bias. Considering the Kilo-Degree Survey, we forecast that ignoring the non-linear halo bias would result in up to 5$σ$ offsets in the recovered cosmological parameters describing structure growth, $S_8$, and the matter… ▽ More We determine the error introduced in a joint halo model analysis of galaxy-galaxy lensing and galaxy clustering observables when adopting the standard approximation of linear halo bias. Considering the Kilo-Degree Survey, we forecast that ignoring the non-linear halo bias would result in up to 5$σ$ offsets in the recovered cosmological parameters describing structure growth, $S_8$, and the matter density parameter, $Ω_{\mathrm{m}}$. We include the scales $10^{-1.3}<r_{\rm{p}} \ / h^{-1}\, \mathrm{Mpc}<10$ in the data vector, and the direction of these offsets are shown to depend on the freedom afforded to the halo model through other nuisance parameters. We conclude that a beyond-linear halo bias correction must therefore be included in future cosmological halo model analyses of large-scale structure observables on non-linear scales. △ Less

Submitted 30 June, 2022; v1 submitted 3 February, 2022; originally announced February 2022.

Comments: v2: Version accepted in MNRAS. Minor clarifications and added robustness test, conclusions unchanged

Report number: YITP-22-07

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

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

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)

The matter density PDF for modified gravity and dark energy with Large Deviations Theory

Authors: Matteo Cataneo, Cora Uhlemann, Christian Arnold, Alex Gough, Baojiu Li, Catherine Heymans

Abstract: We present an analytical description of the probability distribution function (PDF) of the smoothed three-dimensional matter density field for modified gravity and dark energy. Our approach, based on the principles of Large Deviations Theory, is applicable to general extensions of the standard $Λ$CDM cosmology. We show that late-time changes to the law of gravity and background expansion can be in… ▽ More We present an analytical description of the probability distribution function (PDF) of the smoothed three-dimensional matter density field for modified gravity and dark energy. Our approach, based on the principles of Large Deviations Theory, is applicable to general extensions of the standard $Λ$CDM cosmology. We show that late-time changes to the law of gravity and background expansion can be included through Einstein-de Sitter spherical collapse dynamics combined with linear theory calculations and a calibration measurement of the non-linear variance of the smoothed density field from a simple numerical simulation. In a comparison to $N$-body simulations for $f(R)$, DGP and evolving dark energy theories, we find percent level accuracy around the peak of the distribution for predictions in the mildly non-linear regime. A Fisher forecast of an idealised experiment with a Euclid-like survey volume demonstrates the power of combining measurements of the 3D matter PDF with the 3D matter power spectrum. This combination is shown to halve the uncertainty on parameters for an evolving dark energy model, relative to a power spectrum analysis on its own. The PDF is also found to substantially increase the detection significance for small departures from General Relativity, with improvements of up to six times compared to the power spectrum alone. This analysis is therefore very promising for future studies including non-Gaussian statistics, as it has the potential to alleviate the reliance of these analyses on expensive high resolution simulations and emulators. △ Less

Submitted 4 April, 2022; v1 submitted 6 September, 2021; originally announced September 2021.

Comments: 19 pages, 17 figures, 4 tables. Matches MNRAS accepted version. Our code to compute the matter PDF is available at https://github.com/mcataneo/pyLDT-cosmo

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 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

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)

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

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)

Cosmic Shear Cosmology Beyond 2-Point Statistics: A Combined Peak Count and Correlation Function Analysis of DES-Y1

Authors: Joachim Harnois-Déraps, Nicolas Martinet, Tiago Castro, Klaus Dolag, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Qianli Xia

Abstract: We constrain cosmological parameters from a joint cosmic shear analysis of peak-counts and the two-point shear correlation functions, as measured from the Dark Energy Survey (DES-Y1). We find the structure growth parameter $S_8\equiv σ_8\sqrt{Ω_{\rm m}/0.3} = 0.766^{+0.033}_{-0.038}$, which at 4.8% precision, provides one of the tightest constraints on $S_8$ from the DES-Y1 weak lensing data. In o… ▽ More We constrain cosmological parameters from a joint cosmic shear analysis of peak-counts and the two-point shear correlation functions, as measured from the Dark Energy Survey (DES-Y1). We find the structure growth parameter $S_8\equiv σ_8\sqrt{Ω_{\rm m}/0.3} = 0.766^{+0.033}_{-0.038}$, which at 4.8% precision, provides one of the tightest constraints on $S_8$ from the DES-Y1 weak lensing data. In our simulation-based method we determine the expected DES-Y1 peak-count signal for a range of cosmologies sampled in four $w$CDM parameters ($Ω_{\rm m}$, $σ_8$, $h$, $w_0$). We also determine the joint covariance matrix with over 1000 realisations at our fiducial cosmology. With mock DES-Y1 data we calibrate the impact of photometric redshift and shear calibration uncertainty on the peak-count, marginalising over these uncertainties in our cosmological analysis. Using dedicated training samples we show that our measurements are unaffected by mass resolution limits in the simulation, and that our constraints are robust against uncertainty in the effect of baryon feedback. Accurate modelling for the impact of intrinsic alignments on the tomographic peak-count remains a challenge, currently limiting our exploitation of cross-correlated peak counts between high and low redshift bins. We demonstrate that once calibrated, a fully tomographic joint peak-count and correlation functions analysis has the potential to reach a 3% precision on $S_8$ for DES-Y1. Our methodology can be adopted to model any statistic that is sensitive to the non-Gaussian information encoded in the shear field. In order to accelerate the development of these beyond-two-point cosmic shear studies, our simulations are made available to the community upon request. △ Less

Submitted 3 June, 2021; v1 submitted 4 December, 2020; originally announced December 2020.

Comments: 27 pages, 23 figures, updated to MNRAS accepted version

Strong detection of the CMB lensingxgalaxy weak lensingcross-correlation from ACT-DR4,PlanckLegacy and KiDS-1000

Authors: Naomi Clare Robertson, David Alonso, Joachim Harnois-Déraps, Omar Darwish, Arun Kannawad, Alexandra Amon, Marika Asgari, Maciej Bilicki, Erminia Calabrese, Steve K. Choi, Mark J. Devlin, Jo Dunkley, Andrej Dvornik, Thomas Erben, Simone Ferraro, Maria Cristina Fortuna, Benjamin Giblin, Dongwon Han, Catherine Heymans, Hendrik Hildebrandt, J. Colin Hill, Matt Hilton, Shuay-Pwu P. Ho, Henk Hoekstra, Johannes Hubmayr , et al. (26 additional authors not shown)

Abstract: We measure the cross-correlation between galaxy weak lensing data from the Kilo Degree Survey (KiDS-1000, DR4) and cosmic microwave background (CMB) lensing data from the Atacama Cosmology Telescope (ACT, DR4) and the Planck Legacy survey. We use two samples of source galaxies, selected with photometric redshifts, $(0.1<z_{\rm B}<1.2)$ and $(1.2<z_{\rm B}<2)$, which produce a combined detection si… ▽ More We measure the cross-correlation between galaxy weak lensing data from the Kilo Degree Survey (KiDS-1000, DR4) and cosmic microwave background (CMB) lensing data from the Atacama Cosmology Telescope (ACT, DR4) and the Planck Legacy survey. We use two samples of source galaxies, selected with photometric redshifts, $(0.1<z_{\rm B}<1.2)$ and $(1.2<z_{\rm B}<2)$, which produce a combined detection significance of the CMB lensing/weak galaxy lensing cross-spectrum of $7.7σ$. With the lower redshift galaxy sample, for which the cross-correlation is detected at a significance of $5.3σ$, we present joint cosmological constraints on the matter density parameter, $Ω_{\rm m}$, and the matter fluctuation amplitude parameter, $σ_8$, marginalising over three nuisance parameters that model our uncertainty in the redshift and shear calibration, and the intrinsic alignment of galaxies. We find our measurement to be consistent with the best-fitting flat $Λ$CDM cosmological models from both Planck and KiDS-1000. We demonstrate the capacity of CMB-weak lensing cross-correlations to set constraints on either the redshift or shear calibration, by analysing a previously unused high-redshift KiDS galaxy sample $(1.2<z_{\rm B}<2)$, with the cross-correlation detected at a significance of $7σ$. This analysis provides an independent assessment for the accuracy of redshift measurements in a regime that is challenging to calibrate directly owing to known incompleteness in spectroscopic surveys. △ Less

Submitted 23 November, 2020; originally announced November 2020.

Comments: 13 pages, 9 figures, 1 tables, submitted to A&A

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

Discovery of two Einstein crosses from massive post--blue nugget galaxies at z>1 in KiDS

Authors: N. R. Napolitano, R. Li, C. Spiniello, C. Tortora, A. Sergeyev, G. D'Ago, X. Guo, L. Xie, M. Radovich, N. Roy, L. V. E. Koopmans, K. Kuijken, M. Bilicki, T. Erben, F. Getman, C. Heymans, H. Hildebrandt, C. Moya, H. Y. Shan, G. Vernardos, A. H. Wright

Abstract: We report the discovery of two Einstein Crosses (ECs) in the footprint of the Kilo-Degree Survey (KiDS): KIDS J232940-340922 and KIDS J122456+005048. Using integral field spectroscopy from MUSE@VLT, we confirm their gravitational-lens nature. In both cases, the four spectra of the source clearly show a prominence of absorption features, hence revealing an evolved stellar population with little sta… ▽ More We report the discovery of two Einstein Crosses (ECs) in the footprint of the Kilo-Degree Survey (KiDS): KIDS J232940-340922 and KIDS J122456+005048. Using integral field spectroscopy from MUSE@VLT, we confirm their gravitational-lens nature. In both cases, the four spectra of the source clearly show a prominence of absorption features, hence revealing an evolved stellar population with little star formation. The lensing model of the two systems, assuming a singular isothermal ellipsoid (SIE) with external shear, shows that: 1) the two crosses, located at redshift $z=0.38$ and 0.24, have Einstein radius $R_{\rm E}=5.2$ kpc and 5.4 kpc, respectively; 2) their projected dark matter fractions inside the half effective radius are 0.60 and 0.56 (Chabrier IMF); 3) the sources are ultra-compact galaxies, $R_{\rm e}\sim0.9$ kpc (at redshift $z_{\rm s}=1.59$) and $R_{\rm e}\sim0.5$ kpc ($z_{\rm s}=1.10$), respectively. These results are unaffected by the underlying mass density assumption. Due to size, blue color and absorption-dominated spectra, corroborated by low specific star-formation rates derived from optical-NIR spectral energy distribution fitting, we argue that the two lensed sources in these ECs are blue nuggets migrating toward their quenching phase. △ Less

Submitted 19 November, 2020; v1 submitted 18 November, 2020; originally announced November 2020.

Comments: Accepted for publication on APJL

KiDS-1000 Cosmology: constraints beyond flat $Λ$CDM

Authors: Tilman Tröster, Marika Asgari, Chris Blake, Matteo Cataneo, Catherine Heymans, Hendrik Hildebrandt, Benjamin Joachimi, Chieh-An Lin, Ariel G. Sánchez, Angus H. Wright, Maciej Bilicki, Benjamin Bose, Martin Crocce, Andrej Dvornik, Thomas Erben, Benjamin Giblin, Karl Glazebrook, Henk Hoekstra, Shahab Joudaki, Arun Kannawadi, Fabian Köhlinger, Konrad Kuijken, Chris Lidman, Lucas Lombriser, Alexander Mead , et al. (4 additional authors not shown)

Abstract: We present constraints on extensions to the flat $Λ$CDM cosmological model by varying the spatial curvature $Ω_K$, the sum of the neutrino masses $\sum m_ν$, the dark energy equation of state parameter $w$, and the Hu-Sawicki $f(R)$ gravity $f_{R0}$ parameter. With the combined $3\times2$pt measurements of cosmic shear from the Kilo-Degree Survey (KiDS-1000), galaxy clustering from the Baryon Osci… ▽ More We present constraints on extensions to the flat $Λ$CDM cosmological model by varying the spatial curvature $Ω_K$, the sum of the neutrino masses $\sum m_ν$, the dark energy equation of state parameter $w$, and the Hu-Sawicki $f(R)$ gravity $f_{R0}$ parameter. With the combined $3\times2$pt measurements of cosmic shear from the Kilo-Degree Survey (KiDS-1000), galaxy clustering from the Baryon Oscillation Spectroscopic Survey (BOSS), and galaxy-galaxy lensing from the overlap between KiDS-1000, BOSS, and the spectroscopic 2-degree Field Lensing Survey (2dFLenS), we find results that are fully consistent with a flat $Λ$CDM model with $Ω_K=0.011^{+0.054}_{-0.057}$, $\sum m_ν<1.76$ eV (95% CL), and $w=-0.99^{+0.11}_{-0.13}$. The $f_{R0}$ parameter is unconstrained in our fully non-linear $f(R)$ cosmic shear analysis. Considering three different model selection criteria, we find no clear preference for either the fiducial flat $Λ$CDM model or any of the considered extensions. Besides extensions to the flat $Λ$CDM parameter space, we also explore restrictions to common subsets of the flat $Λ$CDM parameter space by fixing the amplitude of the primordial power spectrum to the Planck best-fit value, as well as adding external data from supernovae and lensing of the CMB. Neither the beyond-$Λ$CDM models nor the imposed restrictions explored in this analysis are able to resolve the $\sim 3σ$ tension in $S_8$ between the $3\times2$pt constraints and Planck, with the exception of $w$CDM, where the $S_8$ tension is resolved. The tension in the $w$CDM case persists, however, when considering the joint $S_8$-$w$ parameter space. The joint flat $Λ$CDM CMB lensing and $3\times2$pt analysis is found to yield tight constraints on $Ω_{\rm m}=0.307^{+0.008}_{-0.013}$, $σ_8=0.769^{+0.022}_{-0.010}$, and $S_8=0.779^{+0.013}_{-0.013}$. △ Less

Submitted 2 March, 2021; v1 submitted 30 October, 2020; originally announced October 2020.

Comments: 15 pages, 9 figures, 3 tables, accepted in A&A. This paper concludes the KiDS-1000 series of papers: Heymans, Tröster et al. (arXiv:2007.15632), Asgari et al. (arXiv:2007.15633), Hildebrandt et al. (arXiv:2007.15635), Joachimi et al. (arXiv:2007.01844), and Giblin et al. (arXiv:2007.01845). Data products, likelihoods, and posteriors can be found http://kids.strw.leidenuniv.nl/DR4/lensing.php

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

Photometric selection and redshifts for quasars in the Kilo-Degree Survey Data Release 4

Authors: S. J. Nakoneczny, M. Bilicki, A. Pollo, M. Asgari, A. Dvornik, T. Erben, B. Giblin, C. Heymans, H. Hildebrandt, A. Kannawadi, K. Kuijken, N. R. Napolitano, E. Valentijn

Abstract: We present a catalog of quasars and corresponding redshifts in the Kilo-Degree Survey (KiDS) Data Release 4. We trained machine learning (ML) models, using optical ugri and near-infrared ZYJHK_s bands, on objects known from Sloan Digital Sky Survey (SDSS) spectroscopy. We define inference subsets from the 45 million objects of the KiDS photometric data limited to 9-band detections. We show that pr… ▽ More We present a catalog of quasars and corresponding redshifts in the Kilo-Degree Survey (KiDS) Data Release 4. We trained machine learning (ML) models, using optical ugri and near-infrared ZYJHK_s bands, on objects known from Sloan Digital Sky Survey (SDSS) spectroscopy. We define inference subsets from the 45 million objects of the KiDS photometric data limited to 9-band detections. We show that projections of the high-dimensional feature space can be successfully used to investigate the estimations. The model creation employs two test subsets: randomly selected and the faintest objects, which allows to fit the bias versus variance trade-off. We tested three ML models: random forest (RF), XGBoost (XGB), and artificial neural network (ANN). We find that XGB is the most robust model for classification, while ANN performs the best for combined classification and redshift. The inference results are tested using number counts, Gaia parallaxes, and other quasar catalogs. Based on these tests, we derived the minimum classification probability which provides the best purity versus completeness trade-off: p(QSO_cand) > 0.9 for r < 22 and p(QSO_cand) > 0.98 for 22 < r < 23.5. We find 158,000 quasar candidates in the safe inference subset (r < 22) and an additional 185,000 candidates in the reliable extrapolation regime (22 < r < 23.5). Test-data purity equals 97% and completeness is 94%; the latter drops by 3% in the extrapolation to data fainter by one magnitude than the training set. The photometric redshifts were modeled with Gaussian uncertainties. The redshift error (mean and scatter) equals 0.01 +/- 0.1 in the safe subset and -0.0004 +/- 0.2 in the extrapolation, in a redshift range of 0.14 < z < 3.63. Our success of the extrapolation challenges the way that models are optimized and applied at the faint data end. The catalog is ready for cosmology and active galactic nucleus (AGN) studies. △ Less

Submitted 9 March, 2021; v1 submitted 26 October, 2020; originally announced October 2020.

Comments: We publicly release the catalog at kids.strw.leidenuniv.nl/DR4/quasarcatalog.php , and the code at github.com/snakoneczny/kids-quasars

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

Tightening weak lensing constraints on the ellipticity of galaxy-scale dark matter haloes

Authors: Tim Schrabback, Henk Hoekstra, Ludovic Van Waerbeke, Edo van Uitert, Christos Georgiou, Marika Asgari, Patrick Côté, Jean-Charles Cuillandre, Thomas Erben, Laura Ferrarese, Stephen D. J. Gwyn, Catherine Heymans, Hendrik Hildebrandt, Arun Kannawadi, Konrad Kuijken, Alexie Leauthaud, Martin Makler, Simona Mei, Lance Miller, Anand Raichoor, Peter Schneider, Angus Wright

Abstract: Cosmological simulations predict that galaxies are embedded into triaxial dark matter haloes, which appear approximately elliptical in projection. Weak gravitational lensing allows us to constrain these halo shapes and thereby test the nature of dark matter. Weak lensing has already provided robust detections of the signature of halo flattening at the mass scales of groups and clusters, whereas re… ▽ More Cosmological simulations predict that galaxies are embedded into triaxial dark matter haloes, which appear approximately elliptical in projection. Weak gravitational lensing allows us to constrain these halo shapes and thereby test the nature of dark matter. Weak lensing has already provided robust detections of the signature of halo flattening at the mass scales of groups and clusters, whereas results for galaxies have been somewhat inconclusive. Here we combine data from five surveys (NGVSLenS, KiDS/KV450, CFHTLenS, CS82, and RCSLenS) in order to tighten observational constraints on galaxy-scale halo ellipticity for photometrically selected lens samples. We constrain $f_\rm{h}$, the average ratio between the aligned component of the halo ellipticity and the ellipticity of the light distribution, finding $f_\rm{h}=0.303^{+0.080}_{-0.079}$ for red lenses and $f_\rm{h}=0.217^{+0.160}_{-0.159}$ for blue lenses when assuming elliptical NFW density profiles and a linear scaling between halo ellipticity and galaxy ellipticity. Our constraints for red galaxies constitute the currently most significant ($3.8σ$) systematics-corrected detection of the signature of halo flattening at the mass scale of galaxies. Our results are in good agreement with expectations from the Millennium Simulation that apply the same analysis scheme and incorporate models for galaxy-halo misalignment. Assuming these misalignment models and the analysis assumptions stated above are correct, our measurements imply an average dark matter halo ellipticity for the studied red galaxy samples of $\langle|ε_\rm{h}|\rangle=0.174\pm 0.046$, where $|ε_{h}|=(1-q)/(1+q)$ relates to the ratio $q=b/a$ of the minor and major axes of the projected mass distribution. Similar measurements based on larger upcoming weak lensing data sets can help to calibrate models for intrinsic galaxy alignments. [abridged] △ Less

Submitted 20 November, 2020; v1 submitted 1 October, 2020; originally announced October 2020.

Comments: 23 pages, 11 figures, 2 tables. Accepted for publication in A&A. In this version corrections from the A&A language editor have been applied. Comments welcome

Journal ref: A&A 646, A73 (2021)

HMcode-2020: Improved modelling of non-linear cosmological power spectra with baryonic feedback

Authors: Alexander Mead, Samuel Brieden, Tilman Tröster, Catherine Heymans

Abstract: We present an updated version of the HMcode augmented halo model that can be used to make accurate predictions of the non-linear matter power spectrum over a wide range of cosmologies. Major improvements include modelling of BAO damping in the power spectrum and an updated treatment of massive neutrinos. We fit our model to simulated power spectra and show that we can match the results with an RMS… ▽ More We present an updated version of the HMcode augmented halo model that can be used to make accurate predictions of the non-linear matter power spectrum over a wide range of cosmologies. Major improvements include modelling of BAO damping in the power spectrum and an updated treatment of massive neutrinos. We fit our model to simulated power spectra and show that we can match the results with an RMS error of 2.5 per cent across a range of cosmologies, scales $k < 10\,h\mathrm{Mpc}^{-1}$, and redshifts $z<2$. The error rarely exceeds 5 per cent and never exceeds 16 per cent. The worst-case errors occur at $z\simeq2$, or for cosmologies with unusual dark-energy equations of state. This represents a significant improvement over previous versions of HMcode, and over other popular fitting functions, particularly for massive-neutrino cosmologies with high neutrino mass. We also present a simple halo model that can be used to model the impact of baryonic feedback on the power spectrum. This six-parameter physical model includes gas expulsion by AGN feedback and encapsulates star formation. By comparing this model to data from hydrodynamical simulations we demonstrate that the power spectrum response to feedback is matched at the $<1$ per cent level for $z<1$ and $k<20\,h\mathrm{Mpc}^{-1}$. We also present a single-parameter variant of this model, parametrized in terms of feedback strength, which is only slightly less accurate. We make code available for our non-linear and baryon models at https://github.com/alexander-mead/HMcode and it is also available within CAMB and soon within CLASS. △ Less

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

Comments: 17 pages, 5 figures, 4 appendices; v2 - matches accepted version, new appendix with comparisons between HMcode and 6 different emulators

Clustering of red-sequence galaxies in the fourth data release ofthe Kilo-Degree Survey

Authors: Mohammadjavad Vakili, Henk Hoekstra, Maciej Bilicki, Maria-Cristina Fortuna, Konrad Kuijken, Angus H. Wright, Marika Asgari, Michael Brown, Elisabeth Dombrovskij, Thomas Erben, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Harry Johnston, Shahab Joudaki, Arun Kannawadi

Abstract: We present a sample of luminous red-sequence galaxies to study the large-scale structure in the fourth data release of the Kilo-Degree Survey. The selected galaxies are defined by a red-sequence template, in the form of a data-driven model of the colour-magnitude relation conditioned on redshift. In this work, the red-sequence template is built using the broad-band optical+near infrared photometry… ▽ More We present a sample of luminous red-sequence galaxies to study the large-scale structure in the fourth data release of the Kilo-Degree Survey. The selected galaxies are defined by a red-sequence template, in the form of a data-driven model of the colour-magnitude relation conditioned on redshift. In this work, the red-sequence template is built using the broad-band optical+near infrared photometry of KiDS-VIKING and the overlapping spectroscopic data sets. The selection process involves estimating the red-sequence redshifts, assessing the purity of the sample, and estimating the underlying redshift distributions of redshift bins. After performing the selection, we mitigate the impact of survey properties on the observed number density of galaxies by assigning photometric weights to the galaxies. We measure the angular two-point correlation function of the red galaxies in four redshift bins, and constrain the large scale bias of our red-sequence sample assuming a fixed $Λ$CDM cosmology. We find consistent linear biases for two luminosity-threshold samples (dense and luminous). We find that our constraints are well characterized by the passive evolution model. △ Less

Submitted 30 August, 2020; originally announced August 2020.

Comments: submitted to A&A

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

KiDS-1000 catalogue: Redshift distributions and their calibration

Authors: H. Hildebrandt, J. L. van den Busch, A. H. Wright, C. Blake, B. Joachimi, K. Kuijken, T. Tröster, M. Asgari, M. Bilicki, J. T. A. de Jong, A. Dvornik, T. Erben, F. Getman, B. Giblin, C. Heymans, A. Kannawadi, C. -A. Lin, H. -Y. Shan

Abstract: We present redshift distribution estimates of galaxies selected from the fourth data release of the Kilo-Degree Survey over an area of $\sim1000$ deg$^2$ (KiDS-1000). These redshift distributions represent one of the crucial ingredients for weak gravitational lensing measurements with the KiDS-1000 data. The primary estimate is based on deep spectroscopic reference catalogues that are re-weighted… ▽ More We present redshift distribution estimates of galaxies selected from the fourth data release of the Kilo-Degree Survey over an area of $\sim1000$ deg$^2$ (KiDS-1000). These redshift distributions represent one of the crucial ingredients for weak gravitational lensing measurements with the KiDS-1000 data. The primary estimate is based on deep spectroscopic reference catalogues that are re-weighted with the help of a self-organising map (SOM) to closely resemble the KiDS-1000 sources, split into five tomographic redshift bins in the photometric redshift range $0.1<z_\mathrm{B}\le1.2$. Sources are selected such that they only occupy that volume of nine-dimensional magnitude-space that is also covered by the reference samples (`gold' selection). Residual biases in the mean redshifts determined from this calibration are estimated from mock catalogues to be $\lesssim0.01$ for all five bins with uncertainties of $\sim 0.01$. This primary SOM estimate of the KiDS-1000 redshift distributions is complemented with an independent clustering redshift approach. After validation of the clustering-$z$ on the same mock catalogues and a careful assessment of systematic errors, we find no significant bias of the SOM redshift distributions with respect to the clustering-$z$ measurements. The SOM redshift distributions re-calibrated by the clustering-$z$ represent an alternative calibration of the redshift distributions with only slightly larger uncertainties in the mean redshifts of $\sim 0.01-0.02$ to be used in KiDS-1000 cosmological weak lensing analyses. As this includes the SOM uncertainty, clustering-$z$ are shown to be fully competitive on KiDS-1000 data. △ Less

Submitted 28 January, 2021; v1 submitted 30 July, 2020; originally announced July 2020.

Comments: 15 pages, 7 figures, 3 tables, accepted for publication in A&A. This paper is part of the KiDS-1000 series of papers, accompanying Heymans, Tröster et al. (arXiv:2007.15632), Asgari et al. (arXiv:2007.15633), Joachimi et al. (arXiv:2007.01844), and Giblin et al. (arXiv:2007.01845). Online KiDS-1000 talks and seminars can be viewed at http://kids.strw.leidenuniv.nl/KiDS-1000.php

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

KiDS-1000 Cosmology: Cosmic shear constraints and comparison between two point statistics

Authors: Marika Asgari, Chieh-An Lin, Benjamin Joachimi, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Arun Kannawadi, Benjamin Stölzner, Tilman Tröster, Jan Luca van den Busch, Angus H. Wright, Maciej Bilicki, Chris Blake, Jelte de Jong, Andrej Dvornik, Thomas Erben, Fedor Getman, Henk Hoekstra, Fabian Köhlinger, Konrad Kuijken, Lance Miller, Mario Radovich, Peter Schneider, HuanYuan Shan, Edwin Valentijn

Abstract: We present cosmological constraints from a cosmic shear analysis of the fourth data release of the Kilo-Degree Survey (KiDS-1000), doubling the survey area with nine-band optical and near-infrared photometry with respect to previous KiDS analyses. Adopting a spatially flat $Λ$CDM model, we find $S_8 = σ_8 (Ω_{\rm m}/0.3)^{0.5} = 0.759^{+0.024}_{-0.021}$ for our fiducial analysis, which is in $3σ$… ▽ More We present cosmological constraints from a cosmic shear analysis of the fourth data release of the Kilo-Degree Survey (KiDS-1000), doubling the survey area with nine-band optical and near-infrared photometry with respect to previous KiDS analyses. Adopting a spatially flat $Λ$CDM model, we find $S_8 = σ_8 (Ω_{\rm m}/0.3)^{0.5} = 0.759^{+0.024}_{-0.021}$ for our fiducial analysis, which is in $3σ$ tension with the prediction of the Planck Legacy analysis of the cosmic microwave background. We compare our fiducial COSEBIs (Complete Orthogonal Sets of E/B-Integrals) analysis with complementary analyses of the two-point shear correlation function and band power spectra, finding results to be in excellent agreement. We investigate the sensitivity of all three statistics to a number of measurement, astrophysical, and modelling systematics, finding our $S_8$ constraints to be robust and dominated by statistical errors. Our cosmological analysis of different divisions of the data pass the Bayesian internal consistency tests, with the exception of the second tomographic bin. As this bin encompasses low redshift galaxies, carrying insignificant levels of cosmological information, we find that our results are unchanged by the inclusion or exclusion of this sample. △ Less

Submitted 13 October, 2020; v1 submitted 30 July, 2020; originally announced July 2020.

Comments: 33 pages, 14 figures, 7 tables, Accepted to be published in A&A. This paper is part of the KiDS-1000 series of papers: Heymans, Tröster et al. (arXiv:2007.15632), Hildebrandt et al. (arXiv:2007.15635), Joachimi et al. (arXiv:2007.01844) and Giblin et al. (arXiv:2007.01845). Online KiDS-1000 talks can be viewed at http://kids.strw.leidenuniv.nl/KiDS-1000.php or the KiDS Consortium youtube page

Journal ref: A&A 645, A104 (2021)

KiDS-1000 Cosmology: Multi-probe weak gravitational lensing and spectroscopic galaxy clustering constraints

Authors: Catherine Heymans, Tilman Tröster, Marika Asgari, Chris Blake, Hendrik Hildebrandt, Benjamin Joachimi, Konrad Kuijken, Chieh-An Lin, Ariel G. Sánchez, Jan Luca van den Busch, Angus H. Wright, Alexandra Amon, Maciej Bilicki, Jelte de Jong, Martin Crocce, Andrej Dvornik, Thomas Erben, Maria Cristina Fortuna, Fedor Getman, Benjamin Giblin, Karl Glazebrook, Henk Hoekstra, Shahab Joudaki, Arun Kannawadi, Fabian Köhlinger , et al. (9 additional authors not shown)

Abstract: We present a joint cosmological analysis of weak gravitational lensing observations from the Kilo-Degree Survey (KiDS-1000), with redshift-space galaxy clustering observations from the Baryon Oscillation Spectroscopic Survey (BOSS), and galaxy-galaxy lensing observations from the overlap between KiDS-1000, BOSS and the spectroscopic 2-degree Field Lensing Survey (2dFLenS). This combination of larg… ▽ More We present a joint cosmological analysis of weak gravitational lensing observations from the Kilo-Degree Survey (KiDS-1000), with redshift-space galaxy clustering observations from the Baryon Oscillation Spectroscopic Survey (BOSS), and galaxy-galaxy lensing observations from the overlap between KiDS-1000, BOSS and the spectroscopic 2-degree Field Lensing Survey (2dFLenS). This combination of large-scale structure probes breaks the degeneracies between cosmological parameters for individual observables, resulting in a constraint on the structure growth parameter $S_8=σ_8 \sqrt{Ω_{\rm m}/0.3} = 0.766^{+0.020}_{-0.014}$, that has the same overall precision as that reported by the full-sky cosmic microwave background observations from Planck. The recovered $S_8$ amplitude is low, however, by $8.3 \pm 2.6$ % relative to Planck. This result builds from a series of KiDS-1000 analyses where we validate our methodology with variable depth mock galaxy surveys, our lensing calibration with image simulations and null-tests, and our optical-to-near-infrared redshift calibration with multi-band mock catalogues and a spectroscopic-photometric clustering analysis. The systematic uncertainties identified by these analyses are folded through as nuisance parameters in our cosmological analysis. Inspecting the offset between the marginalised posterior distributions, we find that the $S_8$-difference with Planck is driven by a tension in the matter fluctuation amplitude parameter, $σ_8$. We quantify the level of agreement between the CMB and our large-scale structure constraints using a series of different metrics, finding differences with a significance ranging between $\sim\! 3\,σ$, when considering the offset in $S_{8}$, and $\sim\! 2\,σ$, when considering the full multi-dimensional parameter space. △ Less

Submitted 8 December, 2020; v1 submitted 30 July, 2020; originally announced July 2020.

Comments: 25 pages, 11 figures, 5 tables, A&A accepted, including a new appendix on Intrinsic Alignments. The KiDS-1000 data products are available for download at http://kids.strw.leidenuniv.nl/DR4/lensing.php. This data release includes open source software, the shear-photo-z catalogue, the cosmic shear and 3x2pt data vectors and covariances, and posteriors in the form of Multinest chains

Journal ref: A&A 646, A140 (2021)

Testing KiDS cross-correlation redshifts with simulations

Authors: J. L. van den Busch, H. Hildebrandt, A. H. Wright, C. B. Morrison, C. Blake, B. Joachimi, T. Erben, C. Heymans, K. Kuijken, E. N. Taylor

Abstract: Measuring cosmic shear in wide-field imaging surveys requires accurate knowledge of the redshift distribution of all sources. The clustering-redshift technique exploits the angular cross-correlation of a target galaxy sample with unknown redshifts and a reference sample with known redshifts, and is an attractive alternative to colour-based methods of redshift calibration. We test the performance o… ▽ More Measuring cosmic shear in wide-field imaging surveys requires accurate knowledge of the redshift distribution of all sources. The clustering-redshift technique exploits the angular cross-correlation of a target galaxy sample with unknown redshifts and a reference sample with known redshifts, and is an attractive alternative to colour-based methods of redshift calibration. We test the performance of such clustering redshift measurements using mock catalogues that resemble the Kilo-Degree Survey (KiDS). These mocks are created from the MICE simulation and closely mimic the properties of the KiDS source sample and the overlapping spectroscopic reference samples. We quantify the performance of the clustering redshifts by comparing the cross-correlation results with the true redshift distributions in each of the five KiDS photometric redshift bins. Such a comparison to an informative model is necessary due to the incompleteness of the reference samples at high redshifts. Clustering mean redshifts are unbiased at $|Δz|<0.006$ under these conditions. The redshift evolution of the galaxy bias can be reliably mitigated at this level of precision using auto-correlation measurements and self-consistency relations, and will not become a dominant source of systematic error until the arrival of Stage-IV cosmic shear surveys. Using redshift distributions from a direct colour-based estimate instead of the true redshift distributions as a model for comparison with the clustering redshifts increases the biases in the mean to up to $|Δz|\sim0.04$. This indicates that the interpretation of clustering redshifts in real-world applications will require more sophisticated (parameterised) models of the redshift distribution in the future. If such better models are available, the clustering-redshift technique promises to be a highly complementary alternative to other methods of redshift calibration. △ Less

Submitted 13 April, 2021; v1 submitted 3 July, 2020; originally announced July 2020.

Comments: 21 pages, 18 figures, 10 tables

Journal ref: A&A 642, A200 (2020)

KiDS-1000 catalogue: Weak gravitational lensing shear measurements

Authors: Benjamin Giblin, Catherine Heymans, Marika Asgari, Hendrik Hildebrandt, Henk Hoekstra, Benjamin Joachimi, Arun Kannawadi, Konrad Kuijken, Chieh-An Lin, Lance Miller, Tilman Tröster, Jan Luca van den Busch, Angus H. Wright, Maciej Bilicki, Chris Blake, Jelte de Jong, Andrej Dvornik, Thomas Erben, Fedor Getman, Nicola R. Napolitano, Peter Schneider, HuanYuan Shan

Abstract: We present weak lensing shear catalogues from the fourth data release of the Kilo-Degree Survey, KiDS-1000, spanning 1006 square degrees of deep and high-resolution imaging. Our `gold-sample' of galaxies, with well-calibrated photometric redshift distributions, consists of 21 million galaxies with an effective number density of $6.17$ galaxies per square arcminute. We quantify the accuracy of the… ▽ More We present weak lensing shear catalogues from the fourth data release of the Kilo-Degree Survey, KiDS-1000, spanning 1006 square degrees of deep and high-resolution imaging. Our `gold-sample' of galaxies, with well-calibrated photometric redshift distributions, consists of 21 million galaxies with an effective number density of $6.17$ galaxies per square arcminute. We quantify the accuracy of the spatial, temporal, and flux-dependent point-spread function (PSF) model, verifying that the model meets our requirements to induce less than a $0.1σ$ change in the inferred cosmic shear constraints on the clustering cosmological parameter $S_8 = σ_8\sqrt{Ω_{\rm m}/0.3}$. Through a series of two-point null-tests, we validate the shear estimates, finding no evidence for significant non-lensing B-mode distortions in the data. The PSF residuals are detected in the highest-redshift bins, originating from object selection and/or weight bias. The amplitude is, however, shown to be sufficiently low and within our stringent requirements. With a shear-ratio null-test, we verify the expected redshift scaling of the galaxy-galaxy lensing signal around luminous red galaxies. We conclude that the joint KiDS-1000 shear and photometric redshift calibration is sufficiently robust for combined-probe gravitational lensing and spectroscopic clustering analyses. △ Less

Submitted 8 December, 2020; v1 submitted 3 July, 2020; originally announced July 2020.

Comments: 24 pages, 11 figures, version accepted by A&A. The KiDS-1000 data products are available for download at http://kids.strw.leidenuniv.nl/DR4/lensing.php. This cosmology data release includes open source software, the shear-photo-z catalogue, the cosmic shear and 3x2pt data vectors and covariances, and posteriors in the form of Multinest chains

Journal ref: A&A 645, A105 (2021)

KiDS-1000 Methodology: Modelling and inference for joint weak gravitational lensing and spectroscopic galaxy clustering analysis

Authors: B. Joachimi, C. -A. Lin, M. Asgari, T. Tröster, C. Heymans, H. Hildebrandt, F. Köhlinger, A. G. Sánchez, A. H. Wright, M. Bilicki, C. Blake, J. L. van den Busch, M. Crocce, A. Dvornik, T. Erben, F. Getman, B. Giblin, H. Hoekstra, A. Kannawadi, K. Kuijken, N. R. Napolitano, P. Schneider, R. Scoccimarro, E. Sellentin, H. Y. Shan , et al. (2 additional authors not shown)

Abstract: We present the methodology for a joint cosmological analysis of weak gravitational lensing from the fourth data release of the ESO Kilo-Degree Survey (KiDS-1000) and galaxy clustering from the partially overlapping BOSS and 2dFLenS surveys. Cross-correlations between galaxy positions and ellipticities have been incorporated into the analysis, necessitating a hybrid model of non-linear scales that… ▽ More We present the methodology for a joint cosmological analysis of weak gravitational lensing from the fourth data release of the ESO Kilo-Degree Survey (KiDS-1000) and galaxy clustering from the partially overlapping BOSS and 2dFLenS surveys. Cross-correlations between galaxy positions and ellipticities have been incorporated into the analysis, necessitating a hybrid model of non-linear scales that blends perturbative and non-perturbative approaches, and an assessment of contributions by astrophysical effects. All weak lensing signals are measured consistently via Fourier-space statistics that are insensitive to the survey mask and display low levels of mode mixing. The calibration of photometric redshift distributions and multiplicative gravitational shear bias has been updated, and a more complete tally of residual calibration uncertainties is propagated into the likelihood. A dedicated suite of more than 20000 mocks is used to assess the performance of covariance models and to quantify the impact of survey geometry and spatial variations of survey depth on signals and their errors. The sampling distributions for the likelihood and the $χ^2$ goodness-of-fit statistic have been validated, with proposed changes to the number of degrees of freedom. Standard weak lensing point estimates on $S_8=σ_8\,(Ω_{\rm m}/0.3)^{1/2}$ derived from its marginal posterior are easily misinterpreted to be biased low, and an alternative estimator and associated credible interval have been proposed. Known systematic effects pertaining to weak lensing modelling and inference are shown to bias $S_8$ by no more than 0.1 standard deviations, with the caveat that no conclusive validation data exist for models of intrinsic galaxy alignments. Compared to the previous KiDS analyses, $S_8$ constraints are expected to improve by 20% for weak lensing alone and by 29% for the joint analysis. [abridged] △ Less

Submitted 18 December, 2020; v1 submitted 3 July, 2020; originally announced July 2020.

Comments: 45 pages, 34 figures, 7 tables; minor changes to match version accepted by A&A. The KiDS-1000 data products are available for download at http://kids.strw.leidenuniv.nl/DR4/lensing.php. This data release includes open source software, the shear-photo-z catalogue, the cosmic shear and 3x2pt data vectors and covariances, and posteriors in the form of Multinest chains

Journal ref: A&A 646, A129 (2021)

An adapted filter function for density split statistics in weak lensing

Authors: Pierre Burger, Peter Schneider, Vasiliy Demchenko, Joachim Harnois-Deraps, Catherine Heymans, Hendrik Hildebrandt, Sandra Unruh

Abstract: Context. The density split statistics in weak gravitational lensing analyses probes the correlation between regions of different (foreground) galaxy number densities and their weak lensing signal, measured by the shape distortion of background galaxies. Aims. In this paper, we reconsider density split statistics, by constructing a new angular filter function that is adapted to the expected relatio… ▽ More Context. The density split statistics in weak gravitational lensing analyses probes the correlation between regions of different (foreground) galaxy number densities and their weak lensing signal, measured by the shape distortion of background galaxies. Aims. In this paper, we reconsider density split statistics, by constructing a new angular filter function that is adapted to the expected relation between galaxy number density and shear pattern, in a way that the filter weighting the galaxy number density is matched to the filter that is used to quantify the shear signal. Methods. We use the results of numerical ray-tracing simulations, specifically through Millennium Simulation supplemented by a galaxy distribution based on a semi-analytic model, to construct a matched pair of adapted filter functions for the galaxy density and the tangential shear signal. We compare the performance of our new filter to the previously used top-hat filter, applying both to a different and independent set of numerical simulations (SLICS, cosmo-SLICS). Results. We show that the adapted filter yields a better correlation between the total matter and the galaxy distribution. Furthermore, the adapted filter provides a larger signal-to-noise ratio to constrain the bias between the total matter and the galaxy distribution, and we show that it is, in general, a more sensitive discriminator between different cosmologies, with the exception of cosmologies with very large $σ_8$ values. All analyses lead to the conclusion that our adapted filter should be favored in future density split statistic works. △ Less

Submitted 3 September, 2020; v1 submitted 18 June, 2020; originally announced June 2020.

Comments: 13 pages, 15 figures, accepted for publication by A&A

Journal ref: A&A 642, A161 (2020)

KiDS+GAMA: The weak lensing calibrated stellar-to-halo mass relation of central and satellite galaxies

Authors: Andrej Dvornik, Henk Hoekstra, Konrad Kuijken, Angus H. Wright, Marika Asgari, Maciej Bilicki, Thomas Erben, Benjamin Giblin, Alister W. Graham, Catherine Heymans, Hendrik Hildebrandt, Andrew M. Hopkins, Arun Kannawadi, Chieh-An Lin, Edward N. Taylor, Tilman Tröster

Abstract: We simultaneously present constraints on the stellar-to-halo mass relation for central and satellite galaxies through a weak lensing analysis of spectroscopically classified galaxies. Using overlapping data from the fourth data release of the Kilo-Degree Survey (KiDS), and the Galaxy And Mass Assembly survey (GAMA), we find that satellite galaxies are hosted by halo masses that are… ▽ More We simultaneously present constraints on the stellar-to-halo mass relation for central and satellite galaxies through a weak lensing analysis of spectroscopically classified galaxies. Using overlapping data from the fourth data release of the Kilo-Degree Survey (KiDS), and the Galaxy And Mass Assembly survey (GAMA), we find that satellite galaxies are hosted by halo masses that are $0.53 \pm 0.39$ dex (68\% confidence, $3σ$ detection) smaller than those of central galaxies of the same stellar mass (for a stellar mass of $\log(M_{\star}/M_{\odot}) = 10.6$). This is consistent with galaxy formation models, whereby infalling satellite galaxies are preferentially stripped of their dark matter. We find consistent results with similar uncertainties when comparing constraints from a standard azimuthally averaged galaxy-galaxy lensing analysis and a two-dimensional likelihood analysis of the full shear field. As the latter approach is somewhat biased due to the lens incompleteness and as it does not provide any improvement to the precision when applied to actual data, we conclude that stacked tangential shear measurements are best-suited for studies of the galaxy-halo connection. △ Less

Submitted 18 August, 2020; v1 submitted 18 June, 2020; originally announced June 2020.

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

Journal ref: A&A 642, A83 (2020)

Testing gravity using galaxy-galaxy lensing and clustering amplitudes in KiDS-1000, BOSS and 2dFLenS

Authors: Chris Blake, Alexandra Amon, Marika Asgari, Maciej Bilicki, Andrej Dvornik, Thomas Erben, Benjamin Giblin, Karl Glazebrook, Catherine Heymans, Hendrik Hildebrandt, Benjamin Joachimi, Shahab Joudaki, Arun Kannawadi, Konrad Kuijken, Chris Lidman, David Parkinson, HuanYuan Shan, Tilman Tröster, Jan Luca van den Busch, Christian Wolf, Angus H. Wright

Abstract: The physics of gravity on cosmological scales affects both the rate of assembly of large-scale structure, and the gravitational lensing of background light through this cosmic web. By comparing the amplitude of these different observational signatures, we can construct tests that can distinguish general relativity from its potential modifications. We used the latest weak gravitational lensing data… ▽ More The physics of gravity on cosmological scales affects both the rate of assembly of large-scale structure, and the gravitational lensing of background light through this cosmic web. By comparing the amplitude of these different observational signatures, we can construct tests that can distinguish general relativity from its potential modifications. We used the latest weak gravitational lensing dataset from the Kilo-Degree Survey, KiDS-1000, in conjunction with overlapping galaxy spectroscopic redshift surveys BOSS and 2dFLenS, to perform the most precise existing amplitude-ratio test. We measured the associated E_G statistic with 15-20% errors, in five dz = 0.1 tomographic redshift bins in the range 0.2 < z < 0.7, on projected scales up to 100 Mpc/h. The scale-independence and redshift-dependence of these measurements are consistent with the theoretical expectation of general relativity in a Universe with matter density Omega_m = 0.27 +/- 0.04. We demonstrate that our results are robust against different analysis choices, including schemes for correcting the effects of source photometric redshift errors, and compare the performance of angular and projected galaxy-galaxy lensing statistics. △ Less

Submitted 27 August, 2020; v1 submitted 28 May, 2020; originally announced May 2020.

Comments: 23 pages, 14 figures, version accepted for publication by A&A

Journal ref: A&A 642, A158 (2020)

On the road to percent accuracy IV: ReACT -- computing the non-linear power spectrum beyond $Λ$CDM

Authors: Benjamin Bose, Matteo Cataneo, Tilman Tröster, Qianli Xia, Catherine Heymans, Lucas Lombriser

Abstract: To effectively exploit large-scale structure surveys, we depend on accurate and reliable predictions of non-linear cosmological structure formation. Tools for efficient and comprehensive computational modelling are therefore essential to perform cosmological parameter inference analyses. We present the public software package ReACT, demonstrating its capability for the fast and accurate calculatio… ▽ More To effectively exploit large-scale structure surveys, we depend on accurate and reliable predictions of non-linear cosmological structure formation. Tools for efficient and comprehensive computational modelling are therefore essential to perform cosmological parameter inference analyses. We present the public software package ReACT, demonstrating its capability for the fast and accurate calculation of non-linear power spectra from non-standard physics. We showcase ReACT through a series of analyses on the DGP and $f(R)$ gravity models, adopting LSST-like cosmic shear power spectra. Accurate non-linear modelling with ReACT has the potential to more than double LSST's constraining power on the $f(R)$ parameter, in contrast to an analysis that is limited to the quasi-linear regime. We find that ReACT is sufficiently robust for the inference of consistent constraints on theories beyond $Λ$CDM for current and ongoing surveys. With further improvement, particularly in terms of the accuracy of the non-linear $Λ$CDM power spectrum, ReACT can, in principle, meet the accuracy requirements for future surveys such as Euclid and LSST. △ Less

Submitted 2 September, 2020; v1 submitted 25 May, 2020; originally announced May 2020.

Comments: 14 pages, 5 figures, MNRAS accepted version. ReACT code available here: https://github.com/nebblu/ReACT

KiDS+VIKING-450: Improved cosmological parameter constraints from redshift calibration with self-organising maps

Authors: Angus H. Wright, Hendrik Hildebrandt, Jan Luca van den Busch, Catherine Heymans, Benjamin Joachimi, Arun Kannawadi, Konrad Kuijken

Abstract: We present updated cosmological constraints for the KiDS+VIKING-450 cosmic shear data set (KV450), estimated using redshift distributions and photometric samples defined using self-organising maps (SOMs). Our fiducial analysis finds marginal posterior constraints of $S_8\equivσ_8\sqrt{Ω_{\rm m}/0.3}=0.716^{+0.043}_{-0.038}$; smaller than, but otherwise consistent with, previous work using this dat… ▽ More We present updated cosmological constraints for the KiDS+VIKING-450 cosmic shear data set (KV450), estimated using redshift distributions and photometric samples defined using self-organising maps (SOMs). Our fiducial analysis finds marginal posterior constraints of $S_8\equivσ_8\sqrt{Ω_{\rm m}/0.3}=0.716^{+0.043}_{-0.038}$; smaller than, but otherwise consistent with, previous work using this data set ($|ΔS_8| = 0.023$). We analyse additional samples and redshift distributions constructed in three ways: excluding certain spectroscopic surveys during redshift calibration, excluding lower-confidence spectroscopic redshifts in redshift calibration, and considering only photometric sources which are jointly calibrated by at least three spectroscopic surveys. In all cases, the method utilised here proves robust: we find a maximal deviation from our fiducial analysis of $|ΔS_8| \leq 0.011$ for all samples defined and analysed using our SOM. To demonstrate the reduction in systematic biases found within our analysis, we highlight our results when performing redshift calibration without the DEEP2 spectroscopic data set. In this case we find marginal posterior constraints of $S_8=0.707_{-0.042}^{+0.046}$; a difference with respect to the fiducial that is both significantly smaller than, and in the opposite direction to, the equivalent shift from previous work. These results suggest that our improved cosmological parameter estimates are insensitive to pathological misrepresentation of photometric sources by the spectroscopy used for direct redshift calibration, and therefore that this systematic effect cannot be responsible for the observed difference between $S_8$ estimates made with KV450 and Planck CMB probes. △ Less

Submitted 22 July, 2020; v1 submitted 8 May, 2020; originally announced May 2020.

Comments: 10 pages, 3 figures, 4 appendices, accepted for publication in A&A Letters

Journal ref: A&A 640, L14 (2020)

KiDS+VIKING+GAMA: Testing semi-analytic models of galaxy evolution with galaxy-galaxy-galaxy lensing

Authors: Laila Linke, Patrick Simon, Peter Schneider, Thomas Erben, Daniel J. Farrow, Catherine Heymans, Hendrik Hildebrandt, Andrew M. Hopkins, Arun Kannawadi, Nicola R. Napolitano, Cristóbal Sifón, Angus H. Wright

Abstract: Several semi-analytic models (SAMs) try to explain how galaxies form, evolve and interact inside the dark matter large-scale structure. These SAMs can be tested by comparing their predictions for galaxy-galaxy-galaxy-lensing (G3L), which is weak gravitational lensing around galaxy pairs, with observations. We evaluate the SAMs by Henriques et al. (2015; H15) and by Lagos et al. (2012; L12), implem… ▽ More Several semi-analytic models (SAMs) try to explain how galaxies form, evolve and interact inside the dark matter large-scale structure. These SAMs can be tested by comparing their predictions for galaxy-galaxy-galaxy-lensing (G3L), which is weak gravitational lensing around galaxy pairs, with observations. We evaluate the SAMs by Henriques et al. (2015; H15) and by Lagos et al. (2012; L12), implemented in the Millennium Run, by comparing their predictions for G3L to observations at smaller scales than previous studies and also for pairs of lens galaxies from different populations. We compare the G3L signal predicted by the SAMs to measurements in the overlap of the Galaxy And Mass Assembly survey (GAMA), the Kilo-Degree Survey (KiDS), and the VISTA Kilo-degree Infrared Galaxy survey (VIKING), splitting lens galaxies into two colour and five stellar-mass samples. Using an improved G3L estimator, we measure the three-point correlation of the matter distribution for mixed lens pairs with galaxies from different samples, and unmixed lens pairs with galaxies from the same sample. Predictions by the H15 SAM agree with the observations for all colour-selected and all but one stellar-mass-selected sample with 95% confidence. Deviations occur for lenses with stellar masses below $9.5h^{-2}\mathrm{M}_\odot$ at scales below $0.2h^{-1}\mathrm{Mpc}$. Predictions by the L12 SAM for stellar-mass selected samples and red galaxies are significantly higher than observed, while the predicted signal for blue galaxy pairs is too low. The L12 SAM predicts more pairs of small stellar-mass and red galaxies than the H15 SAM and the observations, as well as fewer pairs of blue galaxies. Likely explanations are different treatments of environmental effects by the SAMs and different models of the initial mass function. We conclude that G3L provides a stringent test for models of galaxy formation and evolution. △ Less

Submitted 12 June, 2020; v1 submitted 5 May, 2020; originally announced May 2020.

Comments: 14 pages, 8 figures, replaced with version accepted to Astronomy & Astrophysics after considering referees comments

Journal ref: A&A 640, A59 (2020)

A hydrodynamical halo model for weak-lensing cross correlations

Authors: A. J. Mead, T. Tröster, C. Heymans, L. Van Waerbeke, I. G. McCarthy

Abstract: On the scale of galactic haloes, the distribution of matter in the cosmos is affected by energetic, non-gravitational processes; so-called baryonic feedback. A lack of knowledge about the details of how feedback processes redistribute matter is a source of uncertainty for weak-lensing surveys, which accurately probe the clustering of matter in the Universe over a wide range of scales. We develop a… ▽ More On the scale of galactic haloes, the distribution of matter in the cosmos is affected by energetic, non-gravitational processes; so-called baryonic feedback. A lack of knowledge about the details of how feedback processes redistribute matter is a source of uncertainty for weak-lensing surveys, which accurately probe the clustering of matter in the Universe over a wide range of scales. We develop a cosmology-dependent model for the matter distribution that simultaneously accounts for the clustering of dark matter, gas and stars. We inform our model by comparing it to power spectra measured from the BAHAMAS suite of hydrodynamical simulations. As well as considering matter power spectra, we also consider spectra involving the electron-pressure field, which directly relates to the thermal Sunyaev-Zel'dovich (tSZ) effect. We fit parameters in our model so that it can simultaneously model both matter and pressure data and such that the distribution of gas as inferred from tSZ has influence on the matter spectrum predicted by our model. We present two variants; one that matches the feedback-induced suppression seen in the matter-matter power spectrum at the per-cent level and a second that matches the matter-matter data slightly less well (~2 per cent), but that is able to simultaneously model the matter-electron pressure spectrum at the ~15 per-cent level. We envisage our models being used to simultaneously learn about cosmological parameters and the strength of baryonic feedback using a combination of tSZ and lensing auto- and cross-correlation data. △ Less

Submitted 3 July, 2020; v1 submitted 30 April, 2020; originally announced May 2020.

Comments: Main paper: 21 pages, 11 figures; Appendix: 16 pages, 13 figures; matches accepted version - referee report addressed, clarifications and extra text compared to v1, improved model, no significant changes of results

Journal ref: A&A 641, A130 (2020)