Deconvolution with Shapelets

We seek to find a shapelet-based scheme for deconvolving galaxy images from the PSF which leads to unbiased shear measurements. Based on the analytic formulation of convolution in shapelet space, we construct a procedure to recover the unconvolved shapelet coefficients under the assumption that the PSF is perfectly known. Using specific simulations, we test this approach and compare it to other published approaches. We show that convolution in shapelet space leads to a shapelet model of order $n_{max}^h = n_{max}^g + n_{max}^f$ with $n_{max}^f$ and $n_{max}^g$ being the maximum orders of the intrinsic galaxy and the PSF models, respectively. Deconvolution is hence a transformation which maps a certain number of convolved coefficients onto a generally smaller number of deconvolved coefficients. By inferring the latter number from data, we construct the maximum-likelihood solution for this transformation and obtain unbiased shear estimates with a remarkable amount of noise reduction compared to established approaches. This finding is particularly valid for complicated PSF models and low $S/N$ images, which renders our approach suitable for typical weak-lensing conditions.Comment: 9 pages, 9 figures, submitted to A&

The Dark UNiverse Explorer (DUNE): Proposal to ESA's Cosmic Vision

The Dark UNiverse Explorer (DUNE) is a wide-field space imager whose primary goal is the study of dark energy and dark matter with unprecedented precision. For this purpose, DUNE is optimised for the measurement of weak gravitational lensing but will also provide complementary measurements of baryonic accoustic oscillations, cluster counts and the Integrated Sachs Wolfe effect. Immediate auxiliary goals concern the evolution of galaxies, to be studied with unequalled statistical power, the detailed structure of the Milky Way and nearby galaxies, and the demographics of Earth-mass planets. DUNE is an Medium-class mission which makes use of readily available components, heritage from other missions, and synergy with ground based facilities to minimise cost and risks. The payload consists of a 1.2m telescope with a combined visible/NIR field-of-view of 1 deg^2. DUNE will carry out an all-sky survey, ranging from 550 to 1600nm, in one visible and three NIR bands which will form a unique legacy for astronomy. DUNE will yield major advances in a broad range of fields in astrophysics including fundamental cosmology, galaxy evolution, and extrasolar planet search. DUNE was recently selected by ESA as one of the mission concepts to be studied in its Cosmic Vision programme.Comment: Accepted in Experimental Astronom

Limitations for shapelet-based weak-lensing measurements

We seek to understand the impact on shape estimators obtained from circular and elliptical shapelet models under two realistic conditions: (a) only a limited number of shapelet modes is available for the model, and (b) the intrinsic galactic shapes are not restricted to shapelet models. We create a set of simplistic simulations, in which the galactic shapes follow a Sersic profile. By varying the Sersic index and applied shear, we quantify the amount of bias on shear estimates which arises from insufficient modeling. Additional complications due to PSF convolution, pixelation and pixel noise are also discussed. Steep and highly elliptical galaxy shapes cannot be accurately modeled within the circular shapelet basis system and are biased towards shallower and less elongated shapes. This problem can be cured partially by allowing elliptical basis functions, but for steep profiles elliptical shapelet models still depend critically on accurate ellipticity priors. As a result, shear estimates are typically biased low. Independently of the particular form of the estimator, the bias depends on the true intrinsic galaxy morphology, but also on the size and shape of the PSF. As long as the issues discussed here are not solved, the shapelet method cannot provide weak-lensing measurements with an accuracy demanded by upcoming missions and surveys, unless one can provide an accurate and reliable calibration, specific for the dataset under investigation.Comment: 8 pages, 5 figures, submitted to A&

Modal decomposition of astronomical images with application to shapelets

The decomposition of an image into a linear combination of digitised basis functions is an everyday task in astronomy. A general method is presented for performing such a decomposition optimally into an arbitrary set of digitised basis functions, which may be linearly dependent, non-orthogonal and incomplete. It is shown that such circumstances may result even from the digitisation of continuous basis functions that are orthogonal and complete. In particular, digitised shapelet basis functions are investigated and are shown to suffer from such difficulties. As a result the standard method of performing shapelet analysis produces unnecessarily inaccurate decompositions. The optimal method presented here is shown to yield more accurate decompositions in all cases.Comment: 12 pages, 17 figures, submitted to MNRA

Prospects for Dark Energy Evolution: a Frequentist Multi-Probe Approach

A major quest in cosmology is the understanding of the nature of dark energy. It is now well known that a combination of cosmological probes is required to break the underlying degeneracies on cosmological parameters. In this paper, we present a method, based on a frequentist approach, to combine probes without any prior constraints, taking full account of the correlations in the parameters. As an application, a combination of current SNIa and CMB data with an evolving dark energy component is first compared to other analyses. We emphasise the consequences of the implementation of the dark energy perturbations on the result for a time varying equation of state. The impact of future weak lensing surveys on the measurement of dark energy evolution is then studied in combination with future measurements of the cosmic microwave background and type Ia supernovae. We present the combined results for future mid-term and long-term surveys and confirm that the combination with weak lensing is very powerful in breaking parameter degeneracies. A second generation of experiment is however required to achieve a 0.1 error on the parameters describing the evolution of dark energy.Comment: Submitted to Astronomy & Astrophysics 14 pages, 8 figure

Effect of Gravitational Lensing on Measurements of the Sunyaev-Zel'dovich Effect

The Sunyaev-Zel'dovich (SZ) effect of a cluster of galaxies is usually measured after background radio sources are removed from the cluster field. Gravitational lensing by the cluster potential leads to a systematic deficit in the residual intensity of unresolved sources behind the cluster core relative to a control field far from the cluster center. As a result, the measured decrement in the Rayleigh-Jeans temperature of the cosmic microwave background is overestimated. We calculate the associated systematic bias which is inevitably introduced into measurements of the Hubble constant using the SZ effect. For the cluster A2218, we find that observations at 15 GHz with a beam radius of 0'.4 and a source removal threshold of 100 microJy underestimate the Hubble constant by 6-10%. If the profile of the gas pressure declines more steeply with radius than that of the dark matter density, then the ratio of lensing to SZ decrements increases towards the outer part of the cluster.Comment: 11 pages, 3 figures, submitted to ApJ

Reliable Shapelet Image Analysis

Aims: We discuss the applicability and reliability of the shapelet technique for scientific image analysis. Methods: We quantify the effects of non-orthogonality of sampled shapelet basis functions and misestimation of shapelet parameters. We perform the shapelet decomposition on artificial galaxy images with underlying shapelet models and galaxy images from the GOODS survey, comparing the publicly available IDL implementation with our new C++ implementation. Results: Non-orthogonality of the sampled basis functions and misestimation of the shapelet parameters can cause substantial misinterpretation of the physical properties of the decomposed objects. Additional constraints, image preprocessing and enhanced precision have to be incorporated in order to achieve reliable decomposition results.Comment: 12 pages, 15 figures, revised version, accepted by A&

A fast empirical method for galaxy shape measurements in weak lensing surveys

We describe a simple and fast method to correct ellipticity measurements of galaxies from the distortion by the instrumental and atmospheric point spread function (PSF), in view of weak lensing shear measurements. The method performs a classification of galaxies and associated PSFs according to measured shape parameters, and corrects the measured galaxy ellipticites by querying a large lookup table (LUT), built by supervised learning. We have applied this new method to the GREAT10 image analysis challenge, and present in this paper a refined solution that obtains the competitive quality factor of Q = 104, without any shear power spectrum denoising or training. Of particular interest is the efficiency of the method, with a processing time below 3 ms per galaxy on an ordinary CPU.Comment: 8 pages, 6 figures. Metric values updated according to the final GREAT10 analysis software (Kitching et al. 2012, MNRAS 423, 3163-3208), no qualitative changes. Associated code available at http://lastro.epfl.ch/megalu

Gravitational Lensing of the X-Ray Background by Clusters of Galaxies

Gravitational lensing by clusters of galaxies affects the cosmic X-ray background (XRB) by altering the observed density and flux distribution of background X-ray sources. At faint detection flux thresholds, the resolved X-ray sources appear brighter and diluted, while the unresolved component of the XRB appears dimmer and more anisotropic, due to lensing. The diffuse X-ray intensity in the outer halos of clusters might be lower than the sky-averaged XRB, after the subtraction of resolved sources. Detection of the lensing signal with a wide-field X-ray telescope could probe the mass distribution of a cluster out to its virialization boundary. In particular, we show that the lensing signature imprinted on the resolved component of the XRB by the cluster A1689, should be difficult but possible to detect out to 8' at the 2-4 sigma level, after 10^6 seconds of observation with the forthcoming AXAF satellite. The lensing signal is fairly insensitive to the lens redshift in the range 0.1<z<0.6. The amplitude of the lensing signal is however sensitive to the faint end slope of the number-flux relation for unresolved X-ray sources, and can thus help constrain models of the XRB. A search for X-ray arcs or arclets could identify the fraction of all faint sources which originate from extended emission of distant galaxies. The probability for a 3 sigma detection of an arclet which is stretched by a factor of about 3 after a 10^6 seconds observation of A1689 with AXAF, is roughly comparable to the fraction of all background X-ray sources that have an intrinsic size of order 1''.Comment: 41 LaTeX pages, 11 postscript figures, 1 table, in AASTeX v4.0 format. To appear in ApJ, April 1, 1997, Vol. 47