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

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 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 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 perform large-scale field theoretical simulations in expanding universe to characterize a network of strings that can form composed bound states. The network consists of two copies of Abelian Higgs strings (which we label p and q, respectively) coupled via a potential term to give pq bound states. The simulations are performed using two different kinds of initial conditions: the first one with a network of p- and q-strings, and the second one with a network of q- and pq-strings. This way, we start from two opposite situations: one with no initial pq-strings, and one with a large initial number of pq-strings. We find that in both cases the system scales, and in both cases the system prefers to have a low fraction of pq-strings. This is somewhat surprising in the case for the second type of conditions, showing that the unzipping mechanism is very efficient. We also find hints that both initial conditions tend to asymptote to a common configuration, though we would need a larger dynamical range to confirm it. The average velocities of the different types of strings in the network have also been explored for the first time.

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

The recent observation that the Cosmic Microwave Background (CMB) may prefer a neutrino excess has triggered a number of works studying this possibility. The effect obtained by the non-interacting massless neutrino excess could be mimicked by some extra radiation component in the early universe, such as a cosmological gravitational wave background. Prompted by the fact that a possible candidate to source those gravitational waves would be cosmic strings, we perform a parameter fitting study with models which considers both cosmic strings and the effective number of neutrinos as free parameters, using CMB and non-CMB data. We find that there is a correlation between cosmic strings and the number of extra relativistic species, and that strings can account for all the extra radiation necessary. In fact, CMB data prefer strings at a 2sigma level, paying the price of a higher extra radiation component. CMB data also give a moderate preference for a model with ns=1. The inclusion of non-CMB data lowers both the preference for strings and for the extra relativistic species.