Joe Zuntz

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

First Season QUIET Observations: Measurements of Cosmic Microwave Background Polarization Power Spectra at 43 GHz in the Multipole Range 25 ≤ ℓ ≤ 475

The Q/U Imaging ExperimenT (QUIET) employs coherent receivers at 43 GHz and 94 GHz, operating on the Chajnantor plateau in the Atacama Desert in Chile, to measure the anisotropy in the polarization of the cosmic microwave background (CMB). QUIET primarily targets the B modes from primordial gravitational waves. The combination of these frequencies gives sensitivity to foreground contributions from diffuse Galactic synchrotron radiation. Between 2008 October and 2010 December, over 10,000 hr of data were collected, first with the 19 element 43 GHz array (3458 hr) and then with the 90 element 94 GHz array. Each array observes the same four fields, selected for low foregrounds, together covering ≈1000 deg^2. This paper reports initial results from the 43 GHz receiver, which has an array sensitivity to CMB fluctuations of 69 μK√s. The data were extensively studied with a large suite of null tests before the power spectra, determined with two independent pipelines, were examined. Analysis choices, including data selection, were modified until the null tests passed. Cross-correlating maps with different telescope pointings is used to eliminate a bias. This paper reports the EE, BB, and EB power spectra in the multipole range l = 25-475. With the exception of the lowest multipole bin for one of the fields, where a polarized foreground, consistent with Galactic synchrotron radiation, is detected with 3σ significance, the E-mode spectrum is consistent with the ΛCDM model, confirming the only previous detection of the first acoustic peak. The B-mode spectrum is consistent with zero, leading to a measurement of the tensor-to-scalar ratio of r = 0.35^(+1.06)_(–0.87). The combination of a new time-stream double-demodulation technique, side-fed Dragonian optics, natural sky rotation, and frequent boresight rotation leads to the lowest level of systematic contamination in the B-mode power so far reported, below the level of r = 0.1.

DESTRUCTION OF MOLECULAR GAS RESERVOIRS IN EARLY-TYPE GALAXIES BY ACTIVE GALACTIC NUCLEUS FEEDBACK

Residual star formation at late times in early-type galaxies and their progenitors must be suppressed in order to explain the population of red, passively evolving systems we see today. Likewise, residual or newly accreted reservoirs of molecular gas that are fueling star formation must be destroyed. This suppression of star formation in early-type galaxies is now commonly attributed to active galactic nucleus (AGN) feedback wherein the reservoir of gas is heated and expelled during a phase of accretion onto the central supermassive black hole. However, direct observational evidence for a link between the destruction of this molecular gas and an AGN phase has been missing so far. We present new mm-wavelength observations from the IRAM 30 m telescope of a sample of low-redshift SDSS early-type galaxies currently undergoing this process of quenching of late-time star formation. Our observations show that the disappearance of the molecular gas coincides within less than 100 Myr with the onset of accretion onto the black hole and is too rapid to be due to star formation alone. Since our sample galaxies are not associated to powerful quasar activity or radio jets, we conclude that low-luminosity AGN episodes are sufficient to suppress residual star formation in early-type galaxies. This suppression mode of AGN feedback is very different from the truncation mode linked to powerful quasar activity during early phases of galaxy formation.

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

Weak gravitational lensing has the potential to constrain cosmological parameters to high nprecision. However, as shown by the Shear Testing Programmes and Gravitational lensing nAccuracy Testing challenges, measuring galaxy shears is a non-trivial task: various methods nintroduce different systematic biases which have to be accounted for. We investigate how npixel noise on the image affects the bias on shear estimates from a maximum likelihood nforward model-fitting approach using a sum of co-elliptical Sersic profiles, in complement nto the theoretical approach of an associated paper. We evaluate the bias using a simple but nrealistic galaxy model and find that the effects of noise alone can cause biases of the order of n1–10 per cent on measured shears, which is significant for current and future lensing surveys. nWe evaluate a simulation-based calibration method to create a bias model as a function of ngalaxy properties and observing conditions. This model is then used to correct the simulated nmeasurements. We demonstrate that, for the simple case in which the correct range of galaxy nmodels is used in the fit, the calibration method can reduce noise bias to the level required for nestimating cosmic shear in upcoming lensing surveys.

THE THIRD GRAVITATIONAL LENSING ACCURACY TESTING (GREAT3) CHALLENGE HANDBOOK

The GRavitational lEnsing Accuracy Testing 3 (GREAT3) challenge is the third in a series of image analysis challenges, with a goal of testing and facilitating the development of methods for analyzing astronomical images that will be used to measure weak gravitational lensing. This measurement requires extremely precise estimation of very small galaxy shape distortions, in the presence of far larger intrinsic galaxy shapes and distortions due to the blurring kernel caused by the atmosphere, telescope optics, and instrumental effects. The GREAT3 challenge is posed to the astronomy, machine learning, and statistics communities, and includes tests of three specific effects that are of immediate relevance to upcoming weak lensing surveys, two of which have never been tested in a community challenge before. These effects include many novel aspects including realistically complex galaxy models based on high-resolution imaging from space; a spatially varying, physically motivated blurring kernel; and a combination of multiple different exposures. To facilitate entry by people new to the field, and for use as a diagnostic tool, the simulation software for the challenge is publicly available, though the exact parameters used for the challenge are blinded. Sample scripts to analyze the challenge data using existing methods will also be provided. See http://great3challenge.info and http://great3.projects.phys.ucl.ac.uk/leaderboard/ for more information.

GREAT3 results - I. Systematic errors in shear estimation and the impact of real galaxy morphology

We present first results from the third GRavitational lEnsing Accuracy Testing (GREAT3) challenge, the third in a sequence of challenges for testing methods of inferring weak gravitational lensing shear distortions from simulated galaxy images. GREAT3 was divided into experiments to test three specific questions, and included simulated space- and ground-based data with constant or cosmologically varying shear fields. The simplest (control) experiment included parametric galaxies with a realistic distribution of signal-to-noise, size, and ellipticity, and a complex point spread function (PSF). The other experiments tested the additional impact of realistic galaxy morphology, multiple exposure imaging, and the uncertainty about a spatially varying PSF; the last two questions will be explored in Paper II. The 24 participating teams competed to estimate lensing shears to within systematic error tolerances for upcoming Stage-IV dark energy surveys, making 1525 submissions overall. GREAT3 saw considerable variety and innovation in the types of methods applied. Several teams now meet or exceed the targets in many of the tests conducted (to within the statistical errors). We conclude that the presence of realistic galaxy morphology in simulations changes shear calibration biases by similar to 1 per cent for a wide range of methods. Other effects such as truncation biases due to finite galaxy postage stamps, and the impact of galaxy type as measured by the S,rsic index, are quantified for the first time. Our results generalize previous studies regarding sensitivities to galaxy size and signal-to-noise, and to PSF properties such as seeing and defocus. Almost all methods results support the simple model in which additive shear biases depend linearly on PSF ellipticity.

Towards a fully consistent parameterization of modified gravity

There is a distinct possibility that current and future cosmological data can be used to constrain Einsteins theory of gravity on the very largest scales. To be able to do this in a model-independent way, it makes sense to work with a general parameterization of modified gravity. Such an approach would be analogous to the Parameterized Post-Newtonian (PPN) approach which is used on the scale of the Solar System. A few such parameterizations have been proposed and preliminary constraints have been obtained. We show that the majority of such parameterizations are only exactly applicable in the quasistatic regime. On larger scales they fail to encapsulate the full behaviour of typical models currently under consideration. We suggest that it may be possible to capture the additions to the `Parameterized Post-Friedmann (PPF) formalism by treating them akin to fluid perturbations.

CosmoSIS: modular cosmological parameter estimation

Cosmological parameter estimation is entering a new era. Large collaborations need to coordinate high-stakes analyses using multiple methods; furthermore such analyses have grown in complexity due to sophisticated models of cosmology and systematic uncertainties. In this paper we argue that modularity is the key to addressing these challenges: calculations should be broken up into interchangeable modular units with inputs and outputs clearly defined. We present a new framework for cosmological parameter estimation, CosmoSIS, designed to connect together, share, and advance development of inference tools across the community. We describe the modules already available in CosmoSIS, including CAMB, Planck, cosmic shear calculations, and a suite of samplers. We illustrate it using demonstration code that you can run out-of-the-box with the installer available at this http URL

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.

Vector field models of modified gravity and the dark sector

We present a comprehensive investigation of cosmological constraints on the class of vector field formulations of modified gravity called generalized Einstein-aether models. Using linear perturbation theory we generate cosmic microwave background and large-scale structure spectra for general parameters of the theory, and then constrain them in various ways. We investigate two parameter regimes: a dark matter candidate where the vector field sources structure formation, and a dark energy candidate where it causes late-time acceleration. We find that the dark matter candidate does not fit the data, and identify five physical problems that can restrict this and other theories of dark matter. The dark energy candidate does fit the data, and we constrain its fundamental parameters; most notably we find that the theorys kinetic index parameter n{sub ae} can differ significantly from its {Lambda}CDM value.