We constrain X-ray spectral shapes for the ensemble of AGN based on the shape of the Cosmic X-ray Background (CXB). Specifically, we rule out regions of X-ray spectral parameter space that do not reproduce the CXB in the energy range 1-100 keV. The key X-ray spectral parameters are the photon index, {\Gamma}; the cutoff energy, Ecutoff; and the reflection scaling factor, R. Assuming each parameter follows a Gaussian distribution, we first explore the parameter space using a Bayesian approach and a fixed X-ray luminosity function (XLF). For {\sigma}_E = 36 keV and {\sigma}_R = 0.14, fixed at the observed values from the Swift-BAT 70-month sample, we allow <R>, <Ecutoff > and <{\Gamma}> to vary subject to reproducing the CXB. We report results for {\sigma}_{\Gamma} = 0.1-0.5. In an alternative approach, we define the parameter distributions, then forward model to fit the CXB by perturbing the XLF using a neural network. This approach allows us to rule out parameter combinations that cannot reproduce the CXB for any XLF. The marginalized conditional probabilities for the four free parameters are: <R> = 0.99^{+0.11}_{-0.26}, <Ecutoff> = 118^{+24}_{-23}, {\sigma}_{\Gamma} = 0.101^{+0.097}_{-0.001} and <{\Gamma}> = 1.9^{+0.08}_{-0.09}. We provide an interactive online tool for users to explore any combination of <Ecutoff>, {\sigma}_E, <{\Gamma}>, {\sigma}_{\Gamma}, <R> and {\sigma}_R including different distributions for each absorption bin, subject to the integral CXB constraint. The distributions observed in many AGN samples can be ruled out by our analysis, meaning these samples can not be representative of the full AGN population. The few samples that fall within the acceptable parameter space are hard X-ray-selected, commensurate with their having fewer selection biases.

MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at 2<z<5. A 6.5-m Magellan telescope will be coupled with DESI spectrographs to achieve multiplexing of 20,000. MegaMapper would be located at Las Campanas Observatory to fully access LSST imaging for target selection.

Tight constraints on the abundance of primordial black holes can be deduced across a vast range of masses, with the exception of those light enough to fully evaporate before nucleosynthesis. This hypothetical population is almost entirely unconstrained, to the point where the early Universe could pass through a matter-dominated phase with primordial black holes as the primary component. The only obvious relic of this phase would be Hawking radiated gravitons which would constitute a stochastic gravitational wave background in the present-day Universe, albeit at frequencies far beyond the scope of any planned detector technology. This paper explores the effects of classical mergers in such a matter dominated phase. For certain ranges of parameters, a significant fraction of the black holes merge, providing an additional, classical source of primordial gravitational waves. The resulting stochastic background typically has a lower amplitude than the Hawking background and lies at less extreme frequencies, but is unlikely to be easily detectable, with a maximal present day density of $\Omega_{GW} \sim 10^{-12}$ and frequencies between $10^{15} - 10^{19}$ Hz. We also asses the impact of radiation accretion on the lifetimes of such primordial black holes and find that it increases the black hole mass by $\sim 14 \%$ and the lifetimes by about $50 \%$. However, this does not qualitatively change any of our conclusions.

21-cm intensity surveys aim to map neutral hydrogen atoms in the universe through hyper-fine emission. Unfortunately, long-wavelength (low-wavenumber) radial modes are highly contaminated by smooth astrophysical foregrounds that are six orders of magnitude brighter than the cosmological signal. This contamination also leaks into higher radial and angular wavenumber modes and forms a foreground wedge. Cosmic tidal reconstruction aims to extract the large-scale signal from anisotropic features in the local small-scale power spectrum through non-linear tidal interactions; losing small-scale modes to foreground wedge will impair its performance. In this paper, we review tidal interaction theory and estimator construction, and derive the theoretical expressions for the reconstructed spectra. We show the reconstruction is robust against peculiar velocities. Removing low line-of-sight $k$ modes, we demonstrate cross-correlation coefficient $r$ is greater than 0.7 on large scales ($k <0.1$ $h/$Mpc) even with a cutoff value $k^c_{\|}=0.1$ $h/$Mpc. Discarding wedge modes yields $0.3< r < 0.5$ and completely removes the dependency on $k^c_{\|}$. Our theoretical predictions agree with these numerical simulations.

The nature of dark energy remains a profound mystery 20 years after the discovery of cosmic acceleration. A very high number density galaxy redshift survey over a wide area (HD GRS Wide) spanning the redshift range of 0.5<z<4 using the same tracer, carried out using massively parallel wide field multi-object slit spectroscopy from space, will provide definitive dark energy measurements with minimal observational systematics by design. The HD GRS Wide will illuminate the nature of dark energy, and lead to revolutionary advances in particle physics and cosmology. It will also trace the cosmic web of dark matter and provide key insight into large-scale structure in the Universe. The required observational facility can be established as part of the probe portfolio by NASA within the next decade.

Upcoming imaging surveys, such as LSST, will provide an unprecedented view of the Universe, but with limited resolution along the line-of-sight. Common ways to increase resolution in the third dimension, and reduce misclassifications, include observing a wider wavelength range and/or combining the broad-band imaging with higher spectral resolution data. The challenge with these approaches is matching the depth of these ancillary data with the original imaging survey. However, while a full 3D map is required for some science, there are many situations where only the statistical distribution of objects (dN/dz) in the line-of-sight direction is needed. In such situations, there is no need to measure the fluxes of individual objects in all of the surveys. Rather a stacking procedure can be used to perform an `ensemble photo-z'. We show how a shallow, higher spectral resolution survey can be used to measure dN/dz for stacks of galaxies which coincide in a deeper, lower resolution survey. The galaxies in the deeper survey do not even need to appear individually in the shallow survey. We give a toy model example to illustrate tradeoffs and considerations for applying this method. This approach will allow deep imaging surveys to leverage the high resolution of spectroscopic and narrow/medium band surveys underway, even when the latter do not have the same reach to high redshift.

Cosmic acceleration is the most surprising cosmological discovery in many decades. Testing and distinguishing among possible explanations requires cosmological measurements of extremely high precision probing the full history of cosmic expansion and structure growth and, ideally, compare and contrast matter and relativistic tracers of the gravity potential. This program is one of the defining objectives of the Wide-Field Infrared Survey Telescope (WFIRST), as set forth in the New Worlds, New Horizons report (NWNH) in 2010. The WFIRST mission has the ability to improve these measurements by 1-2 orders of magnitude compared to the current state of the art, while simultaneously extending their redshift grasp, greatly improving control of systematic effects, and taking a unified approach to multiple probes that provide complementary physical information and cross-checks of cosmological results. We describe in this annual report the activities of the Science Investigation Team (SIT) "Cosmology with the High Latitude Survey (HLS)" during the year 2017. This team was selected by NASA in December 2015 in order to address the stringent challenges of the WFIRST dark energy (DE) program through the Project's formulation phase. This SIT has elected to jointly address Galaxy Redshift Survey, Weak Lensing and Cluster Growth and thus fully embrace the fact that the imaging and spectroscopic elements of the HLS will be realized as an integrated observing program, and they jointly impose requirements on performance and operations. WFIRST is designed to be able to deliver a definitive result on the origin of cosmic acceleration. It is not optimized for Figure of Merit sensitivity but for control of systematic uncertainties and for having multiple techniques each with multiple cross-checks. Our SIT work focuses on understanding the potential systematics in the WFIRST DE measurements.

We measure the growth rate and its evolution using the anisotropic clustering of the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 (DR14) quasar sample, which includes $148\,659$ quasars covering the wide redshift range of $0.8 < z < 2.2$ and a sky area of $2112.90$ $\rm deg^2$. To optimise measurements we deploy a redshift-dependent weighting scheme, which allows us to avoid binning, and perform the data analysis consistently including the redshift evolution across the sample. We perform the analysis in Fourier space, and use the redshift evolving power spectrum multipoles to measure the redshift space distortion parameter $f\sigma_8$ and parameters controlling the anisotropic projection of the cosmological perturbations. We measure $f \sigma_8(z=1.52)=0.43 \pm 0.05 $ and $df\sigma_8/dz (z=1.52)= - 0.16 \pm 0.08$, consistent with the expectation for a $\Lambda$CDM cosmology as constrained by the Planck experiment.

We present an anisotropic analysis of Baryon Acoustic Oscillation (BAO) signal from the SDSS-IV extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 (DR14) quasar sample. The sample consists of 147,000 quasars distributed over a redshift range of $0.8 < z < 2.2$. We apply the redshift weights technique to the clustering of quasars in this sample and achieve a 4.6 per cent measurement of the angular distance measurement $D_M$ at $z = 2.2$ and Hubble parameter $H$ at $z=0.8$. We parameterize the distance-redshift relation, relative to a fiducial model, as a quadratic expansion. The coefficients of this expansion are used to reconstruct the distance-redshift relation and obtain distance and Hubble parameter measurements at all redshifts within the redshift range of the sample. Reporting the result at two characteristic redshifts, we determine $D_M(z=1) = 3405\pm305 \ (r_{\rm d} / r_{\rm d, fid}) \ {\rm Mpc}$, $H(z=1) = 120.7\pm 7.3 \ (r_{\rm d,fid} / r_{\rm d}) \ {\rm km} \ {\rm s}^{-1}{\rm Mpc}^{-1}$ and $D_M(z=2) = 5325\pm249 \ (r_{\rm d} / r_{\rm d, fid}) \ {\rm Mpc}$, $H(z=2) = 189.9\pm 32.9 \ (r_{\rm d,fid} / r_{\rm d}) \ {\rm km} \ {\rm s}^{-1}{\rm Mpc}^{-1}$. These measurements are highly correlated. We assess the outlook of BAO analysis from the final quasar sample by testing the method on a set of mocks that mimic the noise level in the final sample. We demonstrate on these mocks that redshift weighting shrinks the measurement error by over 25 per cent on average. We conclude redshift weighting can bring us closer to the cosmological goal of the final quasar sample.

The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is one of the first of a new generation of galaxy redshift surveys that will cover a large range in redshift with sufficient resolution to measure the baryon acoustic oscillations (BAO) signal. For surveys covering a large redshift range we can no longer ignore cosmological evolution, meaning that either the redshift shells analysed have to be significantly narrower than the survey, or we have to allow for the averaging over evolving quantities. Both of these have the potential to remove signal: analysing small volumes increases the size of the Fourier window function, reducing the large-scale information, while averaging over evolving quantities can, if not performed carefully, remove differential information. It will be important to measure cosmological evolution from these surveys to explore and discriminate between models. We apply a method to optimally extract this differential information to mock catalogues designed to mimic the eBOSS quasar sample. By applying a set of weights to extract redshift space distortion measurements as a function of redshift, we demonstrate an analysis that does not invoke the problems discussed above. We show that our estimator gives unbiased constraints.

We present measurements of the anisotropic galaxy clustering from the Data Release 9 (DR9) CMASS sample of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). We analyze the broad-range shape of the monopole and quadrupole correlation functions to obtain constraints, at the effective redshift $z=0.57$ of the sample, on the Hubble expansion rate $H(z)$, the angular-diameter distance $D_A(z)$, the normalized growth rate $f(z)\sigma_8(z)$, the physical matter density $\Omega_m h^2$, and the biased amplitude of matter fluctuation b\sigma_8(z). We obtain {$H(0.57)$, $D_A(0.57)$, $f(0.57)\sigma_8(0.57)$, $\Omega_m h^2$, $b\sigma_8(0.57)$} = {$87.6_{-6.8}^{+6.7}$, $1396\pm73$, $0.126_{-0.010}^{+0.008}$, $1.19\pm0.14$, $0.428\pm0.066$} and their covariance matrix as well. The parameters which are not well constrained by our of galaxy clustering analysis are marginalized over with wide flat priors. Since no priors from other data sets (i.e., CMB) are adopted and no dark energy models are assumed, our results from BOSS CMASS galaxy clustering alone may be combined with other data sets, i.e. CMB, SNe, lensing or other galaxy clustering data to constrain the parameters of a given cosmological model. We show that the major power on constraining dark energy from the anisotropic galaxy clustering signal, as compared to the angular-averaged one (monopole), arises from including the normalized growth rate $f(z)\sigma_8(z)$. In the case of the wCDM cosmological model our single-probe CMASS constraints, combined with CMB (WMAP9+SPT), yield a value for the dark energy equation of state parameter of $w=-0.90\pm0.11$. Therefore, it is important to include $f(z)\sigma_8(z)$ while investigating the nature of dark energy with current and upcoming large-scale galaxy surveys.

A new generation of interferometric instruments is emerging which aim to use intensity mapping of redshifted $21\,$cm radiation to measure the large-scale structure of the Universe at $z\simeq 1-2$ over wide areas of sky. While these instruments typically have limited angular resolution, they cover huge volumes and thus can be used to provide large samples of rare objects. In this paper we study how well such instruments could find spatially extended large-scale structures, such as cosmic voids, using a matched filter formalism. Such a formalism allows us to work in Fourier space, the natural space for interferometers, and to study the impact of finite $u-v$ coverage, noise and foregrounds on our ability to recover voids. We find that in the absence of foregrounds such instruments would provide enormous catalogs of voids, with high completeness, but that control of foregrounds is key to realizing this goal.

Galaxy clustering on small scales is significantly under-predicted by sub-halo abundance matching (SHAM) models that populate (sub-)haloes with galaxies based on peak halo mass, $M_{\rm peak}$. SHAM models based on the peak maximum circular velocity, $V_{\rm peak}$, have had much better success. The primary reason $M_{\rm peak}$ based models fail is the relatively low abundance of satellite galaxies produced in these models compared to those based on $V_{\rm peak}$. Despite success in predicting clustering, a simple $V_{\rm peak}$ based SHAM model results in predictions for galaxy growth that are at odds with observations. We evaluate three possible remedies that could "save" mass-based SHAM: (1) SHAM models require a significant population of "orphan" galaxies as a result of artificial disruption/merging of sub-haloes in modern high resolution dark matter simulations; (2) satellites must grow significantly after their accretion; and (3) stellar mass is significantly affected by halo assembly history. No solution is entirely satisfactory. However, regardless of the particulars, we show that popular SHAM models based on $M_{\rm peak}$ cannot be complete physical models as presented. Either $V_{\rm peak}$ truly is a better predictor of stellar mass at $z\sim 0$ and it remains to be seen how the correlation between stellar mass and $V_{\rm peak}$ comes about, or SHAM models are missing vital component(s) that significantly affect galaxy clustering.

We present a calibration of halo assembly bias using the Separate Universe technique. Specifically, we measure the response of halo abundances at fixed mass and concentration to the presence of an infinite-wavelength initial perturbation. We develop an analytical framework for describing the concentration dependence of this peak-background split halo bias -- a measure of assembly bias -- relying on the near-Lognormal distribution of halo concentration at fixed halo mass. The combination of this analytical framework and the Separate Universe technique allows us to achieve very high precision in the calibration of the linear assembly bias $b_1$, and qualitatively reproduces known trends such as the monotonic decrease (increase) of $b_1$ with halo concentration at large (small) masses. The same framework extends to the concentration dependence of higher order bias parameters $b_n$, and we present the first calibration of assembly bias in $b_2$. Our calibrations are directly applicable in analytical Halo Model calculations that seek to robustly detect galaxy assembly bias in observational samples. We detect a non-universality in the $b_1 - b_2$ relation arising from assembly bias, and suggest that simultaneous measurements of these bias parameters could be used to both detect the signature of assembly bias as well as mitigate its effects in cosmological analyses.

We present a simple strategy to mitigate the impact of an incomplete spectroscopic redshift galaxy sample as a result of fiber assignment and survey tiling. The method has been designed for the Dark Energy Spectroscopic Instrument (DESI) galaxy survey but may have applications beyond this. We propose a modification to the usual correlation function that nulls the almost purely angular modes affected by survey incompleteness due to fiber assignment. Predictions of this modified statistic can be calculated given a model of the two point correlation function. The new statistic can be computed with a slight modification to the data catalogues input to the standard correlation function code and does not incur any additional computational time. Finally we show that the spherically averaged baryon acoustic oscillation signal is not biased by the new statistic.

We present a Bayesian phase-space reconstruction of the cosmic large-scale matter density and velocity fields from the SDSS-III Baryon Oscillations Spectroscopic Survey Data Release 12 (BOSS DR12) CMASS galaxy clustering catalogue. We rely on a given $\Lambda$CDM cosmology, a mesh resolution in the range of 6-10 $h^{-1}$ Mpc, and a lognormal-Poisson model with a redshift dependent nonlinear bias. The bias parameters are derived from the data and a general renormalised perturbation theory approach. We use combined Gibbs and Hamiltonian sampling, implemented in the \textsc{argo} code, to iteratively reconstruct the dark matter density field and the coherent peculiar velocities of individual galaxies, correcting hereby for coherent redshift space distortions (RSD). Our tests relying on accurate $N$-body based mock galaxy catalogues, show unbiased real space power spectra of the nonlinear density field up to $k\sim0.2\, h$ Mpc$^{-1}$, and vanishing quadrupoles down to $r\sim20\,h^{-1}$ Mpc. We also demonstrate that the nonlinear cosmic web can be obtained from the tidal field tensor based on the Gaussian component of the reconstructed density field. We find that the reconstructed velocities have a statistical correlation coefficient compared to the true velocities of each individual lightcone mock galaxy of $r\sim0.68$ including about 10% of satellite galaxies with virial motions (about $r=0.75$ without satellites). The power spectra of the velocity divergence agree well with theoretical predictions up to $k\sim0.2\,h\,{\rm Mpc}^{-1}$. This work will be especially useful to improve, e.g. BAO reconstructions, kinematic Sunyaev-Zeldovich (kSZ), integrated Sachs-Wolfe (ISW) measurements, or environmental studies.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

We present cosmological results from the final galaxy clustering data set of the Baryon Oscillation Spectroscopic Survey, part of the Sloan Digital Sky Survey III. Our combined galaxy sample comprises 1.2 million massive galaxies over an effective area of 9329 deg^2 and volume of 18.7 Gpc^3, divided into three partially overlapping redshift slices centred at effective redshifts 0.38, 0.51, and 0.61. We measure the angular diameter distance DM and Hubble parameter H from the baryon acoustic oscillation (BAO) method after applying reconstruction to reduce non-linear effects on the BAO feature. Using the anisotropic clustering of the pre-reconstruction density field, we measure the product DM*H from the Alcock-Paczynski (AP) effect and the growth of structure, quantified by f{\sigma}8(z), from redshift-space distortions (RSD). We combine measurements presented in seven companion papers into a set of consensus values and likelihoods, obtaining constraints that are tighter and more robust than those from any one method. Combined with Planck 2015 cosmic microwave background measurements, our distance scale measurements simultaneously imply curvature {\Omega}_K =0.0003+/-0.0026 and a dark energy equation of state parameter w = -1.01+/-0.06, in strong affirmation of the spatially flat cold dark matter model with a cosmological constant ({\Lambda}CDM). Our RSD measurements of f{\sigma}_8, at 6 per cent precision, are similarly consistent with this model. When combined with supernova Ia data, we find H0 = 67.3+/-1.0 km/s/Mpc even for our most general dark energy model, in tension with some direct measurements. Adding extra relativistic species as a degree of freedom loosens the constraint only slightly, to H0 = 67.8+/-1.2 km/s/Mpc. Assuming flat {\Lambda}CDM we find {\Omega}_m = 0.310+/-0.005 and H0 = 67.6+/-0.5 km/s/Mpc, and we find a 95% upper limit of 0.16 eV/c^2 on the neutrino mass sum.

Large redshift surveys capable of measuring the Baryon Acoustic Oscillation (BAO) signal have proven to be an effective way of measuring the distance-redshift relation in cosmology. Building off the work in Zhu et al. (2015), we develop a technique to directly constrain the distance-redshift relation from BAO measurements without splitting the sample into redshift bins. We parametrize the distance-redshift relation, relative to a fiducial model, as a quadratic expansion. We measure its coefficients and reconstruct the distance-redshift relation from the expansion. We apply the redshift weighting technique in Zhu et al. (2015) to the clustering of galaxies from 1000 QuickPM (QPM) mock simulations after reconstruction and achieve a 0.75% measurement of the angular diameter distance $D_A$ at $z=0.64$ and the same precision for Hubble parameter H at $z=0.29$. These QPM mock catalogs are designed to mimic the clustering and noise level of the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12 (DR12). We compress the correlation functions in the redshift direction onto a set of weighted correlation functions. These estimators give unbiased $D_A$ and $H$ measurements at all redshifts within the range of the combined sample. We demonstrate the effectiveness of redshift weighting in improving the distance and Hubble parameter estimates. Instead of measuring at a single 'effective' redshift as in traditional analyses, we report our $D_A$ and $H$ measurements at all redshifts. The measured fractional error of $D_A$ ranges from 1.53% at $z=0.2$ to 0.75% at $z=0.64$. The fractional error of $H$ ranges from 0.75% at $z=0.29$ to 2.45% at $z = 0.7$. Our measurements are consistent with a Fisher forecast to within 10% to 20% depending on the pivot redshift. We further show the results are robust against the choice of fiducial cosmologies, galaxy bias models, and RSD streaming parameters.

This draft report summarizes and details the findings, results, and recommendations derived from the ASCR/HEP Exascale Requirements Review meeting held in June, 2015. The main conclusions are as follows. 1) Larger, more capable computing and data facilities are needed to support HEP science goals in all three frontiers: Energy, Intensity, and Cosmic. The expected scale of the demand at the 2025 timescale is at least two orders of magnitude -- and in some cases greater -- than that available currently. 2) The growth rate of data produced by simulations is overwhelming the current ability, of both facilities and researchers, to store and analyze it. Additional resources and new techniques for data analysis are urgently needed. 3) Data rates and volumes from HEP experimental facilities are also straining the ability to store and analyze large and complex data volumes. Appropriately configured leadership-class facilities can play a transformational role in enabling scientific discovery from these datasets. 4) A close integration of HPC simulation and data analysis will aid greatly in interpreting results from HEP experiments. Such an integration will minimize data movement and facilitate interdependent workflows. 5) Long-range planning between HEP and ASCR will be required to meet HEP's research needs. To best use ASCR HPC resources the experimental HEP program needs a) an established long-term plan for access to ASCR computational and data resources, b) an ability to map workflows onto HPC resources, c) the ability for ASCR facilities to accommodate workflows run by collaborations that can have thousands of individual members, d) to transition codes to the next-generation HPC platforms that will be available at ASCR facilities, e) to build up and train a workforce capable of developing and using simulations and analysis to support HEP scientific research on next-generation systems.

Precision measurements of the large scale structure of the Universe require large numbers of high fidelity mock catalogs to accurately assess, and account for, the presence of systematic effects. We introduce and test a scheme for generating mock catalogs rapidly using suitably derated N-body simulations. Our aim is to reproduce the large scale structure and the gross properties of dark matter halos with high accuracy, while sacrificing the details of the internal structure of the halos. By adjusting global and local time-steps in an N-body code, we demonstrate that we recover halo masses to better than 0.5% and the power spectrum to better than 1% both in real and redshift space for k = 1h/Mpc, while requiring a factor of 4 less CPU time. We also calibrate the redshift spacing of outputs required to generate simulated light cones. We find that outputs separated by every z = 0.05 allow us to interpolate particle positions and velocities to reproduce the real and redshift space power spectra to better than 1% (out to k = 1h/Mpc). We apply these ideas to generate a suite of simulations spanning a range of cosmologies, motivated by the Baryon Oscillation Spectroscopic Survey (BOSS) but broadly applicable to future large scale structure surveys including eBOSS and DESI. As an initial demonstration of the utility of such simulations, we calibrate the shift in the baryonic acoustic oscillation peak position as a function of galaxy bias with higher precision than has been possible so far. This paper also serves to document the simulations, which we make publicly available.

We present distance scale measurements from the baryon acoustic oscillation signal in the CMASS and LOWZ samples from the Data Release 12 of the Baryon Oscillation Spectroscopic Survey (BOSS). The total volume probed is 14.5 Gpc$^3$, a 10 per cent increment from Data Release 11. From an analysis of the spherically averaged correlation function, we infer a distance to $z=0.57$ of $D_V(z)r^{\rm fid}_{\rm d}/r_ {\rm d}=2028\pm21$ Mpc and a distance to $z=0.32$ of $D_V(z)r^{\rm fid}_{\rm d}/r_{\rm d}=1264\pm22$ Mpc assuming a cosmology in which $r^{\rm fid}_{\rm d}=147.10$ Mpc. From the anisotropic analysis, we find an angular diameter distance to $z=0.57$ of $D_{\rm A}(z)r^{\rm fid}_{\rm d}/r_{\rm d}=1401\pm21$ Mpc and a distance to $z=0.32$ of $981\pm20$ Mpc, a 1.5 per cent and 2.0 per cent measurement respectively. The Hubble parameter at $z=0.57$ is $H(z)r_{\rm d}/r^{\rm fid}_{\rm d}=100.3\pm3.7$ km s$^{-1}$ Mpc$^{-1}$ and its value at $z=0.32$ is $79.2\pm5.6$ km s$^{-1}$ Mpc$^{-1}$, a 3.7 per cent and 7.1 per cent measurement respectively. These cosmic distance scale constraints are in excellent agreement with a $\Lambda$CDM model with cosmological parameters released by the recent Planck 2015 results.

Acoustic waves traveling through the early Universe imprint a characteristic scale in the clustering of galaxies, QSOs and inter-galactic gas. This scale can be used as a standard ruler to map the expansion history of the Universe, a technique known as Baryon Acoustic Oscillations (BAO). BAO offer a high-precision, low-systematics means of constraining our cosmological model. The statistical power of BAO measurements can be improved if the `smearing' of the acoustic feature by non-linear structure formation is undone in a process known as reconstruction. In this paper we use low-order Lagrangian perturbation theory to study the ability of $21\,$cm experiments to perform reconstruction and how augmenting these surveys with galaxy redshift surveys at relatively low number densities can improve performance. We find that the critical number density which must be achieved in order to benefit $21\,$cm surveys is set by the linear theory power spectrum near its peak, and corresponds to densities achievable by upcoming surveys of emission line galaxies such as eBOSS and DESI. As part of this work we analyze reconstruction within the framework of Lagrangian perturbation theory with local Lagrangian bias, redshift-space distortions, ${\bf k}$-dependent noise and anisotropic filtering schemes.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

The final step of most large-scale structure analyses involves the comparison of power spectra or correlation functions to theoretical models. It is clear that the theoretical models have parameter dependence, but frequently the measurements and the covariance matrix depend upon some of the parameters as well. We show that a very simple interpolation scheme from an unstructured mesh allows for an efficient way to include this parameter dependence self-consistently in the analysis at modest computational expense. We describe two schemes for covariance matrices. The scheme which uses the geometric structure of such matrices performs roughly twice as well as the simplest scheme, though both perform very well.

We apply a method recently introduced to the statistical literature to directly estimate the precision matrix from an ensemble of samples drawn from a corresponding Gaussian distribution. Motivated by the observation that cosmological precision matrices are often approximately sparse, the method allows one to exploit this sparsity of the precision matrix to more quickly converge to an asymptotic 1/sqrt(Nsim) rate while simultaneously providing an error model for all of the terms. Such an estimate can be used as the starting point for further regularization efforts which can improve upon the 1/sqrt(Nsim) limit above, and incorporating such additional steps is straightforward within this framework. We demonstrate the technique with toy models and with an example motivated by large-scale structure two-point analysis, showing significant improvements in the rate of convergence.For the large-scale structure example we find errors on the precision matrix which are factors of 5 smaller than for the sample precision matrix for thousands of simulations or, alternatively, convergence to the same error level with more than an order of magnitude fewer simulations.

We present a recalibration of the Sloan Digital Sky Survey (SDSS) photometry with new flat fields and zero points derived from Pan-STARRS1 (PS1). Using PSF photometry of 60 million stars with $16 < r < 20$, we derive a model of amplifier gain and flat-field corrections with per-run RMS residuals of 3 millimagnitudes (mmag) in $griz$ bands and 15 mmag in $u$ band. The new photometric zero points are adjusted to leave the median in the Galactic North unchanged for compatibility with previous SDSS work. We also identify transient non-photometric periods in SDSS ("contrails") based on photometric deviations co-temporal in SDSS bands. The recalibrated stellar PSF photometry of SDSS and PS1 has an RMS difference of {9,7,7,8} mmag in $griz$, respectively, when averaged over $15'$ regions.

The Baryon Oscillation Spectroscopic Survey (BOSS), part of the Sloan Digital Sky Survey (SDSS) III project, has provided the largest survey of galaxy redshifts available to date, in terms of both the number of galaxy redshifts measured by a single survey, and the effective cosmological volume covered. Key to analysing the clustering of these data to provide cosmological measurements is understanding the detailed properties of this sample. Potential issues include variations in the target catalogue caused by changes either in the targeting algorithm or properties of the data used, the pattern of spectroscopic observations, the spatial distribution of targets for which redshifts were not obtained, and variations in the target sky density due to observational systematics. We document here the target selection algorithms used to create the galaxy samples that comprise BOSS. We also present the algorithms used to create large scale structure catalogues for the final Data Release (DR12) samples and the associated random catalogues that quantify the survey mask. The algorithms are an evolution of those used by the BOSS team to construct catalogues from earlier data, and have been designed to accurately quantify the galaxy sample. The code used, designated MKSAMPLE, is released with this paper.

We derive constraints on cosmological parameters and tests of dark energy models from the combination of baryon acoustic oscillation (BAO) measurements with cosmic microwave background (CMB) and Type Ia supernova (SN) data. We take advantage of high-precision BAO measurements from galaxy clustering and the Ly-alpha forest (LyaF) in the BOSS survey of SDSS-III. BAO data alone yield a high confidence detection of dark energy, and in combination with the CMB angular acoustic scale they further imply a nearly flat universe. Combining BAO and SN data into an "inverse distance ladder" yields a 1.7% measurement of $H_0=67.3 \pm1.1$ km/s/Mpc. This measurement assumes standard pre-recombination physics but is insensitive to assumptions about dark energy or space curvature, so agreement with CMB-based estimates that assume a flat LCDM cosmology is an important corroboration of this minimal cosmological model. For open LCDM, our BAO+SN+CMB combination yields $\Omega_m=0.301 \pm 0.008$ and curvature $\Omega_k=-0.003 \pm 0.003$. When we allow more general forms of evolving dark energy, the BAO+SN+CMB parameter constraints remain consistent with flat LCDM. While the overall $\chi^2$ of model fits is satisfactory, the LyaF BAO measurements are in moderate (2-2.5 sigma) tension with model predictions. Models with early dark energy that tracks the dominant energy component at high redshifts remain consistent with our constraints. Expansion history alone yields an upper limit of 0.56 eV on the summed mass of neutrino species, improving to 0.26 eV if we include Planck CMB lensing. Standard dark energy models constrained by our data predict a level of matter clustering that is high compared to most, but not all, observational estimates. (Abridged)

We introduce a new method for estimating the covariance matrix for the galaxy correlation function in surveys of large-scale structure. Our method combines simple theoretical results with a realistic characterization of the survey to dramatically reduce noise in the covariance matrix. For example, with an investment of only ~1,000 CPU hours we can produce a model covariance matrix with noise levels that would otherwise require ~35,000 mocks. Non-Gaussian contributions to the model are calibrated against mock catalogs, after which the model covariance is found to be in impressive agreement with the mock covariance matrix. Since calibration of this method requires fewer mocks than brute force approaches, we believe that it could dramatically reduce the number of mocks required to analyse future surveys.

The two-point clustering of dark matter halos is influenced by halo properties besides mass, a phenomenon referred to as halo assembly bias. Using the depth of the gravitational potential well, $V_{\rm max}$, as our secondary halo property, in this paper we present the first study of the scale-dependence assembly bias. In the large-scale linear regime, $r\geq10h^{-1}{\rm Mpc},$ our findings are in keeping with previous results. In particular, at the low-mass end ($M_{\rm vir}<M_{\rm coll}\approx10^{12.5}{\rm M}_{\odot}$), halos with high-$V_{\rm max}$ show stronger large-scale clustering relative to halos with low-$V_{\rm max}$ of the same mass, this trend weakens and reverses for $M_{\rm vir}\geq M_{\rm coll}.$ In the nonlinear regime, assembly bias in low-mass halos exhibits a pronounced scale-dependent "bump" at $500h^{-1}{\rm kpc}-5h^{-1}{\rm Mpc},$ a new result. This feature weakens and eventually vanishes for halos of higher mass. We show that this scale-dependent signature can primarily be attributed to a special subpopulation of ejected halos, defined as present-day host halos that were previously members of a higher-mass halo at some point in their past history. A corollary of our results is that galaxy clustering on scales of $r\sim1-2h^{-1}{\rm Mpc}$ can be impacted by up to $\sim15\%$ by the choice of the halo property used in the halo model, even for stellar mass-limited samples.

The halo occupation distribution (HOD) approach has proven to be an effective method for modeling galaxy clustering and bias. In this approach, galaxies of a given type are probabilistically assigned to individual halos in N-body simulations. In this paper, we present a fast emulator for predicting the fully nonlinear galaxy power spectrum over a range of freely specifiable HOD modeling parameters. The emulator is constructed using results from 100 HOD models run on a large LCDM N-body simulation, with Gaussian Process interpolation applied to a PCA-based representation of the galaxy power spectrum. The total error is currently ~3% (~2% in the simulation and ~1% in the emulation process) from z=1 to z=0, over the considered parameter range. We use the emulator to investigate parametric dependencies in the HOD model, as well as the behavior of galaxy bias as a function of HOD parameters. The emulator is publicly available at http://www.hep.anl.gov/cosmology/CosmicEmu/emu.html.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

We investigate the ability of current implementations of galaxy group finders to recover colour-dependent halo occupation statistics. To test the fidelity of group catalogue inferred statistics, we run three different group finders used in the literature over a mock that includes galaxy colours in a realistic manner. Overall, the resulting mock group catalogues are remarkably similar, and most colour-dependent statistics are recovered with reasonable accuracy. However, it is also clear that certain systematic errors arise as a consequence of correlated errors in group membership determination, central/satellite designation, and halo mass assignment. We introduce a new statistic, the halo transition probability (HTP), which captures the combined impact of all these errors. As a rule of thumb, errors tend to equalize the properties of distinct galaxy populations (i.e. red vs. blue galaxies or centrals vs. satellites), and to result in inferred occupation statistics that are more accurate for red galaxies than for blue galaxies. A statistic that is particularly poorly recovered from the group catalogues is the red fraction of central galaxies as function of halo mass. Group finders do a good job in recovering galactic conformity, but also have a tendency to introduce weak conformity when none is present. We conclude that proper inference of colour-dependent statistics from group catalogues is best achieved using forward modelling (i.e., running group finders over mock data), or by implementing a correction scheme based on the HTP, as long as the latter is not too strongly model-dependent.

The Sloan Digital Sky Survey has surveyed 14,555 square degrees of the sky, and delivered over a trillion pixels of imaging data. We present the large-scale clustering of 1.6 million quasars between z = 0.5 and z = 2.5 that have been classified from this imaging, representing the highest density of quasars ever studied for clustering measurements. This data set spans ~11,000 square degrees and probes a volume of 80(Gpc/h)^3. In principle, such a large volume and medium density of tracers should facilitate high-precision cosmological constraints. We measure the angular clustering of photometrically classified quasars using an optimal quadratic estimator in four redshift slices with an accuracy of ~25% over a bin width of l ~10 - 15 on scales corresponding to matter-radiation equality and larger (l ~ 2 - 30). Observational systematics can strongly bias clustering measurements on large scales, which can mimic cosmologically relevant signals such as deviations from Gaussianity in the spectrum of primordial perturbations. We account for systematics by employing a new method recently proposed by Agarwal et al. (2014) to the clustering of photometrically classified quasars. We carefully apply our methodology to mitigate known observational systematics and further remove angular bins that are contaminated by unknown systematics. Combining quasar data with the photometric luminous red galaxy (LRG) sample of Ross et al. (2011) and Ho et al. (2012), and marginalizing over all bias and shot noise-like parameters, we obtain a constraint on local primordial non-Gaussianity of fNL = -113+/-154 (1\sigma error). [Abridged]

Future baryon acoustic oscillation (BAO) surveys will survey very large volumes, covering wide ranges in redshift. We derive a set of redshift weights to compress the information in the redshift direction to a small number of modes. We suggest that such a compression preserves almost all of the signal for most cosmologies, while giving high signal-to-noise measurements for each combination. We present some toy models and simple worked examples. As an intermediate step, we give a precise meaning to the "effective redshift" of a BAO measurement.

We present a detailed study of how dark matter haloes assemble their mass and grow their (central) potential well. We characterize these via their mass accretion histories (MAHs) and potential well growth histories (PWGHs), which we extract from the Bolshoi simulation and from semi-analytical merger trees supplemented with a method to compute the maximum circular velocity, Vmax, of progenitor haloes. The results of both methods are in excellent agreement, both in terms of the average and the scatter. We show that the MAH and PWGH are tightly correlated, and that growth of the central potential precedes the assembly of mass; the maximum circular velocity is already half the present day value by the time the halo has accreted only 2 percent of its final mass. Finally, we demonstrate that MAHs have a universal form, which we use to develop a new and improved universal model that can be used to compute the average or median MAH and PWGH for a halo of any mass in any LCDM cosmology, without having to run a numerical simulation or a set of halo merger trees.

We present significant improvements in cosmic distance measurements from the WiggleZ Dark Energy Survey, achieved by applying the reconstruction of the baryonic acoustic feature technique. We show using both data and simulations that the reconstruction technique can often be effective despite patchiness of the survey, significant edge effects and shot-noise. We investigate three redshift bins in the redshift range 0.2<$z$<1, and in all three find improvement after reconstruction in the detection of the baryonic acoustic feature and its usage as a standard ruler. We measure model independent distance measures $D_{\mathrm V}(r_{\mathrm s}^\mathrm{fid}/r_{\mathrm s})$ of 1716 $\pm$ 83 Mpc, 2221 $\pm$ 101 Mpc, 2516 $\pm$ 86 Mpc (68% CL) at effective redshifts z = 0.44, 0.6, 0.73, respectively, where $D_{\mathrm V}$ is the volume-average-distance, and $r_{\mathrm s}$ is the sound horizon at the end of the baryon drag epoch. These significantly improved 4.8, 4.5 and 3.4 percent accuracy measurements are equivalent to those expected from surveys with up to 2.5 times the volume of WiggleZ. These measurements are fully consistent with cosmologies allowed by the analyses of the Planck Collaboration and the Sloan Digital Sky Survey.We provide the $D_{\mathrm V}(r_{\mathrm s}^\mathrm{fid}/r_{\mathrm s})$ posterior probability distributions and their covariances. When combining these measurements with temperature fluctuations measurements of Planck, the polarization of WMAP9, and the 6dF Galaxy Survey baryonic acoustic feature, we do not detect deviations from a flat LCDM model. Assuming this model we constrain the current expansion rate to $H_0$ = 67.15 $\pm$ 0.98 kms$^{-1}$Mpc$^{-1}$. Allowing the equation of state of dark energy to vary we obtain $w_\mathrm{DE}$ = -1.080 $\pm$ 0.135. When assuming a curved LCDM model we obtain a curvature value of $\Omega_{\mathrm K}$ = -0.0043 $\pm$ 0.0047.

This document presents the results from the Distances subgroup of the Cosmic Frontier Community Planning Study (Snowmass 2013). We summarize the current state of the field as well as future prospects and challenges. In addition to the established probes using Type IA supernovae and baryon acoustic oscillations, we also consider prospective methods based on clusters, active galactic nuclei, gravitational wave sirens and strong lensing time delays.

We present a one per cent measurement of the cosmic distance scale from the detections of the baryon acoustic oscillations in the clustering of galaxies from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). Our results come from the Data Release 11 (DR11) sample, containing nearly one million galaxies and covering approximately $8\,500$ square degrees and the redshift range $0.2<z<0.7$. We also compare these results with those from the publicly released DR9 and DR10 samples. Assuming a concordance $\Lambda$CDM cosmological model, the DR11 sample covers a volume of 13\,Gpc${}^3$ and is the largest region of the Universe ever surveyed at this density. We measure the correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of over $7\,\sigma$ in both the correlation function and power spectrum. Fitting for the position of the acoustic features measures the distance relative to the sound horizon at the drag epoch, $r_d$, which has a value of $r_{d,{\rm fid}}=149.28\,$Mpc in our fiducial cosmology. We find $D_V=(1264\pm25\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.32$ and $D_V=(2056\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.57$. At 1.0 per cent, this latter measure is the most precise distance constraint ever obtained from a galaxy survey. Separating the clustering along and transverse to the line-of-sight yields measurements at $z=0.57$ of $D_A=(1421\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ and $H=(96.8\pm3.4\,{\rm km/s/Mpc})(r_{d,{\rm fid}}/r_d)$. Our measurements of the distance scale are in good agreement with previous BAO measurements and with the predictions from cosmic microwave background data for a spatially flat cold dark matter model with a cosmological constant.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

We present improved methodology for including covariance matrices in the error budget of Baryon Oscillation Spectroscopic Survey (BOSS) galaxy clustering measurements, revisiting Data Release 9 (DR9) analyses, and describing a method that is used in DR10/11 analyses presented in companion papers. The precise analysis method adopted is becoming increasingly important, due to the precision that BOSS can now reach: even using as many as 600 mock catalogues to estimate covariance of 2-point clustering measurements can still lead to an increase in the errors of ~20%, depending on how the cosmological parameters of interest are measured. In this paper we extend previous work on this contribution to the error budget, deriving formulae for errors measured by integrating over the likelihood, and to the distribution of recovered best-fit parameters fitting the simulations also used to estimate the covariance matrix. Both are situations that previous analyses of BOSS have considered. We apply the formulae derived to Baryon Acoustic Oscillation (BAO) and Redshift-Space Distortion (RSD) measurements from BOSS in our companion papers. To further aid these analyses, we consider the optimum number of bins to use for 2-point measurements using the monopole power spectrum or correlation function for BAO, and the monopole and quadrupole moments of the correlation function for anisotropic-BAO and RSD measurements.

The American Physical Society's Division of Particles and Fields initiated a long-term planning exercise over 2012-13, with the goal of developing the community's long term aspirations. The sub-group "Dark Energy and CMB" prepared a series of papers explaining and highlighting the physics that will be studied with large galaxy surveys and cosmic microwave background experiments. This paper summarizes the findings of the other papers, all of which have been submitted jointly to the arXiv.

We present results on the clustering of 282,068 galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) sample of massive galaxies with redshifts 0.4<z<0.7 which is part of the Sloan Digital Sky Survey III project. Our results cover a large range of scales from ~0.5 to ~90 Mpc/h. We compare these estimates with the expectations of the flat LCDM cosmological model with parameters compatible with WMAP7 data. We use the MultiDark cosmological simulation together with a simple halo abundance matching technique, to estimate galaxy correlation functions, power spectra, abundance of subhaloes and galaxy biases. We find that the LCDM model gives a reasonable description to the observed correlation functions at z~0.5, which is a remarkably good agreement considering that the model, once matched to the observed abundance of BOSS galaxies, does not have any free parameters. However, we find a deviation (>~10%) in the correlation functions for scales less than ~1 Mpc/h and ~10-40 Mpc/h. A more realistic abundance matching model and better statistics from upcoming observations are needed to clarify the situation. We also estimate that about 12% of the "galaxies" in the abundance-matched sample are satellites inhabiting central haloes with mass M>~1e14 M_sun/h. Using the MultiDark simulation we also study the real space halo bias b(r) of the matched catalogue finding that b=2.00+/-0.07 at large scales, consistent with the one obtained using the measured BOSS projected correlation function. Furthermore, the linear large-scale bias depends on the number density n of the abundance-matched sample as b=-0.048-(0.594+/-0.02)*log(n/(h/Mpc)^3). Extrapolating these results to BAO scales we measure a scale-dependent damping of the acoustic signal produced by non-linear evolution that leads to ~2-4% dips at ~3 sigma level for wavenumbers k>~0.1 h/Mpc in the linear large-scale bias.

We measure the luminosity and color dependence and the redshift evolution of galaxy clustering in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Ninth Data Release. We focus on the projected two-point correlation function (2PCF) of subsets of its CMASS sample, which includes about 260,000 galaxies over ~3,300 sq. deg in the redshift range 0.43<z<0.7. To minimize the selection effect on galaxy clustering, we construct well-defined luminosity and color subsamples by carefully accounting for the CMASS galaxy selection cuts. The 2PCF of the whole CMASS sample, if approximated by a power-law, has a correlation length of r_0=7.93\pm0.06Mpc/h and an index of \gamma=1.85\pm0.01. Clear dependences on galaxy luminosity and color are found for the projected 2PCF in all redshift bins, with more luminous and redder galaxies generally exhibiting stronger clustering and steeper 2PCF. The color dependence is also clearly seen for galaxies within the red sequence, consistent with the behavior of SDSS-II main sample galaxies at lower redshifts. At a given luminosity (k+e corrected), no significant evolution of the projected 2PCFs with redshift is detected for red sequence galaxies. We also construct galaxy samples of fixed number density at different redshifts, using redshift-dependent magnitude thresholds. The clustering of these galaxies in the CMASS redshift range is found to be consistent with that predicted by passive evolution. Our measurements of the luminosity and color dependence and redshift evolution of galaxy clustering will allow for detailed modeling of the relation between galaxies and dark matter halos and new constraints on galaxy formation and evolution.

We present measurements of the angular diameter distance to and Hubble parameter at z=0.57 from the measurement of the baryon acoustic peak in the correlation of galaxies from the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey. Our analysis is based on a sample from Data Release 9 of 264,283 galaxies over 3275 square degrees in the redshift range 0.43<z<0.70. We use two different methods to provide robust measurement of the acoustic peak position across and along the line of sight in order to measure the cosmological distance scale. We find D_A(0.57) = 1408 +/- 45 Mpc and H(0.57) = 92.9 +/- 7.8 km/s/Mpc for our fiducial value of the sound horizon. These results from the anisotropic fitting are fully consistent with the analysis of the spherically averaged acoustic peak position presented in Anderson et al, 2012. Our distance measurements are a close match to the predictions of the standard cosmological model featuring a cosmological constant and zero spatial curvature.

We analyze the 2D correlation function of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) CMASS sample of massive galaxies of the ninth data release to measure cosmic expansion H and the angular diameter distance D_A at a mean redshift of <z> = 0.57. We apply, for the first time, a new correlation function technique called clustering wedges. Using a physically motivated model, the anisotropic baryonic acoustic feature in the galaxy sample is detected at a significance level of 4.7 sigma compared to a featureless model. The baryonic acoustic feature is used to obtain model independent constraints cz/H/r_s = 12.28 +- 0.82 (6.7 per-cent accuracy) and D_A/r_s = 9.05 +- 0.27 (3.0 per-cent) with a correlation coefficient of -0.5, where r_s is the sound horizon scale at the end of the baryonic drag era. We conduct thorough tests on the data and 600 simulated realizations, finding robustness of the results regardless of the details of the analysis method. Combining with r_s constraints from the Cosmic Microwave Background we obtain H(0.57) = 90.8 +- 6.2 kms-1Mpc-1 and D_A(0.57) = 1386 +- 45 Mpc. We use simulations to forecast results of the final BOSS CMASS data set. We apply the reconstruction technique on the simulations demonstrating that the sharpening of the anisotropic baryonic acoustic feature should improve the detection as well as tighten constraints of H and D_A by 30 per-cent on average.

We explore the cosmological implications of the clustering wedges, xi_perp(s) and xi_para(s), of the CMASS Data Release 9 (DR9) sample of the Baryon Oscillation Spectroscopic Survey (BOSS). These clustering wedges are defined by averaging the full two-dimensional correlation function, xi(mu,s), over the ranges 0<mu<0.5 and 0.5<mu<1, respectively. These measurements allow us to constrain the parameter combinations D_A(z)/r_s(z_d)=9.03 +- 0.21 and cz/(r_s(z_d)H(z)) = 12.14 +- 0.43 at the mean redsfhit of the sample, z=0.57. We combine the information from the clustering wedges with recent measurements of CMB, BAO and type Ia supernovae to obtain constraints on the cosmological parameters of the standard LCDM model and a number of potential extensions. The information encoded in the clustering wedges is most useful when the dark energy equation of state is allowed to deviate from its standard LCDM value. The combination of all datasets shows no evidence of a deviation from a constant dark energy equation of state, in which case we find w_DE = -1.013 +- 0.064, in complete agreement with a cosmological constant. We explore potential deviations from general relativity by constraining the growth rate f(z)=d ln D(a)/ d ln a, in which case the combination of the CMASS clustering wedges with CMB data implies f(z=0.57)=0.719 +- 0.094, in accordance with the predictions of GR. Our results clearly illustrate the additional constraining power of anisotropic clustering measurements with respect to that of angle-averaged quantities.

We present measurements of the angular diameter distance D_A(z) and the Hubble parameter H(z) at z=0.35 using the anisotropy of the baryon acoustic oscillation (BAO) signal measured in the galaxy clustering distribution of the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) Luminous Red Galaxies (LRG) sample. Our work is the first to apply density-field reconstruction to an anisotropic analysis of the acoustic peak. Reconstruction partially removes the effects of non-linear evolution and redshift-space distortions in order to sharpen the acoustic signal. We present the theoretical framework behind the anisotropic BAO signal and give a detailed account of the fitting model we use to extract this signal from the data. Our method focuses only on the acoustic peak anisotropy, rather than the more model-dependent anisotropic information from the broadband power. We test the robustness of our analysis methods on 160 LasDamas DR7 mock catalogues and find that our models are unbiased at the ~0.2% level in measuring the BAO anisotropy. After reconstruction we measure D_A(z=0.35)=1050+/-38 Mpc and H(z=0.35)=84.4+/-7.0 km/s/Mpc assuming a sound horizon of r_s=152.76 Mpc. Note that these measurements are correlated with a correlation coefficient of 0.58. This represents a factor of 1.4 improvement in the error on D_A relative to the pre-reconstruction case; a factor of 1.2 improvement is seen for H.

We present the measurement of the two-point cross-correlation function (CCF) of 8,198 Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) quasars and 349,608 DR10 CMASS galaxies from the Baryonic Oscillation Spectroscopic Survey (BOSS) at redshift <z>~0.5 (0.3<z<0.9). The cross-correlation function can be reasonably well fit by a power-law model xi_QG(r)=(r/r_0)^(-gamma) on projected scales of r_p=2-25 Mpc/h with r_0=6.61+-0.25 Mpc/h and gamma=1.69+-0.07. We estimate a quasar linear bias of b_Q=1.38+-0.10 at <z>=0.53 from the CCF measurements. This linear bias corresponds to a characteristic host halo mass of ~4x10^12 M_sun/h, compared to ~10^13 M_sun/h characteristic host halo mass for CMASS galaxies. We divide the quasar sample in luminosity and constrain the luminosity dependence of quasar bias to be db_Q/dlogL=0.20+-0.34 or 0.11+-0.32 (depending on different luminosity divisions) for quasar luminosities -23.5>M_i(z=2)>-25.5, implying a weak luminosity dependence of quasar clustering for the bright end of the quasar population at <z>~0.5. We compare our measurements with theoretical predictions, Halo Occupation Distribution (HOD) models and mock catalogs. These comparisons suggest quasars reside in a broad range of host halos, and the host halo mass distributions significantly overlap with each other for quasars at different luminosities, implying a poor correlation between halo mass and instantaneous quasar luminosity. We also find that the quasar HOD parameterization is largely degenerate such that different HODs can reproduce the CCF equally well, but with different outcomes such as the satellite fraction and host halo mass distribution. These results highlight the limitations and ambiguities in modeling the distribution of quasars with the standard HOD approach and the need for additional information in populating quasars in dark matter halos with HOD. [Abridged]

We report on the small scale (0.5<r<40h^-1 Mpc) clustering of 78895 massive (M*~10^11.3M_sun) galaxies at 0.2<z<0.4 from the first two years of data from the Baryon Oscillation Spectroscopic Survey (BOSS), to be released as part of SDSS Data Release 9 (DR9). We describe the sample selection, basic properties of the galaxies, and caveats for working with the data. We calculate the real- and redshift-space two-point correlation functions of these galaxies, fit these measurements using Halo Occupation Distribution (HOD) modeling within dark matter cosmological simulations, and estimate the errors using mock catalogs. These galaxies lie in massive halos, with a mean halo mass of 5.2x10^13 h^-1 M_sun, a large scale bias of ~2.0, and a satellite fraction of 12+/-2%. Thus, these galaxies occupy halos with average masses in between those of the higher redshift BOSS CMASS sample and the original SDSS I/II LRG sample.

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the scale of baryon acoustic oscillations (BAO) in the clustering of matter over a larger volume than the combined efforts of all previous spectroscopic surveys of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7. Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000 quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5. Early results from BOSS include the first detection of the large-scale three-dimensional clustering of the Lyman alpha forest and a strong detection from the Data Release 9 data set of the BAO in the clustering of massive galaxies at an effective redshift z = 0.57. We project that BOSS will yield measurements of the angular diameter distance D_A to an accuracy of 1.0% at redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the same redshifts. Forecasts for Lyman alpha forest constraints predict a measurement of an overall dilation factor that scales the highly degenerate D_A(z) and H^{-1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey is complete. Here, we provide an overview of the selection of spectroscopic targets, planning of observations, and analysis of data and data quality of BOSS.

(abridged) We describe the automated spectral classification, redshift determination, and parameter measurement pipeline in use for the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III (SDSS-III) as of Data Release 9, encompassing 831,000 moderate-resolution optical spectra. We give a review of the algorithms employed, and describe the changes to the pipeline that have been implemented for BOSS relative to previous SDSS-I/II versions, including new sets of stellar, galaxy, and quasar redshift templates. For the color-selected CMASS sample of massive galaxies at redshift 0.4 <~ z <~ 0.8 targeted by BOSS for the purposes of large-scale cosmological measurements, the pipeline achieves an automated classification success rate of 98.7% and confirms 95.4% of unique CMASS targets as galaxies (with the balance being mostly M stars). Based on visual inspections of a subset of BOSS galaxies, we find that ~0.2% of confidently reported CMASS sample classifications and redshifts are incorrect, and ~0.4% of all CMASS spectra are objects unclassified by the current algorithm which are potentially recoverable. The BOSS pipeline confirms that ~51.5% of the quasar targets have quasar spectra, with the balance mainly consisting of stars. Statistical (as opposed to systematic) redshift errors propagated from photon noise are typically a few tens of km/s for both galaxies and quasars, with a significant tail to a few hundreds of km/s for quasars. We test the accuracy of these statistical redshift error estimates using repeat observations, finding them underestimated by a factor of 1.19 to 1.34 for galaxies, and by a factor of 2 for quasars. We assess the impact of sky-subtraction quality, S/N, and other factors on galaxy redshift success. Finally, we document known issues, and describe directions of ongoing development.

In December 2010, NASA created a Science Definition Team (SDT) for WFIRST, the Wide Field Infra-Red Survey Telescope, recommended by the Astro 2010 Decadal Survey as the highest priority for a large space mission. The SDT was chartered to work with the WFIRST Project Office at GSFC and the Program Office at JPL to produce a Design Reference Mission (DRM) for WFIRST. Part of the original charge was to produce an interim design reference mission by mid-2011. That document was delivered to NASA and widely circulated within the astronomical community. In late 2011 the Astrophysics Division augmented its original charge, asking for two design reference missions. The first of these, DRM1, was to be a finalized version of the interim DRM, reducing overall mission costs where possible. The second of these, DRM2, was to identify and eliminate capabilities that overlapped with those of NASA's James Webb Space Telescope (henceforth JWST), ESA's Euclid mission, and the NSF's ground-based Large Synoptic Survey Telescope (henceforth LSST), and again to reduce overall mission cost, while staying faithful to NWNH. This report presents both DRM1 and DRM2.

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

We explore the benefits of using a passively evolving population of galaxies to measure the evolution of the rate of structure growth between z=0.25 and z=0.65 by combining data from the SDSS-I/II and SDSS-III surveys. The large-scale linear bias of a population of dynamically passive galaxies, which we select from both surveys, is easily modeled. Knowing the bias evolution breaks degeneracies inherent to other methodologies, and decreases the uncertainty in measurements of the rate of structure growth and the normalization of the galaxy power-spectrum by up to a factor of two. If we translate our measurements into a constraint on sigma_8(z=0) assuming a concordance cosmological model and General Relativity (GR), we find that using a bias model improves our uncertainty by a factor of nearly 1.5. Our results are consistent with a flat Lambda Cold Dark Matter model and with GR.

We obtain constraints on cosmological parameters from the spherically averaged redshift-space correlation function of the CMASS Data Release 9 (DR9) sample of the Baryonic Oscillation Spectroscopic Survey (BOSS). We combine this information with additional data from recent CMB, SN and BAO measurements. Our results show no significant evidence of deviations from the standard flat-Lambda CDM model, whose basic parameters can be specified by Omega_m = 0.285 +- 0.009, 100 Omega_b = 4.59 +- 0.09, n_s = 0.96 +- 0.009, H_0 = 69.4 +- 0.8 km/s/Mpc and sigma_8 = 0.80 +- 0.02. The CMB+CMASS combination sets tight constraints on the curvature of the Universe, with Omega_k = -0.0043 +- 0.0049, and the tensor-to-scalar amplitude ratio, for which we find r < 0.16 at the 95 per cent confidence level (CL). These data show a clear signature of a deviation from scale-invariance also in the presence of tensor modes, with n_s <1 at the 99.7 per cent CL. We derive constraints on the fraction of massive neutrinos of f_nu < 0.049 (95 per cent CL), implying a limit of sum m_nu < 0.51 eV. We find no signature of a deviation from a cosmological constant from the combination of all datasets, with a constraint of w_DE = -1.033 +- 0.073 when this parameter is assumed time-independent, and no evidence of a departure from this value when it is allowed to evolve as w_DE(a) = w_0 + w_a (1 - a). The achieved accuracy on our cosmological constraints is a clear demonstration of the constraining power of current cosmological observations.

We analyze the density field of galaxies observed by the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) included in the SDSS Data Release Nine (DR9). DR9 includes spectroscopic redshifts for over 400,000 galaxies spread over a footprint of 3,275 deg^2. We identify, characterize, and mitigate the impact of sources of systematic uncertainty on large-scale clustering measurements, both for angular moments of the redshift-space correlation function and the spherically averaged power spectrum, P(k), in order to ensure that robust cosmological constraints will be obtained from these data. A correlation between the projected density of stars and the higher redshift (0.43 < z < 0.7) galaxy sample (the `CMASS' sample) due to imaging systematics imparts a systematic error that is larger than the statistical error of the clustering measurements at scales s > 120h^-1Mpc or k < 0.01hMpc^-1. We find that these errors can be ameliorated by weighting galaxies based on their surface brightness and the local stellar density. We use mock galaxy catalogs that simulate the CMASS selection function to determine that randomly selecting galaxy redshifts in order to simulate the radial selection function of a random sample imparts the least systematic error on correlation function measurements and that this systematic error is negligible for the spherically averaged correlation function. The methods we recommend for the calculation of clustering measurements using the CMASS sample are adopted in companion papers that locate the position of the baryon acoustic oscillation feature (Anderson et al. 2012), constrain cosmological models using the full shape of the correlation function (Sanchez et al. 2012), and measure the rate of structure growth (Reid et al. 2012). (abridged)

We present a fast method of producing mock galaxy catalogues that can be used to compute covariance matrices of large-scale clustering measurements and test the methods of analysis. Our method populates a 2nd-order Lagrangian Perturbation Theory (2LPT) matter field, where we calibrate masses of dark matter halos by detailed comparisons with N-body simulations. We demonstrate the clustering of halos is recovered at ~10 per cent accuracy. We populate halos with mock galaxies using a Halo Occupation Distribution (HOD) prescription, which has been calibrated to reproduce the clustering measurements on scales between 30 and 80 Mpc/h. We compare the sample covariance matrix from our mocks with analytic estimates, and discuss differences. We have used this method to make catalogues corresponding to Data Release 9 of the Baryon Oscillation Spectroscopic Survey (BOSS),producing 600 mock catalogues of the "CMASS" galaxy sample. These mocks enabled detailed tests of methods and errors that formed an integral part of companion analyses of these galaxy data.

We analyze the anisotropic clustering of massive galaxies from the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) sample, which consists of 264,283 galaxies in the redshift range 0.43 < z < 0.7 spanning 3,275 square degrees. Both peculiar velocities and errors in the assumed redshift-distance relation ("Alcock-Paczynski effect") generate correlations between clustering amplitude and orientation with respect to the line-of-sight. Together with the sharp baryon acoustic oscillation (BAO) standard ruler, our measurements of the broadband shape of the monopole and quadrupole correlation functions simultaneously constrain the comoving angular diameter distance (2190 +/- 61 Mpc) to z=0.57, the Hubble expansion rate at z=0.57 (92.4 +/- 4.5 km/s/Mpc), and the growth rate of structure at that same redshift (d sigma8/d ln a = 0.43 +/- 0.069). Our analysis provides the best current direct determination of both DA and H in galaxy clustering data using this technique. If we further assume a LCDM expansion history, our growth constraint tightens to d sigma8/d ln a = 0.415 +/- 0.034. In combination with the cosmic microwave background, our measurements of DA, H, and growth all separately require dark energy at z > 0.57, and when combined imply \Omega_{\Lambda} = 0.74 +/- 0.016, independent of the Universe's evolution at z<0.57. In our companion paper (Samushia et al. prep), we explore further cosmological implications of these observations.

We present measurements of galaxy clustering from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). These use the Data Release 9 (DR9) CMASS sample, which contains 264,283 massive galaxies covering 3275 square degrees with an effective redshift z=0.57 and redshift range 0.43 < z < 0.7. Assuming a concordance Lambda-CDM cosmological model, this sample covers an effective volume of 2.2 Gpc^3, and represents the largest sample of the Universe ever surveyed at this density, n = 3 x 10^-4 h^-3 Mpc^3. We measure the angle-averaged galaxy correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of 5\sigma in both the correlation function and power spectrum. Combining with the SDSS-II Luminous Red Galaxy Sample, the detection significance increases to 6.7\sigma. Fitting for the position of the acoustic features measures the distance to z=0.57 relative to the sound horizon DV /rs = 13.67 +/- 0.22 at z=0.57. Assuming a fiducial sound horizon of 153.19 Mpc, which matches cosmic microwave background constraints, this corresponds to a distance DV(z=0.57) = 2094 +/- 34 Mpc. At 1.7 per cent, this is the most precise distance constraint ever obtained from a galaxy survey. We place this result alongside previous BAO measurements in a cosmological distance ladder and find excellent agreement with the current supernova measurements. We use these distance measurements to constrain various cosmological models, finding continuing support for a flat Universe with a cosmological constant.

We use the 2% distance measurement from our reconstructed baryon acoustic oscillations (BAOs) signature using the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) Luminous Red Galaxies (LRGs) from Padmanabhan et al. (2012) and Xu et al. (2012) combined with cosmic microwave background (CMB) data from Wilkinson Microwave Anisotropy Probe (WMAP7) to measure parameters for various cosmological models. We find a 1.7% measurement of H_0 = 69.8 +/- 1.2 km/s/Mpc and a 5.0% measurement of Omega_m = 0.280 +/- 0.014 for a flat Universe with a cosmological constant. These measurements of H_0 and Omega_m are robust against a range of underlying models for the expansion history. We measure the dark energy equation of state parameter w = -0.97 +/- 0.17, which is consistent with a cosmological constant. If curvature is allowed to vary, we find that the Universe is consistent with a flat geometry (Omega_K = -0.004 +/- 0.005). We also use a combination of the 6 Degree Field Galaxy Survey BAO data, WiggleZ Dark Energy Survey data, Type Ia supernovae (SN) data, and a local measurement of the Hubble constant to explore cosmological models with more parameters. Finally, we explore the effect of varying the energy density of relativistic particles on the measurement of H_0.

We present results from fitting the baryon acoustic oscillation (BAO) signal in the correlation function obtained from the first application of reconstruction to a galaxy redshift survey, namely, the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) luminous red galaxy (LRG) catalogue. We also introduce more careful approaches for deriving a suitable covariance matrix and fitting model for galaxy correlation functions. These all aid in obtaining a more accurate measurement of the acoustic scale and its error. We validate our reconstruction, covariance matrix and fitting techniques on 160 mock catalogues derived from the LasDamas simulations in real and redshift space. We then apply these techniques to the DR7 LRG sample and find that the error on the acoustic scale decreases from ~3.5% before reconstruction to ~1.9% after reconstruction. This factor of 1.8 reduction in the error is equivalent to the effect of tripling the survey volume. We also see an increase in our BAO detection confidence from ~3-sigma to ~4-sigma after reconstruction with our confidence level in measuring the correct acoustic scale increasing from ~3-sigma to ~5-sigma. Using the mean of the acoustic scale probability distributions produced from our fits, we find D_v/r_s = 8.89 +/- 0.31 before reconstruction and 8.88 +/- 0.17 after reconstruction.

We apply the reconstruction technique to the clustering of galaxies from the SDSS DR7 LRG sample, sharpening the baryon acoustic oscillation (BAO) feature and achieving a 1.9% measurement of the distance to z=0.35. This is the first application of reconstruction of the BAO feature in a galaxy redshift survey. We update the reconstruction algorithm of Eisenstein et al, 2007 to account for the effects of survey geometry as well as redshift-space distortions and validate it on 160 LasDamas simulations. We demonstrate that reconstruction sharpens the BAO feature in the angle averaged galaxy correlation function, reducing the nonlinear smoothing scale \Sigma_nl from 8.1 Mpc/h to 4.4 Mpc/h. Reconstruction also significantly reduces the effects of redshift-space distortions at the BAO scale, isotropizing the correlation function. This sharpened BAO feature yields an unbiased distance estimate (< 0.2%) and reduces the scatter from 3.3% to 2.1%. We demonstrate the robustness of these results to the various reconstruction parameters, including the smoothing scale, the galaxy bias and the linear growth rate. Applying this reconstruction algorithm to the SDSS LRG DR7 sample improves the significance of the BAO feature in these data from 3.3 sigma for the unreconstructed correlation function, to 4.2 sigma after reconstruction. We estimate a relative distance scale D_V/r_s to z=0.35 of 8.88+/-0.17, where r_s is the sound horizon and D_V = (D_A^2/H)^{1/3} is a combination of the angular diameter distance D_A and Hubble parameter H. Assuming a sound horizon of 154.25 Mpc, this translates into a distance measurement D_V (z=0.35) = 1.356+/-0.025 Gpc. We find that reconstruction reduces the distance error in the DR7 sample from 3.5% to 1.9%, equivalent to a survey with three times the volume of SDSS.

The Sloan Digital Sky Survey (SDSS) surveyed 14,555 square degrees, and delivered over a trillion pixels of imaging data. We present a study of galaxy clustering using 900,000 luminous galaxies with photometric redshifts, spanning between $z=0.45$ and $z=0.65$, constructed from the SDSS using methods described in Ross et al. (2011). This data-set spans 11,000 square degrees and probes a volume of $3h^{-3} \rm{Gpc}^3$, making it the largest volume ever used for galaxy clustering measurements. We present a novel treatment of the observational systematics and its applications to the clustering signals from the data set. In this paper, we measure the angular clustering using an optimal quadratic estimator at 4 redshift slices with an accuracy of ~15% with bin size of delta_l = 10 on scales of the Baryon Acoustic Oscillations (BAO) (at l~40-400). We derive cosmological constraints using the full-shape of the power-spectra. For a flat Lambda CDM model, when combined with Cosmic Microwave Background Wilkinson Microwave Anisotropy Probe 7 (WMAP7) and H_0 constraints from 600 Cepheids observed by HST, we find \Omega_\Lambda = 0.73 +/- 0.019 and H_0 to be 70.5 +/- 1.6 km/s/Mpc. For an open Lambda CDM model, when combined with WMAP7 + HST, we find $\Omega_K = 0.0035 +/- 0.0054, improved over WMAP7+HST alone by 40%. For a wCDM model, when combined with WMAP7+HST+SN, we find w = -1.071 +/- 0.078, and H_0 to be 71.3 +/- 1.7 km/s/Mpc, which is competitive with the latest large scale structure constraints from large spectroscopic surveys such as SDSS Data Release 7 (DR7) (Reid et al. 2010, Percival et al. 2010, Montesano et al. 2011) and WiggleZ (Blake et al. 2011). The SDSS-III Data Release 8 (SDSS-III DR8) Angular Clustering Data allows a wide range of investigations into the cosmological model, cosmic expansion (via BAO), Gaussianity of initial conditions and neutrino masses. (abridged)

We measure the acoustic scale from the angular power spectra of the Sloan Digital Sky Survey III (SDSS-III) Data Release 8 imaging catalog that includes 872,921 galaxies over ~ 10,000 deg^2 between 0.45<z<0.65. The extensive spectroscopic training set of the Baryon Oscillation Spectroscopic Survey (BOSS) luminous galaxies allows precise estimates of the true redshift distributions of galaxies in our imaging catalog. Utilizing the redshift distribution information, we build templates and fit to the power spectra of the data, which are measured in our companion paper, Ho et al. 2011, to derive the location of Baryon acoustic oscillations (BAO) while marginalizing over many free parameters to exclude nearly all of the non-BAO signal. We derive the ratio of the angular diameter distance to the sound horizon scale D_A/r_s= 9.212 + 0.416 -0.404 at z=0.54, and therefore, D_A= 1411+- 65 Mpc at z=0.54; the result is fairly independent of assumptions on the underlying cosmology. Our measurement of angular diameter distance D_A is 1.4 \sigma higher than what is expected for the concordance LCDM (Komatsu et al. 2011), in accordance to the trend of other spectroscopic BAO measurements for z >~ 0.35. We report constraints on cosmological parameters from our measurement in combination with the WMAP7 data and the previous spectroscopic BAO measurements of SDSS (Percival et al. 2010) and WiggleZ (Blake et al. 2011). We refer to our companion papers (Ho et al. 2011; de Putter et al. 2011) for investigations on information of the full power spectrum.

The discovery of Ultra-Faint Dwarf (UFD) galaxies in the halo of the Milky Way extends the faint end of the galaxy luminosity function to a few hundred solar luminosities. This extremely low luminosity regime poses a significant challenge for the photometric characterization of these systems. We present a suite of simulations aimed at understanding how different observational choices related to the properties of a low luminosity system impact our ability to determine its true structural parameters such as half-light radius and central surface brightness. We focus on estimating half-light radii (on which mass estimates depend linearly) and find that these numbers can have up to 100% uncertainties when relatively shallow photometric surveys, such as SDSS, are used. Our simulations suggest that to recover structural parameters within 10% or better of their true values: (a) the ratio of the field-of-view to the half-light radius of the satellite must be greater than three, (b) the total number of stars, including background objects should be larger than 1000, and (c) the central to background stellar density ratio must be higher than 20. If one or more of these criteria are not met, the accuracy of the resulting structural parameters can be significantly compromised. In the context of future surveys such as LSST, the latter condition will be closely tied to our ability to remove unresolved background galaxies. Assessing the reliability of measured structural parameters will become increasingly critical as the next generation of deep wide-field surveys detects UFDs beyond the reach of current spectroscopic limits.

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

We investigate the effects of potential sources of systematic error on the angular and photometric redshift, z_phot, distributions of a sample of redshift 0.4 < z < 0.7 massive galaxies whose selection matches that of the Baryon Oscillation Spectroscopic Survey (BOSS) constant mass sample. Utilizing over 112,778 BOSS spectra as a training sample, we produce a photometric redshift catalog for the galaxies in the SDSS DR8 imaging area that, after masking, covers nearly one quarter of the sky (9,913 square degrees). We investigate fluctuations in the number density of objects in this sample as a function of Galactic extinction, seeing, stellar density, sky background, airmass, photometric offset, and North/South Galactic hemisphere. We find that the presence of stars of comparable magnitudes to our galaxies (which are not traditionally masked) effectively remove area. Failing to correct for such stars can produce systematic errors on the measured angular auto-correlation function, w, that are larger than its statistical uncertainty. We describe how one can effectively mask for the presence of the stars, without removing any galaxies from the sample, and minimize the systematic error. Additionally, we apply two separate methods that can be used to correct the systematic errors imparted by any parameter that can be turned into a map on the sky. We find that failing to properly account for varying sky background introduces a systematic error on w. We measure w, in four z_phot slices of width 0.05 between 0.45 < z_phot < 0.65 and find that the measurements, after correcting for the systematic effects of stars and sky background, are generally consistent with a generic LambdaCDM model, at scales up to 60 degrees. At scales greater than 3 degrees and z_phot > 0.5, the magnitude of the corrections we apply are greater than the statistical uncertainty in w.

BigBOSS is a Stage IV ground-based dark energy experiment to study baryon acoustic oscillations (BAO) and the growth of structure with a wide-area galaxy and quasar redshift survey over 14,000 square degrees. It has been conditionally accepted by NOAO in response to a call for major new instrumentation and a high-impact science program for the 4-m Mayall telescope at Kitt Peak. The BigBOSS instrument is a robotically-actuated, fiber-fed spectrograph capable of taking 5000 simultaneous spectra over a wavelength range from 340 nm to 1060 nm, with a resolution R = 3000-4800. Using data from imaging surveys that are already underway, spectroscopic targets are selected that trace the underlying dark matter distribution. In particular, targets include luminous red galaxies (LRGs) up to z = 1.0, extending the BOSS LRG survey in both redshift and survey area. To probe the universe out to even higher redshift, BigBOSS will target bright [OII] emission line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are obtained to measure the BAO feature, trace the matter power spectrum at smaller scales, and detect redshift space distortions. BigBOSS will provide additional constraints on early dark energy and on the curvature of the universe by measuring the Ly-alpha forest in the spectra of over 600,000 2.2 < z < 3.5 quasars. BigBOSS galaxy BAO measurements combined with an analysis of the broadband power, including the Ly-alpha forest in BigBOSS quasar spectra, achieves a FOM of 395 with Planck plus Stage III priors. This FOM is based on conservative assumptions for the analysis of broad band power (kmax = 0.15), and could grow to over 600 if current work allows us to push the analysis to higher wave numbers (kmax = 0.3). BigBOSS will also place constraints on theories of modified gravity and inflation, and will measure the sum of neutrino masses to 0.024 eV accuracy.

Strong foreground absorption features from singly-ionized Magnesium (Mg II) are commonly observed in the spectra of quasars and are presumed to probe a wide range of galactic environments. To date, measurements of the average dark matter halo masses of intervening Mg II absorbers by way of large-scale cross-correlations with luminous galaxies have been limited to z<0.7. In this work we cross-correlate 21 strong (W{\lambda}2796>0.6 {\deg}A) Mg II absorption systems detected in quasar spectra from the Sloan Digital Sky Survey Data Release 7 with ~32,000 spectroscopically confirmed galaxies at 0.7<z<1.45 from the DEEP2 galaxy redshift survey. We measure dark matter (DM) halo biases of b_G=1.44\pm0.02 and b_A=1.49\pm0.45 for the DEEP2 galaxies and Mg II absorbers, respectively, indicating that their clustering amplitudes are roughly consistent. Haloes with the bias we measure for the Mg II absorbers have a corresponding mass of 1.8(+4.2/-1.6) \times 10^12h-1M_sun, although the actual mean absorber halo mass will depend on the precise distribution of absorbers within DM haloes. This mass estimate is consistent with observations at z=0.6, suggesting that the halo masses of typical Mg II absorbers do not significantly evolve from z~1. We additionally measure the average W{\lambda}2796>0.6 \AA gas covering fraction to be f =0.5 within 60 h-1kpc around the DEEP2 galaxies, and we find an absence of coincident strong Mg II absorption beyond a projected separation of ~40 h-1kpc. Although the star-forming z>1 DEEP2 galaxies are known to exhibit ubiquitous blueshifted Mg II absorption, we find no direct evidence in our small sample linking W{\lambda}2796>0.6 \AA absorbers to galaxies with ongoing star formation.

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

We measure shifts of the acoustic scale due to nonlinear growth and redshift distortions to a high precision using a very large volume of high-force-resolution simulations. We compare results from various sets of simulations that differ in their force, volume, and mass resolution. We find a consistency within 1.5-sigma for shift values from different simulations and derive shift alpha(z) -1 = (0.300\pm 0.015)% [D(z)/D(0)]^{2} using our fiducial set. We find a strong correlation with a non-unity slope between shifts in real space and in redshift space and a weak correlation between the initial redshift and low redshift. Density-field reconstruction not only removes the mean shifts and reduces errors on the mean, but also tightens the correlations: after reconstruction, we recover a slope of near unity for the correlation between the real and redshift space and restore a strong correlation between the low and the initial redshifts. We derive propagators and mode-coupling terms from our N-body simulations and compared with Zeldovich approximation and the shifts measured from the chi^2 fitting, respectively. We interpret the propagator and the mode-coupling term of a nonlinear density field in the context of an average and a dispersion of its complex Fourier coefficients relative to those of the linear density field; from these two terms, we derive a signal-to-noise ratio of the acoustic peak measurement. We attempt to improve our reconstruction method by implementing 2LPT and iterative operations: we obtain little improvement. The Fisher matrix estimates of uncertainty in the acoustic scale is tested using 5000 (Gpc/h)^3 of cosmological PM simulations from Takahashi et al. (2009). (abridged)

We calculate the real- and redshift-space clustering of massive galaxies at z~0.5 using the first semester of data by the Baryon Oscillation Spectroscopic Survey (BOSS). We study the correlation functions of a sample of 44,000 massive galaxies in the redshift range 0.4<z<0.7. We present a halo-occupation distribution modeling of the clustering results and discuss the implications for the manner in which massive galaxies at z~0.5 occupy dark matter halos. The majority of our galaxies are central galaxies living in halos of mass 10^{13}Msun/h, but 10% are satellites living in halos 10 times more massive. These results are broadly in agreement with earlier investigations of massive galaxies at z~0.5. The inferred large-scale bias (b~2) and relatively high number density (nbar=3e-4 h^3 Mpc^{-3}) imply that BOSS galaxies are excellent tracers of large-scale structure, suggesting BOSS will enable a wide range of investigations on the distance scale, the growth of large-scale structure, massive galaxy evolution and other topics.

We introduce a new statistic omega_l for measuring and analyzing large-scale structure and particularly the baryon acoustic oscillations. omega_l is a band-filtered, configuration space statistic that is easily implemented and has advantages over the traditional power spectrum and correlation function estimators. Unlike these estimators, omega_l can localize most of the acoustic information into a single dip at the acoustic scale while also avoiding sensitivity to the poorly constrained large scale power (i.e., the integral constraint) through the use of a localized and compensated filter. It is also sensitive to anisotropic clustering through pair counting and does not require any binning. We measure the shift in the acoustic peak due to nonlinear effects using the monopole omega_0 derived from subsampled dark matter catalogues as well as from mock galaxy catalogues created via halo occupation distribution (HOD) modeling. All of these are drawn from 44 realizations of 1024^3 particle dark matter simulations in a 1h^{-1}Gpc box at z=1. We compare these shifts with those obtained from the power spectrum and conclude that the results agree. This indicates that any distance measurements obtained from omega_0 and P(k) will be consistent with each other. We also show that it is possible to extract the same amount of acoustic information using either omega_0 or P(k) from equal volume surveys.

We study the profiles of 75 086 elliptical galaxies from the Sloan Digital Sky Survey (SDSS) at both large (50-500 kpc/h) and small (~3 kpc/h) scales. Weak lensing observations in the outskirts of the halo are combined with measurements of the stellar velocity dispersion in the interior regions of the galaxy for stacked galaxy samples. The weak lensing measurements are well characterized by a Navarro, Frenk and White (NFW) profile. The dynamical mass measurements exceed the extrapolated NFW profile even after the estimated stellar masses are subtracted, providing evidence for the modification of the dark matter profile by the baryons. This excess mass is quantitatively consistent with the predictions of the adiabatic contraction (AC) hypothesis. Our finding suggests that the effects of AC during galaxy formation are stable to subsequent bombardment from major and minor mergers. We explore several theoretical and observational systematics and conclude that they cannot account for the inferred mass excess. The most significant source of systematic error is in the IMF, which would have to increase the stellar mass estimates by a factor of two relative to masses from the Kroupa IMF to fully explain the mass excess without AC. Such an increase would create tension with results from SAURON (Cappellari et al. 2006). We demonstrate a connection between the level of contraction of the dark matter halo profile and scatter in the size-luminosity relation, which is a projection of the fundamental plane. Whether or not AC is the mechanism supplying the excess mass, models of galaxy formation and evolution must reconcile the observed halo masses from weak lensing with the comparatively large dynamical masses at the half light radii of the galaxies.

We characterize the nonlinear evolution of the baryon acoustic feature as traced by the dark matter and halos, using a combination of perturbation theory and N-body simulations. We confirm that the acoustic peak traced by the dark matter is both broadened and shifted as structure forms, and that this shift is well described by second-order perturbation theory. These shifts persist for dark matter halos, and are a simple function of halo bias, with the shift (mostly) increasing with increasing bias. Extending our perturbation theory results to halos with simple two parameter bias models (both in Lagrangian and Eulerian space) quantitatively explains the observed shifts. In particular, we demonstrate that there are additional terms that contribute to the shift that are absent for the matter. At z=0 for currently favored cosmologies, the matter shows shifts of ~0.5%, b=1 halos shift the acoustic scale by ~0.2%, while b=2 halos shift it by ~0.5%; these shifts decrease by the square of the growth factor at higher redshifts. These results are easily generalized to galaxies within the halo model, where we show that simple galaxy models show marginally larger shifts than the correspondingly biased halos, due to the contribution of satellites in high mass halos. While our focus here is on real space, our results make specific predictions for redshift space. For currently favored cosmological models, we find that the shifts for halos at z=0 increase by ~0.3%; at high z, they increase by ~0.5% D^2. Our results demonstrate that these theoretical systematics are smaller than the statistical precision of upcoming surveys, even if one ignored the corrections discussed here. Simple modeling, along the lines discussed here, has the potential to reduce these systematics to below the levels of cosmic variance limited surveys.

The baryon acoustic oscillation (BAO) method for constraining the expansion history is adversely affected by non-linear structure formation, which washes out the correlation function peak created at decoupling. To increase the constraining power of low z BAO experiments, it has been proposed that one use the observed distribution of galaxies to "reconstruct'' the acoustic peak. Recently Padmanabhan, White and Cohn provided an analytic formalism for understanding how reconstruction works within the context of Lagrangian perturbation theory. We extend that formalism to include the case of biased tracers of the mass and, because the quantitative validity of LPT is questionable, we investigate reconstruction in N-body simulations. We find that LPT does a good job of explaining the trends seen in simulations for both the mass and for biased tracers and comment upon the implications this has for reconstruction.

We compute in detail the rate at which energy injected by dark matter annihilation heats and ionizes the photon-baryon plasma at z ~ 1000, and provide accurate fitting functions over the relevant redshift range for a broad array of annihilation channels and DM masses. The resulting perturbations to the ionization history can be constrained by measurements of the CMB temperature and polarization angular power spectra. We show that models which fit recently measured excesses in 10-1000 GeV electron and positron cosmic rays are already close to the 95% confidence limits from WMAP. The recently launched Planck satellite will be capable of ruling out a wide range of DM explanations for these excesses. In models of dark matter with Sommerfeld-enhanced annihilation, where sigma v rises with decreasing WIMP velocity until some saturation point, the WMAP5 constraints imply that the enhancement must be close to saturation in the neighborhood of the Earth.

The BigBOSS experiment is a proposed DOE-NSF Stage IV ground-based dark energy experiment to study baryon acoustic oscillations (BAO) and the growth of structure with an all-sky galaxy redshift survey. The project is designed to unlock the mystery of dark energy using existing ground-based facilities operated by NOAO. A new 4000-fiber R=5000 spectrograph covering a 3-degree diameter field will measure BAO and redshift space distortions in the distribution of galaxies and hydrogen gas spanning redshifts from 0.2<z<3.5. The Dark Energy Task Force figure of merit (DETF FoM) for this experiment is expected to be equal to that of a JDEM mission for BAO with the lower risk and cost typical of a ground-based experiment. This project will enable an unprecedented multi-object spectroscopic capability for the U.S. community through an existing NOAO facility. The U.S. community would have access directly to this instrument/telescope combination, as well as access to the legacy archives that will be created by the BAO key project.

Recently a number of analytic prescriptions for computing the non-linear matter power spectrum have appeared in the literature. These typically involve resummation or closure prescriptions which do not have a rigorous error control, thus they must be compared with numerical simulations to assess their range of validity. We present a direct side-by-side comparison of several of these analytic approaches, using a suite of high-resolution N-body simulations as a reference, and discuss some general trends. All of the analytic results correctly predict the behavior of the power spectrum at the onset of non-linearity, and improve upon a pure linear theory description at very large scales. All of these theories fail at sufficiently small scales. At low redshift the dynamic range in scale where perturbation theory is both relevant and reliable can be quite small. We also compute for the first time the 2-loop contribution to standard perturbation theory for CDM models, finding improved agreement with simulations at large redshift. At low redshifts however the 2-loop term is larger than the 1-loop term on quasi-linear scales, indicating a breakdown of the perturbation expansion. Finally, we comment on possible implications of our results for future studies.

Recently Eisenstein and collaborators introduced a method to `reconstruct' the linear power spectrum from a non-linearly evolved galaxy distribution in order to improve precision in measurements of baryon acoustic oscillations. We reformulate this method within the Lagrangian picture of structure formation, to better understand what such a method does, and what the resulting power spectra are. We show that reconstruction does not reproduce the linear density field, at second order. We however show that it does reduce the damping of the oscillations due to non-linear structure formation, explaining the improvements seen in simulations. Our results suggest that the reconstructed power spectrum is potentially better modeled as the sum of three different power spectra, each dominating over different wavelength ranges and with different non-linear damping terms. Finally, we also show that reconstruction reduces the mode-coupling term in the power spectrum, explaining why mis-calibrations of the acoustic scale are reduced when one considers the reconstructed power spectrum.

We show that a suitably defined marked correlation function can be used to break degeneracies in halo-occupation distribution modeling. The statistic can be computed on both 3D and 2D data sets, and should be applicable to all upcoming galaxy surveys. A proof of principle, using mock catalogs created from N-body simulations, is given.

Motivated by the eXciting Dark Matter (XDM) model of Finkbeiner & Weiner, hypothesized to explain the 511 keV signal in the center of the Milky Way, we consider the CMB and 21-cm signatures of models of dark matter with collisional long-lived excited states. We compute the relic excitation fraction from the early universe for a variety of assumptions about the collisional de-excitation cross-section and thermal decoupling. The relic excitation fraction can be as high as 1% for natural regions of parameter space, but could be orders of magnitude smaller. Since the lifetime of the excited state is naturally greater than 10^13s, we discuss the signatures of such relic excitation on cosmic microwave background (CMB) and high-z 21-cm observations. Such models have potentially richer astrophysical signals than the traditional WIMP annihilations and decays, and may have observable consequences for future generations of experiments.

We cross-correlate large scale structure (LSS) observations from a number of surveys with CMB anisotropies from WMAP to investigate the Integrated Sachs-Wolfe (ISW) effect as a function of redshift, covering z~0.1-2.5. Our main goal is to go beyond reporting detections towards developing a reliable likelihood analysis that allows one to determine cosmological constraints from ISW observations. With this in mind we spend a considerable amount of effort in determining the redshift-dependent bias and redshift distribution b(z)*dN/dz of these samples by matching with spectroscopic observations where available, and analyzing auto-power spectra and cross-power spectra between the samples. The data sets we use are 2-Micron All Sky Survey (2MASS) samples, Sloan Digital Sky Survey (SDSS) photometric Luminous Red Galaxies, SDSS photometric quasars and NRAO VLA Sky Survey (NVSS) radio sources. We make a joint analysis of all samples constructing a full covariance matrix, which we subsequently use for cosmological parameter fitting. We report a 3.7 sigma detection of ISW combining all the datasets. We combine the ISW likelihood function with weak lensing of CMB and CMB power spectrum to constrain the equation of state of dark energy and the curvature of the Universe. While ISW does not significantly improve the constraints in the simplest 6-parameter flat Lambda CDM model, it improves constraints on 7-parameter models with curvature by a factor of 3.2 (relative to WMAP alone) to Omega_K=-0.004^{+0.014}_{-0.020}, and with dark energy equation of state by 15% to w=-1.01^{+0.30}_{-0.40}. (Abridged.)

Recent work has shown that the local non-Gaussianity parameter f_NL induces a scale-dependent bias, whose amplitude is growing with scale. Here we first rederive this result within the context of peak-background split formalism and show that it only depends on the assumption of universality of mass function, assuming halo bias only depends on mass. We then use extended Press-Schechter formalism to argue that this assumption may be violated and the scale dependent bias will depend on other properties, such as merging history of halos. In particular, in the limit of recent mergers we find the effect is suppressed. Next we use these predictions in conjunction with a compendium of large scale data to put a limit on the value of f_NL. When combining all data assuming that halo occupation depends only on halo mass, we get a limit of -29 ~ (-65)< f_NL < +70 ~(+93) at 95% (99.7%) confidence. While we use a wide range of datasets, our combined result is dominated by the signal from the SDSS photometric quasar sample. If the latter are modeled as recent mergers then the limits weaken to -31 ~(-96) < f_NL < +70 ~ (+96) . These limits are comparable to the strongest current limits from the WMAP 5 year analysis, with no evidence of a positive signal in f_NL. While the method needs to be thoroughly tested against large scale structure simulations with realistic quasar and galaxy formation models, our results indicate that this is a competitive method relative to CMB and should be further pursued both observationally and theoretically.

We investigate the correlation of gravitational lensing of the cosmic microwave background (CMB) with several tracers of large-scale structure, including luminous red galaxies (LRGs), quasars, and radio sources. The lensing field is reconstructed based on the CMB maps from the Wilkinson Microwave Anisotropy Probe (WMAP) satellite; the LRGs and quasars are observed by the Sloan Digital Sky Survey (SDSS); and the radio sources are observed in the NRAO VLA Sky Survey (NVSS). Combining all three large-scale structure samples, we find evidence for a positive cross-correlation at the $2.5\sigma$ level ($1.8\sigma$ for the SDSS samples and $2.1\sigma$ for NVSS); the cross-correlation amplitude is $1.06\pm 0.42$ times that expected for the WMAP cosmological parameters. Our analysis extends other recent analyses in that we carefully determine bias weighted redshift distribution of the sources, which is needed for a meaningful cosmological interpretation of the detected signal. We investigate contamination of the signal by Galactic emission, extragalactic radio and infrared sources, thermal and kinetic Sunyaev-Zel'dovich effects, and the Rees-Sciama effect, and find all of them to be negligible.

Abridged: We estimate the distances to ~48 million stars detected by the Sloan Digital Sky Survey and map their 3D number density distribution in 100 < D < 20 kpc range over 6,500 deg^2 of sky. The data show strong evidence for a Galaxy consisting of an oblate halo, a disk component, and a number of localized overdensities with exponential disk parameters (bias-corrected for an assumed 35% binary fraction) H_1 = 300 pc, L_1 = 2600 pc, H_2 = 900 pc, L_2 = 3600 pc, and local density normalization of 12%. We find the halo to be oblate, with best-fit axis ratio c/a = 0.64, r^{-2.8} profile, and the local halo-to-thin disk normalization of 0.5%. We estimate the errors of derived model parameters to be no larger than ~20% (disk scales) and ~10% (thick disk normalization). While generally consistent with the above model, the density distribution shows a number of statistically significant localized deviations. We detect two overdensities in the thick disk region at (R, Z) ~ (6.5, 1.5)kpc and (R, Z) ~ (9.5, 0.8) kpc, and a remarkable density enhancement in the halo covering >1000deg^2 of sky towards the constellation of Virgo, at distances of ~6-20 kpc. Compared to a region symmetric with respect to the l=0 line, the Virgo overdensity is responsible for a factor of 2 number density excess and may be a nearby tidal stream or a low-surface brightness dwarf galaxy merging with the Milky Way. After removal of the resolved overdensities, the remaining data are consistent with a smooth density distribution; we detect no evidence of further unresolved clumpy substructure at scales ranging from ~50pc in the disk, to ~1 - 2 kpc in the halo.

We present an analysis of anisotropic baryon acoustic oscillations and elucidate how a mis-estimation of the cosmology, which leads to incorrect values of the angular diameter distance, d_A, and Hubble parameter, H, manifest themselves in changes to the monopole and quadrupole power spectrum of biased tracers of the density field. Previous work has focused on the monopole power spectrum, and shown that the isotropic "dilation" combination d_A^2/H is robustly constrained by an overall shift in the scale of the baryon feature. We extend this by demonstrating that the quadrupole power spectrum is sensitive to an anisotropic "warping" mode d_A H, allowing one to break the degeneracy between d_A and H. We describe a method for measuring this warping, explicitly marginalizing over the form of redshift space distortions. We verify this method on N-body simulations and estimate that d_A H can be measured with a fractional accuracy of ~ 3/sqrt(V) % where the survey volume is estimated in (Gpc/h)^3.

We measure the clustering of a sample of photometrically selected luminous red galaxies around a low redshift (0.2<z<0.6) sample of quasars selected from the Sloan Digital Sky Survey Data Release 5. We make use of a new statistical estimator to obtain precise measurements of the LRG auto-correlations and constrain halo occupation distributions for them. These are used to generate mock catalogs which aid in interpreting our quasar-LRG cross correlation measurements. The cross correlation is well described by a power law with slope 1.8\pm0.1 and r_0=6\pm0.5 h^{-1} Mpc, consistent with observed galaxy correlation functions. We find no evidence for `excess' clustering on 0.1 Mpc scales and demonstrate that this is consistent with the results of Serber et al (2006) and Strand et al (2007), when one accounts for several subtleties in the interpretation of their measurements. Combining the quasar-LRG cross correlation with the LRG auto-correlations, we determine a large-scale quasar bias b_QSO = 1.09\pm0.15 at a median redshift of 0.43, with no observed redshift or luminosity evolution. This corresponds to a mean halo mass <M>~ 10^{12} h^{-1} M_sun, Eddington ratios from 0.01 to 1 and lifetimes less than 10^{7} yr. Using simple models of halo occupation, these correspond to a number density of quasar hosts greater than 10^{-3} h^{3} Mpc^{-3} and stellar masses less than 10^{11} h^{-1} M_sun. The small-scale clustering signal can be interpreted with the aid of our mock LRG catalogs, and depends on the manner in which quasars inhabit halos. We find that our small scale measurements are inconsistent with quasar positions being randomly subsampled from halo centers above a mass threshold, requiring a satellite fraction > 25 per cent.

We present an algorithm to photometrically calibrate wide field optical imaging surveys, that simultaneously solves for the calibration parameters and relative stellar fluxes using overlapping observations. The algorithm decouples the problem of "relative" calibrations, from that of "absolute" calibrations; the absolute calibration is reduced to determining a few numbers for the entire survey. We pay special attention to the spatial structure of the calibration errors, allowing one to isolate particular error modes in downstream analyses. Applying this to the Sloan Digital Sky Survey imaging data, we achieve ~1% relative calibration errors across 8500 sq.deg. in griz; the errors are ~2% for the u band. These errors are dominated by unmodelled atmospheric variations at Apache Point Observatory. These calibrations, dubbed "ubercalibration", are now public with SDSS Data Release 6, and will be a part of subsequent SDSS data releases.

The Sloan Digital Sky Survey (SDSS) and Two Micron All Sky Survey (2MASS) are rich resources for studying stellar astrophysics and the structure and formation history of the Galaxy. As new surveys and instruments adopt similar filter sets, it is increasingly important to understand the properties of the ugrizJHKs stellar locus, both to inform studies of `normal' main sequence stars as well as for robust searches for point sources with unusual colors. Using a sample of ~600,000 point sources detected by SDSS and 2MASS, we tabulate the position and width of the ugrizJHKs stellar locus as a function of g-i color, and provide accurate polynomial fits. We map the Morgan-Keenan spectral type sequence to the median stellar locus by using synthetic photometry of spectral standards and by analyzing 3000 SDSS stellar spectra with a custom spectral typing pipeline. We develop an algorithm to calculate a point source's minimum separation from the stellar locus in a seven-dimensional color space, and use it to robustly identify objects with unusual colors, as well as spurious SDSS/2MASS matches. Analysis of a final catalog of 2117 color outliers identifies 370 white-dwarf/M dwarf (WDMD) pairs, 93 QSOs, and 90 M giant/carbon star candidates, and demonstrates that WDMD pairs and QSOs can be distinguished on the basis of their J-Ks and r-z colors. We also identify a group of objects with correlated offsets in the u-g vs. g-r and g-r vs. r-i color-color spaces, but subsequent follow-up is required to reveal the nature of these objects. Future applications of this algorithm to a matched SDSS-UKIDSS catalog may well identify additional classes of objects with unusual colors by probing new areas of color-magnitude space.

We present measurements of the color and luminosity dependence of galaxy clustering at z~1 in the DEEP2 Galaxy Redshift Survey. Using volume-limited subsamples in bins of both color and luminosity, we find that: 1) The clustering dependence is much stronger with color than with luminosity and is as strong with color at z~1 as is found locally. We find no dependence of the clustering amplitude on color for galaxies on the red sequence, but a significant dependence on color for galaxies within the blue cloud. 2) For galaxies in the range L/L*~0.7-2, a stronger large-scale luminosity dependence is seen for all galaxies than for red and blue galaxies separately. The small-scale clustering amplitude depends significantly on luminosity for blue galaxies, with brighter samples having a stronger rise on scales r_p<0.5 Mpc/h. 3) Redder galaxies exhibit stronger small-scale redshift-space distortions ("fingers of god"), and both red and blue populations show large-scale distortions in xi(r_p,pi) due to coherent infall. 4) While the clustering length, r_0, increases smoothly with galaxy color (in narrow bins), its power-law exponent, gamma, exhibits a sharp jump from the blue cloud to the red sequence. The intermediate color `green' galaxy population likely includes transitional galaxies moving from the blue cloud to the red sequence; on large scales green galaxies are as clustered as red galaxies but show infall kinematics and a small-scale correlation slope akin to the blue galaxy population. 5) We compare our results to a semi-analytic galaxy formation model applied to the Millenium Run simulation. Differences between the data and the model suggest that in the model star formation is shut down too efficiently in satellite galaxies.

We quantify the variability of faint unresolved optical sources using a catalog based on multiple SDSS imaging observations. The catalog covers SDSS Stripe 82, and contains 58 million photometric observations in the SDSS ugriz system for 1.4 million unresolved sources. In each photometric bandpass we compute various low-order lightcurve statistics and use them to select and study variable sources. We find that 2% of unresolved optical sources brighter than g=20.5 appear variable at the 0.05 mag level (rms) simultaneously in the g and r bands. The majority (2/3) of these variable sources are low-redshift (<2) quasars, although they represent only 2% of all sources in the adopted flux-limited sample. We find that at least 90% of quasars are variable at the 0.03 mag level (rms) and confirm that variability is as good a method for finding low-redshift quasars as is the UV excess color selection (at high Galactic latitudes). We analyze the distribution of lightcurve skewness for quasars and find that is centered on zero. We find that about 1/4 of the variable stars are RR Lyrae stars, and that only 0.5% of stars from the main stellar locus are variable at the 0.05 mag level. The distribution of lightcurve skewness in the g-r vs. u-g color-color diagram on the main stellar locus is found to be bimodal (with one mode consistent with Algol-like behavior). Using over six hundred RR Lyrae stars, we demonstrate rich halo substructure out to distances of 100 kpc. We extrapolate these results to expected performance by the Large Synoptic Survey Telescope and estimate that it will obtain well-sampled 2% accurate, multi-color lightcurves for ~2 million low-redshift quasars, and will discover at least 50 million variable stars.

We present a new estimator, omega, of the small scale galaxy correlation function that is robust against the effects of redshift space distortions and large scale structures. The estimator is a weighted integral of the redshift space or angular correlation function and is a convolution of the real space correlation function with a localized filter. This allows a direct comparison with theory, without modeling redshift space distortions and the large scale correlation function. This has a number of advantages over the more traditional w_p estimator, including (i) an insensitivity to large scale structures and the details of the truncation of the line of sight integral, (ii) a compact kernel in xi(r), and (iii) being unbinned. We discuss how this estimator would be used in practice, applying it to a sample of mock galaxies selected from the Millennium simulation.

We compare Stetson's photometric standards with measurements listed in a standard star catalog constructed using repeated SDSS imaging observations. The SDSS catalog includes over 700,000 candidate standard stars from the equatorial stripe 82 (|Dec|<1.266 deg) in the RA range 20h 34' to 4h 00', and with the $r$ band magnitudes in the range 14--21. The distributions of measurements for individual sources demonstrate that the SDSS photometric pipeline correctly estimates random photometric errors, which are below 0.01 mag for stars brighter than (19.5, 20.5, 20.5, 20, 18.5) in ugriz, respectively (about twice as good as for individual SDSS runs). We derive mean photometric transformations between the SDSS gri and the BVRI system using 1165 Stetson stars found in the equatorial stripe 82, and then study the spatial variation of the difference in zeropoints between the two catalogs. Using third order polynomials to describe the color terms, we find that photometric measurements for main-sequence stars can be transformed between the two systems with systematic errors smaller than a few millimagnitudes. The spatial variation of photometric zeropoints in the two catalogs typically does not exceed 0.01 magnitude. Consequently, the SDSS Standard Star Catalog for Stripe 82 can be used to calibrate new data in both the SDSS ugriz and the BVRI systems with a similar accuracy.

We measure the large-scale real-space power spectrum P(k) using luminous red galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS) and use this measurement to sharpen constraints on cosmological parameters from the Wilkinson Microwave Anisotropy Probe (WMAP). We employ a matrix-based power spectrum estimation method using Pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 20 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.01h/Mpc < k < 0.2h/Mpc. Results from the LRG and main galaxy samples are consistent, with the former providing higher signal-to-noise. Our results are robust to omitting angular and radial density fluctuations and are consistent between different parts of the sky. They provide a striking confirmation of the predicted large-scale LCDM power spectrum. Combining only SDSS LRG and WMAP data places robust constraints on many cosmological parameters that complement prior analyses of multiple data sets. The LRGs provide independent cross-checks on Om and the baryon fraction in good agreement with WMAP. Within the context of flat LCDM models, our LRG measurements complement WMAP by sharpening the constraints on the matter density, the neutrino density and the tensor amplitude by about a factor of two, giving Omega_m=0.24+-0.02 (1 sigma), sum m_nu < 0.9 eV (95%) and r<0.3 (95%). Baryon oscillations are clearly detected and provide a robust measurement of the comoving distance to the median survey redshift z=0.35 independent of curvature and dark energy properties. Within the LCDM framework, our power spectrum measurement improves the evidence for spatial flatness, sharpening the curvature constraint Omega_tot=1.05+-0.05 from WMAP alone to Omega_tot=1.003+-0.010. Assuming Omega_tot=1, the equation of state parameter is constrained to w=-0.94+-0.09, indicating the potential for more ambitious future LRG measurements to provide precision tests of the nature of dark energy. All these constraints are essentially independent of scales k>0.1h/Mpc and associated nonlinear complications, yet agree well with more aggressive published analyses where nonlinear modeling is crucial.

Previous calculations of the pregalactic chemistry have found that a small amount of H_2, x[H_2]=n[H_2]/n[H] = 2.6e-6, is produced catalytically through the H^-, H_2^+, and HeH^+ mechanisms. We revisit this standard calculation taking into account the effects of the nonthermal radiation background produced by cosmic hydrogen recombination, which is particularly effective at destroying H^- via photodetachment. We also take into consideration the non-equilibrium level populations of H_2^+, which occur since transitions among the rotational-vibrational levels are slow compared to photodissociation. The new calculation predicts a final H_2 abundance of x[H_2] = 6e-7 for the standard cosmology. This production is due almost entirely to the H^- mechanism, with ~1 per cent coming from HeH^+ and ~0.004 per cent from H_2^+. We evaluate the heating of the diffuse pregalactic gas from the chemical reactions that produce H_2 and from rotational transitions in H_2, and find them to be negligible.

We present the 3D real space clustering power spectrum of a sample of \~600,000 luminous red galaxies (LRGs) measured by the Sloan Digital Sky Survey (SDSS), using photometric redshifts. This sample of galaxies ranges from redshift z=0.2 to 0.6 over 3,528 deg^2 of the sky, probing a volume of 1.5 (Gpc/h)^3, making it the largest volume ever used for galaxy clustering measurements. We measure the angular clustering power spectrum in eight redshift slices and combine these into a high precision 3D real space power spectrum from k=0.005 (h/Mpc) to k=1 (h/Mpc). We detect power on gigaparsec scales, beyond the turnover in the matter power spectrum, on scales significantly larger than those accessible to current spectroscopic redshift surveys. We also find evidence for baryonic oscillations, both in the power spectrum, as well as in fits to the baryon density, at a 2.5 sigma confidence level. The statistical power of these data to constrain cosmology is ~1.7 times better than previous clustering analyses. Varying the matter density and baryon fraction, we find \Omega_M = 0.30 \pm 0.03, and \Omega_b/\Omega_M = 0.18 \pm 0.04, The detection of baryonic oscillations also allows us to measure the comoving distance to z=0.5; we find a best fit distance of 1.73 \pm 0.12 Gpc, corresponding to a 6.5% error on the distance. These results demonstrate the ability to make precise clustering measurements with photometric surveys (abridged).

Weak lensing is emerging as a powerful observational tool to constrain cosmological models, but is at present limited by an incomplete understanding of many sources of systematic error. Many of these errors are multiplicative and depend on the population of background galaxies. We show how the commonly cited geometric test, which is rather insensitive to cosmology, can be used as a ratio test of systematics in the lensing signal at the 1 per cent level. We apply this test to the galaxy-galaxy lensing analysis of the Sloan Digital Sky Survey (SDSS), which at present is the sample with the highest weak lensing signal to noise and has the additional advantage of spectroscopic redshifts for lenses. This allows one to perform meaningful geometric tests of systematics for different subsamples of galaxies at different mean redshifts, such as brighter galaxies, fainter galaxies and high-redshift luminous red galaxies, both with and without photometric redshift estimates. We use overlapping objects between SDSS and the DEEP2 and 2SLAQ spectroscopic surveys to establish accurate calibration of photometric redshifts and to determine the redshift distributions for SDSS. We use these redshift results to compute the projected surface density contrast DeltaSigma around 259 609 spectroscopic galaxies in the SDSS; by measuring DeltaSigma with different source samples we establish consistency of the results at the 10 per cent level (1-sigma). We also use the ratio test to constrain shear calibration biases and other systematics in the SDSS survey data to determine the overall galaxy-galaxy weak lensing signal calibration uncertainty. We find no evidence of any inconsistency among many subsamples of the data.

We discuss the construction of a photometric redshift catalogue of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS), emphasizing the principal steps necessary for constructing such a catalogue -- (i) photometrically selecting the sample, (ii) measuring photometric redshifts and their error distributions, (iii) and estimating the true redshift distribution. We compare two photometric redshift algorithms for these data and find that they give comparable results. Calibrating against the SDSS and SDSS-2dF spectroscopic surveys, we find that the photometric redshift accuracy is $\sigma \sim 0.03$ for redshifts less than 0.55 and worsens at higher redshift ($\sim 0.06$). These errors are caused by photometric scatter, as well as systematic errors in the templates, filter curves, and photometric zeropoints. We also parametrize the photometric redshift error distribution with a sum of Gaussians, and use this model to deconvolve the errors from the measured photometric redshift distribution to estimate the true redshift distribution. We pay special attention to the stability of this deconvolution, regularizing the method with a prior on the smoothness of the true redshift distribution. The methods we develop are applicable to general photometric redshift surveys.

Dark matter (DM) annihilation during hydrogen recombination (z ~ 1000) will alter the recombination history of the Universe, and affect the observed CMB temperature and polarization fluctuations. Unlike other astrophysical probes of DM, this is free of the significant uncertainties in modelling galactic physics, and provides a method to detect and constrain the cosmological abundances of these particles. We parametrize the effect of DM annihilation as an injection of ionizing energy at a rate e_{dm}, and argue that this simple "on the spot'' modification is a good approximation to the complicated interaction of the annihilation products with the photon-electron plasma. Generic models of DM do not change the redshift of recombination, but change the residual ionization after recombination. This broadens the surface of last scattering, suppressing the temperature fluctuations and enhancing the polarization fluctuations. We use the temperature and polarization angular power spectra to measure these deviations from the standard recombination history, and therefore, indirectly probe DM annihilation. (abridged)

We combine the measurements of luminosity dependence of bias with the luminosity dependent weak lensing analysis of dark matter around galaxies to derive the galaxy bias and constrain nonlinear mass and cosmological parameters. We take advantage of theoretical and simulation predictions that predict that while halo bias is rapidly increasing with mass for high masses, it is nearly constant in low mass halos. We use a new weak lensing analysis around the same SDSS galaxies to determine their halo mass probability distribution. These halo mass probability distributions are used to predict the bias for each luminosity subsample and we find an excellent agreement with observed values. The required nonlinear mass suggests slightly lower matter density than usually assumed, Omegam=0.25+/- 0.03 for the simplest models. We combine the bias constraints with those from the WMAP and the SDSS power spectrum analysis to derive new constraints on bias and sigma_8. For the most general parameter space we find sigma_8=0.88+/- 0.06 and b_*=0.99+/- 0.07. In the context of spatially flat models we improve the limit on the neutrino mass for the case of 3 degenerate families from m_nu<0.6eV without bias to m_nu<0.18eV with bias (95% c.l.), which is weakened to m_nu<0.24eV if running is allowed. The corresponding limit for 3 massless + 1 massive neutrino is 1.37eV.

We present a 2.5 sigma detection of the Integrated Sachs-Wolfe (ISW) effect and discuss the constraints it places on cosmological parameters. We cross-correlate microwave temperature maps from the WMAP satellite with a 4000 deg^2 luminous red galaxy (LRG) overdensity map measured by the Sloan Digital Sky Survey. Accurate photometric redshifts allow us to perform a reliable auto-correlation analysis of the LRGs, eliminating the uncertainty in the galaxy bias, and combined with cross correlation signal, constrains cosmological parameters -- in particular, the matter density. We find a 2.5 sigma signal in the Ka, Q, V, and W WMAP bands, after combining the information from multipoles 2 <= l < 400. This is consistent with the expected amplitude of the ISW effect, but requires a lower matter density than is usually assumed: the amplitude, parametrized by the galaxy bias assuming \Omega_M=0.3, \Omega_\Lambda=0.7 and \sigma_8=0.9, is b_g = 4.05 \pm 1.54 for V band, with similar results for the other bands. This should be compared to b_g = 1.82 \pm 0.02 from the auto-correlation analysis. These data provide only a weak confirmation (2.5 sigma) of dark energy, but provide a significant upper limit: \Omega_\Lambda=0.80_{-0.06}^{+0.03} (1 sigma)_{-0.19}^{+0.05} (2 sigma), assuming a cosmology with \Omega_M+\Omega_\Lambda=1, \Omega_b = 0.05, and \sigma_8=0.9, and w=-1. The weak cross-correlation signal rules out low matter density/high dark energy density universes and, in combination with other data, strongly constrains models with w<-1.3. We provide a simple prescription to incorporate these constraints into cosmological parameter estimation methods for (\Omega_M, \sigma_8,w). (abridged)

We summarize the Sloan Digital Sky Survey data acquisition and processing steps, and describe runQA, a pipeline designed for automated data quality assessment. In particular, we show how the position of the stellar locus in color-color diagrams can be used to estimate the accuracy of photometric zeropoint calibration to better than 0.01 mag in 0.03 deg2 patches. Using this method, we estimate that typical photometric zeropoint calibration errors for SDSS imaging data are not larger than ~0.01 mag in the g, r, and i bands, 0.02 mag in the z band, and 0.03 mag in the u band (root-mean-scatter for zeropoint offsets).

The Sloan Digital Sky Survey (SDSS) mosaic camera and telescope have obtained five-band optical-wavelength imaging near the Galactic plane outside of the nominal survey boundaries. These additional data were obtained during commissioning and subsequent testing of the SDSS observing system, and they provide unique wide-area imaging data in regions of high obscuration and star formation, including numerous young stellar objects, Herbig-Haro objects and young star clusters. Because these data are outside the Survey regions in the Galactic caps, they are not part of the standard SDSS data releases. This paper presents imaging data for 832 square degrees of sky (including repeats), in the star-forming regions of Orion, Taurus, and Cygnus. About 470 square degrees are now released to the public, with the remainder to follow at the time of SDSS Data Release 4. The public data in Orion include the star-forming region NGC 2068/NGC 2071/HH24 and a large part of Barnard's loop.

We present the results of a search for gravitational lensing of the cosmic microwave background (CMB) in cross-correlation with the projected density of luminous red galaxies (LRGs). The CMB lensing reconstruction is performed using the first year of Wilkinson Microwave Anisotropy Probe (WMAP) data, and the galaxy maps are obtained using the Sloan Digital Sky Survey (SDSS) imaging data. We find no detection of lensing; our constraint on the galaxy bias derived from the galaxy-convergence cross-spectrum is $b_g=1.81\pm 1.92$ ($1\sigma$, statistical), as compared to the expected result of $b_g\sim 1.8$ for this sample. We discuss possible instrument-related systematic errors and show that the Galactic foregrounds are not important. We do not find any evidence for point source or thermal Sunyaev-Zel'dovich effect contamination.

The Spitzer Space Telescope First Look Survey Infrared Array Camera (IRAC) near and mid-infrared imaging data partially overlaps the Sloan Digital Sky Survey (SDSS), with 313 visually selected (r<17.6 mag) SDSS Main Sample galaxies in the overlap region. The 3.5 and 7.8 um properties of the galaxies are investigated in the context of their visual properties, where the IRAC [3.5] magnitude primarily measures starlight, and the [7.8] magnitude primarily measures PAH emission from the interstellar medium. As expected, we find a strong inverse correlation between [3.5]-[7.8] and visual color; galaxies red in visual colors (`red galaxies') tend to show very little dust and molecular emission (low `PAH-to-star' ratios), and galaxies blue in visual colors (`blue galaxies,' ie, star-forming galaxies) tend to show large PAH-to-star ratios. Red galaxies with high PAH-to-star ratios tend to be edge-on disks reddened by dust lanes. Simple, visually inferred attenuation corrections bring the visual colors of these galaxies in line with those of face-on disks; ie, PAH emission is closely related to attenuation-corrected, optically inferred star-formation rates. Blue galaxies with anomalously low PAH-to-star ratios are all low-luminosity star-forming galaxies. There is some weak evidence in this sample that the deficiency in PAH emission for these low-luminosity galaxies may be related to emission-line metallicity.

Galaxy-galaxy lensing has emerged as a powerful probe of the dark matter halos of galaxies, but is subject to contamination if intrinsically aligned satellites of the lens galaxy are used as part of the source sample. We present a measurement of this intrinsic shear using 200,747 lens galaxies from the Sloan Digital Sky Survey (SDSS) spectroscopic sample and a sample of satellites selected using photometric redshifts. The mean intrinsic shear at transverse separations of 30--446$h^{-1}$ kpc is constrained to be $-0.0062<\Delta\gamma<+0.0066$ (99.9 per cent confidence, including identified systematics), which limits contamination of the galaxy-galaxy lensing signal to at most $\sim 15$ per cent on these scales. We present these limits as a function of transverse separation and lens luminosity. We furthermore investigate shear calibration biases in the SDSS, which can also affect galaxy-galaxy lensing, and conclude that the shear amplitude is calibrated to better than 18 per cent. This includes noise-induced calibration biases in the ellipticity, which are small for the sample considered here, but which can be more important if low signal-to-noise or poorly resolved source galaxies are used.

We study the variation of the dark matter mass fraction of elliptical galaxies as a function of their luminosity, stellar mass, and size using a sample of 29,469 elliptical galaxies culled from the Sloan Digital Sky Survey. We model ellipticals as a stellar Hernquist profile embedded in an adiabatically compressed dark matter halo. This model allows us to estimate a dynamical mass ($M_{dynm}$) at the half-light radius from the velocity dispersion of the spectra, and to compare these to the stellar mass estimates ($M_{*}$) from Kauffmann et al (2003). We find that $M_{*}/L$ is independent of luminosity, while $M_{dynm}/L$ increases with luminosity, implying that the dark matter fraction increases with luminosity. We also observe that at a fixed luminosity or stellar mass, the dark matter fraction increases with increasing galaxy size or, equivalently, increases with decreasing surface brightness: high surface brightness galaxies show almost no evidence for dark matter, while in low surface brightness galaxies, the dark matter exceeds the stellar mass at the half light radius. We relate this to the fundamental plane of elliptical galaxies, suggesting that the tilt of this plane from simple virial predictions is due to the dark matter in galaxies. We find that a simple model where galaxies are embedded in dark matter halos and have a star formation efficiency independent of their surface brightness explains these trends. We estimate the virial mass of ellipticals as being approximately 7-30 times their stellar mass, with the lower limit suggesting almost all of the gas within the virial radius is converted into stars.

We develop a new method to constrain the star formation histories, dust attenuation and stellar masses of galaxies. It is based on two stellar absorption line indices, the 4000 Angstrom break strength and the Balmer absorption line index Hdelta_A. Together, these indices allow us to constrain the mean stellar ages of galaxies and the fractional stellar mass formed in bursts over the past few Gyr. A comparison with broad band photometry then yields estimates of dust attenuation and of stellar mass. We generate a large library of Monte Carlo realizations of different star formation histories, including starbursts of varying strength and a range of metallicities. We use this library to generate median likelihood estimates of burst mass fractions, dust attenuation strengths, stellar masses and stellar mass-to-light ratios for a sample of 122,208 galaxies drawn from the Sloan Digital Sky Survey. The typical 95% confidence range in our estimated stellar masses is +-40%. We study how the stellar mass-to-light ratios of galaxies vary as a function of absolute magnitude, concentration index and photometric pass-band and how dust attenuation varies as a function of absolute magnitude and 4000 Angstrom break strength. We also calculate how the total stellar mass of the present Universe is distributed over galaxies as a function of their mass, size, concentration, colour, burst mass fraction and surface mass density. We find that most of the stellar mass in the local Universe resides in galaxies that have stellar masses \~5\times 10^10 M_sol, half light radii ~3 kpc, and half- light surface mass densities ~10^9 M_sol/kpc^2. The distribution of D(4000) is strongly bimodal, showing a clear division between galaxies dominated by old stellar populations and galaxies with more recent star formation.

We present a survey of the cosmological applications of the next generation of weak lensing surveys, paying special attention to the computational challenges presented by the number of galaxies, $N_{gal} ~$ 10$^{5}$. We focus on optimal methods with no pixelization and derive a multigrid $P^3M$ algorithm that performs the relevant computations in $O(N_{gal} \log N_{gal})$ time. We test the algorithm by studying three applications of weak lensing surveys - convergence map reconstruction, cluster detection and $E$ and $B$ power spectrum estimation using realistic 1 deg^{2} simulations derived from N-body simulations. The map reconstruction is able to reconstruct large scale features without artifacts. Detecting clusters using only weak lensing is difficult because of line of sight contamination and noise, with low completeness if one desires low contamination of the sample. A power spectrum analysis of the convergence field is more promising and we are able to reconstruct the convergence spectrum with no loss of information down to the smallest scales. The numerical methods used here can be applied to other data sets with same $O(N\log N)$ scaling and can be generalised to a sphere.

The combined CfA2/SSRS redshift catalog has recently been improved and made public. We compute the redshift-space power spectrum of this catalog using the regions at -2.5 < delta < 50 north (8h < alpha < 17h) and south (20h < alpha < 3h, b < -13) of the Galactic plane, where it is 98% complete down to Zwicky magnitude 15.5 and contains 13,681 galaxies. Our analysis uses Heavens-Taylor mode expansion, Karhunen-Loeve data compression and the Fisher matrix technique to compute quadratic band power estimates. This allows an exact calculation of window functions, including the integral constraint, in addition to the production of a power spectrum with uncorrelated error bars. Our results with this larger data set agree well with previous studies. We analyze 101Mpc/h and 130Mpc/h volume-limited subsets in addition to the full magnitude-limited sample, and our results are well fit by, e.g., Lambda-CDM models with bias b=1.2, 1.4 and 1.4, respectively. We estimate the effect of extinction using the Schlegel, Finkbeiner & Davis dust map. Our results are exclusively for the redshift space galaxy power spectrum, so they can only be compared with theoretical predictions if appropriate corrections are made for biasing and redshift space distortions.

We study the reliability of dark matter (DM) profile determinations for rich clusters of galaxies based on weak gravitational lensing measurements. We consider three functional forms for the mass distribution: a softened isothermal sphere (SIS), the Navarro, Frenk, & White model (NFW) and the Klypin, Kravstov, Bullock, & Primack model (KKBP). A maximum likelihood parametric method is developed and applied to simulated images of background galaxies (z=2.0) lensed by clusters at z=0.2, 0.5, and 0.8. The images have field size (~10 arcmin) and angular resolution (~0.5 arcsec) characteristic of the new generation of 8-10 meter-class ground-based telescopes. The method recovers the correct virial mass of the DM halo, M_{200}, to within ~10% at all redshifts. The recovered values of \delta_c and r_s are only accurate to ~25%; however, their correlated errors are such that greater accuracy is obtained for the mass. Clusters with an SIS density profile are recognized by the method, while NFW and KKBP clusters are recognized as not being SIS. The method does not discriminate between NFW and KKBP profiles for clusters of a given mass, because they are too similar in their outer regions. Strongly lensed arcs near the cluster centers would be needed to distinguish NFW from KKBP.

We present a state-of-the-art linear redshift distortion analysis of the recently published IRAS Point Source Catalog Redshift Survey (PSCz). The procedure involves linear compression into 4096 Karhunen-Loeve modes culled from a potential pool of about 3 x 10^5 modes, followed by quadratic compression into three separate power spectra, the galaxy-galaxy, galaxy-velocity, and velocity-velocity power spectra. Least squares fitting to the decorrelated power spectra yields a linear redshift distortion parameter beta = Omega_m^0.6/b = 0.41(+0.13,-0.12).