Source author record

Patrick McDonald

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

41works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

41 published item(s)

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

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

Extremizing Temperature Functions of Rods with Robin Boundary Conditions

We compare the solutions of two one-dimensional Poisson problems on an interval with Robin boundary conditions, one with given data, and one where the data has been symmetrized. When the Robin parameter is positive and the symmetrization is symmetric decreasing rearrangement, we prove that the solution to the symmetrized problem has larger increasing convex means. When the Robin parameter equals zero (so that we have Neumann boundary conditions) and the symmetrization is decreasing rearrangement, we similarly show that the solution to the symmetrized problem has larger convex means.

preprint2020arXiv

Massive neutrinos and degeneracies in Lyman-alpha forest simulations

Using a suite of hydrodynamical simulations with cold dark matter, baryons, and neutrinos, we present a detailed study of the effect of massive neutrinos on the 1-D and 3-D flux power spectra of the Lyman-$α$ (Ly$α$) forest. The presence of massive neutrinos in cosmology induces a scale- and time-dependent suppression of structure formation that is strongest on small scales. Measuring this suppression is a key method for inferring neutrino masses from cosmological data, and is one of the main goals of ongoing and future surveys like eBOSS, DES, LSST, Euclid or DESI. The clustering in the Ly$α$ forest traces the quasi-linear power at late times and on small scales. In combination with observations of the cosmic microwave background, the forest therefore provides some of the tightest constraints on the sum of the neutrino masses. However there is a well-known degeneracy between $Σm_ν$ and the amplitude of perturbations in the linear matter power spectrum. We study the corresponding degeneracy in the 1-D flux power spectrum of the Ly$α$ forest, and for the first time also study this degeneracy in the 3-D flux power spectrum. We show that the non-linear effects of massive neutrinos on the Ly$α$ forest, beyond the effect of linear power amplitude suppression, are negligible, and this degeneracy persists in the Ly$α$ forest observables to a high precision. We discuss the implications of this degeneracy for choosing parametrisations of the Ly$α$ forest for cosmological analysis.

preprint2016arXiv

Cherenkov Radiation with Massive, CPT-violating Photons

The source of CPT-violation in the photon sector of the Standard Model Extension arises from a Chern-Simons-like contribution that involves a coupling to a fixed background vector field $k_{AF}^μ$. These Lorentz- and CPT-violating photons have well-known theoretical issues that arise from missing states at low momenta when $k_{AF}^μ$ is timelike. In order to make the theory consistent, a tiny mass for the photon can be introduced, well below current experimental bounds. The implementation of canonical quantization can then be implemented as in the CPT-preserving case by using the Stückelberg mechanism. We explicitly construct a covariant basis of properly-normalized polarization vectors at fixed three-momentum satisfying the momentum space field equations, in terms of which the vector field can be expanded. As an application of the theory, we calculate the Cherenkov radiation rate for the case of purely timelike $k_{AF}^μ$, and find a radiation rate at high energies that has a contribution that does not depend on the mass used to regulate the photons.

preprint2016arXiv

Fast Large Scale Structure Perturbation Theory using 1D FFTs

The usual fluid equations describing the large-scale evolution of mass density in the universe can be written as local in the density, velocity divergence, and velocity potential fields. As a result, the perturbative expansion in small density fluctuations, usually written in terms of convolutions in Fourier space, can be written as a series of products of these fields evaluated at the same location in configuration space. Based on this, we establish a new method to numerically evaluate the 1-loop power spectrum (i.e., Fourier transform of the 2-point correlation function) with one-dimensional Fast Fourier Transforms. This is exact and a few orders of magnitude faster than previously used numerical approaches. Numerical results of the new method are in excellent agreement with the standard quadrature integration method. This fast model evaluation can in principle be extended to higher loop order where existing codes become painfully slow. Our approach follows by writing higher order corrections to the 2-point correlation function as, e.g., the correlation between two second-order fields or the correlation between a linear and a third-order field. These are then decomposed into products of correlations of linear fields and derivatives of linear fields. The method can also be viewed as evaluating three-dimensional Fourier space convolutions using products in configuration space, which may also be useful in other contexts where similar integrals appear.

preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in Fourier-space

We analyse the Baryon Acoustic Oscillation (BAO) signal of the final Baryon Oscillation Spectroscopic Survey (BOSS) data release (DR12). Our analysis is performed in Fourier-space, using the power spectrum monopole and quadrupole. The dataset includes $1\,198\,006$ galaxies over the redshift range $0.2 < z < 0.75$. We divide this dataset into three (overlapping) redshift bins with the effective redshifts $\zeff = 0.38$, $0.51$ and $0.61$. We demonstrate the reliability of our analysis pipeline using N-body simulations as well as $\sim 1000$ MultiDark-Patchy mock catalogues, which mimic the BOSS-DR12 target selection. We apply density field reconstruction to enhance the BAO signal-to-noise ratio. By including the power spectrum quadrupole we can separate the line-of-sight and angular modes, which allows us to constrain the angular diameter distance $D_A(z)$ and the Hubble parameter $H(z)$ separately. We obtain two independent $1.6\%$ and $1.5\%$ constraints on $D_A(z)$ and $2.9\%$ and $2.3\%$ constraints on $H(z)$ for the low ($\zeff=0.38$) and high ($\zeff=0.61$) redshift bin, respectively. We obtain two independent $1\%$ and $0.9\%$ constraints on the angular averaged distance $D_V(z)$, when ignoring the Alcock-Paczynski effect. The detection significance of the BAO signal is of the order of $8σ$ (post-reconstruction) for each of the three redshift bins. Our results are in good agreement with the Planck prediction within $Λ$CDM. 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~\citet{Alam2016} to produce the final cosmological constraints from BOSS.

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.

preprint2015arXiv

Building Unbiased Estimators from Non-Gaussian Likelihoods with Application to Shear Estimation

We develop a general framework for generating estimators of a given quantity which are unbiased to a given order in the difference between the true value of the underlying quantity and the fiducial position in theory space around which we expand the likelihood. We apply this formalism to rederive the optimal quadratic estimator and show how the replacement of the second derivative matrix with the Fisher matrix is a generic way of creating an unbiased estimator (assuming choice of the fiducial model is independent of data). Next apply the approach to estimation of shear lensing, closely following the work of Bernstein and Armstrong (2014). Our first order estimator reduces to their estimator in the limit of zero shear, but it also naturally allows for the case of non-constant shear and the easy calculation of correlation functions or power spectra using standard methods. Both our first-order estimator and Bernstein and Armstrong's estimator exhibit a bias which is quadratic in true shear. Our third-order estimator is, at least in the realm of the toy problem of Bernstein and Armstrong, unbiased to 0.1% in relative shear errors $Δg/|g|$ for shears up to $|g|=0.2$.

preprint2015arXiv

Singular Lorentz-Violating Lagrangians and Associated Finsler Structures

Several lagrangians associated to classical limits of lorenz-violating fermions in the Standard Model extension (SME) have been shown to yield Finsler functions when the theory is expressed in Euclidean space. When spin-couplings are present, the lagrangian can develop singularities that obstruct the construction of a globally defined Legendre transformation, leading to singular Finsler spaces. A specific sector of the SME where such problems arise is studied. It is found that the singular behavior can be eliminated by an appropriate lifting of the problem to an associated algebraic variety. This provides a smooth classical model for the singular problem. In Euclidean space, the procedure involves combining two related singular Finsler functions into a single smooth function with a semi-positive definite quadratic form defined on a desingularized variety.

preprint2014arXiv

DESI and other dark energy experiments in the era of neutrino mass measurements

We present Fisher matrix projections for future cosmological parameter measurements, including neutrino masses, dark energy, curvature, modified gravity, the inflationary perturbation spectrum, non-Gaussianity, and dark radiation. We focus on DESI and generally redshift surveys (BOSS, HETDEX, eBOSS, Euclid, and WFIRST), but also include CMB (Planck) and weak gravitational lensing (DES and LSST) constraints. The goal is to present a consistent set of projections, for concrete experiments, which are otherwise scattered throughout many papers and proposals. We include neutrino mass as a free parameter in most projections, as it will inevitably be relevant -- DESI and other experiments can measure the sum of neutrino masses to ~0.02 eV or better, while the minimum possible sum is ~0.06 eV. We note that the BAO-only use of galaxy clustering is substantially degraded as a dark energy probe in the presence of neutrino mass uncertainty -- using broadband galaxy power is critical, especially pushing it to as small a scale as possible, and big gains are achieved by combining lensing survey constraints with redshift survey constraints. We do not try to be especially innovative, e.g., in careful treatments of potential systematic errors -- these projections are intended as a straightforward baseline for comparison to more detailed analyses.

preprint2014arXiv

Gupta-Bleuler Photon Quantization in the SME

Photon quantization is implemented in the standard model extension (SME) using the Gupta-Bleuler method and BRST concepts. The quantization prescription applies to both the birefringent and non-birefringent CPT-even couplings. A curious incompatibility is found between the presence of the Lorentz-violating terms and the existence of a nontrivial conjugate momentum $Π^0$ yielding problems with covariant quantization procedure. Introduction of a mass regulator term can avoid the vanishing of $Π^0$ and allows for the implementation of a covariant quantization procedure. Field-theoretic calculations involving the SME photons can then be performed using the mass regulator, similar to the conventional procedure used in electrodynamics for infrared-divergence regulation.

preprint2014arXiv

Understanding higher-order nonlocal halo bias at large scales by combining the power spectrum with the bispectrum

Understanding the relation between underlying matter distribution and biased tracers such as galaxy or dark matter halo is essential to extract cosmological information from ongoing or future galaxy redshift surveys. At sufficiently large scales such as the BAO scale, a standard approach for the bias problem on the basis of the perturbation theory (PT) is to assume the `local bias' model in which the density field of biased tracers is deterministically expanded in terms of matter density field at the same position. The higher-order bias parameters are then determined by combining the power spectrum with higher-order statistics such as the bispectrum. As is pointed out by recent studies, however, nonlinear gravitational evolution naturally induces nonlocal bias terms even if initially starting only with purely local bias. As a matter of fact, previous works showed that the second-order nonlocal bias term, which corresponds to the gravitational tidal field, is important to explain the characteristic scale-dependence of the bispectrum. In this paper we extend the nonlocal bias term up to third order, and investigate whether the PT-based model including nonlocal bias terms can simultaneously explain the power spectrum and the bispectrum of simulated halos in $N$-body simulations. We show that the power spectrum, including density and momentum, and the bispectrum between halo and matter in $N$-body simulations can be simultaneously well explained by the model including up to third-order nonlocal bias terms up to k~0.1h/Mpc. Also, the results seem in a good agreement with theoretical predictions of a simple coevolution picture, although the agreement is not perfect. These demonstration clearly shows a failure of the local bias model even at such large scales, and we conclude that nonlocal bias terms should be consistently included in order to model statistics of halos. [abridged]

preprint2014arXiv

What the "simple renormalization group" approach to dark matter clustering really was

McDonald (2007) presented an approach to improving perturbation theory (PT) calculations of the dark matter power spectrum, with a derivation based on the idea of renormalization group flow with time. In spite of a questionable approximation made in deriving it, subsequent comparisons by several groups between the predictions of the resulting equation and N-body simulations showed remarkable improvement relative to "standard" PT (SPT) at similar order. In this brief note I show that the same final equation can be derived cleanly from the point of view not of flowing with time but with non-linear coupling strength, i.e., gradually dialing the coupling from the trivial value zero to the physical one. This understanding makes it clear how to extend the approach to higher order and other statistics. While I do not necessarily think this approach is best among the many, it may be interesting in that it contains a unique way of suppressing UV sensitivity. In passing I remind the reader of references demonstrating that SPT works remarkably well without improvement (except near redshift zero, where, fortunately for SPT, there is very little volume in the Universe).

preprint2013arXiv

Detection of Lyβauto-correlations and Lyα-Lyβcross-correlations in BOSS Data Release 9

The Lyman-$β$ forest refers to a region in the spectra of distant quasars that lies between the rest-frame Lyman-$β$ and Lyman-$γ$ emissions. The forest in this region is dominated by a combination of absorption due to resonant Ly$α$ and Ly$β$ scattering. When considering the 1D Ly$β$ forest in addition to the 1D Ly$α$ forest, the full statistical description of the data requires four 1D power spectra: Ly$α$ and Ly$β$ auto-power spectra and the Ly$α$-Ly$β$ real and imaginary cross-power spectra. We describe how these can be measured using an optimal quadratic estimator that naturally disentangles Ly$α$ and Ly$β$ contributions. Using a sample of approximately 60,000 quasar sight-lines from the BOSS Data Release 9, we make the measurement of the one-dimensional power spectrum of fluctuations due to the Ly$β$ resonant scattering. While we have not corrected our measurements for resolution damping of the power and other systematic effects carefully enough to use them for cosmological constraints, we can robustly conclude the following: i) Ly$β$ power spectrum and Ly$α$-Ly$β$ cross spectra are detected with high statistical significance; ii) the cross-correlation coefficient is $\approx 1$ on large scales; iii) the Ly$β$ measurements are contaminated by the associated OVI absorption, which is analogous to the SiIII contamination of the Ly$α$ forest. Measurements of the Ly$β$ forest will allow extension of the usable path-length for the Ly$α$ measurements while allowing a better understanding of the physics of intergalactic medium and thus more robust cosmological constraints.

preprint2013arXiv

Measurement of Baryon Acoustic Oscillations in the Lyman-alpha Forest Fluctuations in BOSS Data Release 9

We use the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) to detect and measure the position of the Baryonic Acoustic Oscillation (BAO) feature in the three-dimensional correlation function in the Lyman-alpha forest flux fluctuations at a redshift z=2.4. The feature is clearly detected at significance between 3 and 5 sigma (depending on the broadband model and method of error covariance matrix estimation) and is consistent with predictions of the standard LCDM model. We assess the biases in our method, stability of the error covariance matrix and possible systematic effects. We fit the resulting correlation function with several models that decouple the broadband and acoustic scale information. For an isotropic dilation factor, we measure 100x(alpha_iso-1) = -1.6 ^{+2.0+4.3+7.4}_{-2.0-4.1-6.8} (stat.) +/- 1.0 (syst.) (multiple statistical errors denote 1,2 and 3 sigma confidence limits) with respect to the acoustic scale in the fiducial cosmological model (flat LCDM with Omega_m=0.27, h=0.7). When fitting separately for the radial and transversal dilation factors we find marginalised constraints 100x(alpha_par-1) = -1.3 ^{+3.5+7.6 +12.3}_{-3.3-6.7-10.2} (stat.) +/- 2.0 (syst.) and 100x(alpha_perp-1) = -2.2 ^{+7.4+17}_{-7.1-15} +/- 3.0 (syst.). The dilation factor measurements are significantly correlated with cross-correlation coefficient of ~ -0.55. Errors become significantly non-Gaussian for deviations over 3 standard deviations from best fit value. Because of the data cuts and analysis method, these measurements give tighter constraints than a previous BAO analysis of the BOSS DR9 Lyman-alpha forest sample, providing an important consistency test of the standard cosmological model in a new redshift regime.

preprint2013arXiv

The Baryon Acoustic Oscillation Broadband and Broad-beam Array: Design Overview and Sensitivity Forecasts

This work describes a new instrument optimized for a detection of the neutral hydrogen 21cm power spectrum between redshifts of 0.5-1.5: the Baryon Acoustic Oscillation Broadband and Broad-beam (BAOBAB) Array. BAOBAB will build on the efforts of a first generation of 21cm experiments which are targeting a detection of the signal from the Epoch of Reionization at z ~ 10. At z ~ 1, the emission from neutral hydrogen in self-shielded overdense halos also presents an accessible signal, since the dominant, synchrotron foreground emission is considerably fainter than at redshift 10. The principle science driver for these observations are Baryon Acoustic Oscillations in the matter power spectrum which have the potential to act as a standard ruler and constrain the nature of dark energy. BAOBAB will fully correlate dual-polarization antenna tiles over the 600-900MHz band with a frequency resolution of 300 kHz and a system temperature of 50K. The number of antennas will grow in staged deployments, and reconfigurations of the array will allow for both traditional imaging and high power spectrum sensitivity operations. We present calculations of the power spectrum sensitivity for various array sizes, with a 35-element array measuring the cosmic neutral hydrogen fraction as a function of redshift, and a 132-element system detecting the BAO features in the power spectrum, yielding a 1.8% error on the z ~ 1 distance scale, and, in turn, significant improvements to constraints on the dark energy equation of state over an unprecedented range of redshifts from ~0.5-1.5.

preprint2013arXiv

The DESI Experiment, a whitepaper for Snowmass 2013

The Dark Energy Spectroscopic Instrument (DESI) is a massively multiplexed fiber-fed spectrograph that will make the next major advance in dark energy in the timeframe 2018-2022. On the Mayall telescope, DESI will obtain spectra and redshifts for at least 18 million emission-line galaxies, 4 million luminous red galaxies and 3 million quasi-stellar objects, in order to: probe the effects of dark energy on the expansion history using baryon acoustic oscillations (BAO), measure the gravitational growth history through redshift-space distortions, measure the sum of neutrino masses, and investigate the signatures of primordial inflation. The resulting 3-D galaxy maps at z<2 and Lyman-alpha forest at z>2 will make 1%-level measurements of the distance scale in 35 redshift bins, thus providing unprecedented constraints on cosmological models.

preprint2012arXiv

Classical Lagrangians for Momentum Dependent Lorentz Violation

Certain momentum-dependent terms in the fermion sector of the Lorentz-violating Standard Model Extension (SME) yield solvable classical lagrangians of a type not mentioned in the literature. These cases yield new relatively simple examples of Finsler and pseudo-Finsler structures. One of the cases involves antisymmetric $d$-type terms and yields a new example of a relatively simple covariant lagrangian.

preprint2012arXiv

Distribution function approach to redshift space distortions. Part II: N-body simulations

Measurement of redshift-space distortions (RSD) offers an attractive method to directly probe the cosmic growth history of density perturbations. A distribution function approach where RSD can be written as a sum over density weighted velocity moment correlators has recently been developed. We use Nbody simulations to investigate the individual contributions and convergence of this expansion for dark matter. If the series is expanded as a function of powers of mu, cosine of the angle between the Fourier mode and line of sight, there are a finite number of terms contributing at each order. We present these terms and investigate their contribution to the total as a function of wavevector k. For mu^2 the correlation between density and momentum dominates on large scales. Higher order corrections, which act as a Finger-of-God (FoG) term, contribute 1% at k~0.015h/Mpc, 10% at k~0.05h/Mpc at z=0, while for k>0.15h/Mpc they dominate and make the total negative. These higher order terms are dominated by density-energy density correlations which contribute negatively to the power, while the contribution from vorticity part of momentum density auto-correlation is an order of magnitude lower. For mu^4 term the dominant term on large scales is the scalar part of momentum density auto-correlation, while higher order terms dominate for k>0.15h/Mpc. For mu^6 and mu^8 we find it has very little power for k<0.15h/Mpc. We also compare the expansion to the full 2D P^ss(k,mu) as well as to their multipoles. For these statistics an infinite number of terms contribute and we find that the expansion achieves percent level accuracy for kmu<0.15h/Mpc at 6th order, but breaks down on smaller scales because the series is no longer perturbative. We explore resummation of the terms into FoG kernels, which extend the convergence up to a factor of 2 in scale. We find that the FoG kernels are approximately Lorentzian.

preprint2012arXiv

Distribution function approach to redshift space distortions. Part IV: perturbation theory applied to dark matter

We develop a perturbative approach to redshift space distortions (RSD) using the phase space distribution function approach and apply it to the dark matter redshift space power spectrum and its moments. RSD can be written as a sum over density weighted velocity moments correlators, with the lowest order being density, momentum density and stress energy density. We use standard and extended perturbation theory (PT) to determine their auto and cross correlators, comparing them to N-body simulations. We show which of the terms can be modeled well with the standard PT and which need additional terms that include higher order corrections which cannot be modeled in PT. Most of these additional terms are related to the small scale velocity dispersion effects, the so called finger of god (FoG) effects, which affect some, but not all, of the terms in this expansion, and which can be approximately modeled using a simple physically motivated ansatz such as the halo model. We point out that there are several velocity dispersions that enter into the detailed RSD analysis with very different amplitudes, which can be approximately predicted by the halo model. In contrast to previous models our approach systematically includes all of the terms at a given order in PT and provides a physical interpretation for the small scale dispersion values. We investigate RSD power spectrum as a function of μ, the cosine of the angle between the Fourier mode and line of sight, focusing on the lowest order powers of μand multipole moments which dominate the observable RSD power spectrum. Overall we find considerable success in modeling many, but not all, of the terms in this expansion.

preprint2012arXiv

Evidence for Quadratic Tidal Tensor Bias from the Halo Bispectrum

The relation between the clustering properties of luminous matter in the form of galaxies and the underlying dark matter distribution is of fundamental importance for the interpretation of ongoing and upcoming galaxy surveys. The so called local bias model, where galaxy density is a function of local matter density, is frequently discussed as a means to infer the matter power spectrum or correlation function from the measured galaxy correlation. However, gravitational evolution generates a term quadratic in the tidal tensor and thus non-local in the density field, even if this term is absent in the initial conditions (Lagrangian space). Because the term is quadratic, it contributes as a loop correction to the power spectrum, so the standard linear bias picture still applies on large scales, however, it contributes at leading order to the bispectrum for which it is significant on all scales. Such a term could also be present in Lagrangian space if halo formation were influenced by the tidal field. We measure the corresponding coupling strengths from the matter-matter-halo bispectrum in numerical simulations and find a non-vanishing coefficient for the tidal tensor term. We find no scale dependence of the bias parameters up to k=0.1 h/Mpc and that the tidal effect is increasing with halo mass. While the Lagrangian bias picture is a better description of our results than the Eulerian bias picture, our results suggest that there might be a tidal tensor bias already in the initial conditions. We also find that the coefficients of the quadratic density term deviate quite strongly from the theoretical predictions based on the spherical collapse model and a universal mass function. Both quadratic density and tidal tensor bias terms must be included in the modeling of galaxy clustering of current and future surveys if one wants to achieve the high precision cosmology promise of these datasets.

preprint2012arXiv

Non-Gaussian errors of baryonic acoustic oscillations

We revisit the uncertainty in baryon acoustic oscillation (BAO) forecasts and data analyses. In particular, we study how much the uncertainties on both the measured mean dilation scale and the associated error bar are affected by the non-Gaussianity of the non-linear density field. We examine two possible impacts of non-Gaussian analysis: (1) we derive the distance estimators from Gaussian theory, but use 1000 N-Body simulations to measure the actual errors, and compare this to the Gaussian prediction, and (2) we compute new optimal estimators, which requires the inverse of the non-Gaussian covariance matrix of the matter power spectrum. Obtaining an accurate and precise inversion is challenging, and we opted for a noise reduction technique applied on the covariance matrices. By measuring the bootstrap error on the inverted matrix, this work quantifies for the first time the significance of the non-Gaussian error corrections on the BAO dilation scale. We find that the variance (error squared) on distance measurements can deviate by up to 12% between both estimators, an effect that requires a large number of simulations to be resolved. We next apply a reconstruction algorithm to recover some of the BAO signal that had been smeared by non-linear evolution, and we rerun the analysis. We find that after reconstruction, the rms error on the distance measurement improves by a factor of ~1.7 at low redshift (consistent with previous results), and the variance (σ^2) shows a change of up to 18% between optimal and sub-optimal cases (note, however, that these discrepancies may depend in detail on the procedure used to isolate the BAO signal). We finally discuss the impact of this work on current data analyses.

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.

preprint2011arXiv

Distribution function approach to redshift space distortions

We develop a phase space distribution function approach to redshift space distortions (RSD), in which the redshift space density can be written as a sum over velocity moments of the distribution function. These moments are density weighted and their lowest orders are density, momentum density, and stress energy density. The series expansion is convergent on large scales. We perform an expansion of these velocity moments into helicity modes, which are eigenmodes under rotation around the axis of Fourier mode direction, generalizing the scalar, vector, tensor decomposition of perturbations to an arbitrary order. We show that only equal helicity moments correlate and derive the angular dependence of the individual contributions to the redshift space power spectrum in terms of angle mu between wave vector and line of sight. We show that the dominant term of mu^2 dependence on large scales is the cross-correlation between the density and scalar part of momentum density, which can be related to the time derivative of the matter power spectrum. Additional terms contributing and dominating on small scales are the vector part of momentum density-momentum density correlations, the energy density-density correlations, and the scalar part of anisotropic stress density-density correlations. Similarly, we identify 7 terms contributing to mu^4 dependence. Some of the advantages of the distribution function approach are that the series expansion converges on large scales and remains valid in multi-stream situations. We finish with a brief discussion of implications for RSD in galaxies relative to dark matter, highlighting the issue of scale dependent bias of velocity moments correlators.

preprint2011arXiv

Generating mock data sets for large-scale Lyman-α forest correlation measurements

Massive spectroscopic surveys of high-redshift quasars yield large numbers of correlated Lyman α absorption spectra that can be used to measure large-scale structure. Simulations of these surveys are required to accurately interpret the measurements of correlations and correct for systematic errors. An efficient method to generate mock realizations of Lyman α forest surveys is presented which generates a field over the lines of sight to the survey sources only, instead of having to generate it over the entire three-dimensional volume of the survey. The method can be calibrated to reproduce the power spectrum and one-point distribution function of the transmitted flux fraction, as well as the redshift evolution of these quantities, and is easily used for modeling any survey systematic effects. We present an example of how these mock surveys are applied to predict the measurement errors in a survey with similar parameters as the BOSS quasar survey in SDSS-III.

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 Lyman-alpha forest in three dimensions: measurements of large scale flux correlations from BOSS 1st-year data

Using a sample of approximately 14,000 z>2.1 quasars observed in the first year of the Baryon Oscillation Spectroscopic Survey (BOSS), we measure the three-dimensional correlation function of absorption in the Lyman-alpha forest. The angle-averaged correlation function of transmitted flux (F = exp(-tau)) is securely detected out to comoving separations of 60 Mpc/h, the first detection of flux correlations across widely separated sightlines. A quadrupole distortion of the redshift-space correlation function by peculiar velocities, the signature of the gravitational instability origin of structure in the Lyman-alpha forest, is also detected at high significance. We obtain a good fit to the data assuming linear theory redshift-space distortion and linear bias of the transmitted flux, relative to the matter fluctuations of a standard LCDM cosmological model (inflationary cold dark matter with a cosmological constant). At 95% confidence, we find a linear bias parameter 0.16<b<0.24 and redshift-distortion parameter 0.44<beta<1.20, at central redshift z=2.25, with a well constrained combination b(1+β)=0.336 +/- 0.012. The errors on beta are asymmetric, with beta=0 excluded at over 5 sigma confidence level. The value of beta is somewhat low compared to theoretical predictions, and our tests on synthetic data suggest that it is depressed (relative to expectations for the Lyman-alpha forest alone) by the presence of high column density systems and metal line absorption. These results set the stage for cosmological parameter determinations from three-dimensional structure in the Lyman-alpha forest, including anticipated constraints on dark energy from baryon acoustic oscillations.

preprint2010arXiv

Factoring the Dispersion Relation in the Presence of Lorentz Violation

We produce an explicit formula for the dispersion relation for the Dirac Equation in the Standard Model Extension (SME) in the presence of Lorentz violation. Our expression is obtained using a novel techniques which exploit the algebra of quaternions. The dispersion relation is found to conveniently factor in two special cases that each involve a mutually exclusive set of non-vanishing Lorentz-violating parameters. This suggests that a useful approach to studies of Lorentz-violating models is to split the parameter space into two separate pieces, each of which yields a simple, tractable dispersion relation that can be used for analysis.

preprint2010arXiv

Near term measurements with 21 cm intensity mapping: neutral hydrogen fraction and BAO at z<2

It is shown that 21 cm intensity mapping could be used in the near term to make cosmologically useful measurements. Large scale structure could be detected using existing radio telescopes, or using prototypes for dedicated redshift survey telescopes. This would provide a measure of the mean neutral hydrogen density, using redshift space distortions to break the degeneracy with the linear bias. We find that with only 200 hours of observing time on the Green Bank Telescope, the neutral hydrogen density could be measured to 25% precision at redshift 0.54<z<1.09. This compares favourably to current measurements, uses independent techniques, and would settle the controversy over an important parameter which impacts galaxy formation studies. In addition, a 4000 hour survey would allow for the detection of baryon acoustic oscillations, giving a cosmological distance measure at 3.5% precision. These observation time requirements could be greatly reduced with the construction of multiple pixel receivers. Similar results are possible using prototypes for dedicated cylindrical telescopes on month time scales, or SKA pathfinder aperture arrays on day time scales. Such measurements promise to improve our understanding of these quantities while beating a path for future generations of hydrogen surveys.

preprint2010arXiv

Projected Constraints on Modified Gravity Cosmologies from 21 cm Intensity Mapping

We present projected constraints on modified gravity models from the observational technique known as 21 cm intensity mapping, where cosmic structure is detected without resolving individual galaxies. The resulting map is sensitive to both BAO and weak lensing, two of the most powerful cosmological probes. It is found that a 200 m x 200 m cylindrical telescope, sensitive out to z=2.5, would be able to distinguish DGP from most dark energy models, and constrain the Hu & Sawicki f(R) model to |f_{R0}| < 9*10^(-6) at 95% confidence. The latter constraint makes extensive use of the lensing spectrum in the nonlinear regime. These results show that 21 cm intensity mapping is not only sensitive to modifications of the standard model's expansion history, but also to structure growth. This makes intensity mapping a powerful and economical technique, achievable on much shorter time scales than optical experiments that would probe the same era.

preprint2010arXiv

Vector Superfields and Lorentz Violation

We extend Lorentz-violating Supersymmetry models to include vector superfields. The CPT-preserving model generalizes easily, while the obvious attempt at generalizing the CPT-violating model meets serious obstructions. Generalizations of the CPT-preserving but Lorentz-Violating model to higher dimensions are also straightforward. Compactification is used to reduce the six-dimensional theory to an ${\cal{N}} =2$ Lorentz-violating theory in four dimensions, while the ten-dimensional theory is used to produce an ${\cal{N}}=4$ Lorentz-violating theory. This may be useful in future constructions involving Ads/CFT correspondence with Lorentz violation.

preprint2009arXiv

Gravitational redshift and other redshift-space distortions of the imaginary part of the power spectrum

I extend the usual linear-theory formula for large-scale clustering in redshift-space to include gravitational redshift. The extra contribution to the standard galaxy power spectrum is suppressed by k_c^{-2}, where k_c=c k/a H (k is the wavevector, a the expansion factor, and H=\dot{a}/a), and is thus effectively limited to the few largest-scale modes and very difficult to detect; however, a correlation, \propto k_c^{-1}, is generated between the real and imaginary parts of the Fourier space density fields of two different types of galaxy, which would otherwise be zero, i.e., the cross-power spectrum has an imaginary part: P_{ab}(k,μ)/P(k)=(b_a+fμ^2)(b_b+fμ^2) -i(3Ω_m/2)(μ/k_c)(b_a-b_b)+\mathcal{O}(k_c^{-2}), where P(k) is the real-space mass-density power spectrum, b_i are the galaxy biases, μis the cosine of the angle between the wavevector and line of sight, and f=dlnD/dlna (D is the linear growth factor). The total signal-to-noise of measurements of this effect is not dominated by the largest scales -- it converges at k~0.05 h/Mpc. This gravitational redshift result is pedagogically interesting, but naive in that it is gauge dependent and there are other effects of similar form and size, related to the transformation between observable and proper coordinates. I include these effects, which add other contributions to the coefficient of μ/k_c, and add a μ^3/k_c term, but don't qualitatively change the picture. The leading source of noise in the measurement is galaxy shot-noise, not sample variance, so developments that allow higher S/N surveys can make this measurement powerful, although it would otherwise be only marginally detectable in a JDEM-scale survey.

preprint2009arXiv

How to generate a significant effective temperature for cold dark matter, from first principles

I show how to reintroduce velocity dispersion into perturbation theory (PT) calculations of structure in the Universe, i.e., how to go beyond the pressureless fluid approximation, starting from first principles. This addresses a possible deficiency in uses of PT to compute clustering on the weakly non-linear scales that will be critical for probing dark energy. Specifically, I show how to derive a non-negligible value for the (initially tiny) velocity dispersion of dark matter particles, <δv^2>, where δv is the deviation of particle velocities from the local bulk flow. The calculation is essentially a renormalization of the homogeneous (zero order) dispersion by fluctuations 1st order in the initial power spectrum. For power law power spectra with n>-3, the small-scale fluctuations diverge and significant dispersion can be generated from an arbitrarily small starting value -- the dispersion level is set by an equilibrium between fluctuations generating more dispersion and dispersion suppressing fluctuations. For an n=-1.4 power law normalized to match the present non-linear scale, the dispersion would be ~100 km/s. This n corresponds roughly to the slope on the non-linear scale in the real \LambdaCDM Universe, but \LambdaCDM contains much less initial small-scale power -- not enough to bootstrap the small starting dispersion up to a significant value within linear theory (viewed very broadly, structure formation has actually taken place rather suddenly and recently, in spite of the usual "hierarchical" description). The next order PT calculation, which I carry out only at an order of magnitude level, should drive the dispersion up into balance with the growing structure, accounting for small dispersion effects seen recently in simulations.

preprint2008arXiv

How to measure redshift-space distortions without sample variance

We show how to use multiple tracers of large-scale density with different biases to measure the redshift-space distortion parameter beta=f/b=(dlnD/dlna)/b (where D is the growth rate and a the expansion factor), to a much better precision than one could achieve with a single tracer, to an arbitrary precision in the low noise limit. In combination with the power spectrum of the tracers this allows a much more precise measurement of the bias-free velocity divergence power spectrum, f^2 P_m - in fact, in the low noise limit f^2 P_m can be measured as well as would be possible if velocity divergence was observed directly, with rms improvement factor ~[5.2(beta^2+2 beta+2)/beta^2]^0.5 (e.g., ~10 times better than a single tracer for beta=0.4). This would allow a high precision determination of f D as a function of redshift with an error as low as 0.1%. We find up to two orders of magnitude improvement in Figure of Merit for the Dark Energy equation of state relative to Stage II, a factor of several better than other proposed Stage IV Dark Energy surveys. The ratio b_2/b_1 will be determined with an even greater precision than beta, producing, when measured as a function of scale, an exquisitely sensitive probe of the onset of non-linear bias. We also extend in more detail previous work on the use of the same technique to measure non-Gaussianity. Currently planned redshift surveys are typically designed with signal to noise of unity on scales of interest, and are not optimized for this technique. Our results suggest that this strategy may need to be revisited as there are large gains to be achieved from surveys with higher number densities of galaxies.

preprint2006arXiv

Clustering of dark matter tracers: renormalizing the bias parameters

A commonly used perturbative method for computing large-scale clustering of tracers of mass density, like galaxies, is to model the tracer density field as a Taylor series in the local smoothed mass density fluctuations, possibly adding a stochastic component. I suggest a set of parameter redefinitions, eliminating problematic perturbative correction terms, that should represent a modest improvement, at least, to this method. As presented here, my method can be used to compute the power spectrum and bispectrum to 4th order in initial density perturbations, and higher order extensions should be straightforward. While the model is technically unchanged at this order, just reparameterized, the renormalized model is more elegant, and should have better convergence behavior, for three reasons: First, in the usual approach the effects of beyond-linear-order bias parameters can be seen at asymptotically large scales, while after renormalization the linear model is preserved in the large-scale limit, i.e., the effects of higher order bias parameters are restricted to relatively high k. Second, while the standard approach includes smoothing to suppress large perturbative correction terms, resulting in dependence on the arbitrary cutoff scale, no cutoff-sensitive terms appear explicitly after my redefinitions (and, relatedly, my correction terms are less sensitive to high-k, non-linear, power). Third, the 3rd order bias parameter disappears entirely, so my model has one fewer free parameter than usual (this parameter was redundant at the order considered). This model predicts a small modification of the baryonic acoustic oscillation (BAO) signal, in real space, supporting the robustness of BAO as a probe of dark energy, and providing a complete perturbative description over the relevant range of scales.

preprint2004arXiv

Nonrelativisitic Ideal Gasses and Lorentz Violations

We develop statistical mechanics for a nonrelativisitic ideal gas in the presence of Lorentz violating background fields. The analysis is performed using the Standard-Model Extension (SME). We derive the corresponding laws of thermodynamics and find that, to lowest order in Lorentz violation, the scalar thermodynamic variables are corrected by a rotationally invariant combination of the Lorentz terms which can be interpreted in terms of a (frame dependent) effective mass. We find that spin couplings can induce a temperature independent polarization in the gas that is not present in the conventional case.