Showing 1–18 of 18 results for author: Jung, G

Cosmological simulations of scale-dependent primordial non-Gaussianity

Authors: Marco Baldi, Emanuele Fondi, Dionysios Karagiannis, Lauro Moscardini, Andrea Ravenni, William R. Coulton, Gabriel Jung, Michele Liguori, Marco Marinucci, Licia Verde, Francisco Villaescusa-Navarro, Banjamin D. Wandelt

Abstract: We present the results of a set of cosmological N-body simulations with standard $Λ$CDM cosmology but characterized by a scale-dependent primordial non-Gaussianity of the local type featuring a power-law dependence of the $f_{\rm NL}^{\rm loc}(k)$ at large scales followed by a saturation to a constant value at smaller scales where non-linear growth leads to the formation of collapsed cosmic struct… ▽ More We present the results of a set of cosmological N-body simulations with standard $Λ$CDM cosmology but characterized by a scale-dependent primordial non-Gaussianity of the local type featuring a power-law dependence of the $f_{\rm NL}^{\rm loc}(k)$ at large scales followed by a saturation to a constant value at smaller scales where non-linear growth leads to the formation of collapsed cosmic structures. Such models are built to ensure consistency with current Cosmic Microwave Background bounds on primordial non-Gaussianity yet allowing for large effects of the non-Gaussian statistics on the properties of non-linear structure formation. We show the impact of such scale-dependent non-Gaussian scenarios on a wide range of properties of the resulting cosmic structures, such as the non-linear matter power spectrum, the halo and sub-halo mass functions, the concentration-mass relation, the halo and void density profiles, and we highlight for the first time that some of these models might mimic the effects of Warm Dark Matter for several of such observables △ Less

Submitted 11 July, 2024; v1 submitted 9 July, 2024; originally announced July 2024.

Comments: 21 pages, 9 figures, 2 tables; to be submitted to JCAP

Revisiting the large-scale CMB anomalies: The impact of the SZ signal from the Local Universe

Authors: Gabriel Jung, Nabila Aghanim, Jenny Sorce, Benjamin Seidel, Klaus Dolag, Marian Douspis

Abstract: The full sky measurements of the Cosmic Microwave Background (CMB) temperature anisotropies by $\textit{WMAP}$ and $\textit{Planck}$ have highlighted the presence of several unexpected isotropy-breaking features on the largest angular scales. In this work, we investigate the impact of the local large-scale structure on these anomalies through the thermal and kinetic Sunyaev-Zeldovich effects. We… ▽ More The full sky measurements of the Cosmic Microwave Background (CMB) temperature anisotropies by $\textit{WMAP}$ and $\textit{Planck}$ have highlighted the presence of several unexpected isotropy-breaking features on the largest angular scales. In this work, we investigate the impact of the local large-scale structure on these anomalies through the thermal and kinetic Sunyaev-Zeldovich effects. We use a constrained hydrodynamical simulation that reproduces the local Universe in a box of $500\,h^{-1}\,$Mpc to construct full sky maps of the temperature anisotropies produced by these two CMB secondary effects and discuss their statistical properties on large angular scales. We show the significant role played by the Virgo cluster on these scales, and compare it to theoretical predictions and random patches of the universe obtained from the hydrodynamical simulation $\textit{Magneticum}$. We explore three of the main CMB large-scale anomalies -- i.e., lack of correlation, quadrupole-octopole alignment and hemispherical asymmetry -- , both in the latest $\textit{Planck}$ data (PR4), where they are detected at a similar level to the previous releases, and using the simulated secondaries from the local Universe, verifying their negligible impact. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 15 pages, 20 figures

Constraining Primordial Non-Gaussianity from Large Scale Structure with the Wavelet Scattering Transform

Authors: Matteo Peron, Gabriel Jung, Michele Liguori, Massimo Pietroni

Abstract: 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 via a 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 correspondin… ▽ More 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 via a 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 fNL amplitude parameters, 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 fNL's and on cosmological ones. In particular, on fNL_loc for halos, the improvement is about 27%. When B is combined with P, halo constraints from WST are weaker for fNL_loc (at ~ 15% level), but stronger for fNL_eq (~ 25%) and fNL_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 WST in which the origin of the extra information on PNG can be cleanly isolated. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Comments: 25 pages, 13 figures

Quijote-PNG: Optimizing the summary statistics to measure Primordial non-Gaussianity

Authors: Gabriel Jung, Andrea Ravenni, Michele Liguori, Marco Baldi, William R. Coulton, Francisco Villaescusa-Navarro, Benjamin D. Wandelt

Abstract: We apply a suite of different estimators to the Quijote-PNG halo catalogues to find the best approach to constrain Primordial non-Gaussianity (PNG) at non-linear cosmological scales, up to $k_{\rm max} = 0.5 \, h\,{\rm Mpc}^{-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.… ▽ More We apply a suite of different estimators to the Quijote-PNG halo catalogues to find the best approach to constrain Primordial non-Gaussianity (PNG) at non-linear cosmological scales, up to $k_{\rm max} = 0.5 \, h\,{\rm Mpc}^{-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 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 (NN) to perform likelihood-free inference of cosmological and PNG parameters. We assess the performance of different subsets of summary statistics; in the case of $f_\mathrm{NL}^\mathrm{equil}$, 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 analyse the statistics we identified can be implemented either with our ML algorithm or via more traditional estimators, if these are deemed more reliable. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 13 pages, 10 figures

Taming assembly bias for primordial non-Gaussianity

Authors: Emanuele Fondi, Licia Verde, Francisco Villaescusa-Navarro, Marco Baldi, William R. Coulton, Gabriel Jung, Dionysios Karagiannis, Michele Liguori, Andrea Ravenni, Benjamin D. Wandelt

Abstract: Primordial non-Gaussianity of the local type induces a strong scale-dependent bias on the clustering of halos in the late-time Universe. This signature is particularly promising to provide constraints on the non-Gaussianity parameter $f_{\rm NL}$ from galaxy surveys, as the bias amplitude grows with scale and becomes important on large, linear scales. However, there is a well-known degeneracy betw… ▽ More Primordial non-Gaussianity of the local type induces a strong scale-dependent bias on the clustering of halos in the late-time Universe. This signature is particularly promising to provide constraints on the non-Gaussianity parameter $f_{\rm NL}$ from galaxy surveys, as the bias amplitude grows with scale and becomes important on large, linear scales. However, there is a well-known degeneracy between the real prize, the $f_{\rm NL}$ parameter, and the (non-Gaussian) assembly bias i.e., the halo formation history-dependent contribution to the amplitude of the signal, which could seriously compromise the ability of large-scale structure surveys to constrain $f_{\rm NL}$. We show how the assembly bias can be modeled and constrained, thus almost completely recovering the power of galaxy surveys to competitively constrain primordial non-Gaussianity. In particular, studying hydrodynamical simulations, we find that a proxy for the halo properties that determine assembly bias can be constructed from photometric properties of galaxies. Using a prior on the assembly bias guided by this proxy degrades the statistical errors on $f_{\rm NL}$ only mildly compared to an ideal case where the assembly bias is perfectly known. The systematic error on $f_{\rm NL}$ that the proxy induces can be safely kept under control. △ Less

Submitted 2 February, 2024; v1 submitted 16 November, 2023; originally announced November 2023.

Comments: 30 pages, 13 figures. v2: minor updates to match accepted version

Journal ref: JCAP02(2024)048

Quijote-PNG: The Information Content of the Halo Mass Function

Authors: Gabriel Jung, Andrea Ravenni, Marco Baldi, William R. Coulton, Drew Jamieson, Dionysios Karagiannis, Michele Liguori, Helen Shao, Licia Verde, Francisco Villaescusa-Navarro, Benjamin D. Wandelt

Abstract: We study signatures of primordial non-Gaussianity (PNG) in the redshift-space halo field on non-linear scales, using a combination of three summary statistics, namely the halo mass function (HMF), power spectrum, and bispectrum. The choice of adding the HMF to our previous joint analysis of power spectrum and bispectrum is driven by a preliminary field-level analysis, in which we train graph neura… ▽ More We study signatures of primordial non-Gaussianity (PNG) in the redshift-space halo field on non-linear scales, using a combination of three summary statistics, namely the halo mass function (HMF), power spectrum, and bispectrum. The choice of adding the HMF to our previous joint analysis of power spectrum and bispectrum is driven by a preliminary field-level analysis, in which we train graph neural networks on halo catalogues to infer the PNG $f_\mathrm{NL}$ parameter. The covariance matrix and the responses of our summaries to changes in model parameters are extracted from a suite of halo catalogues constructed from the Quijote-PNG N-body simulations. We consider the three main types of PNG: local, equilateral and orthogonal. Adding the HMF to our previous joint analysis of power spectrum and bispectrum produces two main effects. First, it reduces the equilateral $f_\mathrm{NL}$ predicted errors by roughly a factor $2$, while also producing notable, although smaller, improvements for orthogonal PNG. Second, it helps break the degeneracy between the local PNG amplitude, $f_\mathrm{NL}^\mathrm{local}$, and assembly bias, $b_φ$, without relying on any external prior assumption. Our final forecasts for PNG parameters are $Δf_\mathrm{NL}^\mathrm{local} = 40$, $Δf_\mathrm{NL}^\mathrm{equil} = 210$, $Δf_\mathrm{NL}^\mathrm{ortho} = 91$, on a cubic volume of $1 \left(h^{-1}{\rm Gpc}\right)^3$, with a halo number density of $\bar{n}\sim 5.1 \times 10^{-5}~h^3\mathrm{Mpc}^{-3}$, at $z = 1$, and considering scales up to $k_\mathrm{max} = 0.5~h\,\mathrm{Mpc}^{-1}$. △ Less

Submitted 4 February, 2024; v1 submitted 17 May, 2023; originally announced May 2023.

Comments: 17 pages, 11 figures. v3 (minor caption fix)

Journal ref: Astrophys.J. 957 (2023) 1, 50

Quijote-PNG: Quasi-maximum likelihood estimation of Primordial Non-Gaussianity in the non-linear halo density field

Authors: Gabriel Jung, Dionysios Karagiannis, Michele Liguori, Marco Baldi, William R Coulton, Drew Jamieson, Licia Verde, Francisco Villaescusa-Navarro, Benjamin D. Wandelt

Abstract: We study primordial non-Gaussian signatures in the redshift-space halo field on non-linear scales, using a quasi-maximum likelihood estimator based on optimally compressed power spectrum and modal bispectrum statistics. We train and validate the estimator on a suite of halo catalogues constructed from the Quijote-PNG N-body simulations, which we release to accompany this paper. We verify its unbia… ▽ More We study primordial non-Gaussian signatures in the redshift-space halo field on non-linear scales, using a quasi-maximum likelihood estimator based on optimally compressed power spectrum and modal bispectrum statistics. We train and validate the estimator on a suite of halo catalogues constructed from the Quijote-PNG N-body simulations, which we release to accompany this paper. We verify its unbiasedness and near optimality, for the three main types of primordial non-Gaussianity (PNG): local, equilateral, and orthogonal. We compare the modal bispectrum expansion with a $k$-binning approach, showing that the former allows for faster convergence of numerical derivatives in the computation of the score-function, thus leading to better final constraints. We find, in agreement with previous studies, that the local PNG signal in the halo-field is dominated by the scale-dependent bias signature on large scales and saturates at $k \sim 0.2~h\,\mathrm{Mpc}^{-1}$, whereas the small-scale bispectrum is the main source of information for equilateral and orthogonal PNG. Combining power spectrum and bispectrum on non-linear scales plays an important role in breaking degeneracies between cosmological and PNG parameters; such degeneracies remain however strong for equilateral PNG. We forecast that PNG parameters can be constrained with $Δf_\mathrm{NL}^\mathrm{local} = 45$, $Δf_\mathrm{NL}^\mathrm{equil} = 570$, $Δf_\mathrm{NL}^\mathrm{ortho} = 110$, on a cubic volume of $1 \left({ {\rm Gpc}/{ {\rm h}}} \right)^3$, at $z = 1$, considering scales up to $k_\mathrm{max} = 0.5~h\,\mathrm{Mpc}^{-1}$. △ Less

Submitted 18 May, 2023; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: 17 pages, 10 figures. v2: minor updates to match published version

Journal ref: Astrophys.J. 948 (2023) 2, 135

Quijote PNG: The information content of the halo power spectrum and bispectrum

Authors: William R Coulton, Francisco Villaescusa-Navarro, Drew Jamieson, Marco Baldi, Gabriel Jung, Dionysios Karagiannis, Michele Liguori, Licia Verde, Benjamin D. Wandelt

Abstract: We investigate how much can be learnt about four types of primordial non-Gaussianity (PNG) from small-scale measurements of the halo field. Using the QUIJOTE-PNG simulations, we quantify the information content accessible with measurements of the halo power spectrum monopole and quadrupole, the matter power spectrum, the halo-matter cross spectrum and the halo bispectrum monopole. This analysis is… ▽ More We investigate how much can be learnt about four types of primordial non-Gaussianity (PNG) from small-scale measurements of the halo field. Using the QUIJOTE-PNG simulations, we quantify the information content accessible with measurements of the halo power spectrum monopole and quadrupole, the matter power spectrum, the halo-matter cross spectrum and the halo bispectrum monopole. This analysis is the first to include small, non-linear scales, up to $k_\mathrm{max}=0.5 \mathrm{h/Mpc}$, and to explore whether these scales can break degeneracies with cosmological and nuisance parameters making use of thousands of N-body simulations. We perform all the halo measurements in redshift space with a single sample comprised of all halos with mass $>3.2 \times 10^{13}~h^{-1}M_\odot$. For local PNG, measurements of the scale dependent bias effect from the power spectrum using sample variance cancellation provide significantly tighter constraints than measurements of the halo bispectrum. In this case measurements of the small scales add minimal additional constraining power. In contrast, the information on equilateral and orthogonal PNG is primarily accessible through the bispectrum. For these shapes, small scale measurements increase the constraining power of the halo bispectrum by up to $\times4$, though the addition of scales beyond $k\approx 0.3 \mathrm{h/Mpc}$ improves constraints largely through reducing degeneracies between PNG and the other parameters. These degeneracies are even more powerfully mitigated through combining power spectrum and bispectrum measurements. However even with combined measurements and small scale information, equilateral non-Gaussianity remains highly degenerate with $σ_8$ and our bias model. △ Less

Submitted 20 December, 2022; v1 submitted 30 June, 2022; originally announced June 2022.

Comments: Updated to accepted version

Quijote-PNG: Quasi-maximum likelihood estimation of Primordial Non-Gaussianity in the non-linear dark matter density field

Authors: Gabriel Jung, Dionysios Karagiannis, Michele Liguori, Marco Baldi, William R Coulton, Drew Jamieson, Licia Verde, Francisco Villaescusa-Navarro, Benjamin D. Wandelt

Abstract: Future Large Scale Structure surveys are expected to improve over current bounds on primordial non-Gaussianity (PNG), with a significant impact on our understanding of early Universe physics. The level of such improvements will however strongly depend on the extent to which late time non-linearities erase the PNG signal on small scales. In this work, we show how much primordial information remains… ▽ More Future Large Scale Structure surveys are expected to improve over current bounds on primordial non-Gaussianity (PNG), with a significant impact on our understanding of early Universe physics. The level of such improvements will however strongly depend on the extent to which late time non-linearities erase the PNG signal on small scales. In this work, we show how much primordial information remains in the bispectrum of the non-linear dark matter density field by implementing a new, simulation-based, methodology for joint estimation of PNG amplitudes ($f_{\rm NL}$) and standard $Λ$CDM parameters. The estimator is based on optimally compressed statistics, which, for a given input density field, combine power spectrum and modal bispectrum measurements, and numerically evaluate their covariance and their response to changes in cosmological parameters. We train and validate the estimator using a large suite of N-body simulations (QUIJOTE-PNG), including different types of PNG (local, equilateral, orthogonal). We explicitly test the estimator's unbiasedness, optimality and stability with respect to changes in the total number of input realizations. While the dark matter power spectrum itself contains negligible PNG information, as expected, including it as an ancillary statistic increases the PNG information content extracted from the bispectrum by a factor of order $2$. As a result, we prove the capability of our approach to optimally extract PNG information on non-linear scales beyond the perturbative regime, up to $k_{\rm max} = 0.5~h\,{\rm Mpc}^{-1}$, obtaining marginalized $1$-$σ$ bounds of $Δf_{\rm NL}^{\rm local} \sim 16$, $Δf_{\rm NL}^{\rm equil} \sim 77$ and $Δf_{\rm NL}^{\rm ortho} \sim 40$ on a cubic volume of $1~(\mathrm{Gpc}/h)^3$ at $z=1$. At the same time, we discuss the significant information on cosmological parameters contained on these scales. △ Less

Submitted 3 June, 2022; originally announced June 2022.

Comments: 22 pages, 12 figures

Journal ref: Astrophys.J. 940 (2022) 1, 71

Quijote-PNG: Simulations of primordial non-Gaussianity and the information content of the matter field power spectrum and bispectrum

Authors: William R Coulton, Francisco Villaescusa-Navarro, Drew Jamieson, Marco Baldi, Gabriel Jung, Dionysios Karagiannis, Michele Liguori, Licia Verde, Benjamin D. Wandelt

Abstract: Primordial non-Gaussianity (PNG) is one of the most powerful probes of the early Universe and measurements of the large scale structure of the Universe have the potential to transform our understanding of this area. However relating measurements of the late time Universe to the primordial perturbations is challenging due to the non-linear processes that govern the evolution of the Universe. To hel… ▽ More Primordial non-Gaussianity (PNG) is one of the most powerful probes of the early Universe and measurements of the large scale structure of the Universe have the potential to transform our understanding of this area. However relating measurements of the late time Universe to the primordial perturbations is challenging due to the non-linear processes that govern the evolution of the Universe. To help address this issue we release a large suite of N-body simulations containing four types of PNG: \textsc{quijote-png}. These simulations were designed to augment the \textsc{quijote} suite of simulations that explored the impact of various cosmological parameters on large scale structure observables. Using these simulations we investigate how much information on PNG can be extracted by extending power spectrum and bispectrum measurements beyond the perturbative regime at $z=0.0$. This is the first joint analysis of the PNG and cosmological information content accessible with power spectrum and bispectrum measurements of the non-linear scales. We find that the constraining power improves significantly up to $k_\mathrm{max}\approx 0.3 h/{\rm Mpc}$, with diminishing returns beyond as the statistical probes signal-to-noise ratios saturate. This saturation emphasizes the importance of accurately modelling all the contributions to the covariance matrix. Further we find that combining the two probes is a powerful method of breaking the degeneracies with the $Λ$CDM parameters. △ Less

Submitted 26 May, 2023; v1 submitted 3 June, 2022; originally announced June 2022.

Comments: The simulation products are available at https://quijote-simulations.readthedocs.io/en/latest/png.html

Position-Dependent Correlation Function of Weak Lensing Convergence

Authors: D. Munshi, G. Jung, T. D. Kitching, J. McEwen, M. Liguori, T. Namikawa, A. Heavens

Abstract: We provide a systematic study of the position-dependent correlation function in weak lensing convergence maps and its relation to the squeezed limit of the three-point correlation function (3PCF) using state-of-the-art numerical simulations. We relate the position-dependent correlation function to its harmonic counterpart, i.e., the position-dependent power spectrum or equivalently the integrated… ▽ More We provide a systematic study of the position-dependent correlation function in weak lensing convergence maps and its relation to the squeezed limit of the three-point correlation function (3PCF) using state-of-the-art numerical simulations. We relate the position-dependent correlation function to its harmonic counterpart, i.e., the position-dependent power spectrum or equivalently the integrated bispectrum. We use a recently proposed improved fitting function, BiHalofit, for the bispectrum to compute the theoretical predictions as a function of source redshifts. In addition to low redshift results ($z_s=1.0-2.0$), we also provide results for maps inferred from lensing of the cosmic microwave background, i.e., $z_s=1100$. We include a {\em Euclid}-type realistic survey mask and noise. In agreement with the recent studies on the position-dependent power spectrum, we find that the results from simulations are consistent with the theoretical expectations when appropriate corrections are included. Performing a rough estimate, we find that the (S/N) for the detection of the position-dependent correlation function from {\em Euclid}-type mask with $f_{sky}=0.35$, can range between $6-12$ depending on the value of the intrinsic ellipticity distribution parameter $σ_ε = 0.3-1.0$. For reconstructed $κ$ maps using an ideal CMB survey the (S/N) $\approx 1.8$. We also found that a $10\%$ deviation in $σ_8$ can be detected using IB for the optimistic case of $σ_ε=0.3$ with a (S/N) $\approx 5$. The (S/N) for such detection in case of $Ω_M$ is lower. △ Less

Submitted 17 January, 2023; v1 submitted 2 April, 2021; originally announced April 2021.

Comments: 7 pages, 7 figures (PRD in press)

The integrated angular bispectrum of weak lensing

Authors: Gabriel Jung, Toshiya Namikawa, Michele Liguori, Dipak Munshi, Alan Heavens

Abstract: We investigate three-point statistics in weak lensing convergence, through the integrated bispectrum. This statistic involves measuring power spectra in patches, and is thus easy to measure, and avoids the complexity of estimating the very large number of possible bispectrum configurations. The integrated bispectrum principally probes the squeezed limit of the bispectrum. To be useful as a set of… ▽ More We investigate three-point statistics in weak lensing convergence, through the integrated bispectrum. This statistic involves measuring power spectra in patches, and is thus easy to measure, and avoids the complexity of estimating the very large number of possible bispectrum configurations. The integrated bispectrum principally probes the squeezed limit of the bispectrum. To be useful as a set of summary statistics, accurate theoretical predictions of the signal are required, and, assuming Gaussian sampling distributions, the covariance matrix. In this paper, we investigate through simulations how accurate are theoretical formulae for both the integrated bispectrum and its covariance, finding that there a small inaccuracies in the theoretical signal, and more serious deviations in the covariance matrix, which may need to be estimated using simulations. △ Less

Submitted 12 July, 2021; v1 submitted 10 February, 2021; originally announced February 2021.

Comments: 23 pages, 12 figures. v2: minor textual corrections and additions, matches published version

Journal ref: JCAP06(2021)055

The integrated angular bispectrum

Authors: Gabriel Jung, Filippo Oppizzi, Andrea Ravenni, Michele Liguori

Abstract: We study the position-dependent power spectrum and the integrated bispectrum statistic for 2D cosmological fields on the sphere (integrated angular bispectrum). First, we derive a useful, $m$-independent, formula for the full-sky integrated angular bispectrum, based on the construction of azimuthally symmetric patches. We then implement a pipeline for integrated angular bispectrum estimation, incl… ▽ More We study the position-dependent power spectrum and the integrated bispectrum statistic for 2D cosmological fields on the sphere (integrated angular bispectrum). First, we derive a useful, $m$-independent, formula for the full-sky integrated angular bispectrum, based on the construction of azimuthally symmetric patches. We then implement a pipeline for integrated angular bispectrum estimation, including a mean-field correction to account for spurious isotropy-breaking effects in realistic conditions (e.g., inhomogenous noise, sky masking). Finally, we show examples of applications of this estimator to CMB analysis, both using simulations and actual Planck data. Such examples include $f_\mathrm{NL}$ estimation, analyses of non-Gaussianity from secondary anisotropies (ISW-lensing and ISW-tSZ-tSZ bispectra) and studies of non-Gaussian signatures from foreground contamination. △ Less

Submitted 7 April, 2020; originally announced April 2020.

Comments: 38 pages, 14 figures

Journal ref: JCAP 06 (2020) 035

Planck 2018 results. IX. Constraints on primordial non-Gaussianity

Authors: Planck Collaboration, Y. Akrami, F. Arroja, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. R. Bond, J. Borrill, F. R. Bouchet, M. Bucher, C. Burigana, R. C. Butler, E. Calabrese, J. -F. Cardoso, B. Casaponsa, A. Challinor , et al. (135 additional authors not shown)

Abstract: We analyse the Planck full-mission cosmic microwave background (CMB) temperature and E-mode polarization maps to obtain constraints on primordial non-Gaussianity (NG). We compare estimates obtained from separable template-fitting, binned, and modal bispectrum estimators, finding consistent values for the local, equilateral, and orthogonal bispectrum amplitudes. Our combined temperature and polariz… ▽ More We analyse the Planck full-mission cosmic microwave background (CMB) temperature and E-mode polarization maps to obtain constraints on primordial non-Gaussianity (NG). We compare estimates obtained from separable template-fitting, binned, and modal bispectrum estimators, finding consistent values for the local, equilateral, and orthogonal bispectrum amplitudes. Our combined temperature and polarization analysis produces the following results: f_NL^local = -0.9 +\- 5.1; f_NL^equil = -26 +\- 47; and f_NL^ortho = - 38 +\- 24 (68%CL, statistical). These results include the low-multipole (4 <= l < 40) polarization data, not included in our previous analysis, pass an extensive battery of tests, and are stable with respect to our 2015 measurements. Polarization bispectra display a significant improvement in robustness; they can now be used independently to set NG constraints. We consider a large number of additional cases, e.g. scale-dependent feature and resonance bispectra, isocurvature primordial NG, and parity-breaking models, where we also place tight constraints but do not detect any signal. The non-primordial lensing bispectrum is detected with an improved significance compared to 2015, excluding the null hypothesis at 3.5 sigma. We present model-independent reconstructions and analyses of the CMB bispectrum. Our final constraint on the local trispectrum shape is g_NLl^local = (-5.8 +\-6.5) x 10^4 (68%CL, statistical), while constraints for other trispectra are also determined. We constrain the parameter space of different early-Universe scenarios, including general single-field models of inflation, multi-field and axion field parity-breaking models. Our results provide a high-precision test for structure-formation scenarios, in complete agreement with the basic picture of the LambdaCDM cosmology regarding the statistics of the initial conditions (abridged). △ Less

Submitted 14 May, 2019; originally announced May 2019.

Comments: 50 pages, 20 figures

The bispectra of galactic CMB foregrounds and their impact on primordial non-Gaussianity estimation

Authors: Gabriel Jung, Benjamin Racine, Bartjan van Tent

Abstract: We use the binned bispectrum estimator to determine the bispectra of the dust, free-free, synchrotron, and AME galactic foregrounds using maps produced by the Commander component separation method from Planck 2015 data. We find that all of these peak in the squeezed configuration, allowing for potential confusion with in particular the local primordial shape. Applying an additional functionality i… ▽ More We use the binned bispectrum estimator to determine the bispectra of the dust, free-free, synchrotron, and AME galactic foregrounds using maps produced by the Commander component separation method from Planck 2015 data. We find that all of these peak in the squeezed configuration, allowing for potential confusion with in particular the local primordial shape. Applying an additional functionality implemented in the binned bispectrum estimator code, we then use these galactic bispectra as templates in an $f_\mathrm{NL}$ analysis of other maps. After testing and validating the method and code with simulations, we show that we detect the dust in the raw 143 GHz map with the expected amplitude (the other galactic foregrounds are too weak at 143 GHz to be detected) and that no galactic residuals are detected in the cleaned CMB map. We also investigate the effect of the mask on the templates and the effect of the choice of binning on a joint dust-primordial $f_\mathrm{NL}$ analysis. △ Less

Submitted 3 October, 2018; originally announced October 2018.

Comments: 42 pages, 20 figures

Report number: LPT-Orsay-18-87

Journal ref: JCAP11(2018)047

Non-Gaussianity in two-field inflation beyond the slow-roll approximation

Authors: Gabriel Jung, Bartjan van Tent

Abstract: We use the long-wavelength formalism to investigate the level of bispectral non-Gaussianity produced in two-field inflation models with standard kinetic terms. Even though the Planck satellite has so far not detected any primordial non-Gaussianity, it has tightened the constraints significantly, and it is important to better understand what regions of inflation model space have been ruled out, as… ▽ More We use the long-wavelength formalism to investigate the level of bispectral non-Gaussianity produced in two-field inflation models with standard kinetic terms. Even though the Planck satellite has so far not detected any primordial non-Gaussianity, it has tightened the constraints significantly, and it is important to better understand what regions of inflation model space have been ruled out, as well as prepare for the next generation of experiments that might reach the important milestone of Delta f_NL(local) = 1. We derive an alternative formulation of the previously derived integral expression for f_NL, which makes it easier to physically interpret the result and see which types of potentials can produce large non-Gaussianity. We apply this to the case of a sum potential and show that it is very difficult to satisfy simultaneously the conditions for a large f_NL and the observational constraints on the spectral index n_s. In the case of the sum of two monomial potentials and a constant we explicitly show in which small region of parameter space this is possible, and we show how to construct such a model. Finally, the new general expression for f_NL also allows us to prove that for the sum potential the explicit expressions derived within the slow-roll approximation remain valid even when the slow-roll approximation is broken during the turn of the field trajectory (as long as only the epsilon slow-roll parameter remains small). △ Less

Submitted 11 May, 2017; v1 submitted 28 November, 2016; originally announced November 2016.

Comments: 60 pages, 11 figures. v2: Added appendix B on product potentials, minor textual corrections and additions, matches published version

Report number: LPT-Orsay-16-82

Journal ref: JCAP05(2017)019

The Spectrum of Diffuse Cosmic Hard X-Rays Measured with HEAO-1

Authors: D. E. Gruber, J. L. Matteson, L. E. Peterson, G. V. Jung

Abstract: The spectrum of the diffuse isotropic component of cosmic X-rays over the 13-180 keV range was determined by the UCSD/MIT Hard X-Ray and Gamma-Ray instrument (HEAO A4) on the High Energy Astronomical Observatory-1 (HEAO-1). The instrument consists of a complex of actively shielded and collimated scintillation counters, including the Low Energy Detector set from which the data reported here were… ▽ More The spectrum of the diffuse isotropic component of cosmic X-rays over the 13-180 keV range was determined by the UCSD/MIT Hard X-Ray and Gamma-Ray instrument (HEAO A4) on the High Energy Astronomical Observatory-1 (HEAO-1). The instrument consists of a complex of actively shielded and collimated scintillation counters, including the Low Energy Detector set from which the data reported here were obtained. These data join smoothly with the spectrum at lower energies reported by the GSFC HEAO A2 instrument and with that measured to 400 keV by the HEAO A4 Medium Energy Detectors. The HEAO data set also joins the recent results from COMPTEL on the Compton Gamma-Ray Observatory in the 1-10 MeV range, which failed to confirm the existence of an "MeV bump" in this range. Although the spectrum over the entire range 3 keV < E < 100 GeV can be fit by a simple empirical analytic expression, the origin is likely due to a number of distinct source components. The prevailing idea for the origin is that the hard X-ray spectrum is due to X-rays from various AGN components, particularly Seyfert galaxies extending to cosmological distances, and that the low energy gamma-rays may be due to emission from type 1a Supernovae, also integrated to cosmological distances. The higher energy gamma-ray spectrum defined by EGRET, also on the CGRO, may be due to unresolved gamma-ray emitting blazars. Models of production by these source components, extrapolated to the present epoch, must reproduce the observationally derived spectrum. △ Less

Submitted 31 March, 1999; originally announced March 1999.

Comments: LaTeX, 16pp, 3figs

Report number: SP-98-25

OSSE Observations of the Soft Gamma Ray Continuum from the Galactic Plane at Longitude 95 Degrees

Authors: J. G. Skibo, W. N. Johnson, J. D. Kurfess, R. L. Kinzer, G. Jung, J. E. Grove, W. R. Purcell, M. P. Ulmer, N. Gehrels, J. Tueller

Abstract: We present the results of OSSE observations of the soft gamma ray continuum emission from the Galactic plane at longitude 95 degrees. Emission is detected between 50 and 600 keV where the spectrum is fit well by a power law with photon index -2.6+-0.3 and flux (4.0+-0.5) 10^{-2} photons/s/cm^2/rad/MeV at 100 keV. This spectral shape in this range is similar to that found for the continuum emissi… ▽ More We present the results of OSSE observations of the soft gamma ray continuum emission from the Galactic plane at longitude 95 degrees. Emission is detected between 50 and 600 keV where the spectrum is fit well by a power law with photon index -2.6+-0.3 and flux (4.0+-0.5) 10^{-2} photons/s/cm^2/rad/MeV at 100 keV. This spectral shape in this range is similar to that found for the continuum emission from the inner Galaxy but the amplitude is lower by a factor of four. This emission is either due to unresolved and previously unknown point sources or it is of diffuse origin, or a combination of the two. Simultaneous observations with OSSE and smaller field of view instruments operating in the soft gamma ray energy band, such as XTE or SAX, would help resolve this issue. If it is primarily diffuse emission due to nonthermal electron bremsstrahlung, as is the >1 MeV Galactic ridge continuum, then the power in low energy cosmic ray electrons exceeds that of the nuclear component of the cosmic rays by an order of magnitude. This would have profound implications for the origin of cosmic rays and the energetics of the interstellar medium. Alternatively, if the emission is diffuse and thermal, then there must be a component of the interstellar medium at temperatures near 10^9 K. △ Less

Submitted 21 April, 1997; originally announced April 1997.

Comments: 11 pages, Latex, requires AASTEX macros and psfig.tex, 2 postscript figures, Accepted for publication in the Astrophysical Journal Letters