Velocity-dependent One-scale Model Research Articles

Scaling solutions for current-carrying cosmic string networks

Cosmic string networks are the best motivated relics of cosmological phase transitions, being unavoidable in many physically plausible extensions of the Standard Model. Most studies, including those providing constraints from and forecasts of their observational signals, rely on assumptions of featureless networks, neglecting the additional degrees of freedom on the string worldsheet, e.g., charges and currents, which are all but unavoidable in physically realistic models. An extension of the canonical velocity-dependent one-scale model, accounting for all such possible degrees of freedom, has been recently developed. Here we improve its physical interpretation by studying and classifying its possible asymptotic scaling solutions, and in particular how they are affected by the expansion of the Universe and the available energy loss or transfer mechanisms. We find three classes of solutions. For sufficiently fast expansion rates the charges and currents decay and one asymptotes to the Nambu-Goto case, while for slower expansion rates they can dominate the network dynamics. In between the two there is a third regime in which the network, including its charge and current, reaches full scaling. Under specific but plausible assumptions, this intermediate regime corresponds to the matter-dominated era. Our results agree with, and significantly extend, those of previous studies. Published by the American Physical Society 2024

Cosmic F- and D-strings from pure Yang–Mills theory

We discuss the formation of cosmic strings or macroscopic color flux tubes after the deconfinement/confinement phase transition in the pure Yang–Mills theory. Based on holographic dual descriptions, these cosmic strings can be interpreted as fundamental (F-) strings or wrapped D-branes (which we call as D-strings) in the gravity side, depending on the structure of the gauge group. In fact, the reconnection probabilities of the F- and D-strings are suppressed by factors of 1/N2 and e−cN, where c=O(1), in a large-N limit, respectively. Supported by the picture of electric–magnetic duality, we discuss that color flux tubes form after the deconfinement/confinement phase transition, just like the formation of local cosmic strings after spontaneous symmetry breaking in the weak-U(1) gauge theory. We use an extended velocity-dependent one-scale model to describe the dynamics of the string network and calculate the gravitational wave signals from string loops. We also discuss the dependence on the size of produced string loops.

Cosmic strings from pure Yang–Mills theory

We discuss the formation of cosmic strings or macroscopic color flux tubes at the phase transition from the deconfinement to confinement phase in pure Yang--Mills (YM) theory, such as $\mathrm{SU}(N)$, $\mathrm{Sp}(N)$, $\mathrm{SO}(N)$, and $\text{Spin}(N)$, based on the current understanding of theoretical physics. According to the holographic dual descriptions, the cosmic strings are dual to fundamental strings or wrapped D-branes in the gravity side depending on the structure of the gauge group, and the reconnection probability is suppressed by $\mathcal{O}({N}^{\ensuremath{-}2})$ and ${e}^{\ensuremath{-}\mathcal{O}(N)}$, respectively. The pure YM theory thus provides a simple realization of cosmic F- and D-strings without the need for a brane-inflationary scenario or extra dimension. We also review the stability of cosmic strings based on the concept of 1-form symmetry, which further implies the existence of a baryon vertex in some YM theory. We calculate the gravitational wave spectrum that is emitted from the cosmic strings based on an extended velocity-dependent one-scale model and discuss its detectability based on ongoing and planned gravitational-wave experiments. In particular, the SKA and LISA can observe gravitational signals if the confinement scale is higher than $\mathcal{O}({10}^{12})\text{ }\text{ }\mathrm{GeV}$ and $\mathcal{O}({10}^{10})\text{ }\text{ }\mathrm{GeV}$ for $\mathrm{SU}(N)$ with $N=\mathcal{O}(1)$, respectively.

Analytical Scaling Solutions for the Evolution of Cosmic Domain Walls in a Parameter-Free Velocity-Dependent One-Scale Model

We derive an analytical approximation for the linear scaling evolution of the characteristic length L and the root-mean-squared velocity σv of standard frictionless domain wall networks in Friedmann–Lemaître–Robertson–Walker universes with a power law evolution of the scale factor a with the cosmic time t (a∝tλ). This approximation, obtained using a recently proposed parameter-free velocity-dependent one-scale model for domain walls, reproduces well the model predictions for λ close to unity, becoming exact in the λ→1− limit. We use this approximation, in combination with the exact results found for λ=0, to obtain a fit to the model predictions valid for λ∈[0,1] with a maximum error of the order of 1%. This fit is also in good agreement with the results of field theory numerical simulations, especially for λ∈[0.9,1]. Finally, we explicitly show that the phenomenological energy-loss parameter of the original velocity-dependent one-scale model for domain walls vanishes in the λ→1− limit and discuss the implications of this result.

Gravitational waves from global cosmic strings and cosmic archaeology

Global cosmic strings are predicted in many motivated extensions to the Standard Model of particle physics, with close connections to axion dark matter physics. Recent studies suggest that, although subdominant relative to Goldstone emission, gravitational wave (GW) signals from global strings can be detectable with current and planned GW detectors such as LIGO, LISA, DECIGO/BBO, ET/CE and AEDGE/AION, as well as pulsar timing arrays such as PPTA, NANOGrav and SKA. This work is an extensive, updated study on GWs from a global cosmic string network, taking into account of the most recent developments related to the subject. The main analysis is based on the analytical Velocity-dependent One-Scale (VOS) model calibrated with recent simulation results, which provides a generic protocol for such calculations with details given. We also demonstrate how the GW signal can be influenced with variations to the baseline model: this includes considering the uncertainties of model parameters and the potential deviation from the conventional VOS model prediction (i.e. the scaling behavior) as suggested by some of the recent simulation results. Furthermore, we investigated in detail the effect of a non-standard cosmology (e.g. early matter domination or kination) or new particle species on the GW signals from global strings. We demonstrate that the frequency spectrum of GW background from global cosmic strings can be used to probe the cosmic history prior to the Big Bang nucleosynthesis (BBN) (i.e. the primordial dark age) up to a temperature of T ∼ 108 GeV.

High resolution calibration of the cosmic strings velocity dependent one-scale model

The canonical velocity-dependent one-scale (VOS) model for cosmic string evolution must be calibrated using high resolution numerical simulations, We exploit our state of the art graphics processing unit accelerated implementation of the evolution of local Abelian-Higgs string networks to provide a detailed and statistically robust calibration of the VOS model. We rely on the largest set of high resolution simulations carried out to date, with a wide range of cosmological expansion rates, and explore the impact of key numerical parameters, including the dynamic range (comparing box sizes from $1024^3$ to $4096^3$), the lattice spacing, and the choice of numerical estimators for the string velocity. We explore the sensitivity of the VOS model parameters to these numerical parameters, with a particular emphasis on the observationally crucial loop chopping efficiency, and also identify key differences between the equation of state and conjugate momentum estimators for the string velocities, showing that the latter one is more reliable for fast expansion rates (while in Minkowski space the opposite has been previously shown). Finally, we briefly illustrate how our results impact observational constraints on cosmic strings.

Scaling solutions of wiggly cosmic strings

Cosmic strings may have formed in the early universe due to the Kibble mechanism. While string networks are usually modeled as being of Nambu-Goto type, this description is understood to be a convenient approximation, which neglects the typically expected presence of additional degrees of freedom on the string worldsheet. Previous simulations of cosmic strings in expanding universes have established beyond doubt the existence of a significant amount of short-wavelength propagation modes (commonly called wiggles) on the strings, and a wiggly string extension of the canonical velocity-dependent one-scale model has been recently developed. Here we improve the physical interpretation of this model, by studying the possible asymptotic scaling solutions of this model, and in particular how they are affected by the expansion of the universe and the available energy loss or transfer mechanisms -- e.g., the production of loops and wiggles. In addition to the Nambu-Goto solution, to which the wiggly model reduces in the appropriate limit, we find that there are also solutions where the amount of wiggliness can grow as the network evolves or, for specific expansion rates, become a constant. Our results show that full scaling of the network, including the wiggliness, is much more likely in the matter era than in the radiation era, which is in agreement with numerical simulation results.

Generalized velocity-dependent one-scale model for current-carrying strings

We develop an analytic model to quantitatively describe the evolution of superconducting cosmic string networks. Specifically, we extend the velocity-dependent one-scale (VOS) model to incorporate arbitrary currents and charges on cosmic string world sheets under two main assumptions, the validity of which we also discuss. We derive equations that describe the string network evolution in terms of four macroscopic parameters: the mean string separation (or, alternatively, the string correlation length) and the root mean square velocity which are the cornerstones of the VOS model, together with parameters describing the averaged timelike and spacelike current contributions. We show that our extended description reproduces the particular cases of wiggly and chiral cosmic strings, previously studied in the literature. This VOS model enables investigation of the evolution and possible observational signatures of superconducting cosmic string networks for more general equations of state, and these opportunities will be exploited in a companion paper.

Quantifying the effect of cooled initial conditions on cosmic string network evolution

Quantitative studies of the evolution and cosmological consequences of networks of cosmic strings (or other topological defects) require a combination of numerical simulations and analytic modeling with the velocity-dependent one-scale (VOS) model. In previous work, we demonstrated that a GPU-accelerated code for local Abelian-Higgs string networks enables a statistical separation of key dynamical processes affecting the evolution of the string networks and thus a precise calibration of the VOS model. Here we further exploit this code in a detailed study of two important aspects connecting the simulations with the VOS model. First, we study the sensitivity of the model calibration to the presence (or absence) of thermal oscillations due to high gradients in the initial conditions. This is relevant since in some Abelian-Higgs simulations described in the literature a period of artificial (unphysical) dissipation---usually known as cooling---is introduced with the goal of suppressing these oscillations and accelerating the convergence to scaling. We show that a small amount of cooling has no statistically significant impact on the VOS model calibration, while a longer dissipation period does have a noticeable effect. Second, in doing this analysis we also introduce an improved Markov Chain Monte Carlo based pipeline for calibrating the VOS model, Comparison to our previous bootstrap based pipeline shows that the latter accurately determined the best-fit values of the VOS model parameter, but underestimated the uncertainties in some of the parameters. Overall, our analysis shows that the calibration pipeline is robust and can be applied to future much larger field theory simulations.

Comparing parametric and non-parametric velocity-dependent one-scale models for domain wall evolution

We perform a detailed comparison between a recently proposed parameter-free velocity-dependent one-scale model and the standard parametric model for the cosmological evolution of domain wall networks. We find that the latter overestimates the damping of the wall motion due to the Hubble expansion and neglects the direct impact of wall decay on the evolution of the root-mean-square velocity of the network. We show that these effects are significant but may be absorbed into a redefinition of the momentum parameter. We also discuss the implications of these findings for cosmic strings. We compute the energy loss and momentum parameters of the standard parametric model for cosmological domain wall evolution using our {non-parametric} velocity-dependent one-scale model in the context of cosmological models having a power law evolution of the scale factor a with the cosmic time t (a ∝ tλ, 0 < λ < 1), and compare with the results obtained from numerical field theory simulations. We further provide simple linear functions which roughly approximate the dependence of the energy loss and momentum parameters on λ.

Parameter-free velocity-dependent one-scale model for domain walls

We develop a parameter-free velocity-dependent one-scale model for the evolution of the characteristic length $L$ and root-mean-square velocity $\sigma_v$ of standard domain wall networks in homogeneous and isotropic cosmologies. We compare the frictionless scaling solutions predicted by our model, in the context of cosmological models having a power law evolution of the scale factor $a$ as a function of the cosmic time $t$ ($a \propto t^\lambda$, $0< \lambda < 1$), with the corresponding results obtained using field theory numerical simulations. We show that they agree well (within a few $\%$) for root-mean-square velocities $\sigma_v$ smaller than $0.2 \, c$ ($\lambda \ge 0.9$), where $c$ is the speed of light in vacuum, but significant discrepancies occur for larger values of $\sigma_v$ (smaller values of $\lambda$). We identify problems with the determination of $L$ and $\sigma_v$ from numerical field theory simulations which might potentially be responsible for these discrepancies.

Extending and calibrating the velocity dependent one-scale model for cosmic strings with one thousand field theory simulations

Understanding the evolution and cosmological consequences of topological defect networks requires a combination of analytic modeling and numerical simulations. The canonical analytic model for defect network evolution is the Velocity-dependent One-Scale (VOS) model. For the case of cosmic strings, this has so far been calibrated using small numbers of Goto-Nambu and field theory simulations, in the radiation and matter eras, as well as in Minkowski spacetime. But the model is only as good as the available simulations, and it should be extended as further simulations become available. In previous work we presented a General Purpose Graphics Processing Unit implementation of the evolution of cosmological domain wall networks, and used it to obtain an improved VOS model for domain walls. Here we continue this effort, exploiting a more recent analogous code for local Abelian-Higgs string networks. The significant gains in speed afforded by this code enabled us to carry out 1032 field theory simulations of $512^3$ size, with 43 different expansion rates. This detailed exploration of the effects of the expansion rate on the network properties in turn enables a statistical separation of various dynamical processes affecting the evolution of the network. We thus extend and accurately calibrate the VOS model for cosmic strings, including separate terms for energy losses due to loop production and scalar/gauge radiation. By comparing this newly calibrated VOS model with the analogous one for domain walls we quantitatively show that energy loss mechanisms are different for the two types of defects.

Dynamics of junctions and the multitension velocity-dependent one-scale model

The dynamics of string junctions and their influence on the evolution of cosmic superstring networks are studied in full detail. We review kinematic constraints for colliding strings in a Friedmann-Lema\^itre-Robertson-Walker background and obtain the average distribution of possible string configurations after string collisions. The study of small-scale structure enables us to investigate the average growth/reduction rate of string junctions for a given cosmic string network. Incorporating the averaged junction dynamics into the velocity-dependent one-scale model for multi-tension string networks, we improve the semi-analytic description and quantitative understanding of cosmic superstring network evolution.

Scaling properties of cosmological axion strings

There has been recent interest in the evolution and cosmological consequences of global axionic string networks, and in particular in the issue of whether or not these networks reach the scale-invariant scaling solution that is known to exist for the simpler Goto–Nambu and Abelian–Higgs string networks. This is relevant for determining the amount and spectrum of axions they produce. We use the canonical velocity-dependent one-scale model for cosmic defect network evolution to study the evolution of these global networks, confirming the presence of deviations to scale-invariant evolution and in agreement with the most recent numerical simulations. We also quantify the cosmological impact of these corrections and discuss how the model can be used to extrapolate the results of numerical simulations, which have a limited dynamic range, to the full cosmological evolution of the networks, enabling robust predictions of their consequences. Our analysis suggests that around the QCD scale, when the global string network is expected to disappear and produce most of the axions, the number of global strings per Hubble patch should be around ξ∼4.2, but also highlights the need for additional high-resolution numerical simulations.

Effects of biases in domain wall network evolution. II. Quantitative analysis

Domain walls form at phase transitions which break discrete symmetries. In a cosmological context they often overclose the universe (contrary to observational evidence), although one may prevent this by introducing biases or forcing anisotropic evolution of the walls. In a previous work [Correia {\it et al.}, Phys.Rev.D90, 023521 (2014)] we numerically studied the evolution of various types of biased domain wall networks in the early universe, confirming that anisotropic networks ultimately reach scaling while those with a biased potential or biased initial conditions decay. We also found that the analytic decay law obtained by Hindmarsh was in good agreement with simulations of biased potentials, but not of biased initial conditions, and suggested that the difference was related to the Gaussian approximation underlying the analytic law. Here we extend our previous work in several ways. For the cases of biased potential and biased initial conditions we study in detail the field distributions in the simulations, confirming that the validity (or not) of the Gaussian approximation is the key difference between the two cases. For anisotropic walls we carry out a more extensive set of numerical simulations and compare them to the canonical velocity-dependent one-scale model for domain walls, finding that the model accurately predicts the linear scaling regime after isotropization. Overall, our analysis provides a quantitative description of the cosmological evolution of these networks.

Revisiting the velocity-dependent one-scale model for monopoles

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.

Cosmic strings and other topological defects in nonscaling regimes

Cosmic strings are topological defects possibly formed in the early Universe, which may be observable due to their gravitational effects on the cosmic microwave background radiation or gravitational wave experiments. To this effect it is important to quantitatively ascertain the network properties, including their density, velocity or the number of strings present, at the various epochs in the observable Universe. Attempts to estimate these numbers often rely on simplistic approximations for the string parameters, such as assuming that the network is scaling. However, in cosmological models containing realistic amounts of radiation, matter and dark energy a string network is never exactly scaling. Here we use the velocity-dependent one-scale model for the evolution of a string network to better quantify how these networks evolve. In particular we obtain new approximate analytic solutions for the behavior of the network during the radiation-to-matter and matter-to-acceleration transitions (assuming, in the latter case, the canonical $\Lambda$ cold dark matter model), and numerically calculate the relevant quantities for a range of possible dark energy models.

Measuring global monopole velocities, one by one

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.

Stretching and Kibble scaling regimes for Hubble-damped defect networks

The cosmological evolution of topological defect networks can broadly be divided into two stages. At early times they are friction-dominated due to particle scattering and therefore non-relativistic, and may either be conformally stretched or evolve in the Kibble regime. At late times they are relativistic and evolve in the well known linear scaling regime. In this work we show that a sufficiently large Hubble damping (that is a sufficiently fast expansion rate) leads to a linear scaling regime where the network is non-relativistic. This is therefore another realization of a Kibble scaling regime, and also has a conformal stretching regime counterpart which we characterize for the first time. We describe these regimes using analytic arguments in the context of the velocity-dependent one-scale model, and we confirm them using high-resolution $4096^3$ field theory simulations of domain wall networks. We also use these simulations to improve the calibration of this analytic model for the case of domain walls.

Models for small-scale structure on cosmic strings. II. Scaling and its stability

We make use of the formalism described in a previous paper [Martins {\it et al.} Phys. Rev. D90 (2014) 043518] to address general features of wiggly cosmic string evolution. In particular, we highlight the important role played by poorly understood energy loss mechanisms and propose a simple ansatz which tackles this problem in the context of an extended velocity-dependent one-scale model. We find a general procedure to determine all the scaling solutions admitted by a specific string model and study their stability, enabling a detailed comparison with future numerical simulations. A simpler comparison with previous Goto-Nambu simulations supports earlier evidence that scaling is easier to achieve in the matter era than in the radiation era. In addition, we also find that the requirement that a scaling regime be stable seems to notably constrain the allowed range of energy loss parameters.