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 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 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.

With the physical Higgs mass the Standard Model symmetry restoration phase transition is a smooth cross-over. We study the thermodynamics of the cross-over using numerical lattice Monte Carlo simulations of an effective SU(2) X U(1) gauge + Higgs theory, significantly improving on previously published results. We measure the Higgs field expectation value, thermodynamic quantities like pressure, energy density, speed of sound and heat capacity, and screening masses associated with the Higgs and Z fields. While the cross-over is smooth, it is very well defined with a width of only approximately 5 GeV. We measure the cross-over temperature from the maximum of the susceptibility of the Higgs condensate, with the result $T_c = 159.5 \pm 1.5$ GeV. Outside of the narrow cross-over region the perturbative results agree well with non-perturbative ones.

We use large-scale lattice simulations to compute the rate of baryon number violating processes (the sphaleron rate), the Higgs field expectation value, and the critical temperature in the Standard Model across the electroweak phase transition temperature. While there is no true phase transition between the high-temperature symmetric phase and the low-temperature broken phase, the cross-over is sharply defined at $T_c = (159\pm 1)$\,GeV. The sphaleron rate in the symmetric phase ($T> T_c$) is $\Gamma/T^4 = (18\pm 3)\alpha_W^5$, and in the broken phase in the physically interesting temperature range $130\mbox{\,GeV} < T < T_c$ it can be parametrized as $\log(\Gamma/T^4) = (0.83\pm 0.01)T/{\rm GeV} - (147.7\pm 1.9)$. The freeze-out temperature in the early Universe, where the Hubble rate wins over the baryon number violation rate, is $T_* = (131.7\pm 2.3)$\,GeV. These values, beyond being intrinsic properties of the Standard Model, are relevant for e.g. low-scale leptogenesis scenarios.

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.

We propose the liquid-scintillator detector LENA (Low Energy Neutrino Astronomy) as a next-generation neutrino observatory on the scale of 50 kt. The outstanding successes of the Borexino and KamLAND experiments demonstrate the large potential of liquid-scintillator detectors in low-energy neutrino physics. LENA's physics objectives comprise the observation of astrophysical and terrestrial neutrino sources as well as the investigation of neutrino oscillations. In the GeV energy range, the search for proton decay and long-baseline neutrino oscillation experiments complement the low-energy program. Based on the considerable expertise present in European and international research groups, the technical design is sufficiently mature to allow for an early start of detector realization.

We calculate the tension $\sigma$ of the interface between the confined and deconfined phases by the histogram method in SU(3) lattice gauge theory for temporal extents of 4 and 6 using the recent high-statistics data by QCDPAX collaboration. The results are $\sigma/T_c^3 = 0.0292(22)$ and 0.0218(33) for $N_t=4$ and 6, respectively. The ratio $\sigma/T_c^3$ shows a scaling violation similar to that already observed for the latent heat $\latent$. However, we find that the physically interesting dimensionless combinations $(\sigma^{3}/\latent^2 T)^{1/2}$ and $\sigma T/ \latent$ scale within the statistical errors.