We study gravitational wave (GW) production in strongly supercooled cosmological phase transitions, taking particular care of models featuring a complex scalar field with a U$(1)$ symmetric potential. We perform lattice simulations of two-bubble collisions to properly model the scalar field gradients, and compute the GW spectrum sourced by them using thin-wall approximation in many-bubble simulations. We find that in the U$(1)$ symmetric case the low-frequency spectrum is $\propto\omega$ whereas for a real scalar field it is $\propto\omega^3$. In both cases the spectrum decays as $\omega^{-2}$ at high frequencies.

A standard expectation of primordial cosmological inflation is that it dilutes all relics created before its onset to unobservable levels. We present a counterexample to this expectation by demonstrating that a network of cosmic strings diluted by inflation can regrow to a level that is potentially observable today in gravitational waves~(GWs). In contrast to undiluted cosmic strings, whose primary GW signals are typically in the form of a stochastic GW background, the leading signal from a diluted cosmic string network can be distinctive bursts of GWs within the sensitivity reach of current and future GW observatories.

It is commonly assumed that the energy density of the Universe was dominated by radiation between reheating after inflation and the onset of matter domination 54,000 years later. While the abundance of light elements indicates that the Universe was radiation dominated during Big Bang Nucleosynthesis (BBN), there is scant evidence that the Universe was radiation dominated prior to BBN. It is therefore possible that the cosmological history was more complicated, with deviations from the standard radiation domination during the earliest epochs. Indeed, several interesting proposals regarding various topics such as the generation of dark matter, matter-antimatter asymmetry, gravitational waves, primordial black holes, or microhalos during a nonstandard expansion phase have been recently made. In this paper, we review various possible causes and consequences of deviations from radiation domination in the early Universe - taking place either before or after BBN - and the constraints on them, as they have been discussed in the literature during the recent years.

We survey systematically the general parametrisations of particle-physics models for a first-order phase transition in the early universe, including models with polynomial potentials both with and without barriers at zero temperature, and Coleman-Weinberg-like models with potentials that are classically scale-invariant. We distinguish three possibilities for the transition - detonations, deflagrations and hybrids - and consider sound waves and turbulent mechanisms for generating gravitational waves during the transitions in these models, checking in each case the requirement for successful percolation. We argue that in models without a zero-temperature barrier and in scale-invariant models the period during which sound waves generate gravitational waves lasts only for a fraction of a Hubble time after a generic first-order cosmological phase transition, whereas it may last longer in some models with a zero-temperature barrier that feature severe supercooling. We illustrate the implications of these results for future gravitational-wave experiments.

We study strongly supercooled cosmological phase transitions. We perform numerical lattice simulations of two-bubble collisions and demonstrate that, depending on the scalar potential, in the collision the field can either bounce to a false vacuum or remain oscillating around the true vacuum. We study if these cases can be distinguished from their gravitational wave signals and discuss the possibility of black hole formation in the bubble collisions.

We discuss the observability of circular polarisation of the stochastic gravitational-wave background (SGWB) generated by helical turbulence following a first-order cosmological phase transition, using a model that incorporates the effects of both direct and inverse energy cascades. We explore the strength of the gravitational-wave signal and the dependence of its polarisation on the helicity fraction, $\zeta_*$, the strength of the transition, $\alpha$, the bubble size, $R_*$, and the temperature, $T_*$, at which the transition finishes. We calculate the prospective signal-to-noise ratios of the SGWB strength and polarisation signals in the LISA experiment, exploring the parameter space in a way that is minimally sensitive to the underlying particle physics model. We find that discovery of SGWB polarisation is generally more challenging than measuring the total SGWB signal, but would be possible for appropriately strong transitions with large bubble sizes and a substantial polarisation fraction.

We outline the experimental concept and key scientific capabilities of AION (Atom Interferometer Observatory and Network), a proposed UK-based experimental programme using cold strontium atoms to search for ultra-light dark matter, to explore gravitational waves in the mid-frequency range between the peak sensitivities of the LISA and LIGO/Virgo/ KAGRA/INDIGO/Einstein Telescope/Cosmic Explorer experiments, and to probe other frontiers in fundamental physics. AION would complement other planned searches for dark matter, as well as probe mergers involving intermediate mass black holes and explore early universe cosmology. AION would share many technical features with the MAGIS experimental programme in the US, and synergies would flow from operating AION in a network with this experiment, as well as with other atom interferometer experiments such as MIGA, ZAIGA and ELGAR. Operating AION in a network with other gravitational wave detectors such as LIGO, Virgo and LISA would also offer many synergies.

We propose in this White Paper a concept for a space experiment using cold atoms to search for ultra-light dark matter, and to detect gravitational waves in the frequency range between the most sensitive ranges of LISA and the terrestrial LIGO/Virgo/KAGRA/INDIGO experiments. This interdisciplinary experiment, called Atomic Experiment for Dark Matter and Gravity Exploration (AEDGE), will also complement other planned searches for dark matter, and exploit synergies with other gravitational wave detectors. We give examples of the extended range of sensitivity to ultra-light dark matter offered by AEDGE, and how its gravitational-wave measurements could explore the assembly of super-massive black holes, first-order phase transitions in the early universe and cosmic strings. AEDGE will be based upon technologies now being developed for terrestrial experiments using cold atoms, and will benefit from the space experience obtained with, e.g., LISA and cold atom experiments in microgravity. This paper is based on a submission (v1) in response to the Call for White Papers for the Voyage 2050 long-term plan in the ESA Science Programme. ESA limited the number of White Paper authors to 30. However, in this version (v2) we have welcomed as supporting authors participants in the Workshop on Atomic Experiments for Dark Matter and Gravity Exploration held at CERN: ({\tt https://indico.cern.ch/event/830432/}), as well as other interested scientists, and have incorporated additional material.

We study the generation of intergalactic magnetic fields in two models for first-order phase transitions in the early Universe that have been studied previously in connection with the generation of gravitational waves (GWs): the Standard Model supplemented by an $|H|^6$ operator (SM+$H^6$) and a classically scale-invariant model with an extra gauged U(1) $B - L$ symmetry (SM$_{B-L}$). We consider contributions to magnetic fields generated by bubble collisions and by turbulence in the primordial plasma, and we consider the hypotheses that helicity is seeded in the gauge field or kinetically. We study the conditions under which the intergalactic magnetic fields generated may be larger than the lower bounds from blazar observations, and correlate them with the observability of GWs and possible collider signatures. In the SM+$H^6$ model bubble collisions alone cannot yield large enough magnetic fields, whereas turbulence may do so. In the SM$_{B-L}$ model bubble collisions and turbulence may both yield magnetic fields above the blazar bound unless the B$-$L gauge boson is very heavy. In both models there may be observable GW and collider signatures if sufficiently large magnetic fields are generated.

Cosmic string networks offer one of the best prospects for detection of cosmological gravitational waves (GWs). The combined incoherent GW emission of a large number of string loops leads to a stochastic GW background (SGWB), which encodes the properties of the string network. In this paper we analyze the ability of the Laser Interferometer Space Antenna (LISA) to measure this background, considering leading models of the string networks. We find that LISA will be able to probe cosmic strings with tensions $G\mu \gtrsim \mathcal{O}(10^{-17})$, improving by about $6$ orders of magnitude current pulsar timing arrays (PTA) constraints, and potentially $3$ orders of magnitude with respect to expected constraints from next generation PTA observatories. We include in our analysis possible modifications of the SGWB spectrum due to different hypotheses regarding cosmic history and the underlying physics of the string network. These include possible modifications in the SGWB spectrum due to changes in the number of relativistic degrees of freedom in the early Universe, the presence of a non-standard equation of state before the onset of radiation domination, or changes to the network dynamics due to a string inter-commutation probability less than unity. In the event of a detection, LISA's frequency band is well-positioned to probe such cosmic events. Our results constitute a thorough exploration of the cosmic string science that will be accessible to LISA.

We present a formalism based on the functional Schr\"odinger equation to analyse time-dependent tunneling in quantum field theory at the semi-classical level. The full problem is reduced step by step to a finite dimensional quantum mechanical setup and solved using the WKB approximation. As an example, we consider tunneling from a homogeneous oscillating initial state in scalar quantum field theory.

We derive efficiency factors for the production of gravitational waves through bubble collisions and plasma-related sources in strong phase transitions, and find the conditions under which the bubble collisions can contribute significantly to the signal. We use lattice simulations to clarify the dependence of the colliding bubbles on their initial state. We illustrate our findings in two examples, the Standard Model with an extra $|H|^6$ interaction and a classically scale-invariant $U(1)_{\rm B-L}$ extension of the Standard Model. The contribution to the GW spectrum from bubble collisions is found to be negligible in the $|H|^6$ model, whereas it can play an important role in parts of the parameter space in the scale-invariant $U(1)_{\rm B-L}$ model. In both cases the sound-wave period is much shorter than a Hubble time, suggesting a significant amplification of the turbulence-sourced signal. We find, however, that the peak of the plasma-sourced spectrum is still produced by sound waves with the slower-falling turbulence contribution becoming important off-peak.

What is the maximum possible strength of a first-order electroweak phase transition and the resulting gravitational wave (GW) signal? While naively one might expect that supercooling could increase the strength of the transition to very high values, for strong supercooling the Universe is no longer radiation-dominated and the vacuum energy of the unstable minimum of the potential dominates the expansion, which can jeopardize the successful completion of the phase transition. After providing a general treatment for the nucleation, growth and percolation of broken phase bubbles during a first-order phase transition that encompasses the case of significant supercooling, we study the conditions for successful bubble percolation and completion of the electroweak phase transition in theories beyond the Standard Model featuring polynominal potentials. For such theories, these conditions set a lower bound on the temperature of the transition. Since the plasma cannot be significantly diluted, the resulting GW signal originates mostly from sound waves and turbulence in the plasma, rather than bubble collisions. We find the peak frequency of the GW signal from the phase transition to be generically $f \gtrsim 10^{-4}$ Hz. We also study the condition for GW production by sound waves to be long-lasting (GW source active for approximately a Hubble time), showing it is generally not fulfilled in concrete scenarios. Because of this the sound wave GW signal could be weakened, with turbulence setting in earlier, resulting in a smaller overall GW signal as compared to current literature predictions.

We perform a global fit of the extended scalar singlet model with a fermionic dark matter (DM) candidate. Using the most up-to-date results from the $\mathit{Planck}$ measured DM relic density, direct detection limits from the XENON1T (2018) experiment, electroweak precision observables and Higgs searches at colliders, we constrain the 7-dimensional model parameter space. We also find regions in the model parameter space where a successful electroweak baryogenesis (EWBG) can be viable. This allows us to compute the gravitational wave (GW) signals arising from the phase transition, and discuss the potential discovery prospects of the model at current and future GW experiments. Our global fit places a strong upper $\mathit{and}$ lower limit on the second scalar mass, the fermion DM mass and the scalar-fermion DM coupling. In agreement with previous studies, we find that our model can simultaneously yield a strong first-order phase transition and saturate the observed DM abundance. More importantly, the GW spectra of viable points can often be within reach of future GW experiments such as LISA, DECIGO and BBO.

Most cosmological models predict that the universe was hot and dense at the early stages of it's evolution. In this paper we analyse the influence of the thermal bath of Standard Model particles on the dynamics of cosmological Higgs domain walls. This manuscript poses an~extension of our earlier work in which we investigated the evolution of networks of Higgs domain walls neglecting the impact of temperature variation. Using the thermally corrected effective potential of Standard Model we have found that both the position of the local maximum $h_{max}$ separating minima and the width of domain walls strongly depend on temperature $T$. For temperatures higher than $10^{10}\; \textrm{GeV}$ they respectively increase proportionally and decrease inverse proportionally to the increasing temperature. Thus, the energy scale of the problem follows the value of temperature. Our numerical lattice simulations based on the PRS algorithm reveal that Higgs domain walls in the presence of the background thermal bath are highly unstable and decay shortly after formation. Moreover we have found that the fraction of horizons produced by inflation in which Higgs field expectation value is higher then $h_{max}$ needs to be very low in order for the evolution of the~network of the domain walls to end in the electroweak vacuum. This means that Higgs domain walls necessarily were very rare objects and their average energy density was very small. As a result, the domain walls can not significantly effect cosmological observables.

Many motivated extensions of the Standard Model predict the existence of cosmic strings. Gravitational waves originating from the dynamics of the resulting cosmic string network have the ability to probe many otherwise inaccessible properties of the early universe. In this study we show how the spectrum of gravitational waves from a cosmic string network can be used to test the equation of state of the early universe prior to Big Bang Nucleosynthesis (BBN). We also demonstrate that current and planned gravitational wave detectors such as LIGO, LISA, DECIGO/BBO, and ET/CE have the potential to detect signals of a non-standard pre-BBN equation of state and evolution of the early universe (e.g., early non-standard matter domination or kination domination) or new degrees of freedom active in the early universe beyond the sensitivity of terrestrial collider experiments and cosmic microwave background measurements.

We present preliminary results from a global study of the extended scalar singlet model with a fermionic dark matter (DM) candidate. In addition to requiring a successful electroweak baryogenesis, we combine constraints from the DM relic density, direct detection limits from PandaX-II experiment, electroweak precision observables and Higgs searches at colliders. In agreement with previous studies, we find that the model can simultaneously explain (at least a part of) the observed DM abundance and matter-antimatter asymmetry. The viable points often lead to strong gravitational wave (GW) signals that can potentially be probed at future GW experiments.

Reheating after inflation can proceed even if the inflaton couples to Standard Model (SM) particles only gravitationally. However, particle production during the transition between de-Sitter expansion and a decelerating Universe is rather inefficient and the necessity to recover the visible Universe leads to a non-standard cosmological evolution initially dominated by remnants of the inflaton field. We remain agnostic to the specific dynamics of the inflaton field and discuss a generic scenario in which its remnants behave as a perfect fluid with a general barotropic parameter $w$. Using CMB and BBN constraints we derive the allowed range of inflationary scales. We also show that this scenario results in a characteristic primordial Gravitational Wave (GW) spectrum which gives hope for observation in upcoming runs of LIGO as well as in other planned experiments.

Cosmic strings are generic cosmological predictions of many extensions of the Standard Model of particle physics, such as a $U(1)^\prime$ symmetry breaking phase transition in the early universe or remnants of superstring theory. Unlike other topological defects, cosmic strings can reach a scaling regime that maintains a small fixed fraction of the total energy density of the universe from a very early epoch until today. If present, they will oscillate and generate gravitational waves with a frequency spectrum that imprints the dominant sources of total cosmic energy density throughout the history of the universe. We demonstrate that current and future gravitational wave detectors, such as LIGO and LISA, could be capable of measuring the frequency spectrum of gravitational waves from cosmic strings and discerning the energy composition of the universe at times well before primordial nucleosynthesis and the cosmic microwave background where standard cosmology has yet to be tested. This work establishes a benchmark case that gravitational waves may provide an unprecedented, powerful tool for probing the evolutionary history of the very early universe.

Our main interest is the evolution of domain walls of the Higgs field in the early Universe. The aim of this paper is to understand how dynamics of Higgs domain walls could be influenced by yet unknown interactions from beyond the Standard Model. We assume that the Standard Model is valid up to certain, high, energy scale $\Lambda$ and use the framework of the effective field theory to describe physics below that scale. Performing numerical simulations with different values of the scale $\Lambda$ we are able to extend our previous analysis and determine its range of validity. We study domain walls interpolating between the physical electroweak vacuum and the vacuum appearing at very high field strengths. These domain walls could be formed from non-homogeneous configurations of the Higgs field produced by quantum fluctuations during inflation or thermal fluctuations during reheating. Our numerical simulations show that evolution of Higgs domain walls is rather insensitive to interactions beyond the Standard Model as long as masses of new particles are grater than $10^{12}\ \textrm{GeV}$. For lower values of $\Lambda$ the RG improved effective potential is strongly modified at field strengths crucial to the evolution of domain walls. For instance its minima become degenerate for $\Lambda$ around $10^{11}\ \textrm{GeV}$. We find that even in the case when the minima of the potential are nearly degenerate Higgs domain walls decayed shortly after their formation for generic initial conditions. On the other hand, in simulations with specifically chosen initial conditions Higgs domain walls can live longer and enter the scaling regime. We also determine the energy spectrum of gravitational waves produced by decaying domain walls of the Higgs field. For generic initial field configurations the amplitude of the signal is too small to be observed in present and planned detectors.

We analyse a simple extension of the SM with just an additional scalar singlet coupled to the Higgs boson. We discuss the possible probes for electroweak baryogenesis in this model including collider searches, gravitational wave and direct dark matter detection signals. We show that a large portion of the model parameter space exists where the observation of gravitational waves would allow detection while the indirect collider searches would not.

We present measurements of $\rho^0$, $\omega$ and K$^{*0}$ spectra in $\pi^{-} + $C production interactions at 158 GeV/c and $\rho^0$ spectra at 350 GeV/c using the NA61/SHINE spectrometer at the CERN SPS. Spectra are presented as a function of the Feynman's variable $x_\text{F}$ in the range $0 < x_\text{F} < 1$ and $0 < x_\text{F} < 0.5$ for 158 GeV/c and 350 GeV/c respectively. Furthermore, we show comparisons with previous measurements and predictions of several hadronic interaction models. These measurements are essential for a better understanding of hadronic shower development and for improving the modeling of cosmic ray air showers.

We study domain walls which can be created in the Standard Model under the assumption that it is valid up to very high energy scales. We focus on domain walls interpolating between the physical electroweak vacuum and the global minimum appearing at very high field strengths. The creation of the network which ends up in the electroweak vacuum percolating through the Universe is not as difficult to obtain as one may expect, although it requires certain tuning of initial conditions. Our numerical simulations confirm that such domain walls would swiftly decay and thus cannot dominate the Universe. We discuss the possibility of detection of gravitational waves produced in this scenario. We have found that for the standard cosmology the energy density of these gravitational waves is too small to be observed in present and planned detectors.

We consider various models realizing baryogenesis during the electroweak phase transition (EWBG). Our focus is their possible detection in future collider experiments and possible observation of gravitational waves emitted during the phase transition. We also discuss the possibility of a non-standard cosmological history which can facilitate EWBG. We show how acceptable parameter space can be extended due to such a modification and conclude that next generation precision experiments such as the ILC will be able to confirm or falsify many models realizing EWBG. We also show that, in general, collider searches are a more powerful probe than gravitational wave searches. However, observation of a deviation from the SM without any hints of gravitational waves can point to models with modified cosmological history that generically enable EWBG with weaker phase transition and thus, smaller GW signals.

In this paper we investigate the induced inflation with two flat regions: one Starobinsky-like plateau in big field regime and one shorter plateau around the saddle point of the Einstein frame potential. This multi-phase inflationary scenario can be used to solve the problem of classical cosmology. The inflation at the saddle-point plateau is consistent with the data and can have arbitrarily low scale. The results can be useful in the context of the Higgs-Axion relaxation and in a certain limit they are equivalent to the $\alpha$-attractors.

We study the impact on electroweak baryogenesis from a swifter cosmological expansion induced by dark matter. We detail the experimental bounds that one can place on models that realize it, and we investigate the modifications of these bounds that result from a non-standard cosmological history. The modifications can be sizeable if the expansion rate of the Universe increases by several orders of magnitude. We illustrate the impact through the example of scalar field dark matter, which can alter the cosmological history enough to enable a strong-enough first-order phase transition in the Standard Model when it is supplemented by a dimension six operator directly modifying the Higgs boson potential. We show that due to the modified cosmological history, electroweak baryogenesis can be realized, while keeping deviations of the triple Higgs coupling below HL-LHC sensitivies. The required scale of new physics to effectuate a strong-enough first order phase transition can change by as much as twenty percent as the expansion rate increases by six orders of magnitude.

The cosmic evolution of a dark matter model which behaves relativistically in the early Universe is explored. Dark matter is described as a complex scalar field, whose earliest evolution is characterized by a stiff equation of state ($p \simeq \rho$). In this phase, it is the dominant component in the Universe. We present constraints from Big Bang nucleosynthesis and primordial gravity waves from inflation. Also, we study how the associated enhanced expansion rate due to the stiff phase might facilitate a first-order electroweak symmetry breaking phase transition, in light of the recently measured value of the Higgs boson mass.