Cosmological phase transitions in the early Universe may produce relics in the form of a network of cosmic defects. Independently of the order of a phase transition, topology of the defects, and their global or gauge nature, the defects are expected to emit gravitational waves (GWs) as the network energy-momentum tensor adapts itself to maintaining {scaling}. We show that the evolution of any defect network (and for that matter any scaling source) emits a GW background with spectrum $\Omega_{\rm GW} \propto f^3$ for $f \ll f_0$, $\Omega_{\rm GW} \propto 1/f^2$ for $f_0 \lesssim f \lesssim f_{\rm eq}$, and $\Omega_{\rm GW} \propto~const$ (i.e.~exactly scale-invariant) for $f \gg f_{\rm eq}$, where $f_0$ and $ f_{\rm eq}$ denote respectively the frequencies corresponding to the present and matter-radiation equality horizons. This background represents an irreducible emission of GWs from any scaling network of cosmic defects, with its amplitude characterized only by the symmetry breaking scale and the nature of the defects. Using classical lattice simulations we calculate the GW signal emitted by defects created after the breaking of a global symmetry $O(N) \rightarrow O(N-1)$. We obtain the GW spectrum for $N$ between 2 and 20 with two different techniques: integrating over unequal time correlators of the energy momentum tensor, updating our previous work on smaller lattices, and for the first time, comparing the result with the real time evolution of the tensor perturbations sourced by the same defects. Our results validate the equivalence of the two techniques. Using CMB upper bounds on the defects' energy scale, we discuss the difficulty of detecting this GW background in the case of global defects.

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.

In the QCD axion dark matter scenario with post-inflationary Peccei-Quinn symmetry breaking, the number density of axions, and hence the dark matter density, depends on the length of string per unit volume at cosmic time $t$, by convention written $\zeta/t^2$. The expectation has been that the dimensionless parameter $\zeta$ tends to a constant $\zeta_0$, a feature of a string network known as scaling. It has recently been claimed that in larger numerical simulations $\zeta$ shows a logarithmic increase with time, while theoretical modelling suggests an inverse logarithmic correction. Either case would result in a large enhancement of the string density at the QCD transition, and a substantial revision to the axion mass required for the axion to constitute all of the dark matter. With a set of new simulations of global strings we compare the standard scaling (constant-$\zeta$) model to the logarithmic growth and inverse-logarithmic correction models. In the standard scaling model, by fitting to linear growth in the mean string separation $\xi = t/\sqrt{\zeta}$, we find $\zeta_0 = 1.19 \pm 0.20$. We conclude that the apparent corrections to $\zeta$ are artefacts of the initial conditions, rather than a property of the scaling network. The residuals from the constant-$\zeta$ (linear $\xi$) fit also show no evidence for logarithmic growth, restoring confidence that numerical simulations can be simply extrapolated from the Peccei-Quinn symmetry-breaking scale to the QCD scale. Re-analysis of previous work on the axion number density suggests that recent estimates of the axion dark matter mass in the post-inflationary symmetry-breaking scenario we study should be increased by about 50%.

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 calculate gravitational wave power spectra from first order early Universe phase transitions using the Sound Shell Model. The model predicts that the power spectrum depends on the mean bubble separation, the phase transition strength, the phase boundary speed, with the overall frequency scale set by the nucleation temperature. There is also a dependence on the time evolution of the bubble nucleation rate. The gravitational wave peak power and frequency are in good agreement with published numerical simulations, where bubbles are nucleated simultaneously. Agreement is particularly good for detonations, but the total power for deflagrations is predicted higher than numerical simulations show, indicating refinement of the model of the transfer of energy to the fluid is needed for accurate computations. We show how the time-dependence of the bubble nucleation rate affects the shape of the power spectrum: an exponentially rising nucleation rate produces higher amplitude gravitational waves at a longer wavelength than simultaneous nucleation. We present an improved fit for the predicted gravitational wave power spectrum in the form of a double broken power law, where the two breaks in the slope happen at wavenumber corresponding to the mean bubble separation and the thickness of the fluid shell surrounding the expanding bubbles, which in turn is related to the difference of the phase boundary speed from the speed of sound.

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 results from the first simulations of networks of Type I Abelian Higgs cosmic strings to include both matter and radiation eras and Cosmic Microwave Background (CMB) constraints. In Type I strings, the string tension is a slowly decreasing function of the ratio of the scalar and gauge mass-squared, $\beta$. We find that the mean string separation shows no dependence on $\beta$, and that the energy-momentum tensor correlators decrease approximately in proportion to the square of the string tension, with additional O(1) correction factors which asymptote to constants below $\beta \lesssim 0.01$. Strings in models with low self-couplings can therefore satisfy current CMB bounds at higher symmetry-breaking scales. This is particularly relevant for models where the gauge symmetry is broken in a supersymmetric flat direction, for which the effective self-coupling can be extremely small. If our results can be extrapolated to $\beta \simeq 10^{-15}$, even strings formed at $10^{16}$ GeV (approximately the grand unification scale in supersymmetric extensions of the Standard Model) can be compatible with CMB constraints.

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.

Making use of definitive new lattice computations of the Standard Model thermodynamics during the quantum chromodynamic (QCD) phase transition, we calculate the enhancement in the mass distribution of primordial black holes (PBHs) due to the softening of the equation of state. We find that the enhancement peaks at approximately $0.7M_\odot$, with the formation rate increasing by at least two orders of magnitude due to the softening of the equation of state at this time, with a range of approximately $0.3M_\odot<M<1.4M_\odot$ at full width half-maximum. PBH formation is increased by a smaller amount for PBHs with masses spanning a large range, $10^{-3}M_\odot<M_{\rm PBH}<10^{3}M_\odot$, which includes the masses of the BHs that LIGO detected. The most significant source of uncertainty in the number of PBHs formed is now due to unknowns in the formation process, rather than from the phase transition. A near scale-invariant density power spectrum tuned to generate a population with mass and merger rate consistent with that detected by LIGO should also produce a much larger energy density of PBHs with solar mass. The existence of BHs below the Chandresekhar mass limit would be a smoking gun for a primordial origin and they could arguably constitute a significant fraction of the cold dark matter density. They also pose a challenge to inflationary model building which seek to produce the LIGO BHs without overproducing lighter PBHs.

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 investigate cosmic string networks in the Abelian Higgs model using data from a campaign of large-scale numerical simulations on lattices of up to $4096^3$ grid points. We observe scaling or self-similarity of the networks over a wide range of scales, and estimate the asymptotic values of the mean string separation in horizon length units $\dot{\xi}$ and of the mean square string velocity $\bar v^2$ in the continuum and large time limits. The scaling occurs because the strings lose energy into classical radiation of the scalar and gauge fields of the Abelian Higgs model. We quantify the energy loss with a dimensionless radiative efficiency parameter, and show that it does not vary significantly with lattice spacing or string separation. This implies that the radiative energy loss underlying the scaling behaviour is not a lattice artefact, and justifies the extrapolation of measured network properties to large times for computations of cosmological perturbations. We also show that the core growth method, which increases the defect core width with time to extend the dynamic range of simulations, does not introduce significant systematic error. We compare $\dot{\xi}$ and $\bar v^2$ to values measured in simulations using the Nambu-Goto approximation, finding that the latter underestimate the mean string separation by about 25%, and overestimate $\bar v^2$ by about 10%. The scaling of the string separation implies that string loops decay by the emission of massive radiation within a Hubble time in field theory simulations, in contrast to the Nambu-Goto scenario which neglects this energy loss mechanism. String loops surviving for only one Hubble time emit much less gravitational radiation than in the Nambu-Goto scenario, and are consequently subject to much weaker gravitational wave constraints on their tension.

We present the first CMB power spectra from numerical simulations of the global O(N) linear $\sigma$-model with N = 2,3, which have global strings and monopoles as topological defects. In order to compute the CMB power spectra we compute the unequal time correlators (UETCs) of the energy-momentum tensor, showing that they fall off at high wave number faster than naive estimates based on the geometry of the defects, indicating non-trivial (anti-)correlations between the defects and the surrounding Goldstone boson field. We obtain source functions for Einstein-Boltzmann solvers from the UETCs, using a recent method that improves the modelling at the radiation- matter transition. We show that the interpolation function that mimics the transition is similar to other defect models, but not identical, confirming the non-universality of the interpolation function. The CMB power spectra for global strings and monopoles have the same overall shape as those obtained using the non-linear $\sigma$-model approximation, which is well captured by a large-N calculation. However, the amplitudes are larger than the large-N calculation predict, and in the case of global strings much larger: a factor of 20 at the peak. Finally we compare the CMB power spectra with the latest CMB data to put limits on the allowed contribution to the temperature power spectrum at multipole $\ell$ = 10 of 1.7% for global strings and 2.4% for global monopoles. These limits correspond to symmetry-breaking scales of 2.9x1015 GeV (6.3x1014 GeV with the expected logarithmic scaling of the effective string tension between the simulation time and decoupling) and 6.4x1015 GeV respectively. The bound on global strings is a significant one for the ultra-light axion scenario with axion masses ma 10-28 eV. These upper limits indicate that gravitational wave from global topological defects will not be observable at the GW observatory LISA.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 \mu$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations.

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 present cosmic microwave background (CMB) power spectra from recent numerical simulations of cosmic strings in the Abelian Higgs model and compare them to CMB power spectra measured by Planck. We obtain revised constraints on the cosmic string tension parameter $G\mu$. For example, in the $\Lambda$CDM model with the addition of strings and no primordial tensor perturbations, we find $G\mu < 2.0 \times 10^{-7}$ at 95% confidence, about 20% lower than the value obtained from previous simulations, which had 1/64 of the spatial volume. We investigate the source of the difference, showing that the main cause is an improved treatment of the string evolution across the radiation-matter transition. The increased computational volume also makes possible to simulate fully the physical equations of motion, in which the string cores shrink in comoving coordinates. This, and the larger dynamic range, changes the amplitude of the power spectra by only about 10%, demonstrating that field theory simulations of cosmic strings have now reached the required dynamic range for CMB calculations.

We report on the energy-momentum correlators obtained with recent numerical simulations of the Abelian Higgs model, essential for the computation of cosmic microwave background and matter perturbations of cosmic strings. Due to significant improvements both in raw computing power and in our parallel simulation framework, the dynamical range of the simulations has increased four-fold both in space and time, and for the first time we are able to simulate strings with a constant physical width in both the radiation and matter eras. The new simulations improve the accuracy of the measurements of the correlation functions at the horizon scale and confirm the shape around the peak. The normalization is slightly higher in the high wave-number tails, due to a small increase in the string density. We study for the first time the behaviour of the correlators across cosmological transitions, and discover that the correlation functions evolve adiabatically, ie the network adapts quickly to changes in the expansion rate. We propose a new method for constructing source functions for Einstein-Boltzmann integrators, comparing it with two other methods previously used. The new method is more consistent, easier to implement, and significantly more accurate.

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.

I give a short commentary on a seminal article by T W B Kibble in 1976, "Topology of cosmic domains and 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.

A model for the acoustic production of gravitational waves at a first order phase transition is presented. The source of gravitational radiation is the sound waves generated by the explosive growth of bubbles of the stable phase. The model assumes that the sound waves are linear and that their power spectrum is determined by the characteristic form of the sound shell around the expanding bubble. The predicted power spectrum has two length scales, the average bubble separation and the sound shell width when the bubbles collide. The peak of the power spectrum is at wavenumbers set by the sound shell width. For higher wavenumber $k$, the power spectrum decreases as $k^{-3}$. At wavenumbers below the inverse bubble separation, the power spectrum goes as $k^5$. For bubble wall speeds near the speed of sound where these two length scales are distinguished, there is an intermediate $k^{1}$ power law. The detailed dependence of the power spectrum on the wall speed and the other parameters of the phase transition raises the possibility of their constraint or measurement at a future space-based gravitational wave observatory such as eLISA.

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.

latfield2 is a C++ library designed to simplify writing parallel codes for solving partial differen- tial equations, developed for application to classical field theories in particle physics and cosmology. It is a significant rewrite of the latfield framework, moving from a slab domain decomposition to a rod decomposition, where the last two dimension of the lattice are scattered into a two dimensional process grid. Parallelism is implemented using the Message Passing Interface (MPI) standard, and hidden in the basic objects of grid-based simulations: Lattice, Site and Field. It comes with an integrated parallel fast Fourier transform, and I/O server class permitting computation to continue during the writing of large files to disk. latfield2 has been used for production runs on tens of thousands of processor elements, and is expected to be scalable to hundreds of thousands.

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 calculate the cosmic microwave background temperature bispectrum from cosmic strings, for the first time including the contributions from the last scattering surface, using a well-established Gaussian model for the string energy-momentum correlation functions, and a simplified model for the cosmic fluid. We check our approximation for the integrated Sachs-Wolfe (ISW) contribution against the bispectrum obtained from the full sky map of the cosmic string ISW signal used by the Planck team, obtaining good agreement. We validate our model for the last scattering surface contribution by comparing the predicted temperature power spectrum with that obtained from a full Boltzmann code treatment applied to the Unconnected Segment Model of a string network. We find that including the last scattering contribution has only a small impact on the upper limit on the string tension resulting from the bispectrum at Planck resolutions, and argue that the bispectrum is unlikely to be competitive with the power spectrum at any resolution.

We examine the production of dark matter by decaying topological defects in the high mass region $m_{\mathrm{DM}} \gg m_W$ of the Inert Doublet Model, extended with an extra U(1) gauge symmetry. The density of dark matter states (the neutral Higgs states of the inert doublet) is determined by the interplay of the freeze-out mechanism and the additional production of dark matter states from the decays of topological defects, in this case cosmic strings. These decays increase the predicted relic abundance compared to the standard freeze-out only case, and as a consequence the viable parameter space of the Inert Doublet Model can be widened substantially. In particular, for a given dark matter annihilation rate lower dark matter masses become viable. We investigate the allowed mass range taking into account constraints on the energy injection rate from the diffuse $\gamma$-ray background and Big Bang Nucleosynthesis, together with constraints on the dark matter properties coming from direct and indirect detection limits. For the Inert Doublet Model high-mass region, an inert Higgs mass as low as $\sim 200$ GeV is permitted. There is also an upper limit on string mass per unit length, and hence the symmetry breaking scale, from the relic abundance in this scenario. Depending on assumptions made about the string decays, the limits are in the range $10^{12}$ GeV to $10^{13}$ GeV.

In this work we discuss the possibility of cosmic defects being responsible for the B-mode signal measured by the BICEP2 collaboration. We also allow for the presence of other cosmological sources of B-modes such as inflationary gravitational waves and polarized dust foregrounds, which might contribute to or dominate the signal. On the one hand, we find that defects alone give a poor fit to the data points. On the other, we find that defects help to improve the fit at higher multipoles when they are considered alongside inflationary gravitational waves or polarized dust. Finally, we derive new defect constraints from models combining defects and dust. This proceeding is based on previous works [1,2].

We present the first calculation of the bispectrum of the matter perturbations induced by cosmic strings. The calculation is performed in two different ways: the first uses the unequal time correlators (UETCs) of the string network - computed using a Gaussian model previously employed for cosmic string power spectra. The second approach uses the wake model, where string density perturbations are concentrated in sheet-like structures whose surface density grows with time. The qualitative and quantitative agreement of the two gives confidence to the results. An essential ingredient in the UETC approach is the inclusion of compensation factors in the integration with the Green's function of the matter and radiation fluids, and we show that these compensation factors must be included in the wake model also. We also present a comparison of the UETCs computed in the Gaussian model, and those computed in the unconnected segment model (USM) used by the standard cosmic string perturbation package CMBACT. We compare numerical estimates for the bispectrum of cosmic strings to those produced by perturbations from an inflationary era, and discover that, despite the intrinsically non-Gaussian nature of string-induced perturbations, the matter bispectrum is unlikely to produce competitive constraints on a population of cosmic strings.

We analyse the possible contribution of topological defects to cosmic microwave anisotropies, both temperature and polarisation. We allow for the presence of both inflationary scalars and tensors, and of polarised dust foregrounds that may contribute to or dominate the B-mode polarisation signal. We confirm and quantify our previous statements that topological defects on their own are a poor fit to the B-mode signal. However, adding topological defects to a models with a tensor component or a dust component improves the fit around $\ell=200$. Fitting simultaneously to both temperature and polarisation data, we find that textures fit almost as well as tensors ($\Delta\chi^2 = 2.0$), while Abelian Higgs strings are ruled out as the sole source of the B-mode signal at low $\ell$. The 95% confidence upper limits on models combining defects and dust are $G\mu < 2.7\times 10^{-7}$ (Abelian Higgs strings), $G\mu < 9.8\times 10^{-7}$ (semilocal strings) and $G\mu < 7.3\times 10^{-7}$ (textures), a small reduction on the $Planck$ bounds. The most economical fit overall is obtained by the standard $\Lambda$CDM model with a polarised dust component.

We show that the B-mode polarization signal detected at low multipoles by BICEP2 cannot be entirely due to topological defects. This would be incompatible with the high-multipole B-mode polarization data and also with existing temperature anisotropy data. Adding cosmic strings to a model with tensors, we find that B-modes on their own provide a comparable limit on the defects to that already coming from Planck satellite temperature data. We note that strings at this limit give a modest improvement to the best-fit of the B-mode data, at a somewhat lower tensor-to-scalar ratio of $r \simeq 0.15$.

We study dark matter production by decaying topological defects, in particular cosmic strings. In topological defect or "top-down" (TD) scenarios, the dark matter injection rate varies as a power law with time with exponent $p-4$. We find a formula in closed form for the yield for all $p < 3/2$, which accurately reproduces the solution of the Boltzmann equation. We investigate two scenarios ($p=1$, $p=7/6$) motivated by cosmic strings which decay into TeV-scale states with a high branching fraction into dark matter particles. For dark matter models annihilating either by s-wave or p-wave, we find the regions of parameter space where the TD model can account for the dark matter relic density as measured by Planck. We find that topological defects can be the principal source of dark matter, even when the standard freeze-out calculation under-predicts the relic density and hence can lead to potentially large "boost factor" enhancements in the dark matter annihilation rate. We examine dark matter model-independent limits on this scenario arising from unitarity and discuss example model-dependent limits coming from indirect dark matter search experiments. In the four cases studied, the upper bound on $G\mu$ for strings with an appreciable channel into TeV-scale states is significantly more stringent than the current Cosmic Microwave Background limits.

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 demonstrate that any scaling source in the radiation era produces a background of gravitational waves with an exact scale-invariant power spectrum. Cosmic defects, created after a phase transition in the early Universe, are such a scaling source. We emphasise that the result is independent of the topology of the cosmic defects, the order of phase transition, and the nature of the symmetry broken, global or gauged. As an example, using large-scale numerical simulations, we calculate the scale invariant gravitational wave power spectrum generated by the dynamics of a global O(N) scalar theory. The result approaches the large N theoretical prediction as N^(-2), albeit with a large coefficient. The signal from global cosmic strings is O(100) times larger than the large N prediction.

We consider a class of models in which minimal gauged F-term hybrid inflation is coupled renormalisably to the minimal supersymmetric standard model (MSSM), with no extra ingredients; we call this class the "minimal hybrid inflationary supersymmetric standard model" (MHISSM). The singlet inflaton supplies the Higgs mu-term, and allows an exit from inflation to a vacuum characterised by large Higgs vevs. The true ground state is reached after an period of thermal inflation along the Higgs flat direction. The scalar spectral index is reduced from the standard F-term value, to approximately 0.976 in the case where the inflaton potential is dominated by the 1-loop corrections. The reheat temperature following thermal inflation is about 10^9 GeV, solving the gravitino overclosure problem. A Higgs condensate reduces the cosmic string mass per unit length, rendering it compatible with the Cosmic Microwave Background constraints without tuning the inflaton coupling. With the minimal U(1)' gauge symmetry in the inflaton sector, where one of the waterfall fields generates a right-handed neutrino mass, we investigate the Higgs thermal inflation scenario in three popular supersymmetry-breaking schemes: AMSB, GMSB and the CMSSM, focusing on the constraints on the gravitino mass.

We derive the geodesic equation of motion in the presence of weak gravitational fields produced by relativistic sources such as cosmic strings, decomposed into scalar, vector and tensor parts. We find that the vector (gravito-magnetic) force is an important contributor, and for non-relativistic particles we recover the well-known result for the impulse from a moving straight string. Our results can be straightforwardly incorporated into N-body simulations to allow for the presence of cosmic defects or other sources of weak gravitational fields.

We explore the cosmological dynamics of an effective f(R) model constructed from a renormalisation group (RG) improvement of the Einstein--Hilbert action, using the non-perturbative beta functions of the exact renormalisation group equation. The resulting f(R) model has some remarkable properties. It naturally exhibits an unstable de Sitter era in the ultraviolet (UV), dynamically connected to a stable de Sitter era in the IR, via a period of radiation and matter domination, thereby describing a non-singular universe. We find that the UV de Sitter point is one of an infinite set, which make the UV RG fixed point inaccessible to classical cosmological evolution. In the vicinity of the fixed point, the model behaves as R^2 gravity, while it correctly recovers General Relativity at solar system scales. In this simplified model, the fluctuations are too large to be the observed ones, and more ingredients in the action are needed.

We study oscillons, extremely long-lived localized oscillations of a scalar field, with three different potentials: quartic, sine-Gordon model and in a new class of convex potentials. We use an absorbing boundary at the end of the lattice to remove emitted radiation. The energy and the frequency of an oscillon evolve in time and are well fitted by a constant component and a decaying, radiative part obeying a power law as a function time. The power spectra of the emitted radiation show several distinct frequency peaks where oscillons release energy. In two dimensions, and with suitable initial conditions, oscillons do not decay within the range of the simulations, which in quartic theory reach 10^8 time units. While it is known that oscillons in three-dimensional quartic theory and sine-Gordon model decay relatively quickly, we observe a surprising persistence of the oscillons in the convex potential with no sign of demise up to 10^7 time units. This leads us to speculate that an oscillon in such a potential could actually live infinitely long both in two and three dimensions.

We present a significant update of the constraints on the Abelian Higgs cosmic string tension by cosmic microwave background (CMB) data, enabled both by the use of new high-resolution CMB data from suborbital experiments as well as the latest results of the WMAP satellite, and by improved predictions for the impact of Abelian Higgs cosmic strings on the CMB power spectra. The new cosmic string spectra (presented in a previous work) were improved especially for small angular scales, through the use of larger Abelian Higgs string simulations and careful extrapolation. If Abelian Higgs strings are present then we find improved bounds on their contribution to the CMB anisotropies, f10< 0.095, and on their tension, G\mu< 0.57 x 10^-6, both at 95% confidence level using WMAP7 data; and f10 < 0.048 and G\mu < 0.42 x 10^-6 using all the CMB data. We also find that using all the CMB data, a scale invariant initial perturbation spectrum, ns=1, is now disfavoured at 2.4\sigma\ even if strings are present. A Bayesian model selection analysis no longer indicates a preference for strings.

We construct complete sets of (open and closed string) covariant coherent state and mass eigenstate vertex operators in bosonic string theory. By minimally extending the standard definition of coherent states so as to include the string theory requirements, we show that the naive construction of the the closed string coherent states requires the existence of a lightlike compactification of spacetime. When the null winding states in the underlying Hilbert space are projected out the resulting vertex operators satisfy the definition of a coherent state and have a classical interpretation. We present explicitly both the covariant and lightcone gauge realization of the resulting states using the DDF map that relates the two. We also identify the corresponding general lightcone gauge classical solutions around which the quantum states are fluctuating. We go on to show that both the covariant gauge coherent vertex operators, the corresponding lightcone gauge coherent states and the classical solutions all share the same mass and angular momenta and conjecture that the covariant and lightcone gauge states are different manifestations of the same state and share identical interactions. This construction can be used to study the evolution of fundamental cosmic strings as predicted by string theory and may also be useful for other applications where massive string vertex operators are of interest.

We study quantum modifications to cosmology in a Friedmann-Robertson-Walker universe with and without scalar fields by taking the renormalisation group running of gravitational and matter couplings into account. We exploit the Bianchi identity to relate the renormalisation group scale with scale factor and derive the improved cosmological evolution equations. We find two types of cosmological fixed points where the renormalisation group scale either freezes in, or continues to evolve with scale factor. We discuss the implications of each of these, and classify the different cosmological fixed points with and without gravity displaying an asymptotically safe renormalisation group fixed point. We state conditions of existence for an inflating ultraviolet cosmological fixed point for Einstein gravity coupled to a scalar field. We also discuss other fixed point solutions such as "scaling" solutions, or fixed points with equipartition between kinetic and potential energies.

A class of well-motivated models of inflation end by producing cosmic strings. The current status of efforts to calculate and observe the signals from such models are outlined, with a particular emphasis on cosmic strings, and on the Cosmic Microwave Background signal.

COrE (Cosmic Origins Explorer) is a fourth-generation full-sky, microwave-band satellite recently proposed to ESA within Cosmic Vision 2015-2025. COrE will provide maps of the microwave sky in polarization and temperature in 15 frequency bands, ranging from 45 GHz to 795 GHz, with an angular resolution ranging from 23 arcmin (45 GHz) and 1.3 arcmin (795 GHz) and sensitivities roughly 10 to 30 times better than PLANCK (depending on the frequency channel). The COrE mission will lead to breakthrough science in a wide range of areas, ranging from primordial cosmology to galactic and extragalactic science. COrE is designed to detect the primordial gravitational waves generated during the epoch of cosmic inflation at more than $3\sigma $ for $r=(T/S)>=10^{-3}$. It will also measure the CMB gravitational lensing deflection power spectrum to the cosmic variance limit on all linear scales, allowing us to probe absolute neutrino masses better than laboratory experiments and down to plausible values suggested by the neutrino oscillation data. COrE will also search for primordial non-Gaussianity with significant improvements over Planck in its ability to constrain the shape (and amplitude) of non-Gaussianity. In the areas of galactic and extragalactic science, in its highest frequency channels COrE will provide maps of the galactic polarized dust emission allowing us to map the galactic magnetic field in areas of diffuse emission not otherwise accessible to probe the initial conditions for star formation. COrE will also map the galactic synchrotron emission thirty times better than PLANCK. This White Paper reviews the COrE science program, our simulations on foreground subtraction, and the proposed instrumental configuration.

The proposed CMBPol mission will be able to detect the imprint of topological defects on the cosmic microwave background (CMB) provided the contribution is sufficiently strong. We quantify the detection threshold for cosmic strings and for textures, and analyse the satellite's ability to distinguish between these different types of defects. We also assess the level of danger of misidentification of a defect signature as from the wrong defect type or as an effect of primordial gravitational waves. A 0.002 fractional contribution of cosmic strings to the CMB temperature spectrum at multipole ten, and similarly a 0.001 fractional contribution of textures, can be detected and correctly identified at the 3{\sigma} level. We also confirm that a tensor contribution of r = 0.0018 can be detected at over 3{\sigma}, in agreement with the CMBpol mission concept study. These results are supported by a model selection analysis.

We present a significant improvement over our previous calculations of the cosmic string contribution to cosmic microwave background (CMB) power spectra, with particular focus on sub-WMAP angular scales. These smaller scales are relevant for the now-operational Planck satellite and additional sub-orbital CMB projects that have even finer resolutions. We employ larger Abelian Higgs string simulations than before and we additionally model and extrapolate the statistical measures from our simulations to smaller length scales. We then use an efficient means of including the extrapolations into our Einstein-Boltzmann calculations in order to yield accurate results over the multipole range 2 < l < 4000. Our results suggest that power-law behaviour cuts in for l > 3000 in the case of the temperature power spectrum, which then allows cautious extrapolation to even smaller scales. We find that a string contribution to the temperature power spectrum making up 10% of power at l=10 would be larger than the Silk-damped primary adiabatic contribution for l > 3500. Astrophysical contributions such as the Sunyaev-Zeldovich effect also become important at these scales and will reduce the sensitivity to strings, but these are potentially distinguishable by their frequency-dependence.

We provide an analytical expression for the trispectrum of the Cosmic Microwave Background (CMB) temperature anisotropies induced by cosmic strings. Our result is derived for the small angular scales under the assumption that the temperature anisotropy is induced by the Gott-Kaiser-Stebbins effect. The trispectrum is predicted to decay with a non-integer power-law exponent l^(-r) with 6<r<7, depending on the string microstructure, and thus on the string model. For Nambu-Goto strings, this exponent is related to the string mean square velocity and the loop distribution function. We then explore two classes of wavenumber configuration in Fourier space, the kite and trapezium quadrilaterals. The trispectrum can be of any sign and appears to be strongly enhanced for all squeezed quadrilaterals.

The Cosmic Microwave Background (CMB) bispectrum of the temperature anisotropies induced by a network of cosmic strings is derived for small angular scales, under the assumption that the principal cause of temperature fluctuations is the Gott-Kaiser-Stebbins (GKS) effect. We provide analytical expressions for all isosceles triangle configurations in Fourier space. Their overall amplitude is amplified as the inverse cube of the angle and diverges for flat triangles. The isosceles configurations generically lead to a negative bispectrum with a power law decay l^(-6) for large multipole l. However, collapsed triangles are found to be associated with a positive bispectrum whereas the squeezed triangles still exhibit negative values. We then compare our analytical estimates to a direct computation of the bispectrum from a set of 300 statistically independent temperature maps obtained from Nambu-Goto cosmic string simulations in a Friedmann-Lemaitre-Robertson-Walker (FLRW) universe. We find good agreement for the overall amplitude, the power law behaviour and angle dependency of the various triangle configurations. At l~500 the cosmic string GKS effect contributes approximately the same equilateral CMB bispectrum amplitude as an inflationary model with |fNL|~10^3, if the strings contribute about 10% of the temperature power spectrum at l=10. Current bounds on fNL are not derived using cosmic string bispectrum templates, and so our fNL estimate cannot be used to derive bounds on strings. However it does suggest that string bispectrum templates should be included in the search of CMB non-Gaussianities.

We present the first ever calculation of cosmic microwave background CMB anisotropy power spectra from semilocal cosmic strings, obtained via simulations of a classical field theory. Semilocal strings are a type of non-topological defect arising in some models of inflation motivated by fundamental physics, and are thought to relax the constraints on the symmetry breaking scale as compared to models with (topological) cosmic strings. We derive constraints on the model parameters, including the string tension parameter mu, from fits to cosmological data, and find that in this regard BPS semilocal strings resemble global textures more than topological strings. The observed microwave anisotropy at l = 10 is reproduced if Gmu = 5.3x10^{-6} (G is Newton's constant). However as with other defects the spectral shape does not match observations, and in models with inflationary perturbations plus semilocal strings the 95% confidence level upper bound is Gmu<2.0x10^{-6} when CMB data, Hubble Key Project and Big Bang Nucleosynthesis data are used (c.f. Gmu<0.9x10^{-6} for cosmic strings). We additionally carry out a Bayesian model comparison of several models with and without defects, showing models with defects are neither conclusively favoured nor disfavoured at present.

While observations indicate that the predominant source of cosmic inhomogeneities are adiabatic perturbations, there are a variety of candidates to provide auxiliary trace effects, including inflation-generated primordial tensors and cosmic defects which both produce B-mode cosmic microwave background (CMB) polarization. We investigate whether future experiments may suffer confusion as to the true origin of such effects, focusing on the ability of Planck to distinguish tensors from cosmic strings, and show that there is no significant degeneracy.

We perform a multiparameter likelihood analysis to compare measurements of the cosmic microwave background (CMB) power spectra with predictions from models involving cosmic strings. Adding strings to the standard case of a primordial spectrum with power-law tilt n, we find a 2-sigma detection of strings: f_10 = 0.11 +/- 0.05, where f_10 is the fractional contribution made by strings in the temperature power spectrum (at multipole l = 10). CMB data give moderate preference to the model n = 1 with cosmic strings over the standard zero-strings model with variable tilt. When additional non-CMB data are incorporated, the two models become on a par. With variable n and these extra data, we find that f_10 < 0.11, which corresponds to G mu < 0.7x10^-6 (where mu is the string tension and G is the gravitational constant).

We present the first calculation of the possible (local) cosmic string contribution to the cosmic microwave background polarization spectra from simulations of a string network (rather than a stochastic collection of unconnected string segments). We use field theory simulations of the Abelian Higgs model to represent local U(1) strings, including their radiative decay and microphysics. Relative to previous estimates, our calculations show a shift in power to larger angular scales, making the chance of a future cosmic string detection from the B-mode polarization slightly greater. We explore a future ground-based polarization detector, taking the CLOVER project as our example. In the null hypothesis (that cosmic strings make a zero contribution) we find that CLOVER should limit the string tension \mu to G\mu<0.12e-6 (where G is the gravitational constant), above which it is likely that a detection would be possible.

A class of Modified Gravity Models, consisting of inverse powers of linear combination of quadratic curvature invariants, is studied in the full parameter space. We find that singularity-free cosmological solutions, interpolating between an almost-Friedmann universe at Big Bang Nucleosynthesis and an accelerating universe today, exist only in a restricted parameter space. Furthermore, for all parameters of the models, there is an unstable scalar mode of the gravitational field. Therefore we conclude that this class of Modified Gravity Models is not viable.

We present the first field-theoretic calculations of the contribution made by cosmic strings to the temperature power spectrum of the cosmic microwave background (CMB). Unlike previous work, in which strings were modeled as idealized one-dimensional objects, we evolve the simplest example of an underlying field theory containing local U(1) strings, the Abelian Higgs model. Limitations imposed by finite computational volumes are overcome using the scaling property of string networks and a further extrapolation related to the lessening of the string width in comoving coordinates. The strings and their decay products, which are automatically included in the field theory approach, source metric perturbations via their energy-momentum tensor, the unequal-time correlation functions of which are used as input into the CMB calculation phase. These calculations involve the use of a modified version of CMBEASY, with results provided over the full range of relevant scales. We find that the string tension $\mu$ required to normalize to the WMAP 3-year data at multipole $\ell = 10$ is $G\mu = [2.04\pm0.06\textrm{(stat.)}\pm0.12\textrm{(sys.)}] \times 10^{-6}$, where we have quoted statistical and systematic errors separately, and $G$ is Newton's constant. This is a factor 2-3 higher than values in current circulation.

We consider Modified Gravity models involving inverse powers of fourth-order curvature invariants. Using these models' equivalence to the theory of a scalar field coupled to a linear combination of the invariants, we investigate the properties of the propagating modes. Even in the case for which the fourth derivative terms in the field equations vanish, we find that the second derivative terms can give rise to ghosts, instabilities, and superluminal propagation speeds. We establish the conditions which the theories must satisfy in order to avoid these problems in Friedmann backgrounds, and show that the late-time attractor solutions are generically afflicted by superluminally propagating tensor or scalar modes.

We study the evolution of growth and decay laws for the magnetic field coherence length xi, energy E_M and magnetic helicity H in freely decaying 3D MHD turbulence. We show that with certain assumptions, self-similarity of the magnetic power spectrum alone implies that xi~t^{1/2}. This in turn implies that magnetic helicity decays as H~t^{-2s}, where s=(xi_{diff}/xi_H)^2, in terms of xi_{diff}, the diffusion length scale, and xi_H, a length scale defined from the helicity power spectrum. The relative magnetic helicity remains constant, implying that the magnetic energy decays as E_M~t^{-1/2-2s}. The parameter s is inversely proportional to the magnetic Reynolds number Re_M, which is constant in the self-similar regime.

Weakly Interacting Massive Particles (WIMPs) of mass m freeze out at a temperature T_f ~ m/25, i.e. in the range 400 MeV -- 40 GeV for a particle in the typical mass range 10 -- 1000 GeV. The WIMP relic density, which depends on the effective number of relativistic degrees of freedom at T_f, may be measured to better than 1% by Planck, warranting comparable theoretical precision. Recent theoretical and experimental advances in the understanding of high temperature QCD show that the quark gluon plasma departs significantly from ideal behaviour up to temperatures of several GeV, necessitating an improvement of the cosmological equation of state over those currently used. We discuss how this increases the relic density by approximately 1.5 -- 3.5% in benchmark mSUGRA models, with an uncertainly in the QCD corrections of 0.5 -- 1 %. We point out what further work is required to achieve a theoretical accuracy comparable with the expected observational precision, and speculate that the effective number of degrees of freedom at T_f may become measurable in the foreseeable future.

We use the cosmic microwave background angular power spectra to place upper limits on the degree to which global defects may have aided cosmic structure formation. We explore this under the inflationary paradigm, but with the addition of textures resulting from the breaking of a global O(4) symmetry during the early stages of the Universe. As a measure of their contribution, we use the fraction of the temperature power spectrum that is attributed to the defects at a multipole of 10. However, we find a parameter degeneracy enabling a fit to the first-year WMAP data to be made even with a significant defect fraction. This degeneracy involves the baryon fraction and the Hubble constant, plus the normalization and tilt of the primordial power spectrum. Hence, constraints on these cosmological parameters are weakened. Combining the WMAP data with a constraint on the physical baryon fraction from big bang nucleosynthesis calculations and high-redshift deuterium abundance, limits the extent of the degeneracy and gives an upper bound on the defect fraction of 0.13 (95% confidence).

We study the evolution of a stochastic helical magnetic field generated in the early Universe after the electroweak phase transition, using standard magnetohydrodynamics (MHD). We find how the coherence length xi, magnetic energy E_M and magnetic helicity H evolve with time. We show that the self-similarity of the magnetic power spectrum alone implies that xi ~ t^{1/2}. This in turn implies that magnetic helicity decays as H ~ t^{-2s}, and that the magnetic energy decays as E_M ~ t^{-0.5-2s}, where s is inversely proportional to the magnetic Reynolds number Re_M. These laws improve on several previous estimates.

We perform direct numerical simulations of three-dimensional freely decaying magnetohydrodynamic (MHD) turbulence. For helical magnetic fields an inverse cascade effect is observed in which power is transfered from smaller scales to larger scales. The magnetic field reaches a scaling regime with self-similar evolution, and power law behavior at high wavenumbers. We also find power law decay in the magnetic and kinematic energies, and power law growth in the characteristic length scale of the magnetic field.

We have carried out numerical simulations of freely decaying magnetohydrodynamic (MHD) turbulence in three dimensions, which can be applied to the evolution of stochastic magnetic fields in the early Universe. For helical magnetic fields an inverse cascade effect is observed in which magnetic helicity and energy is transfered from smaller scales to larger scales, accompanied by power law growth in the characteristic length scale of the magnetic field. The magnetic field quickly reaches a scaling regime with self-similar evolution, and power law behaviour at high wavenumbers. We also find power law decay in the magnetic and kinematic energies.

There is a reformulation of magnetohydrodynamics in which the fundamental dynamical quantities are the positions and velocities of the lines of magnetic flux in the plasma, which turn out to obey equations of motion very much like ideal strings. We use this approach to study the evolution of a primordial magnetic field generated during the radiation-dominated era in the early Universe. Causality dictates that the field lines form a tangled random network, and the string-like equations of motion, plus the assumption of perfect reconnection, inevitably lead to a self-similar solution for the magnetic field power spectrum. We present the predicted form of the power spectrum, and discuss insights gained from the string approximation, in particular the implications for the existence or not of an inverse cascade.

Recent developments in inflation model building, based on supersymmetry, have produced compelling models in which strings are produced at the end of inflation. In such models the cosmological perturbations are seeded both by the defects and by the quantum fluctuations. We show that such models produce qualitatively new and desirable predictions for CMB anisotropies and the CDM power spectrum. This remark should put an end to the long term animosity between defect and inflationary scenarios of structure formation.

We evaluate the two point functions of the stress energy from the largest string simulations carried out so far. The two point functions are used to calculate the cosmic microwave background (CMB) and cold dark matter (CDM) power spectra from local cosmic string models for structure formation. We find that our spectra differ significantly from those previously calculated for both global and local defects. We find a higher Doppler peak at $l=400-600$ and a less severe bias problem than for global defects. Spectra were obtained for a variety of network energy-decay mechanisms.

In cosmological models where local cosmic strings are formed at the end of a period of inflation, the perturbations are seeded both by the defects and by the quantum fluctuations. In a subset of these models, for example those based on $D$-term inflation, the amplitudes are similar. Using our recent calculations of structure formation with cosmic strings, we point out that in a flat cosmology with zero cosmological constant and 5% baryonic component, strings plus inflation fits the observational data much better than each component individually. The large-angle CMB spectrum is mildly tilted, for Harrison-Zeldovich inflationary fluctuations. It then rises to a thick Doppler bump, covering $\ell=200-600$, modulated by soft secondary undulations. The standard CDM anti-biasing problem is cured, giving place to a slightly biased scenario of galaxy formation.

We compute the power spectra in the cosmic microwave background and cold dark matter (CDM) fluctuations seeded by strings, using the largest string simulations performed so far to evaluate the two-point functions of their stress energy tensor. We find that local strings differ from global defects in that the scalar components of the stress-energy tensor dominate over vector and tensor components. This result has far reaching consequences. We find that cosmic strings exhibit a single Doppler peak of acceptable height at high $\ell$. They also seem to have a less severe bias problem than global defects, although the CDM power spectrum in the ``standard'' cosmology (flat geometry, zero cosmological constant, 5% baryonic component) is the wrong shape to fit large scale structure data.

I report on recent numerical simulations of the simplest field theory with cosmic string solutions, the Abelian Higgs model. We find that random networks of string quickly converge to a scaling solution in which the network scale length $\xi$ increases linearly with time. There are very few loops with sizes less than $\xi$, and the strings are smooth, showing no signs of ``small scale structure''. We claim that particle production is the dominant energy-loss mechanism, not gravitational radiation as previously thought. For strings in Grand Unified Models, stringent constraints can be placed from cosmic ray observations on the string tension $\mu$: we estimate $G\mu < 10^{-9}$, three orders of magnitude lower than the constraint from Cosmic Microwave Background fluctuations.

We present the results of a field theory simulation of networks of strings in the Abelian Higgs model. Starting from a random initial configuration we show that the resulting vortex tangle approaches a self-similar regime in which the length density of lines of zeros of $\phi$ reduces as $t^{-2}$. We demonstrate that the network loses energy directly into scalar and gauge radiation. These results support a recent claim that particle production, and not gravitational radiation, is the dominant energy loss mechanism for cosmic strings. This means that cosmic strings in Grand Unified Theories are severely constrained by high energy cosmic ray fluxes: either they are ruled out, or an implausibly small fraction of their energy ends up in quarks and leptons.

A particular class of variant axion models with two higgs doublets and a singlet is studied. In these models the axion couples either to the $u$-quark or $t$-quark or both, but not to $b$, $c$, $s$, or $d$. When the axion couples to only one quark the models possess the desirable feature of having no domain wall problem, which makes them viable candidates for a cosmological axion string scenario. We calculate the axion couplings to leptons, photons and nucleons, and the astrophysical constraints on the axion decay constant $v_a$ are investigated and compared to the DFSZ axion model. We find that the most restrictive lower bound on $v_a$, that from SN1987a, is lowered by up to a factor of about 30, depending on the model and also the ratio of the vacuum expectation values of the higgs doublets. For scenarios with axionic strings, the allowed window for $v_a$ in the $u$ quark model can be more than two orders of magnitude. For inflationary scenarios, the cosmological upper bound on $v_a/N$, where $N$ is the QCD anomaly factor, is unaffected: however, the variant models have $N$ either 3 or 6 times smaller than the DFSZ model.

The generation of primordial magnetic fields from cosmological phase transitions is discussed, paying particular attention to the electroweak transition and to the various definitions of the `average' field that have been put forward. It is emphasised that only the volume average has dynamical significance as a seed for galactic dynamos. On rather general grounds of causality and energy conservation, it is shown that, in the absence of MHD effects that transfer power in the magnetic field from small to large scales, processes occurring at the electroweak transition cannot generate fields stronger than $10^{-20}$ Gauss on a scale of 0.5 Mpc. However, it is implausible that this upper bound could ever be reached, as it would require all the energy in the Universe to be turned into a magnetic field coherent at the horizon scale. Non-linear MHD effects seem therefore to be necessary if the electroweak transition is to create a primordial seed field.

We present a full characterization of the phase transition in U(1) scalar field theory and of the associated vortex string thermodynamics in 3D. We show that phase transitions in the string densities exist and measure their critical exponents, both for the long string and the short loops. Evidence for a natural separation between these two string populations is presented. In particular our results strongly indicate that an infinite string population will only exist above the critical temperature. Canonical initial conditions for cosmic string evolution are show to correspond to the infinite temperature limit of the theory.

The one-texture universe, introduced by Davis in 1987, is a homogeneous mapping of a scalar field with an $S^3$ vacuum into a closed universe. It has long been known to mathematicians that such solutions, although static, are unstable. We show by explicit construction that there is only one unstable mode, corresponding to collapse of the texture towards a single point, in the case where gravitational backreaction is neglected. We discuss the instability timescale in both static and expanding space-times; in the latter case it is of order of the present age of the universe, suggesting that, though unstable, the one-texture universe could survive to the present. The cosmic microwave background constrains the initial magnitude of this unstable perturbation to be less than of order 10^{-3}.

We examine the scaling properties of an evolving network of strings in Minkowski spacetime and study the evolution of length scales in terms of a 3-scale model proposed by Austin, Copeland and Kibble (ACK). We find good qualitative and some quantitative agreement between the model and our simulations. We also investigate small-scale structure by altering the minimum allowed size for loop production E_c. Certain quantities depend significantly on this parameter: for example the scaling density can vary by a factor of two or more with increasing E_c. Small-scale structure as defined by ACK disappears if no restrictions are placed on loop production, and the fractal dimension of the string changes smoothly from 2 to 1 as the resolution scale is decreased. Loops are nearly all produced at the lattice cut-off. We suggest that the lattice cut-off should be interpreted as corresponding to the string width, and that in a real network loops are actually produced with this size. This leads to a radically different string scenario, with particle production rather than gravitational radiation being the dominant mode of energy dissipation. At the very least, a better understanding of the discretisation effects in all simulations of cosmic strings is called for.

Cosmic strings provide a radically different paradigm for the formation of structure to the prevailing inflationary one. They afford some extra technical complications: for example, the calculation of the power spectrum of matter and radiation perturbations requires the knowledge of the history of the evolution of the defects in the form of two-time correlation functions. We describe some numerical simulations of string networks, designed to measure the two-time correlations during their evolution.

We investigate scaling and correlations of the energy and momentum in an evolving network of cosmic strings in Minkowski space. These quantities are of great interest, as they must be understood before accurate predictions for the power spectra of the perturbations in the matter and radiation in the early Universe can be made. We argue that Minkowski space provides a reasonable approximation to a Friedmann background for string dynamics and we use our results to construct a simple model of the network, in which it is considered to consist of randomly placed segments moving with random velocities. This model works well in accounting for features of the two-time correlation functions, and even better for the power spectra.

A relativistic generalisation of a well-known method for approximating the dynamics of topological defects in condensed matter is constructed, and applied to the evolution of domain walls in a cosmological context. It is shown that there are self-similar ``scaling'' solutions, for which one can in principle calculate many quantities of interest without recourse to numerical simulations. Here, the area density in the scaling regime is calculated in various backgrounds. Remarkably good agreement with numerical simulations is obtained.

Some aspects of the theory of cosmological perturbations from cosmic strings and other topological defects are outlined, with particular reference to a simple example: a spatially flat CDM-dominated universe. The conserved energy-momentum pseudo-tensor is introduced, and the equation for the density perturbation derived from it. It is shown how the scaling hypothesis for defect evolution results in a Harrison-Zel'dovich spectrum for wavelengths well inside the horizon.

The topic of cosmic strings provides a bridge between the physics of the very small and the very large. They are predicted by some unified theories of particle interactions. If they exist, they may help to explain some of the largest-scale structures seen in the Universe today. They are `topological defects' that may have been formed at phase transitions in the very early history of the Universe, analogous to those found in some condensed-matter systems --- vortex lines in liquid helium, flux tubes in type-II superconductors, or disclination lines in liquid crystals. In this review, we describe what they are, why they have been hypothesized and what their cosmological implications would be. The relevant background from the standard models of particle physics and cosmology is described in section 1. In section 2, we review the idea of symmetry breaking in field theories, and show how the defects formed are constrained by the topology of the manifold of degenerate vacuum states. We also discuss the different types of cosmic strings that can appear in different field theories. Section 3 is devoted to the dynamics of cosmic strings, and section 4 to their interaction with other fields. The formation and evolution of cosmic strings in the early Universe is the subject of section 5, while section 6 deals with their observational implications. Finally, the present status of the theory is reviewed in section 7.

We investigate numerically the configurational statistics of strings. The algorithm models an ensemble of global $U(1)$ cosmic strings, or equivalently vortices in superfluid $^4$He. We use a new method which avoids the specification of boundary conditions on the lattice. We therefore do not have the artificial distinction between short and long string loops or a `second phase' in the string network statistics associated with strings winding around a toroidal lattice. Our lattice is also tetrahedral, which avoids ambiguities associated with the cubic lattices of previous work. We find that the percentage of infinite string is somewhat lower than on cubic lattices, 63\% instead of 80\%. We also investigate the Hagedorn transition, at which infinite strings percolate, controlling the string density by rendering one of the equilibrium states more probable. We measure the percolation threshold, the critical exponent associated with the divergence of a suitably defined susceptibility of the string loops, and that associated with the divergence of the correlation length.

The Cosmic Microwave Background (CMB) fluctuations at very small angular scales (less than $10'$) induced by matter sources are computed in a simplified way. The result corrects a previous formula appearing in the literature. The small scale power spectrum from cosmic strings is then calculated by a new analytic method. The result compares extremely well with the spectrum computed by numerical techniques (when the old, incorrect, formula is used). The upper bound on the string parameters derived from OVRO data is re-examined, taking into account the non-Gaussian nature of stringy perturbations on small scales. Assuming a conventional ionization history, the bound is $\gamma G\mu < 11\times 10^{-6}$, where $\gamma^2$ is the number of horizon lengths of string per horizon volume. Current simulations give $\gamma^2 = 31\pm 7$.