Current measurements of the low and high redshift Universe are in tension if we restrict ourselves to the standard six parameter model of flat $\Lambda$CDM. This tension has two parts. First, the Planck satellite data suggest a higher normalization of matter perturbations than local measurements of galaxy clusters. Second, the expansion rate of the Universe today, $H_0$, derived from local distance-redshift measurements is significantly higher than that inferred using the acoustic scale in galaxy surveys and the Planck data as a standard ruler. The addition of a sterile neutrino species changes the acoustic scale and brings the two into agreement; meanwhile, adding mass to the active neutrinos or to a sterile neutrino can suppress the growth of structure, bringing the cluster data into better concordance as well. For our fiducial dataset combination, with statistical errors for clusters, a model with a massive sterile neutrino shows 3.5$\sigma$ evidence for a non-zero mass and an even stronger rejection of the minimal model. A model with massive active neutrinos and a massless sterile neutrino is similarly preferred. An eV-scale sterile neutrino mass -- of interest for short baseline and reactor anomalies -- is well within the allowed range. We caution that 1) unknown astrophysical systematic errors in any of the data sets could weaken this conclusion, but they would need to be several times the known errors to eliminate the tensions entirely; 2) the results we find are at some variance with analyses that do not include cluster measurements; and 3) some tension remains among the datasets even when new neutrino physics is included.

Given recent indications of additional neutrino species and cosmologically significant neutrino masses, we analyze their signatures in the weak lensing shear power spectrum. We find that a shear deficit in the 20-40% range or excess in the 20-60% range cannot be explained by variations in parameters of the flat LambdaCDM model that are allowed by current observations of the expansion history from Type Ia supernovae, baryon acoustic oscillations, and local measures of the Hubble constant H_0, coupled with observations of the cosmic microwave background from WMAP9 and the SPT 2500 square degree survey. Hence such a shear deficit or excess would indicate large masses or extra species, respectively, and we find this to be independent of the flatness assumption. We also discuss the robustness of these predictions to cosmic acceleration physics and the means by which shear degeneracies in joint variation of mass and species can be broken.

We present a method for attaining sub-arcsecond pointing stability during sub- orbital balloon flights, as designed for in the High Altitude Lensing Observatory (HALO) concept. The pointing method presented here has the potential to perform near-space quality optical astronomical imaging at 1-2% of the cost of space-based missions. We also discuss an architecture that can achieve sufficient thermomechanical stability to match the pointing stability. This concept is motivated by advances in the development and testing of Ultra Long Duration Balloon (ULDB) flights which promise to allow observation campaigns lasting more than three months. The design incorporates a multi-stage pointing architecture comprising: a gondola coarse azimuth control system, a multi-axis nested gimbal frame structure with arcsecond stability, a telescope de-rotator to eliminate field rotation, and a fine guidance stage consisting of both a telescope mounted angular rate sensor and guide CCDs in the focal plane to drive a fast-steering mirror. We discuss the results of pointing tests together with a preliminary thermo-mechanical analysis required for sub-arcsecond pointing at high altitude. Possible future applications in the areas of wide-field surveys and exoplanet searches are also discussed.

Any theory invoked to explain cosmic acceleration predicts consistency relations between the expansion history, structure growth, and all related observables. Currently there exist high-quality measurements of the expansion history from Type Ia supernovae, the cosmic microwave background temperature and polarization spectra, and baryon acoustic oscillations. We can use constraints from these datasets to predict what future probes of structure growth should observe. We apply this method to predict what range of cosmic shear power spectra would be expected if we lived in a LambdaCDM universe, with or without spatial curvature, and what results would be inconsistent and therefore falsify the model. Though predictions are relaxed if one allows for an arbitrary quintessence equation of state $-1\le w(z)\le 1$, we find that any observation that rules out LambdaCDM due to excess lensing will also rule out all quintessence models, with or without early dark energy. We further explore how uncertainties in the nonlinear matter power spectrum, e.g. from approximate fitting formulas such as Halofit, warm dark matter, or baryons, impact these limits.

We study the fluctuations in luminosity distances due to gravitational lensing by large scale (> 35 Mpc) structures, specifically voids and sheets. We use a simplified "Swiss cheese" model consisting of a \Lambda -CDM Friedman-Robertson-Walker background in which a number of randomly distributed non-overlapping spherical regions are replaced by mass compensating comoving voids, each with a uniform density interior and a thin shell of matter on the surface. We compute the distribution of magnitude shifts using a variant of the method of Holz & Wald (1998), which includes the effect of lensing shear. The standard deviation of this distribution is ~ 0.027 magnitudes and the mean is ~ 0.003 magnitudes for voids of radius 35 Mpc, sources at redshift z_s=1.0, with the voids chosen so that 90% of the mass is on the shell today. The standard deviation varies from 0.005 to 0.06 magnitudes as we vary the void size, source redshift, and fraction of mass on the shells today. If the shell walls are given a finite thickness of ~ 1 Mpc, the standard deviation is reduced to ~ 0.013 magnitudes. This standard deviation due to voids is a factor ~ 3 smaller than that due to galaxy scale structures. We summarize our results in terms of a fitting formula that is accurate to ~ 20%, and also build a simplified analytic model that reproduces our results to within ~ 30%. Our model also allows us to explore the domain of validity of weak lensing theory for voids. We find that for 35 Mpc voids, corrections to the dispersion due to lens-lens coupling are of order ~ 4%, and corrections to due shear are ~ 3%. Finally, we estimate the bias due to source-lens clustering in our model to be negligible.

We explore the sensitivity of weak gravitational lensing to second-order corrections to the spacetime metric within a cosmological adaptation of the parameterized post-Newtonian framework. Whereas one might expect nonlinearities of the gravitational field to introduce non-Gaussianity into the statistics of the lensing convergence field, we show that such corrections are actually always small within a broad class of scalar-tensor theories of gravity. We show this by first computing the weak lensing convergence within our parameterized framework to second order in the gravitational potential, and then computing the relevant post-Newtonian parameters for scalar-tensor gravity theories. In doing so we show that this potential systematic factor is generically negligible, thus clearing the way for weak lensing to provide a direct tracer of mass on cosmological scales for a wide class of gravity theories despite uncertainties in the precise nature of the departures from general relativity.

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.

We present the 'simage' software suite for the simulation of artificial extragalactic images, based empirically around real observations of the Hubble Ultra Deep Field (UDF). The simulations reproduce galaxies with realistic and complex morphologies via the modeling of UDF galaxies as shapelets. Images can be created in the B, V, i and z bands for both space- and ground-based telescopes and instruments. The simulated images can be produced for any required field size, exposure time, Point Spread Function (PSF), telescope mirror size, pixel resolution, field star density, and a variety of detector noise sources. It has the capability to create images with both a pre-determined number of galaxies or one calibrated to the number counts of pre-existing data sets such as the HST COSMOS survey. In addition, simple options are included to add a known weak gravitational lensing (both shear and flexion) to the simulated images. The software is available in Interactive Data Language (IDL) and can be freely downloaded for scientific, developmental and teaching purposes.

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.

We argued in a previous paper [R. A. Vanderveld et al. 2006, arXiv:astro-ph/0602476] that negative deceleration parameters at the center of symmetry in Lemaitre-Tolman-Bondi cosmological models can only occur if the model is not smooth at the origin. Here we demonstrate explicitly the connection between non-smoothness and the failure of positivity theorems for deceleration. We also address some confusion that has arisen in the literature and respond to some recent criticisms of our arguments.

Recently there have been suggestions that the Type Ia supernova data can be explained using only general relativity and cold dark matter with no dark energy. In "Swiss cheese" models of the Universe, the standard Friedmann-Robertson-Walker picture is modified by the introduction of mass compensating spherical inhomogeneities, typically described by the Lemaitre-Tolman-Bondi metric. If these inhomogeneities correspond to underdense cores surrounded by mass-compensating overdense shells, then they can modify the luminosity distance-redshift relation in a way that can mimic accelerated expansion. It has been argued that this effect could be large enough to explain the supernova data without introducing dark energy or modified gravity. We show that the large apparent acceleration seen in some models can be explained in terms of standard weak field gravitational lensing together with insufficient randomization of void locations. The underdense regions focus the light less than the homogeneous background, thus dimming supernovae in a way that can mimic the effects of acceleration. With insufficient randomization of the spatial location of the voids and of the lines of sight, coherent defocusing can lead to anomalously large demagnification effects. We show that a proper randomization of the voids and lines of sight reduces the effect to the point that it can no longer explain the supernova data.

The fitting of the observed redshifts and magnitudes of type Ia supernovae to what we would see in homogeneous cosmological models has led to constraints on cosmological parameters. However, in doing such fits it is assumed that the sampled supernovae are moving with the Hubble flow, i.e. that their peculiar velocities are zero. In reality, peculiar velocities will modify supernova data in a way that can impact best-fit cosmological parameters. We theoretically quantify this effect in the nonlinear regime with a Monte-Carlo analysis, using data from semi-analytic galaxy catalogs that are built from the Millennium N-body simulation. We find scaling relations for the errors in best-fit parameters resulting solely from peculiar velocities, as a function of the total number of sources in a supernova survey N and its maximum redshift z_max. For low redshift surveys, we find that these errors can be of the same order of magnitude as the errors due to an intrinsic magnitude scatter of 0.1 mag. For a survey with N=2000 and z_max=1.7, we estimate that the expected peculiar velocity-induced errors in the best-fit cosmological constant density and equation of state can be sigma_Lambda~0.009 and sigma_w~0.01, respectively, which are subdominant to the errors due to the intrinsic scatter. We further find that throwing away supernova data below a redshift z~0.01-0.02 can reduce the combined error, due to peculiar velocities and the intrinsic scatter, but by only about 10%.

We calculate the systematic inhomogeneity-induced correction to the cosmological constant that one would infer from an analysis of the luminosities and redshifts of Type Ia supernovae, assuming a homogeneous universe. The calculation entails a post-Newtonian expansion within the framework of second order perturbation theory, wherein we consider the effects of subhorizon density perturbations in a flat, dust dominated universe. Within this formalism, we calculate luminosity distances and redshifts along the past light cone of an observer. The resulting luminosity distance-redshift relation is fit to that of a homogeneous model in order to deduce the best-fit cosmological constant density Omega_Lambda. We find that the luminosity distance-redshift relation is indeed modified, by a small fraction of order 10^{-5}. When fitting this perturbed relation to that of a homogeneous universe, we find that the inferred cosmological constant can be surprisingly large, depending on the range of redshifts sampled. For a sample of supernovae extending from z=0.02 out to z=0.15, we find that Omega_Lambda=0.004. The value of Omega_Lambda has a large variance, and its magnitude tends to get larger for smaller redshifts, implying that precision measurements from nearby supernova data will require taking this effect into account. However, we find that this effect is likely too small to explain the observed value of Omega_Lambda=0.7. There have been previous claims of much larger backreaction effects. By contrast to those calculations, our work is directly related to how observers deduce cosmological parameters from astronomical data.

There has been much debate over whether or not one could explain the observed acceleration of the Universe with inhomogeneous cosmological models, such as the spherically-symmetric Lemaitre-Tolman-Bondi (LTB) models. It has been claimed that the central observer in these models can observe a local acceleration, which would contradict general theorems. We resolve the contradiction by noting that many of the models that have been explored contain a weak singularity at the location of the observer which makes them unphysical. In the absence of this singularity, we show that LTB models must have a positive central deceleration parameter $q_{0}$, in agreement with the general theorems. We also show that it is possible to achieve a negative apparent deceleration parameter at nonzero redshifts in LTB models that do not contain this singularity. However, we find other singularities that tend to arise in LTB models when attempting to match luminosity distance data, and these generally limit the range of redshifts for which these models can mimic observations of an accelerating Universe. Exceptional models do exist that can extend to arbitrarily large redshift without encountering these pathologies, and we show how these may be constructed. These special models exhibit regions with negative effective equation of state parameter, which may fall below negative one, but we have failed to find any singularity-free models that agree with observations. Moreover, models based on dust-filled LTB metrics probably fail to reproduce observed properties of large scale structure.

Spherical gravitational collapse of a cold gas of annihilating particles involves a competition between the free-fall rate $\propto\sqrt{\rho}$ and the (s-wave) annihilation rate $\propto\rho$. Thus, there is a critical density $\rhoann$ above which annihilation proceeds faster than free fall. Gravitational collapse of a cloud of (initial) mass $M$ to a black hole is only possible if $3/32\pi G^3M^2\lesssim\rhoann$, or $M\gtrsim\Mann\equiv (3/32\pi G^3\rhoann)^{1/2}$. For a particle mass $m$ and freeze-out temperature $T_f=m/x_f$, the minimum black hole mass is $\Mann\approx 10^{10}\msun \times(x_f\sqrt{g_\star}/100\omcdm g_{\star S}m({\rm Gev}))$, where $g_{\star S}$ and $g_\star$ are degeneracy factors. The formation of a black hole of initial mass $M_{BH}$ is accompanied by the annihilation of about $M_{ann}$ released in a burst lasting a time $\sim GM_{BH}$ that could reach a total annihilation luminosity $\sim 10^{59} {\rm erg s^{-1}}$. The absence of astronomical observations of such spectacular events suggests either: (i) the branching ratio for CDM annihilation to $e^{+}e^{-}$ pairs or quarks $\lesssim 10^{-10}$, while the branching ratio to $\nu{\bar{\nu}}$ is $\lesssim 10^{-5}$; or (ii) CDM is not made of annihilating particles, but may be in some non-annihilating form, such as axions; or (iii) CDM black holes never form.