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

Semilocal strings -- a particular limit of electroweak strings -- are an interesting example of a stable non-topological defect whose properties resemble those of their topological cousins, the Abrikosov-Nielsen-Olesen vortices. There is, however, one important difference: a network of semilocal strings will contain segments. These are 'dumbbells' whose ends behave almost like global monopoles that are strongly attracted to one another. While closed loops of string will eventually shrink and disappear, the segments can either shrink or grow, and a cosmological network of semilocal strings will reach a scaling regime. We discuss attempts to find a "thermodynamic" description of the cosmological evolution and scaling of a network of semilocal strings, by analogy with well-known descriptions for cosmic strings and for monopoles. We propose a model for the time evolution of an overall lengthscale and typical velocity for the network as well as for its segments, and some supporting (preliminary) numerical evidence.

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

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

We continue a comprehensive numerical study of semilocal string networks and their cosmological evolution. These can be thought of as hybrid networks comprised of (non-topological) string segments, whose core structure is similar to that of Abelian Higgs vortices, and whose ends have long--range interactions and behaviour similar to that of global monopoles. Our study provides further evidence of a linear scaling regime, already reported in previous studies, for the typical length scale and velocity of the network. We introduce a new algorithm to identify the position of the segment cores. This allows us to determine the length and velocity of each individual segment and follow their evolution in time. We study the statistical distribution of segment lengths and velocities for radiation- and matter-dominated evolution in the regime where the strings are stable. Our segment detection algorithm gives higher length values than previous studies based on indirect detection methods. The statistical distribution shows no evidence of (anti)correlation between the speed and the length of the segments.

We present an estimation of the average velocity of a network of global monopoles in a cosmological setting using large numerical simulations. In order to obtain the value of the velocity, we improve some already known methods, and present a new one. This new method estimates individual global monopole velocities in a network, by means of detecting each monopole position in the lattice and following the path described by each one of them. Using our new estimate we can settle an open question previously posed in the literature: velocity-dependent one-scale (VOS) models for global monopoles predict two branches of scaling solutions, one with monopoles moving at subluminal speeds and one with monopoles moving at luminal speeds. Previous attempts to estimate monopole velocities had large uncertainties and were not able to settle that question. Our simulations find no evidence of a luminal branch. We also estimate the values of the parameters of the VOS model. With our new method we can also study the microphysics of the complicated dynamics of individual monopoles. Finally we use our large simulation volume to compare the results from the different estimator methods, as well as to asses the validity of the numerical approximations made.

In this work we calibrate two different analytic models of semilocal strings by constraining the values of their free parameters. In order to do so, we use data obtained from the largest and most accurate field theory simulations of semilocal strings to date, and compare several key properties with the predictions of the models. As this is still work in progress, we present some preliminary results together with descriptions of the methodology we are using in the characterisation of semilocal string networks.

We study field theoretical models for cosmic strings with flat directions in curved space-time. More precisely, we consider minimal models with semilocal, axionic and tachyonic strings, respectively. In flat space-time, the string solutions of these models have a flat direction, i.e., a uniparametric family of configurations with the same energy exists which is associated to a zero mode. We prove that the zero mode survives coupling to gravity, and study the role of the flat direction when coupling the strings to gravity. Even though the total energy of the solution is the same, and thus the global properties of the family of solutions remains unchanged, the energy density, and therefore the gravitational properties, are different. The local structure of the solutions depends strongly on the value of the parameter describing the flat direction; for example, for supermassive strings, the value of the free parameter can determine the size of the universe.