We derive a full analytical approximation to the stochastic gravitational wave background generated by the loops that are produced throughout the cosmological evolution of cosmic string networks. We show that this approximation not only predicts the amplitude of the radiation-era plateau exactly, but also provides a good fit to the high-frequency cut-off and to the low-frequency peak generated by the loops that decay during the matter era, irrespective of cosmic string tension and of the length of loops created. We then find that it provides a good quantitative description of the full stochastic gravitational wave background across the relevant frequency range.

Cosmic string networks offer one of the best prospects for detection of cosmological gravitational waves (GWs). The combined incoherent GW emission of a large number of string loops leads to a stochastic GW background (SGWB), which encodes the properties of the string network. In this paper we analyze the ability of the Laser Interferometer Space Antenna (LISA) to measure this background, considering leading models of the string networks. We find that LISA will be able to probe cosmic strings with tensions $G\mu \gtrsim \mathcal{O}(10^{-17})$, improving by about $6$ orders of magnitude current pulsar timing arrays (PTA) constraints, and potentially $3$ orders of magnitude with respect to expected constraints from next generation PTA observatories. We include in our analysis possible modifications of the SGWB spectrum due to different hypotheses regarding cosmic history and the underlying physics of the string network. These include possible modifications in the SGWB spectrum due to changes in the number of relativistic degrees of freedom in the early Universe, the presence of a non-standard equation of state before the onset of radiation domination, or changes to the network dynamics due to a string inter-commutation probability less than unity. In the event of a detection, LISA's frequency band is well-positioned to probe such cosmic events. Our results constitute a thorough exploration of the cosmic string science that will be accessible to LISA.

A cosmic string network created during an inflationary stage in the early Universe - here defined as i-string network - is expected to enter a transient stretching regime during inflation, in which its characteristic length is stretched to scales much larger than the Hubble radius, before attaining a standard evolution once the network re-enters the Hubble volume after inflation. During the stretching regime, the production of cosmic string loops and consequent emission of gravitational radiation are significantly suppressed. Here, we compute the power spectrum of the stochastic gravitational wave background generated by i-string networks using the velocity-dependent one scale model to describe the network dynamics, and demonstrate that this regime introduces a high-frequency signature on an otherwise standard spectrum of the stochastic gravitational wave background generated by cosmic strings. We argue that, if observed by current or forthcoming experiments, this signature would provide strong evidence for i-strings and, therefore, for (primordial) inflation. We also develop a simple single-parameter algorithm for the computation of the stochastic gravitational wave background generated by i-strings from that of a standard cosmic string network, which may be useful in the determination of the observational constraints to be obtained by current and forthcoming gravitational wave experiments.

We consider a model where particles are described as localized concentrations of energy, with fixed rest mass and structure, which are not significantly affected by their self-induced gravitational field. We show that the volume average of the on-shell matter Lagrangian ${\mathcal L_m}$ describing such particles, in the proper frame, is equal to the volume average of the trace $T$ of the energy-momentum tensor in the same frame, independently of the particle's structure and constitution. Since both ${\mathcal L_m}$ and $T$ are scalars, and thus independent of the reference frame, this result is also applicable to collections of moving particles and, in particular, to those which can be described by a perfect fluid. Our results are expected to be particularly relevant in the case of modified theories of gravity with nonminimal coupling to matter where the matter Lagrangian appears explicitly in the equations of motion of the gravitational and matter fields, such as $f(R,{\mathcal L_m})$ and $f(R,T)$ gravity. In particular, they indicate that, in this context, $f(R,{\mathcal L_m})$ theories may be regarded as a subclass of $f(R,T)$ gravity.

We revisit the physical properties of global and local monopoles and discuss their implications in the dynamics of monopole networks. In particular, we review the Velocity-dependent One-Scale (VOS) model for global and local monopoles and propose physically motivated changes to its equations. We suggest a new form for the acceleration term of the evolution equation of the root-mean-squared velocity and show that, with this change, the VOS model is able to describe the results of radiation and matter era numerical simulations of global monopole networks with a single value of the acceleration parameter $k$, thus resolving the tension previously found in the literature. We also show that the fact that the energy of global monopoles is not localized within their cores affects their dynamics and, thus, the Hubble damping terms in the VOS equations. We study the ultra-relativistic linear scaling regime predicted by the VOS equations and demonstrate that it cannot be attained either on radiation or matter eras and, thus, cannot arise from the cosmological evolution of a global monopole network. We also briefly discuss the implications of our findings for the VOS model for local monopoles.

We consider a model with two real scalar fields which admits phantom domain wall solutions. We investigate the structure and evolution of these phantom domain walls in an expanding homogeneous and isotropic universe. In particular, we show that the increase of the tension of the domain walls with cosmic time, associated to the evolution of the phantom scalar field, is responsible for an additional damping term in their equations of motion. We describe the macroscopic dynamics of phantom domain walls, showing that extended phantom defects whose tension varies on a cosmological timescale cannot be the dark energy.

We explore the dynamics and evolution of the Universe at early and late times, focusing on both dark energy and extended gravity models and their astrophysical and cosmological consequences. Modified theories of gravity not only provide an alternative explanation for the recent expansion history of the universe, but they also offer a paradigm fundamentally distinct from the simplest dark energy models of cosmic acceleration. In this review, we perform a detailed theoretical and phenomenological analysis of different modified gravity models and investigate their consistency. We also consider the cosmological implications of well motivated physical models of the early universe with a particular emphasis on inflation and topological defects. Astrophysical and cosmological tests over a wide range of scales, from the solar system to the observable horizon, severely restrict the allowed models of the Universe. Here, we review several observational probes -- including gravitational lensing, galaxy clusters, cosmic microwave background temperature and polarization, supernova and baryon acoustic oscillations measurements -- and their relevance in constraining our cosmological description of the Universe.

We compute the stochastic gravitational wave background generated by cosmic superstrings using a semi-analytical velocity-dependent model to describe their dynamics. We show that heavier string types may leave distinctive signatures on the stochastic gravitational wave background spectrum within the reach of present and upcoming gravitational wave detectors. We examine the physically motivated scenario in which the physical size of loops is determined by the gravitational backreaction scale and use NANOGRAV data to derive a conservative constraint of $G\mu_F<3.2 \times 10^{-9}$ on the tension of fundamental strings. We demonstrate that approximating the gravitational wave spectrum generated by cosmic superstring networks using the spectrum generated by ordinary cosmic strings with reduced intercommuting probability (which is often done in the literature) leads, in general, to weaker observational constraints on $G\mu_F$. We show that the inclusion of heavier string types is required for a more accurate characterization of the region of the $(g_s,G\mu_F)$ parameter space that may be probed using direct gravitational wave detectors. In particular, we consider the observational constraints that result from NANOGRAV data and show that heavier strings generate a secondary exclusion region of parameter space.

In the letter by Stadnik and Flambaum [Phys. Rev. Lett. 113, 151301 (2014)] it is claimed that topological defects passing through pulsars could be responsible for the observed pulsar glitches. Here, we show that, independently of the detailed network dynamics and defect dimensionality, such proposal is faced with serious difficulties.

We show that simple thermodynamic conditions determine, to a great extent, the equation of state and dynamics of cosmic defects of arbitrary dimensionality. We use these conditions to provide a more direct derivation of the Velocity-dependent One-Scale (VOS) model for the macroscopic dynamics of topological defects of arbitrary dimensionality in a $N+1$-dimensional homogeneous and isotropic universe. We parameterize the modifications to the VOS model associated to the interaction of the topological defects with other fields, including, in particular, a new dynamical degree of freedom associated to the variation of the mass per unit $p$-area of the defects, and compute the corresponding scaling solutions. The observational impact of this new dynamical degree of freedom is also briefly discussed.

We develop a numerical tool for the fast computation of the temperature and polarization power spectra generated by domain wall networks, by extending the publicly available CMBACT code --- that calculates the CMB signatures generated by active sources --- to also describe domain wall networks. In order to achieve this, we adapt the Unconnected Segment model for cosmic strings to also describe domain wall networks, and use it to model the energy-momentum of domain wall networks throughout their cosmological history. We use this new tool to compute and study the TT, EE, TE and BB power spectra generated by standard domain wall networks, and derive a conservative constraint on the energy scale of the domain wall-forming phase transition of $\upeta <0.92\,\,{\rm MeV}$ (which is a slight improvement over the original Zel'dovich bound of $1\,\,{\rm MeV}$).

We consider the possibility that current hints of spatial variations of the fine structure constant at high redshift, based on VLT/UVES and Keck/HIRES observations, could be caused by a biased domain wall network described by a scalar field non-minimally coupled to the electromagnetic field. We show that in order to be responsible for the reported spatial variations of the fine structure constant, the fractional contribution of the domain wall network to the energy density of the Universe should be tightly constrained within the range $10^{-10} < \Omega_{w0} < 10^{-5}$. We also show that the domain wall dynamics should be essentially frictionless, so that its characteristic scale is in the order of the Hubble radius at the present time.

We consider an alternative approach for the computation of the stochastic gravitational wave background generated by small loops produced throughout the cosmological evolution of cosmic string networks and use it to derive an analytical approximation to the corresponding power spectrum. We show that this approximation produces an excellent fit to more elaborate results obtained using the Velocity-dependent One-Scale model to describe cosmic string network dynamics, over a wide frequency range, in the small-loop regime.

We compute the power spectrum of the stochastic gravitational wave background generated by cosmic string networks described by the Velocity-Dependent One-Scale (VOS) model, for a wide range of macroscopic and microscopic parameters. The VOS model --- which has been shown to provide an accurate macroscopic description of the evolution of cosmic string networks --- is used to demonstrate that cosmic string networks are unable to rapidly attain scale-invariant evolution after the transition between the radiation and matter eras. However, in computations of the stochastic gravitational wave background, it is often assumed that the networks experience scale-invariant evolution throughout cosmological history. We demonstrate that this assumption leads to an underestimation of the amplitude and broadness of the peak of the spectrum, that may consequently lead to inaccurate observational constraints on the cosmic string tension.

In this paper we find new scaling laws for the evolution of $p$-brane networks in $N+1$-dimensional Friedmann-Robertson-Walker universes in the weakly-interacting limit, giving particular emphasis to the case of cosmic superstrings ($p=1$) living in a universe with three spatial dimensions (N=3). In particular, we show that, during the radiation era, the root-mean-square velocity is ${\bar v} =1/{\sqrt 2}$ and the characteristic length of non-interacting cosmic string networks scales as $L \propto a^{3/2}$ ($a$ is the scale factor), thus leading to string domination even when gravitational backreaction is taken into account. We demonstrate, however, that a small non-vanishing constant loop chopping efficiency parameter $\tilde c$ leads to a linear scaling solution with constant $L H \ll 1$ ($H$ is the Hubble parameter) and ${\bar v} \sim 1/{\sqrt 2}$ in the radiation era, which may allow for a cosmologically relevant cosmic string role even in the case of light strings. We also determine the impact that the radiation-matter transition has on the dynamics of weakly interacting cosmic superstring networks.

In this paper we derive, directly from the Nambu-Goto action, the relevant components of the acceleration of cosmological featureless $p$-branes, extending previous analysis based on the field theory equations in the thin-brane limit. The component of the acceleration parallel to the velocity is at the core of the velocity-dependent one-scale model for the evolution of $p$-brane networks. We use this model to show that, in a decelerating expanding universe in which the $p$-branes are relevant cosmologically, interactions cannot lead to frustration, except for fine-tuned non-relativistic networks with a dimensionless curvature parameter $k \ll 1$. We discuss the implications of our findings for the cosmological evolution of $p$-brane networks.

We develop a velocity-dependent one-scale model describing p-brane dynamics in flat homogeneous and isotropic backgrounds in a unified framework. We find the corresponding scaling laws in frictionless and friction dominated regimes considering both expanding and collapsing phases.

We develop a velocity-dependent one-scale model for the evolution of domain wall networks in flat expanding or collapsing homogeneous and isotropic universes with an arbitrary number of spatial dimensions, finding the corresponding scaling laws in frictionless and friction dominated regimes. We also determine the allowed range of values of the curvature parameter and the expansion exponent for which a linear scaling solution is possible in the frictionless regime.

The Press-Ryden-Spergel (PRS) algorithm is a modification to the field theory equations of motion, parametrized by two parameters ($\alpha$ and $\beta$), implemented in numerical simulations of cosmological domain wall networks, in order to ensure a fixed comoving resolution. In this paper we explicitly demonstrate that the PRS algorithm provides the correct domain wall dynamics in $N+1$-dimensional Friedmann-Robertson-Walker (FRW) universes if $\alpha+\beta/2=N$, fully validating its use in numerical studies of cosmic domain evolution. We further show that this result is valid for generic thin featureless domain walls, independently of the Lagrangian of the model.

We investigate the potential role of string and monopole-type junctions in the frustration of domain wall networks using a velocity-dependent one-scale model for the characteristic velocity, $v$, and the characteristic length, $L$, of the network. We show that, except for very special network configurations, $v^2 \lsim (HL)^2 \lsim (\rho_\sigma + \rho_\mu)/\rho_m$ where $H$ is the Hubble parameter and $\rho_\sigma$, $\rho_\mu$ and $\rho_m$ are the average density of domain walls, strings and monopole-type junctions. We further show that if domain walls are to provide a significant contribution to the dark energy without generating exceedingly large CMB temperature fluctuations then, at the present time, the network must have a characteristic length $ L_0 \lsim 10 \Omega_{\sigma 0}^{-2/3} {\rm kpc}$ and a characteristic velocity $v_0 \lsim 10^{-5} \Omega_{\sigma 0}^{-2/3}$ where $\Omega_{\sigma 0}=\rho_{\sigma 0}/\rho_{c 0}$ and $\rho_c$ is the critical density. In order to satisfy these constraints with $\Omega_{\sigma 0} \sim 1$, $\rho_{m 0}$ would have to be at least 10 orders of magnitude larger than $\rho_{\sigma 0}$, which would be in complete disagreement with observations. This result provides very strong additional support for the conjecture that no natural frustration mechanism, which could lead to a significant contribution of domain walls to the dark energy budget, exists.

We study the evolution of maximally symmetric $p$-branes with a $S_{p-i}\otimes \mathbbm{R}^i$ topology in flat expanding or collapsing homogeneous and isotropic universes with $N+1$ dimensions (with $N \ge 3$, $p < N$, $0 \le i < p$). We find the corresponding equations of motion and compute new analytical solutions for the trajectories in phase space. For a constant Hubble parameter, $H$, and $i=0$ we show that all initially static solutions with a physical radius below a certain critical value, $r_c^0$, are periodic while those with a larger initial radius become frozen in comoving coordinates at late times. We find a stationary solution with constant velocity and physical radius, $r_c$, and compute the root mean square velocity of the periodic $p$-brane solutions and the corresponding (average) equation of state of the $p$-brane gas. We also investigate the $p$-brane dynamics for $H \neq {\rm constant}$ in models where the evolution of the universe is driven by a perfect fluid with constant equation of state parameter, $w={\cal P}_p/\rho_p$, and show that a critical radius, $r_c$, can still be defined for $ -1 \le w < w_c$ with $w_c=(2-N)/N$. We further show that for $w \sim w_c$ the critical radius is given approximately by $r_c H \propto (w_c-w)^{\gamma_c}$ with $\gamma_c=-1/2$ ($r_c H \to \infty$ when $w \to w_c$). Finally, we discuss the impact that the large scale dynamics of the universe can have on the macroscopic evolution of very small loops.

We study the evolution of biased domain walls in the early universe. We explicitly discuss the roles played by the surface tension and volume pressure in the evolution of the walls, and quantify their effects by looking at the collapse of spherical wall solutions. We then apply our results to a particular mechanism, known as the devaluation scenario, in which the dynamics of biased domain walls was suggested as a possible solution to the cosmological constant problem. Our results indicate that devaluation will in general lead to values of the cosmological constant that differ by several orders of magnitude from the observationally inferred value, $\rho^{1/4}_{vac}\sim10^{-3} \rm eV$. We also argue that the reasons behind this are not specific to a particular realization, and are expected to persist in any scenario of this kind, except if a low energy cut-off on the spectra of vacuum energy densities, of the order of the critical density at the present time, is postulated. This implies that any such scenario will require a fine-tuning similar to the usual one.