Lensing Surveys Research Articles

ALMA 1.3 mm Survey of Lensed Submillimeter Galaxies Selected by Herschel: Discovery of Spatially Extended SMGs and Implications

Abstract We present an ALMA 1.3 mm (Band 6) continuum survey of lensed submillimeter galaxies (SMGs) at z = 1.0 to ∼3.2 with an angular resolution of ∼0.″2. These galaxies were uncovered by the Herschel Lensing Survey and feature exceptionally bright far-infrared continuum emission (S peak ≳ 90 mJy) owing to their lensing magnification. We detect 29 sources in 20 fields of massive galaxy clusters with ALMA. Using both the Spitzer/IRAC (3.6/4.5 μm) and ALMA data, we have successfully modeled the surface brightness profiles of 26 sources in the rest-frame near- and far-infrared. Similar to previous studies, we find the median dust-to-stellar continuum size ratio to be small (R e,dust/R e,star = 0.38 ± 0.14) for the observed SMGs, indicating that star formation is centrally concentrated. This is, however, not the case for two spatially extended main-sequence SMGs with a low surface brightness at 1.3 mm (≲0.1 mJy arcsec−2), in which the star formation is distributed over the entire galaxy (R e,dust/R e,star > 1). As a whole, our SMG sample shows a tight anticorrelation between (R e,dust/R e,star) and far-infrared surface brightness (ΣIR) over a factor of ≃1000 in ΣIR. This indicates that SMGs with less vigorous star formation (i.e., lower ΣIR) lack central starburst and are likely to retain a broader spatial distribution of star formation over the whole galaxies (i.e., larger R e,dust/R e,star). The same trend can be reproduced with cosmological simulations as a result of central starburst and potentially subsequent “inside-out” quenching, which likely accounts for the emergence of compact quiescent galaxies at z ∼ 2.

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 weak lensing surveys (NGVSLenS, KiDS/KV450, CFHTLenS, CS82, and RCSLenS, listed in order of most to least constraining) in order to tighten observational constraints on galaxy-scale halo ellipticity for photometrically selected lens samples. We constrain fh, the average ratio between the aligned component of the halo ellipticity and the ellipticity of the light distribution, finding fh = 0.303−0.079+0.080 for red lens galaxies and fh = 0.217−0.159+0.160 for blue lens galaxies when assuming elliptical Navarro-Frenk-White 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 ⟨|ϵh|⟩ = 0.174 ± 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, which constitute an important source of systematic uncertainty in cosmological weak lensing studies.

The PAU Survey: an improved photo-zsample in the COSMOS field

ABSTRACTWe present – and make publicly available – accurate and precise photometric redshifts in the ACS footprint from the COSMOS field for objects with iAB ≤ 23. The redshifts are computed using a combination of narrow-band photometry from PAUS, a survey with 40 narrow bands spaced at $100\,\mathring{\rm A}$ intervals covering the range from 4500 to $8500\,\mathring{\rm A}$, 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 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 O ii 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 O ii line and Hα, Hβ and $\mathrm{O\,{\small III}}$ are similarly measured and used to generate the emission-line templates. Comparing to public spectroscopic surveys via the quantity Δz ≡ (zphoto − zspec)/(1 + zspec), we find the photometric redshifts to be more precise than previous estimates, with σ68(Δz) ≈ (0.003, 0.009) for galaxies at magnitude iAB ∼ 18 and iAB ∼ 23, respectively, which is three 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 |median(Δz)| ≤ 0.001 at all redshifts and magnitudes considered. Both the added PAUS data and new methodology contribute significantly to the improved results. The catalogue produced with the technique presented here is expected to provide a robust redshift calibration for current and future lensing surveys, and allows one to probe galaxy formation physics in an unexplored luminosity-redshift regime, thanks to its combination of depth, completeness, and excellent redshift precision and accuracy.

Cross-correlating radio continuum surveys and CMB lensing: constraining redshift distributions, galaxy bias, and cosmology

ABSTRACT We measure the harmonic-space auto-power spectrum of the galaxy overdensity in the LOFAR Two-metre Sky Survey (LoTSS) first data release and its cross-correlation with the map of the lensing convergence of the cosmic microwave background (CMB) from the Planck collaboration. We report a ∼5σ detection of the cross-correlation. We show that the combination of the clustering power spectrum and CMB lensing cross-correlation allows us to place constraints on the high-redshift tail of the redshift distribution, one of the largest sources of uncertainty in the use of continuum surveys for cosmology. Our analysis shows a preference for a broader redshift tail than that predicted by the photometric redshifts contained in the LoTSS value-added catalogue, as expected, and more compatible with predictions from simulations and spectroscopic data. Although the ability of CMB lensing to constrain the width and tail of the redshift distribution could also be valuable for the analysis of current and future photometric weak lensing surveys, we show that its performance relies strongly on the redshift evolution of the galaxy bias. Assuming the redshift distribution predicted by the Square Kilometre Array Design simulations, we use our measurements to place constraints on the linear bias of radio galaxies and the amplitude of matter inhomogeneities σ8, finding $\sigma _8=0.69^{+0.14}_{-0.21}$ assuming the galaxy bias scales with the inverse of the linear growth factor, and $\sigma _8=0.79^{+0.17}_{-0.32}$ assuming a constant bias.

Cosmological forecast for non-Gaussian statistics in large-scale weak lensing surveys

Cosmic shear data contains a large amount of cosmological information encapsulated in the non-Gaussian features of the weak lensing mass maps. Weak lensing studies mostly rely on two-point statistics to constrain cosmology from cosmic shear data, that do not capture all of this information. Additional non-Gaussian information can be extracted using non-Gaussian statistics. We compare the constraining power in the Ωm–σ8 plane of three map-based non-Gaussian statistics with the angular power spectrum, namely; peak counts, minimum counts and Minkowski functionals. We further analyze the impact of tomography and systematic effects originating from galaxy intrinsic alignments, multiplicative shear bias and photometric redshift systematics. We forecast the performance of the statistics for a stage-3-like weak lensing survey, spanning an area of 5000 deg2 and restrict ourselves to scales ⩾ 10 arcmin to avoid baryonic effects. The study follows a forward modelling scheme to predict the statistics at different cosmologies based on N-Body simulations. We find, that in our setup, the considered non-Gaussian statistics provide tighter constraints than the angular power spectrum. The peak counts show the greatest potential, increasing the figure-of-Merit (FoM) in the Ωm–σ8 plane by a factor of about 6, while the minimum counts and the Minkowski functionals yield an increase by a factor of about 3 and 2, respectively. A combined analysis using all non-Gaussian statistics in addition to the power spectrum increases the FoM by a factor of 9 and reduces the error on S8 by ≈ 30%. We find that the importance of tomography is diminished when combining non-Gaussian statistics with the angular power spectrum. The non-Gaussian statistics indeed profit less from tomography and the minimum counts and Minkowski functionals add some robustness against galaxy intrinsic alignment in a non-tomographic setting. We further find that a combination of the angular power spectrum and the non-Gaussian statistics allows us to apply conservative scale cuts in the analysis, thus helping to minimize the impact of baryonic and relativistic effects, while conserving the cosmological constraining power. We make the code that was used to conduct this analysis publicly available to simplify performing such analyses in the future.

Minimizing the impact of scale-dependent galaxy bias on the joint cosmological analysis of large-scale structures

ABSTRACT We present a mitigation strategy to reduce the impact of non-linear galaxy bias on the joint ‘3 × 2pt’ cosmological analysis of weak lensing and galaxy surveys. The Ψ-statistics that we adopt are based on Complete Orthogonal Sets of E/B Integrals. As such they are designed to minimize the contributions to the observable from the smallest physical scales where models are highly uncertain. We demonstrate that Ψ-statistics carry the same constraining power as the standard two-point galaxy clustering and galaxy–galaxy lensing statistics, but are significantly less sensitive to scale-dependent galaxy bias. Using two galaxy bias models, motivated by halo model fits to data and simulations, we quantify the error in a standard 3 × 2pt analysis where constant galaxy bias is assumed. Even when adopting conservative angular scale cuts, that degrade the overall cosmological parameter constraints, we find of order 1σ biases for Stage III surveys on the cosmological parameter S8 = σ8(Ωm/0.3)α. This arises from a leakage of the smallest physical scales to all angular scales in the standard two-point correlation functions. In contrast, when analysing Ψ-statistics under the same approximation of constant galaxy bias, we show that the bias on the recovered value for S8 can be decreased by a factor of ∼2, with less conservative scale cuts. Given the challenges in determining accurate galaxy bias models in the highly non-linear regime, we argue that 3 × 2pt analyses should move towards new statistics that are less sensitive to the smallest physical scales.

Photo-z outlier self-calibration in weak lensing surveys

Calibrating photometric redshift errors in weak lensing surveys with external data is extremely challenging. We show that both Gaussian and outlier photo-z parameters can be self-calibrated from the data alone. This comes at no cost for the neutrino masses, curvature and dark energy equation of state w0, but with a 65% degradation when both w0 and wa are varied. We perform a realistic forecast for the Vera Rubin Observatory (VRO) Legacy Survey of Space and Time (LSST) 3× 2 analysis, combining cosmic shear, projected galaxy clustering and galaxy—galaxy lensing. We confirm the importance of marginalizing over photo-z outliers. We examine a subset of internal cross-correlations, dubbed “null correlations”, which are usually ignored in 3× 2 analyses. Despite contributing only ∼ 10% of the total signal-to-noise, these null correlations improve the constraints on photo-z parameters by up to an order of magnitude. Using the same galaxy sample as sources and lenses dramatically improves the photo-z uncertainties too. Together, these methods add robustness to any claim of detected new Physics, and reduce the statistical errors on cosmology by 15% and 10% respectively. Finally, including CMB lensing from an experiment like Simons Observatory or CMB-S4 improves the cosmological and photo-z posterior constraints by about 10%, and further improves the robustness to systematics. To give intuition on the Fisher forecasts, we examine in detail several toy models that explain the origin of the photo-z self-calibration. Our Fisher code LaSSI (Large-Scale Structure Information), which includes the effect of Gaussian and outlier photo-z, shear multiplicative bias, linear galaxy bias, and extensions to LCDM, is publicly available at \\href{https://github.com/EmmanuelSchaan/LaSSI}{https://github.com/EmmanuelSchaan/LaSSI}.

Optimizing tomography for weak gravitational lensing surveys

ABSTRACT The subject of this paper is optimization of weak lensing tomography: we carry out numerical minimization of a measure of total statistical error as a function of the redshifts of the tomographic bin edges by means of a Nelder–Mead algorithm in order to optimize the sensitivity of weak lensing with respect to different optimization targets. Working under the assumption of a Gaussian likelihood for the parameters of a w0wa CDM (cold dark matter) model and using euclid’s conservative survey specifications, we compare an equipopulated, equidistant, and optimized bin setting and find that in general the equipopulated setting is very close to the optimal one, while an equidistant setting is far from optimal and also suffers from the ad hoc choice of a maximum redshift. More importantly, we find that nearly saturated information content can be gained using already few tomographic bins. This is crucial for photometric redshift surveys with large redshift errors. We consider a large range of targets for the optimization process that can be computed from the parameter covariance (or equivalently, from the Fisher matrix), extend these studies to information entropy measures such as the Kullback–Leibler divergence and conclude that in many cases equipopulated binning yields results close to the optimum, which we support by analytical arguments.

A synthetic Roman Space Telescope High-Latitude Imaging Survey: simulation suite and the impact of wavefront errors on weak gravitational lensing

ABSTRACT The Nancy Grace Roman Space Telescope (Roman) mission is expected to launch in the mid-2020s. Its weak lensing program is designed to enable unprecedented systematics control in photometric measurements, including shear recovery, point spread function (PSF) correction, and photometric calibration. This will enable exquisite weak lensing science and allow us to adjust to and reliably contribute to the cosmological landscape after the initial years of observations from other concurrent Stage IV dark energy experiments. This potential requires equally careful planning and requirements validation as the mission prepares to enter its construction phase. We present a suite of image simulations based on galsim that are used to construct a complex, synthetic Roman weak lensing survey that incorporates realistic input galaxies and stars, relevant detector non-idealities, and the current reference 5-yr Roman survey strategy. We present a first study to empirically validate the existing Roman weak lensing requirements flowdown using a suite of 12 matched image simulations, each representing a different perturbation to the wavefront or image motion model. These are chosen to induce a range of potential static and low- and high-frequency time-dependent PSF model errors. We analyse the measured shapes of galaxies from each of these simulations and compare them to a reference, fiducial simulation to infer the response of the shape measurement to each of these modes in the wavefront model. We then compare this to existing analytic flowdown requirements, and find general agreement between the empirically derived response and that predicted by the analytic model.

Sparse Bayesian mass-mapping with uncertainties: Full sky observations on the celestial sphere

ABSTRACT To date weak gravitational lensing surveys have typically been restricted to small fields of view, such that the flat-sky approximation has been sufficiently satisfied. However, with Stage IV surveys (e.g. LSST and Euclid) imminent, extending mass-mapping techniques to the sphere is a fundamental necessity. As such, we extend the sparse hierarchical Bayesian mass-mapping formalism presented in previous work to the spherical sky. For the first time, this allows us to construct maximum a posteriori spherical weak lensing dark-matter mass-maps, with principled Bayesian uncertainties, without imposing or assuming Gaussianty. We solve the spherical mass-mapping inverse problem in the analysis setting adopting a sparsity promoting Laplace-type wavelet prior, though this theoretical framework supports all log-concave posteriors. Our spherical mass-mapping formalism facilitates principled statistical interpretation of reconstructions. We apply our framework to convergence reconstruction on high resolution N-body simulations with pseudo-Euclid masking, polluted with a variety of realistic noise levels, and show a significant increase in reconstruction fidelity compared to standard approaches. Furthermore, we perform the largest joint reconstruction to date of the majority of publicly available shear observational data sets (combining DESY1, KiDS450, and CFHTLens) and find that our formalism recovers a convergence map with significantly enhanced small-scale detail. Within our Bayesian framework we validate, in a statistically rigorous manner, the community’s intuition regarding the need to smooth spherical Kaiser-Squires estimates to provide physically meaningful convergence maps. Such approaches cannot reveal the small-scale physical structures that we recover within our framework.

How accurately can we detect the splashback radius of dark matter haloes and its correlation with accretion rate?

ABSTRACT The splashback radius (Rsp) of dark matter haloes has recently been detected using weak gravitational lensing and cross-correlations with galaxies. However, different methods have been used to measure Rsp and to assess the significance of its detection. In this paper, we use simulations to study the precision and accuracy to which we can detect the splashback radius with 3D density, 3D subhalo, and weak lensing profiles. We study how well various methods and tracers recover Rsp by comparing it with the value measured directly from particle dynamics. We show that estimates of Rsp from density and subhalo profiles correspond to different percentiles of the underlying Rsp distribution of particle orbits. At low accretion rates, a second caustic appears and can bias results. Finally, we show that upcoming lensing surveys may be able to constrain the Rsp–Γdyn relation directly.

Void halo mass function: A promising probe of neutrino mass

Cosmic voids, the underdense regions in the universe, are particularly sensitive to diffuse density components such as cosmic neutrinos. This sensitivity is enhanced by the match between void sizes and the free-streaming scale of massive neutrinos. Using the massive neutrino simulations massivenus, we investigate the effect of neutrino mass on dark matter halos as a function of environment. We find that the halo mass function depends strongly on neutrino mass and that this dependence is more pronounced in voids than in high-density environments. An observational program that measured the characteristic mass of the most massive halos in voids should be able to place novel constraints on the sum of the masses of neutrinos $\ensuremath{\sum}{m}_{\ensuremath{\nu}}$. The neutrino mass effect in the simulations is quite strong: In a ${512}^{3}\text{ }\text{ }{h}^{\ensuremath{-}3}\text{ }{\mathrm{Mpc}}^{3}$ survey, the mean mass of the 1000 most massive halos in the void interiors is $(4.82\ifmmode\pm\else\textpm\fi{}0.11)\ifmmode\times\else\texttimes\fi{}{10}^{12}\text{ }\text{ }{h}^{\ensuremath{-}1}\text{ }{M}_{\ensuremath{\bigodot}}$ for $\ensuremath{\sum}{m}_{\ensuremath{\nu}}=0.6\text{ }\text{ }\mathrm{eV}$ and $(8.21\ifmmode\pm\else\textpm\fi{}0.13)\ifmmode\times\else\texttimes\fi{}{10}^{12}\text{ }\text{ }{h}^{\ensuremath{-}1}\text{ }{M}_{\ensuremath{\bigodot}}$ for $\ensuremath{\sum}{m}_{\ensuremath{\nu}}=0.1\text{ }\text{ }\mathrm{eV}$. Subaru (SuMIRe), Euclid and WFIRST will have both spectroscopic and weak lensing surveys. Covering volumes at least 50 times larger than our simulations, they should be sensitive probes of neutrino mass through void substructure.

Accessing the high- ℓ frontier under the reduced shear approximation with k -cut cosmic shear

The precision of Stage IV cosmic shear surveys will enable us to probe smaller physical scales than ever before, however, model uncertainties from baryonic physics and non-linear structure formation will become a significant concern. The $k$-cut method -- applying a redshift-dependent $\ell$-cut after making the Bernardeau-Nishimichi-Taruya transform -- can reduce sensitivity to baryonic physics; allowing Stage IV surveys to include information from increasingly higher $\ell$-modes. Here we address the question of whether it can also mitigate the impact of making the reduced shear approximation; which is also important in the high-$\kappa$, small-scale regime. The standard procedure for relaxing this approximation requires the repeated evaluation of the convergence bispectrum, and consequently can be prohibitively computationally expensive when included in Monte Carlo analyses. We find that the $k$-cut cosmic shear procedure suppresses the $w_0w_a$CDM cosmological parameter biases expected from the reduced shear approximation for Stage IV experiments, when $\ell$-modes up to $5000$ are probed. The maximum cut required for biases from the reduced shear approximation to be below the threshold of significance is at $k = 5.37 \, h{\rm Mpc}^{-1}$. With this cut, the predicted $1\sigma$ constraints increase, relative to the case where the correction is directly computed, by less than $10\%$ for all parameters. This represents a significant improvement in constraints compared to the more conservative case where only $\ell$-modes up to 1500 are probed, and no $k$-cut is used. We also repeat this analysis for a hypothetical, comparable kinematic weak lensing survey. The key parts of code used for this analysis are made publicly available.

Deep＋ wide lensing surveys can measure the dark matter halos of dwarf galaxies

The advent of new deep+wide photometric lensing surveys will open up the possibility of direct measurements of the dark matter halos of dwarf galaxies. The HSC wide survey will be the first with the statistical capability of measuring the lensing signal with high signal-to-noise at log(M*)=8. At this same mass scale, LSST will have the most overall constraining power with a predicted signal-to-noise for the galaxy-galaxy lensing signal around dwarfs of SN=200. WFIRST and LSST will have the greatest potential to push below the log(M*) = 7 mass scale thanks to the depth of their imaging data. Studies of the dark matter halos of dwarf galaxies at z=0.1 with gravitational lensing are soon within reach. However, further work will be required to develop optimized strategies for extracting dwarfs samples from these surveys, determining redshifts, and accurately measuring lensing on small radial scales. Dwarf lensing will be a new and powerful tool to constrain the halo masses and inner density slopes of dwarf galaxies and to distinguish between baryonic feedback and modified dark matter scenarios.

Fast estimation of aperture mass statistics – I. Aperture mass variance and an application to the CFHTLenS data

ABSTRACT We explore an alternative method to the usual shear correlation function approach for the estimation of aperture mass statistics in weak lensing survey data. Our approach builds on the direct estimator method. In this paper, to test and validate the methodology, we focus on the aperture mass dispersion. After computing the signal and noise for a weighted set of measured ellipticites we show how the direct estimator can be made into a linear order algorithm that enables a fast and efficient computation. We then investigate the applicability of the direct estimator approach in the presence of a real survey mask with holes and chip gaps. For this we use a large ensemble of full ray-tracing mock simulations. By using various weighting schemes for combining information from different apertures we find that inverse variance weighting the individual aperture estimates with an aperture completeness greater than 70 per cent yields an answer that is in close agreement with the correlation function approach. We then apply this approach to the CFHTLenS as a pilot scheme and find that our method recovers to high accuracy the official result for the variance of both the E- and B-mode signal. We then explore the cosmological information content of the direct estimator using the Fisher information approach. We show that there is a only modest loss in cosmological information from the rejection of apertures that are of low completeness. This method unlocks the door to efficient methods for recovering higher order aperture mass statistics in linear order operations.

Higher order spectra of weak lensing convergence maps in parametrized theories of modified gravity

ABSTRACT We compute the low-ℓ limit of the family of higher order spectra for projected (2D) weak lensing convergence maps. In this limit these spectra are computed to an arbitrary order using tree-level perturbative calculations. We use the flat-sky approximation and Eulerian perturbative results based on a generating function approach. We test these results for the lower order members of this family, i.e. the skew- and kurt-spectra against state-of-the-art simulated all-sky weak lensing convergence maps and find our results to be in very good agreement. We also show how these spectra can be computed in the presence of a realistic sky-mask and Gaussian noise. We generalize these results to 3D and compute the equal-time higher order spectra. These results will be valuable in analysing higher order statistics from future all-sky weak lensing surveys such as the Euclid survey at low-ℓ modes. As illustrative examples, we compute these statistics in the context of the Horndeski and beyond Horndeski theories of modified gravity. They will be especially useful in constraining theories such as the Gleyzes–Langlois–Piazza–Vernizzi (GLPV) theories and degenerate higher order scalar-tensor theories as well as the commonly used normal-branch of Dvali–Gabadadze–Porrati model, clustering quintessence models and scenarios with massive neutrinos.

Velocity Dispersions of Massive Quiescent Galaxies from Weak Lensing and Spectroscopy*

Abstract We use MMT spectroscopy and deep Subaru Hyper Suprime-Cam (HSC) imaging to compare the spectroscopic central stellar velocity dispersion of quiescent galaxies with the effective dispersion of the dark matter halo derived from the stacked lensing signal. The spectroscopic survey (the Smithsonian Hectospec Lensing Survey) provides a sample of 4585 quiescent galaxy lenses with measured line-of-sight central stellar velocity dispersion (σ SHELS) that is more than 85% complete for R < 20.6, D n 4000 > 1.5 and M ⋆ > 109.5 M ⊙. The median redshift of the sample of lenses is 0.32. We measure the stacked lensing signal from the HSC deep imaging. The central stellar velocity dispersion is directly proportional to the velocity dispersion derived from the lensing σ Lens, . The independent spectroscopic and weak lensing velocity dispersions probe different scales, ∼3 kpc and ≳100 kpc, respectively, and strongly indicate that the observable central stellar velocity dispersion for quiescent galaxies is a good proxy for the velocity dispersion of the dark matter halo. We thus demonstrate the power of combining high-quality imaging and spectroscopy to shed light on the connection between galaxies and their dark matter halos.

Correlations between triaxial shapes and formation history of dark matter haloes

ABSTRACT The shape of dark matter (DM) haloes plays a critical role in constraining cosmology with upcoming large-scale structure surveys. In this paper, we study the correlations between the triaxial shapes and formation histories in DM haloes in the MultiDark Planck 2 N-body cosmological simulation. We find that halo ellipticity is strongly correlated with halo properties that serve as proxies of halo formation history, such as halo concentration and the normalized peak-centroid offset. These correlations are nearly independent of the halo density peak height. We present a simple model for the correlation between halo ellipticity and concentration using conditional abundance matching, and provide fitting formulae for the multidimensional distributions of triaxial halo shape as a function of halo peak height. We apply our halo shape model to gauge the effects of halo ellipticity and orientation bias on the excess surface mass density profiles in cluster-sized haloes. Our model should be useful for exploring the impact of triaxial halo shape on cosmological constraints in upcoming weak lensing surveys of galaxy clusters.

Molecular gas in CLASH brightest cluster galaxies at z ∼ 0.2 – 0.9

Brightest cluster galaxies (BCGs) are excellent laboratories for the study of galaxy evolution in dense Mpc-scale environments. We used the IRAM-30 m to observe, in CO(1→0), CO(2→1), CO(3→2), or CO(4→3), 18 BCGs at z ∼ 0.2 − 0.9 drawn from the Cluster Lensing And Supernova survey with Hubble (CLASH) survey. Our sample includes RX1532, which is our primary target as it is among the BCGs with the highest star formation rate (SFR ≳100 M⊙ yr−1) in the CLASH sample. We unambiguously detected both CO(1→0) and CO(3→2) in RX1532, yielding a large reservoir of molecular gas, MH2 = (8.7 ± 1.1)×1010 M⊙, and a high level of excitation, r31 = 0.75 ± 0.12. A morphological analysis of the Hubble Space Telescope I-band image of RX1532 reveals the presence of clumpy substructures both within and outside the half-light radius re = (11.6 ± 0.3) kpc, similarly to those found independently both in ultraviolet and in Hα in previous works. We tentatively detected CO(1→0) or CO(2→1) in four other BCGs, with molecular gas reservoirs in the range of MH2 = 2 × 1010 − 11 M⊙. For the remaining 13 BCGs, we set robust upper limits of MH2/M⋆ ≲ 0.1, which are among the lowest molecular-gas-to-stellar-mass ratios found for distant ellipticals and BCGs. In comparison with distant cluster galaxies observed in CO, our study shows that RX1532 (MH2/M⋆ = 0.40 ± 0.05) belongs to the rare population of star-forming and gas-rich BCGs in the distant universe. By using the available X-ray based estimates of the central intra-cluster medium entropy, we show that the detection of large reservoirs of molecular gas MH2 ≳ 1010 M⊙ in distant BCGs is possible when the two conditions are met: (i) high SFR and (ii) low central entropy, which favors the condensation and the inflow of gas onto the BCGs themselves, similarly to what has been previously found for some local BCGs.

Parameter inference for weak lensing using Gaussian Processes and MOPED

ABSTRACT In this paper, we propose a Gaussian Process (GP) emulator for the calculation both of tomographic weak lensing band powers, and of coefficients of summary data massively compressed with the MOPED algorithm. In the former case cosmological parameter inference is accelerated by a factor of ∼10–30 compared with Boltzmann solver class applied to KiDS-450 weak lensing data. Much larger gains of order 103 will come with future data, and MOPED with GPs will be fast enough to permit the Limber approximation to be dropped, with acceleration in this case of ∼105. A potential advantage of GPs is that an error on the emulated function can be computed and this uncertainty incorporated into the likelihood. However, it is known that the GP error can be unreliable when applied to deterministic functions, and we find, using the Kullback–Leibler divergence between the emulator and class likelihoods, and from the uncertainties on the parameters, that agreement is better when the GP uncertainty is not used. In future, weak lensing surveys such as Euclid, and the Legacy Survey of Space and Time, will have up to ∼104 summary statistics, and inference will be correspondingly more challenging. However, since the speed of MOPED is determined not the number of summary data, but by the number of parameters, MOPED analysis scales almost perfectly, provided that a fast way to compute the theoretical MOPED coefficients is available. The GP provides such a fast mechanism.