$\texttt{GENGARS}$: Accurate non-Gaussian initial conditions with arbitrary bispectrum for N-body simulations

Benjamin Dan Wandelt; Emanuele Fondi; Francisco Villaescusa-Navarro; Licia Verde; Marco Baldi; William Coulton

REVIEW 3 major objections 4 minor 65 references

$\texttt{GENGARS}$: Accurate non-Gaussian initial conditions with arbitrary bispectrum for N-body simulations

T0 review · 3 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read GENGARS generates non-Gaussian initial conditions for N-body simulations whose primordial potential has any separable target bispectrum, at a cost set by FFTs rather than brute-force O(N^6) integrals.

desk verdict GENGARS is a practical, well-benchmarked realization of the Wagner-Verde kernel via Schwinger parameterization, but a central identity typo and an overstated 'arbitrary' claim need fixing before publication. read the letter →

arxiv 2508.01855 v1 pith:PVM3ISVO submitted 2025-08-03 astro-ph.CO gr-qc

Emanuele Fondi , Licia Verde , Marco Baldi , William Coulton , Francisco Villaescusa-Navarro , Benjamin Dan Wandelt This is my paper

classification astro-ph.COgr-qc

keywords primordialnon-GaussianitybispectrumN-bodysimulationsinitialconditionsSchwingerparameterizationseparablekernelslarge-scalestructure2LPT-PNG

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper presents GENGARS, a code that generates non-Gaussian initial conditions for N-body simulations whose primordial gravitational potential has a specified separable bispectrum. The central move is to evaluate the reduced-bispectrum kernel previously considered too expensive for high-resolution runs through a Schwinger parameter integral, which turns the kernel's non-separable denominator into a product of one-dimensional factors and makes the whole generation a set of FFT convolutions. For the local, equilateral, and orthogonal templates, the paper shows the generated fields match the target bispectrum while introducing a non-Gaussian correction to the primordial power spectrum that is orders of magnitude smaller than in the standard 2LPT-PNG code, with lower realization-to-realization scatter. This matters because it makes non-standard bispectrum shapes, such as oscillatory feature models, practical to test with simulations, and it removes an unphysical large-scale power excess that the conventional method injects for orthogonal-type shapes.

What carries the argument

The load-bearing object is the reduced-bispectrum kernel $W(k_1,k_2,k_3)=B_\Phi(k_1,k_2,k_3)/\left(2f_{\rm NL}[P_\Phi(k_1)P_\Phi(k_2)+P_\Phi(k_2)P_\Phi(k_3)+P_\Phi(k_1)P_\Phi(k_3)]\right)$, which uniquely fixes the generating kernel and suppresses the squeezed limit that drives spurious power. The Schwinger identity $1/(P_1P_2+P_2P_3+P_1P_3)=\int_0^\infty dt\,\prod_{j=1}^3 e^{-t/P_\Phi(k_j)}/P_\Phi(k_j)$ rewrites the denominator as a separable product, so the non-Gaussian potential becomes a $t$-integral of terms each computable as a convolution via FFTs. Discretizing $t$ with $N_t\sim300$-$400$ steps yields $O(N_t N_i N_g^3 \log N_g)$ scaling instead of the brute-force $O(N_g^6)$.

What would settle it

For a non-separable target shape, compare the truncated Schwinger reconstruction against the brute-force reduced-bispectrum kernel of Eq. (2.11) on the same grid: if the relative difference in the recovered bispectrum or in $P^{\rm NG}_\Phi(k)$ does not shrink as $N_t$ grows from 100 to 1000, the method's accuracy claim is confined to exactly separable templates rather than arbitrary bispectra.

Watch

Extended reading notes

Core claim

The central claim is that the inverse problem of generating a non-Gaussian field with a target bispectrum can be solved accurately and efficiently by the reduced-bispectrum kernel, once that kernel is rendered separable through the Schwinger identity. On the paper's own terms, the method reproduces the mean primordial bispectrum for local, equilateral, and orthogonal templates within the target, while suppressing the spurious non-Gaussian power-spectrum contribution below three orders of magnitude relative to the Gaussian spectrum at $f_{\rm NL}=100$, and to below $0.1\%$ in the evolved large-scale matter power at $z=0$, whereas the comparison code leaves a percent-level artifact for the orthogonal template. The authors also demonstrate the flexibility by generating an oscillatory bispectrum shape without any template-specific derivation.

Load-bearing premise

The whole construction presumes the target bispectrum is exactly separable, and that a few hundred Schwinger integration steps reproduce the true reduced kernel; the paper checks convergence empirically but gives no formal error bound for shapes whose separable decomposition is not exact.

Editorial extensions

If this is right

For any separable bispectrum template, initial conditions can be generated by specifying only the functional form of $B_\Phi$; no per-shape derivation of a generating kernel is needed.
The spurious non-Gaussian contribution to the primordial power spectrum is suppressed relative to 2LPT-PNG, so constraints on the primordial tilt are not contaminated at the percent level for $f_{\rm NL}$ values of order a few hundred.
Realization-to-realization scatter in the measured primordial bispectrum is smaller with GENGARS than with 2LPT-PNG, particularly in squeezed configurations.
At $z=0$, the large-scale matter power spectrum from GENGARS initial conditions shows no percent-level unphysical excess for orthogonal PNG, while the halo mass function is consistent between the two methods.
Computational cost remains higher than 2LPT-PNG (minutes to about an hour on 64 CPUs for $512^3$ grids, versus seconds), so 2LPT-PNG remains preferable for large production suites of standard templates.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

For bispectra that are not exactly separable, the Schwinger-truncated kernel is only an approximation; one could extend the method by numerically factorizing non-separable shapes (for example via eigen-decomposition of the reduced kernel on the triangle grid) and quoting the approximation error.
The paper controls only the three-point function; its lower variance may reflect a reduced six-point function, which could be tested by measuring the connected four-point function (trispectrum) of the generated potentials and comparing with 2LPT-PNG.
The absence of a formal convergence bound for $N_t$ suggests a practical test: for a given shape, one can increase $N_t$ until the measured bispectrum and power-spectrum correction stop changing, and use the residual as a systematic error in simulation analyses.
The implementation could be accelerated with GPUs or by adapting the Orszag 3/2 rule as an optional anti-aliasing choice, which the authors state is planned.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Short version: GENGARS makes the Wagner-Verde reduced-bispectrum kernel computationally practical by rewriting the denominator as a Schwinger integral, and the benchmarks show it reproduces the target bispectrum with less spurious power-spectrum contamination than 2LPT-PNG. It's a genuinely useful tool for the PNG simulation community. But there's a typo in the central identity, and the 'arbitrary bispectrum' framing is too strong—the method is exact only for finite separable sums.

What's new: the specific combination of the reduced-bispectrum kernel (Wagner & Verde 2012) with the Schwinger parameterization (Smith & Zaldarriaga 2011) is not in the earlier literature. The paper demonstrates it on the standard templates and an oscillatory shape, with a clean benchmark design: shared Gaussian fields, 10 realizations, and both initial and evolved statistics. The power-spectrum suppression is consistently demonstrated, and the authors are careful to say that lower variance does not necessarily mean better accuracy because the six-point function is not controlled. That's honest and scientifically sound. The evolved-field comparison with Quijote-PNG is also transparent about unresolved small-scale differences.

The soft spots are addressable. First, Eq. (2.13), the key identity, is misprinted: as written, the integrand has factors of P(k) that should be 1/P(k). The correct form follows from 1/(P1P2+P2P3+P1P3) = ∫ dt [e^{-t/P1}/P1][e^{-t/P2}/P2][e^{-t/P3}/P3], and the typo likely propagates to Eq. (2.16). It needs fixing before the paper can be used as a reference. Second, the title and intro overstate with 'arbitrary bispectrum.' The paper's own Sec. 2.3 restricts to separable bispectra, but the abstract and title should match that scope. For genuinely non-separable shapes, no error control or approximation strategy is provided. Third, the N_t convergence is calibrated empirically without a formal bound—acceptable for a code paper, but it should be stated as a limitation, not a proof of convergence. Finally, the code is promised but not yet public, so independent reproduction isn't possible right now.

Who this is for: anyone running N-body simulations with non-Gaussian initial conditions, particularly for non-standard separable bispectra. The paper deserves a serious referee; with the typo corrected and the language tempered, it will be a solid methods contribution.

Referee Report

3 major / 4 minor

Summary. GENGARS is a code for generating non-Gaussian initial conditions for N-body simulations by implementing the reduced-bispectrum kernel of Wagner & Verde (Eq. 2.11) through a Schwinger parameterization of the denominator (Eqs. 2.13-2.16). For target bispectra that are finite separable sums (Eq. 2.14), the non-Gaussian potential can be computed with FFTs at a cost of O(N_t N_i N_g^3 log N_g). The paper benchmarks against 2LPT-PNG for local, equilateral and orthogonal templates using a shared Gaussian field and 10 realizations, finding unbiased bispectrum recovery with a smaller spurious contribution to the primordial power spectrum and lower realization variance. It demonstrates an oscillatory bispectrum example and compares z=0 matter power spectrum, bispectrum and halo mass function from 10 Quijote-like simulations.

Significance. If the central claims hold, GENGARS is a useful and timely tool: it makes the Wagner-Verde kernel computationally practical, provides a parameter-free kernel choice that reduces a known systematic in 2LPT-PNG initial conditions, and opens the door to simulation studies of feature bispectra. The benchmark design is a clear strength: a shared Gaussian field isolates kernel differences, 10 realizations give error-on-mean estimates, and both initial and evolved statistics are checked. The code is to be made public, which supports reproducibility. The main caveat is that the demonstrated exactness applies to separable bispectra; the 'arbitrary' claims need qualification.

major comments (3)

[Title, Abstract, §2.3, §5] The phrase 'arbitrary bispectrum' in the title and abstract, and the sentence in §5 that GENGARS allows reproducing 'arbitrary shapes without requiring any approximation or analytical treatment beyond specifying the functional form', overstate the scope of the derivation. Exactness requires the target bispectrum to be a finite separable sum, Eq. (2.14); for a non-separable target the user must supply a separable approximation, and the paper gives no prescription, convergence criterion, or error bound for that approximation. The oscillatory template in Eq. (5.1) is itself a finite separable sum, so it does not test the non-separable case. I recommend restricting the claims to separable bispectra throughout (title, abstract, conclusions), or adding an error-controlled approximation layer for non-separable shapes.
[§4.3, Eq. (2.15)] The only numerical approximation in the method is the discretization of the Schwinger t-integral, yet the convergence statement in §4.3 ('N_t ∼ 300–400 steps are sufficient') is not supported by any shown convergence test, quadrature rule, or choice of upper cutoff for the infinite integration range. The comparison described in the text is against the brute-force kernel of Eq. (2.11) for exactly separable templates only, so it does not quantify the error for non-separable targets and does not show how N_t should be chosen for a new shape. Please provide a convergence plot or table (e.g., relative error in the reconstructed bispectrum and in P^NG_Φ as a function of N_t) and specify the integration scheme, so the accuracy claim is reproducible.
[§6.1, §6.2] The claim of 'improved accuracy' of the evolved field is only established on large scales. At small scales the paper states that it is 'not straightforward to assess which implementation is more accurate' and that a 1-loop bispectrum would be needed; the orthogonal matter power spectrum shows scale-dependent differences between the two prescriptions. Since the abstract presents improved accuracy as a headline result, the authors should either restrict that claim to the initial conditions and the linear/large-scale evolved field, or add a quantitative reference calculation that decides which kernel gives the correct small-scale limit.

minor comments (4)

[Eq. (2.16)] The definition of w_i^j(t,k) is typeset ambiguously: it should read w_i^j(t,k) = b_i^j(k) e^{-t/P_Φ(k)} / P_Φ(k), not b_i^j(k) e^{-t/P_Φ(k)} P_Φ(k). Please clarify the notation so that the integral over t reproduces Eq. (2.11).
[§4.1, Fig. 2] The text refers to k^{n_s-4} P_Φ^NG(k), while the figure label uses k^{4-n_s} P_Φ^NG. The latter is the natural choice for isolating the scaling relative to the Gaussian power spectrum; please make text and figure consistent.
[§5, Eq. (5.1)] Please state explicitly that the oscillatory template of Eq. (5.1) is a finite separable sum (the sine of a sum decomposes into products of sines and cosines); otherwise the example may be read as evidence for non-separable shapes.
[§3, Abstract] Minor typos and presentation issues: 'GEnerator of Non-Gaussian ARbitrary Shapes' has an extra space in the expansion, and 'actualisation' in §3 should be 'actualization' or 'updated implementation'. Also, the repository statement would be stronger with a URL or DOI at submission.

Circularity Check

0 steps flagged · score 0.0 of 10

Self-contained construction: the target bispectrum enters through the explicitly defined reduced kernel, and the reported power-spectrum suppression is a derived, non-fitted consequence; no circular step found.

full rationale

GENGARS's central claim is to generate a primordial potential whose bispectrum matches a user-specified separable target. This is an algorithm-construction goal, not a first-principles prediction. The generating kernel is explicitly defined from the target bispectrum in Eq. (2.11), W = B_Phi/(2 f_NL [P1P2 + P2P3 + P1P3]), and inserting this W into Eq. (2.6) makes the leading-order bispectrum equal to the target by algebra. The agreement shown in Figs. 3 and 4 is therefore a numerical consistency check of the implementation and of the Schwinger-parameter discretization, not an independent empirical prediction; this is the intended design, not a hidden circularity. The paper's more substantive claim, the suppression of spurious primordial power-spectrum corrections, is derived from the kernel structure through Eq. (2.10), with the squeezed-limit suppression following from the denominator of Eq. (2.11), and it is checked against the Gaussian baseline and against the independent 2LPT-PNG code without fitting any parameter to the measured P_NG. The cited prior work of Wagner & Verde is background for the kernel choice, but the key relations are restated and used transparently, so the self-citation is not load-bearing. The acknowledged limitations, namely the requirement that the target bispectrum be separable as in Eq. (2.14) and the empirical rather than formal convergence test for N_t (Sec. 4.3), concern approximation error and scope of applicability, not circularity. No equation or fitted parameter is renamed as a prediction, and no result is forced by a self-citation chain. The derivation chain is therefore self-contained and non-circular.

Assumptions & free parameters 1 free parameters · 3 assumptions · 0 invented entities

No new physical entities or fitted model parameters. The only numerical knob is N_t, the number of Schwinger integration steps, which is a convergence parameter rather than a fit. The method assumes separability of the input bispectrum and leading-order PNG, standard in this literature.

free parameters (1)

N_t (Schwinger integration steps) = 300-400
Chosen by testing convergence against brute-force computation of Eq. (2.5) for the reduced kernel; not fitted to match the target bispectrum, but a numerical parameter affecting accuracy.

assumptions (3)

domain assumption Target bispectrum is separable (Eq. 2.14)
The method requires BPhi to be written as a sum of products of one-dimensional functions; many physical shapes must be approximated.
domain assumption Leading-order expansion Phi = PhiG + fNL PhiNG (Eq. 2.5)
The kernel is derived from the tree-level bispectrum; higher-order correlators and fNL^2 corrections are not controlled. This is standard in the field but is an assumption.
standard math Schwinger identity (Eq. 2.12) and its separable application (Eq. 2.13)
The identity 1/f(x) = integral e^{-t f(x)} dt is standard; the specific application uses the substitution that makes the denominator separable. Note: as typeset in the preprint, Eq. (2.13) appears to have an incorrect exponent (e^{-t P} instead of e^{-t/P}), likely a typo; the corrected identity is needed for the derivation to hold.

how reviews work

0 comments

Cite this review

Pith. "Pith review of $\texttt{GENGARS}$: Accurate non-Gaussian initial conditions with arbitrary bispectrum for N-body simulations." pith.science (2026). https://pith.science/paper/PVM3ISVO

@misc{pith2026250801855,
  author       = {Pith},
  title        = {Pith review of: $\textttGENGARS$: Accurate non-Gaussian initial conditions with arbitrary bispectrum for N-body simulations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/PVM3ISVO}},
  note         = {Machine review of arXiv:2508.01855}
}

abstract

Primordial non-Gaussianity is predicted by various inflationary models, and N-body simulations are a crucial tool for studying its imprints on large-scale structure. In this work, we present \texttt{GENGARS} ( GEnerator of Non-Gaussian ARbitrary Shapes), a framework for generating accurate non-Gaussian initial conditions for N-body simulations. It builds upon the formulation introduced by Wagner \& Verde (2012), enabling to generate a primordial gravitational potential with a desired separable bispectrum $B_{\Phi}(k_1,k_2,k_3)$. For the local, equilateral and orthogonal non-Gaussian templates, we benchmark our method against the well-established \texttt{2LPT-PNG} code. We show that \texttt{GENGARS} achieves improved accuracy and lower noise by suppressing spurious contributions to the primordial power spectrum. This paper aims at presenting the method, quantifying its performance and illustrating the benefits and applicable use cases over existing approaches.

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

65 extracted references · 12 canonical work pages

[1]

Maldacena, Non-Gaussian features of primordial fluctuations in single field inflationary models , Journal of High Energy Physics 2003 (2003) 013

J. Maldacena, Non-Gaussian features of primordial fluctuations in single field inflationary models , Journal of High Energy Physics 2003 (2003) 013

work page 2003
[2]

Acquaviva, N

V. Acquaviva, N. Bartolo, S. Matarrese and A. Riotto, Gauge-invariant second-order perturbations and non-Gaussianity from inflation , Nuclear Physics B 667 (2003) 119 [ astro-ph/0209156]

arXiv 2003
[3]

Creminelli and M

P. Creminelli and M. Zaldarriaga, Single field consistency relation for the 3-point function , Journal of Cosmology and Astroparticle Physics 2004 (2004) 006

work page 2004
[4]

Bartolo, E

N. Bartolo, E. Komatsu, S. Matarrese and A. Riotto, Non-Gaussianity from inflation: theory and observations, Phys. Rep. 402 (2004) 103 [ astro-ph/0406398]

arXiv 2004
[5]

Chen, Primordial Non-Gaussianities from Inflation Models , Advances in Astronomy 2010 (2010) 638979 [1002.1416]

X. Chen, Primordial Non-Gaussianities from Inflation Models , Advances in Astronomy 2010 (2010) 638979 [1002.1416]

arXiv 2010
[6]

Meerburg et al., Primordial Non-Gaussianity, 1903.04409

P.D. Meerburg et al., Primordial Non-Gaussianity, 1903.04409. – 16 –

arXiv 1903
[7]

Chen, M.-x

X. Chen, M.-x. Huang, S. Kachru and G. Shiu, Observational signatures and non-Gaussianities of general single-field inflation , JCAP 2007 (2007) 002 [ hep-th/0605045]

arXiv 2007
[8]

Chen and Y

X. Chen and Y. Wang, Large non-Gaussianities with intermediate shapes from quasi-single-field inflation, Phys. Rev. D 81 (2010) 063511 [ 0909.0496]

arXiv 2010

Show all 65 references

[9]

Chen and Y

X. Chen and Y. Wang, Quasi-single field inflation and non-Gaussianities , JCAP 2010 (2010) 027 [0911.3380]

2010 arXiv
[10]

Arkani-Hamed, P

N. Arkani-Hamed, P. Creminelli, S. Mukohyama and M. Zaldarriaga, Ghost inflation , JCAP 2004 (2004) 001 [ hep-th/0312100]

2004 arXiv
[11]

Alishahiha, E

M. Alishahiha, E. Silverstein and D. Tong, DBI in the Sky , Phys. Rev. D 70 (2004) 123505 [hep-th/0404084]

2004 arXiv
[12]

Langlois, S

D. Langlois, S. Renaux-Petel, D.A. Steer and T. Tanaka, Primordial perturbations and non-Gaussianities in DBI and general multifield inflation , Phys. Rev. D 78 (2008) 063523 [ 0806.0336]

2008 arXiv
[13]

Linde and V

A. Linde and V. Mukhanov, Non-Gaussian isocurvature perturbations from inflation , Phys. Rev. D 56 (1997) R535 [ astro-ph/9610219]

1997 arXiv
[14]

Cheung, A.L

C. Cheung, A.L. Fitzpatrick, J. Kaplan, L. Senatore and P. Creminelli, The effective field theory of inflation, Journal of High Energy Physics 2008 (2008) 014 [ 0709.0293]

2008 arXiv
[16]

Arkani-Hamed and J

N. Arkani-Hamed and J. Maldacena, Cosmological collider physics , arXiv e-prints (2015) [arXiv:1503.08043]

2015 arXiv
[17]

Ach´ ucarro et al.,Inflation: Theory and Observations , 2203.08128

A. Ach´ ucarro et al.,Inflation: Theory and Observations , 2203.08128

arXiv
[18]

Planck collaboration, Planck 2018 results. IX. Constraints on primordial non-Gaussianity , Astron. Astrophys. 641 (2020) A9 [ 1905.05697]

2020 arXiv
[19]

Alvarez, T

M. Alvarez, T. Baldauf, J.R. Bond, N. Dalal, R. de Putter, O. Dor´ e et al., Testing Inflation with Large Scale Structure: Connecting Hopes with Reality , arXiv e-prints (2014) arXiv:1412.4671 [ 1412.4671]

2014 arXiv
[20]

Desjacques and U

V. Desjacques and U. Seljak, Primordial non-Gaussianity from the large-scale structure , Classical and Quantum Gravity 27 (2010) 124011 [ 1003.5020]

2010 arXiv
[21]

Sailer, E

N. Sailer, E. Castorina, S. Ferraro and M. White, Cosmology at high redshift — a probe of fundamental physics, JCAP 12 (2021) 049 [ 2106.09713]

2021 arXiv
[22]

Shiveshwarkar, T

C. Shiveshwarkar, T. Brinckmann and M. Loverde, Constraining multi-field inflation using the SPHEREx all-sky survey power spectra , JCAP 05 (2024) 094 [ 2312.15038]

2024 arXiv
[23]

Chaussidon, C

E. Chaussidon, C. Y` eche, A. de Mattia, C. Payerne, P. McDonald, A.J. Ross et al., Constraining primordial non-Gaussianity with DESI 2024 LRG and QSO samples , arXiv e-prints (2024) arXiv:2411.17623 [2411.17623]

2024 arXiv
[25]

Karagiannis, A

D. Karagiannis, A. Lazanu, M. Liguori, A. Raccanelli, N. Bartolo and L. Verde, Constraining Primordial non-Gaussianity with Bispectrum and Power Spectum from Upcoming Optical and Radio Surveys, Monthly Notices of the Royal Astronomical Society 478 (2018) 1341

2018
[26]

Cabass, M.M

G. Cabass, M.M. Ivanov, O.H.E. Philcox, M. Simonovi´ c and M. Zaldarriaga, Constraints on Single-Field Inflation from the BOSS Galaxy Survey , Phys. Rev. Lett. 129 (2022) 021301 [ 2201.07238]

2022 arXiv
[27]

D’Amico, M

G. D’Amico, M. Lewandowski, L. Senatore and P. Zhang, Limits on primordial non-Gaussianities from BOSS galaxy-clustering data , 2201.11518

arXiv
[28]

Cabass, O.H.E

G. Cabass, O.H.E. Philcox, M.M. Ivanov, K. Akitsu, S.-F. Chen, M. Simonovi´ c et al., BOSS Constraints on Massive Particles during Inflation: The Cosmological Collider in Action , 2404.01894

arXiv
[29]

Coulton, F

W.R. Coulton, F. Villaescusa-Navarro, D. Jamieson, M. Baldi, G. Jung, D. Karagiannis et al., Quijote-PNG: Simulations of Primordial Non-Gaussianity and the Information Content of the Matter Field Power Spectrum and Bispectrum , Astrophys. J. 943 (2023) 64 [ 2206.01619]. – 17 –

2023 arXiv
[30]

Coulton, F

W.R. Coulton, F. Villaescusa-Navarro, D. Jamieson, M. Baldi, G. Jung, D. Karagiannis et al., Quijote PNG: The information content of the halo power spectrum and bispectrum , 2206.15450

arXiv
[31]

Andrews, J

A. Andrews, J. Jasche, G. Lavaux and F. Schmidt, Bayesian field-level inference of primordial non-Gaussianity using next-generation galaxy surveys , 2203.08838

arXiv
[32]

G. Jung, A. Ravenni, M. Baldi, W.R. Coulton, D. Jamieson, D. Karagiannis et al., Quijote-PNG: The Information Content of the Halo Mass Function , Astrophys. J. 957 (2023) 50 [ 2305.10597]

2023 arXiv
[33]

G. Jung, A. Ravenni, M. Liguori, M. Baldi, W.R. Coulton, F. Villaescusa-Navarro et al., Quijote-PNG: Optimizing the Summary Statistics to Measure Primordial Non-Gaussianity , Astrophys. J. 976 (2024) 109 [2403.00490]

2024 arXiv
[34]

Baldi, E

M. Baldi, E. Fondi, D. Karagiannis, L. Moscardini, A. Ravenni, W.R. Coulton et al., Cosmological simulations of scale-dependent primordial non-Gaussianity , JCAP 2024 (2024) 053 [ 2407.06641]

2024 arXiv
[35]

Goldstein, O.H.E

S. Goldstein, O.H.E. Philcox, E. Fondi and W.R. Coulton, Wonderings on Wiggly Bispectra: Non-linear Evolution and Reconstruction of Oscillations in the Squeezed Bispectrum , arXiv e-prints (2025) arXiv:2505.13443 [ 2505.13443]

2025 arXiv
[36]

Stahl, B

C. Stahl, B. Famaey, R. Ibata, O. Hahn, N. Martinet and T. Montandon, Scale-dependent local primordial non-Gaussianity as a solution to the S8 tension, 2404.03244

arXiv
[37]

Adame, S

A.G. Adame, S. Avila, V. Gonzalez-Perez, G. Yepes, M. Pellejero, M.S. Wang et al., PNG-UNITsims: Halo clustering response to primordial non-Gaussianities as a function of mass , Astron. Astrophys. 689 (2024) A69 [ 2312.12405]

2024 arXiv
[38]

Peebles, The Large-Scale Structure of the Universe , Princeton University Press (1980)

P.J.E. Peebles, The Large-Scale Structure of the Universe , Princeton University Press (1980)

1980
[39]

Bernardeau, S

F. Bernardeau, S. Colombi, E. Gazta˜ naga and R. Scoccimarro, Large-scale structure of the Universe and cosmological perturbation theory, Phys. Rep. 367 (2002) 1 [ astro-ph/0112551]

2002 arXiv
[40]

Sirko, Initial Conditions to Cosmological N-Body Simulations, or, How to Run an Ensemble of Simulations, Astrophys

E. Sirko, Initial Conditions to Cosmological N-Body Simulations, or, How to Run an Ensemble of Simulations, Astrophys. J. 634 (2005) 728 [ astro-ph/0503106]

2005 arXiv
[41]

Crocce, S

M. Crocce, S. Pueblas and R. Scoccimarro, Transients from initial conditions in cosmological simulations, Monthly Notices of the Royal Astronomical Society 373 (2006) 369 [ astro-ph/0606505]

2006 arXiv
[42]

Scoccimarro, L

R. Scoccimarro, L. Hui, M. Manera and K.C. Chan, Large-scale bias and efficient generation of initial conditions for nonlocal primordial non-Gaussianity , Phys. Rev. D 85 (2012) 083002 [ 1108.5512]

2012 arXiv
[43]

Astrophys

Euclid collaboration, Euclid: The search for primordial features , Astron. Astrophys. 683 (2024) A220 [2309.17287]

2024 arXiv
[44]

Babich, P

D. Babich, P. Creminelli and M. Zaldarriaga, The shape of non-Gaussianities , Journal of Cosmology and Astroparticle Physics 2004 (2004) 009 [ astro-ph/0405356]

2004 arXiv
[45]

Renaux-Petel, Combined local and equilateral non-Gaussianities from multifield DBI inflation , JCAP 2009 (2009) 012 [ 0907.2476]

S. Renaux-Petel, Combined local and equilateral non-Gaussianities from multifield DBI inflation , JCAP 2009 (2009) 012 [ 0907.2476]

2009 arXiv
[46]

Wagner and L

C. Wagner and L. Verde, N-body simulations with generic non-Gaussian initial conditions II: Halo bias , JCAP 03 (2012) 002 [ 1102.3229]

2012 arXiv
[47]

Wagner, L

C. Wagner, L. Verde and L. Boubekeur, N-body simulations with generic non-Gaussian initial conditions I: Power Spectrum and halo mass function , JCAP 1010 (2010) 022 [ 1006.5793]

2010 arXiv
[48]

Schmidt and M

F. Schmidt and M. Kamionkowski, Halo Clustering with Non-Local Non-Gaussianity , Phys. Rev. D 82 (2010) 103002 [ 1008.0638]

2010 arXiv
[49]

Smith and M

K.M. Smith and M. Zaldarriaga, Algorithms for bispectra: forecasting, optimal analysis and simulation , Monthly Notices of the Royal Astronomical Society 417 (2011) 2 [ astro-ph/0612571]

2011 arXiv
[50]

Holman and A.J

R. Holman and A.J. Tolley, Enhanced non-gaussianity from excited initial states , JCAP 05 (2008) 001 [0710.1302]

2008 arXiv
[51]

Senatore, K.M

L. Senatore, K.M. Smith and M. Zaldarriaga, Non-Gaussianities in Single Field Inflation and their Optimal Limits from the WMAP 5-year Data , JCAP 1001 (2010) 028 [ 0905.3746]

2010 arXiv
[52]

Chen, Folded Resonant Non-Gaussianity in General Single Field Inflation , JCAP 12 (2010) 003 [1008.2485]

X. Chen, Folded Resonant Non-Gaussianity in General Single Field Inflation , JCAP 12 (2010) 003 [1008.2485]. – 18 –

2010 arXiv
[53]

LoVerde, A

M. LoVerde, A. Miller, S. Shandera and L. Verde, Effects of Scale-Dependent Non-Gaussianity on Cosmological Structures, JCAP 0804 (2008) 014 [ 0711.4126]

2008 arXiv
[54]

Sefusatti, M

E. Sefusatti, M. Liguori, A.P.S. Yadav, M.G. Jackson and E. Pajer, Constraining Running Non-Gaussianity, JCAP 12 (2009) 022 [ 0906.0232]

2009 arXiv
[55]

Aghanim, Y

Planck Collaboration, N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi et al., Planck 2018 results. VI. Cosmological parameters , ArXiv e-prints 1807.06209 (2018) [ 1807.06209]

2018 arXiv
[56]

McDonald, Primordial non-Gaussianity: large-scale structure signature in the perturbative bias model, Physical Review D 78 (2008) 123519

P. McDonald, Primordial non-Gaussianity: large-scale structure signature in the perturbative bias model, Physical Review D 78 (2008) 123519

2008
[57]

Lewis, A

A. Lewis, A. Challinor and A. Lasenby, Efficient computation of CMB anisotropies in closed FR W models, Astrophys. J. 538 (2000) 473 [ astro-ph/9911177]

2000 arXiv
[58]

Michaux, O

M. Michaux, O. Hahn, C. Rampf and R.E. Angulo, Accurate initial conditions for cosmological N-body simulations: minimizing truncation and discreteness errors , Monthly Notices of the Royal Astronomical Society 500 (2021) 663 [ 2008.09588]

2021 arXiv
[59]

Adame, S

A.G. Adame, S. Avila, V. Gonzalez-Perez, O. Hahn, G. Yepes and M. Manera, Accurate N -body simulations with local Primordial non-Gaussianities: initial conditions and aliasing , arXiv e-prints (2025) arXiv:2506.06200 [ 2506.06200]

2025
[60]

Orszag, On the Elimination of Aliasing in Finite-Difference Schemes by Filtering High-Wavenumber Components, Journal of the Atmospheric Sciences 28 (1971) 1074

S.A. Orszag, On the Elimination of Aliasing in Finite-Difference Schemes by Filtering High-Wavenumber Components, Journal of the Atmospheric Sciences 28 (1971) 1074

1971
[61]

Bouchet, S

F.R. Bouchet, S. Colombi, E. Hivon and R. Juszkiewicz, Perturbative Lagrangian approach to gravitational instability., Astron. Astrophys. 296 (1995) 575 [ astro-ph/9406013]

1995 arXiv
[62]

Springel, The cosmological simulation code gadget-2 , Monthly Notices of the Royal Astronomical Society 364 (2005) 1105 [ astro-ph/0505010]

V. Springel, The cosmological simulation code gadget-2 , Monthly Notices of the Royal Astronomical Society 364 (2005) 1105 [ astro-ph/0505010]

2005 arXiv
[63]

Potter, J

D. Potter, J. Stadel and R. Teyssier, PKDGRA V3: beyond trillion particle cosmological simulations for the next era of galaxy surveys , Computational Astrophysics and Cosmology 4 (2017) 2 [ 1609.08621]

2017 arXiv
[64]

Springel, R

V. Springel, R. Pakmor, O. Zier and M. Reinecke, Simulating cosmic structure formation with the GADGET-4 code, Monthly Notices of the Royal Astronomical Society 506 (2021) 2871 [ 2010.03567]

2021 arXiv
[65]

Feng, M.-Y

Y. Feng, M.-Y. Chu, U. Seljak and P. McDonald, F ASTPM: a new scheme for fast simulations of dark matter and haloes , Monthly Notices of the Royal Astronomical Society 463 (2016) 2273 [ 1603.00476]

2016 arXiv
[66]

Howlett, M

C. Howlett, M. Manera and W.J. Percival, L-PICOLA: A parallel code for fast dark matter simulation , Astronomy and Computing 12 (2015) 109 [ 1506.03737]

2015 arXiv
[67]

Villaescusa-Navarro, C

F. Villaescusa-Navarro, C. Hahn, E. Massara, A. Banerjee, A.M. Delgado, D.K. Ramanah et al., The Quijote simulations , The Astrophysical Journal Supplement Series 250 (2020) 2. – 19 –

2020

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.

[1] [1]

Maldacena, Non-Gaussian features of primordial fluctuations in single field inflationary models , Journal of High Energy Physics 2003 (2003) 013

J. Maldacena, Non-Gaussian features of primordial fluctuations in single field inflationary models , Journal of High Energy Physics 2003 (2003) 013

work page 2003

[2] [2]

Acquaviva, N

V. Acquaviva, N. Bartolo, S. Matarrese and A. Riotto, Gauge-invariant second-order perturbations and non-Gaussianity from inflation , Nuclear Physics B 667 (2003) 119 [ astro-ph/0209156]

arXiv 2003

[3] [3]

Creminelli and M

P. Creminelli and M. Zaldarriaga, Single field consistency relation for the 3-point function , Journal of Cosmology and Astroparticle Physics 2004 (2004) 006

work page 2004

[4] [4]

Bartolo, E

N. Bartolo, E. Komatsu, S. Matarrese and A. Riotto, Non-Gaussianity from inflation: theory and observations, Phys. Rep. 402 (2004) 103 [ astro-ph/0406398]

arXiv 2004

[5] [5]

Chen, Primordial Non-Gaussianities from Inflation Models , Advances in Astronomy 2010 (2010) 638979 [1002.1416]

X. Chen, Primordial Non-Gaussianities from Inflation Models , Advances in Astronomy 2010 (2010) 638979 [1002.1416]

arXiv 2010

[6] [6]

Meerburg et al., Primordial Non-Gaussianity, 1903.04409

P.D. Meerburg et al., Primordial Non-Gaussianity, 1903.04409. – 16 –

arXiv 1903

[7] [7]

Chen, M.-x

X. Chen, M.-x. Huang, S. Kachru and G. Shiu, Observational signatures and non-Gaussianities of general single-field inflation , JCAP 2007 (2007) 002 [ hep-th/0605045]

arXiv 2007

[8] [8]

Chen and Y

X. Chen and Y. Wang, Large non-Gaussianities with intermediate shapes from quasi-single-field inflation, Phys. Rev. D 81 (2010) 063511 [ 0909.0496]

arXiv 2010