Source author record

Peter Schneider

Peter Schneider appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

80works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

80 published item(s)

preprint2025arXiv

KiDS-Legacy: Cosmological constraints from cosmic shear with the complete Kilo-Degree Survey

We present cosmic shear constraints from the completed Kilo-Degree Survey (KiDS), where the cosmological parameter $S_8\equivσ_8\sqrt{Ω_{\rm m}/0.3} = 0.815^{+0.016}_{-0.021}$, is found to be in agreement ($0.73σ$) with results from the Planck Legacy cosmic microwave background experiment. The final KiDS footprint spans $1347$ square degrees of deep nine-band imaging across the optical and near-infrared, along with an extra $23$ square degrees of KiDS-like calibration observations of deep spectroscopic surveys. Improvements in our redshift distribution estimation methodology, combined with our enhanced calibration data and multi-band image simulations, allow us to extend our lensed sample out to a photometric redshift of $z_{\rm B}\leq2.0$. Compared to previous KiDS analyses, the increased survey area and redshift depth results in a $\sim32\%$ improvement in constraining power in terms of $Σ_8\equivσ_8\left(Ω_{\rm m}/0.3\right)^α= 0.821^{+0.014}_{-0.016}$, where $α= 0.58$ has been optimised to match the revised degeneracy direction of $σ_8$ and $Ω_{\rm m}$. We adopt a new physically motivated intrinsic alignment model that depends jointly on the galaxy sample's halo mass and spectral type distributions, and that is informed by previous direct alignment measurements. We also marginalise over our uncertainty on the impact of baryon feedback on the non-linear matter power spectrum. Comparing to previous KiDS analyses, we conclude that the increase seen in $S_8$ primarily results from our improved redshift distribution estimation and calibration, as well as new survey area and improved image reduction. Our companion paper Stölzner et al. (submitted) presents a full suite of internal and external consistency tests, finding the KiDS-Legacy data set to be the most internally robust sample produced by KiDS to date.

preprint2025arXiv

Quantifying the impact of detection bias from blended galaxies on cosmic shear surveys

Increasingly large areas in cosmic shear surveys lead to a reduction of statistical errors, necessitating to control systematic errors increasingly better. One of these systematic effects was initially studied by Hartlap et al. in 2011, namely that image overlap with (bright foreground) galaxies may prevent some distant (source) galaxies to remain undetected. Since this overlap is more likely to occur in regions of high foreground density -- which tend to be the regions in which the shear is largest -- this detection bias would cause an underestimation of the estimated shear correlation function. This detection bias adds to the possible systematic of image blending, where nearby pairs or multiplets of images render shear estimates more uncertain and thus may cause a reduction in their statistical weight. Based on simulations with data from the Kilo-Degree Survey, we study the conditions under which images are not detected. We find an approximate analytic expression for the detection probability in terms of the separation and brightness ratio to the neighbouring galaxies. Applying this fitting formula to weak lensing ray tracing through, and the galaxy distribution in the Millennium Simulation, we estimate that the detection bias alone leads to an underestimate of $S_8=σ_8\sqrt{Ω_\mathrm{m}/0.3}$ by almost 2\% and can therefore not be neglected in current and forthcoming cosmic shear surveys.

preprint2023arXiv

KiDS-1000 Cosmology: Constraints from density split statistics

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}$.

preprint2022arXiv

A revised density split statistic model for general filters

Studying the statistical properties of the large-scale structure in the Universe with weak gravitational lensing is a prime goal of several current and forthcoming galaxy surveys. The power that weak lensing has to constrain cosmological parameters can be enhanced by considering statistics beyond second-order shear correlation functions or power spectra. One such higher-order probe that has proven successful in observational data is the density split statistics (DSS), in which one analyses the mean shear profiles around points that are classified according to their foreground galaxy density. In this paper, we generalise the most accurate DSS model to allow for a broad class of angular filter functions used for the classification of the different local density regions. This approach is motivated by earlier findings showing that an optimised filter can provide tighter constraints on model parameters compared to the standard top-hat case. We build on large deviation theory approaches and approximations thereof to model the matter density PDF, and on perturbative calculations of higher-order moments of the density field. The novel addition relies on the generalisation of these previously employed calculations to allow for general filter functions and is validated on several sets of numerical simulations. The revised model fits well the simulation measurements, with a residual systematic offset that is small compared to the statistical accuracy of current weak lensing surveys. The accuracy of the model is slightly lower for a compensated filter than for a non-negative filter function, and that it increases with the filter size. Using a Fisher matrix approach, we find constraints comparable to the commonly used two-point cosmic shear measures. Hence, our DSS model can be used in competitive analyses of current cosmic shear data, while it may need refinements for forthcoming lensing surveys.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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

preprint2022arXiv

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

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.

preprint2021arXiv

Probing dark energy with tomographic weak-lensing aperture mass statistics

We forecast and optimize the cosmological power of various weak-lensing aperture mass ($M_{\rm ap}$) map statistics for future cosmic shear surveys, including peaks, voids, and the full distribution of pixels (1D $M_{\rm ap}$). These alternative methods probe the non-Gaussian regime of the matter distribution, adding complementary cosmological information to the classical two-point estimators. Based on the SLICS and cosmo-SLICS $N$-body simulations, we build Euclid-like mocks to explore the $S_8 - Ω_{\rm m} - w_0$ parameter space. We develop a new tomographic formalism which exploits the cross-information between redshift slices (cross-$M_{\rm ap}$) in addition to the information from individual slices (auto-$M_{\rm ap}$) probed in the standard approach. Our auto-$M_{\rm ap}$ forecast precision is in good agreement with the recent literature on weak-lensing peak statistics, and is improved by $\sim 50$% when including cross-$M_{\rm ap}$. It is further boosted by the use of 1D $M_{\rm ap}$ that outperforms all other estimators, including the shear two-point correlation function ($γ$-2PCF). When considering all tomographic terms, our uncertainty range on the structure growth parameter $S_8$ is enhanced by $\sim 45$% (almost twice better) when combining 1D $M_{\rm ap}$ and the $γ$-2PCF compared to the $γ$-2PCF alone. We additionally measure the first combined forecasts on the dark energy equation of state $w_0$, finding a factor of three reduction of the statistical error compared to the $γ$-2PCF alone. This demonstrates that the complementary cosmological information explored by non-Gaussian $M_{\rm ap}$ map statistics not only offers the potential to improve the constraints on the recent $σ_8$ - $Ω_{\rm m}$ tension, but also constitutes an avenue to understand the accelerated expansion of our Universe.

preprint2020arXiv

Derived parabolic induction

The classical parabolic induction functor is a fundamental tool on the representation theoretic side of the Langlands program. In this article, we study its derived version. It was shown by the second author that the derived category of smooth $G$-representations over $k$, $G$ a $p$-adic reductive group and $k$ a field of characteristic $p$, is equivalent to the derived category of a certain differential graded $k$-algebra $H_G^\bullet$, whose zeroth cohomology is a classical Hecke algebra. This equivalence predicts the existence of a derived parabolic induction functor on the dg Hecke algebra side, which we construct in this paper. This relies on the theory of six-functor formalisms for differential graded categories developed by O.\ Schnürer. We also discuss the adjoint functors of derived parabolic induction.

preprint2020arXiv

KiDS-1000 catalogue: Weak gravitational lensing shear measurements

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

KiDS+VIKING-450 and DES-Y1 combined: Mitigating baryon feedback uncertainty with COSEBIs

We present cosmological constraints from a joint cosmic shear analysis of the Kilo-Degree Survey (KV450) and the Dark Energy Survey (DES-Y1), conducted using Complete Orthogonal Sets of E/B-Integrals (COSEBIs). With COSEBIs we isolate any B-modes which have a non-cosmic shear origin and demonstrate the robustness of our cosmological E-mode analysis as no significant B-modes are detected. We highlight how COSEBIs are fairly insensitive to the amplitude of the non-linear matter power spectrum at high $k$-scales, mitigating the uncertain impact of baryon feedback in our analysis. COSEBIs, therefore, allow us to utilise additional small-scale information, improving the DES-Y1 joint constraints on $S_8=σ_8(Ω_{\rm m}/0.3)^{0.5}$ and $Ω_{\rm m}$ by $20\%$. Adopting a flat $Λ$CDM model we find $S_8=0.755^{+0.019}_{-0.021}$, which is in $3.2σ$ tension with the Planck Legacy analysis of the cosmic microwave background.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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]

preprint2019arXiv

A gravitational lensing detection of filamentary structures connecting luminous red galaxies

We present a weak lensing detection of filamentary structures in the cosmic web, combining data from the Kilo-Degree Survey, the Red Cluster Sequence Lensing Survey and the Canada-France-Hawaii Telescope Lensing Survey. The line connecting luminous red galaxies with a separation of $3 - 5\, h^{-1}\text{Mpc}$ is chosen as a proxy for the location of filaments. We measure the average weak lensing shear around $\sim$11,000 candidate filaments selected in this way from the Sloan Digital Sky Survey. After nulling the shear induced by the dark matter haloes around each galaxy, we report a $3.4\,σ$ detection of an anisotropic shear signal from the matter that connects them. Adopting a filament density profile, motivated from $N$-body simulations, the average density at the centre of these filamentary structures is found to be $15 \pm 4$ times the critical density.

preprint2019arXiv

Measuring galaxy-galaxy-galaxy-lensing with higher precision and accuracy

Galaxy-galaxy-galaxy lensing (G3L) is a powerful tool for constraining the three-point correlation between the galaxy and matter distribution and thereby models of galaxy evolution. We propose three improvements to current measurements of G3L: (i) a weighting of lens galaxies according to their redshift difference, (ii) adaptive binning of the three-point correlation function, and (iii) accounting for the effect of lens magnification by the cosmic large-scale structure. Improvement (i) is designed to improve the precision of the G3L measurement, whereas improvements (ii) and (iii) remove biases of the estimator. We further show how the G3L signal can be converted from angular into physical scales. The improvements were tested on simple mock data and simulated data based on the Millennium Run with an implemented semi-analytic galaxy model. Our improvements increase the signal-to-noise ratio by 35 % on average at angular scales between 0.1 arcmin and 10 arcmin and physical scales between $0.02$ and $2 \, h^{-1}\, \textrm{Mpc}$. They also remove the bias of the G3L estimator at angular scales below 1 arcmin, which was originally up to 40 %. The signal due to lens magnification is approximately 10 % of the total signal.

preprint2019arXiv

The effects of varying depth in cosmic shear surveys

We present a semi-analytic model for the shear two-point correlation function of a cosmic shear survey with non-uniform depth. Ground-based surveys are subject to depth variations that primarily arise through varying atmospheric conditions. For a survey like the Kilo-Degree Survey (KiDS), we find that the measured depth variation increases the amplitude of the observed shear correlation function at the level of a few percent out to degree-scales, relative to the assumed uniform-depth case. The impact on the inferred cosmological parameters is shown to be insignificant for a KiDS-like survey. For next-generation cosmic shear experiments, however, we conclude that variable depth should be accounted for.

preprint2016arXiv

A canonical torsion theory for pro-p Iwahori-Hecke modules

Let $\mathfrak F$ be a locally compact nonarchimedean field with residue characteristic $p$ and $G$ the group of $\mathfrak{F}$-rational points of a connected split reductive group over $\mathfrak{F}$. We define a torsion pair in the category Mod$(H)$ of modules over the pro-$p$-Iwahori Hecke $k$-algebra $H$ of $G$, where $k$ is an arbitrary field. We prove that, under a certain hypothesis, the torsionfree class embeds fully faithfully into the category Mod${}^I(G)$ of smooth $k$-representations of $G$ generated by their pro-$p$-Iwahori fixed vectors. If the characteristic of $k$ is different from $p$ then this hypothesis is always satisfied and the torsionfree class is the whole category Mod$(H)$. If $k$ contains the residue field of $\mathfrak F$ then we study the case $G = \mathbf{SL_2}(\mathfrak F)$. We show that our hypothesis is satisfied, and we describe explicitly the torsionfree and the torsion classes. If $\mathfrak F\neq \mathbb Q_p$ and $p\neq 2$, then an $H$-module is in the torsion class if and only if it is a union of supersingular finite length submodules; it lies in the torsionfree class if and only if it does not contain any nonzero supersingular finite length module. If $\mathfrak{F} = \mathbb{Q}_p$, the torsionfree class is the whole category Mod$(H)$, and we give a new proof of the fact that Mod$(H)$ is equivalent to Mod${}^I(G)$. These results are based on the computation of the $H$-module structure of certain natural cohomology spaces for the pro-$p$-Iwahori subgroup $I$ of $G$.

preprint2016arXiv

Dependence of GAMA galaxy halo masses on the cosmic web environment from 100 square degrees of KiDS weak lensing data

Galaxies and their dark matter haloes are part of a complex network of mass structures, collectively called the cosmic web. Using the tidal tensor prescription these structures can be classified into four cosmic environments: voids, sheets, filaments and knots. As the cosmic web may influence the formation and evolution of dark matter haloes and the galaxies they host, we aim to study the effect of these cosmic environments on the average mass of galactic haloes. To this end we measure the galaxy-galaxy lensing profile of 91,195 galaxies, within 0.039 < z < 0.263, from the spectroscopic Galaxy And Mass Assembly (GAMA) survey, using ~100 square degrees of overlapping data from the Kilo-Degree Survey (KiDS). In each of the four cosmic environments we model the contributions from group centrals, satellites and neighbouring groups to the stacked galaxy-galaxy lensing profiles. After correcting the lens samples for differences in the stellar mass distribution, we find no dependence of the average halo mass of central galaxies on their cosmic environment. We do find a significant increase in the average contribution of neighbouring groups to the lensing profile in increasingly dense cosmic environments. We show, however, that the observed effect can be entirely attributed to the galaxy density at much smaller scales (within 4 Mpc/h), which is correlated with the density of the cosmic environments. Within our current uncertainties we find no direct dependence of galaxy halo mass on their cosmic environment.

preprint2016arXiv

Generalized shear-ratio tests: A new relation between cosmological distances, and a diagnostic for a redshift-dependent multiplicative bias in shear measurements

We derive a new relation between cosological distances, valid in any (statistically) isotropic space-time and independent of cosmological parameters or even the validity of the field equation of General Relativity. In particular, this relation yields an equation between those distance ratios which are the geometrical factors determining the strength of the gravitational lensing effect of mass concentrations. Considering a combination of weak lensing shear ratios, based on lenses at two different redshifts, and sources at three different redshifts, we derive a relation between shear-ratio tests which must be identically satisfied. A redshift-dependent multiplicative bias in shear estimates will violate this relation, and thus can be probed by this generalized shear-ratio test. Combining the lensing effect for lenses at three different redshifts and three different source redshifts, a relation between shear ratios is derived which must be valid independent of a multiplicative bias. We propose these generalized shear-ratio tests as a diagnostic for the presence of systematics in upcoming weak lensing surveys.

preprint2016arXiv

Revisiting CFHTLenS cosmic shear: Optimal E/B mode decomposition using COSEBIs and compressed COSEBIs

We present a re-analysis of the CFHTLenS weak gravitational lensing survey using Complete Orthogonal Sets of E/B-mode Integrals, known as COSEBIs. COSEBIs provide a complete set of functions to efficiently separate E-modes from B-modes and hence allow for robust and stringent tests for systematic errors in the data. This analysis reveals significant B-modes on large angular scales that were not previously seen using the standard E/B decomposition analyses. We find that the significance of the B-modes is enhanced when the data is split by galaxy type and analysed in tomographic redshift bins. Adding tomographic bins to the analysis increases the number of COSEBIs modes, which results in a less accurate estimation of the covariance matrix from a set of simulations. We therefore also present the first compressed COSEBIs analysis of survey data, where the COSEBIs modes are optimally combined based on their sensitivity to cosmological parameters. In this tomographic CCOSEBIs analysis we find the B-modes to be consistent with zero when the full range of angular scales are considered.

preprint2015arXiv

Coates-Wiles homomorphisms and Iwasawa cohomology for Lubin-Tate extensions

For the $p$-cyclotomic tower of $\mathbb{Q}_p$ Fontaine established a description of local Iwasawa cohomology with coefficients in a local Galois representation $V$ in terms of the $ψ$-operator acting on the attached etale $(φ,Γ)$-module $D(V)$. In this article we generalize Fontaine's result to the case of arbitratry Lubin-Tate towers $L_\infty$ over finite extensions $L$ of $\mathbb{Q}_p$ by using the Kisin-Ren/Fontaine equivalence of categories between Galois representations and $(φ_L,Γ_L)$-module and extending parts of [Herr L.: Sur la cohomologie galoisienne des corps $p$-adiques. Bull. Soc. Math. France 126, 563-600 (1998)], [Scholl A. J.: Higher fields of norms and $(ϕ,Γ)$-modules. Documenta Math.\ 2006, Extra Vol., 685-709]. Moreover, we prove a kind of explicit reciprocity law which calculates the Kummer map over $L_\infty$ for the multiplicative group twisted with the dual of the Tate module $T$ of the Lubin-Tate formal group in terms of Coleman power series and the attached $(φ_L,Γ_L)$-module. The proof is based on a generalized Schmid-Witt residue formula. Finally, we extend the explicit reciprocity law of Bloch and Kato [Bloch S., Kato K.: $L$-functions and Tamagawa numbers of motives. The Grothendieck Festschrift, Vol. I, 333-400, Progress Math., 86, Birkhäuser Boston 1990] Thm. 2.1 to our situation expressing the Bloch-Kato exponential map for $L(χ_{LT}^r)$ in terms of generalized Coates-Wiles homomorphisms, where the Lubin-Tate characater $χ_{LT}$ describes the Galois action on $T.$

preprint2015arXiv

Constrained correlation functions from the Millennium Simulation

Context. In previous work, we developed a quasi-Gaussian approximation for the likelihood of correlation functions, which, in contrast to the usual Gaussian approach, incorporates fundamental mathematical constraints on correlation functions. The analytical computation of these constraints is only feasible in the case of correlation functions of one-dimensional random fields. Aims. In this work, we aim to obtain corresponding constraints in the case of higher-dimensional random fields and test them in a more realistic context. Methods. We develop numerical methods to compute the constraints on correlation functions which are also applicable for two- and three-dimensional fields. In order to test the accuracy of the numerically obtained constraints, we compare them to the analytical results for the one-dimensional case. Finally, we compute correlation functions from the halo catalog of the Millennium Simulation, check whether they obey the constraints, and examine the performance of the transformation used in the construction of the quasi-Gaussian likelihood. Results. We find that our numerical methods of computing the constraints are robust and that the correlation functions measured from the Millennium Simulation obey them. Despite the fact that the measured correlation functions lie well inside the allowed region of parameter space, i.e. far away from the boundaries of the allowed volume defined by the constraints, we find strong indications that the quasi-Gaussian likelihood yields a substantially more accurate description than the Gaussian one.

preprint2015arXiv

Gravitational Lensing Analysis of the Kilo Degree Survey

The Kilo-Degree Survey (KiDS) is a multi-band imaging survey designed for cosmological studies from weak lensing and photometric redshifts. It uses the ESO VLT Survey Telescope with its wide-field camera OmegaCAM. KiDS images are taken in four filters similar to the SDSS ugri bands. The best-seeing time is reserved for deep r-band observations that reach a median 5-sigma limiting AB magnitude of 24.9 with a median seeing that is better than 0.7arcsec. Initial KiDS observations have concentrated on the GAMA regions near the celestial equator, where extensive, highly complete redshift catalogues are available. A total of 109 survey tiles, one square degree each, form the basis of the first set of lensing analyses, which focus on measurements of halo properties of GAMA galaxies. 9 galaxies per square arcminute enter the lensing analysis, for an effective inverse shear variance of 69 per square arcminute. Accounting for the shape measurement weight, the median redshift of the sources is 0.53. KiDS data processing follows two parallel tracks, one optimized for galaxy shape measurement (for weak lensing), and one for accurate matched-aperture photometry in four bands (for photometric redshifts). This technical paper describes how the lensing and photometric redshift catalogues have been produced (including an extensive description of the Gaussian Aperture and Photometry pipeline), summarizes the data quality, and presents extensive tests for systematic errors that might affect the lensing analyses. We also provide first demonstrations of the suitability of the data for cosmological measurements, and explain how the shear catalogues were blinded to prevent confirmation bias in the scientific analyses. The KiDS shear and photometric redshift catalogues, presented in this paper, are released to the community through http://kids.strw.leidenuniv.nl .

preprint2015arXiv

How well can cold-dark-matter substructures account for the observed radio flux-ratio anomalies?

Discrepancies between the observed and model-predicted radio flux ratios are seen in a number of quadruply-lensed quasars. The most favored interpretation of these anomalies is that CDM substructures present in lensing galaxies perturb the lens potentials and alter image magnifications and thus flux ratios. So far no consensus has emerged regarding whether or not the predicted CDM substructure abundance fully accounts for the lensing flux anomaly observations. Accurate modeling relies on a realistic lens sample in terms of both the lens environment and internal structures and substructures. In this paper we construct samples of generalised and specific lens potentials, to which we add (rescaled) subhalo populations from the galaxy-scale Aquarius and the cluster-scale Phoenix simulation suites. We further investigate the lensing effects from subhalos of masses several orders of magnitude below the simulation resolution limit. The resulting flux ratio distributions are compared to the currently best available sample of radio lenses. The observed anomalies in B0128+437, B0712+472 and B1555+375 are more likely to be caused by propagation effects or oversimplified lens modeling, signs of which are already seen in the data. Among the quadruple systems that have closely located image triplets/pairs, the anomalous flux ratios of MG0414+0534 can be reproduced by adding CDM subhalos to its macroscopic lens potential, with a probability of 5%-20%; for B0712+472, B1422+231, B1555+375 and B2045+265, these probabilities are only of a few percent. We hence find that CDM substructures are unlikely to be the whole reason for radio flux anomalies. We discuss other possible effects that might also be at work.

preprint2015arXiv

Rigid character groups, Lubin-Tate theory, and $(φ,Γ)$-modules

The construction of the $p$-adic local Langlands correspondence for $\mathrm{GL}_2(\mathbf{Q}_p)$ uses in an essential way Fontaine's theory of cyclotomic $(φ,Γ)$-modules. Here \emph{cyclotomic} means that $Γ= \mathrm{Gal}(\mathbf{Q}_p(μ_{p^\infty})/\mathbf{Q}_p)$ is the Galois group of the cyclotomic extension of $\mathbf{Q}_p$. In order to generalize the $p$-adic local Langlands correspondence to $\mathrm{GL}_2(L)$, where $L$ is a finite extension of $\mathbf{Q}_p$, it seems necessary to have at our disposal a theory of Lubin-Tate $(φ,Γ)$-modules. Such a generalization has been carried out to some extent, by working over the $p$-adic open unit disk, endowed with the action of the endomorphisms of a Lubin-Tate group. The main idea of our article is to carry out a Lubin-Tate generalization of the theory of cyclotomic $(φ,Γ)$-modules in a different fashion. Instead of the $p$-adic open unit disk, we work over a character variety, that parameterizes the locally $L$-analytic characters on $o_L$. We study $(φ,Γ)$-modules in this setting, and relate some of them to what was known previously.

preprint2015arXiv

The first and second data releases of the Kilo-Degree Survey

The Kilo-Degree Survey (KiDS) is an optical wide-field imaging survey carried out with the VLT Survey Telescope and the OmegaCAM camera. KiDS will image 1500 square degrees in four filters (ugri), and together with its near-infrared counterpart VIKING will produce deep photometry in nine bands. Designed for weak lensing shape and photometric redshift measurements, the core science driver of the survey is mapping the large-scale matter distribution in the Universe back to a redshift of ~0.5. Secondary science cases are manifold, covering topics such as galaxy evolution, Milky Way structure, and the detection of high-redshift clusters and quasars. KiDS is an ESO Public Survey and dedicated to serving the astronomical community with high-quality data products derived from the survey data, as well as with calibration data. Public data releases will be made on a yearly basis, the first two of which are presented here. For a total of 148 survey tiles (~160 sq.deg.) astrometrically and photometrically calibrated, coadded ugri images have been released, accompanied by weight maps, masks, source lists, and a multi-band source catalog. A dedicated pipeline and data management system based on the Astro-WISE software system, combined with newly developed masking and source classification software, is used for the data production of the data products described here. The achieved data quality and early science projects based on the data products in the first two data releases are reviewed in order to validate the survey data. Early scientific results include the detection of nine high-z QSOs, fifteen candidate strong gravitational lenses, high-quality photometric redshifts and galaxy structural parameters for hundreds of thousands of galaxies. (Abridged)

preprint2015arXiv

The masses of satellites in GAMA galaxy groups from 100 square degrees of KiDS weak lensing data

We use the first 100 sq. deg. of overlap between the Kilo-Degree Survey (KiDS) and the Galaxy And Mass Assembly (GAMA) survey to determine the galaxy halo mass of ~10,000 spectroscopically-confirmed satellite galaxies in massive ($M > 10^{13}h^{-1}{\rm M}_\odot$) galaxy groups. Separating the sample as a function of projected distance to the group centre, we jointly model the satellites and their host groups with Navarro-Frenk-White (NFW) density profiles, fully accounting for the data covariance. The probed satellite galaxies in these groups have total masses $\log M_{\rm sub} /(h^{-1}{\rm M}_\odot) \approx 11.7 - 12.2$ consistent across group-centric distance within the errorbars. Given their typical stellar masses, $\log M_{\rm \star,sat}/(h^{-2}{\rm M}_\odot) \sim 10.5$, such total masses imply stellar mass fractions of $M_{\rm \star,sat} /M_{\rm sub} \approx 0.04 h^{-1}$ . The average subhalo hosting these satellite galaxies has a mass $M_{\rm sub} \sim 0.015M_{\rm host}$ independent of host halo mass, in broad agreement with the expectations of structure formation in a $Λ$CDM universe.

preprint2014arXiv

A New Data Compression Method and its Application to Cosmic Shear Analysis

Future large scale cosmological surveys will provide huge data sets whose analysis requires efficient data compression. Calculating accurate covariances is extremely challenging with increasing number of statistics used. Here we introduce a formalism for achieving efficient data compression, based on a local expansion of statistical measures around a fiducial cosmological model. We specifically apply and test this approach for the case of cosmic shear statistics. We demonstrate the performance of our approach, using a Fisher analysis on cosmic shear tomography described in terms of E-/B-mode separating statistics (COSEBIs). We show that our data compression is highly effective in extracting essentially the full cosmological information from a strongly reduced number of observables. Specifically, the number of statistics needed decreases by at least one order of magnitude relative to the COSEBIs, which already compress the data substantially compared to the shear two-point correlation functions. The efficiency appears to be affected only slightly if a highly inaccurate covariance is used for defining the compressed statistics, showing the robustness of the method. We conclude that an efficient data compression is achievable and that the number of compressed statistics depends on the number of model parameters. In addition, we study how well band powers can be obtained from measuring shear correlation functions over a finite interval of separations. We show the strong limitations on the possibility to construct top-hat filters in Fourier space, for which the real-space analog has a finite support, yielding strong bounds on the accuracy of band power estimates. The error on an estimated band-power is larger for a narrower filter and a smaller angular range which for relevant cases can be as large as 10%.

preprint2014arXiv

Can one determine cosmological parameters from multi-plane strong lens systems?

Strong gravitational lensing of sources with different redshifts has been used to determine cosmological distance ratios, which in turn depend on the expansion history. Hence, such systems are viewed as potential tools for constraining cosmological parameters. Here we show that in lens systems with two distinct source redshifts, of which the nearest one contributes to the light deflection towards the more distant one, there exists an invariance transformation which leaves all strong lensing observables unchanged (except the product of time delay and Hubble constant), generalizing the well-known mass-sheet transformation in single plane lens systems. The transformation preserves the relative distribution of mass and light, so that a `mass-follows-light' assumption does not fix the MST. All time delays (from sources on both planes) scale with the same factor -- time-delay ratios are therefore invariant under the MST. Changing cosmological parameters, and thus distance ratios, is essentially equivalent to such a mass-sheet transformation. As an example, we discuss the double source plane system SDSSJ0946+1006, which has been recently studied by Collett and Auger, and show that variations of cosmological parameters within reasonable ranges lead to only a small mass-sheet transformation in both lens planes. Hence, the ability to extract cosmological information from such systems depends heavily on the ability to break the mass-sheet degeneracy.

preprint2014arXiv

Combining Probes of Large-Scale Structure with CosmoLike

Developing accurate analysis techniques to combine various probes of cosmology is essential to tighten constraints on cosmological parameters and to check for inconsistencies in our model of the Universe. In this paper we develop a joint analysis framework for six different second-order statistics calculated from three tracers of the dark matter density field, namely galaxy position, shear, and magnification. We extend a data compression scheme developed in the context of shear-shear statistics (the so-called COSEBIs) to the other five second-order statistics, thereby significantly reducing the number of data points in the joint data vector. We use CosmoLike, a newly developed software framework for joint likelihood analyses, to forecast parameter constraints for the Dark Energy Survey (DES). The simulated MCMCs cover a five dimensional cosmological parameter space comparing the information content of the individual probes to several combined probes (CP) data vectors. Given the significant correlations of these second-order statistics we model all cross terms in the covariance matrix; furthermore we go beyond the Gaussian covariance approximation and use the halo model to include higher order correlations of the density field. We find that adding magnification information (including cross probes with shear and clustering) noticeably increases the information content and that the correct modeling of the covariance (i.e., accounting for non-Gaussianity and cross terms) is essential for accurate likelihood contours from the CP data vector. We also identify several nulltests based on the degeneracy of magnification and shear statistics which can be used to quantify the contamination of data sets by astrophysical systematics and/or calibration issues.

preprint2014arXiv

How well can cold-dark-matter substructures account for the observed lensing flux-ratio anomalies?

Lensing flux-ratio anomalies are most likely caused by gravitational lensing by small-scale dark matter structures. These anomalies offer the prospect of testing a fundamental prediction of the cold dark matter (CDM) cosmological model: the existence of numerous substructures that are too small to host visible galaxies. In two previous studies we found that the number of subhalos in the six high-resolution simulations of CDM galactic halos of the Aquarius project is not sufficient to account for the observed frequency of flux ratio anomalies seen in selected quasars from the CLASS survey. These studies were limited by the small number of halos used, their narrow range of masses (1-2E12 solar masses) and the small range of lens ellipticities considered. We address these shortcomings by investigating the lensing properties of a large sample of halos with a wide range of masses in two sets of high resolution simulations of cosmological volumes and comparing them to a currently best available sample of radio quasars. We find that, as expected, substructures do not change the flux-ratio probability distribution of image pairs and triples with large separations, but they have a significant effect on the distribution at small separations. For such systems, CDM substructures can account for a substantial fraction of the observed flux-ratio anomalies. For large close-pair separation systems, the discrepancies existing between the observed flux ratios and predictions from smooth halo models are attributed to simplifications inherent in these models which do not take account of fine details in the lens mass distributions.

preprint2014arXiv

How well do third-order aperture mass statistics separate E- and B-modes?

With 3rd-order statistics of gravitational shear it will be possible to extract valuable cosmological information from ongoing and future weak lensing surveys which is not contained in standard 2nd-order statistics, due to the non-Gaussianity of the shear field. Aperture mass statistics are an appropriate choice for 3rd-order statistics due to their simple form and their ability to separate E- and B-modes of the shear. However, it has been demonstrated that 2nd-order aperture mass statistics suffer from E-/B-mode mixing because it is impossible to reliably estimate the shapes of close pairs of galaxies. This finding has triggered developments of several new 2nd-order statistical measures for cosmic shear. Whether the same developments are needed for 3rd-order shear statistics is largely determined by how severe this E-/B-mixing is for 3rd-order statistics. We test 3rd-order aperture mass statistics against E-/B-mode mixing, and find that the level of contamination is well-described by a function of $θ/θ_{\rm min}$, with $θ_{\rm min}$ being the cut-off scale. At angular scales of $θ> 10 \;θ_{\rm min}$, the decrease in the E-mode signal due to E-/B-mode mixing is smaller than 1 percent, and the leakage into B-modes is even less. For typical small-scale cut-offs this E-/B-mixing is negligible on scales larger than a few arcminutes. Therefore, 3rd-order aperture mass statistics can safely be used to separate E- and B-modes and infer cosmological information, for ground-based surveys as well as forthcoming space-based surveys such as Euclid.

preprint2014arXiv

Simple characters and coefficient systems on the building

Let F be a non-archimedean local field and G be the group GL(N,F). Let πbe a smooth complex representation of G lying in the Bernstein block B(π) of some simple type in the sense of Bushnell and Kutzko. Refining the approach of the second author and U. Stuhler, we canonically attach to πa subset X_πof the Bruhat-Tits building X of G, as well as a G-equivariant coefficient system C[π] on X_π. Roughly speaking the coefficient system is obtained by taking isotypic components of πaccording to some representations constructed from the Bushnell and Kutzko type of π. We conjecture that when πhas central character, the augmented chain complex associate to C[π] is a projective resolution of πin the category B(π). Moreover we reduce this conjecture to a technical lemma of representation theoretic nature. We prove this lemma when πis an irreducible discrete series of G. We then attach to any irreducible discrete series πof G an explicit pseudo-coefficient f_πand obtain a Lefschetz type formula for the value of the Harish-Chandra character of πat a regular elliptic element. In contrast to that obtained by U. Stuhler and the second author, this formula allows explicit character value computations.

preprint2014arXiv

The 400d Galaxy Cluster Survey weak lensing programme: III: Evidence for consistent WL and X-ray masses at $z\approx 0.5$

Scaling properties of galaxy cluster observables with mass provide central insights into the processes shaping clusters. Calibrating proxies for cluster mass will be crucial to cluster cosmology with upcoming surveys like eROSITA and Euclid. The recent Planck results led to suggestions that X-ray masses might be biased low by $\sim\!40$ %, more than previously considered. We extend the direct calibration of the weak lensing -- X-ray mass scaling towards lower masses (as low as $1\!\times\!10^{14}\,\mathrm{M}_{\odot}$) in a sample representative of the $z\!\sim\!0.4$--$0.5$ population. We investigate the scaling of MMT/Megacam weak lensing (WL) masses for $8$ clusters at $0.39\!\leq\!z\!\leq\!0.80$ as part of the \emph{400d} WL programme with hydrostatic \textit{Chandra} X-ray masses as well as those based on the proxies, e.g. $Y_{\mathrm{X}}\!=\!T_{\mathrm{X}}M_{\mathrm{gas}}$. Overall, we find good agreement between WL and X-ray masses, with different mass bias estimators all consistent with zero. Subdividing the sample, we find the high-mass subsample to show no significant mass bias while for the low-mass subsample, there is a bias towards overestimated X-ray masses at the $\sim\!2σ$ level for some mass proxies. The overall scatter in the mass-mass scaling relations is surprisingly low. Neither observation can be traced back to the parameter settings in the WL analysis. We do not find evidence for a strong ($\sim\!40$ %) underestimate in the X-ray masses, as suggested to reconcile Planck cluster counts and cosmological constraints. For high-mass clusters, our measurements are consistent with studies in the literature. The mass dependent bias, significant at $\sim\!2σ$, may hint at a physically different cluster population (less relaxed clusters with more substructure and mergers); or it may be due to small number statistics.

preprint2013arXiv

A quasi-Gaussian approximation for the probability distribution of correlation functions

Context. Whenever correlation functions are used for inference about cosmological parameters in the context of a Bayesian analysis, the likelihood function of correlation functions needs to be known. Usually, it is approximated as a multivariate Gaussian, though this is not necessarily a good approximation. Aims. We show how to calculate a better approximation for the probability distribution of correlation functions, which we call "quasi-Gaussian". Methods. Using the exact univariate PDF as well as constraints on correlation functions previously derived, we transform the correlation functions to an unconstrained variable for which the Gaussian approximation is well justified. From this Gaussian in the transformed space, we obtain the quasi-Gaussian PDF. The two approximations for the probability distributions are compared to the "true" distribution as obtained from simulations. Additionally, we test how the new approximation performs when used as likelihood in a toy-model Bayesian analysis. Results. The quasi-Gaussian PDF agrees very well with the PDF obtained from simulations; in particular, it provides a significantly better description than a straightforward copula approach. In a simple toy-model likelihood analysis, it yields noticeably different results than the Gaussian likelihood, indicating its possible impact on cosmological parameter estimation.

preprint2013arXiv

Galaxy-galaxy(-galaxy) lensing as a sensitive probe of galaxy evolution

The gravitational lensing effect provides various ways to study the mass environment of galaxies. We investigate how galaxy-galaxy(-galaxy) lensing can be used to test models of galaxy formation and evolution. We consider two semi-analytic galaxy formation models based on the Millennium Run N-body simulation: the Durham model by Bower et al. (2006) and the Garching model by Guo et al. (2011). We generate mock lensing observations for the two models, and then employ Fast Fourier Transform methods to compute second- and third-order aperture statistics in the simulated fields for various galaxy samples. We find that both models predict qualitatively similar aperture signals, but there are large quantitative differences. The Durham model predicts larger amplitudes in general. In both models, red galaxies exhibit stronger aperture signals than blue galaxies. Using these aperture measurements and assuming a linear deterministic bias model, we measure relative bias ratios of red and blue galaxy samples. We find that a linear deterministic bias is insufficient to describe the relative clustering of model galaxies below ten arcmin angular scales. Dividing galaxies into luminosity bins, the aperture signals decrease with decreasing luminosity for brighter galaxies, but increase again for fainter galaxies. This increase is likely an artifact due to too many faint satellite galaxies in massive group and cluster halos predicted by the models. Our study shows that galaxy-galaxy(-galaxy) lensing is a sensitive probe of galaxy evolution.

preprint2013arXiv

Mass-sheet degeneracy, power-law models and external convergence: Impact on the determination of the Hubble constant from gravitational lensing

The light travel time differences in strong gravitational lensing systems allows an independent determination of the Hubble constant. This method has been successfully applied to several lens systems. The formally most precise measurements are, however, in tension with the recent determination of $H_0$ from the Planck satellite for a spatially flat six-parameters $ΛCDM$ cosmology. We reconsider the uncertainties of the method, concerning the mass profile of the lens galaxies, and show that the formal precision relies on the assumption that the mass profile is a perfect power law. Simple analytical arguments and numerical experiments reveal that mass-sheet like transformations yield significant freedom in choosing the mass profile, even when exquisite Einstein rings are observed. Furthermore, the characterization of the environment of the lens does not break that degeneracy which is not physically linked to extrinsic convergence. We present an illustrative example where the multiple imaging properties of a composite (baryons + dark matter) lens can be extremely well reproduced by a power-law model having the same velocity dispersion, but with predictions for the Hubble constant that deviate by $\sim 20%$. Hence we conclude that the impact of degeneracies between parametrized models have been underestimated in current $H_0$ measurements from lensing, and need to be carefully reconsidered.

preprint2013arXiv

Source-position transformation -- an approximate invariance in strong gravitational lensing

The main obstacle for gravitational lensing to determine accurate masses of deflectors, or to determine precise estimates for the Hubble constant, is the degeneracy of lensing observables with respect to the mass-sheet transformation (MST). The MST is a global modification of the mass distribution which leaves all image positions, shapes and flux ratios invariant, but which changes the time delay. Here we show that another global transformation of lensing mass distributions exists which almost leaves image positions and flux ratios invariant, and of which the MST is a special case. Whereas for axi-symmetric lenses this source position transformation exactly reproduces all strong lensing observables, it does so only approximately for more general lens situations. We provide crude estimates for the accuracy with which the transformed mass distribution can reproduce the same image positions as the original lens model, and present an illustrative example of its performance. This new invariance transformation most likely is the reason why the same strong lensing information can be accounted for with rather different mass models.

preprint2013arXiv

The cosmological information of shear peaks: beyond the abundance

We study the cosmological information of weak lensing (WL) peaks, focusing on two other statistics besides their abundance: the stacked tangential-shear profiles and the peak-peak correlation function. We use a large ensemble of simulated WL maps with survey specifications relevant to future missions like Euclid and LSST, to explore the three peak probes. We find that the correlation function of peaks with high signal-to-noise (S/N) measured from fields of size 144 sq. deg. has a maximum of ~0.3 at an angular scale ~10 arcmin. For peaks with smaller S/N, the amplitude of the correlation function decreases, and its maximum occurs on smaller angular scales. We compare the peak observables measured with and without shape noise and find that for S/N~3 only ~5% of the peaks are due to large-scale structures, the rest being generated by shape noise. The covariance matrix of the probes is examined: the correlation function is only weakly covariant on scales < 30 arcmin, and slightly more on larger scales; the shear profiles are very correlated for theta > 2 arcmin, with a correlation coefficient as high as 0.7. Using the Fisher-matrix formalism, we compute the cosmological constraints for {Om_m, sig_8, w, n_s} considering each probe separately, as well as in combination. We find that the correlation function of peaks and shear profiles yield marginalized errors which are larger by a factor of 2-4 for {Om_m, sig_8} than the errors yielded by the peak abundance alone, while the errors for {w, n_s} are similar. By combining the three probes, the marginalized constraints are tightened by a factor of ~2 compared to the peak abundance alone, the least contributor to the error reduction being the correlation function. This work therefore recommends that future WL surveys use shear peaks beyond their abundance in order to constrain the cosmological model.

preprint2013arXiv

Towards an understanding of third-order galaxy-galaxy lensing

Third-order galaxy-galaxy lensing (G3L) is a next generation galaxy-galaxy lensing technique that either measures the excess shear about lens pairs or the excess shear-shear correlations about lenses. It is clear that these statistics assess the three-point correlations between galaxy positions and projected matter density. For future applications of these novel statistics, we aim at a more intuitive understanding of G3L to isolate the main features that possibly can be measured. We construct a toy model ("isolated lens model"; ILM) for the distribution of galaxies and associated matter to determine the measured quantities of the two G3L correlation functions and traditional galaxy-galaxy lensing (GGL) in a simplified context. The ILM presumes single lens galaxies to be embedded inside arbitrary matter haloes that, however, are statistically independent ("isolated") from any other halo or lens position. In the ILM, the average mass-to-galaxy number ratio of clusters of any size cannot change. GGL and galaxy clustering alone cannot distinguish an ILM from any more complex scenario. The lens-lens-shear correlator in combination with second-order statistics enables us to detect deviations from a ILM, though. This can be quantified by a difference signal defined in the paper. We demonstrate with the ILM that this correlator picks up the excess matter distribution about galaxy pairs inside clusters. The lens-shear-shear correlator is sensitive to variations among matter haloes. In principle, it could be devised to constrain the ellipticities of haloes, without the need for luminous tracers, or maybe even random halo substructure. [Abridged]

preprint2012arXiv

A new third-order cosmic shear statistics: Separating E/B-mode correlations on a finite interval

Decomposing the shear signal into E and B-modes properly, i.e. without leakage of B-modes into the E-mode signal and vice versa, has been a long-standing problem in weak gravitational lensing. At the two-point level this problem was resolved by developing the so-called ring statistics, and later the COSEBIs; however, extending these concepts to the three-point level is far from trivial. Currently used methods to decompose three-point shear correlation functions (3PCFs) into E- and B-modes require knowledge of the 3PCF down to arbitrary small scales. This implies that the 3PCF needs to be modeled on scales smaller than the minimum separation of 2 galaxies and subsequently will be biased towards the model, or, in the absence of a model, the statistics is affected by E/B-mode leakage (or mixing). In this paper we derive a new third-order E/B-mode statistic that performs the decomposition using the 3PCF only on a finite interval, and thereby is free of any E/B-mode leakage while at the same time relying solely on information from the data. In addition, we relate this third-order ring statistics to the convergence field, thereby enabling a fast and convenient calculation of this statistic from numerical simulations. We note that our new statistics should be applicable to corresponding E/B-mode separation problems in the CMB polarization field.

preprint2012arXiv

Cosmic Shear Tomography and Efficient Data Compression using COSEBIs

Context. Gravitational lensing is one of the leading tools in understanding the dark side of the Universe. The need for accurate, efficient and effective methods which are able to extract this information along with other cosmological parameters from cosmic shear data is ever growing. COSEBIs, Complete Orthogonal Sets of E-/B-Integrals, is a recently developed statistical measure that encompasses the complete E-/B-mode separable information contained in the shear correlation functions measured on a finite angular range. Aims. The aim of the present work is to test the properties of this newly developed statistics for a higher-dimensional parameter space and to generalize and test it for shear tomography. Methods. We use Fisher analysis to study the effectiveness of COSEBIs. We show our results in terms of figure-of-merit quantities, based on Fisher matrices. Results. We find that a relatively small number of COSEBIs modes is always enough to saturate to the maximum information level. This number is always smaller for 'logarithmic COSEBIs' than for 'linear COSEBIs', and also depends on the number of redshift bins, the number and choice of cosmological parameters, as well as the survey characteristics. Conclusions. COSEBIs provide a very compact way of analyzing cosmic shear data, i.e., all the E-/B-mode separable second-order statistical information in the data is reduced to a small number of COSEBIs modes. Furthermore, with this method the arbitrariness in data binning is no longer an issue since the COSEBIs modes are discrete. Finally, the small number of modes also implies that covariances, and their inverse, are much more conveniently obtainable, e.g., from numerical simulations, than for the shear correlation functions themselves.

preprint2012arXiv

From étale $P_{+}$-representations to $G$-equivariant sheaves on $G/P$

Let $K/\mathbb Q_{p}$ be a finite extension with ring of integers $o$, let $G$ be a connected reductive split $\mathbb Q_{p}$-group of Borel subgroup $P=TN$ and let $α$ be a simple root of $T$ in $N$. We associate to a finitely generated module $D$ over the Fontaine ring over $o $ endowed with a semilinear étale action of the monoid $T_{+} $ (acting on the Fontaine ring via $α$), a $G(\mathbb Q_{p})$-equivariant sheaf of $o$-modules on the compact space $G(\mathbb Q_{p})/P(\mathbb Q_{p})$. Our construction generalizes the representation $D\boxtimes \mathbb P^{1} $ of $ GL(2,\mathbb Q_{p})$ associated by Colmez to a $(φ,Γ)$-module $D$ endowed with a character of $\mathbb Q_{p}^{*}$.

preprint2012arXiv

Measuring primordial non-Gaussianity through weak lensing peak counts

We explore the possibility of detecting primordial non-Gaussianity of the local type using weak lensing peak counts. We measure the peak abundance in sets of simulated weak lensing maps corresponding to three models f_NL={0, +100, -100}. Using survey specifications similar to those of Euclid and without assuming any knowledge of the lens and source redshifts, we find the peak functions of the non-Gaussian models with f_NL=+-100 to differ by up to 15% from the Gaussian peak function at the high-mass end. For the assumed survey parameters, the probability of fitting an f_NL=0 peak function to the f_NL=+-100 peak functions is less than 0.1%. Assuming the other cosmological parameters known, f_NL can be measured with an error Δf_NL ~ 13. It is therefore possible that future weak lensing surveys like Euclid and LSST may detect primordial non-Gaussianity from the abundance of peak counts, and provide complementary information to that obtained from the cosmic microwave background.

preprint2012arXiv

Optimized detection of shear peaks in weak lensing maps

We present a new method to extract cosmological constraints from weak lensing (WL) peak counts, which we denote as `the hierarchical algorithm'. The idea of this method is to combine information from WL maps sequentially smoothed with a series of filters of different size, from the largest down to the smallest, thus increasing the cosmological sensitivity of the resulting peak function. We compare the cosmological constraints resulting from the peak abundance measured in this way and the abundance obtained by using a filter of fixed size, which is the standard practice in WL peak studies. For this purpose, we employ a large set of WL maps generated by ray-tracing through N-body simulations, and the Fisher matrix formalism. We find that if low-S/N peaks are included in the analysis (S/N ~ 3), the hierarchical method yields constraints significantly better than the single-sized filtering. For a large future survey such as Euclid or LSST, combined with information from a CMB experiment like Planck, the results for the hierarchical (single-sized) method are: Δn=0.0039 (0.004); ΔΩm=0.002 (0.0045); Δσ8=0.003 (0.006); Δw=0.019 (0.0525). This forecast is conservative, as we assume no knowledge of the redshifts of the lenses, and consider a single broad bin for the redshifts of the sources. If only peaks with S/N >= 6 are considered, then there is little difference between the results of the two methods. We also examine the statistical properties of the hierarchical peak function: Its covariance matrix has off-diagonal terms for bins with S/N <= 6 and aperture mass of M < 3 x 1e+14 Ms/h, the higher bins being largely uncorrelated and therefore well described by a Poisson distribution.

preprint2012arXiv

Pro-p Iwahori-Hecke algebras are Gorenstein

Let F be a locally compact nonarchimedean field with residue characteristic p and G the group of F-rational points of a connected split reductive group over F. For k an arbitrary field, we study the homological properties of the Iwahori-Hecke k-algebra H' and of the pro-p Iwahori-Hecke k-algebra H of G. We prove that both of these algebras are Gorenstein rings with self-injective dimension bounded above by the rank of G. If G is semisimple, we also show that this upper bound is sharp, that both H and H' are Auslander-Gorenstein and that there is a duality functor on the finitely generated modules of H (respectively H'). We obtain the analogous Gorenstein and Auslander-Gorenstein properties for the graded rings associated to H and H'. When k has characteristic p, we prove that in most cases H and H' have infinite global dimension. In particular, we deduce that the category of smooth k-representations of G=PGL(2,Q_p) generated by their invariant vectors under the pro-p-Iwahori subgroup has infinite global dimension (at least if k is algebraically closed).

preprint2012arXiv

Probing cluster dynamics in RXCJ1504.1-0248 via radial and two-dimensional gas and galaxy properties

We studied one of the most X-ray luminous cluster of galaxies in the REFLEX survey, RXC J1504.1-0248 (hereafter R1504; z=0.2153), using XMM-Newton X-ray imaging spectroscopy, VLT/VIMOS optical spectroscopy, and WFI optical imaging. The mass distributions were determined using both the so-called hydrostatic method with X-ray imaging spectroscopy and the dynamical method with optical spectroscopy, respectively, which yield M^{H.E.}_{500}=(5.81+/-0.49)*1.e14Msun and M^{caustic}_{500}=(4.17+/-0.42)*1e14Msun. According to recent calibrations, the richness derived mass estimates closely agree with the hydrostatic and dynamical mass estimates. The line-of-sight velocities of spectroscopic members reveal a group of galaxies with high-velocities (>1000 km/s) at a projected distance of about r^{H.E.}_{500}=(1.18+/-0.03) Mpc south-east of the cluster centroid, which is also indicated in the X-ray two-dimensional (2-D) temperature, density, entropy, and pressure maps. The dynamical mass estimate is 80% of the hydrostatic mass estimate at r^{H.E.}_{500}. It can be partially explained by the ~20% scatter in the 2-D pressure map that can be propagated into the hydrostatic mass estimate. The uncertainty in the dynamical mass estimate caused by the substructure of the high velocity group is ~14%. The dynamical mass estimate using blue members is 1.23 times that using red members. The global properties of R1504 obey the observed scaling relations of nearby clusters, although its stellar-mass fraction is rather low.

preprint2012arXiv

The 400d Galaxy Cluster Survey weak lensing programme: II: Weak lensing study of seven clusters with MMT/Megacam

Evolution in the mass function of galaxy clusters sensitively traces both the expansion history of the Universe and cosmological structure formation. Robust cluster mass determinations are a key ingredient for a reliable measurement of this evolution, especially at high redshift. Weak gravitational lensing is a promising tool for, on average, unbiased mass estimates. This weak lensing project aims at measuring reliable weak lensing masses for a complete X-ray selected sample of 36 high redshift (0.35<z<0.9) clusters. The goal of this paper is to demonstrate the robustness of the methodology against commonly encountered problems, including pure instrumental effects, the presence of bright (8--9 mag) stars close to the cluster centre, ground based measurements of high-z (z~0.8) clusters, and the presence of massive unrelated structures along the line-sight. We select a subsample of seven clusters observed with MMT/Megacam. Instrumental effects are checked in detail by cross-comparison with an archival CFHT/MegaCam observation. We derive mass estimates for seven clusters by modelling the tangential shear with an NFW profile, in two cases with multiple components to account for projected structures in the line-of-sight. We firmly detect lensing signals from all seven clusters at more than $3.5σ$ and determine their masses, ranging from $10^{14} M_{\odot}$ to $10^{15} M_{\odot}$, despite the presence of nearby bright stars. We retrieve the lensing signal of more than one cluster in the CL 1701+6414 field, while apparently observing CL 1701+6414 through a massive foreground filament. We also find a multi-peaked shear signal in CL 1641+4001. Shear structures measured in the MMT and CFHT images of CL 1701+6414 are highly correlated.

preprint2011arXiv

A bias in cosmic shear from galaxy selection: results from ray-tracing simulations

We identify and study a previously unknown systematic effect on cosmic shear measurements, caused by the selection of galaxies used for shape measurement, in particular the rejection of close (blended) galaxy pairs. We use ray-tracing simulations based on the Millennium Simulation and a semi-analytical model of galaxy formation to create realistic galaxy catalogues. From these, we quantify the bias in the shear correlation functions by comparing measurements made from galaxy catalogues with and without removal of close pairs. A likelihood analysis is used to quantify the resulting shift in estimates of cosmological parameters. The filtering of objects with close neighbours (a) changes the redshift distribution of the galaxies used for correlation function measurements, and (b) correlates the number density of sources in the background with the density field in the foreground. This leads to a scale-dependent bias of the correlation function of several percent, translating into biases of cosmological parameters of similar amplitude. This makes this new systematic effect potentially harmful for upcoming and planned cosmic shear surveys. As a remedy, we propose and test a weighting scheme that can significantly reduce the bias.

preprint2011arXiv

Constrained probability distributions of correlation functions

Context: Two-point correlation functions are used throughout cosmology as a measure for the statistics of random fields. When used in Bayesian parameter estimation, their likelihood function is usually replaced by a Gaussian approximation. However, this has been shown to be insufficient. Aims: For the case of Gaussian random fields, we search for an exact probability distribution of correlation functions, which could improve the accuracy of future data analyses. Methods: We use a fully analytic approach, first expanding the random field in its Fourier modes, and then calculating the characteristic function. Finally, we derive the probability distribution function using integration by residues. We use a numerical implementation of the full analytic formula to discuss the behaviour of this function. Results: We derive the univariate and bivariate probability distribution function of the correlation functions of a Gaussian random field, and outline how higher joint distributions could be calculated. We give the results in the form of mode expansions, but in one special case we also find a closed-form expression. We calculate the moments of the distribution and, in the univariate case, we discuss the Edgeworth expansion approximation. We also comment on the difficulties in a fast and exact numerical implementation of our results, and on possible future applications.

preprint2011arXiv

Cosmic shear covariance: The log-normal approximation

[Abridged] We seek approximations to the cosmic shear covariance that are as easy to use as the common approximations based on normal statistics, but yield more accurate covariance matrices and parameter errors. We derive expressions for the cosmic shear covariance under the assumption that the underlying convergence field follows log-normal statistics. We also derive a simplified version of this log-normal approximation. We use numerical simulations of weak lensing to study how well the normal, log-normal, and simplified log-normal approximations as well as empirical corrections to the normal approximation proposed in the literature reproduce shear covariances for cosmic shear surveys. We find that the normal approximation substantially underestimates the cosmic shear covariances and the inferred parameter confidence regions, in particular for surveys with small fields of view and large galaxy densities, but also for very wide surveys. In contrast, the log-normal approximation yields more realistic covariances and confidence regions, but also requires evaluating slightly more complicated expressions. However, the simplified log-normal approximation, although as simple as the normal approximation, yields confidence regions that are almost as accurate as those obtained from the log-normal approximation. The empirical corrections to the normal approximation do not yield more accurate covariances and confidence regions than the (simplified) log-normal approximation. Moreover, they fail to produce positive-semidefinite data covariance matrices in certain cases, rendering them unusable for parameter estimation. The simplified log-normal approximation should be used in favour of the normal approximation for parameter estimation and parameter error forecasts. Any approximation to the cosmic shear covariance should ensure a positive-semidefinite data covariance matrix.

preprint2011arXiv

K_1 of certain Iwasawa algebras, after Kakde

This paper contains a detailed exposition of the content of section five in Kakde's paper arXiv:1008.0142. We proceed in a slightly more axiomatic way to pin down the exact requirements on the $p$-adic Lie group under consideration. We also make use of our conceptual theory of the completed localization of an Iwasawa algebra as developed in arXiv:0711.2669. This simplifies some of the arguments. Otherwise, with the exception of the notation at certain places, we follow Kakde's paper.

preprint2011arXiv

Star-formation efficiency and metal enrichment of the intracluster medium in local massive clusters of galaxies

We have investigated the baryon-mass content in a subsample of 19 clusters of galaxies extracted from the X-ray flux-limited sample HIFLUGCS according to their positions in the sky. For these clusters, we measured total masses and characteristic radii on the basis of a rich optical spectroscopic data set, the physical properties of the intracluster medium (ICM) using XMM-Newton and ROSAT X-ray data, and total (galaxy) stellar masses utilizing the SDSS DR7 multi-band imaging. The observed (hot) gas-mass fractions are almost constant in this mass range. We confirm that the stellar mass fraction decreases as the total mass increases and shows (20+/-4)% scatter; in addition, we show that it decreases as the central entropy increases. The latter behavior supports a twofold interpretation where heating from merging quenches the star-formation activity of galaxies in massive systems, and feedback from supernovae and/or radio galaxies drives a significant amount of gas to the regions beyond r_{500} or, alternatively, a substantially large amount of intracluster light (ICL) is associated with galaxies in nonrelaxed systems. Furthermore, less massive clusters are confirmed to host less gas per unit total mass; however, they exhibit higher mass fractions in metals, so that their ICM is more metal-rich. This again supports the interpretation that in the potential wells of low-mass systems the star-formation efficiency of galaxies was high or, alternatively, some gas is missing from the hot phase of the ICM. The former hypothesis is preferred as the main driver of the mass-dependent metal enrichment since the total mass-to-optical luminosity ratio increases as the total mass increases.

preprint2011arXiv

The bispectrum covariance beyond Gaussianity: A log-normal approach

To investigate and specify the statistical properties of cosmological fields with particular attention to possible non-Gaussian features, accurate formulae for the bispectrum and the bispectrum covariance are required. The bispectrum is the lowest-order statistic providing an estimate for non-Gaussianities of a distribution, and the bispectrum covariance depicts the errors of the bispectrum measurement and their correlation on different scales. Currently, there do exist fitting formulae for the bispectrum and an analytical expression for the bispectrum covariance, but the former is not very accurate and the latter contains several intricate terms and only one of them can be readily evaluated from the power spectrum of the studied field. Neglecting all higher-order terms results in the Gaussian approximation of the bispectrum covariance. We study the range of validity of this Gaussian approximation for two-dimensional non-Gaussian random fields. For this purpose, we simulate Gaussian and non-Gaussian random fields, the latter represented by log-normal fields and obtained directly from the former by a simple transformation. From the simulated fields, we calculate the power spectra, the bispectra, and the covariance from the sample variance of the bispectra, for different degrees of non-Gaussianity α, which is equivalent to the skewness on a given angular scale θg. We find that the Gaussian approximation provides a good approximation for α<0.6 and a reasonably accurate approximation for α< 1, both on scales >8θg. Using results from cosmic shear simulations, we estimate that the cosmic shear convergence fields are described by α<0.7 at θg~4". We therefore conclude that the Gaussian approximation for the bispectrum covariance is likely to be applicable in ongoing and future cosmic shear studies.

preprint2011arXiv

XMM-Newton/SDSS: star formation efficiency in galaxy clusters and constraints on the matter density parameter

It is believed that the global baryon content of clusters of galaxies is representative of the matter distribution of the universe, and can, therefore, be used to reliably determine the matter density parameter Omega_m. This assumption is challenged by the growing evidence from optical and X-ray observations that the total baryon mass fraction increases towards rich clusters. In this context, we investigate the dependence of stellar, and total baryon mass fractions as a function of mass. To do so, we used a subsample of nineteen clusters extracted from the X-ray flux limited sample HIFLUGCS that have available DR-7 Sloan Digital Sky Survey (SDSS) data. From the optical analysis we derived the stellar masses. Using XMM-Newton we derived the gas masses. Then, adopting a scaling relation we estimate the total masses. Adding the gas and the stellar mass fractions we obtain the total baryonic content that we find to increase with cluster mass, reaching 7-year Wilkinson Microwave Anisotropy Probe (WMAP-7) prediction for clusters with M_500 = 1.6 x 10^{15} M_sun. We observe a decrease of the stellar mass fraction (from 4.5% to ~1.0%) with increasing total mass where our findings for the stellar mass fraction agree with previous studies. This result suggests a difference in the number of stars formed per unit of halo mass, though with a large scatter for low-mass systems. That is, the efficiency of star formation varies on cluster scale that lower mass systems are likely to have higher star formation efficiencies. It follows immediately that the dependence of the stellar mass fraction on total mass results in an increase of the mass-to-light ratio from lower to higher mass systems. We also discuss the consequences of these results in the context of determining the cosmic matter density parameter Omega_m.

preprint2010arXiv

COSEBIs: Extracting the full E-/B-mode information from cosmic shear correlation functions

Cosmic shear is considered one of the most powerful methods for studying the properties of Dark Energy in the Universe. As a standard method, the two-point correlation functions $xi_\pm(theta)$ of the cosmic shear field are used as statistical measures for the shear field. In order to separate the observed shear into E- and B-modes, the latter being most likely produced by remaining systematics in the data set and/or intrinsic alignment effects, several statistics have been defined before. Here we aim at a complete E-/B-mode decomposition of the cosmic shear information contained in the $xi_\pm$ on a finite angular interval. We construct two sets of such E-/B-mode measures, namely Complete Orthogonal Sets of E-/B-mode Integrals (COSEBIs), characterized by weight functions between the $xi_\pm$ and the COSEBIs which are polynomials in $theta$ or polynomials in $ln(theta)$, respectively. Considering the likelihood in cosmological parameter space, constructed from the COSEBIs, we study their information contents. We show that the information grows with the number of COSEBI modes taken into account, and that an asymptotic limit is reached which defines the maximum available information in the E-mode component of the $xi_\pm$. We show that this limit is reached the earlier (i.e., for a smaller number of modes considered) the narrower the angular range is over which $xi_\pm$ are measured, and it is reached much earlier for logarithmic weight functions. For example, for $xi_\pm$ on the interval $1'\le θ\le 400'$, the asymptotic limit for the parameter pair $(Omega_m, sigma_8)$ is reached for $\sim 25$ modes in the linear case, but already for 5 modes in the logarithmic case. The COSEBIs form a natural discrete set of quantities, which we suggest as method of choice in future cosmic shear likelihood analyses.

preprint2010arXiv

Estimate of dark halo ellipticity by lensing flexion

Aims. The predictions of the ellipticity of the dark matter halos from models of structure formation are notoriously difficult to test with observations. A direct measurement would give important constraints on the formation of galaxies, and its effect on the dark matter distribution in their halos. Here we show that galaxy-galaxy flexion provides a direct and potentially powerful method for determining the ellipticity of (an ensemble of) elliptical lenses. Methods. We decompose the spin-1 flexion into a radial and a tangential component. Using the ratio of tangential-to- radial flexion, which is independent of the radial mass profile, the mass ellipticity can be estimated. Results. An estimator for the ellipticity of the mass distribution is derived and tested with simulations. We show that the estimator is slightly biased. We quantify this bias, and provide a method to reduce it. Furthermore, a parametric fitting of the flexion ratio and orientation provides another estimate for the dark halo ellipticity, which is more accurate for individual lenses Overall, galaxy-galaxy flexion appears as a powerful tool for constraining the ellipticity of mass distributions.

preprint2010arXiv

Evidence for the accelerated expansion of the Universe from weak lensing tomography with COSMOS

We present a tomographic cosmological weak lensing analysis of the HST COSMOS Survey. Applying our lensing-optimized data reduction, principal component interpolation for the ACS PSF, and improved modelling of charge-transfer inefficiency, we measure a lensing signal which is consistent with pure gravitational modes and no significant shape systematics. We carefully estimate the statistical uncertainty from simulated COSMOS-like fields obtained from ray-tracing through the Millennium Simulation. We test our pipeline on simulated space-based data, recalibrate non-linear power spectrum corrections using the ray-tracing, employ photometric redshifts to reduce potential contamination by intrinsic galaxy alignments, and marginalize over systematic uncertainties. We find that the lensing signal scales with redshift as expected from General Relativity for a concordance LCDM cosmology, including the full cross-correlations between different redshift bins. For a flat LCDM cosmology, we measure sigma_8(Omega_m/0.3)^0.51=0.75+-0.08 from lensing, in perfect agreement with WMAP-5, yielding joint constraints Omega_m=0.266+0.025-0.023, sigma_8=0.802+0.028-0.029 (all 68% conf.). Dropping the assumption of flatness and using HST Key Project and BBN priors only, we find a negative deceleration parameter q_0 at 94.3% conf. from the tomographic lensing analysis, providing independent evidence for the accelerated expansion of the Universe. For a flat wCDM cosmology and prior w in [-2,0], we obtain w<-0.41 (90% conf.). Our dark energy constraints are still relatively weak solely due to the limited area of COSMOS. However, they provide an important demonstration for the usefulness of tomographic weak lensing measurements from space. (abridged)

preprint2010arXiv

Mass reconstruction by gravitational shear and flexion

Galaxy clusters are considered as excellent probes for cosmology. For that purpose, their mass needs to be measured and their structural properties needs to be understood. We propose a method for galaxy cluster mass reconstruction which combines information from strong lensing, weak lensing shear and flexion. We extend the weak lensing analysis to the inner parts of the cluster and, in particular, improve the resolution of substructure. We use simulations to show that the method recovers the mass density profiles of the cluster. We find that the weak lensing flexion is sensitive to substructure. After combining the flexion data into the joint weak and strong lensing analysis, we can resolve the cluster properties with substructures.

preprint2010arXiv

The 400d Galaxy Cluster Survey Weak Lensing Programme: I: MMT/Megacam Analysis of CL0030+2618 at z=0.50

The mass function of galaxy clusters at high redshifts is a particularly useful probe to learn about the history of structure formation and constrain cosmological parameters. We aim at deriving reliable masses for a high-redshift, high-luminosity sample of clusters of galaxies selected from the 400d survey of X-ray selected clusters. Here, we will focus on a particular object, CL0030+2618 at z=0.50 Using deep imaging in three passbands with the MEGACAM instrument at MMT, we show that MEGACAM is well-suited for measuring gravitational shear. We detect the weak lensing signal of CL0030+2618 at 5.8 sigma significance, using the aperture mass technique. Furthermore, we find significant tangential alignment of galaxies out to ~10 arcmin or >2r_200 distance from the cluster centre. The weak lensing centre of CL0030+2618 agrees with several X-ray measurements and the position of the brightest cluster galaxy. Finally, we infer a weak lensing virial mass of M_200=7.5 10^{14} M_sun for CL0030+2618. Despite complications by a tentative foreground galaxy group in the line of sight, the X-ray and weak lensing estimates for CL0030+2618 are in remarkable agreement. This study paves the way for the largest weak lensing survey of high-redshift galaxy clusters to date.

preprint2009arXiv

Constrained correlation functions

We show that correlation functions have to satisfy contraint relations, owing to the non-negativity of the power spectrum of the underlying random process. Specifically, for any statistically homogeneous and (for more than one spatial dimension) isotropic random field with correlation function $ξ(x)$, we derive inequalities for the correlation coefficients $r_n\equiv ξ(n x)/ξ(0)$ (for integer $n$) of the form $r_{n{\rm l}}\le r_n\le r_{n{\rm u}}$, where the lower and upper bounds on $r_n$ depend on the $r_j$, with $j<n$. Explicit expressions for the bounds are obtained for arbitrary $n$. These constraint equations very significantly limit the set of possible correlation functions. For one particular example of a fiducial cosmic shear survey, we show that the Gaussian likelihood ellipsoid has a significant spill-over into the forbidden region of correlation functions, rendering the resulting best-fitting model parameters and their error region questionable, and indicating the need for a better description of the likelihood function. We conduct some simple numerical experiments which explicitly demonstrate the failure of a Gaussian description for the likelihood of $ξ$. Instead, the shape of the likelihood function of the correlation coefficients appears to follow approximately that of the shape of the bounds on the $r_n$, even if the Gaussian ellipsoid lies well within the allowed region. For more than one spatial dimension of the random field, the explicit expressions of the bounds on the $r_n$ are not optimal. We outline a geometrical method how tighter bounds may be obtained in principle. We illustrate this method for a few simple cases; a more general treatment awaits future work.

preprint2009arXiv

Measuring cosmic shear with the ring statistics

Commonly used methods to decompose E- and B-modes in cosmic shear, namely the aperture mass dispersion and the E/B-mode shear correlation function, suffer from incomplete knowledge of the two-point correlation function (2PCF) on very small and/or very large scales. The ring statistics, the most recently developed cosmic shear measure, improves on this issue and is able to decompose E- and B-modes using a 2PCF measured on a finite interval. First, we improve on the ring statistics' filter function with respect to the signal-to-noise ratio. Second, we examine the ability of the ring statistics to constrain cosmology and compare the results to cosmological constraints obtained with the aperture mass dispersion. Third, we use the ring statistics to measure a cosmic shear signal from CFHTLS (Canada-France-Hawaii Telescope Legacy Survey) data. We consider a scale-dependent filter function for the ring statistics which improves its signal-to-noise ratio. In addition, we show that there exist filter functions which decompose E- and B-modes using a finite range of 2PCFs (EB-statistics) and have higher S/N ratio than the ring statistics. However, we find that data points of the latter are significantly less correlated than data points of the aperture mass dispersion and the EB-statistics. As a consequence the ring statistics is an ideal tool to identify remaining systematics accurately as a function of angular scale. We use the 2PCF of the latest CFHTLS analysis and therefrom calculate the ring statistics and its error bars.

preprint2008arXiv

Dependence of cosmic shear covariances on cosmology - Impact on parameter estimation

In cosmic shear likelihood analyses the covariance is most commonly assumed to be constant in parameter space. Therefore, when calculating the covariance matrix (analytically or from simulations), its underlying cosmology should not influence the likelihood contours. We examine whether the aforementioned assumption holds and quantify how strong cosmic shear covariances vary within a reasonable parameter range. Furthermore, we examine the impact on likelihood contours when assuming different cosmologies in the covariance. We find that covariances vary significantly within the considered parameter range (Omega_m=[0.2;0.4], sigma_8=[0.6;1.0]) and that this has a non-negligible impact on the size of likelihood contours. This impact increases with increasing survey size, increasing number density of source galaxies, decreasing ellipticity noise, and when using non-Gaussian covariances. To improve on the assumption of a constant covariance we present two methods. The adaptive covariance is the most accurate method, but it is computationally expensive. To reduce the computational costs we give a scaling relation for covariances. As a second method we outline the concept of an iterative likelihood analysis. Here, we additionally account for non-Gaussianity using a ray-tracing covariance derived from the Millennium simulation.

preprint2007arXiv

Localisations and Completions of Skew Power Series Rings

This paper is a natural continuation of the study of skew power series rings A initiated in [P. Schneider and O. Venjakob, On the codimension of modules over skew power series rings with applications to Iwasawa algebras, J. Pure Appl. Algebra 204 (2005), 349 - 367.]. We construct skew Laurent series rings B and show the existence of some canonical Ore sets S for the skew power series rings A such that a certain completion of the localisation A_S is isomorphic to B. This is applied to certain Iwasawa algebras. Finally we introduce subrings of overconvergent skew Laurent series rings.

preprint2006arXiv

E- and B-mode mixing from incomplete knowledge of the shear correlation

Aims. We quantify the mixing of the measured cosmic-shear E- and B-modes caused by the lack of shear-correlation measurements on small and large scales, arising from a lack of close projected galaxy pairs and the finite field size, respectively. Methods. We calculate the aperture-mass statistics <M_{ap, perp}^2> and the E-/B-mode shear-correlation functions xi_{E, B +/-} where small- and large-scale cutoffs are taken into account. We assess the deviation of the obtained E-mode to the true E-mode and the introduction of a spurious B-mode. Results. The measured aperture-mass dispersion is underestimated by more than 10% on scales smaller than 12 times the lower cutoff. For a precise measurement of the E- and B-modes at the percent level using a combination of xi_{E, B +} and xi_{E, B -}, a field as large as 7 (2.4) degrees is necessary for ground-based (space-based) observations.

preprint2006arXiv

The ring statistics - how to separate E- and B-modes of cosmic shear correlation functions on a finite interval

Aims. One of the main probes for systematic errors in the cosmic shear signal are the division of the shear field into E- and B-mode shear, where gravitational lensing only produces the former. As shown in a recent note, all currently used E-/B-mode separation methods for the shear correlation functions xi_pm require them to be measured to arbitrarily small and/or large separations which is of course not feasible in practice. Methods. We derive second-order shear statistics which provide a clean separation into E- and B-modes from measurements of xi_pm(theta) over a finite interval only. We call these new statistics the circle and ring statistics, respectively; the latter is obtained by an integral over the former. The mathematical properties of these new shear statistics are obtained, as well as specific expressions for applying them to observed data. Results. It is shown that an E-/B-mode separation can be performed on measurements of xi_pm over a finite interval in angular separation, using the ring statistics. We furthermore generalize this result to derive the most general class of second-order shear statistics which provide a separation of E- and B-mode shear on a given angular interval theta_min <= theta <= theta_max. Our results will be of practical use particularly for future cosmic shear surveys where highly precise measurements of the shear will become available and where control of systematics will be mandatory.

preprint2004arXiv

On the dimension theory of skew power series rings

The first purpose of this paper is to set up a general notion of skew power series rings S over a coefficient ring R, which are then studied by filtered ring techniques. The second subject consists of investigating the class of S-modules which are finitely generated as R-module. In the case that S and R are Auslander regular we show in particular that the codimension of M as S-module is one higher than the codimension of M as R-module. Furthermore its class in the Grothendieck group of S-modules of codimension at most one less vanishes, which is in the spirit of the Gersten conjecture for commutative regular local rings. Finally we apply these results to Iwasawa algebras of p-adic Lie groups.

preprint2000arXiv

Cosmic shear and halo abundances: Analytical versus numerical results

The aperture mass has been shown in a series of recent publications to be a useful quantitative tool for weak lensing studies, ranging from cosmic shear to the detection of a mass-selected sample of dark matter haloes. Quantitative analytical predictions for the aperture mass have been based on a number of simplifying assumptions. In this paper, we test the reliability of these assumptions and the quality of the analytic approximations, using ray-tracing simulations through a cosmological density field generated by very large N-body simulations. We find that those analytic predictions which take into account the non-linear evolution of the matter distribution, such as the dispersion of the aperture mass and the halo abundance, are surprisingly accurately reproduced with our numerical results, whereas the predictions for the skewness, based on quasi-linear theory, are rather imprecise. In particular, we verify numerically that the probability distribution of the aperture mass decreases exponentially for values much larger than the rms. Given the good overall agreement, comparisons between the observed distribution of the aperture mass and the theoretical values provide a powerful tool for testing cosmological models.

preprint1999arXiv

The non-Gaussian tail of cosmic-shear statistics

Due to gravitational instability, an initially Gaussian density field develops non-Gaussian features as the Universe evolves. The most prominent non-Gaussian features are massive haloes, visible as clusters of galaxies. The distortion of high-redshift galaxy images due to the tidal gravitational field of the large-scale matter distribution, called cosmic shear, can be used to investigate the statistical properties of the LSS. In particular, non-Gaussian properties of the LSS will lead to a non-Gaussian distribution of cosmic-shear statistics. The aperture mass ($M_{\rm ap}$) statistics, recently introduced as a measure for cosmic shear, is particularly well suited for measuring these non-Gaussian properties. In this paper we calculate the highly non-Gaussian tail of the aperture mass probability distribution, assuming Press-Schechter theory for the halo abundance and the `universal' density profile of haloes as obtained from numerical simulations. We find that for values of $M_{\rm ap}$ much larger than its dispersion, this probability distribution is closely approximated by an exponential, rather than a Gaussian. We determine the amplitude and shape of this exponential for various cosmological models and aperture sizes, and show that wide-field imaging surveys can be used to distinguish between some of the currently most popular cosmogonies. Our study here is complementary to earlier cosmic-shear investigations which focussed more on two-point statistical properties.

preprint1997arXiv

The cosmological lens equation and the equivalent single-plane gravitational lens

The gravitational lens equation resulting from a single (non-linear) mass concentration (the main lens) plus inhomogeneities of the large-scale structure is shown to be strictly equivalent to the single-plane gravitational lens equation without the cosmological perturbations. The deflection potential (and, by applying the Poisson equation, also the mass distribution) of the equivalent single-plane lens is derived. If the main lens is described by elliptical isopotential curves plus a shear term, the equivalent single-plane lens will be of the same form. Due to the equivalence shown, the determination of the Hubble constant from time delay measurements is affected by the same mass-sheet invariance transformation as for the single-plane lens. If the lens strength is fixed (e.g., by measuring the velocity dispersion of stars in the main lens), the determination of $H_0$ is affected by inhomogeneous matter between us and the lens. The orientation of the mass distribution relative to the image positions is the same for the cosmological lens situation and the single-plane case. In particular this implies that cosmic shear cannot account for a misalignment of the observed galaxy orientation relative to the best-fitting lens model.

preprint1996arXiv

Aperture Multipole Moments from Weak Gravitational Lensing

The projected mass of a gravitational lens inside (circular) apertures can be derived from the measured shear inside an annulus which is caused by the tidal field of the deflecting mass distribution. Here we show that also the multipoles of the two-dimensional mass distribution can be derived from the shear in annuli. We derive several expressions for these mass multipole moments in terms of the shear, which allow large flexibility in the choice of a radial weight function. In contrast to determining multipole moments from weak-lensing mass reconstructions, this approach allows to quantify the signal-to-noise ratio of the multipole moments directly from the observed galaxy ellipticities, and thus to estimate the significance of the multipole detection. Radial weight functions can therefore be chosen such as to optimize the significance of the detection given an assumed radial mass profile. Application of our formulae to numerically simulated clusters demonstrates that the quadrupole moment of realistic cluster models can be detected with high signal-to-noise ratio S/N; in about 85 per cent of the simulated cluster fields S/N >~ 3. We also show that the shear inside a circular annulus determines multipole moments inside and outside the annulus. This is relevant for clusters whose central region is too bright to allow the observation of the shear of background galaxies, or which extend beyond the CCD. We also generalize the aperture mass equation to the case of `radial' weight functions which are constant on arbitrarily-shaped curves which are not necessarily self-similar.

preprint1996arXiv

Deep radio observation of the gravitational lens candidate QSO2345+007

The double QSO2345+007 comprises two optical components separated by 7.1\ arcseconds and is the most prominent `dark matter' gravitational lens candidate. Despite being known for more than a decade, optical spectroscopy and imaging have been unable to determine whether this double QSO is a binary QSO or a gravitational lens system. In this note we report a deep VLA observation of this system, yielding a map with a noise level of 8.5\ts $μ$Jy per beam. We have a $4σ$ detection of a radio source within one arcsecond of the optical position of the brighter A-component of the QSO, but no significant detection of any radio counterpart of the B component. Given that the flux ratio in the optical waveband is $\sim$3--4, the gravitational lens hypothesis for this system would predict a radio flux of image B of $\lesssim 10$\ts $μ$Jy. The non-detection of the B component is thus consistent with, but does not prove, the lens interpretation.

preprint1996arXiv

Detection of (dark) matter concentrations via weak gravitational lensing

The distortion of images of faint background galaxies by (weak) gravitational lensing can be used to measure the mass distribution of the deflector. The image distortions can be used to define a weighted mean of the mass inside a circular aperture, as was first suggested by Kaiser. The aperture mass can be used to {\it detect} dark matter concentrations. Keeping in mind that wide-field cameras will become increasingly available, this method can be used to search for mass concentrations on wide-field images. To do this, the aperture mass measure is generalized to account for different weighting functions. For each such weighting function, a signal-to-noise ratio can be calculated. For an assumed mass profile of the density concentrations, the weighting function can be chosen such as to maximize the resulting signal-to-noise ratio. Numerical simulations which adopt parameters characteristic of 4-m class telescopes are then used to show that dark halos with a velocity dispersion in excess of $\sim 600$\ts km/s can be reliably detected as significant peaks in the signal-to-noise map. The effects of seeing and an anisotropic PSF are then investigated and shown to be less important than might be feared. It is thus suggested that the method of aperture mass measures developed here can be used to obtain a mass-selected sample of dark halos, in contrast to flux-selected samples. Shear fields around high-redshift bright QSOs as detected by Fort et al. provide a first successful application of this strategy. The simplicity of the method allows its routine application to wide-field images of sufficient depth and image quality.

preprint1996arXiv

The lightcurve reconstruction method for measuring the time delay of gravitational lens systems

We propose a new technique to measure the time delay of radio-loud gravitational lens systems, which does not rely on the excessive use of interferometric observations. Instead, the method is based on single-dish flux density monitoring of the (unresolved) lens system's total lightcurve, combined with additional interferometric measurements of the flux density ratio at a few epochs during that monitoring period. The basic idea of the method is to reconstruct the individual image lightcurves from the observed total lightcurve by assuming a range of potential values for the time delay and the magnification ratio of the images. It is then possible to single out the correct reconstruction, and therefore determine the time delay, by checking the consistency of the reconstructed individual lightcurves with the additional interferometric observations. We performed extensive numerical simulations of synthetic lightcurves to investigate the dependence of the performance of this method on various parameters which are involved in the problem. Probably the most promising candidates for applying the method (and also for determining the Hubble constant) are lens systems consisting of multiply imaged compact sources and an Einstein ring, such as B0218+357 from which some of the parameters used for our simulations were adopted.

preprint1995arXiv

Cosmological Applications of Gravitational Lensing

The last decade has seen an enormous increase of activity in the field of gravitational lensing, mainly driven by improvements of observational capabilities. I will review the basics of gravitational lens theory, just enough to understand the rest of this contribution, and will then concentrate on several of the main applications in cosmology. Cluster lensing, and weak lensing, will constitute the main part of this review.

preprint1994arXiv

Light Propagation in Arbitrary Spacetimes and the Gravitational Lens Approximation

We reformulate the transport equation which determines the size, shape and orientation of infinitesimal light beams in arbitrary spacetimes. The behaviour of such light beams near vertices and conjugate points is investigated, with special attention to the singular behaviour of the optical scalars. We then specialize the general transport equation to the case of an approximate metric of an inhomogeneous universe, which is a Friedmann metric `on average' with superposed isolated weak matter inhomogeneities. In a series of well-defined approximations, the equations of gravitational lens theory are derived. Finally, we derive a relative optical focusing equation which describes the focusing of light beams relative to the case that the beam is unaffected by matter inhomogeneities in the universe, from which it follows immediately that no beam can be focused less than one which is unaffected by matter clumps, before it propagates through its first conjugate point.