Lisa Voigt

Thermal conduction and reduced cooling flows in galaxy clusters

Conduction may play an important role in reducing cooling flows in galaxy clusters. We analyse a sample of 16 objects using Chandra data and find that a balance between electron conduction and cooling can exist in the hotter clusters (T 5 keV), provided that the plasma conductivity is close to the unhindered Spitzer value. In the absence of any additional heat sources, a reduced mass inflow must then develop in the cooler objects in the sample. We fit cooling flow models to deprojected spectra 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 measured mass inflow rates are insufficient 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 and find conductivity profiles that follow a strikingly similar temperature dependence to the conductivity values required to prevent cooling. The larger-scale turbulent motions implied by this process are required to have velocities of between 10 and 50 per cent of the speed of sound in the local intracluster gas.

Results of the GREAT08 Challenge: an image analysis competition for cosmological lensing

We present the results of the Gravitational LEnsing Accuracy Testing 2008 (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.

Image analysis for cosmology: results from the GREAT10 Galaxy Challenge

We present the results from the first public blind point-spread function (PSF) reconstruction challenge, the GRavitational lEnsing Accuracy Testing 2010 (GREAT10) Star Challenge. Reconstruction of a spatially varying PSF, sparsely sampled by stars, at non-star positions is a critical part in the image analysis for weak lensing where inaccuracies in the modeled ellipticity e and size R^2 can impact the ability to measure the shapes of galaxies. This is of importance because weak lensing is a particularly sensitive probe of dark energy and can be used to map the mass distribution of large scale structure. Participants in the challenge were presented with 27,500 stars over 1300 images subdivided into 26 sets, where in each set a category change was made in the type or spatial variation of the PSF. Thirty submissions were made by nine teams. The best methods reconstructed the PSF with an accuracy of σ(e) ≈ 2.5 × 10^(–4) and σ(R^2)/R^2 ≈ 7.4 × 10^(–4). For a fixed pixel scale, narrower PSFs were found to be more difficult to model than larger PSFs, and the PSF reconstruction was severely degraded with the inclusion of an atmospheric turbulence model (although this result is likely to be a strong function of the amplitude of the turbulence power spectrum).

Handbook for the GREAT08 Challenge: An image analysis competition for cosmological lensing

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 www.great08challenge.info for the latest information.

Point spread function calibration requirements for dark energy from cosmic shear

Context. The control of systematic effects when measuring background galaxy shapes is one of the main challenges for cosmic shear analyses. Aims. Study the fundamental limitations on shear accuracy due to the measurement of the point spread function (PSF) from the finite number of stars that are available. We translate the accuracy required for cosmological parameter estimation to the minimum number of stars over which the PSF must be calibrated. Methods. We characterise the error made in the shear arising from errors on the PSF. We consider different PSF models, from a simple elliptical Gaussian to various shapelet parametrisations. First we 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 the elliptical Gaussian PSF. Results. 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. This means that these small scales should not be used to constrain cosmology unless the instrument and the observing strategy are optimised to make this variability extremely small. The minimum number of stars varies with the square of the star Signal-to-Noise Ratio, with the complexity of the PSF and with the pixel size. 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. Conclusions. 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. Due to the simplicity of our models these results should be interpreted as optimistic and therefore provide a measure of a systematic “floor’ intrinsic to shape measurements.

Measurement and calibration of noise bias in weak lensing galaxy shape estimation

Weak gravitational lensing has the potential to constrain cosmological parameters to high precision. However, as shown by the Shear Testing Programmes and Gravitational lensing Accuracy Testing challenges, measuring galaxy shears is a non-trivial 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 Sersic profiles, in complement to the theoretical approach of 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 the order of 1–10 per cent 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, for the simple case in which the correct range of galaxy models is used in the fit, the calibration method can reduce noise bias to the level required for estimating cosmic shear in upcoming lensing surveys.

Limitations of model‐fitting methods for lensing shear estimation

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 disc, if these two components have different ellipticities. The limiting accuracy is dependent on the galaxy shape, but we find the most significant biases (similar to 1 per cent of the shear) 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.

IM3SHAPE: A maximum-likelihood galaxy shear measurement code for cosmic gravitational lensing

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 bulgeplus-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 like lihood methods. We assess the bias arising from fitting an incorrect galaxy model using simple n oise-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 obt ain a negligible additive shear bias and sub-two percent level multiplicative bias, which i s 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 futur e experiments. The code is freely available from https://bitbucket.org/joezuntz/im3shap e.

The impact of galaxy colour gradients on cosmic shear measurement

Cosmic shear has been identified as the method with the most potential to constrain dark energy. To capitalize on this potential, it is necessary to measure galaxy shapes with great accuracy, which in turn requires a detailed model for the image blurring by the telescope and atmosphere, 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-centred, 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 the 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 from Simard et al., 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 signal-to-noise ratio of 25 the number of galaxies required in each tomographic redshift bin is of the order of 10^4.

Cosmic shear requirements on the wavelength-dependence of telescope point spread functions

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 introduce 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.