Primordial non-Gaussianity Research Articles

Quijote-PNG: Optimizing the Summary Statistics to Measure Primordial Non-Gaussianity

We apply a suite of different estimators to the Quijote-png halo catalogs to find the best approach to constrain Primordial non-Gaussianity (PNG) at nonlinear cosmological scales, up to kmax=0.5hMpc−1 . The set of summary statistics considered in our analysis includes the power spectrum, bispectrum, halo mass function, marked power spectrum, and marked modal bispectrum. Marked statistics are used here for the first time in the context of the PNG study. We perform a Fisher analysis to estimate their cosmological information content, showing substantial improvements when marked observables are added to the analysis. Starting from these summaries, we train deep neural networks to perform likelihood-free inference of cosmological and PNG parameters. We assess the performance of different subsets of summary statistics; in the case of fNLequil , we find that a combination of the power spectrum and a suitable marked power spectrum outperforms the combination of power spectrum and bispectrum, the baseline statistics usually employed in PNG analysis. A minimal pipeline to analyze the statistics we identified can be implemented either with our ML algorithm or via more traditional estimators, if these are deemed more reliable.

Generalising axion-like particle as the curvaton: sourcing primordial density perturbation and non-Gaussianities

We investigate the non-perturbatively generated axion-like particle (ALP) potential, involving fermions in the dark sector that couple to the ALP, in an early universe cosmological inflationary stage with the ALP being a spectator field. The potential here deviates from the standard cosine nature due to the presence of the two dark sector fermion masses mu\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_u$$\\end{document} and md\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_d$$\\end{document} which couple to the ALP. The ALP is a spectator field during inflation but it starts to oscillate and dominates the energy density of the universe after inflation ends, thereby sourcing isocurvature perturbations, while standard curvature fluctuations from the inflaton are assumed to be sub-dominant. Subsequently the ALP decays converting the isocurvature perturbations to adiabatic perturbations thereby acting as the origin of the primordial density perturbations. We identify the parameter space involving the axion decay constant fϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f_{\\phi }$$\\end{document}, scale of confinement Λ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Lambda $$\\end{document}, ALP mass m and the masses of the fermions, mu\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_u$$\\end{document} and md\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_d$$\\end{document} where it can satisfactorily behave as the curvaton and source the observed primordial density perturbation. We also predict local non-Gaussianity signals for bi-spectrum and tri-spectrum fNL\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f_{NL}$$\\end{document} and gNL\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$g_{NL}$$\\end{document}, as a function of the ratio mu/md\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_u/m_d$$\\end{document}, which are within the allowed range in the latest Planck observations and are detectable with future observations. Particularly we find that the value of fNL\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f_{NL}$$\\end{document} and gNL\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$g_{NL}$$\\end{document} are dependent on the ratio of mu\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_u$$\\end{document} and md\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_d$$\\end{document}: fNL\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f_{NL}$$\\end{document} is more or less positive for all scenarios except mu=md\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_u = m_d$$\\end{document} and gNL\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$g_{NL}$$\\end{document} is always positive irrespective of the ratio between mu\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_u$$\\end{document} and md\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_d$$\\end{document}. The results of our analysis in the limit mu=md\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_u = m_d$$\\end{document} resembles vanilla curvaton scenario while in the limit mu≫md\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_u \\gg m_d$$\\end{document} resembles pure axion cosine potential. This is because the model has the unique feature where the extra matter in the dark sector determines the nature and shape of the potential.

Encyclopædia Inflationaris: Opiparous Edition

The current flow of high-accuracy astrophysical data, among which are the Cosmic Microwave Background (CMB) measurements by the Planck satellite, offers an unprecedented opportunity to constrain the inflationary theory. This is however a challenging project given the size of the inflationary landscape which contains hundreds of different scenarios. Given that there is currently no observational evidence for primordial non-Gaussianities, isocurvature perturbations or any other non-minimal extension of the inflationary paradigm, a reasonable approach is to consider the simplest models first, namely the slow-roll single-field models with minimal kinetic terms. This still leaves us with a very populated landscape, the exploration of which requires new and efficient strategies. It has been customary to tackle this problem by means of approximate model-independent methods while a more ambitious alternative is to study the inflationary scenarios one by one. We have developed the publicly available runtime library ASPIC11https://github.com/cosmicinflation/aspic. to implement this last approach. The ASPIC code provides all routines needed to quickly derive reheating-consistent observable predictions within this class of scenarios. ASPIC has been designed as an evolutive code which presently supports 118 different models. In this paper, for each of the ASPIC models, we present and collect new results in a systematic manner, thereby constituting the first Encyclopædia Inflationaris.

New probe of non-Gaussianities with primordial black hole induced gravitational waves

We propose a new probe of primordial non-Gaussianities (NGs) through the observation of gravitational waves (GWs) induced by ultra-light (MPBH<109g) primordial black holes (PBHs). Interestingly enough, the existence of primordial NG can leave imprints on the clustering properties of PBHs and the spectral shape of induced GW signals. Focusing on a scale-dependent local-type NG, we identify a distinct double-peaked GW energy spectrum that, contingent upon MPBH and the abundance of PBHs at the time of formation, denoted as ΩPBH,f, may fall into the frequency bands of upcoming GW observatories, including LISA, ET, SKA, and BBO. Thus, such a signal can serve as a novel portal for probing primordial NGs. Intriguingly, combining BBN bounds on the GW amplitude, we find for the first time the joint limit on the product of the effective non-linearity parameter for the primordial tri-spectrum, denoted by τ¯NL, and the primordial curvature perturbation power spectrum PR(k), which reads as τ¯NLPR(k)<4×10−20ΩPBH,f−17/9(MPBH104g)−17/9.

PNG-UNITsims: Halo clustering response to primordial non-Gaussianities as a function of mass

This paper presents the PNG-UNITSIMS suite, which includes the largest full N-body simulation to date with local primordial non-Gaussianities (local PNG), the PNG-UNIT. The amplitude of the PNGs is given by f localNL=100. The simulation follows the evolution of 40963 particles in a periodic box with Lbox = 1 h−1 Gpc, resulting in a mass resolution of mp = 1.24 × 109 h−1 M⊙, enough to finely resolve the galaxies targeted by stage-IV spectroscopic surveys. The PNG-UNIT has fixed initial conditions with phases also matching the pre-existing UNIT simulation with Gaussian initial conditions. The fixed and matched initial conditions reduce the simulation uncertainty significantly. In this first study of the PNG-UNITSIMS, we measure the PNG response parameter, p, as a function of the halo mass. halos with masses between 1 × 1012 and 5 × 1013 h−1 M⊙ are well described by the universality relation, given by p = 1. For halos with masses between 2 × 1010 and 1 × 1012 h−1 M⊙ we find that p &lt; 1, at a significance between 1.5 and 3.1σ. Combining all the halos between 2 × 1010 and 5 × 1013 h−1 M⊙, we find p consistent with a value of 0.955 ± 0.013, which is 3σ away from the universality relation. We demonstrate that these findings are robust to mass resolution, scale cuts and uncertainty estimation. We also compare our measurements to separate universe simulations, finding that the PNG-UNITSIMS constraints outperform the former for the setup considered. Using a prior on p as tight as the one reported here for DESI-like forecast can result in fNL constraints comparable to fixing p. At the same time, fixing p to a wrong value (p = 1) may result in up to 2σ biases on fNL.

Cosmology with persistent homology: a Fisher forecast

Persistent homology naturally addresses the multi-scale topological characteristics of the large-scale structure as a distribution of clusters, loops, and voids. We apply this tool to the dark matter halo catalogs from the Quijote simulations, and build a summary statistic for comparison with the joint power spectrum and bispectrum statistic regarding their information content on cosmological parameters and primordial non-Gaussianity. Through a Fisher analysis, we find that constraints from persistent homology are tighter for 8 out of the 10 parameters by margins of 13–50%. The complementarity of the two statistics breaks parameter degeneracies, allowing for a further gain in constraining power when combined. We run a series of consistency checks to consolidate our results, and conclude that our findings motivate incorporating persistent homology into inference pipelines for cosmological survey data.

Primordial magnetic non-Gaussianity with generic vacua and detection prospects in CMB spectral distortions

Primordial magnetic non-Gaussianity with generic vacua and detection prospects in CMB spectral distortions

Optimal constraints on Primordial non-Gaussianity with the eBOSS DR16 quasars in Fourier space

We present constraints on the amplitude of local Primordial Non-Gaussianities (PNG), f NL, using the quasar sample in the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 (DR16). We analyze the power spectrum monopole, testing for the presence of scale-dependent galaxy bias induced by local PNG. Our analysis makes use of optimal redshift weights that maximize the response of the quasar sample to the possible presence of non-zero PNG. We find -4 < f NL < 27 at 68% confidence level, which is among the strongest bounds with Large Scale Structure data. The optimal analysis reduces the error bar by ~10% compared to the standard one, but this improvement is lower than the one expected from previous forecasts. In addition, the larger volume of this dataset, when compared to previous releases of the eBOSS quasar catalog, does not always correspond to a reduction of the final uncertainty on local PNG. This could suggest the presence of still unknown systematic effects in the data.If the quasars have a lower response to local PNG, our optimal constraint becomes -23 < f NL < 21 at 68%, with an improvement of 30% over standardanalyses. We alsoshow how to use the optimal weights to put data-driven priors on the sample's response tolocal PNG.

The renormalization group for large-scale structure: primordial non-Gaussianities

The renormalization group for large-scale structure (RG-LSS) describes the evolution of galaxy bias and stochastic parameters as a function of the cutoff Λ. In this work, we introduce interaction vertices that describe primordial non-Gaussianity into the Wilson-Polchinski framework, thereby extending the free theory to the interacting case.The presence of these interactions forces us to include new operators and bias coefficients to the bias expansion to ensure closure under renormalization.We recover the previously-derived “scale-dependent bias” contributions, as well as a new (subdominant) stochastic contribution.We derive the renormalization group equations governing the RG-LSS for a large class of interactions which account for vertices at linear order in f NL that parametrize interacting scalar and massive spinning fields during inflation.Solving the RG equations, we show the evolution of the non-Gaussian contributions to galaxy clustering as a function of scale.

Constraining primordial non-Gaussianity from large scale structure with the wavelet scattering transform

We investigate the Wavelet Scattering Transform (WST) as a tool for the study of Primordial non-Gaussianity (PNG) in Large Scale Structure (LSS), and compare its performance with that achievable viaa joint analysis with power spectrum and bispectrum (P+B). We consider the three main primordial bispectrum shapes — local, equilateral and orthogonal — and produce Fisher forecast for the corresponding f NL amplitudeparameters, jointly with standard cosmological parameters. We analyze simulations from the publicly available Quijote and Quijote-png N-body suites, studying both the dark matter and halo fields.We find that the WST outperforms the power spectrum alone on all parameters, both on the f NL's and on cosmological ones. In particular, on f NL local for halos, the improvement is about 27%.When B is combined with P, halo constraints from WST are weaker for f NL local (at ∼ 15% level), but stronger for f NL equil (∼ 25%) and f NL ortho (∼ 28%).Our results show that WST, both alone and in combination with P+B, can improve the extraction of information on PNG from LSS data over the one attainable by a standard P+B analysis. Moreover, we identify a class of WSTin which the origin of the extra information on PNG can be cleanly isolated.

Parity violation in primordial tensor non-Gaussianities from matter bounce cosmology

It has been shown that primordial tensor non-Gaussianities from a cubic Weyl action with anon-dynamical coupling are suppressed by the so-called slow-roll parameter in a conventional framework of slow-roll inflation. In this paper, we consider matter bounce cosmology in which the background spacetime is no longer quasi-de Sitter, and hence one might expect that the matter bounce models could predict non-suppressed non-Gaussianities. Nevertheless, we first show that the corresponding non-Gaussian amplitudes from the cubic Weyl term with a non-dynamical coupling are much smaller than those from the conventional slow-roll inflation, in spite of the fact that there is no slow-roll suppression. We then introduce a dynamical coupling that can boost the magnitude of graviton cubic interactions and clarify that there is a parameter region where the tensor non-Gaussianities can be enhanced and can potentially be tested by cosmic microwave background experiments.

Complete analysis of the background and anisotropies of scalar-induced gravitational waves: primordial non-Gaussianity f NL and g NL considered

Investigation of primordial non-Gaussianity holds immense importance in testing the inflation paradigm and shedding light on the physics of the early Universe. In this study, we conduct the complete analysis of scalar-induced gravitational waves (SIGWs) by incorporating the local-type non-Gaussianity f NL and g NL. We develop Feynman-like diagrammatic technique and derive semi-analytic formulas for both the energy-density fraction spectrum and the angular power spectrum. For the energy-density fraction spectrum, we analyze all the relevant Feynman-like diagrams, determining their contributions to the spectrum in an order-by-order fashion. As for the angular power spectrum, our focus lies on the initial inhomogeneities, giving rise to anisotropies in SIGWs, that arise from the coupling between short- and long-wavelength modes due to primordial non-Gaussianity. Our analysis reveals that this spectrum exhibits a typical multipole dependence, characterized by C̃ ℓ ∝ [ℓ(ℓ + 1)]-1, which plays a crucial role in distinguishing between different sources of gravitational waves. Depending on model parameters, significant anisotropies can be achieved. We also show that the degeneracies in model parameters can be broken. The findings of our study underscore the angular power spectrum as a robust probe for investigating primordial non-Gaussianity and the physics of the early Universe. Moreover, our theoretical predictions can be tested using space-borne gravitational-wave detectors and pulsar timing arrays.

Modest set of simulations of local-type primordial non-Gaussianity in the DESI era

Modest set of simulations of local-type primordial non-Gaussianity in the DESI era

Multi-tracing the primordial Universe with future surveys

The fluctuations generated by Inflation are nearly Gaussian in the simplest models, but may be non-Gaussian in more complex models, potentially leading to signatures in the late Universe. In particular, local-type primordial non-Gaussianity induces scale-dependent bias in tracers of the matter distribution. This non-Gaussian imprint in the tracer power spectrum survives at late times on ultra-large scales where nonlinearity is negligible. In order to combat the problem of growing cosmic variance on these scales, we use a multi-tracer analysis that combines different tracers to maximise any primordial signal. Previous work has investigated the combination of a spectroscopic galaxy survey with a 21 cm intensity mapping survey in single-dish mode. We extend this work by considering instead the case where the 21 cm intensity mapping survey is optimised for interferometer mode. As examples, we use two multi-tracer pairs of surveys: one at high redshift (1≤z≤2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$1\\le z\\le 2$$\\end{document}) and one at very high redshift (2≤z≤5\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$2\\le z\\le 5$$\\end{document}). The 21 cm surveys are idealised surveys based on HIRAX and PUMA. We implement foreground-avoidance filters and use detailed models of the interferometer thermal noise. The galaxy surveys are idealised surveys based on Euclid and MegaMapper. Via a simple Fisher forecast, we illustrate the potential of the multi-tracer. Our results show a ∼\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sim $$\\end{document} 20–30% improvement in precision on local primordial non-Gaussianity from the multi-tracer. Furthermore, we investigate the effects on constraints of varying the parameter of non-Gaussian galaxy assembly bias and of varying the parameters of the intensity mapping foreground filters. We find that the non-Gaussian galaxy assembly bias parameter causes a greater change in the constraints on local primordial non-Gaussianity than the foreground filter parameters.

Angular bispectrum and trispectrum of scalar-induced gravitational waves: all contributions from primordial non-Gaussianity f NL and g NL

Studying the primordial non-Gaussianity of inflationary perturbations is crucial for testing the inflation paradigm of the early universe. In this work, we conduct a comprehensive analysis of the angular bispectrum and trispectrum of scalar-induced gravitational waves (SIGWs) in the presence of local-type primordial non-Gaussianity parameterized by f NL and g NL, deriving their semi-analytical formulae for the first time. Our findings indicate that it is the presence of primordial non-Gaussianity that leads to a non-Gaussian SIGW background, suggesting that the angular bispectrum and trispectrum of SIGWs could serve as probes of the primordial non-Gaussianity. Our numerical results further illustrate that f NL and g NL exert significant impacts on the spectral amplitudes, potentially reaching up to 10-5 for the former and 10-8 for the latter. In particular, we demonstrate that the angular bispectrum and trispectrum exhibit characteristic dependence on the angular multipoles and frequency bands. They hold potentials to be measured by gravitational-wave detectors that may advance our understanding of the origin of the universe.

From inflation to dark matter halo profiles: the impact of primordial non-Gaussianities on the central density cusp

It has recently been shown that local primordial non-Gaussianities (PNG) with significant amplitude (|f NL| ∼ 1000), at small (Mpc) scales, can help in forming simulated galaxies with more disky baryonic kinematics than in the Gaussian case, while generating matter power spectra that can differ by up to 20% from the Gaussian case at non-linear scales. Here, we explore in detail the consequences of such small-scale PNG on the dark matter halo profiles. We show in particular that, for negative f NL, dark matter halos formed in collisionless simulations are not always well described by the traditional Navarro-Frenk-White (NFW) profiles, as supported by their sparsity distribution. We conclude that NFW profiles are not as clear attractors for the density profiles of dark matter halos in the presence of PNG than in the case of a Gaussian contrast density field. We show how alternatives to the NFW profile can describe halos both in the Gaussian and non-Gaussian cases. From the combination of our sparsity analysis and the quality of the adjustments of the density profiles with a minimal extension to NFW, we conclude that z = 1 halos carry the most interesting information about PNG.

Analysis of BOSS galaxy data with weighted skew-spectra

We present the first application of the weighted skew-spectra to analyze non-Gaussian information in galaxy survey data. Using the tree-level galaxy skew-spectra together with the one-loop power spectrum multipoles, we analyze the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) galaxy clustering data, and target our search towards the equilateral bispectrum shape of primordial non-Gaussianity. We use the Effective Field Theory model for the galaxy power spectrum and bispectrum, and account for systematic effects, such as the survey geometry. From our likelihood analysis, we find f NL equil = -34+296 -334 at 68% CL, consistent with previous works, while systematic errors from our treatment of the survey geometry lead to an unreliable estimation of f NL ortho. We further constrain the bias and counterterm parameters, while keeping the cosmology fixed to Planck 2018 values. As a check, we also validate our analysis pipeline using the Nseries simulation suite.

BCS in the sky: signatures of inflationary fermion condensation

We consider a Bardeen-Cooper-Schrieffer (BCS)-like model in the inflationary background. We show that with an axial chemical potential, the attractive quartic fermion self-interaction can lead to a BCS-like condensation. In the rigid-de Sitter (dS) limit of inflation where backreaction from the inflaton and graviton is neglected, we perform the first computation of the non-perturbative effective potential that includes the full spacetime curvature effects in the presence of the chemical potential, subject to the mean-field approximation whose validity has been checked via the Ginzburg criterion. The corresponding BCS phase transition is always first-order, when the varying Hubble is interpreted as an effective Gibbons-Hawking temperature of dS spacetime. In the condensed phase, the theory can be understood from UV and IR sides as fermionic and bosonic, respectively. This leads to distinctive signatures in the primordial non-Gaussianity of curvature perturbations. Namely, the oscillatory cosmological collider signal is smoothly turned off at a finite momentum ratio, since different momentum ratios effectively probe different energy scales. In addition, such BCS phase transitions can also source stochastic gravitational waves, which are feasible for future experiments.

CMB bispectrum constraints on DHOST inflation

We present the first direct constraints on a Degenerate Higher Order Scalar Tensor (DHOST) inflation model using the Planck 2018 Cosmic Microwave Background (CMB) results on non-Gaussianities. We identify that the bispectrum consists of a fixed contribution following from the power spectrum and a linear combination of terms depending on five free parameters defining the cubic perturbations to the DHOST model. The former peaks in the squeezed limit, while the latter is maximised in the equilateral limit. We directly confront the model predictions to the CMB bispectrum statistics via the public code CMB-BEST and marginalize over the free parameters. We explicitly show that there are viable DHOST inflationary models satisfying both power spectrum and bispectrum constraints from Planck. However, rather surprisingly, the constraints exclude certain models at the 6σ-level even though they pass the conventional fudge factor tests. In this case and despite having a handful of free parameters, the model's large squeezed bispectrum cannot be cancelled out without introducing a large bispectrum in other limits which are strongly constrained by Planck's non-detection of primordial non-Gaussianity. We emphasize that first-order approximations such as fudge factors, albeit commonly used in the literature, may be misleading and provide weaker constraints. A proper analysis of the constraints from Planck requires a more robust approach, such as the one provided by the CMB-BEST code.

Primordial non-Gaussianities with weak lensing: information on non-linear scales in the Ulagam full-sky simulations

Primordial non-Gaussianities (PNGs) are signatures in the density field that encode particle physics processes from the inflationary epoch. Such signatures have been extensively studied using the Cosmic Microwave Background, through constraining their amplitudes, fX NL, with future improvements expected from large-scale structure surveys; specifically, the galaxy correlation functions. We show that weak lensing fields can be used to achieve competitive and complementary constraints. This is shown via the Ulagam suite of N-body simulations, a subset of which evolves primordial fields with four types of PNGs. We create full-sky lensing maps and estimate the Fisher information from three summary statistics measured on the maps: the moments, the cumulative distribution function, and the 3-point correlation function. We find that the year 10 sample from the Rubin Observatory Legacy Survey of Space and Time (LSST) can constrain PNGs to σ(f NL eq) ≈ 110, σ(f NL or, lss) ≈ 120, σ(f NL loc) ≈ 40. For the former two, this is better than or comparable to expected galaxy clustering-based constraints from the Dark Energy Spectroscopic Instrument (DESI). The PNG information in lensing fields is on non-linear scales and at low redshifts (z ≲ 1.25), with a clear origin in the evolution history of massive halos. The constraining power degrades by ∼60% under scale cuts of ≳ 20 Mpc, showing there is still significant information on scales mostly insensitive to small-scale systematic effects (e.g., baryons). We publicly release the Ulagam suite to enable more survey-focused analyses.