Abstract: 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(GW) proportional to f(3) for f << f(0), Omega(GW) proportional to 1/f(2) for f(0) less than or similar to f less than or similar to feq, and Omega(GW) proportional to const (i.e., exactly scale invariant) for f >> f(eq), where f(0) and f(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) -> 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 cosmic microwave background upper bounds on the defects' energy scale, we discuss the difficulty of detecting this GW background in the case of global defects.

Abstract: We present a comprehensive discussion on lattice techniques for the simulation of scalar and gauge field dynamics in an expanding universe. After reviewing the continuum formulation of scalar and gauge field interactions in Minkowski and FLRW backgrounds, we introduce the basic tools for the discretization of field theories, including lattice gauge invariant techniques. Following, we discuss and classify numerical algorithms, ranging from methods of O(delta t(2)) accuracy like staggered leapfrog and Verlet integration, to Runge-Kutta methods up to O(delta t(4)) accuracy, and the Yoshida and Gauss-Legendre higher-order integrators, accurate up to O(delta t(10)) We adapt these methods for their use in classical lattice simulations of the non-linear dynamics of scalar and gauge fields in an expanding grid in 3+1 dimensions, including the case of 'self-consistent' expansion sourced by the volume average of the fields' energy and pressure densities. We present lattice formulations of canonical cases of: i) Interacting scalar fields, ii) Abelian U(1) gauge theories, and iii) Non-Abelian SU(2) gauge theories. In all three cases we provide symplectic integrators, with accuracy ranging from O(delta t(2)) up to O(delta t(10)) For each algorithm we provide the form of relevant observables, such as energy density components, field spectra and the Hubble constraint. We note that all our algorithms for gauge theories always respect the Gauss constraint to machine precision, including when 'self-consistent' expansion is considered. As a numerical example we analyze the post-inflationary dynamics of an oscillating inflaton charged under SU(2) x U(1). We note that the present manuscript is meant to be part of the theoretical basis for the code CosmoLattice, a multi-purpose MPI-based package for simulating the non-linear evolution of field theories in an expanding universe, publicly available at http://www.cosrnolattice.net.