## Testing the Gaussianity and Statistical Isotropy of the Universe

View this Special IssueReview Article | Open Access

# Cosmic Strings and Their Induced Non-Gaussianities in the Cosmic Microwave Background

**Academic Editor:**Eiichiro Komatsu

#### Abstract

Motivated by the fact that cosmological perturbations of inflationary quantum origin were born Gaussian, the search for non-Gaussianities in the cosmic microwave background (CMB) anisotropies is considered as the privileged probe of nonlinear physics in the early universe. Cosmic strings are active sources of gravitational perturbations and incessantly produce non-Gaussian distortions in the CMB. Even if, on the currently observed angular scales, they can only contribute a small fraction of the CMB angular power spectrum, cosmic strings could actually be the main source of its non-Gaussianities. In this paper, after having reviewed the basic cosmological properties of a string network, we present the signatures Nambu-Goto cosmic strings would induce in various observables ranging from the one-point function of the temperature anisotropies to the bispectrum and trispectrum. It is shown that string imprints are significantly different than those expected from the primordial type of non-Gaussianity and could therefore be easily distinguished.

#### 1. Motivations

The origin of cosmic strings dates back to the discovery that cosmological phase transitions triggered by the spontaneous breakdown of the fundamental interaction symmetries may form topological defects [1–3]. Cosmic strings belong to the class of line-like topological defects, as opposed to point-like monopoles and the membrane shaped domain walls. As shown by Kibble, the appearance of defects in any field theory is related to the topology of the vacuum manifold [3]. If the ground state of a field theory experiences a spontaneous breakdown from a symmetry group to a subgroup , Kibble showed that cosmic strings will be formed if the first homotopy group is nontrivial. In other words, if noncontractile loops can be found in the manifold of equivalent vacua. Similarly, the other homotopy groups and determine the formation of domain walls and monopoles, respectively. Once formed and cooled, these defects cannot be unfolded, precisely due to their nontrivial topological configuration over the vacuum manifold of the theory. This simple statement suggests that cosmic strings, and topological defects in general, are a natural outcome of the unification of the fundamental interactions in the context of Cosmology. As remnants of unified forces, their discovery would be an incredible opportunity to probe extremely high-energy physics with “a telescope”.

In the last thirty years, many works have been devoted to the cosmological consequences, signatures, and searches for topological defects [4–7]. They have pushed cosmic strings to the privileged place to be generically compatible with observations. Indeed, domain walls and monopoles are prone to suffer from the cosmological catastrophe problem; their formation is sufficiently efficient (or their annihilation sufficiently inefficient) to either overclose the universe or spoil the Big-Bang Nucleosynthesis (BBN) predictions [8, 9]. For domain walls, this implies that either they should be extremely light, that is, formed at an energy scale less than a few , or no discrete symmetry should have been broken during the cooling of the universe. There is not so much choice for the monopoles; if interactions were unified, monopoles would have been formed. The homotopy group of with containing the of electroweak interactions is indeed nontrivial.As often with topological defects, sensitivity to the underlying model is such that one can often find a counter-example of any result. Both of these statements, on walls and monopoles, can be evaded in some particular models or with some amount of fine-tuning, as for instance if cosmic strings can be attached to them and catalyse annihilations [10, 11].Cosmic inflation was originally designed to solve the monopole problem. If a phase of accelerated expansion of the universe occurs, then any defects will be diluted enough to no longer have any (dramatic) consequences on cosmology [12–15]. Meanwhile, Inflationary Cosmology solves the flatness and homogeneity problem of the standard Big-Bang model, explains the origin and spectrum of the cosmic microwave background (CMB) anisotropies, as the formation of the large-scale structures [16–18, 17]. Inflation provides a priori an easy solution to the topological defects problem by diluting them to at most one per Hubble radius. However, one has to keep in mind that this mechanism works only if the defects were formed before inflation, and even in that case some may survive [20]. This has to be the case for monopoles and heavy walls, but not for local strings. On the contrary, exhaustive analysis of particle physics motivated inflationary models, embedding the Standard Model , has shown that strings are generically produced at the end of inflation [21]. In this picture, our universe should contain cosmic strings whose properties are closely related to those of the inflation [22–24]. String Theory provides an alternative framework to Field Theories: brane inflationary models propose that the accelerated expansion of the universe is induced by the motion of branes in warped and compact extradimensions [25–28]. Inflation ends when two branes collide and such a mechanism again triggers the formation of one-dimensional cosmological extended objects, dubbed cosmic superstrings [29–32]. These objects may be cosmologically stretched fundamental strings or one-dimensional D-brane [33, 34]. Although cosmic superstrings are of a different nature than their topological analogue, they produce the same gravitational effects and share similar cosmological signatures [35, 36].

Among the expected signatures, cosmic strings induce temperature anisotropies in the CMB with an amplitude typically given by , where is the string energy per unit length and the Newton constant [37]. (To avoid any confusion with Greek tensor indices, we will use the Carter’s notations and for the string energy density and tension [38].) For the Grand Unified Theory (GUT) energy scale, one has , which precisely corresponds to the observed amplitude of the CMB temperature fluctuations [39]. However, the power spectra do not match; topological defects are active sources of gravitational perturbations, that is, they produce perturbations all along the universe history, and cannot produce the characteristic coherent patterns of the acoustic peaks [40–44]. Current CMB data analyses including a string contribution suggest that they can only contribute to at most of the overall anisotropies on the observed angular scales [45, 46]. For Abelian cosmic strings (see Section 2), numerical simulations in Friedmann-Lemaître-Robertson-Walker (FLRW) spacetimes show that this corresponds to an upper two-sigma bound [47]. Direct detection searches provide less stringent limits but are applicable to all cosmic string models: [48–50]. Detecting cosmic strings in the CMB certainly requires one to go further than the power spectrum [51, 52] (see, however, Section 4.5). In fact, strings induce line-like discontinuities in the CMB temperature through the so-called Gott-Kaiser-Stebbins effect, which are intrinsically of non-Gaussian nature [53, 54]. In the inflationary picture, cosmological perturbations find their origin in the quantum fluctuations of the field-metric system, and therefore were born generically Gaussian. Non-Gaussianities can nevertheless appear from non-linear effects during inflation or from couplings to other fields (see the other articles in this issue). These non-Gaussianities are of the primordial type, that is, they exist before the cosmological perturbations reenter the Hubble radius. On the other hand, cosmic strings are a source of non-Gaussianity at all times and, as we will see, produce different signals from the CMB point of view.Notice that second-order perturbations, being non-linear, actively generate non-Gaussianities but at a relatively small amplitude [55–58].

In this paper, we review the non-Gaussian features a cosmological network of cosmic strings produce in the CMB anisotropies. In a first section, we briefly scan various cosmic string models and emphasize their similarities and differences for cosmology. Making observable predictions for cosmic strings faces the problem of understanding their cosmological evolution. Not only one has to solve the local dynamics in curved space, but as extended objects, cosmic strings follow a globally nonlocal evolution: the fate of one string depends on its interactions with the others. The cosmological evolution of a network of cosmic strings is a nontrivial problem which can be overcome by means of numerical simulations. These simulations permit an estimation of the various statistical properties affecting the observational signatures, such as the number of strings per Hubble radius, their shapes, velocities, or the loop density distribution. Latest results in this area, for the Nambu-Goto (NG) type of cosmic strings, are presented in Section 3. Once the statistical properties of a cosmological cosmic strings network are known, it is possible to extract meaningful observables depending only on the unique model parameter . (If no currents are flowing along the string, Lorentz invariance implies that the string tension equals the energy density .) In Section 4, we recap the expected CMB temperature anisotropies induced by cosmic strings, derived from various methods. Particular attention is paid to small angle CMB maps which preserve all of the projected statistical information. We then derive the cosmic string signals expected in various non-Gaussian estimators ranging from the one-point function of the CMB temperature fluctuations to the bispectrum and trispectrum. We conclude in Section 5 and discuss various non-Gaussian aspects which still have to be explored.

#### 2. Cosmic Strings of Various Origins

Cosmic strings of cosmological interest can be of various kinds depending on the microscopic model they stem from. As mentioned in the introduction, they can either be nontrivial stable, or metastable, field configurations or more fundamental objects in String Theory. From a gravitational point of view, they all are, however, line-like energy density and pressure distributions. In the following, we briefly review the different kinds of string having a cosmological interest and we emphasize their similarities and differences.

##### 2.1. Abelian Vortices

The simplest example of cosmic string illustrating the Kibble mechanism is the Abelian Higgs model. The theory is invariant under a local gauge group and the Higgs potential assumes its standard Mexican hat renormalisable formwhere is the self-coupling constant and the vacuum expectation value of the Higgs field . In Minkowski space, the Lagrangian readswhere is the field strength tensor associated with the vector gauge boson and At high enough temperature, loop corrections from the thermal bath restore the symmetry and the effective potential has an overall minimum at [1, 59]. Starting from high enough temperature, one therefore expects the symmetry to be spontaneously broken during the expansion and cooling of the universe. During the phase transition, the Higgs field reaches its new vacuum expectation value . At each spacetime location, the phase will have a given value, all of them being uncorrelated on distances larger than the typical correlation length of the phase transition. As pointed by Kibble, this is at most the horizon size although one expects it to be much smaller [3, 60–63]. As a result, there exists closed paths in space along which varies from to (or a multiple of ). Such phase configurations necessarily encompass a point at which (see Figure 1); the old vacuum has been trapped into a nontrivial configuration of the new vacuum, and this prevents its decay. Such a structure is invariant by translations along the third spatial dimension and is string shaped.

Solitonic solutions of the field equations describing a static straight Abelian string can easily be computed under the Nielsen-Olesen ansatz. The transverse profile of the Higgs and gauge field are assumed to be [64]where stands for a polar coordinate system aligned along the string. The dimensionless radial coordinate has been defined by where, is the mass of the Higgs boson. The integer is the “winding number” and gives the number of times the Higgs winds the potential for one rotation around the string. From (2), the dimensionless equations of motion readwhere is the mass of the vector gauge boson. In Figure 2, we have represented the string solution to these equations in Minkowski spacetime [65, 66]. The boundary conditions are such that the Higgs field vanishes at the center of the string to reach its vacuum expectation value (vev) asymptotically. This typically happens after a length scale given by its Compton wavelength . Similarly, the gauge field boundary conditions are such that it has vanishing derivative in the core and remains finite far from the string. As shown in Figure 2, it actually condenses inside the string with a spatial extension roughly equal to .

The energy content of such a string is given by the stress tensor stemming from the Lagrangian of (2). Along the string worldsheet,which are the only two components which do not vanish after an integration over the transverse coordinates. Integrating the temporal part gives the string energy per unit length , whereas the longitudinal component gives . One finally getswhere is an order unity function at fixed winding number. Increasing the winding number centrifuges the energy density around the core such that is changed in a more complex way [5]. This immediately shows that cosmic strings generically carry an energy density and tension of the order of the symmetry breaking energy scale . Notice that along the string direction the pressure is negative, and we are in presence of a “cosmological constant wire”, as one may expect from a Lorentz invariant vacuum object. Consequently, the trace of the stress tensor vanishes and cosmic strings do not induce any Newtonian gravitational potential. Together with the so-called cosmological scaling behaviour (see below), this is the mere reason why they remain cosmologically acceptable. They do however induce dynamical gravitational effects, the metric far from the string core being Minkowski with a missing angle [67] (see Section 4).

The Abelian string model is intensively used in the literature to explore the string forming phase transition and string interactions [68, 69]. From the Kibble’s argument, one expects the phase of the Higgs field to be random and the resulting string path should be a self-avoiding random walk with a given correlation length [70]. Performing lattice simulations allows to probe in more details the string forming mechanism and gives a more accurate picture of a cosmic string network just after its formation [71–74]. Abelian Higgs simulations are also used to compute the cosmological evolution of such a network [75–77] (see Section 3).

##### 2.2. Other Flux Tubes

*Global String*. The Abelian string provides an explicit example of the formation of line-like topological defects by the spontaneous breakdown of a gauged symmetry. Breaking a global symmetry can also produce topological defects, the so-called global strings. However, in the absence of gauge fields, one can show that global cosmic strings exhibit long-range interactions and Goldstone radiation [78]. Their dynamics can however mimic local strings and being cosmologically acceptable in some regime [79–81].

*Non-Abelian String.* If the broken symmetry group is non-Abelian, the cosmic strings formed during the phase transition exhibit new properties compared to the kind [5]. In particular, the mapping of the Higgs field to the real space can be made along different broken generators . This implies that different type of non-Abelian strings may be formed and will interact with each others according to their respective windings. The classic example being the appearance of a -string from the crossing between a -string and -string [82]. In the cosmological framework, new strings can potentially be formed at each interaction leading to a frustrated intricate configuration [83–86]. Such an outcome depends on the underlying non-Abelian gauge group and Abelian string-like evolution can also be recovered, as for instance in the models [87–91].

*Semilocal String.* String-shaped energy density distribution can also appear even if the vacuum manifold is simply connected. The nontrivial topology argument is indeed only a sufficient condition of defect appearance. The electroweak symmetry breaking scheme enters this class, although the first homotopy group is trivial, semilocal strings can be formed [92, 93]. A simple description of these strings can be obtained by replacing the Higgs field in the Abelian model by a doublet in a global representation [94]. These flux tube configurations are stabilised because they can be energetically favoured for some values of the model parameters, typically for [95]. Let us notice that the currently measured electroweak model parameters do not support stable vortices [94]. Similar vortices could also be formed during a chiral symmetry breaking phase transition [96–103].

*K- and DBI-String.* These are another extensions of the Abelian Higgs string for which the scalar and gauge field kinetic terms are noncanonical, or of the Dirac-Born-Infeld form [104–107]. These strings essentially differ from their Abelian counterparts when the gradient terms are non-vanishing, that is, in the core.

*Current Carrying String.* In minimal extensions of the Abelian Higgs model, one may couple extra-scalar fields to the string forming Higgs field. As shown by Witten, this can lead to the condensation of the extrascalar over the string core [108]. The resulting string is carrying a current that breaks the longitudinal Lorentz invariance; and are no longer degenerated and the string dynamics is affected [38, 65, 109]. One of the most important consequence of these currents is the potential appearance of centrifugally supported loops. If stable, these so-called vortons could efficiently populate the universe and avoiding the overclosure gives strong constraints on the cosmic string energy scale [110, 111]. A similar mechanism works for the fermionic fields which are Yukawa coupled to the string forming Higgs field. They generically produce currents along the string with a discrete mass spectrum, in a way similar to the photon propagation in waveguides [112]. Unless the massive propagation modes are not excited, the resulting loops are however expected to be unstable [113, 114].

##### 2.3. Cosmic Superstrings

Cosmic superstrings are fundamental line-shaped objects that can be formed at the end of brane-inflation (see [34–36, 115, 116] for reviews). The idea that fundamental quantum strings can be stretched to cosmological distances has been mentioned by Witten [117]. If stable, one would expect fundamental strings to be at an energy scale close to the String Theory scale, that is, close to the Planck mass, and this is trivially ruled out by observations. In addition, current CMB constraints tell us that the energy scale of inflation is at most the GUT scale [118], implying that strings formed at a higher energy would have been diluted anyway. The situation changed with the discovery that inflation within String Theory could be a geometrical phenomena induced by the motion of a brane moving in a warped throat, somewhere in the compact manifold of the extra-dimensions [26]. In the KKLMMT model [28, 119], the inflaton is a scalar degree of freedom associated with the position of a D3 brane in a warped throat. Within a ten-dimensional super-gravity ansatz for the metric, in the type IIB String Theory,and the throat is described by the warping function (explicitly, it can be the Klebanov-Strassler conifold [120]). In this system, accelerated expansion of our universe comes from the interaction of this brane with an anti-D3 brane sitting at the bottom of the throat . Current CMB data suggest that inflation preferentially ends by violation of the slow-roll conditions, whereas the system continues to evolve till the two branes collide [121]. The brane interactions at that stage require String Theory calculations and are expected to trigger a reheating era accompanied by a copious production of various D1-branes and fundamental F-strings [29–32]. Since the brane annihilation takes place at the bottom of the throat, due to the warped metric, the cosmic superstring tensions measured by an exterior four-dimensional observer are redshifted by a factor . The resulting effect is to significantly lower the string tension down to acceptable values. In fact, the stability of the produced F-strings and D-strings require additional constraints on the model parameters and the spectrum of superstring tensions depending on the underlying scenario [32]. For instance, in the KKLMMT model, one expects [122].

Cosmic superstrings differ from the Abelian strings in various aspects. In addition to the coexistence of two different types, they can form bound states of F-strings and D-strings. The tension of these -strings depends on , , the binding energy but also on their configuration in the throat [123, 124]. In fact, many of -string properties mimic the non-Abelian type of topological vortex, as the existence of bound states and Y-junctions [125–128]. Such similarities have actually been used to probe the properties of the cosmic superstrings through the more tractable framework of field theory [129–132].

##### 2.4. Infinitely Thin Strings

These are the one-dimensional version of the relativistic point particles. Following Carter macroscopic covariant approach [38, 109, 133, 134], string events can be localised in the four-dimensional spacetime by the so-called embedding functions , where and are a timelike and spacelike internal coordinate of the string worldsheet. Denoting by the four-dimensional metric tensor, one can define the two-dimensional induced metricsuch that the infinitesimal interval between two events reduces to . From its inverse, on can define the first fundamental tensorwhich is nothing but a projector over the string worldsheet. Similarly, is an orthogonal projector and they verifyVariations of the first fundamental form are encoded in the second fundamental tensorwhere a bar quantity stands for the projection of its four-dimensional analogue over the worldsheet, that is, . Integrability imposes = 0 and, by construction, the second fundamental form is, respectively, tangent and orthogonal to the worldsheet on its first and last indices. As a result, contracting the first two tangential indices gives a purely orthogonal vector which measures the string extrinsic curvature [134]

The energy content of a spacetime two-dimensional surface can be characterised by its internal stress energy tensor. Similarly to the cosmological perfect fluid, one may consider a string whose stress-energy tensor is diagonal in a preferred basis. Positivity of the energy conditions ensures that the timelike eigenvalue , while the spacelike eigenvalue should verify [135]. In this frame, represents the energy per unit length of the string and the string tension. Denoting by and the respective timelike and spacelike orthonormal eigenvectors, one haswhereIn the absence of external forces, reparametrisation invariance of the string worldsheet ensures the stress-energy pseudo-conservation from Noether’s theorem [136]As for a cosmological fluid, these equations are not sufficient to close the equations of motion for the string. One has to supplement them by an equation of state of the fluid under scrutiny. The simplest case is the so-called barotropic model for which the equation of state is the relation . One can then introduce the two Legendre conjugated parameterssuch that . Clearly, plays the role of a number density and its Legendre conjugated parameter will therefore be a chemical potential, that is, an effective mass carried per unit number density. Defining their respective worldsheet current density byone can rewrite (16) asContracting (19) with ensures the current conservation along the stringwhile its projection onto the worldsheet gives the momentum transport lawFinally, the orthogonal projection of (19) reduces towhere the string acceleration stands for

As should be clear from (14), the barotropic equation of state breaks Lorentz invariance along the string for . In fact, it describes a wide class of elastic string models [137–139], and as suggested by (20), the scalar current carrying cosmic strings [140]. Conversely, imposing Lorentz invariance along the worldsheet reduces the equation of state to the trivial form (see (14)), which is also the relation found for the Abelian Higgs string. This infinitely thin string is the Nambu-Goto (NG) string and does not possess any internal structure [141]. The associated equations of motion are purely geometrical and do not depend on . From (22), they reduce to the vanishing of the extrinsic curvature vector, that is,which can be rewritten in a coordinate-dependant way by using (13)The connections are for the background spacetime of metric while is the determinant of the induced metric. These equations can also be recovered from the usual NG action with an explicit coordinate system [5]

#### 3. Cosmological Evolution of Nambu-GotoStrings

The previous section shows that the equations of motion of an isolated string depend on the underlying microscopic model. The type of string is more determinant when two strings interact; cosmic superstrings may form bound states, while non-Abelian vortices may weave new vortices from each of their interaction points. Understanding the cosmological evolution of a string network requires one to solve both the local equations of motion for each string and the outcome of their interactions when they meet. Moreover, the evolution of a system of strings starts from an initial configuration which should describe the network configuration just after its formation. Numerical simulations have been used to overcome some of the above-mentioned difficulties and, up to now, Friedmann-Lemaître-Robertson-Walker (FLRW) network simulations have only been performed with Nambu-Goto strings, Abelian Higgs strings and semi-local strings [75–77, 142–149], up to some variations [80, 81, 131, 150]. As a result, extrapolating the following results to other types of string should be made with caution. On the bright side, (22) suggests that as long as the string acceleration remains small compared to , one expects the equations of motion of the string to be close to the NG case (up to the eventual vortons appearance). In the following, we describe the results obtained for NG strings. Some differences exist with the results obtained in the Abelian Higgs string simulations.

Before entering into details, let us summarize two fundamental properties these simulations have revealed. The first is that a cosmic string network avoids cosmological domination by evacuating most of its excess energy through some complex mechanisms, which typically result in transferring energy between the horizon-sized distances and the smaller length scales. For NG simulations, this is the formation of cosmic string loops whereas in Abelian Higgs simulation boson radiation is involved. The second property is that the influence of the initial conditions is expected to disappear on the length scales of astrophysical interests. A network of cosmic strings relaxes towards a cosmological attractor which depends only on the expansion rate; this is the so-called scaling regime.

##### 3.1. Dynamics

The equations of motion for NG strings are the vanishing of the extrinsic curvature vector . In a flat FLRW background, Equation (25) can be simplified with the transverse gauge fixing conditionswith the notation and for the timelike and spacelike string coordinates. Such a choice of coordinates reflects the property that an NG string is Lorentz invariant along the worldsheet; there is no physical longitudinal component of the string velocity. In this gauge, the equations of motion readwhere a “dot” and a “prime” stand, respectively, for differentiation with respect to and . We have also defined the quantityThe conditions in (28) do not completely fix the coordinate degrees of freedom and one can supplement them with the so-called temporal gauge fixing which identifies the timelike coordinate with the background time at the string event: . In the transverse temporal gauge, (28) reads , while (29) simplifies towithand is the conformal Hubble parameter. The vector symbols being understood as three-dimensional spatial vectors. Numerically, it is much more convenient to solve an equivalent set of equations found by Bennett and Bouchet [144]. Defining the new vectors and asevaluated at the new coordinates and , the equations of motion (31) can be recast intoAs an illustrative example, these equations have an exact solution in Minkowski space. Taking , one immediately gets (up to a normalisation constant), and are constant over the characteristics and . Inverting (33) giveswhich describes the propagation of left and right moving string deformations at the speed of light. In the FLRW background, these modes are no longer free moving, but interact through the Hubble term (see (34)). Solving these equations gives the for each strings but does not predict what happens when two strings collide.

##### 3.2. Collisions

In the infinitely thin approach, the outcome of an NG string intersection event cannot be predicted. Stress tensor conservation equations require either that the two strings pass through each others, or they intercommute as sketched in Figure 3. The outcome of a string collision process can only be addressed within the framework of a microscopic model. Numerical simulations of interactions have been performed for a variety of models, and in particular for the Abelian Higgs string as represented in Figure 4. In this case, unless the relative string velocity is close to unity [151, 152], or the strings are almost parallel, string intercommutation generically occurs [153, 154]. Let us notice that for type I Abelian strings (i.e., having ), bound states of the two strings can also be formed at low velocity [155–157].

**(a)**

**(b)**

**(a)**

**(b)**

**(c)**

**(d)**

The situation is not necessarily the same for the other types of string. As already mentioned, non-Abelian strings can weave new strings from their intersection points while current carrying cosmic strings intercommute for bosonic carriers [158]. In the case of cosmic superstrings, it has been shown that they intercommute with a probability depending on the fundamental string coupling, a quantity which can be significantly smaller than unity [159]. In the case of -string collisions, Y-junctions can be formed under some kinematic constraints [126, 157, 160–163]. Concerning NG simulations, string collisions are actually implemented through a phenomenological probability of intercommutation at each intersection event.

##### 3.3. Initial Conditions

Solving the cosmological evolution of a NG string network amounts to solving (34) and (35) along each string, finding all of their intersection points and implementing an intercommutation, or not, with the probability . The network evolution is now uniquely determined once the initial conditions are specified. The simplest way to set initial conditions is through the Vachaspati-Vilenkin (VV) algorithm [164]. Motivated by the Kibble mechanism, one assumes a Higgs field to be uncorrelated above a given correlation length . A cosmic string will cross a given plane if one can find a closed loop along which its phase runs from to a multiple of . On a discrete three-dimensional lattice, of -spacing, it is sufficient to approximate by and randomly choose the phase at each corner from three values , and to decide if a string crosses the associated face. Other symmetry breaking schemes and lattice can be approximated in a similar way [165–169]. In Figure 5, we have shown the initial string network configuration obtained from the VV algorithm. The string paths have been smoothed by replacing the right angles by circles of radius . The initial network configuration obviously depends on the physical parameter , the network correlation length. In FLRW spacetime, there is, however, another physical parameter which has to be specified; the distance to the horizon . From those, the initial string energy density is now uniquely determined. A random transverse velocity field can also be added along each string since one does not expect the strings to be initially at rest in any realistic setup. At this point, let us mention that the numerical implementation of the VV initial conditions introduce two additional purely numerical parameters: the size of the periodic box which contains the simulation, usually normalised to unity in comoving coordinates, and the discretisation step required to represent a string, usually given by , the number of points per correlation length.

##### 3.4. Cosmological Scaling

###### 3.4.1. Long Strings

By switching on the evolution from the initial network, string motion and intersections drastically change the shape of the strings as well as the network aspect (see Figure 6). Naively, without any collisional process, one would expect the string network to dominate the energy density of the universe. In a volume , denoting by the typical correlation length of the network at a given time (initially ), the number of strings should be roughly given by . The resulting energy density should therefore beDue to cosmological expansion one has and . As noted by Kibble, this domination does not occur due to intercommutation processes which allow the formation of loops. In the so-called “one scale model”, Kibble [3] assumes that loops of typical size are formed at a rate equals to (for relativistic speeds, one expects one intercommutation per string during the time ). As a result, during a time interval , the energy density transferred to loops isFrom (37), the energy density of strings which are not loops verifieswhere is a function of the cosmic time. Defining , this equation can be recast intowhere the background cosmological fluid sourcing the universe expansion has an equation of state . The constant solution is an attractor for whichThe energy density associated with strings which are not loops “scales” as matter in the matter era and radiation in the radiation era. Therefore, it is prevented to dominate over the usual cosmological fluids and cannot overclose the universe. Notice however that the total energy density could still dominate the cosmological dynamics if the energy density under the form of loops is not evacuated by some extramechanism. For NG cosmic strings, loops are transformed into radiation due to the emission of gravitational waves [170–173]. Other types of loops may lose energy by different radiative processes, such as particle emission, or even energy leakage into the extra-dimensions in the case of cosmic superstrings [174, 175]. In Figure 6, we have represented an evolved string network at the end of a matter era run. For FLRW simulations within a fixed comoving box with periodic boundary conditions, one cannot evolve the system indefinitely; at some point, periodic boundaries become causally connected. Usually, one stops the run when the distance to the horizon fills the simulation volume, or more rigorously half of it. In this figure, one sees that only a few super-horizon strings remain (black long strings) whereas the box is also filled with a lot of small loops (blue) and a few larger loops having a size typical of distance between two long strings. The latter are freshly formed Kibble loops whereas the existence of the small ones cannot be explained in the framework of the one scale model [144]. Concerning the energy density associated with the super-horizon strings (also called infinite strings), their evolution in the matter and radiation era have been plotted in Figure 7; they “scale” as expected. From [146], one hasAs the behaviour of the energy density associated with long strings suggests, the time evolution drives the string network towards a stable cosmological configuration which does not seem to depend on its initial configuration, at least for the long strings. Figure 6, therefore, displays what a cosmological string network should look like inside a horizon volume, at any time during the matter era. As is clear from Figure 7, the relaxation time required for the energy density of long strings to reach the attractor is small. Concerning the cosmic string loops, their existence and behaviour have been the subject of various claims and analytical works [146, 148, 176–181]. In the following, we present recent results [146, 182] showing that the energy density of loops also reaches a scaling evolution similar to (42).

**(a)**

**(b)**

###### 3.4.2. Loops

As previously mentioned, the small loops observed in NG simulation cannot be explained in the framework of the one scale model. These loops find their origin from the self-intercommutation of strings on length scales typical of their small size. The building of a small scale structure on strings is the outcome of the successive intersection events during which new kinks are produced and propagate along the intercommuted segments. Correlations between the kinks induce, from (33), auto- and cross-correlations between and from which small loops can be produced [176, 179]. In Figure 8, we have plotted the energy density distribution under the form of loops with respect to the conformal time during the radiation and matter era. The simulation performed is one of the largest up to date; the box contains whereas the redshift simulation range reaches almost two orders of magnitude. The loop energy density distribution is defined such that is the energy density carried by all loops having a physical length in the range to . In other words, we measure loop size in unit of the horizon length. (This is the relevant physical length scale of the problem.) A logarithmic binning in of resolution has been used in the range to compute these quantities. From this plot, it is clear that after an overproduction regime characterized by the bump of Figure 8, the energy density distribution of loops of given size relaxes towards a stationary regime in which it scales as . Such an observation implies that, once relaxed, the loop number density distribution is of the formwhere the “scaling function” is found to be well fitted by the power laws [146] withfor the matter and radiation era, respectively. The loop number density distribution, for the matter era run, has been plotted in Figure 9. As Figure 8 already shows, the loop distribution takes more time to reach the scaling regime for the small loops. The relaxation bump is all the more so high and long than is small. In the loop number density distribution, this effect appears as a minimal time decreasing value such that the loop distribution is in scaling at . The redshift range probes by a FLRW string simulation is typically , while for strings formed at the GUT energy scales, one expects a at nucleosynthesis. It is clear that, in the cosmological context, the string network has quite a time to relax; on all of the relevant observable length scales the loop distribution should be in scaling, that is, . Since a power law distribution is scale-free, one concludes that a cosmologically stable string network* does not* exhibit loops of a particular size; this is not surprising since the only length scale involved is the distance to the horizon. These numerical results can be analytically recovered in the framework of the Polchinski-Rocha model [179]. The expected loop number density distribution have been explicitly derived by Rocha in [182] with a predicted power for the matter era and in the radiation era. If not due to statistical errors, these small differences may be explained by the existence of additional fractal microstructure along the strings not considered in the analytical approach (see also Section 4). Of course, the previous statements hold provided the other physical effects which are not included in the simulation do not enter the game, as gravitational radiation and gravitational backreaction. The typical length scales at which they should play a role is typically a multiple of , or some positive power of it [183, 184]. As shown in [182], gravitational radiation indeed cures the energy density divergence that one can extrapolate from Figure 9 when . Let us finally notice that although the long strings are defined by , there is also a small population of Kibble loops. Their typical size being the horizon-sized correlation length of the long string network, they can be defined to be those having , where

**(a)**

**(b)**

**(a)**

**(b)**

##### 3.5. Relaxation Towards Scaling: Memory of The Initial Conditions

Although of less-cosmological interest, the relaxation of the loop energy density distribution towards its scaling regime shows interesting properties which could explain some of the differences observed between Abelian Higgs field simulation and NG simulations. In the left panel of Figure 10, we have plotted the loop energy density distribution for loops smaller than the ones in scaling. At the end of the numerical simulation, these length scales are those having . This plot shows that the formation of the smallest loops is a delayed mechanism which suggests that a cascading process takes place from the initial string network configuration. For , only the increasing part of the relaxation bump appears at the end of the simulation whereas the decaying towards scaling is still visible for the larger . On the right panel, we have plotted the loop number density distribution at the end of the simulation, that is, at the time in the left panel. The loop number density distribution deviates from the asymptotic scaling distribution on two typical length scales.

Firstly, an overdensity compared to scaling is situated at . In other words, there is an overproduction of loops with a typical size equal to the initial correlation length of the string network. Although one expects the system to retain some memory of the initial conditions during the relaxation, it may appear surprising that, in spite of the expansion of the universe, the physical length scale of these loops remains the same. A physical interpretation is that which suggests that, at those small length scales, the system decouples from the Hubble flow. More quantitatively, this effect can be explained in the context of the three scale models [176]. Under some assumptions, the string small-scale correlations can indeed sustain a constant physical length.

Then, there is the overall peak of the loop number density distribution at . Most of the loops which are not in scaling have this size at the end of the simulation. In fact, one can check that these loops start appearing soon after the beginning of the string evolution. This length scale is, again, at a constant physical length and is associated with a purely numerical effect [146]. As discussed, a numerical string is discretised with points. The Bennett-Bouchet code at the basis of the simulations presented here uses an adaptive griding algorithm meaning that loops of any physical size can be formed [144]. The only restriction is that, at a given time, a loop is an object of at least three points. Consequently, when the initial string network starts its evolution, loops smaller than cannot be formed. The existence of a finite numerical resolution therefore adds some unwanted initial correlations of length . Notice that this is not a cutoff but indeed an extra-correlation.

As a basic consequence, one should not trust an NG simulation at those length scales. However, the fact that the initial string network violently relaxes towards scaling by emitting loops at the smallest available correlation lengths has still some physical significance [185]. What happens if we increase ? As discussed in [146], the larger length scales are not affected and only the overall peak is shifted around the new . At these length scales, it is clear that using an NG string to describe a network of topological defects would break down and a reasonable assumption is to assume that the network will now relax by losing energy through the relevant physical mechanism available at those distances. In Abelian Higgs simulation, most of the network energy is emitted through field radiation, up to the point that almost no loop are observed in the simulations [75, 77]. Abelian simulation suffers from low resolution compared to NG ones and this has been a subject of debate to decide whether or not this could explain the absence of loops [76, 77]. The above NG results clearly support that particle and/or gravitational waves; emission is an important mechanism which certainly dominate the relaxation regime. However, when the scaling regime progressively takes place, from large to small length scales, the loop formation mechanism should become dominant. One may speculate that it is not clearly observed in Abelian simulations due to its delayed appearance, as it is actually the case in NG simulations when we compare it to the formation of -sized loops. (A hint for this is the compatibility of the Abelian Higgs string correlators with the Polchinski-Rocha model [77]; this one explaining the NG loop distribution [182].)

#### 4. Cosmic Microwave Background Anisotropies and Non-Gaussianities

At this point, numerical simulations of cosmic string evolution give us the means to derive observable predictions. As should be clear from the previous section, some structures in the numerical simulations are not supposed to be present after a realistic expansion factor of ; these are the structures still in the relaxation regime and such that , at any simulation time. In the following, we will denote by “systematic errors”, the uncertainties inherent to the presence of such nonscaling structures when deriving observable predictions from NG numerical simulations.

##### 4.1. Unequal Time Correlators

The first method used to derive CMB anisotropies has been introduced in [186] and applied in [79, 186–191] for global topological defects and recently in [192, 193] for the Abelian strings.

Cosmic strings are active sources of gravitational perturbations [41, 195] which means that the equations of motion of their induced linear perturbation is of the formwhere is a time differential operator, the perturbation in the energy density, or velocity, and so forth, which is directly related to the CMB temperature anisotropies. Here denotes the source terms, that is, the string stress tensor. From the Green’s function of this equation, one gets, today (at ) and in Fourier spaceThe two-point correlator readsand its determination requires a full-time knowledge of the source term for each mode. Since it is impossible to carry out a simulation over the whole cosmological history, the scaling properties of the cosmic string network can be used to analytically extrapolate the source terms over the required ranges. As shown in [189], as long as the cosmic string network is in a scaling regime, the source terms are the stress tensor components and assume the formNumerical simulations are actually used to determine the scaling functions . In Figure 11, we have represented the temperature anisotropies derived in [194] using such a method from Abelian Higgs string simulations. The current CMB constraint on comes from this power spectrum: at two-sigma, [196].

##### 4.2. Simulated Small Angle Maps

The previous constraint typically corresponds to a string contribution which cannot exceed at the multipole moment . On current observable angular scales, cosmic strings may only be a subdominant fraction of the overall CMB anisotropies. However, string-induced perturbations being non-Gaussian, as opposed to inflationary perturbations of quantum origin, one can go further than deriving the two-point function. Notice that, in principle, the unequal time correlator approach could be used to extrapolate the three- and higher -point function by using the scaling properties of the string network. Another approach is to produce simulated maps of string-induced CMB anisotropies. Again, we face the problem of the small redshift range probed by the numerical simulations. By putting an observer inside the numerical simulation, such maps can only include stringy effects up to a finite redshift, typically . The CMB anisotropies computed in this way are therefore only accurate on large angular scales but can produce full sky maps [198–201]. This limitation can be avoided by stacking maps from different redshifts, an approach outlined in [37] and applied in [197, 202].

Simulations with the observer outside of the numerical box are not well suited for a full-sky map reconstruction, but are perfectly designed for the small angular scales. The reason being that cosmic strings are incessantly sourcing the CMB fluctuations since the last scattering surface, and contrary to the perturbations of inflationary origin, this part cannot be affected by Silk damping. Therefore, at small angular scales, one expects the strings’ signature in the CMB temperature fluctuations to be dominated by their integrated Sachs-Wolfe (ISW) effect from the last scattering surface [202]. In the temporal gauge (), the NG stress tensor derived from (26) readsIn the flat sky approximation, well suited for angles typically smaller than the Hubble angular size at the epoch of interest, Hindmarsh has shown that the ISW temperature anisotropies induced by NG strings can be simplified to [203, 204]where . The wave vector denotes the* transverse* component of the three-dimensional vector with respect to the line of sight , whereas, in the temporal gauge, encodes the string stress tensor distortions of the photon temperature and readsAs can be seen in (51), only the strings that intercept the photon path can imprint their signature in the CMB temperature fluctuations. The previous expression is nothing but the Gott-Kaiser-Stebbins effect in the temporal gauge [54, 203, 205]. As a result, the knowledge of , and therefore of the string trajectories , is only required on our past light cone. In the context of string numerical simulations, the trajectories of all strings are computed during all of the numerical simulation time. Therefore, to compute , one only needs to determine which parts of the string network intercept our past light cone and at what time.

In order to be able to generate a significant amount of maps using such a method, it is more convenient to stack two medium sized simulations along the lines sketched in Figure 12. The first one starts at the last scattering surface and ends at a redshift fixed by the maximum expansion factor achievable in the numerical box. For the simulations we performed, initially (in unit of the fixed comoving simulation volume), which corresponds to and a field of view of (for the current fiducial cosmological parameters [44]). Such a run ends after a 30-fold increase in expansion factor, corresponding to a redshift . We then propagate the photons perturbed by the first run into a second numerical simulation of the same size but starting at . For another 30-fold increase in expansion factor, this run ends at . As can be seen in Figure 12, the second simulation represents a much larger real volume than the first one and therefore subtends a greater angle in the sky. As a result, only the subpart of the second run that matches the angle subtended by the first simulation is actually used. As we will see later on, the CMB temperature maps are weakly sensitive to the string network at low redshifts, simply because there are almost no strings intercepting our past light cone in a recent past, which makes this technique perfectly acceptable. In practice, each of these numerical simulations is started before the redshifts mentioned, in order to give the cosmic string network enough time to relax toward its stable cosmological configuration. As discussed in the previous section, one has to make sure that the structures (strings and loops) we are interested in have indeed reached their scaling behavior during the numerical runs. We switch on the photon propagation inside the runs only after making sure all the large structures (infinite strings and loops) are in their scaling regime. This can be checked by monitoring the evolution of the energy density distributions, and we have chosen to start the photons’ propagation when all loops larger than a third of the horizon size are in scaling. This cutoff is then dynamically pushed to smaller values to include all the loops entering the scaling regime at later times. The cutoff time dependence is simply the function and can be deduced from the loop distribution relaxation times derived in Section 3.4.2. The resulting CMB temperature map is displayed in the left panel of Figure 13 whereas the right panel shows the string paths projected onto our past light cone. Again, these maps are only representative at small angular scales. On larger angles, they represent only the ISW contribution to the overall string anisotropies; for instance, Doppler contributions coming from photon decoupling at last scattering are dominant around (see Figure 11).

**(a)**

**(b)**

The discussion on systematic effects coming from the presence of loops not yet in scaling can be found in [46]. In fact, they have only a small effect. The physical reason being that, due to scaling, the long strings are still the main source of CMB anisotropies even at (reasonably) small angles. Indeed, there are always roughly ten strings per Hubble volume at each time, which means that a patch of is at least crossed by the ten long strings being there at last scattering, plus a few others from lower redshifts.

##### 4.3. Skewness and Kurtosis

The most basic statistical test that can be performed from a set of small angle CMB maps is to plot the one-point function of the temperature anisotropies. As can be seen in Figure 14, the temperature anisotropies induced by cosmic strings are clearly non-Gaussian. From a set of independent CMB maps, one finds the mean sample skewness to be negativewhere the brackets stand for the mean over different realisations while the bar denotes averaging on each map. The variance itself averages toThe quoted errors are statistical and refer to the square root of the variance between the different realisations. Similarly, the mean kurtosis averages toAn analytical approach extending these results to cosmic superstrings can be found in [206]. A simple way to look for strings is to search for large (but rare) temperature fluctuations. Deviations from Gaussianity start to be significant, let us say by a factor of two, only in the tails when the probability distribution becomes typically lower than .

##### 4.4. Real Space Methods

Strings induce step-like discontinuities in the CMB anisotropies and various methods have been designed to probe the non-Gaussianities associated with them.

Multifractal analysis [207] has the advantage of being directly applicable to the time-ordered data retrieved when a CMB telescope scans the sky. As opposed to the fractal dimension of a set which measures how sparse it is, the multifractal spectrum of a measure defined over a set gives how many and which fractal dimensions there are. In the context of cosmic strings, this method has been applied in [208] on one dimensional scan of maps similar to the one in Figure 13, the measure being defined bywhere the integer labels a point along the scan. The multifractal properties of this measure have been shown to be distinctive enough to detect strings, compared to a Gaussian signal, but only when the detector resolution is sufficiently good. One may wonder, under multifractality, how a non-Gaussian string pattern could be distinguished from other non-Gaussian sources. In fact, (56) consists in taking the gradient of the induced CMB fluctuations along the scan. Step-like discontinuities, passed over a gradient filter, become one-dimensional delta functions, and this is a definite string feature that can only be altered by the beam experiment. Denoting by and the horizontal and vertical angular coordinates, the gradient magnitude of the temperature anisotropies is defined byThis definition makes it clear that for a finite temperature step, let us say , being the Heaviside function, the resulting gradient magnitude is a Dirac distribution at the string location. In Figure 15, we have plotted the gradient magnitude of the temperature maps of Figure 13 as well as its convolved version with a Gaussian beam typical of the Planck satellite at . With a finite resolution beam, the discontinuities are now smoothed. Real space methods applied to string are therefore strongly sensitive to the angular resolution. Let us mention that wavelet analysis methods have been also explored in this context [209, 210] or to produce cleaner maps [211].

**(a)**

**(b)**

Directional gradients, obtained by variations with respect to either or , have been discussed in [205] in the context of Minkowski functionals. They are again found to provide a more distinctive non-Gaussian behaviour than the original temperature map for the contour length and genus.

##### 4.5. Temperature Power Spectrum

Moving to Fourier space, the small angle CMB maps also permit a determination of the power spectrum at large multipoles. In Figure 16, we have plotted its mean value over the different maps as well as the one-sigma statistical error around the mean. The overall power at is [46]which is close to the value obtained in Abelian Higgs field simulation (see Figure 11). This is not so surprising since the long strings in both NG and Abelian Higgs simulation have a similar scaling evolution, and as explained above, long strings are the main sources of CMB anisotropies even at the small angles. The power law tail in Figure 16 is the direct consequence of the presence of strings at all times since the last scattering surface: one finds for [46]where only the systematic errors have been reported. Such a power law shows that cosmic strings have to become the dominant primary source of CMB anisotropies at the small angular scales, the fluctuations of inflationary origin being killed by Silk damping at those multipoles. In Figure 17, we have plotted the respective contributions of strings and adiabatic anisotropies of inflationary origin. The cosmological parameters have been set to their fiducial values in the Lambda-Cold-Dark-Matter (LCDM) model and the string energy density is compatible with the current bounds. For the current upper limit on , CMB anisotropies should become dominated by cosmic strings at . An unresolved Sunyaev-Zel’dovich (SZ) component may, however, compromise such a signal. Nevertheless, string induced anisotropies do not depend on the signal frequency whereas the SZ does, and one may hope to disentangle both [46].

##### 4.6. Hindmarsh Approximation

As shown by Hindmarsh in [203], the power law behaviour of the NG string power spectrum at small angles can be analytically recovered. In Fourier space, the power spectrum of the string-induced CMB anisotropies is defined bythe expression of being given by (51). Remembering that the string positions and velocity vectors have to be evaluated on the past light cone, it is more convenient to use the so-called light cone gauge. Instead of identifying the timelike worldsheet coordinate with the background time at the string event, one chooses instead to identify . In this gauge, (51) simplifies towhere we have definedand where the capital indices are two-dimensional. The time parameter then labels the intersections of a set of null hyperplanes with the worldsheet. For our problem, all quantities have to be evaluated at . In a field of view of formal area , one can express the power spectrum asAdding the assumptions that both and obey Gaussian statistics, all of the correlation functions of can now be written in terms of two-point functions only. Using the same notation as in [203], the nonvanishing two-point functions areas well as the quantitiesThe leading terms are given by [203, 212]where we have definedThe correlation length is the projected correlation length on the past light cone, is the mean square projected tangent vector, is the mean square projected velocity and the correlation between projected velocity and curvature. From these assumptions, (63) reduces toWhen gets large, the terms involving the mixed correlator can be shown to be subdominant and only the first term remainswhere . Denoting by the total transverse light-cone gauge length of string in the box of area , one getsAt small angles, the wave number and (74) predicts that . The small difference with (59) is suggestive of a cloud of zero-dimensional objects along the string worldsheet which may be the signature of small loop production in the NG numerical simulations. Let us stress that (74) is not “primordial” but directly approximate the observed angular power spectrum of the CMB temperature anisotropies.

##### 4.7. Bispectrum

###### 4.7.1. Analytical Approach

The success of Hindmarsh approximation to describe the small angular CMB anisotropies power spectrum suggests it can be applied to higher -point functions. In [212], this method was used to derive the bispectrum defined from the three points function byPlugging (61) into the previous expression giveswith , , and . With the Gaussian assumption, the ensemble average of the string observables is lengthy but straightforward and the final result reads [212]The quantities are shorthand for the scalar products . In the same way as for the power spectrum, this expression directly gives the bispectrum of the CMB temperature anisotropies. Its overall dependence varies as . Its sign depends on the sign of defined in (71), and contrary to what one could naively expect : the projected string velocity and curvature vectors are correlated. This can be shown by starting again from the equations of motion (29), but this time, in the light cone gauge. The equation of motion for gives whereas the equation for the transverse components isIn a FLRW background, assuming that is constant, and neglecting higher-order correlations between , , and , we findwhere we have defined , and where is the averaged conformal Hubble parameter. Still assuming that the ensemble is approximately Gaussian in and , the right hand side reduces toThe last term vanishes and the cross correlator is positive; from (71), we deduce that . It is interesting to notice that would vanish in Minkowski spacetime, which can be viewed as a consequence of time reversal invariance. The existence of a cosmic string bispectrum is the consequence of the breaking of the time reversal invariance in a FLRW background.

An illustrative example is to apply (77) to the isosceles triangle configurations in Fourier space such thatwhere denotes the angle between the wavevectors and . The isosceles bispectrum reads

Notice that for , we obtain the peculiar case of an equilateral triangle. In Figure 18, we have plotted the angle dependency of the isosceles bispectrum. These configurations are amplified as in the two flat triangle limits for which either or . Both of these configurations are therefore better suited than the equilateral one to characterize the strings. As suggested by the real space searches, the strings can produce a strong bispectrum signal only if the detector resolution is sufficiently good. Assuming a beam resolution of means that the field of view would contain at maximum roughly Fourier modes. Consequently, the smallest values of achievable would be around radians, with only a few modes saturating this bound.

###### 4.7.2. Numerical Results

The previous analytical results can be compared to the CMB temperature bispectrum derived from the simulated maps of Section 4.2. Numerically, one can use the scale convolution method introduced in [213, 214] and applied to the string bispectrum in [212]. This method relies on the choice of unity window functions in Fourier space peaked around a particular wavenumber . Definingone can construct an estimator of the three point function in Fourier space by remarking thatFor thin enough window functions, remains constant over the window function width and we construct our reduced bispectrum estimator asThe function is the flat sky equivalent of the inverse Wigner-3 symbols and readsFor the window functions such that for , one can approximatefor small enough values of compared to the wavenumber . In this case, (87) can be worked out into

In the left panel of Figure 19, we have plotted the mean string bispectrum and its standard deviation obtained by this method over the string CMB maps. For illustration purpose, this plot is for the isosceles configuration having radians. The right frame of Figure 19 shows the same mean bispectrum but multiplied by , for various small values of . As expected from the analytical results, we recover the behaviour. The wavenumber dependency also matches with the analytical calculations, up to similar slight power differences as we found for the power spectrum. A power law fit against the mean numerical estimator giveswhile the overall amplitude can be evaluated around the minimum variance multipole. At , one getswhich also matches with (83) under some crude estimation of the string parameters [212]. Finally, as suggested by Figure 18, the string bispectrum is mostly negative. Integrated over all possible configurations, one recovers the mean negative sample skewness of (53), thereby explaining its origin as a direct consequence of the breaking of the time reversal symmetry in FLRW spacetimes.

**(a)**

**(b)**

##### 4.8. Trispectrum

The trispectrum of the string induced CMB temperature anisotropies can be derived in a similar way. Starting from the definition of the four-point functionswe define the trispectrum aswith , and . (Notice that our denomination “trispectrum” here contains the unconnected part. This one is however non-vanishing only for parallelogram configurations of the wavevectors.) As shown in [215], the trispectrum and the higher -point functions exhibit unfactorable “flat directions” in the -dimensional space of the integration variables . Physically, it means that the leading order part of the (connected) trispectrum is sensitive to the higher orders of the correlators in (64) to (66). For the trispectrum, the correlator has to be expanded at next-to-leading order, and following the Polchinski-Rocha model [179], we assume a nonanalytical behaviour for at small scalesIn the light-cone gauge, we leave and as undetermined parameters since they cannot be straightforwardly inferred from the numerics performed in the temporal gauge. Nevertheless, because the correlation should be smaller as becomes larger, one should have . In the temporal gauge, is directly related to the power law exponent of the scaling loop distribution functions through . As we are bound to show, the mode dependence of the trispectrum will also be uniquely given by this parameter. Once the tangent vector correlator expressed as in (94), the integrations in (93) can be performed explicitly, except for parallelogram configurations which have to be dealt with as a special case. After some tedious calculations, an interpolating expression for the trispectrum is [215]In this equation, is a number depending only on the parameter while is the trispectrum geometrical factor defined bywhereand stands for cyclic permutations over the indices. The function stands for the normalised incomplete lower gamma function defined byand, finally, has been defined by

As an application, the trispectrum over parallelogram configurations is obtained when and the leading term of (95) simplifies towhere now refers to the parallelogram angle. Under the scaling transformation , the parallelogram trispectrum scales asFor parallelograms, as already mentioned, the trispectrum also gets a contribution from the unconnected part of the four-point function, which is Gaussian and readsUsing (63), ones sees that the unconnected part also behaves as . Therefore, the non-Gaussian contributions for parallelogram configurations remain of the same order of magnitude as the Gaussian ones, with however, and again, an exception in the squeezed limit .

The most interesting situations come when . For these quadrilaterals, the Gaussian contribution vanishes and solely a non-Gaussian statistics can source the trispectrum. At large wavenumber (small angles), one has such that the normalised lower incomplete gamma function in (95) is close to unity:Under the scaling transformation , the geometric factor, and thus (104), scales aswithAs claimed, for NG strings, is directly given by the power law of the loop distribution [146]. Since this exponent is different from the one associated with parallelogram configurations it may actually be used to distinguish a trispectrum sourced by cosmic strings with the one generated by other non-Gaussian effects.

In Figure 20, we have plotted the geometrical factor for the kite quadrilaterals (represented in the same figure), as a function of their opening angles and .

**(a)**

**(b)**

As for the bispectrum, the trispectrum is enhanced on the squeezed configuration obtained when