Source author record

Zachary Slepian

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

15 published item(s)

preprint2022arXiv

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

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

preprint2022arXiv

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

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.

preprint2022arXiv

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

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

preprint2021arXiv

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

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.

preprint2020arXiv

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

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.

preprint2015arXiv

A new look at lines of sight: using Fourier methods for the wide-angle anisotropic 2-point correlation function

The anisotropic 2-point correlation function (2PCF) of galaxies measures pairwise clustering as a function of the pair separation's angle to the line of sight. The latter is often defined as either the angle bisector of the observer-galaxy-pair triangle or the vector from the observer to the separation midpoint. Here we show how to accelerate either of these measurements with Fourier Transforms, using a slight generalization of the Yamamoto et al. (2006) estimator in which each member of the pair is used successively as the line of sight. We also present perturbation theory predictions for our generalized estimator including wide-angle corrections.

preprint2015arXiv

A Simple Analytic Treatment of Linear Growth of Structure with Baryon Acoustic Oscillations

In linear perturbation theory, all information about the growth of structure is contained in the Green's function, or equivalently, transfer function. These functions are generally computed using numerical codes or by phenomenological fitting formula anchored in accurate analytic results in the limits of large and small scale. Here we present a framework for analytically solving all scales, in particular the intermediate scales relevant for the baryon acoustic oscillations (BAO). We solve for the Green's function and transfer function using spherically-averaged overdensities and the approximation that the density of the coupled baryon-photon fluid is constant interior to the sound horizon.

preprint2015arXiv

Accelerating the 2-point and 3-point galaxy correlation functions using Fourier transforms

Though Fourier Transforms (FTs) are a common technique for finding correlation functions, they are not typically used in computations of the anisotropy of the two-point correlation function (2PCF) about the line of sight in wide-angle surveys because the line-of-sight direction is not constant on the Cartesian grid. Here we show how FTs can be used to compute the multipole moments of the anisotropic 2PCF. We also show how FTs can be used to accelerate the 3PCF algorithm of Slepian & Eisenstein (2015). In both cases, these FT methods allow one to avoid the computational cost of pair counting, which scales as the square of the number density of objects in the survey. With the upcoming large datasets of DESI, Euclid, and LSST, FT techniques will therefore offer an important complement to simple pair or triplet counts.

preprint2015arXiv

Computing the Three-Point Correlation Function of Galaxies in $\mathcal{O}(N^2)$ Time

We present an algorithm that computes the multipole coefficients of the galaxy three-point correlation function (3PCF) without explicitly considering triplets of galaxies. Rather, centering on each galaxy in the survey, it expands the radially-binned density field in spherical harmonics and combines these to form the multipoles without ever requiring the relative angle between a pair about the central. This approach scales with number and number density in the same way as the two-point correlation function, allowing runtimes that are comparable, and 500 times faster than a naive triplet count. It is exact in angle and easily handles edge correction. We demonstrate the algorithm on the LasDamas SDSS-DR7 mock catalogs, computing an edge corrected 3PCF out to $90\;{\rm Mpc}/h$ in under an hour on modest computing resources. We expect this algorithm will render it possible to obtain the large-scale 3PCF for upcoming surveys such as Euclid, LSST, and DESI.

preprint2015arXiv

On the signature of the baryon-dark matter relative velocity in the two and three-point galaxy correlation functions

We develop a configuration-space picture of the relative velocity between baryons and dark matter that clearly explains how it can shift the BAO scale in the galaxy-galaxy correlation function. The shift occurs because the relative velocity is non-zero only within the sound horizon and thus adds to the correlation function asymmetrically about the BAO peak. We further show that in configuration space the relative velocity has a localized, distinctive signature in the three-point galaxy correlation function (3PCF). In particular, we find that a multipole decomposition is a favorable way to isolate the relative velocity in the 3PCF, and that there is a strong signature in the l=1 multipole for triangles with 2 sides around the BAO scale. Finally, we investigate a further compression of the 3PCF to a function of only one triangle side that preserves the localized nature of the relative velocity signature while also nicely separating linear from non-linear bias. We expect that this scheme will substantially lessen the computational burden of finding the relative velocity in the 3PCF. The relative velocity's 3PCF signature can be used to correct the shift induced in the galaxy-galaxy correlation function so that no systematic error due to this effect is introduced into the BAO as used for precision cosmology.

preprint2014arXiv

A one-parameter formula for testing slow-roll dark energy: observational prospects

Numerous upcoming observations, such as WFIRST, BOSS, BigBOSS, LSST, Euclid, and Planck, will constrain dark energy (DE)'s equation of state with great precision. They may well find the ratio of pressure to energy density, $w$, is -1, meaning DE is equivalent to a cosmological constant. However, many time-varying DE models have also been proposed. A single parametrization to test a broad class of them and that is itself motivated by a physical picture is therefore desirable. We suggest the simplest model of DE has the same mechanism as inflation, likely a scalar field slowly rolling down its potential. If this is so, DE will have a generic equation of state and the Universe will have a generic dependence of the Hubble constant on redshift independent of the potential's starting value and shape. This equation of state and expression for the Hubble constant offer the desired model-independent but physically motivated parametrization, because they will hold for most of the standard scalar-field models of DE such as quintessence and phantom DE. Up until now two-parameter descriptions of $w$ have been available, but this work finds an additional approximation that leads to a single-parameter model. Using it, we conduct a $χ^2$ analysis and find that experiments in the next seven years should be able to distinguish any of these time-varying DE models on the one hand from a cosmological constant on the other to 73% confidence if $w$ today differs from -1 by 3.5%. In the limit of perfectly accurate measurements of $Ω_m$ and $H_0$, this confidence would rise to 96%. We also include discussion of the current status of DE experiment, a table compiling the techniques each will use, and tables of the precisions of the experiments for which this information was available at the time of publication.

preprint2012arXiv

Ruling Out Bosonic Repulsive Dark Matter in Thermal Equilibrium

Self-interacting dark matter (SIDM), especially bosonic, has been considered a promising candidate to replace cold dark matter (CDM) as it resolves some of the problems associated with CDM. Here, we rule out the possibility that dark matter is a repulsive boson in thermal equilibrium. We develop the model first proposed by Goodman (2000) and derive the equation of state at finite temperature. Isothermal spherical halo models indicate a Bose-Einstein condensed core surrounded by a non-degenerate envelope, with an abrupt density drop marking the boundary between the two phases. Comparing this feature with observed rotation curves constrains the interaction strength of our model's DM particle, and Bullet Cluster measurements constrain the scattering cross section. Both ultimately can be cast as constraints on the particle's mass. We find these two constraints cannot be satisfied simultaneously in any realistic halo model---and hence dark matter cannot be a repulsive boson in thermal equilibrium. It is still left open that DM may be a repulsive boson provided it is not in thermal equilibrium; this requires that the mass of the particle be significantly less than a millivolt.

preprint2012arXiv

The Average Projected Area Theorem - Generalization to Higher Dimensions

In 3-d the average projected area of a convex solid is 1/4 the surface area, as Cauchy showed in the 19th century. In general, the ratio in n dimensions may be obtained from Cauchy's surface area formula, which is in turn a special case of Kubota's theorem. However, while these latter results are well-known to those working in integral geometry or the theory of convex bodies, the results are largely unknown to the physics community---so much so that even the 3-d result is sometimes said to have first been proven by an astronomer in the early 20th century! This is likely because the standard proofs in the mathematical literature are, by and large, couched in terms of concepts that are may not be familiar to many physicists. Therefore, in this work, we present a simple geometrical method of calculating the ratio of average projected area to surface area for convex bodies in arbitrary dimensions. We focus on a pedagogical, physically intuitive treatment that it is hoped will be useful to those in the physics community. We do discuss the mathematical background of the theorem as well, pointing those who may be interested to sources that offer the proofs that are standard in the fields of integral geometry and the theory of convex bodies. We also provide discussion of the applications of the theorem, especially noting that higher-dimensional ratios may be of use for constructing observational tests of string theory. Finally, we examine the limiting behavior of the ratio with the goal of offering intuition on its behavior by pointing out a suggestive connection with a well-known fact in statistics.

preprint2011arXiv

Chance and Chandra (and repulsive dark matter)

A few examples are given of Chandra's work on statistical and stochastic problems that relate to open questions in astrophysics, in particular his theory of dynamical relaxation in systems with inverse-square interparticle forces. The roles of chaos and integrability in this theory require clarification, especially for systems having a dominant central mass. After this prelude, a hypothetical form of repulsive bosonic dark matter is discussed. The repulsion leads to nontrivial thermodynamic behavior, including superfluidity, and would tend to suppress dynamical friction, greatly reducing the drag exerted on rotating galactic bars. However, this form of dark matter can probably be ruled out, at least for parameters that allow halos to reach thermal equilibria within a Hubble time. One combination of the particle mass and interparticle repulsion determines the minimum core radius of dark halos. Bounds on dark-matter collisionality inferred from the Bullet Cluster constrain a second combination. It is possible to satisfy both constraints only for parameters that predict unacceptable rotation curves outside the halo core.

preprint2011arXiv

Dark Energy as Double N-Flation--Observational Predictions

We propose a simple model for dark energy useful for comparison with observations. It is based on the idea that dark energy and inflation should be caused by the same physical process. As motivation, we note that Linde's simple chaotic inflation $V=(1/2)m^{2}ϕ^{2}$ produces values of $n_{s}=0.967$ and $r=0.13$, which are consistent with the WMAP 1-sigma error bars. We therefore propose $V=(1/2)m_{2}^{2}ϕ_{2}^{2}+(1/2)m_{1}^{2}ϕ_{1}^{2}$ with $m_{1}\sim10^{-5}$ and $m_{2}\leq10^{-60}$, where $c=1=\hbar$ and the reduced Planck mass is set to unity. The field $ϕ_{1}$ drives inflation and has damped by now ($ϕ_{1,0}=0$), while $ϕ_{2}$ is currently rolling down its potential to produce dark energy. Using this model, we derive the formula $δw(z)\equiv w(z)+1=δw_{0}(H_0/H(z))^2$ via the slow-roll approximation. Our numerical results from exact and self-consistent solution of the equations of motion for $ϕ_2$ and the Friedmann equations support this formula, and it should hold for any slow-roll dark energy. Our potential can be easily realized in N-flation models with many fields, and is easily falsifiable by upcoming experiments -- for example, if Linde's chaotic inflation is ruled out. But if $r$ values consistent with Linde's chaotic inflation are detected then one should take this model seriously indeed.