BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Oliver H. E. Philcox

Oliver H. E. Philcox contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.COhep-phhep-thastro-ph.IMgr-qcastro-ph.GAphysics.data-anphysics.comp-phcond-mat.dis-nncond-mat.stat-mechhep-exmath-phmath.MPphysics.app-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
14topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

20 published item(s)

preprint2026arXiv

Reanalyzing DESI DR1: 1. $Λ$CDM Constraints from the Power Spectrum and Bispectrum

We present the first independent re-analysis of the galaxy clustering data from DESI Data Release 1, utilizing an effective field theory full-shape model. We analyze the power spectra and bispectra of the public catalogs using a custom-built pipeline based on window-deconvolved quasi-optimal estimators, accounting for a number of systematic effects. Compared to the official collaboration analysis, we add the galaxy power spectrum hexadecapole and the bispectrum monopole, and also introduce a novel stochastic estimator for fiber collisions, which facilitates robust bispectrum analyses. As a first application, we perform a full-shape analysis of the DESI power spectra and bispectra in the context of the standard cosmological model, $Λ$CDM. Using external priors on the physical baryon density and the primordial power spectrum tilt, we constrain the matter density fraction to $Ω_m=0.284\pm 0.011$, the Hubble constant to $H_0=70.7\pm 1.1$ km/s/Mpc, and the mass fluctuation amplitude to $σ_8=0.811\pm 0.030$. The bispectrum sharpens constraints on $σ_8$ and $Ω_m$ by $\approx 10\%$ and shifts $Ω_m$ by $\approx 1σ$ towards the \textit{Planck} $Λ$CDM value. Combining our full-shape likelihood with the official DESI DR2 BAO measurements, cosmological parameters shift further towards the \textit{Planck} values, with $Ω_m=0.296\pm 0.007$, $H_0=68.8\pm 0.6$ km/s/Mpc, $σ_8=0.818\pm 0.029$ (with tighter constraints obtained in joint analyses). Similar results are obtained in a joint analysis with DR1 BAO, accounting for the cross-covariance. Finally, the bispectrum data improves measurements of quadratic bias parameters, which are consistent with predictions from halo occupation distribution models. Our work highlights the importance of higher-order statistics and sets the stage for upcoming full-shape analyses of non-minimal cosmological models.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Squeezing $f_{\rm NL}$ out of the matter bispectrum with consistency relations

We show how consistency relations can be used to robustly extract the amplitude of local primordial non-Gaussianity ($f_{\rm NL}$) from the squeezed limit of the matter bispectrum, well into the non-linear regime. First, we derive a non-perturbative relation between primordial non-Gaussianity and the leading term in the squeezed bispectrum, revising some results present in the literature. This relation is then used to successfully measure $f_{\rm NL}$ from $N$-body simulations. We discuss the dependence of our results on different scale cuts and redshifts. Specifically, the analysis is strongly dependent on the choice of the smallest soft momentum, $q_{\rm min}$, which is the most sensitive to primordial bispectrum contributions, but is largely independent of the choice of the largest hard momentum, $k_{\rm max}$, due to the non-Gaussian nature of the covariance. We also show how the constraints on $f_{\rm NL}$ improve at higher redshift, due to a reduced off-diagonal covariance. In particular, for a simulation with $f_{\rm NL} = 100$ and a volume of $(2.4 \text{ Gpc}/h)^3$, we measure $f_{\rm NL} = 98 \pm 12$ at redshift $z=0$ and $f_{\rm NL} = 97 \pm 8$ at $z=0.97$. Finally, we compare our results with a Fisher forecast, showing that the current version of the analysis is satisfactorily close to the Fisher error. We regard this as a first step towards the realistic application of consistency relations to constrain primordial non-Gaussianity using observations.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

The Disordered Heterogeneous Universe: Galaxy Distribution and Clustering Across Length Scales

Studies of disordered heterogeneous media and galaxy cosmology share a common goal: analyzing the distribution of particles at `microscales' to predict physical properties at `macroscales', whether for a liquid, composite material, or entire Universe. The former theory provides an array of techniques to characterize a wide class of microstructures; in this work, we apply them to the distributions of galaxies. We focus on the lower-order correlation functions, `void' and `particle' nearest-neighbor functions, pair-connectedness functions, percolation properties, and a scalar order metric. Compared to homogeneous Poisson and typical disordered systems, the cosmological simulations exhibit enhanced large-scale clustering and longer tails in the nearest-neighbor functions, due to the presence of quasi-long-range correlations. On large scales, the system appears `hyperuniform', due to primordial density fluctuations, whilst on the smallest scales, the system becomes almost `antihyperuniform', and, via the order metric, is shown to be a highly correlated disordered system. Via a finite scaling analysis, we show that the percolation threshold of the galaxy catalogs is significantly lower than for Poisson realizations; this is consistent with the observation that the galaxy distribution contains larger voids. However, the two sets of simulations share a fractal dimension, implying that they lie in the same universality class. Finally, we consider the ability of large-scale clustering statistics to constrain cosmological parameters using simulation-based inference. Both the nearest-neighbor distribution and pair-connectedness function considerably tighten bounds on the amplitude of cosmological fluctuations at a level equivalent to observing twenty-five times more galaxies. These are a useful alternative to the three-particle correlation, and are computable in much reduced time. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Constraints on Multi-Field Inflation from the BOSS Galaxy Survey

We use redshift-space galaxy clustering data from the BOSS survey to constrain local primordial non-Gaussianity (LPNG). This is of particular importance due to the consistency relations, which imply that a detection of LPNG would rule out all single-field inflationary models. Our constraints are based on the consistently analyzed redshift-space galaxy power spectra and bispectra, extracted from the public BOSS data with optimal window-free estimators. We use a complete perturbation theory model including all one-loop power spectrum corrections generated by LPNG. Our constraint on the amplitude of the local non-Gaussian shape is $f_{\rm NL}^{\rm local}=-33\pm 28$ at 68\%\,CL, yielding no evidence for primordial non-Gaussianity. The addition of the bispectrum tightens the $f_{\rm NL}^{\rm local}$ constraints from BOSS by $20\%$, and allows breaking of degeneracies with non-Gaussian galaxy bias. These results set the stage for the analysis of future surveys, whose larger volumes will yield significantly tighter constraints on LPNG.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Constraints on Single-Field Inflation from the BOSS Galaxy Survey

Non-local primordial non-Gaussianity (NLPNG) is a smoking gun of interactions in single-field inflationary models, and can be written as a combination of the equilateral and orthogonal templates. We present the first constraints on these from the redshift-space galaxy power spectra and bispectra of the Baryon Oscillation Spectroscopic Survey (BOSS) data. These are the first such measurements independent of the cosmic microwave background fluctuations. We perform a consistent analysis that includes all necessary nonlinear corrections generated by NLPNG, and vary all relevant cosmological and nuisance parameters in a global fit to the data. Our conservative analysis yields joint limits on the amplitudes of the equilateral and orthogonal shapes, $f_{\rm NL}^{\rm equil}=940\pm 600$, $f_{\rm NL}^{\rm ortho}= -170\pm 170$ (both at 68\% CL). These can be used to derive constraints on coefficients of the effective single-field inflationary Lagrangian; in particular, we find that the sound speed of inflaton fluctuations has the bound $c_s\geq 0.013$ at 95\% CL. Fixing the quadratic galaxy bias and cosmological parameters, the constraints can be tightened to $f_{\rm NL}^{\rm equil}=260\pm 300$, $f_{\rm NL}^{\rm ortho}= -23\pm 120$ (68\% CL).

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Cosmological constraints without fingers of God

Non-linear redshift-space distortions ("fingers of God") are challenging to model analytically, a fact that limits the applicability of perturbation theory in redshift space as compared to real space. We show how this problem can be mitigated using a new observable, $Q_0$, which can be easily estimated from the redshift space clustering data and is approximately equal to the real space power spectrum. The new statistic does not suffer from fingers of God and can be accurately described with perturbation theory down to $k_{\rm max}\simeq 0.4~h~\text{Mpc}^{-1}$. It can be straightforwardly included in the likelihood at negligible additional computational cost, and yields noticeable improvements on cosmological parameters compared to standard power spectrum multipole analyses. Using both simulations and observational data from the Baryon Oscillation Spectroscopic Survey, we show that improvements vary from $10\%$ to $100\%$ depending on the cosmological parameter considered, the galaxy sample and the survey volume.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Cosmology Intertwined: A Review of the Particle Physics, Astrophysics, and Cosmology Associated with the Cosmological Tensions and Anomalies

In this paper we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant $H_0$, the $σ_8$--$S_8$ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the $5.0\,σ$ tension between the {\it Planck} CMB estimate of the Hubble constant $H_0$ and the SH0ES collaboration measurements. After showing the $H_0$ evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the {\it Planck} CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density $Ω_m$, and the amplitude or rate of the growth of structure ($σ_8,fσ_8$). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the $H_0$--$S_8$ tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals.[Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Cosmology with the Redshift-Space Galaxy Bispectrum Monopole at One-Loop Order

We study the cosmological information content of the redshift-space galaxy bispectrum monopole at one-loop order in perturbation theory. We incorporate all effects necessary for comparison to data: fourth-order galaxy bias, infrared resummation (accounting for the non-linear evolution of baryon acoustic oscillations), ultraviolet counterterms, non-linear redshift-space distortions, stochastic contributions, projection, and binning effects. The model is implemented using FFTLog, and validated with the PT Challenge suite of $N$-body simulations, whose large volume allows for high-precision tests. Focusing on the mass fluctuation amplitude, $σ_8$, and galaxy bias parameters, we find that including one-loop corrections allow us to significantly extend the range of scales over which the bispectrum can be modeled, and greatly tightens constraints on bias parameters. However, this does not lead to noticeable improvements in the $σ_8$ errorbar due to the necessary marginalization over a large number of nuisance parameters with conservative priors. Analyzing a BOSS-volume likelihood, we find that the addition of the one-loop bispectrum may lead to improvements on primordial non-Gaussianity constraints by $\lesssim 30\%$ and on $σ_8$ by $\approx 10\%$, though we caution that this requires pushing the analysis to short scales where the galaxy bias parameters may not be correctly recovered; this may lead to biases in the recovered parameter values. We conclude that restrictive priors from simulations or higher-order statistics such as the bispectrum multipoles will be needed in order to realize the full information content of the galaxy bispectrum.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The BOSS DR12 Full-Shape Cosmology: $Λ$CDM Constraints from the Large-Scale Galaxy Power Spectrum and Bispectrum Monopole

We present a full $Λ$CDM analysis of the BOSS DR12 dataset, including information from the power spectrum multipoles, the real-space power spectrum, the reconstructed power spectrum and the bispectrum monopole. This is the first analysis to feature a complete treatment of the galaxy bispectrum, including a consistent theoretical model and without large-scale cuts. Unlike previous works, the statistics are measured using window-free estimators: this greatly reduces computational costs by removing the need to window-convolve the theory model. Our pipeline is tested using a suite of high-resolution mocks and shown to be robust and precise, with systematic errors far below the statistical thresholds. Inclusion of the bispectrum yields consistent parameter constraints and shrinks the $σ_8$ posterior by $13\%$ to reach $<5\%$ precision; less conservative analysis choices would reduce the error-bars further. Our constraints are broadly consistent with Planck: in particular, we find $H_0 = 69.6^{+1.1}_{-1.3}\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$, $σ_8 = 0.692^{+0.035}_{-0.041}$ and $n_s=0.870^{+0.067}_{-0.064}$, including a BBN prior on the baryon density. When $n_s$ is set by Planck, we find $H_0 = 68.31^{+0.83}_{-0.86}\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$ and $σ_8 = 0.722^{+0.032}_{-0.036}$. Our $S_8$ posterior, $0.751\pm0.039$, is consistent with weak lensing studies, but lower than Planck. Constraints on the higher-order bias parameters are significantly strengthened from the inclusion of the bispectrum, and we find no evidence for deviation from the dark matter halo bias relations. These results represent the most complete full-shape analysis of BOSS DR12 to-date, and the corresponding spectra will enable a variety of beyond-$Λ$CDM analyses, probing phenomena such as the neutrino mass and primordial non-Gaussianity.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Faster Fourier Transform? Computing Small-Scale Power Spectra and Bispectra for Cosmological Simulations in $\mathcal{O}(N^2)$ Time

We present $\mathcal{O}(N^2)$ estimators for the small-scale power spectrum and bispectrum in cosmological simulations. In combination with traditional methods, these allow spectra to be efficiently computed across a vast range of scales, requiring orders of magnitude less computation time than Fast Fourier Transform based approaches alone. These methods are applicable to any tracer; simulation particles, halos or galaxies, and take advantage of the simple geometry of the box and periodicity to remove almost all dependence on large random particle catalogs. By working in configuration-space, both power spectra and bispectra can be computed via a weighted sum of particle pairs up to some radius, which can be reduced at larger $k$, leading to algorithms with decreasing complexity on small scales. These do not suffer from aliasing or shot-noise, allowing spectra to be computed to arbitrarily large wavenumbers. The estimators are rigorously derived and tested against simulations, and their covariances discussed. The accompanying code, HIPSTER, has been publicly released, incorporating these algorithms. Such estimators will be of great use in the analysis of large sets of high-resolution simulations.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Cosmology Without Window Functions: Cubic Estimators for the Galaxy Bispectrum

When analyzing the galaxy bispectrum measured from spectroscopic surveys, it is imperative to account for the effects of non-uniform survey geometry. Conventionally, this is done by convolving the theory model with the the window function; however, the computational expense of this prohibits full exploration of the bispectrum likelihood. In this work, we provide a new class of estimators for the unwindowed bispectrum; a quantity that can be straightforwardly compared to theory. This builds upon the work of Philcox (2021) for the power spectrum, and comprises two parts (both obtained from an Edgeworth expansion): a cubic estimator applied to the data, and a Fisher matrix, which deconvolves the bispectrum components. In the limit of weak non-Gaussianity, the estimator is minimum-variance; furthermore, we give an alternate form based on FKP weights that is close-to-optimal and easy to compute. As a demonstration, we measure the binned bispectrum monopole of a suite of simulations both using conventional estimators and our unwindowed equivalents. Computation times are comparable, except that the unwindowed approach requires a Fisher matrix, computable in an additional $\mathcal{O}(100)$ CPU-hours. Our estimator may be straightforwardly extended to measure redshift-space distortions and the components of the bispectrum in arbitrary separable bases. The techniques of this work will allow the bispectrum to straightforwardly be included in the cosmological analysis of current and upcoming survey data.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Determining the Hubble Constant without the Sound Horizon: Measurements from Galaxy Surveys

Two sources of geometric information are encoded in the galaxy power spectrum: the sound horizon at recombination and the horizon at matter-radiation equality. Analyzing the BOSS DR12 galaxy power spectra using perturbation theory with $Ω_m$ priors from Pantheon supernovae but no priors on $Ω_b$, we obtain constraints on $H_0$ from the second scale, finding $H_0 = 65.1^{+3.0}_{-5.4}\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$; this differs from the best-fit of SH0ES at 95\% confidence. Similar results are obtained if $Ω_m$ is constrained from uncalibrated BAO: $H_0 = 65.6^{+3.4}_{-5.5}\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$. Adding the analogous lensing results from Baxter \& Sherwin 2020, the posterior shifts to $70.6^{+3.7}_{-5.0}\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$. Using mock data, Fisher analyses, and scale-cuts, we demonstrate that our constraints do not receive significant information from the sound horizon scale. Since many models resolve the $H_0$ controversy by adding new physics to alter the sound horizon, our measurements are a consistency test for standard cosmology before recombination. A simple forecast indicates that such constraints could reach $σ_{H_0} \simeq 1.6\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$ in the era of Euclid.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Fewer Mocks and Less Noise: Reducing the Dimensionality of Cosmological Observables with Subspace Projections

Creating accurate and low-noise covariance matrices represents a formidable challenge in modern-day cosmology. We present a formalism to compress arbitrary observables into a small number of bins by projection into a model-specific subspace that minimizes the prior-averaged log-likelihood error. The lower dimensionality leads to a dramatic reduction in covariance matrix noise, significantly reducing the number of mocks that need to be computed. Given a theory model, a set of priors, and a simple model of the covariance, our method works by using singular value decompositions to construct a basis for the observable that is close to Euclidean; by restricting to the first few basis vectors, we can capture almost all the constraining power in a lower-dimensional subspace. Unlike conventional approaches, the method can be tailored for specific analyses and captures non-linearities that are not present in the Fisher matrix, ensuring that the full likelihood can be reproduced. The procedure is validated with full-shape analyses of power spectra from BOSS DR12 mock catalogs, showing that the 96-bin power spectra can be replaced by 12 subspace coefficients without biasing the output cosmology; this allows for accurate parameter inference using only $\sim 100$ mocks. Such decompositions facilitate accurate testing of power spectrum covariances; for the largest BOSS data chunk, we find that: (a) analytic covariances provide accurate models (with or without trispectrum terms); and (b) using the sample covariance from the MultiDark-Patchy mocks incurs a $\sim 0.5σ$ shift in $Ω_m$, unless the subspace projection is applied. The method is easily extended to higher order statistics; the $\sim 2000$-bin bispectrum can be compressed into only $\sim 10$ coefficients, allowing for accurate analyses using few mocks and without having to increase the bin sizes.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Precision analysis of the redshift-space galaxy bispectrum

We study the information content of the angle-averaged (monopole) redshift space galaxy bispectrum. The main novelty of our approach is the use of a systematic tree-level perturbation theory model that includes galaxy bias, IR resummation, and also accounts for nonlinear redshift space distortions, binning, and projection effects. We analyze data from the PT challenge simulations, whose cumulative volume of 566 $h^{-3}$Gpc$^3$ allows for a precise comparison to theoretical predictions. Fitting the power spectrum and bispectrum of our simulated data, and varying all necessary cosmological and nuisance parameters in a consistent Markov chain Monte Carlo analysis, we find that our tree-level bispectrum model is valid up to $k_{\max}=0.08~h{\rm Mpc}^{-1}$ (at $z=0.61$). We also find that inclusion of the bispectrum monopole improves constraints on cosmological parameters by $(5-15)\%$ relative to the power spectrum. The improvement is more significant for the quadratic bias parameters of our simulated galaxies, which we also show to deviate from biases of the host dark matter halos at the $\sim 3σ$ level. Finally, we adjust the covariance and scale cuts to match the volume of the BOSS survey, and estimate that within the minimal $Λ$CDM model the bispectrum data can tighten the constraint on the mass fluctuation amplitude $σ_8$ by roughly $10\%$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Combining Full-Shape and BAO Analyses of Galaxy Power Spectra: A 1.6% CMB-independent constraint on H0

We present cosmological constraints from a joint analysis of the pre- and post-reconstruction galaxy power spectrum multipoles from the final data release of the Baryon Oscillation Spectroscopic Survey (BOSS). Geometric constraints are obtained from the positions of BAO peaks in reconstructed spectra, analyzed in combination with the unreconstructed spectra in a full-shape (FS) likelihood using a joint covariance matrix, giving stronger parameter constraints than FS-only or BAO-only analyses. We introduce a new method for obtaining constraints from reconstructed spectra based on a correlated theoretical error, which is shown to be simple, robust, and applicable to any flavor of density-field reconstruction. Assuming $Λ$CDM with massive neutrinos, we analyze data from two redshift bins $z_\mathrm{eff}=0.38,0.61$ and obtain $1.6\%$ constraints on the Hubble constant $H_0$, using only a single prior on the current baryon density $ω_b$ from Big Bang Nucleosynthesis (BBN) and no knowledge of the power spectrum slope $n_s$. This gives $H_0 = 68.6\pm1.1\,\mathrm{km\,s}^{-1}\mathrm{Mpc}^{-1}$, with the inclusion of BAO data sharpening the measurement by $40\%$, representing one of the strongest current constraints on $H_0$ independent of cosmic microwave background data. Restricting to the best-fit slope $n_s$ from Planck (but without additional priors on the spectral shape), we obtain a $1\%$ $H_0$ measurement of $67.8\pm 0.7\,\mathrm{km\,s}^{-1}\mathrm{Mpc}^{-1}$. We find strong constraints on the cosmological parameters from a joint analysis of the FS, BAO, and Planck data. This sets new bounds on the sum of neutrino masses $\sum m_ν< 0.14\,\mathrm{eV}$ (at $95\%$ confidence) and the effective number of relativistic degrees of freedom $N_\mathrm{eff} = 2.90^{+0.15}_{-0.16}$, though contours are not appreciably narrowed by the inclusion of BAO data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Effective Halo Model: Creating a Physical and Accurate Model of the Matter Power Spectrum and Cluster Counts

We introduce a physically-motivated model of the matter power spectrum, based on the halo model and perturbation theory. This model achieves 1\% accuracy on all $k-$scales between $k=0.02h\,\mathrm{Mpc}^{-1}$ to $k=1h\,\mathrm{Mpc}^{-1}$. Our key ansatz is that the number density of halos depends on the non-linear density contrast filtered on some unknown scale $R$. Using the Effective Field Theory of Large Scale Structure to evaluate the two-halo term, we obtain a model for the power spectrum with only two fitting parameters: $R$ and the effective `sound speed&#39;, which encapsulates small-scale physics. This is tested with two suites of cosmological simulations across a broad range of cosmologies and found to be highly accurate. Due to its physical motivation, the statistics can be easily extended beyond the power spectrum; we additionally derive the one-loop covariance matrices of cluster counts and their combination with the matter power spectrum. This yields a significantly better fit to simulations than previous models, and includes a new model for super-sample effects, which is rigorously tested with separate universe simulations. At low redshift, we find a significant ($\sim 10\%$) exclusion covariance from accounting for the finite size of halos which has not previously been modeled. Such power spectrum and covariance models will enable joint analysis of upcoming large-scale structure surveys, gravitational lensing surveys and cosmic microwave background maps on scales down to the non-linear scale. We provide a publicly released Python code.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

What does the Marked Power Spectrum Measure? Insights from Perturbation Theory

The marked power spectrum is capable of placing far tighter constraints on cosmological parameters (particularly the neutrino mass) than the conventional power spectrum. What new information does it contain beyond conventional statistics? Through the development of a perturbative model, we find that the mark induces a significant coupling between non-Gaussianities, which are usually found on small scales, and large scales, leading to the additional information content. The model is derived in the context of one-loop perturbation theory and validated by comparison to $N$-body simulations across a variety of mark parameters. At moderate redshifts, including for massive neutrino cosmologies, the theory is in good agreement with the simulations. The importance of non-linear gravitational effects on the large-scale spectra complicates the modeling as there is no well-defined convergence radius of the theory at low $z$. Extension to higher perturbative order and biased tracers is possible via a similar approach, and a simple model of the latter is shown to yield promising results. The theory becomes non-perturbative at redshift zero for small smoothing scales, with important contributions from higher-order terms: these will need to be studied before the full power of this tool can be realized.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

RascalC: A Jackknife Approach to Estimating Single and Multi-Tracer Galaxy Covariance Matrices

To make use of clustering statistics from large cosmological surveys, accurate and precise covariance matrices are needed. We present a new code to estimate large scale galaxy two-point correlation function (2PCF) covariances in arbitrary survey geometries that, due to new sampling techniques, runs $\sim 10^4$ times faster than previous codes, computing finely-binned covariance matrices with negligible noise in less than 100 CPU-hours. As in previous works, non-Gaussianity is approximated via a small rescaling of shot-noise in the theoretical model, calibrated by comparing jackknife survey covariances to an associated jackknife model. The flexible code, RascalC, has been publicly released, and automatically takes care of all necessary pre- and post-processing, requiring only a single input dataset (without a prior 2PCF model). Deviations between large scale model covariances from a mock survey and those from a large suite of mocks are found to be be indistinguishable from noise. In addition, the choice of input mock are shown to be irrelevant for desired noise levels below $\sim 10^5$ mocks. Coupled with its generalization to multi-tracer data-sets, this shows the algorithm to be an excellent tool for analysis, reducing the need for large numbers of mock simulations to be computed.

0 citations0 reviewsNo code linkedOpen access