BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Benjamin Joachimi

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

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-ph.IMgr-qchep-ph

Catalog footprint

What is connected

47works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

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

0 citations0 reviewsNo code linkedOpen access
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}$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Explaining deep learning of galaxy morphology with saliency mapping

We successfully demonstrate the use of explainable artificial intelligence (XAI) techniques on astronomical datasets in the context of measuring galactic bar lengths. The method consists of training convolutional neural networks on human classified data from Galaxy Zoo in order to predict general galaxy morphologies, and then using SmoothGrad (a saliency mapping technique) to extract the bar for measurement by a bespoke algorithm. We contrast this to another method of using a convolutional neural network to directly predict galaxy bar lengths. These methods achieved correlation coefficients of 0.76 and 0.59, and root mean squared errors of 1.69 and 2.10 respective to human measurements. We conclude that XAI methods outperform conventional deep learning in this case, which could be reasonably explained by the larger datasets available when training the models. We suggest that our XAI method can be used to extract other galactic features (such as the bulge-to-disk ratio) without needing to collect new datasets or train new models. We also suggest that these techniques can be used to refine deep learning models as well as identify and eliminate bias within training datasets.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Forecasting the potential of weak lensing magnification to enhance LSST large-scale structure analyses

Recent works have shown that weak lensing magnification must be included in upcoming large-scale structure analyses, such as for the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), to avoid biasing the cosmological results. In this work we investigate whether including magnification has a positive impact on the precision of the cosmological constraints, as well as being necessary to avoid bias. We forecast this using an LSST mock catalog and a halo model to calculate the galaxy power spectra. We find that including magnification has little effect on the precision of the cosmological parameter constraints for an LSST galaxy clustering analysis, where the halo model parameters are additionally constrained by the galaxy luminosity function. In particular, we find that for the LSST gold sample ($i < 25.3$) including weak lensing magnification only improves the galaxy clustering constraint on $Ω_{\rm{m}}$ by a factor of 1.03, and when using a very deep LSST mock sample ($i<26.5$) by a factor of 1.3. Since magnification predominantly contributes to the clustering measurement and provides similar information to that of cosmic shear, this improvement would be reduced for a combined galaxy clustering and shear analysis. We also confirm that not modelling weak lensing magnification will catastrophically bias the cosmological results from LSST. Magnification must therefore be included in LSST large-scale structure analyses even though it does not significantly enhance the precision of the cosmological constraints.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Three-point intrinsic alignments of dark matter halos in the IllustrisTNG simulation

We use the IllustrisTNG suite of cosmological simulations to measure intrinsic alignment (IA) bispectra of dark matter subhalos between redshifts 0 and 1. We decompose the intrinsic shear field into E- and B-modes and find that the bispectra $B_{δδ\mathrm{E}}$ and $B_{δ\mathrm{EE}}$, between the matter overdensity field, $δ$, and the E-mode field, are detected with high significance. We also model the IA bispectra analytically using a method consistent with the two-point non-linear alignment model. We use this model and the simulation measurements to infer the intrinsic alignment amplitude $A_\mathrm{IA}$ and find that values of $A_\mathrm{IA}$ obtained from IA power spectra and bispectra agree well at scales up to $k_\mathrm{max}= 2 \, h \mathrm{Mpc}^{-1}$, for example at $z=1$ $A_\mathrm{IA} = 2.13 \pm$ 0.02 from the cross power spectrum between the matter overdensity and E-mode fields and $A_\mathrm{IA} =2.11 \pm$ 0.03 from $B_{δδ\mathrm{E}}$. This demonstrates that a single physically motivated model can jointly model two-point and three-point statistics of intrinsic alignments, thus enabling a cleaner separation between intrinsic alignments and cosmological weak lensing signals.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Group-scale intrinsic galaxy alignments in the Illustris-TNG and MassiveBlack-II simulations

We study the alignments of satellite galaxies, and their anisotropic distribution, with respect to location and orientation of their host central galaxy in MassiveBlack-II and IllustrisTNG simulations. We find that: the shape of the satellite system in halos of mass ($> 10^{13}h^{-1}M_{\odot}$) is well aligned with the shape of the central galaxy at $z=0.06$ with the mean alignment between the major axes being $\sim Δθ= 12^{\circ}$ when compared to a uniform random distribution; that satellite galaxies tend to be anisotropically distributed along the major axis of the central galaxy with a stronger alignment in halos of higher mass or luminosity; and that the satellite distribution is more anisotropic for central galaxies with lower star formation rate, which are spheroidal, and for red central galaxies.Radially we find that satellites tend to be distributed along the major axis of the shape of the stellar component of central galaxies at smaller scales and the dark matter component on larger scales. We find that the dependence of satellite anisotropy on central galaxy properties and the radial distance is similar in both the simulations with a larger amplitude in MassiveBlack-II. The orientation of satellite galaxies tends to point toward the location of the central galaxy at small scales and this correlation decreases with increasing distance, and the amplitude of satellite alignment is higher in high mass halos. However, the projected ellipticities do not exhibit a scale-dependent radial alignment, as has been seen in some observational measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

KiDS+VIKING-450: An internal-consistency test for cosmic shear tomography with a colour-based split of source galaxies

We performed an internal-consistency test of the KiDS+VIKING-450 (KV450) cosmic shear analysis with a colour-based split of source galaxies. Utilising the same measurements and calibrations for both sub-samples, we inspected the characteristics of the shear measurements and the performance of the calibration pipelines. On the modelling side, we examined the observational nuisance parameters, specifically those for the redshift calibration and intrinsic alignments, using a Bayesian analysis with dedicated test parameters. We verified that the current nuisance parameters are sufficient for the KV450 data to capture residual systematics, with slight deviations seen in the second and the third redshift tomographic bins. Our test also showcases the degeneracy between the inferred amplitude of intrinsic alignments and the redshift uncertainties in low redshift tomographic bins. The test is rather insensitive to the background cosmology and, therefore, can be implemented before any cosmological inference is made.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Self-calibration of weak lensing systematic effects using combined two- and three-point statistics

We investigate the prospects for using the weak lensing bispectrum alongside the power spectrum to control systematic uncertainties in a Euclid-like survey. Three systematic effects are considered: the intrinsic alignment of galaxies, uncertainties in the means of tomographic redshift distributions, and multiplicative bias in the measurement of the shear signal. We find that the bispectrum is very effective in mitigating these systematic errors. Varying all three systematics simultaneously, a joint power spectrum and bispectrum analysis reduces the area of credible regions for the cosmological parameters $Ω_\mathrm{m}$ and $σ_8$ by a factor of 90 and for the two parameters of a time-varying dark energy equation of state by a factor of almost 20, compared with the baseline approach of using the power spectrum alone and of imposing priors consistent with the accuracy requirements specified for Euclid. We also demonstrate that including the bispectrum self-calibrates all three systematic effects to the stringent levels required by the forthcoming generation of weak lensing surveys, thereby reducing the need for external calibration data.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The PAU Survey: Intrinsic alignments and clustering of narrow-band photometric galaxies

We present the first measurements of the projected clustering and intrinsic alignments (IA) of galaxies observed by the Physics of the Accelerating Universe Survey (PAUS). With photometry in 40 narrow optical passbands ($450\rm{nm}-850\rm{nm}$), the quality of photometric redshift estimation is $σ_{z} \sim 0.01(1 + z)$ for galaxies in the $19\,\rm{deg}^{2}$ Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) W3 field, allowing us to measure the projected 3D clustering and IA for flux-limited, faint galaxies ($i < 22.5$) out to $z\sim0.8$. To measure two-point statistics, we developed, and tested with mock photometric redshift samples, `cloned' random galaxy catalogues which can reproduce data selection functions in 3D and account for photometric redshift errors. In our fiducial colour-split analysis, we made robust null detections of IA for blue galaxies and tentative detections of radial alignments for red galaxies ($\sim1-3σ$), over scales of $0.1-18\,h^{-1}\rm{Mpc}$. The galaxy clustering correlation functions in the PAUS samples are comparable to their counterparts in a spectroscopic population from the Galaxy and Mass Assembly survey, modulo the impact of photometric redshift uncertainty which tends to flatten the blue galaxy correlation function, whilst steepening that of red galaxies. We investigate the sensitivity of our correlation function measurements to choices in the random catalogue creation and the galaxy pair-binning along the line of sight, in preparation for an optimised analysis over the full PAUS area.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bayesian methods for fitting Baryon Acoustic Oscillations in the Lyman-$α$ forest

We study and compare fitting methods for the Lyman-$α$ (Ly$α$) forest 3D correlation function. We use the nested sampler PolyChord and the community code picca to perform a Bayesian analysis which we compare with previous frequentist analyses. By studying synthetic correlation functions, we find that the frequentist profile likelihood produces results in good agreement with a full Bayesian analysis. On the other hand, Maximum Likelihood Estimation with the Gaussian approximation for the uncertainties is inadequate for current data sets. We compute for the first time the full posterior distribution from the Ly$α$ forest correlation functions measured by the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We highlight the benefits of sampling the full posterior distribution by expanding the baseline analysis to better understand the contamination by Damped Ly$α$ systems (DLAs). We make our improvements and results publicly available as part of the picca package.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

GEOMAX: beyond linear compression for 3pt galaxy clustering statistics

We present the GEOMAX algorithm and its Python implementation for a two-step compression of bispectrum measurements. The first step groups bispectra by the geometric properties of their arguments; the second step then maximises the Fisher information with respect to a chosen set of model parameters in each group. The algorithm only requires the derivatives of the data vector with respect to the parameters and a small number of mock data, producing an effective, non-linear compression. By applying GEOMAX to bispectrum monopole measurements from BOSS DR12 CMASS redshift-space galaxy clustering data, we reduce the $68\%$ credible intervals for the inferred parameters $\left(b_1,b_2,f,σ_8\right)$ by $\left(50.4\%,56.1\%,33.2\%,38.3\%\right)$ with respect to standard MCMC on the full data vector. We run the analysis and comparison between compression methods over one hundred galaxy mocks to test the statistical significance of the improvements. On average GEOMAX performs $\sim15\%$ better than geometrical or maximal linear compression alone and is consistent with being lossless. Given its flexibility, the GEOMAX approach has the potential to optimally exploit three-point statistics of various cosmological probes like weak lensing or line-intensity maps from current and future cosmological data-sets such as DESI, Euclid, PFS and SKA.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Target Neutrino Mass Precision for Determining the Neutrino Hierarchy

Recent works combining neutrino oscillation and cosmological data to determine the neutrino hierarchy found a range of odds in favour of the normal hierarchy. These results arise from differing approaches to incorporating prior knowledge about neutrinos. We develop a hierarchy-agnostic prior and show that the hierarchy cannot be conclusively determined with current data. The determination of the hierarchy is limited by the neutrino mass scale $Σ_ν$ measurement. We obtain a target precision of $σ(Σ_ν) = 0.014$ eV, necessary for conclusively establishing the normal hierarchy with future data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The PAU Survey: An improved photo-$z$ sample in the COSMOS field

We present -- and make publicly available -- accurate and precise photometric redshifts in the ACS footprint from the COSMOS field for objects with $i_{\mathrm{AB}}\leq 23$. The redshifts are computed using a combination of narrow band photometry from PAUS, a survey with 40 narrow bands spaced at $100Å$ intervals covering the range from $4500Å$ to $8500Å$, and 26 broad, intermediate, and narrow bands covering the UV, visible and near infrared spectrum from the COSMOS2015 catalogue. We introduce a new method that models the spectral energy distributions (SEDs) as a linear combination of continuum and emission line templates and computes its Bayes evidence, integrating over the linear combinations. The correlation between the UV luminosity and the OII line is measured using the 66 available bands with the zCOSMOS spectroscopic sample, and used as a prior which constrains the relative flux between continuum and emission line templates. The flux ratios between the OII line and $\mathrm{H}_α$, $\mathrm{H}_β$ and $\mathrm{OIII}$ are similarly measured and used to generate the emission line templates. Comparing to public spectroscopic surveys via the quantity $Δ_z\equiv(z_{\mathrm{photo}}-z_{\mathrm{spec}})/(1+z_{\mathrm{spec}})$, we find the photometric redshifts to be more precise than previous estimates, with $σ_{68}(Δ_z) \approx (0.003, 0.009)$ for galaxies at magnitude $i_{\mathrm{AB}}\sim18$ and $i_{\mathrm{AB}}\sim23$, respectively, which is $3\times$ and $1.66\times$ tighter than COSMOS2015. Additionally, we find the redshifts to be very accurate on average, yielding a median of the $Δ_z$ distribution compatible with $|\mathrm{median}(Δ_z)|\leq0.001$ at all redshifts and magnitudes considered. Both the added PAUS data and new methodology contribute significantly to the improved results.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

CosmoDC2: A Synthetic Sky Catalog for Dark Energy Science with LSST

This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `Outer Rim' run. It covers 440 deg^2 of sky area to a redshift of z=3 and is complete to a magnitude depth of 28 in the r-band. Each galaxy is characterized by a multitude of properties including stellar mass, morphology, spectral energy distributions, broadband filter magnitudes, host halo information and weak lensing shear. The size and complexity of cosmoDC2 requires an efficient catalog generation methodology; our approach is based on a new hybrid technique that combines data-driven empirical approaches with semi-analytic galaxy modeling. A wide range of observation-based validation tests has been implemented to ensure that cosmoDC2 enables the science goals of the planned LSST DESC DC2 analyses. This paper also represents the official release of the cosmoDC2 data set, including an efficient reader that facilitates interaction with the data.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Cosmology from large-scale structure: Constraining $Λ$CDM with BOSS

We reanalyse the anisotropic galaxy clustering measurement from the Baryon Oscillation Spectroscopic Survey (BOSS), demonstrating that using the full shape information provides cosmological constraints that are comparable to other low-redshift probes. We find $Ω_\mathrm{m} = 0.317^{+0.015}_{-0.019}$, $σ_8 = 0.710\pm 0.049$, and $h = 0.704\pm 0.024$ for flat $Λ$CDM cosmologies using uninformative priors on $Ω_\mathrm{c}h^2$, $100θ_\mathrm{MC}$, $\ln 10^{10} A_{s}$, and $n_{s}$, and a prior on $Ω_\mathrm{b}h^2$ that is much wider than current constraints. We quantify the agreement between the Planck 2018 constraints from the cosmic microwave background and BOSS, finding the two data sets to be consistent within a flat $Λ$CDM cosmology using the Bayes factor as well as the prior-insensitive suspiciousness statistic. Combining two low-redshift probes, we jointly analyse the clustering of BOSS galaxies with weak lensing measurements from the Kilo-Degree Survey (KV450). The combination of BOSS and KV450 improves the measurement by up to 45%, constraining $σ_8 = 0.702\pm 0.029$ and $S_8 = σ_8\sqrt{Ω_\mathrm{m}/0.3} = 0.728\pm 0.026$. Over the full 5D parameter space, the odds in favour of single cosmology describing galaxy clustering, lensing, and the cosmic microwave background are $7\pm2$. The suspiciousness statistic signals a $2.1\pm0.3σ$ tension between the combined low-redshift probes and measurements from the cosmic microwave background.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Quantifying Suspiciousness Within Correlated Data Sets

We propose a principled Bayesian method for quantifying tension between correlated datasets with wide uninformative parameter priors. This is achieved by extending the Suspiciousness statistic, which is insensitive to priors. Our method uses global summary statistics, and as such it can be used as a diagnostic for internal consistency. We show how our approach can be combined with methods that use parameter space and data space to identify the existing internal discrepancies. As an example, we use it to test the internal consistency of the KiDS-450 data in 4 photometric redshift bins, and to recover controlled internal discrepancies in simulated KiDS data. We propose this as a diagnostic of internal consistency for present and future cosmological surveys, and as a tension metric for data sets that have non-negligible correlation, such as LSST and Euclid.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

CFHTLenS revisited: assessing concordance with Planck including astrophysical systematics

We investigate the impact of astrophysical systematics on cosmic shear cosmological parameter constraints from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), and the concordance with cosmic microwave background measurements by Planck. We present updated CFHTLenS cosmic shear tomography measurements extended to degree scales using a covariance calibrated by a new suite of N-body simulations. We analyze these measurements with a new model fitting pipeline, accounting for key systematic uncertainties arising from intrinsic galaxy alignments, baryonic effects in the nonlinear matter power spectrum, and photometric redshift uncertainties. We examine the impact of the systematic degrees of freedom on the cosmological parameter constraints, both independently and jointly. When the systematic uncertainties are considered independently, the intrinsic alignment amplitude is the only degree of freedom that is substantially preferred by the data. When the systematic uncertainties are considered jointly, there is no consistently strong preference in favor of the more complex models. We quantify the level of concordance between the CFHTLenS and Planck datasets by employing two distinct data concordance tests, grounded in Bayesian evidence and information theory. We find that the two data concordance tests largely agree with one another, and that the level of concordance between the CFHTLenS and Planck datasets is sensitive to the exact details of the systematic uncertainties included in our analysis, ranging from decisive discordance to substantial concordance as the treatment of the systematic uncertainties becomes more conservative. The least conservative scenario is the one most favored by the cosmic shear data, but it is also the one that shows the greatest degree of discordance with Planck. The data and analysis code are public at https://github.com/sjoudaki/cfhtlens_revisited

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cluster mass profile reconstruction with size and flux magnification on the HST STAGES survey

We present the first measurement of individual cluster mass estimates using weak lensing size and flux magnification. Using data from the HST-STAGES survey of the A901/902 supercluster we detect the four known groups in the supercluster at high significance using magnification alone. We discuss the application of a fully Bayesian inference analysis, and investigate a broad range of potential systematics in the application of the method. We compare our results to a previous weak lensing shear analysis of the same field finding the recovered signal-to-noise of our magnification-only analysis to range from 45% to 110% of the signal-to-noise in the shear-only analysis. On a case-by-case basis we find consistent magnification and shear constraints on cluster virial radius, and finding that for the full sample, magnification constraints to be a factor $0.77 \pm 0.18$ lower than the shear measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Galaxy alignments: An overview

The alignments between galaxies, their underlying matter structures, and the cosmic web constitute vital ingredients for a comprehensive understanding of gravity, the nature of matter, and structure formation in the Universe. We provide an overview on the state of the art in the study of these alignment processes and their observational signatures, aimed at a non-specialist audience. The development of the field over the past one hundred years is briefly reviewed. We also discuss the impact of galaxy alignments on measurements of weak gravitational lensing, and discuss avenues for making theoretical and observational progress over the coming decade.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Galaxy alignments: Observations and impact on cosmology

Galaxy shapes are not randomly oriented, rather they are statistically aligned in a way that can depend on formation environment, history and galaxy type. Studying the alignment of galaxies can therefore deliver important information about the physics of galaxy formation and evolution as well as the growth of structure in the Universe. In this review paper we summarise key measurements of galaxy alignments, divided by galaxy type, scale and environment. We also cover the statistics and formalism necessary to understand the observations in the literature. With the emergence of weak gravitational lensing as a precision probe of cosmology, galaxy alignments have taken on an added importance because they can mimic cosmic shear, the effect of gravitational lensing by large-scale structure on observed galaxy shapes. This makes galaxy alignments, commonly referred to as intrinsic alignments, an important systematic effect in weak lensing studies. We quantify the impact of intrinsic alignments on cosmic shear surveys and finish by reviewing practical mitigation techniques which attempt to remove contamination by intrinsic alignments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Galaxy alignments: Theory, modelling and simulations

The shapes of galaxies are not randomly oriented on the sky. During the galaxy formation and evolution process, environment has a strong influence, as tidal gravitational fields in the large-scale structure tend to align nearby galaxies. Additionally, events such as galaxy mergers affect the relative alignments of both the shapes and angular momenta of galaxies throughout their history. These "intrinsic galaxy alignments" are known to exist, but are still poorly understood. This review will offer a pedagogical introduction to the current theories that describe intrinsic galaxy alignments, including the apparent difference in intrinsic alignment between early- and late-type galaxies and the latest efforts to model them analytically. It will then describe the ongoing efforts to simulate intrinsic alignments using both N-body and hydrodynamic simulations. Due to the relative youth of this field, there is still much to be done to understand intrinsic galaxy alignments and this review summarises the current state of the field, providing a solid basis for future work.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gaussianisation for fast and accurate inference from cosmological data

We present a method to transform multivariate unimodal non-Gaussian posterior probability densities into approximately Gaussian ones via non-linear mappings, such as Box--Cox transformations and generalisations thereof. This permits an analytical reconstruction of the posterior from a point sample, like a Markov chain, and simplifies the subsequent joint analysis with other experiments. This way, a multivariate posterior density can be reported efficiently, by compressing the information contained in MCMC samples. Further, the model evidence integral (i.e. the marginal likelihood) can be computed analytically. This method is analogous to the search for normal parameters in the cosmic microwave background, but is more general. The search for the optimally Gaussianising transformation is performed computationally through a maximum-likelihood formalism; its quality can be judged by how well the credible regions of the posterior are reproduced. We demonstrate that our method outperforms kernel density estimates in this objective. Further, we select marginal posterior samples from Planck data with several distinct strongly non-Gaussian features, and verify the reproduction of the marginal contours. To demonstrate evidence computation, we Gaussianise the joint distribution of data from weak lensing and baryon acoustic oscillations (BAO), for different cosmological models, and find a preference for flat $Λ$CDM. Comparing to values computed with the Savage-Dickey density ratio, and Population Monte Carlo, we find good agreement of our method within the spread of the other two.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A direct measurement of tomographic lensing power spectra from CFHTLenS

We measure the weak gravitational lensing shear power spectra and their cross-power in two photometric redshift bins from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). The measurements are performed directly in multipole space in terms of adjustable band powers. For the extraction of the band powers from the data we have implemented and extended a quadratic estimator, a maximum likelihood method that allows us to readily take into account irregular survey geometries, masks, and varying sampling densities. We find the 68 per cent credible intervals in the $σ_8$-$Ω_{\rm m}$-plane to be marginally consistent with results from $Planck$ for a simple five parameter $Λ$CDM model. For the projected parameter $S_8 \equiv σ_8(Ω_{\rm m}/0.3)^{0.5}$ we obtain a best-fitting value of $S_8 = 0.768_{-0.039}^{+0.045}$. This constraint is consistent with results from other CFHTLenS studies as well as the Dark Energy Survey. Our most conservative model, including modifications to the power spectrum due to baryon feedback and marginalization over photometric redshift errors, yields an upper limit on the total mass of three degenerate massive neutrinos of $Σm_ν< 4.53 \, {\rm eV}$ at 95 per cent credibility, while a Bayesian model comparison does not favour any model extension beyond a simple five parameter $Λ$CDM model. Combining the shear likelihood with $Planck$ breaks the $σ_8$-$Ω_{\rm m}$-degeneracy and yields $σ_8=0.818 \pm 0.013$ and $Ω_{\rm m} = 0.300 \pm 0.011$ which is fully consistent with results from $Planck$ alone.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cross correlation surveys with the Square Kilometre Array

By the time that the first phase of the Square Kilometre Array is deployed it will be able to perform state of the art Large Scale Structure (LSS) as well as Weak Gravitational Lensing (WGL) measurements of the distribution of matter in the Universe. In this chapter we concentrate on the synergies that result from cross-correlating these different SKA data products as well as external correlation with the weak lensing measurements available from CMB missions. We show that the Dark Energy figures of merit obtained individually from WGL/LSS measurements and their independent combination is significantly increased when their full cross-correlations are taken into account. This is due to the increased knowledge of galaxy bias as a function of redshift as well as the extra information from the different cosmological dependences of the cross-correlations. We show that the cross-correlation between a spectroscopic LSS sample and a weak lensing sample with photometric redshifts can calibrate these same photometric redshifts, and their scatter, to high accuracy by modelling them as nuisance parameters and fitting them simultaneously cosmology. Finally we show that Modified Gravity parameters are greatly constrained by this cross-correlations because weak lensing and redshift space distortions (from the LSS survey) break strong degeneracies in common parameterisations of modified gravity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Detection of spatial correlations of fundamental plane residuals, and cosmological implications

The fundamental plane (FP) is a widely used tool to investigate the properties of early-type galaxies, and the tight relation between its parameters has spawned several cosmological applications, including its use as a distance indicator for peculiar velocity surveys and as a means to suppress intrinsic noise in cosmic size magnification measurements. Systematic trends with the large-scale structure across the FP could cause serious biases for these cosmological probes, but may also yield new insights into the early-type population. Here we report the first detection of spatial correlations among offsets in galaxy size from an FP that explicitly accounts for redshift trends, using a sample of about $95,000$ elliptical galaxies from the Sloan Digital Sky Survey. We show that these offsets correlate with the density field out to at least $10h^{-1}$Mpc at $4σ$ significance in a way that cannot be explained by systematic errors in galaxy size estimates. We propose a physical explanation for the correlations by dividing the sample into central, satellite, and field galaxies, identifying trends for each galaxy type separately. Central (satellite) galaxies lie on average above (below) the FP, which we argue could be due to a higher (lower) than average mass-to-light ratio. We fit a simple model to the correlations of FP residuals and use it to predict the impact on peculiar velocity power spectra, finding a contamination larger than $10\,\%$ for $k>0.04\,h/$Mpc. Moreover, cosmic magnification measurements based on an FP could be severely contaminated over a wide range of scales by the intrinsic FP correlations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Enhancing the Cosmic Shear Power Spectrum

Applying a transformation to a non-Gaussian field can enhance the information content of the resulting power spectrum, by reducing the correlations between Fourier modes. In the context of weak gravitational lensing, it has been shown that this gain in information content is significantly compromised by the presence of shape noise. We apply clipping to mock convergence fields, a technique which is known to be robust in the presence of noise and has been successfully applied to galaxy number density fields. When analysed in isolation the resulting convergence power spectrum returns degraded constraints on cosmological parameters. However substantial gains can be achieved by performing a combined analysis of the power spectra derived from both the original and transformed fields. Even in the presence of realistic levels of shape noise, we demonstrate that this approach is capable of reducing the area of likelihood contours within the $Ω_m - σ_8$ plane by more than a factor of three.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

RCSLenS: Testing gravitational physics through the cross-correlation of weak lensing and large-scale structure

The unknown nature of dark energy motivates continued cosmological tests of large-scale gravitational physics. We present a new consistency check based on the relative amplitude of non-relativistic galaxy peculiar motions, measured via redshift-space distortion, and the relativistic deflection of light by those same galaxies traced by galaxy-galaxy lensing. We take advantage of the latest generation of deep, overlapping imaging and spectroscopic datasets, combining the Red Cluster Sequence Lensing Survey (RCSLenS), the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), the WiggleZ Dark Energy Survey and the Baryon Oscillation Spectroscopic Survey (BOSS). We quantify the results using the "gravitational slip" statistic E_G, which we estimate as 0.48 +/- 0.10 at z=0.32 and 0.30 +/- 0.07 at z=0.57, the latter constituting the highest redshift at which this quantity has been determined. These measurements are consistent with the predictions of General Relativity, for a perturbed Friedmann-Robertson-Walker metric in a Universe dominated by a cosmological constant, which are E_G = 0.41 and 0.36 at these respective redshifts. The combination of redshift-space distortion and gravitational lensing data from current and future galaxy surveys will offer increasingly stringent tests of fundamental cosmology.

0 citations0 reviewsNo code linkedOpen access
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)

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Errors on errors - Estimating cosmological parameter covariance

Current and forthcoming cosmological data analyses share the challenge of huge datasets alongside increasingly tight requirements on the precision and accuracy of extracted cosmological parameters. The community is becoming increasingly aware that these requirements not only apply to the central values of parameters but, equally important, also to the error bars. Due to non-linear effects in the astrophysics, the instrument, and the analysis pipeline, data covariance matrices are usually not well known a priori and need to be estimated from the data itself, or from suites of large simulations. In either case, the finite number of realisations available to determine data covariances introduces significant biases and additional variance in the errors on cosmological parameters in a standard likelihood analysis. Here, we review recent work on quantifying these biases and additional variances and discuss approaches to remedy these effects.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Estimating Cosmological Parameter Covariance

We investigate the bias and error in estimates of the cosmological parameter covariance matrix, due to sampling or modelling the data covariance matrix, for likelihood width and peak scatter estimators. We show that these estimators do not coincide unless the data covariance is exactly known. For sampled data covariances, with Gaussian distributed data and parameters, the parameter covariance matrix estimated from the width of the likelihood has a Wishart distribution, from which we derive the mean and covariance. This mean is biased and we propose an unbiased estimator of the parameter covariance matrix. Comparing our analytic results to a numerical Wishart sampler of the data covariance matrix we find excellent agreement. An accurate ansatz for the mean parameter covariance for the peak scatter estimator is found, and we fit its covariance to our numerical analysis. The mean is again biased and we propose an unbiased estimator for the peak parameter covariance. For sampled data covariances the width estimator is more accurate than the peak scatter estimator. We investigate modelling the data covariance, or equivalently data compression, and shown that the peak scatter estimator is less sensitive to biases in the model data covariance matrix than the width estimator, but requires independent realisations of the data to reduce the statistical error. If the model bias on the peak estimator is sufficiently low this is promising, otherwise the sampled width estimator is preferable.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the Probability Distributions of Ellipticity

In this paper we derive an exact full expression for the 2D probability distribution of the ellipticity of an object measured from data, only assuming Gaussian noise in pixel values. This is a generalisation of the probability distribution for the ratio of single random variables, that is well-known, to the multivariate case. This expression is derived within the context of the measurement of weak gravitational lensing from noisy galaxy images. We find that the third flattening, or epsilon-ellipticity, has a biased maximum likelihood but an unbiased mean; and that the third eccentricity, or normalised polarisation chi, has both a biased maximum likelihood and a biased mean. The very fact that the bias in the ellipticity is itself a function of the ellipticity requires an accurate knowledge of the intrinsic ellipticity distribution of the galaxies in order to properly calibrate shear measurements. We use this expression to explore strategies for calibration of biases caused by measurement processes in weak gravitational lensing. We find that upcoming weak lensing surveys like KiDS or DES require calibration fields of order of several square degrees and 1.2 magnitude deeper than the wide survey in order to correct for the noise bias. Future surveys like Euclid will require calibration fields of order 40 square degree and several magnitude deeper than the wide survey. We also investigate the use of the Stokes parameters to estimate the shear as an alternative to the ellipticity. We find that they can provide unbiased shear estimates at the cost of a very large variance in the measurement. The python code used to compute the distributions presented in the paper and to perform the numerical calculations are available on request.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the complementarity of galaxy clustering with cosmic shear and flux magnification

In this paper, we motivate the use of galaxy clustering measurements using photometric redshift information, including a contribution from flux magnification, as a probe of cosmology. We present cosmological forecasts when clustering data alone is used, and when clustering is combined with a cosmic shear analysis. We consider two types of clustering analysis: firstly, clustering with only redshift auto-correlations in tomographic redshift bins; secondly, using all available redshift bin correlations. Finally, we consider how inferred cosmological parameters may be biased using each analysis when flux magnification is neglected. Results are presented for a Stage III ground-based survey, and a Stage IV space-based survey modelled with photometric redshift errors, and values for the slope of the luminosity function inferred from CFHTLenS catalogues. We find that combining clustering information with shear can improve constraints on cosmological parameters, giving an improvement to a Dark Energy Task Force-like figure of merit by a factor of 1.33 when only auto-correlations in redshift are used for the clustering analysis, rising to 1.52 when cross-correlations in redshift are also included. The addition of galaxy-galaxy lensing gives further improvement, with increases in figure of merit by a factor of 2.82 and 3.7 for each type of clustering analysis respectively. The presence of flux magnification in a clustering analysis does not significantly affect the precision of cosmological constraints when combined with cosmic shear and galaxy-galaxy lensing. However if magnification is neglected, inferred cosmological parameter values are biased, with biases in some cosmological parameters larger than statistical errors. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the shear estimation bias induced by the spatial variation of colour across galaxy profiles

The spatial variation of the colour of a galaxy may introduce a bias in the measurement of its shape if the PSF profile depends on wavelength. We study how this bias depends on the properties of the PSF and the galaxies themselves. The bias depends on the scales used to estimate the shape, which may be used to optimise methods to reduce the bias. Here we develop a general approach to quantify the bias. Although applicable to any weak lensing survey, we focus on the implications for the ESA Euclid mission. Based on our study of synthetic galaxies we find that the bias is a few times 10^-3 for a typical galaxy observed by Euclid. Consequently, it cannot be neglected and needs to be accounted for. We demonstrate how one can do so using spatially resolved observations of galaxies in two filters. We show that HST observations in the F606W and F814W filters allow us to model and reduce the bias by an order of magnitude, sufficient to meet Euclid's scientific requirements. The precision of the correction is ultimately determined by the number of galaxies for which spatially-resolved observations in at least two filters are available. We use results from the Millennium Simulation to demonstrate that archival HST data will be sufficient for the tomographic cosmic shear analysis with the Euclid dataset.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Putting the Precision in Precision Cosmology: How accurate should your data covariance matrix be?

Cosmological parameter estimation requires that the likelihood function of the data is accurately known. Assuming that cosmological large-scale structure power spectra data are multivariate Gaussian-distributed, we show the accuracy of parameter estimation is limited by the accuracy of the inverse data covariance matrix - the precision matrix. If the data covariance and precision matrices are estimated by sampling independent realisations of the data, their statistical properties are described by the Wishart and Inverse-Wishart distributions, respectively. Independent of any details of the survey, we show that the fractional error on a parameter variance, or a Figure-of-Merit, is equal to the fractional variance of the precision matrix. In addition, for the only unbiased estimator of the precision matrix, we find that the fractional accuracy of the parameter error depends only on the difference between the number of independent realisations and the number of data points, and so can easily diverge. For a 5% error on a parameter error and N_D << 100 data-points, a minimum of 200 realisations of the survey are needed, with 10% accuracy for the data covariance. If the number of data-points N_D >>100 we need N_S > N_D realisations and a fractional accuracy of <sqrt[2/N_D] in the data covariance. As the number of power spectra data points grows to N_D>10^4 -10^6 this approach will be problematic. We discuss possible ways to relax these conditions: improved theoretical modelling; shrinkage methods; data-compression; simulation and data resampling methods.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Cosmic magnification: nulling the intrinsic clustering signal

We investigate the extent to which the pure magnification effect of gravitational lensing can be extracted from galaxy clustering statistics, by a nulling method which aims to eliminate terms arising from the intrinsic clustering of galaxies. The aim is to leave statistics which are free from the uncertainties of galaxy bias. We find that nulling can be done effectively, leaving data which are relatively insensitive to uncertainties in galaxy bias and its evolution, leading to cosmological parameter estimation which is effectively unbiased. This advantage comes at the expense of increased statistical errors, which are in some cases large, but it offers a robust alternative analysis method to cosmic shear for cosmological imaging surveys designed for weak lensing studies, or to full modelling of the clustering signal including magnification effects.

0 citations0 reviewsNo code linkedOpen access
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)

0 citations0 reviewsNo code linkedOpen access