Source author record

Henk Hoekstra

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

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

92 published item(s)

preprint2025arXiv

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

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

preprint2025arXiv

Quantitative Morphology of Galactic Cirrus in Deep Optical Imaging

Imaging of optical Galactic cirrus, the spatially resolved form of diffuse Galactic light, provides important insights into the properties of the diffuse interstellar medium (ISM) in the Milky Way. While previous investigations have focused mainly on the intensity characteristics of optical cirrus, their morphological properties remain largely unexplored. In this study, we employ several complementary statistical approaches -- local intensity statistics, angular power spectrum / $Δ$-variance analysis, and wavelet scattering transform analysis -- to characterize the morphology of cirrus in deep optical imaging data. We place our investigation of optical cirrus into a multi-wavelength context by comparing the morphology of cirrus seen with the Dragonfly Telephoto Array to that seen with space-based facilities working at longer wavelengths (Herschel 250 $μm$, WISE 12 $μm$, and Planck radiance), as well as with structures seen in the DHIGLS HI column density map. Our statistical methods quantify the similarities and the differences of cirrus morphology in all these datasets. The morphology of cirrus at visible wavelengths resembles that of far-infrared cirrus more closely than that of mid-infrared cirrus; on small scales, anisotropies in the cosmic infrared background and systematics may lead to differences. Across all dust tracers, cirrus morphology can be well described by a power spectrum with a common power-law index $γ\sim-2.9$. We demonstrate quantitatively that optical cirrus exhibits filamentary, coherent structures across a broad range of angular scales. Our results offer promising avenues for linking the analysis of coherent structures in optical cirrus to the underlying physical processes in the ISM that shape them. Furthermore, we demonstrate that these morphological signatures can be leveraged to distinguish and disentangle cirrus from extragalactic light.

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.

preprint2022arXiv

Galaxy cluster aperture masses are more robust to baryonic effects than 3D halo masses

Systematic uncertainties in the mass measurement of galaxy clusters limit the cosmological constraining power of future surveys that will detect more than $10^5$ clusters. Previously, we argued that aperture masses can be inferred more accurately and precisely than 3D masses without loss of cosmological constraining power. Here, we use the Baryons and Haloes of Massive Systems (BAHAMAS) cosmological, hydrodynamical simulations to show that aperture masses are also less sensitive to changes in mass caused by galaxy formation processes. For haloes with $m_\mathrm{200m,dmo} > 10^{14} \, h^{-1} \, \mathrm{M}_\odot$, binned by their 3D halo mass, baryonic physics affects aperture masses and 3D halo masses similarly when measured within apertures similar to the halo virial radius, reaching a maximum reduction of $\approx 3 \, \%$. For lower-mass haloes, $10^{13.5} < m_\mathrm{200m,dmo} / (h^{-1} \, \mathrm{M}_\odot) < 10^{14}$, and aperture sizes $\sim 1 \, h^{-1} \, \mathrm{cMpc}$, representative of weak lensing observations, the aperture mass is consistently reduced less ($\lesssim 5 \, \%$) than the 3D halo mass ($\lesssim 10 \, \%$ for $m_\mathrm{200m}$). The halo mass reduction evolves only slightly, by up to $2$ percentage points, between redshift 0.25 and 1 for both the aperture mass and $m_\mathrm{200m}$. Varying the strength of the simulated feedback so the mean simulated hot gas fraction covers the observed scatter inferred from X-ray observations, we find that the aperture mass is consistently less biased than the 3D halo mass, by up to $2 \, $ percentage points at $m_\mathrm{200m,dmo} = 10^{14} \, h^{-1} \, \mathrm{M}_\odot$. Therefore, cluster aperture mass calibrations provide a fruitful path forward for future cluster surveys to reduce their sensitivity to systematic uncertainties.

preprint2022arXiv

How baryons can significantly bias cluster count cosmology

We quantify two main pathways through which baryonic physics biases cluster count cosmology. We create mock cluster samples that reproduce the baryon content inferred from X-ray observations. We link clusters to their counterparts in a dark matter-only universe, whose abundances can be predicted robustly, by assuming the dark matter density profile is not significantly affected by baryons. We derive weak lensing halo masses and infer the best-fitting cosmological parameters $Ω_\mathrm{m}$, $S_8=σ_8(Ω_\mathrm{m}/0.3)^{0.2}$, and $w_0$ from the mock cluster sample. We find that because of the need to accommodate the change in the density profile due to the ejection of baryons, weak lensing mass calibrations are only unbiased if the concentration is left free when fitting the reduced shear with NFW profiles. However, even unbiased total mass estimates give rise to biased cosmological parameters if the measured mass functions are compared with predictions from dark matter-only simulations. This bias dominates for haloes with $m_\mathrm{500c} < 10^{14.5} \, h^{-1} \, \mathrm{M_\odot}$. For a stage IV-like cluster survey without mass estimation uncertainties, an area $\approx 15000 \, \mathrm{deg^2}$ and a constant mass cut of $m_\mathrm{200m,min} = 10^{14} \, h^{-1} \, \mathrm{M_\odot}$, the biases are $-11 \pm 1 \, \%$ in $Ω_\mathrm{m}$, $-3.29 \pm 0.04 \, \%$ in $S_8$, and $9 \pm 1.5 \, \%$ in $w_0$. The statistical significance of the baryonic bias depends on how accurately the actual uncertainty on individual cluster mass estimates is known. We suggest that rather than the total halo mass, the (re-scaled) dark matter mass inferred from the combination of weak lensing and observations of the hot gas, should be used for cluster count cosmology.

preprint2022arXiv

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

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

preprint2022arXiv

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

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

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.

preprint2022arXiv

Why are we still using 3D masses for cluster cosmology?

The abundance of clusters of galaxies is highly sensitive to the late-time evolution of the matter distribution, since clusters form at the highest density peaks. However, the 3D cluster mass cannot be inferred without deprojecting the observations, introducing model-dependent biases and uncertainties due to the mismatch between the assumed and the true cluster density profile and the neglected matter along the sightline. Since projected aperture masses can be measured directly in simulations and observationally through weak lensing, we argue that they are better suited for cluster cosmology. Using the Mira-Titan suite of gravity-only simulations, we show that aperture masses correlate strongly with 3D halo masses, albeit with large intrinsic scatter due to the varying matter distribution along the sightline. Nonetheless, aperture masses can be measured $\approx 2-3$ times more precisely from observations, since they do not require assumptions about the density profile and are only affected by the shape noise in the weak lensing measurements. We emulate the cosmology dependence of the aperture mass function directly with a Gaussian process. Comparing the cosmology sensitivity of the aperture mass function and the 3D halo mass function for a fixed survey solid angle and redshift interval, we find the aperture mass sensitivity is higher for $Ω_\mathrm{m}$ and $w_a$, similar for $σ_8$, $n_\mathrm{s}$, and $w_0$, and slightly lower for $h$. With a carefully calibrated aperture mass function emulator, cluster cosmology analyses can use cluster aperture masses directly, reducing the sensitivity to model-dependent mass calibration biases and uncertainties.

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.

preprint2021arXiv

Mitigating the effects of undersampling in weak lensing shear estimation with metacalibration

Metacalibration is a state-of-the-art technique for measuring weak gravitational lensing shear from well-sampled galaxy images. We investigate the accuracy of shear measured with metacalibration from fitting elliptical Gaussians to undersampled galaxy images. In this case, metacalibration introduces aliasing effects leading to an ensemble multiplicative shear bias about 0.01 for Euclid, and even larger for the Roman Space Telescope, well exceeding the missions' requirements. We find that this aliasing bias can be mitigated by computing shapes from weighted moments with wider Gaussians as weight functions, thereby trading bias for a slight increase in variance of the measurements. We show that this approach is robust to the point-spread function in consideration and meets the stringent requirements of Euclid for galaxies with moderate to high signal-to-noise ratios. We therefore advocate metacalibration as a viable shear measurement option for weak lensing from upcoming space missions.

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.

preprint2020arXiv

Accounting for object detection bias in weak gravitational lensing studies

Weak lensing by large-scale structure is a powerful probe of cosmology if the apparent alignments in the shapes of distant galaxies can be accurately measured. Most studies have therefore focused on improving the fidelity of the shape measurements themselves, but the preceding step of object detection has been largely ignored. In this paper we study the impact of object detection for a Euclid-like survey and show that it leads to biases that exceed requirements for the next generation of cosmic shear surveys. In realistic scenarios, blending of galaxies is an important source of detection bias. We find that MetaDetection is able to account for blending, leading to average multiplicative biases that meet requirements for Stage IV surveys, provided a sufficiently accurate model for the point spread function is available. Further work is needed to estimate the performance for actual surveys. Combined with sufficiently realistic image simulations, this provides a viable way forward towards accurate shear estimates for Stage IV surveys.

preprint2020arXiv

GAMA+KiDS: Empirical correlations between halo mass and other galaxy properties near the knee of the stellar-to-halo mass relation

We use KiDS weak lensing data to measure variations in mean halo mass as a function of several key galaxy properties (namely: stellar colour, specific star formation rate, Sersic index, and effective radius) for a volume-limited sample of GAMA galaxies in a narrow stellar mass range ($M_* \sim 2$--$5 \times 10^{10}$ Msol). This mass range is particularly interesting, inasmuch as it is where bimodalities in galaxy properties are most pronounced, and near to the break in both the galaxy stellar mass function and the stellar-to-halo mass relation (SHMR). In this narrow mass range, we find that both size and Sersic index are better predictors of halo mass than either colour or SSFR, with the data showing a slight preference for Sersic index. In other words, we find that mean halo mass is more tightly correlated with galaxy structure than either past star formation history or current star formation rate. Our results lead to an approximate lower bound on the dispersion in halo masses among $\log M_* \approx {10.5}$ galaxies: we find that the dispersion is $\gtrsim 0.3$ dex. This would imply either that offsets from the mean SHMR are closely coupled to size/structure, or that the dispersion in the SHMR is larger than past results have suggested. Our results thus provide new empirical constraints on the relationship between stellar and halo mass assembly at this particularly interesting mass range.

preprint2020arXiv

KiDS-1000 catalogue: Weak gravitational lensing shear measurements

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

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

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.

preprint2020arXiv

The PAU Survey: Background light estimation with deep learning techniques

In any imaging survey, measuring accurately the astronomical background light is crucial to obtain good photometry. This paper introduces BKGnet, a deep neural network to predict the background and its associated error. BKGnet has been developed for data from the Physics of the Accelerating Universe Survey (PAUS), an imaging survey using a 40 narrow-band filter camera (PAUCam). Images obtained with PAUCam are affected by scattered light: an optical effect consisting of light multiply that deposits energy in specific detector regions contaminating the science measurements. Fortunately, scattered light is not a random effect, but it can be predicted and corrected for. We have found that BKGnet background predictions are very robust to distorting effects, while still being statistically accurate. On average, the use of BKGnet improves the photometric flux measurements by 7% and up to 20% at the bright end. BKGnet also removes a systematic trend in the background error estimation with magnitude in the i-band that is present with the current PAU data management method. With BKGnet, we reduce the photometric redshift outlier rate

preprint2020arXiv

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

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

preprint2019arXiv

CCCP and MENeaCS: (updated) weak-lensing masses for 100 galaxy clusters

Large area surveys have detected significant samples of galaxy clusters that can be used to constrain cosmological parameters, provided that the masses of the clusters are measured robustly. To improve the calibration of cluster masses using weak gravitational lensing we present new results for 48 clusters at $0.05<z<0.15$, observed as part of the Multi Epoch Nearby Cluster Survey (MENeaCS), and reevaluate the mass estimates for 52 clusters from the Canadian Cluster Comparison Project (CCCP). Updated high-fidelity photometric redshift catalogues of reference deep fields are used in combination with advances in shape measurements and state-of-the-art cluster simulations, yielding an average systematic uncertainty in the lensing signal below 5%, similar to the statistical uncertainty for our cluster sample. We derive a scaling relation with Planck measurements for the full sample and find a bias in the Planck masses of $1-b=0.84 \pm 0.04$. We find no statistically significant trend of the mass bias with redshift or cluster mass, but find that different selections could change the bias by up to 1.5$σ$. We find a gas fraction of $0.139 \pm 0.014$ for 8 relaxed clusters in our sample, which can also be used to infer cosmological parameters.

preprint2019arXiv

The impact of the observed baryon distribution in haloes on the total matter power spectrum

The interpretation of upcoming weak gravitational lensing surveys depends critically on our understanding of the matter power spectrum on scales $k < 10 h/\mathrm{Mpc}$, where baryonic processes are important. We study the impact of galaxy formation processes on the matter power spectrum using a halo model that treats the stars and gas separately from the dark matter distribution. We use empirical constraints from X-ray observations (hot gas) and halo occupation distribution modelling (stars) for the baryons. Since X-ray observations cannot generally measure the hot gas content outside $r_\mathrm{500c}$, we vary the gas density profiles beyond this radius. Compared with dark matter only models, we find a total power suppression of $1 \%$ ($5 \%$) on scales $0.2-1 h/\mathrm{Mpc}$ ($0.5-2h/\mathrm{Mpc}$), where lower baryon fractions result in stronger suppression. We show that groups of galaxies ($10^{13} < m_\mathrm{500c} / (\mathrm{M_\odot}/h) < 10^{14}$) dominate the total power at all scales $k \lesssim 10 h/\mathrm{Mpc}$. We find that a halo mass bias of $30 \%$ (similar to what is expected from the hydrostatic equilibrium assumption) results in an underestimation of the power suppression of up to $4 \%$ at $k = 1 h/\mathrm{Mpc}$, illustrating the importance of measuring accurate halo masses. Contrary to work based on hydrodynamical simulations, our conclusion that baryonic effects can no longer be neglected is not subject to uncertainties associated with our poor understanding of feedback processes. Observationally, probing the outskirts of groups and clusters will provide the tightest constraints on the power suppression for $k \lesssim 1 h/\mathrm{Mpc}$.

preprint2016arXiv

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

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

preprint2016arXiv

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.

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.

preprint2016arXiv

The abundance and spatial distribution of ultra-diffuse galaxies in nearby galaxy clusters

Recent observations have highlighted a significant population of faint but large (r_eff>1.5 kpc) galaxies in the Coma cluster. The origin of these Ultra Diffuse Galaxies (UDGs) remains puzzling, as the interpretation of the observational results has been hindered by the subjective selection of UDGs, and the limited study of only the Coma (and some examples in the Virgo-) cluster. We extend the study of UDGs using 8 clusters in the redshift range 0.044<z<0.063 with deep g- and r-band imaging data taken with MegaCam at the CFHT. We describe an automatic selection pipeline for quantitative identification, tested for completeness using image simulations of these galaxies. We find that the abundance of the UDGs we can detect increases with cluster mass, reaching ~200 in typical haloes of M200~10^15 Msun. The cluster UDGs have colours consistent with the cluster red sequence, and have a steep size distribution that declines as n~r_eff^-3.4. Their radial distribution is significantly steeper than NFW in the outskirts, and is significantly shallower in the inner parts. They follow the same radial distribution as the more massive quiescent galaxies in the clusters, except within the core region of r<0.15XR200 (or <300 kpc). Within this region the number density of UDGs drops and is consistent with zero. These diffuse galaxies can only resist tidal forces down to this cluster-centric distance if they are highly centrally dark-matter dominated. The observation that the radial distribution of more compact dwarf galaxies (r_eff<1.0 kpc) with similar luminosities follows the same distribution as the UDGs, but exist down to a smaller distance of 100kpc from the cluster centres, indicates that they may have similarly massive sub-haloes as the UDGs. Although several scenarios can give rise to the UDG population, our results point to differences in the formation history as the most plausible explanation.

preprint2016arXiv

The evolution in the stellar mass of Brightest Cluster Galaxies over the past 10 billion years

Using a sample of 98 galaxy clusters recently imaged in the near infra-red with the ESO NTT, WIYN and WHT telescopes, supplemented with 33 clusters from the ESO archive, we measure how the stellar mass of the most massive galaxies in the universe, namely Brightest Cluster Galaxies (BCG), increases with time. Most of the BCGs in this new sample lie in the redshift range $0.2<z<0.6$, which has been noted in recent works to mark an epoch over which the growth in the stellar mass of BCGs stalls. From this sample of 132 clusters, we create a subsample of 102 systems that includes only those clusters that have estimates of the cluster mass. We combine the BCGs in this subsample with BCGs from the literature, and find that the growth in stellar mass of BCGs from 10 billion years ago to the present epoch is broadly consistent with recent semi-analytic and semi-empirical models. As in other recent studies, tentative evidence indicates that the stellar mass growth rate of BCGs may be slowing in the past 3.5 billion years. Further work in collecting larger samples, and in better comparing observations with theory using mock images is required if a more detailed comparison between the models and the data is to be made.

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.

preprint2015arXiv

CFHTLenS: Testing the Laws of Gravity with Tomographic Weak Lensing and Redshift Space Distortions

Dark energy may be the first sign of new fundamental physics in the Universe, taking either a physical form or revealing a correction to Einsteinian gravity. Weak gravitational lensing and galaxy peculiar velocities provide complementary probes of General Relativity, and in combination allow us to test modified theories of gravity in a unique way. We perform such an analysis by combining measurements of cosmic shear tomography from the Canada-France Hawaii Telescope Lensing Survey (CFHTLenS) with the growth of structure from the WiggleZ Dark Energy Survey and the Six-degree-Field Galaxy Survey (6dFGS), producing the strongest existing joint constraints on the metric potentials that describe general theories of gravity. For scale-independent modifications to the metric potentials which evolve linearly with the effective dark energy density, we find present-day cosmological deviations in the Newtonian potential and curvature potential from the prediction of General Relativity to be (Delta Psi)/Psi = 0.05 \pm 0.25 and (Delta Phi)/Phi = -0.05 \pm 0.3 respectively (68 per cent CL).

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

Cosmology and fundamental physics with the Euclid satellite

Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

preprint2015arXiv

Evidence for a change in the dominant satellite galaxy quenching mechanism at z=1

We present an analysis of galaxies in groups and clusters at $0.8<z<1.2$, from the GCLASS and GEEC2 spectroscopic surveys. We compute a "conversion fraction" $f_{\rm convert}$ that represents the fraction of galaxies that were prematurely quenched by their environment. For massive galaxies, $M_{\rm star}>10^{10.3}M_\odot$, we find $f_{\rm convert}\sim 0.4$ in the groups and $\sim 0.6$ in the clusters, similar to comparable measurements at $z=0$. This means the time between first accretion into a more massive halo and final star formation quenching is $t_p\sim 2$ Gyr. This is substantially longer than the estimated time required for a galaxy's star formation rate to become zero once it starts to decline, suggesting there is a long delay time during which little differential evolution occurs. In contrast with local observations we find evidence that this delay timescale may depend on stellar mass, with $t_p$ approaching $t_{\rm Hubble}$ for $M_{\rm star}\sim 10^{9.5}M_\odot$. The result suggests that the delay time must not only be much shorter than it is today, but may also depend on stellar mass in a way that is not consistent with a simple evolution in proportion to the dynamical time. Instead, we find the data are well-matched by a model in which the decline in star formation is due to "overconsumption", the exhaustion of a gas reservoir through star formation and expulsion via modest outflows in the absence of cosmological accretion. Dynamical gas removal processes, which are likely dominant in quenching newly accreted satellites today, may play only a secondary role at $z=1$.

preprint2015arXiv

Evidence for the inside-out growth of the stellar mass distribution in galaxy clusters since z~1

We study the radial number density and stellar mass density distributions of satellite galaxies in a sample of 60 massive clusters at 0.04<z<0.26 selected from the Multi-Epoch Nearby Cluster Survey (MENeaCS) and the Canadian Cluster Comparison Project (CCCP). In addition to ~10,000 spectroscopically confirmed member galaxies, we use deep ugri-band imaging to estimate photometric redshifts and stellar masses, and then statistically subtract fore-, and background sources using data from the COSMOS survey. We measure the galaxy number density and stellar mass density distributions in logarithmically spaced bins over 2 orders of magnitude in radial distance from the BCGs. For projected distances in the range 0.1<R/R200<2.0, we find that the stellar mass distribution is well-described by an NFW profile with a concentration of c=2.03+/-0.20. However, at smaller radii we measure a significant excess in the stellar mass in satellite galaxies of about $10^{11}$ Msun per cluster, compared to these NFW profiles. We do obtain good fits to generalized NFW profiles with free inner slopes, and to Einasto profiles. To examine how clusters assemble their stellar mass component over cosmic time, we compare this local sample to the GCLASS cluster sample at z~1, which represents the approximate progenitor sample of the low-z clusters. This allows for a direct comparison, which suggests that the central parts (R<0.4 Mpc) of the stellar mass distributions of satellites in local galaxy clusters are already in place at z~1, and contain sufficient excess material for further BCG growth. Evolving towards z=0, clusters appear to assemble their stellar mass primarily onto the outskirts, making them grow in an inside-out fashion.

preprint2015arXiv

Evolution of the luminosity-to-halo mass relation of LRGs from a combined SDSS-DR10+RCS2 analysis

We study the evolution of the luminosity-to-halo mass relation of Luminous Red Galaxies (LRGs). We select a sample of 52 000 LOWZ and CMASS LRGs from the Baryon Oscillation Spectroscopic Survey (BOSS) SDSS-DR10 in the ~450 deg^2 that overlaps with imaging data from the second Red-sequence Cluster Survey (RCS2), group them into bins of absolute magnitude and redshift and measure their weak lensing signals. The source redshift distribution has a median of 0.7, which allows us to study the lensing signal as a function of lens redshift. We interpret the lensing signal using a halo model, from which we obtain the halo masses as well as the normalisations of the mass-concentration relations. We find that the concentration of haloes that host LRGs is consistent with dark matter only simulations once we allow for miscentering or satellites in the modelling. The slope of the luminosity-to-halo mass relation has a typical value of 1.4 and does not change with redshift, but we do find evidence for a change in amplitude: the average halo mass of LOWZ galaxies increases by 25_{-14}^{+16} % between z=0.36 and 0.22 to an average value of 6.43+/-0.52 x 10^13 h70^-1 Msun. If we extend the redshift range using the CMASS galaxies and assume that they are the progenitors of the LOWZ sample, we find that the average mass of LRGs increases by 80^{+39}_{-28} % between z=0.6 and 0.2

preprint2015arXiv

Gravitational Lensing Analysis of the Kilo Degree Survey

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

preprint2015arXiv

Intrinsic alignments of galaxies in the EAGLE and cosmo-OWLS simulations

We report results for the alignments of galaxies in the EAGLE and cosmo-OWLS simulations as a function of galaxy separation and halo mass. The combination of these hydro-cosmological simulations enables us to span four orders of magnitude in halo mass ($10.7<log_{10}(M_{200}/[h^{-1}M_\odot])<15$) and a large range of separations ($-1<log_{10}(r/[h^{-1}Mpc])< 2$). We focus on two classes of alignments: the orientations of galaxies with respect to either the directions to, or the orientations of, surrounding galaxies. We find that the strength of the alignment is a strongly decreasing function of the distance between galaxies. The orientation-direction alignment can remain significant up to ~100 Mpc, for galaxies hosted by the most massive haloes in our simulations. Galaxies hosted by more massive subhaloes show stronger alignment. At a fixed halo mass, more aspherical or prolate galaxies exhibit stronger alignments. The spatial distribution of satellites is anisotropic and significantly aligned with the major axis of the main host halo. The major axis of satellite galaxies, when all stars are considered, are preferentially aligned towards the centre of the main host halo. The predicted projected direction-orientation alignment, $ε_{g+}(r_{p})$, is in broad agreement with recent observations when only stars within the typical observable extent of a galaxy are used to define galaxy orientations. We find that the orientation-orientation alignment is weaker than the orientation-direction alignment on all scales. Overall, the strength of galaxy alignments depends strongly on the subset of stars that are used to measure the orientations of galaxies and it is always weaker than the alignment of the dark matter haloes. Thus, alignment models that use halo orientation as a direct proxy for galaxy orientation will overestimate the impact of intrinsic alignments on weak lensing analyses.

preprint2015arXiv

Mapping the Galaxy Color-Redshift Relation: Optimal Photometric Redshift Calibration Strategies for Cosmology Surveys

Calibrating the photometric redshifts of >10^9 galaxies for upcoming weak lensing cosmology experiments is a major challenge for the astrophysics community. The path to obtaining the required spectroscopic redshifts for training and calibration is daunting, given the anticipated depths of the surveys and the difficulty in obtaining secure redshifts for some faint galaxy populations. Here we present an analysis of the problem based on the self-organizing map, a method of mapping the distribution of data in a high-dimensional space and projecting it onto a lower-dimensional representation. We apply this method to existing photometric data from the COSMOS survey selected to approximate the anticipated Euclid weak lensing sample, enabling us to robustly map the empirical distribution of galaxies in the multidimensional color space defined by the expected Euclid filters. Mapping this multicolor distribution lets us determine where - in galaxy color space - redshifts from current spectroscopic surveys exist and where they are systematically missing. Crucially, the method lets us determine whether a spectroscopic training sample is representative of the full photometric space occupied by the galaxies in a survey. We explore optimal sampling techniques and estimate the additional spectroscopy needed to map out the color-redshift relation, finding that sampling the galaxy distribution in color space in a systematic way can efficiently meet the calibration requirements. While the analysis presented here focuses on the Euclid survey, similar analysis can be applied to other surveys facing the same calibration challenge, such as DES, LSST, and WFIRST.

preprint2015arXiv

The Canadian Cluster Comparison Project: detailed study of systematics and updated weak lensing masses

Masses of clusters of galaxies from weak gravitational lensing analyses of ever larger samples are increasingly used as the reference to which baryonic scaling relations are compared. In this paper we revisit the analysis of a sample of 50 clusters studied as part of the Canadian Cluster Comparison Project. We examine the key sources of systematic error in cluster masses. We quantify the robustness of our shape measurements and calibrate our algorithm empirically using extensive image simulations. The source redshift distribution is revised using the latest state-of-the-art photometric redshift catalogs that include new deep near-infrared observations. Nonetheless we find that the uncertainty in the determination of photometric redshifts is the largest source of systematic error for our mass estimates. We use our updated masses to determine b, the bias in the hydrostatic mass, for the clusters detected by Planck. Our results suggest 1-b=0.76+-0.05(stat)}+-0.06(syst)}, which does not resolve the tension with the measurements from the primary cosmic microwave background.

preprint2015arXiv

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

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

preprint2015arXiv

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

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

preprint2015arXiv

Weak-lensing-inferred scaling relations of galaxy clusters in the RCS2: mass-richness, mass-concentration, mass-bias, and more

We study a sample of ~10^4 galaxy clusters in the redshift range 0.2<z<0.8 with masses M_200 > 5x10^13 h_70^-1 M_sun, discovered in the second Red-sequence Cluster Survey (RCS2). The depth and excellent image quality of the RCS2 enable us to detect the cluster-mass cross-correlation up to z~0.7. To obtain cluster masses, concentrations and halo biases, we fit a cluster halo model simultaneously to the lensing signal and to the projected density profile of red-sequence cluster members, as the latter provides tight constraints on the cluster miscentring distribution. We parametrise the mass-richness relation as M_200 = A x (N_200/20)^alpha, and find A = (15.0 +- 0.8) x 10^13 h_70^-1 M_sun and alpha = 0.73 +- 0.07 at low redshift (0.2<z<0.35). At intermediate redshift (0.35<z<0.55), we find a higher normalisation, which points at a fractional increase of the richness towards lower redshift caused by the build-up of the red-sequence. The miscentring distribution is well constrained. Only ~30% of our BCGs coincide with the peak of the dark matter distribution. The distribution of the remaining BCGs are modelled with a 2D-Gaussian, whose width increases from 0.2 to 0.4 h_70^-1 Mpc towards higher masses; the ratio of width and r_200 is constant with mass and has an average value of 0.44 +- 0.01. The mass-concentration and mass-bias relation agree fairly well with literature results at low redshift, but have a higher normalisation at higher redshifts, which may be due to selection and projection effects. The concentration of the satellite distribution decreases with mass and is correlated with the concentration of the halo.

preprint2014arXiv

CFHTLenS: A Weak Lensing Shear Analysis of the 3D-Matched-Filter Galaxy Clusters

We present the cluster mass-richness scaling relation calibrated by a weak lensing analysis of >18000 galaxy cluster candidates in the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). Detected using the 3D-Matched-Filter cluster-finder of Milkeraitis et al., these cluster candidates span a wide range of masses, from the small group scale up to $\sim10^{15} M_{\odot}$, and redshifts 0.2 $\lesssim z\lesssim$ 0.9. The total significance of the stacked shear measurement amounts to 54$σ$. We compare cluster masses determined using weak lensing shear and magnification, finding the measurements in individual richness bins to yield 1$σ$ compatibility, but with magnification estimates biased low. This first direct mass comparison yields important insights for improving the systematics handling of future lensing magnification work. In addition, we confirm analyses that suggest cluster miscentring has an important effect on the observed 3D-MF halo profiles, and we quantify this by fitting for projected cluster centroid offsets, which are typically $\sim$ 0.4 arcmin. We bin the cluster candidates as a function of redshift, finding similar cluster masses and richness across the full range up to $z \sim$ 0.9. We measure the 3D-MF mass-richness scaling relation $M_{200} = M_0 (N_{200} / 20)^β$. We find a normalization $M_0 \sim (2.7^{+0.5}_{-0.4}) \times 10^{13} M_{\odot}$, and a logarithmic slope of $β\sim 1.4 \pm 0.1$, both of which are in 1$σ$ agreement with results from the magnification analysis. We find no evidence for a redshift-dependence of the normalization. The CFHTLenS 3D-MF cluster catalogue is now available at cfhtlens.org.

preprint2014arXiv

CFHTLenS: Co-evolution of galaxies and their dark matter haloes

Galaxy-galaxy weak lensing is a direct probe of the mean matter distribution around galaxies. The depth and sky coverage of the CFHT Legacy Survey yield statistically significant galaxy halo mass measurements over a much wider range of stellar masses ($10^{8.75}$ to $10^{11.3} M_{\odot}$) and redshifts ($0.2 < z < 0.8$) than previous weak lensing studies. At redshift $z \sim 0.5$, the stellar-to-halo mass ratio (SHMR) reaches a maximum of $4.0\pm0.2$ percent as a function of halo mass at $\sim 10^{12.25} M_{\odot}$. We find, for the first time from weak lensing alone, evidence for significant evolution in the SHMR: the peak ratio falls as a function of cosmic time from $4.5 \pm 0.3$ percent at $z \sim 0.7$ to $3.4 \pm 0.2$ percent at $z \sim 0.3$, and shifts to lower stellar mass haloes. These evolutionary trends are dominated by red galaxies, and are consistent with a model in which the stellar mass above which star formation is quenched "downsizes" with cosmic time. In contrast, the SHMR of blue, star-forming galaxies is well-fit by a power law that does not evolve with time. This suggests that blue galaxies form stars at a rate that is balanced with their dark matter accretion in such a way that they evolve along the SHMR locus. The redshift dependence of the SHMR can be used to constrain the evolution of the galaxy population over cosmic time.

preprint2014arXiv

CFHTLenS: Cosmological constraints from a combination of cosmic shear two-point and three-point correlations

Higher-order, non-Gaussian aspects of the large-scale structure carry valuable information on structure formation and cosmology, which is complementary to second-order statistics. In this work we measure second- and third-order weak-lensing aperture-mass moments from CFHTLenS and combine those with CMB anisotropy probes. The third moment is measured with a significance of $2σ$. The combined constraint on $Σ_8 = σ_8 (Ω_{\rm m}/0.27)^α$ is improved by 10%, in comparison to the second-order only, and the allowed ranges for $Ω_{\rm m}$ and $σ_8$ are substantially reduced. Including general triangles of the lensing bispectrum yields tighter constraints compared to probing mainly equilateral triangles. Second- and third-order CFHTLenS lensing measurements improve Planck CMB constraints on $Ω_{\rm m}$ and $σ_8$ by 26% for flat $Λ$CDM. For a model with free curvature, the joint CFHTLenS-Planck result is $Ω_{\rm m} = 0.28 \pm 0.02$ (68% confidence), which is an improvement of 43% compared to Planck alone. We test how our results are potentially subject to three astrophysical sources of contamination: source-lens clustering, the intrinsic alignment of galaxy shapes, and baryonic effects. We explore future limitations of the cosmological use of third-order weak lensing, such as the nonlinear model and the Gaussianity of the likelihood function.

preprint2014arXiv

Constraints on the Alignment of Galaxies in Galaxy Clusters from $\sim$14,000 Spectroscopic Members

Torques acting on galaxies lead to physical alignments, but the resulting ellipticity correlations are difficult to predict. As they constitute a major contaminant for cosmic shear studies, it is important to constrain the intrinsic alignment signal observationally. We measured the alignments of satellite galaxies within 90 massive galaxy clusters in the redshift range 0.05<z<0.55 and quantified their impact on the cosmic shear signal. We combined a sample of 38,104 galaxies with spectroscopic redshifts with high-quality data from the Canada-France-Hawaii Telescope. We used phase-space information to select 14,576 cluster members, 14,250 of which have shape measurements and measured three different types of alignment: the radial alignment of satellite galaxies toward the brightest cluster galaxies (BCGs), the common orientations of satellite galaxies and BCGs, and the radial alignments of satellites with each other. Residual systematic effects are much smaller than the statistical uncertainties. We detect no galaxy alignment of any kind out to at least 3 r200. The signal is consistent with zero for both blue and red galaxies, bright and faint ones, and also for subsamples of clusters based on redshift, dynamical mass, and dynamical state. These conclusions are unchanged if we expand the sample with bright cluster members from the red sequence. We augment our constraints with those from the literature to estimate the importance of the intrinsic alignments of satellites compared to those of central galaxies, for which the alignments are described by the linear alignment model. Comparison of the alignment signals to the expected uncertainties of current surveys such as the Kilo-Degree Survey suggests that the linear alignment model is an adequate treatment of intrinsic alignments, but it is not clear whether this will be the case for larger surveys.

preprint2014arXiv

Hubble Space Telescope/Advanced Camera for Surveys Confirmation of the Dark Substructure in A520

We present the results from a weak gravitational lensing study of the merging cluster A520 based on the analysis of Hubble Space Telescope/Advanced Camera for Surveys (ACS) data. The excellent data quality allows us to reach a mean number density of source galaxies of ~109 per sq. arcmin, which improves both resolution and significance of the mass reconstruction compared to a previous study based on Wide Field Planetary Camera 2 (WFPC2) images. We take care in removing instrumental effects such as the trailing of charge due to radiation damage of the ACS detector and the position-dependent point spread function (PSF). This new ACS analysis confirms the previous claims that a substantial amount of dark mass is present between two luminous subclusters. We examine the distribution of cluster galaxies and observe very little light at this location. We find that the centroid of the dark peak in the current ACS analysis is offset to the southwest by ~1 arcmin with respect to the centroid from the WFPC2 analysis. Interestingly, this new centroid is in better spatial agreement with the location where the X-ray emission is strongest, and the mass-to-light ratio estimated with this centroid is much higher 813+-78 M_sun/L_Rsun than the previous value; the aperture mass based on the WFPC2 centroid provides a slightly lower, but consistent mass. Although we cannot provide a definite explanation for the presence of the dark peak, we discuss a revised scenario, wherein dark matter with a more conventional range sigma_DM/m_DM < 1 cm^2/g of self-interacting cross-section can lead to the detection of this dark substructure. If supported by detailed numerical simulations, this hypothesis opens up the possibility that the A520 system can be used to establish a lower limit of the self-interacting cross-section of dark matter.

preprint2014arXiv

MC^2: Constraining the Dark Matter Distribution of the Violent Merging Galaxy Cluster CIZA J2242.8+5301 by Piercing through the Milky Way

CIZA J2242.8+5301 is a merging system with a prominent (~2 Mpc long) radio relic, which together with the morphology of the X-ray emission provides strong evidence for a violent collision along the N-S axis. We present our constraints on the dark matter distribution of this unusual system using Subaru and CFHT imaging data. Measuring a high S/N lensing signal from this cluster is potentially a challenging task because of its proximity to the Milky Way plane (|b|~5 deg). We overcome this challenge with careful observation planning and systematics control, which enables us to successfully map the dark matter distribution of the cluster with high fidelity. The resulting mass map shows that the mass distribution is highly elongated along the N-S merger axis inferred from the orientation of the radio relics. Based on our mass reconstruction, we identify two sub-clusters, which coincide with the cluster galaxy distributions. We determine their masses using MCMC analysis by simultaneously fitting two NFW halos without fixing their centroids. The resulting masses of the northern and southern systems are $M_{200}=11.0_{-3.2}^{+3.7}\times10^{14} M_{\sun}$ and $9.8_{-2.5}^{+3.8}\times10^{14} M_{\sun}$, respectively, indicating that we are witnessing a post-collision of two giant systems of nearly equal mass. When the mass and galaxy centroids are compared in detail, we detect ~ 1' (~190 kpc) offsets in both northern and southern sub-clusters. We find that the galaxy luminosity-mass offset for the northern clump is statistically significant at the ~2 sigma level whereas the detection is only marginal for the southern sub-cluster in part because of a relatively large mass centroid error. We conclude that it is yet premature to uniquely attribute the galaxy-mass misalignment to SIDM and discuss caveats.

preprint2014arXiv

Scaling Relations and X-ray Properties of Moderate-Luminosity Galaxy Clusters from 0.3 < z < 0.6 with XMM-Newton

We present new X-ray temperatures and improved X-ray luminosity estimates for 15 new and archival XMM-Newton observations of galaxy clusters at intermediate redshift with mass and luminosities near the galaxy group/cluster division (M2500 < $2.4\times 10^{14} h_{70}^{-1} M_\odot$, L < $2\times 10^{44}$ erg $s^{-1}$, 0.3< z < 0.6). These clusters have weak-lensing mass measurements based on Hubble Space Telescope observations of clusters representative of an X-ray selected sample (the ROSAT 160SD survey). The angular resolution of XMM-Newton allows us to disentangle the emission of these galaxy clusters from nearby point sources, which significantly contaminated previous X-ray luminosity estimates for six of the fifteen clusters. We extend cluster scaling relations between X-ray luminosity, temperature, and weak-lensing mass for low-mass, X-ray-selected clusters out to redshift $\approx$0.45. These relations are important for cosmology and the astrophysics of feedback in galaxy groups and clusters. Our joint analysis with a sample of 50 clusters in a similar redshift range but with larger masses (M500 < $21.9 \times 10^{14} M_\odot$, $0.15 \leq z \leq 0.55$) from the Canadian Cluster Comparison Project finds that within r2500, $M \propto L^{0.44 +/- 0.05}$, $T \propto L^{0.23 +/- 0.02}$, and $M \propto T^{1.9 +/- 0.2}$. The estimated intrinsic scatter in the M-L relation for the combined sample is reduced to $σ_{log(M|L)}=0.10$, from $σ_{log(M|L)}=0.26$ with the original ROSAT measurements. We also find an intrinsic scatter for the T-L relation, $σ_{log(T|L)}=0.07 +/- 0.01$.

preprint2014arXiv

The Phase Space and Stellar Populations of Cluster Galaxies at z ~ 1: Simultaneous Constraints on the Location and Timescale of Satellite Quenching

We investigate the velocity vs. position phase space of z ~ 1 cluster galaxies using a set of 424 spectroscopic redshifts in 9 clusters drawn from the GCLASS survey. Dividing the galaxy population into three categories: quiescent, star-forming, and poststarburst, we find that these populations have distinct distributions in phase space. Most striking are the poststarburst galaxies, which are commonly found at small clustercentric radii with high clustercentric velocities, and appear to trace a coherent ``ring" in phase space. Using several zoom simulations of clusters we show that the coherent distribution of the poststarbursts can be reasonably well-reproduced using a simple quenching scenario. Specifically, the phase space is best reproduced if satellite quenching occurs on a rapid timescale (0.1 < tau_{Q} < 0.5 Gyr) after galaxies make their first passage of R ~ 0.5R_{200}, a process that takes a total time of ~ 1 Gyr after first infall. We compare this quenching timescale to the timescale implied by the stellar populations of the poststarburst galaxies and find that the poststarburst spectra are well-fit by a rapid quenching (tau_{Q} = 0.4^{+0.3}_{-0.4} Gyr) of a typical star-forming galaxy. The similarity between the quenching timescales derived from these independent indicators is a strong consistency check of the quenching model. Given that the model implies satellite quenching is rapid, and occurs well within R_{200}, this would suggest that ram-pressure stripping of either the hot or cold gas component of galaxies are the most plausible candidates for the physical mechanism. The high cold gas consumption rates at z ~ 1 make it difficult to determine if hot or cold gas stripping is dominant; however, measurements of the redshift evolution of the satellite quenching timescale and location may be capable of distinguishing between the two.

preprint2014arXiv

The rise and fall of star-formation in $\bf z\sim0.2$ merging galaxy clusters

CIZA J2242.8+5301 (`Sausage') and 1RXS J0603.3+4213 (`Toothbrush') are two low-redshift ($z\sim0.2$), massive ($\sim2\times10^{15}M_\odot$), post-core passage merging clusters, which host shock waves traced by diffuse radio emission. To study their star-formation properties, we uniformly survey the `Sausage' and `Toothbrush' clusters in broad and narrow band filters and select a sample of $201$ and $463$ line emitters, down to a rest-frame equivalent width ($13$Å). We robustly separate between H$α$ and higher redshift emitters using a combination of optical multi-band (B, g, V, r, i, z) and spectroscopic data. We build H$α$ luminosity functions for the entire cluster region, near the shock fronts, and away from the shock fronts and find striking differences between the two clusters. In the dynamically younger, $1$ Gyr old `Sausage' cluster we find numerous ($59$) H$α$ emitters above a star-formation rate (SFR) of $0.17$ M$_{\sun}$ yr$^{-1}$ surprisingly located in close proximity to the shock fronts, embedded in very hot intra-cluster medium plasma. The SFR density for the cluster population is at least at the level of typical galaxies at $z\sim2$. Down to the same star-formation rate, the possibly dynamically more evolved `Toothbrush' cluster has only $9$ H$α$ galaxies. The cluster H$α$ galaxies fall on the SFR-stellar mass relation $z\sim0.2$ for the field. However, the `Sausage' cluster has an H$α$ emitter density $>20$ times that of blank fields. If the shock passes through gas-rich cluster galaxies, the compressed gas could collapse into dense clouds and excite star-formation for a few $100$ Myr. This process ultimately leads to a rapid consumption of the molecular gas, accelerating the transformation of gas-rich field spirals into cluster S0s or ellipticals.

preprint2013arXiv

A Census of Stellar Mass in 10 Massive Haloes at z~1 from the GCLASS Survey

We study the stellar mass content of massive haloes in the redshift range 0.86<z<1.34, by measuring: (1) The stellar mass in the central galaxy versus total dynamical halo mass. (2) The total stellar mass (including satellites) versus total halo mass. (3) The radial stellar mass and number density profiles for the ensemble halo. We use a Ks-band selected catalogue for the 10 clusters in the Gemini Cluster Astrophysics Spectroscopic Survey (GCLASS), with photometric redshifts and stellar masses measured from 11-band SED fitting. Combining the photometric catalogues with the deep spectroscopic component of GCLASS, we correct the cluster galaxy sample for interlopers. We also perform a dynamical analysis of the cluster galaxies to estimate the halo mass M200 for each cluster based on a measurement of its velocity dispersion. (1) We find that the central galaxy stellar mass fraction decreases with total halo mass, and that this is in reasonable quantitative agreement with measurements from abundance matching studies at z~1. (2) The total stellar mass fractions of these systems decrease with halo mass, indicating that lower mass systems are more efficient at transforming baryons into stars. We find the total stellar mass to be a good proxy for total halo mass, with a small intrinsic scatter. When we compare these results with literature measurements, we find that the stellar mass fraction at fixed halo mass shows no significant evolution in the range 0<z<1. (3) We measure a relatively high NFW concentration parameter c ~ 7 for the stellar mass distribution in these clusters. A simple model shows that the stellar mass content of GCLASS can evolve in typical distributions observed at lower redshifts if the clusters primarily accrete stellar mass onto the outskirts. [Abridged]

preprint2013arXiv

CFHTLenS tomographic weak lensing cosmological parameter constraints: Mitigating the impact of intrinsic galaxy alignments

We present a finely-binned tomographic weak lensing analysis of the Canada-France-Hawaii Telescope Lensing Survey, CFHTLenS, mitigating contamination to the signal from the presence of intrinsic galaxy alignments via the simultaneous fit of a cosmological model and an intrinsic alignment model. CFHTLenS spans 154 square degrees in five optical bands, with accurate shear and photometric redshifts for a galaxy sample with a median redshift of zm =0.70. We estimate the 21 sets of cosmic shear correlation functions associated with six redshift bins, each spanning the angular range of 1.5<theta<35 arcmin. We combine this CFHTLenS data with auxiliary cosmological probes: the cosmic microwave background with data from WMAP7, baryon acoustic oscillations with data from BOSS, and a prior on the Hubble constant from the HST distance ladder. This leads to constraints on the normalisation of the matter power spectrum sigma_8 = 0.799 +/- 0.015 and the matter density parameter Omega_m = 0.271 +/- 0.010 for a flat Lambda CDM cosmology. For a flat wCDM cosmology we constrain the dark energy equation of state parameter w = -1.02 +/- 0.09. We also provide constraints for curved Lambda CDM and wCDM cosmologies. We find the intrinsic alignment contamination to be galaxy-type dependent with a significant intrinsic alignment signal found for early-type galaxies, in contrast to the late-type galaxy sample for which the intrinsic alignment signal is found to be consistent with zero.

preprint2013arXiv

CFHTLenS: Combined probe cosmological model comparison using 2D weak gravitational lensing

We present cosmological constraints from 2D weak gravitational lensing by the large-scale structure in the Canada-France Hawaii Telescope Lensing Survey (CFHTLenS) which spans 154 square degrees in five optical bands. Using accurate photometric redshifts and measured shapes for 4.2 million galaxies between redshifts of 0.2 and 1.3, we compute the 2D cosmic shear correlation function over angular scales ranging between 0.8 and 350 arcmin. Using non-linear models of the dark-matter power spectrum, we constrain cosmological parameters by exploring the parameter space with Population Monte Carlo sampling. The best constraints from lensing alone are obtained for the small-scale density-fluctuations amplitude sigma_8 scaled with the total matter density Omega_m. For a flat LambdaCDM model we obtain sigma_8(Omega_m/0.27)^0.6 = 0.79+-0.03. We combine the CFHTLenS data with WMAP7, BOSS and an HST distance-ladder prior on the Hubble constant to get joint constraints. For a flat LambdaCDM model, we find Omega_m = 0.283+-0.010 and sigma_8 = 0.813+-0.014. In the case of a curved wCDM universe, we obtain Omega_m = 0.27+-0.03, sigma_8 = 0.83+-0.04, w_0 = -1.10+-0.15 and Omega_K = 0.006+0.006-0.004. We calculate the Bayesian evidence to compare flat and curved LambdaCDM and dark-energy CDM models. From the combination of all four probes, we find models with curvature to be at moderately disfavoured with respect to the flat case. A simple dark-energy model is indistinguishable from LambdaCDM. Our results therefore do not necessitate any deviations from the standard cosmological model.

preprint2013arXiv

CFHTLenS: Mapping the Large Scale Structure with Gravitational Lensing

We present a quantitative analysis of the largest contiguous maps of projected mass density obtained from gravitational lensing shear. We use data from the 154 deg2 covered by the Canada-France-Hawaii Telescope Lensing Survey. Our study is the first attempt to quantitatively characterize the scientific value of lensing maps, which could serve in the future as a complementary approach to the study of the dark universe with gravitational lensing. We show that mass maps contain unique cosmological information beyond that of traditional two-points statistical analysis techniques. Using a series of numerical simulations, we first show how, reproducing the CFHTLenS observing conditions, gravitational lensing inversion provides a reliable estimate of the projected matter distribution of large scale structure. We validate our analysis by quantifying the robustness of the maps with various statistical estimators. We then apply the same process to the CFHTLenS data. We find that the 2-points correlation function of the projected mass is consistent with the cosmological analysis performed on the shear correlation function discussed in the CFHTLenS companion papers. The maps also lead to a significant measurement of the third order moment of the projected mass, which is in agreement with analytic predictions, and to a marginal detection of the fourth order moment. Tests for residual systematics are found to be consistent with zero for the statistical estimators we used. A new approach for the comparison of the reconstructed mass map to that predicted from the galaxy distribution reveals the existence of giant voids in the dark matter maps as large as 3 degrees on the sky. Our analysis shows that lensing mass maps can be used for new techniques such as peak statistics and the morphological analysis of the projected dark matter distribution.

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

CFHTLenS: The relation between galaxy dark matter haloes and baryons from weak gravitational lensing

We present a study of the relation between dark matter halo mass and the baryonic content of host galaxies, quantified via luminosity and stellar mass. Our investigation uses 154 deg2 of Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) lensing and photometric data, obtained from the CFHT Legacy Survey. We employ a galaxy-galaxy lensing halo model which allows us to constrain the halo mass and the satellite fraction. Our analysis is limited to lenses at redshifts between 0.2 and 0.4. We express the relationship between halo mass and baryonic observable as a power law. For the luminosity-halo mass relation we find a slope of 1.32+/-0.06 and a normalisation of 1.19+0.06-0.07x10^13 h70^-1 Msun for red galaxies, while for blue galaxies the best-fit slope is 1.09+0.20-0.13 and the normalisation is 0.18+0.04-0.05x10^13 h70^-1 Msun. Similarly, we find a best-fit slope of 1.36+0.06-0.07 and a normalisation of 1.43+0.11-0.08x10^13 h70^-1 Msun for the stellar mass-halo mass relation of red galaxies, while for blue galaxies the corresponding values are 0.98+0.08-0.07 and 0.84+0.20-0.16x10^13 h70^-1 Msun. For red lenses, the fraction which are satellites tends to decrease with luminosity and stellar mass, with the sample being nearly all satellites for a stellar mass of 2x10^9 h70^-2 Msun. The satellite fractions are generally close to zero for blue lenses, irrespective of luminosity or stellar mass. This, together with the shallower relation between halo mass and baryonic tracer, is a direct confirmation from galaxy-galaxy lensing that blue galaxies reside in less clustered environments than red galaxies. We also find that the halo model, while matching the lensing signal around red lenses well, is prone to over-predicting the large-scale signal for faint and less massive blue lenses. This could be a further indication that these galaxies tend to be more isolated than assumed. [abridged]

preprint2013arXiv

Defining a weak lensing experiment in space

This paper describes the definition of a typical next-generation space-based weak gravitational lensing experiment. We first adopt a set of top-level science requirements from the literature, based on the scale and depth of the galaxy sample, and the avoidance of systematic effects in the measurements which would bias the derived shear values. We then identify and categorise the contributing factors to the systematic effects, combining them with the correct weighting, in such a way as to fit within the top-level requirements. We present techniques which permit the performance to be evaluated and explore the limits at which the contributing factors can be managed. Besides the modelling biases resulting from the use of weighted moments, the main contributing factors are the reconstruction of the instrument point spread function (PSF), which is derived from the stellar images on the image, and the correction of the charge transfer inefficiency (CTI) in the CCD detectors caused by radiation damage.

preprint2013arXiv

Describing galaxy weak lensing measurements from tenths to tens of Mpc and up to z~0.6 with a single model

The clustering of galaxies and the matter distribution around them can be described using the halo model complemented with a realistic description of the way galaxies populate dark matter haloes. This has been used successfully to describe statistical properties of samples of galaxies at z<0.2. Without adjusting any model parameters, we compare the predicted weak lensing signal induced by Luminous Red Galaxies to measurements from SDSS DR7 on much larger scales (up to ~90 h_{70}^{-1} Mpc) and at higher redshift (z~0.4). We find excellent agreement, suggesting that the model captures the main properties of the galaxy-dark matter connection. To extend the comparison to lenses at even higher redshifts we complement the SDSS data with shape measurements from the deeper RCS2, resulting in precise lensing measurements for lenses up to z~0.6. These measurements are also well described using the same model. Considering solely these weak lensing measurements, we robustly assess that, up to z~0.6, the number of central galaxies as a function of halo mass is well described by a log-normal distribution with scatter $σ_{\log L_{\rm c}}=0.146\pm0.011$, in agreement with previous independent studies at lower redshift. Our results demonstrate the value of complementing the information about the properties of the (lens) galaxies provided by SDSS with deeper, high-quality imaging data.

preprint2013arXiv

Discovery of a Rich Cluster at z = 1.63 using the Rest-Frame 1.6um "Stellar Bump Sequence" Method

We present a new two-color algorithm, the "Stellar Bump Sequence" (SBS), that is optimized for robustly identifying candidate high-redshift galaxy clusters in combined wide-field optical and mid-infrared (MIR) data. The SBS algorithm is a fusion of the well-tested cluster red-sequence method of Gladders & Yee (2000) with the MIR 3.6um - 4.5um cluster detection method developed by Papovich (2008). As with the cluster red-sequence method, the SBS identifies candidate overdensities within 3.6um - 4.5um color slices, which are the equivalent of a rest-frame 1.6um stellar bump "red-sequence". In addition to employing the MIR colors of galaxies, the SBS algorithm incorporates an optical/MIR (z' - 3.6um) color cut. This cut effectively eliminates foreground 0.2 < z < 0.4 galaxies which have 3.6um - 4.5um colors that are similarly red as z > 1.0 galaxies and add noise when searching for high-redshift galaxy overdensities. We demonstrate using the z ~ 1 GCLASS cluster sample that similar to the red sequence, the stellar bump sequence appears to be a ubiquitous feature of high-redshift clusters, and that within that sample the color of the stellar bump sequence increases monotonically with redshift and provides photometric redshifts accurate to dz = 0.05. We apply the SBS method in the XMM-LSS SWIRE field and show that it robustly recovers the majority of confirmed optical, MIR, and X-ray-selected clusters at z > 1.0 in that field. Lastly, we present confirmation of SpARCS J022427-032354 at z = 1.63, a new cluster detected with the method and confirmed with 12 high-confidence spectroscopic redshifts obtained using FORS2 on the VLT. We conclude with a discussion of future prospects for using the algorithm.

preprint2013arXiv

Effect of baryonic feedback on two- and three-point shear statistics: prospects for detection and improved modelling

Accurate knowledge of the effect of feedback from galaxy formation on the matter distribution is a key requirement for future weak lensing experiments. Recent studies using hydrodynamic simulations have shown that different baryonic feedback scenarios lead to significantly different two-point shear statistics. In this paper we extend earlier work to three-point shear statistics. We show that, relative to the predictions of dark matter only models, the amplitude of the signal can be reduced by as much as 30-40% on scales of a few arcminutes. We find that baryonic feedback may affect two- and three-point shear statistics differently and demonstrate that this can be used to assess the fidelity of various feedback models. In particular, upcoming surveys such as Euclid might be able to discriminate between different feedback models by measuring both second- and third-order statistics. Because it will likely remain impossible to predict baryonic feedback with high accuracy from first principles, we argue in favour of phenomenological models that can capture the relevant effects of baryonic feedback processes in addition to changes in cosmology. We construct such a model by modifying the dark matter-only halo model to characterise the generic effects of energetic feedback using a small number of parameters. We use this model to perform a likelihood analysis in a simplified case in which two- and three-point shear statistics are measured between 0.5 and 20 arcmin and in which the amplitude of fluctuations, sigma8, the matter density parameter, Om, and the dark energy parameter, w0, are the only unknown free parameters. We demonstrate that for weak lensing surveys such as Euclid, marginalising over the feedbac parameters describing the effects of baryonic processes, such as outflows driven by feedback from star formation and AGN, may be able to mitigate the bias affecting Om, sigma8 and w0.

preprint2013arXiv

Masses of galaxy clusters from gravitational lensing

Despite consistent progress in numerical simulations, the observable properties of galaxy clusters are difficult to predict ab initio. It is therefore important to compare both theoretical and observational results to a direct measure of the cluster mass. This can be done by measuring the gravitational lensing effects caused by the bending of light by the cluster mass distribution. In this review we discuss how this phenomenon can be used to determine cluster masses and study the mass distribution itself. As sample sizes increase, the accuracy of the weak lensing mass estimates needs to improve accordingly. We discuss the main practical aspects of these measurements. We review a number of applications and highlight some recent results.

preprint2013arXiv

Search for Rayleigh scattering in the atmosphere of GJ1214b

We investigate the atmosphere of GJ1214b, a transiting super-Earth planet with a low mean density, by measuring its transit depth as a function of wavelength in the blue optical portion of the spectrum. It is thought that this planet is either a mini-Neptune, consisting of a rocky core with a thick, hydrogen-rich atmosphere, or a planet with a composition dominated by water. Most observations favor a water-dominated atmosphere with a small scale-height, however, some observations indicate that GJ1214b could have an extended atmosphere with a cloud layer muting the molecular features. In an atmosphere with a large scale-height, Rayleigh scattering at blue wavelengths is likely to cause a measurable increase in the apparent size of the planet towards the blue. We observed the transit of GJ1214b in the B-band with the FOcal Reducing Spectrograph (FORS) at the Very Large Telescope (VLT) and in the g-band with both ACAM on the William Hershel Telescope (WHT) and the Wide Field Camera (WFC) at the Isaac Newton Telescope (INT). We find a planet-to-star radius ratio in the B-band of 0.1162+/-0.0017, and in the g-band 0.1180+/-0.0009 and 0.1174+/-0.0017 for the WHT & INT observations respectively. These optical data do not show significant deviations from previous measurements at longer wavelengths. In fact, a flat transmission spectrum across all wavelengths best describes the combined observations. When atmospheric models are considered a small scale-height water-dominated model fits the data best.

preprint2013arXiv

The environmental dependence of the stellar mass function at z~1: Comparing cluster and field between the GCLASS and UltraVISTA surveys

We present the stellar mass functions (SMFs) of star-forming and quiescent galaxies from observations of 10 rich clusters in the Gemini Cluster Astrophysics Spectroscopic Survey (GCLASS) in the redshift range 0.86<z<1.34. We compare our results with field measurements at similar redshifts using data from a Ks-band selected catalogue of the COSMOS/UltraVISTA field. We construct a Ks-band selected multi-colour catalogue for the clusters in 11 photometric bands covering u-8um, and estimate photometric redshifts and stellar masses using SED fitting techniques. To correct for interlopers in our cluster sample, we use the deep spectroscopic component of GCLASS, which contains spectra for 1282 identified cluster and field galaxies taken with Gemini/GMOS. Both the photometric and spectroscopic samples are sufficiently deep that we can probe the SMF down to masses of 10^10 Msun. We distinguish between star-forming and quiescent galaxies using the rest-frame U-V versus V-J diagram, and find that the best-fitting Schechter parameters alpha and M* are similar within the uncertainties for these galaxy types within the different environments. However, there is a significant difference in the shape and normalisation of the total SMF between the clusters and the field sample. This difference in the total SMF is primarily a reflection of the increased fraction of quiescent galaxies in high-density environments. We apply a simple quenching model that includes components of mass- and environment-driven quenching, and find that in this picture 45% of the star-forming galaxies, which normally would be forming stars in the field, are quenched by the cluster. If galaxies in clusters and the field quench their star formation via different mechanisms, these processes have to conspire in such a way that the shapes of the quiescent and star-forming SMF remain similar in these different environments.

preprint2013arXiv

Weak Lensing Calibrated M-T Scaling Relation of Galaxy Groups in the COSMOS Field

The scaling between X-ray observables and mass for galaxy clusters and groups is instrumental for cluster based cosmology and an important probe for the thermodynamics of the intracluster gas. We calibrate a scaling relation between the weak lensing mass and X-ray spectroscopic temperature for 10 galaxy groups in the COSMOS field, combined with 55 higher mass clusters from the literature. The COSMOS data includes HST imaging and redshift measurements of 46 source galaxies per square arcmin, enabling us to perform unique weak lensing measurements of low mass systems. Our sample extends the mass range of the lensing calibrated M-T relation an order of magnitude lower than any previous study, resulting in a power-law slope of 1.48$^{+0.13}_{-0.09}$. The slope is consistent with the self-similar model, predictions from simulations, and observations of clusters. However, X-ray observations relying on mass measurements derived under the assumption of hydrostatic equilibrium have indicated that masses at group scales are lower than expected. Both simulations and observations suggest that hydrostatic mass measurements can be biased low. Our external weak lensing masses provides the first observational support for hydrostatic mass bias at group level, showing an increasing bias with decreasing temperature and reaching a level of 30-50% at 1 keV.

preprint2012arXiv

A Giant Metrewave Radio Telescope/Chandra view of IRAS 09104+4109: A type 2 QSO in a cooling flow

IRAS 09104+4109 is a rare example of a dust enshrouded type 2 QSO in the centre of a cool-core galaxy cluster. Previous observations of this z=0.44 system showed that as well as powering the hyper-luminous infrared emission of the cluster-central galaxy, the QSO is associated with a double-lobed radio source. However, the steep radio spectral index and misalignment between the jets and ionised optical emission suggested that the orientation of the QSO had recently changed. We use a combination of new, multi-band Giant Metrewave Radio Telescope observations and archival radio data to confirm that the jets are no longer powered by the QSO, and estimate their age to be 120-160 Myr. This is in agreement with the ~70-200 Myr age previously estimated for star-formation in the galaxy. Previously unpublished Very Long Baseline Array data reveal a 200 pc scale double radio source in the galaxy core which is more closely aligned with the current QSO axis and may represent a more recent period of jet activity. These results suggest that the realignment of the QSO, the cessation of jet activity, and the onset of rapid star-formation may have been caused by a gas-rich galaxy merger. A Chandra X-ray observation confirms the presence of cavities associated with the radio jets, and we estimate the energy required to inflate them to be ~7.7x10^60 erg. The mechanical power of the jets is sufficient to balance radiative cooling in the cluster, provided they are efficiently coupled to the intra-cluster medium (ICM). We find no evidence of direct radiative heating and conclude that the QSO either lacks the radiative luminosity to heat the ICM, or that it requires longer than 100-200 Myr to significantly impact its environment. [Abridged]

preprint2012arXiv

CFHTLenS tomographic weak lensing: Quantifying accurate redshift distributions

The Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) comprises deep multi-colour (u*g'r'i'z') photometry spanning 154 square degrees, with accurate photometric redshifts and shape measurements. We demonstrate that the redshift probability distribution function summed over galaxies provides an accurate representation of the galaxy redshift distribution accounting for random and catastrophic errors for galaxies with best fitting photometric redshifts z_p < 1.3. We present cosmological constraints using tomographic weak gravitational lensing by large-scale structure. We use two broad redshift bins 0.5 < z_p <= 0.85 and 0.85 < z_p <= 1.3 free of intrinsic alignment contamination, and measure the shear correlation function on angular scales in the range ~1-40 arcmin. We show that the problematic redshift scaling of the shear signal, found in previous CFHTLS data analyses, does not afflict the CFHTLenS data. For a flat Lambda-CDM model and a fixed matter density Omega_m=0.27, we find the normalisation of the matter power spectrum sigma_8=0.771 \pm 0.041. When combined with cosmic microwave background data (WMAP7), baryon acoustic oscillation data (BOSS), and a prior on the Hubble constant from the HST distance ladder, we find that CFHTLenS improves the precision of the fully marginalised parameter estimates by an average factor of 1.5-2. Combining our results with the above cosmological probes, we find Omega_m=0.2762 \pm 0.0074 and sigma_8=0.802 \pm 0.013.

preprint2012arXiv

CFHTLenS: The Canada-France-Hawaii Telescope Lensing Survey

We present the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) that accurately determines a weak gravitational lensing signal from the full 154 square degrees of deep multi-colour data obtained by the CFHT Legacy Survey. Weak gravitational lensing by large-scale structure is widely recognised as one of the most powerful but technically challenging probes of cosmology. We outline the CFHTLenS analysis pipeline, describing how and why every step of the chain from the raw pixel data to the lensing shear and photometric redshift measurement has been revised and improved compared to previous analyses of a subset of the same data. We present a novel method to identify data which contributes a non-negligible contamination to our sample and quantify the required level of calibration for the survey. Through a series of cosmology-insensitive tests we demonstrate the robustness of the resulting cosmic shear signal, presenting a science-ready shear and photometric redshift catalogue for future exploitation.

preprint2012arXiv

Constraints on the shapes of galaxy dark matter haloes from weak gravitational lensing

We study the shapes of galaxy dark matter haloes by measuring the anisotropy of the weak gravitational lensing signal around galaxies in the second Red-sequence Cluster Survey (RCS2). We determine the average shear anisotropy within the virial radius for three lens samples: all galaxies with 19<m_r'<21.5, and the `red' and `blue' samples, whose lensing signals are dominated by massive low-redshift early-type and late-type galaxies, respectively. To study the environmental dependence of the lensing signal, we separate each lens sample into an isolated and clustered part and analyse them separately. We also measure the azimuthal dependence of the distribution of physically associated galaxies around the lens samples. We find that these satellites preferentially reside near the major axis of the lenses, and constrain the angle between the major axis of the lens and the average location of the satellites to <theta>=43.7 deg +/- 0.3 deg for the `all' lenses, <theta>=41.7 deg +/- 0.5 deg for the `red' lenses and <theta>=42.0 deg +/- 1.4 deg for the `blue' lenses. For the `all' sample, we find that the anisotropy of the galaxy-mass cross-correlation function <f-f_45>=0.23 +/- 0.12, providing weak support for the view that the average galaxy is embedded in, and preferentially aligned with, a triaxial dark matter halo. Assuming an elliptical Navarro-Frenk-White (NFW) profile, we find that the ratio of the dark matter halo ellipticity and the galaxy ellipticity f_h=e_h/e_g=1.50+1.03-1.01, which for a mean lens ellipticity of 0.25 corresponds to a projected halo ellipticity of e_h=0.38+0.26-0.25 if the halo and the lens are perfectly aligned. For isolated galaxies of the `all' sample, the average shear anisotropy increases to <f-f_45>=0.51+0.26-0.25 and f_h=4.73+2.17-2.05, whilst for clustered galaxies the signal is consistent with zero. (abridged)

preprint2012arXiv

On Combining Galaxy Clustering and Weak Lensing to Unveil Galaxy Biasing via the Halo Model

We formulate the concept of non-linear and stochastic galaxy biasing in the framework of halo occupation statistics. Using two-point statistics in projection, we define the galaxy bias function, b_g(r_p), and the galaxy-dark matter cross-correlation function, R_{gm}(r_p), where r_p is the projected distance. We use the analytical halo model to predict how the scale dependence of b_g and R_{gm}, over the range 0.1 Mpc/h < r_p < 30 Mpc/h, depends on the non-linearity and stochasticity in halo occupation models. In particular we quantify the effect due to the presence of central galaxies, the assumption for the radial distribution of satellite galaxies, the richness of the halo, and the Poisson character of the probability to have a certain number of satellite galaxies in a halo of a certain mass. Overall, brighter galaxies reveal a stronger scale dependence, and out to a larger radius. In real-space, we find that galaxy bias becomes scale independent, with R_{gm} = 1, for radii r > 1 - 5 Mpc/h, depending on luminosity. However, galaxy bias is scale-dependent out to much larger radii when one uses the projected quantities defined in this paper. These projected bias functions have the advantage that they are more easily accessible observationally and that their scale dependence carries a wealth of information regarding the properties of galaxy biasing. To observationally constrain the parameters of the halo occupation statistics and to unveil the origin of galaxy biasing we propose the use of the bias function Gamma_{gm}(r_p)=b_g(r_p)/R_{gm}(r_p). This function is obtained via a combination of weak gravitational lensing and galaxy clustering, and it can be measured using existing and forthcoming imaging and spectroscopic galaxy surveys.

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.

preprint2012arXiv

Origins of weak lensing systematics, and requirements on future instrumentation (or knowledge of instrumentation)

The first half of this paper explores the origin of systematic biases in the measurement of weak gravitational lensing. Compared to previous work, we expand the investigation of PSF instability and fold in for the first time the effects of non-idealities in electronic imaging detectors and imperfect galaxy shape measurement algorithms. Together, these now explain the additive A(l) and multiplicative M(l) systematics typically reported in current lensing measurements. We find that overall performance is driven by a product of a telescope/camera's *absolute performance*, and our *knowledge about its performance*. The second half of this paper propagates any residual shear measurement biases through to their effect on cosmological parameter constraints. Fully exploiting the statistical power of Stage IV weak lensing surveys will require additive biases A<1.8e-12 and multiplicative biases M<4.0e-3. These can be allocated between individual budgets in hardware, calibration data and software, using results from the first half of the paper. If instrumentation is stable and well-calibrated, we find extant shear measurement software from GREAT10 already meet requirements on galaxies detected at S/N=40. Averaging over a population of galaxies with a realistic distribution of sizes, it also meets requirements for a 2D cosmic shear analysis from space. If used on fainter galaxies or for 3D cosmic shear tomography, existing algorithms would need calibration on simulations to avoid introducing bias at a level similar to the statistical error. Requirements on hardware and calibration data are discussed in more detail in a companion paper. Our analysis is intentionally general, but is specifically being used to drive the hardware and ground segment performance budget for the design of the European Space Agency's recently-selected Euclid mission.

preprint2012arXiv

Quantifying the effect of baryon physics on weak lensing tomography

We use matter power spectra from cosmological hydrodynamic simulations to quantify the effect of baryon physics on the weak gravitational lensing shear signal. The simulations consider a number of processes, such as radiative cooling, star formation, supernovae and feedback from active galactic nuclei (AGN). Van Daalen et al. (2011) used the same simulations to show that baryon physics, in particular the strong feedback that is required to solve the overcooling problem, modifies the matter power spectrum on scales relevant for cosmological weak lensing studies. As a result, the use of power spectra from dark matter simulations can lead to significant biases in the inferred cosmological parameters. We show that the typical biases are much larger than the precision with which future missions aim to constrain the dark energy equation of state, w_0. For instance, the simulation with AGN feedback, which reproduces X-ray and optical properties of groups of galaxies, gives rise to a ~40% bias in w_0. We demonstrate that the modification of the power spectrum is dominated by groups and clusters of galaxies, the effect of which can be modelled. We consider an approach based on the popular halo model and show that simple modifications can capture the main features of baryonic feedback. Despite its simplicity, we find that our model, when calibrated on the simulations, is able to reduce the bias in w_0 to a level comparable to the size of the statistical uncertainties for a Euclid-like mission. While observations of the gas and stellar fractions as a function of halo mass can be used to calibrate the model, hydrodynamic simulations will likely still be needed to extend the observed scaling relations down to halo masses of 10 ^12 M_sun/h.

preprint2012arXiv

Stellar mass versus velocity dispersion as tracer of the lensing signal around bulge-dominated galaxies

We present the results of a weak gravitational lensing analysis to determine whether the stellar mass or the velocity dispersion is more closely related to the amplitude of the lensing signal around galaxies - and hence to the projected distribution of dark matter. The lensing signal on scales smaller than the virial radius corresponds most closely to the lensing velocity dispersion in the case of a singular isothermal profile, but is on larger scales also sensitive to the clustering of the haloes. We select over 4000 lens galaxies at a redshift z<0.2 with concentrated (or bulge-dominated) surface brightness profiles from the ~300 square degree overlap between the Red-sequence Cluster Survey 2 (RCS2) and the data release 7 (DR7) of the Sloan Digital Sky Survey (SDSS). We consider both the spectroscopic velocity dispersion and a model velocity dispersion (a combination of the stellar mass, the size and the Sersic index of a galaxy). Comparing the model and spectroscopic velocity dispersion we find that they correlate well for galaxies with concentrated brightness profiles. We find that the stellar mass and the spectroscopic velocity dispersion trace the amplitude of the lensing signal on small scales equally well. The model velocity dispersion, however, does significantly worse. A possible explanation is that the halo properties that determine the small-scale lensing signal - mainly the total mass - also depend on the structural parameters of galaxies, such as the effective radius and Sersic index, but we lack data for a definitive conclusion.

preprint2012arXiv

The Canadian Cluster Comparison Project: weak lensing masses and SZ scaling relations

The Canadian Cluster Comparison Project is a comprehensive multi-wavelength survey targeting 50 massive X-ray selected clusters of galaxies to examine baryonic tracers of cluster mass and to probe the cluster-to-cluster variation in the thermal properties of the hot intracluster medium. In this paper we present the weak lensing masses, based on the analysis of deep wide-field imaging data obtained using the Canada-France-Hawaii-Telescope. The final sample includes two additional clusters that were located in the field-of-view. We take these masses as our reference for the comparison of cluster properties at other wavelengths. In this paper we limit the comparison to published measurements of the Sunyaev-Zel'dovich effect. We find that this signal correlates well with the projected lensing mass, with an intrinsic scatter of 12\pm5% at ~r_2500, demonstrating it is an excellent proxy for cluster mass.

preprint2012arXiv

The Environmental Dependence of the Incidence of Galactic Tidal Features

In a sample of 54 galaxy clusters (0.04<z<0.15) containing 3551 early-type galaxies suitable for study, we identify those with tidal features both interactively and automatically. We find that ~3% have tidal features that can be detected with data that reaches a 3-sigma sensitivity limit of 26.5 mag arcsec^-2. Regardless of the method used to classify tidal features, or the fidelity imposed on such classifications, we find a deficit of tidally disturbed galaxies with decreasing clustercentric radius that is most pronounced inside of ~0.5R_200. We cannot distinguish whether the trend arises from an increasing likelihood of recent mergers with increasing clustercentric radius or a decrease in the lifetime of tidal features with decreasing clustercentric radius. We find no evidence for a relationship between local density and the incidence of tidal features, but our local density measure has large uncertainties. We find interesting behavior in the rate of tidal features among cluster early-types as a function of clustercentric radius and expect such results to provide constraints on the effect of the cluster environment on the structure of galaxy halos, the build-up of the red sequence of galaxies, and the origin of the intracluster stellar population.

preprint2011arXiv

Galaxy-galaxy lensing constraints on the relation between baryons and dark matter in galaxies in the Red Sequence Cluster Survey 2

We present the results of a study of weak gravitational lensing by galaxies using imaging data that were obtained as part of the second Red Sequence Cluster Survey (RCS2). In order to compare to the baryonic properties of the lenses we focus here on the ~300 square degrees that overlap with the DR7 of the SDSS. The depth and image quality of the RCS2 enables us to significantly improve upon earlier work for luminous galaxies at z>=0.3. Comparison with dynamical masses from the SDSS shows a good correlation with the lensing mass for early-type galaxies. For low luminosity (stellar mass) early-type galaxies we find a satellite fraction of ~40% which rapidly decreases to <10% with increasing luminosity (stellar mass). The satellite fraction of the late-types has a value in the range 0-15%. We find that early-types in the range 10^10<L_r<10^11.5 Lsun have virial masses that are about five times higher than those of late-type galaxies and that the mass scales as M_200 \propto L^2.34 +0.09 \ -0.16. We also measure the virial mass-to-light ratio, and find for L_200<10^11 Lsun a value of M_200/L_200=42+-10 for early-types, which increases for higher luminosities to values that are consistent with those observed for groups and clusters of galaxies. For late-type galaxies we find a lower value of M_200/L_200=17+-9. Our measurements also show that early- and late-type galaxies have comparable halo masses for stellar masses M_*<10^11 Msun, whereas the virial masses of early-type galaxies are higher for higher stellar masses. Finally, we determine the efficiency with which baryons have been converted into stars. Our results for early-type galaxies suggest a variation in efficiency with a minimum of ~10% for a stellar mass M_*,200=10^12 Msun. The results for the late-type galaxies are not well constrained, but do suggest a larger value. (abridged)

preprint2011arXiv

Lensing Magnification: A novel method to weigh high-redshift clusters and its application to SpARCS

We introduce a novel method to measure the masses of galaxy clusters at high redshift selected from optical and IR Spitzer data via the red-sequence technique. Lyman-break galaxies are used as a well understood, high-redshift background sample allowing mass measurements of lenses at unprecedented high redshifts using weak lensing magnification. By stacking a significant number of clusters at different redshifts with average masses of ~1-3x10^14M_sun, as estimated from their richness, we can calibrate the normalisation of the mass-richness relation. With the current data set (area: 6 deg^2) we detect a magnification signal at the >3-sigma level. There is good agreement between the masses estimated from the richness of the clusters and the average masses estimated from magnification, albeit with large uncertainties. We perform tests that suggest the absence of strong systematic effects and support the robustness of the measurement. This method - when applied to larger data sets in the future - will yield an accurate calibration of the mass-observable relations at z>~1 which will represent an invaluable input for cosmological studies using the galaxy cluster mass function and astrophysical studies of cluster formation. Furthermore this method will probably be the least expensive way to measure masses of large numbers of z>1 clusters detected in future IR-imaging surveys.

preprint2011arXiv

The Gemini Cluster Astrophysics Spectroscopic Survey (GCLASS): The Role of Environment and Self-Regulation in Galaxy Evolution at z ~ 1

We evaluate the effects of environment and stellar mass on galaxy properties at 0.85 < z < 1.20 using a 3.6um-selected spectroscopic sample of 797 cluster and field galaxies drawn from the GCLASS survey. We confirm that for galaxies with LogM* > 9.3 the well-known correlations between environment and properties such as star-forming fraction (f_SF), SFR, SSFR, D(4000), and color are already in place at z ~ 1. We separate the effects of environment and stellar mass on galaxies by comparing the properties of star-forming and quiescent galaxies at fixed environment, and fixed stellar mass. The SSFR of star-forming galaxies at fixed environment is correlated with stellar mass; however, at fixed stellar mass it is independent of environment. The same trend exists for the D(4000) measures of both the star-forming and quiescent galaxies and shows that their properties are determined primarily by their stellar mass, not by their environment. Instead, it appears that environment's primary role is to control the fraction of star-forming galaxies. Using the spectra we identify candidate poststarburst galaxies and find that those with 9.3 < LogM* < 10.7 are 3.1 +/- 1.1 times more common in high-density regions compared to low-density regions. The clear association of poststarbursts with high-density regions as well as the lack of a correlation between the SSFRs and D(4000)s of star-forming galaxies with their environment suggests that at z ~ 1 the environmental-quenching timescale must be rapid. Lastly, we construct a simple quenching model which demonstrates that the lack of a correlation between the D(4000) of quiescent galaxies and their environment results naturally if self quenching dominates over environmental quenching at z > 1, or if the evolution of the self-quenching rate mirrors the evolution of the environmental-quenching rate at z > 1, regardless of which dominates.

preprint2011arXiv

The impact of high spatial frequency atmospheric distortions on weak lensing measurements

High precision cosmology with weak gravitational lensing requires a precise measure of the Point Spread Function across the imaging data where the accuracy to which high spatial frequency variation can be modelled is limited by the stellar number density across the field. We analyse dense stellar fields imaged at the Canada-France-Hawaii Telescope to quantify the degree of high spatial frequency variation in ground-based imaging Point Spread Functions and compare our results to models of atmospheric turbulence. The data shows an anisotropic turbulence pattern with an orientation independent of the wind direction and wind speed. We find the amplitude of the high spatial frequencies to decrease with increasing exposure time as $t^{-1/2}$, and find a negligibly small atmospheric contribution to the Point Spread Function ellipticity variation for exposure times $t>180$ seconds. For future surveys analysing shorter exposure data, this anisotropic turbulence will need to be taken into account as the amplitude of the correlated atmospheric distortions becomes comparable to a cosmological lensing signal on scales less than $\sim 10$ arcminutes. This effect could be mitigated, however, by correlating galaxy shear measured on exposures imaged with a time separation greater than 50 seconds, for which we find the spatial turbulence patterns to be uncorrelated.

preprint2011arXiv

The Multi-Epoch Nearby Cluster Survey: type Ia supernova rate measurement in z~0.1 clusters and the late-time delay time distribution

We describe the Multi-Epoch Nearby Cluster Survey (MENeaCS), designed to measure the cluster Type Ia supernova (SN Ia) rate in a sample of 57 X-ray selected galaxy clusters, with redshifts of 0.05 < z < 0.15. Utilizing our real time analysis pipeline, we spectroscopically confirmed twenty-three cluster SN Ia, four of which were intracluster events. Using our deep CFHT/Megacam imaging, we measured total stellar luminosities in each of our galaxy clusters, and we performed detailed supernova detection efficiency simulations. Bringing these ingredients together, we measure an overall cluster SN Ia rate within R_{200} (1 Mpc) of 0.042^{+0.012}_{-0.010}^{+0.010}_{-0.008} SNuM (0.049^{+0.016}_{-0.014}^{+0.005}_{-0.004} SNuM) and a SN Ia rate within red sequence galaxies of 0.041^{+0.015}_{-0.015}^{+0.005}_{-0.010} SNuM (0.041^{+0.019}_{-0.015}^{+0.005}_{-0.004} SNuM). The red sequence SN Ia rate is consistent with published rates in early type/elliptical galaxies in the `field'. Using our red sequence SN Ia rate, and other cluster SNe measurements in early type galaxies up to $z\sim1$, we derive the late time (>2 Gyr) delay time distribution (DTD) of SN Ia assuming a cluster early type galaxy star formation epoch of z_f=3. Assuming a power law form for the DTD, Ψ(t)\propto t^s, we find s=-1.62\pm0.54. This result is consistent with predictions for the double degenerate SN Ia progenitor scenario (s\sim-1), and is also in line with recent calculations for the double detonation explosion mechanism (s\sim-2). The most recent calculations of the single degenerate scenario delay time distribution predicts an order of magnitude drop off in SN Ia rate \sim 6-7 Gyr after stellar formation, and the observed cluster rates cannot rule this out.

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

Spectroscopic Confirmation of Three Red-Sequence Selected Galaxy Clusters at z=0.87, 1.16 and 1.21 from the SpARCS Survey

The Spitzer Adaptation of the Red-sequence Cluster Survey (SpARCS) is a z'-passband imaging survey of the 50 deg^2 Spitzer SWIRE Legacy fields, designed with the primary aim of creating the first large, homogeneously selected sample of massive clusters at z>1. SpARCS uses an infrared adaptation of the two-filter cluster red-sequence technique. In this paper we report Keck/LRIS spectroscopic confirmation of two new exceptionally rich galaxy clusters, SpARCS J161315+564930 at z=0.871+/-0.002, with 14 high-confidence members and a rest-frame velocity dispersion of sigma_v= 1230+/-320 km s^-1, and SpARCS J161641+554513 at z=1.161+/-0.003, with seven high-confidence members (including one AGN) and a rest-frame velocity dispersion of sigma_v=950+/-330 km s^-1. We also report confirmation of a third new system, SpARCS J161037+552417 at z=1.210+/-0.002, with seven high-confidence members and a rest-frame velocity dispersion of sigma_v=410+/-300 km s^-1. These three new spectroscopically confirmed clusters further demonstrate the efficiency and effectiveness of two-filter imaging for detecting bona fide galaxy clusters at high redshift. We conclude by demonstrating that prospects are good for the current generation of surveys aiming to estimate cluster redshifts and masses at z>~1 directly from optical-infrared imaging.

preprint2010arXiv

Structure Formation by Fifth Force: Power Spectrum from N-Body Simulations

We lay out the framework to numerically study nonlinear structure formation in the context of scalar-field-coupled cold dark matter models (phiCDM models) where the scalar field phi serves as dynamical dark energy. Adopting parameters for the scalar field which yield a realistic CMB spectrum, we generate the initial conditions for our Nbody simulations, which follow the spatial distributions of dark matter and the scalar field, by solving their equations of motion using the multilevel adaptive grid technique. We show that the spatial configuration of the scalar field tracks both the voids and clusters of dark matter. The phiCDM model differs from standard LCDM at small scales with observable modifications of, e.g., the mass function of halos as well as the matter power spectrum. Nevertheless, the predictions of both models for the Hubble expansion and the CMB spectrum are virtually indistinguishable. Hence, galaxy cluster counts and weak lensing observations, which probe structure formation at small scales, are needed to falsify this class of models.

preprint2010arXiv

Substructure lensing in galaxy clusters as a constraint on low-mass sterile neutrinos in tensor-vector-scalar theory: The straight arc of Abell 2390

Certain covariant theories of the modified Newtonian dynamics paradigm seem to require an additional hot dark matter (HDM) component - in the form of either heavy ordinary neutrinos or more recently light sterile neutrinos (SNs) with a mass around 11eV - to be relieved of problems ranging from cosmological scales down to intermediate ones relevant for galaxy clusters. Here we suggest using gravitational lensing by galaxy clusters to test such a marriage of neutrino HDM and modified gravity, adopting the framework of tensor-vector-scalar theory (TeVeS). Unlike conventional cold dark matter (CDM), such HDM is subject to strong phase-space constraints, which allows one to check cluster lens models inferred within the modified framework for consistency. Since the considered HDM particles cannot collapse into arbitrarily dense clumps and only form structures well above the galactic scale, systems which indicate the need for dark substructure are of particular interest. As a first example, we study the cluster lens Abell 2390 and its impressive straight arc with the help of numerical simulations. Based on our results, we outline a general and systematic approach to model cluster lenses in TeVeS which significantly reduces the calculation complexity. We further consider a simple bimodal lens configuration, capable of producing the straight arc, to demonstrate our approach. We find that such a model is marginally consistent with the hypothesis of 11eV SNs. Future work including more detailed and realistic lens models may further constrain the necessary SN distribution and help to conclusively assess this point. Cluster lenses could therefore provide an interesting discriminator between CDM and such modified gravity scenarios supplemented by SNs or other choices of HDM.

preprint2010arXiv

The mass-L_x relation for moderate luminosity X-ray clusters

We present measurements of the masses of a sample of 25 moderate X-ray luminosity clusters of galaxies from the 160 square degree ROSAT survey. The masses were obtained from a weak lensing analysis of deep F814W images obtained using the Advanced Camera for Surveys (ACS). We present an accurate empirical correction for the effect of charge transfer (in)efficiency on the shapes of faint galaxies. A significant lensing signal is detected around most of the clusters. The lensing mass correlates tightly with the cluster richness. We measured the intrinsic scatter in the scaling relation between M_2500 and L_X and find the best fit power law slope and normalisation to be alpha=0.68+-0.07 and M_X=(1.2+-0.12)10^14M_sun (for L_X=2x10^44 erg/s). These results agree well with a number of recent studies, but the normalisation is lower compared to the study of Rykoff et al. (2008b). One explanation for this difference may be the fact that (sub)structures projected along the line-of-sight boost both the galaxy counts and the lensing mass. Such superpositions lead to an increased mass at a given L_X when clusters are binned by richness.

preprint2009arXiv

Large-scale Structure and Dynamics of the Most X-ray Luminous Galaxy Cluster Known -- RX J1347-1145

We present photometric, spectroscopic and weak lensing analysis of the large-scale structure and dynamics of the most X-ray luminous galaxy cluster known, RX J1347-1145, at z=0.451. We spectroscopically confirmed 47 new members with LDSS3 on the Magellan telescope. Together with previously known members, we measure a new velocity dispersion of 1163(+/-97) km/s. The mass inferred from our velocity dispersion is M200 = 1.16^{+0.32}_{-0.27}x10^{15} solar mass, with r200=1.85Mpc, under the assumption of a singular isothermal sphere. We also present a weak lensing analysis using deep CFHT data on this cluster, and find a deprojected mass of 1.47^{+0.46}_{-0.43}x10^{15} solar mass within r200, in excellent agreement with our dynamical estimate. Thus, our new dynamical mass estimate is consistent with that from weak lensing and X-ray studies in the literature, resolving a previously claimed discrepancy. We photometrically detect and spectroscopically confirm another massive cluster with sigma=780(+/-100) km/s and M200=3.4^{+1.4}_{-1.1}x10^{14} solar mass ~7Mpc south-west of RX J1347-1145, which we refer to as RXJ1347-SW. Our spectroscopic survey reveals a possible excess of galaxies in velocity space in the region between RX J1347-1145 and RXJ1347-SW; comparing with simulations, this excess appears consistent with that expected from a large filamentary structure traced by galaxies connecting these two clusters.

preprint2009arXiv

Sersiclets - A Matched Filter Extension of Shapelets for Weak Lensing Studies

The precision study of dark matter using weak lensing by large scale structure is strongly constrained by the accuracy with which one can measure galaxy shapes. Several methods have been devised but none have demonstrated the ability to reach the level of precision required by future weak lensing surveys. In this Letter we explore new avenues to the existing Shapelets approach, combining a priori knowledge of the galaxy profile with the power of orthogonal basis function decomposition. This Letter discusses the new issues raised by this matched filter approach and proposes promising alternatives to shape measurement techniques. In particular it appears that the use of a matched filter (e.g. Sersic profile) restricted to elliptical radial fitting functions resolves several well known Shapelet issues.

preprint2008arXiv

Spectroscopic Confirmation of a Massive Red-Sequence-Selected Galaxy Cluster at z = 1.34 in the SpARCS-South Cluster Survey

The Spitzer Adaptation of the Red-sequence Cluster Survey (SpARCS) is a z'-passband imaging survey, consisting of deep (z' ~ 24 AB) observations made from both hemispheres using the CFHT 3.6m and CTIO 4m telescopes. The survey was designed with the primary aim of detecting galaxy clusters at z >~ 1. In tandem with pre-existing 3.6um observations from the Spitzer Space Telescope SWIRE Legacy Survey, SpARCS detects clusters using an infrared adaptation of the two-filter red-sequence cluster technique. The total effective area of the SpARCS cluster survey is 41.9 deg^2. In this paper, we provide an overview of the 13.6 deg^2 Southern CTIO/MOSAICII observations. The 28.3 deg^2 Northern CFHT/MegaCam observations are summarized in a companion paper by Muzzin et al. (2008). In this paper, we also report spectroscopic confirmation of SpARCS J003550-431224, a very rich galaxy cluster at z = 1.335, discovered in the ELAIS-S1 field. To date, this is the highest spectroscopically confirmed redshift for a galaxy cluster discovered using the red-sequence technique. Based on nine confirmed members, SpARCS J003550-431224 has a preliminary velocity dispersion of 1050 +/- 230 km/s. With its proven capability for efficient cluster detection, SpARCS is a demonstration that we have entered an era of large, homogeneously-selected z > 1 cluster surveys.

preprint2008arXiv

Spectroscopic Confirmation of Two Massive Red-Sequence-Selected Galaxy Clusters at z ~ 1.2 in the SpARCS-North Cluster Survey

The Spitzer Adaptation of the Red-sequence Cluster Survey (SpARCS) is a deep z'-band imaging survey covering the Spitzer SWIRE Legacy fields designed to create the first large homogeneously-selected sample of massive clusters at z > 1 using an infrared adaptation of the cluster red-sequence method. We present an overview of the northern component of the survey which has been observed with CFHT/MegaCam and covers 28.3 deg^2. The southern component of the survey was observed with CTIO/MOSAICII, covers 13.6 deg^2, and is summarized in a companion paper by Wilson et al. (2008). We also present spectroscopic confirmation of two rich cluster candidates at z ~ 1.2. Based on Nod-and-Shuffle spectroscopy from GMOS-N on Gemini there are 17 and 28 confirmed cluster members in SpARCS J163435+402151 and SpARCS J163852+403843 which have spectroscopic redshifts of 1.1798 and 1.1963, respectively. The clusters have velocity dispersions of 490 +/- 140 km/s and 650 +/- 160 km/s, respectively which imply masses (M200) of (1.0 +/- 0.9) x 10^{14} M_{solar} and (2.4 +/- 1.8) x 10^{14} M_{solar}. Confirmation of these candidates as bona fide massive clusters demonstrates that two-filter imaging is an effective, yet observationally efficient, method for selecting clusters at z > 1.

preprint2008arXiv

The X-ray Cluster Normalization of the Matter Power Spectrum

The number density of galaxy clusters provides tight statistical constraints on the matter fluctuation power spectrum normalization, traditionally phrased in terms of sigma_8, the root mean square mass fluctuation in spheres with radius 8 h^-1 Mpc. We present constraints on sigma_8 and the total matter density Omega_m0 from local cluster counts as a function of X-ray temperature, taking care to incorporate and minimize systematic errors that plagued previous work with this method. In particular, we present new determinations of the cluster luminosity - temperature and mass - temperature relations, including their intrinsic scatter, and a determination of the Jenkins mass function parameters for the same mass definition as the mass - temperature calibration. Marginalizing over the 12 uninteresting parameters associated with this method, we find that the local cluster temperature function implies sigma_8 (Omega_m0/0.32)^alpha = 0.86+/-0.04 with alpha = 0.30 (0.41) for Omega_m0 < 0.32 (Omega_mo > 0.32) (68% confidence for two parameters). This result agrees with a wide range of recent independent determinations, and we find no evidence of any additional sources of systematic error for the X-ray cluster temperature function determination of the matter power spectrum normalization. The joint WMAP5 + cluster constraints are: Omega_m0 = 0.30+0.03/-0.02 and sigma_8 = 0.85+0.04/-0.02 (68% confidence for two parameters).

preprint2008arXiv

Weak Gravitational Lensing and its Cosmological Applications

Weak gravitational lensing is a unique probe of the dark side of the universe: it provides a direct way to map the distribution of dark matter around galaxies, clusters of galaxies and on cosmological scales. Furthermore, the measurement of lensing induced distortions of the shapes of distant galaxies is a powerful probe of dark energy. In this review we describe how lensing measurements are made and interpreted. We discuss various systematic effects that can hamper progress and how they may be overcome. We review some of the recent results in weak lensing by galaxies, galaxy clusters and cosmic shear and discuss the prospects for dark energy measurements from planned surveys.