Source author record

Gary M. Bernstein

Gary M. Bernstein appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2021arXiv

A search of the full six years of the Dark Energy Survey for outer Solar System objects

We present the results of a search for outer Solar System objects in the full six years of data (Y6) from the Dark Energy Survey (DES). The DES covered a contiguous $5000$ deg$^2$ of the southern sky with $\approx 80,000$ $3$ deg$^2$ exposures in the $grizY$ optical/IR filters between 2013 and 2019. This search yielded 815 trans-Neptunian objects (TNOs), one Centaur and one Oort cloud comet, with 461 objects reported for the first time in this paper. We present methodology that builds upon our previous search carried out on the first four years of data. Here, all DES images were reprocessed with an improved detection pipeline that leads to an average completeness gain of 0.47 mag per exposure, as well as an improved transient catalog production and optimized algorithms for linkage of detections into orbits. All objects were verified by visual inspection and by computing the sub-threshold significance, the total signal-to-noise ratio in the stack of images in which the object's presence is indicated by the orbit fit, but no detection was reported. This yields a highly pure catalog of TNOs complete to $r \approx 23.8$ mag and distances $29<d<2500$ au. The Y6 TNOs have minimum (median) of 7 (12) distinct nights' detections and arcs of 1.1 (4.2) years, and will have $grizY$ magnitudes available in a further publication. We present software for simulating our observational biases that enable comparisons of population models to our detections. Initial inferences demonstrating the statistical power of the DES catalog are: the data are inconsistent with the CFEPS-L7 model for the classical Kuiper Belt; the 16 ``extreme'' TNOs ($a>150$ au, $q>30$ au) are consistent with the null hypothesis of azimuthal isotropy; and non-resonant TNOs with $q>38$ au, $a>50$ au show a highly significant tendency to be sunward of the major mean motion resonances, whereas this tendency is not present for $q<38$ au.

preprint2020arXiv

Propagating sample variance uncertainties in redshift calibration: simulations, theory and application to the COSMOS2015 data

Cosmological analyses of galaxy surveys rely on knowledge of the redshift distribution of their galaxy sample. This is usually derived from a spectroscopic and/or many-band photometric calibrator survey of a small patch of sky. The uncertainties in the redshift distribution of the calibrator sample include a contribution from shot noise, or Poisson sampling errors, but, given the small volume they probe, they are dominated by sample variance introduced by large-scale structures. Redshift uncertainties have been shown to constitute one of the leading contributions to systematic uncertainties in cosmological inferences from weak lensing and galaxy clustering, and hence they must be propagated through the analyses. In this work, we study the effects of sample variance on small-area redshift surveys, from theory to simulations to the COSMOS2015 data set. We present a three-step Dirichlet method of resampling a given survey-based redshift calibration distribution to enable the propagation of both shot noise and sample variance uncertainties. The method can accommodate different levels of prior confidence on different redshift sources. This method can be applied to any calibration sample with known redshifts and phenotypes (i.e. cells in a self-organizing map, or some other way of discretizing photometric space), and provides a simple way of propagating prior redshift uncertainties into cosmological analyses. As a worked example, we apply the full scheme to the COSMOS2015 data set, for which we also present a new, principled SOM algorithm designed to handle noisy photometric data. We make available a catalog of the resulting resamplings of the COSMOS2015 galaxies.

preprint2020arXiv

Testing the isotropy of the Dark Energy Survey's extreme trans-Neptunian objects

We test whether the population of "extreme" trans-Neptunian objects (eTNOs) detected in the Y4 Dark Energy Survey (DES) data exhibit azimuthal asymmetries which might be evidence of gravitational perturbations from an unseen super-Earth in a distant orbit. By rotating the orbits of the detected eTNOs, we construct a synthetic population which, when subject to the DES selection function, reproduces the detected distribution of eTNOs in the orbital elements $a,e,$ and $i$ as well as absolute magnitude $H$, but has uniform distributions in mean anomaly $M$, longitude of ascending node $Ω,$ and argument of perihelion $ω.$ We then compare the detected distributions in each of $Ω, ω,$ and $\varpi\equivΩ+ω$ to those expected from the isotropic population, using Kuiper's variant of the Kolmogorov-Smirnov test. The three angles are tested for each of 4 definitions of the eTNO population, choosing among $a>(150,250)$ AU and perihelion $q>(30,37)$ AU. These choices yield 3--7 eTNOs in the DES Y4 sample. Among the twelve total tests, two have the likelihood of drawing the observed angles from the isotropic population at $p<0.05.$ The 3 detections at $a>250, q>37$ AU, and the 4 detections at $a>250, q>30$ AU, have $Ω$ distribution with $p=0.03$ of coming from the isotropic construction, but this is not strong evidence of anisotropy given the 12 different tests. The DES data taken on their own are thus consistent with azimuthal isotropy and do not require a "Planet 9" hypothesis. The limited sky coverage and object count mean, however, that the DES data by no means falsify this hypothesis.

preprint2020arXiv

Trans-Neptunian objects found in the first four years of the Dark Energy Survey

We present a catalog of 316 trans-Neptunian bodies detected by the Dark Energy Survey (DES). These objects include 245 discoveries by DES (139 not previously published) detected in $\approx 60,000$ exposures from the first four seasons of the survey ("Y4" data). The survey covers a contiguous 5000 deg$^2$ of the southern sky in the $grizY$ optical/NIR filter set, with a typical TNO in this part of the sky being targeted by $25-30$ Y4 exposures. We describe the processes for detection of transient sources and the linkage into TNO orbits, which are made challenging by the absence of the few-hour repeat observations employed by TNO-optimized surveys. We also describe the procedures for determining detection efficiencies vs. magnitude and estimating rates of false-positive linkages. This work presents all TNOs which were detected on $\ge 6$ unique nights in the Y4 data and pass a "sub-threshold confirmation" test wherein we demand the the object be detectable in a stack of the individual images in which the orbit indicates an object should be present, but was not detected. This eliminates false positives and yields TNO detections complete to $r\lesssim 23.3$ mag with virtually no dependence on orbital properties for bound TNOs at distance $30\,{\rm AU}<d<2500\,{\rm AU}.$ The final DES TNO catalog is expected to yield $>0.3$ mag more depth, and arcs of $>4$ years for nearly all detections.

preprint2018arXiv

Evidence for Color Dichotomy in the Primordial Neptunian Trojan Population

In the current model of early Solar System evolution, the stable members of the Jovian and Neptunian Trojan populations were captured into resonance from the leftover reservoir of planetesimals during the outward migration of the giant planets. As a result, both Jovian and Neptunian Trojans share a common origin with the primordial disk population, whose other surviving members constitute today's trans-Neptunian object (TNO) populations. The cold classical TNOs are ultra-red, while the dynamically excited "hot" population of TNOs contains a mixture of ultra-red and blue objects. In contrast, Jovian and Neptunian Trojans are observed to be blue. While the absence of ultra-red Jovian Trojans can be readily explained by the sublimation of volatile material from their surfaces due to the high flux of solar radiation at 5AU, the lack of ultra-red Neptunian Trojans presents both a puzzle and a challenge to formation models. In this work we report the discovery by the Dark Energy Survey (DES) of two new dynamically stable L4 Neptunian Trojans,2013 VX30 and 2014 UU240, both with inclinations i >30 degrees, making them the highest-inclination known stable Neptunian Trojans. We have measured the colors of these and three other dynamically stable Neptunian Trojans previously observed by DES, and find that 2013 VX30 is ultra-red, the first such Neptunian Trojan in its class. As such, 2013 VX30 may be a "missing link" between the Trojan and TNO populations. Using a simulation of the DES TNO detection efficiency, we find that there are 162 +/- 73 Trojans with Hr < 10 at the L4 Lagrange point of Neptune. Moreover, the blue-to-red Neptunian Trojan population ratio should be higher than 17:1. Based on this result, we discuss the possible origin of the ultra-red Neptunian Trojan population and its implications for the formation history of Neptunian Trojans.

preprint2016arXiv

An accurate and practical method for inference of weak gravitational lensing from galaxy images

We demonstrate highly accurate recovery of weak gravitational lensing shear using an implementation of the Bayesian Fourier Domain (BFD) method proposed by Bernstein & Armstrong (2014, BA14), extended to correct for selection biases. The BFD formalism is rigorously correct for Nyquist-sampled, background-limited, uncrowded image of background galaxies. BFD does not assign shapes to galaxies, instead compressing the pixel data D into a vector of moments M, such that we have an analytic expression for the probability P(M|g) of obtaining the observations with gravitational lensing distortion g along the line of sight. We implement an algorithm for conducting BFD's integrations over the population of unlensed source galaxies which measures ~10 galaxies/second/core with good scaling properties. Initial tests of this code on ~10^9 simulated lensed galaxy images recover the simulated shear to a fractional accuracy of m=0.0021+-0.0004, substantially more accurate than has been demonstrated previously for any generally applicable method. Deep sky exposures generate a sufficiently accurate approximation to the noiseless, unlensed galaxy population distribution assumed as input to BFD. Potential extensions of the method include simultaneous measurement of magnification and shear; multiple-exposure, multi-band observations; and joint inference of photometric redshifts and lensing tomography.

preprint2015arXiv

Characterization and correction of charge-induced pixel shifts in DECam

Interaction of charges in CCDs with the already accumulated charge distribution causes both a flux dependence of the point-spread function (an increase of observed size with flux, also known as the brighter/fatter effect) and pixel-to-pixel correlations of the Poissonian noise in flat fields. We describe these effects in the Dark Energy Camera (DECam) with charge dependent shifts of effective pixel borders, i.e. the Antilogus et al. (2014) model, which we fit to measurements of flat-field Poissonian noise correlations. The latter fall off approximately as a power-law r^-2.5 with pixel separation r, are isotropic except for an asymmetry in the direct neighbors along rows and columns, are stable in time, and are weakly dependent on wavelength. They show variations from chip to chip at the 20% level that correlate with the silicon resistivity. The charge shifts predicted by the model cause biased shape measurements, primarily due to their effect on bright stars, at levels exceeding weak lensing science requirements. We measure the flux dependence of star images and show that the effect can be mitigated by applying the reverse charge shifts at the pixel level during image processing. Differences in stellar size, however, remain significant due to residuals at larger distance from the centroid.

preprint2015arXiv

GalSim: The modular galaxy image simulation toolkit

GALSIM is a collaborative, open-source project aimed at providing an image simulation tool of enduring benefit to the astronomical community. It provides a software library for generating images of astronomical objects such as stars and galaxies in a variety of ways, efficiently handling image transformations and operations such as convolution and rendering at high precision. We describe the GALSIM software and its capabilities, including necessary theoretical background. We demonstrate that the performance of GALSIM meets the stringent requirements of high precision image analysis applications such as weak gravitational lensing, for current datasets and for the Stage IV dark energy surveys of the Large Synoptic Survey Telescope, ESA's Euclid mission, and NASA's WFIRST-AFTA mission. The GALSIM project repository is public and includes the full code history, all open and closed issues, installation instructions, documentation, and wiki pages (including a Frequently Asked Questions section). The GALSIM repository can be found at https://github.com/GalSim-developers/GalSim .

preprint2014arXiv

Resampling images in Fourier domain

When simulating sky images, one often takes a galaxy image $F(x)$ defined by a set of pixelized samples and an interpolation kernel, and then wants to produce a new sampled image representing this galaxy as it would appear with a different point-spread function, a rotation, shearing, or magnification, and/or a different pixel scale. These operations are sometimes only possible, or most efficiently executed, as resamplings of the Fourier transform $\tilde F(u)$ of the image onto a $u$-space grid that differs from the one produced by a discrete Fourier transform (DFT) of the samples. In some applications it is essential that the resampled image be accurate to better than 1 part in $10^3$, so in this paper we first use standard Fourier techniques to show that Fourier-domain interpolation with a wrapped sinc function yields the exact value of $\tilde F(u)$ in terms of the input samples and kernel. This operation scales with image dimension as $N^4$ and can be prohibitively slow, so we next investigate the errors accrued from approximating the sinc function with a compact kernel. We show that these approximations produce a multiplicative error plus a pair of ghost images (in each dimension) in the simulated image. Standard Lanczos or cubic interpolators, when applied in Fourier domain, produce unacceptable artifacts. We find that errors $<1$ part in $10^3$ can be obtained by (1) 4-fold zero-padding of the original image before executing the $x\rightarrow u$ DFT, followed by (2) resampling to the desired $u$ grid using a 6-point, piecewise-quintic interpolant that we design expressly to minimize the ghosts, then (3) executing the DFT back to $x$ domain.

preprint2013arXiv

Bayesian Lensing Shear Measurement

We derive an estimator of weak gravitational lensing shear from background galaxy images that avoids noise-induced biases through a rigorous Bayesian treatment of the measurement. The derived shear estimator disposes with the assignment of ellipticities to individual galaxies that is typical of previous approaches to galaxy lensing. Shear estimates from the mean of the Bayesian posterior are unbiased in the limit of large number of background galaxies, regardless of the noise level on individual galaxies. The Bayesian formalism requires a prior describing the (noiseless) distribution of the target galaxy population over some parameter space; this prior can be constructed from low-noise images of a subsample of the target population, attainable from long integrations of a fraction of the survey field. We find two ways to combine this exact treatment of noise with rigorous treatment of the effects of the instrumental point-spread function and sampling. The Bayesian model fitting (BMF) method assigns a likelihood of the pixel data to galaxy models (e.g. Sersic ellipses). The Bayesian Fourier domain (BFD) method compresses the pixel data to a small set of weighted moments calculated after PSF correction in Fourier space. A numerical test using a simplified model of a biased galaxy measurement process demonstrates that the Bayesian formalism recovers applied shears to $<1$ part in $10^3$ accuracy and provides accurate uncertainty estimates. BFD is the first shear measurement algorithm that is model-free and requires no approximations or ad hoc assumptions in correcting for the effects of PSF, noise, or sampling on the galaxy images. These algorithms are good candidates for attaining the part-per-thousand shear inference required for hemisphere-scale weak gravitational lensing surveys. (abridged)

preprint2012arXiv

Atmospheric dispersion effects in weak lensing measurements

The wavelength dependence of atmospheric refraction causes elongation of finite-bandwidth images along the elevation vector, which produces spurious signals in weak gravitational lensing shear measurements unless this atmospheric dispersion is calibrated and removed to high precision. Because astrometric solutions and PSF characteristics are typically calibrated from stellar images, differences between the reference stars' spectra and the galaxies' spectra will leave residual errors in both the astrometric positions (dr) and in the second moment (width) of the wavelength-averaged PSF (dv) for galaxies.We estimate the level of dv that will induce spurious weak lensing signals in PSF-corrected galaxy shapes that exceed the statistical errors of the DES and the LSST cosmic-shear experiments. We also estimate the dr signals that will produce unacceptable spurious distortions after stacking of exposures taken at different airmasses and hour angles. We also calculate the errors in the griz bands, and find that dispersion systematics, uncorrected, are up to 6 and 2 times larger in g and r bands,respectively, than the requirements for the DES error budget, but can be safely ignored in i and z bands. For the LSST requirements, the factors are about 30, 10, and 3 in g, r, and i bands,respectively. We find that a simple correction linear in galaxy color is accurate enough to reduce dispersion shear systematics to insignificant levels in the r band for DES and i band for LSST,but still as much as 5 times than the requirements for LSST r-band observations. More complex corrections will likely be able to reduce the systematic cosmic-shear errors below statistical errors for LSST r band. But g-band effects remain large enough that it seems likely that induced systematics will dominate the statistical errors of both surveys, and cosmic-shear measurements should rely on the redder bands.

preprint2011arXiv

Cosmology without cosmic variance

We examine the improvements in constraints on the linear growth factor G and its derivative f=d ln G / dln a that are available from the combination of a large-scale galaxy redshift survey with a weak gravitational lensing survey of background sources. In the linear perturbation theory limit, the bias-modulation method of McDonald & Seljak allows one to distinguish the real-space galaxy clustering from the peculiar velocity signal in each Fourier mode. The ratio of lensing signal to galaxy clustering in transverse modes yields the bias factor b of each galaxy subset (as per Pen 2004), hence calibrating the conversion from galaxy real-space density to matter density in every mode. In combination these techniques permit measure of the growth rate f in each Fourier mode. This yields a measure of the growth rate free of sample variance, i.e. the uncertainty in f can be reduced without bound by increasing the number of redshifts within a finite volume. In practice, the gain from the absence of sample variance is bounded by the limited range of bias modulation among dark-matter halos. Nonetheless, the addition of background weak lensing data to a redshift survey increases information on G and f by an amount equivalent to a 10-fold increase in the volume of a standard redshift-space distortion measurement---if the lensing signal can be measured to sub-percent accuracy. This argues that a combined lensing and redshift survey over a common low-redshift volume is a more powerful test of general relativity than an isolated redshift survey over larger volume at high redshift. An example case is that a survey of ~10^6 redshifts over half the sky in the redshift range $z=0.5\pm 0.05$ can determine the growth exponent γfor the model $f=Ω_m^γ$ to an accuracy of $\pm 0.015$, using only modes with k<0.1h/Mpc, but only if a weak lensing survey is conducted in concert. [Abridged]

preprint2011arXiv

Lossy compression of weak lensing data

Future orbiting observatories will survey large areas of sky in order to constrain the physics of dark matter and dark energy using weak gravitational lensing and other methods. Lossy compression of the resultant data will improve the cost and feasibility of transmitting the images through the space communication network. We evaluate the consequences of the lossy compression algorithm of Bernstein et al. (2010) for the high-precision measurement of weak-lensing galaxy ellipticities. This square-root algorithm compresses each pixel independently, and the information discarded is by construction less than the Poisson error from photon shot noise. For simulated space-based images (without cosmic rays) digitized to the typical 16 bits per pixel, application of the lossy compression followed by image-wise lossless compression yields images with only 2.4 bits per pixel, a factor of 6.7 compression. We demonstrate that this compression introduces no bias in the sky background. The compression introduces a small amount of additional digitization noise to the images, and we demonstrate a corresponding small increase in ellipticity measurement noise. The ellipticity measurement method is biased by the addition of noise, so the additional digitization noise is expected to induce a multiplicative bias on the galaxies' measured ellipticities. After correcting for this known noise-induced bias, we find a residual multiplicative ellipticity bias of m ~ -4x10^{-4}. This bias is small when compared to the many other issues that precision weak lensing surveys must confront, and furthermore we expect it to be reduced further with better calibration of ellipticity measurement methods.

preprint2010arXiv

Improved Constraints on the Gravitational Lens Q0957+561. II. Strong Lensing

We present a detailed strong lensing analysis of an HST/ACS legacy dataset for the first gravitational lens, Q0957+561. With deep imaging we identify 24 new strongly lensed features, which we use to constrain mass models. We model the stellar component of the lens galaxy using the observed luminosity distribution, and the dark matter halo using several different density profiles. We draw on the weak lensing analysis by Nakajima et al. (2009) to constrain the mass sheet and environmental terms in the lens potential. Adopting the well-measured time delay, we find H_0 = 85 (+14/-13) km/s/Mpc (68% CL) using lensing constraints alone. The principal uncertainties in H_0 are tied to the stellar mass-to-light ratio (a variant of the radial profile degeneracy in lens models). Adding constraints from stellar population synthesis models, we obtain H_0 = 79.3 (+6.7/-8.5) km/s/Mpc (68% CL). We infer that the lens galaxy has a rising rotation curve and a dark matter distribution with an inner core. Intriguingly, we find the quasar flux ratios predicted by our models to be inconsistent with existing radio measurements, suggesting the presence of substructure in the lens.

preprint2010arXiv

Noise and bias in square-root compression schemes

We investigate data compression schemes for proposed all-sky diffraction-limited visible/NIR sky surveys aimed at the dark energy problem. We show that lossy square-root compression to 1 bit of noise per pixel, followed by standard lossless compression algorithms, reduces the images to 2.5-4 bits per pixel, depending primarily upon the level of cosmic-ray contamination of the images. Compression to this level adds noise equivalent to <10% penalty in observing time. We derive an analytic correction to flux biases inherent to the square-root compression scheme. Numerical tests on simple galaxy models confirm that galaxy fluxes and shapes are measured with systematic biases <~10^{-4} induced by the compression scheme, well below the requirements of supernova and weak gravitational lensing dark-energy experiments. An accompanying paper (Vanderveld 2009) bounds the shape biases using realistic simulated images of the high-Galactic-latitude sky. The square-root preprocessing step has advantages over simple (linear) decimation when there are many bright objects or cosmic rays in the field, or when the background level will vary.

preprint2010arXiv

Shape measurement biases from underfitting and ellipticity gradients

Precision weak gravitational lensing experiments require measurements of galaxy shapes accurate to <1 part in 1000. We investigate measurement biases, noted by Voigt and Bridle (2009) and Melchior et al. (2009), that are common to shape measurement methodologies that rely upon fitting elliptical-isophote galaxy models to observed data. The first bias arises when the true galaxy shapes do not match the models being fit. We show that this "underfitting bias" is due, at root, to these methods' attempts to use information at high spatial frequencies that has been destroyed by the convolution with the point-spread function (PSF) and/or by sampling. We propose a new shape-measurement technique that is explicitly confined to observable regions of k-space. A second bias arises for galaxies whose ellipticity varies with radius. For most shape-measurement methods, such galaxies are subject to "ellipticity gradient bias". We show how to reduce such biases by factors of 20--100 within the new shape-measurement method. The resulting shear estimator has multiplicative errors <1 part in 1000 for high-S/N images, even for highly asymmetric galaxies. Without any training or recalibration, the new method obtains Q=3000 in the GREAT08 Challenge of blind shear reconstruction on low-noise galaxies, several times better than any previous method.

preprint2009arXiv

The cosmology dependence of weak lensing cluster counts

We present the main results of a numerical study of weak lensing cluster counting. We examine the scaling with cosmology of the projected-density-peak mass function. Our main conclusion is that the projected-peak and the three-dimensional mass functions scale with cosmology in an astonishingly close way. This means that, despite being derived from a two-dimensional field, the weak lensing cluster abundance can be used to constrain cosmology in the same way as the three-dimensional mass function probed by other types of surveys.

preprint2009arXiv

The impact of correlated projections on weak lensing cluster counts

Large-scale structure projections are an obstacle in converting the shear signal of clusters detected in weak-lensing maps into virial masses. However, this step is not necessary for constraining cosmology with the shear-peak abundance, if we are able to predict its amplitude. We generate a large ensemble of N-body simulations spanning four cosmological models, with total volume V~1 (Gpc/h)^3 per model. Variations to the matter density parameter and amplitude of fluctuations are considered. We measure the abundance of peaks in the mass density projected in ~100 Mpc/h slabs to determine the impact of structures spatially correlated with the simulation clusters, identified by the 3D friends-of-friends algorithm. The halo model shows that the choice of the smoothing filter for the density field is important in reducing the contribution of correlated projections to individual halo masses. Such contributions are less than 2% in the case of the optimal, compensated filter used throughout this analysis. We measure the change in the mass of peaks when projected in slabs of various thicknesses. Peaks in slabs of 26 Mpc/h and 102 Mpc/h suffer an average mass change of less than 2% compared to their mass in slabs of 51 Mpc/h. We then explore the cosmology dependence of the projected-peak mass function, and find that, for a wide range of slab thicknesses (<500 Mpc/h), it scales with cosmology in exactly the same way as the 3D friends-of-friends mass function and the Sheth-Tormen formula. This extends the earlier result of Marian et al. (2009). Finally, we show that for all cosmological models considered, the low and intermediate mass bins of the peak abundance can be described using a modified Sheth-Tormen functional form to within 10%-20% accuracy.

preprint2008arXiv

Comprehensive Two-Point Analyses of Weak Gravitational Lensing Surveys

We present a framework for analyzing weak gravitational lensing survey data, including lensing and source-density observables, plus spectroscopic redshift calibration data. All two-point observables are predicted in terms of parameters of a perturbed Robertson-Walker metric, making the framework independent of the models for gravity, dark energy, or galaxy properties. For Gaussian fluctuations the 2-point model determines the survey likelihood function and allows Fisher-matrix forecasting. The framework includes nuisance terms for the major systematic errors: shear measurement errors, magnification bias and redshift calibration errors, intrinsic galaxy alignments, and inaccurate theoretical predictions. We propose flexible parameterizations of the many nuisance parameters related to galaxy bias and intrinsic alignment. For the first time we can integrate many different observables and systematic errors into a single analysis. As a first application of this framework, we demonstrate that: uncertainties in power-spectrum theory cause very minor degradation to cosmological information content; nearly all useful information (excepting baryon oscillations) is extracted with ~3 bins per decade of angular scale; and the rate at which galaxy bias varies with redshift substantially influences the strength of cosmological inference. The framework will permit careful study of the interplay between numerous observables, systematic errors, and spectroscopic calibration data for large weak-lensing surveys.

preprint2008arXiv

Multipole Formulae for Gravitational Lensing Shear and Flexion

The gravitational lensing equations for convergence, potential, shear, and flexion are simple in polar coordinates and separate under a multipole expansion once the shear and flexion spinors are rotated into a ``tangential'' basis. We use this to investigate whether the useful monopole aperture-mass shear formulae generalize to all multipoles and to flexions. We re-derive the result of Schneider and Bartelmann that the shear multipole m at radius R is completely determined by the mass multipole at R, plus specific moments Q^m_in and Q^m_out of the mass multipoles internal and external, respectively, to R. The m>=0 multipoles are independent of Q_out. But in contrast to the monopole, the m<0 multipoles are independent of Q_in. These internal and external mass moments can be determined by shear (and/or flexion) data on the complementary portion of the plane, which has practical implications for lens modelling. We find that the ease of E/B separation in the monopole aperture moments does {\em not} generalize to m!=0: the internal monopole moment is the {\em only} non-local E/B discriminant available from lensing observations. We have also not found practical local E/B discriminants beyond the monopole, though they could exist. We show also that the use of weak-lensing data to constrain a constant shear term near a strong-lensing system is impractical without strong prior constraints on the neighboring mass distribution.