Weak lensing, which is the deflection of light by matter along the line of sight, has proven to be an efficient method for constraining models of structure formation and reveal the nature of dark energy. So far, most weak-lensing studies have focused on the shear field that can be measured directly from the ellipticity of background galaxies. However, within the context of forthcoming full-sky weak-lensing surveys such as Euclid, convergence maps (mass maps) offer an important advantage over shear fields in terms of cosmological exploitation. While it carry the same information, the lensing signal is more compressed in the convergence maps than in the shear field. This simplifies otherwise computationally expensive analyses, for instance, non-Gaussianity studies. However, the inversion of the non-local shear field requires accurate control of systematic effects caused by holes in the data field, field borders, shape noise, and the fact that the shear is not a direct observable (reduced shear). We present the two mass-inversion methods that are included in the official Euclid data-processing pipeline: the standard Kaiser & Squires method (KS), and a new mass-inversion method (KS+) that aims to reduce the information loss during the mass inversion. This new method is based on the KS method and includes corrections for mass-mapping systematic effects. The results of the KS+ method are compared to the original implementation of the KS method in its simplest form, using the Euclid Flagship mock galaxy catalogue. In particular, we estimate the quality of the reconstruction by comparing the two-point correlation functions and third- and fourth- order moments obtained from shear and convergence maps, and we analyse each systematic effect independently and simultaneously...

The data from the Euclid mission will enable the measurement of the photometric redshifts, angular positions, and weak lensing shapes for over a billion galaxies. This large dataset will allow for cosmological analyses using the angular clustering of galaxies and cosmic shear. The cross-correlation (XC) between these probes can tighten constraints and it is therefore important to quantify their impact for Euclid. In this study we carefully quantify 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 at understanding the amount of additional information that XC can provide for parameters encoding systematic effects, such as galaxy bias or intrinsic alignments (IA). We follow the formalism presented in Euclid Collaboration: Blanchard et al. (2019) and make use of the codes validated therein. We show that XC improves the dark energy Figure of Merit (FoM) by a factor $\sim 5$, whilst it also reduces the uncertainties on galaxy bias by $\sim 17\%$ and the uncertainties on IA by a factor $\sim 4$. We observe that the role of XC on the final parameter constraints is qualitatively the same irrespective of the galaxy bias model used. We also show that XC can help in distinguishing between different IA models, and that if IA terms are neglected then this can lead to significant biases on the cosmological parameters. We find that the XC terms are necessary to extract the full information content from the data in future analyses. They help in better constraining the cosmological model, and lead to a better understanding of the systematic effects that contaminate these probes. Furthermore, we find that XC helps in constraining the mean of the photometric-redshift distributions, but it requires a more precise knowledge of this mean in order not to degrade the final FoM. [Abridged]

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.

The Euclid mission will observe well over a billion galaxies out to $z\sim6$ and beyond. This will offer an unrivalled opportunity to investigate several key questions for understanding galaxy formation and evolution. The first step for many of these studies will be the selection of a sample of quiescent and star-forming galaxies, as is often done in the literature by using well known colour techniques such as the `UVJ' diagram. However, given the limited number of filters available for the Euclid telescope, the recovery of such rest-frame colours will be challenging. We therefore investigate the use of observed Euclid colours, on their own and together with ground-based u-band observations, for selecting quiescent and star-forming galaxies. The most efficient colour combination, among the ones tested in this work, consists of the (u-VIS) and (VIS-J) colours. We find that this combination allows users to select a sample of quiescent galaxies complete to above $\sim70\%$ and with less than 15$\%$ contamination at redshifts in the range $0.75<z<1$. For galaxies at high-z or without the u-band complementary observations, the (VIS-Y) and (J-H) colours represent a valid alternative, with $>65\%$ completeness level and contamination below 20$\%$ at $1<z<2$ for finding quiescent galaxies. In comparison, the sample of quiescent galaxies selected with the traditional UVJ technique is only $\sim20\%$ complete at $z<3$, when recovering the rest-frame colours using mock Euclid observations. This shows that our new methodology is the most suitable one when only Euclid bands, along with u-band imaging, are available.

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

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. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate 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 that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such 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.

The Weak Equivalence Principle (WEP), stating that two bodies of different compositions and/or mass fall at the same rate in a gravitational field (universality of free fall), is at the very foundation of General Relativity. The MICROSCOPE mission aims to test its validity to a precision of $10^{-15}$, two orders of magnitude better than current on-ground tests, by using two masses of different compositions (titanium and platinum alloys) on a quasi-circular trajectory around the Earth. This is realised by measuring the accelerations inferred from the forces required to maintain the two masses exactly in the same orbit. Any significant difference between the measured accelerations, occurring at a defined frequency, would correspond to the detection of a violation of the WEP, or to the discovery of a tiny new type of force added to gravity. MICROSCOPE's first results show no hint for such a difference, expressed in terms of E\"otv\"os parameter $\delta(Ti,Pt)=[-1\pm{}9{\rm (stat)}\pm{}9{\rm (syst)}] \times{}10^{-15}$ (both 1$\sigma$ uncertainties) for a titanium and platinum pair of materials. This result was obtained on a session with 120 orbital revolutions representing 7\% of the current available data acquired during the whole mission. The quadratic combination of 1$\sigma$ uncertainties leads to a current limit on $\delta$ of about $1.3\times{}10^{-14}$.

The Canada-France Imaging Survey (CFIS) will map the northern high Galactic latitude sky in the $u$-band ("CFIS-u", 10,000$\, {\rm deg^2}$) and in the $r$-band ("CFIS-r", 5,000$\, {\rm deg^2}$), enabling a host of stand-alone science investigations, and providing some of the ground-based data necessary for photometric redshift determination for the Euclid mission. In this first contribution we present the $u$-band component of the survey, describe the observational strategy, and discuss some first highlight results, based on approximately one third of the final area. We show that the Galactic anticenter structure is distributed continuously along the line of sight, out to beyond 20 kpc, and possesses a metallicity distribution that is essentially identical to that of the outer disk sampled by APOGEE. This suggests that it is probably a buckled disk of old metal-rich stars, rather than a stream or a flare. We also discuss the future potential for CFIS-u in discovering star-forming dwarf galaxies around the Local Group, the characterization of the white dwarf and blue straggler population of the Milky Way, as well as its sensitivity to low-surface brightness structures in external galaxies.

We present the chemical distribution of the Milky Way, based on 2,900$\, {\rm deg^2}$ of $u$-band photometry taken as part of the Canada-France Imaging Survey. When complete, this survey will cover 10,000$\, {\rm deg^2}$ of the Northern sky. By combing the CFHT $u$-band photometry together with SDSS and Pan-STARRS $g,r,$ and $i$, we demonstrate that we are able to measure reliably the metallicities of individual stars to $\sim 0.2$ dex, and hence additionally obtain good photometric distance estimates. This survey thus permits the measurement of metallicities and distances of the dominant main-sequence population out to approximately 30 kpc, and provides much higher number of stars at large extraplanar distances than have been available from previous surveys. We develop a non-parametric distance-metallicity decomposition algorithm and apply it to the sky at $30\deg < |b| < 70\deg$ and to the North Galactic Cap. We find that the metallicity-distance distribution is well-represented by three populations whose metallicity distributions do not vary significantly with vertical height above the disk. As traced in main-sequence stars, the stellar halo component shows a vertical density profile that is close to exponential, with a scale height of around 3 kpc. This may indicate that the inner halo was formed partly from disk stars ejected in an ancient minor merger.

The MICROSCOPE space mission, launched on April 25, 2016, aims to test the weak equivalence principle (WEP) with a 10^-15 precision. To reach this performance requires an accurate and robust data analysis method, especially since the possible WEP violation signal will be dominated by a strongly colored noise. An important complication is brought by the fact that some values will be missing -therefore, the measured time series will not be strictly regularly sampled. Those missing values induce a spectral leakage that significantly increases the noise in Fourier space, where the WEP violation signal is looked for, thereby complicating scientific returns. Recently, we developed an inpainting algorithm to correct the MICROSCOPE data for missing values. This code has been integrated in the official MICROSCOPE data processing pipeline because it enables us to significantly measure an equivalence principle violation (EPV) signal in a model-independent way, in the inertial satellite configuration. In this work, we present several improvements to the method that may allow us now to reach the MICROSCOPE requirements for both inertial and spin satellite configurations. The main improvement has been obtained using a prior on the power spectrum of the colored-noise that can be directly derived from the incomplete data. We show that after reconstructing missing values with this new algorithm, a least-squares fit may allow us to significantly measure an EPV signal with a 0.96x10^-15 precision in the inertial mode and 1.2x10^-15 precision in the spin mode. Although, the inpainting method presented in this paper has been optimized to the MICROSCOPE data, it remains sufficiently general to be used in the general context of missing data in time series dominated by an unknown colored-noise. The improved inpainting software, called ICON, is freely available at http://www.cosmostat.org/software/icon.

This is the third in a series of papers that develop a new and flexible model to predict weak-lensing (WL) peak counts, which have been shown to be a very valuable non-Gaussian probe of cosmology. In this paper, we compare the cosmological information extracted from WL peak counts using different filtering techniques of the galaxy shear data, including linear filtering with a Gaussian and two compensated filters (the starlet wavelet and the aperture mass), and the nonlinear filtering method MRLens. We present improvements to our model that account for realistic survey conditions, which are masks, shear-to-convergence transformations, and non-constant noise. We create simulated peak counts from our stochastic model, from which we obtain constraints on the matter density $\Omega_\mathrm{m}$, the power spectrum normalisation $\sigma_8$, and the dark-energy parameter $w_0$. We use two methods for parameter inference, a copula likelihood, and approximate Bayesian computation (ABC). We measure the contour width in the $\Omega_\mathrm{m}$-$\sigma_8$ degeneracy direction and the figure of merit to compare parameter constraints from different filtering techniques. We find that starlet filtering outperforms the Gaussian kernel, and that including peak counts from different smoothing scales helps to lift parameter degeneracies. Peak counts from different smoothing scales with a compensated filter show very little cross-correlation, and adding information from different scales can therefore strongly enhance the available information. Measuring peak counts separately from different scales yields tighter constraints than using a combined peak histogram from a single map that includes multiscale information. Our results suggest that a compensated filter function with counts included separately from different smoothing scales yields the tightest constraints on cosmological parameters from WL peaks.

We propose a new mass-mapping algorithm, specifically designed to recover small-scale information from a combination of gravitational shear and flexion. Including flexion allows us to supplement the shear on small scales in order to increase the sensitivity to substructures and the overall resolution of the convergence map without relying on strong lensing constraints. In order to preserve all available small scale information, we avoid any binning of the irregularly sampled input shear and flexion fields and treat the mass-mapping problem as a general ill-posed inverse problem, regularised using a robust multi-scale wavelet sparsity prior. The resulting algorithm incorporates redshift, reduced shear, and reduced flexion measurements for individual galaxies and is made highly efficient by the use of fast Fourier estimators. We test our reconstruction method on a set of realistic weak lensing simulations corresponding to typical HST/ACS cluster observations and demonstrate our ability to recover substructures with the inclusion of flexion which are lost if only shear information is used. In particular, we can detect substructures at the 15$^{\prime \prime}$ scale well outside of the critical region of the clusters. In addition, flexion also helps to constrain the shape of the central regions of the main dark matter halos. Our mass-mapping software, called Glimpse2D, is made freely available at http://www.cosmostat.org/software/glimpse .

Missing data are a common problem in experimental and observational physics. They can be caused by various sources, either an instrument's saturation, or a contamination from an external event, or a data loss. In particular, they can have a disastrous effect when one is seeking to characterize a colored-noise-dominated signal in Fourier space, since they create a spectral leakage that can artificially increase the noise. It is therefore important to either take them into account or to correct for them prior to e.g. a Least-Square fit of the signal to be characterized. In this paper, we present an application of the {\it inpainting} algorithm to mock MICROSCOPE data; {\it inpainting} is based on a sparsity assumption, and has already been used in various astrophysical contexts; MICROSCOPE is a French Space Agency mission, whose launch is expected in 2016, that aims to test the Weak Equivalence Principle down to the $10^{-15}$ level. We then explore the {\it inpainting} dependence on the number of gaps and the total fraction of missing values. We show that, in a worst-case scenario, after reconstructing missing values with {\it inpainting}, a Least-Square fit may allow us to significantly measure a $1.1\times10^{-15}$ Equivalence Principle violation signal, which is sufficiently close to the MICROSCOPE requirements to implement {\it inpainting} in the official MICROSCOPE data processing and analysis pipeline. Together with the previously published KARMA method, {\it inpainting} will then allow us to independently characterize and cross-check an Equivalence Principle violation signal detection down to the $10^{-15}$ level.

In asteroseismology, the observed time series often suffers from incomplete time coverage due to gaps. The presence of periodic gaps may generate spurious peaks in the power spectrum that limit the analysis of the data. Various methods have been developed to deal with gaps in time series data. However, it is still important to improve these methods to be able to extract all the possible information contained in the data. In this paper, we propose a new approach to handle the problem, the so-called inpainting method. This technique, based on a sparsity prior, enables to judiciously fill-in the gaps in the data, preserving the asteroseismic signal, as far as possible. The impact of the observational window function is reduced and the interpretation of the power spectrum is simplified. This method is applied both on ground and space-based data. It appears that the inpainting technique improves the oscillation modes detection and estimation. Additionally, it can be used to study very long time series of many stars because its computation is very fast. For a time series of 50 days of CoRoT-like data, it allows a speed-up factor of 1000, if compared to methods of the same accuracy.

The NASA Kepler mission has observed more than 190,000 stars in the constellations of Cygnus and Lyra. Around 4 years of almost continuous ultra high-precision photometry have been obtained reaching a duty cycle higher than 90% for many of these stars. However, almost regular gaps due to nominal operations are present in the light curves at different time scales. In this paper we want to highlight the impact of those regular gaps in asteroseismic analyses and we try to find a method that minimizes their effect in the frequency domain. To do so, we isolate the two main time scales of quasi regular gaps in the data. We then interpolate the gaps and we compare the power density spectra of four different stars: two red giants at different stages of their evolution, a young F-type star, and a classical pulsator in the instability strip. The spectra obtained after filling the gaps in the selected solar-like stars show a net reduction in the overall background level, as well as a change in the background parameters. The inferred convective properties could change as much as 200% in the selected example, introducing a bias in the p-mode frequency of maximum power. When global asteroseismic scaling relations are used, this bias can lead up to a variation in the surface gravity of 0.05 dex. Finally, the oscillation spectrum in the classical pulsator is cleaner compared to the original one.

We present an algorithm using Principal Component Analysis (PCA) to subtract galaxies from imaging data, and also two algorithms to find strong, galaxy-scale gravitational lenses in the resulting residual image. The combined method is optimized to find full or partial Einstein rings. Starting from a pre-selection of potential massive galaxies, we first perform a PCA to build a set of basis vectors. The galaxy images are reconstructed using the PCA basis and subtracted from the data. We then filter the residual image with two different methods. The first uses a curvelet (curved wavelets) filter of the residual images to enhance any curved/ring feature. The resulting image is transformed in polar coordinates, centered on the lens galaxy center. In these coordinates, a ring is turned into a line, allowing us to detect very faint rings by taking advantage of the integrated signal-to-noise in the ring (a line in polar coordinates). The second way of analysing the PCA-subtracted images identifies structures in the residual images and assesses whether they are lensed images according to their orientation, multiplicity and elongation. We apply the two methods to a sample of simulated Einstein rings, as they would be observed with the ESA Euclid satellite in the VIS band. The polar coordinates transform allows us to reach a completeness of 90% and a purity of 86%, as soon as the signal-to-noise integrated in the ring is higher than 30, and almost independent of the size of the Einstein ring. Finally, we show with real data that our PCA-based galaxy subtraction scheme performs better than traditional subtraction based on model fitting to the data. Our algorithm can be developed and improved further using machine learning and dictionary learning methods, which would extend the capabilities of the method to more complex and diverse galaxy shapes.

In this paper, we compare three methods to reconstruct galaxy cluster density fields with weak lensing data. The first method called FLens integrates an inpainting concept to invert the shear field with possible gaps, and a multi-scale entropy denoising procedure to remove the noise contained in the final reconstruction, that arises mostly from the random intrinsic shape of the galaxies. The second and third methods are based on a model of the density field made of a multi-scale grid of radial basis functions. In one case, the model parameters are computed with a linear inversion involving a singular value decomposition. In the other case, the model parameters are estimated using a Bayesian MCMC optimization implemented in the lensing software Lenstool. Methods are compared on simulated data with varying galaxy density fields. We pay particular attention to the errors estimated with resampling. We find the multi-scale grid model optimized with MCMC to provide the best results, but at high computational cost, especially when considering resampling. The SVD method is much faster but yields noisy maps, although this can be mitigated with resampling. The FLens method is a good compromise with fast computation, high signal to noise reconstruction, but lower resolution maps. All three methods are applied to the MACS J0717+3745 galaxy cluster field, and reveal the filamentary structure discovered in Jauzac et al. 2012. We conclude that sensitive priors can help to get high signal to noise, and unbiased reconstructions.

The next generation of galaxy surveys, aiming to observe millions of galaxies, are expensive both in time and cost. This raises questions regarding the optimal investment of this time and money for future surveys. In a previous work, it was shown that a sparse sampling strategy could be a powerful substitute for the contiguous observations. However, in this previous paper a regular sparse sampling was investigated, where the sparse observed patches were regularly distributed on the sky. The regularity of the mask introduces a periodic pattern in the window function, which induces periodic correlations at specific scales. In this paper, we use the Bayesian experimental design to investigate a random sparse sampling, where the observed patches are randomly distributed over the total sparsely sampled area. We find that, as there is no preferred scale in the window function, the induced correlation is evenly distributed amongst all scales. This could be desirable if we are interested in specific scales in the galaxy power spectrum, such as the Baryonic Acoustic Oscillation (BAO) scales. However, for constraining the overall galaxy power spectrum and the cosmological parameters, there is no preference over regular or random sampling. Hence any approach that is practically more suitable can be chosen and we can relax the regular-grid condition for the distribution of the observed patches.

This paper describes the definition of a typical next-generation space-based weak gravitational lensing experiment. We first adopt a set of top-level science requirements from the literature, based on the scale and depth of the galaxy sample, and the avoidance of systematic effects in the measurements which would bias the derived shear values. We then identify and categorise the contributing factors to the systematic effects, combining them with the correct weighting, in such a way as to fit within the top-level requirements. We present techniques which permit the performance to be evaluated and explore the limits at which the contributing factors can be managed. Besides the modelling biases resulting from the use of weighted moments, the main contributing factors are the reconstruction of the instrument point spread function (PSF), which is derived from the stellar images on the image, and the correction of the charge transfer inefficiency (CTI) in the CCD detectors caused by radiation damage.

The first half of this paper explores the origin of systematic biases in the measurement of weak gravitational lensing. Compared to previous work, we expand the investigation of PSF instability and fold in for the first time the effects of non-idealities in electronic imaging detectors and imperfect galaxy shape measurement algorithms. Together, these now explain the additive A(l) and multiplicative M(l) systematics typically reported in current lensing measurements. We find that overall performance is driven by a product of a telescope/camera's *absolute performance*, and our *knowledge about its performance*. The second half of this paper propagates any residual shear measurement biases through to their effect on cosmological parameter constraints. Fully exploiting the statistical power of Stage IV weak lensing surveys will require additive biases A<1.8e-12 and multiplicative biases M<4.0e-3. These can be allocated between individual budgets in hardware, calibration data and software, using results from the first half of the paper. If instrumentation is stable and well-calibrated, we find extant shear measurement software from GREAT10 already meet requirements on galaxies detected at S/N=40. Averaging over a population of galaxies with a realistic distribution of sizes, it also meets requirements for a 2D cosmic shear analysis from space. If used on fainter galaxies or for 3D cosmic shear tomography, existing algorithms would need calibration on simulations to avoid introducing bias at a level similar to the statistical error. Requirements on hardware and calibration data are discussed in more detail in a companion paper. Our analysis is intentionally general, but is specifically being used to drive the hardware and ground segment performance budget for the design of the European Space Agency's recently-selected Euclid mission.

The aperture mass statistic is a common tool used in weak lensing studies. By convolving lensing maps with a filter function of a specific scale, chosen to be larger than the scale on which the noise is dominant, the lensing signal may be boosted with respect to the noise. This allows for detection of structures at increased fidelity. Furthermore, higher-order statistics of the aperture mass (such as its skewness or kurtosis), or counting of the peaks seen in the resulting aperture mass maps, provide a convenient and effective method to constrain the cosmological parameters. In this paper, we more fully explore the formalism underlying the aperture mass statistic. We demonstrate that the aperture mass statistic is formally identical to a wavelet transform at a specific scale. Further, we show that the filter functions most frequently used in aperture mass studies are not ideal, being non-local in both real and Fourier space. In contrast, the wavelet formalism offers a number of wavelet functions that are localized both in real and Fourier space, yet similar to the 'optimal' aperture mass filters commonly adopted. Additionally, for a number of wavelet functions, such as the starlet wavelet, very fast algorithms exist to compute the wavelet transform. This offers significant advantages over the usual aperture mass algorithm when it comes to image processing time, demonstrating speed-up factors of ~ 5 - 1200 for aperture radii in the range 2 to 64 pixels on an image of 1024 x 1024 pixels.

Weak gravitational lensing has become a common tool to constrain the cosmological model. The majority of the methods to derive constraints on cosmological parameters use second-order statistics of the cosmic shear. Despite their success, second-order statistics are not optimal and degeneracies between some parameters remain. Tighter constraints can be obtained if second-order statistics are combined with a statistic that is efficient to capture non-Gaussianity. In this paper, we search for such a statistical tool and we show that there is additional information to be extracted from statistical analysis of the convergence maps beyond what can be obtained from statistical analysis of the shear field. For this purpose, we have carried out a large number of cosmological simulations along the {\sigma}8-{\Omega}m degeneracy, and we have considered three different statistics commonly used for non-Gaussian features characterization: skewness, kurtosis and peak count. To be able to investigate non-Gaussianity directly in the shear field we have used the aperture mass definition of these three statistics for different scales. Then, the results have been compared with the results obtained with the same statistics estimated in the convergence maps at the same scales. First, we show that shear statistics give similar constraints to those given by convergence statistics, if the same scale is considered. In addition, we find that the peak count statistic is the best to capture non-Gaussianities in the weak lensing field and to break the {\sigma}8-{\Omega}m degeneracy. We show that this statistical analysis should be conducted in the convergence maps: first, because there exist fast algorithms to compute the convergence map for different scales, and secondly because it offers the opportunity to denoise the reconstructed convergence map, which improves non-Gaussian features extraction.

To derive the convergence field from the gravitational shear (gamma) of the background galaxy images, the classical methods require a convolution of the shear to be performed over the entire sky, usually expressed thanks to the Fast Fourier transform (FFT). However, it is not optimal for an imperfect geometry survey. Furthermore, FFT implicitly uses periodic conditions that introduce errors to the reconstruction. A method has been proposed that relies on computation of an intermediate field u that combines the derivatives of gamma and on convolution with a Green kernel. In this paper, we study the wavelet Helmholtz decomposition as a new approach to reconstructing the dark matter mass map. We show that a link exists between the Helmholtz decomposition and the E/B mode separation. We introduce a new wavelet construction, that has a property that gives us more flexibility in handling the border problem, and we propose a new method of reconstructing the dark matter mass map in the wavelet space. A set of experiments based on noise-free images illustrates that this Wavelet Helmholtz decomposition reconstructs the borders better than all other existing methods.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

The measurements of the statistical properties of the Cosmic Microwave Background (CMB) fluctuations enable us to probe the physics of the very early Universe especially at the epoch of inflation. A particular interest lays on the detection of the non-Gaussianity of the CMB as it can constrain the current proposed models of inflation and structure formation, or possibly point out new models. The current approach to measure the degree of non-Gaussianity of the CMB is to estimate a single parameter which is highly model-dependent. The bispectrum is a natural and widely studied tool for measuring the non-Gaussianity in a model-independent way. This paper sets the grounds for a full CMB bispectrum estimator based on the decomposition of the sphere onto projected patches. The mean bispectrum estimated this way can be calculated quickly and is model-independent. This approach is very flexible, allowing exclusion of some patches in the processing or consideration of just a specific region of the sphere.

We present high spatial resolution optical imaging and polarization observations of the PSR B0540-69.3 and its highly dynamical pulsar wind nebula (PWN) performed with HST, and compare them with X-ray data obtained with the Chandra X-ray Observatory. We have studied the bright region southwest of the pulsar where a bright "blob" is seen in 1999. We show that it may be a result of local energy deposition around 1999, and that the emission from this then faded away. Polarization data from 2007 show that the polarization properties show dramatic spatial variations at the 1999 blob position arguing for a local process. Several other positions along the pulsar-"blob" orientation show similar changes in polarization, indicating previous recent local energy depositions. In X-rays, the spectrum steepens away from the "blob" position, faster orthogonal to the pulsar-"blob" direction than along this axis of orientation. This could indicate that the pulsar-"blob" orientation is an axis along where energy in the PWN is mainly injected, and that this is then mediated to the filaments in the PWN by shocks. We highlight this by constructing an [S II]-to-[O III]-ratio map. We argue, through modeling, that the high [S II]/[O III] ratio is not due to time-dependent photoionization caused by possible rapid Xray emission variations in the "blob" region. We have also created a multiwavelength energy spectrum for the "blob" position showing that one can, to within 2sigma, connect the optical and X-ray emission by a single power law. We obtain best power-law fits for the X-ray spectrum if we include "extra" oxygen, in addition to the oxygen column density in the interstellar gas of the Large Magellanic Cloud and the Milky Way. This oxygen is most naturally explained by the oxygen-rich ejecta of the supernova remnant. The oxygen needed likely places the progenitor mass in the 20 - 25 Msun range.

Weak gravitational lensing has proven to be a powerful tool to map directly the distribution of dark matter in the Universe. The technique, currently used, relies on the accurate measurement of the gravitational shear that corresponds to the first-order distortion of the background galaxy images. More recently, a new technique has been introduced that relies on the accurate measurement of the gravitational flexion that corresponds to the second-order distortion of the background galaxy images. This technique should probe structures on smaller scales than that of a shear analysis. The goal of this paper is to compare the ability of shear and flexion to reconstruct the dark matter distribution by taking into account the dispersion in shear and flexion measurements. Our results show that the flexion is less sensitive than shear for constructing the convergence maps on scales that are physically feasible for mapping, meaning that flexion alone not be used to do convergence map reconstruction, even on small scales.

In Helio- and asteroseismology, it is important to have continuous, uninterrupted, data sets. However, seismic observations usually contain gaps and we need to take them into account. In particular, if the gaps are not randomly distributed, they will produce a peak and a series of harmonics in the periodogram that will destroy the stellar information. An interpolation of the data can be good solution for this problem. In this paper we have studied an interpolation method based on the so-called 'inpainting' algorithms. To check the algorithm, we used both VIRGO and CoRoT satellite data to which we applied a realistic artificial window of a real CoRoT observing run to introduce gaps. Next we compared the results with the original, non-windowed data. Therefore, we were able to optimize the algorithm by minimizing the difference between the power spectrum density of the data with gaps and the complete time series. In general, we find that the power spectrum of the inpainted time series is very similar to the original, unperturbed one. Seismic inferences obtained after interpolating the data are the same as in the original case.

This paper reviews statistical methods recently developed to reconstruct and analyze dark matter mass maps from weak lensing observations. The field of weak lensing is motivated by the observations made in the last decades showing that the visible matter represents only about 4-5% of the Universe, the rest being dark. The Universe is now thought to be mostly composed by an invisible, pressureless matter -potentially relic from higher energy theories- called "dark matter" (20-21%) and by an even more mysterious term, described in Einstein equations as a vacuum energy density, called "dark energy" (70%). This "dark" Universe is not well described or even understood, so this point could be the next breakthrough in cosmology. Weak gravitational lensing is believed to be the most promising tool to understand the nature of dark matter and to constrain the cosmological model used to describe the Universe. Gravitational lensing is the process in which light from distant galaxies is bent by the gravity of intervening mass in the Universe as it travels towards us. This bending causes the image of background galaxies to appear slightly distorted and can be used to extract significant results for cosmology. Future weak lensing surveys are already planned in order to cover a large fraction of the sky with large accuracy. However this increased accuracy also places greater demands on the methods used to extract the available information. In this paper, we will first describe the important steps of the weak lensing processing to reconstruct the dark matter distribution from shear estimation. Then we will discuss the problem of statistical estimation in order to set constraints on the cosmological model. We review the methods which are currently used especially new methods based on sparsity.

The GRavitational lEnsing Accuracy Testing 2008 (GREAT08) Challenge focuses on a problem that is of crucial importance for future observations in cosmology. The shapes of distant galaxies can be used to determine the properties of dark energy and the nature of gravity, because light from those galaxies is bent by gravity from the intervening dark matter. The observed galaxy images appear distorted, although only slightly, and their shapes must be precisely disentangled from the effects of pixelisation, convolution and noise. The worldwide gravitational lensing community has made significant progress in techniques to measure these distortions via the Shear TEsting Program (STEP). Via STEP, we have run challenges within our own community, and come to recognise that this particular image analysis problem is ideally matched to experts in statistical inference, inverse problems and computational learning. Thus, in order to continue the progress seen in recent years, we are seeking an infusion of new ideas from these communities. This document details the GREAT08 Challenge for potential participants. Please visit http://www.great08challenge.info for the latest information.

Weak gravitational lensing provides a unique method to map directly the dark matter in the Universe. The majority of lensing analyses uses the two-point statistics of the cosmic shear field to constrain the cosmological model yielding degeneracies, such as that between sigma_8 and Omega_M respectively the r.m.s. of the mass fluctuations at a scale of 8 Mpc/h and the matter density parameter both at z = 0. However, the two-point statistics only measure the Gaussian properties of the field and the weak lensing field is non-Gaussian. It has been shown that the estimation of non-Gaussian statistics on weak lensing data can improve the constraints on cosmological parameters. In this paper, we systematically compare a wide range of non-Gaussian estimators in order to determine which one provides tighter constraints on the cosmological parameters. These statistical methods include skewness, kurtosis and the Higher Criticism test in several sparse representations such as wavelet and curvelet; as well as the bispectrum, peak counting and a new introduced statistic called Wavelet Peak Counting (WPC). Comparisons based on sparse representations show that the wavelet transform is the most sensitive to non-Gaussian cosmological structures. It appears also that the best statistic for non-Gaussian characterization in weak lensing mass maps is the WPC. Finally, we show that the sigma_8 -Omega_m degeneracy could be even better broken if the WPC estimation is performed on weak lensing mass maps filtered by the wavelet method, MRLens.

We study the correlation between the locations of galaxy-galaxy strong lensing candidates and tracers of large-scale structure from both weak lensing or X-ray emission. The COSMOS survey is a unique data set, combining deep, high resolution and contiguous imaging in which strong lenses have been discovered, plus unparalleled multiwavelength coverage. To help interpret the COSMOS data, we have also produced mock COSMOS strong and weak lensing observations, based on ray-tracing through the Millenium simulation. In agreement with the simulations, we find that strongly lensed images with the largest angular separations are found in the densest regions of the COSMOS field. This is explained by a prevalence among the lens population in dense environments of elliptical galaxies with high total-to-stellar mass ratios, which can deflect light through larger angles. However, we also find that the overall fraction of elliptical galaxies with strong gravitational lensing is independent of the local mass density; this observation is not true of the simulations, which predict an increasing fraction of strong lenses in dense environments. The discrepancy may be a real effect, but could also be explained by various limitations of our analysis. For example, our visual search of strong lens systems could be incomplete and suffer from selection bias; the luminosity function of elliptical galaxies may differ between our real and simulated data; or the simplifying assumptions and approximations used in our lensing simulations may be inadequate. Work is therefore ongoing. Automated searches for strong lens systems will be particularly important in better constraining the selection function.

With increasingly large data sets, weak lensing measurements are able to measure cosmological parameters with ever greater precision. However this increased accuracy also places greater demands on the statistical tools used to extract the available information. To date, the majority of lensing analyses use the two point-statistics of the cosmic shear field. These can either be studied directly using the two-point correlation function, or in Fourier space, using the power spectrum. But analyzing weak lensing data inevitably involves the masking out of regions or example to remove bright stars from the field. Masking out the stars is common practice but the gaps in the data need proper handling. In this paper, we show how an inpainting technique allows us to properly fill in these gaps with only $N \log N$ operations, leading to a new image from which we can compute straight forwardly and with a very good accuracy both the pow er spectrum and the bispectrum. We propose then a new method to compute the bispectrum with a polar FFT algorithm, which has the main advantage of avoiding any interpolation in the Fourier domain. Finally we propose a new method for dark matter mass map reconstruction from shear observations which integrates this new inpainting concept. A range of examples based on 3D N-body simulations illustrates the results.

We have performed a 70 billion dark-matter particles N-body simulation in a 2 $h^{-1}$ Gpc periodic box, using the concordance, cosmological model as favored by the latest WMAP3 results. We have computed a full-sky convergence map with a resolution of $\Delta \theta \simeq 0.74$ arcmin$^{2}$, spanning 4 orders of magnitude in angular dynamical range. Using various high-order statistics on a realistic cut sky, we have characterized the transition from the linear to the nonlinear regime at $\ell \simeq 1000$ and shown that realistic galactic masking affects high-order moments only below $\ell < 200$. Each domain (Gaussian and non-Gaussian) spans 2 decades in angular scale. This map is therefore an ideal tool for testing map-making algorithms on the sphere. As a first step in addressing the full map reconstruction problem, we have benchmarked in this paper two denoising methods: 1) Wiener filtering applied to the Spherical Harmonics decomposition of the map and 2) a new method, called MRLens, based on the modification of the Maximum Entropy Method on a Wavelet decomposition. While the latter is optimal on large spatial scales, where the signal is Gaussian, MRLens outperforms the Wiener method on small spatial scales, where the signal is highly non-Gaussian. The simulated full-sky convergence map is freely available to the community to help the development of new map-making algorithms dedicated to the next generation of weak-lensing surveys.

Ordinary baryonic particles (such as protons and neutrons) account for only one-sixth of the total matter in the Universe. The remainder is a mysterious "dark matter" component, which does not interact via electromagnetism and thus neither emits nor reflects light. As dark matter cannot be seen directly using traditional observations, very little is currently known about its properties. It does interact via gravity, and is most effectively probed through gravitational lensing: the deflection of light from distant galaxies by the gravitational attraction of foreground mass concentrations. This is a purely geometrical effect that is free of astrophysical assumptions and sensitive to all matter -- whether baryonic or dark. Here we show high fidelity maps of the large-scale distribution of dark matter, resolved in both angle and depth. We find a loose network of filaments, growing over time, which intersect in massive structures at the locations of clusters of galaxies. Our results are consistent with predictions of gravitationally induced structure formation, in which the initial, smooth distribution of dark matter collapses into filaments then into clusters, forming a gravitational scaffold into which gas can accumulate, and stars can be built.

Context: The detection and identification of oscillation modes (in terms of their $\ell$, $m$ and successive $n$) is a great challenge for present and future asteroseismic space missions. The "peak tagging" is an important step in the analysis of these data to provide estimations of stellar oscillation mode parameters, i.e., frequencies, rotation rates, and further studies on the stellar structure. Aims: To increase the signal-to-noise ratio of the asteroseismic spectra computed from time series representative of MOST and CoRoT observations (30- and 150-day observations). Methods: We apply the curvelet transform -- a recent image processing technique which looks for curved patterns -- to echelle diagrams built using asteroseismic power spectra. In this diagram the eigenfrequencies appear as smooth continuous ridges. To test the method we use Monte Carlo simulations of several sun-like stars with different combinations of rotation rates, rotation-axis inclination and signal-to-noise ratios. Results: The filtered diagrams enhance the contrast between the ridges of the modes and the background allowing a better tagging of the modes and a better extraction of some stellar parameters. Monte Carlo simulations have also shown that the region where modes can be detected is enlarged at lower and higher frequencies compared to the raw spectra. Even more, the extraction of the mean rotational splitting from modes at low frequency can be done more easily than using the raw spectrum.

This paper presents a new method for the reconstruction of weak lensing mass maps. It uses the multiscale entropy concept, which is based on wavelets, and the False Discovery Rate which allows us to derive robust detection levels in wavelet space. We show that this new restoration approach outperforms several standard techniques currently used for weak shear mass reconstruction. This method can also be used to separate E and B modes in the shear field, and thus test for the presence of residual systematic effects. We concentrate on large blind cosmic shear surveys, and illustrate our results using simulated shear maps derived from N-Body Lambda-CDM simulations with added noise corresponding to both ground-based and space-based observations.

We present a new method for the reconstruction of Sunyaev-Zel'dovich (SZ) galaxy clusters in future SZ-survey experiments using multiband bolometer cameras such as Olimpo, APEX, or Planck. Our goal is to optimise SZ-Cluster extraction from our observed noisy maps. We wish to emphasize that none of the algorithms used in the detection chain is tuned on prior knowledge on the SZ -Cluster signal, or other astrophysical sources (Optical Spectrum, Noise Covariance Matrix, or covariance of SZ Cluster wavelet coefficients). First, a blind separation of the different astrophysical components which contribute to the observations is conducted using an Independent Component Analysis (ICA) method. Then, a recent non linear filtering technique in the wavelet domain, based on multiscale entropy and the False Discovery Rate (FDR) method, is used to detect and reconstruct the galaxy clusters. Finally, we use the Source Extractor software to identify the detected clusters. The proposed method was applied on realistic simulations of observations. As for global detection efficiency, this new method is impressive as it provides comparable results to Pierpaoli et al. method being however a blind algorithm. Preprint with full resolution figures is available at the URL: w10-dapnia.saclay.cea.fr/Phocea/Vie_des_labos/Ast/ast_visu.php?id_ast=728