BTZ dynamics and chaos

We find an effective action for gravitational interactions with scalars in $AdS_3$ to first order in $G_N$ at the conformal boundary. This action can be understood as an action for the Brown-Henneaux modes and is given by the square-root product of right and left moving Schwarzian derivatives for conformal transformations of the boundary. We thus reproduce the result $\lambda_L=\frac{2\pi}{\beta}$ for OTOC computed first in arXiv:1412.6087 for a Schwarzchild black hole in $AdS_3$. Applying the same procedure to rotating BTZ we find the Lyapunov index to be $\lambda_L=\frac{2\pi}{\beta_+}>\frac{2\pi}{\beta}$ where $\beta_+=\beta(1-\mu_L)$, with $\mu_L=\frac{r_-}{r_+}$ being the chemical potential for angular momentum. We thus comment on a possible modification to a part of the proof given in arXiv:1503.01409 to accommodate this result.


Introduction
A very intriguing phenomena of strongly coupled thermal systems is chaos. In a classical sense, phase space trajectories which differ in their initial values by a small amount tend to grow exponentially far apart at later times. A good quantum mechanical analogue of this would be the time scale when C(t) = [W (t), V (0)] 2 β (1.1) becomes equal to 2 W W β V V β . Here, W &V are simple Hermitian operators with O(1) degrees of freedom. The exponential increase in the time of C(t) ≈ e λ L t can be considered as the Lypunov index generically associated with chaotic systems [1]. The chaotic behaviour of thermal large N CFTs is related to the chaotic behaviour of black holes via the gauge-gravity duality, the latter are conjectured to be the fastest "scramblers" of information [2]. Shenker and Stanford [3] first computed the out-of -time-ordered (otoc) term in (1.1) holographically using the eikonal approximation. In [3] the next order in probe approximation in G N was computed for a 2 → 2 scattering for 2 minimally coupled scalars in AdS d Schwarzchild interacting with each other only via gravity. The Lyapunov index thus obtained was λ L = 2π/β, β being the inverse temperature of the AdS d Schwarzchild. This lead Maldacena, Shenker and Stanford [4] to propose a bound on the Lyapunov index of large N thermal QFTs to be λ L ≤ 2π/β by using generic arguments of unitarity and analyticity of Whightman functions on the complex plane. It was assumed that holographic CFTs saturate this bound as evidenced by [3,5].
Further interest in chaotic systems was heightened by the study of the SYK [6]and SYKlike models initiated first in [7] and [8]; in [8] Maldacena and Stanford found that the otoc for the fermions has λ L = 2π/β, this computation was done in the strong coupling (zero temperature) limit of the SYK model where the model is conformal. In order to compute the leading contribution to the 4pt. function they had to break the conformal invariance at zero temperature. The modes which are responsible for maximizing chaos where shown to be the modes related by diffeomorphism which now have an action due to breaking of conformal invariance. Their effective action was computed and found to be the Schwarzian derivative of reparametrizations of the thermal circle. Many interesting properties of the SYK model have since been uncovered [9][10][11][12][13][14][15][16] to quote a few.
This lead to investigations to ascertain the bulk degrees of freedom which are responsible for similar chaotic behaviour. The dynamics of near extremal black holes was found to be captured by a 2d dilaton-gravity theory of Jackiw [33] and Teitelboim [34] in [35]. The nAdS 2 dynamics of Jackiw-Teitelboim (JT) action is essentially dictated by its asymptotic symmetries since the theory possesses zero propagating degrees of freedom. The effective action for these modes was captured by the Schwarzian action for the AdS 2 boundary diffeomorphisms [36,37], similar ideas where pursued in [38]. Explicit computations on near extremal RN AdS 4 black holes [39] corroborated this understanding from a higher dimensional perspective.
The dynamics of the 2d gravity theory reproducing the Schwarzian effective action have since been studied [40,41]. Apart from the JT action, the Polyakov action with a cosmological constant was also studied and found to describe the AdS 2 bulk dynamics dual to the soft modes of SYK [42]. This was done by analysing the action for the co-adjoint orbits of the Virasoro group which can be thought to describe the soft modes. Different aspects of AdS 2 gravity were also covered in [38,[43][44][45][46][47][48][49][50][51]. The effect of AdS 2 arising in rotating horizons was also studied in [52], here a large N SYK like system was modelled to mimic the near horizon near extremal symmetries of Kerr-Neuman blackholes in AdS 4 . For past works the reader may refer to [53,54] and references therein.
It is also worth noting that a theory of open string governed by the Nambu-Goto action probing an AdS Schwarzchild geometry also exhibits maximal chaos [67]. In such a system the scrambling time is governed by the string tension. A Schwarzian effective action has also been uncovered for such systems as being comprised of the reparametrizations of the world sheet [68,69].
It would be an interesting question to ask if such modes can be found in thermal large N CFTs such that their effective actions govern the chaotic behaviour of the system, like in the SYK model studied in [8]. The present holographic understanding of this phenomenon allows one to visualize these modes close to extremality in the near horizon region for at least non-rotating geometries. The near extremal geometries possess a near horizon AdS 2 throat, and bulk scattering of the form studied in [3] excite graviton modes which can be described by a JT theory confined to this throat region [4]. It is worth noting that the holographic computations in [3,5] do not assume extremality. It would therefore be worthwhile to understand how these modes behave away from extremality and also in the entirety of a black hole in AdS.
To this end we address a simpler problem of that in AdS 3 which like in the dilaton-gravity theory in AdS 2 has only boundary degrees of freedom. In fact in AdS 3 these have been well studied and are called the Brown-Henneaux modes [70] which are in one-to-one correspondence with 2d infinite conformal symmetries of the boundary CFT 2 . In section 2 using the gauge gravity prescription we first formally equate the computation of eikonal scattering in the bulk done in [3] to computing correlators in the boundary CFT upto linear order in G N . This we do by computing the effective action for conformal transformations on the boundary obtained from the bulk on-shell path integral. We then (section 3) compute the effective action for the Brown-Henneaux modes about a rotating BTZ and find it to be the product of square-root Schwarzian derivatives, each for left and right moving conformal transformations of the boundary.
In section 4 we proceed to compute the correction to the 4pt function of 2 boundary operators -computed in the probe approximation; to linear order in G N . We thus reproduce the answer of [3] for the non-rotating BTZ case of λ L = 2π/β. The similar procedure when used for the rotating BTZ yields λ L = 2π/β + where β ± = β(1 ∓ µ L ) with µ L = r − /r + being the chemical potential associated with angular momentum. We thus find that for the rotating BTZ λ L = 2π/β + > 2π/β.
We end with section 5 with some conclusions and discuss the possible implications of the result. In particular we point out a possible modification of a part of the proof given in [4] so as to allow for a modified bound in the presence of a chemical potential for angular momentum.

Bulk Computation
In this section we heuristically equate the eikonal approximate calculation of [3] to the one generally done while introducing the AdS/CFT correspondence i.e. equating the bulk on-shell (small G N ) path-integral to the generating function of boundary correlators: where φ 0 is the boundary value of the scalar field in the bulk which sources a scalar operator O in the boundary CFT. Like in [3] we will consider 2 minimally coupled scalar fields in the bulk with masses m 1 & m 2 1 , with no interaction terms between them Here we have not written down the boundary terms for the actions which make the variational problem well defined and render the on-shell action finite.
We will concern ourselves with the computation of the bulk 4pt. function φ 1 φ 1 φ 2 φ 2 , which is equal to the boundary 4pt. function where each of the bulk coordinates is taken to the boundary 2 . Using the bulk path integral expression, the 4pt function is given by Here for simplicity of notation we have not mentioned the space time dependence. Since we would be concerned with the limit in which classical gravity dominates we would be interested in the saddle point evaluation of the above path-integral. Further one usually considers the probe approximation in which the scalars act as probes for a given metric which satisfies vacuum Einstein's equation with cosmological constant Λ. We will denote this solution as g µν . In this limit since the matter action is quadratic the answer is readily computed from the expression

5)
1 There is summation in i and the space time integrals are suppressed for brevity.
where l is the AdS radius.
hereφ i solves the Klein-Gordon equation in the backgroundḡ µν 3 . Let us consider the effect of first order back reaction in orders of G N by considering where g µν =ḡ µν + h µν and T µν is determined entirely in terms ofφ i s. Therefore we can rewrite (2.4) as where φ 1 φ 1 ḡ+h denotes the 2pt. function evaluated in background metricḡ µν + h µν for an h µν determined by (2.6) and constrained by boundary conditions on the bulk metric. Note, that here we have only considered diagrams in which the graviton lines attach to different scalar legs, corrections to scalar propagator due to gravitons attaching to the same scalar propagator have been ignored.
In the eikonal approximation i.e. in the limit when the scalar field momenta are taken to be large and light-like, (2.7) reduces to [71] where we have assumed that the T µν due to the matter fields are like shockwaves with lightlike momenta. In this approximation the incoming and the out-going momenta are the same and we get the effect of infinite graviton ladder exchanges between the 2 scalar propagators. This precisely the origin of the e iδ(s) factor in [5]. Here h µν is the response to the shockwave generated by high momentum scalar propagators in T µν . This approximation is justified in the shock wave analysis of [5] since the scattering is arranged to have a maximum contribution from the near horizon region. By the time any scalar perturbation reaches the horizon area it would be blue-shifted exponentially, thus giving the first leading contribution to the correction to the 4pt function. The eikonal approximation used in [3] effectively means that the contribution comes from the bifurcate horizon. One can then use bulk to boundary propagators to compute the correlation function on the boundary.
Let us now try to rephrase the same computation in AdS 3 . We write the probe approximate 4pt function as here we have haven't included S grav [ḡ] since it would a constant. The correction it would receive from the gravity path integral would be captured in (2.10) We could now consider the above path integral being dominated by the gravity saddle by considering the small G N limit. Therefore we would only be seeking contributions from geometries satisfying vacuum Einstein's equations. This way we should be able to reproduce the leading correction in the small G N limit. In the above metric path integral only boundary degrees of freedom contribute in AdS 3 . These are in one-to one correspondence with the boundary conformal transformations 4 . Further O 1 O 1 g ≈ O 1 O 1 ḡ+h would simply correspond to the change in the 2pt. functions due to conformal transformations. Therefore atleast in AdS 3 we must be able to capture the effect of chaos via (2.10).

AdS story from the boundary
In this section we will derive the effective action for the soft modes by calculating the gravity bulk on-shell action about an arbitrary BTZ geometry. As justified in the previous section this amounts to finding the effective action for the diffeomorphisms in the bulk respecting Dirichlet boundary conditions. It is known that such configurations in the bulk are dual to states in the boundary CFT which are created by the action of 2 copies of commuting Virasoro algebra. We begin by writing the most general bulk configuration with the boundary metric being flat, we begin with a Lorenztian metric in the Fefferman-Graham gauge [72,73].
in the Fefferman-Graham gauge with T ±± = (r + ± r − ) 2 . It is worthwhile to notice that the radial coordinate in (3.1) sees the horizon at r h = r 2 The bulk action with relevant boundary(counter) terms is [73,75] 16πG where Λ = −1/l 2 and l is the length of AdS 3 . The 1/l term in the boundary action is used to make the on-shell action finite. Therefore the on-shell value of (3.3) for arbitrary metrics given by (3.1) is Since gravity in 3-dim is non-dynamical in the bulk, all solutions to the bulk action (3.3) can be obtained from one another via diffeomorphisms. In fact the Fefferman-Graham theorem [73] allows one to express any solution in the form (3.1) as a diffeomorphisms of the other and we thus only need to concern ourselves with diffeomorphisms that preserve this form to generate all solutions. We would therefore like to know the action (3.4) associated with such a diffeomorphism, having started from a particular solution (for eg.) of the form of the BTZ metric (3.2). To this end we would like to know how T ±± depend on these diffeomorphisms.
We notice that for infinitesimal diffeomorphisms which maintain the form of (3.1), the change in T ±± is given by [70,76] δT We also note that change in a Schwarzian derivative {T (u), u} due to a diffeo u → u + ǫ(u) is: One can find the full non-linear completion of (3.6) [77] which takes a Poincáre AdS 3 (with T ±± = 0) to (3.1). Under such a diffeomorphism the stress tensor is proportional to the Schwarzian for boundary conformal transformations. Comparing (3.5) and (3.7) we can deduce that under where for infinitesimal diffeomorphisms X ± ≡ x ± + ξ ± (0) . One notices that for X ± = x ± =⇒ T ±± = 0. This value can be shifted by defining where the choice of X ± ′2 makes sure that the linear in T ±± terms in (3.5) remain the same.
Here, L ± define the charge of the BTZ metric about which the change in the parameters T ±± is measured 7 .
Therefore the on-shell action in (3.4) is where L ± decides which bulk configuration one mesures the change from. The above action is defined on the boundary of AdS 3 , the integral in the AdS radial direction receives contribution from r = r h and the boundary r = ∞. The divergent contribution from the boundary at r = ∞ is cancelled by the holographic boundary counter-terms, the only finite contribution comes from the horizon. Under infinitesimal diffeomorphisms X ± → x ± + ǫ ± (x ± ) the quadratic action takes the form where we have ignored boundary terms. Since we would be interested in computing OTOC's later we Euclideanize the above action The quadratic action divides itself into left and right sector. The action (3.10) evidently has the symmetries of the Schwarzian derivative, the infinitesimal versions of which are manifested in (3.12). We would correspondingly have got the above action by working in the bulk in a Euclidean setting to begin with. This would have invariably required us to have the angular momentum associated with the BTZ metric to be imaginary so as to have a real Euclidean metric.
One can in principle derive the one-loop exact action for the boundary gravitons in AdS 3 as was done recently by [78] using the Chern-Simmons prescription. The authors obtain a theory of re-parametrizations which encodes loop contributions of boundary gravitons in a perturbation in 1/c 8 .

Propagators
We are now in the Euclidean setting where the time τ = −it is along the imaginary direction while the space-like coordinate φ is real. The x ± coordinates of Euclideanised BTZ metric can be regarded to have complex periodicities [79] x We regard the integral in (3.12) to be in one such periodic interval, therefore the above action splits into two 1-dim actions where α ± = −2πl . For convenience we have defined and we analyse the propagator for ǫ + , for which we evaluate the Green's function G + for the operator Here, we observe that G + would depend onz by a function of the ratioz/β + . The zero β + }. We will compute the G + first for the Schwarzchild case and try and generalize for the rotating BTZ case.
For AdS 3 Schwarzchild, β + = β − = β ∈ R. Assuming G + be a function ofz/β we solve (3.16) for real values ofz/β i.e. G + (τ /β) and then using Schwarz's theorem analytically continue it for arbitrary complex values ofz/β [80]. This allows us to express both sides of the (3.16) as a discrete sum, thus where the prime on the sum denotes n / ∈ {−1, 0, 1}. Doing the relevant Matsubara summation and analytically continuing in the complexz/β plane yields Here, The last 2 terms 9 in (3.18) are comprised of the zero modes we neglected in the sum and would drop out of any computation which respects the bulk isometries. An identical expression would exist for G − in terms of z/β. 9 Explicitly: a = 1 + π 2 6 + 7π 4

Rotating BTZ propagators
One could extend the above method naively to rotating BTZ case. This would require solving the (3.16) for the real values ofz/β + where Ω = µβ. Same holds true for z/β − ∈ R. It is quite clear from the outset that one could define Euclidean coordinates {τ ,φ}. (3.20) Which is a conformal transformation on the boundary metric: We do not however perform such a transformation on the propagator, we merely use (3.20) for making the coordinated dependence look simple. Therefore solving (3.16) for real values ofz/β and analytically continuing we get Note that the conformal transformation (3.20) isn't one of the SL(2, R) zero modes.

4pt correlator
In this section we will use the propagators obtained in the last section to compute the next to leading order in G N corrections to the 4pt. function. We consider the first the leading contribution to the Euclidean 4pt. function of four scalars [79] where z 12 = z 1 − z 2 and V 1 ≡ V (z 1 , z 1 ). We would be interested in seeing how they would depend on the bulk on-shell metrics in the path integral (2.10). As explained before, these would correspond to computing the change in (4.1) due to conformal transformations parametrized by ǫ It can be seen that B above is invariant under the SL(2, R) zero modes of The correction to (4.1) is obtained by Wick contracting the ǫ ± s with each other using the propagator (3.18) and it's complex conjugate. We will first analyse this around AdS 3 Schwarzchild and then in a generic rotating BTZ background.

AdS 3 Schwarzchild
For the case of AdS 3 Schwarzchild we have β ± = β. The reality condition of (3.19) implies φ = 0, i.e. we compute the propagators G ± along the τ real line and then analytically continue. Here to proceed we first have to order the Euclidean times for the operators in question and then use the appropriate propagator value for Wick contraction. We then add arbitrary Lorenztian time arguments corresponding to each operator and then read off the answer. The Euclidean answer to the expression looks like  Similarly we introduce Lorentzian time coordinates for each of the operators i.e. z → τ + iφ − it, we fix τ 1 = β, τ 3 = β/4, τ 2 = β/2, τ 4 = 3β/4 and further fix t Here one clearly sees the exponential behaviour of the correlator for both the boundary null coordinates t ± φ as β (4.6) Thus we see that the Schwarzian action (3.10) associated with Brown-Henneaux modes (3.6) are responsible for the maximal Lypunov index at least in non-rotating BTZ. Note that since we have been cavalier about causality in our propagators (3.22) we would not reproduce the (t − |φ|) behaviour in the exponent like [3] 10 .

Rotating BTZ
Let's do the similar exercise for rotating BTZ where β ± are complex parameters. Here in order to use the propagator in (3.18) we would have to use a shifted Euclidean time in (3.20) τ = (τ + µφ)/(1 + µ 2 ) for ordering the different Euclidean times. i.e. for time ordered correlator we arrangeτ 1 >τ 2 >τ 3 >τ 4 while for out-of-time-orderedτ 1 >τ 3 >τ 2 >τ 4 . We would then fix theτ on a circle of period β. We will return later to the detail of how specifying the Euclidean time ordering in shifted Euclidean time does not effect the end result.
Having done this (6.1) shows a polynomial growth in time.
where we define the Lorenztian angular velocity as µ L = iµ = r − /r + as it would has risen in a Lorenztian bulk geometry, thus yielding 11 A some what similar conclusion was reached in [81] with mutual information computed between the left and right intervals of | TFD corresponding to a rotating BTZ. This was computed both by computing mutual information via Rényi entropy in the 2D CFT with a chemical potential µ perturbed by a heavy operator, and from the bulk by employing the Ryu-Takayanagi prescription of minimal surfaces in a shock wave background. In [81] the symmetry between β ± is broken by the spatial arrangement of the heavy operator relative to the entangling interval in question. This arrangement is such that only one of the modes with a smaller temperature effects the entangling region for positive times. For a different spatial arrangement one may seem to find that the scrambling time computed is governed by the higher of the 2 temperatures.
In the above computations we had chosen to order the Euclidean times in the shifted coordinate frame ofτ = (τ + µφ)/(1 + µ 2 ) instead of ordering their τ coordinates. The reason why we still end up computing the desired correlator in the desired frame is as follows: The behaviour of the 4pt function is encapsulated in the sum conformal blocks. This in the limit of large central charge for holographic CFTs can be seen to be captured by only the Virasoro conformal block of the identity operator [82]. The conformal blocks are given by hypergeometric functions F (a, b, c, η) 12 that are known to have a branch point at η = 1 where η is the conformal cross ratio 13 . It was shown in [82] that the behaviour of η for large Lorentzian times depends on the ordering of the Euclidean times 14 . For out-of-time-ordering of the Euclidean times, η circles the branch point η = 1 in the complex plane; while this is not the case when the Euclidean times are time-ordered. The authors of [82] then use the monodromy of F (a, b, c, η) to compute it's behaviour for small η thus obtaining the Lyapunov exponent for 2dim holographic CFTs. One can therefore use this behaviour of the cross ratio to discern if one is indeed computing time-ordered or out-of-time-ordered correlators.
We define the cross ratio η as and similarly for the barred coordinates which are complex conjugates. η can be analytically continued to include Lorentzian times by τ → τ − it. The tilde coordinates are related to un-tilde coordinates by (3.20) i.e.
In all the above cases we chooseφ 1,2 = 0,t 3,4 = 0 &t 1,2 =t,φ 3,4 =φ, therefore for fixed values of Euclidean timesτ i , η(t,φ) is a complex function. It was demonstrated in [82] that for |t| > |φ| > 0, the out-of-time-ordered cross-ratio η oto (t,φ) circles the branch point η = 1 as one variest; while no such behaviour is found in the case for the time-ordered cross-ratio η to (t,φ). In order to know how the shifted Euclidean time ordering effects the answer in the un-tilde frame we would like to see how η (andη) behaves as a function of {t, φ} for the respective orderings inτ i .
The cross-ratio η for arbitraryτ i looks like where in the second line above we have performed the Lorentzian coordinate transformation (4.8). A similar expression holds forη with β − → β + and t + φ → t − φ. Here β± = β(1 ∓ µ L ). We plot the second expression in (4.12) for different values of µ L ∈ {0, 1}, τ i = β − (i − 1) 3β 4 + (i − 1)ǫ β 4 and β = 2π with ǫ ∈ {0, 1} for the time-ordered case. For the out-of-time-ordered case we exchangeτ 2 ↔τ 3 as compared to the time-ordered case. We find the same behaviour as seen in [82] for different values of chemical potential µ L and ǫ. Here ǫ varies the adjacent Euclidean time differences on the thermal circle linearly from β/4 to zero. Therefore even though we ordered shifted Euclidean times we end up computing the desired correlator in the un-tilde Lorentzian frame spanned by {t, φ} 15 .
One could further ask the question whether a purely CFT computation of the form of [82] can recover(or corroborate) the results obtained for the rotating BTZ. This would require using the hypergeometric function F (a, b, c, η) in the large central charge limit to incorporate not only temperature but also µ L due to rotation. In [83] it was shown that effect of introducing heavy operators at time-like infinities in 2dim CFT corresponds to evaluating the global conformal blocks on a new 2dim background metric obtained by mapping the flat Euclidean space to a torus by 16 ω = e −2πiz/β − ,ω = e 2πiz/β + . where {ω,ω} are the Euclidean flat-space light-cone coordinates. We leave this as an exercise for the near future.

Results and Discussions
The bulk understanding of how the Lyapunov index is 2π/β has as of yet always relied upon near the extremal property of an AdS black hole exhibiting an AdS 2 throat. In some sense the back reaction of the scalars in the bulk gets the most contribution from this region. Any deviation in the AdS 2 geometry is captured in an action like the Jackiw-Teitielboim thus giving rise to a Schwarzian action. What we have demonstrated here -at least in AdS 3 ; is that the Schwarzian arises even when one is far away from extremality. Moreover we find this as an effective action at the boundary of AdS 3 rather than at some screen in the interior. It would be interesting to investigate how such an action can be arrived at for black holes in AdS d>3 , this would indeed give some understanding of the soft modes in higher dimensional large N CFTs.
In the probe approximation there is an inherent conformal symmetry in the bulk emanating from the asymptotic symmetries of AdS 3 . This would correspond to the 2 copies of Virasoro algebra in the boundary CFT. Any arbitrary solution to the Einstein's equation with matter would not have such a symmetry. Expanding perturbatively about the probe approximation in orders of G N (i.e. back reaction) breaks this symmetry spontaneously. The action (3.10) therefore can be seen as the action cost associated with conformal transformations at the boundary of AdS 3 when one tries to go away from the probe approximation to linear order in G N .
Extremality can be reached in the simplest possible manner by turning on charges for the AdS black hole, the top down understanding of [37] in such a setting was explained in [39] in AdS 4 . Here the authors studied a probe uncharged scalar in the bulk thus having no dynamics for the gauge field. It would be interesting to analyse how the near horizon picture in [37] is reached for rotating geometries close to extremality in AdS d>3 . In [84] 5d rotating Kerr geometries were analysed close to their near extremal limit in the near horizon throat region. There the authors have discovered a generalized JT action consisting of a dilaton and an additional scalar.
In (4.7) we take the view that β ± are complex to begin with. The complex value of β ± = β∓iΩ is required to make sense of the Euclidean BTZ metric as a real quantity. Further the iǫ 17 prescription that we use requires us to first compute a Euclidean correlator and then analytically continue it to desired Lorentzian times. This is similar to the technique employed in [81] for computing the mutual entanglement from the 2D CFT, as computing the Rényi entropy involves analytically continuing the Euclidean 4pt correlator ψσσψ † 18 to obtain the Lorentzian answer. This requires making the left and right moving temperatures real: β ± = β ∓ Ω 19 when all Euclidean times have been put to zero. (4.7) would therefore yield a growth in the scrambling time that would be greater than the Lyapunov index of 2π/β due to the left moving (anti-holomorphic) modes for x + i.e. 2π β−Ω .
The presence of rotation in the bulk implies a CFT with a chemical potential corresponding to angular momentum. 1d SYK and gauged-SYK models have been studied in the presence of a chemical potential [21]. Here the Lypunov index computed was found to be bounded by 2π/β. This also bodes well with the intuition that holding other charges fixed makes the system less chaotic. However the chemical potential present in such cases were associated to an internal symmetry and not a space-time symmetry.
The analysis of section (4.2) for the case of rotating BTZ seems to yield a result in contradiction with the mathematical proof in [4]. The proof in section (4.1) of [4] is basically based on the maximal modulus theorem for a bounded holomorphic function. Here we try to give a simple understanding as to how one may try to reconcile the result of section (4.2) of this paper with that of [4]. 17 Here Euclidean time is τ instead of ǫ. 18 ψ is the heavy operator generating the shock-wave in the BTZ while σ is the twist operator. 19 In our analysis β ± ∼ 1 √ are associated with x ± respectively.
For the case of 2d CFTs one could argue a change in the Lyapunov index using the fact that the function of the form of V 1 V 2 W 3 W 4 V 1 V 2 W 3 W 4 of Virasoro primaries must be a scalar under conformal transformations.
For the specific case of 2dim CFT dual to the Schwarzchild-AdS 3 this function was shown [3] to have a behaviour of The finite diffeomorphism which takes Schwarzchild-AdS 3 to a generic rotating BTZ blackhole corresponds to a conformal transformation on the boundary of the form (4.8) where β ± = (1 ∓ µ L ) and δx is the difference in any two of the coordinates in the correlator. This would imply studying a holographic CFT with a chemical potential corresponding to rotation turned on. The proof in section (4.1) of [4] relies crucially upon mapping the half strip of width β/2 in Euclidean time τ to a disk via a conformal map .

(5.3)
This map has a periodicity under τ → τ + β which is exhibited by It was argued that this half strip is the minimum area of analyticity a function of the form Smaller the width of the strip, greater would be the bound on chaos. However conformal transformation of the form (5.2) suggests that there exist 2dim holographic CFTs for has a thinner region of analyticity of width β − /2 = β(1 − µ L )/2 controlled by the chemical potential µ L . This suggests that repeating the analysis of [4] for systems at finite temperature and chemical potential might in certain cases yield a smaller region of analyticity in the Euclidean time direction.
One could very well have guessed such a result simply by observing that the effective action (3.10) doesn't mix the left and right movers, each of which have different inverse temperatures i.e. β. The conformal blocks used in [82] to study chaos in 2dim holographic CFTs also exhibit such a holomorphic and anti-holomorphic factorization. From such a point of view it seems reasonable to expect that there would be 2 Lyapunov indices governed by the 2 temperatures; this point of view does not seem to be at odds with the statement of maximal chaos first proposed in [4]. Therefore the maximal growth would be governed by the smaller of the two i.e. min[β + , β−], while the surface gravity of the bulk would be related to the average β + +β − 2 . However one can easily see that this would also lead to extremal BTZ metric exhibiting chaotic behaviour with the Lyapunov index being one of the 2 temperatures which is non-vanishing. This would precisely be the temperature which would be greater than that given by the surface gravity for a generic rotating BTZ. This raises an important question in higher dimensional Kerr-AdS black holes of the possibility of extremal chaos. It would be very interesting to see how these considerations would have to be modified when analysing rotating geometries in AdS d>3 as unlike AdS 3 the bulk degrees of freedom of the metric would also participate in the dynamics. It is also worth noting that bulk explanations of the dynamics of chaotic behaviour using the 2dim Jackiw-Teitelboim model in the near horizon region of near extremal black holes [37] cannot be used to explain such extremal chaos. Therefore it would be interesting to find a similar JT like model(if possible) in the near horizon region of extremal and/or near extremal region to explain such a behaviour. This would provide a useful tool in analysing such behaviour in higher dimensional Kerr-AdS black holes.
The result of section (4.2) is also validated by the analysis of mutual information for late times computed in the BTZ geometry subjected to a shock wave [81]. Here the author found the Lyapunov index to be related to the smaller of the two temperatures i.e. λ L = 2π β − in the conventions of this paper. The mutual information in the |T F D state corresponding to an eternal BTZ subjected to a shock wave is computed in [81] on the boundary by taking the limit of the Rényi entropy, and in the bulk by employing the Ryu-Takayanagi prescription of minimal area. However the spatial arrangement of the heavy operator in [81] w.r.t. the entangling region under consideration only sees the effect of one of the modes.
We have analysed the bulk on-shell action (3.10) upto quadratic order in fluctuations about a generic BTZ metric. It is in this approximation that the left and the right sectors seem to decouple and the scrambling times for the left and the right sectors are governed by their respective temperatures. The full non-linear action (3.10) if expanded beyond quadratic order in fluctuations would indeed have terms which mix the left and the right moving sectors. It would be interesting to study how the three temperatures β & β ± participate in the process of reaching a thermal equilibrium at late times when the boundary CFT is excited. This would require a more detailed analysis for longer time periods than the scrambling times of the left and right sectors 20 .
The techniques used in this paper seem to be too well suited for AdS 3 . As mentioned before that generalizing this to higher dimensional AdS black holes would be interesting, it would nonetheless be easier to analyse the rotating BTZ along the lines of [3] by computing bulk eikonal scattering which seems to be an analysis suited for all dimensions 21 .
To conclude, this work also suggests that if the Lypunov index associated with rotating AdS black holes in Einstein-Hilbert theory have maximal chaos, then for large N thermal CFTs with chemical potential associated with angular momentum the chaos bound is λ L = 2π β(1−µ L ) . It would be interesting to find a more thorough generalization of the proof in [4] for large-N CFTs with chemical potential associated with angular momentum.

Acknowledgements
The author is indebted to Gautam Mandal for fruitful discussions throughout the duration of this work. The author also benefited by earlier work with Gautam Mandal, Pranjal Nayak, Nemani Suryanarayana and Spenta Wadia. The author also acknowledges Adwait Gaikwad, Anurag Kaushal, R. Loganayagam, Pranjal Nayak, Ronak Soni, Shiraz Minwalla & M. V. Vishal for their numerous discussions during the completion of this work. The author also acknowledges Abhijit Gadde and Justin David for giving inputs on this work. The author also found the academic environment of the Indian Strings Meet 2018 held in IISER Thiruvanathapuram conducive towards the completion of this work.