Cosmic shear is a powerful probe of cosmological parameters, but its potential can be fully utilised only if galaxy shapes are measured with great accuracy. Two major effects have been identified which are likely to account for most of the bias for maximum likelihood methods in recent shear measurement challenges. Model bias occurs when the true galaxy shape is not well represented by the fitted model. Noise bias occurs due to the non-linear relationship between image pixels and galaxy shape. In this paper we investigate the potential interplay between these two effects when an imperfect model is used in the presence of high noise. We present analytical expressions for this bias, which depends on the residual difference between the model and real data. They can lead to biases not accounted for in previous calibration schemes. By measuring the model bias, noise bias and their interaction, we provide a complete statistical framework for measuring galaxy shapes with model fitting methods from GRavitational lEnsing Accuracy Testing (GREAT) like images. We demonstrate the noise and model interaction bias using a simple toy model, which indicates that this effect can potentially be significant. Using real galaxy images from the Cosmological Evolution Survey (COSMOS) we quantify the strength of the model bias, noise bias and their interaction. We find that the interaction term is often a similar size to the model bias term, and is smaller than the requirements of the current and shortly upcoming galaxy surveys.

In this paper we present results from the weak lensing shape measurement GRavitational lEnsing Accuracy Testing 2010 (GREAT10) Galaxy Challenge. This marks an order of magnitude step change in the level of scrutiny employed in weak lensing shape measurement analysis. We provide descriptions of each method tested and include 10 evaluation metrics over 24 simulation branches. GREAT10 was the first shape measurement challenge to include variable fields; both the shear field and the Point Spread Function (PSF) vary across the images in a realistic manner. The variable fields enable a variety of metrics that are inaccessible to constant shear simulations including a direct measure of the impact of shape measurement inaccuracies, and the impact of PSF size and ellipticity, on the shear power spectrum. To assess the impact of shape measurement bias for cosmic shear we present a general pseudo-Cl formalism, that propagates spatially varying systematics in cosmic shear through to power spectrum estimates. We also show how one-point estimators of bias can be extracted from variable shear simulations. The GREAT10 Galaxy Challenge received 95 submissions and saw a factor of 3 improvement in the accuracy achieved by shape measurement methods. The best methods achieve sub-percent average biases. We find a strong dependence in accuracy as a function of signal-to-noise, and indications of a weak dependence on galaxy type and size. Some requirements for the most ambitious cosmic shear experiments are met above a signal-to-noise ratio of 20. These results have the caveat that the simulated PSF was a ground-based PSF. Our results are a snapshot of the accuracy of current shape measurement methods and are a benchmark upon which improvement can continue. This provides a foundation for a better understanding of the strengths and limitations of shape measurement methods.

We present and describe im3shape, a new publicly available galaxy shape measurement code for weak gravitational lensing shear. im3shape performs a maximum likelihood fit of a bulge-plus-disc galaxy model to noisy images, incorporating an applied point spread function. We detail challenges faced and choices made in its design and implementation, and then discuss various limitations that affect this and other maximum likelihood methods. We assess the bias arising from fitting an incorrect galaxy model using simple noise-free images and find that it should not be a concern for current cosmic shear surveys. We test im3shape on the GREAT08 Challenge image simulations, and meet the requirements for upcoming cosmic shear surveys in the case that the simulations are encompassed by the fitted model, using a simple correction for image noise bias. For the fiducial branch of GREAT08 we obtain a negligible additive shear bias and sub-two percent level multiplicative bias, which is suitable for analysis of current surveys. We fall short of the sub-percent level requirement for upcoming surveys, which we attribute to a combination of noise bias and the mis-match between our galaxy model and the model used in the GREAT08 simulations. We meet the requirements for current surveys across all branches of GREAT08, except those with small or high noise galaxies, which we would cut from our analysis. Using the GREAT08 metric we we obtain a score of Q=717 for the usable branches, relative to the goal of Q=1000 for future experiments. The code is freely available from https://bitbucket.org/joezuntz/im3shape

Weak gravitational lensing has the potential to constrain cosmological parameters to high precision. However, as shown by the Shear TEsting Programmes (STEP) and GRavitational lEnsing Accuracy Testing (GREAT) Challenges, measuring galaxy shears is a nontrivial task: various methods introduce different systematic biases which have to be accounted for. We investigate how pixel noise on the image affects the bias on shear estimates from a Maximum-Likelihood forward model-fitting approach using a sum of co-elliptical S\'{e}rsic profiles, in complement to the theoretical approach of an an associated paper. We evaluate the bias using a simple but realistic galaxy model and find that the effects of noise alone can cause biases of order 1-10% on measured shears, which is significant for current and future lensing surveys. We evaluate a simulation-based calibration method to create a bias model as a function of galaxy properties and observing conditions. This model is then used to correct the simulated measurements. We demonstrate that this method can effectively reduce noise bias so that shear measurement reaches the level of accuracy required for estimating cosmic shear in upcoming lensing surveys.

GRavitational lEnsing Accuracy Testing 2010 (GREAT10) is a public image analysis challenge aimed at the development of algorithms to analyze astronomical images. Specifically, the challenge is to measure varying image distortions in the presence of a variable convolution kernel, pixelization and noise. This is the second in a series of challenges set to the astronomy, computer science and statistics communities, providing a structured environment in which methods can be improved and tested in preparation for planned astronomical surveys. GREAT10 extends upon previous work by introducing variable fields into the challenge. The "Galaxy Challenge" involves the precise measurement of galaxy shape distortions, quantified locally by two parameters called shear, in the presence of a known convolution kernel. Crucially, the convolution kernel and the simulated gravitational lensing shape distortion both now vary as a function of position within the images, as is the case for real data. In addition, we introduce the "Star Challenge" that concerns the reconstruction of a variable convolution kernel, similar to that in a typical astronomical observation. This document details the GREAT10 Challenge for potential participants. Continually updated information is also available from http://www.greatchallenges.info.

Cosmic shear has been identified as the method with the most potential to constrain dark energy. To capitalise on this potential it is necessary to measure galaxy shapes with great accuracy, which in turn requires a detailed model for the image blurring, the Point Spread Function (PSF). In general the PSF varies with wavelength and therefore the PSF integrated over an observing filter depends on the spectrum of the object. For a typical galaxy the spectrum varies across the galaxy image, thus the PSF depends on the position within the image. We estimate the bias on the shear due to such colour gradients by modelling galaxies using two co-centered, co-elliptical Sersic profiles, each with a different spectrum. We estimate the effect of ignoring colour gradients and find the shear bias from a single galaxy can be very large depending on the properties of the galaxy. We find that halving the filter width reduces the shear bias by a factor of about 5. We show that, to first order, tomographic cosmic shear two point statistics depend on the mean shear bias over the galaxy population at a given redshift. For a single broad filter, and averaging over a small galaxy catalogue, we find a mean shear bias which is subdominant to the predicted statistical errors for future cosmic shear surveys. However, the true mean shear bias may exceed the statistical errors, depending on how accurately the catalogue represents the observed distribution of galaxies in the cosmic shear survey. We then investigate the bias on the shear for two-filter imaging and find that the bias is reduced by at least an order of magnitude. Lastly, we find that it is possible to calibrate galaxies for which colour gradients were ignored using two-filter imaging of a fair sample of noisy galaxies, if the galaxy model is known. For a S/N of 25 the number of galaxies required in each tomographic redshift bin is of order 1e4.

Cosmic shear requires high precision measurement of galaxy shapes in the presence of the observational Point Spread Function (PSF) that smears out the image. The PSF must therefore be known for each galaxy to a high accuracy. However, for several reasons, the PSF is usually wavelength dependent, therefore the differences between the spectral energy distribution of the observed objects introduces further complexity. In this paper we investigate the effect of the wavelength-dependence of the PSF, focusing on instruments in which the PSF size is dominated by the diffraction-limit of the telescope and which use broad-band filters for shape measurement. We first calculate biases on cosmological parameter estimation from cosmic shear when the stellar PSF is used uncorrected. Using realistic galaxy and star spectral energy distributions and populations and a simple three-component circular PSF we find that the colour-dependence must be taken into account for the next generation of telescopes. We then consider two different methods for removing the effect (i) the use of stars of the same colour as the galaxies and (ii) estimation of the galaxy spectral energy distribution using multiple colours and using a telescope model for the PSF. We find that both of these methods correct the effect to levels below the tolerances required for per-cent level measurements of dark energy parameters. Comparison of the two methods favours the template-fitting method because its efficiency is less dependent on galaxy redshift than the broad-band colour method and takes full advantage of deeper photometry.

We present the results of the GREAT08 Challenge, a blind analysis challenge to infer weak gravitational lensing shear distortions from images. The primary goal was to stimulate new ideas by presenting the problem to researchers outside the shear measurement community. Six GREAT08 Team methods were presented at the launch of the Challenge and five additional groups submitted results during the 6 month competition. Participants analyzed 30 million simulated galaxies with a range in signal to noise ratio, point-spread function ellipticity, galaxy size, and galaxy type. The large quantity of simulations allowed shear measurement methods to be assessed at a level of accuracy suitable for currently planned future cosmic shear observations for the first time. Different methods perform well in different parts of simulation parameter space and come close to the target level of accuracy in several of these. A number of fresh ideas have emerged as a result of the Challenge including a re-examination of the process of combining information from different galaxies, which reduces the dependence on realistic galaxy modelling. The image simulations will become increasingly sophisticated in future GREAT challenges, meanwhile the GREAT08 simulations remain as a benchmark for additional developments in shear measurement algorithms.

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.

Gravitational lensing shear has the potential to be the most powerful tool for constraining the nature of dark energy. However, accurate measurement of galaxy shear is crucial and has been shown to be non-trivial by the Shear TEsting Programme. Here we demonstrate a fundamental limit to the accuracy achievable by model-fitting techniques, if oversimplistic models are used. We show that even if galaxies have elliptical isophotes, model-fitting methods which assume elliptical isophotes can have significant biases if they use the wrong profile. We use noise-free simulations to show that on allowing sufficient flexibility in the profile the biases can be made negligible. This is no longer the case if elliptical isophote models are used to fit galaxies made up of a bulge plus a disk, if these two components have different ellipticities. The limiting accuracy is dependent on the galaxy shape but we find the most significant biases for simple spiral-like galaxies. The implications for a given cosmic shear survey will depend on the actual distribution of galaxy morphologies in the universe, taking into account the survey selection function and the point spread function. However our results suggest that the impact on cosmic shear results from current and near future surveys may be negligible. Meanwhile, these results should encourage the development of existing approaches which are less sensitive to morphology, as well as methods which use priors on galaxy shapes learnt from deep surveys.

The control of systematic effects when measuring galaxy shapes is one of the main challenges for cosmic shear analyses. In this context, we study the fundamental limitations on shear accuracy due to the measurement of the Point Spread Function (PSF) from the finite number of stars. In order to do that, we translate the accuracy required for cosmological parameter estimation to the minimum number of stars over which the PSF must be calibrated. We first derive our results analytically in the case of infinitely small pixels (i.e. infinitely high resolution). Then image simulations are used to validate these results and investigate the effect of finite pixel size in the case of an elliptical gaussian PSF. Our results are expressed in terms of the minimum number of stars required to calibrate the PSF in order to ensure that systematic errors are smaller than statistical errors when estimating the cosmological parameters. On scales smaller than the area containing this minimum number of stars, there is not enough information to model the PSF. In the case of an elliptical gaussian PSF and in the absence of dithering, 2 pixels per PSF Full Width at Half Maximum (FWHM) implies a 20% increase of the minimum number of stars compared to the ideal case of infinitely small pixels; 0.9 pixels per PSF FWHM implies a factor 100 increase. In the case of a good resolution and a typical Signal-to-Noise Ratio distribution of stars, we find that current surveys need the PSF to be calibrated over a few stars, which may explain residual systematics on scales smaller than a few arcmins. Future all-sky cosmic shear surveys require the PSF to be calibrated over a region containing about 50 stars.

Accurate measurements of the mass distribution in galaxy and cluster halos are essential to test the cold dark matter (CDM) paradigm. The cosmological model predicts a universal shape for the density profile in all halos, independent of halo mass. Its profile has a `cuspy' centre, with no evidence for the constant density core. In this paper we carry out a careful analysis of twelve galaxy clusters, using Chandra data to compute the mass distribution in each system under the assumption of hydrostatic equilibrium. Due to their low concentration, clusters provide ideal objects for studying the central cusps in dark matter halos. The majority of the systems are consistent with the CDM model, but 4 objects exhibit flat inner density profiles. We suggest that the flat inner profile found for these clusters is due to an underestimation of the mass in the cluster centre (rather than any problem with the CDM model), since these objects also have a centrally peaked gas mass fraction. We discuss possible causes for erroneously low mass measurements in the cores of some systems.

Conduction may play an important role in reducing cooling flows in galaxy clusters. We analyse a sample of sixteen objects using Chandra data and find that a balance between conduction and cooling can exist in the hotter clusters (T > 5 keV), provided the plasma conductivity is close to the unhindered Spitzer value. In the absence of any additional heat sources, a reduced mass inflow must develop in the cooler objects in the sample. We fit cooling flow models to deprojected data and compare the spectral mass deposition rates found to the values required to account for the excess luminosity, assuming Spitzer-rate heat transfer over the observed temperature gradients. The mass inflow rates found are lower than is necessary to maintain energy balance in at least five clusters. However, emission from cooling gas may be partially absorbed. We also compute the flux supplied by turbulent heat transport (Cho et al. 2003) and find conductivity profiles which follow a strikingly similar temperature dependence to the conductivity values required to prevent cooling. Finally, we show that the cluster radio luminosities vary by over five orders of magnitude in objects with X-ray luminosities differing by no more than a factor of a few. This suggests that there is unlikely to be a straightforward correlation between the mechanical power provided by the radio lobes and the rate of energy loss in cooling flow clusters.

On the basis of the universal gas fraction in clusters of galaxies, we estimate that the effective thermal conductivity required to balance radiative cooling in the cores, where the gas temperature is 3-10keV, is about one tenth of the Spitzer rate. This confirms that thermal conduction can be important for the energy balance provided that it is not highly suppressed by magnetic fields in the gas. We determine the global effective conductivity in a sample of 29 clusters using published X-ray data on the inferred cooling rates and show that most lie between one and one tenth of the Spitzer rate. More work on the profiles in cooling flow clusters is required to test the conduction hypothesis further. We examine the possibility that conduction operates during galaxy formation, and show that it provides a simple explanation for the upper-mass cutoff in galaxy masses.

Chandra and XMM-Newton observations have confirmed the presence of large temperature gradients within the cores of many relaxed clusters of galaxies. Here we investigate whether thermal conduction operating over those gradients can supply sufficient heat to offset radiative cooling. Narayan & Medvedev (2001) and Gruzinov (2002) have noted, using published results on cluster temperatures, that conduction within a factor of a few of the Spitzer rate is sufficient to balance bremsstrahlung cooling. From a detailed study of the temperature and emission measure profiles of Abell 2199 and Abell 1835, we find that the heat flux required by conduction is consistent with or below the rate predicted by Spitzer in the outer regions of the core. Conduction may therefore explain the lack of observational evidence for large mass cooling rates inferred from arguments based simply on radiative cooling, provided that conductivity is suppressed by no more than a factor of three below the full Spitzer rate. To stem cooling in the cluster centre, however, would necessitate conductivity values at least a factor of two larger than the Spitzer values, which we consider implausible. This may provide an explanation for the observed star formation and optical nebulosities in cluster cores. The solution is likely to be time dependent. We briefly discuss the possible origin of the cooler gas and the implications for massive galaxies.