Numerical Simulations of Early Universe Sources of Gravitational Waves

28 Jul 2025, 09:00 → 15 Aug 2025, 19:00 Europe/Stockholm

Albano Building 3

Alberto Roper Pol (University of Geneva)

Venue

Nordita, Stockholm, Sweden

---

Scope

We are happy to announce the three-week program on numerical simulations of early Universe sources of gravitational waves at Nordita, Stockholm, that will take place from July 28 to August 15 of 2025.

The main objectives of the program are to review the different possible sources contributing to a cosmological background of gravitational waves, the state-of-the-art numerical codes that are available and used by the community, to discuss numerical techniques for these studies, and to discuss what are the priorities that we should address as a community in the coming years.

For this purpose, the program is divided into three weeks, covering the following potential sources of GWs in the early Universe: 

Week 1 (Jul 28-Aug 1): Inflation, (p)reheating, and primordial black holes

Week 2 (Aug 4-8): First order phase transitions and primordial turbulence

Week 3 (Aug 11-15): Topological defects: cosmic strings, domain walls, and others

Registration is now open and will close on March 15th.

To submit an abstract for a contributed talk, please provide the details while registering and then send an abstract in the call for abstracts section by the deadline, March 15th.

---

Confirmed Invited speakers

• Pierre Auclair (Louvain University)
• Jorge Baeza-Ballesteros (DESY)
• Jose Juan Blanco-Pillado (University of the Basque Country)
• Axel Brandenburg (Nordita)
• Malte Buschmann (University of Amsterdam)
• Angelo Caravano (IAP, Paris)
• Jose Ricardo Correia (University of Helsinki)
• Emanuela Dimastrogiovanni (University of Groningen)
• Matteo Fasiello (IFT)
• Marco Gorghetto (DESY)
• Mark Hindmarsh (University of Helsinki)
• Oksana Iarygina (Nordita)
• Ryusuke Jinno (Kobe University)
• Tina Kahniashvili (Carnegie Mellon University)
• Marek Lewicki (University of Warsaw)
• Joanes Lizarraga (University of the Basque Country)
• Swagat Mishra (University of Nottingham)
• Ilia Musco (University of Nova Gorica)
• Gerasimos Rigopoulos (Newcastle University)
• Henrique Rubira (LMU/Cambridge)
• Kari Rummukainen (University of Helsinki)
• Philipp Schicho (University of Geneva)
• Lara Sousa (University of Porto)
• Francisco Torrenti (University of Barcelona)
• Tanmay Vachaspati (Arizona State University)
• Jorinde van de Vis (CERN)
• Masahide Yamaguchi (Institute for Basic Science, Daejeon)

---

Accommodation

We will provide accommodation in Stockholm for some accepted applicants, prioritizing early-career researchers and based on travel fund needs. Please indicate if you would like Nordita to cover your accommodation in the registration form and we will confirm accommodation and travel support shortly after the registration deadline.

We have reserved some apartments at Biz Apartment Gärdet at a reduced price. If you would like to reserve one of these apartments and cover the costs on your own, please indicate so in the registration form.

---

Application/Registration

Registration is open and will close on March 15th. Registrants will receive an on-site/remote participation confirmation from the organizers shortly after this date.

---

Organizers

Chiara Caprini (CERN and University of Geneva)

Amelia Drew (University of Cambridge)

Daniel Figueroa (IFIC, Valencia)

Alberto Roper Pol (University of Geneva)

David Weir (University of Helsinki)

---

Code of conduct

All participants are expected to uphold our Code of Conduct, which aims to foster a harassment-free environment where everyone feels comfortable engaging in both scientific and social activities.

---

Sponsored by:

Code_of_Conduct_Nordita.pdf

Monday 28 July

1 Lectures on Inflation. Part I

Speaker: Matteo Fasiello

2 Studying Inhomogeneous Inflation with Numerical Relativity

Cosmic inflation is the leading paradigm for describing the early universe, addressing fundamental issues such as the horizon and flatness problems. However, a key unresolved question is the nature of its initial conditions. In this talk, I will explore how numerical relativity helps study inflationary spacetimes with inhomogeneous initial conditions, particularly in the presence of strong gravitational effects from large inhomogeneities. Full numerical simulations allow us to map out the phase space of initial conditions that lead to sufficient duration of slow roll inflation versus those that do not. The results strongly depend on the inflationary model, with a rule of thumb that the models with near- or super-Planckian characteristic scales are more robust to matter and geometric inhomogeneities than those with sub-Planckian scales. We mainly focus on the study of α-attractor models and our simulation results allow us to establish a lower bound on the tensor-to-scalar ratio r.

Speaker: Panagiotis Giannadakis (King's College London)

11:00 Coffee break

3 TBA

Speaker: Emanuela Dimastrogiovanni

12:30 Lunch break

4 AlbaNova Colloquium: Numerical Simulations of Early Universe Sources of Gravitational Waves

Speaker: Alberto Roper Pol (University of Geneva)

Tuesday 29 July

5 Lectures on Inflation. Part II

Speaker: Matteo Fasiello

6 What is the maximum temperature ever reached in the universe?

Gravitational waves are naturally sourced by hydrodynamical fluctuations in a thermal medium, as the one that filled the universe before recombination. Since the corresponding gravitational wave spectrum is expected to show rapid growth at high frequencies, f ∼ 1...1000 Hz, unprecedented prospects to detect these signals may be offered by the proposed Einstein Telescope. While the spectral shape is well understood, the peak amplitude is set by the plasma temperature at emission. We present a model-independent numerical method to estimate an upper bound for the maximal temperature reached after inflation and discuss the validity of the formalism across possible scale hierarchies.

Speaker: Simona Procacci (U. Geneva)

11:00 Coffee break

7 Towards Preheating after Inflation: Inflation Fragmentation, Oscillon Formation and Decay

The transition from cosmic inflation to the hot Big Bang, known as reheating, remains a key open question in cosmology. During its early stage, called preheating, the inflaton field decays explosively via parametric resonance into lighter bosonic offspring fields. However, when these external couplings are weak, strong self-interaction (cohesive force) causes the oscillating inflaton condensate to fragment, forming extremely long-lived scalar-field lumps known as oscillons. We investigate the conditions for oscillon formation during preheating, particularly in the presence of external couplings, within the class of inflationary potentials favored by the latest CMB observations. Using high-resolution (3+1)-dimensional lattice simulations on CosmoLattice platform, we systematically map the parameter space that supports oscillon formation. Our results suggest that preheating may have proceeded through both oscillon decay and the conventional decay of the inflaton condensate, offering new insights into the early reheating dynamics.

Speaker: Swagat Saurav Mishra (University of Nottingham, UK)

12:30 Lunch break

8 Stochastic Inflation in (Numerical) General Relativity

Speaker: Gerasimos Rigopoulos

Wednesday 30 July

9 Lectures on Inflation. Part III

Speaker: Matteo Fasiello

10 Primordial gravitational waves from from fully relativistic inflation simulations

Speaker: Paul Shellard

11:00 Coffee break

11 Backreaction and cosmic butterflies: what simulations can tell us about inflation

Speaker: Angelo Caravano

12:30 Lunch break

12 GW background signals from Inflation: lattice calculation

Speaker: Joanes Lizarraga

Thursday 31 July

13 Gravitational waves from axion inflation with non-Abelian gauge fields

Speaker: Oksana Iarygina (Stockholm University, Nordita)

14 Stochastic simulation of reheating and/or warm inflation

The late stage of the reheating process may be captured by a two-component approach, in which a self-interacting plasma has already attained local equilibrium, while the inflaton field is still far from equilibrium. This should be particularly suitable if the plasma contains non-Abelian gauge bosons, which are believed to equilibrate fast. We describe the foundations of such an approach, which can in principle be studied both in a linear and non-linear regime. Recent progress towards a gauge-invariant numerical implementation of the linear regime is summarized, and steps towards determining the curvature and tensor power spectra are outlined.

Speaker: Mikko Laine

11:00 Coffee break

15 Simulations of inflationary magnetogenesis

Speaker: Axel Brandenburg (Stockholm University, Nordita)

12:30 Lunch break

Friday 1 August

16 Gravitational waves from cosmic domain walls

Speaker: Francisco Torrenti

11:00 Coffee break

17 Numerical Simulations and Primordial Blak Holes

Speaker: ILIA Musco (Sapienza University of Rome, INFN)

12:30 Lunch break

Panel discussion week 1

Conveners: Chiara Caprini, Daniel Figueroa

Monday 4 August

18 TBA

Speaker: Marek Lewicki

19 Bubble Nucleation and Gravitational Waves from Strongly Coupled QFT's

Speaker: Nicklas Ramberg (SISSA Trieste)

11:00 Coffee break

20 Completing the perturbative program for cosmological phase transitions

Speaker: Philipp Schicho

12:30 Lunch break

Tuesday 5 August

21 The bubble wall velocity in first order phase transitions

Speaker: Jorinde van de Vis

22 Simulating cosmic bubbles on the lattice and in the lab

Bubble nucleation plays a pivotal role in many models of particle physics and the early Universe, and is a promising potential source of cosmological gravitational waves. However, we lack a satisfying theoretical understanding of this process, with existing approaches working only in imaginary (Euclidean) time, and relying on assumptions that have yet to be empirically tested. A promising route forward is to use cold-atom systems which undergo first-order phase transitions that are analogous to vacuum decay. In this talk, I will present recent theoretical work to understand this analogy using semiclassical lattice simulations, and will discuss possibilities and challenges for realising these analogues in the laboratory.

Speaker: Alex Jenkins (University of Cambridge)

11:00 Coffee break

23 TBA

Speaker: Ryusuke Jinno (DESY)

12:30 Lunch break

24 Origin of chiral magnetic effect, production of turbulence and generation of large-scale magnetic fields

In the standard model of particle physics, the chiral anomaly can occur in relativistic plasmas and plays an important role in the early Universe, proto-neutron stars, heavy-ion collisions, and quantum materials. It gives rise to a chiral magnetic effect if the number densities of left- and right-handed electrically charged fermions are
unequal. At high energies, the dynamics of a plasma with charged fermions can be described in terms of chiral magnetohydrodynamics.
We show that a chiral magnetic effect can result just from spatial fluctuations of the chemical potential, causing a chiral dynamo instability, magnetically driven turbulence, and ultimately a
generation of large-scale magnetic field through the magnetic alpha effect. This have consequences for the dynamics of certain high-energy plasmas, such as the early Universe. We discuss how the chiral magnetic effect can be a source for gravitational waves.

Authors: I. Rogachevskii , J. Schober, A. Brandenburg

Speaker: Prof. Igor Rogachevskii (Ben-Gurion University of the Negev)

Wednesday 6 August

25 Gas from first order phase transitions

Speaker: Mark Hindmarsh (University of Sussex)

26 Cold baryogenesis revisited

The matter-antimatter asymmetry of the Universe represents one of the main open questions in particle physics and cosmology. In this talk, we will present a novel realization of cold baryogenesis (a mechanism involving the formation and decay of topological defects associated with the gauge group of the Standard Model known as SU(2) textures) that relies on the out-of-equilibrium dynamics during a strong first order electroweak phase transition. By performing extensive lattice simulations of the Higgs doublet and gauge field dynamics, we evaluate the related Chern-Simons number production as well as the rate of baryon number violation, as a function of the parameters of the phase transition and the shape of the Higgs potential. We finally provide an estimate for the total baryon asymmetry generated this way.

Speaker: Simone Blasi (DESY Hamburg)

11:00 Coffee break

27 Phase transitions on the lattice

Speaker: Kari Rummukainen (University of Oulu)

12:30 Lunch break

28 Langer’s nucleation rate on the lattice

First-order phase transitions in the early universe provide a possible mechanism for producing observable gravitational waves. Predicting the gravitational wave spectrum requires accurate nucleation rate computations, as the rate determines the transition temperature and duration. These computations often rely on Langer's nucleation rate formula, which has long resisted validation in numerical simulations. Here, we present work that, for the first time, demonstrates agreement between Langer's formula and lattice simulations. Our findings clarify the conditions necessary for successful lattice simulations of nucleation and reveal insights into nucleation processes and the limitations of Langer's formula.

Speaker: Joonas Hirvonen (University of Nottingham)

Thursday 7 August

29 Higgsless simulations: Gravitation waves from sound waves

Speaker: Henrique Rubira

30 Gravitational waves from decaying sources: strong phase transitions

Speaker: Isak Stomberg (DESY)

11:00 Coffee break

31 What do we learn from pulsar timing data: testing early Universe physical processes

Speaker: Tina Kahniashvili (Carnegie Mellon University (USA) & Ilia State University (Georgia))

12:30 Lunch break

32 TBA

Speaker: Kenneth Marschall

Friday 8 August

33 Gravitational Waves from First-Order Phase Transitions

Speaker: Antonino Midiri

34 Gravitational wave production: the interplay between vortical and compressional motions.

Speaker: Madeline Salome

35 GWs from generalized fluid perturbations during first-order phase transitions

Given their weak interaction with different degrees of freedom, gravitational waves (GWs) offer fresh opportunities to probe the earliest moments of the Universe, and physics beyond the Standard Model. For instance, a first-order phase transition (FOPT) in the primordial plasma at the electroweak scale could emit a recognizable signal in the frequency range of the upcoming detectors like LISA. During a FOPT, bubbles of the stable phase nucleate, expand, and collide, generating perturbations that result in GWs. Thus, to model the consequent GW spectrum realistically, it is crucial to understand the fluid perturbations of the primordial plasma. Given the high energy and the relativistic speeds involved in strong enough phase transitions, the shape and amplitude of the power spectrum could be significantly affected by nonlinearities and the generated turbulence. While these topics have been generally studied in fluid dynamics for decades, their impact on GW spectra from FOPTs remains unclear. We investigate the production of vorticity and turbulence in the relativistic regime beyond linear theory. Using a semi-analytical approach, we report estimates of the GW spectrum, as well as the time scales and the strength of vorticity production when the initial field is purely compressional, as in the case of sound waves. This analysis helps constrain the applicability of models that approximate the fluid perturbations induced by bubble collisions as a superposition of the longitudinal modes, such as sound waves. Furthermore, it lays the foundation for more detailed numerical studies of the relativistic regime in the future.

Speaker: Deepen Garg (University of Bonn)

11:00 Coffee break

36 Gravitational Effects on Sound Waves: A Perturbative Approach for Large Bubbles in Cosmological First-Order Phase Transitions

If a cosmological first-order phase transition takes place over the cosmological time scale, gravity must affect the profile of sound waves. To investigate such a regime beyond the self-similarity, we combine a hydrodynamic scheme in the presence of gravity with a fluid computation scheme under energy injection from the bubble wall.
In this talk, I present the results of our (1+1)d hydrodynamic simulations around an expanding bubble in the cosmological background. We find that gravitational effects lead to a thinner fluid shell and reduce the kinetic energy budget, in qualitative agreement with previous discussions on late-time fluid behavior in an expanding universe. Moreover, our simulations reveal the development of sub-structures in the fluid profile for accelerating bubble walls. We also discuss the potential impact of these effects on the broadening of the SGWB spectral plateau.

Speaker: Jun'ya Kume (University of Padova, INFN Padova)

37 Generating gravitational wave spectra from equations of state for phase transitions in the early universe

The Sound Shell Model provides a computationally efficient way of calculating the gravitational wave spectra of first-order cosmological phase transitions, reproducing the results of lattice simulations for intermediate strength transitions. The Sound Shell Model has been encapsulated in the Python-based simulation framework PTtools, which enables easy generation of the gravitational wave spectra. Using PTtools one can simulate hundreds of spectra in a matter of minutes on a laptop. This enables charting the full parameter space.

The vast majority of existing simulations have been based on the bag model equation of state, which assumes the ultrarelativistic sound speed $c_s = \frac{1}{\sqrt{3}}$. As the latest addition, PTtools has been extended to include support for arbitrary equations of state and therefore for a temperature- and phase-dependent sound speed $c_s(T,\phi)$, extending the simulations beyond the ultrarelativistic assumption of the bag model.

In addition to PTtools, we have developed the web-based plotting utility
PTPlot, which makes generating the gravitational wave spectra possible through a web browser. Integrating PTtools with PTPlot enables easy access to the full Sound Shell Model for the research community.

Speaker: Mr Mika Mäki (University of Helsinki)

12:30 Lunch break

Panel discussion week 2

Conveners: Alberto Roper Pol (University of Geneva), David Weir (University of Helsinki)

Monday 11 August

38 Lectures on topological defects. Part I

Speaker: Tanmay Vachaspati (Arizona State University)

39 Gravitational Signatures of Domain Walls

In this talk, we present the results of simulations that explore the gravitational wave spectra produced by domain walls in a $\mathbf{Z}_2$ model, using the publicly available $CosmoLattice$ code. Focusing on the approach to scaling, we investigate the impact of various initial fluctuations and mass-to-Hubble ratios. We demonstrate that the Velocity-Dependent One-Scale (VOS) model accurately describes the evolution toward scaling after just a few e-folds, irrespective of the initial fluctuation conditions. Using a $2048^3$ grid, we compute the gravitational wave spectra for a domain wall system in an expanding Universe with different equations of state. Additionally, we conduct semi-analytical studies of the unequal time correlator (UTC) to refine our understanding of gravitational wave production and the dynamics of domain walls.

Speaker: Aäron Rase (Vrije Universiteit Brussel)

11:00 Coffee break

40 TBA

Speaker: Jose Juan Blanco-Pillado

12:30 Lunch break

Tuesday 12 August

41 Lectures on topological defects. Part II

Speaker: Tanmay Vachaspati (Arizona State University)

42 Simulating primordial black hole formation from domain wall collapse

Speaker: Matthew Elley

11:00 Coffee break

43 TBA

Speaker: Malte Buschmann

12:30 Lunch break

44 TBA

Speaker: Masahide Yamaguchi

Wednesday 13 August

45 Lectures on topological defects. Part III

Speaker: Tanmay Vachaspati (Arizona State University)

46 Gravitational Wave Stairway from Topological Defects.

Speaker: Nicklas Ramberg (SISSA Trieste)

11:00 Coffee break

47 TBA

Speaker: Lara Sousa

12:30 Lunch break

Thursday 14 August

48 Gravitational waves from axions

Speaker: Marco Gorghetto (DESY Hamburg)

11:00 Coffee break

49 TBA

Speaker: José Correia (University of Helsinki)

12:30 Lunch break

Friday 15 August

50 Gravitational wave and particle emission from cosmic string loops

Speaker: Jorge Baeza-Ballesteros (University of Valencia/IFIC)

11:00 Coffee break

51 Cosmic string loop fragmentation

Speaker: Pierre Auclair (APC)

12:30 Lunch break

Panel discussion week 3

Conveners: Amelia Drew (ICTP, Trieste), Daniel Figueroa