Alignment of Quasar Polarizations on Large Scales Explained by Warped Cosmic Strings . PART II : The Second Order Contribution

We find an azimuthal-angle dependent approximate wave like solution to second order on a warped five-dimensional manifold with a self-gravitating U(1) scalar gauge field (cosmic string) on the brane using the multiple-scale method. The spectrum of the several orders of approximation show maxima of the energy distribution dependent on the azimuthal-angle and the winding numbers n of the subsequent orders of scalar field. This breakup of the quantized flux quanta does not lead to instability of the asymptotic wavelike solution, due to the suppression of the n-dependency in the energy momentum tensor components by the warp factor. This effect is triggered by the contribution of the five dimensional Weyl tensor on the brane. This contribution can be understood as dark energy and can trigger the self-acceleration of the universe without the need of a cosmological constant. There is a striking relation between the symmetry breaking of the Higgs field described by the winding number and the SO(2) breaking of the axially symmetric configuration into a discrete subgroup of rotations about 180 . The discrete sequence of non-axially symmetric deviations, cancelled by the emission of gravitational waves in order to restore the SO(2) symmetry, triggers the pressure zz T for discrete values of the azimuthal-angle. There can be a possible relation between the recently discovered angle-preferences of polarization axes of quasars on large scales and our theoretical predicted angle-dependency and can be an evidence for the existence of cosmic strings. The discovery of the increase of polarization rate in smaller subgroups of the several largequasar groups (LQGs), the red shift dependency and the relative orientation of the spin axes with respect to the major axes of their host LQGs, point at a fractional azimuthal structure, were also found in our cosmic string model. This peculiar discontinuous large scale structure, i.e., polarizations directions of multiples of, for example, π 2 or π 4 , can be explained by the spectrum of How to cite this paper: Slagter, R.J. (2017) Alignment of Quasar Polarizations on Large Scales Explained by Warped Cosmic Strings. PART II: The Second Order Contribution. Journal of Modern Physics, 8, 163-180. https://doi.org/10.4236/jmp.2017.82015 Received: January 17, 2017 Accepted: February 14, 2017 Published: February 17, 2017 Copyright © 2017 by author and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access

quantized flux quanta does not lead to instability of the asymptotic wavelike solution, due to the suppression of the n-dependency in the energy momentum tensor components by the warp factor.This effect is triggered by the contribution of the five dimensional Weyl tensor on the brane.This contribution can be understood as dark energy and can trigger the self-acceleration of the universe without the need of a cosmological constant.There is a striking relation between the symmetry breaking of the Higgs field described by the winding number and the SO(2) breaking of the axially symmetric configuration into a discrete subgroup of rotations about 180  .The discrete sequence of non-axially symmetric deviations, cancelled by the emission of gravitational waves in order to restore the SO(2) symmetry, triggers the pressure zz T for discrete values of the azimuthal-angle.There can be a possible relation between the recently discovered angle-preferences of polarization axes of quasars on large scales and our theoretical predicted angle-dependency and can be an evidence for the existence of cosmic strings.The discovery of the increase of polarization rate in smaller subgroups of the several largequasar groups (LQGs), the red shift dependency and the relative orientation of the spin axes with respect to the major axes of their host LQGs, point at a fractional azimuthal structure, were also found in our cosmic string model.This peculiar discontinuous large scale structure, i.e., polarizations directions of multiples of, for example, π 2 or π 4 , can be explained by the spectrum of

Introduction
General relativity theory (GRT) is by far the most successful theory constructed by theoretical physicists.Its predictive power is impressive.Famous empirical confirmed examples are the Kerr black hole and the emission of gravitational waves by merging black holes.There are, however, predicted phenomena not yet detected by observations.An example is the axially symmetric (spinning) compact object, i.e., the Papapetrou or Lewis-van Stockum solution.Another well studied object is the self-gravitating cosmic string solution.For an overview, see the work of Vilenkin and Shellard [1] and Anderson [2].Cosmic strings are U (1) scalar gauge vortex solutions in general relativity in the framework of GUT's.
This U(1) scalar gauge field with a "Mexican hat" potential has lived up its reputation in the theory of superconductivity, where vortex lines occur as topological defects and in the standard model of particle physics.In cosmology it could trigger the inflationary period of expansion and could solve the horizon and flatness problem.It came as a surprise that the relativistic string-like vortex solution of Nielsen and Olesen [3] can be found in GRT [4] [5].However, general relativistic cosmic string is still not found directly or indirectly by observations.The interest in cosmic strings faded away when one found inconsistencies with the power spectrum of the CMB: Cosmic strings cannot provide satisfactory explanation for the magnitude of the initial density perturbations from which galaxies and clusters grew.It turns out that the upper bound of the mass per unit length . Further, the special pattern of a lensing effect of cosmic strings is not found yet.Studies of the radiative effects of cosmic strings embedded in a FLRW spacetime show that cylindrical gravitational radiation is rapidly damped and is negligible in any physical regime [6].These string-cosmology spacetimes essentially look like a scaled version of a string in a vacuum spacetime.There is, however, another possibility to test the existence of cosmic strings.The recently discovered alignment of quasar polarizations on very large scales [7] [8] could be explained by considering cosmic strings on a warped brane world spacetime [9] [10].It was realized recently that cosmic strings could be produced within the framework of superstring theory inspired cosmological models and a revival of cosmic strings occurred.These so-called cosmic superstrings can play the role of cosmic strings in the framework of string theory or M-theory, i.e., brane world models.Supersymmetric GUT's can even demand the existence of cosmic strings.Super-massive strings with an energy density of 1 Gµ  are interesting [11] because their gravitational impact will be much stronger than GUT strings.They could be produced when the universe underwent phase transitions at energies much higher than the GUT scale.Brane world cosmological models were first proposed by Arkani-Hamed, Dimopoulos and Dvali (ADD) [12] [13] and Antoniadis et al. [14] and extended by Randall and Sundrum (RS) [15] [16].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 electro-weak scale, (hierarchy problem) 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 warp factor.The weakness of gravity in these models are fundamental and the Planck energy could be of ( ) that may be accessible by LHC.It is possible that effective 4D Kaluza-Klein (KK)-modes are obtained from the perturbative 5D graviton.These KK-modes will be massive from the brane viewpoint.Further, one usually considers a fine-tuning between the tension on the brane (4D cosmological constant) and the 5D tension in order to ensure a zero effective cosmological constant.For an overview of brane world models, see, for example, Roy [17] [18] [19] and Shiromizu et al. [20] and references therein.Cosmic strings could have tremendous mass in the bulk, while their warped manifestations in the brane show consistency with the observed bound of by the warp factor [21], even if its value was at the Planck scale.Wavelike disturbances triggered by the huge mass of the cosmic string in the bulk, could have indeed observational effects in the brane.Evidence of these objects would give us information at very high energies in the early stages of the universe.Maybe they might actually provide us the best observational window upon fundamental string theory.One reason for this conjecture is the possibility of the extension of fundamental strings into the bulk in warped spacetimes.
In this research we extend our previous work on the relation of the alignment of quasar polarizations with warped cosmic strings [10].Recent research in this field [22] shows peculiar correlation of quasar polarization orientations in LQGs.The alignments can be parallel or perpendicular to each other.Ever other discretizations are not excluded.We shall see that this behaviour fits in our brane world model.However, there are some peculiar issues to be addressed.A network of long cosmic strings can be characterised by a single scaling length, the persistence length of the inter-string distance.Numerical simulations have shown that the evolution of the network must be scale-invariant.It is a great challenge to prove that there is a comparable scale-invariance in the alignments of the polarization vectors of quasars.By the discovery of high-redshift supermassive black holes, 6 z > , one could extend the investigation on quasar alignment to this era.There is already a tiny indication of z-dependency in the observed quasar polarization orientations [7].Further, the formation of smallscale closed loops of cosmic strings can disturb the energy balance of long strings.In our model this problem is overcome by the warp factor: it makes the cosmic string more massive during the evolution.
In Section 2 we outline the multiple-scale method on a warped brane world spacetime.In Section 3 we calculate the metric perturbations to second order.In Section 4 we derive the matter field equations to second order and indicate the possible relation with axially symmetric instabilities caused by radiationreactions.In Section 5 we discuss the possible connection of the warp factor with conformal invariance.In the appendices we collected all the relevant equations in order to keep the main text readable and self-contained.

The Multiple-Scale Approximation on a Warped Brane World Spacetime
We will investigate vortex-like solution on a warped five-dimensional Friedmann-Lematre-Robertson-Walker (FLRW) model in cylindrical coordinates [9] [10]: with  a warpfactor dependent of , r t and the bulk dimension y.The selfgravitating scalar-gauge field, parameterized as resides on the brane.η is the vacuum expectation value of the scalar field, n the winding number and  the gauge coupling constant.The winding number (number of jumps in phase of the scalar field when one goes around the flux tube) is related to the quantized flux in the Ginsberg Landau theory of superconductivity (Abrikosov vortices) and the discrete values of the topological charge in the sin-Gordon theory.The exact solution of  [9] follows from the 5D Einstein equation ( ) M is the fundamental 5D Planck mass.The scalar-gauge field equations become [4] ( ) ( ) ∇ the covariant derivative with respect to 4 g µν and V the potential of the Abelian Higgs model.The star represents the complex conjugated.F µν is the Maxwell tensor.The modified Einstein eq- uations become [20]  and n µ the unit vector normal to the brane.We will consider here eff 0 Λ = , so we are dealing with the RS-fine tuning condition [15].The last two terms on the righthand side of Equation ( 5) represent the quadratic contribution of the energy-momentum tensor and the electric part of the five dimensional Weyl tensor respectively.All standard model fields are bound to the brane, although gravity may propagate into the extra dimension.Is is obvious, that the cosmic string can build up a huge mass 1 Gµ  by the warp factor and can induce massive KK-modes felt on the brane.The warp factor causes perturbations to be damped as they move away from the brane, so gravity looks four dimensional, at least perturbatively, to a brane world observer.Brane world models can also explain the acceleration of the universe without the need of a cosmological constant [18].Disturbances can survive the natural damping by expansion of the universe due to the warp factor.This effect was also found numerically [9].
Here we will consider the modified cosmic string features on the warped spacetime Equation ( 1) and use the multiple-scale approximation [23] [24] [25] in order to find partial differential equations (PDE's) for the perturbations to second order of the metric and matter fields.This extends our previous research [10].
Let us expand the metric field and the scalar-gauge fields in the multiple-scale scheme ( , , ( , , , , with g µν the background metric and , Φ A µ the background scalar and gauge fields.For the scalar field we take different winding numbers, so different magnetic flux quantization for the background field and higher order perturbations.We define , e , , , e , , e in Z t r ϕ ξ Ξ = . So we break up the original vortex with winding number n in our case in three strings with winding numbers 1 2 , n n and 3 n .One can prove [26] that this breakup remains stable if the gauge to scalar mass is > 1.In our case stability will be guaranteed by inverse powers of the warp factor.Further, , which will fulfil the highest order perturbation equation of the gauge field as we shall see.Rapid perturbations occur in the direction of the wave vector l x The Einstein equations in subsequent orders of approximation become We will consider here the equations up to order ( ) ω .We used the notation for the several terms in the expansion of a tensor, vector or scalar: The contribution from the bulk space, µν  , must be calculated with the 5D Riemann tensor.If we consider 0 l l µ µ = , i.e., the Eikonal equation (which follows from the ( ) 1 ω − scalar equation), then one obtains from Equation ( 8) a set of restrictions on h µν , such as the "gauge condition" . Then we let survive 11 13 14 44 , , , h h h h and 55 h as independent first order perturbations of the metric.

The Metric Perturbations up to Second Order
The PDE's for the background fields 1 , , W ψ γ and the first order perturbations h µν can be derived from Equation ( 9)and were found in [10].By integrating Equation ( 10) with respect to ξ , one obtains with τ the period of the rapid variation.The last two terms on the right hand side of Equation ( 11) can be interpreted as "back-reaction" contribution of the KK-modes of the perturbative 5D graviton and can act as a cosmological constant.These equations can be used to eliminate the background fields from the first order PDE's.One then obtains propagation equations for the first order perturbations.For completeness we collected them in appendix A. From these equations one observes that to first order there is an interaction between the high-frequency perturbations from the bulk, the matter fields on the brane and the evolution of ij h  , also found in the numerical solution [9].The bulk contribution 55 h  is amplified by (also found in [10]).
Let us now investigate the higher order equations in ω , i.e., Equation (10), which will provide us first order equations of t k µν ∂  and second order eq- uations for tt h µν ∂ . With the help of an algebraic manipulation program, one obtains, for example, the equations for  ( )  , also a ( )  -term in connection with the second order perturbation Z  , amplified by the warp factor.In the equation for the first order counterpart equation, i.e., Equation (24), there is only the ( ) , then the maxima in ϕ of these two terms belonging to the perturbations of first and second order respectively, are out-of-phase.In the next section this will also become clear by considering the energy-current components of the energy momentum tensor.From Equation (33) and Equation (26) we can obtain a second order PDE for 55 h if we impose constraint conditions on 55 k  or integrate Equation (33) with respect to ξ .This can also be done for the other components.So we can construct with the multiple-scale method uniformly valid wavelike approximations to solutions of perturbation problems without without "resonance" interactions between the consecutive orders of perturbation theory.This result is related to the Cauchy problem.In any field theory where there is a gauge freedom (as, for example, in GRT and Maxwell theory), one has to specify gauge conditions in order to determine the dynamical evolution for some initial set of Cauchy data.In Maxwell theory one usually choose the Lorentz gauge.In GRT one has constraint equations because the system is over determined.These constraints are usually PDE's of first order.In our approximation scheme we have Equation ( 8) which leads to conditions on h µν  .The first order equations Equations ( 23)-( 27) can also be considered as constraint equations for the second order wave equations for h µν , as example, in the case of 55 h .So we can construct a dynamical evolution of the system of equations which fulfil the Cauchy data.
We shall see in the next sections how this result can be applied in context with the recently observed alignment of the polarization vectors of quasars over large distances and the discreteness in the azimuthal dependency of the polarization axes.

The Matter Field Equations and the Energy-Momentum Tensor Expansion
From Equation (4) we obtain, after substituting the expansions of Equation ( 6), from the ( ) where we have integrated the equation with respect to ξ .On the right hand side we see again the high-frequency contribution to the field equation.In our case, this back reaction term turns out to be zero.So the first order equation is just the unperturbed equation for X .See Equation (28) of Appendix A. This also holds for the equation for A µ .See Equation (29) of appendix A. If we substitute back the integrated equations into the original equations, we then obtain the first order perturbations (for ). See Equation (30)-Equation (32) of Appendix A. Further, we used 0 l A α α = , otherwise the real and imaginary parts of Ψ  interact as the propagation progresses.Again, there appears a ϕ -dependent term in the propagation equation for 0 B , amplified by 1 W .This deviation from axially symmetry was also found by [23].More insight in this ϕ -dependency can be obtained by studying the second-order matter field equations.One obtains for Z  and C  again first order differential equations.See Appendix B. We used the fact that the complex conjugate of the full complex second order scalar equation also must be satisfied.The appearance of the terms ( ) ( ) ( ) where  is an expression in backgrounds fields, 11 h and 44 h .We have again two periodic functions ( ) ( )  with fre- quency difference of a factor two and where one of the functions is amplified by W .It will be necessary to study these equations numerically in order to compare the amplitudes of these two periodic functions with those of the azimuthal-angle dependent maxima of the quasar polarization alignment.
The second order equation for 0 B can be obtained from the sum of the t- and r-components of the second order gauge field equations.
We can calculate the three first terms of the energy momentum tensor n n P − + .In the next order ( ) tt T there will be terms proportional with higher orders of ( ) n n P − + .
These higher order terms will be suppressed by the warp factor, so the vortex will not become unstable as is the case when one breakup the vortex string in multiple flux [27].
The most interesting behavior arises in the angular component There is a relation between the phase freedom e inϕ of our scalar field and the secular instability of an initially quasi-stationary axially symmetric configuration caused by radiative reaction.The small non-axially symmetric deformations turn out to be of the form e imϕ with m an integer [28].This broken symmetry, described by the inverse of the angular momentum J, is comparable with the symmetry breaking of the Higgs field considered in our model.An axially symmetric systems is invariant under rotations in two dimensions, the SO(2) group.The breaking of this symmetry can be expressed in the equatorial eccentricity.The particular orientation of the ellipsoid in the (x-y) frame can be expressed through the azimuthal angle ϕ .This discrete change into non-axially symmetry must be cancelled by emission of gravitational energy (and is amplified in our model by the 5D contribution), otherwise we are saddled with a helical time coordinate, t t Jϕ → + and must give up Lorentz invariance.This is clear from the fact that our metric will then possess a t g ϕ term.The angular momentum in (x-y) plane is determined by the currents of the momentum density, x T x T ρ νµ ν ρµ ∼ − and can be calculated in our case with the off-diagonal components of T µν of Appendix C. For example ( )

Quasar Polarization Alignment and Scale-Invariance
In order to explain the recently found large-scale alignments of the polarization vectors of quasars in LQGs at cosmological redshifts 1.5 z  by cosmic strings, it would be desirable to find a kind of scale (conformal)-invariance, because it is conjectured that cosmic string networks evolve scale-invariant just after the radiation dominated era of our universe.The intricate features of the polarization axes alignments of the quasars in LQGs show at least a kind of coevolution at very large scales [22], so a study of conformal symmetry of our model could be of interest.
Gravity theory invariant under → Ω is local conformal invariant and must be spontaneously broken because our world appears not to be scale invariant [29] [30].Let us rewrite our spacetime of Equation ( 1) with where we renamed the warp factor as Ω .If we consider the ( ) , r t -dependent part of Ω and consider (brane) case of the metric Equation (1), we then obtain for the Ricci scalar ( ) The Ricci scalar transforms under → Ω as [31] ( ) So for conformal invariancy of 5 R , the second term on the right hand side of Equation ( 17) must vanish.For 5 ˆ0 R = we then find 1 0, with constraint equations 18) is just the equation of a vibrating circular drum.The general solution for the boundary conditions scalar gauge field in the brane, one can try to formulate again the conformal invariance.This is a peculiar issue in theoretical physics till now.Einstein equations and the scalar equation (Klein-Gordon equation)) are not conformally invariant.One has to modify Einstein's equations to make it conformally invariant and make the energy momentum tensor traceless [31] [32].Our Ωfield can play a crucial role in this context if one introduces an unavoidable dilaton field.

Conclusions
It is found on a five dimensional warped brane world spacetime, using a multiple-scale approximation scheme, that to second order the metric and scalar gauge field show a spectrum of azimuthal-angle dependent wavelike modes with extremal values dependent of the winding numbers of the background, first and second order perturbations of the scalar field.
In four dimensional models, this local field theory admits vortex-like behavior and is a generalization of the Ginzburg-Landau theory of superconductivity.A lattice of Abrikosov vortices can be formed, carrying a quantized flux dependent of the winding number or "vortex charge" n .Vortices with 1 n > are unlikely, since the energy is reduced if they split up into single vortices.However, in the general relativistic case, gravity comes into play and the energy of the configuration of the vortices must be calculated covariantly by means of the energy momentum tensor.This general relativistic vortex solution (cosmic string) can build up a huge mass per unit length in the bulk and can induce massive Kaluza-Klein modes felt on the brane, where the standard model fields reside.
Disturbances don't fade away during the expansion of the universe due to the warp factor.
The jump in the phase of the scalar field is related to the secular instability of the initially stationary axially symmetric configuration caused by the radiation reaction.The breaking of the axially symmetry, described by the inverse of the angular momentum, is e imϕ  (m an integer and ϕ the azimuthal angle), comparable with the symmetry breaking of the scalar field.The recovery of the SO(2) symmetry from the equatorial eccentricity is triggered by the emission of gravitational waves.Our model can be used to explain the mysterious alignment of quasar axes with the large-scale structure of our universe and can serve as evidence for the existence of cosmic strings.The found fractional azimuthal structure is comparable with the angle-preferences of the polarization axes of quasars.
There is a strong evidence of scaling (conformal) behaviour of long nonintercommuting cosmic strings networks during the radiation-dominated era.
High or low initial string densities tend toward a fixed scaling value.However, standard cosmology, constraints such strings to be very light and will fade away (or disappear by the forming of closed loops).In our model they can survive by the warp factor.This fact makes the comparison with the alignment of quasars possible.So it would be desirable to have also a scale-invariant alignment ) ( ) ( ) The background matter fields ( , X P ) become [9]   ( ) and perturbation equations ( ) These background equations don't contain perturbations terms, due to our (simplified) gauge conditions.
with  a function of the background fields and the fields ij h .The other components of the metric perturbations are obtained in the same way.The second order equations for Z  and C  are ( )

C. The Energy Momentum Components
For the several orders of the energy-momentum tensor components we find 4 0, )

4 Λ
the brane tension and

4 G
µν the Einstein tensor calculated on the brane metric

7 )
Substituting the expansions into the effective Einstein and matter field equations, one obtains a set of PDE's for , , g h k µν µν µν and the scalar gauge fields , , , , A B µ µ Φ Ψ Ξ and C µ .The perturbations can be ϕ -dependent.

2 1 W 1 W
. It is a reflection of the massive KK modes felt on the brane.The most interesting equation is the differential equation for 14 h  , i.e., the ( ) , t ϕ component, Equation (24).It triggers the ϕ - , can have extremal values on [ ] 0, π , if we choose, for example, then have the term cos 2ϕ , which has two extremal values on [ ] otherwise Ψ must be zero.From the first order equation of the scalar field we obtain ( ) ( )

with 0 J , 0 Y
Bessel functions and n C coefficients dependent of ( ) f r .These solutions represent for suitable boundary conditions, the standing normal modes of the brane in the vacuum case.In general, Ω can also depend on the azimuthal angle.This dependency is found, in our non-vacuum situation, in the preceding sections in the multiple-scale approximation.So one could conclude that the warp factor in the vacuum case fulfils a scalar wave equations representing fluctuations of the brane in the ground state.It represents the amount of local "stretching" of the 4D geometry.In the non-vacuum case, with the U(1) in Equation (35) and Equation (36) respectively.They can be isolated in order to get a wave equation for the first order perturbations.
next manuscript to SCIRP and we will provide best service for you:Accepting pre-submission inquiries through Email, Facebook, LinkedIn, Twitter, etc.A wide selection of journals (inclusive of 9 subjects, more than 200 journals) Providing 24-hour high-quality service User-friendly online submission system Fair and swift peer-review system Efficient typesetting and proofreading procedure Display of the result of downloads and visits, as well as the number of cited articles Maximum dissemination of your research work Submit your manuscript at: http://papersubmission.scirp.org/Or contact jmp@scirp.org 4T µν .See Appendix C. In 4T ϕϕ , i.e., can change sign dynamically.A numerical solution can give a decisive answer. ) ) )