## A Stochasticity Threshold in Holography and the Instability of AdS

## Pallab Basu, a Chethan Krishnan, b Ayush Saurabh a,b,c

a International Center for Theoretical Sciences,

Indian Institute of Science Campus,

Bangalore - 560012,

India

b Center for High Energy Physics,

Indian Institute of Science,

Bangalore - 560012,

India

c School of M. A. C. E., The University of Manchester,

Manchester - M13 9PL, United Kingdom

E-mail: pallabbasu@gmail.com , chethan.krishnan@gmail.com ,

ayushsaurabh@hotmail.com

Abstract: We give strong numerical evidence that a self-interacting probe scalar field in AdS, with only a few modes turned on initially, will undergo fast thermalization only if it is above a certain energetic threshold. Below the threshold the energy stays close to constant in a few modes for a very long time instead of cascading quickly. This indicates the existance of a Strong Stochasticity Threshold (SST) in holography. The idea of SST is familiar from certain statistical mechanical systems, and we suggest that it exists also in AdS gravity. This would naturally reconcile the generic non-linear instability of AdS observed by Bizon and Rostworowski, with the Fermi-Pasta-Ulam-Tsingou-like quasi-periodocity noticed recently for some classes of initial conditions. We show that our simple set-up captures many of the relevant features of the full gravity-scalar system.

Keywords: Gravitational Collapse in AdS/CFT, AdS Instability

## Contents

|   1 | Introduction and Conclusion                |   1 |
|-----|--------------------------------------------|-----|
| 2   | Self-interacting Scalar in AdS             |   3 |
| 3   | Holographic Strong Stochasticity Threshold |   5 |
| 3.1 | Multi-scale Analysis                       |   5 |
| 3.2 | Truncated Model                            |   6 |
| 3.3 | Full theory                                |   8 |
| 4   | Discussion                                 |  10 |

## 1 Introduction and Conclusion

In this paper, we will investigate the thermalization properties of a self-interacting probe scalar field in AdS. Our goal is to get some intuition about holographic thermalization in AdS gravity. We will find that for the questions we are interested in, the scalar probe has remarkable similarities to the full gravity system.

The set up we will consider in detail is as follows. We will launch the system with only a few modes excited, and then trace its evolution as the energy migrates to other modes. We find that the qualitative behavior of the system is drastically different above and below a certain threshold energy, E c .

We study this problem, first by restricting our attention to a finite-mode truncation of the system and keeping track of the (suitably time-averaged) energy per mode as the system evolves. We find that the energy per mode remains roughly constant when the system is kicked off (with a small number of modes excited) with E &lt; E T c . But when E &gt; E T c we find that there is fast thermalization, and that the energy in each mode in the truncated system approaches the same value. The transition from quasi-periodicity to fast thermalization is a very sharp one in the truncated case, even though the threshold energy and the thermalization times are dependent on the initial modes we start with.

We also investigate the self-interacting scalar probe, without truncating it to a finite number of modes. Note that now, since all infinite number of modes are present, the average energy per mode in the final thermal state should be zero if one starts with a finite energy. We find that a threshold energy, ∼ E U c , exists also in this un-truncated system, above which the thermalization time falls quickly by many orders of magnitude. The behavior of the full un-truncated system is less sharp

that the truncated case, reminiscent more of a crossover than a phase transition, but the phenomenon is again qualitatively very robust. In both the truncated and untruncated cases, we find that the thermalization timescales do not exhibit sharp changes as we further increase the energy above their respective E c .

The above paragraphs summarize our main results in this paper: when the initial modes are of small enough amplitude, quasi-periodic behavior persists and ergodicity does not set in for a long time, for the probe scalar in AdS. Now we turn to a brief summary of the relevance of these results as well as a discussion of some auxiliary observations. In the concluding section, we will suggest that a similar statement also holds for dynamical gravity in AdS.

Our primary motivation for studying this system is as an instructive toy model for holographic thermalization and for questions about the stability of AdS. AdS spacetimes with a reflecting boundary condition have been conjectured to be unstable at the non-linear level [1, 2]. It was argued that generic perturbations evolve into black holes after multiple reflections at the boundary [1]. This should not be thought of as an instability of the AdS ground state, but rather as the statement that generic excited states (even low-lying ones) eventually thermalize. This question is obviously of great physical interest not just from the viewpoint of understanding dynamical black holes and horizon formation (and the various related paradoxes), but also from the perspective of understanding thermalization in the dual strongly coupled field theory.

Even though this non-linear turbulent instability is believed to be generic, it has been noted recently by [3] that when only certain low-lying modes are excited, and the Einstein-(massless-)scalar system in AdS 1 is allowed to evolve, the system in fact does not thermalize for long periods. While the existence of the non-linear instabiltiy in the system for low amplitude coherent initial data is robust, the claims regarding the existance of the instability when only a few low-lying modes are excited are challenged by the results of [3]. They found instead that the evolution is quasiperiodic, and suggested that this behavior is analogous to the evolution of the famous Fermi-Pasta-Ulam-Tsingou (FPUT) lattice. One version of the FPUT paradox is the observation that when a self-interacting fluctuating chain with the Hamiltonian 2

$$H ( q, p ) = \sum _ { i = 1 } ^ { N } \frac { 1 } { 2 } p _ { i } ^ { 2 } + \frac { 1 } { 2 } ( q _ { i + 1 } - q _ { i } ) ^ { 2 } + \frac { \mu } { 4 } ( q _ { i + 1 } - q _ { i } ) ^ { 4 }$$

is allowed to evolve, the system does not thermalize for a very long time when only low lying modes are excited initially. Instead it exhibits quasi-periodic behavior. N here is the number of lattice sites.

1 This was the system where [1] had noted evidence for the non-linear instability of AdS. But, see also [4]. A very recent analytic approach to secular term resummation in the AdS instability problem can be found in [5].

2 This Hamiltonian is referred to as the FPUT β -system.

However, it is known in statistical mechanics that similar slow thermalization behavior is in fact not tied only to the FPUT lattice, but that other non-linear systems also exhibit similar behavior [6, 7] for low amplitudes. One such example is the quartic self-intreacting scalar on a lattice:

$$H ( q, p ) = \sum _ { i = 1 } ^ { N } \frac { 1 } { 2 } p _ { i } ^ { 2 } + \frac { 1 } { 2 } ( q _ { i + 1 } - q _ { i } ) ^ { 2 } + \frac { \mu } { 4 } q _ { i } ^ { 4 }$$

This is the massless version of a case considered in [6]. A bounded lattice with Dirichlet boundary conditions is analogous to the AdS box with reflecting boundary conditions, so one might wonder whether a self-interacting scalar field in the AdS geometry captures any of these features. As a by-product of our investigation, we will find indeed that it does: we find that the solutions are quasiperiodic for a long time, and we will also be able to reproduce the qualitative features of the relevant curves presented in [3] with our simple probe scalar system.

The conclusion of our work is that we find that when the initial modes are of small enough amplitude, quasi-periodic behavior persists and ergodicity does not set in . 3 These observations are in alignment with the idea that for many non-linear systems, there is a Strong Stochastisicty Threshold (SST) and thermalization happens only if there is enough energy (amplitude) in the modes. This remark has been made in various forms before in the statistical physics literature and some relevant papers are [6, 7, 9]. If the strong stochaticity threshold result holds not just for the probe scalar, but also for the backreacting gravity system, the observations made in [1] and [3] as well as the results of this paper can be seen to be consistent with this same general phenomenon, namely the existence of the SST in holographic systems; see the discussion in our section 3. Indeed, this can be seen as the holographic manifestation of the Strong Stochastic Threshold (HSST): for small enough energies in the low-lying modes, one has quasi-periodic behavior and for high enough energy the behavior changes qualitatively.

In the rest of the paper, we first introduce the systems under consideration and write down the equations of motion (section 2), describe our numerical set up and present the various numerical results that justify the claims in the abstract and introduction (section 3) and conclude with a discussion of the prospects, discussions and open questions (section 4).

## 2 Self-interacting Scalar in AdS

In this paper we will propagate a probe scalar field with the action

$$S = \int d ^ { 4 } x \sqrt { - g } \left ( \frac { 1 } { 2 } \nabla _ { \mu } \phi \nabla ^ { \mu } \phi + V ( \phi ) \right )$$

3 At least during the time period of the simulation! It is expected that given enough time these systems will [8] ultimately thermalize.

on a non-dynamical AdS, with reflecting boundary conditions. The dynamics of this system has previously been studied in [10], as a holographic model for quantum quench. We will take the potential to be that of a massless scalar m = 0 in what follows, but with a quartic coupling

$$V ( \phi ) = \frac { \lambda } { 4! } \phi ^ { 4 }.$$

We work with global AdS in the coordinates:

$$d s ^ { 2 } = \sec ^ { 2 } x ( - d t ^ { 2 } + d x ^ { 2 } + \sin ^ { 2 } x d \Omega ^ { 2 } ).$$

Specifically, we look for spherically symmetric solutions of the non-linear wave equation arising from the Lagrangian. The non-trivial (radial) part of the equation reduces to

$$\phi ^ { ( 2, 0 ) } + \Box _ { s } \phi \equiv \phi ^ { ( 2, 0 ) } - \phi ^ { ( 0, 2 ) } - \frac { 2 } { \sin x \cos x } \phi ^ { ( 0, 1 ) } = - \frac { \lambda } { 6 \cos ^ { 2 } ( x ) } \phi ^ { 3 } \, \quad \quad ( 2. 4 )$$

We will find it useful to have solutions of the eigenvalue equation glyph[square] s e j = ω e 2 j j for the spatial Laplacian glyph[square] s :

$$e _ { j } ( x ) = 4 \sqrt { \frac { ( j + 1 ) ( j + 2 ) } { \pi } } \cos ^ { 3 } ( x ) \, _ { 2 } F _ { 1 } \left ( - j, j + 3 ; \frac { 3 } { 2 } ; \sin ^ { 2 } ( x ) \right ), \\ \omega _ { j } ^ { 2 } = ( 2 j + 3 ) ^ { 2 }.$$

We will evolve the wave equation starting with various initial conditions, and we will be interested in keeping track of the total energy as well as the energy per mode. To avoid overall energy non-conservation and losses due to numerical artifacts, we will use a symplectic Hamiltonian formalism for dealing with the equations of motion numerically. The Hamiltonian density takes the form

$$\mathcal { H } = 4 \pi l ^ { 2 } \tan ^ { 2 } x \left ( \frac { l ^ { 2 } V ( \phi ) } { \cos ^ { 2 } ( x ) } + \frac { 1 } { 2 } \phi ^ { ( 0, 1 ) } ( t, x ) ^ { 2 } \right ) + \frac { \Pi ( t, x ) ^ { 2 } } { 8 \pi l ^ { 2 } \tan ^ { 2 } ( x ) },$$

where we have temporarily re-instated the AdS scale l for clarity, and Π is the canonical momentum. We have done the integration over the angles of the sphere in AdS already: in other words, this Hamiltonian density when integrated over the x gives us the full Hamiltonian of the system. In terms of the mode sum, this takes the form

$$H = \frac { 1 } { 2 } \sum _ { i = 0 } ^ { \infty } \left ( \dot { A } _ { i } ( t ) ^ { 2 } + \omega _ { i } ^ { 2 } A _ { i } ( t ) ^ { 2 } \right ) + \sum _ { i, j, k, l = 0 } ^ { \infty } \mathcal { C } _ { i j k l } A _ { i } ( t ) A _ { j } ( t ) A _ { k } ( t ) A _ { l } ( t ) \quad \ \ ( 2. 8 )$$

where the C ijkl coefficients can be straightforwardly computed.

## 3 Holographic Strong Stochasticity Threshold

## 3.1 Multi-scale Analysis

At low energy (i.e. A 0 glyph[lessmuch] 1) resonant transfer between modes is the most important mechanism of energy transfer [9]. The naive perturbation theory in λ breaks down at larget due to resonances. However in often-used multi-time analysis, one separates the slow and fast parts of the degrees of freedom by defining A t i ( ) ≈ B τ e i ( ) iωt + ¯ B τ e i ( ) -iωt , where τ ≡ glyph[epsilon1] 2 t is the slow-moving time. If we take a t -average of the Lagrangian over a fundamental period and assume τ dependent terms to be constant. We get,

$$L _ { M S } = i \sum \omega _ { i } ( B _ { i } \dot { \bar { B } } _ { i } - B _ { i } \dot { \bar { B } } _ { i } ) + \sum \mathcal { C } _ { i j k l } B _ { i } ( \tau ) B _ { j } ( \tau ) B _ { k } ( \tau ) B _ { l } ( \tau ), \quad ( 3. 1 )$$

with ω i + ω j -ω k -ω l = 0 in the second term. Here, in the interaction part only t -independent terms satisfying the resonance condition survive. One can systematically include corrections by introducing other time scales in (3.1), but even the zero'th order theory gives us information about how energy is transferred across resonances.

One interesting observation is that other than the usual time translation invariance, (3.1) has a set of enhanced symmetries. From the dilatation symmetry of the form B k ( τ ) → αB k ( 1 α 2 τ ), one may naively presume that the thermalization time goes inversely with the square of initial amplitude. However (3.1) has some other enhanced symmetries that are not possessed by the full system, so it seems possible that these other symmetries might forbid the thermalization of the theory in the multi-scale approximation. Here we list the conserved charges and the corresponding symmetries:

$$Q _ { 0 } = \sum B _ { k } \bar { B } _ { k }, \ \text{ Symmetry} \, \colon B _ { k } \to e ^ { i \theta } B _ { k },$$

$$Q _ { 1 } = \sum k B _ { k } \bar { B } _ { k }, \ \text{ Symmetry} \, \colon B _ { k } \to e ^ { i k \theta } B _ { k },$$

$$E = \sum _ { \omega _ { i } + \omega _ { j } - \omega _ { k } - \omega _ { l } = 0 } \mathcal { C } _ { i j k l } B _ { i } \bar { ( } \tau ) B _ { j } \bar { ( } \tau ) B _ { k } ( \tau ) B _ { l } ( \tau ), \ \text{ Symmetry} \, \colon t \to t + \alpha. \quad \ \ ( 3. 4 )$$

A closely related fact is that the multi-scale theory has a set of quasi-periodic solutions given by B k ( ) = t α k exp( -i ( β t 0 +( i -1) β t 1 )). Plugging this into the multiscale equation of motion we can solve for β , β 0 1 , α k , as also observed in [3]. The enhanced symmetries and the accompanying slow thermalization are consistent with the previous observation [11] that the thermalization times are of the order of e -1 /A δ in some statistical systems, where δ is some positive constant.

To provide evidence that our truncated system can capture the relevant features of the gravity system, we first note that the scalar evolution is quasi-periodic for a long time, as already emphasized. This was the case in [3] as well. We also present plots of some of the quantities. These have qualitative similarities with the results in [3]. In particular our figures 1(a), 1(b) (which were made using a Two-Time Formalism), should be compared to figures 3 and 5 in [3].

GLYPH&lt;160&gt;

(a) The upper envelope of ˙ ( φ x = 0) , is qualita2 tively identical to figure 3 in [3]. The fine structure here in the envelope is due to a finer smoothing scheme.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0624v1-figure-3f5239ee-13ee-4d4e-887b-f82c33f7487d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0624v1-figure-91c735d0-4fcb-47f5-bafd-a05ae967a263.jpg)

(b) Running time-average energy per mode ¯ E J vs t for equal-energy two-mode initial data. This is qualitatively identical to figure 5 in [3]. The curves correspond to j = 0 1 2 3 top to bottom. , , ,

Figure 1 : Both these plots have been made with initial data that contained equal energy in the j = 0 1 modes. ,

## 3.2 Truncated Model

As a preliminary exercise, we will first study a system where we truncate the full system and restrict the modes to only the first 9 energy eigenstates: ie., the sum in (2.8) will be from 0 to 8 . 4 We integrate the resulting Hamilton's equations of motion with an initial condition A 1 (0) = A 0 while setting the rest of the initial A i 's and all the conjugate momenta to zero. To do this, without incurring fictitious losses in the total energy due to numerics, we resort to a sympletic integration of the equations that preserves the phase space volume.

Due to the presence of non-linearity, the dynamics of the system is not exactly integrable. However, this fact does not tell us much about the approach of the system towards thermalization as it evolves. In fact, thermalization is a somewhat tricky idea to precisely capture, so we will use the following approach to define it. A thermal system shows equipartition of energy among modes, i.e. the virial theorem holds. For the kinetic part it means that the time average

$$\langle \dot { A } _ { i } ( t ) ^ { 2 } \rangle _ { t } \sim 1 / t \int _ { 0 } ^ { t } \dot { A } _ { i } ( \tau ) ^ { 2 } d \tau$$

for larget is independent of the mode number. So we will use this time-avegared kinetic energy per mode as a quantity that keeps track of the approach to thermalization.

We have investigated the full time evolution of the truncated theory. When the initial amplitude is small, within the scope of our numerics we do not see any tendency of the theory to equilibrate even after long times of integration. In fact, we

4 We also study the case where we restrict to 17 modes, but the results are qualitatively similar.

Figure 2 : Plots of time average of kinetic energies (in the truncated model) for modes- 0 1 2 3 from the top. , , ,

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0624v1-figure-527599f9-6ce8-43d1-9253-0a9f490a6086.jpg)

always find that the average kinetic energy of high enough modes are exponentially suppressed, i.e. log( 〈 ˙ A n ( ) t 2 〉 t ) ∼ -O n ( ) for largen (see Fig 3).

Figure 3 : Mode number vs lograthim of 〈 ˙ A n ( ) t 2 〉 t for large t .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0624v1-figure-42ba1fec-4942-493e-9f58-e4fb7a8ed416.jpg)

On the other hand, according to the theory of [9] (see also [6]) one expects that at higher energy direct energy transfer between modes (outside resonant couplings) would be important. For example in our case, we find that for mode j = 0, if we start with an initial value of A &gt; 0 0 707(approx), the theory quickly thermalizes. We . have considered various other modes and combinations of modes as well: a similar story (but with other values of the threshold amplitude/energy) holds also in these cases. This is evidence for a strong stochasticity threshold (SST).

In the mode-truncated system that we discussed in this subsection, the SST is a very sharp feature and the transition is quite drastic (almost as if it were a phase transition). When we keep all the modes and do not do the truncation, we will again see that there is a big qualitative difference in thermalization times as we cross the threshold energy, but the transition is less sharp and more reminiscent of a crossover rather than a phase transition. This is what we discuss in the next subsection.

Figure 4 : Plots of time averages of kinetic energy with time, for modes- 6(top) and 7(bottom).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0624v1-figure-30a84a1c-81bc-41d7-a10a-d59202087964.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0624v1-figure-bbfcf581-5736-4310-973c-aec57d23e1ac.jpg)

## 3.3 Full theory

Here we look at the time evolution of full scalar equations in AdS numerically. To understand the time evolution numerically, we discretize the hamiltonian on a lattice with N = 256 ( and also 512 points) and time evolve the Hamilton's EOM's using a simplectic Runge-Kutta algorithm.

We start with various initial conditions by turning on various combination of modes. We will confine ourselves to low lying modes and highest mode we would go up to is ten. Keeping the ratio of amplitudes of modes fixed, we change the overall amplitude. We always see a SST, below which the theory does not have effective thermalization and over which the theory quickly thermalizes.

We here describe the main story with only the eighth mode turned on initially. We look at the time evolution of the system for a time scale of the order t = 1 × 10 8 (to compare ω 0 = 3 and T 0 = 2 π ω 0 ≈ 2 094). . Technically thermalization time may be tricky to define. Here we would not go into the technicalities. For our purpose we will pick two modes (taken to be 7 , 8 in the plots) and look at their average kinetic energy (averaging started at t ≈ 2 75 . × 10 6 ) to have a better idea of thermalization. For a

Kinavg

Figure 5 : Plots of time average kinetic energy with time, for the modes- 0, 2, 4, and 6 (from the top at left).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0624v1-figure-f6581f61-41f2-4ddc-bfb1-7306d0d2009b.jpg)

small initial amplitude the theory does not show a tendency to thermalize within the integration time scale Fig 4(a). The situation changes after quickly after A 8 (0) ≈ 1 08 . (Figure-4(c)). Between Fig-4(a) and Fig-4(d), we clearly see a crossover behavior and a few percent increase of initial energy leads to a huge decrease in thermalization time.

It is also interesting to note how thermalization happens. For that we made detailed plots of average kinetic energy of many modes. For low enough initial amplitude, initially for some time the system is trapped in a non-thermal steady state. For example, in Fig-5(b), before the average kinetic energies thermilizes, they level out like in Fig-5(a). A vestige of this phenomenon can also be seen in Fig-5(c). One possible explanation is that system was trapped in valley in the phase space with a narrow escape, and when the phase space point reaches that escape location, the kinetic energy shoots up. This phenomenon is not observed for A 8 (0) &gt; 1 3 (see . Fig- 5(d)). Also thermalization behavior of the system is only weakly dependent on the initial energy after the thermalization threshold. It is believed that the dominant

channel of energy transfer between modes changes from resonance transfers to nonresonant transfer after the thermalization threshold [6].

Another noteworthy feature of the curves is the upward rise in the time-averaged kinetic energy per mode close to the onset of thermalization. We do not believe that this is a numerical artifact , but we also do not understand the physical significance 5 of this. But it is tempting to speculate that this is due to an escape from a trapped region in phase space as we mentioned in the previous paragraph.

To capture some of the features of the threshold, we also show the spatial profiles of the scalar field after a large number of time steps, above and below the SST: the drastic transition of the waveform from quasi-periodic to chaotic is evident (see Fig 6), as expected of a system that is close to thermalization. It is also noteworthy that unlike in the case of the mode-truncated system, after thermalization, the system relaxes to zero average energy per mode as it should. We have checked in specific cases that the average energy per mode comes down to very low values (1% of original). But the computational time required to evolve the system for long times is large, so we have not checked it for all cases.

## 4 Discussion

In this paper, we have demonstrated the existance of a Holographic Strong Styochasticity Threshold for thermalization of an interacting probe scalar field. We conjecture that a similar threshold exists also in dynamical gravity in AdS. Since it is only in the limit of small 6 glyph[epsilon1] that the TTF analysis of [3] is expected to hold, and since this corresponds precisely to the case where the amplitudes are forced to be small [3], and since the instability identified by Bizon and Rostworoswki is expcted to be generic, we believe that an SST-like mechanism analogous to the one we found for the probe scalar should operate also for gravity. In a thought experiment we may consider turning on gravity at a late time: if enough energy is transferred to localized modes in the position space, even if solely by scalar self-interaction, then eventually a black hole would form (see e.g. [12] for connection of classical chaos and black hole formation).

Our discussion opens up various intriguing questions, and we conclude this paper by mentioning two of them, which are currently under investigation.

- · Is the total energy in the modes that is the relevant quntity for deciding the onset of thermalization or is it the maximum energy in any one mode? More pertinently, does the number of modes that one turns on affect the threshold

5 The timescale of the rise is a fraction of the total time-range, and is very many orders of magnitude bigger than the integration time-step. Besides, we have checked that the total energy in the system in our sympletic integrator is constant to within negligible numerical errors.

6 This glyph[epsilon1] is the (square-root) of the ratio of the two 'times' in the TTF formalism, see [3].

Figure 6 : Radial profile of the scalar field at t ≈ 6 283 . × 10 . 5

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0624v1-figure-a3825187-a892-4efa-9379-f422f3f9358c.jpg)

at which thermalization sets in? How do things change when we slowly start increasing the number of modes. The results of [1] suggest that for smooth profiles, the instability is fairly generic. It will be interesting to understand the connection between the existance of this energy threshold and the threshold of black hole formation, a la Choptuik. Naively one might have thought that Coptuik's result is the holographic analog of the SST, but in AdS for coherent initial data, the observation of [1] is that black holes form for arbitrarily low amplitudes.

- · Is there something that can be said about the system analytically? Work along this lines is currently under progress. See [13].

## Acknowledgments

PB thanks Abhishek Dhar, Sumit Das and Oscar Dias for discussions; and YITP, Kyoto for hospitality during part of this work. CK thanks Ben Craps, Alex Buchel,

Umut Gursoy, Sachin Vaidya and Leopoldo Pando Zayas for discussions/correspondence; and the Spinoza Institute (Utrecht) and ETH (Zurich) for hospitality during part of this work.

## References

- [1] P. Bizon and A. Rostworowski, On weakly turbulent instability of anti-de Sitter space , Phys.Rev.Lett. 107 (2011) 031102, [ arXiv:1104.3702 ].
- [2] P. Bizon, Is AdS stable? , Gen.Rel.Grav. 46 (2014) 1724, [ arXiv:1312.5544 ].
- [3] V. Balasubramanian, A. Buchel, S. R. Green, L. Lehner, and S. L. Liebling, Holographic Thermalization, stability of AdS, and the Fermi-Pasta-Ulam-Tsingou paradox , arXiv:1403.6471 .
- [4] O. J. Dias, G. T. Horowitz, D. Marolf, and J. E. Santos, On the Nonlinear Stability of Asymptotically Anti-de Sitter Solutions , Class.Quant.Grav. 29 (2012) 235019, [ arXiv:1208.5772 ].
- [5] B. Craps, O. Evnin, and J. Vanhoof, Renormalization group, secular term resummation and AdS (in)stability , arXiv:1407.6273 .
- [6] M. Pettini and M. Cerruti-Sola, Strong stochasticity threshold in nonlinear large hamiltonian systems: Effect on mixing times , Phys. Rev. A 44 (Jul, 1991) 975-987.
- [7] K. Yoshimura, Strong stochasticity threshold in some anharmonic lattices , Phys. D 104 (May, 1997) 148-162.
- [8] G. Gallavotti, The Fermi-Pasta-Ulam Problem: A Status Report , .
- [9] F. Izrailev and B. Chirikov, Statistical properties of a nonlinear string , Sov. Phys. Dokl 11 (1966), no. 1 30-32.
- [10] P. Basu, D. Das, S. R. Das, and T. Nishioka, Quantum Quench Across a Zero Temperature Holographic Superfluid Transition , JHEP 1303 (2013) 146, [ arXiv:1211.7076 ].
- [11] M. Pettini and M. Landolfi, Relaxation properties and ergodicity breaking in nonlinear hamiltonian dynamics , Phys. Rev. A 41 (Jan, 1990) 768-783.
- [12] L. A. Pando Zayas, A Quantum Rosetta Stone for the Information Paradox , arXiv:1405.3655 .
- [13] A. Ponno, L. Galgani, and F. Guerra, Analytical estimate of stochasticity thresholds in Fermi-Pasta-Ulam and φ 4 models , Phys.Rev. E61 (2000) 7081-7086.
Pallab Basu Chethan Krishnan Ayush Saurabh
A Stochasticity Threshold in Holography and and the Instability of AdS
We give strong numerical evidence that a self-interacting probe scalar field
in AdS, with only a few modes turned on initially, will undergo fast
thermalization only if it is above a certain energetic threshold. Below the
threshold the energy stays close to constant in a few modes for a very long
time instead of cascading quickly. This indicates the existance of a Strong
Stochasticity Threshold (SST) in holography. The idea of SST is familiar from
certain statistical mechanical systems, and we suggest that it exists also in
AdS gravity. This would naturally reconcile the generic non-linear instability
of AdS observed by Bizon and Rostworowski, with the
Fermi-Pasta-Ulam-Tsingou-like quasi-periodocity noticed recently for some
classes of initial conditions. We show that our simple set-up captures many of
the relevant features of the full gravity-scalar system.