arxiv: v1 [hep-th] 28 Aug 2016

Transcription

1 Quantum cat map dynamics on AdS 2 arxiv: v1 [hep-th] 28 Aug 2016 Minos Axenides 1, Emmanuel Floratos 1,2 and Stam Nicolis 3 1 NCSR Demokritos, Institute of Nuclear and Particle Physics Aghia Paraskevi, Attiki, Greece 2 Physics Department, University of Athens, Zografou University Campus Athens, Greece 3 CNRS Laboratoire de Mathématiques et Physique Théorique (UMR 7350) Fédération de Recherche Denis Poisson (FR 2964) Département de Physique, Université François Rabelais de Tours Parc Grandmont, Tours, France axenides@inp.demokritos.gr, mflorato@phys.uoa.gr, stam.nicolis@lmpt.univ-tours.fr Abstract We present a toy model for the chaotic unitary scattering of single particle wave packets on the radial AdS 2 geometry of extremal BH horizons. Based on our recent work for the discretization of the AdS 2 space-time, which describes a finite and random geometry, by modular arithmetic, we investigate the validity of the eigenstate thermalization hypothesis (ETH), as well as that of the fast scrambling time bound conjecture (STB), for an observer with time evolution operator the quantum Arnol d cat map (QACM). We find that the QACM, while possessing a linear spectrum, has eigenstates, which can be expressed in closed form, are found to be random and to satisfy the assumptions of the ETH.The implications are that the dynamics is described by a chaotic, unitary, single particle S-matrix, which completely delocalizes and randomizes initial gaussian wave packets. Applying results obtained by Dyson and Falk for the periods of the Arnol d Cat Map(ACM),which are related to its mixing time, we also find that the thermalization of infalling wave packets in this particular model is exponentially fast saturating the STB, under the constraint that the finite dimension of the single particle Hilbert space takes values in the set of Fibonacci integers.

2 Contents 1 Introduction 1 2 Modular discretization: from the torus to AdS The Arnol d cat map and Fibonacci chaos on AdS 2 [N] 6 4 Chaotic eigenstates of the QACM on AdS ETH and the scrambling time bound for the QACM 13 6 Summary and conclusions 18 A Coset geometry and light cones on AdS 2 19 B The Heisenberg Weyl group and the Weil representation of PSL 2 [p] 22 C The construction of the QACM eigenstates 26 1 Introduction A very interesting revival of the old relation between the near horizon shock wave BH geometries with gravitational memory effects and the information paradox has recently appeared [1, 2] and references therein). It seems possible in principle, that the horizon region of BH could form a random basis of purely geometrical data of all of its past and recent history, through the t Hooft mechanism of permanent space-time displacements caused by high energy scattering events of infalling wave packets [1]. In the language of refs. [2], such data can be identified with the soft hair of the BH, whose origin is the infinite number of conservation laws, described by the BMS group. This new reincarnation of the t Hooft-Susskind horizon holography [3] gives hope, that the emitted Hawing radiation will be influenced by the shock-wave spacetime geometry structure and it will encode, in its angular and time correlations, past and recent BH memories. Although a realistic calculation with a truly chaotic horizon region, in general, still isn t possible, some progress could be made using simple mathematical toy models, which can describe how the information carried by infalling wave packets, could be encoded in the Hawking radiation, through a chaotic single particle S-matrix [4, 5, 6, 7, 8, 9]. Proposals for a chaotic, discretized, dynamics for the microscopic degrees of freedom of the stretched horizon have been discussed for quite some time in the literature [10, 11, 12, 13, 14, 15, 16]. 1

3 Recently, we constructed a consistent, modular, discretization, AdS 2 [p], of the AdS 2 near horizon geometry of nearly extremal black holes, for every prime integer p [17]. AdS 2 [p] has a random structure due to the modular arithmetic, which can capture relevant properties of the shock wave geometry, that includes permanent random spacetime displacements. These, supposedly, encode memories of the past scattering events that formed the BH horizon. As explained in ref. [17], by explicit calculation, this discretization has the merit of realizing the holographic correspondence between the bulk, AdS 2 [p] and its boundary RP 1 [p], the discrete projective line. The reason this discrete holography exists at all is that it is possible to realize the action of the finite symmetry group PSL 2 [p] as an isometry group of the bulk and, also, as the (Möbius) conformal group on the boundary. This discrete geometry, also, provides a natural framework for describing the discrete single particle dynamics, via observers, with time evolution operators that are elements of the isometry group of AdS 2 [p], consistent with known results of the superconformal mechanics of extremal black holes [18]. In the present work we shall take the classical evolution map to be the well known Arnol d cat map (ACM) [19], for an observer on a nice slice [17].This map is appropriate for describing the single particle dynamics, since it belongs to the discrete isometry group of AdS 2 [p], which is PSL 2 [p]. The ACM can be quantized, i.e. it is possible to construct explicitly a p p unitary evolution operator, called the quantum Arnol d cat map (QACM). This definition uses the Weil representations of SL 2 [p] and especially those that correspond to its projective action, by PSL 2 [p] on AdS 2 [p]. This construction extends the results for the case of the discrete torus [20, 21, 22, 23, 24]. We find a linear spectrum and the exact analytic eigenstates of the QACM, which are random in a specific way, that is, the squares of the absolute values of the amplitudes (probabilities) are drawn from a (discrete) Gaussian distribution, while their phases have a flat distribution. These properties provide support, for the Eigenstate Thermalizaton Hypothesis [25], as applied to the closed quantum system of observer+bh horizon. The implications are that the dynamics of infalling localized wave packets is described by a chaotic unitary S-matrix, which completely delocalizes and randomizes initial Gaussian wave packets. Moreover, by identifying the scrambling time of the dynamical system with its mixing time [19], we are able to evaluate it analytically, using results of Dyson and Falk for the periods of the QACM [26]. For chaotic dynamical systems the mixing time has a very precise meaning as the time required for an initial distribution function to become flat (i.e. to thermalize). In the case of the QACM the periods, also, provide the degeneracies of the energy levels and thereby imply the existence of non trivial conservation laws, which are determined explicitly. These degeneracies present an obstacle for realizing quantum ergodicity for all large values of p. This problem is known as the problem of quantum unique ergodicity 2

4 and has been discussed in the quantum chaos literature [27]. Finally we find that the scrambling time bound [28] is saturated for the QACM under the condition that the dimension of the single particle Hilbert space of states takes values in the sequence of Fibonacci integers. An introduction to the requisite tools from arithmetic geometry and computational number theory can be found in ref. [29]. Next we proceed with the plan of the paper. In section 2 We present the discretization of AdS 2 through the hyperbolic lattice of its integral points and the corresponding lattice of causal diamonds.we find the natural action on this geometry of the discrete isometry group PSL 2 [Z] and of its mod p reduction, PSL 2 [p]. This discretization procedure is in strict analogy to that of the torus geometry but of course the action of the isometry group is projective and not linear as in the Torus case. In section 3 we recall the properties of the ACM, its relation to the Fibonacci sequence and its periods mod p. We study its group of symmetries, inside PSL 2 [Z], i.e. the set of elements of PSL 2 [Z], that commute with it. These properties are used to describe the dynamics of a nice slice observer with time evolution defined by the succesive applications of the ACM for motions on the AdS 2 [p] space-time. In section. 4 we construct explicitly the exact Quantum Arnol d cat map, using the metaplectic (or Weil) representation of SL 2 [p] which is reducible and splits into (p + 1)/2 and (p 1)/2 dimensional irreducible ones.one of these two unitary irreps, depending on the form of the prime number p, is also a representation of PSL 2 [p]. We determine analytically the spectrum and the eigenstates of the QACM in the appropriate representation of PSL 2 [p] and we compute their degeneracies. We obtain the interesting result that, while the spectrum is linear, the eigenstates are chaotic in the specific way, that is, the squares of the absolute values of the amplitudes (probabilities) are drawn from a (discrete) Gaussian distribution, while their phases have a flat distribution. These results support the validity of the eigenstate thermalization hypothesis (ETH) for quantum ergodicity or unitary thermalization. We discuss also relevant problems for the unique ergodicity hypothesis pertaining to this map. In section 5 we use the results obtained in section 4 to study the spectrum of scrambling times for various values of p and we find that for Fibonacci integer values the scrambling time bound of Susskind and Sekino is saturated. Finally in section 6 we discuss our results and open problems for future work. In the appendix A we recall our construction for the continuous geometry of AdS 2 as a ruling surface which facilitates the visualization of the modular discretization as well as the coset nature of the geometry. In appendix B we collect all the necessary material for the detailed construction of the Weil representation of SL 2 (p) and PSL 2 (p) and in appendix C we present the technical details for the analytic construction of the eigenstates and eigenvalues of the QACM. 3

5 2 Modular discretization: from the torus to AdS 2 In this section we shall review the modular discretization of AdS 2, introduced in [17]. The coset geometry and the light cone structure of AdS 2 as well as its global cover(the Einstein strip) for the standard continuous case is reviewed in appendix A. We define the modular discretization by replacing the set of real numbers, R, by the set of integers modulo N.The so obtained coset finite geometry AdS 2 [N] = SL 2 (Z N )/SO(1, 1, Z N ) is a finite random set of points in the embedding Minkowski spacetime due to the modular arithmetic. This type of geometries are studied in mathematics under the name of arithmetic geometries [29]. This discretization is a possible toy model for describing nonlocality, chaos and quantum information processing in the vicinity of the BH horizon. The discretization parameter N acts as an infrared cutoff and 1/N as an ultrarviolet cutoff for the bulk.the basic scale here is the AdS 2 radius which we set equal to one. We have shown that this model of the bulk spacetime realizes the holographic correspondence with the boundary and the UV/IR connection holds [17]. The infrared cutoff is needed in order to get a compact phase space for the dynamical system. This is an important requirement for chaotic systems to be mixing [19]. In the next section we shall introduce an extension of the Arnol d cat map (which is usually considered to act on a two dimensional toroidal phase space) for the dynamics in the radial and time AdS 2 geometry of the nearly extremal BH. We proceed in analogy with the dynamics of the ACM on the torus. For the study of the chaotic and mixing properties of the ACM on the torus, we first compactify the standard planar phase space R 2 by passing to the toroidal phase space T 2 = R 2 /Z 2. Correspondingly for the AdS 2 geometry, we introduce the hyperbolic lattice of its integral points, which can be generated by the adjoint action of the modular group SL 2 (Z) on an appropriate set of initial points. The vertices of the lattice in AdS 2 are all the integer solutions of the equation x x 2 1 x 2 2 = 1. It is known that the number of such points grows as the area of AdS 2, when an infrared cutoff, L, on x 2 ( x 2 L), goes to infinity. This is property is important for the implementation of holography since it implies that most of the area of AdS 2 is concentrated on its boundary. For the ACM it is known that all of the toroidal phase space is a chaotic attractor. The unstable periodic orbits leading to this attractor have rational coordinates with the same denominator N for each orbit. To study these periodic orbits one has to restrict the action of ACM on a Z N Z N discrete torus.the corresponding period T [N] of ACM is a random function of N and its values are the periods mod N of the Fibonacci sequence [26], as we shall describe in the next section. 4

6 Although for the recovery of truly chaotic orbits one has to send N to infinity,with a finite resolution bigger than 1/N, the ACM produces chaos already at moderate values of N. Our proposal for the chaotic dynamics on AdS 2 [N] essentially boils down to the requirement of studying periodic orbits, leading to chaotic mixing, by passing to the mod N projection of SL(2, Z) i.e the group SL(2, Z N ). In other words we pass to the quotient SL 2 (Z)/Γ N, where Γ N is the principal modular subgroup of the modular group. All this can be simply realized by selecting all the integral points of AdS 2 and then taking their mod N projection. The set of points of the finite geometry of AdS 2 [N] is, by definition, the set of integer solutions mod N of the equation x x 2 1 x mod N (2.1) This can be parametrized in analogy with the ruling coordinates in the continuum case (Appendix A) as follows: x 0 (a b µ) mod N x 1 (b + a µ) mod N x 2 µ mod N (2.2) where a 2 + b 2 1 mod N and a, b, µ {0, 1, 2,..., N 1}. The random set of points in the embedding Minkowski spacetime M 2,1 thus produced are to be considered as space-time defects and they may be useful for describing the permanent infalling shock wave effects on the BH horizon s past history through the t Hooft scattering mechanism. Thus, the discretized, spatial part consists of N points and the Hilbert space of single particle states has dimension N. The microscopic degrees of freedom of the near horizon discrete geometry, have a Hilbert space of dimension 2 N (e.g. in the case of spin 1/2 states in each spatial point). From this we conclude that the entropy, S BH, of such configurations, is proportional to log 2 N = N. To derive the Horizon BH entropy from its radial, near horizon, geometry is quite subtle and requires a special treatment of an IR cutoff [30]. which, in the case at hand, is N. The discussion of the group theoretic properties of this discrete geometry is facilitated, if we restrict N to be a prime integer, p. The extension for arbitrary, odd, integer values of N is easily realized by using appropriate factorization theorems [24]. The finite geometry, AdS 2 [p], has as isometry group the finite modular group, PSL 2 [p]. This group is obtained as the reduction mod p, of all elements of P SL(2, Z). The kernel of this homomorphism is the principal congruent subgroup, Γ p. The order of PSL 2 [p] is p(p 2 1)/2 and the order of its dilatation subgroup is (p 1)/2, thus, the number of points of AdS 2 [p] is p(p + 1). 5

7 Numerical experiments suggest the following recursion relation for the number of points of AdS 2 [p k ], Sol(p k ), Sol(p k ) = p 2(k 1) Sol(p) Sol(p k ) = p 2k 1 (p + 1) (2.3) where Sol(p) = p(p + 1) and k = 1, 2,... for any prime integer p. Therefore, by the prime factorization theorem, we can determine the number of points for any, odd, integer N. For N = 2 n we find Sol(2) = 4, Sol(4) = 24, and Sol(2 k ) = 4Sol(2 k 1 ), for k 3. We remark that N = 4 is an exception. The solution is Sol(2 k ) = 2 2k+1, for k 3. Therefore we may deduce the expression for the number of points, for any integer N by prime factorization. From these results we deduce that, for large N, the number of solutions, mod N, scales like the area, i.e. N 2. So most of the points of AdS 2 [N] are close to its boundary and holography is possible in this case too. We close this section by some observations about the construction of the AdS 2 integral lattice using a specific base of generators of the modular group SL(2, Z). This specific base has the interesting property to organize the integral points into a lattice of causal diamonds. The base generators, R and L, of SL(2, Z), satisfy the braid relation and can be chosen to be L = ( ) RLR = LRL (2.4) and R = ( ) (2.5) Starting from the points (1, 0, 0) and (0, 1, 0), the corresponding SL(2, Z) orbits define vertices of light cone diamonds of the continuous AdS 2. These properties hold obviously also for the modn projection SL 2 [N]. The details will be presented in forthcoming work. 3 The Arnol d cat map and Fibonacci chaos on AdS 2 [N] The motion of p branes, near the horizon of large extremal supersymmetric BH, at the probe approximation, has been studied several years ago and it has been shown that the induced worldvolume dynamics is governed by superconformal Hamiltonians and the corresponding superconformal quantum mechanics [18]. This means that the Hamiltonian is an element of the isometry group of the AdS 2 space time P SL(2, R). In analogy, for our discrete finite geometry AdS 2 [N] we may choose as the one time step evolution map an element A of its isometry group, SL(2, Z N ). 6

8 The discrete time evolution of the motion of a particle on AdS 2, for an initial point, X 0 AdS 2, is given by the Weyl transformation (cf. Appendix A) X n+1 = AX n A 1 = A n X 0 A n (3.1) where n = 0, 1, 2,... labels the (stroboscopic) time of the ACM observer. We point out that this action differs from the action of ACM on the torus, which is linear in A [19]. Indeed, the action of the matrices A and A produce the same orbit. The action of SL(2, Z N )) on AdS 2 thus is identified with that of its projective counterpart P SL(2, Z N ). This will play an important role in choosing the Hilbert space of single particle states, to carry the unitary irreps of this group. We make a specific choice, introduced in ref. [17] of an observer, whose motion is along a nice slice, described by the Arnol d cat map(acm): ( ) 1 1 A = (3.2) 1 2 The map corresponds to successive kicks, forwards and backwards along the light cone of AdS 2, since ( ) ( ) A = = LR 1 (3.3) This is an element of P SL(2, Z N ). We choose to use this particular map, for the following reasons: The ACM has been thoroughly studied for its area preserving action on the classical toroidal phase space and it is known to possess ergodicity,exponentially fast mixing and an infinite number of unstable periodic orbits.the important property of mixing which is the technical definition of scrambling assumes that the phase space is compact. The interesting feature of these toroidal periodic orbits is that they are all known. Indeed all the rational points of the torus, with common denominator N, are mapped onto themselves under the action of the ACM, which then has a definite period, T [N], depending randomly on N [26]. These are the only unstable periodic orbits of ACM and, as N goes to infinity, we cover all the rational points of the torus and finally we reach the continuum and the, truly, chaotic orbits. The ACM has been studied intensively also as a toy model for semiclassical quantum chaos on the toroidal phase space [20],although the degeneracies in its spectrum impose additional constraints on its quantum ergodic properties [27]. 7

9 Here we shall study the classical and quantum motion of particles under the Arnol d cat map on the discretized AdS 2 [N] geometry. As we explained in the previous section this is the discretized deformation of the radial and time geometry AdS 2 of the near horizon region of extremal BH s. This motion is the longitutinal motion of probes and it differs from the motion along the horizon, which is the two dimensional sphere. On the other hand, physically, the scrambling of information on the horizon happens at the same time as the longitudinal (radial) scrambling [31]. We are interested in finding the shortest possible periods mod N, which will determine the fastest scrambling time of radially infalling pure states on the BH horizon, under the action of the quantum cat map (QACM). Moreover, we would like to understand how the quantum chaotic properties of this map are realized. For that we have first to understand the classical ACM dynamics on AdS 2 [N]. An important property of ACM is that it is known to generate the sequence of Fibonacci numbers, defined by ( ) ) ( ) fn 1 = (3.4) We observe that therefore fn f n+1 A = ( ( ) = ( f n ( ) A n f2n 1 f = 2n f 2n f 2n+1 ) 2 (3.5) (3.6) with n = 1, 2, 3,... Falk and Dyson [26] studied the periods, T (N), of the iteration ( ) A n f2n 1 f mod N = 2n mod N (3.7) f 2n f 2n+1 for various classes of the integers N. T (N) is the smallest, positive, integer, such that ( ) 1 0 A T (N) = mod N (3.8) 0 1 Thus T (N) is, also, the period of the Fibonacci sequence mod N, which is known to be a random function of N cf. fig. 1. We define the scrambling time, t scrambling for the cat map to be equal to the mixing time, t mixing of the discrete dynamical system on AdS 2 [N]. Since T (N) is the period of A mod N, the maximum available time for scrambling, is proportional to T (N)/2, since the mixing time is the time for any initial distribution to become almost flat. t mixing T (N) 2 8 (3.9)

10 Figure 1: The period of the Arnol d cat map as a function of the order of the prime N, for the first 200 primes. for large N. By [26], if N = f 2k, then T (N) = 2k. Therefore, t mixing = k. Now we recall that the solution of the Fibonacci recurrence is given by f n = 1 ( λ n 1 ) (3.10) 5 λ n with λ = (3.11) known as the Golden Ratio. For n = 2k >> 1, f 2k = N (1/ 5) exp(2k log λ), so T (N) ln N. These orbits mod f 2k are short orbits and, in order to get mixing we have to take, correspondingly large, values of k. At the quantum level, the role of short orbits has been connected with that of scars [27]. We remark that the restriction to prime values of N for which 5 is a quadratic residue mod N will make possible the analytic construction of eigenstates and eigenvalues of the QACM. It appears that, up to now, explicit expressions for the eigenstates and eigenvalues of the QACM are not known for generic N [20]. The deterministic, chaotic, orbits of ACM on AdS 2 [N] can be obtained as follows: Choosing as initial point X 0 = σ 3 we obtain, from eq. (3.1), ( ) f2n 1 f X n = 2n+1 + f2n 2 2f 2n f 2n 1 2f 2n f 2n+1 f 2n 1 f 2n+1 f2n 2 mod N (3.12) 9

11 Such short periods of the ACM imply, also, the existence of, non trivial, conservation laws, that is, elements of SL(2, Z N ), that commute with it. These form an abelian group, the commutant, G(A).For prime values of N, it is cyclic, i.e. there exists a primitive element, whose powers generate all the others. Among the elements of this group, obviously, are the powers of A mod N; the non trivial conservation laws are described by the complement thereof. The general element of G(A) has the form ( ) k l C(k, l) = (3.13) l k + l with k, l integers, satisfying the constraint k 2 + kl l 2 = 1. This can be cast in the form of Pell s equation x 2 5y 2 = 1 (3.14) with x = k + (l/2) and y = l/2, in which case l must be even. The trivial conservation laws are given by the Fibonacci numbers, k = f 2n 1, l = f 2n, for all n; in this case, integer solutions of Pell s equation correspond to n = 3m, with m = 1, 2,... [29]. For prime values of N the period of the ACM divides the period of the commutant. If the two periods are equal, the ACM is a primitive element of G(A) and there aren t any non-trivial conservation laws. If they re not, then the ACM is a power of the primitive element of the commutant. This power determines the degeneracies of the quantum ACM, as we shall see in the next section. 4 Chaotic eigenstates of the QACM on AdS 2 Quantum mechanics for a local observer in AdS 2 is defined, once a choice of time evolution has been made. As discussed in the previous section the isometry group, SL(2, R), can be used to identify the time evolution for an observer with an element of this group. Since we have discretized the geometry, locally, the canonical variables, for any observer, are the exponentials of position and momentum operators, that define the generators of the finite Heisenberg Weyl group, HW N (see appendix B). The classical isometry of the discretized geometry is P SL(2, Z N ). Having chosen as time evolution map, A, the Arnol d cat map, we have, in fact, specified the observer and its time is defined as the number of iterations of this map. We can now construct the corresponding quantum evolution map, QACM, choosing, for simplicity, N = p prime. This will be, also, the dimension of the single particle Hilbert space of this observer. The unique Hilbert space, for all observers, is defined by the irrep of the Heisenberg Weyl group, as we discussed in ref. [17]. Every observer in the bulk can reconstruct the algebra of his/her observables from those of the conformal field theory (for AdS 2, this is conformal quantum mechanics) on the boundary using bulk to boundary Green functions. 10

12 To construct the unitary (quantum) evolution operator, U(A), corresponding to the, classical, Arnol d cat map, A, we shall use the Weil representation of SL 2 [p]. The detailed construction of U(A) is given, for completeness, in appendix B. This representation, by construction, is the direct sum of two, irreducible, representations of PSL 2 [p], of dimensions (p + 1)/2 and (p 1)/2. Since the action of A on AdS 2 does not distinguish the action of A from that of A, which is, also, an element of SL 2 [p], it realizes a projective action which is one way, how quantum mechanics on AdS 2 differs from the torus, T 2. Therefore it is necessary to choose one of these two representations, thereby imposing the constraint that it is, also, a representation of PSL 2 [p]. One important property of the quantization procedure is that, for any two elements, A 1, A 2 PSL 2 [p], U(A 1 A 2 ) = U(A 1 )U(A 2 ). This implies the interesting fact that, to calculate the quantum evolution, at time n = 1, 2, 3,..., it suffices to compute U(A n ), which is equal to [U(A)] n realizing a a very big simplification in the calculation of time correlation functions. Therefore, the period of the quantum map is equal to that of the classical map and this determines the degeneracies of the spectra and the conservation laws. In the following we sketch the main steps of the construction of the eigenstates and eigenvalues of the QACM for the technical details cf. appendix B. The basic idea comes from the observation that the classical ACM can be diagonalized over the finite field F p = {0, 1, 2,..., p 1}, if 5 is a quadratic residue mod p. So to avoid unnecessary technical complications, we choose the prime p to be of the form 4k 1, since, in that case, if 5 is a quadratic residue mod p and it is easy to construct 5 mod p. (If 5 isn t a quadratic residue mod p, we must work in the corresponding quadratic extension.) We can then check that a 5 k mod p satisfies a 2 5 mod p. The eigenvalues of A are then, λ ± 3 ± a mod p (4.1) 2 Moreover, there is an element, R SL 2 [p], that diagonalizes A, A = RD A R 1 (4.2) where D A = diag (λ +, λ ). Using the explicit formulae in appendix B, we deduce that and that U(D A ) is the circulant matrix U(A) = U(R)U(D A )U(R) (4.3) l U(D A ) k = δ λ+ k,l = δ k,λ l k, l = 0, 1,..., p 1 (4.4) The eigenstates of U(D A ) are the multiplicative characters, π 0, π 1,..., π p 1, of F p = {1, 2,..., p 1}, given by the expressions k π 0 = δ k,0 k = 0, 1, 2,..., p 1 (4.5) 11

13 0 π n = 0 k π n = e 2πin p 1 Indg(k) p 1 (4.6) where k, n = 1, 2,..., p 1 and Ind g (k) is the discrete logarithm of k with respect to the base g, where g is a primitive element of F p; i.e. g Indg(k) k mod p (4.7) It follows that Ind g (k l) = Ind g (k) + Ind g (l). Having determined the eigenstates, let us now provide the expressions of the eigenvalues. We remark that U(D A ) π 0 = π 0 U(D A ) π n = e 2πi p 1 nindg(λ +) π n (4.8) and can read off the eigenvalues of U(D A ) and, thus, of U(A). The eigenvectors of U(A), ψ n, are given by ψ n = U(R) π n (4.9) This calculation becomes effective using the explicit form of U(R), derived from the Weil representation (cf. appendix B). The period of QACM is the period, T (p), of the ACM and is, also, the order of the element(s), λ ± ; since λ ± are integers in F p, this order divides p 1, the order of F p. So there exists an integer, τ p, such that p 1 = τ p T (p). Therefore, Ind g (λ + ) = τ p. This is, precisely, the degeneracy of the eigenvalues of QACM, which are phases, e iεn. From the above we obtain ε n = 2π p 1 τ pn (4.10) with n = 0, 1, 2,..., p 1. Since n labels the eigenstates, too and τ p /(p 1) = 1/T (p) we can determine the degenerate eigenstates. With these tools we can write explicit expressions for the eigenstates, ψ n, p 1 k ψ 0 = k U(R) 0 p 1 k ψ n = k U(R) l l π n 1 k ψ n = k U(R) l e 2πi p 1 nindg(l) (4.11) p 1 l=0 where k = 0, 1, 2,..., p 1, which will help to understand their chaotic properties. The degeneracies of the spectrum imply the existence of non trivial conservation laws, that reduce the size of the attractor. As discussed in the last part of section 3, we can determine explicitly, depending on p, the commutant of A, G(A). This group is cyclic since we have chosen prime values for p,and, if its order is different from the period of ACM, there is a unique element, B, which generates G(A). The corresponding quantum oeprator, U(B), generates the quantum conservation laws. 12 l=1

14 5 ETH and the scrambling time bound for the QACM Recently there has been a lot of activity around the question of the thermodynamics of closed quantum systems [25]. An important role in this question has been assigned to the specific mechanisms of thermalization of various subsystems. A particularly interesting proposal is the Eigenstate Thermalization Hypothesis(ETH). Here we shall not describe the general set up of the thermodynamic statements of this proposal, but rather, we shall show that, in our particular chaotic, quantum, model of single particle scattering in the near horizon region of an extremal black hole, (some of) the basic ingredients of the ETH scenario can be unambiguously identified. The ETH assumes that the closed quantum system has a complete set of chaotic states in the specific sense that their probabilities are sampled from a Gaussian pdf, while their phases are sampled from a flat pdf. Under these assumptions, it can be shown, quite generally, that any initial pure state of a, typical, subsystem, will evolve to a mixed state, described by the equilibrium thermal density matrix, with an effective temperature defined by the average energy of the system. The time required for, unitary, equilibration of the subsystem, is called the scrambling time and it is, obviously, interesting to study how it is bounded, for various physical systems [28]. This is, obviously, a very interesting framework, to discuss the black hole information paradox, where it has been recently conjectured that such a bound exists, is proportional to (β/2π) log S, with S the entropy and β the (inverse) temperature and that black holes saturate it. In the AdS/CFT approach to this problem, we have the tools to study thermodynamics of gravitational backgrounds through the thermodynamics of the boundary, non gravitational, conformal field theory. This is a consistent description of the thermodynamics of local, gravitational, observers, for which the observables are, indeed, defined, unambiguously, on the boundary. These are the sources for the boundary conformal field theory. Chaos is realized within the ETH, assuming that the dynamics of the closed quantum system is ergodic and mixing. This can be shown, using a random matrix description for the dynamics [25]. For the case of the black hole, chaos can be described by shock wave geometries in the near horizon region. For the thermodynamics of the black hole we expect to obtain a random ensemble of geometries and to determine its measure. Assuming ergodicity and a unitary, gravitational, dynamics, this is equivalent to picking out a typical random geometry, as a gravitational background. Our approach is to use a particular arithmetic, namely modular, discretization of the geometry, which is consistent with unitarity and holography. Such a discretization has the merit of being, automatically, random, at the classical and quantum levels, due to the properties of modular arithmetic, as discussed in the previous sections, so it is reasonable to assume that it is a typical random geometry. 13

15 In the following we shall present arguments that support the statement that the QACM eigenfunctions do satisfy the assumptions of the ETH. The chaotic properties of the eigenstates can be traced back to the chaotic character of the discrete logarithm, Ind g (l). The definition of chaos we shall adopt, which is the only one consistent with computational and algorithmic complexity, is that of algorithmic chaos. The effective computation of the discrete logarithm is a classic example of a non compressible algorithm, i.e. that cannot be done in polynomial time, with respect to the number of the input bits of l. On the other hand, using quantum algorithms, Shor and others have shown that it can be reduced to polynomial complexity. These considerations are consistent with the results of an old but very interesting paper of Ford et al. in ref. [20] that shows that the complexity of the QACM is (log N) 2, in contrast with the classical one which is N. This paper created a lot of discussion in the quantum chaos community (cf. the paper by Berry in Les Houches 1989 [20]). The explicit expressions for the states ψ n are sums of random phases, with fixed, complex, amplitudes. This leads, using the large number theorems, for p, to Gaussian distributions of the state components. In the next section, we will provide numerical evidence for this claim and shall discuss some of the consequences regarding the randomness of the matrix QACM itself. As discussed in section 3 the mixing time for the classical ACM scales as the logarithm of the discretization parameter N whenever N takes values in the Fibonacci sequence. The time required for unitary thermalization of a wavepacket (the scrambling time) is identified here with the mixing time (t scrambling = t mixing ). This is so, because the period of the classical and the quantum ACMs coincide as a result of the construction of U(A). We find thus that the scrambling time of the QACM is proportional to log N when the dimension of the single particle Hilbert space,n, takes values in the sequence of Fibonacci integers. We recall also that the entropy S of the AdS 2 is proportional to N, the spatial extent of the geometry. This leads to the saturation of the scrambling time bound, t scrambling log S, of Hayden-Preskill and Sekino-Susskind. The prefactor, which would be the inverse of the temperature, here is to be replaced by the effective temperature of the closed quantum system according to the ETH scenario. For, while the Hawking temperature of the extremal black hole is zero, the chaotic dynamics of the extremal black hole microstates defines a consistently closed system, since the extremal black hole doesn t radiate. This point deserves a fuller analysis, that will be reported in future work. We may define finally the single particle scattering S matrix as the evolution operator, evaluated at half the period of the QACM. A consequence of the chaotic character of the QACM eigenstates is that this matrix is random and completely delocalizes and scrambles initial Gaussian wavepackets. Closing this section we shall present numerical support of our arguments for the chaotic nature of the eigenstates of QACM (cf. fig. 2 for the ground state for p = 461). These 14

16 results were obtained by using Mathematica codes. For p = is a quadratic residue; but 461 isn t of the form 4k 1; however it is possible to obtain the eigenvalues and eigenstates of the QACM numerically. In addition, for p = 461, there exist two invariant subspaces, of dimensionalities (p + 1)/2 = 231 and (p 1)/2 = 230. The, projective, representation, that is appropriate for AdS 2, as discussed above, is the second one. To highlight the symmetry of the ground state, we haven t projected onto the half length, but display the full length. From this figure it s possible to deduce that the probability Figure 2: Symmetric ground state (in the subspace of dimensionality (p 1)/2) of the quantum Arnol d cat map, for p = 461. distribution function (PDF) for the values of the amplitude squared is, indeed, Gaussian, and the phases are uniformly distributed cf. fig. 3 as well as for the chaotic evolution in time of an initial Gaussian wave packet. This is illustrated in the following figures 4, for p = 461, for which the period of the QACM is found to be equal to 23. These results highlight that, as p grows, the chaotic behavior of the eigenstates is enhanced for subsequences of primes, for which the degeneracy of the spectrum of the QACM remains small equivalently, the ratio of the period of the QACM to p 1 decreases. For Fibonacci dimensions of the Hilbert space of states, the period grows as log p, which saturates the STB. So, for large p, an exponentially small part of the space con- 15

17 Figure 3: The bin values for the squared amplitude and of the phases of the ground state, for p =

18 Figure 4: Evolution of a Gaussian wavepacket, for t = 1, 2, 3, 4. For t = 3 we remark the appearence of oscillations, precursors to chaotic behavior. For 11 > t 4, we observe complete delocalization, but remark that patterns still persist. 17

19 tributes to the chaotic dynamics, the rest are copies of it. This means that there is a large number of islands of chaos. This holds for the dynamics of the probes of the radial and temporal, AdS 2, geometry of the near horizon region; far from the horizon, at distances large compared to log p, in units of the AdS 2 radius, the behavior becomes regular. The detailed crossover remains to be elucidated. Of course this numerical analysis is suggestive and will be completed in future work; but the big picture it defines is expected to be valid. 6 Summary and conclusions In this work we proposed a toy model for the chaotic scattering of single particle wave packets in the radial AdS 2 space time geometry of extremal (or nearly) extremal BHs. In recent discussions of the chaotic scattering the focus has been on the dynamics of the microscopic degrees of freedon on the streched horizon, but it is evident that although the dynamics of longitudinal and transverse scrambling will be different, the time duration will be the same [31]. We were able to discretize the coset structure of the geometry of this space time, by introducing a, modular invariant, infrared and ultraviolet cutoff. We obtained an AdS 2 /CFT1 holography and we provided the eigenstates and eigenvalues for the Quantum Chaotic Arnold catmap,as well as the single particle S-matrix. These eigenstates are chaotic in the sense of the eigenstate thermalization Hypothesis and create mixing and chaos for any infalling Gaussian wave packet. An interesting property of this model is,is that we can fix the dimension of the single particle Hilbert space of states, so that to saturate the scrambling time bound of Hayden-Preskill, Sekino and Susskind for an observer with time evolution defined by the QACM. Our results provide a toy model mechanism to explain, how the incoming information of a pure state can be scattered back as thermal radiation described by a density matrix (through the ETH scenario),while at the same time preserving unitarity in the single particle Hilbert space. For future research along these lines,we think it would be interesting to extend this toy model to the construction of the many particle or field theoretic chaotic scattering S-matrix, on the AdS 2 geometry and also investigate in detail how the ETH works in this case as well as to study the corresponding scrambling time. Finally describing the geometry with finite dimensional, N = p n, Hilbert qudit spaces,we provided a framework of contact with the complexity theory of quantum algorithms and quantum circuits since the finite unitary scattering matrix can be written, as we shall show in a future work, as a tensor product of elementary qudit gates [32]. 18

20 Acknowledgements: EGF and SN acknowledge the warm hospitality of the CERN Theory Division and in particular the organizers of the 2016 CERN Winter School for stimulating exchanges. SN acknowledges the warm hospitality of the Institute of Nuclear and Particle Physics of the NCSR Demokritos. The research of EGF was partly implemented under the ARISTEIA-I action (Code no. 1612, D.654) and title Holographic Hydrodynamics of the operational programme education and lifelong learning and is co-funded by the European Social Fund (ESF) and National Resources. The research of MA was supported in part by the grant MIS ORASY(NSRF ACTION, KRIPIS) of the European Regional Development Fund. A Coset geometry and light cones on AdS 2 Iin the near horizon region of spherically symmetric 4d extremal black holes the geometry is known to be of the form AdS 2 S 2, where the AdS 2 = SL(2, R)/SO(1, 1, R), factor describes the geometry of the radial and time coordinates and S 2 is the horizon surface ref(sen). The AdS 2 spacetime, is a one-sheeted hyperboloid defined through its global embedding in Minkowski spacetime with one space and two time like dimensions by the equation x x 2 1 x 2 2 = 1 (A.1) The boundaries of AdS 2 consist of two time like disconnected circles, where AdS 2 approaches, asymptotically, the light cone of M 1,2, x x 2 1 x 2 2 = 0 (A.2) AdS 2 is, at the same time, the homogeneous space, SO(1, 2)/SO(1, 1). This case, is special, in that SO(1, 2) has a double cover, SL(2, R), so we have AdS 2 = SL(2, R)/SO(1, 1). In order to establish our notation and conventions, we proceed with the Weyl construction of the double covering group, SL(2, R). To every point, x µ AdS 2, µ = 0, 1, 2, we assign the traceless, real, 2 2 matrix ( M(x) x 0 x 1 + x 2 x 1 x 2 x 0 Its determinant is, det M(x) = x 2 0 x x 2 2 = 1. The action of any A SL(2, R) on AdS 2 is defined through the mapping ) (A.3) M(x ) = AM(x)A 1 (A.4) 19

21 This induces an SO(1, 2) transformation on (x µ ) µ=0,1,2, x L(A)x (A.5) Choosing as the origin of coordinates, the base point p (1, 0, 0), its stability group SO(1, 1), is the group of Lorentz transformations in the x 0 = 0 plane of M 1,2 or equivalently, the scaling subgroup, D, of SL(2, R) ( ) λ 0 D S(λ) 0 λ 1 (A.6) for λ R. For this choice of the stability point, we define the coset, h A, by decomposing A as A = h A S(λ A ) (A.7) Thus, we associate uniquely to every point x AdS 2 the corresponding coset representative h A (x). We introduce now, the global coordinate system,defined by the straight lines that generate AdS 2 and for which it can be checked easily that they form its complete set of light cones. Consider the two lines, l ± (p), passing through the point p M 1,2, orthogonal to the x 0 axis and at angles ±π/4 to the x 1 = 0 plane. They are defined by the intersection of AdS 2 and the plane x 0 = 1 cf. fig. 5. The coordinates of any point, q + l + (p), q l (p) are given as, (1, µ ±, ±µ ± ), µ ± R correspondingly. We can parametrize any point x µ, of AdS 2, by the intersection of the local light cone lines, l ± (x), with coordinates µ ± and φ ± through the relations x 0 = cos φ ± µ ± sin φ ± x 1 = sin φ ± + µ ± cos φ ± x 2 = ±µ ± (A.8) These can be inverted as follows: e iφ ± = x 0 + ix 1 1 ± ix 2 µ ± = ±x 2 (A.9) The geometric meaning of the coordinates φ and µ is that µ parametrizes the x 2, space like, coordinate and, thus, µ ± 2 parametrizes the light cone lines l± (x). The angle φ ± is the azimuthal angle of the intersection of l ± (x) with the plane (x 0, x 1 ). From eq. (A.9), by re-expressing numerator and denominator in polar coordinates, we find φ = τ σ (A.10) 20

22 where τ and σ are the arguments of the complex numbers x 0 + ix 1 and 1 + ix 2. The corresponding coset parametrization (group coset motion which brings the origin to the point x) is: h(µ ±, φ ± ) = R(φ ± )T ± (µ ± ) (A.11) where and R(φ) = ( cos φ/2 sin φ/2 sin φ/2 cos φ/2 T + (µ) = [T ( µ)] T = ( 1 µ 0 1 ) ) (A.12) (A.13) It is easy to see also, that T ± (µ ± ), acting on the base point X(p), generate the light cone l ± (p), so we identify these one parameter subgroups with the light cones at p. Figure 5: The light cone of AdS 2 at p = (1, 0, 0). 21

23 In the literature the study of fields on AdS 2 requires an extension, to the universal covering of this spacetime, ÃdS 2, together with appropriate boundary conditions, in order to avoid closed time like geodesics and reflection of waves from the boundary. This extension can be parametrized using as time coordinate the azimuthal angle τ, by extending its range, ( π, π) to ( π + 2πk, π + 2π(k + 1)), k = ±1, ±2,... and by the space coordinate σ ( π/2, π/2), defined in eqs. (A.10). The extension of the range of τ parametrizes the infinitely sheeted Riemann surface of the function log( ), used in deriving eq. (A.10). These parametrizations induce specific metrics on AdS 2 and ÃdS 2, namely: ds 2 = (1 + µ 2 )dφ 2 + 2dφdµ (A.14) and ds 2 = 1 [ dτ 2 + dσ 2] (A.15) cos 2 σ The latter metric, with τ σ + φ R describes the Einstein strip, that has two disconnected boundaries at σ ±π/2. The Einstein strip is the union of an infinite number of copies of AdS 2. It is interesting to note that the coset structure of AdS 2 can be elevated to ÃdS 2 by using the universal covering group of SL(2, R), which has been explicitly constructed in ref. [33]. This point deserves a discussion in its own right and will be reported in future work. B The Heisenberg Weyl group and the Weil representation of PSL 2 [p] Detailed references to this and the following appendix can be found in refs. [17] and [34]. The Finite Heisenberg group HW p, is defined as the set of 3 3 matrices of the form g(r, s, t) = r 1 0 (B.1) t s 1 where r, s, t belong to Z p (integers modulo p), where the multiplication of two elements is carried modulo p. When p is a prime integer there is a unique p-dimensional unitary irreducible and faithful representation of this group, given by the following matrices J r,s,t = ω t P r Q s (B.2) where ω = e 2πi/p, i.e. the p th primitive root of unity and the matrices P, Q are defined as P kl = δ k 1,l Q kl = ω k δ kl (B.3) 22

24 where k, l = 0,..., p 1. It is to be observed that, if ω is replaced with ω k, for k = 1, 2,..., p 1 all the relations above remain intact. Since p is prime all the resulting representations are p-dimensional and inequivalent. The matrices P, Q satisfy the fundamental Heisenberg commutation relation of Quantum Mechanics in an exponentiated form Q P = ω P Q (B.4) In the above, Q represents the position operator on the circle Z p of the p roots of unity and P the corresponding momentum operator. These two operators are related by the diagonalising unitary matrix F of P, QF = F P so F is the celebrated Discrete Fourier Transform matrix F kl = 1 p ω kl, with k, l = 0,..., p 1 (B.5) (B.6) An important subset of HW p consists of the magnetic translations J r,s = ω rs/2 P r Q s (B.7) with r, s = 0,..., p 1. These matrices are unitary (J r,s = J r, s ) and traceless, and they form a basis for the Lie algebra of SL(p, C). They satisfy the important relation J r,s J r,s = ω(r s rs )/2 J r+r,s+s (B.8) This relation implies that the magnetic translations form a projective representation of the translation group Z p Z p. The factor of 1/2 in the exponent of (B.8) must be taken modulo p. The SL 2 (p) appears here as the automorphism group of magnetic translations and this defines the Weil metaplectic representation. If we consider the action of an element ( ) a b A = (B.9) c d on the coordinates (r, s) of the periodic torus Z p Z p, this induces a unitary automorphism U(A) on the magnetic translations, since the representation of Heisenberg group is unitary and irreducible, U(A)J r,s U (A) = J r,s (B.10) where (r, s ) are given by ( a b (r, s ) = (r, s) c d 23 ) (B.11)

25 This relation determines U(A) up to a phase and in the case of A SL 2 [p], the phase can be fixed to give an exact (and not projective) unitary representation of SL 2 [p]. The detailed formula of U(A) has been given by Balian and Itzykson [35]. Depending on the specific values of the a, b, c, d parameters of the matrix A, we distinguish the following cases: δ 0 : δ = 0, b 0 : δ = b = 0, c 0 : δ = b = 0 = c = 0 : U(A) = σ(1)σ(δ) p U(A) = σ( 2b) p U(A) = σ(2c) p U(1) = I r,s ω br2 +(d a)rs cs 2 2δ s ω s 2 2b J s(a 1)/b,s r ω r 2 2c P r where δ = 2 a d and σ(a) is the quadratic Gauss sum given by J r,s (B.12) σ(a) = 1 p 1 ω ak2 = (a p) p k=0 { 1 for p = 4k + 1 i for p = 4k 1 (B.13) while the Legendre symbol takes the values (a p) = ±1 depending on whether a is or is not a square modulo p. It is possible to perform explicitly the above Gaussian sums noticing that (J r,s ) k,l = δ r,k l ω k+l 2 s (B.14) where all indices take the values k, l, r, s = 0,..., p 1. This has been done in [23, 24]. In the case δ = 2 a d 0 mod p and c 0 mod p, the result is U(A) k,l = ( 2c p) { } 1 ω ak2 2kl+dl 2 2c (B.15) p i If c 0 mod p, then we transform the matrix A to one with c 0 mod p. The cases δ 0 mod p can be worked out easily using the expressions of the matrix elements of J r,s, given in (B.2). It is interesting to notice that redefining ω to become ω k for k = 1, 2,..., p 1, the matrix U(A) transforms to the matrix U(A k ), where A k is the 2 2 matrix A k = ( ) a bk c/k d, which belongs to the same conjugacy class with A as long as k is a quadratic residue. If k = p 1 we pass from the representation U(A) to the complex conjugate one U(A). The Weyl representation presented above, provides the interesting result that the unitary matrix corresponding to the SL 2 (p) element a = ( ) is -up to a phase- the Discrete Finite Fourier Transform (B.6) U(a) = ( 1) k+1 i n F where n = 0 for p = 4k + 1 and n = 1 for p = 4k 1. 24