Search

We propose that observations of supermassive galaxies contain cosmological statistical constraining power similar to conventional cluster cosmology, and we provide promising indications that the associated systematic errors are comparably easier to control. We consider a fiducial spectroscopic and stellar mass complete sample of galaxies drawn from the Dark Energy Spectroscopic Instrument (DESI) and forecast how constraints on Ωm–σ8 from this sample will compare with those from number counts of clusters based on richness λ. At fixed number density, we find that massive galaxies offer similar constraints to galaxy clusters. However, a mass-complete galaxy sample from DESI has the potential to probe lower halo masses than standard optical cluster samples (which are typically limited to λ ≳ 20 and M halo ≳ 1013.5 M⊙  h −1); additionally, it is straightforward to cleanly measure projected galaxy clustering w p for such a DESI sample, which we show can substantially improve the constraining power on Ωm. We also compare the constraining power of M *-limited samples to those from larger but mass-incomplete samples [e.g. the DESI Bright Galaxy Survey (BGS) sample]; relative to a lower number density M *-limited samples, we find that a BGS-like sample improves statistical constraints by 60 per cent for Ωm and 40 per cent for σ8, but this uses small-scale information that will be harder to model for BGS. Our initial assessment of the systematics associated with supermassive galaxy cosmology yields promising results. The proposed samples have a ∼10 per cent satellite fraction, but we show that cosmological constraints may be robust to the impact of satellites. These findings motivate future work to realize the potential of supermassive galaxies to probe lower halo masses than richness-based clusters and to potentially avoid persistent systematics associated with optical cluster finding. [ABSTRACT FROM AUTHOR]

The overlap of galaxy surveys and cosmic microwave background (CMB) experiments presents an ideal opportunity for joint cosmological data set analyses. In this paper we develop a halo model-based method for the first joint analysis combining these two experiments using 10 correlated two-point functions (10 × 2pt) derived from galaxy position, galaxy shear, CMB lensing convergence, and Compton- y fields. We explore this method using the Vera Rubin Observatory Legacy Survey of Space and Time (LSST) and the Simons Observatory (SO) as examples. We find such LSS × CMB joint analyses lead to significant improvement in Figure-of-Merit of Ωm and S 8 over the constraints from using LSS-only probes within Λ cold dark matter (ΛCDM) model. We identify that the shear– y and y–y correlations are the most valuable additions when thermal Sunyaev–Zel'dolvich (tSZ) is included. We further identify the dominant sources of halo model uncertainties in the small-scale modelling, and investigate the impact of halo self-calibration due to the inclusion of small-scale tSZ information. [ABSTRACT FROM AUTHOR]

Kinematic lensing (KL) is a new cosmological measurement technique that combines traditional weak lensing (WL) shape measurements of disc galaxies with their kinematic information. Using the Tully–Fisher relation, KL breaks the degeneracy between intrinsic and observed ellipticity and significantly reduces the impact of multiple systematics that are present in traditional WL. We explore the performance of KL given the instrument capabilities of the Roman Space Telescope , assuming overlap of the High Latitude Imaging Survey (HLIS) and the High Latitude Spectroscopy Survey (HLSS) over 2000 deg2. Our KL suitable galaxy sample has a number density of n gal = 4 arcmin−1 with an estimated shape noise level of σϵ = 0.035. We quantify the cosmological constraining power on Ωm– S 8 and wp – wa by running simulated likelihood analyses that account for redshift and shear calibration uncertainties, intrinsic alignment, and baryonic feedback. Compared to a traditional WL survey, we find that KL significantly improves the constraining power on Ωm– S 8 (FoMKL = 1.70FoMWL) and wp – wa (FoMKL = 3.65FoMWL). We also explore a 'narrow tomography KL survey' using 30 instead of the default 10 tomographic bins; however, we find no meaningful enhancement to the figure of merit even when assuming a significant time dependence in our fiducial dark energy input scenarios. [ABSTRACT FROM AUTHOR]

We simulate the scientific performance of the Nancy Grace Roman Space Telescope High Latitude Survey (HLS) on dark energy and modified gravity. The 1.6-yr HLS Reference survey is currently envisioned to image 2000 deg2 in multiple bands to a depth of ∼26.5 in Y ,  J ,  H and to cover the same area with slit-less spectroscopy beyond z  = 3. The combination of deep, multiband photometry and deep spectroscopy will allow scientists to measure the growth and geometry of the Universe through a variety of cosmological probes (e.g. weak lensing, galaxy clusters, galaxy clustering, BAO, Type Ia supernova) and, equally, it will allow an exquisite control of observational and astrophysical systematic effects. In this paper, we explore multiprobe strategies that can be implemented, given the telescope's instrument capabilities. We model cosmological probes individually and jointly and account for correlated systematics and statistical uncertainties due to the higher order moments of the density field. We explore different levels of observational systematics for the HLS survey (photo -z and shear calibration) and ultimately run a joint likelihood analysis in N- dim parameter space. We find that the HLS reference survey alone can achieve a standard dark energy FoM of >300 when including all probes. This assumes no information from external data sets, we assume a flat universe however, and includes realistic assumptions for systematics. Our study of the HLS reference survey should be seen as part of a future community-driven effort to simulate and optimize the science return of the Roman Space Telescope. [ABSTRACT FROM AUTHOR]

We explore synergies between the Nancy Grace Roman Space Telescope and the Vera Rubin Observatory's Legacy Survey of Space and Time (LSST). Specifically, we consider scenarios where the currently envisioned survey strategy for the Roman Space Telescope 's High Latitude Survey (HLS reference), i.e. 2000 deg2 in four narrow photometric bands is altered in favour of a strategy of rapid coverage of the LSST area (to full LSST depth) in one band. We find that in only five months, a survey in the W -band can cover the full LSST survey area providing high-resolution imaging for >95 per cent of the LSST Year 10 gold galaxy sample. We explore a second, more ambitious scenario where the Roman Space Telescope spends 1.5 yr covering the LSST area. For this second scenario, we quantify the constraining power on dark energy equation-of-state parameters from a joint weak lensing and galaxy clustering analysis. Our survey simulations are based on the Roman Space Telescope exposure-time calculator and redshift distributions from the CANDELS catalogue. Our statistical uncertainties account for higher order correlations of the density field, and we include a wide range of systematic effects, such as uncertainties in shape and redshift measurements, and modelling uncertainties of astrophysical systematics, such as galaxy bias, intrinsic galaxy alignment, and baryonic physics. We find a significant increase in constraining power for the joint LSST + HLS wide survey compared to LSST Y10 (FoMHLSwide = 2.4 FoMLSST) and compared to LSST + HLS (FoMHLSwide = 5.5 FoMHLSref). [ABSTRACT FROM AUTHOR]

We present a method of combining cluster abundances and large-scale two-point correlations, namely galaxy clustering, galaxy–cluster cross-correlations, cluster autocorrelations, and cluster lensing. This data vector yields comparable cosmological constraints to traditional analyses that rely on small-scale cluster lensing for mass calibration. We use cosmological survey simulations designed to resemble the Dark Energy Survey Year 1 (DES-Y1) data to validate the analytical covariance matrix and the parameter inferences. The posterior distribution from the analysis of simulations is statistically consistent with the absence of systematic biases detectable at the precision of the DES-Y1 experiment. We compare the χ2 values in simulations to their expectation and find no significant difference. The robustness of our results against a variety of systematic effects is verified using a simulated likelihood analysis of DES-Y1-like data vectors. This work presents the first-ever end-to-end validation of a cluster abundance cosmological analysis on galaxy catalogue level simulations. [ABSTRACT FROM AUTHOR]

Accurate covariance matrices for two-point functions are critical for inferring cosmological parameters in likelihood analyses of large-scale structure surveys. Among various approaches to obtaining the covariance, analytic computation is much faster and less noisy than estimation from data or simulations. However, the transform of covariances from Fourier space to real space involves integrals with two Bessel integrals, which are numerically slow and easily affected by numerical uncertainties. Inaccurate covariances may lead to significant errors in the inference of the cosmological parameters. In this paper, we introduce a 2D-FFTLog algorithm for efficient, accurate, and numerically stable computation of non-Gaussian real-space covariances for both 3D and projected statistics. The 2D-FFTLog algorithm is easily extended to perform real-space bin-averaging. We apply the algorithm to the covariances for galaxy clustering and weak lensing for a Dark Energy Survey Year 3-like and a Rubin Observatory's Legacy Survey of Space and Time Year 1-like survey, and demonstrate that for both surveys, our algorithm can produce numerically stable angular bin-averaged covariances with the flat sky approximation, which are sufficiently accurate for inferring cosmological parameters. The code CosmoCov for computing the real-space covariances with or without the flat-sky approximation is released along with this paper. [ABSTRACT FROM AUTHOR]

The tangential shear signal receives contributions from physical scales in the galaxy–matter correlation function well below the transverse scale at which it is measured. Since small scales are difficult to model, this non-locality has generally required stringent scale cuts or new statistics for cosmological analyses. Using the fact that uncertainty in these contributions corresponds to an uncertainty in the enclosed projected mass around the lens, we provide an analytic marginalization scheme to account for this. Our approach enables the inclusion of measurements on smaller scales without requiring numerical sampling over extra free parameters. We extend the analytic marginalization formalism to retain cosmographic (shear-ratio) information from small-scale measurements that would otherwise be removed due to modelling uncertainties, again without requiring the addition of extra sampling parameters. We test the methodology using simulated likelihood analysis of a Dark Energy Survey Year 5-like galaxy–galaxy lensing and galaxy clustering datavector. We demonstrate that we can remove parameter biases due to the presence of an unmodelled one-halo contamination of the galaxy–galaxy lensing signal, and use the shear-ratio information on small scales to improve cosmological parameter constraints. [ABSTRACT FROM AUTHOR]

ABSTRACT We examine sources of scatter in scaling relations between galaxy cluster mass and thermal Sunyaev-Zel'dovich (SZ) effect using cluster samples extracted from cosmological hydrodynamical simulations. Overall, the scatter of the mass-SZ scaling relation is well correlated with the scatter in the mass-halo concentration relation with more concentrated haloes having stronger integrated SZ signals at fixed mass. Additional sources of intrinsic scatter are projection effects from correlated structures, which cause the distribution of scatter to deviate from lognormality and skew it towards higher inferred masses, and the dynamical state of clusters. We study the evolution of merging clusters based on simulations of 39 clusters and their cosmological environment with high time resolution. This sample enables us to study for the first time the detailed evolution of merging clusters around the scaling relation for a cosmologically representative distribution of merger parameters. Major mergers cause an asymmetric scatter such that the inferred mass of merging systems is biased low. We find mergers to be the dominant source of bias towards low inferred masses: over 50 per cent of outliers on this side of the scaling relation underwent a major merger within the last Gyr. As the fraction of dynamically disturbed clusters increases with redshift, our analysis indicates that mergers cause a redshift-dependent bias in scaling relations. Furthermore, we find the SZ morphology of massive clusters to be well correlated with the clusters' dynamical state, suggesting that morphology may be used to constrain merger fractions and identify merger-induced outliers of the scaling relation. [ABSTRACT FROM AUTHOR]