Search

The masses of galaxy clusters can be measured using data obtained exclusively from wide photometric surveys in one of two ways: directly from the amplitude of the weak lensing signal or, indirectly, through the use of scaling relations calibrated using binned lensing measurements. In this paper, we build on a recently proposed idea and implement an alternative method based on the radial profile of the satellite distribution. This technique relies on splashback, a feature associated with the apocentre of recently accreted galaxies that offers a clear window into the phase-space structure of clusters without the use of velocity information. We carry out this dynamical measurement using the stacked satellite distribution around a sample of luminous red galaxies in the fourth data release of the Kilo-Degree Survey and validate our results using abundance-matching and lensing masses. To illustrate the power of this measurement, we combine dynamical and lensing mass estimates to robustly constrain scalar–tensor theories of gravity at cluster scales. Our results exclude departures from General Relativity of the order of unity. We conclude the paper by discussing the implications for future data sets. Because splashback mass measurements scale only with the survey volume, stage-IV photometric surveys are well-positioned to use splashback to provide high-redshift cluster masses.

Around $10^5$ strongly lensed galaxies are expected to be discovered with Euclid and the LSST. Utilising these large samples to study the inner structure of lens galaxies requires source redshifts, to turn lens models into mass measurements. However, obtaining spectroscopic source redshifts for large lens samples is prohibitive with the capacity of spectroscopic facilities. Alternatively, we study the possibility of obtaining source photometric redshifts (photo-zs) for large lens samples. Our strategy consists of deblending the source and lens light by simultaneously modelling the lens and background source in all available photometric bands, and feeding the derived source colours to a template-fitting photo-z algorithm. We describe the lens and source light with a Sersic profile, and the lens mass with a Singular Isothermal Ellipsoid. We test our approach on a simulated and a real sample of lenses, both in broad-band photometry of the Hyper Suprime-Cam survey. We identify the deviations of the lens light from a Sersic profile and the contrast between the lens and source image as the main drivers of the source colour measurement error. We split the real sample based on the ratio $\Lambda$ of the lens to source surface brightness measured at the image locations. In the $\Lambda, Comment: 15 pages, 13 figures, Accepted for Publication in A&A

The diffuse light within galaxy groups and clusters provides valuable insight into the growth of massive cosmic structures. Groups are particularly interesting in this context, because they represent the link between galactic haloes and massive clusters. However, low surface brightness makes their diffuse light extremely challenging to detect individually. Stacking many groups is a promising alternative, but its physical interpretation is complicated by possible systematic variations of diffuse light profiles with other group properties. Another issue is the often ambiguous choice of group centre. We explore these challenges using mock observations for 497 galaxy groups and clusters with halo masses from $~ 10^{12} \textrm{M}_{\odot}$ to $1.5 \times 10^{15}\textrm{M}_{\odot}$ at redshift $0.1$ from the Hydrangea cosmological hydrodynamic simulations. In 18 per cent of groups with at least five galaxies above $10^{9} \textrm{M}_{\odot}$ in stellar mass, the $r$-band brightest galaxy is not the one at the centre of the gravitational potential; line-of-sight projections account for half of these cases. Miscentring does not significantly affect the ensemble average mass density profile or the surface brightness profile for our sample: even within ambiguously centred haloes, different centring choices lead to only a 1 per cent change in the total fraction of diffuse intra-group light, $f_{\textrm{IGL}}$. We find strong correlations of $f_{\textrm{IGL}}$ with the luminosity of the central group galaxy and halo mass. Stacking groups in narrow bins of central galaxy luminosity will therefore make the physical interpretation of the signal more straightforward than combining systems across a wide range of mass., 18 pages, 12 figures, Accepted for publication in MNRAS

Expanding from previous work we present weak lensing measurements for a total sample of 30 distant ($z_\mathrm{median}=0.93$) massive galaxy clusters from the South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) Survey, measuring galaxy shapes in Hubble Space Telescope (HST) Advanced Camera for Surveys images. We remove cluster members and preferentially select $z\gtrsim 1.4$ background galaxies via $V-I$ colour, employing deep photometry from VLT/FORS2 and Gemini-South/GMOS. We apply revised calibrations for the weak lensing shape measurements and the source redshift distribution to estimate the cluster masses. In combination with earlier Magellan/Megacam results for lower-redshifts clusters we infer refined constraints on the scaling relation between the SZ detection significance and the cluster mass, in particular regarding its redshift evolution. The mass scale inferred from the weak lensing data is lower by a factor $0.76^{+0.10}_{-0.14}$ (at our pivot redshift $z=0.6$) compared to what would be needed to reconcile a flat Planck $\nu\Lambda$CDM cosmology (in which the sum of the neutrino masses is a free parameter) with the observed SPT-SZ cluster counts. In order to sensitively test the level of (dis-)agreement between SPT clusters and Planck, further expanded weak lensing follow-up samples are needed., Comment: 32 pages, 20 figures, 12 tables, matches the version accepted for publication in MNRAS

Galaxy clusters are excellent probes to study the effect of environment on galaxy formation and evolution. Along with high-quality observational data, accurate cosmological simulations are required to improve our understanding of galaxy evolution in these systems. In this work, we compare state-of-the-art observational data of massive galaxy clusters ($>10^{14} \textrm{M}_{\odot}$) at different redshifts ($010^{14} \textrm{M}_{\odot}$ at $z=0$). We compare three fundamental observables of galaxy clusters: the total stellar mass to halo mass ratio, the stellar mass function (SMF), and the radial mass density profile of the cluster galaxies. In the first two of these, the simulations agree well with the observations, albeit with a slightly too high abundance of $M_\star \lesssim 10^{10} \textrm{M}_{\odot}$ galaxies at $z \gtrsim 1$. The NFW concentrations of cluster galaxies increase with redshift, in contrast to the decreasing dark matter halo concentrations. This previously observed behaviour is therefore due to a qualitatively different assembly of the smooth DM halo compared to the satellite population. Quantitatively, we however find a discrepancy in that the simulations predict higher stellar concentrations than observed at lower redshifts ($z, Comment: 14 pages, 9 figures (excluding appendices), Accepted for publication in MNRAS

This report is the result of a joint discussion between the Rubin andEuclid scientific communities.The workpresented in this report was focused on designing and recommending an initial set of Derived Data products (DDPs) that could realize the science goals enabled by joint processing. All interested Rubin and Euclid data rights holders were invited to contribute via an online discussion forum and a series of virtual meetings. Strong interest in enhancing science with joint DDPs emerged from across a wide range of astrophysical domains: Solar System, the Galaxy, the Local Volume, from the nearby to the primaeval Universe, and cosmology.&nbsp

We present SKiLLS, a suite of multi-band image simulations for the weak lensing analysis of the complete Kilo-Degree Survey (KiDS), dubbed KiDS-Legacy analysis. The resulting catalogues enable joint shear and redshift calibration, enhancing the realism and hence accuracy over previous efforts. To create a large volume of simulated galaxies with faithful properties and to a sufficient depth, we integrated cosmological simulations with high-quality imaging observations. We also improved the realism of simulated images by allowing the point spread function (PSF) to differ between CCD images, including stellar density variations and varying noise levels between pointings. Using realistic variable shear fields, we accounted for the impact of blended systems at different redshifts. Although the overall correction is minor, we found a clear redshift-bias correlation in the blending-only variable shear simulations, indicating the non-trivial impact of this higher-order blending effect. We also explored the impact of the PSF modelling errors and found a small yet noticeable effect on the shear bias. Finally, we conducted a series of sensitivity tests, including changing the input galaxy properties. We conclude that our fiducial shape measurement algorithm, lensfit, is robust within the requirements of lensing analyses with KiDS. As for future weak lensing surveys with tighter requirements, we suggest further investments in understanding the impact of blends at different redshifts, improving the PSF modelling algorithm and developing the shape measurement method to be less sensitive to the galaxy properties., Comment: 28 pages, 31 figures, 2 tables, minor revisions to match the final accepted version

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 \xi_\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 \xi_\pm^E/B(\vt) and their ambiguous components \xi_\pm^amb(\vt), that can be derived from the measured \xi_\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 \xi_+=\xi_+^E+\xi_+^B+\xi_+^amb and \xi_-=\xi_-^E-\xi_-^B+\xi_-^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 \xi_\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., Comment: 24 pages, A&A, in press

Stanford, S. A., et al., The Complete Calibration of the Color–Redshift Relation (C3R2) survey is obtaining spectroscopic redshifts in order to map the relation between galaxy color and redshift to a depth of i ∼ 24.5 (AB). The primary goal is to enable sufficiently accurate photometric redshifts for Stage iv dark energy projects, particularly Euclid and the Nancy Grace Roman Space Telescope (Roman), which are designed to constrain cosmological parameters through weak lensing. We present 676 new high-confidence spectroscopic redshifts obtained by the C3R2 survey in the 2017B–2019B semesters using the DEIMOS, LRIS, and MOSFIRE multiobject spectrographs on the Keck telescopes. Combined with the 4454 redshifts previously published by this project, the C3R2 survey has now obtained and published 5130 high-quality galaxy spectra and redshifts. If we restrict consideration to only the 0.2 < zp < 2.6 range of interest for the Euclid cosmological goals, then with the current data release, C3R2 has increased the spectroscopic redshift coverage of the Euclid color space from 51% (as reported by Masters et al.) to the current 91%. Once completed and combined with extensive data collected by other spectroscopic surveys, C3R2 should provide the spectroscopic calibration set needed to enable photometric redshifts to meet the cosmology requirements for Euclid, and make significant headway toward solving the problem for Roman.

We present measurements of the radial gravitational acceleration around isolated galaxies, comparing the expected gravitational acceleration given the baryonic matter with the observed gravitational acceleration, using weak lensing measurements from the fourth data release of the Kilo-Degree Survey. These measurements extend the radial acceleration relation (RAR) by 2 decades into the low-acceleration regime beyond the outskirts of the observable galaxy. We compare our RAR measurements to the predictions of two modified gravity (MG) theories: MOND and Verlinde's emergent gravity. We find that the measured RAR agrees well with the MG predictions. In addition, we find a difference of at least $6\sigma$ between the RARs of early- and late-type galaxies (split by S\'{e}rsic index and $u-r$ colour) with the same stellar mass. Current MG theories involve a gravity modification that is independent of other galaxy properties, which would be unable to explain this behaviour. The difference might be explained if only the early-type galaxies have significant ($M_{gas} \approx M_*$) circumgalactic gaseous haloes. The observed behaviour is also expected in $\Lambda$CDM models where the galaxy-to-halo mass relation depends on the galaxy formation history. We find that MICE, a $\Lambda$CDM simulation with hybrid halo occupation distribution modelling and abundance matching, reproduces the observed RAR but significantly differs from BAHAMAS, a hydrodynamical cosmological galaxy formation simulation. Our results are sensitive to the amount of circumgalactic gas; current observational constraints indicate that the resulting corrections are likely moderate. Measurements of the lensing RAR with future cosmological surveys will be able to further distinguish between MG and $\Lambda$CDM models if systematic uncertainties in the baryonic mass distribution around galaxies are reduced., Comment: 31 pages, 15 figures, published in A&A

21 pags., 12 figs., 4 tabs., The See Change survey was designed to make z > 1 cosmological measurements by efficiently discovering high-redshift Type Ia supernovae (SNe Ia) and improving cluster mass measurements through weak lensing. This survey observed twelve galaxy clusters with the Hubble Space Telescope (HST) spanning the redshift range z = 1.13-1.75, discovering 57 likely transients and 27 likely SNe Ia at z ~ 0.8-2.3. As in similar previous surveys, this proved to be a highly efficient use of HST for supernova observations; the See Change survey additionally tested the feasibility of maintaining, or further increasing, the efficiency at yet higher redshifts, where we have less detailed information on the expected cluster masses and star formation rates. We find that the resulting number of SNe Ia per orbit is a factor of ~8 higher than for a field search, and 45% of our orbits contained an active SN Ia within 22 rest-frame days of peak, with one of the clusters by itself yielding 6 of the SNe Ia. We present the survey design, pipeline, and supernova discoveries. Novel features include fully blinded supernova searches, the first random forest candidate classifier for undersampled IR data (with a 50% detection threshold within 0.05 mag of human searchers), real-time forward-modeling photometry of candidates, and semi-automated photometric classifications and follow-up forecasts. We also describe the spectroscopic follow-up, instrumental in measuring host galaxy redshifts. The cosmology analysis of our sample will be presented in a companion paper. * Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by AURA, Inc., under NASA contract NAS 5-26555, under programs 13677, 14327., Based in part on observations obtained at the Gemini Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science Foundation (United States), National Research Council (Canada), CONICYT (Chile), Ministerio de Ciencia, Tecnología e Innovación Productiva (Argentina), Ministério da Ciência, Tecnologia e Inovação (Brazil), and Korea Astronomy and Space Science Institute (Republic of Korea). Based on observations collected at the Table 4. European Organisation for Astronomical Research in the Southern Hemisphere under ESO program(s) 294.A-5025(A), 095.A-0830(A, B, C), 096.A-0926(B, C), 097.A-0442(A, B, C), and 0100.A-0851(A). G.W. acknowledges support from the National Science Foundation through grant AST-1517863, by HST program numbers GO-13677/14327.01 and GO-15294, and by grant No. 80NSSC17K0019 issued through the NASA Astrophysics Data Analysis Program (ADAP). G.A. and R.G. acknowledge support from HST program NASA HST-GO14163.002-A, and by grant No. NNH16AC25I issued through NASA ADAP. M. J. J. acknowledges support for the current research from the National Research Foundation (NRF) of Korea under the programs 2017R1A2B2004644 and 2020R1A4A2002885. H. Hildebrandt is supported by a Heisenberg grant of the Deutsche Forschungsgemeinschaft (Hi 1495/5-1) as well as an ERC Consolidator Grant (No. 770935). Support for program numbers GO-13677/14327.01 and GO-15294 was provided by NASA through a grant from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Incorporated, under NASA contract NAS5-26555. This work was also partially supported by the Office of Science, Office of High Energy Physics, of the U.S. Department of Energy, under contract no. DE-AC02-05CH11231. Facilities: Hubble Space Telescope, Keck:I (LRIS), Keck:I (MOSFIRE), VLT, Gemini, GTC, Subaru. Software: astropy (Astropy Collaboration 2013), DrizzlePac (Gonzaga et al. 2012), iPython (Pérez & Granger 2007), Matplotlib (Hunter 2007), Numpy (van der Walt et al. 2011), SEP (Barbary 2016), scikit-learn (Pedregosa et al. 2011), SciPy

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. 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 S8 = σ8? (ωm/0.3) = 0.766-0.014+0.020, which has the same overall precision as that reported by the full-sky cosmic microwave background observations from Planck. The recovered S8 amplitude is low, however, by 8.3? ±? 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 S8-difference with Planck is driven by a tension in the matter fluctuation amplitude parameter, σ8. We quantify the level of agreement between the cosmic microwave background and our large-scale structure constraints using a series of different metrics, finding differences with a significance ranging between ∼3σ, when considering the offset in S8, and ∼2σ, when considering the full multi-dimensional parameter space., Heymans, Catherine, et al.

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., Comment: 12 pages, 8 figures, 3 tables. Revised to match version published in A&A

Intrinsic alignments (IAs) of galaxies are an important contaminant for cosmic shear studies, but the modelling is complicated by the dependence of the signal on the source galaxy sample. In this paper, we use the halo model formalism to capture this diversity and examine its implications for Stage-III and Stage-IV cosmic shear surveys. We account for the different IA signatures at large and small scales, as well for the different contributions from central/satellite and red/blue galaxies, and we use realistic mocks to account for the characteristics of the galaxy populations as a function of redshift. We inform our model using the most recent observational findings: we include a luminosity dependence at both large and small scales and a radial dependence of the signal within the halo. We predict the impact of the total IA signal on the lensing angular power spectra, including the current uncertainties from the IA best-fits to illustrate the range of possible impact on the lensing signal: the lack of constraints for fainter galaxies is the main source of uncertainty for our predictions of the IA signal. We investigate how well effective models with limited degrees of freedom can account for the complexity of the IA signal. Although these lead to negligible biases for Stage-III surveys, we find that, for Stage-IV surveys, it is essential to at least include an additional parameter to capture the redshift dependence., 21 pages, 15 figure, accepted for publication in MNRAS

We present constraints on extensions to the standard cosmological model of a spatially flat Universe governed by general relativity, a cosmological constant (Λ), and cold dark matter (CDM) by varying the spatial curvature ωK, the sum of the neutrino masses χ mν, the dark energy equation of state parameter w, and the Hu-Sawicki f(R) gravity fR0 parameter. With the combined 3 × 2 pt measurements of cosmic shear from the Kilo-Degree Survey (KiDS-1000), galaxy clustering from the Baryon Oscillation Spectroscopic Survey (BOSS), and galaxy-galaxy lensing from the overlap between KiDS-1000, BOSS, and the spectroscopic 2-degree Field Lensing Survey, we find results that are fully consistent with a flat ΛCDM model with ωK = 0.011-0.057+0.054, χ mν, Tröster, T., et al.

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 f(h), the average ratio between the aligned component of the halo ellipticity and the ellipticity of the light distribution, finding, f(h) = 0.303(-0.079)(+0.080), f h, = 0 ., 303, 0.079, + 0.080, for red lens galaxies and, f(h) = 0.217(-0.159)(+0.160), 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 sigma) 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 |E-h|=0.174 +/- 0.046, where |E-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.

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., version accepted by A&A; 24 pages; included reference to companion paper Kannawadi et al. arXiv:2010.04164

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 $\Lambda$CDM model, we find $S_8 = \sigma_8 (\Omega_{\rm m}/0.3)^{0.5} = 0.759^{+0.024}_{-0.021}$ for our fiducial analysis, which is in $3\sigma$ 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., Comment: 33 pages, 14 figures, 7 tables, Accepted to be published in A&A. This paper is part of the KiDS-1000 series of papers: Heymans, Tr\"oster et al. (arXiv:2007.15632), Hildebrandt et al. (arXiv:2007.15635), Joachimi et al. (arXiv:2007.01844) and Giblin et al. (arXiv:2007.01845). Online KiDS-1000 talks can be viewed at http://kids.strw.leidenuniv.nl/KiDS-1000.php or the KiDS Consortium youtube page

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. We study the performance of a fully data-driven approach, based on MetaDetection, focusing on the more realistic case of observations with an anisotropic PSF. Under the assumption that PSF anisotropy is the only source of additive shear bias, we show how unbiased shear estimates can be obtained from the observed data alone. To do so, we exploit the finding that the multiplicative shear bias obtained with MetaDetection is nearly insensitive to the PSF ellipticity. In practice, this assumption can be validated by comparing the empirical corrections obtained from observations to those from simulated data. We show that our data-driven approach meets the stringent requirements for upcoming space and ground-based surveys, although further optimisation is possible., 11 pages, accepted for publication in Astronomy and Astrophysics

We measure the cross-correlation between galaxy weak lensing data from the Kilo Degree Survey (KiDS-1000, DR4) and cosmic microwave background (CMB) lensing data from the Atacama Cosmology Telescope (ACT, DR4) and the Planck Legacy survey. We use two samples of source galaxies, selected with photometric redshifts, $(0.1, Comment: 13 pages, 9 figures, 1 tables, submitted to A&A

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 $\Omega_\mathrm{m}$, $S_8=\sigma_8(\Omega_\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 $\Omega_\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., Comment: 17 pages, 11 figures. Updated wrong definition of spherical overdensity mass m200m, our calculations used the correct definition, as can be seen from the publicly-available code at https://github.com/StijnDebackere/lensing_haloes/blob/afe4f8699003/lensing_haloes/results.py#L311-L321 . Mock data, including cosmology MAPs, available at https://doi.org/10.5281/zenodo.4469436

The outskirts of accreting dark matter haloes exhibit a sudden drop in density delimiting their multi-stream region. Due to the dynamics of accretion, the location of this physically motivated edge strongly correlates with the halo growth rate. Using hydrodynamical zoom-in simulations of high-mass clusters, we explore this definition in realistic simulations and find an explicit connection between this feature in the dark matter and galaxy profiles. We also show that the depth of the splashback feature correlates well with the direction of filaments and, surprisingly, the orientation of the brightest cluster galaxy. Our findings suggest that galaxy profiles and weak-lensing masses can define an observationally viable mass-size scaling relation for galaxy clusters, which can be used to extract cosmological information., Comment: 9 pages, 9 figures. Fig. 6 and 9 updated with correct labels

We present measurements of $f_h$, the ratio of the aligned components of the projected halo and galaxy ellipticities, for a sample of central galaxies using weak gravitational lensing data from the Kilo-Degree Survey (KiDS). Using a lens galaxy shape estimation that is more sensitive to outer galaxy regions, we find $f_{\rm h}=0.50\pm0.20$ for our full sample and $f_{\rm h}=0.55\pm0.19$ for an intrinsically red (and therefore higher stellar-mass) sub-sample, rejecting the hypothesis of round halos and/or galaxies being un-aligned with their parent halo at $2.5\sigma$ and $2.9\sigma$, respectively. We quantify the 93.4% purity of our central galaxy sample using numerical simulations and overlapping spectroscopy from the Galaxy and Mass Assembly survey. This purity ensures that the interpretation of our measurements is not complicated by the presence of a significant fraction of satellite galaxies. Restricting our central galaxy ellipticity measurement to the inner isophotes, we find $f_{\rm h}=0.34\pm0.17$ for our red sub-sample, suggesting that the outer galaxy regions are more aligned with their dark matter halos compared to the inner regions. Our results are in agreement with previous studies and suggest that lower mass halos are rounder and/or less aligned with their host galaxy than samples of more massive galaxies, studied in galaxy groups and clusters., Comment: Accepted for publication in A&A

We use stellar and dynamical mass profiles, combined with a stellar population analysis, of 32 brightest cluster galaxies (BCGs) at redshifts of 0.05 $\leq z \leq$ 0.30, to place constraints on their stellar Initial Mass Function (IMF). We measure the spatially-resolved stellar population properties of the BCGs, and use it to derive their stellar mass-to-light ratios ($\Upsilon_{\star \rm POP}$). We find young stellar populations ($, Comment: 15 pages, 11 figures, accepted for publication in MNRAS

Aims. The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. The estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on various individual methodologies and numerical implementations, which were developed for different observational probes and for the combination thereof. In this paper we present validated forecasts, which combine both theoretical and observational ingredients for different cosmological probes. This work is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. Methods. We describe in detail the methods adopted for Fisher matrix forecasts, which were applied to galaxy clustering, weak lensing, and the combination thereof. We estimated the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations, some of which are made publicly available, in different programming languages, together with a reference training-set of input and output matrices for a set of specific models. These can be used by the reader to validate their own implementations if required. Results. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setting, for example flat or non-flat spatial cosmologies, or different cuts at non-linear scales. The numerical implementations are now reliable for these settings. We present the results for an optimistic and a pessimistic choice for these types of settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy figure of merit by at least a factor of three., 0, info:eu-repo/semantics/published

We investigate the stellar and dynamical mass profiles in the centres of 25 brightest cluster galaxies (BCGs) at redshifts of 0.05 $\leq z \leq$ 0.30. Our spectroscopy enables us to robustly measure the Gauss-Hermite higher order velocity moments $h_{3}$ and $h_{4}$, which we compare to measurements for massive early-type galaxies, and central group galaxies. We measure positive central values for $h_{4}$ for all the BCGs. We derive the stellar mass-to-light ratio ($\Upsilon_{\star \rm DYN}$), and velocity anisotropy ($\beta$) based on a Multi-Gaussian Expansion (MGE) and axisymmetric Jeans Anisotropic Methods (JAM, cylindrically- and spherically-aligned). We explicitly include a dark matter halo mass component, which is constrained by weak gravitational lensing measurements for these clusters. We find a strong correlation between anisotropy and velocity dispersion profile slope, with rising velocity dispersion profiles corresponding to tangential anisotropy and decreasing velocity dispersion profiles corresponding to radial anisotropy. The rising velocity dispersion profiles can also indicate a significant contribution from the intracluster light (ICL) to the total light (in projection) in the centre of the galaxy. For a small number of BCGs with rising velocity dispersion profiles, a variable stellar mass-to-light ratio can also account for the profile shape, instead of tangential anisotropy or a significant ICL contribution. We note that, for some BCGs, a variable $\beta_{z}(r)$ (from radial to tangential anisotropy) can improve the model fit to the observed kinematic profiles. The observed diversity in these properties illustrates that BCGs are not the homogeneous class of objects they are often assumed to be., Comment: 25 pages, accepted for publication in MNRAS

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). The images obtained with PAUCam are affected by scattered light: an optical effect consisting of light multiply reflected that deposits energy in specific detector regions affecting 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 per cent and up to 20 per cent 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 by 35 per cent for the best 20 per cent galaxies selected with a photometric quality parameter., Funding for PAUS has been provided by Durham University (via the ERC StG DEGAS-259586), ETH Zurich, Leiden University (via ERC StG ADULT-279396 and Netherlands Organisation for Scientific Research (NWO) Vici grant 639.043.512) and University College London. The PAUS participants from Spanish institutions are partially supported by MINECO under grants CSD2007-00060, AYA2015-71825, ESP2015-88861, FPA2015-68048, SEV-2016-0588, SEV-2016-0597, and MDM-2015-0509, some of which include ERDF funds from the European Union. IEEC and IFAE are partially funded by the CERCA program of the Generalitat de Catalunya. The PAU data center is hosted by the Port d’Informacio Cientifica (PIC), maintained through a collaboration of CIEMAT and IFAE, with additional support from Universitat Autonoma de Barcelona and ERDF. CosmoHub has been developed by PIC and was partially funded by the ‘Plan Estatal de Investigacion Cientifica y Tecnica y de Innovacion program of the Spanish government. We gratefully acknowledge the support of NVIDIA Corporation with the donation of the Titan V GPU used for this research. This project has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 776247. AA is supported by a Royal Society Wolfson Fellowship. MS has been supported by the National Science Centre (grant UMO2016/23/N/ST9/02963).

We present a tomographic cosmic shear analysis of the Kilo-Degree Survey (KiDS) combined with the VISTA Kilo-Degree Infrared Galaxy Survey (VIKING). This is the first time that a full optical to near-infrared data set has been used for a wide-field cosmological weak lensing experiment. This unprecedented data, spanning $450~$deg$^2$, allows us to improve significantly the estimation of photometric redshifts, such that we are able to include robustly higher-redshift sources for the lensing measurement, and - most importantly - solidify our knowledge of the redshift distributions of the sources. Based on a flat $\Lambda$CDM model we find $S_8\equiv\sigma_8\sqrt{\Omega_{\rm m}/0.3}=0.737_{-0.036}^{+0.040}$ in a blind analysis from cosmic shear alone. The tension between KiDS cosmic shear and the Planck-Legacy CMB measurements remains in this systematically more robust analysis, with $S_8$ differing by $2.3\sigma$. This result is insensitive to changes in the priors on nuisance parameters for intrinsic alignment, baryon feedback, and neutrino mass. KiDS shear measurements are calibrated with a new, more realistic set of image simulations and no significant B-modes are detected in the survey, indicating that systematic errors are under control. When calibrating our redshift distributions by assuming the 30-band COSMOS-2015 photometric redshifts are correct (following the Dark Energy Survey and the Hyper Suprime-Cam Survey), we find the tension with Planck is alleviated. The robust determination of source redshift distributions remains one of the most challenging aspects for future cosmic shear surveys., Comment: 31 pages, 14 figures, accepted for publication by A&A; data products available at http://kids.strw.leidenuniv.nl

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, Comment: 20 pages, 9 figures, 4 tables, submitted to MNRAS

[Context]The data from the Euclid mission will enable the measurement of the angular positions and weak lensing shapes of over a billion galaxies, with their photometric redshifts obtained together with ground-based observations. This large dataset, with well-controlled systematic effects, will allow for cosmological analyses using the angular clustering of galaxies (GCph) and cosmic shear (WL). For Euclid, these two cosmological probes will not be independent because they will probe the same volume of the Universe. The cross-correlation (XC) between these probes can tighten constraints and is therefore important to quantify their impact for Euclid., [Aims] In this study, we therefore extend the recently published Euclid forecasts by carefully quantifying the impact of XC not only on the final parameter constraints for different cosmological models, but also on the nuisance parameters. In particular, we aim to decipher the amount of additional information that XC can provide for parameters encoding systematic effects, such as galaxy bias, intrinsic alignments (IAs), and knowledge of the redshift distributions., [Methods] We follow the Fisher matrix formalism and make use of previously validated codes. We also investigate a different galaxy bias model, which was obtained from the Flagship simulation, and additional photometric-redshift uncertainties; we also elucidate the impact of including the XC terms on constraining these latter., [Results] Starting with a baseline model, we show that the XC terms reduce the uncertainties on galaxy bias by ∼17% and the uncertainties on IA by a factor of about four. The XC terms also help in constraining the γ parameter for minimal modified gravity models. Concerning galaxy bias, we observe that the role of the XC terms on the final parameter constraints is qualitatively the same irrespective of the specific galaxy-bias model used. For IA, we show that the XC terms can help in distinguishing between different models, and that if IA terms are neglected then this can lead to significant biases on the cosmological parameters. Finally, we show that the XC terms can lead to a better determination of the mean of the photometric galaxy distributions., [Conclusions] We find that the XC between GCph and WL within the Euclid survey is necessary to extract the full information content from the data in future analyses. These terms help in better constraining the cosmological model, and also lead to a better understanding of the systematic effects that contaminate these probes. Furthermore, we find that XC significantly helps in constraining the mean of the photometric-redshift distributions, but, at the same time, it requires more precise knowledge of this mean with respect to single probes in order not to degrade the final "figure of merit"., Tutusaus acknowledges support from the Spanish Ministry of Science, Innovation and Universities through grant ESP2017-89838-C3-1-R, and the H2020 programme of the European Commission through grant 776247. M. Martinelli has received the support of a fellowship from “la Caixa” Foundation (ID 100010434), with fellowship code LCF/BQ/PI19/11690015, and the support of the Spanish Agen-cia Estatal de Investigacion through the grant “IFT Centro de Excelencia Severo Ochoa SEV-2016-0597”. S. Camera acknowledges support from the “Departments of Excellence 2018−2022” Grant awarded by the Italian Ministry of Education, University and Research (miur) L. 232/2016. S. Camera is supported by miur through Rita Levi Montalcini project “prometheus – Probing and Relating Observables with Multi-wavelength Experiments To Help Enlightening the Universe’s Structure”. F. Lacasa acknowledges support by the Swiss National Science Foundation. S. Ilić acknowledges financial support from the European Research Council under the European Union’s Seventh Framework Programme (FP7/2007–2013)/ERC Grant Agreement No. 617656 “Theories and Models of the Dark Sector: Dark Matter, Dark Energy and Gravity”. V. Yankelevich acknowledges funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 769130). P. Ntelis is funded by “Centre National d’Études Spatiales” (CNES), for the Euclid project. The Euclid Consortium acknowledges the European Space Agency and the support of a number of agencies and institutes that have supported the development of Euclid, in particular the Academy of Finland, the Agenzia Spaziale Italiana, the Belgian Science Policy, the Canadian Euclid Consortium, the Centre National d’Etudes Spatiales, the Deutsches Zentrum für Luft-und Raumfahrt, the Danish Space Research Institute, the Fundação para a Ciência e a Tecnologia, the Ministerio de Economia y Compet-itividad, the National Aeronautics and Space Administration, the Netherlandse Onderzoekschool Voor Astronomie, the Norwegian Space Agency, the Romanian Space Agency, the State Secretariat for Education, Research and Innovation (SERI) at the Swiss Space Office (SSO), and the United Kingdom Space Agency. A detailed complete list is available on the Euclid website (http:// www.euclid-ec.org).

Alarcon, A., et al., We 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., Argonne National Laboratory’s work was supported by the US Department of Energy, Office of High Energy Physics. Argonne, a U.S. Department of Energy Office of Science Laboratory, is operated by UChicago Argonne LLC under contract no. DE-AC02-06CH11357. The PAU Survey is partially supported by MINECO under grantsCSD2007-00060, AYA2015-71825, ESP2017-89838, PGC2018-094773, PGC2018-102021, SEV-2016-0588, SEV-2016-0597, MDM-2015-0509 and Juan de la Cierva fellowship and LACEGAL and EWC Marie Sklodowska-Curie grant No 734374 and no. 776247 with ERDF funds from the EU Horizon 2020 Programme, some of which include ERDF funds from the European Union. IEEC and IFAE are partially funded by the CERCA and Beatriu de Pinos program of the Generalitat de Catalunya. Funding for PAUS has also been provided by Durham University (via the ERC StG DEGAS-259586), ETH Zurich, Leiden University (via ERC StG ADULT-279396 and Netherlands Organisation for Scientific Research (NWO) Vici grant 639.043.512), University College London and from the European Union’s Horizon 2020 research and innovation programme under the grant agreement No 776247 EWC. The PAU data centre is hosted by the Port d’Informació Científica (PIC), maintained through a collaboration of CIEMAT and IFAE, with additional support from Universitat Autònoma de Barcelona and ERDF. We acknowledge the PIC services department team for their support and fruitful discussions. This project has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 776247. H. Hildebrandt is supported by a Heisenberg grant of the Deutsche Forschungsgemeinschaft (Hi 1495/5-1) as well as an ERC Consolidator Grant (No. 770935). H. Hoekstra acknowledges support from Vici grant 639.043.512 from the Netherlands Organization for Scientific Research (NWO). GM acknowledges support from ST/P006744/1. CMB acknowledges support from ST/P000541/1 and ST/T000244/1. PN acknowledges support from ST/P000541/1 and ST/T000244/1. This project has received funding from the European Union’s Horizon 2020 research and innovation programme under the Maria Skłodowska-Curie (grant agreement No 754510), the National Science Centre of Poland (grant UMO-2016/23/N/ST9/02963) and by the Spanish Ministry of Science and Innovation through Juan de la Cierva-formacion program (reference FJC2018-038792-I).

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., Comment: 17 pages, 12 figures, 3 tables; matches the published version in MNRAS

WL measurements have well-known shear estimation biases, which can be partially corrected for with the use of image simulations. We present an analysis of simulated images that mimic HST/ACS observations of high-redshift galaxy clusters, including cluster specific issues such as non-weak shear and increased blending. Our synthetic galaxies have been generated to match the observed properties of the background-selected samples in the real images. First, we used simulations with galaxies on a grid to determine a revised signal-to-noise-dependent correction for multiplicative shear measurement bias, and to quantify the sensitivity of our bias calibration to mismatches of galaxy or PSF properties between the real data and the simulations. We studied the impact of increased blending and light contamination from cluster and foreground galaxies, finding it negligible for $z>0.7$ clusters, whereas there is an effect at the $\sim 1\%$ level for lower redshift clusters. Finally, we studied the impact of fainter neighbours and selection bias mimicking the positions and magnitudes of galaxies in CANDELS data. The initial SExtractor object detection causes a selection bias of $-0.028 \pm 0.002$, reduced to $-0.016 \pm 0.002$ by further cuts. We compared our CANDELS-based estimate to a grid-based analysis, with added clustered galaxies reaching even fainter magnitudes, yielding a refined estimate of $\sim -0.013$. Our pipeline is calibrated to an accuracy of $\sim 0.015$, which is fully sufficient for current and near-future weak lensing studies of high-redshift clusters. As an application, we used it for a refined analysis of three highly relaxed clusters from the SPT-SZ survey, including measurements down to $r>200$ kpc. Compared to previously employed scales, this tightens the cluster mass constraints by a factor 1.38 on average., 21 pages, 19 figures. Matches the version accepted for publication in A&A

Stage IV weak lensing experiments will offer more than an order of magnitude leap in precision. We must therefore ensure that our analyses remain accurate in this new era. Accordingly, previously ignored systematic effects must be addressed. In this work, we evaluate the impact of the reduced shear approximation and magnification bias, on the information obtained from the angular power spectrum. To first-order, the statistics of reduced shear, a combination of shear and convergence, are taken to be equal to those of shear. However, this approximation can induce a bias in the cosmological parameters that can no longer be neglected. A separate bias arises from the statistics of shear being altered by the preferential selection of galaxies and the dilution of their surface densities, in high-magnification regions. The corrections for these systematic effects take similar forms, allowing them to be treated together. We calculated the impact of neglecting these effects on the cosmological parameters that would be determined from Euclid, using cosmic shear tomography. To do so, we employed the Fisher matrix formalism, and included the impact of the super-sample covariance. We also demonstrate how the reduced shear correction can be calculated using a lognormal field forward modelling approach. These effects cause significant biases in Omega_m, sigma_8, n_s, Omega_DE, w_0, and w_a of -0.53 sigma, 0.43 sigma, -0.34 sigma, 1.36 sigma, -0.68 sigma, and 1.21 sigma, respectively. We then show that these lensing biases interact with another systematic: the intrinsic alignment of galaxies. Accordingly, we develop the formalism for an intrinsic alignment-enhanced lensing bias correction. Applying this to Euclid, we find that the additional terms introduced by this correction are sub-dominant., Comment: 16 pages, 6 figures, submitted to Astronomy & Astrophysics on 16/12/2019, accepted on 04/03/2020. SSC Fisher procedure corrected

We present a new measurement of $E_{\rm G}$, which combines measurements of weak gravitational lensing, real-space galaxy clustering and redshift space distortions. This statistic was proposed as a consistency test of General Relativity (GR) that is insensitive to linear, deterministic galaxy bias and the matter clustering amplitude. We combine deep imaging data from KiDS with overlapping spectroscopy from 2dFLenS, BOSS DR12 and GAMA and find $E_{\rm G}(\overline{z}=0.267)=0.43 \pm 0.13$ (GAMA), $E_{\rm G}(\overline{z}=0.305)=0.27 \pm 0.08$ (LOWZ+2dFLOZ) and $E_{\rm G}(\overline{z}=0.554)=0.26 \pm 0.07$ (CMASS+2dFHIZ). We demonstrate that the existing tension in the value of the matter density parameter hinders the robustness of this statistic as solely a test of GR. We find that our $E_{\rm G}$ measurements, as well as existing ones in the literature, favour a lower matter density cosmology than the Cosmic Microwave Background. For a flat $\Lambda$CDM Universe and assuming GR, we find $\Omega_{\rm m}(z=0)=0.25\pm0.03$. With this paper we publicly release the 2dFLenS dataset at: \url{http://2dflens.swin.edu.au}., Comment: 17 pages, 8 figures

The convolution of galaxy images by the point-spread function (PSF) is the dominant source of bias for weak gravitational lensing studies, and an accurate estimate of the PSF is required to obtain unbiased shape measurements. The PSF estimate for a galaxy depends on its spectral energy distribution (SED), because the instrumental PSF is generally a function of the wavelength. In this paper we explore various approaches to determine the resulting `effective' PSF using broad-band data. Considering the Euclid mission as a reference, we find that standard SED template fitting methods result in biases that depend on source redshift, although this may be remedied if the algorithms can be optimised for this purpose. Using a machine-learning algorithm we show that, at least in principle, the required accuracy can be achieved with the current survey parameters. It is also possible to account for the correlations between photometric redshift and PSF estimates that arise from the use of the same photometry. We explore the impact of errors in photometric calibration, errors in the assumed wavelength dependence of the PSF model and limitations of the adopted template libraries. Our results indicate that the required accuracy for Euclid can be achieved using the data that are planned to determine photometric redshifts.

We present measurements of the weak gravitational lensing shear power spectrum based on $450$ sq. deg. of imaging data from the Kilo Degree Survey. We employ a quadratic estimator in two and three redshift bins and extract band powers of redshift auto-correlation and cross-correlation spectra in the multipole range $76 \leq \ell \leq 1310$. The cosmological interpretation of the measured shear power spectra is performed in a Bayesian framework assuming a $\Lambda$CDM model with spatially flat geometry, while accounting for small residual uncertainties in the shear calibration and redshift distributions as well as marginalising over intrinsic alignments, baryon feedback and an excess-noise power model. Moreover, massive neutrinos are included in the modelling. The cosmological main result is expressed in terms of the parameter combination $S_8 \equiv \sigma_8 \sqrt{\Omega_{\rm m}/0.3}$ yielding $S_8 = \ 0.651 \pm 0.058$ (3 z-bins), confirming the recently reported tension in this parameter with constraints from Planck at $3.2\sigma$ (3 z-bins). We cross-check the results of the 3 z-bin analysis with the weaker constraints from the 2 z-bin analysis and find them to be consistent. The high-level data products of this analysis, such as the band power measurements, covariance matrices, redshift distributions, and likelihood evaluation chains are available at http://kids.strw.leidenuniv.nl/, Comment: 23 pages, 18 figures, 5 tables; results unchanged, version accepted for publication by MNRAS. Data products available at http://kids.strw.leidenuniv.nl/

We measure the cross-correlation between Fermi-LAT gamma-ray photons and over 1000 deg$^2$ of weak lensing data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), the Red Cluster Sequence Lensing Survey (RCSLenS), and the Kilo Degree Survey (KiDS). We present the first measurement of tomographic weak lensing cross-correlations and the first application of spectral binning to cross-correlations between gamma rays and weak lensing. The measurements are performed using an angular power spectrum estimator while the covariance is estimated using an analytical prescription. We verify the accuracy of our covariance estimate by comparing it to two internal covariance estimators. Based on the non-detection of a cross-correlation signal, we derive constraints on weakly interacting massive particle (WIMP) dark matter. We compute exclusion limits on the dark matter annihilation cross-section $\langle\sigma_\rm{ann} v \rangle$, decay rate $\Gamma_\rm{dec}$, and particle mass $m_\rm{DM}$. We find that in the absence of a cross-correlation signal, tomography does not significantly improve the constraining power of the analysis. Assuming a strong contribution to the gamma-ray flux due to small-scale clustering of dark matter and accounting for known astrophysical sources of gamma rays, we exclude the thermal relic cross-section for particle masses of $m_\rm{DM}\lesssim 20$ GeV., Comment: 19 pages, 17 figures, accepted in MNRAS

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}$., 21 pages, 19 figures, accepted for publication in MNRAS

We investigate the disruption of group and cluster satellite galaxies with total mass (dark matter plus baryons) above 10^10 M_sun in the Hydrangea simulations, a suite of 24 high-resolution cosmological hydrodynamical zoom-in simulations based on the EAGLE model. The simulations predict that ~50 per cent of satellites survive to redshift z = 0, with higher survival fractions in massive clusters than in groups and only small differences between baryonic and pure N-body simulations. For clusters, up to 90 per cent of galaxy disruption occurs in lower-mass sub-groups (i.e., during pre-processing); 96 per cent of satellites in massive clusters that were accreted at z < 2 and have not been pre-processed survive. Of those satellites that are disrupted, only a few per cent merge with other satellites, even in low-mass groups. The survival fraction changes rapidly from less than 10 per cent of those accreted at high z to more than 90 per cent at low z. This shift, which reflects faster disruption of satellites accreted at higher z, happens at lower z for more massive galaxies and those accreted onto less massive haloes. The disruption of satellite galaxies is found to correlate only weakly with their pre-accretion baryon content, star formation rate, and size, so that surviving galaxies are nearly unbiased in these properties. These results suggest that satellite disruption in massive haloes is uncommon, and that it is predominantly the result of gravitational rather than baryonic processes., 26 pages, 24 figures (16 pages, 14 figures without appendices); accepted by MNRAS (a few extra references added)

Exploiting the full statistical power of future cosmic shear surveys will necessitate improvements to the accuracy with which the gravitational lensing signal is measured. We present a framework for calibrating shear with image simulations that demonstrates the importance of including realistic correlations between galaxy morphology, size and more importantly, photometric redshifts. This realism is essential so that selection and shape measurement biases can be calibrated accurately for a tomographic cosmic shear analysis. We emulate Kilo-Degree Survey (KiDS) observations of the COSMOS field using morphological information from {\it Hubble} Space Telescope imaging, faithfully reproducing the measured galaxy properties from KiDS observations of the same field. We calibrate our shear measurements from lensfit, and find through a range of sensitivity tests that lensfit is robust and unbiased within the allowed 2 per cent tolerance of our study. Our results show that the calibration has to be performed by selecting the tomographic samples in the simulations, consistent with the actual cosmic shear analysis, because the joint distributions of galaxy properties are found to vary with redshift. Ignoring this redshift variation could result in misestimating the shear bias by an amount that exceeds the allowed tolerance. To improve the calibration for future cosmic shear analyses, it will be essential to also correctly account for the measurement of photometric redshifts, which requires simulating multi-band observations., 31 pages, 17 figures and 2 tables. Accepted for publication in A&A. Matches the published version

In this paper we derive a full expression for the propagation of multiplicative and additive shape measurement biases into the cosmic shear power spectrum. In doing so we identify several new terms that are associated with selection effects, as well as cross-correlation terms between the multiplicative and additive biases and the shear field. The computation of the resulting bias in the shear power spectrum scales as the fifth power of the maximum multipole considered. Consequently the calculation is unfeasible for large l-modes, and the only tractable way to assess the full impact of shape measurement biases on cosmic shear power spectrum is through forward modelling of the effects. To linear order in bias parameters the shear power spectrum is only affected by the mean of the multiplicative bias field over a survey and the cross correlation between the additive bias field and the shear field. If the mean multiplicative bias is zero then second order convolutive terms are expected to be orders of magnitude smaller., 10 pages, accepted to the Open Journal of Astrophysics

In this paper we address the problem of finding optimal cosmic shear tomographic bins. We generalise the definition of a cosmic shear tomographic bin to be a set of commonly labelled voxels in photometric colour space; rather than bins defined directly in redshift. We explore this approach by using a self-organising map to define the multi-dimensional colour space, and a we define a 'label space' of connected regions on the self-organising map using overlapping elliptical disks. This allows us to then find optimal labelling schemes by searching the label space. We use a metric that is the signal-to-noise ratio of a dark energy equation of state measurement, and in this case we find that for up to five tomographic bins the optimal colour-space labelling is an approximation of an equally-spaced binning in redshift; that is in all cases the best configuration. We also show that such a redefinition is more robust to photometric redshift outliers than a standard tomographic bin selection., Comment: 10 pages, submitted to PRD, comments welcome

A number of radio observations have revealed the presence of large synchrotron-emitting sources associated with the intra-cluster medium. There is strong observational evidence that the emitting particles have been (re-)accelerated by shocks and turbulence generated during merger events. The particles that are accelerated are thought to have higher initial energies than those in the thermal pool but the origin of such mildly relativistic particles remains uncertain and needs to be further investigated. The galaxy cluster Abell 1914 is a massive galaxy cluster in which X-ray observations show clear evidence of merging activity. We carried out radio observations of this cluster with the LOw Frequency ARay (LOFAR) at 150 MHz and the Giant Metrewave Radio Telescope (GMRT) at 610 MHz. We also analysed Very Large Array (VLA) 1.4 GHz data, archival GMRT 325 MHz data, CFHT weak lensing data and Chandra observations. Our analysis shows that the ultra-steep spectrum source (4C38.39; $\alpha \lesssim -2$), previously thought to be part of a radio halo, is a distinct source with properties that are consistent with revived fossil plasma sources. Finally, we detect some diffuse emission to the west of the source 4C38.39 that could belong to a radio halo., Comment: 8 figures, 2 tables and 12 pages. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available on lofar-surveys.org upon publication of the journal

The outer regions of galaxies are more susceptible to the tidal interactions that lead to intrinsic alignments of galaxies. The resulting alignment signal may therefore depend on the passband if the colours of galaxies vary spatially. To quantify this, we measured the shapes of galaxies with spectroscopic redshifts from the GAMA survey using deep gri imaging data from the KiloDegree Survey. The performance of the moment-based shape measurement algorithm DEIMOS was assessed using dedicated image simulations, which showed that the ellipticities could be determined with an accuracy better than 1% in all bands. Additional tests for potential systematic errors did not reveal any issues. We measure a significant difference of the alignment signal between the g,r and i-band observations. This difference exceeds the amplitude of the linear alignment model on scales below 2 Mpc/h. Separating the sample into central/satellite and red/blue galaxies, we find that that the difference is dominated by red satellite galaxies., 16 pages, 13 figures, accepted, to appear in A&A