Recognition: unknown
Gravitational wave emission from nonspherical collapse in an early matter-dominated era using N-body simulations
Pith reviewed 2026-05-08 17:36 UTC · model grok-4.3
The pith
Reliable predictions of gravitational waves from early nonspherical collapse require full numerical tracking of nonlinear dynamics.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using N-body simulations of nonspherical collapse initialized through a Zel'dovich deformation of a homogeneous sphere and evolved in an Einstein-de Sitter background, we find that a reliable prediction of the emitted gravitational waves requires a fully numerical treatment of the nonlinear collapse dynamics. In particular, fitting-based procedures and Zel'dovich-based estimates fail to capture the post-shell-crossing evolution and can over- or under-estimate the emitted power of the gravitational waves. After averaging over realizations weighted by the Doroshkevich and BBKS distributions, the two spectra have similar shapes and remain within the same overall order of magnitude at the peak.
What carries the argument
the semirelativistic N-body simulation that numerically evolves the quadrupole moment of the collapsing patch to compute the gravitational wave emission
If this is right
- Averaged spectra from Doroshkevich and BBKS peak distributions have similar shapes and lie within the same order of magnitude in amplitude.
- The dominant contribution arises from peaks of height around three standard deviations.
- Larger variance in the density field significantly enhances the overall signal strength.
- Varying the horizon mass and reheating temperature maps the present-day spectrum across the sensitivity bands of pulsar timing arrays to very high-frequency detectors.
- Gravitational waves from such collapses can serve as a probe of the pre-BBN thermal history.
Where Pith is reading between the lines
- Analytical estimates of gravitational wave backgrounds from early matter-dominated eras may need systematic revision once nonlinear collapse is included.
- The framework could be applied to collapse in other expanding backgrounds or with different initial perturbation statistics.
- A detected stochastic background with the reported spectral shape would constrain the duration of any early matter-dominated phase.
- The similarity between results from two different peak statistics indicates robustness to the precise choice of weighting.
Load-bearing premise
The collapsing patch is initialized through a Zel'dovich deformation of a homogeneous sphere and evolved in an Einstein-de Sitter background.
What would settle it
A calculation that continues the Zel'dovich approximation past shell crossing and obtains a gravitational wave power spectrum that agrees with the full N-body result within the reported uncertainties would falsify the claim that numerical treatment is required.
read the original abstract
We study the dynamics of the collapse of a nonspherical overdense patch during an early matter-dominated era and the associated production of gravitational waves (GWs) using a semirelativistic N-body framework that we develop. The collapsing patch is initialized through a Zel'dovich deformation of a homogeneous sphere and evolved in an Einstein--de Sitter background, while the emitted signal is computed directly from the numerical quadrupole evolution. We show that a reliable prediction of the signal requires a fully numerical treatment of the nonlinear collapse dynamics. In particular, fitting-based procedures and Zel'dovich-based estimates fail to capture the post-shell-crossing evolution and can over/under-estimate the emitted power of the GWs. After averaging over realizations weighted by the Doroshkevich and BBKS (peak theory) distributions, we find that the two spectra have similar shapes and remain within the same overall order of magnitude at the peak amplitude, although the BBKS result is systematically smaller. The dominant contribution arises from peaks of relatively modest height, around $\nu \simeq 3$, while a larger variance significantly enhances the signal. Finally, by varying the horizon mass and reheating temperature, we map the present-day GW spectra to the sensitivity bands of different classes of detectors. In this way, the signal can populate a broad range of frequencies, from pulsar timing arrays to very high-frequency experiments, showing that GWs from nonspherical collapse can provide a probe of the pre-BBN thermal history.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a semirelativistic N-body framework to simulate the collapse of nonspherical overdense patches initialized via Zel'dovich deformation of a homogeneous sphere in an Einstein-de Sitter background during an early matter-dominated era. Gravitational wave emission is extracted directly from the time-dependent quadrupole moment. The central claim is that only this fully numerical treatment reliably captures post-shell-crossing dynamics, while fitting-based procedures and Zel'dovich estimates fail and produce over- or under-estimates of GW power. After averaging realizations weighted by the Doroshkevich and BBKS peak distributions, the resulting spectra have similar shapes and remain within the same order of magnitude at peak amplitude (with BBKS systematically smaller), dominated by modest peaks around ν ≃ 3. Varying horizon mass and reheating temperature maps the present-day spectra across detector frequency bands.
Significance. If the numerical results hold, the work establishes that nonlinear collapse dynamics must be treated fully numerically for accurate GW predictions in matter-dominated eras, providing a benchmark that improves on prior approximations. Strengths include the direct quadrupole evolution, explicit parameter variation to connect to observables, and averaging over realistic peak-theory distributions. This could position GWs from such collapses as a probe of pre-BBN thermal history across pulsar timing arrays to high-frequency detectors.
major comments (2)
- [Description of the semirelativistic N-body framework and GW signal extraction] The claim that numerical results are required because approximations fail post-shell-crossing depends on the semirelativistic quadrupole formula remaining accurate when local densities and velocities become large after shell-crossing. The manuscript provides no quantitative error budget or comparison against fully relativistic simulations to show that neglected relativistic corrections or gauge effects do not alter the radiated power at the level of the reported differences between methods.
- [Averaging and statistical results] In the section on averaging over realizations weighted by the Doroshkevich and BBKS distributions: the exact sampling procedure, number of realizations per distribution, and weighting implementation are not specified in sufficient detail to assess the robustness of the conclusion that the two spectra have similar shapes while the BBKS amplitude is systematically smaller.
minor comments (1)
- [Abstract] The abstract uses 'Einstein--de Sitter' with a double dash; ensure consistent typesetting of the background metric name throughout the manuscript.
Simulated Author's Rebuttal
We thank the referee for the thorough review and valuable feedback on our manuscript. We address the major comments point by point below, indicating where revisions will be made to improve clarity and address concerns.
read point-by-point responses
-
Referee: The claim that numerical results are required because approximations fail post-shell-crossing depends on the semirelativistic quadrupole formula remaining accurate when local densities and velocities become large after shell-crossing. The manuscript provides no quantitative error budget or comparison against fully relativistic simulations to show that neglected relativistic corrections or gauge effects do not alter the radiated power at the level of the reported differences between methods.
Authors: We acknowledge the validity of this concern. Our framework is semirelativistic by design, and while we believe the quadrupole approximation holds because the collapse remains non-relativistic (velocities << c in the early matter-dominated era), we did not provide a quantitative error estimate or direct comparison to full GR simulations. Such a comparison is outside the current scope. In the revised version, we will expand the discussion to include an assessment of the approximation's limitations, referencing studies on the quadrupole formula's accuracy in similar systems, and note that reported differences between methods are likely larger than expected relativistic corrections. revision: partial
-
Referee: In the section on averaging over realizations weighted by the Doroshkevich and BBKS distributions: the exact sampling procedure, number of realizations per distribution, and weighting implementation are not specified in sufficient detail to assess the robustness of the conclusion that the two spectra have similar shapes while the BBKS amplitude is systematically smaller.
Authors: We appreciate this observation and agree that more details are needed for reproducibility. In the revised manuscript, we will add text detailing the sampling procedure, the number of realizations used for each distribution, the method of sampling from the probability density functions, and the exact weighting implementation used to compute the averaged spectra. This will allow proper assessment of the robustness of our conclusions. revision: yes
- Providing a quantitative error budget via comparison to fully relativistic simulations, as this would require new, computationally intensive work beyond the scope of the current study.
Circularity Check
Numerical quadrupole evolution is independent of fitted or self-referential inputs
full rationale
The derivation proceeds from explicit N-body integration of the quadrupole moment in an EdS background, initialized via Zel'dovich deformation but evolved fully numerically to extract GW power. Comparisons to Zel'dovich-based estimates and fitting procedures are performed against external approximations, not by construction. No parameter is fitted to a data subset and then relabeled as a prediction; no uniqueness theorem or ansatz is imported via self-citation to force the central result; the averaging over Doroshkevich/BBKS distributions is a post-processing step using standard peak theory, not a reduction of the emitted signal itself. The chain therefore contains independent dynamical content.
Axiom & Free-Parameter Ledger
free parameters (2)
- horizon mass
- reheating temperature
axioms (3)
- domain assumption Einstein-de Sitter background for the collapsing patch evolution
- domain assumption Zel'dovich deformation of a homogeneous sphere for initial conditions
- standard math Direct computation of GW signal from numerical quadrupole evolution
Forward citations
Cited by 1 Pith paper
-
Reviving primordial black hole formation in slow first-order phase transitions
Primordial black hole formation from slow first-order phase transitions remains viable if slow reheating after supercooling creates an early matter-dominated era allowing small overdensities to grow, producing black h...
Reference graph
Works this paper leans on
-
[1]
R. Allahverdi, M.A. Amin, A. Berlin et al.,The first three seconds: a review of possible expansion histories of the early universe,2006.16182
-
[2]
Carr and S.W
B.J. Carr and S.W. Hawking,Black holes in the early universe,Monthly Notices of the Royal Astronomical Society168(1974) 399 [https://academic.oup.com/mnras/article-pdf/168/2/399/8079885/mnras168-0399.pdf]
1974
- [3]
-
[4]
Primordial Black Holes as a dark matter candidate
A.M. Green and B.J. Kavanagh,Primordial Black Holes as a dark matter candidate,J. Phys. G48(2021) 043001 [2007.10722]
- [5]
-
[6]
Halo Models of Large Scale Structure
A. Cooray and R.K. Sheth,Halo Models of Large Scale Structure,Phys. Rept.372(2002) 1 [astro-ph/0206508]
work page Pith review arXiv 2002
-
[7]
A.L. Erickcek and K. Sigurdson,Reheating effects in the matter power spectrum and implications for substructure,Phys. Rev. D84(2011) 083503 [1106.0536]
-
[8]
Goodman and J
J. Goodman and J. Binney,Gravitational collapse of pressure- less inhomogeneous spheroids,Monthly Notices of the Royal Astronomical Society203(1983) 265 [https://academic.oup.com/mnras/article-pdf/203/2/265/3876955/mnras203-0265.pdf]
1983
-
[9]
Thorne,Multipole expansions of gravitational radiation,Rev
K.S. Thorne,Multipole expansions of gravitational radiation,Rev. Mod. Phys.52(1980) 299
1980
-
[10]
K. Jedamzik, M. Lemoine and J. Martin,Generation of gravitational waves during early structure formation between cosmic inflation and reheating,JCAP04(2010) 021 [1002.3278]
-
[11]
Nakama,Stochastic gravitational waves associated with primordial black holes formed during an early matter era,Phys
T. Nakama,Stochastic gravitational waves associated with primordial black holes formed during an early matter era,Phys. Rev. D101(2020) 063519
2020
-
[12]
B. Eggemeier, J.C. Niemeyer, K. Jedamzik and R. Easther,Stochastic gravitational waves from postinflationary structure formation,Phys. Rev. D107(2023) 043503 [2212.00425]
- [13]
-
[14]
B. Eggemeier, P. Hayman, J.C. Niemeyer and R. Easther,Postinflationary structure formation boosted by parametric self-resonance,Phys. Rev. D109(2024) 043521 [2311.08780]
-
[16]
N. Fernandez, J.W. Foster, B. Lillard and J. Shelton,Stochastic gravitational waves from early structure formation,Phys. Rev. Lett.133(2024) 111002 [2312.12499]
-
[17]
X.-X. Zeng, Z. Ning, Z.-Y. Yuwen, S.-J. Wang, H. Deng and R.-G. Cai,Relic gravitational waves from primordial gravitational collapses,2504.11275
work page internal anchor Pith review Pith/arXiv arXiv
-
[18]
I. Dalianis and C. Kouvaris,Gravitational waves from density perturbations in an early matter domination era,JCAP07(2021) 046 [2012.09255]
-
[19]
I. Dalianis and C. Kouvaris,Gravitational waves from collapse of pressureless matter in the early universe,JCAP10(2024) 006 [2403.15126]
-
[20]
H. Assadullahi and D. Wands,Gravitational waves from an early matter era,Phys. Rev. D79 (2009) 083511 [0901.0989]
-
[21]
K. Kohri and T. Terada,Semianalytic calculation of gravitational wave spectrum nonlinearly induced from primordial curvature perturbations,Phys. Rev. D97(2018) 123532 [1804.08577]
-
[22]
K. Inomata, K. Kohri, T. Nakama and T. Terada,Gravitational Waves Induced by Scalar Perturbations during a Gradual Transition from an Early Matter Era to the Radiation Era, JCAP10(2019) 071 [1904.12878]
-
[23]
K. Inomata, K. Kohri, T. Nakama and T. Terada,Enhancement of Gravitational Waves Induced by Scalar Perturbations due to a Sudden Transition from an Early Matter Era to the Radiation Era,Phys. Rev. D100(2019) 043532 [1904.12879]
- [24]
- [25]
-
[26]
The poltergeist mechanism – Enhancement of scalar-induced gravitational waves with early matter-dominated era
K. Inomata, K. Kohri and T. Terada, “The poltergeist mechanism – Enhancement of scalar-induced gravitational waves with early matter-dominated era.” 11, 2025
2025
-
[27]
C. Rampf,Cosmological vlasov–poisson equations for dark matter: Recent developments and connections to selected plasma problems,Rev. Mod. Plasma Phys.5(2021) 10 [2110.06265]
-
[28]
Rampf et al.,Shell-crossing in quasi-one-dimensional flow,Mon
C. Rampf et al.,Shell-crossing in quasi-one-dimensional flow,Mon. Not. Roy. Astron. Soc. 471(2017) 671
2017
-
[29]
Adamek et al.,Distance-redshift relation in plane symmetric universes,Phys
J. Adamek et al.,Distance-redshift relation in plane symmetric universes,Phys. Rev. D89 (2014) 063543
2014
-
[30]
S. Pueblas and R. Scoccimarro,Generation of vorticity and velocity dispersion by orbit crossing,Phys. Rev. D80(2009) 043504 [0809.4606]
-
[31]
Doroshkevich,Spatial structure of perturbations and origin of galactic rotation in fluctuation theory,Astrophysics6(1970) 320
A.G. Doroshkevich,Spatial structure of perturbations and origin of galactic rotation in fluctuation theory,Astrophysics6(1970) 320
1970
-
[32]
Bardeen, J.R
J.M. Bardeen, J.R. Bond, N. Kaiser and A.S. Szalay,The statistics of peaks of gaussian random fields,Astrophys. J.304(1986) 15
1986
-
[33]
Hockney and J.W
R.W. Hockney and J.W. Eastwood,Computer Simulation Using Particles, McGraw-Hill, New York (1981)
1981
-
[34]
V. Springel, N. Yoshida and S.D.M. White,GADGET: A Code for collisionless and gasdynamical cosmological simulations,New Astron.6(2001) 79 [astro-ph/0003162]. – 54 –
-
[35]
D.J. Price and J.J. Monaghan,An energy-conserving formalism for adaptive gravitational force softening in SPH and N-body codes,Mon. Not. Roy. Astron. Soc.374(2007) 1347 [astro-ph/0610872]
-
[36]
Monaghan and J.C
J.J. Monaghan and J.C. Lattanzio,A refined particle method for astrophysical problems, A&A 149(1985) 135
1985
-
[37]
Lebovitz,The virial tensor and its application to self-gravitating fluids,Astrophysical Journal134(1961) 500
N.R. Lebovitz,The virial tensor and its application to self-gravitating fluids,Astrophysical Journal134(1961) 500
1961
-
[38]
Layzer,A preface to cosmogony
D. Layzer,A preface to cosmogony. i. the energy equation and the virial theorem for cosmic distributions,Astrophysical Journal138(1963) 174
1963
-
[39]
Zel’dovich,An approximate theory for large density perturbations,Astron
Y.B. Zel’dovich,An approximate theory for large density perturbations,Astron. Astrophys.5 (1970) 84
1970
-
[40]
Peters and J
P.C. Peters and J. Mathews,Gravitational radiation from point masses in a keplerian orbit, Phys. Rev.131(1963) 435
1963
-
[41]
K. Saikawa and S. Shirai,Primordial gravitational waves, precisely: The role of thermodynamics in the Standard Model,JCAP05(2018) 035 [1803.01038]
- [42]
-
[43]
Albert escriva github repository
A. Escriva, “Albert escriva github repository.”https://github.com/albert-escriva, 2026. Accessed: 2026-04-20
2026
-
[44]
N. Yoshida, N. Sugiyama and L. Hernquist,The evolution of baryon density fluctuations in multicomponent cosmological simulations,Monthly Notices of the Royal Astronomical Society 344(2003) 481 [astro-ph/0305210]
-
[45]
C.C. Lin, L. Mestel and F.H. Shu,The Gravitational Collapse of a Uniform Spheroid., ApJ 142(1965) 1431
1965
-
[46]
A. Escriv` a and C.-M. Yoo,Simulations of ellipsoidal primordial black hole formation,Phys. Rev. D112(2025) 083518 [2410.03452]
-
[47]
Schmitz, JHEP01, 097 (2021), arXiv:2002.04615 [hep-ph]
K. Schmitz,New Sensitivity Curves for Gravitational-Wave Signals from Cosmological Phase Transitions,JHEP01(2021) 097 [2002.04615]
- [48]
-
[49]
Khlopov and A.G
M.Y. Khlopov and A.G. Polnarev,PRIMORDIAL BLACK HOLES AS A COSMOLOGICAL TEST OF GRAND UNIFICATION,Phys. Lett. B97(1980) 383
1980
-
[50]
Primordial black hole formation in the matter-dominated phase of the Universe,
T. Harada, C.-M. Yoo, K. Kohri, K.-i. Nakao and S. Jhingan,Primordial black hole formation in the matter-dominated phase of the Universe,Astrophys. J.833(2016) 61 [1609.01588]
-
[51]
Spins of primordial black holes formed in the matter-dominated phase of the Universe,
T. Harada, C.-M. Yoo, K. Kohri and K.-I. Nakao,Spins of primordial black holes formed in the matter-dominated phase of the Universe,Phys. Rev. D96(2017) 083517 [1707.03595]
- [52]
- [53]
-
[54]
I. Dalianis,Constraints on the curvature power spectrum from primordial black hole evaporation,JCAP08(2019) 032 [1812.09807]
-
[55]
M. Kawasaki, K. Kohri and N. Sugiyama,Cosmological constraints on late time entropy production,Phys. Rev. Lett.82(1999) 4168 [astro-ph/9811437]. – 55 –
-
[56]
M. Kawasaki, K. Kohri and N. Sugiyama,MeV scale reheating temperature and thermalization of neutrino background,Phys. Rev. D62(2000) 023506 [astro-ph/0002127]
-
[57]
K. Ichikawa, M. Kawasaki and F. Takahashi,The Oscillation effects on thermalization of the neutrinos in the Universe with low reheating temperature,Phys. Rev. D72(2005) 043522 [astro-ph/0505395]
-
[58]
P.F. de Salas, M. Lattanzi, G. Mangano, G. Miele, S. Pastor and O. Pisanti,Bounds on very low reheating scenarios after Planck,Phys. Rev. D92(2015) 123534 [1511.00672]
-
[59]
T. Hasegawa, N. Hiroshima, K. Kohri, R.S.L. Hansen, T. Tram and S. Hannestad,MeV-scale reheating temperature and thermalization of oscillating neutrinos by radiative and hadronic decays of massive particles,JCAP12(2019) 012 [1908.10189]
-
[60]
N. Barbieri, T. Brinckmann, S. Gariazzo, M. Lattanzi, S. Pastor and O. Pisanti,Current Constraints on Cosmological Scenarios with Very Low Reheating Temperatures,Phys. Rev. Lett.135(2025) 181003 [2501.01369]. [61]NANOGravcollaboration,The NANOGrav 15 yr Data Set: Evidence for a Gravitational-wave Background,Astrophys. J. Lett.951(2023) L8 [2306.16213]. [6...
-
[61]
Search for an isotropic gravitational-wave background with the Parkes Pulsar Timing Array
D.J. Reardon et al.,Search for an Isotropic Gravitational-wave Background with the Parkes Pulsar Timing Array,Astrophys. J. Lett.951(2023) L6 [2306.16215]
work page internal anchor Pith review arXiv 2023
-
[62]
Janssenet al., PoSAASKA14, 037 (2015), arXiv:1501.00127 [astro-ph.IM]
G. Janssen et al.,Gravitational wave astronomy with the SKA,PoSAASKA14(2015) 037 [1501.00127]. [66]EPTA, InPTA:collaboration,The second data release from the European Pulsar Timing Array - III. Search for gravitational wave signals,Astron. Astrophys.678(2023) A50 [2306.16214]
-
[63]
H. Xu et al.,Searching for the Nano-Hertz Stochastic Gravitational Wave Background with the Chinese Pulsar Timing Array Data Release I,Res. Astron. Astrophys.23(2023) 075024 [2306.16216]
work page internal anchor Pith review arXiv 2023
- [64]
-
[65]
J. Crowder and N.J. Cornish,Beyond LISA: Exploring future gravitational wave missions, Phys. Rev. D72(2005) 083005 [gr-qc/0506015]
-
[66]
Laser Interferometer Space Antenna
P. Amaro-Seoane, H. Audley, S. Babak, J. Baker, E. Barausse, P. Bender et al.,Laser Interferometer Space Antenna,arXiv e-prints(2017) arXiv:1702.00786 [1702.00786]
work page internal anchor Pith review arXiv 2017
- [67]
-
[68]
Z. Luo, Y. Wang, Y. Wu, W. Hu and G. Jin,The Taiji program: A concise overview,PTEP 2021(2021) 05A108
2021
-
[69]
Z.-C. Liang, Y.-M. Hu, Y. Jiang, J. Cheng, J.-d. Zhang and J. Mei,Science with the TianQin Observatory: Preliminary results on stochastic gravitational-wave background,Phys. Rev. D 105(2022) 022001 [2107.08643]
-
[70]
Sesana et al.,Unveiling the gravitational universe at µ-Hz frequencies,Exper
A. Sesana et al.,Unveiling the gravitational universe atµ-Hz frequencies,Exper. Astron.51 (2021) 1333 [1908.11391]. – 56 – [75]KAGRA, Virgo, LIGO Scientificcollaboration,Upper limits on the isotropic gravitational-wave background from Advanced LIGO and Advanced Virgo’s third observing run, Phys. Rev. D104(2021) 022004 [2101.12130]. [76]LIGO Scientific, VI...
-
[71]
Barsotti, L
L. Barsotti, L. McCuller, M. Evans and P. Fritschel,The A+ design curve, Tech. Rep. LIGO-T1800042-v5, LIGO Laboratory (Mar., 2018)
2018
-
[72]
Sathyaprakash et al.,Scientific Objectives of Einstein Telescope,Class
B. Sathyaprakash et al.,Scientific Objectives of Einstein Telescope,Class. Quant. Grav.29 (2012) 124013 [1206.0331]
-
[73]
Cosmic Explorer: The U.S. Contribution to Gravitational-Wave Astronomy beyond LIGO
D. Reitze et al.,Cosmic Explorer: The U.S. Contribution to Gravitational-Wave Astronomy beyond LIGO,Bull. Am. Astron. Soc.51(2019) 035 [1907.04833]
work page internal anchor Pith review arXiv 2019
- [74]
- [75]
-
[76]
Y. Jiang and T. Suyama,Spectrum of high-frequency gravitational waves from graviton bremsstrahlung by the decay of inflaton: case with polynomial potential,JCAP02(2025) 041 [2410.11175]. – 57 –
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.