Marek Lewicki

We investigate how perturbative reheating after inflation modifies the primordial gravitational wave (GW) spectrum generated by cosmological phase transitions. Within a specific inflationary setup, we show that the thermodynamic quantities that control the phase transition depend on the effective equation of state of the cosmological background, which is itself set by the form of the inflaton potential during reheating. Assuming reheating proceeds via perturbative dissipation of the inflaton condensate into boson or fermion pairs, we find that phase transitions taking place in this epoch generally produce GW signals that are systematically suppressed compared with the standard radiation-dominated scenario. We also identify characteristic spectral features that may arise in this case, which could serve as distinctive signatures of the modified expansion history during reheating.

We study gravitational waves from supercooled cosmological first-order phase transitions. If such a transition is followed by inefficient reheating, the evolution history of the universe is modified by a period of early matter domination. This leaves an imprint on the predicted gravitational-wave spectra. Using Fisher analysis we show the parameter space in reach of upcoming gravitational wave observatories where reheating can be probed due to its impact on the stochastic background produced by the transition. We use both the simplified geometric parametrisation and the thermodynamical one explicitly including the decay rate of the field undergoing the transition as a parameter determining the spectrum. We show the expansion history following the transition can be probed provided the transition is very strong which is naturally realised in classically scale invariant models generically predicting supercooling. Moreover, in such a scenario the decay rate of the scalar undergoing the phase transition, a parameter most likely inaccessible to accelerators, can be determined through the spectrum analysis.

We investigate the hydrodynamic solutions for expanding bubbles in cosmological first-order phase transitions going beyond local thermal equilibrium approximation. Under the assumption of a tangensoidal field profile, we supplement the matching conditions with the entropy produced due to the interaction of the bubble wall with ambient plasma. This allows us to analytically compute the corresponding fluid profiles and find the bubble-wall velocity. We show that due to the entropy production, two stable solutions corresponding to a deflagration or hybrid and a detonation can coexist. Finally, we use numerical real-time simulations of bubble growth to show that in such cases it is typically the faster detonation solution which is realised. This effect can be explained in terms of the fluid profile not being fully formed into the predicted steady-state solution as the wall accelerates past this slower solution.

We initiate the LISA template databank for stochastic gravitational wave backgrounds sourced by cosmic strings. We include two templates, an analytical template, which enables more flexible searches, and a numerical template derived directly from large Nambu-Goto simulations of string networks. Using searches based on these templates, we forecast the parameter space within the reach of the experiment and the precision with which their parameters will be reconstructed, provided a signal is observed. The reconstruction permits probing the Hubble expansion and new relativistic DoF in the early universe. We quantify the impact that astrophysical foregrounds can have on these searches. Finally, we discuss the impact that these observations would have on our understanding of the fundamental models behind the string networks. Overall, we prove that LISA has great potential for probing cosmic string models and may reach tensions as low as Gμ = 10^-16 – 10^-17, which translates into energy scales of the order 10¹¹ GeV.

We study the evolution of growing vacuum bubbles. The bubble walls interact with the surrounding fluid and may, consequently, reach a terminal velocity. If the mean free path of the particles in the fluid is much shorter than the bubble wall thickness, the fluid is locally in thermal equilibrium and the wall’s terminal velocity can be determined by entropy conservation. On the other hand, if local thermal equilibrium inside the wall cannot be maintained, the wall velocity can be estimated from the pressure impacted by ballistic particle dynamics at the wall. We find that the latter case leads to slightly slower bubble walls. Expectedly, we find the largest differences in the terminal velocity when the fluid is entirely ballistic. This observation indicates that the non-equilibrium effects inside walls are relevant. To study bubble evolution, we perform hydrodynamic lattice simulations in the case of local thermal equilibrium and N-body simulations in the ballistic case to investigate the dynamical effects during expansion. Both simulations show that even if a stationary solution exists in theory it may not be reached depending on the dynamics of the accelerating bubble walls.

Slow first-order phase transitions generate large inhomogeneities that can lead to the formation of primordial black holes. We show that the gravitational wave spectrum then consists of a primary component sourced by bubble collisions and a secondary one induced by large perturbations. The latter gives the dominant peak if β/H_{0}<12, impacting, in particular, the interpretation of the recent pulsar timing array data. The gravitational wave signal associated with a particular primordial black hole population is stronger than in typical scenarios because of a negative non-Gaussianity of the perturbations, and it has a distinguishable shape with two peaks.

We study the evolution of growing vacuum bubbles. The bubble walls interact with the surrounding fluid and may, consequently, reach a terminal velocity. If the mean free path of the particles in the fluid is much shorter than the bubble wall thickness, the fluid is locally in thermal equilibrium and the wall's terminal velocity can be determined by entropy conservation. On the other hand, if local thermal equilibrium inside the wall cannot be maintained, the wall velocity can be estimated from the pressure impacted by ballistic particle dynamics at the wall. We find that the latter case leads to slightly slower bubble walls. Expectedly, we find the largest differences in the terminal velocity when the fluid is entirely ballistic. This observation indicates that the non-equilibrium effects inside walls are relevant. To study bubble evolution, we perform hydrodynamic lattice simulations in the case of local thermal equilibrium and $N$-body simulations in the ballistic case to investigate the dynamical effects during expansion. Both simulations show that even if a stationary solution exists in theory it may not be reached depending on the dynamics of the accelerating bubble walls.