-
Euclid preparation. LensMC, weak lensing cosmic shear measurement with forward modelling and Markov Chain Monte Carlo sampling
Authors:
Euclid Collaboration,
G. Congedo,
L. Miller,
A. N. Taylor,
N. Cross,
C. A. J. Duncan,
T. Kitching,
N. Martinet,
S. Matthew,
T. Schrabback,
M. Tewes,
N. Welikala,
N. Aghanim,
A. Amara,
S. Andreon,
N. Auricchio,
M. Baldi,
S. Bardelli,
R. Bender,
C. Bodendorf,
D. Bonino,
E. Branchini,
M. Brescia,
J. Brinchmann,
S. Camera
, et al. (217 additional authors not shown)
Abstract:
LensMC is a weak lensing shear measurement method developed for Euclid and Stage-IV surveys. It is based on forward modelling to deal with convolution by a point spread function with comparable size to many galaxies; sampling the posterior distribution of galaxy parameters via Markov Chain Monte Carlo; and marginalisation over nuisance parameters for each of the 1.5 billion galaxies observed by Eu…
▽ More
LensMC is a weak lensing shear measurement method developed for Euclid and Stage-IV surveys. It is based on forward modelling to deal with convolution by a point spread function with comparable size to many galaxies; sampling the posterior distribution of galaxy parameters via Markov Chain Monte Carlo; and marginalisation over nuisance parameters for each of the 1.5 billion galaxies observed by Euclid. The scientific performance is quantified through high-fidelity images based on the Euclid Flagship simulations and emulation of the Euclid VIS images; realistic clustering with a mean surface number density of 250 arcmin$^{-2}$ ($I_{\rm E}<29.5$) for galaxies, and 6 arcmin$^{-2}$ ($I_{\rm E}<26$) for stars; and a diffraction-limited chromatic point spread function with a full width at half maximum of $0.^{\!\prime\prime}2$ and spatial variation across the field of view. Objects are measured with a density of 90 arcmin$^{-2}$ ($I_{\rm E}<26.5$) in 4500 deg$^2$. The total shear bias is broken down into measurement (our main focus here) and selection effects (which will be addressed elsewhere). We find: measurement multiplicative and additive biases of $m_1=(-3.6\pm0.2)\times10^{-3}$, $m_2=(-4.3\pm0.2)\times10^{-3}$, $c_1=(-1.78\pm0.03)\times10^{-4}$, $c_2=(0.09\pm0.03)\times10^{-4}$; a large detection bias with a multiplicative component of $1.2\times10^{-2}$ and an additive component of $-3\times10^{-4}$; and a measurement PSF leakage of $α_1=(-9\pm3)\times10^{-4}$ and $α_2=(2\pm3)\times10^{-4}$. When model bias is suppressed, the obtained measurement biases are close to Euclid requirement and largely dominated by undetected faint galaxies ($-5\times10^{-3}$). Although significant, model bias will be straightforward to calibrate given the weak sensitivity. LensMC is publicly available at https://gitlab.com/gcongedo/LensMC
△ Less
Submitted 13 August, 2024; v1 submitted 1 May, 2024;
originally announced May 2024.
-
Power asymmetry in WMAP and Planck temperature sky maps as measured by a local variance estimator
Authors:
Y. Akrami,
Y. Fantaye,
A. Shafieloo,
H. K. Eriksen,
F. K. Hansen,
A. J. Banday,
K. M. Górski
Abstract:
We revisit the question of hemispherical power asymmetry in the WMAP and Planck temperature sky maps by measuring the local variance over the sky and on disks of various sizes. For the 2013 Planck sky map we find that none of the 1000 available isotropic Planck "Full Focal Plane" simulations have a larger variance asymmetry than that estimated from the data, suggesting the presence of an anisotrop…
▽ More
We revisit the question of hemispherical power asymmetry in the WMAP and Planck temperature sky maps by measuring the local variance over the sky and on disks of various sizes. For the 2013 Planck sky map we find that none of the 1000 available isotropic Planck "Full Focal Plane" simulations have a larger variance asymmetry than that estimated from the data, suggesting the presence of an anisotropic signature formally significant at least at the $3.3σ$ level. For the WMAP 9-year data we find that 5 out of 1000 simulations have a larger asymmetry. The preferred direction for the asymmetry from the Planck data is $(l,b)=(212^\circ,-13^\circ)$, in good agreement with previous reports of the same hemispherical power asymmetry.
△ Less
Submitted 24 March, 2014; v1 submitted 4 February, 2014;
originally announced February 2014.
-
Supersymmetry vis-à-vis Observation: Dark Matter Constraints, Global Fits and Statistical Issues
Authors:
Yashar Akrami
Abstract:
Weak-scale supersymmetry is one of the most favoured theories beyond the Standard Model of particle physics that elegantly solves various theoretical and observational problems in both particle physics and cosmology. In this thesis, I describe the theoretical foundations of supersymmetry, issues that it can address and concrete supersymmetric models that are widely used in phenomenological studies…
▽ More
Weak-scale supersymmetry is one of the most favoured theories beyond the Standard Model of particle physics that elegantly solves various theoretical and observational problems in both particle physics and cosmology. In this thesis, I describe the theoretical foundations of supersymmetry, issues that it can address and concrete supersymmetric models that are widely used in phenomenological studies. I discuss how the predictions of supersymmetric models may be compared with observational data from both colliders and cosmology. I show why constraints on supersymmetric parameters by direct and indirect searches of particle dark matter are of particular interest in this respect. Gamma-ray observations of astrophysical sources, in particular dwarf spheroidal galaxies, by the Fermi satellite, and recording nuclear recoil events and energies by future ton-scale direct detection experiments are shown to provide powerful tools in searches for supersymmetric dark matter and estimating supersymmetric parameters. I discuss some major statistical issues in supersymmetric global fits to experimental data. In particular, I further demonstrate that existing advanced scanning techniques may fail in correctly mapping the statistical properties of the parameter spaces even for the simplest supersymmetric models. Complementary scanning methods based on Genetic Algorithms are proposed.
△ Less
Submitted 2 November, 2011;
originally announced November 2011.
-
Statistical coverage for supersymmetric parameter estimation: a case study with direct detection of dark matter
Authors:
Yashar Akrami,
Christopher Savage,
Pat Scott,
Jan Conrad,
Joakim Edsjö
Abstract:
Models of weak-scale supersymmetry offer viable dark matter (DM) candidates. Their parameter spaces are however rather large and complex, such that pinning down the actual parameter values from experimental data can depend strongly on the employed statistical framework and scanning algorithm. In frequentist parameter estimation, a central requirement for properly constructed confidence intervals i…
▽ More
Models of weak-scale supersymmetry offer viable dark matter (DM) candidates. Their parameter spaces are however rather large and complex, such that pinning down the actual parameter values from experimental data can depend strongly on the employed statistical framework and scanning algorithm. In frequentist parameter estimation, a central requirement for properly constructed confidence intervals is that they cover true parameter values, preferably at exactly the stated confidence level when experiments are repeated infinitely many times. Since most widely-used scanning techniques are optimised for Bayesian statistics, one needs to assess their abilities in providing correct confidence intervals in terms of the statistical coverage. Here we investigate this for the Constrained Minimal Supersymmetric Standard Model (CMSSM) when only constrained by data from direct searches for dark matter. We construct confidence intervals from one-dimensional profile likelihoods and study the coverage by generating several pseudo-experiments for a few benchmark sets of pseudo-true parameters. We use nested sampling to scan the parameter space and evaluate the coverage for the benchmarks when either flat or logarithmic priors are imposed on gaugino and scalar mass parameters. The sampling algorithm has been used in the configuration usually adopted for exploration of the Bayesian posterior. We observe both under- and over-coverage, which in some cases vary quite dramatically when benchmarks or priors are modified. We show how most of the variation can be explained as the impact of explicit priors as well as sampling effects, where the latter are indirectly imposed by physicality conditions. For comparison, we also evaluate the coverage for Bayesian credible intervals, and observe significant under-coverage in those cases.
△ Less
Submitted 13 July, 2011; v1 submitted 18 November, 2010;
originally announced November 2010.