Conformal Flow on S 3 and Weak Field Integrability in AdS 4

Transcription

1 Commun. Math. Phys. 353, ) Digital Object Identifier DOI) /s Communications in Mathematical Physics Conformal Flow on S 3 and Weak Field Integrability in AdS 4 Piotr Bizoń 1,4, Ben Craps 2, Oleg Evnin 2,3, Dominika Hunik 1, Vincent Luyten 2, Maciej Maliborski 4 1 Institute of Physics, Jagiellonian University, Kraków, Poland. bizon@th.if.uj.edu.pl; dominika.hunik@uj.edu.pl 2 Theoretische Natuurkunde, Vrije Universiteit Brussel and The International Solvay Institutes, Brussels, Belgium. Ben.Craps@vub.ac.be; Vincent.Luyten@vub.ac.be 3 Department of Physics, Faculty of Science, Chulalongkorn University, Bangkok, Thailand. oleg.evnin@gmail.com 4 Max-Planck-Institut für Gravitationsphysik, Albert-Einstein-Institut, Golm, Germany. maciej.maliborski@aei.mpg.de Received: 30 August 2016 / Accepted: 20 March 2017 Published online: 13 May 2017 The Authors) This article is an open access publication Abstract: We consider the conformally invariant cubic wave equation on the Einstein cylinder R S 3 for small rotationally symmetric initial data. This simple equation captures many key challenges of nonlinear wave dynamics in confining geometries, while a conformal transformation relates it to a self-interacting conformally coupled scalar in four-dimensional anti-de Sitter spacetime AdS 4 ) and connects it to various questions of AdS stability. We construct an effective infinite-dimensional time-averaged dynamical system accurately approximating the original equation in the weak field regime. It turns out that this effective system, which we call the conformal flow, exhibits some remarkable features, such as low-dimensional invariant subspaces, a wealth of stationary states for which energy does not flow between the modes), as well as solutions with nontrivial exactly periodic energy flows. Based on these observations and close parallels to the cubic Szegő equation, which was shown by Gérard and Grellier to be Lax-integrable, it is tempting to conjecture that the conformal flow and the corresponding weak field dynamics in AdS 4 are integrable as well. 1. Introduction Propagation of nonlinear waves in confining geometries presents significant challenges because the key mechanism stabilizing the evolution of waves on unbounded domains, which is the dispersion of energy by radiation, is missing in confined settings. Consequently, an arbitrarily small perturbation of a ground state can lead to complicated long-time behavior. The central physical problem in this context is that of energy transfer, namely, understanding how the energy injected into the system gets distributed over the degrees of freedom in the course of its evolution. This problem has been studied in the physics literature from a statistical viewpoint under the name of wave weak) turbulence [1,2], but it has only recently emerged as an active research topic in dispersive PDE theory. In this deterministic mathematical approach, the energy cascades from low to high modes are quantified by the growth of higher Sobolev norms of solutions, and

2 1180 P. Bizoń, B. Craps, O. Evnin, D. Hunik, V. Luyten, M. Maliborski the main question is whether these norms can become unbounded in finite or infinite time. The past few years have witnessed significant progress in the understanding of this issue in the context of nonlinear Schrödinger equations and their variations) on tori [3 5]. However, almost nothing is known, as far as we can tell, about manifestations of these phenomena in other interesting evolution equations of mathematical physics. Our attention in this article will be focused on the weak field dynamics for the conformally invariant cubic wave equation on the three-sphere. This very natural geometrical model describes relativistic scalar waves propagating on a compact manifold and interacting with the background geometry and among themselves. The spherical geometry prevents the waves from dispersing to infinity and thus ensures that the nonlinear self-interactions remain important for all times, inducing complicated energy transfer patterns. Another important circumstance is that all normal mode frequencies of the corresponding linearized theory are integer in appropriate units. Such a fully resonant spectrum ensures that nonlinearities can produce significant effects over long times for arbitrarily small perturbation amplitudes, generating highly nontrivial weak field dynamics. The weak field dynamics of the system we are considering is closely related to investigations of nonlinear stability of anti-de Sitter AdS) spacetime initiated in [6] see [7,8] for brief reviews and further references). A conformal transformation relates our considerations to the dynamics of the cubic Klein Gordon equation with mass m 2 = 2 conformally coupled self-interacting probe scalar) on the four-dimensional AdS space. In contrast to the bulk of AdS stability research, we do not consider the backreaction of perturbations on the AdS geometry. However, a number of features of the weak field dynamics remain unaffected by our simplification. The weak field dynamics of a massless self-interacting probe scalar in AdS was introduced as a toy model for the AdS stability problem in [9].) Our main tool for analyzing the weak field dynamics, given a fully resonant spectrum of the linearized system, is the time averaging method that goes back to Bogoliubov and Krylov and is described in many books on perturbation theory e.g., [10]). For a contemporary treatment adapted to studies of nonlinear PDEs, see [11]. Applied to the PDE we are considering, the time averaging method produces a simplified infinite-dimensional dynamical system, which we shall refer to as the conformal flow. By standard theorems underlying time averaging, this effective system in PDE theory sometimes called the resonant system) accurately approximates the dynamics of the original PDE for small amplitude fields of order Oε) on long time-scales of order O1/ε 2 ), and in particular it accurately captures the energy transfer patterns on such time scales. Time averaging and related techniques were introduced to the AdS stability problem closely related to our current studies in [12 14].) The conformal flow displays a number of highly special features not apparently present in the original PDE and suggestive of integrability. This includes additional conserved quantities, low-dimensional invariant subspaces, a variety of stationary states without any energy transfer, as well as solutions with exactly periodic energy transfer patterns. Our original PDE can be viewed as an infinite system of oscillators with integer frequencies and a quartic potential. We do not see any indications that it might be integrable; only the time-averaged system describing its weak field dynamics displays signs of integrability. While extensive searches for quartic integrable systems i.e., integrable mechanical analogs of our original PDE) have been conducted in the past see [15]for a review and references to earlier work and [16,17] for some later contributions), we are

3 Conformal Flow on S 3 and Weak Field Integrability in AdS not aware of any broad scans for systems displaying weak field integrability of the type suggested by our results. The conformal flow introduced in this paper bears a strong resemblance to the cubic Szegő equation that has been designed and studied in a series of papers by Gérard and Grellier, of which [18 21] are particularly relevant for our purposes. This equation also emerges as the resonant approximation of a nonlocal transport equation, called the cubic half-wave equation [19], which is a special case of one-dimensional models of wave turbulence introduced in [22]. In [18 21] and related publications, a number of remarkable properties were established for the cubic Szegő equation, in particular, a Lax pair structure, existence of finite-dimensional invariant subspaces, as well as some weak turbulence phenomena. The bulk of our effort will be directed at recovering analogs of a subset of these results for the conformal flow, a system with a natural geometric origin. We furthermore believe that the analogy goes deeper than what we are able to explicitly demonstrate in this article. There are also close similarities between the conformal flow and resonant approximations for the cubic Schrödinger equation, either with harmonic trapping [23] or in the infinite volume limit on the 2-dimensional torus [24,25]. Another system that shares some properties with the conformal flow and the cubic Szegő equation is the Lowest Landau Level evolution equation, which appears in the studies of rapidly rotating Bose- Einstein condensates see [26] and references therein). We furthermore mention the curious case of the Charney Hasegawa Mima CHM) equation explored in [11]. Also in that case, effective integrable dynamics emerges in the time-averaged approximation which the authors call the effective equation ). However, a substantial difference is that the linearized spectrum of the CHM equation is not fully resonant and splits into disconnected resonant clusters. In such a situation, the resulting effective time-averaged system literally separates into independent finite-dimensional integrable subsystems. This simplification does not happen in our case because our linearized spectrum is fully resonant. We comment on how our work is related to extensive investigations of the last decade in the area of AdS/CFT Integrability see [27] for a review). One of the central ingredients of that line of research is integrability of sigma models describing string worldsheets) on target spaces involving AdS factors. While we are also presently talking about integrability emerging due to special features of the AdS geometry, there are also essential differences. We are considering a field system evolving in AdS, rather than a sigma model whose fields are maps from a two-dimensional surface into AdS). We are furthermore not talking about exact integrability of our AdS field system, but rather about the integrability of the effective weak field theory conformal flow) emerging from it within the standard time-averaged approximation. Our exposition is organized as follows. In Sect. 2, we formulate our scalar field dynamics on a three-sphere, describe its connection to AdS stability problems and apply time averaging to derive the corresponding effective system, which we name the conformal flow. In Sect. 3, we give a very elementary pragmatic summary of some properties of the cubic Szegő equation relevant in our context, aimed mostly at the physics audience. In Sect. 4, we construct a three-dimensional invariant subspace of the conformal flow and discuss some complex plane properties of the generating functions that encode the conformal flow amplitudes. In Sect. 5, we analyze stationary states of the conformal flow, in which the initial conditions are adjusted to ensure that no energy transfer occurs. In Sect. 6, we exhibit explicit solutions with periodic energy flows on the three-dimensional invariant subspace of the conformal flow and discuss their properties. Finally, in Sect. 7,

4 1182 P. Bizoń, B. Craps, O. Evnin, D. Hunik, V. Luyten, M. Maliborski we give a summary of what we have practically demonstrated and an outlook on more ambitious directions, keeping the likelihood of Lax-integrability in mind. 2. Conformally Coupled Scalar Field on the Einstein Cylinder and in AdS 4 and Its Weak Field Dynamics As a simple model of confining geometry we consider the Einstein cylinder, which is the globally hyperbolic spacetime M, g) with topology R S 3 and metric g = dt 2 + r 2 dx 2 +sin 2 xdω 2), 1) where x [0,π], dω 2 is the round metric on the unit 2-sphere and r is the radius of S 3. This spacetime has a constant scalar curvature Rg) = 6/r 2. As a model of nonlinear dynamics on M, we consider the semilinear wave equation for the real scalar field φ : M R g 16 Rg) ) φ λφ 3 = 0, 2) where g := g μν μ ν is the wave operator associated with g and λ is a constant. For concreteness, we assume that λ>0 which corresponds to a defocusing nonlinearity); however, in the small data regime that we focus upon here, all the results below hold true in the focusing case λ<0aswell. Due to the identity 2 g 1 ) 6 R 2 g) 1 φ) = 3 g 16 ) Rg) φ, 3) Eq. 2) is conformally invariant. After rescaling t t/r and φ r λφ,ittakesthe dimensionless form φ tt S 3φ + φ + φ 3 = 0, 4) which can be interpreted as the cubic Klein Gordon equation on the unit S 3 with a unit mass. For simplicity, we shall restrict our analysis to rotationally symmetric fields, which depend only on t, x). Substituting φt, x) = vt, x)/ sinx) into 4) we obtain a nonlinear string equation v tt v xx + sin 2 x = 0, 5) with Dirichlet boundary conditions vt, 0) = vt, π) = 0, which are enforced by regularity of φt, x) on M. Remark 1. We could have arrived at the same equation considering the conformally coupled self-interacting scalar 2) in the 4-dimensional anti-de Sitter AdS 4 ) spacetime with the metric g = 1 cos 2 dt 2 + dx 2 +sin 2 xdω 2), 6) x where x varies between 0 and π/2. Indeed, g = cos 2 x g, hence it follows from 3) that the field redefinition φt, x) = φt, x) cos x converts Eq. 2) from the Einstein cylinder to AdS. In order to define the evolution, one has to impose a boundary condition at x = π/2. Here, we shall impose the Dirichlet condition φπ/2) = 0 on the equator, v3

5 Conformal Flow on S 3 and Weak Field Integrability in AdS corresponding to the reflecting boundary conditions most commonly used in the AdS research context. This model is just a subsector of the model on the full Einstein cylinder, which can be implemented by imposing the reflection symmetry φt, x) = φt,π x). We emphasize this connection because problems involving nonlinear dynamics of small AdS perturbations have received a significant amount of attention in recent years, and our objectives here have much in common with that body of work. Our goal is to understand the evolution of small smooth initial data v0, x), t v0, x). Decomposing the solution into a Fourier series, vt, x) = c n t) sinn +1)x, 7) n=0 we get from 5) an infinite system of coupled oscillators d 2 c n dt 2 + n +1)2 c n = jkl S jkln c j c k c l, 8) with the interaction coefficients S jkln = 2 π dx π 0 sin 2 sin j +1)x sink +1)x sinl +1)x sinn +1)x. 9) x To factor out fast linear oscillations in 8), we change the variables using variation of constants c n = β n e in+1)t + β n e in+1)t, 10) dc n = in +1) β n e in+1)t β n e in+1)t). 11) dt This transforms the system 8)into 2in +1) dβ n dt = jkl S jkln c j c k c l e in+1)t, 12) where each c j in the sum is given by 10). Thus, each term in the sum has a factor e i t, where = n +1) ± j +1) ± k +1) ± l +1), with all the three plus-minus signs independent. The terms with = 0 correspond to resonant interactions, while those with = 0 are non-resonant. Passing to slow time τ = ε 2 t and rescaling β n t) = εα n τ), we see that for ε going to zero the non-resonant terms e i τ/ε2 are highly oscillatory, and therefore expected to be negligible in a sense we are about to specify). Keeping only the resonant terms in 12) which is equivalent to time-averaging), we obtain an infinite autonomous dynamical system 2in +1) dα n dτ = 3 S jkln α j α k ᾱ l, 13) jkl where the summation runs over the set of indices { j, k, l} for which = 0. Evaluating the integrals for the S-coefficients, one finds that this set reduces to { jkl j + k l = n} and for such resonant combinations of indices, S jkln = min j, k, l, n) )

6 1184 P. Bizoń, B. Craps, O. Evnin, D. Hunik, V. Luyten, M. Maliborski Note that terms with n = j + k + l + 2 could have been present in principle, but the corresponding S-coefficients vanish, which can be verified by direct calculation. This is directly parallel to the selection rules that have been extensively discussed in the AdS stability literature [13,14,28,29]. It is part of the standard lore in nonlinear perturbation theory that solutions to 12) starting from small initial data of size Oε) are well approximated by solutions of 13) on timescales of order Oε 2 ). More precisely, if β n t) and α n τ) are solutions of 12) and 13), respectively, and β n 0) = εα n 0) for each n, then β n t) εα n τt)) Oε 2 ) for t Oε 2 ). In other words, on this timescale the dynamics of solutions of Eq. 5) is dominated by resonant interactions. Straightforward proofs normally phrased for systems with a finite number of degrees of freedom) can be found in textbooks on nonlinear perturbation theory, e.g., [10]. For convenience, we shall henceforth rescale τ to remove the numerical factors 2 and 3 from13). Renaming so rescaled τ back to t and using =d/dt, we finally arrive at the following system, which we call the conformal flow in +1) α n = [minn, j, k, n + j k) +1]ᾱ j α k α n+ j k. 15) Remark 2. Note that the AdS formulation described in Remark 1 simply corresponds to setting all even-numbered modes to 0, while only keeping odd-numbered modes α 2m+1 with m running from 0 to infinity. It is easy to see that if one implements this constraint in 15) and expresses everything through the new mode counting index m, one gets back Eq. 15). Thus, the conformal flow is equally applicable to the dynamics in R S 3 and in AdS 4. The conformal flow 15) is Hamiltonian with H = [minn, j, k, n + j k) +1]ᾱ n ᾱ j α k α n+ j k 16) n=0 and symplectic form n 2in +1) dᾱ n dα n : in +1) α n = 1 H. 17) 2 ᾱ n It enjoys the following one-parameter groups of symmetries where λ, θ are real parameters): Scaling: α n t) λα n λ 2 t), 18) Global phase shift: α n t) e iθ α n t), 19) Mode-dependent phase shift: α n t) e inθ α n t). 20) The latter two symmetries respect the Hamiltonian structure and give rise to two conserved quantities in addition to the Hamiltonian itself): Charge: Q = n +1) α n 2, 21) n=0

7 Conformal Flow on S 3 and Weak Field Integrability in AdS Linear energy: E = n +1) 2 α n 2. 22) We emphasize that these two conservation laws have no obvious counterparts in the original wave equation 4). Their prototypes were presented in the AdS stability context in [9,14]. n=0 3. An Elementary Introduction to the Cubic Szegő Equation The main content of our current treatment will be in presenting a range of remarkably simple dynamical behaviors of the conformal flow system 15). In anticipation of these results, it is useful to contemplate for a moment the apparent special features of Eq. 15) itself. One can trivially rewrite 15)as i α n = minn, j, k, n + j k) +1 ᾱ j α k α n+ j k. 23) n +1 One notable property of this representation is that many coefficients of the trilinear form on the right hand side are simply 1. Indeed, this happens for the contribution of any triplet of modes with all frequencies higher than the recepient mode n. It turns out useful to consider the following simpler equation, in which all numerical coefficients are set to 1, i α n = ᾱ j α k α n+ j k. 24) We shall see that, while the systems 23) and 24) are distinct, their dynamics is qualitatively similar, with a number of parallels between the emerging algebraic structures. Actually, the simplified system 24) is the Fourier representation of a simple paradifferential equation that has been studied before. To show this, let ut, e iθ ) be a complex function on the circle whose nonnegative Fourier coefficients are given by α n t) and negative Fourier coefficients vanish, i.e. ut, e iθ ) = α n t)e inθ, 25) n=0 and assume that n=0 α n 2 <. Mathematically, this means that u belongs to the Hardy space on the circle L 2 + S1 ) L 2 S 1 ) [30]. Furthermore, we define an orthogonal projector : L 2 S 1 ) L 2 + S1 ) called the Szegő projector) that filters out negative frequencies ) α n e inθ := α n e inθ. 26) n= n=0 Then, it is easy to see that the system 24) is equivalent to i t u = u 2 u). 27)

8 1186 P. Bizoń, B. Craps, O. Evnin, D. Hunik, V. Luyten, M. Maliborski This equation was introduced by Gérard and Grellier in [18] under the name of the cubic Szegő equation. They showed that this equation has a remarkably deep structure, including a Lax pair, finite-dimensional invariant subspaces and weakly turbulent solutions. We shall demonstrate below that the conformal flow exhibits at least some of this structure in a modified form), and conjecture that the similarity goes even further. In view of these parallels, we shall present below a very elementary summary of the features of the cubic Szegő equation that are important to us. Readers interested in more details are referred to the original publications [18 21]. The cubic Szegő equation is Hamiltonian with H = n=0 ᾱ n ᾱ j α k α n+ j k = 1 2π 0 2π u 4 dθ 28) and symplectic form n 2idᾱ n dα n, and has the same symmetries 18 20) as15). Hence, it possesses two conserved quantities analogous to 21) and 22), in addition to the Hamiltonian itself: M = P = n=0 n=0 α n 2 = 1 2π n α n 2 = 1 2π 2π 0 2π 0 u 2 dθ, 29) i θ u)ūdθ. 30) Borrowing terminology from studies of the nonlinear Schrödinger equation, we shall refer to these charges as mass and momentum, respectively. Note that M + P = u 2 H 1/2, where H s L 2 + is the Sobolev space on the circle equipped with the norm u 2 H s = n +1) 2s α n t) 2. 31) n=0 This fact can be used to control the growth in time of higher Sobolev norms and thereby prove that the Cauchy problem for the Szegő equation is globally well posed for smooth initial data in H s with s > 1/2 [18]. Any function ue iθ ) of the form 25)inL 2 + S1 ) can be identified with a holomorphic function uz) inside the unit disk z < 1 in the complex plane [30], hence we can write ut, z) = α n t)z n. 32) n=0 In terms of ut, z), the cubic Szegő equation 27) takes the form i t ut, z) = 1 dw 2πi w z ũt,w)ut,w)2, 33) where w =1 ũt, z) = ᾱ n z n. 34) n=0

9 Conformal Flow on S 3 and Weak Field Integrability in AdS General properties of the Cauchy-type integral appearing in 33) suggest the existence of meromorphic solutions for u with a finite number of isolated poles moving in the complex plane outside the unit disk. Indeed, it was shown in [18] that the cubic Szegő equation admits finite-dimensional invariant subspaces given by rational functions of arbitrarily high degree. Moreover, due to the Lax-integrability and, consequently, infinitely many conserved quantities), the Szegő flow restricted to such subspaces is completely integrable in the sense of Liouville. The simplest invariant manifold with nontrivial dynamics is parametrized by the following single pole ansatz bt) + at)z ut, z) =, p < 1. 35) 1 pt)z We shall now demonstrate explicitly by elementary means that the subspace given by such functions is indeed dynamically invariant under the flow, and analyze its evolution. The Fourier coefficients corresponding to 35) are α 0 = b and α n = a + bp)p n 1 for n 1. 36) Expressed through the parameters of our ansatz, the mass and momentum take the form M = b 2 a + bp 2 a + bp 2 +, P = 1 p 2 1 p 2 ) 2. 37) Substitution of 35) into33) produces exactly three distinct dependences on z on both sides, namely, 1 pz) 2, 1 pz) 1 and a z-independent term, resulting in the following three equations for at), bt), and pt) iȧ = Ma, 38) iḃ = M + P)b + Pa p, 39) i ṗ = Mp + a b, 40) where we used the conservations laws 37) to simplify some terms. These equations describe a three-dimensional Liouville-integrable system which can be easily solved in [20] an explicit formula for the general solution was derived using harmonic analysis tools). One first integrates 38) to find a, whereupon the remaining equations become linear. The solution corresponding to initial conditions a0) = a 0, b0) = b 0, p0) = p 0 for simplicity assumed to be real) is at) = a 0 e imt, 41) bt) = b 0 cos ωt i b ) 0M + P) +2a 0 p 0 P sin ωt e 2 i M+P)t, 42) 2ω pt) = p 0 cos ωt i p ) 0M + P) +2a 0 b 0 sin ωt e 2 i M P)t, 43) 2ω where ω = 2 1 M + P) 2 4Pa0 2 and the expressions for M and P in terms of a 0, b 0, p 0 ) are given in 37). Thus, the solution is quasiperiodic for all initial conditions, while the corresponding mass spectrum, α 0 2 = M P1 p 2 ) and α n 2 = P1 p 2 ) 2 p 2n 1) for n 1, 44)

10 1188 P. Bizoń, B. Craps, O. Evnin, D. Hunik, V. Luyten, M. Maliborski is expressible through p 2 and the conserved quantities, and hence exactly periodic. In the special case M = P = a0 2,wehaveω = 0 and then at) = a 0 e imt, bt) = a 0 p 0 e imt, pt) = p 0. 45) The evolution of two-mode initial data a 0 = 1, b 0 = 2ε, p 0 = 0 was used in [18] to illustrate a very interesting instability phenomenon. For these data, 43) gives i pt) = sinωt) 1+ε 2 e 2iε2t, ω = 2ε 1+ε 2, 46) hence, in the regime of ε 0, pt n ) 1 ε 2 /2 for a sequence of times t n = nπ 2ω, meaning that the energy goes periodically to arbitrarily high modes. Put differently, although all solutions are quasiperiodic, their radius of analyticity is not uniformly bounded from below. Note that 46) implies instability of the one-mode stationary state ut, z) = a 0 e imt z. As mentioned above, one can consider more general meromorphic solutions with an arbitrary number of time-dependent simple poles outside the unit disk ut, z) = N k=1 b k t) 1 p k t)z, p k < 1, N 2. 47) One can convince oneself that the ansatz is consistent in the sense that substituting it in 33) results only in z-dependences of the form 1 p k z) 1 and 1 p k z) 2 on both sides, producing 2N ordinary differential equations for the 2N functions b k and p k.itis possible to consistently impose p N = 0, producing an ansatz with one parameter less. The single pole ansatz 35) is precisely of such a form. The equations of motion within the subspace defined by 47) were written down explicitly in [18]. As in the simplest case of the three-dimensional invariant subspace discussed above, the dynamics on all finite-dimensional invariant subspaces is Liouville-integrable and bounded, and hence quasiperiodic. Among the quasiperiodic solutions there exist special solutions with time-independent amplitudes α n. Such solutions, which we call stationary states, were classified in [18] in the case of finite dimensional invariant manifolds. They are given by either finite Blaschke products or ut, z) = Nω = ut, z) = ce i c 2 t N k=1 p k z 1 p k z, p k < 1, c C, 48) cz l 1 p N z N e iλt, pt) = p0)e iωt c 2, λ = 1 p 2N ) 2, c 2, 49) 1 p 2N where N 1 and l N 1 are nonegative integers. It was recently shown in [21] that outside the finite-dimensional invariant subspaces there exists a dense set of smooth solutions whose radius of analyticity tends to zero for a sequence of times t n. Consequently, the Sobolev norms with s > 1/2 are

11 Conformal Flow on S 3 and Weak Field Integrability in AdS unbounded. Such a weakly turbulent behavior is somewhat surprising in a completely integrable model. The coexistence of integrability and turbulence is possible because the infinitely many conserved quantities are too weak to control higher regularity properties of solutions. 4. Three-Dimensional Invariant Subspaces of the Conformal Flow We now return to the conformal flow 15) and attempt to treat it along the lines applied above to the cubic Szegő equation. While the treatment is closely parallel, we shall be able to immediately recover only a part of the results available for the cubic Szegő equation. Using the generating function ut, z) and its conjugate ũt, z) defined as in 32) and 34), we find that the system 15) is equivalent to the integro-differential equation i t z zu) = 1 2πi w =1 ) dw wuw) zuz) 2 w ũw). 50) w z The following summation formula has been used while obtaining this complex plane representation: [minn, j, k, n + j k) +1]ρ j θ k = = 1 θ n+1 1 θ)1 ρ)1 θρ) 1+θ + + θ n 1 ρ)1 θρ). 51) Derivations of 50) and 51) and further comments are given in Appendix A. Although Eq. 50) looks like a somewhat more elaborate version of 33), we do not see an immediate way to characterize its meromorphic solutions and defer it to future work. Nonetheless, the lowest-dimensional nontrivial invariant subspace is easy to construct, as we shall explicitly demonstrate now. The ansatz relevant for the three-dimensional invariant subspace, analogous to 36), is α n = b + an)p n, 52) where b, a, p are complex-valued functions of time. The corresponding generating function has the following combination of poles outside the unit disk ut, z) = bt) at) 1 pt)z + at) bt) + at) bt))pt)z = 1 pt)z) 2 1 pt)z) 2. 53) We note that the two-mode initial data α n = 0forn 2) are accommodated within this ansatz as a special limiting case p 0 with b and a b)p finite. All of our general statements about the dynamics on the three-dimensional invariant subspace apply to solutions starting with such initial configurations. While we could have used the complex plane representation 50) to establish the validity of our ansatz, it is instructive to apply brute force summations in this case. Substitution of 52) into15) yields

12 1190 P. Bizoń, B. Craps, O. Evnin, D. Hunik, V. Luyten, M. Maliborski i1+n) ḃ + ȧn + nb + an) ṗ ) p = [minn, j, k, n + j k) +1] b + āj)b + ak)b + an + j k)) p 2 j. Note that the p n factor has consistently cancelled on the two sides. It remains to show that different n-dependences on both sides can be matched and produce a sufficiently small number of equations. All summations in 54) can be performed by applying the relation [minn, j, k, n + j k) +1] j K k L p 2 j = ) K L 1+θ + + θ n ρ ρ θ θ ) 1 ρ)1 θρ) 54), 55) θ=1 ρ= p 2 which follows from 51). While we give explicit expressions for the sums involved in 54) in Appendix B, the only thing one needs to know about these sums to establish the validity of our ansatz is that they are all polynomials of degree L in n, times n +1), the latter factor coming from Faulhaber s sums l + + n l, with l L, originating from 1 + θ + + θ n in 55). Therefore, explicit counting tells us that, upon substituting the summation formulas of Appendix B in 54), the factor n +1) will cancel on the two sides, leaving behind a statement that two quadratic polynomials in n equal each other. Matching the coefficients of these polynomials produces three ordinary differential equations for three functions b, a, p, confirming the validity of our ansatz. The explicit equations for b, a and p are given by i ṗ 1+y) 2 = p 6 iȧ 1+y) 2 = a 6 ) 2y a 2 + ba, 56) ) 5 b y 2 +4y) a 2 + 6y 1) ba +10yāb, 57) iḃ 1+y) 2 = b b 2 + 6y 2 +2y) a 2 +2ybā ) ) + a 2y b 2 + 4y +2)y 2 a 2 + y 2 ba, 58) where we have introduced the following notation, which will turn out to be useful later: y = p 2 1 p 2. 59) The conservation laws 21) and 22) take the form Q = 1+y) 2 b 2 +4y Re ba) +2y3y +1) a 2), 60) E = 1+y) 2 1+2y) b 2 +4y3y +2) Re ba) +4y6y 2 +6y +1) a 2). 61)

13 Conformal Flow on S 3 and Weak Field Integrability in AdS We remark that if ut, z) is a solution of the conformal flow, so is z N ut, z N+1 ) for any non-negative integer N. This automatically generates an infinite number of other three-dimensional invariant subspaces of the conformal flow, in which only subsets of modes are activated. We shall return to equations 56 58) in sect. 6, and analyze the dynamics within the three-dimensional dynamically invariant subspace 52) explicitly. Before proceeding in that direction, we shall discuss special solutions for which α n are time-independent. 5. Stationary States The conformal flow 15) admits solutions of the form α n t) = A n e iλ nt, 62) where the frequencies λ n and complex amplitudes A n are time-independent, and λ n are linear in n, that is λ n = λ nω for some real λ and ω. For such solutions there is no energy transfer between the modes, hence we call them stationary states. The cubic Szegő equation admits a variety of stationary states [18]; such solutions have also been considered in the context of AdS stability [12,31,32], where they were referred to as quasiperiodic solutions. Substituting 62) into15) we get a nonlinear eigenvalue problem n +1)λ nω)a n = [minn, j, k, n + j k) +1] Ā j A k A n+ j k. 63) The simplest solutions of this algebraic system, easily seen by inspection, A n = c δ Nn, λ = c 2, ω = 0forc C, give the one-mode stationary states for any non-negative integer N α n = c δ Nn e i c 2t. 64) We will see shortly that these trivial solutions are the endpoints of two-parameter families of stationary solutions. We note in passing that they are seeds for time-periodic solutions of the original equation 5) whose construction will be described elsewhere along the lines of [33]. Within the ansatz 52), stationary states take the form bt) = b0)e iλt, at) = a0)e iλt, pt) = p0)e iωt. 65) Plugging this into the system 56 58), we obtain a system of algebraic equations that can be solved explicitly. The case a0) = 0 yields For this solution Q = reads bt) = c exp i c 2 t 1 p 2 ) 2 ). 66) c 2 1 p 2 ) 2, hence λ = Q. The corresponding generating function ut, z) = c 1 pz e iλt. 67)

14 1192 P. Bizoń, B. Craps, O. Evnin, D. Hunik, V. Luyten, M. Maliborski For nonzero a we get the following two-parameter family of stationary states with ω = 0 b0) = 2c p 2, a0) = c1 p 2 ), λ = c 2 p 2 1 p 2 ) 2 68) For this solution Q = 2 c 2 p 2 1 p 2 ) 2, E = 4 c 2 p 2 1+ p 2 ) 1 p 2 ) 3, 69) hence λ = 2 1 Q. The corresponding generating function reads ut, z) = c 2 p p 2) pz 1 pz) 2 e iλt. 70) In addition, for p p := , there is a pair of two-parameter families of stationary states with nonzero ω. The range of p is restricted by the condition that κ := p 4 14 p should be real, so we have a tangential bifurcation at p. These two solutions labelled by ±) aregivenby b ± 0) = c1+5 p 2 ± κ), a ± 0) = 2c1 p 2 ), 71) ω ± = c 2 1+ p 2 ± κ 3 1 p 2, 72) λ ± = 2 c p 2 ) 3 1 p 2 ± 3+4 p 2 )κ 1 p 2 ) 2, 73) and their conserved quantities are The corresponding generating function reads Q ± = 6 7 λ ± + ω ± ), E ± = 6ω ±. 74) ut, z) = c 1+5 p 2 ± κ) p 2 ± κ)pz 1 pz) 2 e iλ ±t, 75) where pt) = p0)e iω ±t. Outside the three-dimensional invariant subspace given by 52) there exist many other stationary states. For instance, we verified that, in close parallel to the cubic Szegő equation, any finite Blaschke product ut, z) = ce i c 2 t N k=1 p k z 1 p k z yields a stationary state. There are also stationary states where only every Nth mode is activated, for instance ut, z) = 76) czn 1 c 2 1 p N z N e iλt, λ = 1 p 2N ) 2. 77)

15 Conformal Flow on S 3 and Weak Field Integrability in AdS Remark 3. For all the above stationary states with ω = 0wehaveλ = Q N+1, where N is the number of zeroes of the generating function counted with multiplicity). It would be interesting to find a reason that underlies this quantization structure. It is well known that stationary states of Hamiltonian systems admit a variational characterization. In our case, it follows from the Hamilton equation 17) that stationary states 62) are the critical points of the functional K := 1 H λq + ωe Q), 78) 2 hence stationary states with ω = 0areextremaofH for fixed Q, while stationary states with nonzero ω are extrema of H for fixed Q and E. This fact is very helpful in determining stability properties of stationary states, as will be described elsewhere. 6. Dynamics on the Three-Dimensional Invariant Subspace We now return to equations 56 58) describing the dynamics on the three-dimensional invariant subspace of the conformal flow and demonstrate how to solve them explicitly. Once again, the situation is closely parallel to the three-dimensional invariant subspace of the cubic Szegő equation with its underlying one-dimensional periodic motion. From 56), we obtain for y defined by 59) This equation and 57) imply that ẏ = 1 3 y1+y)3 Im ba). 79) S = a 2 y1+y) 3 80) is conserved. This quantity is related to the conformal flow Hamiltonian 16)by H = Q 2 2S 2. 81) Equations 60), 61) and 80) can be resolved to express b 2, a 2 and Re ba) through Q, E, S and y as follows b 2 2Q E +3yQ +2S) = 1+y) 3, a 2 = E Q 2S 2yQ +6S) Re ba) = 4y1+y) 3. From these relations, together with 79), one obtains ) y y + ẏ 2 = Q2 +12S 2 36 EQ +2S) Q 2 +12S 2 S y1+y) 3, 82) ) E Q 2S)2 4Q 2 +12S 2. 83) ) Note that this equation has the algebraic form of energy conservation for an ordinary one-dimensional harmonic oscillator. The solution of 83) reads yt) = B + A sin t + ψ), = 1 Q 2 +12S 2) 1/2, 84) 6

16 1194 P. Bizoń, B. Craps, O. Evnin, D. Hunik, V. Luyten, M. Maliborski where B and A are constants depending on E, Q, and S; ψ is a phase that must be determined from the entire set of initial conditions. By Eq. 82), this exactly periodic motion is transferred to b 2, a 2 and Re ba), and hence to the mode energy spectrum α n 2 = b + na 2 p 2n = b 2 +2nRe ba) + n 2 a 2) p 2n. 85) One will thus observe exact returns of the energy spectrum to the initial configurations for all solutions within the three-dimensional invariant subspace of the conformal flow. The turning points of the periodic motion of y described by 84), y ± = B ± A, are given by the roots of the quadratic polynomial appearing on the right-hand side of 83) the special case A = 0 corresponds to stationary states). In terms of the energy spectrum 85), these turning points provide lower and upper bounds for the inverse and direct cascades of energy, respectively. One of the key questions in this context is how large can y grow starting from some small y. To answer this question, we consider the ratio 1+y + 1+y + ) 2 = 1+y 1+y )1+y + ) 1+y + + y ) 2. 86) 1+y + y + + y y + An advantage of the last representation is that it only contains combinations of y + and y directly expressible through the coefficients of the quadratic polynomial on the righthand side of 83). One thus gets 1+y + + y ) 2 4E 2 Q +2S) 2 4Q +2S)2 = 1+y + y + + y y + Q 2 +12S 2 )E + Q +2S) 2 Q 2 +12S 2. 87) We furthermore notice that Q = 1+y) 2 b +2ya 2 +2S 2S. 88) Therefore, dividing the numerator and denominator of 87) byq 2 and replacing the numerator by its maximum and denominator by its minimum, we obtain a simple uniform upper bound which very likely can be tightened with extra work) 1+y + 1+y ) This proves that the transfer of energy to high frequencies or, equivalently, the growth of the higher Sobolev norms 31)) is uniformly bounded. We recall from Sect. 3 that there is no such bound for the Szegő flow on the three-dimensional invariant subspace 35). In this sense, the conformal flow is much less turbulent than the Szegő flow. Heuristically, this is not very surprising in view of the fact that the interaction coefficients in 23) minn, j, k, n + j k) +1 n +1 90) decay when the recipient mode number n is large and at least one of the source mode numbers is small, undermining the efficiency of energy transfer from low to high frequencies. This is in contrast to the Szegő flow where all coupling coefficients are equal to 1.

17 Conformal Flow on S 3 and Weak Field Integrability in AdS Summary and Open Questions Starting with a naturally defined geometric PDE 4) on a three-sphere describing a self-interacting conformally coupled scalar field, we have considered its effective timeaveraged weak field dynamics, arriving at the conformal flow 15). The conformal flow is both structurally similar and displays a number of dynamical parallels to the previously known cubic Szegő equation 24). Building on these analogies, we have revealed the dynamics on the three-dimensional invariant subspace of the conformal flow described by the ansatz 52). Within this subspace, the dynamics is Liouville-integrable and bounded, and hence quasiperiodic, with exactly periodic energy flows 85) displaying an alternating sequence of direct and inverse cascades. Unlike the cubic Szegő case, the conformal flow dynamics within the three-dimensional invariant subspace does not display turbulent behaviors not even the weak weak turbulence of [3]). Other similar three-dimensional subspaces can be immediately constructed, as follows from simple properties of the conformal flow equation. There are furthermore direct parallels between the structure of stationary states of the Szegő and conformal flows these states are special solutions for which no energy transfer between the modes occurs). Of these, the Blaschke product states 76) are particularly intriguing, as they clearly demonstrate that the parallels between the Szegő and conformal flows extend beyond the three-dimensional invariant subspaces that have been our main focus here. We conclude with a list of open questions that we leave for future investigations: Is the conformal flow 15) Lax-integrable? Is there a way to re-express the conformal flow through a projector similar to the Szegő projector? Apart from the conceptual importance, such a projector should make calculations easier. The complex plane representation of the conformal flow 50) can be seen as a first step in this direction. Are there higher dimensional invariant subspaces? The three-dimensional subspace we have considered, with its double-pole generating function 53), and parallels to the cubic Szegő case strongly suggest that more elaborate meromorphic ansätze should work; however, we have not been able to find them. Are there weakly turbulent solutions? We have explicitly demonstrated that no turbulent behavior occurs within the three-dimensional invariant subspace 52). This does not in principle exclude turbulence for other initial conditions; however, we find it unlikely because the conformal flow generally appears to display less energy spread than the Szegő flow. Complete classification of stationary states and analysis of their stability properties would contribute to the overall picture of the conformal flow dynamics. To what extent can the properties of the conformal flow, both explicitly demonstrated and putative, be structurally stable with respect to variations of the flow equation? Do these dynamical patterns have applications to the weak field limit of other related systems? Can the Szegő and conformal flows be just two representative members of a large hierarchy of equations? To what extent can the highly structured dynamics of the conformal flow be transferred to the original conformally invariant cubic wave equation 5) on the threesphere? Standard results on time averaging guarantee that in the weak field regime the conformal flow accurately approximates our original wave equation for a long but limited time. Can our findings have further implications on longer time scales? Acknowledgements. P.B. gratefully acknowledges the stimulating atmosphere at the trimestre Ondes Non Linéaires in IHES, where part of this work was done; he is particularly indebted to Patrick Gérard for

18 1196 P. Bizoń, B. Craps, O. Evnin, D. Hunik, V. Luyten, M. Maliborski helpful discussions. The work of P.B., D.H., and M.M. was supported by the Polish National Science Centre grant no. DEC-2012/06/A/ST2/ P.B. and M.M. gratefully acknowledge the support of the Alexander von Humboldt Foundation. D.H. acknowledges a scholarship of Marian Smoluchowski Research Consortium Matter Energy Future from KNOW funding. The work of B.C. and V.L. is supported in part by the Belgian Federal Science Policy Office through the Interuniversity Attraction Pole P7/37, by FWO-Vlaanderen through projects G020714N and G044016N, and by Vrije Universiteit Brussel VUB) through the Strategic Research Program High-Energy Physics. V.L. is supported by a PhD fellowship from the VUB Research Council. B.C. thanks the organizers of the NumHol2016 meeting at Santiago de Compostela and of the Nordita program Black Holes and Emergent Spacetime for hospitality while this work was in progress. The work of O.E. is funded under the CUniverse research promotion project by Chulalongkorn University grant reference CUAASC). O.E. would like to thank Latham Boyle, Stephen Green and especially Luis Lehner for stimulating discussions and warm hospitality during his visit to the Perimeter Institute Waterloo, Canada). The part of this research conducted during that visit was supported by the Perimeter Institute for Theoretical Physics. Research at Perimeter Institute is supported by the Government of Canada through the Department of Innovation, Science and Economic Development and by the Province of Ontario through the Ministry of Research and Innovation.) Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original authors) and the source, provide a link to the Creative Commons license, and indicate if changes were made. A. Complex Plane Representation for the Conformal Flow We first introduce the generating functions u and ũ according to 32) and34). It is true for any contour enclosing the origin but not enclosing the singularities of u, and any contour enclosing all singularities of ũ that α m = 1 2πi dz z uz)z m, ᾱ m = 1 2πi dz z ũz)zm. 91) Substituting these expressions to the right-hand side of conformal flow 15), one gets 1 ds dv dw in +1) α n = 2πi) 3 ũs)uv)uw) s v wn+1 s ) j w ) k [minn, j, k, n + j k) +1]. 92) w v The summation can be performed using the following master formula [minn, j, k, n + j k) +1]ρ j θ k = 1 θ n+1 1+θ + + θn = 1 θ)1 ρ)1 θρ) 1 ρ)1 θρ). 93) This summation formula can be derived by many methods, for example, by noticing that 2minn, j, k, n + j k) = n + j k n k j, 94) and performing brute force resummations of geometric series. The best way to understand it, however, is essentially combinatorial in nature. The right-hand side of 93) can be rewritten as 1+θ + + θ n ) ρ m 1+θ + + θ m). 95) m=0 The coefficient of ρ j θ k in this expression is just the coefficient of θ k in the product 1+θ + + θ n ) 1+θ + + θ j ). 96)

19 Conformal Flow on S 3 and Weak Field Integrability in AdS Counting the number of relevant pairings of powers of θ in the two polynomials of the product reproduces the coefficient on the left-hand side of 93), thus verifying the summation formula. Using 93), one can rewrite 92)as in +1) α n = 1 2πi) 3 ds s dv v w n+1) v n+1) dw ũs)uv)uw) 1 s/v)1 s/w)1 w/v). 97) One must have s < v and s < w in order for the sum to converge.) Multiplying by z n with z < v and z < w ) and summing over n, one obtains after elementary simplification i t z zu) = 1 2πi ds [ 1 s ũs) 2πi ] uv) 2 dv. 98) 1 s/v)v z) Since by construction, the v-integration contour must not enclose any singularities of u, the integral inside the square brackets is simply given by a sum of two residues, resulting in the complex plane representation we have quoted in the main text: i t z zu) = 1 2πi ds sus) zuz) s ũs) s z ) 2. 99) B. A Few Summation Formulas We assemble below the explicit summation formulas necessary for deriving 56 58) from 54). A crucial feature of these expressions in our context is that all of them are proportional to n +1). [minn, j, k, n + j k) +1] p 2 j = [minn, j, k, n + j k) +1] j p 2 j = n +1 1 p 2 ) 2, 100) 2n +1) p 2 1 p 2 ) 3, 101) [minn, j, k, n + j k) +1]k p 2 j nn +1) n +1) p 2 = 21 p 2 + ) 2 1 p 2 ) 3, 102) [minn, j, k, n + j k) +1] j 2 p 2 j 2n +1) p 2 6n +1) p 4 = 1 p 2 + ) 3 1 p 2 ) 4, 103) [minn, j, k, n + j k) +1]k 2 p 2 j nn +1)2n +1) = 61 p 2 ) 2 + n +1)2 p 2 2n +1) p 4 1 p 2 + ) 3 1 p 2 ) 4, 104) [minn, j, k, n + j k) +1] jk p 2 j = n +1)2 p 2 3n +1) p 4 1 p + ) 3 1 p 2 ) 4, 105) [minn, j, k, n + j k) +1] j 2 k p 2 j = n +1)2 p 2 3n +1)n +3) p 4 12n +1) p 6 1 p 2 + ) 3 1 p 2 ) p 2 ) 5, 106) [minn, j, k, n + j k) +1] jk 2 p 2 j

20 1198 P. Bizoń, B. Craps, O. Evnin, D. Hunik, V. Luyten, M. Maliborski = n +1)2n2 +4n +3) p 2 n +1)3n +7) p 4 8n +1) p 6 31 p 2 ) p 2 ) p 2 ) ) References 1. Zakharov, V.E., L vov, V.S., Falkovich, G.: Kolmogorov Spectra of Turbulence I: Wave Turbulence. Springer, Berlin 1992) 2. Nazarenko, S.: Wave Turbulence. Springer, Berlin 2011) 3. Colliander, J., Keel, M., Staffilani, G., Takaoka, H., Tao, T.: Transfer of energy to high frequencies in the cubic defocusing nonlinear Schrödinger equation. Invent. Math. 181, ). arxiv: [math.ap] 4. Hani, Z.: Long-time instability and unbounded Sobolev orbits for some periodic nonlinear Schrödinger equations. Arch. Ration. Mech. Anal. 211, ). arxiv: [math.ap] 5. Guardia, M., Kaloshin, V.: Growth of Sobolev norms in the cubic defocusing nonlinear Schrödinger equation. J. Eur. Math. Soc. 17, ). arxiv: [math.ap] 6. Bizoń, P., Rostworowski, A.: On weakly turbulent instability of anti-de Sitter space. Phys. Rev. Lett. 107, ). arxiv: [gr-qc] 7. Bizoń, P.: Is AdS stable? Gen. Rel. Grav. 46, ). arxiv: [gr-qc] 8. Craps, B., Evnin, O.: AdS in)stability: an analytic approach. Fortsch. Phys. 64, ). arxiv: [gr-qc] 9. Basu, P., Krishnan, C., Saurabh, A.: A stochasticity threshold in holography and the instability of AdS. Int. J. Mod. Phys. A 30, ). arxiv: [hep-th] 10. Murdock, J.A.: Perturbations: Theory and Methods. SIAM, Philadelphia 1987) 11. Kuksin, S., Maiocchi, A.: The effective equation method. In: New Approaches to Nonlinear Waves. Springer, Berlin 2016). arxiv: [math-ph] 12. Balasubramanian, V., Buchel, A., Green, S. R., Lehner, L., Liebling, S. L.: Holographic thermalization, stability of anti-de Sitter space, and the Fermi-Pasta-Ulam paradox. Phys. Rev. Lett. 113, ). arxiv: [hep-th] 13. Craps, B., Evnin, O., Vanhoof, J.: Renormalization group, secular term resummation and AdS in)stability. JHEP ) 48. arxiv: [gr-qc] 14. Craps, B., Evnin, O., Vanhoof, J.: Renormalization, averaging, conservation laws and AdS in)stability. JHEP ) 108. arxiv: [gr-qc] 15. Lakshmanan, M., Sahadevan, R.: Painlevé analysis, Lie symmetries and integrability of coupled nonlinear oscillators of polynomial type. Phys. Rep. 224, ) 16. Bruschi, M., Calogero, F.: Integrable systems of quartic oscillators. Phys. Lett. A 273, ) 17. Bruschi, M., Calogero, F.: Integrable systems of quartic oscillators II. Phys. Lett. A 327, ). arxiv:nlin/ [nlin.si] 18. Gérard, P., Grellier, S.: The cubic Szegő equation. Ann. Scient. Éc. Norm. Sup. 43, ). arxiv: [math.cv] 19. Gérard, P., Grellier, S.: Effective integrable dynamics for a certain nonlinear wave equation. Anal. PDE 5, ). arxiv: [math.ap] 20. Gérard, P., Grellier, S.: An explicit formula for the cubic Szegő equation. Trans. Am. Math. Soc. 367, ). arxiv: [math.ap] 21. Gérard, P., Grellier, S.: The cubic Szegő equation and Hankel operators. arxiv: [math.ap] 22. Majda, A.J., McLaughlin, D.W., Tabak, E.G.: A one-dimensional model for dispersive wave turbulence. J. Nonlinear Sci. 7, ) 23. Hani, Z., Thomann, L.: Asymptotic behavior of the nonlinear Schrödinger equation with harmonic trapping. Commun. Pur. App. Math. 69, ). arxiv: [math.ap] 24. Faou, E., Germain, P., Hani, Z.: The weakly nonlinear large box limit of the 2D cubic nonlinear Schrödinger equation. J. Am. Math. Soc. 29, ). arxiv: [math.ap] 25. Germain, P., Hani, Z., Thomann, L.: On the continuous resonant equation for NLS: I. Deterministic analysis. J. Math. Pure Appl. 105, ). arxiv: [math.ap] 26. Germain, P., Thomann, L.: On the high frequency limit of the LLL equation. Quart. Appl. Math. 74, ). arxiv: [math.ap] 27. Beisert, N. et al.: Review of AdS/CFT integrability: an overview. Lett. Math. Phys. 99, ). arxiv: [hep-th] 28. Yang, I-S.: Missing top of the AdS resonance structure. Phys. Rev. D 91, ). arxiv: [hep-th] 29. Evnin, O., Nivesvivat, R.: AdS perturbations, isometries, selection rules and the Higgs oscillator. JHEP 1601, ). arxiv: [hep-th] 30. Garnett, J.B.: Bounded Analytic Functions. Springer, Berlin 2007)