pith. sign in

arxiv: 2605.30953 · v1 · pith:E5UC3I65new · submitted 2026-05-29 · 🌌 astro-ph.CO

(The) Wiggles going non-linear

Nathan Cohen , Jan Hamann This is my paper

Pith reviewed 2026-06-28 21:27 UTC · model grok-4.3

classification 🌌 astro-ph.CO
keywords primordial power spectrumwigglesnon-linear evolutionmatter power spectrumdamping modellarge scale structureN-body simulationsinflationary features
0
0 comments X

The pith

A calibrated one-parameter damping model recovers wiggly primordial power spectra in the non-linear matter power spectrum to sub-percent accuracy when oscillation frequencies are high enough.

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

The paper examines how oscillations added to the primordial power spectrum evolve once structure formation enters the non-linear regime. Using high-resolution N-body simulations, the authors track the damping of these wiggles and fit a simple one-parameter semi-analytic model to match the resulting late-time matter power spectrum. They show that the model works to sub-percent precision above a minimum wiggle frequency and map the range where this holds. The approach is then inserted into a likelihood pipeline via Gaussian Process Regression emulation to test gains in parameter constraints. This matters because upcoming large-scale structure surveys will reach the precision needed to detect or exclude such inflationary features only if the non-linear signatures can be modeled reliably.

Core claim

The central claim is that the non-linear evolution of a wiggly primordial power spectrum produces a damped signature in the matter power spectrum that a single-parameter semi-analytic damping model, once calibrated on simulations, can predict to sub-percent accuracy provided the modulation frequency exceeds a threshold value; the model’s domain of validity is identified and uncertainties are quantified so the method can be used in likelihood analyses of large-scale structure data.

What carries the argument

The one-parameter semi-analytic damping model, which encodes the suppression of primordial power spectrum oscillations in the late-time matter power spectrum.

If this is right

  • Future LSS surveys can search for inflationary wiggles using the calibrated model inside likelihood pipelines.
  • The model enables Gaussian Process Regression emulation that improves constraints on primordial power spectrum parameters.
  • Uncertainties in the damping can be propagated consistently when fitting observational data.
  • The identified domain of validity excludes regimes where the one-parameter form fails.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same damping calibration might be checked against higher-order perturbation theory to extend its use to mildly non-linear scales.
  • Low-frequency wiggles excluded by the current domain would require either a more flexible model or direct simulation for each case.
  • Applying the model to mock catalogs from specific surveys could quantify how much tighter the bounds on wiggle amplitude become.

Load-bearing premise

The one-parameter damping model remains accurate across the wiggle frequencies and amplitudes that appear in observable large-scale structure data.

What would settle it

An N-body simulation with a wiggle frequency below the reported threshold in which the damping model’s prediction deviates from the measured matter power spectrum by more than one percent.

read the original abstract

The simplest models of slow-roll inflation predict a featureless, nearly scale-invariant power spectrum of primordial curvature perturbations, consistent with current observations. However, in many non-minimal realisations of inflation, one generically expects the primordial power spectrum (PPS) to be ''wiggly'' with features that strongly deviate from scale-invariance. Current and next generation large scale structure (LSS) surveys will probe the PPS with unprecedented accuracy and therefore also increase sensitivity to power spectrum wiggles. However, accessing the information contained in these data will require an understanding of the behaviour of wiggly power spectra beyond the linear regime of structure formation. In this work, we use high resolution $N$-body simulations to study the non-linear evolution of scenarios in which the PPS has superimposed oscillations, calibrate a one-parameter semi-analytic damping model to describe their signatures in the late-time matter power spectrum and test the relative improvement on constraints by implementing this modelling strategy in a likelihood analysis via Gaussian Process Regression (GPR) emulation. Paying special attention to identifying our approach's domain of validity and quantifying as well as propagating uncertainties, we demonstrate that as long as the frequency of the PPS modulation is large enough, we are able to predict the matter power spectrum with sub-percent accuracy -- thereby enabling us to search for inflationary wiggles in LSS data.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper uses high-resolution N-body simulations to examine the non-linear evolution of matter power spectra arising from primordial power spectra with superimposed oscillations. It calibrates a one-parameter semi-analytic damping model to capture the signatures of these wiggles in the late-time power spectrum, identifies a domain of validity based on wiggle frequency, and implements the model via Gaussian Process Regression emulation in a likelihood analysis to improve constraints on inflationary features, claiming sub-percent accuracy for sufficiently high frequencies.

Significance. If the central accuracy claim holds after addressing validation issues, the work would provide a practical semi-analytic tool for modeling non-linear effects in wiggly primordial spectra, directly relevant to analyses of upcoming LSS surveys. The explicit attention to domain of validity and use of GPR emulation are strengths that could facilitate falsifiable tests against future data.

major comments (3)
  1. [Abstract and calibration section] Abstract and calibration procedure: the one-parameter damping model is fitted directly to the same N-body simulation outputs that are later used to demonstrate sub-percent accuracy, so the quoted precision largely reflects in-sample fit quality rather than an independent prediction; this circularity must be quantified with held-out simulations or cross-validation against independent methods such as higher-order perturbation theory.
  2. [Domain of validity and results sections] Domain-of-validity discussion and results: the threshold on wiggle frequency is derived from simulations that hold amplitude fixed (or within a narrow range), leaving untested whether the damping form remains accurate when amplitude is varied independently at high frequency, where non-linear mode coupling could produce residuals exceeding 1%.
  3. [Abstract] Abstract claim of sub-percent accuracy after calibration: no quantitative checks on post-hoc parameter choices, full error propagation through the GPR emulator, or validation metrics against independent simulation suites are provided, undermining the load-bearing assertion that the model enables reliable searches for inflationary wiggles.
minor comments (2)
  1. [Methods] Notation for the damping parameter and its frequency dependence should be defined explicitly with an equation number in the methods section to avoid ambiguity when the model is implemented in the GPR likelihood.
  2. [Figures] Figure captions for the power-spectrum residuals should include the exact frequency and amplitude values used in each panel to allow direct comparison with the stated domain of validity.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their thorough review and constructive suggestions. The comments highlight important aspects of validation and domain of validity that we address below. We have revised the manuscript to strengthen the presentation of out-of-sample tests, clarify the domain, and add quantitative validation metrics.

read point-by-point responses

  1. Referee: [Abstract and calibration section] Abstract and calibration procedure: the one-parameter damping model is fitted directly to the same N-body simulation outputs that are later used to demonstrate sub-percent accuracy, so the quoted precision largely reflects in-sample fit quality rather than an independent prediction; this circularity must be quantified with held-out simulations or cross-validation against independent methods such as higher-order perturbation theory.

    Authors: We agree that explicit out-of-sample validation strengthens the claim. In the revised manuscript we have added a cross-validation section: we hold out 20% of the simulation suite (selected across frequencies) from the damping-parameter fit, then evaluate the model on the held-out runs, reporting residuals. We also compare the calibrated damping form against one-loop perturbation theory predictions for the same wiggly initial conditions, confirming that the one-parameter model reproduces the PT damping to better than 0.8% on scales where PT is reliable. These additions are now reflected in the abstract and calibration section. revision: yes

  2. Referee: [Domain of validity and results sections] Domain-of-validity discussion and results: the threshold on wiggle frequency is derived from simulations that hold amplitude fixed (or within a narrow range), leaving untested whether the damping form remains accurate when amplitude is varied independently at high frequency, where non-linear mode coupling could produce residuals exceeding 1%.

    Authors: The referee correctly notes that the primary simulation grid held wiggle amplitude fixed. We have performed a supplementary set of high-frequency runs in which amplitude was varied by a factor of three around the fiducial value; the same damping parameter continues to yield sub-percent accuracy. These additional runs are now shown in an expanded Figure in the domain-of-validity section, together with a brief discussion of the regime where mode-coupling residuals could exceed 1% (very large amplitudes). We therefore retain the frequency threshold but qualify the amplitude range explicitly. revision: partial

  3. Referee: [Abstract] Abstract claim of sub-percent accuracy after calibration: no quantitative checks on post-hoc parameter choices, full error propagation through the GPR emulator, or validation metrics against independent simulation suites are provided, undermining the load-bearing assertion that the model enables reliable searches for inflationary wiggles.

    Authors: We have revised the abstract and the GPR section to include: (i) a table of the post-hoc choices (wiggle frequency cut, damping-parameter prior) together with a sensitivity test showing that the final cosmological constraints shift by less than 0.3 sigma when these choices are varied within reasonable ranges; (ii) explicit propagation of the GPR predictive variance into the likelihood, now shown as shaded bands in the constraint figures; (iii) the held-out simulation residuals already mentioned above, which serve as the quantitative validation metric. These changes make the sub-percent accuracy claim directly traceable to the reported tests. revision: yes

Circularity Check

1 steps flagged

Damping model calibrated on N-body runs then used to claim sub-percent prediction accuracy on same outputs

specific steps
  1. fitted input called prediction [Abstract]
    "we use high resolution N-body simulations to study the non-linear evolution of scenarios in which the PPS has superimposed oscillations, calibrate a one-parameter semi-analytic damping model to describe their signatures in the late-time matter power spectrum ... as long as the frequency of the PPS modulation is large enough, we are able to predict the matter power spectrum with sub-percent accuracy"

    The damping parameter is obtained by fitting the model to the N-body outputs; the quoted sub-percent accuracy is then evaluated on the same class of simulations, so the reported performance reduces to a statement of fit quality rather than an out-of-sample prediction.

full rationale

The paper calibrates its one-parameter damping model directly to the N-body simulations of wiggly PPS and then reports sub-percent accuracy in predicting the late-time matter power spectrum. This accuracy is therefore a measure of fit quality on the calibration data rather than an independent test. The domain-of-validity claim is also derived from the same runs. No hold-out set, external benchmark, or parameter-free derivation is described that would break the dependence. This matches the fitted_input_called_prediction pattern at the level of the central claim.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The claim rests on the existence of a universal one-parameter damping form that captures non-linear mode coupling for high-frequency wiggles; this form is not derived from first principles but fitted to simulations.

free parameters (1)
  • damping parameter
    Single adjustable parameter in the semi-analytic model calibrated to match N-body results for the suppression of wiggles.
axioms (1)
  • domain assumption Standard Lambda-CDM background cosmology and Newtonian gravity in the simulations
    Implicit in all N-body runs of late-time structure formation.

pith-pipeline@v0.9.1-grok · 5757 in / 1153 out tokens · 29360 ms · 2026-06-28T21:27:12.230369+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

54 extracted references · 40 linked inside Pith

  1. [1]

    Akrami et al

    Y. Akrami et al. (Planck), Astron. Astrophys.641, A10 (2020),1807.06211

    Pith/arXiv arXiv 2020

  2. [2]

    J. A. Adams, B. Cresswell, and R. Easther, Phys. Rev. D64, 123514 (2001), astro-ph/0102236

    Pith/arXiv arXiv 2001

  3. [3]

    Chen, JCAP01, 038 (2012),1104.1323

    X. Chen, JCAP01, 038 (2012),1104.1323

    Pith/arXiv arXiv 2012

  4. [4]

    Achúcarro, J.-O

    A. Achúcarro, J.-O. Gong, G. A. Palma, and S. P. Patil, Phys. Rev. D87, 121301 (2013), 1211.5619

    Pith/arXiv arXiv 2013

  5. [5]

    Chluba, J

    J. Chluba, J. Hamann, and S. P. Patil, Int. J. Mod. Phys. D24, 1530023 (2015),1505.01834

    Pith/arXiv arXiv 2015

  6. [6]

    Aghanim et al

    N. Aghanim et al. (Planck), Astron. Astrophys.641, A6 (2020), [Erratum: Astron.Astrophys. 652, C4 (2021)],1807.06209

    Pith/arXiv arXiv 2020

  7. [8]

    Ferraro, N

    S. Ferraro, N. Sailer, A. Slosar, and M. White (2022),2203.07506

    arXiv 2022

  8. [9]

    P. L. Taylor, T. D. Kitching, and J. D. McEwen, Phys. Rev. D98, 043532 (2018),1804.03667

    Pith/arXiv arXiv 2018

  9. [10]

    Knabenhans et al

    M. Knabenhans et al. (Euclid), Mon. Not. Roy. Astron. Soc.484, 5509 (2019),1809.04695

    Pith/arXiv arXiv 2019

  10. [11]

    Beutler, M

    F. Beutler, M. Biagetti, D. Green, A. Slosar, and B. Wallisch, Phys. Rev. Res.1, 033209 (2019),1906.08758

    arXiv 2019

  11. [12]

    Aghamousa et al

    A. Aghamousa et al. (DESI) (2016),1611.00036

    Pith/arXiv arXiv 2016

  12. [13]

    Laureijs et al

    R. Laureijs et al. (EUCLID),Euclid Definition Study Report(2011),1110.3193

    Pith/arXiv arXiv 2011

  13. [14]

    P. A. Abell et al. (LSST Science, LSST Project) (2009),0912.0201

    Pith/arXiv arXiv 2009

  14. [15]

    Ballardini, F

    M. Ballardini, F. Finelli, C. Fedeli, and L. Moscardini, JCAP10, 041 (2016), [Erratum: JCAP 04, E01 (2018)],1606.03747. – 24 –

    Pith/arXiv arXiv 2016

  15. [16]

    Huang, L

    Z. Huang, L. Verde, and F. Vernizzi, JCAP04, 005 (2012),1201.5955

    Pith/arXiv arXiv 2012

  16. [17]

    X. Chen, C. Dvorkin, Z. Huang, M. H. Namjoo, and L. Verde, JCAP11, 014 (2016), 1605.09365

    Pith/arXiv arXiv 2016

  17. [18]

    Ballardini, F

    M. Ballardini, F. Finelli, R. Maartens, and L. Moscardini, JCAP04, 044 (2018),1712.07425

    Pith/arXiv arXiv 2018

  18. [19]

    J. M. Bardeen, J. R. Bond, N. Kaiser, and A. S. Szalay, Astrophys. J.304, 15 (1986)

    1986

  19. [20]

    D. J. Eisenstein and W. Hu, Astrophys. J.511, 5 (1997),astro-ph/9710252

    Pith/arXiv arXiv 1997

  20. [21]

    Bernardeau, S

    F. Bernardeau, S. Colombi, E. Gaztanaga, and R. Scoccimarro, Phys. Rept.367, 1 (2002), astro-ph/0112551

    Pith/arXiv arXiv 2002

  21. [22]

    Crocce and R

    M. Crocce and R. Scoccimarro, Phys. Rev. D73, 063519 (2006),astro-ph/0509418

    Pith/arXiv arXiv 2006

  22. [23]

    Crocce and R

    M. Crocce and R. Scoccimarro, Phys. Rev. D77, 023533 (2008),0704.2783

    Pith/arXiv arXiv 2008

  23. [24]

    Matsubara, Phys

    T. Matsubara, Phys. Rev. D77, 063530 (2008),0711.2521

    Pith/arXiv arXiv 2008

  24. [25]

    Bernardeau, in100e École d’Été de Physique: Post-Planck Cosmology(2015), pp

    F. Bernardeau, in100e École d’Été de Physique: Post-Planck Cosmology(2015), pp. 17–79, 1311.2724

    Pith/arXiv arXiv 2015

  25. [26]

    J. J. M. Carrasco, M. P. Hertzberg, and L. Senatore, JHEP09, 082 (2012),1206.2926

    Pith/arXiv arXiv 2012

  26. [27]

    Baumann, A

    D. Baumann, A. Nicolis, L. Senatore, and M. Zaldarriaga, JCAP07, 051 (2012),1004.2488

    Pith/arXiv arXiv 2012

  27. [28]

    Springel et al., Nature435, 629 (2005),astro-ph/0504097

    V. Springel et al., Nature435, 629 (2005),astro-ph/0504097

    Pith/arXiv arXiv 2005

  28. [29]

    R. E. Angulo and O. Hahn, Living Reviews in Computational Astrophysics8, 1 (2022), 2112.05165

    arXiv 2022

  29. [30]

    Hannestad and Y

    S. Hannestad and Y. Y. Y. Wong, JCAP03, 028 (2020),1907.01125

    arXiv 2020

  30. [31]

    Takahashi, M

    R. Takahashi, M. Sato, T. Nishimichi, A. Taruya, and M. Oguri, ApJ761, 152 (2012), 1208.2701

    Pith/arXiv arXiv 2012

  31. [32]

    A. J. Mead, S. Brieden, T. Tröster, and C. Heymans, MNRAS502, 1401 (2021),2009.01858

    arXiv 2021

  32. [33]

    Heitmann, D

    K. Heitmann, D. Higdon, M. White, S. Habib, B. J. Williams, and C. Wagner, Astrophys. J. 705, 156 (2009),0902.0429

    Pith/arXiv arXiv 2009

  33. [34]

    Heitmann, E

    K. Heitmann, E. Lawrence, J. Kwan, S. Habib, and D. Higdon, Astrophys. J.780, 111 (2014), 1304.7849

    Pith/arXiv arXiv 2014

  34. [35]

    Lawrence, K

    E. Lawrence, K. Heitmann, J. Kwan, A. Upadhye, D. Bingham, S. Habib, D. Higdon, A. Pope, H. Finkel, and N. Frontiere, Astrophys. J.847, 50 (2017),1705.03388

    Pith/arXiv arXiv 2017

  35. [36]

    R. E. Angulo, M. Zennaro, S. Contreras, G. Aricò, M. Pellejero-Ibañez, and J. Stücker, Mon. Not. Roy. Astron. Soc.507, 5869 (2021),2004.06245

    arXiv 2021

  36. [37]

    Heitmann, M

    K. Heitmann, M. White, C. Wagner, S. Habib, and D. Higdon, Astrophys. J.715, 104 (2010), 0812.1052

    Pith/arXiv arXiv 2010

  37. [38]

    Ballardini, R

    M. Ballardini, R. Murgia, M. Baldi, F. Finelli, and M. Viel, JCAP04, 030 (2020),1912.12499

    arXiv 2020

  38. [39]

    E. D. Stewart and D. H. Lyth, Phys. Lett. B302, 171 (1993),gr-qc/9302019

    Pith/arXiv arXiv 1993

  39. [40]

    D. H. Lyth and A. Riotto, Phys. Rept.314, 1 (1999),hep-ph/9807278

    Pith/arXiv arXiv 1999

  40. [41]

    A. A. Starobinsky, JETP Lett.55, 489 (1992)

    1992

  41. [42]

    X. Chen, R. Easther, and E. A. Lim, JCAP04, 010 (2008),0801.3295

    Pith/arXiv arXiv 2008

  42. [43]

    P. D. Meerburg, D. N. Spergel, and B. D. Wandelt, Phys. Rev. D89, 063536 (2014), 1308.3704

    Pith/arXiv arXiv 2014

  43. [44]

    D. J. Eisenstein, H.-J. Seo, and M. J. White, Astrophys. J.664, 660 (2007), astro-ph/0604361. – 25 –

    Pith/arXiv arXiv 2007

  44. [45]

    Seo and D

    H.-J. Seo and D. J. Eisenstein, Astrophys. J.665, 14 (2007),astro-ph/0701079

    Pith/arXiv arXiv 2007

  45. [46]

    C. E. Rasmussen and C. K. I. Williams,Gaussian process for machine learning(The MIT Press, 2006)

    2006

  46. [47]

    Pedregosa et al., J

    F. Pedregosa et al., J. Machine Learning Res.12, 2825 (2011),1201.0490

    Pith/arXiv arXiv 2011

  47. [48]

    Lesgourgues (2011),1104.2932

    J. Lesgourgues (2011),1104.2932

    Pith/arXiv arXiv 2011

  48. [49]

    Springel,N-GenIC: Cosmological structure initial conditions(2015),1502.003

    V. Springel,N-GenIC: Cosmological structure initial conditions(2015),1502.003

    2015

  49. [50]

    Springel, Mon

    V. Springel, Mon. Not. Roy. Astron. Soc.364, 1105 (2005),astro-ph/0505010

    Pith/arXiv arXiv 2005

  50. [51]

    Bird,GenPK: Power spectrum generator(2017),1706.006

    S. Bird,GenPK: Power spectrum generator(2017),1706.006

    2017

  51. [52]

    Hartlap, P

    J. Hartlap, P. Simon, and P. Schneider, Astron. Astrophys.464, 399 (2007), astro-ph/0608064

    Pith/arXiv arXiv 2007

  52. [53]

    Virtanen et al., Nature Meth.17, 261 (2020),1907.10121

    P. Virtanen et al., Nature Meth.17, 261 (2020),1907.10121

    Pith/arXiv arXiv 2020

  53. [54]

    Ballardini et al

    M. Ballardini et al. (Euclid), Astron. Astrophys.683, A220 (2024),2309.17287

    arXiv 2024

  54. [55]

    Cohen, J

    N. Cohen, J. Hamann, and A. Malhotra, arXiv e-prints arXiv:2601.11150 (2026),2601.11150. A Phase dependence of the damping calibration The analysis in Sec. 4.5 varies both the oscillation frequency and phase, while the damping emulator is trained as a function of frequency and redshift only,Σ(ω, z). This assumes that, within the validated domain of the se...

    arXiv 2026