Tauberian-Cardy formula with spin

We prove a $2$ dimensional Tauberian theorem in context of $2$ dimensional conformal field theory. The asymptotic density of states with conformal weight $(h,\bar{h})\to (\infty,\infty)$ for any arbitrary spin is derived using the theorem. We further rigorously show that the error term is controlled by the twist parameter and insensitive to spin. The sensitivity of the leading piece towards spin is discussed. We identify a universal piece in microcanonical entropy when the averaging window is large. An asymptotic spectral gap on $(h,\bar{h})$ plane, hence the asymptotic twist gap is derived. We prove an universal inequality stating that in a compact unitary $2$D CFT without any conserved current $Ag\leq \frac{\pi(c-1)r^2}{24}$ is satisfied, where $g$ is the twist gap over vacuum and $A$ is the minimal"areal gap", generalizing the minimal gap in dimension to $(h',\bar{h}')$ plane and $r=\frac{4\sqrt{3}}{\pi}\simeq 2.21$. We investigate density of states in the regime where spin is parametrically larger than twist with both going to infinity. Moreover, the large central charge regime is studied. We also probe finite twist, large spin behavior of density of states.

The Cardy formula [1] for the asymptotic density of states has recently been rigorously derived with an estimate for the error term in [2,3]. A natural question is to ask whether one can generalize the formalism so as to make it sensitive to the spin or equivalently to the conformal weights h,h separately. This necessitates working out a 2 dimensional Tauberian theorem, which we achieve here. The motivations for investigating Cardy formula on (h ,h ) plane are several. First of all, the notion of infinity on a 2d plane is richer than ∆ → ∞ limit. We will see that the finer details of the Cardy formula actually depend on how infinity is approached unless one makes extra assumption about the spectrum. Furthermore, there have been interesting developments in the direction of lightcone bootstrap in recent times [4][5][6][7][8], our analysis puts some of these results on rigorous footing. Another amazing feature is the ability to investigate the "areal" notion of spectral gap. If we probe the (h ,h ) plane with circular areas of radius R, centered at (h,h), then we find the optimal value of R which guarantees that the area contains at least one state. Again unless we put in extra assumption, the value of R depends on how infinity is approached and thus showing a richer asymptotic behavior. If we one assumes existence of twist gap, it turns out that the twist gap is complementary to asymptotic spectral gap in some sense, which we will make precise in due course.
The naive Cardy like analysis provides us with an expression for the asymptotic density of states where h andh are of the same order. One can re-express this as a function of dimension ∆ and spin J with ∆ J. Now a natural question is to ask whether the result is valid when ∆ and J is not of the same order. For example, in the large charge expansion literature [9][10][11][12][13][14][15][16][17][18], the regime where J ∆ 1/n with n > 1 is being probed. It turns out that only a part of the answer coming from the naive Cardy like analysis is meaningful while the rest of it is comparable to the error term. We emphasize that the analysis is only possible because now we have a rigorous estimate of the error term due to the Tauberian theorem that we prove in this paper.
With our rigorous treatment, it is possible to address issues regarding whether we can trust the naive Cardy formula when h andh are not of the same order. It turns out that the answer to this question is intimately connected with the existence of twist gap. We show that we can trust the naive Cardy formula for all the operators when max(h,h) = min(h,h) Υ with 1 ≤ Υ < 2. It is also shown that with the assumption of twist gap, the validity of Cardy formula for primaries for c > 1 CFTs does not require any restriction on Υ.
The another motivation for taking up a rigorous study of Cardy formula is to be able to probe the large central charge (c) sector, to be specific, to derive the density of states when h/c,h/c are finite but c is very large. This part is in the spirit of result derived in [19]. A nice feature that reveals itself through the rigorous treatment is a curious connection between validity of Cardy regime and the twist gap above the vacuum. These features are important in the context of holography.
The plan of the paper is to quote the main results here in the beginning and discuss its consequences in terms of CFT data, such that the current section can be thought of as mostly self contained. The next section §2 gives some intuitive understanding of the technical stuff that follows. From §3 onwards, we plunge into technical proofs with a healthy relaxing intermission in §6, where we numerically verify our results. For readers going for a really quick ride, we have highlighted the main equations and results in what follows.

Integrated density of states
We prove a 2 dimensional Tauberian theorem in context of 2 dimensional conformal field theory. The asymptotic density of states with conformal weight (h,h) → (∞, ∞) is derived using the theorem. We find that the error term is controlled by the twist parameter. We note that as (h,h) → (∞, ∞), the twist also goes to ∞. We remark that the regime of validity depends on whether we put in the assumption of having a twist gap.
Definition: by finite twist gap, we mean there exists a number τ * > 0 such that there is no operator with twist τ ∈ (0, τ * ) and there are finite number of zero twist operators 1 with dimension less than c/12.
We make two remarks: a) the fact that there are finite number of zero twist operators with dimension less than c/12 is always true since there are finite number of operators with dimension less than c/12 for finite central charge, b) Not having any operator with twist τ ∈ (0, τ * ) disallows having 0 as twist accumulation point.

Main theorems on integrated density of states
No assumption on twist gap: We show that for finite central charge c, the number of states with conformal weights less than or equal to some specified large conformal weight h,h is given by: (1) h,h→∞  Assuming a twist gap: It turns out that if we assume a finite twist gap, we can trust eq. (1.1) even when h andh are not of the same order but h =h υ with 1/2 < υ < 2. In such a scenario, the error term becomes O(τ Υ 4 −1/2 ), where Υ = max (υ, 1/υ). The Υ characterizes how h andh are of different order asymptotically in a symmetrized fashion, for example, if we approach the infinity along the curve h =h 1.1 orh = h 1.1 , we have Υ = 1.1. Thus our error estimation is symmetric if we reflect the line of approach to infinity about h =h line.
We have for 1 ≤ Υ < 2, (1. 2) The eq. (1.1) and eq. (1.2) are two of the central results obtained in this paper. If h andh are not of the same order, we basically probe the large spin sector of density of states 3 , to be precise, the regime where spin is parametrically larger than the twist but both goes to infinity.
The basic structure of both the eq. (1.1) and eq. (1.2) is that they have leading 2 If we say f = O(1), we mean |f | < M for a fixed positive number M . 3 A cautionary remark is that here in this paper unless otherwise mentioned, the twist is NOT kept finite while taking this limit. This can be contrasted to the scenario in the usual large spin expansion [21], where one keeps the twist finite. exponential piece multiplied with a subleading polynomial suppression. The error term is then further suppressed by a polynomial piece. Now if Υ ≥ 2, one can see the error term in (1.2) is not really suppressed, hence is not in fact an error term. Thus we can not trust the polynomially suppressed terms. In this regime, we are able to show that We further remark that for CFTs where the partition function nicely factorizes into holomorphic and antiholomorphic pieces, the leading result directly follows from the analogous result for large ∆ = h +h, proven in [2], nonetheless the error term in analogues of eq. (1.1) and eq. (1.2) goes like O(h −1/2 ), hence, in such a case, we have more control over the approximation. Below we will digress a bit and touch upon some of the interesting results that can be extracted from the above before coming back to summazing our main results in the next subsection §1.2.
Rich structure of asymptotic approach: The integrated density of states show distinct leading behavior depending on how the asymptotic infinity is approached.
In [2], it has been shown that as ∆ → ∞, we have We remark that in the asymptotic limit, both F MZ (∆ → ∞) and F (h → ∞,h → ∞) count the total number of operators. But these functions approach infinity in a different manner (see the figure 1). To be concrete, let us choose h =h = ∆/2, thus we have  Figure 1: Approaching to infinity on (h ,h ) plane: The blue lines denote how the number of operators with (h ,h ) is counted such that h +h ≤ ∆ and then we let ∆ → ∞, this is given by the function F MZ (∆), originally calculated in [2]. On the other hand, the black lines denote how the number of operator is counted upto some value of h =h = ∆/2 i.e. h ≤ ∆/2,h ≤ ∆/2 and then we let ∆ → ∞. This approach to infinity is captured by the function F (∆/2, ∆/2), calculated in this paper. We see that the square of size ∆/2 with one vertex at origin and another one at (∆/2, ∆/2) are always contained within the rightangled triangular region, created by h axis,h axis and h +h = ∆ line. This is consistent with the observation that leading behavior of F (∆/2, ∆/2) is suppressed compared to F MZ (∆). i.e.
We see that the square S of size ∆/2 with one vertex at origin and another one at (∆/2, ∆/2) is always contained within the rightangled triangular region T, created by h axis,h axis and h +h = ∆ line. This is consistent with the observation that leading behavior of F (∆/2, ∆/2) is suppressed compared to F MZ (∆). In fact, one can similarly study the distribution of the operators in rectangular (or square) areas such that the rectangle is contained within T, and one vertex is on the line h +h = ∆ (see the figure 2). This study reveals that the among such areas, the square S contains the most number of operators while any other rectangular region contains fewer number of operators, in fact the number is exponentially suppressed compared to that of the square S .
Spin sensitivity of the asymptotics: One can make a detailed analysis of spin sensitivity of the above result, which we expound on §5.2.  Figure 2: (h ,h ) plane : asymptotically, the rectangular region (blue shaded) contains exponentially less number of operators compared to square region (red shaded). They are contained within the rightangled triangle, created by h axis,h axis and h +h = ∆ line. Here ∆ = 12.
Windowed entropy with respect to h andh: An immediate consequence of the eq. (1.1) is the expression for "windowed" entropy S δ,δ . The windowed entropy is defined as logarithm of number of states within a rectangular window of side length 2δ and 2δ, centered at (h,h). This is analogous to entropy defined as in microcanonical ensemble by proper "binning", where the bin size is dictated by δ,δ. As we take h → ∞,h → ∞, we can keep the bin size δ,δ order one or let them scale like h α and h α respectively. We find that where for 3/8 < α ≤ 1/2, we have, : where the functions s ± (δ,δ) are determined in the section §3, in particular, we have s ± ≡ exp(c ± ), and c ± is given by (3.25). We remark that when the bin size is large, there is a universal correction to Cardy formula given by the sinhyperbolic functions. This is analogous to what is found in [2] from the analysis sensitive to dimension only.
Windowed entropy with respect to ∆ + J: One can define a microcanonical entropy with respect to ∆ + J = 2 max{h,h} (name this parameter κ) as The asymptotic behavior of S κ δ is given by where for large enough bin size (δ κ α ) we have

c > 1 CFTs-results specific for primaries
One can make the results in the previous subsection specific to Virasoro primaries only, in fact do better. This boils down essentially repeating the argument presented in §3, §4 and §5 with minor modification. The idea of extending the argument from §3, §4 and §5 to this case is similar in spirit and practice to how [2] obtained the specific results for primary using methods suitable to study all the operators. The details can be found in §3, specifically eq. (3.37) onwards. Without much ado, here is the result: for finite central charge c, we find the integrated density of states specific for primaries behave like (from now on, we will be using the superscript "Vir" to denote the result specific for primaries): (1.12) The "windowed" entropy (we have considered bin of size 2δ by 2δ just like what we did for the analysis of all the operators) for Virasoro primaries is given by (1.13) where for 1/8 < α ≤ 1/2, we have : where the functions s ± (δ,δ) are the same functions that appear in the analysis for all the operators.
Large spin, large twist sector for primaries: If we assume a finite twist gap (as defined in {1.1}), the result given in eq. (1.12) is true irrespective of whether h andh are of the order one or not. Thus unlike the case for all the operators, here we can trust the polynomially suppressed correction for all values of υ, where h =h υ .

Large spin, finite twist sector
The large spin, finite twist sector is not entirely asymptotic regime since the quantity knows about low lying spectrum in one of the weights. It turns out we can only put an upper bound in this case. There is an O(1) error in the estimation. While for the upper bound this does not cause any trouble, for the lower bound, it makes thing tricky. In particular, the lower bound on the density of states, appropriately integrated, contains an exponential piece as expected from extended Cardy formula [4,5,7,8] but it comes with a multiplicative order one number, which can become negative unless proven otherwise.

Analysis for all the operators:
In what follows, we will keep h finite and let h → ∞, the windowed entropy S ft δ,δ is found to be bounded above by where M is an order one number. Here S ft δ,δ and S ft h,δ,δ are defined as The number M is given (or estimated) by where χh is the character for the conserved current with weight (h, 0),h ≥ 0 (including the Identity) and c ± is an order one h,h independent number, defined in §7. M is a finite number as the absolute value of the sum overh is bounded above by the partition function evaluated at β =β = 2π, which is a finite number.
Analysis for primaries with/without conserved currents: The above result can also be made specific to primaries: and M Vir is an order one number, given by where the zero twist primaries have weight (h, 0),h ≥ 0 (including the Identity) and c + is an order one h,h independent number, defined in §7. M Vir is a finite number since the sum inside the log is convergent. This happens because the absolute value of the sum is bounded by the partition function evaluated at β =β = 2π, which is a finite number.
Analysis for primaries for CFT with no nontrivial conserved current: If we assume that there is no nontrivial conserved current i.e the only zero twist primary is the Identity and there is a finite twist gap (the finite twist gap as defined in {1.1}, combined with the absence of nontrivial conserved current implies the usual twist gap condition used in the literature, for example in [20]). we show that 19) where M is an order one h independent number. If we assume that (h + δ − c−1 24 ) is a very small number compared to 1 c−1 , this matches with the leading result appeared in the lightcone bootstrap program [4][5][6][7][8] i.e.
(1.20) whereM = M + log(4π 2 ). We remark that the limit is very subtle here. There are several scales. The scale set byh is the largest one and we are seeking an asymptotic behavior inh. Then there are two fixed parameters h and c. We are probing the regime where (h + δ − c−1 24 ) is a very small number compared to 1 c−1 . The details of the calculation can be found at the end of §7.

Asymptotic spectral gap
The idea about deriving an upper bound on spectral gap comes from binning the states. If we make the bin size very small, we can not prove a positive lower bound on the number of states in that bin, because the bin might not have any state at all. As we increase the bin size, the chances are more that we find such positive lower bound. If we find a positive lower bound for a specific bin size centered at some large h,h; that would immediately imply existence of an upper bound on the asymptotic spectral gap.

Probing spectral gap via Circle of order one area
With/without twist gap: Here we do not put any assumption on twist gap.
Let us consider a square S of side 4 √ 3γ π + g centered at (h,h) on (h ,h ) plane. Here g can be any arbitrarily small positive number. In the limit h → ∞,h → ∞ we have where the asymptotic region is reached along a curve for which max(h,h) Thus the spectral gap along this curve is bounded above by a circle of radius γr √ 2 and the best possible value of r that we find is r = 4 √ 3 π , this being the circle circumscribing the square.
An immediate corollary is that the asymptotic twist gap is upper bounded by 42γ. For c > 1, the argument can be made specific for primaries, hence the asymptotic gap. This in some sense complements the bound on primary twist gap over the vacuum 4 [8,20].
We suspect that either by suitable choice of function or by the better estimate of heavy sector of the partition function, r and/or length of a side of the bounding square can be made to 1. If this can be done, then the bound becomes optimal for γ = 1, (assuming we always consider circular/square region) since tensoring chiral Monster CFT with antichiral Monster CFT saturates the bound. One can see the saturation by circumscribing a square of unit length by a circle of radius 1 √ 2 on (h ,h ) plane [See the fig. 3]. Nonetheless, the optimality along a curve for γ = 1 is not guaranteed. An immediate corollary of finding such a circle is that the asymptotic twist gap is upper bounded by 2 √ 2 along h =h curve. The same bound holds for asymptotic primary twist gap. We remark that in terms of twist, the above gap might not be optimal, since if we tensor chiral Monster CFT with anti-chiral monster CFT, the asymptotic twist gap is 2. If one can find a bounding square of side length given by 1, that would reproduce the optimal twist gap 2.
The above result and the conjectures can not be applied to a scenario, where infinity is approached along a curve where h is of widely different order compared tō h, in particular, say, if we want to approach the infinity along the curve h υ =h with υ = 1. To circumnavigate this issue, we assume existence of twist gap g. We remind the readers that by finite twist gap, we mean there exists a number τ * > 0 such that there is no operator with twist τ ∈ (0, τ * ) and there are finite number of zero twist operators with dimension less than c/12, and g ≥ τ * . Moreover, assuming existence of g helps us to get rid of dependence on γ. One can always circumscribe a square of unit length by a circle of radius 1 √ 2 . Thus any circle centered at (h,h) and of radius 1 √ 2 + g with g > 0 would contain at least an operator. This is consistent with our result that asymptotically any circle of radius γr √ 2 + g with γ ≥ 1 will contain at least one operator where γ is defined as the fourth root of asymptotic ratio of max(h,h) and min(h,h). We show r = 4 √ 3 π > 1 and suspect that r can be made to 1. Here we have shown three such examples of containment. For each choice, two circles with different radius have been drawn to show that a circle of γ √ 2 + g will contain at least one operator, as long as g > 0.
On the diagonal we have γ = 1. The red dots denote the contained operators. The black circle shows that one can not ensure containment with circle of radius less than 1 √ 2 , showing optimality along the h =h curve if one can show r = 1. We remark that as we have only shown that r = 4 √ 3 π > 1 (the dashed black circle), we have not yet reached the optimal bound. The blue circle suggests that the optimal bound should be 1 √ 2 even for the case where h =h. Assuming a twist gap makes this bound insensitive to the line of approach towards infinity, but sensitive to the gap. See fig. 4.
CFT with twist gap g: Now we assume that the CFT has a twist gap as defined in {1.1}.
For a CFT with twist gap g (as defined in {1.1}) and central charge c > 1 , one can have a bounding circle C specific to primaries having a radius σr irrespective of how infinity is approached, such that the bounding circle contains at least one operator.
Thus for such a scenario there exists h * andh * , two order one numbers such that Again this is obtained by circumscribing the appropriate bounding square [See the fig. 4]. The superscript "Vir" on ρ Vir denotes that it is density of primary operators as opposed to all the operators. In a compact unitary 2D CFT without any conserved current, one can use the upper bound of twist gap due to Hartman, appearing in [20], to deduce Now we will explain the sense in which the minimal gap is complementary to twist gap: Suppose we consider a 2D compact unitary CFT with twist gap g such that it does not have any zero twist primaries (conserved currents) except the Identity: if asymptotically there exists a circle of minimal area a A on (h ,h ) plane such that it does not contain any operator, we immediately deduce the following inequality a One can imagine that operators having conformal weight (h,h) are denoted by the point (h,h) on (h ,h ) plane. We name this set to be S. Now consider the set is the Euclidean distance on the plane between the points a and b. Existence of a circle with minimal area means Inf S d > 0. Asymptotically minimal area means that we look at the plane for h > h * &h >h * and construct the set S d and consider its infimum.
This can be thought of as an upper bound on twist gap if the minimal areal gap A is known (note that minimal areal gap obtained from the full spectra has to be less Figure 4: Assuming twist gap g: operator spectrum of chiral Monster CFT tensored with its antichiral avatar on (h ,h ) plane: each vertex in the lattice represents the presence of operators. One can always circumscribe a square of unit length by a circle of radius 1 √ 2 . Thus any circle centered at (h,h) and of radius κ √ 2 + g with g > 0, κ ≥ 1 would contain at least an operator. This is consistent with our result that asymptotically any square of side length 4 √ 3 π , hence any circle of radius σr √ 2 + g will contain at least one operator, where σ = max Again we suspect that r can be made to 1. In the Monster example g = 4, hence σ = 1.
than or equal to the asymptotic minimal gap). If one can make r = 1 and show that The similar analysis can be done for all the operators assuming a twist gap. The only difference is that σ would be given by σ = max 1, c 12g . We elucidate on these bounds in the §3.

Probing spectral gap via Strips
Instead of squares , we can think of covering the (h ,h ) plane via strips of finite width and ask what is the minimum width of the strip that guarantees existence of at least one state in the strip. We can consider three kind of strips (see figure 5) : (1.28) Figure 5: We consider three kind of strips: the red vertical one is H 2 (h), the blue horizontal one is H 1 (h) and the black one is H 3 (∆). In each cases, we see that there is a minimum width of the strip such that the strips contain at least one operator. If we tensor chiral Monster CFT with its anti-chiral avatar, we can not make the width of the horizontal and vertical strip less than 1. For g > c/12, which is the case for this tensored CFT, the bound we find using the function given in [3] is in fact 1, thus optimality is achieved for g > c/12. Nonetheless, it is not clear whether the bound (which is √ 2) found for the case g ≤ c/12 is optimal or not. A technical remark is that if it is not optimal, that would hint that one could have estimated the heavy part of the partition function in a better way.
Thus the asymptotic spectral gap is bounded above by 1.
In this work we show that The same result holds true for the H 2 strip with h replaced byh, where H 1 , H 2 are defined in (1.26). This comes from putting a positive lower bound on the right hand side of the following inequality: We achieve this as a corollary of the lemma proven in §4. This shows that Asymptotically on the (h ,h ) plane, if the width of the horizontal or vertical strip is bigger than Θ(g − c/12) + √ 2Θ(c/12 − g), it contains at least one Virasoro primary. This is in fact the optimal bound for g > c/12 and saturated by tensoring chiral Monster CFT with its anti-chiral avatar.

Analysis at large central charge
We consider the c → ∞ limit and parametrize the conformal weights in following way: Let us define * = min ( ,¯ ) and * = max ( ,¯ ).
With/without twist gap: We show that 5 For min (h/c − 1/24), (h/c − 1/24) = * > 1 6 , the microcanonical entropy for order one window δ,δ O(1) is given by (1.34) Assuming finite twist gap: If we assume a finite twist gap g (as defined in {1.1}), then The regime of validity of the Cardy result as in (1.33) can be extended to * > We remark that if one assumes that g > c 12 , then it is possible to achieve the regime of validity of the Cardy fromula for * > 1 24 . This does not follow automatically from the above eq. (1.35); rather one has to go through the analysis from the scratch using a different separation of light sector and heavy sector, where the light sector contains the Identity alone and the heavy sector has the states with h,h > c/24. This result is consistent with the fact that one can show Cardy like behavior for * > 1 24 if tensor product of two copies of extremal CFT with c = 24k is studied in the k → ∞ limit.

Relevant recent work:
There has been recent surge in analyzing the asymptotics of CFT data on a rigorous footing. The results have been obtained [22,23] using techniques borrowed from a part of mathematics literature, which goes by the name of Tauberian theorems. The appendix C of [24] emphasizes the importance of Ingham theorem [25] in analyzing Cardy's result [1] for the asymptotic density of states in 2D CFT. Subsequently, the complex Tauberian theorems, as appeared originally in [26] is utilized in the work of [27]. A complete rigorous treatment of Cardy formula appeared in the work [2], where they figured out the density of states in ∆ → ∞ limit with a rigorous optimal estimate of the error term. The improvement of the result along with a proof of the conjecture made in [2] has been put forward in [3]. An rigorous analysis of the asymptotics of three point coefficients [28] appeared recently [29] where the main challenge was to circumnavigate the negativity issue for the analysis of three point coefficients.

Set up
We consider a 2D CFT with spectrum of operators having conformal weights (h ,h ). We assign different real temperatures β,β to the left-moving and the right-moving sectors respectively. The partition function Z(β,β) is given by: The modular invariance of the partition function yields: We further define the following measure Our goal is to estimate the integral of the measure dF over different regions.

A semi technical glimpse of the subtleties
One of the key step in deriving the Cardy formula is the intuitive understanding that at high temperature, the partition function is dominated by the heavy states, thus doing an inverse Laplace transform of the high temperature behavior of the partition function should produce the asymptotic density of states. Schematically, Produces asymptotic density of states + Error , (2.4)
= Error in asymptotic density of states . (2.5) The underlying assumption while doing the above is that the inverse Laplace transform of the error term is bounded as well, thus producing an error compared to the leading behavior of the asymptotic density of states. The Tauberian formalism justifies this step by carefully estimating the error terms. The way it works is following: one bounds the number of states within an order an window centered at some heavy ∆, from above and below by some convolution ( ) of partition function at high temperature (β) and bandlimited function φ ± , schematically this looks like Intuitively, at this stage, we know that heavy states contribute to this partition function at high temperature. Now we implement modular transformation, the par-tition function at high temperature becomes partition function at low temperature. Schematically we have At low temperature, low lying states contribute the most i.e Z 4π 2 β → ∞ is dominated by low lying states. So, following [19], we separate the low lying states from heavy states; the low lying states constitute the "light" part, while the "heavy" part is complement of that. The "light" part contains finite number of operators at finite central charge. So we do the inverse Laplace transform of this "Light" part to get the leading answer ρ * . Thus ρ * reproduces the Z 4π 2 β → ∞ behavior upon doing Laplace transformation. We are still left with the "heavy" part contribution of Z 4π 2 β+it . This part can be shown to produce a subleading correction to the leading piece of asymptotic density of states, thus justifying (2.5) using the bound proven in [19]. This requires relating the "heavy" part at temperature β to the "light" part at temperature β . The upshot of this discussion is that we have a full control of the error term and its inverse Laplace transformation. In practice, we only consider the Identity operator among all the operators in the "light" region. So one might worry about the error coming from that but since there are finite number of operators in the "Light" region, one can do inverse Laplace transformation term by term and show that each of them is exponentially suppressed and hence the finite sum of them. We emphasize the "finiteness" of finite sum is really very important for this and this is precisely why we need to treat heavy part separately 6 . In fact, we remind the readers that in the large central charge limit, we have infinite number of operators even in the "light" sector, hence we need an extra assumption of sparseness, as done in [19].
The immediate generalization of this technique used in the [2] has obstacles because of the cross terms present in the analysis, for example, say contribution from the states where h is large buth is not that large. The most obvious way to generalize the argument is to use the generalized HKS [19] cut: dividing the (h ,h ) plane into two regions, where the "light" region (call it L) contains all the operators with one of h orh being less than c/24 while the heavy region is the complement of them. It is possible to make a similar statement about this "heavy" region, relating it to the "light" part using HKS like argument. Nonetheless, one then stumbles upon the issue of defining ρ * , which is supposed to reproduce the leading contribution to the partition function at high temperature (β → 0, β → ∞), to be precise, the light part of the partition function at temperature β . Now the issue is that there are infinite number of operators in this region L. Unlike the case in [2], we just can not take the Identity operator to prove that this produces the leading behavior and say the rest are suppressed. In particular, the previous argument of term by term exponential suppression fails because there are infinite of them. The take home message is that it is not a priori clear whether just considering the vacuum to calculate ρ * is good enough, because infinite number of other operators might conspire to spoil the "leading contribution", even if each one of them is exponentially suppressed. We reemphasize that [2] did not face this problem, since in their case, the light region was ∆ < c/12 and the region has finite number of operators and everything is under control, so in principle their ρ * was defined having contribution from all those states with ∆ < c/12 and in practice derivable from the vacuum. To circumnavigate this problem, in §3, we use the original HKS cut i.e. we define the "light" region to be the one where h +h < c/12 and the "heavy" region is the complement of the light region. See fig. 6. The immediate cost for doing this is that we can make comment  Figure 6: The (h ,h ) plane: the light red shaded region contributes to Z L and given by the states with h +h < c/12. The complement of this contributes to Z H . The naive HKS cut makes the union of red and blue shaded region to be the "Light" region. But this brings in the problem of having infinite number of operators in the light region. This cut is applicable to the analysis sensitive to primaries as well with a shifted central charge c → c − 1. The picture on the right emphasizes that one needs to treat the zero twist operators separately for some part of the analysis.
only when h andh is of the same order asymptotically. The details depend on how infinity is approached. But this is expected because intuitively the infinity can be reached several ways on a plane. This issue persists for the analysis sensitive to primary as well. It turns out that one can bypass this restriction of h andh being of the same order, if one assumes a twist gap (as defined in {1.1} we require finite number of conserved currents with dimension less than c/12 to be precise). In that scenario, the results in §3 holds true even if h andh are not of the same order asymptotically. The estimation of the "heavy" region becomes really involved in this case and carried out in details in §3. One has to further separate out the "heavy" zero twist operators and estimate their contribution separately (see the right fig. 6).
We further point out that estimating the integrated density of states require us to prove another lemma, which is special to 2 D Tauberian theorem. We achieve this in §4. The lemma is then fed into the main proof in §5. The result for the integrated density of states also requires an estimate of number of states within an order one area on (h ,h ) plane, which is achieved in §3. In §6, we verify our results using 2D Ising model. The large spin, finite twist is discussed in §7. The large central charge regime is discussed in §8. We conclude with a list of open of problems along with a brief discussion.

O(1) rectangular window
In this section, we study the number of states lying on an order one rectangular/square area of (h ,h ) plane, centered at some large (h,h) and having sides of length 2δ and 2δ. We are interested in h → ∞,h → ∞ limit with δ,δ being fixed order one numbers. We do it in two ways, in the first subsection, we do it generically without any assumption on twist gap while in the next subsection, we assume existence of a twist gap. The assumption of the twist gap (as defined in {1.1}) facilitates studying the regime where h is not of the same order ash.

Generic analysis: with/without Twist gap
Following [2], let us choose functions Φ ± (h,h) such that where Θ h,h,δ,δ (h ,h ) is the indicator function of the rectangle and defined as In principle, one can choose the energy window for the microcanonical ensemble to be of a different shape, for example, a circle. But for now, we consider it to be a rectangle. Now, from eq. (3.1) we obtain 3) Multiplying the above by the density of states ρ(h ,h ) and integrating yields the following inequality: At this point, we use the Fourier transformΦ ± (t,t), defined as This facilitates us to rewrite the inequality in (3.4) as where we have Next we need to split Z = Z L + Z H . In [19], for the mixed temperature analysis, the light sector is chosen to be {(h,h)|h < c 24 orh < c 24 } and the heavy sector is the complement of that. Here, we choose a different light sector {(h,h)|h +h < c 12 } which has finite size (for large central charge, this is not true and one needs to have an extra sparseness condition on the low lying spectra [2,19]). So at least, in principle, we can choose ρ * (h,h) such that ρ * reproduce the contributions from all operators in the light sector. In practice, we take ρ * such that it reproduces only the vacuum state contribution, that is, which has the asymptotic behavior Then we can write e −(β+it+β+it)c/24 Z L 4π 2 β+it , 4π 2 β+it = L ρ * ,L (β + it,β + it). So we get (3.10) Now we can estimate the contribution from the heavy sector using the HKS bound. We choose Φ ± such thatΦ ± (t,t) have finite support [−Λ ± , Λ ± ] for t andt. A possible choice can be made via modifying the functions appearing in [2,3] a little bit. To be concrete, let us make the following choices: where we have We remark that for Φ − , the locus 1 − (a/δ) 2 − (b/δ) 2 = 0 is inside the rectangle region, hence it is a valid choice conforming to the inequality (3.1). Here a = h − h and b =h −h . At this point, our aim is to show that is sub-leading. We will make use of the most basic HKS bound [19] for ∆ in a clever way. We remark that by requiring the saddle of the light sector is located at (h,h), we find β = π c 6h << 1,β = π c 6h << 1, so some terms such as βδ can be dropped from the bounds as it goes to 0 for large h,h. Then using the fact that Φ ± is a bandlimited function, we have Now Z H has contribution from states where either h orh is greater than c/24. Since the contributing states have h +h > c/12, both can not be less than or equal to c 24 . We illustrate the case for h > c/24.
where β * is defined as β * = min β,β , hence we have (for small enough β andβ) Thus we have where the last Z H,∆ = ∆>c/12 e −β(∆−c/12) is the heavy contribution from the original HKS bound for ∆. We will be showing that the above term is subleading. There are two pieces, one with e π 2 c 6β and another with e π 2 c 6β . Let us illustrate the subleading nature of the term with e π 2 c 6β . The other term can be treated similarly. In β,β → 0 limit, we have where τ 2 γ 4 h * (clearly, γ > 1). This can be achieved by choosing Then at large h,h, we have the following bound, starting from (3.10) (the second term i.e. the term with the absolute value in (3.10) has already been shown to be subleading and we rewrite the first term as an integral over h ,h below): We can evaluate this integral by saddle point approximation, where c ± is defined as With the previous choice of Φ ± , we have where we must optimize over 0 < Λ ± < √ 2π γ to get the tightest bound while keeping δ,δ arbitrary.
The condition for the lower bound to be positive is . (3.26) The allowed region increases as we increase Λ − . And the minimum area such that there has to be at least one operator is given by 4δδ = 48γ 2 π 2 = 4.86γ 2 at δ =δ = 2 √ 3γ π . This analysis can be made sensitive to primaries only, thus gives an asymptotic gap between primaries. We suspect that the above does not give the tightest bound for spectral gap (this intuition is coming from the similar analysis done for the spectral gap in ∆, appearing in [2,3])! Next we can keep δ,δ arbitrary and optimize over 0 < Λ ± < √ 2π γ to get the tightest bound.

Analysis of O(1) window assuming a twist gap
We have seen that the result proven in the previous subsection holds true when h andh are of the same order asymptotically. When we make the analysis sensitive to primary, this feature persists. Nonetheless, we can circumnavigate this issue by assuming an existence of twist gap (as defined in {1.1}). One can also do this for the analysis sensitive to all the operators. The only catch is that one has to separately treat the zero twist operators with dimension greater than c/12. We revisit the analysis of suppression of heavy region in the light of the above discussion. At first, we take up the analysis for all the operators. We go back to the eq. (3.16) and redo the first part of the analysis. We separate out Z H into two pieces, one with zero twist heavy operators, we name it Z (0) H , while the other one contains all the heavy operators with non-zero twist, we name it Z (τ ) H . We start with analyzing Z (τ ) H , assuming a twist gap. We have two scenarios.
(3.30) Here going from the first inequality to the second one, we used 2h > g, 2h > g; going from the second one to the third (last) one, we used g ≤ c/12. The rest of the analysis goes in a similar manner, as done subsequently after (3.16) and we deduce: . (3.31) Scenario II: If g ≥ c 12 , we have h > c/24 andh > c/24, leading to (3.32) The rest of the analysis goes in a similar manner and we deduce: (3.34) The subscript on ∆ in the second line denotes the actual conformal weights of the operator and in the penultimate line, we have extended the sum to all the heavy operators. Now one can see the zero twist heavy sector is suppressed as long as we choose Λ ± < √ 2π. Thus, combining everything, we have where ζ 2 = c 12g . We can combine the above to write for all g, The above has immediate implication in terms of asymptotic gap for all the operators, in particular, the gap does not depend on the term γ anymore. Nonetheless, as we already know existence of descendants asymptotically, the asymptotic gap for all the operators is not so illuminating, so we will not illustrate upon this. Rather we come back to this when we make our analysis sensitive to primaries and in that scenario, the result about asymptotic gap is indeed illuminating.
Analysis for primaries: twist gap complementary to asymptotic spectral gap: The analysis for primaries proceeds in a similar manner. We will be estimating the following object exp S Vir where ρ Vir (h ,h ) is the density of primaries. Instead of the partition function we consider the following object (for c > 1, the expansion of this object is universal) (3.38) Under modular transformation, we have Then we define the crossing ρ Vir (3.40) Explicitly, ρ Vir * would be given by the following function: (3.41) The analysis pertaining to the estimation of the heavy part presented before for the analysis of all the operators can be used as a stepping stone for a similar analysis for primaries for c > 1 CFTs. We again use bandlimited functions and we deduce that the support Λ ± has to satisfy 7 : The leading answer comes out to be where c ± is defines as in the Eq. (3.24).
Asymptotic gap: Now we come back to our discussion of asymptotic gap of primaries. We use the function given in Eq. (3.12) but now with constraint as given in Eq. (3.42). Thus the asymptotic binding square will have length 4 √ 3σ π and the binding circle would have radius rσ √ 2 + g with g > 0, and σ, r are given by If we consider tensoring the chiral Monster CFT with its antichiral avatar, we find g = 4 and our result predicts that the asymptotic spectral gap involves a circle of radius 2 √ 6 π irrespective of how infinity is approached. This is above the suspected optimal value 1 (see fig. 4). In a unitary compact CFT without conserved currents, there is a bound on twist gap [20]: In that scenario, we have As a result, we deduce the universal inequality satisfied by the "areal" spectral 7 Should we not assume twist gap, we would have Λ ± < √ 2π γ , just like the analysis for all the operators without assuming twist gap. gap A and twist gap g: where we have shown r = 4 √ 3 π 2.21 > 1 and we suspect that it can be made to 1.

Lemma: density of states on strip of order one width
In this section, we prove a lemma which is going to play a pivotal role in the next section, where we are going to prove an asymptotic result for the integrated density of states i.e number of states upto a large (h,h) threshold. This also helps us to derive the asymptotic spectral gap via strip like regions as defined in {1.4.2}. We start by defining the following functions The aim of this section is to prove the following lemma: Let us focus on the quantity Q, the argument for P follows in a similar manner. In order to estimate Q, we write down the master inequality: where we have used (4.6) Now we separate the contribution to Z 4π 2 β+ıt , 4π 2 β into two pieces Z L (contribution from the "light" sector) and Z H (contribution from the "heavy" sector). The inequality in (4.4) can be written as For notational simplicity, let us name the terms The idea is to show that I 2 ± is subleading with respect to I 1 ± . The argument closely follows the argument presented in §3. Let us concentrate on I 2 ± first. We have We notice that (for the heavy sector, h +h > c/12, thus one of them has to be greater than c/24) . h > c/24 (4.11) Thus we have (4.12) The above analysis is analogous to the one presented in §3. Now we choose , (4.13) and use the HKS bound [19] to estimate Z H,∆ . This leads to the following inequality: (4.14) To estimate I 1 ± , we consider ρ * (h,h), the crossing kernel, defined as Here ρ * (h,h) = ρ * (h)ρ * (h), and ρ * (h), ρ * (h) are given by (4.16) Hence, we have where we have used separability of ρ * in the variable h andh and have defined Comparing the inequalities (4.17) and (4.14), we see that by choosing γΛ ± < √ 2π, with max(h,h) = γ 4 min(h,h); one can make I 2 ± subleading, consequently in the h,h → ∞ limit we have By symmetry we obtain (4.20) The best possible value of c ± (c ± ) can be obtained from [3]. For the verification purpose, here we choose the function given in [2] for estimating Q The above function yields almost the same bound as found in [2], except for the fact that we have to take care of the constraint Λ < √ 2π γ . In particular we find the following bounding function s ± (δ) = log c ± : where a * = 3.38 [2]. We verify the eq. (4.19) in the §6 using the above values of c ± (in particular, see the fig. 11). Similar verification can be done for the eq. (4.20) as well.
The analysis with the assumption of twist gap g proceeds as in the end of §3. We do not repeat the analysis here. We just state the result. In that scenario, one obtains . (4.24) If we make it specific for primaries, c → c − 1 and we have where the second equality follows only if there is no conserved currents because of the bound on twist gap.
We can use this lemma to prove the result about asymptotic spectral gap in terms of strips, as mentioned in {1.4.2}. For this purpose, we use the magic function introduced in [3]. We can let γ = 1 in the above analysis. We have to keep in mind that now the support of the Fourier transform of φ − satisfies Λ < √ 2π γ , thus the minimal value of δ comes out to be 1 √ 2 in stead of 1/2 as in [3]. If g > c/12, then it is possible to redo the heavy estimation using a different kind of separation between light and heavy operators, leading to the constraint Λ < 2π γ , thereby making the value of δ being 1/2. This is in fact the optimal bound for g > c/12 and saturated by tensoring chiral Monster CFT with its antichiral avatar.

The main 2D Tauberian theorem
We prove in this section where τ is the twist of the state with h,h and given by τ = 2min{h,h} and h =h υ with 1/2 < υ < 2 and Υ = max (υ, 1/υ). When Υ = 1, this reduces to the eq (1.1).
In order to prove this we define Since the leading term is already produced by ρ * (h,h), our job is to show that  In order to prove the above, we proceed as in [2] and introduce the following kernel: Here we have done slight abuse of notation. It is implicitly assumed that the function G(ν) depends onβ instead of β. Now it can be shown that [2]: where G ± is defined exactly like in [2] for both the variable ν andν. At this point, we use the kernel given in (5.7) and integrate it against δρ(∆). This yields us the following equation (5.10) Most of the terms can be estimated using techniques from [2]. The new players in the game are the cross terms, for example the term: (5.11) Below we will illustrate how to handle these terms. We remark that this is what requires us to prove the lemma in the previous section. For concreteness, consider the following term θ(−ν)θ(ν)G + (ν) and analyze it carefully. The analysis for the other terms in Z goes exactly in the same manner. In what follows, we will be using the inequalities for β > 0: (5.14) The inequalities (5.13) and (5.14) have been derived in the appendix of [2].
Consider the term For the term with ρ * (h ,h ), the estimation procedure mimics the one presented in the section 5 of [2]. In particular, we have The h integral is done via saddle point method, note it is important to have the factor min (1, (h − h ) −2 ) for the validity of saddle point approximation. This is why, theh integral can not be done using saddle, hence does not produce any polynomial suppression inh. Now consider the term and we divide it into three pieces The estimate for a 1 and a 3 again proceeds like in [2] and we obtain The estimation of the term a 2 would require the lemma from the previous section §4. We subdivide a 2 into three different parts: We have already estimated a 22 in the previous section §4. This is basically order one window. Since, the estimation of a 21 and a 23 proceeds in similar manner, we would demonstrate the estimation of the term a 21 .  The different colors denote the different methods of treating them. As we will see, the region with red dot requires us to further subdivide it. We remark that the corner black ones can be colored blue as well, nonetheless in the main text, we have done the estimation in the "blue" way. Basically, we have to estimate i S i where S i is given by The appearance of 4 is due to the fact that there are 4 terms in the integrand defining S i and each is suppressed in a same manner. Here Now for R i with i = 4, 5, 6 we have (5.30) For R 4 and R 6 we observe the following (5.31) For analyzing the region R 5 , we are required to subdivide it into 9 regions again where each of the region is Cartesian product of order one interval in h andh . Now one can use the lemma proven in the section §3 to show that The estimation for S (2) can be done by saddle point method: We are yet to estimate the following term where y = β orβ. Essentially this is exactly the same argument as in §3. Then by choosing Λ < √ 2π (when h h , otherwise we need to choose γΛ < √ 2π, see the discussion in §3) and using the HKS [19] argument, one can show that the above term is exponentially suppressed compared to the leading answer coming from ρ * (h,h). When h is not of the order ofh, we need to assume existence of twist gap (as defined in {1.1}) and proceed like we did in §3. This concludes our analysis and hence the proof of the main theorem.

Sensitivity of Asymptotics towards spin J
The asymptotic formula given in eq. (1.1) and derived above can be rewritten in terms of dimension ∆ = h +h and spin J = |h −h|: . It turns out that even when ∆ and J is not of the same order, we can do order by order correction to this integrated density of states by spin.
First of all, when J is of order one, we should just ignore J dependence of the eq. (5.40). In fact, J = ∆ 1/n with n > 4/3, one can ignore J dependence. Thus in this regime, we have When 16/15 < n < 8/7 (n can not less than one because of unitarity bound), we have (5.43) We remark that not all the term in the exponential are meaningful, we have to do an J/∆ expansion and the only meaningful terms are of the form exp ∆ − with < 1/4 (since exp ∆ − 1 + ∆ − becomes comparable to error term for < 1/4).
For example, when n = 15 14 , it is meaningful to keep the following terms only: One can generalize the equation (5.43) for J = ∆ 1/n with 2 m+3 2 m+3 −1 < n < 2 m+2 2 m+2 −1 ≤ 8/7 where m ≥ 1 is a fixed integer. We define f (m) = 2 m 2 m −1 and we have : 6 Verification: 2D Ising model We verify the bounds on O(1) correction to entropy associated with order one window centered at some large h andh, proven in §3 as shown in the figure 8, 9 and 10. In §3, we have obtained two kinds of bounds, one without assuming any twist gap while the other one assumes existence of a twist gap. For 2D Ising model, it turns out that the bound coming from assuming a twist gap (indeed 2D Ising model has a twist gap) is stronger than the one without using the information about twist gap. Thus in the figures below, we verify the bounds that uses the information about twist gap. In fact, as we have mentioned in §3, the use of twist gap actually enables us to probe regions where h andh are not of the same order. The figure 10 elucidates such a scenario.
We verify the lemma proven in §4 as shown in the figure 11 and figure 12 using 2 dimensional Ising CFT. The partition function for 2 dimensional Ising CFT is given by We verify the main theorem proven in §5 in fig. 13 and 14. We plot the total number of states upto some (h,h) as a function of h. We have considered two different cases where the asymptote is approached along different curves. We compare it against the asymptotic formula we derive in §5.

Finite twist-Large spin
In this section, our aim is to estimate the finite twist and large spin sector of the density of states. Without loss of generality, we will assume h is finite andh → ∞. We will probe the following quantity U h (h) in the limith → ∞:  Intuitively, it is clear that one can not make both β andβ to approach zero while looking at the partition function, as this would probe the regime h,h → ∞. This suggests that we should letβ → 0 and keep β fixed. We can prove an upper bound on density of states in this sector. Let us write down the following inequality: At this point we choose Again we separate the partition function into two pieces; the light sector, where the contribution comes from two kinds of states: a) all the states with conformal weight (h , 0) with h ≥ 0 and b) the states such that h +h < c 12 , and the heavy sector, which is defined to be the complement of this. We remark the heavy sector does not contain any operator with conformal weight (h , 0). In the usual asymptotic analysis as done previously, the operators with (h , 0) such that h > c 12 is put into the heavy sector, but here we can not do that since β is finite, there is no separation between light or heavy in the sense of HKS [19] i.e. the operators with (h , 0) contribute on equal footing as the operator (0, 0). The upshot of the above discussion is that one can define a kernel ρ ft * ("ft" stands for finite twist) such that where the sum on the right hand side is over all the states with weight (h i , 0) and we have used β = 4π 2 β = 2π. As a result ρ ft * is given by Once we have defined ρ * , we follow the usual method and our next aim is to show that the heavy part has a subleading contribution. The estimation of the heavy part is done depending on whether both h andh is greater than c/24 or one of them is less than c/24. In the later case, we have to assume existence of twist gap. Last but not the least, we have to estimate the contribution from the states with (0,h ) withh > c/12. Following the methods in §3, we can show that they are indeed subleading. The leading answer is then given by Theh integral can be done using saddle point approximation and we obtain where χ 0 (q) is the vacuum character and q = e −β where we have assumed that there is no nontrivial conserved current. Here c + is defined as One can extend this argument to a scenario where we have nontrivial conserved currents, then we would have a sum over characters for all the conserved currents This sum overh is convergent as the absolute value of the sum is bounded above by partition function evaluated at β =β = 2π. Similar result applies to h,h getting swapped. Thus for finite h andh → ∞ limit we have where χh is the character for the conserved current with weight (h, 0). The similar result specific for the primaries with finite h andh → ∞ limit, would read: where the zero twist primaries have weight (h, 0). Here again, the sum overh is convergent as the absolute value of the sum is bounded above by partition function evaluated at β =β = 2π.
CFT without nontrivial zero twist primary: In case the CFT does not have any conserved current, we do not need to worry about (h , 0) operators anymore and we can do much better as it is possible to choose β = 2π and still define ρ ft-Vir * This parallels the analysis for Virasoro primary in §3 (see (3.41)).
In the present case, using the kernel ρ ft-Vir * Now by appropriately choosing β we can recover the "square-root" edge present in the analysis in [7]. The square root edge in h dependence of the density of states should produce a factor of (h − (c − 1)/24) 3/2 . In particular, we choose β = 1 (h− c− 1 24 ) and let h − c 24 to be very small 8 The above result is consistent with the leading result as reported in [4,5,7,8].

Holographic CFTs
Holographic CFTs are the ones characterized by a sparse low lying spectrum and large central charge. The sparseness condition is first derived in [19], then rederived in [2], where it emerges naturally out of the Tauberian formalism. In the context of asymptotic behavior of OPE coefficients in large c CFTs, a stronger sparseness condition appears as elucidated in [29,30]. In this section we will be exploring such CFTs with large central charge and a low lying sparse spectra. In particular, we derive an expression for the density of states in the limit h,h ∼ c → ∞. Following As done in §3 we separate out the contribution to the partition function into two pieces: the light Z L , and the heavy Z H .Here we choose We are required to show that Z H term is sub-leading. Z H term gets contribution whenever h orh is greater than c/24 and ∆ > c/12. It can be shown that Z H contribution is sub-leading (the method is exactly similar as in §3, one has to be careful about e −βc/24−βc/24 factor, since c is not finite in the analysis.). The result of the analysis is summarized below: where γ 4 = * * ≥ 1, * = max ( ,¯ ), * = min ( ,¯ ). The above requires that 1 − γ 2 12 * > 0, that is, * * = ¯ > 1 12 2 .
In fact, it turns out that we will be requiring much more stronger condition on ,¯ : * > where the crossing kernel of the operator with conformal dimension (h L ,h L ) is given by we find that A ∼ O(ρ * (h,h)). (8.18) where the condition β,β > 2π translates to ,¯ > 1/24. Below we will see that we need to have a more stronger condition i.e ,¯ > 1/6. We are left to investigate the region where either of the h L orh L is greater than c/24. This is basically given by the brown region in fig. 15. Without loss of generality, let us look at the region where h L > c/24. Now the crossing would be given by  4. It would be nice to improve on the value of r and possibly prove that r = 1 (the parameter appearing in the asymptotic "areal" spectral gap) either by some suitable choice of magic functions or by better estimate of the heavy sector of the partition function. The naive generalization from [3] would not suffice. So one needs to be more creative. And this might shed light on the twist gap and provide a way to expound on the proposed gap in [8].
Our work should be thought of a part of modular bootstrap program [20,24,28,[33][34][35][38][39][40][41][42][43][44]. On a more general ground, it would be interesting to see whether Tauberian theorems and/or Modular bootstrap program can say anything about the chaotic, irrational CFTs. An approach borrowing ideas from Tauberian techniques and that of extremal functions appearing in [45][46][47][48] might be useful in this regard. Furthermore, for holographic CFTs, we can only achieve a reduced regime of validity of Cardy formula compared to what is reported in [19]. It might be possible to improve our result. Albeit, we remark that that if the twist gap is greater than c/12, it is possible to achieve the regime of validity of Cardy formula as predicted in [19]. We hope to come back to these problems in future.