Multi-channel hunt for the Big Bang's echo: lenses, gas, and X-rays versus primordial non-Gaussianity ✛ ✦
Context
Primordial non-Gaussianity is a deviation of the initial density field's statistics from a Gaussian distribution. It arises in various inflation scenarios and is sensitive to particle physics at energies around $10^{15}$ GeV, making it a unique window into fundamental theory. Traditional methods using cosmic microwave background and galaxy clustering have already achieved high precision, but they make little use of nonlinear structures. Meanwhile, massive galaxy clusters are especially sensitive to the non-Gaussian tails of the density distribution. Thanks to the expansion of the universe, described by Hubble's law, new sky surveys like the Legacy Survey of Space and Time (LSST) at the Vera Rubin Observatory, the Simons Observatory, and eROSITA will for the first time study these structures using three observables at once: weak gravitational lensing, thermal Sunyaev-Zeldovich effect, and X-ray emission.
Methods
The authors built synthetic sky maps using N-body simulations from the Ulagam suite, where the initial conditions were modified to reproduce six different models of primordial non-Gaussianity (local, equilateral, quasi-single field, and others). Using numerical simulation and the semi-analytic BaryonForge model, they accounted for the effects of dark matter and gas on the observables: weak lensing was modeled via the density integral along the line of sight, clusters produced a microwave signal through photon scattering off hot electrons, and the X-ray flux was calculated from thermal emission of the gas taking metallicity into account. Realistic noise components and correlated background sources, including active galactic nuclei, were also added. To compress the information, second- and third-order moments of the smoothed fields were computed, and then the Fisher method was used to estimate the potential accuracy of future measurements.
Results
The analysis showed that the weak lensing data from LSST (14,000 sq. degrees) alone can provide constraints on the non-Gaussianity amplitude $f_{NL}$ comparable to the best results from the Planck satellite, and for some models—twice as good. However, the need to simultaneously fit cosmology and baryonic physics parameters reduces accuracy by about 2–3 times. Maps of the Sunyaev-Zeldovich effect and X-ray emission individually give weaker constraints, but using them jointly with lensing breaks degeneracies between parameters. As a result, combining all three fields provides a 2-fold improvement in precision compared to lensing alone, even after marginalizing over astrophysical uncertainties. A particularly strong improvement comes from cross-calibrating the gas ejection radius around clusters. Important contributions came from pioneers: Fritz Zwicky first pointed out the role of dark matter in clusters, Rashid Sunyaev predicted the scattering of CMB photons off hot gas, and Alan Guth laid the foundations of inflationary theory.
Implications
The results demonstrate that future surveys can turn constraining primordial non-Gaussianity into a multi-instrumental challenge, much like today's combination of cosmic microwave background, supernova, and baryon acoustic oscillation data for the standard cosmological model. This paves the way for detailed testing of inflation scenarios involving multiple fields and heavy particles that were previously nearly inaccessible. Moreover, joint analysis automatically reduces systematic errors associated with accounting for ordinary baryonic matter.
Future development
In the coming years, precision will improve with higher resolution maps (NSIDE 4096) and the inclusion of fourth-order moments and topological statistics. Additional observables such as the kinematic Sunyaev-Zeldovich effect and fast radio bursts will help pin down the baryonic physics model even more tightly. New simulations with better resolution will refine the contribution of low-mass clusters, and machine learning will efficiently extract information from the full variety of maps.
Impact
The work will impact cosmology, hot gas astrophysics, and particle physics, linking the observed structure of the universe to processes at Planck-scale energies.
Next steps
Next steps include applying the developed pipeline to real data from DES, DECADE, and eROSITA DR1, as well as integrating it into the analysis of future projects like SPHEREx and Athena.
Key open problems
The research directly addresses the problem of the nature of inflation and physics beyond the Standard Model, including mass hierarchy and dark energy, since the shape of non-Gaussianity may point to new fields and interactions in the early universe.
The knowledge machine's view
This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.
What the work does
The authors generated synthetic sky maps from N-body simulations of Ulagam with modified initial conditions for six models of primordial non-Gaussianity. Using the semi-analytic model BaryonForge, they simulated weak lensing, the Sunyaev-Zel'dovich effect, and X-ray emission, adding realistic noise and backgrounds. They then computed second- and third-order moments of the smoothed fields and used the Fisher method to estimate expected constraints on f_NL for future surveys.
Where its strength is
The strength of this work lies in the combined use of three complementary observables to break degeneracies between cosmological and baryonic parameters. This significantly improves constraints on primordial non-Gaussianity compared to using lensing alone.
Where it could go next
- It might be interesting to see how the non-Gaussianity constraints from these multi-wavelength maps compare to limits obtainable from gravitational-wave observations of bubble collisions in multi-field inflation scenarios.based on: Giant Bubbles Beyond the Horizon: Their Collisions Are Heard in Gravitational Waves arXiv:2607.24734v1
- A nearby work connects multi-field inflation models with the birth of the universe from quantum tunneling to predictions for gravitational waves; perhaps similar scenarios could produce distinctive non-Gaussianity signatures that could be searched for in your maps.based on: Uneven Tunnel Chain: How Gravity Rewrote the History of the Early Universe arXiv:2607.27193v1
- Another adjacent work considers vacuum decay as a mechanism for the birth of the universe; such a scenario could lead to specific non-Gaussianity that could be tested with your approach combining lensing, SZ, and X-ray.based on: Dark Start: How Vacuum Decay Birthed the Big Bang arXiv:2606.06587v1
Worth trying
- Try applying the same moment analysis to maps generated from simulations with initial conditions corresponding to bubble or tunneling scenarios, to assess how strongly they differ from the considered non-Gaussianity models.
- Investigate how adding information about cluster shapes (e.g., via SZ and X-ray profiles) could improve constraints on non-Gaussianity within your method.
- Explore whether correlations between maps on large scales can be used to isolate the contribution from superhorizon modes, which may be linked to non-minimal inflationary scenarios.
Why this may matter
Extending the analysis to a broader class of inflationary models that predict non-Gaussianity could help narrow down theoretical possibilities and bring us closer to understanding physics at ultra-high energies inaccessible to laboratory experiments.
Nearby works in our archive
- Ancient Spacetime Ripples Reshape the Face of Galaxiessimilarity 0.68
- The Cosmic Recipe Gets More Complicatedsimilarity 0.67
- Giant Bubbles Beyond the Horizon: Their Collisions Are Heard in Gravitational Wavessimilarity 0.66
- The Cosmic Loom: The Puzzle of Missing Antimatter and Inflationsimilarity 0.65
- Tiny Black Holes Don’t Warp the Cosmossimilarity 0.65
- A hiccup in the early Universe gave rise to a gravitational humsimilarity 0.65
- Uneven Tunnel Chain: How Gravity Rewrote the History of the Early Universesimilarity 0.64
- Dark Start: How Vacuum Decay Birthed the Big Bangsimilarity 0.63
The closest work in the archive sits at 0.68 by meaning — the usual level for our corpus.
Original abstract
as written by the authors arXiv:2607.06692
Primordial non-Gaussianities (PNGs) are a broad class of features in the initial density field that are connected to the particle physics of the early Universe. Measuring the amplitude of these features directly constrains fundamental physics from these earliest epochs and lends insight into energy scales that cannot be probed with terrestrial experiments. Using a new class of simulation methods, we propagate these signatures to their impact on the formation of non-linear structure and quantify the constraining power in non-Gaussian summary statistics of weak lensing, thermal Sunyaev Zeldovich (tSZ), and X-ray surveys. We use semi-analytic baryon models that consistently include astrophysical effects across all these observables, and use foreground modeling approaches that explicitly fold in correlations between the various components. We find that the tSZ and X-ray fields have significant information about PNGs, and additionally can help self-calibrate a broad set of nuisance parameters/models by breaking parameter degeneracies. Using the second and third moments of the lensing, tSZ, and X-ray fields, we find a factor of 2 improvement in PNG constraints relative to using lensing alone. Larger improvements are expected when including more scales and other complementary summary statistics. Our multi-wavelength map maker can be found at https://github.com/DhayaaAnbajagane/Vaanam. The simulations and software pipelines used in this analysis are publicly available.
Although the energy of inflation can reach $10^{15}$ GeV, a trillion times higher than the Large Hadron Collider, its imprint on galaxy distribution is today measured to within a few percent. It's like reconstructing the shape of a thrown pebble from the ripples on a pond's surface.
Key numbers
- constraint improvement: 2x when adding tSZ and X-ray
- inflation energy: ~1015 GeV
- LSST survey area: 14,000 sq. degrees
- number of simulations: over 2000 full-sky maps
- map resolution: NSIDE=1024 (pixels ~3.2')