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

Based on our recent work on the discretization of the radial AdS2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {AdS}_2$$\end{document} 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 as 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][3][4][5][6].
It seems possible in principle that the horizon region of the 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][2][3]. a e-mail: axenides@inp.demokritos.gr b e-mail: mflorato@phys.uoa.gr c e-mail: stam.nicolis@lmpt.univ-tours.fr In the language of Refs. [4][5][6], 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 a result of integrating the whole BMS soft graviton cloud, accompanying the hard Hawking quanta [7,8].
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 [22][23][24][25][26][27].
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 well known that the near-horizon geometry is described by a metric that factorizes into a radialtemporal 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 so-called modular discretization, AdS 2 [ p], for every prime integer p [28][29][30][31]. 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 Refs. [28][29][30][31], 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 [32].
In the present work we specify the infalling, accelerating, observer by the well-known Arnol'd cat map (ACM) [33][34][35]. 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 Weyl 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 [36][37][38][39][40][41][42][43].
An introduction to the requisite tools from arithmetic geometry and computational number theory can be found in Refs. [44][45][46].
Next we present the plan of the paper.
In Sect. 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 Sect. 3 we construct explicitly the exact quantum Arnol'd cat map, using the metaplectic (or Weyl) 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 [47][48][49][50][51][52][53][54].
In Sect. 4 we review the Eigenstate Thermalization Hypothesis and we stress that its premises can be checked to hold within our framework. We use the results obtained in Sect. 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 Sect. 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 Weyl 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.
We define the modular discretization by replacing the set of real numbers, R, by the set of integers modulo N . The coset finite geometry AdS 2 [N ] = SL 2 (Z N )/SO(1, 1, Z N ) so obtained is a discrete deformation of its continuous counterpart, AdS 2 [R] = SL 2 (R)/SO(1, 1, R). AdS 2 [N ] is a finite and random set of points in the embedding Minkowski space-time. When N is prime this is an "arithmetic geometry" in the mathematical literature [44][45][46], 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. 1is 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.
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 [43].
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 PSL(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 AdS 2 [ p k ], Sol( p k ): where Sol( p) = p( p + 1) and k = 1, 2, . . . for any prime integer p.
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 [55].
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 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 Refs. [28][29][30][31], 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 it 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 [33][34][35]. 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 [36][37][38][39], although the degeneracies in its spectrum impose additional constraints on its quantum ergodic properties [56][57][58][59][60].
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 [61][62][63].
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 [64,65].
An important property of ACM is that it is well 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

Falk and Dyson [66] 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 well 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 [67][68][69][70][71][72] 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 to the mixing time, t mixing of the dynamical system, here the cat map, on AdS 2 [N ] [33][34][35]. 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 [66], 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" [56][57][58][59][60].
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, 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, . . . [44][45][46]73].
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 are no non-trivial conservation laws at all. If they are not equal, 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 the initial point 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 we note 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.

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 PSL(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 Refs. [28][29][30][31]. 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 Weyl 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) n -realizing 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 [36][37][38][39]; 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 is not 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 where D A = diag (λ + , λ − ). We can deduce that and that U (D A ) is the circulant matrix l|U (D A )|k = δ λ + k,l = δ k,λ − l k, l = 0, 1, . . . , p − 1.
(3.4) 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.
It follows that Ind g (k · l) = Ind g (k) + Ind g (l).
Having determined the eigenstates of U (D A ), let us now provide the expressions of the eigenvalues. We remark that This calculation becomes effective using the explicit form of U (R), derived from the Weyl representation (cf. Appendix A).
The period of QACM is the period, T ( p), of the ACM; and it is also the order of the element(s), λ ± . Since the λ ± 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 Sect. 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 operator, U (B), generates the quantum conservation laws.
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) [47][48][49][50][51][52][53][54]: 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 to be 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 [47][48][49][50][51][52][53][54] 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 singleparticle 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 [67][68][69][70][71][72].
ETH, therefore, is a very interesting framework within which to discuss these issues.
It has recently been 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, nongravitational, 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 [47][48][49][50][51][52][53][54]. 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 a 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 [44][45][46]73]. On the other hand, using quantum algorithms, Shor and others have shown that it can be reduced to polynomial complexity [74] These considerations are consistent with the results of an old but very interesting paper of Ford et al. in Refs. [36][37][38][39], which 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 [36][37][38][39]). 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 Sect. 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 singleparticle 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 does not radiate. This point deserves a fuller analysis, which will be reported in future work.
We may define finally the single-particle scattering Smatrix 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 is not 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, which is appropriate for AdS 2 , as discussed above, is the second one. To highlight the symmetry of the ground state, we have not projected onto the half length, but we display the full length. From this figure it is 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.  Fig. 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 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 space-time 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.

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 freedom on the stretched horizon. It is evident that, although the dynamics of longitudinal and transverse scrambling will be different, the time duration will be the same [64,65].
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 Arnol'd 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 singleparticle Hilbert space of states, so as 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 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 frame-work of contact with the complexity theory of quantum algorithms and quantum circuits for the AdS/CFT correspondence [75], since the finite unitary scattering matrix can be written, as we shall show in future work, as a tensor product of elementary qudit gates.

A The Weyl 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. [28][29][30][31] and [76].
The finite Heisenberg-Weyl group, H W 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 out 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: where ω = e 2πi/ p , i.e. the pth 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 diagonalizing unitary matrix F of P, so F is the celebrated discrete Fourier transform matrix, The detailed formula of U (A) has been given by Balian and Itzykson [77]. Depending on the specific values of the a, b, c, d parameters of the matrix A, we distinguish the following cases: 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 where all indices take the values k, l, r, s = 0, . . . , p − 1. This has been done in [42,43]. 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 Weyl representation presented above provides the interesting result that the unitary matrix corresponding to the 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: 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 Sect. 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: and their corresponding eigenstates, |ψ n ± |ψ n ± = O p 1 2 (I ± S) |ψ n , (A. 27) respectively. It can be checked that the V ± (A) have the same period as U (A).
These two irreducible representations of SL2[ p] are both appropriate for the torus, T 2 but they are not, however, for realizing the projective action of SL 2 [ p] on AdS 2 , as discussed in Sect. 3. We must choose that one of the two which is, also, an (irreducible) representation of PSL 2 [ p] [78].
We can obtain irreducible representations of PSL 2  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, the 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 2dimensional 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 − .