Search

The last few years have seen the emergence of a wide array of novel techniques for analyzing high-precision data from upcoming galaxy surveys, which aim to extend the statistical analysis of galaxy clustering data beyond the linear regime and the canonical two-point (2pt) statistics. We test and benchmark some of these new techniques in a community data challenge "Beyond-2pt", initiated during the Aspen 2022 Summer Program "Large-Scale Structure Cosmology beyond 2-Point Statistics," whose first round of results we present here. The challenge dataset consists of high-precision mock galaxy catalogs for clustering in real space, redshift space, and on a light cone. Participants in the challenge have developed end-to-end pipelines to analyze mock catalogs and extract unknown ("masked") cosmological parameters of the underlying $\Lambda$CDM models with their methods. The methods represented are density-split clustering, nearest neighbor statistics, BACCO power spectrum emulator, void statistics, LEFTfield field-level inference using effective field theory (EFT), and joint power spectrum and bispectrum analyses using both EFT and simulation-based inference. In this work, we review the results of the challenge, focusing on problems solved, lessons learned, and future research needed to perfect the emerging beyond-2pt approaches. The unbiased parameter recovery demonstrated in this challenge by multiple statistics and the associated modeling and inference frameworks supports the credibility of cosmology constraints from these methods. The challenge data set is publicly available and we welcome future submissions from methods that are not yet represented., Comment: New submissions welcome! Challenge data available at https://github.com/ANSalcedo/Beyond2ptMock

Massive particles leave imprints on primordial non-Gaussianity via couplings to the inflaton, even despite their exponential dilution during inflation: practically, the Universe acts as a Cosmological Collider. We present the first dedicated search for spin-zero particles using BOSS redshift-space galaxy power spectrum and bispectrum multipoles, as well as Planck CMB non-Gaussianity data. We demonstrate that some Cosmological Collider models are well approximated by the standard equilateral and orthogonal parametrization; assuming negligible inflaton self-interactions, this facilitates us translating Planck non-Gaussianity constraints into bounds on Collider models. Many models have signatures that are not degenerate with equilateral and orthogonal non-Gaussianity and thus require dedicated searches. Here, we constrain such models using BOSS three-dimensional redshift-space galaxy clustering data, focusing on spin-zero particles in the principal series and constraining their couplings to the inflaton at varying speed and mass, marginalizing over the unknown inflaton self-interactions. This is made possible through an improvement in Cosmological Bootstrap techniques and the combination of perturbation theory and halo occupation distribution models for galaxy clustering. Our work sets the standard for inflationary spectroscopy with cosmological observations, providing the ultimate link between physics on the largest and smallest scales., Comment: 35 pages, 15 figures, 6 tables

We study how well perturbative forward modeling can constrain cosmological parameters compared to conventional analyses. We exploit the fact that in perturbation theory the field-level posterior can be computed analytically in the limit of small noise. In the idealized case where the only relevant parameter for the nonlinear evolution is the nonlinear scale, we argue that information content in this posterior is the same as in the $n$-point correlation functions computed at the same perturbative order. In the real universe other parameters can be important, and there are possibly enhanced effects due to nonlinear interactions of long and short wavelength fluctuations that can either degrade the signal or increase covariance matrices. We identify several different parameters that control these enhancements and show that for some shapes of the linear power spectrum they can be large. This leads to degradation of constraints in the standard analyses, even though the effects are not dramatic for a $\Lambda$CDM-like cosmology. The aforementioned long-short couplings do not affect the field-level inference which remains optimal. Finally, we show how in these examples calculation of the perturbative posterior motivates new estimators that are easier to implement in practice than the full forward modelling but lead to nearly optimal constraints on cosmological parameters., Comment: 33 pages, 1 figure, matches published version

We present a framework for self-consistent cosmological analyses of the full-shape anisotropic bispectrum, including the quadrupole $(\ell=2)$ and hexadecapole $(\ell=4)$ moments. This features a novel window-free algorithm for extracting the latter quantities from data, derived using a maximum-likelihood prescription. Furthermore, we introduce a theoretical model for the bispectrum multipoles (which does not introduce new free parameters), and test both aspects of the pipeline on several high-fidelity mocks, including the PT Challenge suite of gigantic cumulative volume. This establishes that the systematic error is significantly below the statistical threshold, both for the measurement and modeling. As a realistic example, we extract the large-scale bispectrum multipoles from BOSS DR12 and analyze them in combination with the power spectrum data. Assuming a minimal $\Lambda$CDM model, with a BBN prior on the baryon density and a \textit{Planck} prior on $n_s$, we can extract the remaining cosmological parameters directly from the clustering data. The inclusion of the unwindowed higher-order $(\ell>0)$ large-scale bispectrum multipoles is found to moderately improve one-dimensional cosmological parameter posteriors (at the $5\%-10\%$ level), though these multipoles are detected only in three out of four BOSS data segments at $\approx 5\sigma$. Combining information from the power spectrum and bispectrum multipoles, the real space power spectrum, and the post-reconstructed BAO data, we find $H_0 = 68.2\pm 0.8~\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$, $\Omega_m =0.33\pm 0.01$ and $\sigma_8 = 0.736\pm 0.033$ (the tightest yet found in perturbative full-shape analyses). Our estimate of the growth parameter $S_8=0.77\pm 0.04$ agrees with both weak lensing and CMB results., Comment: 40 pages, 8 figures, 6 tables; estimators and data are publicly available at https://github.com/oliverphilcox/Spectra-Without-Windows

We search for ultra-light axions as dark matter (DM) and dark energy particle candidates, for axion masses $10^{-32}\,\mathrm{eV} \leq m_\mathrm{a} \leq 10^{-24}\,\mathrm{eV}$, by a joint analysis of cosmic microwave background (CMB) and galaxy clustering data -- and consider if axions can resolve the tension in inferred values of the matter clustering parameter $S_8$. We give legacy constraints from Planck 2018 CMB data, improving 2015 limits on the axion density $\Omega_\mathrm{a} h^2$ by up to a factor of three; CMB data from the Atacama Cosmology Telescope and the South Pole Telescope marginally weaken Planck bounds at $m_\mathrm{a} = 10^{-25}\,\mathrm{eV}$, owing to lower (and theoretically-consistent) gravitational lensing signals. We jointly infer, from Planck CMB and full-shape galaxy power spectrum and bispectrum data from the Baryon Oscillation Spectroscopic Survey (BOSS), that axions are, today, $< 10\%$ of the DM for $m_\mathrm{a} \leq 10^{-26}\,\mathrm{eV}$ and $< 1\%$ for $10^{-30}\,\mathrm{eV} \leq m_\mathrm{a} \leq 10^{-28}\,\mathrm{eV}$. BOSS data strengthen limits, in particular at higher $m_\mathrm{a}$ by probing high-wavenumber modes ($k < 0.4 h\,\mathrm{Mpc}^{-1}$). BOSS alone finds a preference for axions at $2.7 \sigma$, for $m_\mathrm{a} = 10^{-26}\,\mathrm{eV}$, but Planck disfavours this result. Nonetheless, axions in a window $10^{-28}\,\mathrm{eV} \leq m_\mathrm{a} \leq 10^{-25}\,\mathrm{eV}$ can improve consistency between CMB and galaxy clustering data, e.g., reducing the $S_8$ discrepancy from $2.7 \sigma$ to $1.6 \sigma$, since these axions suppress structure growth at the $8 h^{-1}\,\mathrm{Mpc}$ scales to which $S_8$ is sensitive. We expect improved constraints with upcoming high-resolution CMB and galaxy lensing and future galaxy clustering data, where we will further assess if axions can restore cosmic concordance., Comment: 53 pages, 22 figures. Minor changes to match version accepted for publication in JCAP

We forecast the constraints on single-field inflation from the bispectrum of future high-redshift surveys such as MegaMapper. Considering non-local primordial non-Gaussianity (NLPNG), we find that current methods will yield constraints of order $\sigma(f_{\rm NL}^{\rm eq})\approx 23$, $\sigma(f_{\rm NL}^{\rm orth})\approx 12$ in a joint power-spectrum and bispectrum analysis, varying both nuisance parameters and cosmology, including a conservative range of scales. Fixing cosmological parameters and quadratic bias parameter relations, the limits tighten significantly to $\sigma(f_{\rm NL}^{\rm eq})\approx 17$, $\sigma(f_{\rm NL}^{\rm orth})\approx 8$. These compare favorably with the forecasted bounds from CMB-S4: $\sigma(f_{\rm NL}^{\rm eq})\approx 21$, $\sigma(f_{\rm NL}^{\rm orth})\approx 9$, with a combined constraint of $\sigma(f_{\rm NL}^{\rm eq})\approx 14$, $\sigma(f_{\rm NL}^{\rm orth})\approx 7$; this weakens only slightly if one instead combines with data from the Simons Observatory. We additionally perform a range of Fisher analyses for the error, forecasting the dependence on nuisance parameter marginalization, scale cuts, and survey strategy. Lack of knowledge of bias and counterterm parameters is found to significantly limit the information content; this could be ameliorated by tight simulation-based priors on the nuisance parameters. The error-bars decrease significantly as the number of observed galaxies and survey depth is increased: as expected, deep dense surveys are the most constraining, though it will be difficult to reach $\sigma(f_{\rm NL})\approx 1$ with current methods. The NLPNG constraints will tighten further with improved theoretical models (incorporating higher-loop corrections and improved understanding of nuisance parameters), as well as the inclusion of additional higher-order statistics., Comment: 6 pages, 3 figures, accepted by Phys. Lett. B

Could new physics break the mirror symmetry of the Universe? Utilizing recent measurements of the parity-odd four-point correlation function of BOSS galaxies, we probe the physics of inflation by placing constraints on the amplitude of a number of parity-violating models. Within canonical models of (single-field, slow-roll) inflation, no parity-asymmetry can occur; however, it has recently been shown that breaking of the standard assumptions can lead to parity violation within the Effective Field Theory of Inflation (EFTI). In particular, we consider the Ghost Condensate and Cosmological Collider scenarios -- the former for the leading and subleading operators in the EFTI and the latter for different values of mass and speed of an exchanged spin-$1$ particle -- for a total of $18$ models. Each instance yields a definite prediction for the inflationary trispectrum, which we convert to a late-time galaxy correlator prediction (through a highly non-trivial calculation) and constrain using the observed data. We find no evidence for inflationary parity-violation (with each of the $18$ models having significances below $2\sigma$), and place the first constraints on the relevant coupling strengths, at a level comparable with the theoretical perturbativity bounds. This is also the first time Cosmological Collider signatures have directly been searched for in observational data. We further show that possible secondary parity-violating signatures in galaxy clustering can be systematically described within the Effective Field Theory of Large-Scale Structure. We argue that these late-time contributions are subdominant compared to the primordial parity-odd signal for a vast region of parameter space. In summary, the results of this paper disfavor the notion that the recent hints of parity-violation observed in the distribution of galaxies are due to new physics., Comment: 20+9 spooky pages, 4 figures, accepted by Phys. Rev. D. Code available at https://github.com/oliverphilcox/Parity-Odd-4PCF

We summarize recent progress in the development, application, and understanding of effective field theories and highlight promising directions for future research. This Report is prepared as the TF02 "Effective Field Theory" topical group summary for the Theory Frontier as part of the Snowmass 2021 process., Comment: 12 pages

We derive a set of no-go theorems and yes-go examples for the parity-odd primordial trispectrum of curvature perturbations. We work at tree-level in the decoupling limit of the Effective Field Theory of Inflation and assume scale invariance and a Bunch-Davies vacuum. We show that the parity-odd scalar trispectrum vanishes in the presence of any number of scalar fields with arbitrary mass and any parity-odd scalar correlator vanishes in the presence of any number of spinning fields with massless de Sitter mode functions, in agreement with the findings of Liu, Tong, Wang and Xianyu [1]. The same is true for correlators with an odd number of conformally-coupled external fields. We derive these results using both the (boostless) cosmological bootstrap, in particular the Cosmological Optical Theorem, and explicit perturbative calculations. We then discuss a series of yes-go examples by relaxing the above assumptions one at the time. In particular, we provide explicit results for the parity-odd trispectrum for (i) violations of scale invariance in single-clock inflation, (ii) the modified dispersion relation of the ghost condensate (non-Bunch-Davies vacuum), and (iii) interactions with massive spinning fields. Our results establish the parity-odd trispectrum as an exceptionally sensitive probe of new physics beyond vanilla inflation., Comment: 38 pages, 5 figures. v2 fixes typo in Eq. (2.8)

We summarize progress made in theoretical astrophysics and cosmology over the past decade and areas of interest for the coming decade. This Report is prepared as the TF09 "Astrophysics and Cosmology" topical group summary for the Theory Frontier as part of the Snowmass 2021 process., Comment: 57 pages

We derive parity-even graviton bispectra in the Effective Field Theory of Inflation (EFToI) to all orders in derivatives. Working in perturbation theory, we construct all cubic interactions that can contribute to tree-level graviton bispectra, showing that they all come from EFToI operators containing two or three powers of the extrinsic curvature and its covariant derivatives: all other operators can be removed by field redefinitions or start at higher-order in perturbations. For operators cubic in the extrinsic curvature, where the single-clock consistency relations are satisfied without a correction to the graviton two-point function, we use the Manifestly Local Test (MLT) to efficiently extract the effects of evolving graviton fluctuations to the end of inflation. Despite the somewhat complicated nature of the bulk interactions, the final boundary correlators take a very compact form. For operators quadratic in the extrinsic curvature, the leading order bispectra are a sum of contact and single exchange diagrams, which are tied together by spatial diffeomorphisms, and to all orders in derivatives we derive these bispectra by computing the necessary bulk time integrals. For single exchange diagrams we exploit factorisation properties of the bulk-bulk propagator for massless gravitons and write the result as a finite sum over residues. Perhaps surprisingly, we show these single exchange contributions have only total-energy poles and also satisfy the MLT.

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, $\sigma_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 $\sigma_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 $\sigma_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., Comment: 51 pages, 6 figures, one-loop redshift-space galaxy bispectrum code is publicly available at https://github.com/oliverphilcox/OneLoopBispectrum

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., Comment: 20 pages, 5 figures

Small fluctuations around homogeneous and isotropic expanding backgrounds are the main object of study in cosmology. Their origin and evolution is sensitive to the physical processes that happen during inflation and in the late Universe. As such, they hold the key to answering many of the major open questions in cosmology. Given a large separation of relevant scales in many examples of interest, the most natural description of these fluctuations is formulated in terms of effective field theories. This was the main avenue for many of the important modern developments in theoretical cosmology, which provided a unifying framework for a plethora of cosmological models and made a clear connection between the fundamental cosmological parameters and observables. In this review we summarize these results in the context of effective field theories of inflation, large-scale structure, and dark energy., Comment: Contribution to Snowmass 2021

Cosmic inflation provides a window to the highest energy densities accessible in nature, far beyond those achievable in any realistic terrestrial experiment. Theoretical insights into the inflationary era and its observational probes may therefore shed unique light on the physical laws underlying our universe. This white paper describes our current theoretical understanding of the inflationary era, with a focus on the statistical properties of primordial fluctuations. In particular, we survey observational targets for three important signatures of inflation: primordial gravitational waves, primordial non-Gaussianity and primordial features. With the requisite advancements in analysis techniques, the tremendous increase in the raw sensitivities of upcoming and planned surveys will translate to leaps in our understanding of the inflationary paradigm and could open new frontiers for cosmology and particle physics. The combination of future theoretical and observational developments therefore offer the potential for a dramatic discovery about the nature of cosmic acceleration in the very early universe and physics on the smallest scales., Comment: Contribution to Snowmass 2021; 103 pages, 8 figures, 378 endorsers; v2: minor changes

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)., Comment: 6+7 pages, 2+3 figures, 0+1 table, code available at github.com/Michalychforever/CLASS-PT, with custom MontePython likelihoods available at github.com/oliverphilcox/full_shape_likelihoods