Source author record

Martin White

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

ResearcherUnclaimed source record

astro-ph.CO hep-ph astro-ph.IM astro-ph.GA astro-ph hep-ex astro-ph.HE physics.data-an Methodology

Catalog footprint

What is connected

150works

9topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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

Building this map preview

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

preprint2023arXiv

Second Data Release of the COSMOS Lyman-alpha Mapping and Tomographic Observation: The First 3D Maps of the Detailed Cosmic Web at 2.05<z<2.55

We present the second data release of the COSMOS Lyman-Alpha Mapping And Tomography Observations (CLAMATO) Survey conducted with the LRIS spectrograph on the Keck-I telescope. This project used Lyman-alpha forest absorption in the spectra of faint star forming galaxies and quasars at z ~ 2-3 to trace neutral hydrogen in the intergalactic medium. In particular, we use 320 objects over a footprint of ~0.2 deg^2 to reconstruct the absorption field at 2.05 < z < 2.55 at ~2 h^{-1}Mpc resolution. We apply a Wiener filtering technique to the observed data to reconstruct three dimensional maps of the field over a volume of 4.1 x 10^5 comoving cubic Mpc. In addition to the filtered flux maps, for the first time we infer the underlying dark matter field through a forward modeling framework from a joint likelihood of galaxy and Lyman-alpha forest data, finding clear examples of the detailed cosmic web consisting of cosmic voids, sheets, filaments, and nodes. In addition to traditional figures, we present a number of interactive three dimensional models to allow exploration of the data and qualitative comparisons to known galaxy surveys. We find that our inferred over-densities are consistent with those found from galaxy fields. Our reduced spectra, extracted Lyman-alpha forest pixel data, and reconstructed tomographic maps are available publicly at https://doi.org/10.5281/zenodo.7524313

preprint2022arXiv

A new analysis of galaxy 2-point functions in the BOSS survey, including full-shape information and post-reconstruction BAO

We present a new method for consistent, joint analysis of the pre- and post-reconstruction two-point functions of the BOSS survey. The post-reconstruction correlation function is used to accurately measure the distance-redshift relation and expansion history, while the pre-reconstruction power spectrum multipoles constrain the broad-band shape and the rate-of-growth of large-scale structure. Our technique uses Lagrangian perturbation theory to self-consistently work at the level of two-point functions, i.e.\ directly with the measured data, without approximating the constraints with summary statistics normalized by the drag scale. Combining galaxies across the full redshift range and both hemispheres we constrain $Ω_m=0.303 \pm 0.0082$, $H_0=69.23 \pm 0.77$ and $σ_8=0.733 \pm 0.047$ within the context of $Λ$CDM. These constraints are in good agreement both with the Planck primary CMB anisotropy data and recent cosmic shear surveys.

preprint2022arXiv

A Spectroscopic Road Map for Cosmic Frontier: DESI, DESI-II, Stage-5

In this white paper, we present an experimental road map for spectroscopic experiments beyond DESI. DESI will be a transformative cosmological survey in the 2020s, mapping 40 million galaxies and quasars and capturing a significant fraction of the available linear modes up to z=1.2. DESI-II will pilot observations of galaxies both at much higher densities and extending to higher redshifts. A Stage-5 experiment would build out those high-density and high-redshift observations, mapping hundreds of millions of stars and galaxies in three dimensions, to address the problems of inflation, dark energy, light relativistic species, and dark matter. These spectroscopic data will also complement the next generation of weak lensing, line intensity mapping and CMB experiments and allow them to reach their full potential.

preprint2022arXiv

Cosmological Analysis of Three-Dimensional BOSS Galaxy Clustering and Planck CMB Lensing Cross Correlations via Lagrangian Perturbation Theory

We present a formalism for jointly fitting pre- and post-reconstruction redshift-space clustering (RSD) and baryon acoustic oscillations (BAO) plus gravitational lensing (of the CMB) that works directly with the observed 2-point statistics. The formalism is based upon (effective) Lagrangian perturbation theory and a Lagrangian bias expansion, which models RSD, BAO and galaxy-lensing cross correlations within a consistent dynamical framework. As an example we present an analysis of clustering measured by the Baryon Oscillation Spectroscopic Survey in combination with CMB lensing measured by Planck. The post-reconstruction BAO strongly constrains the distance-redshift relation, the full-shape redshift-space clustering constrains the matter density and growth rate, and CMB lensing constrains the clustering amplitude. Using only the redshift space data we obtain $Ω_\mathrm{m} = 0.303\pm 0.008$, $H_0 = 69.21\pm 0.78$ and $σ_8 = 0.743\pm 0.043$. The addition of lensing information, even when restricted to the Northern Galactic Cap, improves constraints to $Ω_m = 0.300 \pm 0.008$, $H_0 = 69.21 \pm 0.77$ and $σ_8 = 0.707 \pm 0.035$, in tension with CMB and cosmic shear constraints. The combination of $Ω_m$ and $H_0$ are consistent with Planck, though their constraints derive mostly from redshift-space clustering. The low $σ_8$ value are driven by cross correlations with CMB lensing in the low redshift bin ($z\simeq 0.38$) and at large angular scales, which show a $20\%$ deficit compared to expectations from galaxy clustering alone. We conduct several systematics tests on the data and find none that could fully explain these tensions.

preprint2022arXiv

Cosmological constraints from the tomographic cross-correlation of DESI Luminous Red Galaxies and Planck CMB lensing

We use luminous red galaxies selected from the imaging surveys that are being used for targeting by the Dark Energy Spectroscopic Instrument (DESI) in combination with CMB lensing maps from the Planck collaboration to probe the amplitude of large-scale structure over $0.4\le z\le 1$. Our galaxy sample, with an angular number density of approximately $500\,\mathrm{deg}^{-2}$ over 18,000 sq.deg., is divided into 4 tomographic bins by photometric redshift and the redshift distributions are calibrated using spectroscopy from DESI. We fit the galaxy autospectra and galaxy-convergence cross-spectra using models based on cosmological perturbation theory, restricting to large scales that are expected to be well described by such models. Within the context of $Λ$CDM, combining all 4 samples and using priors on the background cosmology from supernova and baryon acoustic oscillation measurements, we find $S_8=σ_8(Ω_m/0.3)^{0.5}=0.73\pm 0.03$. This result is lower than the prediction of the $Λ$CDM model conditioned on the Planck data. Our data prefer a slower growth of structure at low redshift than the model predictions, though at only modest significance.

preprint2022arXiv

Cosmology at high redshift -- a probe of fundamental physics

An observational program focused on the high redshift ($2<z<6$) Universe has the opportunity to dramatically improve over upcoming LSS and CMB surveys on measurements of both the standard cosmological model and its extensions. Using a Fisher matrix formalism that builds upon recent advances in Lagrangian perturbation theory, we forecast constraints for future spectroscopic and 21-cm surveys on the standard cosmological model, curvature, neutrino mass, relativistic species, primordial features, primordial non-Gaussianity, dynamical dark energy, and gravitational slip. We compare these constraints with those achievable by current or near-future surveys such as DESI and Euclid, all under the same forecasting formalism, and compare our formalism with traditional linear methods. Our Python code FishLSS $-$ used to calculate the Fisher information of the full shape power spectrum, CMB lensing, the cross-correlation of CMB lensing with galaxies, and combinations thereof $-$ is publicly available.

preprint2022arXiv

Global fit of 2HDM with future collider results

In this work, we summarize a global fit study of Type-II two Higgs doublet models (2HDM), and explore the impact of future SM-like Higgs and Z-pole precision measurements on the allowed parameter space. The work is based on the study results of a global fit of 2HDMs with the tool GAMBIT, utilising various current constraints including theoretical constraints (unitarity, perturbativity and vacuum stability), Higgs searches at colliders, electroweak physics and flavour constraints. We further investigate the ability of future facilities, such as the HL-LHC, CEPC, ILC and FCC-ee to explore the 2HDM parameter space.

preprint2022arXiv

Priors on red galaxy stochasticity from hybrid effective field theory

We investigate the stochastic properties of typical red galaxy samples in a controlled numerical environment. We use Halo Occupation Distribution (HOD) modelling to create mock realizations of three separate bright red galaxy samples consistent with datasets used for clustering and lensing analyses in modern galaxy surveys. Second-order Hybrid Effective Field Theory (HEFT) is used as a field-level forward model to describe the full statistical distribution of these tracer samples, and their stochastic power spectra are directly measured and compared to the Poisson shot-noise prediction. While all of the galaxy samples we consider are hosted within haloes with sub-Poisson stochasticity, we observe that the galaxy samples themselves possess stochasticities that range from sub-Poisson to super-Poisson, in agreement with predictions from the halo model. As an application of our methodology, we place priors on the expected degree of non-Poisson stochasticity in cosmological analyses using such samples. We expect these priors will be useful in reducing the complexity of the full parameter space for future analyses using second-order Lagrangian bias models. More generally, the techniques outlined here present the first application of hybrid EFT methods to characterize models of the galaxy--halo connection at the field level, revealing new connections between once-disparate modelling frameworks.

preprint2022arXiv

Simple and statistically sound recommendations for analysing physical theories

Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Our examples can be reproduced with the code publicly available at https://doi.org/10.5281/zenodo.4322283.

preprint2022arXiv

Snowmass 2021 CMB-S4 White Paper

This Snowmass 2021 White Paper describes the Cosmic Microwave Background Stage 4 project CMB-S4, which is designed to cross critical thresholds in our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. We provide an overview of the science case, the technical design, and project plan.

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Cosmology and Fundamental Physics from the three-dimensional Large Scale Structure

Advances in experimental techniques make it possible to map the high redshift Universe in three dimensions at high fidelity in the near future. This will increase the observed volume by many-fold, while providing unprecedented access to very large scales, which hold key information about primordial physics. Recently developed theoretical techniques, together with the smaller size of non-linearities at high redshift, allow the reconstruction of an order of magnitude more "primordial modes", and should improve our understanding of the early Universe through measurements of primordial non-Gaussianity and features in the primordial power spectrum. In addition to probing the first epoch of accelerated expansion, such measurements can probe the Dark Energy density in the dark matter domination era, tightly constraining broad classes of dynamical Dark Energy models. The shape of the matter power spectrum itself has the potential to detect sub-percent fractional amounts of Early Dark Energy to $z \sim 10^5$, probing Dark Energy all the way to when the Universe was only a few years old. The precision of these measurements, combined with CMB observations, also has the promise of greatly improving our constraints on the effective number of relativistic species, the masses of neutrinos, the amount of spatial curvature and the gravitational slip. Studies of linear or quasi-linear large-scale structure with redshift surveys and the CMB currently provide our tightest constraints on cosmology and fundamental physics. Pushing the redshift and volume frontier will provide guaranteed, significant improvements in the state-of-the-art in a manner that is easy to forecast and optimize.

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: High Density Galaxy Clustering in the Regime of Cosmic Acceleration

Joint studies of imaging and spectroscopic samples, informed by theory and simulations, offer the potential for comprehensive tests of the cosmological model over redshifts z<1.5. Spectroscopic galaxy samples at these redshifts can be increased beyond the planned Dark Energy Spectroscopic Instrument (DESI) program by at least an order of magnitude, thus offering significantly more constraining power for these joint studies. Spectroscopic observations of these galaxies in the latter half of the 2020's and beyond would leverage the theory and simulation effort in this regime. In turn, these high density observations will allow enhanced tests of dark energy, physics beyond the standard model, and neutrino masses that will greatly exceed what is currently possible. Here, we present a coordinated program of simulations, theoretical modeling, and future spectroscopy that would enable precise cosmological studies in the accelerating epoch where the effects of dark energy are most apparent.

preprint2022arXiv

Snowmass2021: Opportunities from Cross-survey Analyses of Static Probes

Cosmological data in the next decade will be characterized by high-precision, multi-wavelength measurements of thousands of square degrees of the same patches of sky. By performing multi-survey analyses that harness the correlated nature of these datasets, we will gain access to new science, and increase the precision and robustness of science being pursued by each individual survey. However, effective application of such analyses requires a qualitatively new level of investment in cross-survey infrastructure, including simulations, associated modeling, coordination of data sharing, and survey strategy. The scientific gains from this new level of investment are multiplicative, as the benefits can be reaped by even present-day instruments, and can be applied to new instruments as they come online.

preprint2022arXiv

The DESI $N$-body Simulation Project -- II. Suppressing sample variance with fast simulations

Dark Energy Spectroscopic Instrument (DESI) will construct a large and precise three-dimensional map of our Universe. The survey effective volume reaches $\sim20\Gpchcube$. It is a great challenge to prepare high-resolution simulations with a much larger volume for validating the DESI analysis pipelines. \textsc{AbacusSummit} is a suite of high-resolution dark-matter-only simulations designed for this purpose, with $200\Gpchcube$ (10 times DESI volume) for the base cosmology. However, further efforts need to be done to provide a more precise analysis of the data and to cover also other cosmologies. Recently, the CARPool method was proposed to use paired accurate and approximate simulations to achieve high statistical precision with a limited number of high-resolution simulations. Relying on this technique, we propose to use fast quasi-$N$-body solvers combined with accurate simulations to produce accurate summary statistics. This enables us to obtain 100 times smaller variance than the expected DESI statistical variance at the scales we are interested in, e.g. $k < 0.3\hMpc$ for the halo power spectrum. In addition, it can significantly suppress the sample variance of the halo bispectrum. We further generalize the method for other cosmologies with only one realization in \textsc{AbacusSummit} suite to extend the effective volume $\sim 20$ times. In summary, our proposed strategy of combining high-fidelity simulations with fast approximate gravity solvers and a series of variance suppression techniques sets the path for a robust cosmological analysis of galaxy survey data.

preprint2022arXiv

The MegaMapper: A Stage-5 Spectroscopic Instrument Concept for the Study of Inflation and Dark Energy

In this white paper, we present the MegaMapper concept. The MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at $2<z<5$. In order to achieve path-breaking results with a mid-scale investment, the MegaMapper combines existing technologies for critical path elements and pushes innovative development in other design areas. To this aim, we envision a 6.5-m Magellan-like telescope, with a newly designed wide field, coupled with DESI spectrographs, and small-pitch robots to achieve multiplexing of at least 26,000. This will match the expected achievable target density in the redshift range of interest and provide a 10x capability over the existing state-of the art, without a 10x increase in project budget.

preprint2021arXiv

Convergent Bayesian Global Fits of 4D Composite Higgs Models

Models in which the Higgs boson is a composite pseudo-Nambu-Goldstone boson offer attractive solutions to the Higgs mass naturalness problem. We consider three such models based on the minimal $SO(5) \rightarrow SO(4)$ symmetry breaking pattern, and perform convergent global fits on the models under a Bayesian framework in order to find the regions of their parameter spaces that best fit a wide range of constraints, including recent Higgs measurements. We use a novel technique to analyse the fine-tuning of the models, quantifying the tuning as the Kullback-Leibler divergence from the prior to the posterior probability on the parameter space. Each model is found to be able to satisfy all constraints at the $3σ$ level simultaneously. As a by-product of the fits, we analyse the collider phenomenology of our models in these viable regions. In two of the three models, we find that the $g g \rightarrow H \rightarrow γγ$ cross section is less than ${\sim}90$% that predicted by the SM, which is already in slight tension with experiment and could potentially be ruled out in the future high-luminosity run of the LHC. In addition, the lightest fermions $F$ arising from the new strong dynamics in these models are seen in general to lie above ${\sim}1.1$ TeV, with the $F \rightarrow tW^{+}$ and $F \rightarrow \bar{b}W^{+}$ decays offering particularly promising channels for probing these models in future collider searches.

preprint2021arXiv

CosmoBit: A GAMBIT module for computing cosmological observables and likelihoods

We introduce $\sf{CosmoBit}$, a module within the open-source $\sf{GAMBIT}$ software framework for exploring connections between cosmology and particle physics with joint global fits. $\sf{CosmoBit}$ provides a flexible framework for studying various scenarios beyond $Λ$CDM, such as models of inflation, modifications of the effective number of relativistic degrees of freedom, exotic energy injection from annihilating or decaying dark matter, and variations of the properties of elementary particles such as neutrino masses and the lifetime of the neutron. Many observables and likelihoods in $\sf{CosmoBit}$ are computed via interfaces to $\sf{AlterBBN}$, $\sf{CLASS}$, $\sf{DarkAges}$, $\sf{MontePython}$, $\sf{MultiModeCode}$, and $\sf{plc}$. This makes it possible to apply a wide range of constraints from large-scale structure, Type Ia supernovae, Big Bang Nucleosynthesis and the cosmic microwave background. Parameter scans can be performed using the many different statistical sampling algorithms available within the $\sf{GAMBIT}$ framework, and results can be combined with calculations from other $\sf{GAMBIT}$ modules focused on particle physics and dark matter. We include extensive validation plots and a first application to scenarios with non-standard relativistic degrees of freedom and neutrino temperature, showing that the corresponding constraint on the sum of neutrino masses is much weaker than in the standard scenario.

preprint2021arXiv

Does SUSY have friends? A new approach for LHC event analysis

We present a novel technique for the analysis of proton-proton collision events from the ATLAS and CMS experiments at the Large Hadron Collider. For a given final state and choice of kinematic variables, we build a graph network in which the individual events appear as weighted nodes, with edges between events defined by their distance in kinematic space. We then show that it is possible to calculate local metrics of the network that serve as event-by-event variables for separating signal and background processes, and we evaluate these for a number of different networks that are derived from different distance metrics. Using a supersymmetric electroweakino and stop production as examples, we construct prototype analyses that take account of the fact that the number of simulated Monte Carlo events used in an LHC analysis may differ from the number of events expected in the LHC dataset, allowing an accurate background estimate for a particle search at the LHC to be derived. For the electroweakino example, we show that the use of network variables outperforms both cut-and-count analyses that use the original variables and a boosted decision tree trained on the original variables. The stop example, deliberately chosen to be difficult to exclude due its kinematic similarity with the top background, demonstrates that network variables are not automatically sensitive to BSM physics. Nevertheless, we identify local network metrics that show promise if their robustness under certain assumptions of node-weighted networks can be confirmed.

preprint2021arXiv

Neural Network Acceleration of Large-scale Structure Theory Calculations

We make use of neural networks to accelerate the calculation of power spectra required for the analysis of galaxy clustering and weak gravitational lensing data. For modern perturbation theory codes, evaluation time for a single cosmology and redshift can take on the order of two seconds. In combination with the comparable time required to compute linear predictions using a Boltzmann solver, these calculations are the bottleneck for many contemporary large-scale structure analyses. In this work, we construct neural network-based surrogate models for Lagrangian perturbation theory (LPT) predictions of matter power spectra, real and redshift space galaxy power spectra, and galaxy--matter cross power spectra that attain $\sim 0.1\%$ (at one sigma) accuracy over a broad range of scales in a $w$CDM parameter space. The neural network surrogates can be evaluated in approximately one millisecond, a factor of 1000 times faster than the full Boltzmann code and LPT computations. In a simulated full-shape redshift space galaxy power spectrum analysis, we demonstrate that the posteriors obtained using our surrogates are accurate compared to those obtained using the full LPT model. We make our surrogate models public at https://github.com/sfschen/EmulateLSS, so that others may take advantage of the speed gains they provide to enable rapid iteration on analysis settings, something that is essential in complex contemporary large-scale structure analyses.

preprint2021arXiv

Thermal WIMPs and the Scale of New Physics: Global Fits of Dirac Dark Matter Effective Field Theories

We assess the status of a wide class of WIMP dark matter (DM) models in light of the latest experimental results using the global fitting framework $\textsf{GAMBIT}$. We perform a global analysis of effective field theory (EFT) operators describing the interactions between a gauge-singlet Dirac fermion and the Standard Model quarks, the gluons and the photon. In this bottom-up approach, we simultaneously vary the coefficients of 14 such operators up to dimension 7, along with the DM mass, the scale of new physics and several nuisance parameters. Our likelihood functions include the latest data from $\mathit{Planck}$, direct and indirect detection experiments, and the LHC. For DM masses below 100 GeV, we find that it is impossible to satisfy all constraints simultaneously while maintaining EFT validity at LHC energies. For new physics scales around 1 TeV, our results are influenced by several small excesses in the LHC data and depend on the prescription that we adopt to ensure EFT validity. Furthermore, we find large regions of viable parameter space where the EFT is valid and the relic density can be reproduced, implying that WIMPs can still account for the DM of the universe while being consistent with the latest data.

preprint2020arXiv

Consistent Modeling of Velocity Statistics and Redshift-Space Distortions in One-Loop Perturbation Theory

The peculiar velocities of biased tracers of the cosmic density field contain important information about the growth of large scale structure and generate anisotropy in the observed clustering of galaxies. Using N-body data, we show that velocity expansions for halo redshift-space power spectra are converged at the percent-level at perturbative scales for most line-of-sight angles $μ$ when the first three pairwise velocity moments are included, and that the third moment is well-approximated by a counterterm-like contribution. We compute these pairwise-velocity statistics in Fourier space using both Eulerian and Lagrangian one-loop perturbation theory using a cubic bias scheme and a complete set of counterterms and stochastic contributions. We compare the models and show that our models fit both real-space velocity statistics and redshift-space power spectra for both halos and a mock sample of galaxies at sub-percent level on perturbative scales using consistent sets of parameters, making them appealing choices for the upcoming era of spectroscopic, peculiar-velocity and kSZ surveys.

preprint2020arXiv

Exploring the 2HDM with Global Fits in $\texttt{GAMBIT}$

In this work, we present preliminary results of a global fit of the type-II two-Higgs-doublet model (2HDM) with the tool $\texttt{GAMBIT}$. Our study includes various constraints, including the theoretical constraints (unitarity, perturbativity and vacuum stability), Higgs searches at colliders, electroweak physics and flavour constraints. With the latest experimental results, our results not only confirm past studies but also go further in probing the model. We find, for example, that the measurements of $B\rightarrow K^*μ^+μ^-$ angular observables cannot be explained in the type-II 2HDM.

preprint2020arXiv

GAMBIT and its Application in the Search for Physics Beyond the Standard Model

The Global and Modular Beyond-Standard Model Inference Tool (GAMBIT) is an open source software framework for performing global statistical fits of particle physics models, using a wide range of particle and astroparticle data. In this review, we describe the design principles of the package, the statistical and sampling frameworks, the experimental data included, and the first two years of physics results generated with it. This includes supersymmetric models, axion theories, Higgs portal dark matter scenarios and an extension of the Standard Model to include right-handed neutrinos. Owing to the broad spectrum of physics scenarios tackled by the GAMBIT community, this also serves as a convenient, self-contained review of the current experimental and theoretical status of the most popular models of dark matter.

preprint2020arXiv

Global analyses of Higgs portal singlet dark matter models using GAMBIT

We present global analyses of effective Higgs portal dark matter models in the frequentist and Bayesian statistical frameworks. Complementing earlier studies of the scalar Higgs portal, we use GAMBIT to determine the preferred mass and coupling ranges for models with vector, Majorana and Dirac fermion dark matter. We also assess the relative plausibility of all four models using Bayesian model comparison. Our analysis includes up-to-date likelihood functions for the dark matter relic density, invisible Higgs decays, and direct and indirect searches for weakly-interacting dark matter including the latest XENON1T data. We also account for important uncertainties arising from the local density and velocity distribution of dark matter, nuclear matrix elements relevant to direct detection, and Standard Model masses and couplings. In all Higgs portal models, we find parameter regions that can explain all of dark matter and give a good fit to all data. The case of vector dark matter requires the most tuning and is therefore slightly disfavoured from a Bayesian point of view. In the case of fermionic dark matter, we find a strong preference for including a CP-violating phase that allows suppression of constraints from direct detection experiments, with odds in favour of CP violation of the order of 100:1. Finally, we present DDCalc 2.0.0, a tool for calculating direct detection observables and likelihoods for arbitrary non-relativistic effective operators.

preprint2020arXiv

Imaging Systematics and Clustering of DESI Main Targets

We evaluate the impact of imaging systematics on the clustering of luminous red galaxies (LRG), emission-line galaxies (ELG) and quasars (QSO) targeted for the upcoming Dark Energy Spectroscopic Instrument (DESI) survey. Using Data Release 7 of the DECam Legacy Survey, we study the effects of astrophysical foregrounds, stellar contamination, differences between north galactic cap and south galactic cap measurements, and variations in imaging depth, stellar density, galactic extinction, seeing, airmass, sky brightness, and exposure time before presenting survey masks and weights to mitigate these effects. With our sanitized samples in hand, we conduct a preliminary analysis of the clustering amplitude and evolution of the DESI main targets. From measurements of the angular correlation functions, we determine power law fits $r_0 = 7.78 \pm 0.26$ $h^{-1}$Mpc, $γ= 1.98 \pm 0.02$ for LRGs and $r_0 = 5.45 \pm 0.1$ $h^{-1}$Mpc, $γ= 1.54 \pm 0.01$ for ELGs. Additionally, from the angular power spectra, we measure the linear biases and model the scale dependent biases in the weakly nonlinear regime. Both sets of clustering measurements show good agreement with survey requirements for LRGs and ELGs, attesting that these samples will enable DESI to achieve precise cosmological constraints. We also present clustering as a function of magnitude, use cross-correlations with external spectroscopy to infer $dN/dz$ and measure clustering as a function of luminosity, and probe higher order clustering statistics through counts-in-cells moments.

preprint2020arXiv

Minimal 4D Composite Higgs Models Under Current LHC Constraints

We present preliminary results of the first convergent global fits of several minimal composite Higgs models. Our fits are performed using the differential evolution optimisation package $\texttt{Diver}$. A variety of physical constraints are taken into account, including a wide range of exclusion bounds on heavy resonance production from Run 2 of the LHC. As a by-product of the fits, we analyse the collider phenomenology of the lightest new up-type and down-type resonances in the viable regions of our models, finding some low-mass resonances that can be probed in future collider searches.

preprint2020arXiv

Model-independent approach for incorporating interference effects in collider searches for new resonances

The presence of large-mass resonances in the data collected at the Large Hadron Collider would provide direct evidence of physics beyond the Standard Model. A key challenge in current resonance searches at the LHC is the modelling of signal--background interference effects, which can severely distort the shape of the reconstructed invariant mass distribution relative to the case where there is no interference. Such effects are strongly dependent on the beyond the Standard Model theory that must be considered as unknown if one aims to minimise any theoretical bias on the search results. In this paper, we describe a procedure which employs a physically-motivated, model-independent template functional form that can be used to model interference effects, both for the characterisation of positive discoveries, and in the presentation of null results. We illustrate the approach with the example of a scalar resonance decaying into a pair of photons.

preprint2020arXiv

Packed Ultra-wideband Mapping Array (PUMA): Astro2020 RFI Response

The Packed Ultra-wideband Mapping Array (PUMA) is a proposed low-resolution transit interferometric radio telescope operating over the frequency range 200 - 1100MHz. Its rich science portfolio will include measuring structure in the universe from redshift z = 0.3 to 6 using 21cm intensity mapping, detecting one million fast radio bursts, and monitoring thousands of pulsars. It will allow PUMA to advance science in three different areas of physics (the physics of dark energy, the physics of cosmic inflation and time-domain astrophysics). This document is a response to a request for information (RFI) by the Panel on Radio, Millimeter, and Submillimeter Observations from the Ground (RMS) of the Decadal Survey on Astronomy and Astrophysics 2020. We present the science case of PUMA, the development path and major risks to the project.

preprint2020arXiv

Reconstruction with velocities

Reconstruction is becoming a crucial procedure of galaxy clustering analysis for future spectroscopic redshift surveys to obtain subpercent level measurement of the baryon acoustic oscillation scale. Most reconstruction algorithms rely on an estimation of the displacement field from the observed galaxy distribution. However, the displacement reconstruction degrades near the survey boundary due to incomplete data and the boundary effects extend to $\sim100\,\mathrm{Mpc}/h$ within the interior of the survey volume. We study the possibility of using radial velocities measured from the cosmic microwave background observation through the kinematic Sunyaev-Zel'dovich effect to improve performance near the boundary. We find that the boundary effect can be reduced to $\sim30-40\,\mathrm{Mpc}/h$ with the velocity information from Simons Observatory. This is especially helpful for dense low redshift surveys where the volume is relatively small and a large fraction of total volume is affected by the boundary.

preprint2020arXiv

Self-consistent approach for measuring the energy spectra and composition of cosmic rays and determining the properties of hadronic interactions at high energy

Air showers, produced by the interaction of energetic cosmic rays with the atmosphere, are an excellent alternative to study particle physics at energies beyond any human-made particle accelerator. For that, it is necessary to identify first the mass composition of the primary cosmic ray (and its energy). None of the existing high energy interaction models have been able to reproduce coherently all air shower observables over the entire energy and zenith angle phase space. This is despite having tried all possible combinations for the cosmic ray mass composition. This proposal outlines a self-consistent strategy to study high energy particle interactions and identify the energy spectra and mass composition of cosmic rays. This strategy involves the participation of different particle accelerators and astrophysics experiments. This is important to cover the entire cosmic ray energy range and a larger phase-space of shower observables to probe the high energy interaction models.

preprint2020arXiv

Simulations and symmetries

We investigate the range of applicability of a model for the real-space power spectrum based on N-body dynamics and a (quadratic) Lagrangian bias expansion. This combination uses the highly accurate particle displacements that can be efficiently achieved by modern N-body methods with a symmetries-based bias expansion which describes the clustering of any tracer on large scales. We show that at low redshifts, and for moderately biased tracers, the substitution of N-body-determined dynamics improves over an equivalent model using perturbation theory by more than a factor of two in scale, while at high redshifts and for highly biased tracers the gains are more modest. This hybrid approach lends itself well to emulation. By removing the need to identify halos and subhalos, and by not requiring any galaxy-formation-related parameters to be included, the emulation task is significantly simplified at the cost of modeling a more limited range in scale.

preprint2020arXiv

TARDIS Paper I: A Constrained Reconstruction Approach to Modeling the z~2.5 Cosmic Web Probed by Lyman-alpha Forest Tomography

Recent Lyman-$α$ forest tomography measurements of the intergalactic medium (IGM) have revealed a wealth of cosmic structures at high redshift ($z\sim 2.5$). In this work, we present the Tomographic Absorption Reconstruction and Density Inference Scheme (TARDIS), a new chrono-cosmographic analysis tool for understanding the formation and evolution of these observed structures. We use maximum likelihood techniques with a fast non-linear gravitational model to reconstruct the initial density field of the observed regions. We find that TARDIS allows accurate reconstruction of smaller scale structures than standard Wiener filtering techniques. Applying this technique to mock Lyman-$α$ forest data sets that simulate ongoing and future surveys such as CLAMATO, Subaru-PFS or the ELTs, we are able to infer the underlying matter density field at observed redshift and classify the cosmic web structures. We find good agreement with the underlying truth both in the characteristic eigenvalues and eigenvectors of the pseudo-deformation tensor, with the eigenvalues inferred from 30m-class telescopes correlated at $r=0.95$ relative to the truth. As an output of this method, we are able to further evolve the inferred structures to late time ($z=0$), and also track the trajectories of coeval $z=2.5$ galaxies to their $z=0$ cosmic web environments.

preprint2020arXiv

The Reconstructed Power Spectrum in the Zeldovich Approximation

Density-reconstruction sharpens the baryon acoustic oscillations signal by undoing some of the smoothing incurred by nonlinear structure formation. In this paper we present an analytical model for reconstruction based on the Zeldovich approximation, which for the first time includes a complete set of counterterms and bias terms up to quadratic order and can fit real and redshift-space data pre- and post-reconstruction data in both Fourier and configuration space over a wide range of scales. We compare our model to n-body data at $z = 0$ from the {\tt DarkSky} simulation \cite{skillman14}, finding sub-percent agreement in both real space and in the redshift-space power spectrum monopole out to $k = 0.4\ h$ Mpc$^{-1}$, and out to $k = 0.2\ h$ Mpc$^{-1}$ in the quadrupole, with comparable agreement in configuration space. We compare our model with several popular existing alternatives, updating existing theoretical results for exponential damping in wiggle/no-wiggle splits of the BAO signal and discuss the usually-ignored effect of higher bias contributions on the reconstructed signal. In the appendices, we re-derive the former within our formalism, present exploratory results on higher-order corrections due to nonlinearities inherent to reconstruction, and present numerical techniques with which to calculate the redshift-space power spectrum of biased tracers within the Zeldovich approximation.

preprint2019arXiv

Axion global fits with Peccei-Quinn symmetry breaking before inflation using GAMBIT

We present global fits of cosmologically stable axion-like particle and QCD axion models in the mass range 0.1 neV to 10 eV. We focus on the case where the Peccei-Quinn symmetry is broken before the end of inflation, such that the initial value of the axion field can be considered to be homogeneous throughout the visible Universe. We include detailed likelihood functions from light-shining-through-wall experiments, haloscopes, helioscopes, the axion relic density, horizontal branch stars, supernova 1987A, white dwarf cooling, and gamma-ray observations. We carry out both frequentist and Bayesian analyses, with and without the inclusion of white dwarf cooling. We explore the degree of fine-tuning present in different models and identify parameter regions where it is possible for QCD axion models to account for both the dark matter in the Universe and the cooling hints, comparing them to specific DFSZ- and KSVZ-type models. We find the most credible parameter regions, allowing us to set (prior-dependent) upper and lower bounds on the axion mass. Our analysis also suggests that QCD axions in this scenario most probably make up a non-negligible but sub-dominant component of the dark matter in the Universe.

preprint2019arXiv

ColliderBit: a GAMBIT module for the calculation of high-energy collider observables and likelihoods

We describe ColliderBit, a new code for the calculation of high energy collider observables in theories of physics beyond the Standard Model (BSM). ColliderBit features a generic interface to BSM models, a unique parallelised Monte Carlo event generation scheme suitable for large-scale supercomputer applications, and a number of LHC analyses, covering a reasonable range of the BSM signatures currently sought by ATLAS and CMS. ColliderBit also calculates likelihoods for Higgs sector observables, and LEP searches for BSM particles. These features are provided by a combination of new code unique to ColliderBit, and interfaces to existing state-of-the-art public codes. ColliderBit is both an important part of the GAMBIT framework for BSM inference, and a standalone tool for efficiently applying collider constraints to theories of new physics.

preprint2019arXiv

Combined collider constraints on neutralinos and charginos

Searches for supersymmetric electroweakinos have entered a crucial phase, as the integrated luminosity of the Large Hadron Collider is now high enough to compensate for their weak production cross-sections. Working in a framework where the neutralinos and charginos are the only light sparticles in the Minimal Supersymmetric Standard Model, we use gambit to perform a detailed likelihood analysis of the electroweakino sector. We focus on the impacts of recent ATLAS and CMS searches with 36 fb$^{-1}$ of 13 TeV proton-proton collision data. We also include constraints from LEP and invisible decays of the $Z$ and Higgs bosons. Under the background-only hypothesis, we show that current LHC searches do not robustly exclude any range of neutralino or chargino masses. However, a pattern of excesses in several LHC analyses points towards a possible signal, with neutralino masses of $(m_{\tildeχ_1^0}, m_{\tildeχ_2^0}, m_{\tildeχ_3^0}, m_{\tildeχ_4^0})$ = (8-155, 103-260, 130-473, 219-502) GeV and chargino masses of $(m_{\tildeχ_1^{\pm}}, m_{\tildeχ_2^{\pm}})$ = (104-259, 224-507) GeV at the 95% confidence level. The lightest neutralino is mostly bino, with a possible modest Higgsino or wino component. We find that this excess has a combined local significance of $3.3σ$, subject to a number of cautions. If one includes LHC searches for charginos and neutralinos conducted with 8 TeV proton-proton collision data, the local significance is lowered to 2.9$σ$. We briefly consider the implications for dark matter, finding that the correct relic density can be obtained through the Higgs-funnel and $Z$-funnel mechanisms, even assuming that all other sparticles are decoupled. All samples, gambit input files and best-fit models from this study are available on Zenodo.

preprint2019arXiv

CosmoDC2: A Synthetic Sky Catalog for Dark Energy Science with LSST

This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `Outer Rim' run. It covers 440 deg^2 of sky area to a redshift of z=3 and is complete to a magnitude depth of 28 in the r-band. Each galaxy is characterized by a multitude of properties including stellar mass, morphology, spectral energy distributions, broadband filter magnitudes, host halo information and weak lensing shear. The size and complexity of cosmoDC2 requires an efficient catalog generation methodology; our approach is based on a new hybrid technique that combines data-driven empirical approaches with semi-analytic galaxy modeling. A wide range of observation-based validation tests has been implemented to ensure that cosmoDC2 enables the science goals of the planned LSST DESC DC2 analyses. This paper also represents the official release of the cosmoDC2 data set, including an efficient reader that facilitates interaction with the data.

preprint2019arXiv

Higgs Quark Flavor Violation: Simplified Models and Status of General Two-Higgs-Doublet Model

We study quark flavor violating interactions mediated by the Higgs boson $h$. We consider observables involving a third generation quark, of both the up and the down quark sectors, like $h\rightarrow bs$ and $t\rightarrow ch$. Using an effective field theory approach we systematically list all the possible tree-level ultraviolet completions, which comprise models with vector-like quarks and/or extra scalars. We provide upper bounds on the flavor violating transitions allowed by current limits stemming from low energy processes, such as meson mixing and $b \rightarrow s γ$. We find that scenarios with vector-like quarks always have very suppressed flavor-violating transitions, while a general two Higgs doublet model may have a sizeable rate. To study the latter case in detail, we perform a full numerical simulation taking into account all relevant theoretical and phenomenological constraints. Our results show that ${\rm BR}(t\rightarrow ch)$ [${\rm BR}(h\rightarrow bs)$] are still allowed at the subpercent [percent] level, which are being [may be] explored at the LHC [future colliders]. Finally, we have found that the mild mass-splitting discrepancy with respect to the SM in the $B_s$ meson system can be accommodated in the Two-Higgs-Doublet Model. If confirmed, it yields the prediction ${\rm BR}(h\rightarrow bs)\simeq 10^{-4}$, if the new contribution to the mass-splitting is dominated by tree-level Higgs boson exchange.

preprint2019arXiv

SpecBit, DecayBit and PrecisionBit: GAMBIT modules for computing mass spectra, particle decay rates and precision observables

We present the GAMBIT modules SpecBit, DecayBit and PrecisionBit. Together they provide a new framework for linking publicly available spectrum generators, decay codes and other precision observable calculations in a physically and statistically consistent manner. This allows users to automatically run various combinations of existing codes as if they are a single package. The modular design allows software packages fulfilling the same role to be exchanged freely at runtime, with the results presented in a common format that can be easily passed to downstream dark matter, collider and flavour codes. These modules constitute an essential part of the broader GAMBIT framework, a major new software package for performing global fits. In this paper we present the observable calculations, data, and likelihood functions implemented in the three modules, as well as the conventions and assumptions used in interfacing them with external codes. We also present 3-BIT-HIT, a command-line utility for computing mass spectra, couplings, decays and precision observables in the MSSM, which shows how the three modules can be easily used independently of GAMBIT.

preprint2019arXiv

Wide angle effects for peculiar velocities

The line-of-sight peculiar velocities of galaxies contribute to their observed redshifts, breaking the translational invariance of galaxy clustering down to a rotational invariance around the observer. This becomes important when the line-of-sight direction varies significantly across a survey, leading to what are known as `wide angle' effects in redshift space distortions. Wide-angle effects will also be present in measurements of the momentum field, i.e. the galaxy density-weighted velocity field, in upcoming peculiar velocity surveys. In this work we study how wide-angle effects modify the predicted correlation function and power spectrum for momentum statistics, both in auto-correlation and in cross-correlation with the density field. Using both linear theory and the Zeldovich approximation, we find that deviations from the plane-parallel limit are large and could become important in data analysis for low redshift surveys. We point out that even multipoles in the cross-correlation between density and momentum are non-zero regardless of the choice of line of sight, and therefore contain new cosmological information that could be exploited. We discuss configuration-space, Fourier-space and spherical analyses, providing exact expressions in each case rather than relying on an expansion in small angles. We hope these expressions will be of use in the analysis of upcoming surveys for redshift-space distortions and peculiar velocities.

preprint2018arXiv

A global fit of the MSSM with GAMBIT

We study the seven-dimensional Minimal Supersymmetric Standard Model (MSSM7) with the new GAMBIT software framework, with all parameters defined at the weak scale. Our analysis significantly extends previous weak-scale, phenomenological MSSM fits, by adding more and newer experimental analyses, improving the accuracy and detail of theoretical predictions, including dominant uncertainties from the Standard Model, the Galactic dark matter halo and the quark content of the nucleon, and employing novel and highly-efficient statistical sampling methods to scan the parameter space. We find regions of the MSSM7 that exhibit co-annihilation of neutralinos with charginos, stops and sbottoms, as well as models that undergo resonant annihilation via both light and heavy Higgs funnels. We find high-likelihood models with light charginos, stops and sbottoms that have the potential to be within the future reach of the LHC. Large parts of our preferred parameter regions will also be accessible to the next generation of direct and indirect dark matter searches, making prospects for discovery in the near future rather good.

preprint2018arXiv

DarkBit: A GAMBIT module for computing dark matter observables and likelihoods

We introduce DarkBit, an advanced software code for computing dark matter constraints on various extensions to the Standard Model of particle physics, comprising both new native code and interfaces to external packages. This release includes a dedicated signal yield calculator for gamma-ray observations, which significantly extends current tools by implementing a cascade decay Monte Carlo, as well as a dedicated likelihood calculator for current and future experiments (gamlike). This provides a general solution for studying complex particle physics models that predict dark matter annihilation to a multitude of final states. We also supply a direct detection package that models a large range of direct detection experiments (DDcalc), and provides the corresponding likelihoods for arbitrary combinations of spin-independent and spin-dependent scattering processes. Finally, we provide custom relic density routines along with interfaces to DarkSUSY, micrOMEGAs, and the neutrino telescope likelihood package nuLike. DarkBit is written in the framework of the Global And Modular Beyond the Standard Model Inference Tool (GAMBIT), providing seamless integration into a comprehensive statistical fitting framework that allows users to explore new models with both particle and astrophysics constraints, and a consistent treatment of systematic uncertainties. In this paper we describe its main functionality, provide a guide to getting started quickly, and show illustrative examples for results obtained with DarkBit (both as a standalone tool and as a GAMBIT module). This includes a quantitative comparison between two of the main dark matter codes (DarkSUSY and micrOMEGAs), and application of DarkBit's advanced direct and indirect detection routines to a simple effective dark matter model.

preprint2018arXiv

Global fits of GUT-scale SUSY models with GAMBIT

We present the most comprehensive global fits to date of three supersymmetric models motivated by grand unification: the Constrained Minimal Supersymmetric Standard Model (CMSSM), and its Non-Universal Higgs Mass generalisations NUHM1 and NUHM2. We include likelihoods from a number of direct and indirect dark matter searches, a large collection of electroweak precision and flavour observables, direct searches for supersymmetry at LEP and Runs I and II of the LHC, and constraints from Higgs observables. Our analysis improves on existing results not only in terms of the number of included observables, but also in the level of detail with which we treat them, our sampling techniques for scanning the parameter space, and our treatment of nuisance parameters. We show that stau co-annihilation is now ruled out in the CMSSM at more than 95\% confidence. Stop co-annihilation turns out to be one of the most promising mechanisms for achieving an appropriate relic density of dark matter in all three models, whilst avoiding all other constraints. We find high-likelihood regions of parameter space featuring light stops and charginos, making them potentially detectable in the near future at the LHC. We also show that tonne-scale direct detection will play a largely complementary role, probing large parts of the remaining viable parameter space, including essentially all models with multi-TeV neutralinos.

preprint2018arXiv

Status of the scalar singlet dark matter model

One of the simplest viable models for dark matter is an additional neutral scalar, stabilised by a $\mathbb{Z}_2$ symmetry. Using the GAMBIT package and combining results from four independent samplers, we present Bayesian and frequentist global fits of this model. We vary the singlet mass and coupling along with 13 nuisance parameters, including nuclear uncertainties relevant for direct detection, the local dark matter density, and selected quark masses and couplings. We include the dark matter relic density measured by Planck, direct searches with LUX, PandaX, SuperCDMS and XENON100, limits on invisible Higgs decays from the Large Hadron Collider, searches for high-energy neutrinos from dark matter annihilation in the Sun with IceCube, and searches for gamma rays from annihilation in dwarf galaxies with the Fermi-LAT. Viable solutions remain at couplings of order unity, for singlet masses between the Higgs mass and about 300 GeV, and at masses above $\sim$1 TeV. Only in the latter case can the scalar singlet constitute all of dark matter. Frequentist analysis shows that the low-mass resonance region, where the singlet is about half the mass of the Higgs, can also account for all of dark matter, and remains viable. However, Bayesian considerations show this region to be rather fine-tuned.

preprint2017arXiv

FlavBit: A GAMBIT module for computing flavour observables and likelihoods

Flavour physics observables are excellent probes of new physics up to very high energy scales. Here we present FlavBit, the dedicated flavour physics module of the global-fitting package GAMBIT. FlavBit includes custom implementations of various likelihood routines for a wide range of flavour observables, including detailed uncertainties and correlations associated with LHCb measurements of rare, leptonic and semileptonic decays of B and D mesons, kaons and pions. It provides a generalised interface to external theory codes such as SuperIso, allowing users to calculate flavour observables in and beyond the Standard Model, and then test them in detail against all relevant experimental data. We describe FlavBit and its constituent physics in some detail, then give examples from supersymmetry and effective field theory illustrating how it can be used both as a standalone library for flavour physics, and within GAMBIT.

preprint2017arXiv

GAMBIT: The Global and Modular Beyond-the-Standard-Model Inference Tool

We describe the open-source global fitting package GAMBIT: the Global And Modular Beyond-the-Standard-Model Inference Tool. GAMBIT combines extensive calculations of observables and likelihoods in particle and astroparticle physics with a hierarchical model database, advanced tools for automatically building analyses of essentially any model, a flexible and powerful system for interfacing to external codes, a suite of different statistical methods and parameter scanning algorithms, and a host of other utilities designed to make scans faster, safer and more easily-extendible than in the past. Here we give a detailed description of the framework, its design and motivation, and the current models and other specific components presently implemented in GAMBIT. Accompanying papers deal with individual modules and present first GAMBIT results. GAMBIT can be downloaded from gambit.hepforge.org.

preprint2016arXiv

Combined analysis of effective Higgs portal dark matter models

We combine and extend the analyses of effective scalar, vector, Majorana and Dirac fermion Higgs portal models of dark matter (DM), in which DM couples to the Standard Model (SM) Higgs boson via an operator of the form $\mathcal{O}_{\textrm{DM}}\, H^\dagger H$. For the fermion models, we take an admixture of scalar $\overlineψ ψ$ and pseudoscalar $\overlineψ iγ_5 ψ$ interaction terms. For each model, we apply constraints on the parameter space based on the Planck measured DM relic density and the LHC limits on the Higgs invisible branching ratio. For the first time, we perform a consistent study of the indirect detection prospects for these models based on the WMAP7/Planck observations of the cosmic microwave background, a combined analysis of 15 dwarf spheroidal galaxies by Fermi-LAT and the upcoming Cherenkov Telescope Array (CTA). We also perform a correct treatment of the momentum-dependent direct search cross section that arises from the pseudoscalar interaction term in the fermionic DM theories. We find, in line with previous studies, that current and future direct search experiments such as LUX and XENON1T can exclude much of the parameter space, and we demonstrate that a joint observation in both indirect and direct searches is possible for high mass weakly interacting massive particles. In the case of a pure pseudoscalar interaction of a fermionic DM candidate, future gamma-ray searches are the only class of experiment capable of probing the high mass range of the theory.

preprint2016arXiv

Combining galaxy and 21cm surveys

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.

preprint2016arXiv

Connecting massive galaxies to dark matter halos in BOSS - I. Is galaxy color a stochastic process in high-mass halos?

We use subhalo abundance matching (SHAM) to model the stellar mass function (SMF) and clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) "CMASS" sample at $z\sim0.5$. We introduce a novel method which accounts for the stellar mass incompleteness of CMASS as a function of redshift, and produce CMASS mock catalogs which include selection effects, reproduce the overall SMF, the projected two-point correlation function $w_{\rm p}$, the CMASS $dn/dz$, and are made publicly available. We study the effects of assembly bias above collapse mass in the context of "age matching" and show that these effects are markedly different compared to the ones explored by Hearin et al. (2013) at lower stellar masses. We construct two models, one in which galaxy color is stochastic ("AbM" model) as well as a model which contains assembly bias effects ("AgM" model). By confronting the redshift dependent clustering of CMASS with the predictions from our model, we argue that that galaxy colors are not a stochastic process in high-mass halos. Our results suggest that the colors of galaxies in high-mass halos are determined by other halo properties besides halo peak velocity and that assembly bias effects play an important role in determining the clustering properties of this sample.

preprint2016arXiv

Cosmological perturbation theory in 1+1 dimensions

Many recent studies have highlighted certain failures of the standard Eulerian-space cosmological perturbation theory (SPT). Its problems include (1) not capturing large-scale bulk flows [leading to an O(1) error in the 1-loop SPT prediction for the baryon acoustic peak in the correlation function], (2) assuming that the Universe behaves as a pressureless, inviscid fluid, and (3) treating fluctuations on scales that are non-perturbative as if they were. Recent studies have highlighted the successes of perturbation theory in Lagrangian space or theories that solve equations for the effective dynamics of smoothed fields. Both approaches mitigate some or all of the aforementioned issues with SPT. We discuss these physical developments by specializing to the simplified 1D case of gravitationally interacting sheets, which allows us to substantially reduces the analytic overhead and still (as we show) maintain many of the same behaviors as in 3D. In 1D, linear-order Lagrangian perturbation theory ("the Zeldovich approximation") is exact up to shell crossing, and we prove that n^{th}-order Eulerian perturbation theory converges to the Zeldovich approximation as n goes to infinity. In no 1D cosmology that we consider (including a CDM-like case and power-law models) do these theories describe accurately the matter power spectrum on any mildly nonlinear scale. We find that theories based on effective equations are much more successful at describing the dynamics. Finally, we discuss many topics that have recently appeared in the perturbation theory literature such as beat coupling, the shift and smearing of the baryon acoustic oscillation feature, and the advantages of Fourier versus configuration space. Our simplified 1D case serves as an intuitive review of these perturbation theory results.

preprint2016arXiv

Gamma rays from muons from WIMPs: Implementation of radiative muon decays for dark matter analyses

Dark matter searches in gamma ray final states often make use of the fact that photons can be produced from final state muons. Modern Monte Carlo generators and DM codes include the effects of final state radiation from muons produced in the dark matter annihilation process itself, but neglect the O(1%) radiative correction that arises from the subsequent muon decay. After implementing this correction we demonstrate the effect that it can have on dark matter phenomenology by considering the case of dark matter annihilation to four muons via scalar mediator production. We first show that the AMS-02 positron excess can no longer easily be made consistent with this final state once the Fermi-LAT dwarf limits are calculated with the inclusion of radiative muon decays, and we next show that the Fermi-LAT galactic centre gamma excess can be improved with this final state after inclusion of the same effect. We provide code and tables for the implementation of this effect in the popular dark matter code micrOMEGAs, providing a solution for any model producing final state muons.

preprint2016arXiv

Improving photometric redshifts with Ly$α$ tomography

Forming a three dimensional view of the Universe is a long-standing goal of astronomical observations, and one that becomes increasingly difficult at high redshift. In this paper we discuss how tomography of the intergalactic medium (IGM) at $z\simeq 2.5$ can be used to estimate the redshifts of massive galaxies in a large volume of the Universe based on spectra of galaxies in their background. Our method is based on the fact that hierarchical structure formation leads to a strong dependence of the halo density on large-scale environment. A map of the latter can thus be used to refine our knowledge of the redshifts of halos and the galaxies and AGN which they host. We show that tomographic maps of the IGM at a resolution of $2.5\,h^{-1}$Mpc can determine the redshifts of more than 90 per cent of massive galaxies with redshift uncertainty $Δz/(1+z)=0.01$. Higher resolution maps allow such redshift estimation for lower mass galaxies and halos.

preprint2016arXiv

Including parameter dependence in the data and covariance for cosmological inference

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.

preprint2016arXiv

On using angular cross-correlations to determine source redshift distributions

We investigate how well the redshift distribution of a population of extragalactic objects can be reconstructed using angular cross-correlations with a sample whose redshifts are known. We derive the minimum variance quadratic estimator, which has simple analytic representations in very applicable limits and is significantly more sensitive than earlier proposed estimation procedures. This estimator is straightforward to apply to observations, it robustly finds the likelihood maximum, and it conveniently selects angular scales at which fluctuations are well approximated as independent between redshift bins and at which linear theory applies. We find that the linear bias times number of objects in a redshift bin generally can be constrained with cross-correlations to fractional error (10^2 n/N)^1/2, where N is the total number of spectra per dz and n is the number of redshift bins spanned by the bulk of the unknown population. The error is often independent of the sky area and sampling fraction. Furthermore, we find that sub-percent measurements of the angular source density per unit redshift, dN/dz, are in principle possible, although cosmic magnification needs to be accounted for at fractional errors of <~ 10 per cent. We discuss how the sensitivity to dN/dz changes as a function of photometric and spectroscopic depth and how to optimize the survey strategy to constrain dN/dz. We also quantify how well cross-correlations of photometric redshift bins can be used to self-calibrate a photometric redshift sample. Simple formulae that can be quickly applied to gauge the utility of cross correlating different samples are given.

preprint2016arXiv

Redshift Weights for Baryon Acoustic Oscillations : Application to Mock Galaxy Catalogs

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.

preprint2016arXiv

Revealing the z~2.5 Cosmic Web With 3D Lyman-Alpha Forest Tomography: A Deformation Tensor Approach

Studies of cosmological objects should take into account their positions within the cosmic web of large-scale structure. Unfortunately, the cosmic web has only been extensively mapped at low-redshifts ($z<1$), using galaxy redshifts as tracers of the underlying density field. At $z>1$, the required galaxy densities are inaccessible for the foreseeable future, but 3D reconstructions of Lyman-$α$ forest absorption in closely-separated background QSOs and star-forming galaxies already offer a detailed window into $z\sim2-3$ large-scale structure. We quantify the utility of such maps for studying the cosmic web by using realistic $z=2.5$ Ly$α$ forest simulations matched to observational properties of upcoming surveys. A deformation tensor-based analysis is used to classify voids, sheets, filaments and nodes in the flux, which is compared to those determined from the underlying dark matter field. We find an extremely good correspondence, with $70\%$ of the volume in the flux maps correctly classified relative to the dark matter web, and $99\%$ classified to within 1 eigenvalue. This compares favorably to the performance of galaxy-based classifiers with even the highest galaxy densities at low-redshift. We find that narrow survey geometries can degrade the cosmic web recovery unless the survey is $\gtrsim 60\,h^{-1}\,\mathrm{Mpc}$ or $\gtrsim 1\,\mathrm{deg}$ on the sky. We also examine halo abundances as a function of the cosmic web, and find a clear dependence as a function of flux overdensity, but little explicit dependence on the cosmic web. These methods will provide a new window on cosmological environments of galaxies at this very special time in galaxy formation, "high noon", and on overall properties of cosmological structures at this epoch.

preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Observational systematics and baryon acoustic oscillations in the correlation function

We present baryon acoustic oscillation (BAO) scale measurements determined from the clustering of 1.2 million massive galaxies with redshifts 0.2 < z < 0.75 distributed over 9300 square degrees, as quantified by their redshift-space correlation function. In order to facilitate these measurements, we define, describe, and motivate the selection function for galaxies in the final data release (DR12) of the SDSS III Baryon Oscillation Spectroscopic Survey (BOSS). This includes the observational footprint, masks for image quality and Galactic extinction, and weights to account for density relationships intrinsic to the imaging and spectroscopic portions of the survey. We simulate the observed systematic trends in mock galaxy samples and demonstrate that they impart no bias on baryon acoustic oscillation (BAO) scale measurements and have a minor impact on the recovered statistical uncertainty. We measure transverse and radial BAO distance measurements in 0.2 < z < 0.5, 0.5 < z < 0.75, and (overlapping) 0.4 < z < 0.6 redshift bins. In each redshift bin, we obtain a precision that is 2.7 per cent or better on the radial distance and 1.6 per cent or better on the transverse distance. The combination of the redshift bins represents 1.8 per cent precision on the radial distance and 1.1 per cent precision on the transverse distance. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

preprint2016arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the correlation function of LOWZ and CMASS galaxies in Data Release 12

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 $Λ$CDM model with cosmological parameters released by the recent Planck 2015 results.

preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

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-$α$ 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.

preprint2016arXiv

The DESI Experiment Part II: Instrument Design

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= λ/Δλ$ 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.

preprint2016arXiv

The Gaussian streaming model and Lagrangian effective field theory

We update the ingredients of the Gaussian streaming model (GSM) for the redshift-space clustering of biased tracers using the techniques of Lagrangian perturbation theory, effective field theory (EFT) and a generalized Lagrangian bias expansion. After relating the GSM to the cumulant expansion, we present new results for the real-space correlation function, mean pairwise velocity and pairwise velocity dispersion including counter terms from EFT and bias terms through third order in the linear density, its leading derivatives and its shear up to second order. We discuss the connection to the Gaussian peaks formalism. We compare the ingredients of the GSM to a suite of large N-body simulations, and show the performance of the theory on the low order multipoles of the redshift-space correlation function and power spectrum. We highlight the importance of a general biasing scheme, which we find to be as important as higher-order corrections due to non-linear evolution for the halos we consider on the scales of interest to us.

preprint2016arXiv

The SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Overview and Early Data

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.

preprint2016arXiv

The XMM Cluster Survey: The Halo Occupation Number of BOSS galaxies in X-ray clusters

We present a direct measurement of the mean halo occupation distribution (HOD) of galaxies taken from the eleventh data release (DR11) of the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey (BOSS). The HOD of BOSS low-redshift (LOWZ: $0.2 < z < 0.4$) and Constant-Mass (CMASS: $0.43 <z <0.7$) galaxies is inferred via their association with the dark-matter halos of 174 X-ray-selected galaxy clusters drawn from the XMM Cluster Survey (XCS). Halo masses are determined for each galaxy cluster based on X-ray temperature measurements, and range between ${\rm log_{10}} (M_{180}/M_{\odot}) = 13-15$. Our directly measured HODs are consistent with the HOD-model fits inferred via the galaxy-clustering analyses of Parejko et al. for the BOSS LOWZ sample and White et al. for the BOSS CMASS sample. Under the simplifying assumption that the other parameters that describe the HOD hold the values measured by these authors, we have determined a best-fit alpha-index of 0.91$\pm$0.08 and $1.27^{+0.03}_{-0.04}$ for the CMASS and LOWZ HOD, respectively. These alpha-index values are consistent with those measured by White et al. and Parejko et al. In summary, our study provides independent support for the HOD models assumed during the development of the BOSS mock-galaxy catalogues that have subsequently been used to derive BOSS cosmological constraints.

preprint2015arXiv

A Lagrangian effective field theory

We have continued the development of Lagrangian, cosmological perturbation theory for the low-order correlators of the matter density field. We provide a new route to understanding how the effective field theory (EFT) of large-scale structure can be formulated in the Lagrandian framework and a new resummation scheme, comparing our results to earlier work and to a series of high-resolution N-body simulations in both Fourier and configuration space. The `new' terms arising from EFT serve to tame the dependence of perturbation theory on small-scale physics and improve agreement with simulations (though with an additional free parameter). We find that all of our models fare well on scales larger than about two to three times the non-linear scale, but fail as the non-linear scale is approached. This is slightly less reach than has been seen previously. At low redshift the Lagrangian model fares as well as EFT in its Eulerian formulation, but at higher $z$ the Eulerian EFT fits the data to smaller scales than resummed, Lagrangian EFT. All the perturbative models fare better than linear theory.

preprint2015arXiv

Clustering of intermediate redshift quasars using the final SDSS III-BOSS sample

We measure the two-point clustering of spectroscopically confirmed quasars from the final sample of the Baryon Oscillation Spectroscopic Survey (BOSS) on comoving scales of 4 < s < 22 Mpc/h. The sample covers 6950 deg^2 (~ 19 (Gpc/h)^3) and, over the redshift range 2.2 < z < 2.8, contains 55,826 homogeneously selected quasars, which is twice as many as in any similar work. We deduce b_Q = 3.54 +/- 0.10 ; the most precise measurement of quasar bias to date at these redshifts. This corresponds to a host halo mass of ~ 2 x 10^12 ~ M_sun/h with an implied quasar duty cycle of ~1 percent. The real-space projected correlation function is well-fit by a power law of index -2 and correlation length r0 = (8.12 +/- 0.22), Mpc/h over scales of 4 < rp < 25 ~ Mpc/h. To better study the evolution of quasar clustering at moderate redshift, we extend the redshift range of our study to z ~ 3.4 and measure the bias and correlation length of three subsamples over 2.2 < z < 3.4. We find no significant evolution of r0 or bias over this range, implying that the host halo mass of quasars decreases somewhat with increasing redshift. We find quasar clustering remains similar over a decade in luminosity, contradicting a scenario in which quasar luminosity is monotonically related to halo mass at z ~ 2.5. Our results are broadly consistent with previous BOSS measurements, but they yield more precise constraints based upon a larger and more uniform data set.

preprint2015arXiv

Collider constraints on tuning in composite Higgs models

Two potential sources of tuning exist in composite Higgs models: one comes from keeping the Higgs VEV below the compositeness scale and one comes from keeping the Higgs light after constraints on the top partner masses are applied. We construct a measure that determines whether these tunings are independent or not and combines them appropriately. We perform a comprehensive scan of the parameter space for three explicit models and report the minimum tuning values compatible with existing collider constraints. Tuning values are given as functions of resonance masses and deviations to the Higgs couplings so the effect of future constraints can easily be quantified. The current minimum tuning in the minimal model is 2.5-5% and will be decreased to around 0.8-3.3% if no top partners are observed over the lifetime of the LHC.

preprint2015arXiv

Dependence of the Cosmic Microwave Background Lensing Power Spectrum on the Matter Density

The anisotropies in the cosmic microwave background (CMB) provide our best laboratory for testing models of the formation and evolution of large-scale structure. The rich features in the cosmic microwave background anisotropy spectrum, in combination with highly precise observations and theoretical predictions, also allow us to simultaneously constrain a number of cosmological parameters. As observations have progressed, measurements at smaller angular scales have provided increasing leverage. These smaller angular scales provide sensitive measures of the matter density through the effect of gravitational lensing. In this work we provide an analytic explanation of the manner in which the lensing of CMB anisotropies depends on the matter density, finding that the dominant effect comes from the shape of the matter power spectrum set by the decay of small-scale potentials between horizon crossing and matter-radiation equality.

preprint2015arXiv

Estimating sparse precision matrices

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.

preprint2015arXiv

Interpreting the Fermi-LAT gamma ray excess in the simplified framework

We test the plausibility that a Majorana fermion dark matter candidate with a scalar mediator explains the gamma ray excess from the Galactic center. Assuming that the mediator couples to all third generation fermions we calculate observables for dark matter abundance and scattering on nuclei, gamma, positron, and anti-proton cosmic ray fluxes, radio emission from dark matter annihilation, and the effect of dark matter annihilations on the CMB. After discarding the controversial radio observation the rest of the data prefers a dark matter (mediator) mass in the 10-100 (3-1000) GeV region and weakly correlated couplings to bottom quarks and tau leptons with values of 10^{-3}-1 at the 68\% credibility level.

preprint2015arXiv

LUX likelihood and limits on spin-independent and spin-dependent WIMP couplings with LUXCalc

We present LUXCalc, a new utility for calculating likelihoods and deriving WIMP-nucleon coupling limits from the recent results of the LUX direct search dark matter experiment. After a brief review of WIMP-nucleon scattering, we derive LUX limits on the spin-dependent WIMP-nucleon couplings over a broad range of WIMP masses, under standard assumptions on the relevant astrophysical parameters. We find that, under these and other common assumptions, LUX excludes the entire spin-dependent parameter space consistent with a dark matter interpretation of DAMA's anomalous signal, the first time a single experiment has been able to do so. We also revisit the case of spin-independent couplings, and demonstrate good agreement between our results and the published LUX results. Finally, we derive constraints on the parameters of an effective dark matter theory in which a spin-1 mediator interacts with a fermionic WIMP and Standard Model fermions via axial-vector couplings. A detailed appendix describes the use of LUXCalc with standard codes to place constraints on generic dark matter theories.

preprint2015arXiv

Reconstruction within the Zeldovich approximation

The Zeldovich approximation, 1st order Lagrangian perturbation theory, provides a good description of the clustering of matter and galaxies on large scales. The acoustic feature in the large-scale correlation function of galaxies imprinted by sound waves in the early Universe has been successfully used as a `standard ruler' to constrain the expansion history of the Universe. The standard ruler can be improved if a process known as density field reconstruction is employed. In this paper we develop the Zeldovich formalism to compute the correlation function of biased tracers in both real- and redshift-space using the simplest reconstruction algorithm with a Gaussian kernel and compare to N-body simulations. The model qualitatively describes the effects of reconstruction on the simulations, though its quantitative success depends upon how redshift-space distortions are handled in the reconstruction algorithm.

preprint2015arXiv

SDSS-III Baryon Oscillation Spectroscopic Survey Data Release 12: galaxy target selection and large scale structure catalogues

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.

preprint2015arXiv

Shadow of a Colossus: A z=2.45 Galaxy Protocluster Detected in 3D Ly-a Forest Tomographic Mapping of the COSMOS Field

Using moderate-resolution optical spectra from 58 background Lyman-break galaxies and quasars at $z\sim 2.3-3$ within a $11.5'\times13.5'$ area of the COSMOS field ($\sim 1200\,\mathrm{deg}^2$ projected area density or $\sim 2.4\,h^{-1}\,\mathrm{Mpc}$ mean transverse separation), we reconstruct a 3D tomographic map of the foreground Ly$α$ forest absorption at $2.2<z<2.5$ with an effective smoothing scale of $σ_{3d}\approx3.5\,h^{-1}\,\mathrm{Mpc}$ comoving. Comparing with 61 coeval galaxies with spectroscopic redshifts in the same volume, we find that the galaxy positions are clearly biased towards regions with enhanced IGM absorption in the tomographic map. We find an extended IGM overdensity with deep absorption troughs at $z=2.45$ associated with a recently-discovered galaxy protocluster at the same redshift. Based on simulations matched to our data, we estimate the enclosed dark matter mass within this IGM overdensity to be $M_{\rm dm} (z=2.45) = (9\pm4)\times 10^{13}\,h^{-1}\,\mathrm{M_\odot}$, and argue based on this mass and absorption strength that it will form at least one $z\sim0$ galaxy cluster with $M(z=0) = (3\pm 2) \times 10^{14}\,h^{-1}\mathrm{M_\odot}$, although its elongated nature suggests that it will likely collapse into two separate clusters. We also point out a compact overdensity of six MOSDEF galaxies at $z=2.30$ within a $r\sim 1\,h^{-1}\,\mathrm{Mpc}$ radius and $Δz\sim 0.006$, which does not appear to have a large associated IGM overdensity. These results demonstrate the potential of Ly$α$ forest tomography on larger volumes to study galaxy properties as a function of environment, as well as revealing the large-scale IGM overdensities associated with protoclusters and other features of large-scale structure.

preprint2015arXiv

The Eleventh and Twelfth Data Releases of the Sloan Digital Sky Survey: Final Data from SDSS-III

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.

preprint2014arXiv

A 2.5% measurement of the growth rate from small-scale redshift space clustering of SDSS-III CMASS galaxies

We perform the first fit to the anisotropic clustering of SDSS-III CMASS DR10 galaxies on scales of ~ 0.8 - 32 Mpc/h. A standard halo occupation distribution model evaluated near the best fit Planck LCDM cosmology provides a good fit to the observed anisotropic clustering, and implies a normalization for the peculiar velocity field of M ~ 2 x 10^13 Msun/h halos of f*sigma8(z=0.57) = 0.450 +/- 0.011. Since this constraint includes both quasi-linear and non-linear scales, it should severely constrain modified gravity models that enhance pairwise infall velocities on these scales. Though model dependent, our measurement represents a factor of 2.5 improvement in precision over the analysis of DR11 on large scales, f*sigma8(z=0.57) = 0.447 +/- 0.028, and is the tightest single constraint on the growth rate of cosmic structure to date. Our measurement is consistent with the Planck LCDM prediction of 0.480 +/- 0.010 at the ~1.9 sigma level. Assuming a halo mass function evaluated at the best fit Planck cosmology, we also find that 10% of CMASS galaxies are satellites in halos of mass M ~ 6 x 10^13 Msun/h. While none of our tests and model generalizations indicate systematic errors due to an insufficiently detailed model of the galaxy-halo connection, the precision of these first results warrant further investigation into the modeling uncertainties and degeneracies with cosmological parameters.

preprint2014arXiv

Comparing Simple Quasar Demographics Models

This paper explores several simple model variations for the connections among quasars, galaxies, and dark matter halos for redshifts 1 < z < 6. A key component of these models is that we enforce a self-consistent black hole (BH) history by tracking both BH mass and BH growth rate at all redshifts. We connect objects across redshift with a simple constant-number-density procedure, and choose a fiducial model with a relationship between BH and galaxy growth rates that is linear and evolves in a simple way with redshift. Within this fiducial model, we find the quasar luminosity function (QLF) by calculating an "intrinsic" luminosity based on either the BH mass or BH growth rate, and then choosing a model of quasar variability with either a lognormal or truncated power-law distribution of instantaneous luminosities. This gives four model variations, which we fit to the observed QLF at each redshift. With the best-fit models in hand, we undertake a detailed comparison of the four fiducial models, and explore changes to our fiducial model of the BH-galaxy relationship. Each model variation can successfully fit the observed QLF, the shape of which is generally set by the "intrinsic" luminosity at the faint end and by the scatter due to variability at the bright end. We focus on accounting for the reasons that physically different models can make such similar predictions, and on identifying what observational data or physical arguments are most essential in breaking the degeneracies among models.

preprint2014arXiv

Cross-correlation of cosmic far-infrared background anisotropies with large scale structures

We measure the cross-power spectra between luminous red galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS)-III Data Release Eight (DR8) and cosmic infrared background (CIB) anisotropies from Planck and data from the Improved Reprocessing (IRIS) of the Infrared Astronomical Satellite (IRAS) at 353, 545, 857, and 3000 GHz, corresponding to 850, 550, 350 and 100 micron, respectively, in the multipole range 100<l<1000. Using approximately 6.5 10^5 photometrically determined LRGs in 7760 deg^2 of the northern hemisphere in the redshift range 0.45 < z < 0.65, we model the far-infrared background (FIRB) anisotropies with an extended version of the halo model. With these methods, we confirm the basic picture obtained from recent analyses of FIRB anisotropies with Herschel and Planck, that the most efficient halo mass at hosting star forming galaxies is log(M_ eff/M_\odot)=12.84+/-0.15. We estimate the percentage of FIRB anisotropies correlated with LRGs as approximately 11.8 %, 3.9 %, 1.8 %, and 1.0 % of the total at 3000, 857, 545, and 353 GHz, respectively. At redshift z~0.55, the bias of FIRB galaxies with respect to the dark matter density field has the value b_{FIRB}~1.45, and the mean dust temperature of FIRB galaxies is T_d=26 K. Finally, we discuss the impact of present and upcoming cross-correlations with far-infrared background anisotropies on the determination of the global star formation history and the link between galaxies and dark matter.

preprint2014arXiv

Lyman-alpha Forest Tomography from Background Galaxies: The First Megaparsec-Resolution Large-Scale Structure Map at z>2

We present the first observations of foreground Lyman-$α$ forest absorption from high-redshift galaxies, targeting 24 star-forming galaxies (SFGs) with $z\sim 2.3-2.8$ within a $5' \times 15'$ region of the COSMOS field. The transverse sightline separation is $\sim 2\,h^{-1}\mathrm{Mpc}$ comoving, allowing us to create a tomographic reconstruction of the 3D Ly$α$ forest absorption field over the redshift range $2.20\leq z\leq 2.45$. The resulting map covers $6\,h^{-1}\mathrm{Mpc} \times 14\,h^{-1}\mathrm{Mpc}$ in the transverse plane and $230\,h^{-1}\mathrm{Mpc}$ along the line-of-sight with a spatial resolution of $\approx 3.5\,h^{-1}\mathrm{Mpc}$, and is the first high-fidelity map of large-scale structure on $\sim\mathrm{Mpc}$ scales at $z>2$. Our map reveals significant structures with $\gtrsim 10\,h^{-1}\mathrm{Mpc}$ extent, including several spanning the entire transverse breadth, providing qualitative evidence for the filamentary structures predicted to exist in the high-redshift cosmic web. Simulated reconstructions with the same sightline sampling, spectral resolution, and signal-to-noise ratio recover the salient structures present in the underlying 3D absorption fields. Using data from other surveys, we identified 18 galaxies with known redshifts coeval with our map volume enabling a direct comparison to our tomographic map. This shows that galaxies preferentially occupy high-density regions, in qualitative agreement with the same comparison applied to simulations. Our results establishes the feasibility of the CLAMATO survey, which aims to obtain Ly$α$ forest spectra for $\sim 1000$ SFGs over $\sim 1 \,\mathrm{deg}^2$ of the COSMOS field, in order to map out IGM large-scale structure at $\langle z \rangle \sim 2.3$ over a large volume $(100\,h^{-1}\mathrm{Mpc})^3$.

preprint2014arXiv

Observational Requirements for Lyman-alpha Forest Tomographic Mapping of Large-Scale Structure at z ~ 2

The z > 2 Lyman-alpha (Lya) forest traces the underlying dark-matter distribution on large scales and, given sufficient sightlines, can be used to create 3D maps of large-scale structure. We examine the observational requirements to construct such maps and estimate the signal-to-noise as a function of exposure time and sightline density. Sightline densities at z = 2.25 are n_los = [360, 1200,3300] deg^{-2} at limiting magnitudes of g =[24.0, 24.5,25.0], resulting in transverse sightline separations of d_perp = [3.6, 1.9, 1.2] h^{-1} Mpc, which roughly sets the reconstruction scale. We simulate these reconstructions using mock spectra with realistic noise properties, and find that spectra with S/N = 4 per angstrom can be used to generate maps that clearly trace the underlying dark-matter at overdensities of rho/<rho> ~ 1. For the VLT/VIMOS spectrograph, exposure times t_exp = [4, 6, 10] hrs are sufficient for maps with spatial resolution epsilon_3d = [5.0, 3.2, 2.3] h^{-1} Mpc. Assuming ~ 250 h^{-1} Mpc is probed along the line-of-sight, 1 deg^2 of survey area would cover a comoving volume of ~ 10^6 h^{-3} Mpc^3 at <z>=2.3, enabling efficient mapping of large volumes with 8-10m telescopes. These maps could be used to study galaxy environments, detect proto-clusters, and study the topology of large-scale structure at high-z.

preprint2014arXiv

On Estimating Lyman-alpha Forest Correlations between Multiple Sightlines

The next frontier of Lyman-alpha forest studies is the reconstruction of 3D correlations from a dense sample of background sources. The measurement of 3D correlations has the potential to improve constraints on fundamental cosmological parameters, ionizing background models, and the reionization history. This study addresses the sensitivity of spectroscopic surveys to 3D correlations in the Lyman-alpha forest. We show that the sensitivity of a survey to this signal can be quantified by just a single number, a noise-weighted number density of sources on the sky. We investigate how the sensitivity of a spectroscopic quasar (or galaxy) survey scales as a function of its depth, area, and redshift. We propose a simple method for weighting sightlines with varying S/N levels to estimate the correlation function, and we show that this estimator generally performs nearly as well as the minimum variance quadratic estimator. In addition, we show that the sensitivity of a quasar survey to the flux correlation function is generally maximized if it observes each field just long enough to achieve S/N ~ 2 in a 1 A pixel on an L_* quasar while acquiring spectra for all quasars with L > L_*: Little is gained by integrating longer on the same targets or by including fainter quasars. We quantify how these considerations relate to constraints on the angular diameter distance, the curvature of space-time, and the reionization history.

preprint2014arXiv

Optimal Redshift Weighting For Baryon Acoustic Oscillations

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.

preprint2014arXiv

Protocluster Discovery in Tomographic Ly$α$ Forest Flux Maps

We present a new method of finding protoclusters using tomographic maps of Ly$α$ Forest flux. We review our method of creating tomographic flux maps and discuss our new high performance implementation, which makes large reconstructions computationally feasible. Using a large N-body simulation, we illustrate how protoclusters create large-scale flux decrements, roughly 10 $h^{-1}$Mpc across, and how we can use this signal to find them in flux maps. We test the performance of our protocluster finding method by running it on the ideal, noiseless map and tomographic reconstructions from mock surveys, and comparing to the halo catalog. Using the noiseless map, we find protocluster candidates with about 90% purity, and recover about 75% of the protoclusters that form massive clusters ($> 3 \times 10^{14} \, h^{-1} M_{\odot}$). We construct mock surveys similar to the ongoing COSMOS Lyman-Alpha Mapping And Tomography Observations (CLAMATO) survey. While the existing data has an average sightline separation of 2.3 $h^{-1}$Mpc, we test separations of 2 - 6 $h^{-1}$Mpc to see what can be tolerated for our application. Using reconstructed maps from small separation mock surveys, the protocluster candidate purity and completeness are very close what was found in the noiseless case. As the sightline separation increases, the purity and completeness decrease, although they remain much higher than we initially expected. We extended our test cases to mock surveys with an average separation of 15 $h^{-1}$Mpc, meant to reproduce high source density areas of the BOSS survey. We find that even with such a large sightline separation, the method can still be used to find some of the largest protoclusters.

preprint2014arXiv

Tests of redshift-space distortions models in configuration space for the analysis of the BOSS final data release

Observations of redshift-space distortions in spectroscopic galaxy surveys offer an attractive method for observing the build-up of cosmological structure, which depends both on the expansion rate of the Universe and our theory of gravity. In preparation for analysis of redshift-space distortions from the Baryon Oscillation Spectroscopic Survey (BOSS) final data release we compare a number of analytic and phenomenological `streaming' models, specified in configuration space, to mock catalogs derived in different ways from several N-body simulations. The galaxies in each mock catalog have properties similar to those of the higher redshift galaxies measured by BOSS but differ in the details of how small-scale velocities and halo occupancy are determined. We find that all of the analytic models fit the simulations over a limited range of scales while failing at small scales. We discuss which models are most robust and on which scales they return reliable estimates of the rate of growth of structure: we find that models based on some form of resummation can fit our N-body data for BOSS-like galaxies above $30\,h^{-1}$Mpc well enough to return unbiased parameter estimates.

preprint2014arXiv

The Clustering of Galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS): measuring growth rate and geometry with anisotropic clustering

We use the observed anisotropic clustering of galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 11 CMASS sample to measure the linear growth rate of structure, the Hubble expansion rate and the comoving distance scale. Our sample covers 8498 ${\rm deg}^2$ and encloses an effective volume of 6.0 ${\rm Gpc}^3$ at an effective redshift of $\bar{z} = 0.57$. We find $fσ_8 = 0.441 \pm 0.044$, $H = 93.1 \pm 3.0\ {\mathrm{km}\ \mathrm{s}^{-1} \mathrm{Mpc}^{-1}}$ and $D_{\rm A} = 1380 \pm 23\ {\rm Mpc}$ when fitting the growth and expansion rate simultaneously. When we fix the background expansion to the one predicted by spatially-flat $Λ$CDM model in agreement with recent Planck results, we find $fσ_8 = 0.447 \pm 0.028$ (6 per cent accuracy). While our measurements are generally consistent with the predictions of $Λ$CDM and General Relativity, they mildly favor models in which the strength of gravitational interactions is weaker than what is predicted by General Relativity. Combining our measurements with recent cosmic microwave background data results in tight constraints on basic cosmological parameters and deviations from the standard cosmological model. Separately varying these parameters, we find $w = -0.983 \pm 0.075$ (8 per cent accuracy) and $γ= 0.69 \pm 0.11$ (16 per cent accuracy) for the effective equation of state of dark energy and the growth rate index, respectively. Both constraints are in good agreement with the standard model values of $w=-1$ and $γ= 0.554$.

preprint2014arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the Data Release 10 and 11 galaxy samples

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 $Λ$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\,σ$ 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.

preprint2014arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: galaxy clustering measurements in the low redshift sample of Data Release 11

We present the distance measurement to z = 0.32 using the 11th data release of the Sloan Digital Sky Survey-III Baryon Acoustic Oscillation Survey (BOSS). We use 313,780 galaxies of the low-redshift (LOWZ) sample over 7,341 square-degrees to compute $D_V = (1264 \pm 25)(r_d/r_{d,fid})$ - a sub 2% measurement - using the baryon acoustic feature measured in the galaxy two-point correlation function and power-spectrum. We compare our results to those obtained in DR10. We study observational systematics in the LOWZ sample and quantify potential effects due to photometric offsets between the northern and southern Galactic caps. We find the sample to be robust to all systematic effects found to impact on the targeting of higher-redshift BOSS galaxies, and that the observed north-south tensions can be explained by either limitations in photometric calibration or by sample variance, and have no impact on our final result. Our measurement, combined with the baryonic acoustic scale at z = 0.57, is used in Anderson et al. (2013a) to constrain cosmological parameters.

preprint2014arXiv

The Tenth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Apache Point Observatory Galactic Evolution Experiment

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.

preprint2014arXiv

The Zeldovich approximation

This year marks the 100th anniversary of the birth of Yakov Zel'dovich. Amongst his many legacies is the Zel'dovich approximation for the growth of large-scale structure, which remains one of the most successful and insightful analytic models of structure formation. We use the Zel'dovich approximation to compute the two-point function of the matter and biased tracers, and compare to the results of N-body simulations and other Lagrangian perturbation theories. We show that Lagrangian perturbation theories converge well and that the Zel'dovich approximation provides a good fit to the N-body results except for the quadrupole moment of the halo correlation function. We extend the calculation of halo bias to 3rd order and also consider non-local biasing schemes, none of which remove the discrepancy. We argue that a part of the discrepancy owes to an incorrect prediction of inter-halo velocity correlations. We use the Zel'dovich approximation to compute the ingredients of the Gaussian streaming model and show that this hybrid method provides a good fit to clustering of halos in redshift space down to scales of tens of Mpc.

preprint2014arXiv

Using large galaxy surveys to distinguish z~0.5 quiescent galaxy models

One of the most striking properties of galaxies is the bimodality in their star-formation rates. A major puzzle is why any given galaxy is star-forming or quiescent, and a wide range of physical mechanisms have been proposed as solutions. We consider how observations, such as might be available in upcoming large galaxy surveys, might distinguish different galaxy quenching scenarios. To do this, we combine an N-body simulation and multiple prescriptions from the literature to create several quiescent galaxy mock catalogues. Each prescription uses a different set of galaxy properties (such as history, environment, centrality) to assign individual simulation galaxies as quiescent. We find how and how much the resulting quiescent galaxy distributions differ from each other, both intrinsically and observationally. In addition to tracing observational consequences of different quenching mechanisms, our results indicate which sorts of quenching models might be most readily disentangled by upcoming observations and which combinations of observational quantities might provide the most discriminatory power. Our observational measures are auto, cross, and marked correlation functions, pro- jected density distributions, and group multiplicity functions, which rely upon galaxy positions, stellar masses and of course quiescence. Although degeneracies between models are present for individual observations, using multiple observations in concert allows us to distinguish between all ten models we consider. In addition to identifying intrinsic and observational consequences of quiescence prescriptions and testing these quiescence models against each other and observations, these methods can also be used to validate colors (or other history and environment dependent properties) in simulated mock catalogues.

preprint2014arXiv

Velocity Bias from the Small Scale Clustering of SDSS-III BOSS Galaxies

We present the measurements and modelling of the projected and redshift-space clustering of CMASS galaxies in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Data Release 11. For a volume-limited luminous red galaxy sample in the redshift range of $0.48<z<0.55$, we perform halo occupation distribution modelling of the small- and intermediate-scale ($0.1$--$60h^{-1}{\rm {Mpc}}$) projected and redshift-space two-point correlation functions, with an accurate model built on high resolution $N$-body simulations. To interpret the measured redshift-space distortions, the distribution of galaxy velocities must differ from that of the dark matter inside haloes of $\sim 10^{13}$--$10^{14}h^{-1}{\rm M_{\odot}}$, i.e. the data require the existence of galaxy velocity bias. Most notably, central galaxies on average are not at rest with respect to the core of their host haloes (defined by the inner 25% of particles around the halo potential minimum), but rather move around it with a 1D velocity dispersion of $0.22^{+0.03}_{-0.04}$ times that of the dark matter, implying a spatial offset from the centre at the level of $\lesssim$1% of the halo virial radius. The luminous satellite galaxies move more slowly than the dark matter, with velocities $0.86^{+0.08}_{-0.03}$ times those of the dark matter, which suggests that the velocity and spatial distributions of these satellites cannot both be unbiased. The constraints mainly arise from the Fingers-of-God effect at nonlinear scales and the smoothing to the Kaiser effect in the translinear regime; the robustness of the results is demonstrated by a variety of tests. In addition, no clear evidence is found for a strong luminosity dependence of the velocity bias. We discuss the implications of the existence of galaxy velocity bias for investigations of galaxy formation and cosmology.

preprint2013arXiv

A New Method to Directly Measure the Jeans Scale of the Intergalactic Medium Using Close Quasar Pairs

Although the baryons in the intergalactic medium (IGM) trace dark matter on Mpc scales, small-scale (~100 kpc) fluctuations are suppressed by pressure support, analogous to the classical Jeans argument. This Jeans filtering scale has fundamental cosmological implications: it provides a thermal record of heat injected by UV photons during reionization events, determines the clumpiness of the IGM, and sets the minimum mass scale for gravitational collapse, a key quantity in galaxy formation. Unfortunately, it is extremely challenging to measure via the standard analysis of purely longitudinal Lyman-alpha forest spectra, because the thermal Doppler broadening of absorption lines is highly degenerate with Jeans pressure smoothing. In this work we show that the Jeans scale can be directly measured by characterizing the coherence of correlated Lyman-alpha absorption in quasar pairs with separations small enough to resolve it. We present a novel technique for this purpose, based on the probability distribution function (PDF) of phase angle differences of homologous longitudinal Fourier modes in close quasar pair spectra. A Bayesian formalism is introduced based on the phase angle PDF, and MCMC techniques are used to characterize the precision of a future Jeans scale measurement, and explore degeneracies with other thermal parameters governing the IGM. A semi-analytical model of the IGM is used to generate a grid of 500 thermal models from a dark matter simulation. Our full parameter study indicates that a realistic sample of only 20 close quasar pair spectra can pinpoint the Jeans scale to ~ 5% precision, independent of the parameters governing the temperature-density relation of the IGM. We show that this new method is insensitive to a battery of systematics such as continuum fitting errors, imprecise knowledge of the noise and spectral resolution, and metal-line absorption.

preprint2013arXiv

An analytic model for redshift-space distortions

Understanding the formation and evolution of large-scale structure is a central problem in cosmology and enables precise tests of General Relativity on cosmological scales and constraints on dark energy. An essential ingredient is an accurate description of the pairwise velocities of biased tracers of the matter field. In this paper we compute the first and second moments of the pairwise velocity distribution by extending the Convolution Lagrangian Perturbation theory (CLPT) formalism of Carlson et al. (2012). Our predictions outperform standard perturbation theory calculations in many cases when compared to statistics measured in N-body simulations. We combine the CLPT predictions of real-space clustering and velocity statistics in the Gaussian streaming model of Reid & White (2011) to obtain predictions for the monopole and quadrupole correlation functions accurate to 2 and 4 per cent respectively down to <25Mpc/h for halos hosting the massive galaxies observed by SDSS-III BOSS. We also discuss contours of the 2D correlation function and clustering "wedges". We generalize the scheme to cross-correlation functions.

preprint2013arXiv

Baryon Acoustic Oscillations in the Ly-α forest of BOSS quasars

We report a detection of the baryon acoustic oscillation (BAO) feature in the three-dimensional correlation function of the transmitted flux fraction in the \Lya forest of high-redshift quasars. The study uses 48,640 quasars in the redshift range $2.1\le z \le 3.5$ from the Baryon Oscillation Spectroscopic Survey (BOSS) of the third generation of the Sloan Digital Sky Survey (SDSS-III). At a mean redshift $z=2.3$, we measure the monopole and quadrupole components of the correlation function for separations in the range $20\hMpc<r<200\hMpc$. A peak in the correlation function is seen at a separation equal to $(1.01\pm0.03)$ times the distance expected for the BAO peak within a concordance $Λ$CDM cosmology. This first detection of the BAO peak at high redshift, when the universe was strongly matter dominated, results in constraints on the angular diameter distance $\da$ and the expansion rate $H$ at $z=2.3$ that, combined with priors on $H_0$ and the baryon density, require the existence of dark energy. Combined with constraints derived from Cosmic Microwave Background (CMB) observations, this result implies $H(z=2.3)=(224\pm8){\rm km\,s^{-1}Mpc^{-1}}$, indicating that the time derivative of the cosmological scale parameter $\dot{a}=H(z=2.3)/(1+z)$ is significantly greater than that measured with BAO at $z\sim0.5$. This demonstrates that the expansion was decelerating in the range $0.7<z<2.3$, as expected from the matter domination during this epoch. Combined with measurements of $H_0$, one sees the pattern of deceleration followed by acceleration characteristic of a dark-energy dominated universe.

preprint2013arXiv

Dark Energy: A Short Review

The accelerating expansion of the universe is the most surprising cosmological discovery in many decades. In this short review, we briefly summarize theories for the origin of cosmic acceleration and the observational methods being used to test these theories. We then discuss the current observational state of the field, with constraints from the cosmic microwave background (CMB), baryon acoustic oscillations (BAO), Type Ia supernovae (SN), direct measurements of the Hubble constant ($H_0$), and measurements of galaxy and matter clustering. Assuming a flat universe and dark energy with a constant equation-of-state parameter $w = P/ρ$, the combination of Planck CMB temperature anisotropies, WMAP CMB polarization, the Union2.1 SN compilation, and a compilation of BAO measurements yields $w = -1.10^{+0.08}_{-0.07}$, consistent with a cosmological constant ($w=-1$). However, with these constraints the cosmological constant model predicts a value of $H_0$ that is lower than several of the leading recent estimates, and it predicts a parameter combination $σ_8(Ω_m)^{0.5}$ that is higher than many estimates from weak gravitational lensing, galaxy clusters, and redshift-space distortions. Individually these tensions are only significant at the ~$2σ$ level, but they arise in multiple data sets with independent statistics and distinct sources of systematic uncertainty. The tensions are stronger with Planck CMB data than they were with WMAP because of the smaller statistical errors and the higher central value of $Ω_m.$ With the improved measurements expected from the next generation of data sets, these tensions may diminish, or they may sharpen in a way that points towards a more complete physical understanding of cosmic acceleration.

preprint2013arXiv

Mock galaxy catalogs using the quick particle mesh method

Sophisticated analysis of modern large-scale structure surveys requires mock catalogs. Mock catalogs are used to optimize survey design, test reduction and analysis pipelines, make theoretical predictions for basic observables and propagate errors through complex analysis chains. We present a new method, which we call "quick particle mesh", for generating many large-volume, approximate mock catalogs at low computational cost. The method is based on using rapid, low-resolution particle mesh simulations that accurately reproduce the large-scale dark matter density field. Particles are sampled from the density field based on their local density such that they have N-point statistics nearly equivalent to the halos resolved in high-resolution simulations, creating a set of mock halos that can be populated using halo occupation methods to create galaxy mocks for a variety of possible target classes.

preprint2013arXiv

Should we still believe in constrained supersymmetry?

We calculate Bayes factors to quantify how the feasibility of the constrained minimal supersymmetric standard model (CMSSM) has changed in the light of a series of observations. This is done in the Bayesian spirit where probability reflects a degree of belief in a proposition and Bayes' theorem tells us how to update it after acquiring new information. Our experimental baseline is the approximate knowledge that was available before LEP, and our comparison model is the Standard Model with a simple dark matter candidate. To quantify the amount by which experiments have altered our relative belief in the CMSSM since the baseline data we compute the Bayes factors that arise from learning in sequence the LEP Higgs constraints, the XENON100 dark matter constraints, the 2011 LHC supersymmetry search results, and the early 2012 LHC Higgs search results. We find that LEP and the LHC strongly shatter our trust in the CMSSM (with $M_0$ and $M_{1/2}$ below 2 TeV), reducing its posterior odds by a factor of approximately two orders of magnitude. This reduction is largely due to substantial Occam factors induced by the LEP and LHC Higgs searches.

preprint2013arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological implications of the full shape of the clustering wedges in the data release 10 and 11 galaxy samples

We explore the cosmological implications of the angle-averaged correlation function, xi(s), and the clustering wedges, xi_perp(s) and xi_para(s), of the LOWZ and CMASS galaxy samples from Data Release 10 and 11 of the SDSS-III Baryon Oscillation Spectroscopic Survey. Our results show no significant evidence for a deviation from the standard LCDM model. The combination of the information from our clustering measurements with recent data from the cosmic microwave background is sufficient to constrain the curvature of the Universe to Omega_k = 0.0010 +- 0.0029, the total neutrino mass to Sum m_nu < 0.23 eV (95% confidence level), the effective number of relativistic species to N_eff=3.31 +- 0.27, and the dark energy equation of state to w_DE = -1.051 +- 0.076. These limits are further improved by adding information from type Ia supernovae and baryon acoustic oscillations from other samples. In particular, this data set combination is completely consistent with a time-independent dark energy equation of state, in which case we find w_DE=-1.024 +- 0.052. We explore the constraints on the growth-rate of cosmic structures assuming f(z)=Omega_m(z)^gamma and obtain gamma=0.69 +- 0.15, in agreement with the predictions from general relativity of gamma=0.55.

preprint2013arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Luminosity and Color Dependence and Redshift Evolution

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 γ=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.

preprint2013arXiv

The Coyote Universe III: Simulation Suite and Precision Emulator for the Nonlinear Matter Power Spectrum

Many of the most exciting questions in astrophysics and cosmology, including the majority of observational probes of dark energy, rely on an understanding of the nonlinear regime of structure formation. In order to fully exploit the information available from this regime and to extract cosmological constraints, accurate theoretical predictions are needed. Currently such predictions can only be obtained from costly, precision numerical simulations. This paper is the third in a series aimed at constructing an accurate calibration of the nonlinear mass power spectrum on Mpc scales for a wide range of currently viable cosmological models, including dark energy. The first two papers addressed the numerical challenges, and the scheme by which an interpolator was built from a carefully chosen set of cosmological models. In this paper we introduce the "Coyote Univers"' simulation suite which comprises nearly 1,000 N-body simulations at different force and mass resolutions, spanning 38 wCDM cosmologies. This large simulation suite enables us to construct a prediction scheme, or emulator, for the nonlinear matter power spectrum accurate at the percent level out to k~1 h/Mpc. We describe the construction of the emulator, explain the tests performed to ensure its accuracy, and discuss how the central ideas may be extended to a wider range of cosmological models and applications. A power spectrum emulator code is released publicly as part of this paper.

preprint2013arXiv

The large-scale Quasar-Lyman α Forest Cross-Correlation from BOSS

We measure the large-scale cross-correlation of quasars with the Lyman α forest absorption in redshift space, using ~ 60000 quasar spectra from Data Release 9 (DR9) of the Baryon Oscillation Spectroscopic Survey (BOSS). The cross-correlation is detected over a wide range of scales, up to comoving separations r of 80 Mpc/h. For r > 15 Mpc/h, we show that the cross-correlation is well fitted by the linear theory prediction for the mean overdensity around a quasar host halo in the standard ΛCDM model, with the redshift distortions indicative of gravitational evolution detected at high confidence. Using previous determinations of the Lyman α forest bias factor obtained from the Lyman α autocorrelation, we infer the quasar bias factor to be b_q = 3.64^+0.13_-0.15 at a mean redshift z=2.38, in agreement with previous measurements from the quasar auto-correlation. We also obtain a new estimate of the Lyman α forest redshift distortion factor, β_F = 1.1 +/- 0.15, slightly larger than but consistent with the previous measurement from the Lyman α forest autocorrelation. The simple linear model we use fails at separations r < 15 Mpc/h, and we show that this may reasonably be due to the enhanced ionization due to radiation from the quasars. We also provide the expected correction that the mass overdensity around the quasar implies for measurements of the ionizing radiation background from the line-of-sight proximity effect.

preprint2013arXiv

The z=5 Quasar Luminosity Function from SDSS Stripe 82

We present a measurement of the Type I quasar luminosity function at z=5 using a large sample of spectroscopically confirmed quasars selected from optical imaging data. We measure the bright end (M_1450<-26) with Sloan Digital Sky Survey (SDSS) data covering ~6000 deg^2, then extend to lower luminosities (M_1450<-24) with newly discovered, faint z~5 quasars selected from 235 deg^2 of deep, coadded imaging in the SDSS Stripe 82 region (the celestial equator in the Southern Galactic Cap). The faint sample includes 14 quasars with spectra obtained as ancillary science targets in the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), and 59 quasars observed at the MMT and Magellan telescopes. We construct a well-defined sample of 4.7<z<5.1 quasars that is highly complete, with 73 spectroscopic identifications out of 92 candidates. Our color selection method is also highly efficient: of the 73 spectra obtained, 71 are high redshift quasars. These observations reach below the break in the luminosity function (M_1450* ~ -27). The bright end slope is steep (beta <~ -4), with a constraint of beta < -3.1 at 95% confidence. The break luminosity appears to evolve strongly at high redshift, providing an explanation for the flattening of the bright end slope reported previously. We find a factor of ~2 greater decrease in the number density of luminous quasars (M_1450<-26) from z=5 to z=6 than from z=4 to z=5, suggesting a more rapid decline in quasar activity at high redshift than found in previous surveys. Our model for the quasar luminosity function predicts that quasars generate ~30% of the ionizing photons required to keep the universe ionized at z=5.

preprint2012arXiv

A Simple Model for Quasar Demographics

We present a simple model for the relationship between quasars, galaxies, and dark matter halos from 0.5<z<6. In the model, black hole (BH) mass is linearly related to galaxy mass, and galaxies are connected to dark matter halos via empirically constrained relations. A simple "scattered" light bulb model for quasars is adopted, wherein BHs shine at a fixed fraction of the Eddington luminosity during accretion episodes, and Eddington ratios are drawn from a lognormal distribution that is redshift-independent. This model has two free, physically meaningful parameters at each redshift: the normalization of the Mbh-Mgal relation and the quasar duty cycle; these parameters are fit to the observed quasar luminosity function (LF) over the interval 0.5<z<6. This simple model provides an excellent fit to the LF at all epochs, and also successfully predicts the observed projected two-point correlation of quasars from 0.5<z<2.5. It is significant that a single quasar duty cycle at each redshift is capable of reproducing the extant observations. The data are therefore consistent with a scenario wherein quasars are equally likely to exist in galaxies, and therefore dark matter halos, over a wide range in masses. The knee in the quasar LF is a reflection of the knee in the stellar mass-halo mass relation. Future constraints on the quasar LF and quasar clustering at high redshift will provide strong constraints on the model. In the model, the autocorrelation function of quasars becomes a strong function of luminosity only at the very highest luminosities, and will be difficult to observe because such quasars are so rare. Cross-correlation techniques may provide useful constraints on the bias of such rare objects.

preprint2012arXiv

Acoustic scale from the angular power spectra of SDSS-III DR8 photometric luminous galaxies

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

preprint2012arXiv

Clustering of Sloan Digital Sky Survey III Photometric Luminous Galaxies: The Measurement, Systematics and Cosmological Implications

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 Ω_Λ= 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 $Ω_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)

preprint2012arXiv

Convolution Lagrangian perturbation theory for biased tracers

We present a new formulation of Lagrangian perturbation theory which allows accurate predictions of the real- and redshift-space correlation functions of the mass field and dark matter halos. Our formulation involves a non-perturbative resummation of Lagrangian perturbation theory and indeed can be viewed as a partial resummation of the formalism of Matsubara (2008a,b) in which we keep exponentiated all of the terms which tend to a constant at large separation. One of the key features of our method is that we naturally recover the Zel'dovich approximation as the lowest order of our expansion for the matter correlation function. We compare our results against a suite of N-body simulations and obtain good agreement for the correlation functions in real-space and for the monopole correlation function in redshift space. The agreement becomes worse for higher multipole moments of the redshift-space, halo correlation function. Our formalism naturally includes non-linear bias and explains the strong bias-dependence of the multipole moments of the redshift-space correlation function seen in N-body simulations.

preprint2012arXiv

Cross-Correlation of SDSS DR7 Quasars and DR10 BOSS Galaxies: The Weak Luminosity Dependence of Quasar Clustering at z~0.5

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]

preprint2012arXiv

Natural gauge mediation with a bino NLSP at the LHC

Natural models of supersymmetry with a gravitino LSP provide distinctive signatures at the LHC. For a neutralino NLSP, sparticles can decay to two high energy photons plus missing energy. We use the ATLAS diphoton search with 4.8 fb^{-1} of data to place limits in both the stop-gluino and neutralino-chargino mass planes for this scenario. If the neutralino is heavier than 50 GeV, the lightest stop must be heavier than 580 GeV, the gluino heavier than 1100 GeV and charginos must be heavier than approximately 300-470 GeV. This provides the first nontrivial constraints in natural gauge mediation models with a neutralino NLSP decaying to photons, and implies a fine tuning of at least a few percent in such models.

preprint2012arXiv

New Neutrino Mass Bounds from Sloan Digital Sky Survey III Data Release 8 Photometric Luminous Galaxies

We present neutrino mass bounds using 900,000 luminous galaxies with photometric redshifts measured from Sloan Digital Sky Survey III Data Release Eight (SDSS DR8). The galaxies have photometric redshifts between $z = 0.45$ and $z = 0.65$, and cover 10,000 square degrees and thus probe a volume of 3$h^{-3}$Gpc$^3$, enabling tight constraints to be derived on the amount of dark matter in the form of massive neutrinos. A new bound on the sum of neutrino masses $\sum m_ν< 0.26$ eV, at 95% confidence level (CL), is obtained after combining our sample of galaxies, which we call "CMASS", with WMAP 7 year Cosmic Microwave Background (CMB) data and the most recent measurement of the Hubble parameter from the Hubble Space Telescope (HST). This constraint is obtained with a conservative multipole range choice of $30 < \ell < 200$ in order to minimize non-linearities, and a free bias parameter in each of the four redshift bins. We study the impact of assuming this linear galaxy bias model using mock catalogs, and find that this model causes a small ($\sim 1-1.5 σ$) bias in $Ω_{\rm DM} h^2$. For this reason, we also quote neutrino bounds based on a conservative galaxy bias model containing additional, shot noise-like free parameters. In this conservative case, the bounds are significantly weakened, e.g. $\sum m_ν< 0.36$ eV (95% confidence level) for WMAP+HST+CMASS ($\ell_{\rm max}=200$). We also study the dependence of the neutrino bound on multipole range ($\ell_{\rm max}=150$ vs $\ell_{\rm max}=200$) and on which combination of data sets is included as a prior. The addition of supernova and/or Baryon Acoustic Oscillation data does not significantly improve the neutrino mass bound once the HST prior is included. [abridged]

preprint2012arXiv

The Baryon Oscillation Spectroscopic Survey of SDSS-III

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.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: a large sample of mock galaxy catalogues

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.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Analysis of potential systematics

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)

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the Data Release 9 Spectroscopic Galaxy Sample

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σin both the correlation function and power spectrum. Combining with the SDSS-II Luminous Red Galaxy Sample, the detection significance increases to 6.7σ. 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.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological implications of the large-scale two-point correlation function

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.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: measurements of the growth of structure and expansion rate at z=0.57 from anisotropic clustering

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 Ω_Λ = 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.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: the low redshift sample

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.

preprint2012arXiv

The Clustering of Galaxies in the SDSS-III DR9 Baryon Oscillation Spectroscopic Survey: Testing Deviations from $Λ$ and General Relativity using anisotropic clustering of galaxies

We use the joint measurement of geometry and growth from anisotropic galaxy clustering in the Baryon Oscillation Spectroscopic Survey Data Release 9 CMASS sample reported by Reid et al. to constrain dark energy properties and possible deviations from the General Relativity. Assuming GR and taking a prior on the linear matter power spectrum at high redshift from the cosmic microwave background (CMB), anisotropic clustering of the CMASS DR9 galaxies alone constrains $Ω_{\rm m} = 0.308 \pm 0.022$ and $100Ω_{\rm k} = 5.9 \pm 4.8$ for $w = -1$, or $w = -0.91 \pm 0.12$ for $Ω_k = 0$. When combined with the full CMB likelihood, the addition of the anisotropic clustering measurements to the spherically-averaged BAO location increases the constraining power on dark energy by a factor of 4 in a flat CDM cosmology with constant dark energy equation of state $w$ (giving $w = -0.87 \pm 0.05$). This impressive gain depends on our measurement of both the growth of structure and Alcock-Paczynski effect, and is not realised when marginalising over the amplitude of redshift space distortions. Combining with both the CMB and Supernovae Type Ia (SNeIa), we find $Ω_{\rm m} = 0.281 \pm 0.014$ and $1000Ω_{\rm k}=-9.2\pm5.0$ for $w = -1$, or $w_0 = -1.13 \pm 0.12$ and $w_{\rm a}=0.65 \pm 0.36$ assuming $Ω_k = 0$. Finally, when a $Λ$CDM background expansion is assumed, the combination of our estimate of the growth rate with previous growth measurements provides tight constraints on the parameters describing possible deviations from GR giving $γ= 0.64 \pm 0.05$. For one parameter extensions of the flat $Λ$CDM model, we find a $\sim 2σ$ preference either for $w > -1$ or slower growth than in GR. However, the data is fully consistent with the concordance model, and the evidence for these additional parameters is weaker than $2σ$.

preprint2012arXiv

The clustering of intermediate redshift quasars as measured by the Baryon Oscillation Spectroscopic Survey

We measure the quasar two-point correlation function over the redshift range 2.2<z<2.8 using data from the Baryon Oscillation Spectroscopic Survey. We use a homogeneous subset of the data consisting of 27,129 quasars with spectroscopic redshifts---by far the largest such sample used for clustering measurements at these redshifts to date. The sample covers 3,600 square degrees, corresponding to a comoving volume of 9.7(Gpc/h)^3 assuming a fiducial LambdaCDM cosmology, and it has a median absolute i-band magnitude of -26, k-corrected to z=2. After accounting for redshift errors we find that the redshift space correlation function is fit well by a power-law of slope -2 and amplitude s_0=(9.7\pm 0.5)Mpc/h over the range 3<s<25Mpc/h. The projected correlation function, which integrates out the effects of peculiar velocities and redshift errors, is fit well by a power-law of slope -1 and r_0=(8.4\pm 0.6)Mpc/h over the range 4<R<16Mpc/h. There is no evidence for strong luminosity or redshift dependence to the clustering amplitude, in part because of the limited dynamic range in our sample. Our results are consistent with, but more precise than, previous measurements at similar redshifts. Our measurement of the quasar clustering amplitude implies a bias factor of b~3.5 for our quasar sample. We compare the data to models to constrain the manner in which quasars occupy dark matter halos at z~2.4 and infer that such quasars inhabit halos with a characteristic mass of <M>~10^{12}Msun/h with a duty cycle for the quasar activity of 1 per cent.

preprint2012arXiv

The large-scale cross-correlation of Damped Lyman Alpha Systems with the Lyman Alpha Forest: First Measurements from BOSS

We present the first measurement of the large-scale cross-correlation of Lyman alpha forest absorption and Damped Lyman alpha systems (DLA), using the 9th Data Release of the Baryon Oscillation Spectroscopic Survey (BOSS). The cross-correlation is clearly detected on scales up to 40 Mpc/h and is well fitted by the linear theory prediction of the standard Cold Dark Matter model of structure formation with the expected redshift distortions, confirming its origin in the gravitational evolution of structure. The amplitude of the DLA-Lyman alpha cross-correlation depends on only one free parameter, the bias factor of the DLA systems, once the Lyman alpha forest bias factors are known from independent Lyman alpha forest correlation measurements. We measure the DLA bias factor to be b_D = (2.17 +/- 0.20) beta_F^{0.22}, where the Lyman alpha forest redshift distortion parameter beta_F is expected to be above unity. This bias factor implies a typical host halo mass for DLAs that is much larger than expected in present DLA models, and is reproduced if the DLA cross section scales with halo mass as M_h^alpha, with alpha= 1.1 +/- 0.1 for beta_F=1. Matching the observed DLA bias factor and rate of incidence requires that atomic gas remains extended in massive halos over larger areas than predicted in present simulations of galaxy formation, with typical DLA proper sizes larger than 20 kpc in host halos of masses ~ 10^12 solar masses. We infer that typical galaxies at z ~ 2 to 3 are surrounded by systems of atomic clouds that are much more extended than the luminous parts of galaxies and contain ~ 10% of the baryons in the host halo.

preprint2012arXiv

The Ninth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Baryon Oscillation Spectroscopic Survey

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.

preprint2012arXiv

The SDSS-III Baryon Oscillation Spectroscopic Survey: The Quasar Luminosity Function from Data Release Nine

We present a new measurement of the optical Quasar Luminosity Function (QLF), using data from the Sloan Digital Sky Survey-III: Baryon Oscillation Spectroscopic Survey (SDSS-III: BOSS). From the SDSS-III Data Release Nine (DR9), we select a uniform sample of 22,301 i<=21.8 quasars over an area of 2236 sq. deg with confirmed spectroscopic redshifts between 2.2<z<3.5, filling in a key part of the luminosity-redshift plane for optical quasar studies. We derive the completeness of the survey through simulated quasar photometry, and check this completeness estimate using a sample of quasars selected by their photometric variability within the BOSS footprint. We investigate the level of systematics associated with our quasar sample using the simulations, in the process generating color-redshift relations and a new quasar k-correction. We probe the faint end of the QLF to M_i(z=2.2) = -24.5 and see a clear break in the QLF at all redshifts up to z=3.5. We find that a log-linear relation (in log[Phi*] - M*) for a luminosity and density evolution (LEDE) model adequately describes our data within the range 2.2<z<3.5; across this interval the break luminosity increases by a factor of ~2.3 while Phi* declines by a factor of ~6. At z<2.2 our data is reasonably well fit by a pure luminosity evolution (PLE) model. We see only a weak signature of "AGN downsizing", in line with recent studies of the hard X-ray luminosity function. We compare our measured QLF to a number of theoretical models and find that models making a variety of assumptions about quasar triggering and halo occupation can fit our data over a wide range of redshifts and luminosities.

preprint2012arXiv

The Sloan Digital Sky Survey quasar catalog: ninth data release

We present the Data Release 9 Quasar (DR9Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the survey, are spectrocopically confirmed as quasars via visual inspection, have luminosities Mi[z=2]<-20.5 (in a $Λ$CDM cosmology with H0 = 70 km/s/Mpc, $Ω_{\rm M}$ = 0.3, and $Ω_Λ$ = 0.7) and either display at least one emission line with full width at half maximum (FWHM) larger than 500 km/s or, if not, have interesting/complex absorption features. It includes as well, known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. This catalog contains 87,822 quasars (78,086 are new discoveries) detected over 3,275 deg$^{2}$ with robust identification and redshift measured by a combination of principal component eigenspectra newly derived from a training set of 8,632 spectra from SDSS-DR7. The number of quasars with $z>2.15$ (61,931) is ~2.8 times larger than the number of z>2.15 quasars previously known. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 7,533 broad absorption line quasars and gives their characteristics. For each object the catalog presents five-band (u,g,r,i,z) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys.

preprint2011arXiv

Ameliorating Systematic Uncertainties in the Angular Clustering of Galaxies: A Study using SDSS-III

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.

preprint2011arXiv

Mass Function Predictions Beyond LCDM

The mass distribution of halos, as specified by the halo mass function, is a key input for several cosmological probes. The sizes of $N$-body simulations are now such that, for the most part, results need no longer be statistics-limited, but are still subject to various systematic uncertainties. We investigate and discuss some of the reasons for these differences. Quantifying error sources and compensating for them as appropriate, we carry out a high-statistics study of dark matter halos from 67 $N$-body simulations to investigate the mass function and its evolution for a reference $Λ$CDM cosmology and for a set of $w$CDM cosmologies. For the reference $Λ$CDM cosmology (close to WMAP5), we quantify the breaking of universality in the form of the mass function as a function of redshift, finding an evolution of as much as 10% away from the universal form between redshifts $z=0$ and $z=2$. For cosmologies very close to this reference we provide a fitting formula to our results for the (evolving) $Λ$CDM mass function over a mass range of $6\cdot 10^{11}-3\cdot 10^{15}$ M$_{\odot}$ to an estimated accuracy of about 2%. The set of $w$CDM cosmologies is taken from the Coyote Universe simulation suite. The mass functions from this suite (which includes a $Λ$CDM cosmology and others with $w\simeq-1$) are described by the fitting formula for the reference $Λ$CDM case at an accuracy level of 10%, but with clear systematic deviations. We argue that, as a consequence, fitting formulae based on a universal form for the mass function may have limited utility in high precision cosmological applications.

preprint2011arXiv

SDSS-III: Massive Spectroscopic Surveys of the Distant Universe, the Milky Way Galaxy, and Extra-Solar Planetary Systems

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)

preprint2011arXiv

The Eighth Data Release of the Sloan Digital Sky Survey: First Data from SDSS-III

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.

preprint2011arXiv

The SDSS-III Baryon Oscillation Spectroscopic Survey: Quasar Target Selection for Data Release Nine

The SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), a five-year spectroscopic survey of 10,000 deg^2, achieved first light in late 2009. One of the key goals of BOSS is to measure the signature of baryon acoustic oscillations in the distribution of Ly-alpha absorption from the spectra of a sample of ~150,000 z>2.2 quasars. Along with measuring the angular diameter distance at z\approx2.5, BOSS will provide the first direct measurement of the expansion rate of the Universe at z > 2. One of the biggest challenges in achieving this goal is an efficient target selection algorithm for quasars over 2.2 < z < 3.5, where their colors overlap those of stars. During the first year of the BOSS survey, quasar target selection methods were developed and tested to meet the requirement of delivering at least 15 quasars deg^-2 in this redshift range, out of 40 targets deg^-2. To achieve these surface densities, the magnitude limit of the quasar targets was set at g <= 22.0 or r<=21.85. While detection of the BAO signature in the Ly-alpha absorption in quasar spectra does not require a uniform target selection, many other astrophysical studies do. We therefore defined a uniformly-selected subsample of 20 targets deg^-2, for which the selection efficiency is just over 50%. This "CORE" subsample will be fixed for Years Two through Five of the survey. In this paper we describe the evolution and implementation of the BOSS quasar target selection algorithms during the first two years of BOSS operations. We analyze the spectra obtained during the first year. 11,263 new z>2.2 quasars were spectroscopically confirmed by BOSS. Our current algorithms select an average of 15 z > 2.2 quasars deg^-2 from 40 targets deg^-2 using single-epoch SDSS imaging. Multi-epoch optical data and data at other wavelengths can further improve the efficiency and completeness of BOSS quasar target selection. [Abridged]

preprint2011arXiv

Towards an accurate model of the redshift space clustering of halos in the quasilinear regime

Observations of redshift-space distortions in spectroscopic galaxy surveys offer an attractive method for measuring the build-up of cosmological structure, which depends both on the expansion rate of the Universe and our theory of gravity. Galaxies occupy dark matter halos, whose redshift space clustering has a complex dependence on bias that cannot be inferred from the behavior of matter. We identify two distinct corrections on quasilinear scales (~ 30-80 Mpc/h): the non-linear mapping between real and redshift space positions, and the non-linear suppression of power in the velocity divergence field. We model the first non-perturbatively using the scale-dependent Gaussian streaming model, which we show is accurate at the <0.5 (2) per cent level in transforming real space clustering and velocity statistics into redshift space on scales s>10 (s>25) Mpc/h for the monopole (quadrupole) halo correlation functions. We use perturbation theory to predict the real space pairwise halo velocity statistics. Our fully analytic model is accurate at the 2 per cent level only on scales s > 40 Mpc/h. Recent models that neglect the corrections from the bispectrum and higher order terms from the non-linear real-to-redshift space mapping will not have the accuracy required for current and future observational analyses. Finally, we note that our simulation results confirm the essential but non-trivial assumption that on large scales, the bias inferred from real space clustering of halos is the same one that determines their pairwise infall velocity amplitude at the per cent level.

preprint2010arXiv

Cluster Galaxy Dynamics and the Effects of Large Scale Environment

We use a high-resolution N-body simulation to study how the influence of large-scale structure in and around clusters causes correlated signals in different physical probes and discuss some implications this has for multi-physics probes of clusters. We pay particular attention to velocity dispersions, matching galaxies to subhalos which are explicitly tracked in the simulation. We find that not only do halos persist as subhalos when they fall into a larger host, groups of subhalos retain their identity for long periods within larger host halos. The highly anisotropic nature of infall into massive clusters, and their triaxiality, translates into an anisotropic velocity ellipsoid: line-of-sight galaxy velocity dispersions for any individual halo show large variance depending on viewing angle. The orientation of the velocity ellipsoid is correlated with the large-scale structure, and thus velocity outliers correlate with outliers caused by projection in other probes. We quantify this orientation uncertainty and give illustrative examples. Such a large variance suggests that velocity dispersion estimators will work better in an ensemble sense than for any individual cluster, which may inform strategies for obtaining redshifts of cluster members. We similarly find that the ability of substructure indicators to find kinematic substructures is highly viewing angle dependent. While groups of subhalos which merge with a larger host halo can retain their identity for many Gyr, they are only sporadically picked up by substructure indicators. (Abridged)

preprint2010arXiv

Embedding realistic surveys in simulations through volume remapping

Connecting cosmological simulations to real-world observational programs is often complicated by a mismatch in geometry: while surveys often cover highly irregular cosmological volumes, simulations are customarily performed in a periodic cube. We describe a technique to remap this cube into elongated box-like shapes that are more useful for many applications. The remappings are one-to-one, volume-preserving, keep local structures intact, and involve minimal computational overhead.

preprint2010arXiv

High-precision predictions for the acoustic scale in the non-linear regime

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)

preprint2010arXiv

Particle mesh simulations of the Lyman-alpha forest and the signature of Baryon Acoustic Oscillations in the intergalactic medium

We present a set of ultra-large particle-mesh simulations of the LyA forest targeted at understanding the imprint of baryon acoustic oscillations (BAO) in the inter-galactic medium. We use 9 dark matter only simulations which can, for the first time, simultaneously resolve the Jeans scale of the intergalactic gas while covering the large volumes required to adequately sample the acoustic feature. Mock absorption spectra are generated using the fluctuating Gunn-Peterson approximation which have approximately correct flux probability density functions (PDFs) and small-scale power spectra. On larger scales there is clear evidence in the redshift space correlation function for an acoustic feature, which matches a linear theory template with constant bias. These spectra, which we make publicly available, can be used to test pipelines, plan future experiments and model various physical effects. As an illustration we discuss the basic properties of the acoustic signal in the forest, the scaling of errors with noise and source number density, modified statistics to treat mean flux evolution and misestimation, and non-gravitational sources such as fluctuations in the photo-ionizing background and temperature fluctuations due to HeII reionization.

preprint2010arXiv

Shot noise and reconstruction of the acoustic peak

We study the effect of noise in the density field, such as would arise from a finite number density of tracers, on reconstruction of the acoustic peak within the context of Lagrangian perturbation theory. Reconstruction performs better when the density field is determined from denser tracers, but the gains saturate at n~1e-4(h/Mpc)^3. For low density tracers it is best to use a large smoothing scale to define the shifts, but the optimum is very broad.

preprint2010arXiv

The clustering of massive galaxies at z~0.5 from the first semester of BOSS data

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.

preprint2010arXiv

The Coyote Universe I: Precision Determination of the Nonlinear Matter Power Spectrum

Near-future cosmological observations targeted at investigations of dark energy pose stringent requirements on the accuracy of theoretical predictions for the clustering of matter. Currently, N-body simulations comprise the only viable approach to this problem. In this paper we demonstrate that N-body simulations can indeed be sufficiently controlled to fulfill these requirements for the needs of ongoing and near-future weak lensing surveys. By performing a large suite of cosmological simulation comparison and convergence tests we show that results for the nonlinear matter power spectrum can be obtained at 1% accuracy out to k~1 h/Mpc. The key components of these high accuracy simulations are: precise initial conditions, very large simulation volumes, sufficient mass resolution, and accurate time stepping. This paper is the first in a series of three, with the final aim to provide a high-accuracy prediction scheme for the nonlinear matter power spectrum.

preprint2010arXiv

The Coyote Universe II: Cosmological Models and Precision Emulation of the Nonlinear Matter Power Spectrum

The power spectrum of density fluctuations is a foundational source of cosmological information. Precision cosmological probes targeted primarily at investigations of dark energy require accurate theoretical determinations of the power spectrum in the nonlinear regime. To exploit the observational power of future cosmological surveys, accuracy demands on the theory are at the one percent level or better. Numerical simulations are currently the only way to produce sufficiently error-controlled predictions for the power spectrum. The very high computational cost of (precision) N-body simulations is a major obstacle to obtaining predictions in the nonlinear regime, while scanning over cosmological parameters. Near-future observations, however, are likely to provide a meaningful constraint only on constant dark energy equation of state 'wCDM' cosmologies. In this paper we demonstrate that a limited set of only 37 cosmological models -- the "Coyote Universe" suite -- can be used to predict the nonlinear matter power spectrum at the required accuracy over a prior parameter range set by cosmic microwave background observations. This paper is the second in a series of three, with the final aim to provide a high-accuracy prediction scheme for the nonlinear matter power spectrum for wCDM cosmologies.

preprint2010arXiv

What determines satellite galaxy disruption?

In hierarchical structure formation, dark matter halos that merge with larger halos can persist as subhalos. These subhalos are likely hosts of visible galaxies. While the dense halo environment rapidly strips subhalos of their dark mass, the compact luminous material can remain intact for some time, making the correspondence of galaxies with severely stripped subhalos unclear. Many galaxy evolution models assume that satellite galaxies eventually merge with their central galaxy, but this ignores the possibility of satellite tidal disruption. We use a high-resolution $N$-body simulation of cosmological volume to explore satellite galaxy merging/disruption criteria based on dark matter subhalo dynamics. We explore the impact that satellite merging/disruption has on the Halo Occupation Distribution and radial profile of the remnants. Using abundance matching to assign stellar mass/luminosity to subhalos, we compare with observed galaxy clustering, satellite fractions, cluster satellite luminosity functions, finding that subhalos reproduce well these observables on all scales. Our results imply that satellite subhalos corresponding to >0.2L* galaxies must be well-resolved down to 1-3% of their mass at infall to robustly trace the galaxy population. We also explore a simple analytic model based on dynamical friction for satellite galaxy infall, finding good agreement with our subhalo catalog and observations.

preprint2009arXiv

Incorporating Photometric Redshift Probability Density Information into Real-Space Clustering Measurements

The use of photometric redshifts in cosmology is increasing. Often, however these photo-zs are treated like spectroscopic observations, in that the peak of the photometric redshift, rather than the full probability density function (PDF), is used. This overlooks useful information inherent in the full PDF. We introduce a new real-space estimator for one of the most used cosmological statistics, the 2-point correlation function, that weights by the PDF of individual photometric objects in a manner that is optimal when Poisson statistics dominate. As our estimator does not bin based on the PDF peak it substantially enhances the clustering signal by usefully incorporating information from all photometric objects that overlap the redshift bin of interest. As a real-world application, we measure QSO clustering in the Sloan Digital Sky Survey (SDSS). We find that our simplest binned estimator improves the clustering signal by a factor equivalent to increasing the survey size by a factor of 2-3. We also introduce a new implementation that fully weights between pairs of objects in constructing the cross-correlation and find that this pair-weighted estimator improves clustering signal in a manner equivalent to increasing the survey size by a factor of 4-5. Our technique uses spectroscopic data to anchor the distance scale and it will be particularly useful where spectroscopic data (e.g, from BOSS) overlaps deeper photometry (e.g.,from Pan-STARRS, DES or the LSST). We additionally provide simple, informative expressions to determine when our estimator will be competitive with the autocorrelation of spectroscopic objects. Although we use QSOs as an example population, our estimator can and should be applied to any clustering estimate that uses photometric objects.

preprint2009arXiv

Reconstructing baryon oscillations

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.

preprint2008arXiv

Forecasting Cosmological Constraints from Redshift Surveys

preprint2008arXiv

Halo Assembly Bias in Hierarchical Structure Formation

We investigate the origin of halo assembly bias, the dependence of halo clustering on assembly history. We relate halo assembly to peak properties measured in the Lagrangian space of the initial linear Gaussian random density field, and show how these same Lagrangian properties determine large-scale bias. We focus on the two regimes where assembly bias has been observed to be significant: at masses very large and very small compared to the nonlinear mass scale. At high masses, we show that assembly bias is expected from the statistics of the peaks of Gaussian random fluctuations, and we show that the extent of assembly bias found in N-body simulations of rare halos is in excellent agreement with our theoretical prediction. At low masses, we argue that assembly bias largely arises from a sub-population of low mass halos whose mass accretion has ceased. Due to their arrested development, these halos naturally become unbiased, in contrast to their anti-biased peers. We show that a simple toy model incorporating these effects can roughly reproduce the bias trends found in N-body simulations.

preprint2008arXiv

Red Galaxy Growth and the Halo Occupation Distribution

We have traced the past 7 Gyr of red galaxy stellar mass growth within dark matter halos. We have determined the halo occupation distribution, which describes how galaxies reside within dark matter halos, using the observed luminosity function and clustering of 40,696 0.2<z<1.0 red galaxies in Bootes. Half of 10^{11.9} Msun/h halos host a red central galaxy, and this fraction increases with increasing halo mass. We do not observe any evolution of the relationship between red galaxy stellar mass and host halo mass, although we expect both galaxy stellar masses and halo masses to evolve over cosmic time. We find that the stellar mass contained within the red population has doubled since z=1, with the stellar mass within red satellite galaxies tripling over this redshift range. In cluster mass halos most of the stellar mass resides within satellite galaxies and the intra-cluster light, with a minority of the stellar mass residing within central galaxies. The stellar masses of the most luminous red central galaxies are proportional to halo mass to the power of a third. We thus conclude that halo mergers do not always lead to rapid growth of central galaxies. While very massive halos often double in mass over the past 7 Gyr, the stellar masses of their central galaxies typically grow by only 30%.

preprint2004arXiv

Baryons and weak lensing power spectra

We provide a simple toy model to elucidate the effects of baryonic cooling on predictions of the convergence power spectrum in weak gravitational lensing. Our model suggests that the effects of baryonic cooling on dark matter halo profiles can alter the dark matter power spectrum at k~10 h/Mpc and lead to percent level changes in the convergence angular power spectrum at wavenumbers of several thousand. These effects would be measurable with planned weak lensing surveys and will impact constraints on cosmological parameters.

preprint2002arXiv

High redshift galaxies and the Lyman-alpha forest in a CDM universe

We use a cosmological hydrodynamic simulation of a cold dark matter universe to investigate theoretically the relationship between high redshift galaxies and the Lyman=alpha forest at redshift z=3. Galaxies in the simulation are surrounded by halos of hot gas, which nevertheless contain enough neutral hydrogen to cause a Ly-alpha flux decrement, its strength increasing with galaxy mass. A comparison with recent observational data by Adelberger et. al on the Ly-alpha forest around galaxies reveals that actual galaxies may have systematically less Ly-alpha absorption within 1 Mpc of them than our simulated galaxies. In order to investigate this possibility, we add several simple prescriptions for galaxy feedback on the IGM to the evolved simulation. These include the effect of photoionizing background radiation coming from galactic sources, galactic winds whose only effect is to deposit thermal energy into the IGM, and another, kinetic model for winds, which are assumed to evacuate cavities in the IGM around galaxies. We find that only the latter is able to produce a large effect, enough to match the tentative observational data, given the energy available from star formation in the simulated galaxies. Another intriguing possibility is that a selection effect is responsible, so that galaxies with low Ly-alpha absorption are preferentially included in the sample. This is also viable, but predicts very different galaxy properties (including clustering) than the other scenarios.

preprint1998arXiv

Baryonic Signatures in Large-Scale Structure

We investigate the consequences of a non-negligible baryon fraction for models of structure formation in Cold Dark Matter dominated cosmologies, emphasizing in particular the existence of oscillations in the present-day matter power spectrum. These oscillations are the remnants of acoustic oscillations in the photon-baryon fluid before last scattering. For acceptable values of the cosmological and baryon densities, the oscillations modulate the power by up to 10%, with a `period' in spatial wavenumber which is close to Delta k approximately 0.05/ Mpc. We study the effects of nonlinear evolution on these features, and show that they are erased for k > 0.2 h/ Mpc. At larger scales, the features evolve as expected from second-order perturbation theory: the visibility of the oscillations is affected only weakly by nonlinear evolution. No realistic CDM parameter combination is able to account for the claimed feature near k = 0.1 h/ Mpc in the APM power spectrum, or the excess power at 100 Mpc/h wavelengths quoted by several recent surveys. Thus baryonic oscillations are not predicted to dominate existing measurements of clustering. We examine several effects which may mask the features which are predicted, and conclude that future galaxy surveys may be able to detect the oscillatory features in the power spectrum provided baryons comprise more than 15% of the total density, but that it will be a technically challenging achievement.

preprint1996arXiv

Cold Dark Matter Models with a Cosmological Constant

We use linear and quasi-linear perturbation theory to analyse cold dark matter models of structure formation in spatially flat models with a cosmological constant. Both a tilted spectrum of density perturbations and a significant gravitational wave contribution to the microwave anisotropy are allowed as possibilities. We provide normalizations of the models to microwave anisotropies, as given by the four-year {\it COBE} observations, and show how all the normalization information for such models, including tilt, can be condensed into a single fitting function which is independent of the value of the Hubble parameter. We then discuss a wide variety of other types of observations. We find that a very wide parameter space is available for these models, provided $Ω_0$ is greater than about 0.3, and that large-scale structure observations show no preference for any particular value of $Ω_0$ in the range 0.3 to 1.

preprint1996arXiv

Cold dark matter models with high baryon content

Recent results have suggested that the density of baryons in the Universe, OmegaB, is much more uncertain than previously thought, and may be significantly higher. We demonstrate that a higher OmegaB increases the viability of critical-density cold dark matter (CDM) models. High baryon fraction offers the twin benefits of boosting the first peak in the microwave anisotropy power spectrum and of suppressing short-scale power in the matter power spectrum. These enable viable CDM models to have a larger Hubble constant than otherwise possible. We carry out a general exploration of high OmegaB CDM models, varying the Hubble constant h and the spectral index n. We confront a variety of observational constraints and discuss specific predictions. Although some observational evidence may favour baryon fractions as high as 20 per cent, we find that values around 10 to 15 per cent provide a reasonable fit to a wide range of data. We suggest that models with OmegaB in this range, with h about 0.5 and n about 0.8, are currently the best critical-density CDM models.

preprint1995arXiv

Cold Dark Matter Resuscitated?

The Cold Dark Matter (CDM) model has an elegant simplicitly which makes it very predictive, but when its parameters are fixed at their `canonical' values its predictions are in conflict with observational data. There is, however, much leeway in the initial conditions within the CDM framework. We advocate a re-examination of the CDM model, taking into account modest variation of parameters from their canonical values. We find that CDM models with $n=0.8$--0.9 and $h=0.45$--0.50 can fit the available data. Our ``best fit'' CDM model has $n=0.9$, $h=0.45$ and $C_2^{T}/C_2^{S}=0.7$. We discuss the current state of observations which could definitely rule out this model.

preprint1993arXiv

Cosmic Variance in CMB Anisotropies: From $1^{\circ}$ to COBE

Cosmic Microwave Background (CMB) anisotropies that result from quantum fluctuations during inflation are explored and the impact of their ``cosmic variance'' on the ability to use existing data to probe inflationary models is studied. We calculate the rms temperature fluctuation, and its cosmic variance, for a number of experiments and for models with primordial power spectra which range from $n={1\over 2}$ to $1$. We find: (1) cosmic variance obscures the information which can be extracted, so a comparison of the rms temperature fluctuation on small scales with the COBE result can fix $n$ to only $\approx\pm 0.2$ at best; (2) measurements of the rms fluctuation on $1^{\circ}$ scales may not allow one to unambiguously infer the tensor contribution to the COBE anisotropy; (3) comparison of this contribution with the predictions of inflation are ambiguous if the quadrupole anisotropy alone is utilized. We discuss means for minimizing the uncertainty due to cosmic variance in comparisons between experiments.

preprint1993arXiv

Decuplet Contributions to Hyperon Axial Vector Form Factors

We consider the predictions of chiral perturbation theory for $SU(3)$ breaking in the axial vector form factor $g_1$ measured in semileptonic hyperon decays. We confirm that if only octet baryon intermediate states are included, the non-analytic corrections are $\sim 100\%$. These corrections are dominated by an $SU(3)$-symmetric wavefunction renormalization, which explains the fact that the ``corrected'' predictions still fit the data well. We argue that the large corrections are nonetheless strong evidence that this chiral expansion is breaking down. Following a recent suggestion of Jenkins and Manohar, we then include contributions from decuplet baryon intermediate states. Unlike these authors, we do not neglect the octet--decuplet mass difference $Δ$. We find that the effects of $Δ\ne 0$ significantly change the pattern of corrections: we still find that the decuplet corrections can cancel the large octet contributions in a non-trivial way, but the corrections no longer favor the $SU(6)$ values of the axial couplings. We also argue that $D$ and $F$ axial couplings cannot be reliably extracted from calculations which keep only the non-analytic corrections.

preprint1992arXiv

Contribution of Long Wavelength Gravitational Waves to the CMB Anisotropy

We present an in depth discussion of the production of gravitational waves from an inflationary phase that could have occurred in the early universe, giving derivations for the resulting spectrum and energy density. We also consider the large-scale anisotropy in the cosmic microwave background radiation coming from these waves. Assuming that the observed quadrupole anisotropy comes mostly from gravitational waves (consistent with the predictions of a flat spectrum of scalar density perturbations and the measured dipole anisotropy) we describe in detail how to derive a value for the scale of inflation of $(1.5-5)\times 10^{16}$GeV, which is at a particularly interesting scale for particle physics. This upper limit corresponds to a 95\% confidence level upper limit on the scale of inflation assuming only that the quadrupole anisotropy from gravitational waves is not cancelled by another source. Direct detection of gravitational waves produced by inflation near this scale will have to wait for the next generation of detectors.

Martin White

What is connected

Connect this record

See the researcher in context

Building this map preview

150 published item(s)

Second Data Release of the COSMOS Lyman-alpha Mapping and Tomographic Observation: The First 3D Maps of the Detailed Cosmic Web at 2.05<z<2.55

A new analysis of galaxy 2-point functions in the BOSS survey, including full-shape information and post-reconstruction BAO

A Spectroscopic Road Map for Cosmic Frontier: DESI, DESI-II, Stage-5

Cosmological Analysis of Three-Dimensional BOSS Galaxy Clustering and Planck CMB Lensing Cross Correlations via Lagrangian Perturbation Theory

Cosmological constraints from the tomographic cross-correlation of DESI Luminous Red Galaxies and Planck CMB lensing

Cosmology at high redshift -- a probe of fundamental physics

Global fit of 2HDM with future collider results

Priors on red galaxy stochasticity from hybrid effective field theory

Simple and statistically sound recommendations for analysing physical theories

Snowmass 2021 CMB-S4 White Paper

Snowmass2021 Cosmic Frontier White Paper: Cosmology and Fundamental Physics from the three-dimensional Large Scale Structure

Snowmass2021 Cosmic Frontier White Paper: High Density Galaxy Clustering in the Regime of Cosmic Acceleration

Snowmass2021: Opportunities from Cross-survey Analyses of Static Probes

The DESI $N$-body Simulation Project -- II. Suppressing sample variance with fast simulations

The MegaMapper: A Stage-5 Spectroscopic Instrument Concept for the Study of Inflation and Dark Energy

Convergent Bayesian Global Fits of 4D Composite Higgs Models

CosmoBit: A GAMBIT module for computing cosmological observables and likelihoods

Does SUSY have friends? A new approach for LHC event analysis

Neural Network Acceleration of Large-scale Structure Theory Calculations

Thermal WIMPs and the Scale of New Physics: Global Fits of Dirac Dark Matter Effective Field Theories

Consistent Modeling of Velocity Statistics and Redshift-Space Distortions in One-Loop Perturbation Theory

Exploring the 2HDM with Global Fits in $\texttt{GAMBIT}$

GAMBIT and its Application in the Search for Physics Beyond the Standard Model

Global analyses of Higgs portal singlet dark matter models using GAMBIT

Imaging Systematics and Clustering of DESI Main Targets

Minimal 4D Composite Higgs Models Under Current LHC Constraints

Model-independent approach for incorporating interference effects in collider searches for new resonances

Packed Ultra-wideband Mapping Array (PUMA): Astro2020 RFI Response

Reconstruction with velocities

Self-consistent approach for measuring the energy spectra and composition of cosmic rays and determining the properties of hadronic interactions at high energy

Simulations and symmetries

TARDIS Paper I: A Constrained Reconstruction Approach to Modeling the z~2.5 Cosmic Web Probed by Lyman-alpha Forest Tomography

The Reconstructed Power Spectrum in the Zeldovich Approximation

Axion global fits with Peccei-Quinn symmetry breaking before inflation using GAMBIT

ColliderBit: a GAMBIT module for the calculation of high-energy collider observables and likelihoods

Combined collider constraints on neutralinos and charginos

CosmoDC2: A Synthetic Sky Catalog for Dark Energy Science with LSST

Higgs Quark Flavor Violation: Simplified Models and Status of General Two-Higgs-Doublet Model

SpecBit, DecayBit and PrecisionBit: GAMBIT modules for computing mass spectra, particle decay rates and precision observables

Wide angle effects for peculiar velocities

A global fit of the MSSM with GAMBIT

DarkBit: A GAMBIT module for computing dark matter observables and likelihoods

Global fits of GUT-scale SUSY models with GAMBIT

Status of the scalar singlet dark matter model

FlavBit: A GAMBIT module for computing flavour observables and likelihoods

GAMBIT: The Global and Modular Beyond-the-Standard-Model Inference Tool

Combined analysis of effective Higgs portal dark matter models

Combining galaxy and 21cm surveys

Connecting massive galaxies to dark matter halos in BOSS - I. Is galaxy color a stochastic process in high-mass halos?

Cosmological perturbation theory in 1+1 dimensions

Gamma rays from muons from WIMPs: Implementation of radiative muon decays for dark matter analyses

Improving photometric redshifts with Ly$α$ tomography

Including parameter dependence in the data and covariance for cosmological inference

On using angular cross-correlations to determine source redshift distributions

Redshift Weights for Baryon Acoustic Oscillations : Application to Mock Galaxy Catalogs

Revealing the z~2.5 Cosmic Web With 3D Lyman-Alpha Forest Tomography: A Deformation Tensor Approach

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Observational systematics and baryon acoustic oscillations in the correlation function

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the correlation function of LOWZ and CMASS galaxies in Data Release 12

The DESI Experiment Part I: Science,Targeting, and Survey Design

The DESI Experiment Part II: Instrument Design

The Gaussian streaming model and Lagrangian effective field theory

The SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Overview and Early Data

The XMM Cluster Survey: The Halo Occupation Number of BOSS galaxies in X-ray clusters

A Lagrangian effective field theory

Clustering of intermediate redshift quasars using the final SDSS III-BOSS sample

Collider constraints on tuning in composite Higgs models

Dependence of the Cosmic Microwave Background Lensing Power Spectrum on the Matter Density

Estimating sparse precision matrices

Interpreting the Fermi-LAT gamma ray excess in the simplified framework

LUX likelihood and limits on spin-independent and spin-dependent WIMP couplings with LUXCalc

Reconstruction within the Zeldovich approximation

SDSS-III Baryon Oscillation Spectroscopic Survey Data Release 12: galaxy target selection and large scale structure catalogues

Shadow of a Colossus: A z=2.45 Galaxy Protocluster Detected in 3D Ly-a Forest Tomographic Mapping of the COSMOS Field

The Eleventh and Twelfth Data Releases of the Sloan Digital Sky Survey: Final Data from SDSS-III