Vacuum block thermalization in semi-classical 2d CFT

The universal nature of black hole collapse in asymptotically $AdS_3$ gravitational theories suggests that its holographic dual process, thermalization, should similarly be fixed by the universal features of 2d CFT with large central charge $c$. It is known that non-equilibrium states with scaling dimensions of order $c$ can be sorted into states that eventually thermalize and those that fail to do so. By proving an equivalence between bounded Virasoro coadjoint orbits and certain (in)stability intervals of Hill's equation it is shown that a state that fails to thermalize can be promoted to a thermalizing state by preparing the system beforehand with an energy greater than an appropriate threshold energy. It is generally a difficult problem to ascertain whether a state will thermalize or not. As partial progress to this problem a set of lower bounds are presented for the treshold energy, which can alternatively be interpreted as criteria for thermalization.


Introduction
Strongly coupled quantum field theories provide examples of very interesting but extremely complex chaotic phenomena. Unfortunately the amount of precise statements that can be made about systems at strong coupling is very limited, and these typically only apply to very specific models. One family of strongly coupled theories where one can make rather powerful statements is 2d conformal field theory. In particular, if one narrows it down further and considers 2d CFT in the regime of large central charge, in which case the amount of accessible information grows even further.
Better yet a lot of these statements are often universal, depending only on conformal symmetry, modular invariance or the presence of the identity operator in the operator spectrum, while otherwise completely independent of the microscopic details of the theory. Examples include the famous Cardy formula [8], [9] and the universal ground state entanglement entropy [10] [11]. A particularly fruitful common feature of 2d CFT is the universal contribution of the vacuum representation to the Hilbert space and its state-operator dual, the identity operator. In the 't Hooft limit in higher dimensions the presence of the identity operator block ensures large-N factorization, in 2d dimensions the extended vacuum representation under the much larger Virasoro algebra causes the identity block to capture substantially more non-trivial information of the field interactions, in particular this has been exploited in for instance [3] [12].
Of course one of the most, if not the most, exciting developments in theoretical physics of the last decades is the AdS/CFT correspondence. Its most pressing contributions are in the form of implications to quantum gravity. Rightfully so, most of the AdS/CFT literature is devoted to formulating (or possibly defining) quantum gravity in terms of conformal field theory. But the flipside of the conjecture provides us with a lot of intuition on strongly coupled conformal field theory with large numbers of degrees of freedom in terms of semiclassical gravity. For one we expect the gravitational sector of the AdS 3 bulk dual to exhibit a phase transition in the form of black hole collapse. In the boundary CFT we hence expect that non-equilibrium states evolve at asymptotically late Lorentzian times to either evolve to an equilibrium state that is either in an integrable phase or a thermal phase.
The key word in the CFT description of black hole collapse is thermalization, it is typically said that a state thermalizes if at asymptotic late times a certain class of sufficiently 'simple' correlators on the state approach the expectation values that would have been obtained if they had been computed on a thermal state instead. The subject of this paper will be to consider the identity block contribution to a two-point function of light probe operators on a heavy non-equilibrium state, specifically at leading order in a large central charge expansion.

Summary of the central message
The distinguishing function determining whether a state thermalizes or not is the stress tensor expectation value of the state. A family of stress tensor expectation values which are connected to each other through conformal transformations form a Virasoro coadjoint orbit, all of these distinct orbits have been classified [5,6,13,32]. While the role of these orbits within the context of 3d gravity with a negative cosmological constant is known [7,[15][16][17][18][19][20]29], in [2] the authors established from the conformal field theory perspective that the orbit in which the stress tensor expectation value of a state is contained determines its ultimate fate.
In this paper it will be demonstrated that given a non-equilibrium state that would keep the system in an integrable phase, there exists a state associated energy scale such that one can trigger thermalization. Specifically if we prepare the system in advance in an energy eigenstate with an energy above this scale and release the non-equilibrium state on top of this eigenstate the system will thermalize. Furthermore it will be shown that unitarity combined with Virasoro representation theory ensures that there can only be at most one such transition scale per state.
While the gravitational bulk picture of this phenomena is very intuitive (see figure 1) a pure conformal field theory description of this transition for CFTs with large numbers of Figure 1: The bulk gravitational picture is particularly intuitive. If a non-equilibrium collapse state is prepared that will collapse into a non-thermal geometry (e.g. a conical defect geometry), then if the experiment were to be repeated with an initial conical defect with mass m in the asymptotic past, then there exists a critical value of m such that the conical defect is pushed over the critical opening angle and a BTZ black hole is formed instead.
degrees of freedom appears to be lacking in the literature, this paper attempts to fill that gap.
on a more pragmatic note, this approach naturally presents a set of lower bounds on the transition energy in terms of the other Virasoro conserved charges of the state. These provide us with a set of necessary (but not sufficient) conditions that a state needs to satisfy in order to thermalize. It is an open problem to find an analytically accessible diagnostic deciding whether a state thermalizes or not [20], in principal the answer is given by a path-ordered integral but it is generally not feasible for these to be analytically computed. The necessary conditions in terms of Virasoro charges gives a partial answer to this question in terms of much simpler contour integrals instead. Alternatively, if the Virasoro conserved charges are interpreted as thermodynamic equation of state variables, then the derived lower bounds can be thought of as bounding the equation of state in the thermalizing phase. This paper has been organized in the following way, section 2 provides a relatively lowtechnical review of the identity block contribution to probe two-point functions on heavy states and the role of Virasoro coadjoint orbits. Section 3 discusses some of the general features of Hill's equation, the class of ODEs that naturally occur in these two-point functions and presents an equivalence between a subspace of Hill's equations and the space of bounded Virasoro orbits. In section 4 a set of lower bounds are derived for the transition energy from an elliptic to hyperbolic state. Furthermore two appendices have been added, the first of which contains a more technical review of the classification of Virasoro coadjoint orbits. The second appendix provides some of the more technical aspects of the proof for the argued equivalence in section 3. For a reader who is interested in an overview of the results but wishes to avoid the technical derivations reading section 2 is in principal sufficient.

Black hole collapse and the Virasoro identity block
One particular feature of asymptotically AdS 3 gravity is that its black hole solutions, the BTZ black holes, possess a mass gap. As a result not all collapsing states can form black holes. Naively one could expect that the total energy of a state could act as a diagnostic to determine whether a collapse process' final state consists of a black hole. Interestingly this does not quite work due to the rich boundary dynamics attributed to the Brown-Henneaux conserved charges [30]. Simply put, arbitrarily large amounts of energy can be stored in the Virasoro boundary charges without contributing to a potential final state black hole mass, in fact one could remove or increase this energy by means of a conformal transformation [20]. As an analogy, one can consider a bound system of particles, one could increase the energy of this system by arbitrarily large amounts by boosting it to higher and higher velocities but it will not affect the Schwarzschild radius of the bound system.
By the holographic nature of asymptotic AdS 3 gravity it is expected that there exists a CFT analogue to this dichotomy between states |V that collapse to a black hole and states that fail to form black hole states. The holographic dictionary states that the analogue would be that some initial states will equilibriate to states that are indistinguishable from being thermal [21] while other states remain in an integrable phase. The class of 2d CFTs under consideration will be very broad, mirroring the universality of black hole thermodynamics in gravitational theories. The focus will be on 2d CFTs with very large central charge c, due to the semi-classical nature of black hole thermodynamics but mostly due to immense analytical simplification it provides. Secondly we will assume a sparse operator spectrum at low scaling dimensions, this allows one to find a kinematical regime where the dominant contribution to correlation function is due to the vacuum Virasoro block, as far as we know there exist no examples of theories with gravitational duals that violate this assumption.
Similarly we will want to examine the late time behaviour of a set of states that is as broad as possible. One area where the large c approximation particularly shines is the calculation of correlation functions of operators with scaling dimensions that are small with respect to the central charge. Unfortunately in order to overcome the Casimir energy of the Lorentzian cylinder |E C | = c/24 one needs to consider states |V with scaling dimensions H V proportional to c (typically called heavy states). Put explicitly we will want to consider states of the generic form with heavy operators O ik and possibly some sufficiently smooth smearing function g(z). The insertion locations z ik are not completely free since we will want to compute correlators involving this state where there exists an OPE channel where the identity block contribution dominates. One way to achieve this is to restrain z ik to a thin annulus whose exterior boundary is given by the unit circle. The heaviness of the state provides an additional complication that was tackled by the authors of [3], [4] in what they call the uniformization approach. They applied there strategy to states generated by acting on the origin with a primary operator, in [2] the authors applied this method to a wider class of heavy states, and showed that while the resulting differential equations typically can not be solved analytically some qualitative features can be derived. It was for instance shown that if the heavy state thermalizes at late time that the temperature matches the temperature of the bulk dual black hole in the Chern-Simons formulation of AdS 3 gravity. Some of the results of [2] are briefly covered in section 2.1.
It will be shown that given an integrable heavy state increasing the energy of this state in a controlled way will cause it to cross a threshold that will trigger thermalization at late Lorentzian time. A state that thermalizes at late times will be referred to as a hyperbolic state, conversely states that fail to thermalize will be called elliptic states, this terminology makes particular sense in the following discussion but conveniently (though not coincidentally) this convention was already standardized in [16].
We will see that the existence of a single such threshold is fixed by a combination of representation theory of the Virasoro algebra and classic properties of linear ordinary differential equations.

probe two-point function
To probe the behavior of the final state of an out-of-equilibrium state |V consider the latetime limit of a two-point function of light operators on the state the coordinates z are on the radial plane. The operator Q h is light in the sense that h c as such the amount by which they affect the dynamics of |V surmounts to a 1/c correction, for this reason they form appropriate probes to measure physics of the state |V , in order for the the correlator to be radially ordered we restrict the probe coordinates z 1 , z 2 to the unit circle. Note that while this correlator is radially ordered on the Euclidean radial plane the Lorentzian correlator of interest is out of time order, this has to be taken into account when analytically continuing to the Lorentzian cylinder. The direct interest of this note is the late Lorentzian time limit, in [2] it was shown that after mapping the radial plane to the Euclidean cylinder z = exp(τ + iφ) and analytically continuing the Euclidean time coordinate τ time-evolving in Lorentzian time corresponds to letting the probe observers undergo counter-clockwise circles on the radial plane. The non-trivial multi-valuedness of the correlator A(z 1 , z 2 ) ensures that the result is not necessarily periodic in time.
The correlator will be simplified in two directions, first we will only look at the leading order in a 1/c expansion secondly we will only consider the identity block exchange, i.e. we only consider the exchange of identity operators when performing OPE contractions between the operators of the state (2.1) and the probe operators. Unlike in higher dimensional CFT where the vacuum block contribution merely computes the factorized correlator in 2d CFT the vacuum block contribution to a correlation function captures highly non-trivial information between the interaction of the state onto the probes. One of the Virasoro descendants of the identity operator is the stress tensor, the AdS/CFT dictionary states that its dual field is the graviton field, therefore the bulk interpretation is that Virasoro identity block resums all graviton exchanges between the probes and the geometry generated by the heavy state.
It was found through means closely related to the uniformization of punctured Riemann surfaces [3][23] that the correlator A(z 1 , z 2 ) could be written in the form where ψ(z) satisfies the Fuchsian ODE and where the expectation value T (z) is defined through The freedom in picking a particular solution ψ(z) of (2.4) the reflects the freedom in conjoining a uniformizing coordinate system with a global conformal transformation without spoiling its uniformization property. As a complex function T (z) cannot be entirely generic, due to the fact that the in-and out-states are identical up to Hermitian conjugation, taking the complex conjugate of (2.5) shows that it is subject to the following Schwarz reflection property this for one implies that z 2 T (z) will be real-valued along the unit circle which will be relevant when we consider the associated Hill system in section 3. Schwarz reflection also restrict the monodromy matrix around the unit circle M ψ of a basis of solutions of (2.4). There has to exist a basis of solutions such that M ψ is contained within SU (1, 1) [22], which restricts the eigenvalues of M ψ in any basis to consist of both purely real numbers or both pure phases. Additionally the determinant constraint fixes the eigenvalues to be each others inverse (µ, 1/µ) The late Lorentzian time limit corresponds to letting both the coordinates z 1 and z 2 perform a large number of cycles around the unit circle. If the eigenvalues are given by real numbers (µ, 1/µ) the probe correlator is proportional to identifying n i with a Lorentzian time through t i = 2πn i we find that (2.7) resembles a twopoint function evaluated on a thermal state with a temperature T V given by This concludes a condensed version of some of the results in [2].

Temperature as a conformal invariant and Virasoro orbits
The arguments reviewed above demonstrate that whether a heavy state |V thermalizes and what temperature it thermalizes to is controlled by the eigenvalues of the monodromy matrix M ψ of the independent solutions of the linear ODE (2.4). It can be shown that these eigenvalues are independent of conformal frame, consider the single-valued holomorphic mapping w(z), under this transformation the stress tensor T (z) transforms as it can easily be verified that if ψ(z) solves the ODE (2.4) that after the substitution (2.9) the ODE (2.4) will be solved by Therefore the solution φ(w) under the curve w(e iθ ) inherits the monodromy of the solution ψ(z) under traversing the unit circle 1 . Hence the matrix M ψ forms a conformal invariant of the stress tensor. A set of functions T (z) that are connected to each other through the action (2.9) is given by a Virasoro coadjoint orbit, see appendix A for a short review based on [5], [6], [31], for an additional review see [7]. All disconnected Virasoro coadjoint orbits have been classified, for instance the set of all hyperbolic orbits is designated by B n (b). Here the index n indicates the number of zeroes of the independent solutions of (2.4) in the diagonal basis, the argument b uniquely parametrizes a certain representative element of the orbit. While all the various Virasoro coadjoint orbits have been classified (see table 1) not all of them are physically relevant. If we define the holomorphic contribution to the energy the usual way as the zeromode of the stress tensor it is known that all orbits contain elements with arbitrarily large energies. This is not a problem but it does complicate the process of sorting T (z) into an orbit as it almost completely disconnects the temperature and the energy as quantities. What is problematic is that some orbits possess an energy that is unbounded from below, as one of the assumptions of unitary conformal field theory is that the vacuum state is the lowest energy state in any conformal frame this leads to a contradiction, hence we can dismiss these orbits as being unphysical. Each Virasoro orbit that possesses a lower bound on its energy also possesses a primary representative, meaning that for any T (z) contained within that orbit there exists a singlevalued transformation that takes it to a stress tensor of the form T (z) = H/z 2 . This is exactly the form of the stress tensor expectation value of a state created by a primary operator 1 There is a small subtlety present concerning winding number, one could construct a map w(z) such that a closed curve in w-space has to be traversed multiple times before z(w) performs a single S 1 circle in z-space.
In [5], [6] this issue is bypassed by considering only mappings S 1 → S 1 with winding number 0, in this paper the issue will largely be disregarded.

Orbit
Stabilizer SL(2, R) class Bounded energy Constant representative Table 1: Summarized form of the various features of the Virasoro orbits as classified in [5] and reviewed in appendix A.
with scaling weight H inserted at the origin, O H (0)|0 . They are related to the constant representatives of [6], [5] through the coordinate transformation z = e iθ , in both cases they are the representatives of an orbit whose energy saturates the lower bound. The fact that each physical orbit contains such a primary representatives has some implications with respect to the eigenstate thermalization hypothesis [21], it implies that any generic heavy state is sharply dominated by a single conformal family, see [2] for a short discussion on this point. The set of distinct physical orbits are dominated by the hyperbolic B 0 (b) orbits and the elliptic C ν orbits with 0 < ν < 1. Besides these the physical spectrum consists of one parabolic orbit P + 0 and an exceptional orbit E 1 respectively corresponding to the BTZ cross-over point and the vacuum. Each of these orbits possesses a primary representatives of the form From the AdS/CFT dictionary mass ←→ scaling dimension it can be seen that it is consistent to attribute to a hyperbolic state a bulk interpretation as being conformally related to a BTZ black hole in the center of AdS. Similarly we can think of elliptic orbits as being conformally equivalent to sourcing a conical defect geometry and of the parabolic orbit as the minimal mass BTZ geometry. As these are the lowest energy representatives they solidify the CFT intuition that |V can be thought of as being dominated by a single descendant state of Virasoro representation whose highest weight state has a scaling dimension given by the numerical factor in front of the primary representative stress tensor.
One would rightly think that it would be very difficult to couple a generic state |V to the appropriate primary representative since they are related to each other by the solutions of extremely complex differential equations [5]. Fortunately there exists a large body of known theorems about ODEs closely related to (2.4). These theorems will allow us to prove that we can change the primary representative of a state by acting on a state with a primary operator at the origin, |V H = O H (0)|V . Given a set of stress tensor expectation values we can think of this family as a curve through the space of all possible stress tensor expectation values parametrized by H. In section 3 it will be proven that this curve intersects a bounded Virasoro coadjoint orbit only once, by showing that the constant representative of the orbit increases monotonically as a function of H. In fact this tells us that given an elliptic state |V there exists a certain state-specific critical scaling dimension H * such that |V H * falls within the class of hyperbolic states. The bulk interpretation is that without any matter in the center of AdS an elliptic orbit T (z) would source a geometry that would evolve to a non-thermal geometry (e.g. a conical defect), but if there was a point source of matter in the center of AdS before releasing |V then the total scattering mass might be sufficient to trigger black hole collapse.
Finding the exact value for the energy neccessary to trigger such a transition is a very difficult problem but it is possible to determine a set of bounds on the transition energy in terms of the Virasoro conserved charges of the stress tensor. In section 4 it will be derived that for any integer n = 0 the transition value of L 0 designated L |V is bounded from below by From this one can easily interpret this instead as an infinite set of neccessary conditions for a stress tensor T (z) to be contained within a hyperbolic orbit to emphasize the meaning of these inequalities, all of these need to hold in order for a state |V associated to T (z) to thermalize, if just one of them fails it implies that T (z) is contained within an elliptic orbit and hence |V fails to thermalize.

From Fuchs to Hill's equation and its implications
In order to determine the late-time fate of a state |V one needs to know the monodromy of the solutions to (2.4) around the unit circle. This suggests that for the purpose of determining the monodromy, or (as we shall see) equivalently the orbit of T (z), knowing the solutions on the full complex plane is largely superfluous. Therefore it makes sense to restrict ODE to the unit circle, as a result one obtains a real-valued ODE of Hill class. The class of Hill's equations are of the following general form where the coefficient function Q(x) is a periodic function satisfying There exists a very substantial library of known theorems on the qualitative behavior of this large class of ODEs [1] which will be exploited in the upcoming sections. Turning our attention back to the Fuchs equation where now T (z) is complex meromorphic function subject to the reflection constraint as mentioned before this directly implies that the product z 2 T (z)| |z|=1 is real-valued. To rewrite this ODE into Hill form first map the unit circle to the interval [0, π] by means of z = e 2ix , in which case the ODE transforms to This can (almost) be brought into the canonical form above by means of the transformation In order to fully connect to the canonical form (3.1) the zero-mode needs to be subtracted out of the coefficient function Hence we obtain the Hill system (3.1) whose elements are described in table 2.
There are a few interesting things here to note. First that this is exclusively a well-defined Hill system if the coefficient function is guaranteed to be real for all real values of x. This in turn is guaranteed by the reflection symmetry of the stress-tensor expectation value T (z) Secondly the number λ which plays a very important role in various stability theorems of the Hill equation, has a clear physical interpretation, being for one related to the holomorphic Euclidean energy L 0 . Additionally if the total energy of the state is converted into a black hole (i.e. no boundary gravitons) λ satisfies |λ| ∝ T BH [2] [3]. In the next section we will further expand on the CFT role of λ. Table 2: The constituents of the Hill system in terms of CFT data.

The off-set parameter as a CFT quantity
In a large amount of the theorems we will quote a deciding role will be played by the off-set parameter λ. Within the CFT context the parameter λ is not so esoteric as will be shown. Consider a state |V and its stress tensor expectation value T (z) From the state |V we can construct a family of states by acting with primary operators whose holomorphic stress tensor expectation value is given by Both the Fuchs and Hill equation take as input the expectation value of the stress-tensor on the Euclidean radial plane, hence in this specific context no analytic continuation to Lorentzian time is required, as a result the apparent breakdown of radial ordering is of no concern. While no general proof is known it is generally believed that in the large c domain correlators take the general form e − c 6 f (z i ) [10,[24][25][26], this combined with the holomorphic Ward identities fixes T H (z) to the form demanding that both T (z) and T H (z) satisfy the reflection symmetry (2.6) sets the coefficient This demonstrates that acting on a state with a primary operator has the effect of shifting the L 0 mode of the stress tensor without affecting the other modes. From the table 2 we can as a result read off that acting with a primary operator on the state has the effect of changing the relevant Hill equation to one with the exact same potential Q(x) but with a lower off-set parameter λ.
Acting with a primary operator on the state has a clear bulk interpretation, depending on the scaling dimension H it corresponds to creating in the asymptotic past either a conical defect or a black hole in the center of AdS 3 . In fact one could feasibly construct an experiment in which the value of λ is varied. First prepare an ensemble of systems in various energy eigenstates then at a fixed reference create the collapse state |V on top of the eigenstate by acting with the operators that create |V out of the vacuum.

Floquet's theorem and its implications
The main advantage of rewriting the Fuchs equation into Hill form is that it provides access to the general theorems that apply to Hill's equation. Some of these classic results will be quoted without proof in this section, these have all been taken from [1] which in addition contains all the missing proofs. Some definitions are required, given a Hill system (3.1) and a set of normalized solutions y 1 (x) and y 2 (x) such that y 1 (0) = 1, y 1 (0) = 0, y 2 (0) = 0, y 2 (0) = 1. (3.13) We can define a characteristic equation ρ 2 − y 1 (π) + y 2 (π) ρ + 1 = 0, (3.14) whose roots we denote by where p i (x) are periodic functions with period π. If ρ 1 = ρ 2 then Hill's equation has a periodic solution with either period π (when ρ 1 = ρ 2 = 1) or period 2π (when ρ 1 = ρ 2 = −1). Let p(x) denote this periodic solution and let y(x) be another linearly independent solution. Then

17)
with a constant θ. The case θ = 0 corresponds to y 1 (π) + y 2 (π) = ±2, y 2 (π) = 0, y 1 (π) = 0. (3.18) The Floquet basis of solutions is equivalent to the eigenbasis of solutions discussed in section 2, this theorem is essentially a slightly more powerful version of the statement that the monodromy matrix around the unit circle of the solutions of the Fuchs equation is restricted to the set of SU (1, 1) matrices. It shares the same conclusion that, depending on wether (3.14) has real solutions, the ODE solutions can broadly be divided into a class of solutions that grows/shrinks exponentially and a class that stays bounded at all times. The former will be referred to as unstable solutions and the latter as stable solutions. This brings us to the second theorem, which tells us how the unstable and stable solutions are distributed over parameter space The solutions of (3.1) are stable in the intervals At the endpoints of these intervals the solutions are generally unstable. This is always true for λ = λ 0 . The solutions are stable for λ = λ 2n+1 or λ = λ 2n+2 if and only if λ 2n+1 = λ 2n+2 , and they are stable for λ = λ 2n+1 or λ = λ 2n+2 if and only if λ 2n+1 = λ 2n+2 .
If we promote a basis of solutions to functions of x and λ as y i (x, λ) then the characteristic values solve the equations y 1 (π, λ n ) + y 2 (π, λ n ) = 2, (3.20) y 1 (π, λ n ) + y 2 (π, λ n ) = −2. (3.21) Due to the fact that λ = 1 − 24L 0 /c, the oscillation theorem is consistent with the physical intuition that if the energy is taken to infinity while keeping all the other Virasoro conserved charges fixed the state equilibriates to an unstable (black hole) phase with an arbitrarily high temperature. If we now think of the energy as a free parameter in the spirit of primary operator insertions presented in the last section then we can see that lowering the energy causes the solutions to settle down into a stable phase at some point λ 0 . In principle the oscillation theorem states that instabilities can reoccur if the energy is lowered further but this possiblity will be shown to be excluded. In order to connect the stability intervals of the oscillation theorem to the Virasoro coadjoint orbits we need one more theorem Theorem 3. Either all nontrivial real solutions of (3.1) have only a finite number of zeroes, or all real solutions have infinitely many zeroes. Let λ 0 be the smallest value of λ for which (3.1) has a periodic solution. Then for λ ≤ λ 0 all nontrivial real solutions have only a finite number of zeroes, but for λ > λ 0 , every real solution has infinitely many zeroes.
This theorem hints at the fact that the zeroth instability interval is special since it is the only one where solutions possess a finite number of zeroes.

Equivalence to Virasoro coadjoint orbits
At this point two ways of classifying stress tensors have been discussed, one can either sort T (z) into one of the Virasoro coadjoint orbits or one can try to find out in what (in)stability interval its associated Hill system falls. In this section it will be shown that both approaches are equivalent.
The basis of solutions to the Fuchs equation that diagonalizes the monodromy matrix around the unit circle corresponds to the Floquet basis of solutions of the Hill system. From this we can determine that the hyperbolic class of Fuchs equations, those with eigenvalues |µ i | = 1 correspond to unstable solutions of the Hill system. The set of stress tensors T (z) that lead to hyperbolic Fuchs equations is given by the set of hyperbolic Virasoro coadjoint orbits B n (b). The integer n counts the number of zeroes of the solutions of the Fuchs equation on the unit circle [5] in the diagonal basis, hence for the Hill system the number n corresponds to the number of zeroes on the interval (0, π) of the Floquet basis.
Theorem 3 states that the only instability interval whose solutions can have a finite number of zeroes is the zeroth instability interval λ ∈ (−∞, λ 0 ). The solutions associated to the B 0 (b) orbit stress tensors possess no zeroes in the Floquet basis, therefore they have to be contained within the zeroth instability interval. Conversely, all orbits B n>0 (b) lead to Hill solutions that possess at least one zero on the interval (0, π). because of the real-valuedness of the Floquet factors ρ i this zero has to be attributed to the periodic parts p i (x), hence on the entire real line these solutions have an infinite number of zeroes and they have to be contained within the higher instability intervals.
On physical grounds all orbits B n>0 (b) can be dismissed as they correspond to states which are conformally related to states with an energy lower than the vacuum. This implies that all instability intervals other than the zeroth interval can be dismissed as being unphysical.
From this we can easily establish that all stability intervals other than the first first are similarly unphysical. Simply take such a hypothetical state |V such that its associated Hill system falls in one of the higher stability intervals. Then we can construct a new state |V H = O H (x)|V by acting with a primary operator. As derived in section 3.1 this will only affect the Hill system by lowering the value of λ. Therefore if we are in a higher stability interval we can pick a scaling dimension H such that |V H is contained within a higher instability interval. As a result we conclude that O H (0)|V is an unphysical state and cannot be contained within the Hilbert of our CFT which implies that |V cannot be contained within the physical Hilbert space either.
To complete the equivalence one needs to show the converse statement, that every point on zeroth instability interval and first stability interval respectively correspond to a unique B 0 (b) orbit or C 0<ν<1 orbit. The proofs for these statements are quite a bit more technical and as a result have been relegated to appendix B. The underlying idea is straightforward though, there exists a simple bijective mapping of the Floquet factor α in e iαx and the constant representative of the orbit, in appendix B it is shown that on the zeroth instability interval and the first stability interval the Floquet factor is a monotonic function of λ, which interpolates every single possible Floquet factor. As a result every single constant representative is hit exactly once on the interval λ ∈ (−∞, λ 1 ).
In short we conclude that every physical bounded Virasoro coadjoint orbit has to correspond to a point on the zeroth instability interval or the first instability interval, and conversely that every point on these two (in)stability intevals in a one-to-one manner exhausts the full set of B 0 (b), C 0<ν<1 and the P + 0 orbits. This completes the proof for the proposed equivalence between physical orbits and the first two (in)stability intervals. As a direct corollary this provides a CFT version of the fact that there can only be one non-thermal to thermal gravitational phase transition in the bulk.

Putting upper bounds on the transition energy
Virasoro representation theory restricts the set of physically relevant characteristic values to λ 0 (other than λ 1 which will not be discussed). When λ 0 is converted into an energy scale L |V through λ 0 = 1 − 24 c L |V it corresponds to the energy at which the family of states crosses over to a state that eventually thermalizes at late time. Alternatively one can think of the energy scale L |V as a diagnostic tool that lets one easily determine whether a state |V will thermalize or remain in an integrable phase. Unfortunately determining λ 0 analytically is an extremely difficult task, either given by the solution to a path-ordered integral [2], or by the solution of an infinite determinant problem [1]. In this section we will instead derive a set of upper bounds for the transition energy L |V in terms of the Virasoro conserved charges. For this purpose consider the following eigenvalue problem which can be constructed out of Hill's equation d 2 dx 2 + λ + Q(x) y n (x) = σ n y n (x), y n (x + π) = y n (x), (4.2) from the oscillation theorem we know that if y(x) is constrained to have periodic boundary conditions y(x + π) = y(x) that the eigenvalues σ n have to be related to the characteristic values through σ n = λ−λ n . If we are given a stress tensor in the thermal phase we know from the arguments above that λ < λ 0 . Hence in the thermal phase all eigenvalues are negative and is a negative definite operator. Assuming that the set of periodic solutions y n (x) forms a complete set on the space of periodic functions we can derive the following inequality which should hold for any real-valued periodic function w(x). We can think of (−)H(x) as a Hamiltonian of some quantum system in which case estimating λ 0 corresponds to estimating the ground state energy. Note that the inequality (4.4) suggests that one possible strategy to sort T (z) into an orbit would be to check whether there exists a periodic function w(x) such that the inequality (4.4) fails to hold. Mind though that if T (z) is contained within the elliptic phase that H(x) will only have one positive eigenvalue, hence trial-and-error is extremely unlikely to produce a counter-example. By exploiting the fact that the spectrum of H(x) is bounded from above a more general set of bounds can be derived. This can be done by applying the variational method to estimate the ground state energy of a Hamiltonian. Since σ 0 is the largest eigenvalue of H(x) the following inequality has to hold for any mormalized periodic function w(x) Take the following normalized set of test functions w n (x) = 2 π sin(nx).
To obtain this expression the identities Q + λ = 1 − 24 c e 4ix T (x) and L n = 1 2πi dz z n+1 T (z) have been used. Since the largest eigenvalue of H(x) is given by σ 0 we can establish the following general upper bound This provides us with an infinite set of lower bounds on the black hole transition energy. Therefore if L 0 is smaller than the left hand-side for any value of n we know for a fact that the energy of the state |V is below the transition energy and we can conclude that the stress tensor has to be contained within an elliptic orbit, or to phrase it quantitatively, if there exists any n such that then T (z) has to be some element of one of the C ν elliptic orbits.
By considering the different set of test functions w n (x) = 2 π cos(nx) a similar set of lower bounds can be obtained 2 The interpretation of an orbit as a classical Virasoro representation suggests that the L n charges are due to lowering operators and as such the dependence on these numbers is in a sense redundant, the reflection condition (3.4) imposes that L n = L * −n this allows us to eliminate the lowering operator charges. Furthermore the two sets of lower bounds (4.9) and (4.10) can be combined to give This provides the following general necessary conditions for whether a state will thermalize L 0 > c 24 (1 − n 2 ) + |Re(L −n )|, (4.12) these are necessary in the sense that the state |V will exclusively thermalize if all of these condtions hold. The major downside is of course that these are only necessary conditions not sufficient conditions. The CFT interpretation is that |V is generically not a primary state but a descendent. Only the contribution to the energy from the primary contributes to the constant representative. These criteria tell us whether, after filtering out the energy contribution from a particular Virasoro raising mode, the energy is still sufficient to overcome the holomorphic Casimir energy equal to c/24. Alternatively one can think of the inequalities (4.12) as bounding an equation of state in the thermalizing phase if the Virasoro conserved charges are interpreted as thermodynamic equation of state variables.

Discussion
The universality of black hole formation and black hole thermodynamics of gravity theories on asymptotic AdS 3 suggests that the dual description of thermalization should equally be fixed by the universal features of 2d conformal field theory. This paper set out to demonstrate that a similar transition to a thermal state is fixed by a combination of Virasoro representation theory, properties of linear differential equations and the universal vacuum block contribution to correlators. As a consequence a set of upper bounds were derived that firstly let one, in some cases, determine if a given heavy state is elliptic (i.e. non-thermal) and secondly provides lower bounds on the deficient energy required to trigger thermalization. To probe thermalization it is typical to compute correlators of a set of simple operators and see if, to leading order in a statistical mechanics large system limit, they match thermal expectation values. This paper specifically considers two-point functions of probe operators, these operators have small scaling dimensions compared to the system. Therefore they can be interpreted as corresponding to low-energy excitations in radial quantization. In the case considered the large central charge expansion plays the role of the large system limit. In principal the number of probe operators can be scaled up relatively straightforwardly as long as the resulting product operator is parametrically light compared to the central charge. Physically this correponds to the fact that a sufficiently complicated quantum measurement can resolve the microstate of the ensemble which in the process destroys the statistical mechanics interpretation. Quantitatively the Virasoro Ward identity states that one would have to sum up all light contributions to the stress tensor expecation value if there are order O(c) of these O c 0 contributions the large c expansion will break down.
Intuitively black holes formed by dynamical collapse would be an ideal arena to study the information paradox. In fact since the leading order large central charge calculation reproduces the correlator decay discussed in Maldacena's formulation of the information paradox [33], it is tempting to compute perturbative corrections. Especially since the authors of [4] showed that these corrections can be computed systematically. Unfortunately almost exactly the same authors demonstrated that these perturbative corrections can not resolve the information paradox [34], see also [35] for a recent discussion. Perturbative 1/c corrections are not entirely devoid of practical meaning though, and are in fact very efficient for the purpose of calculating gravitational loop corrections. The first 1/c corrections of the probe part contributes the Eikonal resummation of ladder diagrams of graviton exchanges between the probe fields and higher order corrections resum loop corrections on top of that.
As is well-known, the Virasoro algebra is the central extension of the Witt algebra. As such we have to motivate a particular way to extend this algebra. The elements are easy, the basic elements of Witt algebra were the tangent vector fields on the circle parametrized by functions f (θ). To extend the elements we simply append to each function f (θ) a number to create a doublet (f (θ), t), these will be the elements of the extended algebra. The complication is to find a new bracket on this extended algebra, i.e given two doublets (f (θ), The standard CFT approach is to define the bracket through the 2d conformal Ward identity From this Ward identity we can associate the following charge operator to the infinitessimal conformal transformation (z) we associate the following Lie bracket to these conserved charges We know the effect of the infinitessimal transformation on the stress tensor hence we can work out the right-hand side of (A.6) Note that the boundary term coming from the integration by parts vanishes due to the assumed single-valuedness of T (z) and i (z). This candidate Lie bracket does not look antisymmetric yet, this can be fixed by performing two integrations by parts on the anomalous term This commutator algebra on the conformal charges of 2d CFT motivates the following bracket on the extended Witt algebra by transforming back to our coordinate system through z = e iθ . After some algebra we find Write the space (f (θ)∂ θ , a) in terms of the following basis It is then easy to check that we find the usual Virasoro algebra (in a slightly unconventional form) The additional generator Z commutes with all other generators, hence this is a proper central extension of the Witt algebra.

A.1 Coadjoint action
In the last subsection it was reviewed how the Witt algebra can be viewed as the Lie bracket of vector field on S 1 and the Virasoro algebra as a central extensions of this algebra. To define the adjoint and coadjoint action on these vector fields we need to define an inner product. The vector fields f (θ)∂ θ and quadratic differential fields b(θ)dθ 2 possess a natural inner product b(θ), f (θ) = The adjoint representation is defined through the action of the vector fields onto themselves i.e. 16) from this definition we can construct the adjoint action through the bracket (A.10). Inserting this into the inner product tells us the effect of the adjoint action Given the adjoint action we define the coadjoint action as the unique operator on the dual space such that its effect on inner products is the same as that of the adjoint action 3 , specifically Since the right hand-side is known we can compute what the effect of the coadjoint action is on a quadratic differential pair (b, a) (A.19) after performing some integrations by parts we obtain from which we can read off the effect of the coadjoint action This expression looks unfamiliar but if we apply the change of variables z = e iθ on the right hand-side of (A.20) we obtain g f (A. 22) hence in z-coordinates the coadjoint action takes the form of the variation of the 2d CFT stress-tensor under infinitessimal conformal transformations generated by g. This justifies the description maintained throughout the body of the text, that the sets of all stress-energy tensors that are connected to each other through single-valued conformal transformations is isomorphic to the coadjoint orbits generated from a quadratic differential field on the extended space by acting with the coadjoint action.

A.2 Classifying orbits
In the previous subsection it was shown that due to the transformation rule of the stress tensor it naturally lives inside a Virasoro coadjoint orbit, i.e. the space of all quadratic differentials that are continuously connected to some reference quadratic differential field. These coadjoint orbits possess a natural manifold structure [31], where we can think of the coadjoint action parametrized by a pair (f, x) as (almost) the tangent space at the reference point (b, a). This is not entirely accurate, there are possibly also pairs (h, y) that leave the reference point (b, a) invariant, these have to be modded out. Define the central extension of the algebra of vector fields on S 1 as dif f (S 1 ) and G b as the the algebra spanned by the vector fields (h, y) such that ad * (h,y) (b, a) = 0 specifically In that case the tangent space of the coadjoint orbit at the reference point (b, a) is given by dif f (S 1 )/H b . This tangent space can be integrated to a full manifold designated by Dif f (S 1 )/H b , these will be identified as the coadjoint orbits. As a result the coadjoint orbits are classified by their reference point and their stabilizer algebra. Turning back to the stabilizer algebra, since the coadjoint action is given by we know that given a fixed point (b, a) the defining equation for the stabilizer vectors h is 25) or after a change of variables z = exp(iθ) As a first step towards the well-studied [5] connection between Virasoro coadjoint orbits and the differential equations classified as Hill's equation we note that given a basis of solutions ψ 1 (z) and ψ 2 (z) to the linear ODE the independent solutions of (A.26) will be given by There is one additional complication, the condition (A.26) is a local condition, there is no guarantee that the local solutions (A.28) respect the periodicity condition of the circle. In a particularly impressive proof in [6] it is shown that the number of global stabilizing vector fields of a given reference point is always either 1 or 3.

Orbit Stabilizer algebra generators
Dif f (S 1 )/S 1 l 0 Dif f (S 1 )/P SL (n) (2,R) l 0 , l n , l −n Dif f (S 1 )/T n,∆ f (θ), where f is a vector field with 2n simple roots Dif f (S 1 )/T n,± f (θ), where f is a vector field with n double roots Table 3: The classification of Virasoro coadjoint orbits purely in terms of group manifolds, i.e. without their reference points The classes of solutions to (A.26) have been classified [6], the resulting (either 1 or 3 dimensional) stabilizer groups that can be uplifted from the solutions are typically listed as [5] S 1 , , P SL (n) (2, R), T n,∆ ,T n,± . (A.29) Respectively, S 1 forms the set of rigid rotations of the circle, P SL (n) (2, R) is the group generated from the subalgebra of dif f (S 1 ) spanned by (L 0 , L −n , L n ). T n,∆ andT n,± are one-dimensional groups whose Lie algebras consist of vector fields on S 1 with respectively 2n simple zeros or n double zeros, since the number of zeros is an orbit invariant The list of these group manifolds that form the various Virasoro coadjoint orbits has been summarized in table 3 [14].

A.3 Hill's equation, Virasoro coadjoint orbits and SL(2, R) conjugacy classes
So far this concludes the classification as presented in [6], the group manifolds formed by the coadjoint orbits form somewhat intangible structures. In this section we will follow [5] and link the Virasoro coadjoint orbits to something that is of direct physical interest for our purposes; the conjugacy classes of the monodromy matrices of the solutions ψ(z) of around the unit circle. The information of the coadjoint orbit enters the ODE through T (z) which is related to the earlier quadratic differential fields b(θ)dθ 2 through T (z) = b(x(z)) with z = e iθ . In order to directly connect to application in 2d CFT described in section 2 we identify T (z) with the stress tensor expectation value of the CFT used throughout the body of the text, in that vein c designates the central charge of theory. Being a second order linear differential equation, the solution space of (A.30) is spanned by a basis of solutions of the form ψ ψ ψ(z) = (ψ 1 (z), ψ 2 (z)) T , the monodromy matrix ψ ψ ψ(z) around the unit circle M ψ ∈ SL(2, R) is given by ψ ψ ψ e 2πi z = M ψ ψ ψ ψ(z) (A.31) The infinitessimal transformation under a vector field (f, y) was given in (A.22) it is known that this expression integrates to the transformation rule The set of eigenvalues of an element of SL(2, R) are invariant under conjugation, the det(M ψ ) = 1 constraint restricts the two eigenvalues of M ψ to be each others inverse, Floquet's theorem further restricts the relevant classes to either be purely real or pure phase with cross-over points where both eigenvalues are either 1 or both -1. The orbits with eigenvalues falling in these categories are respectively designated as hyperbolic, elliptic and parabolic or exceptional orbits. By constructing explicit examples of stress tensors T (z) whose associated solutions have a monodromy within these conjugacy classes the authors of [5] list the following classes of Virasoro orbits B n (b), C(ν), P ± n , E n (A.38) these respectively correspond hyperbolic, elliptic, parabolic and exceptional orbits, but they can have different stabilizer algebras corresponding to the affixes n and ± the real arguments b and ν parametrize the set of representative stress tensors.

A.4 Relevant qualitative features of the individual orbits
As we shall see not all the orbits mentioned in the previous section are of direct physical interest though. To each stress tensor we an associate an energy expectation value through the conserved charge L 0 while within each orbit this number is unbounded from above, there are only a few orbits for which this number bounded from below. Since we are considering unitary CFTs the energy expectation value on the entire spectrum of states is bounded from below by the energy expectation value of the vacuum state (in any conformal frame). Therefore one can dismiss any orbit without a lower bound as incapable of being associated to a physical state of the CFT spectrum. The full list of orbits which possess a lower bound on their energy is given by In addition some of these orbits have the special property that they contain constant representatives 4 these orbits and their representatives are given by Orbit Stabilizer SL(2, R) class Bounded energy Constant representative P SL (n) (2, R) exceptional no no Table 4: Relevant features of the Virasoro coadjoint orbits as classified in [5].
The representative functions T (z) have the form of the stress-tensor expectation values on a state given by conformal primary inserted at the origin acting on the vacuum with a scaling dimension given by the numerical prefactor of T (z). This also serves as a first indicator that something seems to change at the BTZ mass threshold c/24. Some of the details of the classification have been summarized in table 4.

B Floquet factor as a monotonic function of λ
The number λ contained within Hill's equation possesses within our context the interpretation as being related to the total energy of our state on the radial plane through λ = 1 − 24 c L 0 . It would hence seem sensible that decreasing λ would increase the mass of the conformal primary associated to an orbit. In this section we will prove that this is the case, in more practical terms, this establishes that as we decrease λ from λ 1 to negative infinity we sweep every physically allowed orbit and we reach every orbit exactly once. The scaling dimension H of the primary representative of an orbit T (z) = H/z 2 is simply related to the Floquet factor α in (3.16) through α = 1 − 24 H c , (B.1) therefore we will argue the equivalent statement that varying λ from λ 1 to negative infinity monotonically interpolates every possible Floquet factor. The proof is somewhat technical, mostly due to the fact that any proof has to naturally exclude the instability intervals that are not the zeroth interval due to the clear fact that the Floquet factor cannot be monotonic within these intervals. The proof for the first stability interval has already been done in the past, in [1] it is proven that Hill's discriminant is monotonically decreasing within this interval.
To prove the monotonicity of the Floquet factor in the zeroth instability interval consider the Floquet basis of solutions y 1 (x) = e −βx p 1 (x), y 2 (x) = e βx p 2 (x) with p i (x + π) = p i ( The same consideration that applied to y i (x) apply to v i (x), i.e. the functions v i (x) can be taken to be strictly greater than zero for all x. The solutions of Hill's equation should continuously flow into each other as we vary λ. Consider the solutionỹ 2 (x) to (B.5) obtained from the exponentially increasing y 2 (x) by continuously shifting λ → λ − . As a solution it should be expressible as a linear combination of Floquet solutions y 2 (x) = av 1 (x) + bv 2 (x) = ae −γx q 1 (x) + be γx q 2 (x), (B.7) since y 2 (x) is positive everywhere we can choose an such that the solutionỹ 2 (x) is positive everywhere as well, this implies that both a, b > 0. Now consider the ratioỹ 2 (x)/y 2 (x) in the limit x → ∞ lim x→∞ỹ 2 (x) y 2 (x) = lim x→∞ ae −γx q 1 (x) + be γx q 2 (x) e βx p 2 (x) = lim x→∞ be (γ−β)x q 2 (x) p 2 (x) . (B.8) Since q 2 (x) is a periodic function and hence bounded and since p 2 (x) can vanish nowhere we can conclude that ifỹ 2 (x)/y 2 (x) blows up in the limit x → ∞ this implies that γ > β.
We will now construct the solutionỹ 2 (x) for small by perturbing around the solution y 2 (x) i.e.ỹ 2 (x) = y 2 (x) + z(x). (B.9) To linear order in we find that z(x) satisfies the ODE z (x) + (Q(x) + λ)z(x) = y 2 (x), (B.10) by means of variation of parameters we find that the relevant particular solution for z(x) is given by z(x) = y 2 (x) x 0 y 1 (x )y 2 (x )dx − y 1 (x) x 0 y 2 (x ) 2 dx = e βx p 2 (x) x 0 p 1 (x )p 2 (x )dx − e −βx p 1 (x) x 0 e 2βx p 2 2 (x )dx . (B.11) Due to the periodic nature of p 2 (x) it is bounded from above, take M to be an upper bound such that for all x p 2 (x) < M then z(x) > e βx p 2 (x) The negative term in parenthesis is bounded whereas the positive term diverges linearly with x in the limit x → ∞, again due to the positive definite and periodic nature of p i (x). Looking back at the ratioỹ 2 (x)/y 2 (x) in the large x limit this ratio diverges to positive infinity in the limit x → ∞, from which we can establish that γ > β. This completes the proof that the Floquet factor is a stricly decreasing function of λ as long as λ < λ 0 .