Source author record

Stefan Hilbert

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

34works
6topics
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

34 published item(s)

preprint2020arXiv

CFHTLenS: Galaxy bias as function of scale, stellar mass, and colour. Conflicts with predictions by semi-analytic models

Galaxy models predict a tight relation between the clustering of galaxies and dark matter on cosmological scales, but predictions differ notably in the details. We used this opportunity and tested two semi-analytic models by the Munich and Durham groups with data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). For the test we measured the scale-dependent galaxy bias factor $b(k)$ and correlation factor $r(k)$ from linear to non-linear scales of $k\approx10\,h\,\rm Mpc^{-1}$ at two redshifts $\bar{z}=0.35,0.51$ for galaxies with stellar mass between $5\times10^9$ and $3\times10^{11}\,h_{\rm 70}^{-2}\,{\rm M}_\odot$. Our improved gravitational lensing technique accounts for the intrinsic alignment of sources and the magnification of lens galaxies for better constraints for the galaxy-matter correlation $r(k)$. Galaxy bias in CFHTLenS increases with $k$ and stellar mass, it is colour-dependent, revealing the individual footprints of galaxy types. Despite a reasonable model agreement for the relative change with both scale and galaxy properties, there is a clear conflict for $b(k)$ with no model preference: the model galaxies are too weakly clustered. This may flag a model problem at $z\gtrsim0.3$ for all stellar masses. As in the models, however, there is a high correlation $r(k)$ between matter and galaxy density on all scales, and galaxy bias is typically consistent with a deterministic bias on linear scales. Only our blue and low-mass galaxies of about $7\times10^9\,h_{\rm 70}^{-2}\,{\rm M}_\odot$ at $\bar{z}=0.51$ show, contrary to the models, a weak tendency towards a stochastic bias on linear scales where $r_{\rm ls}=0.75\pm0.14\,{\rm(stat.)}\pm0.06\,{\rm(sys.)}$. This result is of interest for cosmological probes, such as $E_{\rm G}$, that rely on a deterministic galaxy bias.

preprint2020arXiv

Cosmic dissonance: new physics or systematics behind a short sound horizon?

Persistent tension between low-redshift observations and the Cosmic Microwave Background radiation (CMB), in terms of two fundamental distance scales set by the sound horizon $r_d$ and the Hubble constant $H_0$, suggests new physics beyond the Standard Model or residual systematics. We examine recently updated distance calibrations from Cepheids, gravitational lensing time-delay observations, and the Tip of the Red Giant Branch. Calibrating the Baryon Acoustic Oscillations (BAO) and Type Ia supernovae with combinations of the distance indicators, we obtain a joint and self-consistent measurement of $H_0$ and $r_d$ at low redshift, independent of cosmological models and CMB inference. In an attempt to alleviate the tension between late-time and CMB-based measurements, we consider four extensions of the standard $Λ$CDM model. The sound horizon from our different measurements is $r_d=(137\pm3^{stat.}\pm2^{syst.})$~Mpc. Depending on the adopted distance indicators, the $combined$ tension in $H_0$ and $r_d$ ranges between 2.3 and 5.1 $σ$. We find that modifications of $Λ$CDM that change the physics after recombination fail to solve the problem, for the reason that they only resolve the tension in $H_0$, while the tension in $r_d$ remains unchanged. Pre-recombination extensions (with early dark energy or the effective number of neutrinos $\rm{N}_{\rm{eff}}=3.24 \pm 0.16$) are allowed by the data, unless the calibration from Cepheids is included. Results from time-delay lenses are consistent with those from distance-ladder calibrations and point to a discrepancy between absolute distance scales measured from the CMB (assuming the standard cosmological model) and late-time observations. New proposals to resolve this tension should be examined with respect to reconciling not only the Hubble constant but also the sound horizon derived from the CMB and other cosmological probes.

preprint2020arXiv

The Accuracy of Weak Lensing Simulations

We investigate the accuracy of weak lensing simulations by comparing the results of five independently developed lensing simulation codes run on the same input $N$-body simulation. Our comparison focuses on the lensing convergence maps produced by the codes, and in particular on the corresponding PDFs, power spectra and peak counts. We find that the convergence power spectra of the lensing codes agree to $\lesssim 2\%$ out to scales $\ell \approx 4000$. For lensing peak counts, the agreement is better than $5\%$ for peaks with signal-to-noise $\lesssim 6$. We also discuss the systematic errors due to the Born approximation, line-of-sight discretization, particle noise and smoothing. The lensing codes tested deal in markedly different ways with these effects, but they nonetheless display a satisfactory level of agreement. Our results thus suggest that systematic errors due to the operation of existing lensing codes should be small. Moreover their impact on the convergence power spectra for a lensing simulation can be predicted given its numerical details, which may then serve as a validation test.

preprint2020arXiv

Time Delay Lens Modeling Challenge: I. Experimental Design

Strong gravitational lenses with measured time delay are a powerful tool to measure cosmological parameters, especially the Hubble constant ($H_0$). Recent studies show that by combining just three multiply-imaged AGN systems, one can determine $H_0$ to 2.4% precision. Furthermore, the number of time-delay lens systems is growing rapidly, enabling the determination of $H_0$ to 1% precision in the near future. However, as the precision increases it is important to ensure that systematic errors and biases remain subdominant. For this purpose, challenges with simulated datasets are a key component in this process. Following the experience of the past challenge on time delay, where it was shown that time delays can indeed be measured precisely and accurately at the sub-percent level, we now present the "Time Delay Lens Modeling Challenge" (TDLMC). The goal of this challenge is to assess the present capabilities of lens modeling codes and assumptions and test the level of accuracy of inferred cosmological parameters given realistic mock datasets. We invite scientists to model a set of simulated HST observations of 50 mock lens systems. The systems are organized in rungs, with the complexity and realism increasing going up the ladder. The goal of the challenge is to infer $H_0$ for each rung, given the HST images, the time delay, and stellar velocity dispersion of the deflector for a fixed background cosmology. The TDLMC challenge starts with the mock data release on 2018 January 8th. The deadline for blind submission is different for each rung. The deadline for Rung 0-1 is 2018 September 8; the deadline for Rung 2 is 2019 April 8 and the one for Rung 3 is 2019 September 8. This first paper gives an overview of the challenge including the data design, and a set of metrics to quantify the modeling performance and challenge details.

preprint2019arXiv

H0LiCOW XII. Lens mass model of WFI2033-4723 and blind measurement of its time-delay distance and $H_0$

We present the lens mass model of the quadruply-imaged gravitationally lensed quasar WFI2033-4723, and perform a blind cosmographical analysis based on this system. Our analysis combines (1) time-delay measurements from 14 years of data obtained by the COSmological MOnitoring of GRAvItational Lenses (COSMOGRAIL) collaboration, (2) high-resolution $\textit{Hubble Space Telescope}$ imaging, (3) a measurement of the velocity dispersion of the lens galaxy based on ESO-MUSE data, and (4) multi-band, wide-field imaging and spectroscopy characterizing the lens environment. We account for all known sources of systematics, including the influence of nearby perturbers and complex line-of-sight structure, as well as the parametrization of the light and mass profiles of the lensing galaxy. After unblinding, we determine the effective time-delay distance to be $4784_{-248}^{+399}~\mathrm{Mpc}$, an average precision of $6.6\%$. This translates to a Hubble constant $H_{0} = 71.6_{-4.9}^{+3.8}~\mathrm{km~s^{-1}~Mpc^{-1}}$, assuming a flat $Λ$CDM cosmology with a uniform prior on $Ω_\mathrm{m}$ in the range [0.05, 0.5]. This work is part of the $H_0$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration, and the full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper (H0LiCOW XIII).

preprint2019arXiv

H0LiCOW XIII. A 2.4% measurement of $H_{0}$ from lensed quasars: $5.3σ$ tension between early and late-Universe probes

We present a measurement of the Hubble constant ($H_{0}$) and other cosmological parameters from a joint analysis of six gravitationally lensed quasars with measured time delays. All lenses except the first are analyzed blindly with respect to the cosmological parameters. In a flat $Λ$CDM cosmology, we find $H_{0} = 73.3_{-1.8}^{+1.7}$, a 2.4% precision measurement, in agreement with local measurements of $H_{0}$ from type Ia supernovae calibrated by the distance ladder, but in $3.1σ$ tension with $Planck$ observations of the cosmic microwave background (CMB). This method is completely independent of both the supernovae and CMB analyses. A combination of time-delay cosmography and the distance ladder results is in $5.3σ$ tension with $Planck$ CMB determinations of $H_{0}$ in flat $Λ$CDM. We compute Bayes factors to verify that all lenses give statistically consistent results, showing that we are not underestimating our uncertainties and are able to control our systematics. We explore extensions to flat $Λ$CDM using constraints from time-delay cosmography alone, as well as combinations with other cosmological probes, including CMB observations from $Planck$, baryon acoustic oscillations, and type Ia supernovae. Time-delay cosmography improves the precision of the other probes, demonstrating the strong complementarity. Allowing for spatial curvature does not resolve the tension with $Planck$. Using the distance constraints from time-delay cosmography to anchor the type Ia supernova distance scale, we reduce the sensitivity of our $H_0$ inference to cosmological model assumptions. For six different cosmological models, our combined inference on $H_{0}$ ranges from $\sim73$-$78~\mathrm{km~s^{-1}~Mpc^{-1}}$, which is consistent with the local distance ladder constraints.

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.

preprint2016arXiv

Cosmic Shear Results from the Deep Lens Survey - II: Full Cosmological Parameter Constraints from Tomography

We present a tomographic cosmic shear study from the Deep Lens Survey (DLS), which, providing a limiting magnitude r_{lim}~27 (5 sigma), is designed as a pre-cursor Large Synoptic Survey Telescope (LSST) survey with an emphasis on depth. Using five tomographic redshift bins, we study their auto- and cross-correlations to constrain cosmological parameters. We use a luminosity-dependent nonlinear model to account for the astrophysical systematics originating from intrinsic alignments of galaxy shapes. We find that the cosmological leverage of the DLS is among the highest among existing >10 sq. deg cosmic shear surveys. Combining the DLS tomography with the 9-year results of the Wilkinson Microwave Anisotropy Probe (WMAP9) gives Omega_m=0.293_{-0.014}^{+0.012}, sigma_8=0.833_{-0.018}^{+0.011}, H_0=68.6_{-1.2}^{+1.4} km/s/Mpc, and Omega_b=0.0475+-0.0012 for LCDM, reducing the uncertainties of the WMAP9-only constraints by ~50%. When we do not assume flatness for LCDM, we obtain the curvature constraint Omega_k=-0.010_{-0.015}^{+0.013} from the DLS+WMAP9 combination, which however is not well constrained when WMAP9 is used alone. The dark energy equation of state parameter w is tightly constrained when Baryonic Acoustic Oscillation (BAO) data are added, yielding w=-1.02_{-0.09}^{+0.10} with the DLS+WMAP9+BAO joint probe. The addition of supernova constraints further tightens the parameter to w=-1.03+-0.03. Our joint constraints are fully consistent with the final Planck results and also the predictions of a LCDM universe.

preprint2015arXiv

CFHTLenS: Weak lensing constraints on the ellipticity of galaxy-scale matter haloes and the galaxy-halo misalignment

We present weak lensing constraints on the ellipticity of galaxy-scale matter haloes and the galaxy-halo misalignment. Using data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), we measure the weighted-average ratio of the aligned projected ellipticity components of galaxy matter haloes and their embedded galaxies, $f_\mathrm{h}$, split by galaxy type. We then compare our observations to measurements taken from the Millennium Simulation, assuming different models of galaxy-halo misalignment. Using the Millennium Simulation we verify that the statistical estimator used removes contamination from cosmic shear. We also detect an additional signal in the simulation, which we interpret as the impact of intrinsic shape-shear alignments between the lenses and their large-scale structure environment. These alignments are likely to have caused some of the previous observational constraints on $f_\mathrm{h}$ to be biased high. From CFHTLenS we find $f_\mathrm{h}=-0.04 \pm 0.25$ for early-type galaxies, which is consistent with current models for the galaxy-halo misalignment predicting $f_\mathrm{h}\simeq 0.20$. For late-type galaxies we measure $f_\mathrm{h}=0.69_{-0.36}^{+0.37}$ from CFHTLenS. This can be compared to the simulated results which yield $f_\mathrm{h}\simeq 0.02$ for misaligned late-type models.

preprint2015arXiv

Meaning of temperature in different thermostatistical ensembles

Depending on the exact experimental conditions, the thermodynamic properties of physical systems can be related to one or more thermostatistical ensembles. Here, we survey the notion of thermodynamic temperature in different statistical ensembles, focusing in particular on subtleties that arise when ensembles become non-equivalent. The 'mother' of all ensembles, the microcanonical ensemble, uses entropy and internal energy (the most fundamental, dynamically conserved quantity) to derive temperature as a secondary thermodynamic variable. Over the past century, some confusion has been caused by the fact that several competing microcanonical entropy definitions are used in the literature, most commonly the volume and surface entropies introduced by Gibbs. It can be proved, however, that only the volume entropy satisfies exactly the traditional form of the laws of thermodynamics for a broad class of physical systems, including all standard classical Hamiltonian systems, regardless of their size. This mathematically rigorous fact implies that negative 'absolute' temperatures and Carnot efficiencies $>1$ are not achievable within a standard thermodynamical framework. As an important offspring of microcanonical thermostatistics, we shall briefly consider the canonical ensemble and comment on the validity of the Boltzmann weight factor. We conclude by addressing open mathematical problems that arise for systems with discrete energy spectrum.

preprint2014arXiv

Cosmological constraints from the CFHTLenS shear measurements using a new, accurate and flexible way of predicting nonlinear mass clustering

We explore the cosmological constraints from cosmic shear using a new way of modelling the non-linear matter correlation functions. The new formalism extends the method of Angulo & White (2010), which manipulates outputs of $N$-body simulations to represent the three-dimensional nonlinear mass distribution in different cosmological scenarios. We show that predictions from our approach for shear two-point correlations at $1$ to $300$ arcmin separations are accurate at the $\sim10$\% level, even for extreme changes in cosmology. For moderate changes, with target cosmologies similar to that preferred by analyses of recent Planck data, the accuracy is close to $\sim5$\%. We combine this approach with a MonteCarlo Markov Chain sampler to explore constraints on a $Λ$CDM model from the shear correlation functions measured in the Canada-France Hawaii Telescope Lensing Survey (CFHTLenS). We obtain constraints on the parameter combination $σ_8 (Ω_m/0.27)^{0.6} = 0.801 \pm 0.028$. Combined with results from CMB data, we obtain marginalised constraints on $σ_8 = 0.81 \pm 0.01$ and $Ω_m = 0.29 \pm 0.01$. These results are fully compatible with previous analyses, which supports the validity of our approach. We discuss the advantages of our method and the potential it offers, including a path to incorporate in detail the effects of baryons, among others effects, in future high-precision cosmological analyses.

preprint2014arXiv

Thermodynamic laws in isolated systems

The recent experimental realization of exotic matter states in isolated quantum systems and the ensuing controversy about the existence of negative absolute temperatures demand a careful analysis of the conceptual foundations underlying microcanonical thermostatistics. Here, we provide a detailed comparison of the most commonly considered microcanonical entropy definitions, focussing specifically on whether they satisfy or violate the zeroth, first and second law of thermodynamics. Our analysis shows that, for a broad class of systems that includes all standard classical Hamiltonian systems, only the Gibbs volume entropy fulfills all three laws simultaneously. To avoid ambiguities, the discussion is restricted to exact results and analytically tractable examples.

preprint2013arXiv

CFHTLenS: The Environmental Dependence of Galaxy Halo Masses from Weak Lensing

We use weak gravitational lensing to analyse the dark matter halos around satellite galaxies in galaxy groups in the CFHTLenS dataset. This dataset is derived from the CFHTLS-Wide survey, and encompasses 154 sq. deg of high-quality shape data. Using the photometric redshifts, we divide the sample of lens galaxies with stellar masses in the range 10^9 Msun to 10^10.5 Msun into those likely to lie in high-density environments (HDE) and those likely to lie in low-density environments (LDE). Through comparison with galaxy catalogues extracted from the Millennium Simulation, we show that the sample of HDE galaxies should primarily (~61%) consist of satellite galaxies in groups, while the sample of LDE galaxies should consist of mostly (~87%) non-satellite (field and central) galaxies. Comparing the lensing signals around samples of HDE and LDE galaxies matched in stellar mass, the lensing signal around HDE galaxies clearly shows a positive contribution from their host groups on their lensing signals at radii of ~500--1000 kpc, the typical separation between satellites and group centres. More importantly, the subhalos of HDE galaxies are less massive than those around LDE galaxies by a factor 0.65 +/- 0.12, significant at the 2.9 sigma level. A natural explanation is that the halos of satellite galaxies are stripped through tidal effects in the group environment. Our results are consistent with a typical tidal truncation radius of ~40 kpc.

preprint2013arXiv

Cosmic shear results from the deep lens survey - I: Joint constraints on omega_m and sigma_8 with a two-dimensional analysis

We present a cosmic shear study from the Deep Lens Survey (DLS), a deep BVRz multi-band imaging survey of five 4 sq. degree fields with two National Optical Astronomy Observatory (NOAO) 4-meter telescopes at Kitt Peak and Cerro Tololo. For both telescopes, the change of the point-spread-function (PSF) shape across the focal plane is complicated, and the exposure-to-exposure variation of this position-dependent PSF change is significant. We overcome this challenge by modeling the PSF separately for individual exposures and CCDs with principal component analysis (PCA). We find that stacking these PSFs reproduces the final PSF pattern on the mosaic image with high fidelity, and the method successfully separates PSF-induced systematics from gravitational lensing effects. We calibrate our shears and estimate the errors, utilizing an image simulator, which generates sheared ground-based galaxy images from deep Hubble Space Telescope archival data with a realistic atmospheric turbulence model. For cosmological parameter constraints, we marginalize over shear calibration error, photometric redshift uncertainty, and the Hubble constant. We use cosmology-dependent covariances for the Markov Chain Monte Carlo analysis and find that the role of this varying covariance is critical in our parameter estimation. Our current non-tomographic analysis alone constrains the Omega_M-sigma_8 likelihood contour tightly, providing a joint constraint of Omega_M=0.262+-0.051 and sigma_8=0.868+-0.071. We expect that a future DLS weak-lensing tomographic study will further tighten these constraints because explicit treatment of the redshift dependence of cosmic shear more efficiently breaks the Omega_M-sigma_8 degeneracy. Combining the current results with the Wilkinson Microwave Anisotropy Probe 7-year (WMAP7) likelihood data, we obtain Omega_M=0.278+-0.018 and sigma_8=0.815+-0.020.

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

Improving the precision of time-delay cosmography with observations of galaxies along the line of sight

In order to use strong gravitational lens time delays to measure precise and accurate cosmological parameters the effects of mass along the line of sight must be taken into account. We present a method to achieve this by constraining the probability distribution function of the effective line of sight convergence k_ext. The method is based on matching the observed overdensity in the weighted number of galaxies to that found in mock catalogs with k_ext obtained by ray-tracing through structure formation simulations. We explore weighting schemes based on projected distance, mass, luminosity, and redshift. This additional information reduces the uncertainty of k_ext from sigma_k $0.06 to ~0.04 for very overdense lines of sight like that of the system B1608+656. For more common lines of sight, sigma_k is reduced to ~<0.03, corresponding to an uncertainty of ~<3% on distance. This uncertainty has comparable effects on cosmological parameters to that arising from the mass model of the deflector and its immediate environment. Photometric redshifts based on g, r, i and K photometries are sufficient to constrain k_ext almost as well as with spectroscopic redshifts. As an illustration, we apply our method to the system B1608+656. Our most reliable k_ext estimator gives sigma_k=0.047 down from 0.065 using only galaxy counts. Although deeper multi-band observations of the field of B1608+656 are necessary to obtain a more precise estimate, we conclude that griK photometry, in addition to spectroscopy to characterize the immediate environment, is an effective way to increase the precision of time-delay cosmography.

preprint2013arXiv

Inconsistent thermostatistics and negative absolute temperatures

A considerable body of experimental and theoretical work claims the existence of negative absolute temperatures in spin systems and ultra-cold quantum gases. Here, we clarify that such findings can be attributed to the use of a popular yet inconsistent entropy definition, which violates fundamental thermodynamic relations and fails to produce sensible results for simple analytically tractable classical and quantum systems. Within a mathematically consistent thermodynamic formalism, based on an entropy concept originally derived by Gibbs, absolute temperature remains positive even for systems with bounded spectrum. We address spurious arguments against the Gibbs formalism and comment briefly on heat engines with efficiencies greater than one.

preprint2013arXiv

Noiseless Gravitational Lensing Simulations

The microphysical properties of the DM particle can, in principle, be constrained by the properties and abundance of substructures in DM halos, as measured through strong gravitational lensing. Unfortunately, there is a lack of accurate theoretical predictions for the lensing signal of substructures, mainly because of the discreteness noise inherent to N-body simulations. Here we present Recursive-TCM, a method that is able to provide lensing predictions with an arbitrarily low discreteness noise, without any free parameters or smoothing scale. This solution is based on a novel way of interpreting the results of N-body simulations, where particles simply trace the evolution and distortion of Lagrangian phase-space volume elements. We discuss the advantages of this method over the widely used cloud-in-cells and adaptive-kernel smoothing density estimators. Applying the new method to a cluster-sized DM halo simulated in warm and cold DM scenarios, we show how the expected differences in their substructure population translate into differences in the convergence and magnification maps. We anticipate that our method will provide the high-precision theoretical predictions required to interpret and fully exploit strong gravitational lensing observations.

preprint2013arXiv

Reconstructing the Lensing Mass in the Universe from Photometric Catalogue Data

High precision cosmological distance measurements towards individual objects such as time delay gravitational lenses or type Ia supernovae are affected by weak lensing perturbations by galaxies and groups along the line of sight. In time delay gravitational lenses, "external convergence," kappa, can dominate the uncertainty in the inferred distances and hence cosmological parameters. In this paper we attempt to reconstruct kappa, due to line of sight structure, using a simple halo model. We use mock catalogues from the Millennium Simulation, and calibrate and compare our reconstructed P(kappa) to ray-traced kappa "truth" values; taking into account realistic observational uncertainties. We find that the reconstruction of kappa provides an improvement in precision of ~50% over galaxy number counts. We find that the lowest-kappa lines of sight have the best constrained P(kappa). In anticipation of large future samples of lenses, we find that selecting the third of the systems with the highest precision kappa estimates gives a sample of unbiased time delay distance measurements with just ~1% uncertainty due to line of sight external convergence effects. Photometric data are sufficient to pre-select the best-constrained lines of sight, and can be done before investment in light-curve monitoring. Conversely, we show that selecting lines of sight with high external shear could, with the reconstruction model presented, induce biases of up to 1% in time delay distance. We find that a major potential source of systematic error is uncertainty in the high mass end of the stellar mass-halo mass relation; this could introduce ~2% biases on the time-delay distance if completely ignored. We suggest areas for the improvement of this general analysis framework (including more sophisticated treatment of high mass structures) that should allow yet more accurate cosmological inferences to be made.

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.

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

Measuring primordial non-Gaussianity with weak-lensing surveys

We study the ability of future weak lensing (WL) surveys to constrain primordial non-Gaussianity of the local type. We use a large ensemble of simulated WL maps with survey specifications relevant to Euclid and LSST. The simulations assume Cold Dark Matter cosmologies that vary certain parameters around fiducial values: the non-Gaussianity parameter f_NL, the matter density parameter Omega_m, the amplitude of the matter power spectrum sigma_8, the spectral index of the primordial power spectrum n_s, and the dark-energy equation-of-state parameter w_0. We assess the sensitivity of the cosmic shear correlation functions, the third-order aperture mass statistics, and the abundance of shear peaks to these parameters. We find that each of the considered probes provides unmarginalized constraints of Delta f_NL ~ 20 on f_NL. Marginalized constraints from any individual WL probe are much weaker due to strong correlations between parameters. However, the parameter errors can be substantially reduced by combining information from different WL probes. Combining all WL probes yields the following marginal (68% confidence level) uncertainties: Delta f_NL ~ 50, Delta Omega_m ~ 0.002, Delta sigma_8 ~ 0.004, Delta n_s ~ 0.007, and Delta w_0 ~ 0.03. We examine the bias induced by neglecting f_NL on the constraints on the other parameters. We find sigma_8 and w_0 to be the most affected. Moreover, neglecting non-Gaussianity leads to a severe underestimation of the uncertainties in the other cosmological parameters.

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

Probing Satellite Halos with Weak Gravitational Lensing

We demonstrate the possibility of detecting tidal stripping of dark matter subhalos within galaxy groups using weak gravitational lensing. We have run ray-tracing simulations on galaxy catalogues from the Millennium Simulation to generate mock shape catalogues. The ray-tracing catalogues assume a halo model for galaxies and groups, using various models with different distributions of mass between galaxy and group halos to simulate different stages of group evolution. Using these mock catalogues, we forecast the lensing signals that will be detected around galaxy groups and satellite galaxies, as well as test two different methods for isolating the satellites' lensing signals. A key challenge is to determine the accuracy to which group centres can be identified. We show that with current and ongoing surveys, it will possible to detect stripping in groups of mass 10^12--10^15 Msun.

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

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

Lensed arc statistics: comparison of Millennium-simulation galaxy clusters to Hubble Space Telescope observations of an X-ray selected sample

It has been debated for a decade whether there is a large overabundance of strongly lensed arcs in galaxy clusters, compared to expectations from LambdaCDM cosmology. We perform ray tracing through the most massive halos of the Millennium simulation at several redshifts in their evolution, using the Hubble Ultra Deep Field as a source image, to produce realistic simulated lensed images. We compare the lensed arc statistics measured from the simulations to those of a sample of 45 X-ray selected clusters, observed with the Hubble Space Telescope, that we have analysed in Horesh et al. (2010). The observations and the simulations are matched in cluster masses, redshifts, observational effects, and the algorithmic arc detection and selection. At z=0.6 there are too few massive-enough clusters in the Millennium volume for a proper statistical comparison with the observations. At redshifts 0.3<z<0.5, however, we have large numbers of simulated and observed clusters, and the latter are an unbiased selection from a complete sample. For these redshifts, we find excellent agreement between the observed and simulated arc statistics, in terms of the mean number of arcs per cluster, the distribution of number of arcs per cluster, and the angular separation distribution. At z ~ 0.2 some conflict remains, with real clusters being ~3 times more efficient arc producers than their simulated counterparts. This may arise due to selection biases in the observed subsample at this redshift, to some mismatch in masses between the observed and simulated clusters, or to physical effects that arise at low redshift and enhance the lensing efficiency, but which are not represented by the simulations.

preprint2010arXiv

Effects of large-scale structure on the accuracy of weak lensing mass measurements

Weak gravitational lensing has become an important method to determine the masses of galaxy clusters. The intrinsic shapes of the galaxies are a dominant source of uncertainty, but there are other limitations to the precision that can be achieved. In this paper we revisit a typically ignored source of uncertainty: structure along the line-of sight. Using results from the Millennium Simulation we confirm the validity of analytical calculations that have shown that such random projections are particularly important for studies of the cluster density profile. In general the contribution of large-scale structure to the total error budget is comparable to the statistical errors. We find that the precision of the mass measurement can be improved only slightly by modelling the large-scale structure using readily available data.

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

Weak lensing from space: first cosmological constraints from three-point shear statistics

We use weak lensing data from the Hubble Space Telescope COSMOS survey to measure the second- and third-moments of the cosmic shear field, estimated from about 450,000 galaxies with average redshift <z> ~ 1.3. We measure two- and three-point shear statistics using a tree-code, dividing the signal in E, B and mixed components. We present a detection of the third-order moment of the aperture mass statistic and verify that the measurement is robust against systematic errors caused by point spread function (PSF) residuals and by the intrinsic alignments between galaxies. The amplitude of the measured three-point cosmic shear signal is in very good agreement with the predictions for a WMAP7 best-fit model, whereas the amplitudes of potential systematics are consistent with zero. We make use of three sets of large Lambda CDM simulations to test the accuracy of the cosmological predictions and to estimate the influence of the cosmology-dependent covariance. We perform a likelihood analysis using the measurement and find that the Omega_m-sigma_8 degeneracy direction is well fitted by the relation: sigma_8 (Omega_m/0.30)^(0.49)=0.78+0.11/-0.26. We present the first measurement of a more generalised three-point shear statistic and find a very good agreement with the WMAP7 best-fit cosmology. The cosmological interpretation of this measurement gives sigma_8 (Omega_m/0.30)^(0.46)=0.69 +0.08/-0.14. Furthermore, the combined likelihood analysis of this measurement with the measurement of the second order moment of the aperture mass improves the accuracy of the cosmological constraints, showing the high potential of this combination of measurements to infer cosmological constraints.

preprint2009arXiv

Abundances, masses, and weak-lensing mass profiles of galaxy clusters as a function of richness and luminosity in LambdaCDM cosmologies

We test the concordance LCDM cosmology by comparing predictions for the mean properties of galaxy clusters to observations. We use high-resolution N-body simulations of cosmic structure formation and semi-analytic models of galaxy formation to compute the abundance, mean density profile, and mass of galaxy clusters as a function of richness and luminosity, and we compare these predictions to observations of clusters in the Sloan Digital Sky Survey (SDSS) maxBCG catalogue. We discuss the scatter in the mass-richness relation, the reconstruction of the cluster mass function from the mass-richness relation, and fits to the weak-lensing cluster mass profiles. The impact of cosmological parameters on the predictions is investigated by comparing results from galaxy models based on the Millennium Simulation (MS) and another WMAP1 simulation to those from a WMAP3 simulation. We find that the simulated weak-lensing mass profiles and the observed profiles of the SDSS maxBCG clusters agree well in shape and amplitude. The mass-richness relations in the simulations are close to the observed relation, with differences lesssim 30%. The MS and WMAP1 simulations yield cluster abundances similar to those observed, whereas abundances in the WMAP3 simulation are 2-3 times lower. The differences in cluster abundance, mass, and density amplitude between the simulations and the observations can be attributed to differences in the underlying cosmological parameters, in particular the power spectrum normalisation sigma_8. Better agreement between predictions and observations should be reached with a normalisation $0.722<sigma_8<0.9$ (probably closer to the upper value), i.e. between the values underlying the two simulation sets.

preprint2009arXiv

Cluster strong lensing in the Millennium Simulation: The effect of galaxies and structures along the line-of-sight

We use ray-tracing through the Millennium simulation to study how secondary matter structures along the line-of-sight and the stellar mass in galaxies affect strong cluster lensing, in particular the cross-section for giant arcs. Furthermore, we investigate the distribution of the cluster Einstein radii and the radial distribution of giant arcs. We find that additional structures along the line-of-sight increase the strong-lensing optical depth by ~10-25%, while strong-lensing cross-sections of individual clusters are frequently boosted by as much as ~50%. The enhancement is mainly due to structures that are not correlated with the lens. Cluster galaxies increase the strong-lensing optical depth by up to a factor of 2, while interloping galaxies are not significant. We conclude that these effects need to be taken into account for predictions of the giant arc abundance, but they are not large enough to fully account for the reported discrepancy between predicted and observed abundances. Furthermore, we find that Einstein radii defined via the area enclosed by the critical curve are 10-30% larger than those defined via radial surface mass density profiles. The contributions of radial and tangential arcs to the radial distribution of arcs can be clearly distinguished. The radial distribution of tangential arcs is very broad and extends out to several Einstein radii. Thus, individual arcs are not well suited for constraining Einstein radii.