A New Fate of a Warped 5D FRW Model with a U(1) Scalar Gauge Field
Abstract
If we live on the weak brane with zero effective cosmological constant in a warped 5D bulk spacetime, gravitational waves and brane fluctuations can be generated by a part of the 5D Weyl tensor and carries information of the gravitational field outside the brane. We consider on a cylindrical symmetric warped FLRW background a U(1) selfgravitating scalar field coupled to a gauge field without bulk matter. It turns out that brane fluctuations can be formed dynamically, due to the modified energymomentum tensor components of the scalargauge field (“cosmic string”). As a result, we find that the latetime behavior could be significantly deviate from the standard evolution of the universe. The effect is triggered by the timedependent warpfactor with two branches of the form
( with constants) and the modified brane equations comparable with a dark energy effect. This is a braneworld mechanism, not present in standard 4D FLRW, where the large disturbances are rapidly damped as the expansion proceed. Because gravity can propagate in the bulk, the cosmic string can build up a huge angle deficit (or mass per unit length) by the warpfactor and can induce massive KKmodes felt on the brane. Disturbances in the spatial components of the stressenergy tensor cause cylindrical symmetric waves, amplified due to the presence of the bulk space and warpfactor. They could survive the natural damping due to the expansion of the universe. It turns out that one of the metric components becomes singular at the moment the warp factor develops a extremum. This behavior could have influence on the possibility of a transition from acceleration to deceleration or vice versa.

06.20.Jr, 95.30.Dr, 95.30.Sf, 98.62.Ra, 98.80.k, 98.80.Es, 98.80.Jk
1 Introduction
It is conjectured that the expansion of our universe is accelerating. Recent observations provide strong evidence of this acceleration. The explanation of this remarkable phenomenon is rather difficult: one needs a dark energy field with an effectively negative pressure, . From several independent observational data, one finds and , where and stand for the energy densities of matter and dark energy respectively with respect to the critical ones [1, 2, 3, 4, 5]. We should live now in the cosmological constant dominated era (and approximately flat, ), while at earlier times there was a radiation and matter dominated era.
There are some fundamental questions. First, will the acceleration last forever. Secondly, why there is the huge discrepancy between in Planck units and vacuum energy density which is times greater than the value we need to accelerate the expansion of the universe. Thirdly, the value must be incredibly finetuned, . So it would be a logical step to try to explain the latetime acceleration without the need for dark energy [6]. It is also conjectured that one needs an inflaton field in the very early stage of our universe, to solve the flatness and horizon problem in the standard model of cosmology. This is the inflationary cold dark matter model with cosmological constant (CDM). It could also predict the existence of fluctuations we observe in the CMB shortly before the end of inflation. The inflatonfield could be the wellknown scalarHiggs field with the mexicanhat potential. This model has lived up to its reputation. It was successful in the explanation of superconductivity, i.e., the GinsburgLandau theory, in the standard model of particle physics, in the general relativistic solution of the selfgravitating NielsenOlesen vortex (cosmic string) [7, 8, 9] and could play a fundamental role in warped spacetimes [10, 11, 12]. Cosmic strings are U(1) scalargauge vortex solutions in general relativity (GR)[9, 13, 14]. It is conjectured that in any field theory which admits cosmic strings, a network of strings inevitably forms at some point during the early universe. However it is doubtful if they persist to the present time. Evidence of these objects would give us information at very high energies in the early stages of the universe. From recent observations by COBE, WMAP and Planck satellites, it was concluded that cosmic strings cannot provide satisfactory explanation for the magnitude of the initial density perturbations from which galaxies and clusters grew. The interest in cosmic strings faded away, mainly because of the inconsistencies with the power spectrum of CMB.
Cosmological cosmic strings can also be investigated on a FriedmannLemaitreRobertsonWalker (FLRW) background. However, the stringcosmology spacetime essentially looks like a scaled version of a string in a vacuum spacetime and the corrections in the field equations due to cylindrical gravitational radiation are rapidly damped and are negligible in any physical regime [15]. The reason of this result comes from the notion that the string radius is much smaller than the Hubble radius in the postinflationary era. The ratio , where is the factored out overall scale factor in the metric, is then of the order . The resulting field equations look like the generalized NielsenOlesen vortex system. Only in the preinflationary epoch, radiative corrections could be very large. While the inflaton field plays a crucial role in the early stage of our universe, it could play a comparable role at much later times, if we modify gravity by considering warped braneworld models. It could be possible that there exists a correlation between the accelerating universe and large extra dimensions in braneworld models.
When it was realized that cosmic strings could be produced within the framework of superstring theory inspired cosmological models, a revival of cosmic strings occurred. These socalled cosmic superstrings can play the role of cosmic strings in the framework of string theory or Mtheory, i.e., braneworld models. Supersymmetric GUT’s can even demand the existence of cosmic strings. Physicists speculate that extra spatial dimensions could exist in addition to our ordinary 4dimensional spacetime. The idea that spacetime could have more than four dimensions was first proposed by Kaluza and Klein (KK) in the early 20 century [16, 17]. These theories can be used to explain several of the shortcomings of the Standard Model, i.e., the unknown origin of dark energy and dark matter and the weakness of gravity (hierarchy problem). In these models, the weakness of gravity might be fundamental. One might naively imagine that these extra dimensions must be very small, i.e., curled up and never observable. Supermassive strings with , could be produced when the universe underwent phase transitions at energies much higher than the GUT scale. Recently there is growing interest in the socalled braneworld models, first proposed by ArkaniHamed, Dimopoulos and Dvali (ADD) [18, 19] and I. Antoniadis et al. [20], which was extended by Randall and Sundrum (RS) [21]. In these models, the extra dimension can be very large compared to the ones predicted in string theory, i.e., of order of millimeters. The difference with the standard superstring model is that the compactification rely on the curvature of the bulk. The huge discrepancy between the electroweak scale, and the gravitational mass scale will be suppressed by the volume of the extra dimension , or the curvature in that region. This effect can also be achieved in the RS models by a warpfactor. In the RS2 model, there are two branes, the visible and the gravity brane at . The branes have equal and opposite tensions . At low energy, a negative bulk cosmological constant will prevent gravity to leak into the extra dimension, , with the corresponding energy scale and the compactification radius. The squeeze the gravitational field closer to the brane at . In the RS–1 model, one pushes the negative tension brane . If one finetunes the , then this ensures a zero effective cosmological constant on the brane. The infinite extra dimension makes, however, a finite contribution to the 5D volume due to the warpfactor. Because of the finite separation of the branes in the RS2 model, one obtains socalled effective 4D modes (KKmodes) of the perturbative 5D graviton on the 4D brane. These KKmodes will be massive from the brane viewpoint. In the RS1 model, the discrete spectrum disappears and will form a continuous spectrum [22, 23, 38].
It will be clear that compact objects, such as black holes and cosmic strings, could have tremendous mass in the bulk, while their warped manifestations in the brane can be consistent with observations. So, braneworld models could overcome the observational bounds one encounters in cosmic string models. could be warped down to GUT scale, even if its value was at the Planck scale. Although static solutions of the U(1) gauge string on a warped spacetime show significant deviation from the classical solution in 4D [25], one is interested in the dynamical evolution of the effective brane equations. Wavelike disturbances triggered by the huge mass of the cosmic string in the bulk, could have observational effects in the brane. One conjectures that these disturbances could act as an effective dark energy field. The question is if the pulselike cylindrical waves have the desired asymptotic behavior at null infinity.
Here, we will investigate the latetime evolution of a warped 5D FLRW model when a U(1) scalar field coupled to a gauge field is present. In section II, we outline the model and in section III, we present some numerical solutions of the model. Finally, we have summarized our results in section IV.
2 The Field Equations
In order to embed a cosmological cosmic string in a FLRW model, one has to modify the metric. The spacetime of a cosmic string is invariant under boosts in the direction, whereas the FLRW spacetime is not. Gregory[15] pointed out that in the static case, if the radius of the string core is small compared to the Hubble radius , the embedding is easily done. For a nonstatic radiating cylindrically symmetric spacetime, one can analyze the behavior of gravitational waves generated by the local strings on the expansion of the universe if we consider the zerothickness limit of nonsingular spacetimes containing a cylindrical distribution of stressenergy embedded on a cosmological background.
The importance of cylindrical symmetric gravitational waves was first noticed by Beck [27] and rediscovered by Einstein and Rosen (ER)[28]. A substantial part of the present knowledge on gravitational waves originated from their wavelike solutions of the fully nonlinear gravitational field equations. However, these original ER cylindrical wave solutions of the Einstein’s vacuum field equations inhabit a universe that is conical rather than flat at infinite cylindrical radius[29]. To overcome this problem, one could impose artificially nonconicality constraints[30, 31] in order to obtain a Kasnerlike spacetime. In a warped 5dimensional braneworld model with a U(1) gauge cosmic string, these constraints are superfluous because the effective brane spacetime is nonconical[10]. Much insight into wavelike solutions on a cylindrically spacetime can also be obtained from well known stationary axially symmetric solutions, like the Weyl, Papapetrou, Lewisvan Stockum, and Kerr solution, through a complex transformation [32]. The cylindrical symmetric metric is (with two killing vectors , both spacelike and hypersurface orthogonal)
(1) 
Without any matter source, one then obtains a wave equation for , which is decoupled from . Once is solved, can be solved by quadratures and the behavior of and at null infinity is well known, i.e., the EinsteinRosen pulse waves. One can construct incoming and outgoing radiation solutions [33, 34] in terms of the retarded and advanced coordinates and . For a nonvacuum situation, the equations will not decouple.
It was found by Gregory [15] that a U(1) cosmic string can be embedded into a flat FLRW spacetime along the polar axis, when the Hubble radius is much larger than the stringcore. The approximated spacetime is
(2) 
with the mass per unit length of the string, or angle deficit and and when respectively. One can match this spacetime on the wellknown FLRW spacetime by a transformation from to by defining and , with when or 1 respectively [35]. One should also like to impose matching conditions across the boundary of the interior and exterior spacetime (although for radiating strings, this will be difficult). For latetime solutions, it is not necessary to consider the matching conditions between the interior string solutions and the exterior FLRW spacetime. It turns out that when the width of the strings is smaller than the Hubble radius, the disturbances are neglectable [15].
On a warped spacetime, where the warpfactor has a significant effect on the cosmic string solution [10, 11], disturbances can be significant. In the bulk, the mass of the cosmic string will be huge, inducing a back reaction on the brane equations.
2.1 The Bulk Equations
Following Shiromizu et al. [38], we have 5dimensional Einstein equations with a bulk cosmological constant and brane energymomentum as source
(3) 
with , the brane tension, , , and the unit normal to the brane. The is the fundamental 5D Planck mass, which is much smaller than the effective Planck mass on the brane, GeV. We consider here the matter field confined to the brane, i.e., the U(1) scalargauge field, written, in polar coordinates , in the form [9]
(4) 
with , the vacuum expectation value of the scalar field, the coupling constant and the radial coordinate. Further, we take for the potential of the scalar field the wellknown “Mexicanhat” form . Let us consider the cylindrically symmetric warped spacetime
(5) 
with be the warpfactor and be the bulk space coordinate.
From the 5D equations we obtain from the component of the 5D Einstein equations (a dot means )
(6) 
If we write as particular solution of Eq. (6) and substitute this expression into the 5D Einstein equations, we can then separate the ydependent parts containing :
(7) 
with a simplified solution (for the moment, )
(8) 
We see that only a negative bulk cosmological constant make sense. For one obtains a solution with positive bulk cosmological constant. The yindependent part of the 5D Einstein equation for becomes
(9) 
with a prime . A typical solution is
(10) 
with some constants. So we have two branches, i.e., the plus and minus sign in Eq. (10). Figure 1 represents the typical plots of for different sets of the constants and . In general, can possess a saddlepoint or extremal values, as can be seen from the timedependent part of of Figure 1. We shall see in section (2.2) that the appearance of a minimum ( or maximum in case of the other branch) will make the equation for divergent. This could have a significant influence on a transition from acceleration to deceleration or vice versa. We shall also see in Section 2.2 that these solutions for the warpfactor are consistent with the supplementary equations following from the 5D Einstein and the Bianchi equations. It turns out that the equations for can not be isolated from the effective 4D brane equations, indicating that is a really a warpfactoreffect.
2.2 The Effective Brane Equations
The field equations induced on the brane can be derived using the Gauss–Codazzi equations together with the Israel–Darmois junction conditions at the brane and the symmetry [36]. The modified Einstein equations become
(11) 
where and is the vacuum energy in the brane (brane tension). The latter equality sign is a consequence of the relation between the 4 and 5dimensional Planck mass in the braneworld approach, [37, 38]. If in addition the brane tension is related to the 5dimensional coupling constant and the cosmological constant by , then and we are dealing with the RSfine tuning condition[21]
The first correction term is the quadratic term in the energymomentum tensor arising from the extrinsic curvature terms in the projected Einstein tensor
(12) 
The second correction term is given by
(13) 
and is a part of the 5D Weyl tensor and carries information of the gravitational field outside the brane and is constrained by the motion of the matter on the brane, i.e., the Codazzi equation
(14) 
Further, we have for the extrinsic curvature from the junction conditions
(15) 
and the 4D contracted Bianchi equations
(16) 
From the 5D Einstein and Bianchi equations one obtains supplementary equations
(17) 
(18)  
(19) 
(20) 
with the Lie derivative with respect to and the “magnetic” part of the bulk Weyl tensor, . From these supplementary equations one can also separate the equations in Eq. (7) and in Eq. (9) for and which proves that the supplementary equations are consistent with the bulk equations. From the 4D scalargauge field equations we obtain
(21) 
(22) 
which are consistent with those obtained by Gregory [15]. From the Einstein equations we obtain
(23)  
(24)  
(25)  
(26) 
(27)  
(28)  
(29) 
In the 4D case[15], assuming , one can factor out a timedependent “scale”factor in the 4D FLRW spacetime in the metric and terms in the field equations of the form , could safely omitted. The overall average cosmological expansion is not affected by the string. In our case, this is not allowed. One proves[15] that in the 4D case the ratio . In our warped 5D case one cannot compare with . It is determined by Eq.(10) and will depend on the parameters of the cosmic string, for example the mass per unit length . By this warpfactor , one obtains a huge angle deficit, or equivalently, a tremendous mass per unit length , while the warped manifestation in the brane will be warped down to GUT scale, consistent with observations.
The time dependent part of the warpfactor causes disturbances of the order much larger than the expected values in the 4D case. We observe that the several terms on the righthand side of the metric field equation Eq. (29) are not influenced by the sign of the two branches of the warpfactor. On the righthand side of Eq. (26) we have in the denominator the overall term . So for an extremum of , the solution diverges. This could have a significant influence on a transition from acceleration to a deceleration or vice versa. A reflecting point [ see figure 1] does not cause this problem.
3 Numerical Solution
We solved numerically the set of PDE’s and plotted in figure 2 a typical solution for the choice of .
We used valueboundary conditions at and Neumann boundary conditions at . For the initial values we choose: (WeberWheeler pulse wave), , and . It is evident that the emission of retarded wave energy of the scalargauge field, triggered by , changes the evolution of the induced brane metric. This is a nonlocal effect from the bulk. The relevant components of the energymomentum tensor are and :
(30) 
(31) 
We observe that there are two terms which will become dominant for different epochs of cosmic time, i.e., the potential term , acting as a tension and the gauge field term . They both contain the term in the numerator and denominator respectively. The same arguments hold for the quadratic term . The pure gauge field term in and will oscillate in time but can change sign in , so alternate between pressure (plus sign) and tension (minus sign). In the static case, and will only contain the term and is positive, meaning that work is required to squeeze the lines of the magnetic field toward the axis in the radial direction. In the case of a negative sign one would need work in order to stretch the lines of the magnetic field. From the numerical solution of figure 2 we observe that can become negative. In the 4D case[9], this change of sign of is initiated by the change of ( the gauge to scalar mass ratio) from a value to a value . Here we find a dynamical change. The same arguments hold for the pure scalar term . However, it doesn’t contain the warpfactor.
In figure 2 we also plotted the Cenergy (abbreviation for “cylindrical” energy introduced by Thorne[39]), the BrownYork quasilocal Cenergy [41] and the advanced and retarded Cenergy flow along the two nulldirections, () and (). The Cenergy was originally introduced by Thorne as a very useful tool in analyzing the dynamics of a wholecylindersymmetric system (i.e., cylindrical symmetry with hypersurface orthogonal Killing fields) and to describe the nature of EinsteinRosen gravitational waves and electromagnetic waves. It was useful to demonstrate the resistance of magneticfield lines to cylindrically gravitational collapse and to discuss the dynamics of Melvin’s magnetic universe when large radial perturbations are introduced[40]. Evidence for a bona fide “energy” status of is the fact that its coordinate gradient is related to a flux vector field which is covariantly conserved[41]. It was found that the BrownYork quasilocal Cenergy description is a better measure of the Cenergy than the original one of Thorne. is derived by means of a Hamiltonian reduction of the field equations, wherein the energy arises naturally from the Hamiltonian constraint, while Thorne’s arguments were merely heuristic. Another argument for the choice of is the fact that the asymptotic metric in the vacuum case can be written in the form
(32) 
with the asymptotic value for . This can be seen by the rescaling , and in the asymptotic form of metric of (5). This metric is exactly the metric exterior to a cosmic string with mass per unit length . This is only valid in the weakfield limit . Because doesn’t decouple from the equation for , we have two degrees of freedom, i.e., and , for the gravitational waves (in the EinsteinRosen case there is one degree of freedom, because can be obtained by quadrature when is known). Further we have the warpfactor as extra metric component. It turns out that the wavelike behavior depends critically on the ratio of the scalar to gauge masses, . We took . For one obtains the features of a global string, which is in general singular[11]. From the plots of figure 2 we can draw some global conclusions concerning regular behavior, the interaction of de cosmic string with the gravitational waves and the deviation from the 4D FLRW counterpart model ( with the same matter field and cylindrical symmetry). The Cenergy is radiated along the null directions in agreement with other research on cylindrical symmetric spacetimes. Although our numerical code runs out of accuracy for larger r and t values, we observe a falloff of the WeberWheeler pulse . The induced disturbances in the scalar field X travel backwards closer to the axes of the cosmic string and the scalar field returns to its original static form with its core closer to r=0.
Some remarks must be made with regards to the asymptotic behavior of the metric Eq.(32). If we define , then the metric Eq.(32) becomes Minkowski. However, now takes the values . So if , we obtain a angle deficit and the spacetime is conical. Because is related to the metric component [41], we can plot this angle deficit residue and compare its behavior with the 4D counterpart model of Gregory[15]. In the static 4D case, so in the absence of radiation, the angle deficit for isolated strings at large radial distances is related to the mass per unit length of the cosmic string, i.e., . In the weak field approximation one has . So for GUT strings () we expect that the asymptotic behavior is practically Minkowski minus a wedge. The angle deficit will also depend on . A typical value for the angle deficit is for [13]. If one admits radiative effects on the 4D timedependent spacetime, one proves[15] that the stringcosmology spacetime essentially looks like a scaled version of a cosmic string solution, not desirable (see Eq. (32)) because we live now in a flat FLRW spacetime. In deriving this result, one makes the simplification that the gravitational wave energy decouples from the cosmic string energy. The energy of the string is then exponentially damped and one finds that asymptotically remains constant for all time. In our 5D warped model, we made no simplifications, because we solved the full set of coupled wave equations numerically. Moreover, we have the contribution from the bulk as explained in section (2.2). We observe from figure 2 that increases with time. One can prove[14], that for increasing , the angle deficit becomes and one deals with a Kasnerlike spacetime. So we conclude that in our warped cosmological model we don’t deal with a scaled version of the standard cosmic string. The notion of a conical spacetime fades away on our effective spacetime. For values of above the GUT scale, gravitational effects become more relevant and the weakfield approximation can no longer be applied. These socalled supermassive cosmic strings possess in 4D, however, a curvature singularity[14]. Moreover, they probably would have formed before an inflationary era and subsequently inflated away. In our warped 5D spacetime these supermassive cosmic strings arise in a different manner, i.e., by the warpfactor. We can also compare the plots of and of figure 2, in order to draw some conclusions about the asymptotic conical behavior. The metric component approaches in the classical 4D cosmic string solution to the value , so in the weak field approximation one expects asymptotically that approaches a constant value and will show a small deviation from . This constant value is dependent of the parameters and of the model. However not in our case. We observe that the warped counterpart metric component behavior is different: it increases for larger values of r and t and seems to be insensitive for the values and . So we observe no conical behavior. This absence of the angle deficit in a 5 dimensional warped cosmic string model was also found earlier[10]. If one investigate the brane evolution via the modified second Friedmann equations, then on the righthand side new terms will appear which contain , the quadratic contribution of the energymomentum tensor ( see Eq.(26) and (29)) and a bulk radiation term containing the constants and ( see Eq.(7) and Eq.(10)). From the behavior of of figure 1 and the occurrence of in the denominator of the evolution equation Eq.(26)of one could conclude that an ever lasting exponential expansion is preferable over the other possibilities. The dependency of our model on the several (integration) constants is currently under investigation by the authors.
4 Conclusion
The latetime epoch of accelerated expansion can be explained within the framework of the 4D standard model of cosmology by hypothesizing a dark energy component to the total energy density of the universe with negative pressure equation of state , where the equation of state parameter w lies near 1. Dark energy can possibly be understood as arising from vacuum energy, but this explanation introduces the cosmological constant problem. Alternatively, it is possible that there is no dark energy, but instead a lowenergy/large scale ( i.e., “infrared”) modification to general relativity that accounts for latetime acceleration, schematically, . The higher dimensional graviton has massive 4D modes felt on the brane (KK modes). Gravity leakage from the 4D brane can cause this selfacceleration without the need to introduce dark energy.
We find a significant new latetime behavior of the warped 5D FLRW cosmological model, when a selfgravitating U(1) scalar field coupled to a gauge field (cosmic string) is present on the brane and the effective cosmological constant on the brane is zero (RS finetuning). Compared with the 4D counterpart model, the cylindrically symmetric disturbances have a significant impact on the latetime expansion rate of the universe. This is caused by the fact that one has no bound on the mass per unit length of the cosmic string, in contrast to the 4D model of Gregory [15]. We find an exact solution for the warpfactor, . The timedependent part, which plays the role of a “scale” factor, is a monotone exponential increasing function with an inflection point or has a extremum. In the latter case, the metric component becomes singular at the moment the warp factor develops a extremum. This behavior could have influence on the possibility of a transition from acceleration to deceleration or vice versa. In the first case the accelerated expansion can be accomplished without the need of an effective brane cosmological constant as dark energy. The disturbances of the scalar and gauge fields interact with gravitational waves and could act as an effective dark energy field and trigger brane fluctuations. The scalar field sheds off wave energy and the oscillations slowly decay and eventually the variables return to their static values.
The recently found largescale alignment of quasar polarizations in our Universe[26] could be explained by our model. The modified energy momentum tensor component induces traveling waves in the zdirection on the brane.
References
 Spergel, D. N. et al.: FirstYear WMAP Observations: determination of cosmological parameters. Astrophys. J. Suppl. 148, 175 (2003) [arXiv: astroph/0302209]
 Peiris, H. V. et al.: FirstYear WMAP observations: Implications for inflation. Astrophys. J. Suppl. 148, 213 (2003) [arXiv: astroph/0302225]
 Spergel, D. N. et al.: ThreeYear WMAP observations: Implications for cosmology. Astrophys. J. Suppl. 170 377 (2007) [arXiv: astroph/0603449]
 Perlmutter, S.: et al.: Measurements of and from 42 HighRedshift Supernovae Astrophys. J. 517, 565 (1999) [arXiv: astroph/9812133]
 Riess, A. G. et al.: New Hubble Space Telescope Discoveries of Type Ia Supernovae at : narrowing constraints on the early behavior of dark energy. Astrophys. J. 659, 98 (2007) [arXiv: astroph/0611572]
 Maartens, R.: Dark Energy and Dark Gravity J. Phys. Conf. Ser. 68, 012046 (2007), Dark energy from braneworld gravity. Lect. Notes. Phys. 720, 323 (2007) [arXiv: astroph/0602415v1]
 Felsager, B.: Geometry, Particles and Fields. Odense University Press, Odnese (1981)
 Nielsen, H. B., Olesen, P.: Vortexline models for dual strings. Nucl. Phys. B 61, 45 (1973)
 Garfinkle, D.: General relativistic strings. Phys. Rev. D 32, 1323 (1985)
 Slagter, R. J., Masselink, D.: Warped selfgravitating U(1) gauge cosmic strings in 5D. Int. J. Mod. Phys. D 21, 1250060 (2012) [arXiv:grqc/1202.4487]
 Slagter, R. J.: Time evolution of a warped cosmic string. Int. J. Mod. Phys. D 10, 1237 (2014)
 Slagter, R. J.: Dark energy generated by warped cosmic strings. arXiv: grqc/14077505. (2014)
 LagunaCastillo, P., E. A. Matzner, E. A.: Coupled field solutions for U(1) gauge cosmic strings. Phys. Rev. D 36, 3663 (1987)
 LagunaCastillo, P., Garfinkle, D.: Spacetime of supermassive U(1) gauge cosmic strings. Phys. Rev. D 40, 1011 (1989)
 Gregory, R.: Cosmological cosmic strings. Phys. Rev. 39, 2108 (1989)
 Kaluza, T.: Zum unitts problem in der physik. Sitzungsber. K. Preuss. Akad. Wiss. Berlin, 966 (1921)
 Klein, O.: Quantentheorie und fnfdimensionale relativittstheorie. Zeitsch. Phys. 37, 895 (1926)
 ArkaniHamed, N, Dimopoulos, S., Dvali, G.: The hierarchy problem and new dimensions at a millimeter. Phys. Lett. B 429, 263 (1998) (arXiv: hepph/9803315)
 ArkaniHamed, N., Dimopoulos, S., Dvali, G.: Phenomenology, astrophysics and cosmology of theories with submillimeter dimensions and TeV scale quantum gravity. Phys. Rev. D 59, 086004 (1999) (arXiv: hepph/9807344)
 Antoniadis, I., ArkaniHamed, N., Dimopoulos, S., Dvali, G.: New dimensions at a millimeter to a Fermi and superstrings at a TeV. Phys. Lett. B 436, 257 (1998) [arXiv: hepph/9804398]
 Randall, L., Sundrum, R.: Large mass hierarchy from a small extra dimension. Phys. Rev. Lett. 83, 3370 (1999) [arXiv: hepph/9905221]; An alternative to compactification. Phys. Rev. Lett. 83, 4690 (1999) [arXiv: hepth/9906064]
 Maartens, R.: Braneworld cosmological perturbations: a covariant approach. Prog. Theor. Phys. Suppl. 148, 213 (2003) [arXiv: grqc/0312059]
 Maartens, R., Koyama, K.: Braneworld gravity. Liv. Rev. Relt 13, 5 (2010) [arXiv: grqc/10043962v2]
 Shiromizu, T., Maeda, K., Sasaki, M.: The Einstein equations on the 3brane world. Phys. Rev. D 62, 024012 (2000)
 Slagter, R. J.: Cosmic strings in brane world models. In: Dabrowski, M. P., et al. (ed.) Proceedings of the conference Multiverse and Fundamental Cosmology, pp.3942. Amercan Inst. Phys., NewYork. (2012)
 Hutsemekers, D., Braibant, L., Pelgrims, V., Sluse, D.: Alignment of quasar polarizations with largescale structures. Astron. Astrophys. 572, A 18 (2014) [arXiv: astroph/1409.6098]
 Beck, G.: Zur theorie binarer gravitationsfelder. Z. Physik 33, 713 (1925)
 Einstein, A., Rosen, N. JK.: On gravitational waves. J. Franklin Inst. 223, 43 (1937)
 Weber, J., Wheeler, J. A.: Reality of the cylindrical gravitational waves of Einstein and Rosen. Rev. Mod. Phys. 29, 509 (1957)
 Gowdy, R. A., Edmonds, B. D.: Cylindrical gravitational waves in expanding universes: models for waves from compact sources. Phys. Rev. D75,084011 (2007)
 Gowdy, R. A.: Cylindrical gravitational waves in expanding universes: explicit pulse solutions.arXiv: grqc/07100987v1 (2007)
 H. Stephani, H., Kramer, D., MacCallum, M., Hoenselaers, C., Herlt, E.: In: Landshoff, P. V., et al.(ed.) Exact Solutions of Einstein’s Field Equations. Cambridge University Press, Cambridge, (2009)
 Stachel, J. J.: Cylindrical gravitational waves. J. Math. Phys. 7, 1321 (1968)
 Ashtekar, A., Bicak, J., Schmidt, B. G.: Behavior of EinsteinRosen waves at null infinity. Phys. Rev. D 55, 687 (1997) [arXiv: grqc/9608041v1]
 Anderson, M. R.: In: Foster, B., et al. (ed.) The Mathematical Theory of Cosmic Strings. Institute of Physics Publishing, Bristol, (2003)
 Sasaki, M., Shiromizu, T., Maeda, K.,: Gravity, stability and energy conservation on the RandallSundrum brane world. Phys. Rev. D 62, 024008 (2000) [arXiv: hepth/9912233]
 Cartier, C., Durrer, R., Ruser, M.: On graviton production in braneworld cosmology. Phys. Rev. D 72, 104018 (2005)[arXiv:hepth/0507006v1]
 Shiromizu, T., Maeda, K., Sasaki, M.,: The Einstein equations on the 3brane world. Phys. Rev. D 62, 024012 (2000) [arXiv: grqc/9910076]
 Thorne, K. S.: Energy of infinitely long, cylindrically symmetric systems in general relativity. Phys. Rev 138, B251 (1965)
 Thorne, K. S.: The absolute stability of Melvin’s magnetic universe. Phys. Rev 139, B244 (1965)
 Goncalves, M. C. V.: Unpolarized radiative cylindrical spacetimes: Trapped surfaces and quasilocal energy. Class. Quant. Grav. 20, 3749 (2003) [arXiv: grqc/0212125]
 Slagter, R. J.: Nonlinear gravitational waves as dark energy in warped spacetimes. To appear in the proceedings of the 14th Marcel Grossmann Meeting, Rome (2015)