Showing 1–50 of 66 results for author: Slepian, Z

Parity-Odd Power Spectra: Concise Statistics for Cosmological Parity Violation

Authors: Drew Jamieson, Angelo Caravano, Jiamin Hou, Zachary Slepian, Eiichiro Komatsu

Abstract: We introduce the Parity-Odd Power (POP) spectra, a novel set of observables for probing parity violation in cosmological $N$-point statistics. POP spectra are derived from composite fields obtained by applying nonlinear transformations, involving also gradients, curls, and filtering functions, to a scalar field. This compresses the parity-odd trispectrum into a power spectrum. These new statistics… ▽ More We introduce the Parity-Odd Power (POP) spectra, a novel set of observables for probing parity violation in cosmological $N$-point statistics. POP spectra are derived from composite fields obtained by applying nonlinear transformations, involving also gradients, curls, and filtering functions, to a scalar field. This compresses the parity-odd trispectrum into a power spectrum. These new statistics offer several advantages: they are computationally fast to construct, estimating their covariance is less demanding compared to estimating that of the full parity-odd trispectrum, and they are simple to model theoretically. We measure the POP spectra on simulations of a scalar field with a specific parity-odd trispectrum shape. We compare these measurements to semi-analytic theoretical calculations and find agreement. We also explore extensions and generalizations of these parity-odd observables. △ Less

Submitted 21 June, 2024; originally announced June 2024.

Comments: 18 pages, 5 figures

On a Generating Function for the Isotropic Basis Functions and Other Connected Results

Authors: Zachary Slepian, Jessica Chellino, Jiamin Hou

Abstract: Recently isotropic basis functions of $N$ unit vector arguments were presented; these are of significant use in measuring the N-Point Correlation Functions (NPCFs) of galaxy clustering. Here we develop the generating function for these basis functions -- i.e. that function which, expanded in a power series, has as its angular part the isotropic functions. We show that this can be developed using b… ▽ More Recently isotropic basis functions of $N$ unit vector arguments were presented; these are of significant use in measuring the N-Point Correlation Functions (NPCFs) of galaxy clustering. Here we develop the generating function for these basis functions -- i.e. that function which, expanded in a power series, has as its angular part the isotropic functions. We show that this can be developed using basic properties of the plane wave. A main use of the generating function is as an efficient route to obtaining the Cartesian basis expressions for the isotropic functions. We show that the methods here enable computing difficult overlap integrals of multiple spherical Bessel functions, and we also give related expansions of the Dirac Delta function into the isotropic basis. Finally, we outline how the Cartesian expressions for the isotropic basis functions might be used to enable a faster NPCF algorithm on the CPU. △ Less

Submitted 20 April, 2024; originally announced June 2024.

Comments: 28 pages, no figures, comments welcome

A Model for the Redshift-Space Galaxy 4-Point Correlation Function

Authors: William Ortolá Leonard, Zachary Slepian, Jiamin Hou

Abstract: The field of cosmology is entering an epoch of unparalleled wealth of observational data thanks to galaxy surveys such as DESI, Euclid, and Roman. Therefore, it is essential to have a firm theoretical basis that allows the effective analysis of the data. With this purpose, we compute the nonlinear, gravitationally-induced connected galaxy 4-point correlation function (4PCF) at the tree level in St… ▽ More The field of cosmology is entering an epoch of unparalleled wealth of observational data thanks to galaxy surveys such as DESI, Euclid, and Roman. Therefore, it is essential to have a firm theoretical basis that allows the effective analysis of the data. With this purpose, we compute the nonlinear, gravitationally-induced connected galaxy 4-point correlation function (4PCF) at the tree level in Standard Perturbation Theory (SPT), including redshift-space distortions (RSD). We begin from the trispectrum and take its inverse Fourier transform into configuration space, exploiting the isotropic basis functions of Cahn & Slepian (2023). We ultimately reduce the configuration-space expression to low-dimensional radial integrals of the power spectrum. This model will enable the use of the BAO feature in the connected 4PCF to sharpen our constraints on the expansion history of the Universe. It will also offer an additional avenue for determining the galaxy bias parameters, and thus tighten our cosmological constraints by breaking degeneracies. Survey geometry can be corrected in the 4PCF, and many systematics are localized, which is an advantage over data analysis with the trispectrum. Finally, this work is a first step in using the parity-even 4PCF to calibrate out any possible systematics in the parity-odd 4PCF; comparing our model to the measurement will offer a consistency test that can serve as a "canary in the coal mine" for systematics in the odd sector. △ Less

Submitted 21 March, 2024; v1 submitted 23 February, 2024; originally announced February 2024.

Comments: 100 pages, 36 figures

Triple-Spherical Bessel Function Integrals with Exponential and Gaussian Damping: Towards an Analytic N-Point Correlation Function Covariance Model

Authors: Jessica Chellino, Zachary Slepian

Abstract: Spherical Bessel functions appear commonly in many areas of physics wherein there is both translation and rotation invariance, and often integrals over products of several arise. Thus, analytic evaluation of such integrals with different weighting functions (which appear as toy models of a given physical observable, such as the galaxy power spectrum) is useful. Here we present a generalization of… ▽ More Spherical Bessel functions appear commonly in many areas of physics wherein there is both translation and rotation invariance, and often integrals over products of several arise. Thus, analytic evaluation of such integrals with different weighting functions (which appear as toy models of a given physical observable, such as the galaxy power spectrum) is useful. Here we present a generalization of a recursion-based method for evaluating such integrals. It gives relatively simple closed-form results in terms of Legendre functions (for the exponentially-damped case) and Gamma, incomplete Gamma functions, and hypergeometric functions (for the Gaussian-damped case). We also present a new, non-recursive method to evaluate integrals of products of spherical Bessel functions with Gaussian damping in terms of incomplete Gamma functions and hypergeometric functions. △ Less

Submitted 21 December, 2023; v1 submitted 3 August, 2023; originally announced August 2023.

Comments: 29 pages

Journal ref: Proc. R. Soc. A. 479 (2023)

NANCY: Next-generation All-sky Near-infrared Community surveY

Authors: Jiwon Jesse Han, Arjun Dey, Adrian M. Price-Whelan, Joan Najita, Edward F. Schlafly, Andrew Saydjari, Risa H. Wechsler, Ana Bonaca, David J Schlegel, Charlie Conroy, Anand Raichoor, Alex Drlica-Wagner, Juna A. Kollmeier, Sergey E. Koposov, Gurtina Besla, Hans-Walter Rix, Alyssa Goodman, Douglas Finkbeiner, Abhijeet Anand, Matthew Ashby, Benedict Bahr-Kalus, Rachel Beaton, Jayashree Behera, Eric F. Bell, Eric C Bellm , et al. (184 additional authors not shown)

Abstract: The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GAL… ▽ More The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GALAH, 4MOST, WEAVE, MOONS, PFS, UVEX, NEO Surveyor, etc.). Roman can uniquely provide uniform high-spatial-resolution (~0.1 arcsec) imaging over the entire sky, vastly expanding the science reach and precision of all of these near-term and future surveys. This imaging will not only enhance other surveys, but also facilitate completely new science. By imaging the full sky over two epochs, Roman can measure the proper motions for stars across the entire Milky Way, probing 100 times fainter than Gaia out to the very edge of the Galaxy. Here, we propose NANCY: a completely public, all-sky survey that will create a high-value legacy dataset benefiting innumerable ongoing and forthcoming studies of the universe. NANCY is a pure expression of Roman's potential: it images the entire sky, at high spatial resolution, in a broad infrared bandpass that collects as many photons as possible. The majority of all ongoing astronomical surveys would benefit from incorporating observations of NANCY into their analyses, whether these surveys focus on nearby stars, the Milky Way, near-field cosmology, or the broader universe. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: Submitted to the call for white papers for the Roman Core Community Survey (June 16th, 2023), and to the Bulletin of the AAS

GTC Follow-up Observations of Very Metal-Poor Star Candidates from DESI

Authors: Carlos Allende Prieto, David S. Aguado, Jonay I. González Hernández, Rafael Rebolo, Joan Najita, Christopher J. Manser, Constance Rockosi, Zachary Slepian, Mar Mezcua, Monica Valluri, Rana Ezzeddine, Sergey E. Koposov, Andrew P. Cooper, Arjun Dey, Boris T. Gänsicke, Ting S. Li, Katia Cunha, Siwei Zou, Jessica Nicole Aguilar, Steven Ahlen, David Brooks, Todd Claybaugh, Shaun Cole, Sarah Eftekharzadeh, Kevin Fanning , et al. (26 additional authors not shown)

Abstract: The observations from the Dark Energy Spectroscopic Instrument (DESI) will significantly increase the numbers of known extremely metal-poor stars by a factor of ~ 10, improving the sample statistics to study the early chemical evolution of the Milky Way and the nature of the first stars. In this paper we report high signal-to-noise follow-up observations of 9 metal-poor stars identified during the… ▽ More The observations from the Dark Energy Spectroscopic Instrument (DESI) will significantly increase the numbers of known extremely metal-poor stars by a factor of ~ 10, improving the sample statistics to study the early chemical evolution of the Milky Way and the nature of the first stars. In this paper we report high signal-to-noise follow-up observations of 9 metal-poor stars identified during the DESI commissioning with the Optical System for Imaging and low-Intermediate-Resolution Integrated Spectroscopy (OSIRIS) instrument on the 10.4m Gran Telescopio Canarias (GTC). The analysis of the data using a well-vetted methodology confirms the quality of the DESI spectra and the performance of the pipelines developed for the data reduction and analysis of DESI data. △ Less

Submitted 27 October, 2023; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 14 pages, 4 figures, accepted for publication in ApJ, data available from https://doi.org/10.5281/zenodo.8363303

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

Algorithm to Produce a Density Field with Given Two, Three, and Four-Point Correlation Functions

Authors: Zachary Slepian

Abstract: Here we show how to produce a 3D density field with a given set of higher-order correlation functions. Our algorithm enables producing any desired two-point, three-point, and four-point functions, including odd-parity for the latter. We note that this algorithm produces the desired correlations about a set of "primary" points, matched to how the spherical-harmonic-based algorithms ENCORE and CADEN… ▽ More Here we show how to produce a 3D density field with a given set of higher-order correlation functions. Our algorithm enables producing any desired two-point, three-point, and four-point functions, including odd-parity for the latter. We note that this algorithm produces the desired correlations about a set of "primary" points, matched to how the spherical-harmonic-based algorithms ENCORE and CADENZA measure them. These "primary points" must be used as those around which the correlation functions are measured. We also outline how to generalize the algorithm to i) N-point correlations with $N>4$, ii) dimensions other than 3, and iii) beyond scalar quantities. This algorithm should find use in verifying analysis pipelines for higher-order statistics in upcoming galaxy redshift surveys such as DESI, Euclid, Roman, and Spherex, as well as intensity mapping △ Less

Submitted 30 March, 2024; v1 submitted 29 May, 2023; originally announced June 2023.

Comments: 10 pages, accepted Royal Astronomical Society Techniques and Instruments (RASTI)

SARABANDE: 3/4 Point Correlation Functions with Fast Fourier Transforms

Authors: James Sunseri, Zachary Slepian, Stephen Portillo, Jiamin Hou, Sule Kahraman, Douglas P. Finkbeiner

Abstract: We present a new $\texttt{python}$ package SARABANDE for measuring 3 & 4 Point Correlation Functions (3/4 PCFs) in $\mathcal{O}(N_{\rm g} \log N_{\rm g})$ time using Fast Fourier Transforms (FFTs), with $N_{\rm g}$ the number of grid points used for the FFT. SARABANDE can measure both projected and full 3 and 4 PCFs on gridded 2D and 3D datasets. The general technique is to generate suitable angul… ▽ More We present a new $\texttt{python}$ package SARABANDE for measuring 3 & 4 Point Correlation Functions (3/4 PCFs) in $\mathcal{O}(N_{\rm g} \log N_{\rm g})$ time using Fast Fourier Transforms (FFTs), with $N_{\rm g}$ the number of grid points used for the FFT. SARABANDE can measure both projected and full 3 and 4 PCFs on gridded 2D and 3D datasets. The general technique is to generate suitable angular basis functions on an underlying grid, radially bin these to create kernels, and convolve these kernels with the original gridded data to obtain expansion coefficients about every point simultaneously. These coefficients are then combined to give us the 3/4 PCF as expanded in our basis. We apply SARABANDE to simulations of the Interstellar Medium (ISM) to show the results and scaling of calculating both the full and projected 3/4 PCFs. △ Less

Submitted 25 October, 2022; v1 submitted 18 October, 2022; originally announced October 2022.

Comments: 16 Pages, 8 Figures, 8 Algorithms, 1 code package

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

Measurement of Parity-Odd Modes in the Large-Scale 4-Point Correlation Function of SDSS BOSS DR12 CMASS and LOWZ Galaxies

Authors: Jiamin Hou, Zachary Slepian, Robert N. Cahn

Abstract: A tetrahedron is the simplest shape that cannot be rotated into its mirror image in 3D. The 4-Point Correlation Function (4PCF), which quantifies excess clustering of quartets of galaxies over random, is the lowest-order statistic sensitive to parity violation. Each galaxy defines one vertex of the tetrahedron. Parity-odd modes of the 4PCF probe an imbalance between tetrahedra and their mirror ima… ▽ More A tetrahedron is the simplest shape that cannot be rotated into its mirror image in 3D. The 4-Point Correlation Function (4PCF), which quantifies excess clustering of quartets of galaxies over random, is the lowest-order statistic sensitive to parity violation. Each galaxy defines one vertex of the tetrahedron. Parity-odd modes of the 4PCF probe an imbalance between tetrahedra and their mirror images. We measure these modes from the largest currently available spectroscopic samples, the 280,067 Luminous Red Galaxies (LRGs) of the Baryon Oscillation Spectroscopic Survey (BOSS) DR12 LOWZ ($\bar{z} = 0.32$) and the 803,112 LRGS of BOSS DR12 CMASS ($\bar{z} = 0.57$). In LOWZ we find $3.1σ$ evidence for a non-zero parity-odd 4PCF, and in CMASS we detect a parity-odd 4PCF at $7.1σ$. Gravitational evolution alone does not produce this effect; parity-breaking in LSS, if cosmological in origin, must stem from the epoch of inflation. We have explored many sources of systematic error and found none that can produce a spurious parity-odd \textit{signal} sufficient to explain our result. Underestimation of the \textit{noise} could also lead to a spurious detection. Our reported significances presume that the mock catalogs used to calculate the covariance sufficiently capture the covariance of the true data. We have performed numerous tests to explore this issue. The odd-parity 4PCF opens a new avenue for probing new forces during the epoch of inflation with 3D LSS; such exploration is timely given large upcoming spectroscopic samples such as DESI and Euclid. △ Less

Submitted 23 June, 2023; v1 submitted 7 June, 2022; originally announced June 2022.

Comments: 45 pages, 38 figures, accepted by MNRAS

A Two-Fluid Treatment of the Effect of Neutrinos on the Matter Density

Authors: Farshad Kamalinejad, Zachary Slepian

Abstract: We study the linear effects of massive neutrinos on the matter density fluctuations on both large and small scales by solving the fluid equations that are obtained from Boltzmann hierarchy. These equations set the evolution of density and velocity perturbations of cold dark matter and neutrinos after the latter become non-relativistic. Due to the free-streaming effects, the massive neutrinos suppr… ▽ More We study the linear effects of massive neutrinos on the matter density fluctuations on both large and small scales by solving the fluid equations that are obtained from Boltzmann hierarchy. These equations set the evolution of density and velocity perturbations of cold dark matter and neutrinos after the latter become non-relativistic. Due to the free-streaming effects, the massive neutrinos suppress the total matter power spectrum on small scales while leaving the large scales untouched. Approximating the neutrino distribution function as peaked at one value gives an equation of state that can then be used in an analytic solution of the fluid equations. We then test the validity of fluid approximation with full Boltzmann solvers such as \textsc{class}. We also comment on the effects of dark energy on our results. △ Less

Submitted 24 March, 2022; originally announced March 2022.

On a General Method for Resolving Integrals of Multiple Spherical Bessel Functions Against Power Laws into Distributions

Authors: Kiersten Meigs, Zachary Slepian

Abstract: We here present a method of performing integrals of products of spherical Bessel functions (SBFs) weighted by a power-law. Our method, which begins with double-SBF integrals, exploits a differential operator $\hat{D}$ defined via Bessel's differential equation. Application of this operator raises the power-law in steps of two. We also here display a suitable base integral expression to which this… ▽ More We here present a method of performing integrals of products of spherical Bessel functions (SBFs) weighted by a power-law. Our method, which begins with double-SBF integrals, exploits a differential operator $\hat{D}$ defined via Bessel's differential equation. Application of this operator raises the power-law in steps of two. We also here display a suitable base integral expression to which this operator can be applied for both even and odd cases. We test our method by showing that it reproduces previously-known solutions. Importantly, it also goes beyond them, offering solutions in terms of singular distributions, Heaviside functions, and Gauss's hypergeometric,$\;_2{\rm F}_1$ for $all$ double-SBF integrals with positive semi-definite integer power-law weight. We then show how our method for double-SBF integrals enables evaluating $arbitrary$ triple-SBF overlap integrals, going beyond the cases currently in the literature. This in turn enables reduction of arbitrary quadruple, quintuple, and sextuple-SBF integrals and beyond into tractable forms. △ Less

Submitted 14 December, 2021; originally announced December 2021.

Comments: 20 pages, 1 figure, submitted

Accelerating BAO Scale Fitting Using Taylor Series

Authors: Matthew Hansen, Alex Krolewski, Zachary Slepian

Abstract: The Universe is currently undergoing accelerated expansion driven by dark energy. Dark energy's essential nature remains mysterious: one means of revealing it is by measuring the Universe's size at different redshifts. This may be done using the Baryon Acoustic Oscillation (BAO) feature, a standard ruler in the galaxy 2-Point Correlation Function (2PCF). In order to measure the distance scale, one… ▽ More The Universe is currently undergoing accelerated expansion driven by dark energy. Dark energy's essential nature remains mysterious: one means of revealing it is by measuring the Universe's size at different redshifts. This may be done using the Baryon Acoustic Oscillation (BAO) feature, a standard ruler in the galaxy 2-Point Correlation Function (2PCF). In order to measure the distance scale, one dilates and contracts a template for the 2PCF in a fiducial cosmology, using a scaling factor $α$. The standard method for finding the best-fit $α$ is to compute the likelihood over a grid of roughly 100 values of it. This approach is slow; in this work, we propose a significantly faster way. Our method writes the 2PCF as a polynomial in $α$ by Taylor-expanding it about $α= 1$, exploiting that we know the fiducial cosmology sufficiently well that $α$ is within a few percent of unity. The likelihood resulting from this expansion may then be analytically solved for the best-fit $α$. Our method is 48-85$\times$ faster than a directly comparable approach in which we numerically minimize $α$, and $\sim$$12,000 \times$ faster than the standard iterative method. Our work will be highly enabling for upcoming large-scale structure redshift surveys such as that by Dark Energy Spectroscopic Instrument (DESI). △ Less

Submitted 29 November, 2022; v1 submitted 13 December, 2021; originally announced December 2021.

Comments: 13 pages, 6 figures, accepted MNRAS

A Test for Cosmological Parity Violation Using the 3D Distribution of Galaxies

Authors: Robert N. Cahn, Zachary Slepian, Jiamin Hou

Abstract: We show that the galaxy 4-Point Correlation Function (4PCF) can test for cosmological parity violation. The detection of cosmological parity violation would reflect previously unknown forces present at the earliest moments of the Universe. Recent developments both in rapidly evaluating galaxy $N$-Point Correlation Functions (NPCFs) and in determining the corresponding covariance matrices make the… ▽ More We show that the galaxy 4-Point Correlation Function (4PCF) can test for cosmological parity violation. The detection of cosmological parity violation would reflect previously unknown forces present at the earliest moments of the Universe. Recent developments both in rapidly evaluating galaxy $N$-Point Correlation Functions (NPCFs) and in determining the corresponding covariance matrices make the search for parity violation in the 4PCF possible in current and upcoming surveys such as those undertaken by Dark Energy Spectroscopic Instrument (DESI), the $Euclid$ satellite, and the Vera C. Rubin Observatory (VRO). △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: 5 pages

Analytic Gaussian Covariance Matrices for Galaxy $N$-Point Correlation Functions

Authors: Jiamin Hou, Robert N. Cahn, Oliver H. E. Philcox, Zachary Slepian

Abstract: We derive analytic covariance matrices for the $N$-Point Correlation Functions (NPCFs) of galaxies in the Gaussian limit. Our results are given for arbitrary $N$ and projected onto the isotropic basis functions of Cahn & Slepian (2020), recently shown to facilitate efficient NPCF estimation. A numerical implementation of the 4PCF covariance is compared to the sample covariance obtained from a set… ▽ More We derive analytic covariance matrices for the $N$-Point Correlation Functions (NPCFs) of galaxies in the Gaussian limit. Our results are given for arbitrary $N$ and projected onto the isotropic basis functions of Cahn & Slepian (2020), recently shown to facilitate efficient NPCF estimation. A numerical implementation of the 4PCF covariance is compared to the sample covariance obtained from a set of lognormal simulations, Quijote dark matter halo catalogues, and MultiDark-Patchy galaxy mocks, with the latter including realistic survey geometry. The analytic formalism gives reasonable predictions for the covariances estimated from mock simulations with a periodic-box geometry. Furthermore, fitting for an effective volume and number density by maximizing a likelihood based on Kullback-Leibler divergence is shown to partially compensate for the effects of a non-uniform window function. △ Less

Submitted 3 August, 2021; originally announced August 2021.

Comments: 31 pages, 17 figures

A First Detection of the Connected 4-Point Correlation Function of Galaxies Using the BOSS CMASS Sample

Authors: Oliver H. E. Philcox, Jiamin Hou, Zachary Slepian

Abstract: We present an $8.1σ$ detection of the non-Gaussian 4-Point Correlation Function (4PCF) using a sample of $N_{\rm g} \approx 8\times 10^5$ galaxies from the BOSS CMASS dataset. Our measurement uses the $\mathcal{O}(N_{\rm g}^2)$ NPCF estimator of Philcox et al. (2021), including a new modification to subtract the disconnected 4PCF contribution (arising from the product of two 2PCFs) at the estimato… ▽ More We present an $8.1σ$ detection of the non-Gaussian 4-Point Correlation Function (4PCF) using a sample of $N_{\rm g} \approx 8\times 10^5$ galaxies from the BOSS CMASS dataset. Our measurement uses the $\mathcal{O}(N_{\rm g}^2)$ NPCF estimator of Philcox et al. (2021), including a new modification to subtract the disconnected 4PCF contribution (arising from the product of two 2PCFs) at the estimator level. This approach is unlike previous work and ensures that our signal is a robust detection of gravitationally-induced non-Gaussianity. The estimator is validated with a suite of lognormal simulations, and the analytic form of the disconnected contribution is discussed. Due to the high dimensionality of the 4PCF, data compression is required; we use a signal-to-noise-based scheme calibrated from theoretical covariance matrices to restrict to $\sim$ $100$ basis vectors. The compression has minimal impact on the detection significance and facilitates traditional $χ^2$-like analyses using a suite of mock catalogs. The significance is stable with respect to different treatments of noise in the sample covariance (arising from the limited number of mocks), but decreases to $4.7σ$ when a minimum galaxy separation of $14 h^{-1}\mathrm{Mpc}$ is enforced on the 4PCF tetrahedra (such that the statistic can be modelled more easily). The detectability of the 4PCF in the quasi-linear regime implies that it will become a useful tool in constraining cosmological and galaxy formation parameters from upcoming spectroscopic surveys. △ Less

Submitted 3 August, 2021; originally announced August 2021.

Comments: 26 pages, 9 figures, submitted to Phys. Rev. D. Code available at https://github.com/oliverphilcox/Parity-Even-4PCF

Clustering in Massive Neutrino Cosmologies via Eulerian Perturbation Theory

Authors: Alejandro Aviles, Arka Banerjee, Gustavo Niz, Zachary Slepian

Abstract: We introduce an Eulerian Perturbation Theory to study the clustering of tracers for cosmologies in the presence of massive neutrinos. Our approach is based on mapping recently-obtained Lagrangian Perturbation Theory results to the Eulerian framework. We add Effective Field Theory counterterms, IR-resummations and a biasing scheme to compute the one-loop redshift-space power spectrum. To assess our… ▽ More We introduce an Eulerian Perturbation Theory to study the clustering of tracers for cosmologies in the presence of massive neutrinos. Our approach is based on mapping recently-obtained Lagrangian Perturbation Theory results to the Eulerian framework. We add Effective Field Theory counterterms, IR-resummations and a biasing scheme to compute the one-loop redshift-space power spectrum. To assess our predictions, we compare the power spectrum multipoles against synthetic halo catalogues from the Quijote simulations, finding excellent agreement on scales $k\lesssim 0.25 \,h \text{Mpc}^{-1}$. One can obtain the same fitting accuracy using higher wave-numbers, but then the theory fails to give a correct estimation of the linear bias parameter. We further discuss the implications for the tree-level bispectrum. Finally, calculating loop corrections is computationally costly, hence we derive an accurate approximation wherein we retain only the main features of the kernels, as produced by changes to the growth rate. As a result, we show how FFTLog methods can be used to further accelerate the loop computations with these reduced kernels. △ Less

Submitted 11 November, 2021; v1 submitted 25 June, 2021; originally announced June 2021.

Comments: 42 pages, 10 figures

Report number: FERMILAB-PUB-21-280-T

Journal ref: JCAP11(2021)028

Efficient Computation of $N$-point Correlation Functions in $D$ Dimensions

Authors: Oliver H. E. Philcox, Zachary Slepian

Abstract: We present efficient algorithms for computing the $N$-point correlation functions (NPCFs) of random fields in arbitrary $D$-dimensional homogeneous and isotropic spaces. Such statistics appear throughout the physical sciences, and provide a natural tool to describe stochastic processes. algorithms for computing the NPCF components have $\mathcal{O}(n^N)$ complexity (for a data set containing $n$ p… ▽ More We present efficient algorithms for computing the $N$-point correlation functions (NPCFs) of random fields in arbitrary $D$-dimensional homogeneous and isotropic spaces. Such statistics appear throughout the physical sciences, and provide a natural tool to describe stochastic processes. algorithms for computing the NPCF components have $\mathcal{O}(n^N)$ complexity (for a data set containing $n$ particles); their application is thus computationally infeasible unless $N$ is small. By projecting the statistic onto a suitably-defined angular basis, we show that the estimators can be written in a separable form, with complexity $\mathcal{O}(n^2)$, or $\mathcal{O}(n_{\rm g}\log n_{\rm g})$ if evaluated using a Fast Fourier Transform on a grid of size $n_{\rm g}$. Our decomposition is built upon the $D$-dimensional hyperspherical harmonics; these form a complete basis on the $(D-1)$-sphere and are intrinsically related to angular momentum operators. Concatenation of $(N-1)$ such harmonics gives states of definite combined angular momentum, forming a natural separable basis for the NPCF. As $N$ and $D$ grow, the number of basis components quickly becomes large, providing a practical limitation to this (and all other) approaches: however, the dimensionality is greatly reduced in the presence of symmetries; for example, isotropic correlation functions require only states of zero combined angular momentum. We provide a \textsc{Julia} package implementing our estimators, and show how they can be applied to a variety of scenarios within cosmology and fluid dynamics. The efficiency of such estimators will allow higher-order correlators to become a standard tool in the analysis of random fields. △ Less

Submitted 2 February, 2022; v1 submitted 18 June, 2021; originally announced June 2021.

Comments: 12 pages, 3 figures, accepted by PNAS. Code available at https://github.com/oliverphilcox/NPCFs.jl

ENCORE: An $\mathcal{O}(N_{\rm g}^2)$ Estimator for Galaxy $N$-Point Correlation Functions

Authors: Oliver H. E. Philcox, Zachary Slepian, Jiamin Hou, Craig Warner, Robert N. Cahn, Daniel J. Eisenstein

Abstract: We present a new algorithm for efficiently computing the $N$-point correlation functions (NPCFs) of a 3D density field for arbitrary $N$. This can be applied both to a discrete spectroscopic galaxy survey and a continuous field. By expanding the statistics in a separable basis of isotropic functions built from spherical harmonics, the NPCFs can be estimated by counting pairs of particles in space,… ▽ More We present a new algorithm for efficiently computing the $N$-point correlation functions (NPCFs) of a 3D density field for arbitrary $N$. This can be applied both to a discrete spectroscopic galaxy survey and a continuous field. By expanding the statistics in a separable basis of isotropic functions built from spherical harmonics, the NPCFs can be estimated by counting pairs of particles in space, leading to an algorithm with complexity $\mathcal{O}(N_{\rm g}^2)$ for $N_{\rm g}$ particles, or $\mathcal{O}\left(N_\mathrm{FFT}\log N_\mathrm{FFT}\right)$ when using a Fast Fourier Transform with $N_\mathrm{FFT}$ grid-points. In practice, the rate-limiting step for $N>3$ will often be the summation of the histogrammed spherical harmonic coefficients, particularly if the number of radial and angular bins is large. In this case, the algorithm scales linearly with $N_{\rm g}$. The approach is implemented in the ENCORE code, which can compute the 3PCF, 4PCF, 5PCF, and 6PCF of a BOSS-like galaxy survey in $\sim$ $100$ CPU-hours, including the corrections necessary for non-uniform survey geometries. We discuss the implementation in depth, along with its GPU acceleration, and provide practical demonstration on realistic galaxy catalogs. Our approach can be straightforwardly applied to current and future datasets to unlock the potential of constraining cosmology from the higher-point functions. △ Less

Submitted 13 October, 2021; v1 submitted 18 May, 2021; originally announced May 2021.

Comments: 25 pages, 6 figures, accepted by MNRAS. Code available at https://github.com/oliverphilcox/encore

An Exact Integral-to-Sum Relation for Products of Bessel Functions

Authors: Oliver H. E. Philcox, Zachary Slepian

Abstract: A useful identity relating the infinite sum of two Bessel functions to their infinite integral was discovered in Dominici et al. (2012). Here, we extend this result to products of $N$ Bessel functions, and show it can be straightforwardly proven using the Abel-Plana theorem, or the Poisson summation formula. For $N=2$, the proof is much simpler than that of Dominici et al., and significantly enlar… ▽ More A useful identity relating the infinite sum of two Bessel functions to their infinite integral was discovered in Dominici et al. (2012). Here, we extend this result to products of $N$ Bessel functions, and show it can be straightforwardly proven using the Abel-Plana theorem, or the Poisson summation formula. For $N=2$, the proof is much simpler than that of Dominici et al., and significantly enlarges the range of validity. △ Less

Submitted 5 May, 2021; v1 submitted 20 April, 2021; originally announced April 2021.

Comments: 15 pages, 1 figure. Submitted to Proc. Roy. Soc. A

Kepler's Goat Herd: An Exact Solution to Kepler's Equation for Elliptical Orbits

Authors: Oliver H. E. Philcox, Jeremy Goodman, Zachary Slepian

Abstract: A fundamental relation in celestial mechanics is Kepler's equation, linking an orbit's mean anomaly to its eccentric anomaly and eccentricity. Being transcendental, the equation cannot be directly solved for eccentric anomaly by conventional treatments; much work has been devoted to approximate methods. Here, we give an explicit integral solution, utilizing methods recently applied to the "geometr… ▽ More A fundamental relation in celestial mechanics is Kepler's equation, linking an orbit's mean anomaly to its eccentric anomaly and eccentricity. Being transcendental, the equation cannot be directly solved for eccentric anomaly by conventional treatments; much work has been devoted to approximate methods. Here, we give an explicit integral solution, utilizing methods recently applied to the "geometric goat problem" and to the dynamics of spherical collapse. The solution is given as a ratio of contour integrals; these can be efficiently computed via numerical integration for arbitrary eccentricities. The method is found to be highly accurate in practice, with our C++ implementation outperforming conventional root-finding and series approaches by a factor greater than two. △ Less

Submitted 13 May, 2021; v1 submitted 29 March, 2021; originally announced March 2021.

Comments: 7 pages, 4 figures, accepted by MNRAS. Code available at https://github.com/oliverphilcox/Keplers-Goat-Herd

A Uniform Spherical Goat (Problem): Explicit Solution for Homologous Collapse's Radial Evolution in Time

Authors: Zachary Slepian, Oliver H. E. Philcox

Abstract: The homologous collapse from rest of a uniform density sphere under its self gravity is a well-known toy model for the formation dynamics of astronomical objects ranging from stars to galaxies. Equally well-known is that the evolution of the radius with time cannot be explicitly obtained because of the transcendental nature of the differential equation solution. Rather, both radius and time are wr… ▽ More The homologous collapse from rest of a uniform density sphere under its self gravity is a well-known toy model for the formation dynamics of astronomical objects ranging from stars to galaxies. Equally well-known is that the evolution of the radius with time cannot be explicitly obtained because of the transcendental nature of the differential equation solution. Rather, both radius and time are written parametrically in terms of the development angle $θ$. We here present an explicit integral solution for radius as a function of time, exploiting methods from complex analysis recently applied to the mathematically-similar 'geometric goat problem'. Our solution can be efficiently evaluated using a Fast Fourier Transform and allows for arbitrary sampling in time, with a simple Python implementation that is $\sim$$100\times$ faster than using numerical root-finding to achieve arbitrary sampling. Our explicit solution is advantageous relative to the usual approach of first generating a uniform grid in $θ$, since this latter results in a non-uniform radial or time sampling, less useful for applications such as generation of sub-grid physics models. △ Less

Submitted 20 March, 2021; v1 submitted 17 March, 2021; originally announced March 2021.

Comments: 4 pages, 3 figures, submitted to MNRAS. Python implementation available at https://gist.github.com/oliverphilcox/559f086f1bf63b23d55c508b2f47bad3

Beyond Yamamoto: Anisotropic Power Spectra and Correlation Functions with Pairwise Lines-of-Sight

Authors: Oliver H. E. Philcox, Zachary Slepian

Abstract: Conventional estimators of the anisotropic power spectrum and two-point correlation function (2PCF) adopt the `Yamamoto approximation', fixing the line-of-sight of a pair of galaxies to that of just one of its members. Whilst this is accurate only to first-order in the characteristic opening angle $θ_\max$, it allows for efficient implementation via Fast Fourier Transforms (FFTs). This work presen… ▽ More Conventional estimators of the anisotropic power spectrum and two-point correlation function (2PCF) adopt the `Yamamoto approximation', fixing the line-of-sight of a pair of galaxies to that of just one of its members. Whilst this is accurate only to first-order in the characteristic opening angle $θ_\max$, it allows for efficient implementation via Fast Fourier Transforms (FFTs). This work presents practical algorithms for computing the power spectrum and 2PCF multipoles using pairwise lines-of-sight, adopting either the galaxy midpoint or angle bisector definitions. Using newly derived infinite series expansions for spherical harmonics and Legendre polynomials, we construct estimators accurate to arbitrary order in $θ_\max$, though note that the midpoint and bisector formalisms themselves differ at fourth order. Each estimator can be straightforwardly implemented using FFTs, requiring only modest additional computational cost relative to the Yamamoto approximation. We demonstrate the algorithms by applying them to a set of realistic mock galaxy catalogs, and find both procedures produce comparable results for the 2PCF, with a slight preference for the bisector power spectrum algorithm, albeit at the cost of greater memory usage. Such estimators provide a useful method to reduce wide-angle systematics for future surveys. △ Less

Submitted 30 April, 2021; v1 submitted 16 February, 2021; originally announced February 2021.

Comments: 18+9 pages, 5 figures, accepted by Phys. Rev. D. Code available at http://github.com/oliverphilcox/BeyondYamamoto

Journal ref: Phys. Rev. D 103, 123509 (2021)

Information Content of Higher-Order Galaxy Correlation Functions

Authors: Lado Samushia, Zachary Slepian, Francisco Villaescusa-Navarro

Abstract: The shapes of galaxy N-point correlation functions can be used as standard rulers to constrain the distance-redshift relationship and thence the expansion rate of the Universe. The cosmological density fields traced by late-time galaxy formation are initially nearly Gaussian, and hence all the cosmological information can be extracted from their 2-Point Correlation Function (2PCF) or its Fourier-s… ▽ More The shapes of galaxy N-point correlation functions can be used as standard rulers to constrain the distance-redshift relationship and thence the expansion rate of the Universe. The cosmological density fields traced by late-time galaxy formation are initially nearly Gaussian, and hence all the cosmological information can be extracted from their 2-Point Correlation Function (2PCF) or its Fourier-space analog the power spectrum. Subsequent nonlinear evolution under gravity, as well as halo and then galaxy formation, generate higher-order correlation functions. Since the mapping of the initial to the final density field is, on large scales, invertible, it is often claimed that the information content of the initial field's power spectrum is equal to that of all the higher-order functions of the final, nonlinear field. This claim implies that reconstruction of the initial density field from the nonlinear field renders analysis of higher-order correlation functions of the latter superfluous. We here show that this claim is false when the N-point functions are used as standard rulers. Constraints available from joint analysis of the galaxy power spectrum and bispectrum (Fourier-space analog of the 3-Point Correlation Function) can, in some cases, exceed those offered by the initial power spectrum even when the reconstruction is perfect. We provide a mathematical justification for this claim and also demonstrate it using a large suite of N-body simulations. In particular, we show that for the z = 0 real-space matter field in the limit of vanishing shot noise, taking modes up to k_max = 0.2 h/Mpc, using the bispectrum alone offers a factor of two reduction in the variance on the cosmic distance scale relative to that available from the power spectrum. △ Less

Submitted 2 February, 2021; originally announced February 2021.

Comments: 14 pages, 5 figures. Summitted to the Monthly Notices of the Royal Astronomical Society

Testing the theory of gravity with DESI: estimators, predictions and simulation requirements

Authors: Shadab Alam, Christian Arnold, Alejandro Aviles, Rachel Bean, Yan-Chuan Cai, Marius Cautun, Jorge L. Cervantes-Cota, Carolina Cuesta-Lazaro, N. Chandrachani Devi, Alexander Eggemeier, Sebastien Fromenteau, Alma X. Gonzalez-Morales, Vitali Halenka, Jian-hua He, Wojciech A. Hellwing, Cesar Hernandez-Aguayo, Mustapha Ishak, Kazuya Koyama, Baojiu Li, Axel de la Macorra, Jennifer Menesses Rizo, Christopher Miller, Eva-Maria Mueller, Gustavo Niz, Pierros Ntelis , et al. (11 additional authors not shown)

Abstract: Shortly after its discovery, General Relativity (GR) was applied to predict the behavior of our Universe on the largest scales, and later became the foundation of modern cosmology. Its validity has been verified on a range of scales and environments from the Solar system to merging black holes. However, experimental confirmations of GR on cosmological scales have so far lacked the accuracy one wou… ▽ More Shortly after its discovery, General Relativity (GR) was applied to predict the behavior of our Universe on the largest scales, and later became the foundation of modern cosmology. Its validity has been verified on a range of scales and environments from the Solar system to merging black holes. However, experimental confirmations of GR on cosmological scales have so far lacked the accuracy one would hope for -- its applications on those scales being largely based on extrapolation and its validity sometimes questioned in the shadow of the unexpected cosmic acceleration. Future astronomical instruments surveying the distribution and evolution of galaxies over substantial portions of the observable Universe, such as the Dark Energy Spectroscopic Instrument (DESI), will be able to measure the fingerprints of gravity and their statistical power will allow strong constraints on alternatives to GR. In this paper, based on a set of $N$-body simulations and mock galaxy catalogs, we study the predictions of a number of traditional and novel estimators beyond linear redshift distortions in two well-studied modified gravity models, chameleon $f(R)$ gravity and a braneworld model, and the potential of testing these deviations from GR using DESI. These estimators employ a wide array of statistical properties of the galaxy and the underlying dark matter field, including two-point and higher-order statistics, environmental dependence, redshift space distortions and weak lensing. We find that they hold promising power for testing GR to unprecedented precision. The major future challenge is to make realistic, simulation-based mock galaxy catalogs for both GR and alternative models to fully exploit the statistic power of the DESI survey and to better understand the impact of key systematic effects. Using these, we identify future simulation and analysis needs for gravity tests using DESI. △ Less

Submitted 8 October, 2021; v1 submitted 11 November, 2020; originally announced November 2020.

Comments: 69 pages, 36 figures, 6 tables; accepted for publication in JCAP. The modified gravity HOD galaxy catalogues used in this paper can be downloaded from http://icc.dur.ac.uk/data/, by searching for the title of this paper from the dataset table

Improving the Line of Sight for the Anisotropic 3-Point Correlation Function of Galaxies: Centroid and Unit-Vector-Average Methods Scaling as $\mathcal{O}(N^2)$

Authors: Karolina Garcia, Zachary Slepian

Abstract: The 3-Point Correlation Function (3PCF), which measures correlations between triplets of galaxies encodes information about peculiar velocities, which distort the observed positions of galaxies along the line of sight away from their true positions. To access this information, we must track the 3PCF's dependence not only on each triangle's shape, but also on its orientation with respect to the lin… ▽ More The 3-Point Correlation Function (3PCF), which measures correlations between triplets of galaxies encodes information about peculiar velocities, which distort the observed positions of galaxies along the line of sight away from their true positions. To access this information, we must track the 3PCF's dependence not only on each triangle's shape, but also on its orientation with respect to the line of sight. Consequently, different choices for the line of sight will affect the measured 3PCF. Up to now, the line of sight has been taken as the direction to a single triplet member (STM), but which triplet member is used impacts the 3PCF by ~20% of the statistical error for a BOSS-like survey. For DESI (2019-24), which is 5X more precise, this would translate to 100% of the statistical error, increasing the total error bar by 40%. We here propose a new method that is fully symmetric between the triplet members, and uses either the average of the three galaxy position vectors (which we show points to the triangle centroid), or the average of their unit (direction) vectors. Naively, these approaches would seem to require triplet counting, scaling as $N^3$, with $N$ the number of objects in the survey. By harnessing the solid harmonic shift theorem, we here show how these methods can be evaluated scaling as $N^2$. We expect that they can be used to make a robust, systematics-free measurement of the anisotropic 3PCF of upcoming redshift surveys such as DESI. So doing will in turn open an additional channel to constrain the growth rate of structure and thereby learn the matter density as well as test the theory of gravity. △ Less

Submitted 6 November, 2020; originally announced November 2020.

Comments: 16 pages, submitted MNRAS

A Non-Degenerate Neutrino Mass Signature in the Galaxy Bispectrum

Authors: Farshad Kamalinejad, Zachary Slepian

Abstract: In the Standard Model, neutrinos are massless, yet oscillation experiments show in fact they do have a small mass. Currently only the differences of the masses' squares are known, and an upper bound on the sum. However, upcoming surveys of the Universe's large-scale structure (LSS) can probe the neutrino mass by exposing how neutrinos modulate galaxy clustering. But these measurements are challeng… ▽ More In the Standard Model, neutrinos are massless, yet oscillation experiments show in fact they do have a small mass. Currently only the differences of the masses' squares are known, and an upper bound on the sum. However, upcoming surveys of the Universe's large-scale structure (LSS) can probe the neutrino mass by exposing how neutrinos modulate galaxy clustering. But these measurements are challenging: in looking at the clustering of galaxy pairs, the effect of neutrinos is degenerate with galaxy formation, the details of which are unknown. Marginalizing over them degrades the constraints. Here we show that using correlations of galaxy triplets---the 3-Point Correlation Function or its Fourier-space analog the bispectrum---can break the degeneracy between galaxy formation physics (known as biasing) and the neutrino mass. Specifically, we find a signature of neutrinos in the bispectrum's dipole moment (with respect to triangle opening angle) that is roughly orthogonal to the contribution of galaxy biases. This signature was missed in previous works by failing to account for how neutrinos alter mode-coupling between perturbations on different scales. Our proposed signature will contribute to upcoming LSS surveys' such as DESI making a robust detection of the neutrino mass. We estimate that it can offer several-$σ$ evidence for non-zero $m_ν$ with DESI from the bispectrum alone, and that this is independent from information in the galaxy power spectrum. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Comments: 9 pages including references, submitted

Isotropic N-Point Basis Functions and Their Properties

Authors: Robert N. Cahn, Zachary Slepian

Abstract: Isotropic functions of positions $\hat{\bf r}_1, \hat{\bf r}_2,\ldots, \hat{\bf r}_N$, i.e. functions invariant under simultaneous rotations of all the coordinates, are conveniently formed using spherical harmonics and Clebsch-Gordan coefficients. An orthonormal basis of such functions provides a formalism suitable for analyzing isotropic distributions such as those that arise in cosmology, for in… ▽ More Isotropic functions of positions $\hat{\bf r}_1, \hat{\bf r}_2,\ldots, \hat{\bf r}_N$, i.e. functions invariant under simultaneous rotations of all the coordinates, are conveniently formed using spherical harmonics and Clebsch-Gordan coefficients. An orthonormal basis of such functions provides a formalism suitable for analyzing isotropic distributions such as those that arise in cosmology, for instance in the clustering of galaxies as revealed by large-scale structure surveys. The algebraic properties of the basis functions are conveniently expressed in terms of 6-$j$ and 9-$j$ symbols. The calculation of relations among the basis functions is facilitated by "Yutsis" diagrams for the addition and recoupling of angular momenta. △ Less

Submitted 27 October, 2020; originally announced October 2020.

Comments: 57 pages, submitted to Journal of Mathematical Physics, comments welcome

Classification of Magnetohydrodynamic Simulations using Wavelet Scattering Transforms

Authors: Andrew K. Saydjari, Stephen K. N. Portillo, Zachary Slepian, Sule Kahraman, Blakesley Burkhart, Douglas P. Finkbeiner

Abstract: The complex interplay of magnetohydrodynamics, gravity, and supersonic turbulence in the interstellar medium (ISM) introduces non-Gaussian structure that can complicate comparison between theory and observation. We show that the Wavelet Scattering Transform (WST), in combination with linear discriminant analysis (LDA), is sensitive to non-Gaussian structure in 2D ISM dust maps. WST-LDA classifies… ▽ More The complex interplay of magnetohydrodynamics, gravity, and supersonic turbulence in the interstellar medium (ISM) introduces non-Gaussian structure that can complicate comparison between theory and observation. We show that the Wavelet Scattering Transform (WST), in combination with linear discriminant analysis (LDA), is sensitive to non-Gaussian structure in 2D ISM dust maps. WST-LDA classifies magnetohydrodynamic (MHD) turbulence simulations with up to a 97\% true positive rate in our testbed of 8 simulations with varying sonic and Alfvénic Mach numbers. We present a side-by-side comparison with two other methods for non-Gaussian characterization, the Reduced Wavelet Scattering Transform (RWST) and the 3-Point Correlation Function (3PCF). We also demonstrate the 3D-WST-LDA and apply it to classification of density fields in position-position-velocity (PPV) space, where density correlations can be studied using velocity coherence as a proxy. WST-LDA is robust to common observational artifacts, such as striping and missing data, while also sensitive enough to extract the net magnetic field direction for sub-Alfvénic turbulent density fields. We include a brief analysis of the effect of point spread functions and image pixelization on 2D-WST-LDA applied to density fields, which informs the future goal of applying WST-LDA to 2D or 3D all-sky dust maps to extract hydrodynamic parameters of interest. △ Less

Submitted 22 October, 2020; originally announced October 2020.

Comments: 21 pages, 13 figures, submitted to ApJ

The Catalogue for Astrophysical Turbulence Simulations (CATS)

Authors: B. Burkhart, S. Appel, S. Bialy, J. Cho, A. J. Christensen, D. Collins, C. Federrath, D. Fielding, D. Finkbeiner, A. S. Hill, J. C. Ibanez-Mejia, M. R. Krumholz, A. Lazarian, M. Li, P. Mocz, M. -M. Mac Low, J. Naiman, S. K. N. Portillo, B. Shane, Z. Slepian, Y. Yuan

Abstract: Turbulence is a key process in many fields of astrophysics. Advances in numerical simulations of fluids over the last several decades have revolutionized our understanding of turbulence and related processes such as star formation and cosmic ray propagation. However, data from numerical simulations of astrophysical turbulence are often not made public. We introduce a new simulation-oriented databa… ▽ More Turbulence is a key process in many fields of astrophysics. Advances in numerical simulations of fluids over the last several decades have revolutionized our understanding of turbulence and related processes such as star formation and cosmic ray propagation. However, data from numerical simulations of astrophysical turbulence are often not made public. We introduce a new simulation-oriented database for the astronomical community: The Catalogue for Astrophysical Turbulence Simulations (CATS), located at www.mhdturbulence.com. CATS includes magnetohydrodynamic (MHD) turbulent box simulation data products generated by the public codes athena++, arepo, enzo, and flash. CATS also includes several synthetic observational data sets, such as turbulent HI data cubes. We also include measured power spectra and 3-point correlation functions from some of these data. We discuss the importance of open source statistical and visualization tools for the analysis of turbulence simulations such as those found in CATS. △ Less

Submitted 21 October, 2020; originally announced October 2020.

Comments: Accepted to ApJ

Rotation method for accelerating multiple-spherical Bessel function integrals against a numerical source function

Authors: Zachary Slepian, Yin Li, Marcel Schmittfull, Zvonimir Vlah

Abstract: A common problem in cosmology is to integrate the product of two or more spherical Bessel functions (sBFs) with different configuration-space arguments against the power spectrum or its square, weighted by powers of wavenumber. Naively computing them scales as $N_{\rm g}^{p+1}$ with $p$ the number of configuration space arguments and $N_{\rm g}$ the grid size, and they cannot be done with Fast Fou… ▽ More A common problem in cosmology is to integrate the product of two or more spherical Bessel functions (sBFs) with different configuration-space arguments against the power spectrum or its square, weighted by powers of wavenumber. Naively computing them scales as $N_{\rm g}^{p+1}$ with $p$ the number of configuration space arguments and $N_{\rm g}$ the grid size, and they cannot be done with Fast Fourier Transforms (FFTs). Here we show that by rewriting the sBFs as sums of products of sine and cosine and then using the product to sum identities, these integrals can then be performed using 1-D FFTs with $N_{\rm g} \log N_{\rm g}$ scaling. This "rotation" method has the potential to accelerate significantly a number of calculations in cosmology, such as perturbation theory predictions of loop integrals, higher order correlation functions, and analytic templates for correlation function covariance matrices. We implement this approach numerically both in a free-standing, publicly-available \textsc{Python} code and within the larger, publicly-available package \texttt{mcfit}. The rotation method evaluated with direct integrations already offers a factor of 6-10$\times$ speed-up over the naive approach in our test cases. Using FFTs, which the rotation method enables, then further improves this to a speed-up of $\sim$$1000-3000\times$ over the naive approach. The rotation method should be useful in light of upcoming large datasets such as DESI or LSST. In analysing these datasets recomputation of these integrals a substantial number of times, for instance to update perturbation theory predictions or covariance matrices as the input linear power spectrum is changed, will be one piece in a Monte Carlo Markov Chain cosmological parameter search: thus the overall savings from our method should be significant. △ Less

Submitted 29 November, 2019; originally announced December 2019.

Comments: 21 pages, 14 figures, submitted MNRAS; code available at https://github.com/eelregit/sbf_rotation

Accelerating Computation of the Nonlinear Mass by an Order of Magnitude

Authors: Alex Krolewski, Zachary Slepian

Abstract: The nonlinear mass is a characteristic scale in halo formation that has wide-ranging applications across cosmology. Naively, computing it requires repeated numerical integration to calculate the variance of the power spectrum on different scales and determine which scales exceed the threshold for nonlinear collapse. We accelerate this calculation by working in configuration space and approximating… ▽ More The nonlinear mass is a characteristic scale in halo formation that has wide-ranging applications across cosmology. Naively, computing it requires repeated numerical integration to calculate the variance of the power spectrum on different scales and determine which scales exceed the threshold for nonlinear collapse. We accelerate this calculation by working in configuration space and approximating the correlation function as a polynomial at $r <= 5$ $h^{-1}$ Mpc. This enables an analytic rather than numerical solution, accurate across a variety of cosmologies to 0.1$-$1% (depending on redshift) and 10$-$20 times faster than the naive numerical method. We also present a further acceleration (40$-$80 times faster than the naive method) in which we determine the polynomial coefficients using a Taylor expansion in the cosmological parameters rather than re-fitting a polynomial to the correlation function. Our acceleration greatly reduces the cost of repeated calculation of the nonlinear mass. This will be useful for MCMC analyses to constrain cosmological parameters from the highly nonlinear regime, e.g. with data from upcoming surveys. △ Less

Submitted 1 November, 2019; originally announced November 2019.

Comments: 8 pages, 5 figures. Comments welcome! Code available at https://github.com/akrolewski/NonlinearMassFaster

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.

Astrobites as a Community-led Model for Education, Science Communication, and Accessibility in Astrophysics

Authors: Gourav Khullar, Susanna Kohler, Tarini Konchady, Mike Foley, Amber L. Hornsby, Mithi A. de los Reyes, Nora Elisa Chisari, V. Ashley Villar, Kaitlyn Shin, Caitlin Doughty, Nora Shipp, Joanna Ramasawmy, Zephyr Penoyre, Tim Lichtenberg, Kate Storey-Fisher, Oliver Hall, Briley Lewis, Aaron B. Pearlman, Alejandro Cárdenas-Avendaño, Joanna S. Bridge, Elena González-Egea, Vatsal Panwar, Zachary Slepian, Mara Zimmerman

Abstract: Support for early career astronomers who are just beginning to explore astronomy research is imperative to increase retention of diverse practitioners in the field. Since 2010, Astrobites has played an instrumental role in engaging members of the community -- particularly undergraduate and graduate students -- in research. In this white paper, the Astrobites collaboration outlines our multi-facete… ▽ More Support for early career astronomers who are just beginning to explore astronomy research is imperative to increase retention of diverse practitioners in the field. Since 2010, Astrobites has played an instrumental role in engaging members of the community -- particularly undergraduate and graduate students -- in research. In this white paper, the Astrobites collaboration outlines our multi-faceted online education platform that both eases the transition into astronomy research and promotes inclusive professional development opportunities. We additionally offer recommendations for how the astronomy community can reduce barriers to entry to astronomy research in the coming decade. △ Less

Submitted 22 July, 2019; originally announced July 2019.

Comments: 14 pages, 5 figures (including cover page and list of endorsers); Submitted as an Astro2020 Decadal Survey State of the Profession White Paper. Access the Astrobites in the Classroom lesson plans here: https://tinyurl.com/y6kfszbn

Astro2020 Project White Paper: The Cosmic Accelerometer

Authors: Stephen S. Eikenberry, Anthony Gonzalez, Jeremy Darling, Jochen Liske, Zachary Slepian, Guido Mueller, John Conklin, Paul Fulda, Claudia Mendes de Oliveira, Misty Bentz, Sarik Jeram, Chenxing Dong, Amanda Townsend, Lilianne Mariko Izuti Nakazono, Robert Quimby, William Welsh, Joseph Harrington, Nicholas Law

Abstract: We propose an experiment, the Cosmic Accelerometer, designed to yield velocity precision of $\leq 1$ cm/s with measurement stability over years to decades. The first-phase Cosmic Accelerometer, which is at the scale of the Astro2020 Small programs, will be ideal for precision radial velocity measurements of terrestrial exoplanets in the Habitable Zone of Sun-like stars. At the same time, this expe… ▽ More We propose an experiment, the Cosmic Accelerometer, designed to yield velocity precision of $\leq 1$ cm/s with measurement stability over years to decades. The first-phase Cosmic Accelerometer, which is at the scale of the Astro2020 Small programs, will be ideal for precision radial velocity measurements of terrestrial exoplanets in the Habitable Zone of Sun-like stars. At the same time, this experiment will serve as the technical pathfinder and facility core for a second-phase larger facility at the Medium scale, which can provide a significant detection of cosmological redshift drift on a 6-year timescale. This larger facility will naturally provide further detection/study of Earth twin planet systems as part of its external calibration process. This experiment is fundamentally enabled by a novel low-cost telescope technology called PolyOculus, which harnesses recent advances in commercial off the shelf equipment (telescopes, CCD cameras, and control computers) combined with a novel optical architecture to produce telescope collecting areas equivalent to standard telescopes with large mirror diameters. Combining a PolyOculus array with an actively-stabilized high-precision radial velocity spectrograph provides a unique facility with novel calibration features to achieve the performance requirements for the Cosmic Accelerometer. △ Less

Submitted 18 July, 2019; originally announced July 2019.

Comments: 13 pages, 2 figures; Submitted as an Astro 2020 Decadal Survey Project White Paper

Astro2020 White Paper: A Direct Measure of Cosmic Acceleration

Authors: Stephen Eikenberry, Anthony Gonzalez, Jeremy Darling, Zachary Slepian, Guido Mueller, John Conklin, Paul Fulda, Sarik Jeram, Chenxing Dong, Amanda Townsend, Manunya Likamonsavad

Abstract: Nearly a century after the discovery that we live in an expanding Universe, and two decades after the discovery of accelerating cosmic expansion, there remains no direct detection of this acceleration via redshift drift - a change in the cosmological expansion velocity versus time. Because cosmological redshift drift directly determines the Hubble parameter H(z), it is arguably the cleanest possib… ▽ More Nearly a century after the discovery that we live in an expanding Universe, and two decades after the discovery of accelerating cosmic expansion, there remains no direct detection of this acceleration via redshift drift - a change in the cosmological expansion velocity versus time. Because cosmological redshift drift directly determines the Hubble parameter H(z), it is arguably the cleanest possible measurement of the expansion history, and has the potential to constrain dark energy models (e.g. Kim et al. 2015). The challenge is that the signal is small - the best observational constraint presently has an uncertainty several orders of magnitude larger than the expected signal (Darling 2012). Nonetheless, direct detection of redshift drift is becoming feasible, with upcoming facilities such as the ESO-ELT and SKA projecting possible detection within two to three decades. This timescale is uncomfortably long given the potential of this cosmological test. With dedicated experiments it should be possible to rapidly accelerate progress and detect redshift drift with only a five-year observational baseline. Such a facility would also be ideal for precision radial velocity measurements of exoplanets, which could be obtained as a byproduct of the ongoing calibration measurements for the experiment. △ Less

Submitted 30 March, 2019; originally announced April 2019.

Comments: White paper submitted to the Astro2020 Decadal Survey. 6 pages

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

Automatic Kalman-Filter-based Wavelet Shrinkage Denoising of 1D Stellar Spectra

Authors: Sankalp Gilda, Zachary Slepian

Abstract: We propose a non-parametric method to denoise 1D stellar spectra based on wavelet shrinkage followed by adaptive Kalman thresholding. Wavelet shrinkage denoising involves applying the Discrete Wavelet Transform (DWT) to the input signal, `shrinking' certain frequency components in the transform domain, and then applying inverse DWT to the reduced components. The performance of this procedure is in… ▽ More We propose a non-parametric method to denoise 1D stellar spectra based on wavelet shrinkage followed by adaptive Kalman thresholding. Wavelet shrinkage denoising involves applying the Discrete Wavelet Transform (DWT) to the input signal, `shrinking' certain frequency components in the transform domain, and then applying inverse DWT to the reduced components. The performance of this procedure is influenced by the choice of base wavelet, the number of decomposition levels, and the thresholding function. Typically, these parameters are chosen by `trial and error', which can be strongly dependent on the properties of the data being denoised. We here introduce an adaptive Kalman-filter-based thresholding method that eliminates the need for choosing the number of decomposition levels. We use the `Haar' wavelet basis, which we found to be the best-suited for 1D stellar spectra. We introduce various levels of Poisson noise into synthetic PHOENIX spectra, and test the performance of several common denoising methods against our own. It proves superior in terms of noise suppression and peak shape preservation. We expect it may also be of use in automatically and accurately filtering low signal-to-noise galaxy and quasar spectra obtained from surveys such as SDSS, Gaia, LSST, PESSTO, VANDELS, LEGA-C, and DESI. △ Less

Submitted 2 July, 2020; v1 submitted 12 March, 2019; originally announced March 2019.

Comments: 21 pages, 2 appendices, 17 figures. Published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society 490.4 (2019): 5249-5269

On Decoupling the Integrals of Cosmological Perturbation Theory

Authors: Zachary Slepian

Abstract: Perturbation theory (PT) is often used to model statistical observables capturing the translation and rotation-invariant information in cosmological density fields. PT produces higher-order corrections by integration over linear statistics of the density fields weighted by kernels resulting from recursive solution of the fluid equations. These integrals quickly become high-dimensional and naively… ▽ More Perturbation theory (PT) is often used to model statistical observables capturing the translation and rotation-invariant information in cosmological density fields. PT produces higher-order corrections by integration over linear statistics of the density fields weighted by kernels resulting from recursive solution of the fluid equations. These integrals quickly become high-dimensional and naively require increasing computational resources the higher the order of the corrections. Here we show how to decouple the integrands that often produce this issue, enabling PT corrections to be computed as a sum of products of independent 1-D integrals. Our approach is related to a commonly used method for calculating multi-loop Feynman integrals in Quantum Field Theory, the Gegenbauer Polynomial $x$-Space Technique (GPxT). We explicitly reduce the three terms entering the 2-loop power spectrum, formally requiring 9-D integrations, to sums over successive 1-D radial integrals. These 1-D integrals can further be performed as convolutions, rendering the scaling of this method $N_{\rm g} \log N_{\rm g}$ with $N_{\rm g}$ the number of grid points used for each Fast Fourier Transform. This method should be highly enabling for upcoming large-scale structure redshift surveys where model predictions at an enormous number of cosmological parameter combinations will be required by Monte Carlo Markov Chain searches for the best-fit values. △ Less

Submitted 6 December, 2018; originally announced December 2018.

Comments: 25 pages, 3 figures, submitted MNRAS, comments welcome

A Physical Picture of Bispectrum Baryon Acoustic Oscillations in the Interferometric Basis

Authors: Hillary L. Child, Zachary Slepian, Masahiro Takada

Abstract: We present a picture of the matter bispectrum in a novel "interferometric" basis designed to highlight interference of the baryon acoustic oscillations (BAO) in the power spectra composing it. Triangles where constructive interference amplifies BAO provide stronger cosmic distance constraints than triangles with destructive interference. We show that the amplitude of the BAO feature in the full cy… ▽ More We present a picture of the matter bispectrum in a novel "interferometric" basis designed to highlight interference of the baryon acoustic oscillations (BAO) in the power spectra composing it. Triangles where constructive interference amplifies BAO provide stronger cosmic distance constraints than triangles with destructive interference. We show that the amplitude of the BAO feature in the full cyclically summed bispectrum can be decomposed into simpler contributions from single terms or pairs of terms in the perturbation theory bispectrum, and that across large swathes of our parameter space the full BAO amplitude is described well by the amplitude of BAO in a single term. The dominant term is determined largely by the $F^{(2)}$ kernel of Eulerian standard perturbation theory. We present a simple physical picture of the BAO amplitude in each term; the BAO signal is strongest in triangle configurations where two wavenumbers differ by a multiple of the BAO fundamental wavelength. △ Less

Submitted 29 November, 2018; originally announced November 2018.

Comments: 21 pages, 15 figures

Bispectrum as Baryon Acoustic Oscillation Interferometer

Authors: Hillary L. Child, Masahiro Takada, Takahiro Nishimichi, Tomomi Sunayama, Zachary Slepian, Salman Habib, Katrin Heitmann

Abstract: The galaxy bispectrum, measuring excess clustering of galaxy triplets, offers a probe of dark energy via baryon acoustic oscillations (BAOs). However up to now it has been severely underused due to the combinatorically explosive number of triangles. Here we exploit interference in the bispectrum to identify triangles that amplify BAOs. This approach reduces the computational cost of estimating cov… ▽ More The galaxy bispectrum, measuring excess clustering of galaxy triplets, offers a probe of dark energy via baryon acoustic oscillations (BAOs). However up to now it has been severely underused due to the combinatorically explosive number of triangles. Here we exploit interference in the bispectrum to identify triangles that amplify BAOs. This approach reduces the computational cost of estimating covariance matrices, offers an improvement in BAO constraints equivalent to lengthening BOSS by 30%, and simplifies adding bispectrum BAO information to future large-scale redshift survey analyses. △ Less

Submitted 20 December, 2018; v1 submitted 28 June, 2018; originally announced June 2018.

Comments: 6 pages, 3 figures; revised to match published version

Journal ref: Phys. Rev. D 98, 123521 (2018)

Overview of the DESI Legacy Imaging Surveys

Authors: Arjun Dey, David J. Schlegel, Dustin Lang, Robert Blum, Kaylan Burleigh, Xiaohui Fan, Joseph R. Findlay, Doug Finkbeiner, David Herrera, Stephanie Juneau, Martin Landriau, Michael Levi, Ian McGreer, Aaron Meisner, Adam D. Myers, John Moustakas, Peter Nugent, Anna Patej, Edward F. Schlafly, Alistair R. Walker, Francisco Valdes, Benjamin A. Weaver, Christophe Yeche Hu Zou, Xu Zhou, Behzad Abareshi , et al. (135 additional authors not shown)

Abstract: The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerr… ▽ More The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project. △ Less

Submitted 19 February, 2019; v1 submitted 23 April, 2018; originally announced April 2018.

Comments: 47 pages, 18 figures; accepted for publication in the Astronomical Journal

nbodykit: an open-source, massively parallel toolkit for large-scale structure

Authors: Nick Hand, Yu Feng, Florian Beutler, Yin Li, Chirag Modi, Uros Seljak, Zachary Slepian

Abstract: We present nbodykit, an open-source, massively parallel Python toolkit for analyzing large-scale structure (LSS) data. Using Python bindings of the Message Passing Interface (MPI), we provide parallel implementations of many commonly used algorithms in LSS. nbodykit is both an interactive and scalable piece of scientific software, performing well in a supercomputing environment while still taking… ▽ More We present nbodykit, an open-source, massively parallel Python toolkit for analyzing large-scale structure (LSS) data. Using Python bindings of the Message Passing Interface (MPI), we provide parallel implementations of many commonly used algorithms in LSS. nbodykit is both an interactive and scalable piece of scientific software, performing well in a supercomputing environment while still taking advantage of the interactive tools provided by the Python ecosystem. Existing functionality includes estimators of the power spectrum, 2 and 3-point correlation functions, a Friends-of-Friends grouping algorithm, mock catalog creation via the halo occupation distribution technique, and approximate N-body simulations via the FastPM scheme. The package also provides a set of distributed data containers, insulated from the algorithms themselves, that enable nbodykit to provide a unified treatment of both simulation and observational data sets. nbodykit can be easily deployed in a high performance computing environment, overcoming some of the traditional difficulties of using Python on supercomputers. We provide performance benchmarks illustrating the scalability of the software. The modular, component-based approach of nbodykit allows researchers to easily build complex applications using its tools. The package is extensively documented at http://nbodykit.readthedocs.io, which also includes an interactive set of example recipes for new users to explore. As open-source software, we hope nbodykit provides a common framework for the community to use and develop in confronting the analysis challenges of future LSS surveys. △ Less

Submitted 15 December, 2017; originally announced December 2017.

Comments: 18 pages, 7 figures. Feedback very welcome. Code available at https://github.com/bccp/nbodykit and for documentation, see http://nbodykit.readthedocs.io

Developing the 3-Point Correlation Function For the Turbulent Interstellar Medium

Authors: Stephen K. N. Portillo, Zachary Slepian, Blakesley Burkhart, Sule Kahraman, Douglas P. Finkbeiner

Abstract: We present the first application of the angle-dependent 3-Point Correlation Function (3PCF) to the density fields magnetohydrodynamic (MHD) turbulence simulations intended to model interstellar (ISM) turbulence. Previous work has demonstrated that the angle-averaged bispectrum, the 3PCF's Fourier-space analog, is sensitive to the sonic and Alfvénic Mach numbers of turbulence. Here we show that int… ▽ More We present the first application of the angle-dependent 3-Point Correlation Function (3PCF) to the density fields magnetohydrodynamic (MHD) turbulence simulations intended to model interstellar (ISM) turbulence. Previous work has demonstrated that the angle-averaged bispectrum, the 3PCF's Fourier-space analog, is sensitive to the sonic and Alfvénic Mach numbers of turbulence. Here we show that introducing angular information via multipole moments with respect to the triangle opening angle offers considerable additional discriminatory power on these parameters. We exploit a fast, order $N_{\rm g} \log N_{\rm g}$ ($N_{\rm g}$ the number of grid cells used for a Fourier Transform) 3PCF algorithm to study a suite of MHD turbulence simulations with 10 different combinations of sonic and Alfvénic Mach numbers over a range from sub to super-sonic and sub to super-Alfvénic. The 3PCF algorithm's speed for the first time enables full quantification of the time-variation of our signal: we study 9 timeslices for each condition, demonstrating that the 3PCF is sufficiently time-stable to be used as an ISM diagnostic. In future, applying this framework to 3-D dust maps will enable better treatment of dust as a cosmological foreground as well as reveal conditions in the ISM that shape star formation. △ Less

Submitted 1 October, 2018; v1 submitted 27 November, 2017; originally announced November 2017.

Comments: 19 pages, 16 figures; version accepted by ApJ

Journal ref: The Astrophysical Journal, Volume 862, Issue 2, article id. 119, 15 pp. (2018)

Too hot to handle? Analytic solutions for massive neutrino or warm dark matter cosmologies

Authors: Zachary Slepian, Stephen KN Portillo

Abstract: We obtain novel closed form solutions to the Friedmann equation for cosmological models containing a component whose equation of state is that of radiation $(w=1/3)$ at early times and that of cold pressureless matter $(w=0)$ at late times. The equation of state smoothly transitions from the early to late-time behavior and exactly describes the evolution of a species with a Dirac Delta function di… ▽ More We obtain novel closed form solutions to the Friedmann equation for cosmological models containing a component whose equation of state is that of radiation $(w=1/3)$ at early times and that of cold pressureless matter $(w=0)$ at late times. The equation of state smoothly transitions from the early to late-time behavior and exactly describes the evolution of a species with a Dirac Delta function distribution in momentum magnitudes $|\vec{p}_0|$ (i.e. all particles have the same $|\vec{p}_0|$). Such a component, here termed "hot matter", is an approximate model for both neutrinos and warm dark matter. We consider it alone and in combination with cold matter and with radiation, also obtaining closed-form solutions for the growth of super-horizon perturbations in each case. The idealized model recovers $t(a)$ to better than $1.5\%$ accuracy for all $a$ relative to a Fermi-Dirac distribution (as describes neutrinos). We conclude by adding the second moment of the distribution to our exact solution and then generalizing to include all moments of an arbitrary momentum distribution in a closed form solution. △ Less

Submitted 4 October, 2017; originally announced October 2017.

Comments: 13 pages, 7 figures, MNRAS submitted

A Practical Computational Method for the Anisotropic Redshift-Space 3-Point Correlation Function

Authors: Zachary Slepian, Daniel J. Eisenstein

Abstract: We present an algorithm enabling computation of the anisotropic redshift-space galaxy 3-point correlation function (3PCF) scaling as $N^2$, with $N$ the number of galaxies. Our previous work showed how to compute the isotropic 3PCF with this scaling by expanding the radially-binned density field around each galaxy in the survey into spherical harmonics and combining these coefficients to form mult… ▽ More We present an algorithm enabling computation of the anisotropic redshift-space galaxy 3-point correlation function (3PCF) scaling as $N^2$, with $N$ the number of galaxies. Our previous work showed how to compute the isotropic 3PCF with this scaling by expanding the radially-binned density field around each galaxy in the survey into spherical harmonics and combining these coefficients to form multipole moments. The $N^2$ scaling occurred because this approach never explicitly required the relative angle between a galaxy pair about the primary galaxy. Here we generalize this work, demonstrating that in the presence of azimuthally-symmetric anisotropy produced by redshift-space distortions (RSD) the 3PCF can be described by two triangle side lengths, two independent total angular momenta, and a spin. This basis for the anisotropic 3PCF allows its computation with negligible additional work over the isotropic 3PCF. We also present the covariance matrix of the anisotropic 3PCF measured in this basis. Our algorithm tracks the full 5-D redshift-space 3PCF, uses an accurate line of sight to each triplet, is exact in angle, and easily handles edge correction. It will enable use of the anisotropic large-scale 3PCF as a probe of RSD in current and upcoming large-scale redshift surveys. △ Less

Submitted 28 September, 2017; originally announced September 2017.

Comments: 17 pages, 2 figures, MNRAS submitted

Galactos: Computing the Anisotropic 3-Point Correlation Function for 2 Billion Galaxies

Authors: Brian Friesen, Md. Mostofa Ali Patwary, Brian Austin, Nadathur Satish, Zachary Slepian, Narayanan Sundaram, Deborah Bard, Daniel J Eisenstein, Jack Deslippe, Pradeep Dubey, Prabhat

Abstract: The nature of dark energy and the complete theory of gravity are two central questions currently facing cosmology. A vital tool for addressing them is the 3-point correlation function (3PCF), which probes deviations from a spatially random distribution of galaxies. However, the 3PCF's formidable computational expense has prevented its application to astronomical surveys comprising millions to bill… ▽ More The nature of dark energy and the complete theory of gravity are two central questions currently facing cosmology. A vital tool for addressing them is the 3-point correlation function (3PCF), which probes deviations from a spatially random distribution of galaxies. However, the 3PCF's formidable computational expense has prevented its application to astronomical surveys comprising millions to billions of galaxies. We present Galactos, a high-performance implementation of a novel, O(N^2) algorithm that uses a load-balanced k-d tree and spherical harmonic expansions to compute the anisotropic 3PCF. Our implementation is optimized for the Intel Xeon Phi architecture, exploiting SIMD parallelism, instruction and thread concurrency, and significant L1 and L2 cache reuse, reaching 39% of peak performance on a single node. Galactos scales to the full Cori system, achieving 9.8PF (peak) and 5.06PF (sustained) across 9636 nodes, making the 3PCF easily computable for all galaxies in the observable universe. △ Less

Submitted 31 August, 2017; originally announced September 2017.

Comments: 11 pages, 7 figures, accepted to SuperComputing 2017