Showing 1–50 of 134 results for author: McDonald, P

Mitigation of DESI fiber assignment incompleteness effect on two-point clustering with small angular scale truncated estimators

Authors: M. Pinon, A. de Mattia, P. McDonald, E. Burtin, V. Ruhlmann-Kleider, M. White, D. Bianchi, A. J. Ross, J. Aguilar, S. Ahlen, D. Brooks, E. Chaussidon, T. Claybaugh, S. Cole, A. de la Macorra, B. Dey, P. Doel, K. Fanning, J. E. Forero-Romero, E. Gaztañaga, S. Gontcho A Gontcho, C. Howlett, D. Kirkby, T. Kisner, A. Kremin , et al. (27 additional authors not shown)

Abstract: We present a method to mitigate the effects of fiber assignment incompleteness in two-point power spectrum and correlation function measurements from galaxy spectroscopic surveys, by truncating small angular scales from estimators. We derive the corresponding modified correlation function and power spectrum windows to account for the small angular scale truncation in the theory prediction. We vali… ▽ More We present a method to mitigate the effects of fiber assignment incompleteness in two-point power spectrum and correlation function measurements from galaxy spectroscopic surveys, by truncating small angular scales from estimators. We derive the corresponding modified correlation function and power spectrum windows to account for the small angular scale truncation in the theory prediction. We validate this approach on simulations reproducing the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) with and without fiber assignment. We show that we recover unbiased cosmological constraints using small angular scale truncated estimators from simulations with fiber assignment incompleteness, with respect to standard estimators from complete simulations. Additionally, we present an approach to remove the sensitivity of the fits to high $k$ modes in the theoretical power spectrum, by applying a transformation to the data vector and window matrix. We find that our method efficiently mitigates the effect of fiber assignment incompleteness in two-point correlation function and power spectrum measurements, at low computational cost and with little statistical loss. △ Less

Submitted 7 June, 2024; originally announced June 2024.

Comments: 32 pages, 22 figures

The Construction of Large-scale Structure Catalogs for the Dark Energy Spectroscopic Instrument

Authors: A. J. Ross, J. Aguilar, S. Ahlen, S. Alam, A. Anand, S. Bailey, D. Bianchi, S. Brieden, D. Brooks, E. Burtin, A. Carnero Rosell, E. Chaussidon, T. Claybaugh, S. Cole, K. Dawson, A. de la Macorra, A. de Mattia, Arjun Dey, Biprateep Dey, P. Doel, K. Fanning, S. Ferraro, J. Ereza, A. Font-Ribera, J. E. Forero-Romero , et al. (59 additional authors not shown)

Abstract: We present the technical details on how large-scale structure (LSS) catalogs are constructed from redshifts measured from spectra observed by the Dark Energy Spectroscopic Instrument (DESI). The LSS catalogs provide the information needed to determine the relative number density of DESI tracers as a function of redshift and celestial coordinates and, e.g., determine clustering statistics. We produ… ▽ More We present the technical details on how large-scale structure (LSS) catalogs are constructed from redshifts measured from spectra observed by the Dark Energy Spectroscopic Instrument (DESI). The LSS catalogs provide the information needed to determine the relative number density of DESI tracers as a function of redshift and celestial coordinates and, e.g., determine clustering statistics. We produce catalogs that are weighted subsamples of the observed data, each matched to a weighted `random' catalog that forms an unclustered sampling of the probability density that DESI could have observed those data at each location. Precise knowledge of the DESI observing history and associated hardware performance allows for a determination of the DESI footprint and the number of times DESI has covered it at sub-arcsecond level precision. This enables the completeness of any DESI sample to be modeled at this same resolution. The pipeline developed to create LSS catalogs has been designed to easily allow robustness tests and enable future improvements. We describe how it allows ongoing work improving the match between galaxy and random catalogs, such as including further information when assigning redshifts to randoms, accounting for fluctuations in target density, accounting for variation in the redshift success rate, and accommodating blinding schemes. △ Less

Submitted 26 May, 2024; originally announced May 2024.

Comments: Supporting publication of DESI 2024II: Sample definitions, characteristics, and two-point clustering statistics

HOD-Dependent Systematics in Emission Line Galaxies for the DESI 2024 BAO analysis

Authors: C. Garcia-Quintero, J. Mena-Fernández, A. Rocher, S. Yuan, B. Hadzhiyska, O. Alves, M. Rashkovetskyi, H. Seo, N. Padmanabhan, S. Nadathur, C. Howlett, M. Ishak, L. Medina-Varela, P. McDonald, A. J. Ross, Y. Xie, X. Chen, A. Bera, J. Aguilar, S. Ahlen, U. Andrade, S. BenZvi, D. Brooks, E. Burtin, S. Chen , et al. (51 additional authors not shown)

Abstract: The Dark Energy Spectroscopic Instrument (DESI) will provide precise measurements of Baryon Acoustic Oscillations (BAO) to constrain the expansion history of the Universe and set stringent constraints on dark energy. Therefore, precise control of the global error budget due to various systematic effects is required for the DESI 2024 BAO analysis. In this work, we focus on the robustness of the BAO… ▽ More The Dark Energy Spectroscopic Instrument (DESI) will provide precise measurements of Baryon Acoustic Oscillations (BAO) to constrain the expansion history of the Universe and set stringent constraints on dark energy. Therefore, precise control of the global error budget due to various systematic effects is required for the DESI 2024 BAO analysis. In this work, we focus on the robustness of the BAO analysis against the Halo Occupation Distribution (HOD) modeling for the Emission Line Galaxy (ELG) tracer. Based on a common dark matter simulation, our analysis relies on HOD mocks tuned to early DESI data, namely the One-Percent survey data. To build the mocks, we use several HOD models for the ELG tracer as well as extensions to the baseline HOD models. Among these extensions, we consider distinct recipes for galactic conformity and assembly bias. We perform two independent analyses in the Fourier space and in the configuration space. We recover the BAO signal from two-point measurements after performing reconstruction on our mocks. Additionally, we also apply the control variates technique to reduce sample variance noise. Our BAO analysis can recover the isotropic BAO parameter $α_\text{iso}$ within 0.1\% and the Alcock Paczynski parameter $α_\text{AP}$ within 0.3\%. Overall, we find that our systematic error due to the HOD dependence is below 0.17\%, with the Fourier space analysis being more robust against the HOD systematics. We conclude that our analysis pipeline is robust enough against the HOD systematics for the ELG tracer in the DESI 2024 BAO analysis. △ Less

Submitted 12 April, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Comments: Supporting publication of DESI 2024 III: Baryon Acoustic Oscillations from Galaxies and Quasars

DESI 2024 VI: Cosmological Constraints from the Measurements of Baryon Acoustic Oscillations

Authors: DESI Collaboration, A. G. Adame, J. Aguilar, S. Ahlen, S. Alam, D. M. Alexander, M. Alvarez, O. Alves, A. Anand, U. Andrade, E. Armengaud, S. Avila, A. Aviles, H. Awan, B. Bahr-Kalus, S. Bailey, C. Baltay, A. Bault, J. Behera, S. BenZvi, A. Bera, F. Beutler, D. Bianchi, C. Blake, R. Blum , et al. (178 additional authors not shown)

Abstract: We present cosmological results from the measurement of baryon acoustic oscillations (BAO) in galaxy, quasar and Lyman-$α$ forest tracers from the first year of observations from the Dark Energy Spectroscopic Instrument (DESI), to be released in the DESI Data Release 1. DESI BAO provide robust measurements of the transverse comoving distance and Hubble rate, or their combination, relative to the s… ▽ More We present cosmological results from the measurement of baryon acoustic oscillations (BAO) in galaxy, quasar and Lyman-$α$ forest tracers from the first year of observations from the Dark Energy Spectroscopic Instrument (DESI), to be released in the DESI Data Release 1. DESI BAO provide robust measurements of the transverse comoving distance and Hubble rate, or their combination, relative to the sound horizon, in seven redshift bins from over 6 million extragalactic objects in the redshift range $0.1<z<4.2$. DESI BAO data alone are consistent with the standard flat $Λ$CDM cosmological model with a matter density $Ω_\mathrm{m}=0.295\pm 0.015$. Paired with a BBN prior and the robustly measured acoustic angular scale from the CMB, DESI requires $H_0=(68.52\pm0.62)$ km/s/Mpc. In conjunction with CMB anisotropies from Planck and CMB lensing data from Planck and ACT, we find $Ω_\mathrm{m}=0.307\pm 0.005$ and $H_0=(67.97\pm0.38)$ km/s/Mpc. Extending the baseline model with a constant dark energy equation of state parameter $w$, DESI BAO alone require $w=-0.99^{+0.15}_{-0.13}$. In models with a time-varying dark energy equation of state parametrized by $w_0$ and $w_a$, combinations of DESI with CMB or with SN~Ia individually prefer $w_0>-1$ and $w_a<0$. This preference is 2.6$σ$ for the DESI+CMB combination, and persists or grows when SN~Ia are added in, giving results discrepant with the $Λ$CDM model at the $2.5σ$, $3.5σ$ or $3.9σ$ levels for the addition of Pantheon+, Union3, or DES-SN5YR datasets respectively. For the flat $Λ$CDM model with the sum of neutrino mass $\sum m_ν$ free, combining the DESI and CMB data yields an upper limit $\sum m_ν< 0.072$ $(0.113)$ eV at 95% confidence for a $\sum m_ν>0$ $(\sum m_ν>0.059)$ eV prior. These neutrino-mass constraints are substantially relaxed in models beyond $Λ$CDM. [Abridged.] △ Less

Submitted 24 April, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Comments: This DESI Collaboration Key Publication is part of the 2024 publication series using the first year of observations (see https://data.desi.lbl.gov/doc/papers). Typos corrected and a new figure and discussion added to Appendix A

DESI 2024 IV: Baryon Acoustic Oscillations from the Lyman Alpha Forest

Authors: DESI Collaboration, A. G. Adame, J. Aguilar, S. Ahlen, S. Alam, D. M. Alexander, M. Alvarez, O. Alves, A. Anand, U. Andrade, E. Armengaud, S. Avila, A. Aviles, H. Awan, S. Bailey, C. Baltay, A. Bault, J. Bautista, J. Behera, S. BenZvi, F. Beutler, D. Bianchi, C. Blake, R. Blum, S. Brieden , et al. (174 additional authors not shown)

Abstract: We present the measurement of Baryon Acoustic Oscillations (BAO) from the Lyman-$α$ (Ly$α$) forest of high-redshift quasars with the first-year dataset of the Dark Energy Spectroscopic Instrument (DESI). Our analysis uses over $420\,000$ Ly$α$ forest spectra and their correlation with the spatial distribution of more than $700\,000$ quasars. An essential facet of this work is the development of a… ▽ More We present the measurement of Baryon Acoustic Oscillations (BAO) from the Lyman-$α$ (Ly$α$) forest of high-redshift quasars with the first-year dataset of the Dark Energy Spectroscopic Instrument (DESI). Our analysis uses over $420\,000$ Ly$α$ forest spectra and their correlation with the spatial distribution of more than $700\,000$ quasars. An essential facet of this work is the development of a new analysis methodology on a blinded dataset. We conducted rigorous tests using synthetic data to ensure the reliability of our methodology and findings before unblinding. Additionally, we conducted multiple data splits to assess the consistency of the results and scrutinized various analysis approaches to confirm their robustness. For a given value of the sound horizon ($r_d$), we measure the expansion at $z_{\rm eff}=2.33$ with 2\% precision, $H(z_{\rm eff}) = (239.2 \pm 4.8) (147.09~{\rm Mpc} /r_d)$ km/s/Mpc. Similarly, we present a 2.4\% measurement of the transverse comoving distance to the same redshift, $D_M(z_{\rm eff}) = (5.84 \pm 0.14) (r_d/147.09~{\rm Mpc})$ Gpc. Together with other DESI BAO measurements at lower redshifts, these results are used in a companion paper to constrain cosmological parameters. △ Less

Submitted 12 April, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Comments: This DESI Collaboration Key Publication is part of the 2024 publication series using the first year of observations (see https://data.desi.lbl.gov/doc/papers)

DESI 2024 III: Baryon Acoustic Oscillations from Galaxies and Quasars

Authors: DESI Collaboration, A. G. Adame, J. Aguilar, S. Ahlen, S. Alam, D. M. Alexander, M. Alvarez, O. Alves, A. Anand, U. Andrade, E. Armengaud, S. Avila, A. Aviles, H. Awan, S. Bailey, C. Baltay, A. Bault, J. Behera, S. BenZvi, F. Beutler, D. Bianchi, C. Blake, R. Blum, S. Brieden, A. Brodzeller , et al. (171 additional authors not shown)

Abstract: We present the DESI 2024 galaxy and quasar baryon acoustic oscillations (BAO) measurements using over 5.7 million unique galaxy and quasar redshifts in the range 0.1<z<2.1. Divided by tracer type, we utilize 300,017 galaxies from the magnitude-limited Bright Galaxy Survey with 0.1<z<0.4, 2,138,600 Luminous Red Galaxies with 0.4<z<1.1, 2,432,022 Emission Line Galaxies with 0.8<z<1.6, and 856,652 qu… ▽ More We present the DESI 2024 galaxy and quasar baryon acoustic oscillations (BAO) measurements using over 5.7 million unique galaxy and quasar redshifts in the range 0.1<z<2.1. Divided by tracer type, we utilize 300,017 galaxies from the magnitude-limited Bright Galaxy Survey with 0.1<z<0.4, 2,138,600 Luminous Red Galaxies with 0.4<z<1.1, 2,432,022 Emission Line Galaxies with 0.8<z<1.6, and 856,652 quasars with 0.8<z<2.1, over a ~7,500 square degree footprint. The analysis was blinded at the catalog-level to avoid confirmation bias. All fiducial choices of the BAO fitting and reconstruction methodology, as well as the size of the systematic errors, were determined on the basis of the tests with mock catalogs and the blinded data catalogs. We present several improvements to the BAO analysis pipeline, including enhancing the BAO fitting and reconstruction methods in a more physically-motivated direction, and also present results using combinations of tracers. We present a re-analysis of SDSS BOSS and eBOSS results applying the improved DESI methodology and find scatter consistent with the level of the quoted SDSS theoretical systematic uncertainties. With the total effective survey volume of ~ 18 Gpc$^3$, the combined precision of the BAO measurements across the six different redshift bins is ~0.52%, marking a 1.2-fold improvement over the previous state-of-the-art results using only first-year data. We detect the BAO in all of these six redshift bins. The highest significance of BAO detection is $9.1σ$ at the effective redshift of 0.93, with a constraint of 0.86% placed on the BAO scale. We find our measurements are systematically larger than the prediction of Planck-2018 LCDM model at z<0.8. We translate the results into transverse comoving distance and radial Hubble distance measurements, which are used to constrain cosmological models in our companion paper [abridged]. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: This DESI Collaboration Key Publication is part of the 2024 publication series using the first year of observations (see https://data.desi.lbl.gov/doc/papers)

The 3D Lyman-$α$ Forest Power Spectrum from eBOSS DR16

Authors: Roger de Belsunce, Oliver H. E. Philcox, Vid Irsic, Pat McDonald, Julien Guy, Nathalie Palanque-Delabrouille

Abstract: We measure the three-dimensional power spectrum (P3D) of the transmitted flux in the Lyman-$α$ (Ly-$α$) forest using the complete extended Baryon Oscillation Spectroscopic Survey data release 16 (eBOSS DR16). This sample consists of 205,012 quasar spectra in the redshift range 2 <= z <= 4 at an effective redshift z=2.334. We propose a pair-count spectral estimator in configuration space, weighting… ▽ More We measure the three-dimensional power spectrum (P3D) of the transmitted flux in the Lyman-$α$ (Ly-$α$) forest using the complete extended Baryon Oscillation Spectroscopic Survey data release 16 (eBOSS DR16). This sample consists of 205,012 quasar spectra in the redshift range 2 <= z <= 4 at an effective redshift z=2.334. We propose a pair-count spectral estimator in configuration space, weighting each pair by exp(ikr), for wave vector k and pixel pair separation r, effectively measuring the anisotropic power spectrum without the need for fast Fourier transforms. This accounts for the window matrix in a tractable way, avoiding artifacts found in Fourier-transform based power spectrum estimators due to the sparse sampling transverse to the line-of-sight of Ly-$α$ skewers. We extensively test our pipeline on two sets of mocks: (i) idealized Gaussian random fields with a sparse sampling of Ly-$α$ skewers, and (ii) log-normal LyaCoLoRe mocks including realistic noise levels, the eBOSS survey geometry and contaminants. On eBOSS DR16 data, the Kaiser formula with a non-linear correction term obtained from hydrodynamic simulations yields a good fit to the power spectrum data in the range 0.02 <= k <= 0.35 h/Mpc at the 1-2 sigma level with a covariance matrix derived from LyaCoLoRe mocks. We demonstrate a promising new approach for full-shape cosmological analyses of Ly-$α$ forest data from cosmological surveys such as eBOSS, the currently observing Dark Energy Spectroscopic Instrument and future surveys such as the Prime Focus Spectrograph, WEAVE-QSO and 4MOST. △ Less

Submitted 13 March, 2024; originally announced March 2024.

Comments: 14 pages, 11 figures

Baryon Acoustic Oscillation Theory and Modelling Systematics for the DESI 2024 results

Authors: Shi-Fan Chen, Cullan Howlett, Martin White, Patrick McDonald, Ashley J. Ross, Hee-Jong Seo, Nikhil Padmanabhan, J. Aguilar, S. Ahlen, S. Alam, O. Alves, R. Blum, D. Brooks, X. Chen, S. Cole, T. M. Davis, K. Dawson, A. de la Macorra, Arjun Dey, Z. Ding, P. Doel, S. Ferraro, A. Font-Ribera, D. Forero-Sánchez, J. E. Forero-Romero , et al. (33 additional authors not shown)

Abstract: This paper provides a comprehensive overview of how fitting of Baryon Acoustic Oscillations (BAO) is carried out within the upcoming Dark Energy Spectroscopic Instrument's (DESI) 2024 results using its DR1 dataset, and the associated systematic error budget from theory and modelling of the BAO. We derive new results showing how non-linearities in the clustering of galaxies can cause potential bias… ▽ More This paper provides a comprehensive overview of how fitting of Baryon Acoustic Oscillations (BAO) is carried out within the upcoming Dark Energy Spectroscopic Instrument's (DESI) 2024 results using its DR1 dataset, and the associated systematic error budget from theory and modelling of the BAO. We derive new results showing how non-linearities in the clustering of galaxies can cause potential biases in measurements of the isotropic ($α_{\mathrm{iso}}$) and anisotropic ($α_{\mathrm{ap}}$) BAO distance scales, and how these can be effectively removed with an appropriate choice of reconstruction algorithm. We then demonstrate how theory leads to a clear choice for how to model the BAO and develop, implement and validate a new model for the remaining smooth-broadband (i.e., without BAO) component of the galaxy clustering. Finally, we explore the impact of all remaining modelling choices on the BAO constraints from DESI using a suite of high-precision simulations, arriving at a set of best-practices for DESI BAO fits, and an associated theory and modelling systematic error. Overall, our results demonstrate the remarkable robustness of the BAO to all our modelling choices and motivate a combined theory and modelling systematic error contribution to the post-reconstruction DESI BAO measurements of no more than $0.1\%$ ($0.2\%$) for its isotropic (anisotropic) distance measurements. We expect the theory and best-practices laid out to here to be applicable to other BAO experiments in the era of DESI and beyond. △ Less

Submitted 21 February, 2024; originally announced February 2024.

Comments: 29 pages, 18 figures, 1 table, submitted to MNRAS

Homological properties of the relative Frobenius morphism

Authors: Peter M. McDonald

Abstract: This work concerns maps of commutative noetherian local rings containing a field of positive characteristic. Given such a map $\varphi$ of finite flat dimension, the results relate homological properties of the relative Frobenius of $\varphi$ to those of the fibers of $\varphi$. The focus is on the complete intersection property and the Gorenstein property. This work concerns maps of commutative noetherian local rings containing a field of positive characteristic. Given such a map $\varphi$ of finite flat dimension, the results relate homological properties of the relative Frobenius of $\varphi$ to those of the fibers of $\varphi$. The focus is on the complete intersection property and the Gorenstein property. △ Less

Submitted 3 January, 2024; originally announced January 2024.

Comments: 14 pages, comments welcome!

MSC Class: 13A35; 13D05; 13B10; 13H10

Mitigating the noise of DESI mocks using analytic control variates

Authors: Boryana Hadzhiyska, Martin J. White, Xinyi Chen, Lehman H. Garrison, Joseph DeRose, Nikhil Padmanabhan, Cristhian Garcia-Quintero, Juan Mena-Fernández, Shi-Fan Chen, Hee-Jong Seo, Patrick McDonald, Jessica Aguilar, Steven Ahlen, David Brooks, Todd Claybaugh, Axel de la Macorra, Peter Doel, Andreu Font-Ribera, Jaime E. Forero-Romero, Satya Gontcho A Gontcho, Klaus Honscheid, Anthony Kremin, Martin Landriau, Marc Manera, Ramon Miquel , et al. (8 additional authors not shown)

Abstract: In order to address fundamental questions related to the expansion history of the Universe and its primordial nature with the next generation of galaxy experiments, we need to model reliably large-scale structure observables such as the correlation function and the power spectrum. Cosmological $N$-body simulations provide a reference through which we can test our models, but their output suffers f… ▽ More In order to address fundamental questions related to the expansion history of the Universe and its primordial nature with the next generation of galaxy experiments, we need to model reliably large-scale structure observables such as the correlation function and the power spectrum. Cosmological $N$-body simulations provide a reference through which we can test our models, but their output suffers from sample variance on large scales. Fortunately, this is the regime where accurate analytic approximations exist. To reduce the variance, which is key to making optimal use of these simulations, we can leverage the accuracy and precision of such analytic descriptions using Control Variates (CV). The power of control variates stems from utilizing inexpensive but highly correlated surrogates of the statistics one wishes to measure. The stronger the correlation between the surrogate and the statistic of interest, the larger the variance reduction delivered by the method. We apply two control variate formulations to mock catalogs generated in anticipation of upcoming data from the Dark Energy Spectroscopic Instrument (DESI) to test the robustness of its analysis pipeline. Our CV-reduced measurements offer a factor of 5-10 improvement in the measurement error compared with the raw measurements. We explore the relevant properties of the galaxy samples that dictate this reduction and comment on the improvements we find on some of the derived quantities relevant to Baryon Acoustic Oscillation (BAO) analysis. We also provide an optimized package for computing the power spectra and other two-point statistics of an arbitrary galaxy catalog as well as a pipeline for obtaining CV-reduced measurements on any of the AbacusSummit cubic box outputs. We make our scripts publicly available and report a speed improvement of $\sim$10 for a grid size of $N_{\rm mesh} = 256^3$ compared with \texttt{nbodykit}. △ Less

Submitted 16 October, 2023; v1 submitted 23 August, 2023; originally announced August 2023.

Comments: 17 pages, 9 figures, public package (for power spectrum and control variates estimation)

Multiplier ideals and klt singularities via (derived) splittings

Authors: Peter M. McDonald

Abstract: Let $X$ be a normal, excellent, noetherian scheme over $\operatorname{Spec}\mathbb{Q}$ with a dualizing complex. In this note, we find an alternate characterization of the multiplier ideal of $X$, as defined by de Fernex-Hacon, by considering maps $π_*ω_Y\to\mathcal{O}_X$ where $π:Y\to X$ ranges over all regular alterations. As a corollary to this result, we give a derived splinter characterizatio… ▽ More Let $X$ be a normal, excellent, noetherian scheme over $\operatorname{Spec}\mathbb{Q}$ with a dualizing complex. In this note, we find an alternate characterization of the multiplier ideal of $X$, as defined by de Fernex-Hacon, by considering maps $π_*ω_Y\to\mathcal{O}_X$ where $π:Y\to X$ ranges over all regular alterations. As a corollary to this result, we give a derived splinter characterization of klt singularities, akin to the characterization of rational singularities given by Kovács and Bhatt. We also give an analogous description of the test ideal in characteristic $p>2$ as a corollary to a result of Epstein-Schwede. △ Less

Submitted 22 February, 2024; v1 submitted 15 July, 2023; originally announced July 2023.

Comments: 15 pages, comments welcome! Result generalized to excellent schemes with dualizing complexes

MSC Class: 14F18; 14B05; 14E15; 13A35

The Early Data Release of the Dark Energy Spectroscopic Instrument

Authors: DESI Collaboration, A. G. Adame, J. Aguilar, S. Ahlen, S. Alam, G. Aldering, D. M. Alexander, R. Alfarsy, C. Allende Prieto, M. Alvarez, O. Alves, A. Anand, F. Andrade-Oliveira, E. Armengaud, J. Asorey, S. Avila, A. Aviles, S. Bailey, A. Balaguera-Antolínez, O. Ballester, C. Baltay, A. Bault, J. Bautista, J. Behera, S. F. Beltran , et al. (240 additional authors not shown)

Abstract: The Dark Energy Spectroscopic Instrument (DESI) completed its five-month Survey Validation in May 2021. Spectra of stellar and extragalactic targets from Survey Validation constitute the first major data sample from the DESI survey. This paper describes the public release of those spectra, the catalogs of derived properties, and the intermediate data products. In total, the public release includes… ▽ More The Dark Energy Spectroscopic Instrument (DESI) completed its five-month Survey Validation in May 2021. Spectra of stellar and extragalactic targets from Survey Validation constitute the first major data sample from the DESI survey. This paper describes the public release of those spectra, the catalogs of derived properties, and the intermediate data products. In total, the public release includes good-quality spectral information from 466,447 objects targeted as part of the Milky Way Survey, 428,758 as part of the Bright Galaxy Survey, 227,318 as part of the Luminous Red Galaxy sample, 437,664 as part of the Emission Line Galaxy sample, and 76,079 as part of the Quasar sample. In addition, the release includes spectral information from 137,148 objects that expand the scope beyond the primary samples as part of a series of secondary programs. Here, we describe the spectral data, data quality, data products, Large-Scale Structure science catalogs, access to the data, and references that provide relevant background to using these spectra. △ Less

Submitted 15 June, 2023; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 43 pages, 7 figures, 17 tables, submitted to AJ, DESI EDR references added

Validation of the Scientific Program for the Dark Energy Spectroscopic Instrument

Authors: DESI Collaboration, A. G. Adame, J. Aguilar, S. Ahlen, S. Alam, G. Aldering, D. M. Alexander, R. Alfarsy, C. Allende Prieto, M. Alvarez, O. Alves, A. Anand, F. Andrade-Oliveira, E. Armengaud, J. Asorey, S. Avila, A. Aviles, S. Bailey, A. Balaguera-Antolínez, O. Ballester, C. Baltay, A. Bault, J. Bautista, J. Behera, S. F. Beltran , et al. (239 additional authors not shown)

Abstract: The Dark Energy Spectroscopic Instrument (DESI) was designed to conduct a survey covering 14,000 deg$^2$ over five years to constrain the cosmic expansion history through precise measurements of Baryon Acoustic Oscillations (BAO). The scientific program for DESI was evaluated during a five month Survey Validation (SV) campaign before beginning full operations. This program produced deep spectra of… ▽ More The Dark Energy Spectroscopic Instrument (DESI) was designed to conduct a survey covering 14,000 deg$^2$ over five years to constrain the cosmic expansion history through precise measurements of Baryon Acoustic Oscillations (BAO). The scientific program for DESI was evaluated during a five month Survey Validation (SV) campaign before beginning full operations. This program produced deep spectra of tens of thousands of objects from each of the stellar (MWS), bright galaxy (BGS), luminous red galaxy (LRG), emission line galaxy (ELG), and quasar target classes. These SV spectra were used to optimize redshift distributions, characterize exposure times, determine calibration procedures, and assess observational overheads for the five-year program. In this paper, we present the final target selection algorithms, redshift distributions, and projected cosmology constraints resulting from those studies. We also present a `One-Percent survey' conducted at the conclusion of Survey Validation covering 140 deg$^2$ using the final target selection algorithms with exposures of a depth typical of the main survey. The Survey Validation indicates that DESI will be able to complete the full 14,000 deg$^2$ program with spectroscopically-confirmed targets from the MWS, BGS, LRG, ELG, and quasar programs with total sample sizes of 7.2, 13.8, 7.46, 15.7, and 2.87 million, respectively. These samples will allow exploration of the Milky Way halo, clustering on all scales, and BAO measurements with a statistical precision of 0.28% over the redshift interval $z<1.1$, 0.39% over the redshift interval $1.1<z<1.9$, and 0.46% over the redshift interval $1.9<z<3.5$. △ Less

Submitted 12 January, 2024; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 42 pages, 18 figures, accepted by AJ

First Detection of the BAO Signal from Early DESI Data

Authors: Jeongin Moon, David Valcin, Michael Rashkovetskyi, Christoph Saulder, Jessica Nicole Aguilar, Steven Ahlen, Shadab Alam, Stephen Bailey, Charles Baltay, Robert Blum, David Brooks, Etienne Burtin, Edmond Chaussidon, Kyle Dawson, Axel de la Macorra, Arnaud de Mattia, Govinda Dhungana, Daniel Eisenstein, Brenna Flaugher, Andreu Font-Ribera, Jaime E. Forero-Romero, Cristhian Garcia-Quintero, Satya Gontcho A Gontcho, Julien Guy, Malik Muhammad Sikandar Hanif , et al. (43 additional authors not shown)

Abstract: We present the first detection of the baryon acoustic oscillations (BAO) signal obtained using unblinded data collected during the initial two months of operations of the Stage-IV ground-based Dark Energy Spectroscopic Instrument (DESI). From a selected sample of 261,291 Luminous Red Galaxies spanning the redshift interval 0.4 < z < 1.1 and covering 1651 square degrees with a 57.9% completeness le… ▽ More We present the first detection of the baryon acoustic oscillations (BAO) signal obtained using unblinded data collected during the initial two months of operations of the Stage-IV ground-based Dark Energy Spectroscopic Instrument (DESI). From a selected sample of 261,291 Luminous Red Galaxies spanning the redshift interval 0.4 < z < 1.1 and covering 1651 square degrees with a 57.9% completeness level, we report a ~5 sigma level BAO detection and the measurement of the BAO location at a precision of 1.7%. Using a Bright Galaxy Sample of 109,523 galaxies in the redshift range 0.1 < z < 0.5, over 3677 square degrees with a 50.0% completeness, we also detect the BAO feature at ~3 sigma significance with a 2.6% precision. These first BAO measurements represent an important milestone, acting as a quality control on the optimal performance of the complex robotically-actuated, fiber-fed DESI spectrograph, as well as an early validation of the DESI spectroscopic pipeline and data management system. Based on these first promising results, we forecast that DESI is on target to achieve a high-significance BAO detection at sub-percent precision with the completed 5-year survey data, meeting the top-level science requirements on BAO measurements. This exquisite level of precision will set new standards in cosmology and confirm DESI as the most competitive BAO experiment for the remainder of this decade. △ Less

Submitted 19 October, 2023; v1 submitted 17 April, 2023; originally announced April 2023.

Comments: 17 pages, 10 figures, 4 tables. MNRAS, 525, 5406. Author accepted manuscript matching the published version

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 525, Issue 4, November 2023, Pages 5406-5422

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

Authors: David J. Schlegel, Juna A. Kollmeier, Greg Aldering, Stephen Bailey, Charles Baltay, Christopher Bebek, Segev BenZvi, Robert Besuner, Guillermo Blanc, Adam S. Bolton, Ana Bonaca, Mohamed Bouri, David Brooks, Elizabeth Buckley-Geer, Zheng Cai, Jeffrey Crane, Regina Demina, Joseph DeRose, Arjun Dey, Peter Doel, Xiaohui Fan, Simone Ferraro, Douglas Finkbeiner, Andreu Font-Ribera, Satya Gontcho A Gontcho , et al. (64 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 9 September, 2022; originally announced September 2022.

Comments: Contributed White Paper to Snowmass 2021. arXiv admin note: substantial text overlap with arXiv:1907.11171. text overlap with arXiv:2209.03585

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

Authors: David J. Schlegel, Simone Ferraro, Greg Aldering, Charles Baltay, Segev BenZvi, Robert Besuner, Guillermo A. Blanc, Adam S. Bolton, Ana Bonaca, David Brooks, Elizabeth Buckley-Geer, Zheng Cai, Joseph DeRose, Arjun Dey, Peter Doel, Alex Drlica-Wagner, Xiaohui Fan, Gaston Gutierrez, Daniel Green, Julien Guy, Dragan Huterer, Leopoldo Infante, Patrick Jelinsky, Dionysios Karagiannis, Stephen M. Kent , et al. (40 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 8 September, 2022; originally announced September 2022.

Comments: Contribution to Snowmass 2021

Overview of the Instrumentation for the Dark Energy Spectroscopic Instrument

Authors: B. Abareshi, J. Aguilar, S. Ahlen, Shadab Alam, David M. Alexander, R. Alfarsy, L. Allen, C. Allende Prieto, O. Alves, J. Ameel, E. Armengaud, J. Asorey, Alejandro Aviles, S. Bailey, A. Balaguera-Antolínez, O. Ballester, C. Baltay, A. Bault, S. F. Beltran, B. Benavides, S. BenZvi, A. Berti, R. Besuner, Florian Beutler, D. Bianchi , et al. (242 additional authors not shown)

Abstract: The Dark Energy Spectroscopic Instrument (DESI) has embarked on an ambitious five-year survey to explore the nature of dark energy with spectroscopy of 40 million galaxies and quasars. DESI will determine precise redshifts and employ the Baryon Acoustic Oscillation method to measure distances from the nearby universe to z > 3.5, as well as measure the growth of structure and probe potential modifi… ▽ More The Dark Energy Spectroscopic Instrument (DESI) has embarked on an ambitious five-year survey to explore the nature of dark energy with spectroscopy of 40 million galaxies and quasars. DESI will determine precise redshifts and employ the Baryon Acoustic Oscillation method to measure distances from the nearby universe to z > 3.5, as well as measure the growth of structure and probe potential modifications to general relativity. In this paper we describe the significant instrumentation we developed for the DESI survey. The new instrumentation includes a wide-field, 3.2-deg diameter prime-focus corrector that focuses the light onto 5020 robotic fiber positioners on the 0.812 m diameter, aspheric focal surface. The positioners and their fibers are divided among ten wedge-shaped petals. Each petal is connected to one of ten spectrographs via a contiguous, high-efficiency, nearly 50 m fiber cable bundle. The ten spectrographs each use a pair of dichroics to split the light into three channels that together record the light from 360 - 980 nm with a resolution of 2000 to 5000. We describe the science requirements, technical requirements on the instrumentation, and management of the project. DESI was installed at the 4-m Mayall telescope at Kitt Peak, and we also describe the facility upgrades to prepare for DESI and the installation and functional verification process. DESI has achieved all of its performance goals, and the DESI survey began in May 2021. Some performance highlights include RMS positioner accuracy better than 0.1", SNR per \sqrtÅ > 0.5 for a z > 2 quasar with flux 0.28e-17 erg/s/cm^2/A at 380 nm in 4000s, and median SNR = 7 of the [OII] doublet at 8e-17 erg/s/cm^2 in a 1000s exposure for emission line galaxies at z = 1.4 - 1.6. We conclude with highlights from the on-sky validation and commissioning of the instrument, key successes, and lessons learned. (abridged) △ Less

Submitted 22 May, 2022; originally announced May 2022.

Comments: 78 pages, 32 figures, submitted to AJ

Nuclear spin relaxation in aqueous paramagnetic ion solutions

Authors: David A. Faux, Örs Istók, Arifah A. Rahaman, Peter J. McDonald, Eoin McKiernan, Dermot F. Brougham

Abstract: An angular time-dependent probability density function describing Brownian or anomalous rotational dynamics of fixed-length atom-to-atom vectors is presented. The probability density function, which fully incorporates angular boundary conditions, is applied to aqueous ion complexes. The rotational dynamics of ion-$^1$H vectors are shown by molecular dynamics (MD) simulation to be Brownian. A Brown… ▽ More An angular time-dependent probability density function describing Brownian or anomalous rotational dynamics of fixed-length atom-to-atom vectors is presented. The probability density function, which fully incorporates angular boundary conditions, is applied to aqueous ion complexes. The rotational dynamics of ion-$^1$H vectors are shown by molecular dynamics (MD) simulation to be Brownian. A Brownian shell model is presented which yields a closed form expression for the frequency-dependent nuclear-magnetic-resonance spin-lattice relaxation rate $T_1^{-1}(ω)$ based on a distance parameter and time constant. Appropriate combinations of shell and/or continuum models are shown to provide excellent fully-quantitative fits to experimental $T_1^{-1}(ω)$ dispersion curves from aqueous manganese(II), iron(III) and copper(II) chloride solutions. The distance parameters and time constants obtained from the fits are in good agreement with independent experimental and MD data in the literature. The Brownian shell model is a significant enhancement to existing particle-particle models that describe the rotational correlation function as a single exponential and are unable to provide the correct distance dependence for a shell of $^1$H spin density preventing a match to experiment without an arbitrary scaling factor. △ Less

Submitted 12 September, 2021; originally announced September 2021.

Comments: 39 pages, 7 figures, includes an appendix

Unified galaxy power spectrum measurements from 6dFGS, BOSS, and eBOSS

Authors: Florian Beutler, Patrick McDonald

Abstract: We make use of recent developments in the analysis of galaxy redshift surveys to present an easy to use matrix-based analysis framework for the galaxy power spectrum multipoles, including wide-angle effects and the survey window function. We employ this framework to derive the deconvolved power spectrum multipoles of 6dFGS DR3, BOSS DR12 and the eBOSS DR16 quasar sample. As an alternative to the s… ▽ More We make use of recent developments in the analysis of galaxy redshift surveys to present an easy to use matrix-based analysis framework for the galaxy power spectrum multipoles, including wide-angle effects and the survey window function. We employ this framework to derive the deconvolved power spectrum multipoles of 6dFGS DR3, BOSS DR12 and the eBOSS DR16 quasar sample. As an alternative to the standard analysis, the deconvolved power spectrum multipoles can be used to perform a data analysis agnostic of survey specific aspects, like the window function. We show that in the case of the BOSS dataset, the Baryon Acoustic Oscillation (BAO) analysis using the deconvolved power spectra results in the same likelihood as the standard analysis. To facilitate the analysis based on both the convolved and deconvolved power spectrum measurements, we provide the window function matrices, wide-angle matrices, covariance matrices and the power spectrum multipole measurements for the datasets mentioned above. Together with this paper we publish a \code{Python}-based toolbox to calculate the different analysis components. The appendix contains a detailed user guide with examples for how a cosmological analysis of these datasets could be implemented. We hope that our work makes the analysis of galaxy survey datasets more accessible to the wider cosmology community. △ Less

Submitted 27 October, 2021; v1 submitted 11 June, 2021; originally announced June 2021.

Comments: 39 pages, 16 figures, one table

Torsional rigidity and isospectral planar sets

Authors: Joseph Comer, Patrick McDonald

Abstract: We prove that a certain pair of isospectral planar sets are distinguished by torsional rigidity. We prove that a certain pair of isospectral planar sets are distinguished by torsional rigidity. △ Less

Submitted 16 May, 2021; originally announced May 2021.

Comments: 10 pages, 1 figure

MSC Class: 58J53 (primary) 35J25 (secondary)

Water as a Levy rotor

Authors: David A. Faux, Arifah A. Rahaman, Peter J. McDonald

Abstract: A probability density function describing the angular evolution of a fixed-length atom-atom vector as a Lévy rotor is derived containing just two dynamical parameters: the Lévy parameter $α$ and a rotational time constant $τ$. A Lévy parameter $α\!<\!2$ signals anomalous (non-Brownian) motion. A molecular dynamics simulation of water at 298\,K validates the probability density function for the int… ▽ More A probability density function describing the angular evolution of a fixed-length atom-atom vector as a Lévy rotor is derived containing just two dynamical parameters: the Lévy parameter $α$ and a rotational time constant $τ$. A Lévy parameter $α\!<\!2$ signals anomalous (non-Brownian) motion. A molecular dynamics simulation of water at 298\,K validates the probability density function for the intra-molecular $^1$H--$^1$H dynamics of water. The rotational dynamics of water is found to be approximately Brownian at sub-picosecond time intervals but becomes increasingly anomalous at longer times due to hydrogen-bond breaking and reforming. The rotational time constant lies in the range $8 \! < \! τ\! < \! 11$\,ps. The Lévy rotor model is used to estimate the intra-molecular contribution to the longitudinal nuclear-magnetic-resonance relaxation rate $R_{1,{\rm intra}}$ due to dipolar $^1$H--$^1$H interactions. It is found that $R_{1,{\rm intra}}$ contributes $65\,\pm 7$\% to the overall relaxation rate of water at room temperature. △ Less

Submitted 14 September, 2021; v1 submitted 26 April, 2021; originally announced April 2021.

Comments: 4 pages, 3 figures, for the main paper. Supplementary material included and in the right place

Extremizing Temperature Functions of Rods with Robin Boundary Conditions

Authors: Jeffrey J. Langford, Patrick McDonald

Abstract: 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 z… ▽ More 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. △ Less

Submitted 23 January, 2021; originally announced January 2021.

Comments: 18 pages

MSC Class: Primary 34B08; Secondary 34C10

An emulator for the Lyman-$α$ forest in beyond-$Λ$CDM cosmologies

Authors: Christian Pedersen, Andreu Font-Ribera, Keir K. Rogers, Patrick McDonald, Hiranya V. Peiris, Andrew Pontzen, Anže Slosar

Abstract: Interpreting observations of the Lyman-$α$ forest flux power spectrum requires interpolation between a small number of expensive simulations. We present a Gaussian process emulator modelling the 1D flux power spectrum as a function of the amplitude and slope of the small-scale linear matter power spectrum, and the state of the intergalactic medium at the epoch of interest ($2 < z < 4$). This param… ▽ More Interpreting observations of the Lyman-$α$ forest flux power spectrum requires interpolation between a small number of expensive simulations. We present a Gaussian process emulator modelling the 1D flux power spectrum as a function of the amplitude and slope of the small-scale linear matter power spectrum, and the state of the intergalactic medium at the epoch of interest ($2 < z < 4$). This parameterisation enables the prediction of the flux power spectrum in extended cosmological models that are not explicitly included in the training set, eliminating the need to construct bespoke emulators for a number of extensions to $Λ$CDM. Our emulator is appropriate for cosmologies in which the linear matter power spectrum is described to percent level accuracy by just an amplitude and slope across the epoch of interest, and in the regime probed by eBOSS/DESI data. We demonstrate this for massive neutrino cosmologies, where the emulator is able to predict the flux power spectrum in a $Σm_ν=0.3$ eV neutrino cosmology to sub-percent accuracy, without including massive neutrinos in the training simulations. Further parameters would be required to describe models with sharp features in the linear power, such as warm or light axion dark matter. This work will facilitate the combination of upcoming DESI data with observations of the cosmic microwave background, to obtain constraints on neutrino mass and other extensions to $Λ$CDM cosmology. △ Less

Submitted 14 May, 2021; v1 submitted 30 November, 2020; originally announced November 2020.

Comments: 26 pages, 11 figures. Minor changes to match version published in JCAP

Journal ref: JCAP 05 (2021) 033

International observational campaign of the 2014 eclipse of EE Cep

Authors: D. Pieńkowski, C. Gałan, T. Tomov, K. Gazeas, P. Wychudzki, M. Mikołajewski, D. Kubicki, B. Staels, S. Zoła, P. Pakońska, B. Dȩbski, T. Kundera, W. Ogłoza, M. Dróżdż, A. Baran, M. Winiarski, M. Siwak, D. Dimitrov, D. Kjurkchieva, D. Marchev, A. Armiński, I. Miller, Z. Kołaczkowski, D. Moździerski, E. Zahajkiewicz , et al. (44 additional authors not shown)

Abstract: Context. EE Cep is one of few eclipsing binary systems with a dark, dusty disk around an invisible object similar to ε Aur. The system is characterized by grey and asymmetric eclipses every 5.6 yr, with a significant variation in their photometric depth, ranging from ~ 0 m .5 to ~ 2 m .0. Aims. The main aim of the observational campaign of the EE Cep eclipse in 2014 was to test the model of disk p… ▽ More Context. EE Cep is one of few eclipsing binary systems with a dark, dusty disk around an invisible object similar to ε Aur. The system is characterized by grey and asymmetric eclipses every 5.6 yr, with a significant variation in their photometric depth, ranging from ~ 0 m .5 to ~ 2 m .0. Aims. The main aim of the observational campaign of the EE Cep eclipse in 2014 was to test the model of disk precession (Galan et al. 2012). We expected that this eclipse would be one of the deepest with a depth of ~ 2 m .0. Methods. We collected multicolor observations from almost 30 instruments located in Europe and North America. This photometric data covers 243 nights during and around the eclipse. We also analyse the low- and high-resolution spectra from several instruments. Results. The eclipse was shallow with a depth of 0 m .71 in V-band. The multicolor photometry illustrates small color changes during the eclipse with a total amplitude of order ~ +0 m . 15 in B-I color index. The linear ephemeris for this system is updated by including new times of minima, measured from the three most recent eclipses at epochs E = 9, 10 and 11. New spectroscopic observations were acquired, covering orbital phases around the eclipse, which were not observed in the past and increased the data sample, filling some gaps and giving a better insight into the evolution of the H α and NaI spectral line profiles during the primary eclipse. Conclusions. The eclipse of EE Cep in 2014 was shallower than expected 0 m .71 instead of ~ 2 m . 0. This means that our model of disk precession needs revision. △ Less

Submitted 16 January, 2020; originally announced January 2020.

Comments: 10 pages, 13 figures, 29 tables in appendix, submited to Astronomy & Astrophysics

Massive neutrinos and degeneracies in Lyman-alpha forest simulations

Authors: Christian Pedersen, Andreu Font-Ribera, Thomas D. Kitching, Patrick McDonald, Simeon Bird, Anže Slosar, Keir K. Rogers, Andrew Pontzen

Abstract: 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 suppres… ▽ More 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. △ Less

Submitted 17 March, 2020; v1 submitted 21 November, 2019; originally announced November 2019.

Comments: 18 pages, 7 figures

Astro2020 APC White Paper: The MegaMapper: a z > 2 spectroscopic instrument for the study of Inflation and Dark Energy

Authors: David J. Schlegel, Juna A. Kollmeier, Greg Aldering, Stephen Bailey, Charles Baltay, Christopher Bebek, Segev BenZvi, Robert Besuner, Guillermo Blanc, Adam S. Bolton, Mohamed Bouri, David Brooks, Elizabeth Buckley-Geer, Zheng Cai, Jeffrey Crane, Arjun Dey, Peter Doel, Xiaohui Fan, Simone Ferraro, Andreu Font-Ribera, Gaston Gutierrez, Julien Guy, Henry Heetderks, Dragan Huterer, Leopoldo Infante , et al. (52 additional authors not shown)

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

Submitted 25 July, 2019; originally announced July 2019.

Improved renormalization group computation of likelihood functions for cosmological data sets

Authors: Patrick McDonald

Abstract: Evaluation of likelihood functions for cosmological large scale structure data sets (including CMB, galaxy redshift surveys, etc.) naturally involves marginalization, i.e., integration, over an unknown underlying random signal field. Recently, I showed how a renormalization group method can be used to carry out this integration efficiently by first integrating out the smallest scale structure, i.e… ▽ More Evaluation of likelihood functions for cosmological large scale structure data sets (including CMB, galaxy redshift surveys, etc.) naturally involves marginalization, i.e., integration, over an unknown underlying random signal field. Recently, I showed how a renormalization group method can be used to carry out this integration efficiently by first integrating out the smallest scale structure, i.e., localized structure on the scale of differences between nearby data cells, then combining adjacent cells in a coarse graining step, then repeating this process over and over until all scales have been integrated. Here I extend the formulation in several ways in order to reduce the prefactor on the method's linear scaling with data set size. The key improvement is showing how to integrate out the difference between specific adjacent cells before summing them in the coarse graining step, compared to the original formulation in which small-scale fluctuations were integrated more generally. I suggest some other improvements in details of the scheme, including showing how to perform the integration around a maximum likelihood estimate for the underlying random field. In the end, an accurate likelihood computation for a million-cell Gaussian test data set runs in two minutes on my laptop, with room for further optimization and straightforward parallelization. △ Less

Submitted 21 June, 2019; originally announced June 2019.

Comments: 10 pg, 2 fig

Journal ref: Phys. Rev. D 100, 043511 (2019)

Inflation and Dark Energy from spectroscopy at $z > 2$

Authors: Simone Ferraro, Michael J. Wilson, Muntazir Abidi, David Alonso, Behzad Ansarinejad, Robert Armstrong, Jacobo Asorey, Arturo Avelino, Carlo Baccigalupi, Kevin Bandura, Nicholas Battaglia, Chetan Bavdhankar, José Luis Bernal, Florian Beutler, Matteo Biagetti, Guillermo A. Blanc, Jonathan Blazek, Adam S. Bolton, Julian Borrill, Brenda Frye, Elizabeth Buckley-Geer, Philip Bull, Cliff Burgess, Christian T. Byrnes, Zheng Cai , et al. (118 additional authors not shown)

Abstract: The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at… ▽ More The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $σ(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space. △ Less

Submitted 21 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey

Testing Gravity Using Type Ia Supernovae Discovered by Next-Generation Wide-Field Imaging Surveys

Authors: A. G. Kim, G. Aldering, P. Antilogus, A. Bahmanyar, S. BenZvi, H. Courtois, T. Davis, H. Feldman, S. Ferraro, S. Gontcho A Gontcho, O. Graur, R. Graziani, J. Guy, C. Harper, R. Hložek, C. Howlett, D. Huterer, C. Ju, P. -F. Leget, E. V. Linder, P. McDonald, J. Nordin, P. Nugent, S. Perlmutter, N. Regnault , et al. (7 additional authors not shown)

Abstract: In the upcoming decade cadenced wide-field imaging surveys will increase the number of identified $z<0.3$ Type~Ia supernovae (SNe~Ia) from the hundreds to the hundreds of thousands. The increase in the number density and solid-angle coverage of SNe~Ia, in parallel with improvements in the standardization of their absolute magnitudes, now make them competitive probes of the growth of structure and… ▽ More In the upcoming decade cadenced wide-field imaging surveys will increase the number of identified $z<0.3$ Type~Ia supernovae (SNe~Ia) from the hundreds to the hundreds of thousands. The increase in the number density and solid-angle coverage of SNe~Ia, in parallel with improvements in the standardization of their absolute magnitudes, now make them competitive probes of the growth of structure and hence of gravity. The peculiar velocity power spectrum is sensitive to the growth index $γ$, which captures the effect of gravity on the linear growth of structure through the relation $f=Ω_M^γ$. We present the first projections for the precision in $γ$ for a range of realistic SN peculiar-velocity survey scenarios. In the next decade the peculiar velocities of SNe~Ia in the local $z<0.3$ Universe will provide a measure of $γ$ to $\pm 0.01$ precision that can definitively distinguish between General Relativity and leading models of alternative gravity. △ Less

Submitted 18 March, 2019; originally announced March 2019.

Comments: Decadal Survey White Paper submission

Renormalization group computation of likelihood functions for cosmological data sets

Authors: Patrick McDonald

Abstract: I show how a renormalization group (RG) method can be used to incrementally integrate the information in cosmological large-scale structure data sets (including CMB, galaxy redshift surveys, etc.). I show numerical tests for Gaussian fields, where the method allows arbitrarily close to exact computation of the likelihood function in order $\sim N$ time, even for problems with no symmetry, compared… ▽ More I show how a renormalization group (RG) method can be used to incrementally integrate the information in cosmological large-scale structure data sets (including CMB, galaxy redshift surveys, etc.). I show numerical tests for Gaussian fields, where the method allows arbitrarily close to exact computation of the likelihood function in order $\sim N$ time, even for problems with no symmetry, compared to $N^3$ for brute force linear algebra (where $N$ is the number of data points -- to be fair, methods already exist to solve the Gaussian problem in at worst $N \log N$ time, and this method will not necessarily be faster in practice). The method requires no sampling or other Monte Carlo (random) element. Non-linearity/non-Gaussianity can be accounted for to the extent that terms generated by integrating out small scale modes can be projected onto a sufficient basis, e.g., at least in the sufficiently perturbative regime. The formulas to evaluate are straightforward and require no understanding of quantum field theory, but this paper may also serve as a pedagogical introduction to Wilsonian RG for astronomers. △ Less

Submitted 19 October, 2018; originally announced October 2018.

Comments: 13 pg, 4 fig

Journal ref: Phys. Rev. D 99, 043538 (2019)

How to estimate the 3D power spectrum of the Lyman-$α$ forest

Authors: Andreu Font-Ribera, Patrick McDonald, Anže Slosar

Abstract: We derive and numerically implement an algorithm for estimating the 3D power spectrum of the Lyman-$α$ (Ly-$α$) forest flux fluctuations. The algorithm exploits the unique geometry of Ly-$α$ forest data to efficiently measure the cross-spectrum between lines of sight as a function of parallel wavenumber, transverse separation and redshift. The key to fast evaluation is to approximate the global co… ▽ More We derive and numerically implement an algorithm for estimating the 3D power spectrum of the Lyman-$α$ (Ly-$α$) forest flux fluctuations. The algorithm exploits the unique geometry of Ly-$α$ forest data to efficiently measure the cross-spectrum between lines of sight as a function of parallel wavenumber, transverse separation and redshift. The key to fast evaluation is to approximate the global covariance matrix as block-diagonal, where only pixels from the same spectrum are correlated. We then compute the eigenvectors of the derivative of the signal covariance with respect to cross-spectrum parameters, and project the inverse-covariance-weighted spectra onto them. This acts much like a radial Fourier transform over redshift windows. The resulting cross-spectrum inference is then converted into our final product, an approximation of the likelihood for the 3D power spectrum expressed as second order Taylor expansion around a fiducial model. We demonstrate the accuracy and scalability of the algorithm and comment on possible extensions. Our algorithm will allow efficient analysis of the upcoming Dark Energy Spectroscopic Instrument dataset. △ Less

Submitted 29 December, 2017; v1 submitted 30 October, 2017; originally announced October 2017.

Comments: 29 pages, many figures. Minor changes in v2, accepted in JCAP

Large-scale structure perturbation theory without losing stream crossing

Authors: Patrick McDonald, Zvonimir Vlah

Abstract: We suggest an approach to perturbative calculations of large-scale clustering in the Universe that includes from the start the stream crossing (multiple velocities for mass elements at a single position) that is lost in traditional calculations. Starting from a functional integral over displacement, the perturbative series expansion is in deviations from (truncated) Zel'dovich evolution, with term… ▽ More We suggest an approach to perturbative calculations of large-scale clustering in the Universe that includes from the start the stream crossing (multiple velocities for mass elements at a single position) that is lost in traditional calculations. Starting from a functional integral over displacement, the perturbative series expansion is in deviations from (truncated) Zel'dovich evolution, with terms that can be computed exactly even for stream-crossed displacements. We evaluate the one-loop formulas for displacement and density power spectra numerically in 1D, finding dramatic improvement in agreement with N-body simulations compared to the Zel'dovich power spectrum (which is exact in 1D up to stream crossing). Beyond 1D, our approach could represent an improvement over previous expansions even aside from the inclusion of stream crossing, but we have not investigated this numerically. In the process we show how to achieve effective-theory-like regulation of small-scale fluctuations without free parameters. △ Less

Submitted 5 January, 2018; v1 submitted 8 September, 2017; originally announced September 2017.

Comments: added pedagogical explanation of key math trick in appendix

Journal ref: Phys. Rev. D 97, 023508 (2018)

Comparison Results, Exit Time Moments, And Eigenvalues On Riemannian Manifolds With A Lower Ricci Curvature Bound

Authors: Don Colladay, Jeffrey J. Langford, Patrick McDonald

Abstract: We study the relationship between the geometry of smoothly bounded domains in complete Riemannian manifolds and the associated sequence of $L^1$-norms of exit time moments for Brownian motion. We establish bounds for Dirichlet eigenvalues and, for closed manifolds, we establish a comparison result for elements of the moment sequence. We study the relationship between the geometry of smoothly bounded domains in complete Riemannian manifolds and the associated sequence of $L^1$-norms of exit time moments for Brownian motion. We establish bounds for Dirichlet eigenvalues and, for closed manifolds, we establish a comparison result for elements of the moment sequence. △ Less

Submitted 13 June, 2017; v1 submitted 6 June, 2017; originally announced June 2017.

Comments: 21 pages. Added results for the case of equality

MSC Class: 58J65; 58J50; 35P15

Generalizing the Minkowski Question Mark Function to a Family of Multidimensional Continued Fractions

Authors: Thomas Garrity, Peter McDonald

Abstract: The Minkowski question mark function, maping the unit interval to itself, is a continuous, strictly increasing, one-to-one and onto function that has derivative zero almost everywhere. Key to these facts are the basic properties of continued fractions. Thus the question mark function is a naturally occurring number theoretic singular function. This paper generalizes the question mark function to t… ▽ More The Minkowski question mark function, maping the unit interval to itself, is a continuous, strictly increasing, one-to-one and onto function that has derivative zero almost everywhere. Key to these facts are the basic properties of continued fractions. Thus the question mark function is a naturally occurring number theoretic singular function. This paper generalizes the question mark function to the 216 triangle partition (TRIP) maps. These are multidimensional continued fractions which generate a family of almost all known multidimensional continued fractions. We show for each TRIP map that there is a natural candidate for its analog of the Minkowski question mark function. We then show that the analog is singular for 96 of the TRIP maps and show that 60 more are singular under an assumption of ergodicity. △ Less

Submitted 21 February, 2017; v1 submitted 31 January, 2017; originally announced January 2017.

Comments: 40 pages, font problems fixed

The DESI Experiment Part II: Instrument Design

Authors: DESI Collaboration, Amir Aghamousa, Jessica Aguilar, Steve Ahlen, Shadab Alam, Lori E. Allen, Carlos Allende Prieto, James Annis, Stephen Bailey, Christophe Balland, Otger Ballester, Charles Baltay, Lucas Beaufore, Chris Bebek, Timothy C. Beers, Eric F. Bell, José Luis Bernal, Robert Besuner, Florian Beutler, Chris Blake, Hannes Bleuler, Michael Blomqvist, Robert Blum, Adam S. Bolton, Cesar Briceno , et al. (268 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 13 December, 2016; v1 submitted 31 October, 2016; originally announced November 2016.

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

Authors: DESI Collaboration, Amir Aghamousa, Jessica Aguilar, Steve Ahlen, Shadab Alam, Lori E. Allen, Carlos Allende Prieto, James Annis, Stephen Bailey, Christophe Balland, Otger Ballester, Charles Baltay, Lucas Beaufore, Chris Bebek, Timothy C. Beers, Eric F. Bell, José Luis Bernal, Robert Besuner, Florian Beutler, Chris Blake, Hannes Bleuler, Michael Blomqvist, Robert Blum, Adam S. Bolton, Cesar Briceno , et al. (268 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 13 December, 2016; v1 submitted 31 October, 2016; originally announced November 2016.

Covariant Quantization of CPT-violating Photons

Authors: D. Colladay, P. McDonald, J. P. Noordmans, R. Potting

Abstract: We perform the covariant canonical quantization of the CPT- and Lorentz-symmetry-violating photon sector of the minimal Standard-Model Extension, which contains a general (timelike, lightlike, or spacelike) fixed background tensor $k_{AF}^μ$. Well-known stability issues, arising from complex-valued energy states, are solved by introducing a small photon mass, orders of magnitude below current expe… ▽ More We perform the covariant canonical quantization of the CPT- and Lorentz-symmetry-violating photon sector of the minimal Standard-Model Extension, which contains a general (timelike, lightlike, or spacelike) fixed background tensor $k_{AF}^μ$. Well-known stability issues, arising from complex-valued energy states, are solved by introducing a small photon mass, orders of magnitude below current experimental bounds. We explicitly construct a covariant basis of polarization vectors, in which the photon field can be expanded. We proceed to derive the Feynman propagator and show that the theory is microcausal. Despite the occurrence of negative energies and vacuum-Cherenkov radiation, we do not find any runaway stability issues, because the energy remains bounded from below. An important observation is that the ordering of the roots of the dispersion relations is the same in any observer frame, which allows for a frame-independent condition that selects the correct branch of the dispersion relation. This turns out to be critical for the consistency of the quantization. To our knowledge, this is the first system for which quantization has consistently been performed, in spite of the fact that the theory contains negative energies in some observer frames. △ Less

Submitted 1 October, 2016; originally announced October 2016.

Comments: 38 pages, 2 figures

Journal ref: Phys. Rev. D 95, 025025 (2017)

Exit Time Moments and Eigenvalue Estimates

Authors: Emily B. Dryden, Jeffrey J. Langford, Patrick McDonald

Abstract: We study estimates involving the principal Dirichlet eigenvalue associated to a smoothly bounded domain in a complete Riemannian manifold and L1-norms of exit time moments of Brownian motion. Our results generalize a classical inequality of Polya. We study estimates involving the principal Dirichlet eigenvalue associated to a smoothly bounded domain in a complete Riemannian manifold and L1-norms of exit time moments of Brownian motion. Our results generalize a classical inequality of Polya. △ Less

Submitted 2 January, 2017; v1 submitted 31 August, 2016; originally announced August 2016.

Comments: 15 pages, added appendix on variational quotients, added lower bounds, improved exposition

MSC Class: 35P15

Torsional Rigidity, Isospectrality and Quantum Graphs

Authors: Don Colladay, Leon Kaganovskiy, Patrick McDonald

Abstract: We study torsional rigidity for graph and quantum graph analogs of well-known pairs of isospectral non-isometric planar domains. We prove that such isospectral pairs are distinguished by torsional rigidity. We study torsional rigidity for graph and quantum graph analogs of well-known pairs of isospectral non-isometric planar domains. We prove that such isospectral pairs are distinguished by torsional rigidity. △ Less

Submitted 19 July, 2016; originally announced July 2016.

Comments: 15 pages, 5 figures

MSC Class: 58J53

Heat content determines planar triangles

Authors: Reed Meyerson, Patrick McDonald

Abstract: We prove that heat content determines planar triangles. We prove that heat content determines planar triangles. △ Less

Submitted 12 July, 2016; originally announced July 2016.

Comments: 9 pages, 1 figure

MSC Class: 58J50

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

Authors: Florian Beutler, Hee-Jong Seo, Ashley J. Ross, Patrick McDonald, Shun Saito, Adam S. Bolton, Joel R. Brownstein, Chia-Hsun Chuang, Antonio J. Cuesta, Daniel J. Eisenstein, Andreu Font-Ribera, Jan Niklas Grieb, Nick Hand, Francisco-Shu Kitaura, Chirag Modi, Robert C. Nichol, Will J. Percival, Francisco Prada, Sergio Rodriguez-Torres, Natalie A. Roe, Nicholas P. Ross, Salvador Salazar-Albornoz, Ariel G. Sánchez, Donald P. Schneider, Anže Slosar , et al. (4 additional authors not shown)

Abstract: 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 re… ▽ More 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. △ Less

Submitted 11 July, 2016; originally announced July 2016.

Explicit calculation of nuclear magnetic resonance relaxation rates in small pores to elucidate molecular scale fluid dynamics

Authors: David A Faux, Peter J McDonald

Abstract: A model linking the molecular-scale dynamics of fluids confined to nano-pores to nuclear magnetic resonance (NMR) relaxation rates is proposed. The model is fit to experimental NMR dispersions for water and oil in an oil shale assuming that each fluid is characterised by three time constants and Lévy statistics. Results yield meaningful and consistent intra-pore dynamical time constants, insight i… ▽ More A model linking the molecular-scale dynamics of fluids confined to nano-pores to nuclear magnetic resonance (NMR) relaxation rates is proposed. The model is fit to experimental NMR dispersions for water and oil in an oil shale assuming that each fluid is characterised by three time constants and Lévy statistics. Results yield meaningful and consistent intra-pore dynamical time constants, insight into diffusion mechanisms and pore morphology. The model is applicable to a wide range of porous systems and advances NMR dispersion as a powerful tool for measuring nano-porous fluid properties. △ Less

Submitted 29 June, 2016; originally announced June 2016.

Comments: 4 pages, 4 figures, submitted to Phys. Rev. Lett. and under review

Journal ref: Phys. Rev. E 95, 033117 (2017)

The Clustering of Galaxies in the Completed SDSS-III Baryon Oscillation Spectroscopic Survey: Cosmic Flows and Cosmic Web from Luminous Red Galaxies

Authors: Metin Ata, Francisco-Shu Kitaura, Chia-Hsun Chuang, Sergio Rodríguez-Torres, Raul E. Angulo, Simone Ferraro, Hector Gil-Marín, Patrick McDonald, Carlos Hernández Monteagudo, Volker Müller, Gustavo Yepes, Mathieu Autefage, Falk Baumgarten, Florian Beutler, Joel R. Brownstein, Angela Burden, Daniel J. Eisenstein, Hong Guo, Shirley Ho, Cameron McBride, Mark Neyrinck, Matthew D. Olmstead, Nikhil Padmanabhan, Will J. Percival, Francisco Prada , et al. (9 additional authors not shown)

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

Submitted 17 January, 2017; v1 submitted 31 May, 2016; originally announced May 2016.

Comments: 21 pages, 17 figures submitted to MNRAS

Linear redshift space distortions for cosmic voids based on galaxies in redshift space

Authors: Chia-Hsun Chuang, Francisco-Shu Kitaura, Yu Liang, Andreu Font-Ribera, Cheng Zhao, Patrick McDonald, Charling Tao

Abstract: Cosmic voids found in galaxy surveys are defined based on the galaxy distribution in redshift space. We show that the large scale distribution of voids in redshift space traces the fluctuations in the dark matter density field δ(k) (in Fourier space with μbeing the line of sight projected k-vector): δ_v^s(k) = (1 + β_v μ^2) b^s_v δ(k), with a beta factor that will be in general different than the… ▽ More Cosmic voids found in galaxy surveys are defined based on the galaxy distribution in redshift space. We show that the large scale distribution of voids in redshift space traces the fluctuations in the dark matter density field δ(k) (in Fourier space with μbeing the line of sight projected k-vector): δ_v^s(k) = (1 + β_v μ^2) b^s_v δ(k), with a beta factor that will be in general different than the one describing the distribution of galaxies. Only in case voids could be assumed to be quasi-local transformations of the linear (Gaussian) galaxy redshift space field, one gets equal beta factors β_v=β_g=f/b_g with f being the growth rate, and b_g, b^s_v being the galaxy and void bias on large scales defined in redshift space. Indeed, in our mock void catalogs we measure void beta factors being in good agreement with the galaxy one. Further work needs to be done to confirm the level of accuracy of the beta factor equality between voids and galaxies, but in general the void beta factor needs to be considered as a free parameter for RSD studies. △ Less

Submitted 3 August, 2017; v1 submitted 17 May, 2016; originally announced May 2016.

Comments: 10 pages, 9 figures; matches the version accepted by PRD

Journal ref: Phys. Rev. D 95, 063528 (2017)

Completely Controlling the Dimensions of Formal Fiber Rings at Prime Ideals of Small Height

Authors: Sarah M. Fleming, Lena Ji, S. Loepp, Peter M. McDonald, Nina Pande, David Schwein

Abstract: Let $T$ be a complete equicharacteristic local (Noetherian) UFD of dimension $3$ or greater. Assuming that $|T| = |T/m|$, where $m$ is the maximal ideal of $T$, we construct a local UFD $A$ whose completion is $T$ and whose formal fibers at height one prime ideals have prescribed dimension between zero and the dimension of the generic formal fiber. If, in addition, $T$ is regular and has character… ▽ More Let $T$ be a complete equicharacteristic local (Noetherian) UFD of dimension $3$ or greater. Assuming that $|T| = |T/m|$, where $m$ is the maximal ideal of $T$, we construct a local UFD $A$ whose completion is $T$ and whose formal fibers at height one prime ideals have prescribed dimension between zero and the dimension of the generic formal fiber. If, in addition, $T$ is regular and has characteristic zero, we can construct $A$ to be excellent. △ Less

Submitted 14 April, 2016; originally announced April 2016.

Comments: 20 pages

Fast Large Scale Structure Perturbation Theory using 1D FFTs

Authors: Marcel Schmittfull, Zvonimir Vlah, Patrick McDonald

Abstract: 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 locati… ▽ More 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. △ Less

Submitted 5 May, 2016; v1 submitted 14 March, 2016; originally announced March 2016.

Comments: 10+4 pages, published version

Journal ref: Phys. Rev. D 93, 103528 (2016)

FastPM: a new scheme for fast simulations of dark matter and halos

Authors: Yu Feng, Man-Yat Chu, Uros Seljak, Patrick McDonald

Abstract: We introduce FastPM, a highly-scalable approximated particle mesh N-body solver, which implements the particle mesh (PM) scheme enforcing correct linear displacement (1LPT) evolution via modified kick and drift factors. Employing a 2-dimensional domain decomposing scheme, FastPM scales extremely well with a very large number of CPUs. In contrast to COmoving-LAgrangian (COLA) approach, we do not re… ▽ More We introduce FastPM, a highly-scalable approximated particle mesh N-body solver, which implements the particle mesh (PM) scheme enforcing correct linear displacement (1LPT) evolution via modified kick and drift factors. Employing a 2-dimensional domain decomposing scheme, FastPM scales extremely well with a very large number of CPUs. In contrast to COmoving-LAgrangian (COLA) approach, we do not require to split the force or track separately the 2LPT solution, reducing the code complexity and memory requirements. We compare FastPM with different number of steps ($N_s$) and force resolution factor ($B$) against 3 benchmarks: halo mass function from Friends of Friends halo finder, halo and dark matter power spectrum, and cross correlation coefficient (or stochasticity), relative to a high resolution TreePM simulation. We show that the modified time stepping scheme reduces the halo stochasticity when compared to COLA with the same number of steps and force resolution. While increasing $N_s$ and $B$ improves the transfer function and cross correlation coefficient, for many applications FastPM achieves sufficient accuracy at low $N_s$ and $B$. For example, $N_s=10$ and $B=2$ simulation provides a substantial saving (a factor of 10) of computing time relative to $N_s=40$, $B=3$ simulation, yet the halo benchmarks are very similar at $z=0$. We find that for abundance matched halos the stochasticity remains low even for $N_s=5$. FastPM compares well against less expensive schemes, being only 7 (4) times more expensive than 2LPT initial condition generator for $N_s=10$ ($N_s=5$). Some of the applications where FastPM can be useful are generating a large number of mocks, producing non-linear statistics where one varies a large number of nuisance or cosmological parameters, or serving as part of an initial conditions solver. △ Less

Submitted 12 September, 2017; v1 submitted 1 March, 2016; originally announced March 2016.

Comments: 14 pages, 13 figures. The development is hosted at https://github.com/rainwoodman/fastpm . Major revision addressing referee report from MNRAS; Update code repository URL in comments

Cherenkov Radiation with Massive, CPT-violating Photons

Authors: Don Colladay, Patrick McDonald, Robertus Potting

Abstract: 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… ▽ More 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. △ Less

Submitted 23 March, 2016; v1 submitted 1 March, 2016; originally announced March 2016.

Journal ref: Phys. Rev. D 93, 125007 (2016)

Controlling the Dimensions of Formal Fibers of a Unique Factorization Domain at the Height One Prime Ideals

Authors: Sarah M. Fleming, Lena Ji, S. Loepp, Peter M. McDonald, Nina Pande, David Schwein

Abstract: Let T be a complete local (Noetherian) equidimensional ring with maximal ideal m such that the Krull dimension of T is at least two and the depth of T is at least two. Suppose that no integer of T is a zerodivisor and that |T|=|T/m|. Let d and t be integers such that 1 $\leq$ d $\leq$ dimT-1, 0 $\leq$ t $\leq$ dimT - 1, and d - 1 $\leq$ t. Assume that, for every p in AssT, ht(p) $\leq$ d-1 and tha… ▽ More Let T be a complete local (Noetherian) equidimensional ring with maximal ideal m such that the Krull dimension of T is at least two and the depth of T is at least two. Suppose that no integer of T is a zerodivisor and that |T|=|T/m|. Let d and t be integers such that 1 $\leq$ d $\leq$ dimT-1, 0 $\leq$ t $\leq$ dimT - 1, and d - 1 $\leq$ t. Assume that, for every p in AssT, ht(p) $\leq$ d-1 and that if z is a regular element of T and Q is in Ass(T/zT), then ht(Q) $\leq$ d. We construct a local unique factorization domain A such that the completion of A is T and such that the dimension of the formal fiber ring at every height one prime ideal of A is d - 1 and the dimension of the formal fiber ring of A at (0) is t. △ Less

Submitted 26 January, 2016; originally announced January 2016.

Comments: 15 pages

MSC Class: 13J10

Singular Lorentz-Violating Lagrangians and Associated Finsler Structures

Authors: Don Colladay, Patrick McDonald

Abstract: 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… ▽ More 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. △ Less

Submitted 3 July, 2015; originally announced July 2015.

Comments: 11 pages