No Arabic abstract
Aims. One of the main probes for systematic errors in the cosmic shear signal are the division of the shear field into E- and B-mode shear, where gravitational lensing only produces the former. As shown in a recent note, all currently used E-/B-mode separation methods for the shear correlation functions xi_pm require them to be measured to arbitrarily small and/or large separations which is of course not feasible in practice. Methods. We derive second-order shear statistics which provide a clean separation into E- and B-modes from measurements of xi_pm(theta) over a finite interval only. We call these new statistics the circle and ring statistics, respectively; the latter is obtained by an integral over the former. The mathematical properties of these new shear statistics are obtained, as well as specific expressions for applying them to observed data. Results. It is shown that an E-/B-mode separation can be performed on measurements of xi_pm over a finite interval in angular separation, using the ring statistics. We furthermore generalize this result to derive the most general class of second-order shear statistics which provide a separation of E- and B-mode shear on a given angular interval theta_min <= theta <= theta_max. Our results will be of practical use particularly for future cosmic shear surveys where highly precise measurements of the shear will become available and where control of systematics will be mandatory.
Cosmic shear is considered one of the most powerful methods for studying the properties of Dark Energy in the Universe. As a standard method, the two-point correlation functions $xi_pm(theta)$ of the cosmic shear field are used as statistical measures for the shear field. In order to separate the observed shear into E- and B-modes, the latter being most likely produced by remaining systematics in the data set and/or intrinsic alignment effects, several statistics have been defined before. Here we aim at a complete E-/B-mode decomposition of the cosmic shear information contained in the $xi_pm$ on a finite angular interval. We construct two sets of such E-/B-mode measures, namely Complete Orthogonal Sets of E-/B-mode Integrals (COSEBIs), characterized by weight functions between the $xi_pm$ and the COSEBIs which are polynomials in $theta$ or polynomials in $ln(theta)$, respectively. Considering the likelihood in cosmological parameter space, constructed from the COSEBIs, we study their information contents. We show that the information grows with the number of COSEBI modes taken into account, and that an asymptotic limit is reached which defines the maximum available information in the E-mode component of the $xi_pm$. We show that this limit is reached the earlier (i.e., for a smaller number of modes considered) the narrower the angular range is over which $xi_pm$ are measured, and it is reached much earlier for logarithmic weight functions. For example, for $xi_pm$ on the interval $1le thetale 400$, the asymptotic limit for the parameter pair $(Omega_m, sigma_8)$ is reached for $sim 25$ modes in the linear case, but already for 5 modes in the logarithmic case. The COSEBIs form a natural discrete set of quantities, which we suggest as method of choice in future cosmic shear likelihood analyses.
We present new methods for mapping the curl-free (E-mode) and divergence-free (B-mode) components of spin 2 signals using spin directional wavelets. Our methods are equally applicable to measurements of the polarisation of the cosmic microwave background (CMB) and the shear of galaxy shapes due to weak gravitational lensing. We derive pseudo and pure wavelet estimators, where E-B mixing arising due to incomplete sky coverage is suppressed in wavelet space using scale- and orientation-dependent masking and weighting schemes. In the case of the pure estimator, ambiguous modes (which have vanishing curl and divergence simultaneously on the incomplete sky) are also cancelled. On simulations, we demonstrate the improvement (i.e., reduction in leakage) provided by our wavelet space estimators over standard harmonic space approaches. Our new methods can be directly interfaced in a coherent and computationally-efficient manner with component separation or feature extraction techniques that also exploit wavelets.
Aims. We quantify the mixing of the measured cosmic-shear E- and B-modes caused by the lack of shear-correlation measurements on small and large scales, arising from a lack of close projected galaxy pairs and the finite field size, respectively. Methods. We calculate the aperture-mass statistics <M_{ap, perp}^2> and the E-/B-mode shear-correlation functions xi_{E, B +/-} where small- and large-scale cutoffs are taken into account. We assess the deviation of the obtained E-mode to the true E-mode and the introduction of a spurious B-mode. Results. The measured aperture-mass dispersion is underestimated by more than 10% on scales smaller than 12 times the lower cutoff. For a precise measurement of the E- and B-modes at the percent level using a combination of xi_{E, B +} and xi_{E, B -}, a field as large as 7 (2.4) degrees is necessary for ground-based (space-based) observations.
We explore the stability of the variance and skewness of the cosmic gravitational convergence field, using two different approaches: first we simulate a whole MEGACAM survey (100 sq. degrees). The reconstructed mass map, obtained from a shear map, shows that the state-of-the-art data analysis methods can measure weak-lensing statistics at angular scales ranging from 2.5 to 25. We looked also at the influence of a varying signal-to-noise ratio over the shear map (due to local variations of source density) on the mass reconstruction, by means of Monte-Carlo simulation. The effect at small scales can easily be corrected-for in most of the relevant cases. These results enhance the confidence in the capability of future large surveys to measure accurately cosmologically interesting quantities.
A systematic analysis of the structure of single-baryon correlation functions calculated with lattice QCD is performed, with a particular focus on characterizing the structure of the noise associated with quantum fluctuations. The signal-to-noise problem in these correlation functions is shown, as long suspected, to result from a sign problem. The log-magnitude and complex phase are found to be approximately described by normal and wrapped normal distributions respectively. Properties of circular statistics are used to understand the emergence of a large time noise region where standard energy measurements are unreliable. Power-law tails in the distribution of baryon correlation functions, associated with stable distributions and Levy flights, are found to play a central role in their time evolution. A new method of analyzing correlation functions is considered for which the signal-to-noise ratio of energy measurements is constant, rather than exponentially degrading, with increasing source-sink separation time. This new method includes an additional systematic uncertainty that can be removed by performing an extrapolation, and the signal-to-noise problem re-emerges in the statistics of this extrapolation. It is demonstrated that this new method allows accurate results for the nucleon mass to be extracted from the large-time noise region inaccessible to standard methods. The observations presented here are expected to apply to quantum Monte Carlo calculations more generally. Similar methods to those introduced here may lead to practical improvements in analysis of noisier systems.