Non-Integrability of Strings in $AdS_{6}\times S^{2}\times\Sigma$ Background and its 5D Holographic Duals

In this manuscript we study Liouvillian non-integrability of strings in $AdS_{6}\times S^{2}\times\Sigma$ background. We consider soliton strings and look for simple solutions in order to reduce the equations to only one linear second order differential equation called Normal Variation Equation (NVE). We study truncations in $\eta$ and $\sigma$ variables showing their applicability or not to catch (non) integrability of models. With this technique we are able to study many recent cases considered in the literature: the abelian and non-abelian T-duals, the $(p,q)$-5-brane system, the T$_{N}$, $+_{MN}$ theories and the $\tilde{T}_{N,P}$ and $+_{P,N}$ quivers. We show that all of them are not integrable. Finally, we consider the general case at the boundary $\sigma=\sigma_0$ for large $\sigma_0$ and show that we can get general conclusions about integrability. For example, beyond the above quivers, we show generically that long quivers are not integrable. In order to establish the results, we numerically study the string dynamical system seeking by chaotic behaviour. Such a characteristic gives one more piece of evidence for non-integrability of the background studied.

In a recent paper, in the context of the AdS/CF T correspondence [1], a precise relation between type II string theory and d = 5 conformal field theories was described [2].Given the AdS space-time string background, when the BPS equations ensuring half-SUSY and Bianchi identity are imposed, a potential function obeying a linear PDE (a Laplace equation) is found with proper boundary conditions.This potential contains all information about the background and, because of the quantization of Page charges, one of the boundary conditions implies that it is written in terms of a Rank function, a convex polygonal function, linear by pieces, with integer values at integer points.Then, this Rank Function is put in contact with a quiver gauge theory: the ranks of the gauge and flavour groups are organized in it.Because the Laplace equation is central in this description, it is referred as an electrostatic description.In the research for half-BPS backgrounds of the type AdS D ×S 2 ×Σ 8−D the same steps can be built for several dimensions D [3] (see references therein).After the field theory string background pair is identified, various tests and predictions of the correspondence follows in a simpler way.In this work we want to take profit of this formalism in order to study Liouvillian (non) integrability of strings in the AdS 6 × S 2 × Σ background.
As the correspondence signals, several characteristics should match in the string side and in gauge theory side.An important characteristic is just the integrability of the model: it helps us in using powerful techniques to study conjectured relations nonperturbatively.In the case of the AdS 5 × S 5 it is already known the existence of integrability structures behind it.The study of integrability, from a classical perspective, can be made by following some specific strategies.One of them is by searching for existence/non-existence of chaos in the associated dynamical system given by the string model: by finding chaotic behaviour (studying, for example, the Lyapunov exponents and Poincare section's structure), integrability is excluded.There is an good number of references of chaos research in string related topics [22][23][24][25][26][27][28][29][30][31][32][33][34][35][36][38][39][40].Another path is to look for a Lax pair formulation of the string model.Lax pairs, if they exist for a classical system, can be used to generate a tower of integrals of motion and this will give support to integrability.In the case of string theory, this approach is very successful when the string background is coset-like, as in the case of AdS 5 × S 5 [41][42][43] and the simple case of R × S 3 [44].But in fact, there is no general guide to find a Lax pair formulation.
The last path we would like to cite is an analytical method to discuss the existence of integrability in a dynamical system.The idea is to find a string soliton and show that the dynamics of such an object is (non) integrable in some sense.This method has been recently used to study integrability of a lot of string backgrounds, their respective duals and in other models [35,36,38,[45][46][47][48].Given a system of differential equations, the analysis of the variational equation around a particular solution can show its (non) integrability.In other words, if a nonlinear system admits first integrals, the variational equation will admit it too.Disproving this for a given class of functions will imply in the non integrability of the initial nonlinear system.The mathematical establishment of integrability through the normal variation equation (NVE) has been made by some tests that were improved along the years.First, there is Ziglin's theorem relating the existence of first integral of motion with monodromy matrices around the straight line solution, the basis to linearize the system of differential equations [49,50].After that, techniques of differential Galois' theory applied to the NVE equation were introduced [51][52][53].In this work we make use of the improvement made by Kovacic [54].It gives, through an specific algorithm, an answer to the existence of integrability: once the NVE is written in a linear form with polynomial coefficients, it suffices to check a group of criteria.In fact, Kovacic provided a way to construct the solutions.In the case of string models, it basically consists in the following: first we find the equations of motion for the l degrees of freedom of a proposed string soliton.Next, we find simple solutions for (l − 1) of these equations which are replaced in the last one.They give us the normal variation equation (NVE).It is a linear second order differential equation given by z + Bz + Az = 0.
With the equation above at hand, we can use the Kovacic's criteria to seek if a Liouvillian solution do exist.As will be explained, the functions A, B and its derivatives determine the existence of a closed form of Liouvillian solutions.We should point out that the non-integrability in some of these backgrounds has been studied in Ref. [55].The authors use the assumption that they allow for a regular critical point for the warp factor, which is necessary to have a GKP string solution.With this assumption they were able to show non-integrability for the background of Ref. [56].However that background does not include D7 branes and further generalizations [8,9,12,19].Namely, the only explicit examples without D7 branes considered here are the T N and + M N .For these, our results enforce the findings of Ref. [55] about non-integrability.In this paper our results include D7 branes and therefore are more general: we are able to find non integrability for a larger class of backgrounds.
In this work we study analytical and numerical (non) integrability of strings in AdS 6 ×S 2 ×Σ background in direct connection to its 5D Holographic Duals, as suggested by the electrostatic method which is reviewed in section 2 .In section 3 we study the string dynamics in the given background and write the NVE that will be basis for the conclusions.In sections 4 and 5 we apply Kovacic's criteria to several potentials, running for regions where σ = σ 0 and η = η 0 , including those supporting quiver gauge models.In section 6 we discuss the potential for a general case, studying how to see and avoid logarithmic terms at the final U -function in order to apply Kovacic's criteria.In section 7 we supplement the results coming from the analytical process by finding chaotic behaviour in the dynamical system.We compute string trajectories, power spectra, Lyapunov exponents and Poincaré sections in order to establish the results.Finally, we present our conclusions.

The Electrostatic Description
In this section we quickly review the correspondence between strings in AdS 6 × S 2 × Σ background and d = 5 SCFT through the electrostatic viewpoint.

The Background.
The type IIB background that is important in this work is described in Ref. [2].The full configuration consists of the metric, the dilaton, B 2 , C 2 and C 0 -fields in the NS and Ramond sectors respectively.In string frame it is given by In the equations above the range of η and σ are the intervals [0, P ] and −∞ < σ < ∞ respectively.We also use the following parametrization The background depends only of one potential function V (σ, η), which solves a linear partial differential equation given by Next, we define to arrive at a Laplace equation given by The boundary conditions are Due to the Laplace equation above and boundary conditions, the authors of Ref. [2] called this an "electrostatic description" of the system.The solution of the equation above is given by where the rank function is obtained as R(η) is interpreted as a charge density at σ = 0, extended along 0 ≤ η ≤ P (P ∈ Z).Because of the quantization of Page charges, it is fixed to be a convex piece-wise linear function.It is this function that will make possible the precise connection between the string model and the quiver gauge models in the correspondence.

The Holographic Duals.
Now we review the holographic duals of the background just described.For this we follow the lines of Refs.[2] and [13].In order to find a general solution, the quantization of the Page charges is necessary.This and the boundary conditions (5) enforce the Rank function R(η) to be given by Depending on the choice of the Rank function R, the number of D7, D5 and N S5 branes can be determined.This also fix the strong coupling CFT to which a linear quiver theory flows.The gauge group of the linear quiver is given by Π P −1 i=1 SU (N i ), with each gauge connected by bifundamental hypermultiplets.Finally, some of the color groups can have SU (N f ) flavor groups.However, the relation to the holographic dual is trustable only in the limit of very large P .The only exception to the general solution above are the abelian and nonabelian T-duals, in which the boundary conditions (5) are not satisfied and the dual CFTs are not well known.In the rest of this section we give the details of some particular solutions we will consider later.

Abelian and non-Abelian T-Duals.
The first AdS 6 solutions to Type IIB supergravity were first constructed in [20,21] by acting with Abelian and non-Abelian T-duality on the Brandhuber-Oz solution to massive IIA.Later, it was shown to fit in the construction of D'Hoker, Gutperle and Uhlemann [12] and therefore also in the electrostatic description of [2].Despite of being the first IIB constructions, the dual theories are not yet well known (see [12] and references therein).
The abelian case includes D7/O7-branes and provides an example of a Riemann surface with the topology of an annulus [12].The electrostatic description has a potential given by The non-Abelian case arises from the upper half-plane [21].In the electrostatic description it has potential given by These potentials can be obtained as partial polynomial expansions for σ ≥ 0 from the full electrostatic description [2].We should point that they do not satisfy the boundary conditions (5).

2.2.2
The T N and + M N Theories.
The electrostatic description for these models were given in Ref. [2].The T N theory was first studied in Refs.[15][16][17].It is the strongly-coupled UV fixed point of the linear quiver gauge with junctions of N D5, N N S5 and N (1, 1) 5-branes.Its electrostatic description has potential given by The + M N theory was already studied in Ref. [14] and is defined on the intersection of N D5 and M N S5-branes.In field theory it can be defined as the UV fixed point of the linear quiver gauge theory.In this case, the electrostatic description has potential given by In Fig. 1 we give a brane diagram2 in order to visualize the quivers above.We should point that both the theories above have N f = 2N .For more examples with this number of flavors see Ref. [13].

2.2.3
The TN,P and + P,N Theories.
These theories were both studied in Ref. [2].The rank function associated to TN,P theory is given by The number of D7-branes is given by P N .The number of D5-branes is determined by the value of R(η) at the points η = 1, 2, 3, etc.With this we get a total of N P (P − 1)/2 D5-branes.Finally, we also have P NS5-branes.In fig. 2 we give this quiver 3 .N 2N 3N . . .PN (P-1)N Figure 2: The TN,P quiver, with P N D7-branes, N P (P − 1)/2 D5-branes and P NS5-branes.
The potential can be computed by replacing (12) in Eq. ( 6) and is given by Since the theory is trustable only for large P , we must expand the potential.The expression was found in Ref. [2] and is given by For + P,N theory the associated rank function is defined as In this case we have N D7-branes at η = 1, N D7-branes at η = P − 1.We also have N D5-branes at each point η = 1, 2, ..., P − 1.With this we get a total of N (P − 1) D5-branes.In fig. 3  The potential can be found by replacing (15) in Eq. ( 6) and is given by We should point that, as before, the above theories are trustable only in the large P limit and we will see that this will give us the very simple potential In the next section we study the dynamics of strings in these backgrounds in order to seek for integrability in an analytic way.
3 Dynamics of Strings in AdS 6 × S 2 × Σ.Now we consider the dynamics of strings in the background (1) that is the subject of this work.The string action is given by supplemented by the Virasoro constraints Our soliton is a string at the center of the AdS space, which rotates and wraps on the following coordinates (τ and σ are the world-sheet coordinates) Here κ is an integer number that indicates how many times the string wraps the corresponding direction.It is straightforward to check that this truncation is consistent, since they are a solution of the generic equations of motion.For this, the coordinates above must be subject to the equations of motion as described below.We get the effective Lagrangian and The equations of motion can be obtained from the Lagrangian above and are given by In the first of the equations above, E is a constant of integration and has been used in the last three equations.It is easy to verify that the derivative of the Virasoro constraints (22) vanishes if Eqs. (23)(24)(25)(26) are used.Since Eqs.(23)(24)(25)(26) define the τ evolution of the string configuration, we can study its (non) integrability.In the next, we consider the possibility of finding simple solutions and study these aspects for the configuration (20).

Finding Simple Solutions.
The first step we take is to look for some simple solutions of the eom (23)(24)(25)(26).
As cited in the introduction, the general procedure is to find a solution to Eqs. ( 25) and ( 26) which must be replaced in the NVE of Eq. ( 24).However, we just need to solve (25) or (26) and use the constraint (22).We will see that, with this at hand, we obtain general conclusions without choosing any specific form of the background.In the next sections we will apply this idea to many cases.First we note that χ = χ = χ = 0 is a solution to the second equation in (24).Replacing this in the other equations we get The equations above can be further simplified.By using equation ( 23), the constraint can be written as With this we get Finally, we fluctuate χ by χ = 0 + z(τ ) in equation ( 24) to get the NVE The coefficients A and B depend on σ, η.Therefore, in principle, we should solve for σ, η in such a way that the NVE becomes a linear second order differential equation.However, if we choose the simple solution σ = σ 0 = constant or η = η 0 = constant, we see that a linear equation can be obtained.We analyse now both cases.
We note that we can have a simple solution of Eq. ( 28) given by With this the NVE (30) and η equations are simplified to and η = 1 2 Despite the simplification, the coefficients of the NVE are yet η dependent.The general procedure is to choose a specific background and solve for η.However, depending on the background, solve equation ( 33) and determine η can be a very difficult task.Besides this, we would need to choose a specific background, which would spoil the generality of this study.
In order to solve that, we remember that the constraint ( 27) reduces to Now, by using the result above in Eq. (33) we see that the equation of motion for η can be written as From equations (34) and (35) we see that η and η depends only on η(τ ).This suggests that we use τ = τ (η) and we arrive at a new NVE given by Now we use ( 34) and ( 35) to get where the quantities above must be taken at σ = σ 0 and z means derivative with respect to η.Therefore, there is no need to solve the equation for η in order to obtain the final NVE.
In the next step, we consider that a simple solution of Eq. ( 29) can be found and is given by In this case the constraint ( 27) simplifies and we get From the equations above we could determine σ.However, as in the σ = σ 0 case this will not be necessary.By using ( 37) and ( 38), the σ equation ( 28) becomes From equations ( 38) and (39) we see that σ and σ depend only on σ(τ ).This suggests that we use the parameter τ = τ (σ) and the NVE, equation (30), becomes Now, by using (38) and (39) we finally write where the prime means a derivative with σ and the quantities above must by taken at η = η 0 .Again, we point that there is no need to solve the σ equation in order to obtain the second order linear differential equation.The results of the last subsections show that the behaviour of is crucial in order to discover what is the consistent truncation of our system of equations.In the next sections we will apply the results above and analyse some specific backgrounds.Later we generalize our results.

The Kovacic's Criteria of Liouvillian Integrability.
As shown in the last subsections, we can find consistent truncations of string equations in order to get the NVE as a homogeneous second order linear equation.With this at hand, we can study Liouvillian integrability.Interestingly, Kovacic provided not only an algorithm to find the solutions, but also a set of necessary, but not sufficient, conditions to analyse if an equation is Liouvillian integrable.Consider the general NVE in Eq. (30).First, we transform it to a first order differential equation in the following way where With the potential U (x), Kovacic set some general conditions for integrability.First of all, his criteria is valid only if U (x) is a fractional polynomial.If this is the case, the conditions are the following: • Case 1: every pole of U (x) has order 1 or has even order.The order of the function U (x) at infinity is either even or greater than 2.
• Case 2: U (x) has either one pole of order 2, or poles of odd-order greater than 2 .
• Case 3: the order of the poles of U (x) does not exceed 2, and its order at infinity is at least 2.
If none of the conditions above are satisfied, the analytic solution (if it exists), is non-Liouvillian.With the help of the criteria above, we will study integrability of the string dynamic equations in the rest of our manuscript.The interested reader can find more detailed explanations and specific examples of this method in Refs.[36,38,45]. 4 The Case σ = σ 0 .
In this section study the possibility of obtaining a consistent truncation given by σ = σ 0 .We apply the results of the last section to some simple specific backgrounds and, later, we consider the general case.We first consider the Tduals, the (p, q)-5-brane system and finally we analyse the possibility σ 0 = 0 in a general way.

Type IIA Abelian T-dual.
We now study the abelian T-dual of the D4/D8 system in massive Type IIA theory [12].In the electrostatic description the potential is given by Eq. ( 8) and therefore we have As said before, in order that σ = σ 0 be a solution of Eq. ( 28) we must obey for some σ 0 .From the last expression we see that the boundaries σ 0 = 0, σ 0 = ∞ are not solutions to the equation above.However, we can try b 1 + 4b 4 σ 3 = 0. Using this in Eq.( 36) we get where we have used a ≡ b 1 /(4b 4 ).Since a > 0, we get complex coefficients and σ = σ 0 is not a consistent truncation of the equations of motion.

Type IIA non-Abelian T-dual.
The next simple background is given by the non-Abelian T-dual of the type IIA D4/D8 system.The electrostatic description is given by ( 9), namely and 1 For σ = σ 0 to be a solution we need that Just as with the abelian case, σ = 0 and σ = ∞ are not solutions.However, we can try 16a 4 σ 3 + a 1 = 0.In this case we get, from Eq. ( 36), that where we have used a ≡ a 1 /(16a 4 ).Since a > 0 we get complex coefficients σ = σ 0 is not a consistent truncation of the eom.

4.3
The TN,P and + P,N Theories.
As said before, the theory above is trustable only for large P .The TN,P solution is studied in Ref. [2] with potential given by Eq. ( 14).Now we expand it around σ = 0 to get In Ref. [57] the author gives this potential up to order σ.However we introduced terms of order σ 2 since the background functions depends on second derivative of σ.With the potential above we get and therefore σ = 0 is a good truncation.With this we get Then the potential becomes .
This last potential is not polynomial and we cannot use the Kovacic's criteria.However we can change the η variable to y = √ η, compute a new NVE and write the following Ũ (y): Now we have it in a polynomial form and, even after that procedure, this last potential does not obey any of the Kovacic's criteria: we conclude that this model is not integrable.Next we consider the + P,N theory.This solution, studied in Ref. [2], has potential given by Eq. (17).By expanding it in σ we get From the result above, we have Hence the truncation σ = 0 is consistent.We can also obtain the U function.
Nevertheless, we see that we get a logarithm dependence in η.This potential, thus, cannot be written as a polynomial function and we cannot apply the Kovacic's criteria.We point that the cases above were studied in Ref. [57].The author used a further truncation η = 0.However, if we use η = σ = χ = 0 in (22), we find ṫf 1 = 0. Then we get ṫ = 0 or f 1 = 0.The first choice is not consistent and the second give us a null background making the equations of motion trivially null.Therefore, with this choice it is not possible to conclude about integrability.From the viewpoint given in our work there is no need to fix η = 0 to make the truncation consistent.Despite this, due to the shape of the potential, the Kovacic's criteria were applied only for the TN,P .
In the previous section we learnt that σ = σ 0 sometimes is a good solution for the string's equations of motion, sometimes not.In this section we apply the results of section 3 to search the other possibility of obtaining a simple solution for η = η 0 .We show that this works for several backgrounds: the abelian and non-abelian T-duals, the T N , + M N , TN,P and + P,N theories.Finally, we consider the general case, which includes any long quiver.

Type IIA Abelian T-dual
As we already know, in this case the potential is given by Eq. ( 8) and therefore we have 1 We should point out that the equation above is valid for any value of η.Then, the coefficients of the NVE (40) becomes where we have used κ = E.The coefficients are already rational functions and can be analysed with the Kovacic's criteria.We get a potential (42) given by In the expression above we have used b 1 = 2ab 4 in order to simplify it.By analysing the U −function, we see that it does not satisfy all the three possible necessary Kovacic's conditions described above.The solution to the equation should then be non-Liouvillian.Then η = η 0 give us a consistent truncation in order to study integrability.

Type IIA non-Abelian T-dual
For this case we have the potential given by Eq. ( 9) and As in the abelian case, again we have the equation above valid for any value of η and a simple solution to ( 29) is given by η = η 0 .Therefore, the coefficients of the NVE (( 40)) become , and where we have used κ = E and a = a 1 /(8a 4 ).We see that the coefficients above are not rational functions.We could try to solve this with a change of coordinate but, since the term in the square root is cubic, it is not possible to find an inverse.But we can note that all terms with square roots cancel at η 0 = 0. Therefore we choose this particular value to get The coefficients above are already rational functions and we find from Eq.( 42) that 2 Analysing the U −function, we see that it does not satisfy all the three possible necessary Kovacic's conditions and the solution must be non-Liouvillian.However, this case is less general if compared to the abelian one, since we had to choose η 0 = 0. We can also conclude that η = 0 provides a consistent solution that allow us to analyse the integrability of the system.

T N Theory
This solution is studied in Ref. [13].Its electrostatic description has potential given by (10) [36] Since we look for a solution with η = 0, we can expand the potential above to obtain with , and With the quantities above we get where With this, it is clear that in the limit η = 0, Eq. ( 52) is null.Therefore, it is a consistent truncation in order to study the NVE (40).The same kind of procedure will be repeated in all the next cases.The coefficients are given by + 9N e .
The potential (42) will be given by We cannot apply the Kovacic's criteria to the case above since the coefficients are not fractional polynomials.However, for very large σ it reduces to where O(e −cσ ) means terms which decrease exponentially.Therefore, in this region, the potential has the desired shape.Nevertheless, the order of the fractional polynomial is negative and does not satisfy all Kovacic's conditions.Since we are analysing only one particular point, this is only a hint and we cannot conclude that the solution is not integrable.We need more regions to enforce this.

+ M N Theory
This solution is studied in Ref. [13], with electrostatic description given in Ref. [36].The potential is given by ( 11) Now we expand in η as in the last subsection to get , and With the expressions above we find that with Again we see that the last expression is null for η = 0. Therefore, this is a good truncation of our system and we can study the NVE (40) .
Finally, the potential ( 42) is given by .
Again, we have that the potential above is not a fractional polynomial.However, it is simple to show that for very large σ we have where O(e −cσ ) means terms which decrease exponentially.Therefore, in this region, we can test the Kovacic's criteria by analysing the U −function.We see that it does not satisfy all the conditions.As before, this is a hint and we cannot conclude that the solution is not integrable.

5.5
The TN,P and + P,N Theories.
In this subsection we analyse the integrability of the TN,P and + P,N theories.These theories are trustable only in the large P limit (so we must expand them in 1/P ).The TN,P solution is studied in Ref. [2] with potential given by Eq. ( 14).Now we must expand it around η = 0 to get With the above potential we get where we have used a 2 = 24P 2 log(2).Therefore η = 0 is a good truncation.With this we find .
where b 2 = 8P 2 log(2).Therefore the potential becomes Just like the σ = 0 truncation, we find a polynomial function and we conclude, by using the Kovacic's criteria, that it is not integrable.
Next we consider the + P,N theory.This solution, studied in Ref. [2], has potential given by Eq. 17.By expanding the potential around η = 0 we get Therefore, η = 0 is a good truncation.We compute .
For the U function we find However, just as the σ = 0 case, we have a logarithm and therefore we can not obtain a polynomial expansion for the U function.The conclusion is the same: we can not apply the Kovacic's criteria for the + P,N theory.

General Observations about Expansions and
Integrability.
In this section we study some aspects of Polylogarithmic expansions in order to understand how to apply (or not) the Kovacic's criteria. where This express the fact that the Polylogarithm Li s (z) is singular close to z = 1.The emergence of this logarithm is problematic since they are not polynomial and we can not apply Kovacic's criteria.The only particular case in which this does not happens is when s = P − 1.This is due to the fact that e − π(−is+iη) P = −e − π(i+iη) P .In this case, the Polylogarithm is polynomial around z = −1.The only case in which we have an s = P − 1 contribution is the TN,P theory, where we have c s = N P δ s,P −1 .Any other case will contain s = 1 and we necessarily have logarithm contributions an in the + P,N case considered before.This explains the fact that, for the quivers, we have got a polynomial U function only for the TN,P case.We also point that we did not expand in 1/P .Therefore the above result is valid for any quiver, such as the + M N .

The Truncation at σ = 0.
We consider now general aspects of σ = 0 truncation.We stress that the Type IIA Abelian and non-Abelian T-duals are not in the class discussed in the last subsection.However, as we saw before, the truncation η = 0 is consistent for these cases and we were able to study integrability.In this case we must expand the potential (55) in σ to get We find that and therefore ∂ η f 3 is linear in σ.We conclude that, for any theory described by the potential (55), σ = 0 is a good truncation.We also point that we did not expand in 1/P .Then the above result is valid for any quiver, such as the + M N .When analysing integrability, for the cases above, we arrive at the same problem as in the η = 0 case.This is the reason why, also in σ = 0, we were not able to apply the Kovacic's criteria.Another important point is that the Type IIA Abelian and non-Abelian T-duals are not in the above class.As we saw before, the truncation σ = σ 0 is not consistent for these cases and we must use the truncation η = 0.

The Truncation at Large σ 0
In the last subsections we saw that, despite of being consistent, the truncations η = 0 and σ = 0 are not good to study integrability of long quivers.The reason is the appearance of logarithm dependences in the U function and this spoils the applicability of the Kovacic's criteria.However, if we consider σ = σ 0 with σ 0 large, we can circumvent that difficulty.
In order to analyse the long quiver with some generality we remember that the general potential is given by [2] V (σ, η) = where c s does depends on P .In the limit of large σ/P limit, only the k = 1 term survive and we have Now we consider the long quiver, with large P >> 1 where we have defined From the expression above we find which is null in the large σ limit.Therefore, we can apply the method if we consider σ 0 , P very large.As explained before, we use the NVE (36) with Therefore the U function is given by By analysing this last U , we see that it does not satisfy all the three necessary Kovacic's conditions.The solution to the equation should then be non-Liouvillian.The limit above includes the (p, q)-5-Brane background studied in Ref. [4] and also any long quiver.It is important to note that f (P ), given in Eq. ( 58), determines what is the specific long quiver.However, f (P ) does not appear in all the expressions and particularly in the U function.Therefore, we conclude that strings propagating in any long quiver, and (p, q)-5-Branes are not integrable.

Numerical Analysis.
In this section we change the focus and turn to a numerical analysis of the string dynamical system in the AdS 6 background.The main goal here is to find signals of chaotic motion due to the non-linearity of the eom for the quivers given.In finding chaotic behaviour, we give another piece of evidence for non-integrability of the models studied in this paper.Basically, chaos can be understood by the high sensitivity of the dynamical system to the initial conditions: given two trajectories with different (but close) initial conditions, they will drastically differ from each other after time evolves.The way the system behaves shows some characteristics that can be measured.For example, the distance between trajectories with time can be computed by the Lyapunov exponents, as we will see.Another interesting characteristic is related to the periodicity of the system.In chaotic motions periodicity is lost, and this can be seem by computing the power spectra associated to the trajectories.Despite this, by studying the Poincaré sections we can understand the behaviour of trajectories in the phase space in order to see how periodicity is lost if we, for example, increase the energy of the string soliton.All of these can be better understood in the nice reference [36] (see also [37]).In the discussion below we explain in more details these aspects.We show, as examples, results related to the backgrounds associated to the + P N and TN,P quivers described respectively by the potentials ( 17) and (14).
First of all, it is important to remember that the string soliton is placed at the center of the AdS 6 part of the background.In this way, its motion is restricted to the S 2 × Σ part.In other words, we should study the behaviour with time of the coordinates χ(t), σ(t) and η(t).Furthermore, due to the need of large P for a trustable quiver field theory, we let the η coordinate runs to a large interval but still being bounded.In another words, the range of η and σ are the intervals [0, P ] and −∞ < σ < ∞ respectively.This is closely related to the discussion in the work done in [36] where the coordinates evolved are bounded.Given these, we turn now attention to the study of trajectories and some of its characteristics.With this in mind we made a numerical evolution of the hamiltonian equations of motion coming from eq. ( 21).

String Trajectories and Power Spectra.
Related to the trajectories, the first characteristic we call attention is to their almost periodic oscillations for low energies.In figure (4a) we show a plot of trajectories with time (the horizontal line) for the TN,P quiver with E 11.We represent in red, green and blue respectively the σ, η and χ coordinates.Despite some small oscillations, the general behaviour is almost periodic (the real non-periodicity will be shown in the subsequent discussion).If we increase the energy, then we see the trajectories in figure (4b) with more oscillations: despite the fact that the time axes don't show the same regions, we clearly see an increasing in frequencies, for example, if we compare the plots between t = 0 and t = 200.Another interesting point is the change of the form of oscillations with different energies.In a periodic system we would wait change of frequencies and amplitudes, but no change of, for example, sinusoidal description of them.In order to make the numerics we chose P = 100 in this work.The initial conditions taken for the energies are, for E = 11.04,χ(0) = 0.1, p χ (0) = 0.006, σ(0) = 0.1, p σ (0) = 0.0001, η(0) = 0.03, p η (0) = 0.006.For E = 330.41,χ(0) = 1, p χ (0) = 2, σ(0) = 7, p σ (0) = 2, η(0) = 4, p η (0) = 5.For the +P N quiver we ran the initial conditions χ(0) = 3, p χ (0) = 5, σ(0) = 3, p σ (0) = 5, η(0) = 2, p η (0) = 5 to get E = 267.18.In this last case, we again can see basically the same trajectory behaviour for a high energetic string already present in the TN,P case.A better way to see a chaotic signal is to plot the behaviour of the two compact coordinates χ and η (χ is a bounded coordinate in S 2 and η is bounded too, even for higher values, due to the quiver potential definition).We can see in figures 6 and 7 that characteristics for different energies for the TN,P quiver.We clearly see the string bouncing back   again and again at η = 100.For a high value of the energy, the trajectories in fact become more and more complex.
In order to catch non-periodicity characteristics for the string trajectories presented above, we can compute their power spectra.By taking Fourier transforms of the coordinates, the periodic, quasi-periodic or chaotic behaviour will appear.When a signal is perfectly periodic with a frequency ω its Fourier spectrum will show a vertical line at the characteristic frequency of the system.If the signal is chaotic, what we expect for high energies, the power spectra presents a noisy band of frequencies.In the figures 8a and 8b we see exactly that behaviour for the quivers + P N and TN,P .For the + P N case we set the string energy to E 534 in which the initial conditions are given by χ(0) = 3, p χ (0) = 10, σ(0) = 3, p σ (0) = 10, η(0) = 2 and p η (0) = 10.In this way we see clearly its noisy power spectra due to chaotic behaviour.

Lyapunov Exponents.
With the trajectories at hand, we can now search for more concrete chaos signals.We can estimate the Lyapunov exponent related to the dynamical system of the string soliton.For dynamical systems with a non-zero Lyapunov, the time evolution associated to two close trajectories in the phase space sensitively de-      we obtain several values for χ(0) such that the Virasoro constraint is obeyed for a given value of energy.For the quiver TN,P , for E 330, we chose σ(0) = 7 and computed random values for p σ (0) ∈ [0, 5], p χ (0) ∈ [0, 10], η(0) ∈ [0, 20] and p η (0) ∈ [0, 10].We plot the points (η, p η ) every time χ(t) = 0 in figures 11 and 12, including sections for lower energies.
In these pictures we can see how the organization of points changes if we change the string energy.For lower energies, E 11 for TN,P and E 5 for + P N we can see a more or less organized pattern of quasi periodic trajectories.The absence of KAM circular curves shows up for the higher energies pointed above for both quivers.In another words, the distribution of points does not show any organization and, in fact is lost.And these just add to the conclusion of existence of non-periodic trajectories or, in other words, to the existence of chaos.

Conclusions.
In this manuscript we studied non-integrability of strings in AdS 6 × S 2 × Σ.By using the electrostatic method [2], we considered a very large class of theories: the long quiver, the abelian and non-abelian T-duals, the (p, q)-five-brane system, the T N and + M N theories.In the electrostatic method it is claimed that general dual backgrounds can be written in terms of a only function given by Eq. ( 6) This method gives a clear way to link the string theory to its dual gauge theory.Of course, if the AdS/CFT conjecture is proved correct, results related to both sides of the duality will follow as consequence.With this consideration, the idea of this paper, following previous studies, is to point that the   holographic duals should be (non) integrable if the string equations are (non) integrable.For this, we first carefully studied the general dynamics of strings in the background (61) by regarding a consistent truncation of the string equations in the supergravity background.The usual procedure of the literature is to find simple solutions for (l − 1) of the equations, which can be replaced in the last one to obtain the NVE.However we shown that one of these equations can be replaced by the Virasoro constraint and this provides a simplification of the procedure.In this case we have three equations for σ, η and χ.Previous results of the literature argue that we must solve the equations for η and σ and replace this in the variation of χ in order to find an homogeneous second order linear equations [36,38,45].However, depending on the background, this can be a difficult task.By using the Virasoro constraint we shown that we just have to find a simple solution for η or σ.The fact is that the Virasoro constraint fixes the value of the Hamiltonian of the 1d system (21) to zero, and the constraint H = 0 will always imply one of the equations of motion.We also found a general condition to discover how simple this solution can be.We found that η = η 0 or σ = σ 0 , respectively, are consistent truncations if We first applied the conditions above to look for the simple solution σ = σ 0 .We applied it to the abelian and non-abelian T-duals, the TN,P and + P,N cases.We found that for the abelian and non-abelian T-duals cases σ = σ 0 is not a consistent truncation.For TN,P we found σ = 0 a good consistent truncation, providing the applicability of Kovacic's criteria: in this case the equations of motion are non-integrable.For the + P,N we found again σ = 0 as a good truncation, but we found logarithmic dependences on the coefficients which turn into non-applicability of the Kovacic's criteria.In summary, σ = 0 is not always a good truncation, as suggested in previous studies [36,38,45].Turning to the η = η 0 possibility of truncation we found that all cases can be studied.More specifically, η = η 0 is a good truncation for the abelian T-dual case, η = 0 is good for the non-abelian T dual, T N , + M N , TN,P and + P,N cases.For the non-abelian T-dual case we directly applied the Kovacic's criteria finding it a non-integrable model.For the T N and + M N cases we found coefficients in non-fractional polynomial form, which do not makes possible the application of the criteria.However, we turn to the large σ limit in which we can apply them: in this region, the criteria tell us that these models are non-integrable.For the TN,P case, as in σ = 0 truncation, we found it again non-integrable, as expected.Nevertheless, in the + P,N we found again logarithmic dependence inside the coefficients which turn into non-applicability of the Kovacic's criteria.
Next, we considered the general case, with arbitrary potential given by (61).We found, for σ = σ 0 and for η = η 0 , a potential which is, in general, not a fractional polynomial and therefore the Kovacic's criteria cannot be applied.Nevertheless, we found the origin of logarithm contributions to the coefficients entering in the NVE equations by studying expansions of Polylogarithmic functions appearing in the general potential.In order to circumvent these contribu-tions we shown that we can go to the large σ and P limits, finding well behaved coefficients which, consequently, gives us a nice U − function to apply the criteria.In this case, by analysing the pole structure, we were able to show that it does not satisfy the Kovacic's conditions and, therefore, it is not integrable.It is interesting that in these limits we found an universal potential for long quivers.With universal we mean that f (P ), given in Eq. ( 58), determines what is the specific long quiver and it does not enter in the final expressions for the U function.This limit describes the long quivers and we conclude that all long quivers are not integrable.This includes not just the ones studied here, but all its versions, as described in in Ref. [13].This includes the + N,M , T N , Y N , + N ,T 2K,K,2 , T N,K,j and + N,M,j theories [11,[14][15][16][17]19].
Finally, we made a numerical treatment of the evolution of the string dynamical system.The main goal was to give more evidence of non-integrability by searching for chaotic behaviour hidden in the string eom.We plotted the trajectories for quivers + P N and TN,P cases, observing the increase of path's complexity for coordinates evolving in time.We pointed the sign of chaotic trajectories by the plots of the bounded coordinates η and χ, showing its behaviour for higher energetic strings.We added to this the power spectra of trajectories by showing in fact the non-existence of periodic motion for the quivers treated.After this we made computations of Lyapunov exponents for the + P N and TN,P cases.We concluded that the large time behaviour of trajectories contributes to non-zero values of Lyapunov exponents which signals for chaos underlined in the system.For last, we obtained its Poincaré sections.These sections points to the non-periodicity of trajectories in phase space as we increase the string energy: the more it grows, the more non-periodic is the behaviour of trajectories.This is another signal of chaos.Therefore, we add these to the non-integrability box of the AdS 6 background.
After all discussion, the main conclusion of this paper is that, for a very large class of models, we show analytically and numerically that the string equations of motion in the AdS 6 × S 2 × Σ background are not integrable.An interesting point to be studied is the near Penrose limit behaviour of these models.In this limit the string model is solvable and quantizable.However, it is interesting to see what happens when the metric is perturbed and look for, as an example, chaotic motion and non-integrability appearing again.This question is left for a future work.

Figure 1 :
Figure 1: Brane diagrams for the a) T N and b) + N,M theories.The (p, q)-5branes are represented by straight lines and 7-branes represented by filled black dots.The angles are determined by the (p, q) charges and the 5d SCFTs are realized by intersections at a point.

Figure 4 :
Figure 4: String trajectories evolving in time in the case of the TN,P quiver.In red, green and blue we show respectively the σ, η and χ coordinates.It is clear how complex they become as we increase the energy from E 11 to E 330.
Plot of χ coordinate for E 267.

Figure 5 :
Figure 5: Plots with an example of trajectories for the + P N quiver.In this case we can see basically the same kind of behaviour present in the TN,P quiver.

Figure 6 :
Figure 6: In this plot we show the behavior of the compact coordinates χ and η (with cos χ rescaled by 10) with E 11 for the TN,P quiver.

Figure 7 :
Figure 7:  In this plot we show the behaviour of the compact coordinates χ and η (with cos χ rescaled by 10) with E 330 for the TN,P quiver in order to see more complex trajectories.

Figure 8 :Figure 10 :
Figure8: Plots of power spectra for high energetic strings in the quivers + P N and TN,P .In blue we show the spectra for the χ coordinate and, in orange, for η.There is no principal frequency appearing in any case.The spectra are dominated by noise.