In a vacuum first order phase transition gravitational waves are generated from collision of bubbles of the true vacuum. The spectrum from such collisions takes the form of a broken power law. We consider a toy model for such a phase transition where the dynamics of the scalar field depends on a single parameter $\bar{\lambda}$, which controls how thin the bubble wall is at nucleation and how close to degenerate the vacua are relative to the barrier. We extend on our previous work by performing a series of simulations with a range of $\bar{\lambda}$. The peak of the gravitational wave power spectrum varies by up to a factor of $1.3$, which is probably an unobservable effect. We find that the ultraviolet (UV) power law in the gravitational wave spectrum becomes steeper as $\bar{\lambda} \rightarrow 0$, varying between $k^{-1.4}$ and $k^{-2.3}$ for the $\bar{\lambda}$ considered. This provides some evidence that the form of the underlying effective potential of a vacuum first order phase transition could be determined from the gravitational wave spectrum it produces.

We present results from large-scale numerical simulations of a first order thermal phase transition in the early universe, in order to explore the shape of the acoustic gravitational wave and the velocity power spectra. We compare the results with the predictions of the recently proposed sound shell model. For the gravitational wave power spectrum, we find that the predicted $k^{-3}$ behaviour, where $k$ is the wavenumber, emerges clearly for detonations. The power spectra from deflagrations show similar features, but exhibit a steeper high-$k$ decay and an extra feature not accounted for in the model. There are two independent length scales: the mean bubble separation and the thickness of the sound shell around the expanding bubble of the low temperature phase. It is the sound shell thickness which sets the position of the peak of the power spectrum. The low wavenumber behaviour of the velocity power spectrum is consistent with a causal $k^{3}$, except for the thinnest sound shell, where it is steeper. We present parameters for a simple broken power law fit to the gravitational wave power spectrum for wall speeds well away from the speed of sound where this form can be usefully applied. We examine the prospects for the detection, showing that a LISA-like mission has the sensitivity to detect a gravitational wave signal from sound waves with an RMS fluid velocity of about $0.05c$, produced from bubbles with a mean separation of about $10^{-2}$ of the Hubble radius. The shape of the gravitational wave power spectrum depends on the bubble wall speed, and it may be possible to estimate the wall speed, and constrain other phase transition parameters, with an accurate measurement of a stochastic gravitational wave background.

We investigate the potential for observing gravitational waves from cosmological phase transitions with LISA in light of recent theoretical and experimental developments. Our analysis is based on current state-of-the-art simulations of sound waves in the cosmic fluid after the phase transition completes. We discuss the various sources of gravitational radiation, the underlying parameters describing the phase transition and a variety of viable particle physics models in this context, clarifying common misconceptions that appear in the literature and identifying open questions requiring future study. We also present a web-based tool, PTPlot, that allows users to obtain up-to-date detection prospects for a given set of phase transition parameters at LISA.

We have performed the first 3-dimensional simulations of strong first-order thermal phase transitions in the early Universe. For deflagrations, we find that the rotational component of the fluid velocity increases as the transition strength is increased. For detonations, however, the rotational velocity component remains constant and small. We also find that the efficiency with which kinetic energy is transferred to the fluid falls below theoretical expectations as we increase the transition strength. The probable origin of the kinetic energy deficit is the formation of reheated droplets of the metastable phase during the collision, slowing the bubble walls. The rate of increase in the gravitational wave energy density for deflagrations in strong transitions is suppressed compared to that predicted in earlier work. This is largely accounted for by the reduction in kinetic energy. Current modelling therefore substantially overestimates the gravitational wave signal for strong transitions with deflagrations, in the most extreme case by a factor of $10^{3}$. Detonations, however, are less affected.

We present the first end-to-end nonperturbative analysis of the gravitational wave power spectrum from a thermal first-order electroweak phase transition (EWPT), using the framework of dimensionally reduced effective field theory and pre-existing nonperturbative simulation results. We are able to show that a first-order EWPT in any beyond the Standard Model (BSM) scenario that can be described by a Standard Model-like effective theory at long distances will produce gravitational wave signatures too weak to be observed at existing and planned detectors. This implies that colliders are likely to provide the best chance of exploring the phase structure of such theories, while transitions strong enough to be detected at gravitational wave experiments require either previously neglected higher-dimension operators or light BSM fields to be included in the dimensionally reduced effective theory and therefore necessitate dedicated nonperturbative studies. As a concrete application, we analyze the real singlet-extended Standard Model and identify regions of parameter space with single-step first-order transitions, comparing our findings to those obtained using a fully perturbative method. We discuss the prospects for exploring the electroweak phase diagram in this model at collider and gravitational wave experiments in light of our nonperturbative results.

Models of symmetry breaking in the early universe can produce networks of cosmic strings threading 't Hooft-Polyakov monopoles. In certain cases there is a larger global symmetry group and the monopoles split into so-called semipoles. These networks are all known as cosmic necklaces. We carry out large-scale field theory simulations of the simplest model containing these objects, confirming that the energy density of networks of cosmic necklaces approaches scaling, i.e. that it remains a constant fraction of the background energy density. The number of monopoles per unit comoving string length is constant, meaning that the density fraction of monopoles decreases with time. Where the necklaces carry semipoles rather than monopoles, we perform the first simulations large enough to demonstrate that they also maintain a constant number per unit comoving string length. We also compare our results to a number of analytical models of cosmic necklaces, finding that none explains our results. We put forward evidence that annihilation of poles on the strings is controlled by a diffusive process, a possibility not considered before. The observational constraints derived in our previous work for necklaces with monopoles can now be safely applied to those with semipoles as well.

We conduct large scale numerical simulations of gravitational wave production at a first order vacuum phase transition. We find a power law for the gravitational wave power spectrum at high wavenumber which falls off as $k^{-1.5}$ rather than the $k^{-1}$ produced by the envelope approximation. The peak of the power spectrum is shifted to slightly lower wave numbers from that of the envelope approximation. The envelope approximation reproduces our results for the peak power less well, agreeing only to within an order of magnitude. After the bubbles finish colliding the scalar field oscillates around the true vacuum. An additional feature is produced in the UV of the gravitational wave power spectrum, and this continues to grow linearly until the end of our simulation. The additional feature peaks at a length scale close to the bubble wall thickness and is shown to have a negligible contribution to the energy in gravitational waves, providing the scalar field mass is much smaller than the Planck mass.

We review the production of gravitational waves by an electroweak first order phase transition. The resulting signal is a good candidate for detection at next-generation gravitational wave detectors, such as LISA. Detection of such a source of gravitational waves could yield information about physics beyond the Standard Model that is complementary to that accessible to current and near-future collider experiments. We summarise efforts to simulate and model the phase transition and the resulting production of gravitational waves.

We perform the first numerical simulations of necklaces in a non-Abelian gauge theory. Necklaces are composite classical solutions which can be interpreted as monopoles trapped on strings, rather generic structures in a Grand Unified Theory. We generate necklaces from random initial conditions, modelling a phase transition in the early Universe, and study the evolution. For all cases, we find that the necklace system shows scaling behaviour similar to that of a network of ordinary cosmic strings. Furthermore, our simulations indicate that comoving distance between the monopoles or semipoles along the string asymptotes to a constant value at late times. This means that while the monopole-to-string energy density ratio decreases as the inverse of the scale factor, a horizon-size length of string has a large number of monopoles, significantly affecting the dynamics of string loops. We argue that gravitational wave bounds from millisecond pulsar timing on the string tension in the Nambu-Goto scenario are greatly relaxed.

We study the properties of classical vortex solutions in a non-Abelian gauge theory. A system of two adjoint Higgs fields breaks the SU(2) gauge symmetry to $Z_2$, producing 't Hooft-Polyakov monopoles trapped on cosmic strings, termed beads; there are two charges of monopole and two degenerate string solutions. The strings break an accidental discrete $Z_2$ symmetry of the theory, explaining the degeneracy of the ground state. Further symmetries of the model, not previously appreciated, emerge when the masses of the two adjoint Higgs fields are degenerate. The breaking of the enlarged discrete symmetry gives rise to additional string solutions and splits the monopoles into four types of `semipole': kink solutions that interpolate between the string solutions, classified by a complex gauge invariant magnetic flux and a $Z_4$ charge. At special values of the Higgs self-couplings, the accidental symmetry broken by the string is continuous, giving rise to supercurrents on the strings. The SU(2) theory can be embedded in a wide class of Grand Unified Theories, including SO(10). We argue that semipoles and supercurrents are generic on GUT strings.

In real-time lattice simulations of cosmic strings in the Abelian Higgs model, the broken translational invariance introduces lattice artefacts; relativistic strings therefore decelerate and radiate. We introduce two different methods to construct a moving string on the lattice, and study in detail the lattice effects on moving strings. We find that there are two types of lattice artefact: there is an effective maximum speed with which a moving string can be placed on the lattice, and a moving string also slows down, with the deceleration approximately proportional to the exponential of the velocity. To mitigate this, we introduce and study an improved discretisation, based on the tree-level L\"{u}scher-Weisz action, which is found to reduce the deceleration by an order of magnitude, and to increase the string speed limit by an amount equivalent to halving the lattice spacing. The improved algorithm is expected to be very useful for 3D simulations of cosmic strings in the early universe, where one wishes to simulate as large a volume as possible.

We study the envelope approximation and its applicability to first-order phase transitions in the early universe. We demonstrate that the power laws seen in previous studies exist independently of the nucleation rate. We also compare the envelope approximation prediction to results from large-scale phase transition simulations. For phase transitions where the contribution to gravitational waves from scalar fields dominates over that from the coupled plasma of light particles, the envelope approximation is in agreement, giving a power spectrum of the same form and order of magnitude. In all other cases the form and amplitude of the gravitational wave power spectrum is markedly different and new techniques are required.

We investigate the potential for the eLISA space-based interferometer to detect the stochastic gravitational wave background produced by strong first-order cosmological phase transitions. We discuss the resulting contributions from bubble collisions, magnetohydrodynamic turbulence, and sound waves to the stochastic background, and estimate the total corresponding signal predicted in gravitational waves. The projected sensitivity of eLISA to cosmological phase transitions is computed in a model-independent way for various detector designs and configurations. By applying these results to several specific models, we demonstrate that eLISA is able to probe many well-motivated scenarios beyond the Standard Model of particle physics predicting strong first-order cosmological phase transitions in the early Universe.

We study the resonant decay of the primordial Standard Model Higgs condensate after inflation into $SU(2)$ gauge bosons on the lattice. We find that the non-Abelian interactions between the gauge bosons quickly extend the momentum distribution towards high values, efficiently destroying the condensate after the onset of backreaction. For the inflationary scale $H = 10^8$ GeV, we find that 90% of the Higgs condensate has decayed after $n \sim 10$ oscillation cycles. This differs significantly from the Abelian case where, given the same coupling strengths, most of the condensate would persist after the resonance.

We present details of numerical simulations of the gravitational radiation produced by a first order thermal phase transition in the early universe. We confirm that the dominant source of gravitational waves is sound waves generated by the expanding bubbles of the low-temperature phase. We demonstrate that the sound waves have a power spectrum with a power-law form between the scales set by the average bubble separation (which sets the length scale of the fluid flow $L_\text{f}$) and the bubble wall width. The sound waves generate gravitational waves whose power spectrum also has a power-law form, at a rate proportional to $L_\text{f}$ and the square of the fluid kinetic energy density. We identify a dimensionless parameter $\tilde\Omega_\text{GW}$ characterising the efficiency of this "acoustic" gravitational wave production whose value is $8\pi\tilde\Omega_\text{GW} \simeq 0.8 \pm 0.1$ across all our simulations. We compare the acoustic gravitational waves with the standard prediction from the envelope approximation. Not only is the power spectrum steeper (apart from an initial transient) but the gravitational wave energy density is generically larger by the ratio of the Hubble time to the phase transition duration, which can be 2 orders of magnitude or more in a typical first order electroweak phase transition.

First order phase transitions in the early Universe generate gravitational waves, which may be observable in future space-based gravitational wave observatiories, e.g. the European eLISA satellite constellation. The gravitational waves provide an unprecedented direct view of the Universe at the time of their creation. We study the generation of the gravitational waves during a first order phase transition using large-scale simulations of a model consisting of relativistic fluid and an order parameter field. We observe that the dominant source of gravitational waves is the sound generated by the transition, resulting in considerably stronger radiation than earlier calculations have indicated.

We report on the first 3-dimensional numerical simulations of first-order phase transitions in the early universe to include the cosmic fluid as well as the scalar field order parameter. We calculate the gravitational wave (GW) spectrum resulting from the nucleation, expansion and collision of bubbles of the low-temperature phase, for phase transition strengths and bubble wall velocities covering many cases of interest. We find that the compression waves in the fluid continue to be a source of GWs long after the bubbles have merged, a new effect not taken properly into account in previous modelling of the GW source. For a wide range of models the main source of the GWs produced by a phase transition is therefore the sound the bubbles make.

In this paper we use spontaneous flux production in annular superconductors to shed light on the Kibble-Zurek scenario. In particular, we examine the effects of finite size and external fields, neither of which is directly amenable to the KZ analysis. Supported by 1D and 3D simulations, the properties of a superconducting ring are seen to be well represented by analytic Gaussian approximations which encode the KZ scales indirectly. Experimental results for annuli in the presence of external fields corroborate these findings.