Analytical marginalization over photometric redshift uncertainties in cosmic shear analyses
Monthly Notices of the Royal Astronomical Society Oxford University Press 522:4 (2023) 5037-5048
Abstract:
As the statistical power of imaging surveys grows, it is crucial to account for all systematic uncertainties. This is normally done by constructing a model of these uncertainties and then marginalizing over the additional model parameters. The resulting high dimensionality of the total parameter spaces makes inferring the cosmological parameters significantly more costly using traditional Monte Carlo sampling methods. A particularly relevant example is the redshift distribution, p(z ), of the source samples, which may require tens of parameters to describe fully. However, relatively tight priors can be usually placed on these parameters through calibration of the associated systematics. In this paper, we show, quantitatively, that a linearization of the theoretical prediction with respect to these calibrated systematic parameters allows us to analytically marginalize over these extra parameters, leading to a factor of ∼30 reduction in the time needed for parameter inference, while accurately recovering the same posterior distributions for the cosmological parameters that would be obtained through a full numerical marginalization over 160 p(z ) parameters. We demonstrate that this is feasible not only with current data and current achievable calibration priors but also for future Stage-IV data sets.On the functional form of the radial acceleration relation
(2023)
Combining cosmic shear data with correlated photo-z uncertainties: constraints from DESY1 and HSC-DR1
Journal of Cosmology and Astroparticle Physics IOP Publishing (2023)
Abstract:
An accurate calibration of the source redshift distribution p(z) is a key aspect in the analysis of cosmic shear data. This, one way or another, requires the use of spectroscopic or high-quality photometric samples. However, the difficulty to obtain colour-complete spectroscopic samples matching the depth of weak lensing catalogs means that the analyses of different cosmic shear datasets often use the same samples for redshift calibration. This introduces a source of statistical and systematic uncertainty that is highly correlated across different weak lensing datasets, and which must be accurately characterised and propagated in order to obtain robust cosmological constraints from their combination. In this paper we introduce a method to quantify and propagate the uncertainties on the source redshift distribution in two different surveys sharing the same calibrating sample. The method is based on an approximate analytical marginalisation of the p(z) statistical uncertainties and the correlated marginalisation of residual systematics. We apply this method to the combined analysis of cosmic shear data from the DESY1 data release and the HSC-DR1 data, using the COSMOS 30-band catalog as a common redshift calibration sample. We find that, although there is significant correlation in the uncertainties on the redshift distributions of both samples, this does not change the final constraints on cosmological parameters significantly. The same is true also for the impact of residual systematic uncertainties from the errors in the COSMOS 30-band photometric redshifts. Additionally, we show that these effects will still be negligible in Stage-IV datasets. Finally, the combination of DESY1 and HSC-DR1 allows us to constrain the “clumpiness” parameter to S8 = 0.768+0.021 −0.017. This corresponds to a ∼ √ 2 improvement in uncertainties with respect to either DES or HSC alone.Black hole merger simulations in wave dark matter environments
Physical Review D American Physical Society 107:2 (2023) 024035
Abstract:
The interaction of binary black hole mergers with their environments can be studied using numerical relativity simulations. These start only a short finite time before merger, at which point appropriate initial conditions must be imposed. A key task is therefore to identify the configuration that is appropriate for the binary and its environment at this stage of the evolution. In this work we study the behavior of wave dark matter around equal mass black hole binaries, finding that there is a preferred, quasistationary profile that persists and grows over multiple orbits, in contrast to heavier mass dark matter where any overdensity tends to be dispersed by the binary motion. While different initial configurations converge to the preferred quasistationary one after several orbits, unwanted transient oscillations are generated in the process, which may have an impact on the signal in short simulation runs. We also point out that naively superimposing the matter onto a circular binary results in artificially eccentric orbits due to the matter backreaction, which is an effect of the initial conditions and not a signature of dark matter. We discuss the further work required so that comparison of waveforms obtained with environments to vacuum cases can be done in a meaningful way.Where is the ringdown: reconstructing quasinormal modes from dispersive waves
Physical Review D American Physical Society 106 (2022) 104002