The quantum cat map on the modular discretization of extremal black hole horizons

Based on our recent work on the discretization of the radial AdS$_2$ geometry of extremal BH horizons,we present a toy model for the chaotic unitary evolution of infalling single particle wave packets. We construct explicitly the eigenstates and eigenvalues for the single particle dynamics for an observer falling into the BH horizon, with time evolution operator the quantum Arnol'd cat map (QACM). Using these results we investigate the validity of the eigenstate thermalization hypothesis (ETH), as well as that of the fast scrambling time bound (STB). We find that the QACM, while possessing a linear spectrum, has eigenstates, which are random and satisfy the assumptions of the ETH. We also find that the thermalization of infalling wave packets in this particular model is exponentially fast, thereby saturating the STB, under the constraint that the finite dimension of the single--particle Hilbert space takes values in the set of Fibonacci integers.


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].
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 provides a new framework to study mechanisms by which past and recent memories of the shock-wave spacetime geometry structure are encoded in the angular and time correlations of the emitted Hawing radiation.
Quite recently the non-unitarity of the Hawking radiation has been interpreted as the result of integrating all the BMS soft graviton cloud, accompanying the hard Hawking quanta [3].
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, is scrambled, 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].
So it is imperative to separate the issue of the chaotic dynamics and geometry of the near horizon region from any issues regarding Hawking radiation.
It is possible to do so, when studying probes of extremal black hole horizons; for, in that case, since the Hawking temperature vanishes, while the entropy does not, there is no Hawking radiation.
In this case, it is known that the near horizon geometry is described by a metric, that factorizes into a radial-temporal part that can be identified as an AdS 2 manifold, while the angular part describes the charges.
It is possible to construct a model for random AdS 2 geometries, inspired by the shock wave geometries, necessary for the chaotic dynamics of the probes, by introducing a consistent discretization of the AdS 2 near horizon geometry of nearly extremal black holes. This is the socalled the modular discretization, AdS 2 [p], for every prime integer p [15]. In this framework, the entropy of the black hole, is identified with the Kolmogorov-Sinai entropy of the deterministic, chaotic, dynamics of the geometric, microscopic, degrees of freedom, defining the near horizon geometry.
AdS 2 [p] is a specific discrete deformation of its continuous counterpart.It has a random structure due to the modular arithmetic. As explained in ref. [15], by explicit calculation, this specific discretization is chosen among many possible discretizations because it provides a way of constructing an 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 discrete and finite symmetry group of AdS 2 [p], which is PSL 2 [p]: it acts as an isometry group of the bulk and as the (Möbius) conformal group on the boundary.
This discrete geometry provides also a natural framework for describing the single particle dynamics, via observers, with time evolution operators that are elements of the isometry group. This is a discrete analog of the superconformal quantum mechanics of probes near the horizons of large extremal black holes [16].
In the present work we specify the infalling, accelerating, observer by the well known Arnol'd cat map (ACM) [17]. This map defines a, particular, observer, ACM, who, by performing single-particle scattering on the horizon of the black hole, can probe the randomness of the geometry and it is consistent with the isometries of the background since it belongs to the discrete isometry group of AdS 2 [p]. The discreteness of the geometry implies that the global coordinates of AdS 2 [p] are discrete. On the other hand, the time, measured by the ACM observer, is the iteration step of the corresponding map.
In this work we study the quantum dynamics of the probe in this discrete, background geometry.
It is important to stress that both, probe and background geometry, have a finite dimensional space of states. What we study is how superpositions of the states of the probe evolve, when the background geometry is found in any given, fixed state.
The ACM can be "quantized", i.e. it is possible to define 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 [18,19,20,21,22].
An introduction to the requisite tools from arithmetic geometry and computational number theory can be found in ref. [23].
Next we proceed with the plan of the paper: In section 2 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 [p], i.e. the set of elements of PSL 2 [p], that commute with it.
In section. 3 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]. This particular representation defines, for every prime, p, the Hilbert space of states of the infalling wavepackets on AdS 2 [p], while the other one is appropriate for the case of the torus.
We determine analytically the spectrum and the eigenstates of the QACM and we compute their degeneracies.
We find the interesting result that, while the spectrum is linear, the eigenstates are chaotic in a very 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 are known to be the premises for the eigenstate thermalization hypothesis (ETH) for quantum ergodicity or unitary thermalization [24].
In section 4 we review the Eigenstate Thermalization Hypothesis and we stress that its premises can be checked to hold within ou our framework. We use the results obtained in section 3 to study the spectrum of scrambling times and we find that for Fibonacci integer values the scrambling time bound of Susskind and Sekino is saturated.
Finally in section 5 we discuss our results and open problems for future work.
In appendix A we collect all the necessary material for the detailed construction of the Weil representation of SL 2 (p) and PSL 2 (p) and we present the technical details for the analytic construction of the eigenstates and eigenvalues of the QACM.

The Arnol'd cat map and Fibonacci chaos on AdS [N ]
We review the description of the modular discretization of AdS 2 [15].
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 is a finite and random set of points in the embedding Minkowski spacetime. When N is prime this is an "arithmetic geometry" in the mathematical literature [23], that is, a geometry over a finite field.
This discretization can be used to describe nonlocality, chaos and quantum information processing in the vicinity of the BH horizon, as well as defining a discrete version of the AdS 2 /CFT 1 holography.
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 This set-an example of which is shown in fig. 1-is constructed by noting that, for any triplet of integers, (k, l, m), that satisfy eq. (2.1) mod N , there exists an integer M ≡ 1 mod N , such that the triple of rational points (k/M, l/M, m/M ) ≡ (x 0 , x 1 , x 2 ), satisfies the equation it defines rational points of the continuum AdS 2 manifold. If we fix an "infrared cutoff", L, for x 2 , |x 2 | ≤ L, by increasing the denominator M , we can obtain in this way a rational approximant to the continuum AdS 2 geometry.
A straightforward prescription for constructing all the solutions of eq. (2.1), i.e. the points of AdS 2 [N ], is as follows: where a 2 + b 2 ≡ 1 mod N and a, b, µ ∈ {0, 1, 2, . . . , N − 1} [15]. Thus, the discretized, spatial part, along x 2 , consists of N points and the Hilbert space of single-particle states has dimension N . The global, AdS 2 , time, is parametrized by the points of the discrete circle, a 2 + b 2 ≡ 1 mod N . The proper time of the ACM observer is identified with the number of iterations of the ACM mod N . Due to the mod N arithmetic, the global, AdS 2 time and the proper time of the ACM are periodic.
If we assign, to each spatial point, a two-state system, the microscopic degrees of freedom of the near-horizon discrete geometry, have a Hilbert space of dimension ∝ 2 N . From this we conclude that the entropy, S BH , of such configurations, is proportional to log 2 N = N log 2, which is the number of spatial points.
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 [22].
The finite geometry, AdS 2 [p], has as isometry group the finite projective 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).
It is easy to find the number of points of AdS 2 [N ], for any integer N . Numerical experiments suggest the following recursion relation for the number of points of where Sol(p) = p(p + 1) and k = 1, 2, . . . for any prime integer p.
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 [28].
Next we discuss the discrete time evolution of the motion of a particle on AdS 2 [N ]. To every point x µ ∈ AdS 2 [N ], where µ = 0, 1, 2, can be assigned a traceless, 2 × 2, matrix X X ≡ x 0 The discrete time evolution, for an observer, defined by its evolution matrix, A ∈ PSL 2 [N ], is given by the recursion relation where n = 0, 1, 2, . . . labels the (stroboscopic) time of the observer and X 0 is the initial point of the trajectory. We make a specific choice, introduced in ref. [15], of an observer, described by the Arnol'd cat map(ACM): The map corresponds to successive kicks, forwards and backwards along the light cone of AdS 2 [N ], since and is an element of PSL 2 (Z N ). We point out that this action differs from the action of ACM on the torus, which is linear in A [17]. 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 ACM, acting on AdS 2 [N ], according to (2.5), induces a discrete Lorentz transformation. The mod N prescription guarantees mixing and, thus, ergodicity of the dynamics.
• The mod N prescription provides chaotic orbits on AdS 2 [N ] but, since ACM has a finite period, T (N ), depending, randomly, on N , all of these orbits are, also, periodic. Their chaotic nature can be seen up to evolution time ≤ T (N )/2.
• The ACM has been studied intensively also as a toy model for semiclassical quantum chaos on the toroidal phase space [18], although the degeneracies in its spectrum impose additional constraints on its quantum ergodic properties [26].
Here we extend the study to the classical and quantum motion of particles under the Arnol'd cat map on the discretized AdS 2 [N ] geometry.
As we pointed out above, this is a discretized deformation of the continuous, radial and time, geometry of the near horizon region of extremal BHs [29].
The motion we study is the longitudinal motion of probes and it differs from the motion along the horizon of the black hole, 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].
An important property of ACM is that it is known to generate the sequence of Fibonacci numbers, f n , with n = 1, 2, . . ., defined by where f 0 = 0 and f 1 = 1. We observe that with n = 1, 2, 3, . . . Falk and Dyson [25] studied the periods, T (N ), of the iteration for various classes of the integers N . T (N ) is the smallest, positive, integer, such that 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. 2. We now turn to the discussion of the scrambling time, t scrambling , of the black hole horizon geometry, which has been introduced in [27] as the time necessary for uniform spreading of the distribution of the microscopic degrees of freedom, near the black hole horizon, under an external perturbation. Here we use the probe approximation for an infalling, single-particle, wavepacket as such a perturbation and we assume that the scrambling time of the black hole horizon is the same as that of the wavepacket. For a Gaussian wavepacket the scrambling time is the time necessary for its uniform spreading along the horizon. Moreover, we focus on the radial dependence of the spreading, for which the scrambling time is the same as that of the transverse spreading.
More technically, the scrambling time, t scrambling , as defined above, is identical with the mixing time, t mixing of the dynamical system, here the cat map, on AdS 2 [N ] [17]. Since T (N ) is the period of A mod N , the maximum available time for scrambling is proportional to T (N )/2, therefore From the work of Falk and Dyson [25], if N = f 2k , then T (N ) = 2k. Therefore, t scrambling = t mixing ≤ k. We recall that the solution of the Fibonacci recurrence is given by These orbits mod f 2k are "short" orbits and, in order to get mixing, we have to take "large" values of k. At the quantum level, the role of "short" orbits has been connected with that of "scars" [26].
Such "short" periods of the ACM imply 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 with k, l integers, satisfying the constraint k 2 + kl − l 2 = 1. This can be cast in the form of Pell's equation .17) 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, . . . [23].
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. The deterministic, chaotic, orbits of ACM on AdS 2 [N ] can be obtained as follows: If we take as initial point X 0 ≡ (x 0 , x 1 , x 2 ), we find the corresponding sequence, {X n }, Writing out the results we find that (x 2 ) is given by the action of elements L n ∈ SO(2, 1), with integer coefficients mod N , that act on the initial point (x 0 , x 1 , x 2 ), for every time step n: It is noteworthy that L n has fractional coefficients, which means, even without applying the mod operation, that the corresponding combinations of integer coordinates on the hyperboloid are even. Also, that, for even N , these matrices must be defined separately.
The relevance of these expressions is that they highlight the, classical, deterministic, chaotic dynamics of the ACM on AdS 2 [N ].
Indeed, the motion of a particle under ACM, along the spatial direction, x 2 , of AdS 2 [N ] can be seen to be fully chaotic and mixing, cf. fig. 3. mod 683 as a function of the time step n. The period is found to be equal to 684 = N + 1.

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 operator 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 .
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. [15]. 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.
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 PSL 2 [p].
The detailed construction of U (A) is given, for completeness, in appendix A, for the group SL 2 [p]. This representation, by construction, is the direct sum of two, irreducible, representations, 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, it realizes a projective action-which is 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, . 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)] nrealizing 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.
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 [18]; so our results, for these primes, are new.
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 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 check, in this case, that a ≡ 5 k mod p satisfies a 2 ≡ 5 mod p. The eigenvalues of A are then, Moreover, there is an element, R ∈ SL 2 [p], that diagonalizes A, where D A = diag (λ + , λ − ). We can deduce that 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 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.
Having determined the eigenstates of U (D A ), let us now provide the expressions of the eigenvalues. We remark that p−1 nIndg(λ + ) |π n (3.8) and can read off the eigenvalues of U (A). The eigenvectors of U (A), |ψ n , are given by |ψ n = U (R)|π n (3.9) This calculation becomes effective using the explicit form of U (R), derived from the Weil representation (cf. appendix A.) 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 τ p n (3.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 , 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 2, 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.

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 [24].
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) [24]: The time average of any observable of a subsystem of a closed quantum system, which is, initially, in a pure state, for large times, converges to the thermal average of the observable, along with exponentially small corrections, O(e −S ), where S is the entropy, defined by the thermal density matrix of the system. Since the total system is closed, it can be taken in a pure state and the temperature in the thermal density matrix is an effective temperature, defined by the energy average in the initial state of the total system. It has been shown [24] that a way to realize this hypothesis is to assume 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.
These are the basic premises for thermalization to be possible. 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.
As we shall show below, the system of the black hole near horizon geometry, including the infalling wavepacket as a subsystem, has a complete set of chaotic eigenstates with, precisely, these properties.
An important quantity that describes thermalization is the time required for thermalization. Recent studies have highlighted the relevance of the Hamiltonian dynamics of integrable and chaotic systems for determining the thermalization time.
On the other hand, for unitary, thermalization of the subsystem, it is interesting to study how-and if-the time required, is bounded, from below, as a function of the entropy, for various physical systems [27].
It has been conjectured that black hole horizons, considered, along with their probes, as closed quantum systems, are the fastest scramblers [27].
ETH, therefore, is a very interesting framework within which to discuss these issues. It has been recently conjectured that such a bound exists and that it 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 [24]. 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 to this end is to use a particular arithmetic, namely modular, discretization of the geometry, which, while random, is consistent with unitarity and holography.
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 [23]. On the other hand, using quantum algorithms, Shor and others have shown that it can be reduced to polynomial complexity [30] These considerations are consistent with the results of an old but very interesting paper of Ford et al. in ref. [18] 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 [18]).
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 2 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 mixing = 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, which depends on the complete set of its chaotic eigenstates, 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. 4 for the ground state for p = 461). These results were obtained by using Mathematica codes. For p = 461, 5 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 distribution function (PDF) for the values of the amplitude squared is, indeed, Gaussian, and the phases are uniformly distributed cf. fig. 5.
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 spacetime contributes to the chaotic dynamics, the rest are copies of it. This means that there is a large number of "islands of chaos". ETH thus holds within each such island separately.
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. Figure 4: The squared amplitude, P n ≡ |a n | 2 , vs. the component label, 0 ≤ n ≤ p − 1, in the symmetric ground state (in the subspace of dimensionality (p − 1)/2) of the quantum Arnol'd cat map, for p = 461.

Summary and conclusions
In this work we proposed a toy model for the chaotic scattering of single particle wave packets in the modular discretization of 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. 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 /CFT 1 holography and we provided the eigenstates and eigenvalues for the quantum chaotic Arnold cat map, as well as the single particle S-matrix.
These eigenstates are chaotic in the sense of the eigenstate thermalization hypothesis. They create mixing and chaos for any infalling Gaussian wave packet. An interesting property of this model 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 Figure 5: The binned histograms for the squared amplitude, |a n | 2 , and for the (relative) phase of the amplitude in the ground state, for p = 461.

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 modular AdS 2 [N ] 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 multi-qudit spaces, we provided a framework of contact with the complexity theory of quantum algorithms and quantum circuits for the AdS/CFT correspondance [32], 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.
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-448332-ORASY(NSRF 2007-13 ACTION, KRIPIS) of the European Regional Development Fund.
A The Weil representation of PSL 2 [p] and the construction of the QACM eigenstates Detailed references to this and the following appendix can be found in refs. [15] and [34]. The finite Heisenberg-Weyl group HW p , is defined as the set of 3 × 3 matrices of the form 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 (A. 2) where ω = e 2πi/p , i.e. the p th primitive root of unity and the matrices P, Q are defined as 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 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 , so F is the celebrated Discrete Fourier Transform matrix An important subset of HW p consists of the magnetic translations J r,s = ω rs/2 P r Q s (A.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 (A.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 (A.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 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: where δ = 2 − a − d and σ(a) is the quadratic Gauss sum given by 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 (A.14) where all indices take the values k, l, r, s = 0, . . . , p − 1. This has been done in [21,22]. In the case δ = 2 − a − d = 0 mod p and c = 0 mod p, the result is 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 (A.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 Weil representation presented above, provides the interesting result that the unitary matrix corresponding to the SL 2 (p) element a = 0 −1 To obtain the block diagonal form of the above matrices U ± (A), we rotate with the orthogonal matrix of the eigenvectors of S. This p-dimensional orthogonal matrix, dubbed here O p , can be obtained in a maximally symmetric form (along the diagonal as well as along the anti-diagonal) using the eigenvectors of S in the following order: In the first (p + 1)/2 columns we put the eigenvectors of S of eigenvalue equal to 1, and in the next (p − 1)/2 columns the eigenvectors of eigenvalue equal to −1 in the specific order given below: (e 0 ) k = δ k0 , (e + j ) k = 1 √ 2 (δ k,j + δ k,−j ), j = 1, . . . , p−1 2 (e − j ) k = 1 √ 2 (δ k,j − δ k,−j ), j = p+1 2 , . . . , p (A. 19) where k = 0, . . . , p − 1. Different orderings of eigenvectors may lead to different forms of the matrices U ± (A). The so obtained orthogonal matrix O p has the property In the following using the above material we shall provide the details of the construction of the eigenstates of the QACM.
As we discussed in section 3, the first step consists in diagonalizing the ACM in SL 2 [p] and this can be done by an element R ∈ SL 2 [p] given by where a is the square root of 5 mod p. The matrix elements of U (A) and U (R) can be constructed explicitly, for any prime p, in particular for p ≡ 3 mod 4: where k, l = 0, 1, 2, . . . , p − 1 and ω = exp(2πi/p). Plugging these expressions in eq. (3.11) we obtain the explicit forms of the QACM eigenstates, k|ψ n , k, n = 0, 1, 2, . . . , p − 1: We can project U (A) and its eigenstates onto the two irreducible subspaces, of dimension (p ± 1)/2, obtaining the block-diagonal forms V ± (A) .26) and their corresponding eigenstates, |ψ n ± |ψ n ± = O p 1 2 (I ± S) |ψ n (A. 27) respectively. It can be checked that V ± (A) have the same period as U (A). These two irreducible representations of SL2[p] are both appropriate for the torus, T 2 , not, however, for realizing the projective action of SL 2 [p] on AdS 2 , as discussed in section 3. We must choose that one of the two, which is, also, a (irreducible) representation of PSL 2 [p] [33]. where n = 0 for p = 4k + 1 and n = 1 for p = 4k − 1. For PSL 2 [p] the element a 2 = −I is identified with the identity matrix, I. Therefore, we should choose, among the two irreducible representations of SL 2 [p], of dimension (p + 1)/2 and (p − 1)/2, that one, for which U (a 2 ) = I.
We can easily check that this happens for the p+1 2 dimensional representation, when p ≡ 1 mod 4, and for the p−1 2 −dimensional one, when p ≡ 3 mod 4. In our construction we found it simpler to work with primes of the latter form, therefore the representation is that defined by V − (A), with eigenstates, |ψ n − .