Integrability of the Einstein-nonlinear $SU(2)$ $\sigma$-model in a nontrivial topological sector

The integrability of the\ $\Lambda-$Einstein-nonlinear $SU(2)$ $\sigma$-model with nonvanishing cosmological charge is studied. We apply the method of singularity analysis of differential equations and we show that the equations for the gravitational field are integrable. The first few terms of the solution are presented.

Nonlinear σ−models are important theoretical models in Physics for the properties that they provide which are of special interest [1,2]. The Einstein nonlinear σ−models, in which the total Action Integral is the sum of the Einstein-Hilbert Action Integral and the Action Integral which corresponds to the nonlinear matter source, have provided different kinds of solutions for the gravitational equations. Specifically it has been shown that there exist black hole solutions with a regular event horizon which asymptotically approach the Schwarzschild spacetime, in the context of the Einstein-Skyrme Model, which violates the "no hair" conjecture for black holes (see for instance [3] and references therein).
The purpose of this work is to study the integrability of the field equations of the Einstein-nonlinear SU (2) σ−model, which have been studied previously in [4]. We do this by using the method of singularity analysis of differential equations 1 . The application of singularity analysis in gravitational theories has been applied by many researchers in the past, for instance in the case of the Mixmaster Universe (Bianchi IX) [5][6][7], in scalar field cosmology [8] and in modified theories of gravity [9][10][11].
Consider a Riemannian manifold M with metric g µν of Lorentzian signature. The action integral of the field equations for the Einstein-nonlinear SU (2) σ−model in a four-dimensional manifold is given by where S EH is the Einstein-Hilbert action with the cosmological constant, i.e., S EH = dx 4 (R − 2Λ), and S (σ) is the action integral of the nonlinear sigma model [15] where U (x ν ) is the SU (2)-valued scalar and K is a positive constant. The physical implication of the action, (2), is that it describes the dynamics of low energy pions. The gravitational field equations are derived by variation of the action integral (1) with respect to the metric tensor g µν . This leads to the following set of equations, in which the left hand side of (3) corresponds to the Einstein-Hilbert action where G µν is the Einstein tensor and Λ is the cosmological constant. The right hand side of (3) is that of the nonlinear σ-model and provides the matter source. The explicit form of the energy-momentum tensor is Furthermore variation with respect to the scalar-valued U (x ν ), in (1) leads to the constraint equation U −1 U ;µ ;ν g µν = 0, while the latter can follow from the application of the Bianchi identity in (3), that is, T µν ;ν = 0. By following the Ansatz, which was proposed in [16][17][18][19] and its generalizations [20][21][22][23], for the parametrization of the SU (2) algebra in [4] Ayón-Beato, Canfora and Zanelli found that for the four-dimensional spacetime, the conservation equation, T µν ;ν = 0, is satisfied always and the energy-momentum tensor, (4), is expressed only in terms of the fields F , ρ and N . Specifically, the matter field equations U −1 U ;µ ;ν g µν = 0, are identically satisfied so that one has only to deal with the Einstein's field equations while the solution provides a nontrivial topological sector [4].
In the minisuperspace approach the gravitational field equations (3) can arise from the Euler-Lagrange equations of the singular Lagrangian, where the equation gives the constraint equation or the G r r − T r r = 0 (8) component of (3). In (6) we can see the dynamical terms which corresponds to the R (3) curvature term of (5) of the cosmological constant and of the σ−model. Because the field equations are singular there could exist a nonlocal conservation law which is generated by the conformal Killing vectors of the minisuperspace, for details see [26,27]. In our consideration, as the minisuperspace of (6) has dimension two, there exists an infinite number of conformal Killing vectors and an infinite number of nonlocal conservation laws. Hence with the use of a nonlocal conservation law the two second-order differential equations and the first-order differential equation, which describe the gravitational field equations, can be reduced to the second-order nonautonomous equation [25], where y = y (r) = ρF, N (r) is in which and the new gauge has been selected to be so that ρ = r . The term V ef f in (10) includes all the potential terms of (6) so that the Lagrangian (6) is that of geodesic equations in a two-dimensional manifold. The nonautonomous equation (9) can always be written in the form of an autonomous third-order differential equation by introducing the new variables r → Y (x) and y → Y ,x . The third-order equation is for which reduction with the autonomous symmetry, ∂ x , leads to the original equation, (9). We found that (13), except the autonomous one, does not admit any other point symmetry vector for any value of the cosmological constant. That is an interesting result because it indicates that there exists a unique relation among equations (9) and (13). On the other hand, it has been found in [25] that equation (9) admits a rescaling symmetry when the cosmological constant is zero. The application of this symmetry vector reduced equation (9) to a first-order Abel equation. However, for nonvanishing cosmological constant only special solutions of the form y (r) = σ 0 r 3 have been derived, where σ 0 = −1, K 4 . Those special solutions correspond to specific initial conditions which describe the asymptotic behaviour of the general evolution of the system.
Below we assume the case of nonvanishing cosmological constant and we perform the singularity analysis. Note that, if (13) is integrable, then (9) is also integrable which means that the gravitational field equations for the Einsteinnonlinear SU (2) σ−model are also integrable, that is, the dynamical system which follows from the action integral, (1), is integrable.
We define the new variable Φ (x) = Y (x) 2 and we search for power-law solutions Φ (x) = αχ p in (13) from where we have the following possible sets (χ = x − x 0 ), where x 0 is the location of the putative movable singularity) and We consider the values of (14) from which we can see that these are the power-law solutions of (9) which we described above. The next step, in order to test if (13) passes the singularity test, is to determine the resonances. Let α = − 2 K , which is the case in which the solution leads to a Lorentzian signature spacetime, for details see [25]. Then by substituting Φ (χ) = αχ −1 + γχ −1+s (16) into (13) and taking the terms linear in γ, we have s 2s 2 − s − 3 = 0, which gives the triple solution Here we remark that the resonances are the same and for α = 1 2 , from where we can say that (13) passes the singularity test and it is integrable.
As far as the solution is concerned, we can write it in a series form in which from s 2 , we have that the powers of χ in the series increase by 1 2 . Therefore the solution is where m 3 and m 0 are arbitrary constants 2 and m 1 , m 2 , m I have to be determined. In particular they are functions of m 0 , m 3 , K and Λ.
Hence we substitute (18) into (13) and find for the coefficient constants that m 1 = 0, and from which we can see that for m 0 = − 2 K it follows that m 2 = 0and m 4 = 0. For values of χ such that χ −1 >> χ 2 the solution takes the following form Now the spacetime, (5), has the following form Furthermore from the second dominant behavior p = − 1 4 we find the resonances which provides us with the second solution which is a right and left Laurent expansion. The free parameters of solution (24) are the n −1 , and n 4 . One important issue of the general solution of (22) is the number of constants which are two from the solution (18); the {m 0 , m 3 } , {n −1 , n 4 } and the constant K which corresponds to the matter source. One would expect three integration constants for the system (3). The latter is hidden in χ, as for the singularity analysis we move to the complex plane in which χ = x − x 0 and x 0 is the position of the singularity. However, that is not an essential constant because it can always be absorbed with the transformation of the coordinate, x, which means that at the end there are only two constants.
On the other hand, as we discussed above, the asymptotic behaviour of the general solution (18) is Φ (χ) ≃ χ −1 , which is the dominant term around the singularity of (13). If we start far from that solution, that is, from different initial conditions, then in the asymptotic behaviour of the solution only the terms of lower powers of χ contribute in the solution. Hence the solution can be described well from the first terms of the Laurent series. We demonstrate that in fig.  1 where the numerical solution of (13) (using Mathematica's NDSolve routine) is given and compared with analytical solutions from the Laurent expansion (18) where we considered the first four terms of (18), (m 0 − m 3 ) (Analytic Sol.1), the six first terms (m 0 − m 5 ) (Analytic Sol..2) and the seven first terms (m 0 − m 6 ) (Analytic Sol..3). Figure  2 presents the absolute error between the analytical solutions and the numerical solution. We observe that those specific initial conditions, where the singularity x 0 is at x 0 ≃ 2.10, are very close in the region of the singularity the approximation works well. Of course eventually the error becomes big and other terms of the Laurent expansion have to be considered.
We remark that in the case of vanishing cosmological constant the singularity analysis provides that the resonances depend on the parameter K. Last but not least an important observation which we can extract from the singularity analysis about the stability of the leading order term is that the solution is unstable. The reason for that is the existence of the terms given by the right Laurent expansion which dominate as far as we are moving from the singularity, see also [28] and references therein. of Conicyt. PGL Leach thanks the Instituto de Ciencias Físicas y Matemáticas of the UACh for the hospitality provided while this work carried out and acknowledges the National Research Foundation of South Africa and the University of KwaZulu-Natal for financial support.