Derivation of the two Schwarzians effective action for the Sachdev–Ye-Kitaev spectral form factor

The Sachdev–Ye-Kitaev model spectral form factor exhibits absence of information loss, in the form of a ramp and a plateau that are typical in random matrix theory. In a large N collective fields description, the ramp was reproduced by Saad et al. (A semiclassical ramp in SYK and in gravity, arXiv:1806.06840) by replica symmetry breaking saddles. We derive a two sides Schwarzians effective action for fluctuations around the ramp critical saddles, by computing responses to a smeared version of the two replica kinetic kernel. Our result confirms [1], where the form of the action was heuristically guessed by indirect arguments supported by numerical evidences.


Introduction
The Bekenstein-Hawking [2,3] black hole entropy formula in conjunction with the holographic principle [4,5] suggest that a region of space surrounded by a boundary surface of a e-mail: matteo.cardella@unimi.it (corresponding author) finite area, should be described by a finite dimensional Hilbert space. It is not known how such an holographic description works, except then in few notable cases. In the context of AdS/CFT [6], despite a large amounts of results, several relevant questions related to black holes cannot be cast in quantitative terms, due to the complexity of the boundary theory at finite temperature.
A manifestation of the black hole information problem in holographic models [7] appears in the behavior of correlation functions of boundary operators at very large separation times. Finiteness of the entropy demands that a correlation function cannot decay to zero at large times, since this would violate quantum mechanics in the form of quantum information loss. For the sake of illustration, let us consider a thermal two point function for a boundary operatorÔ(t) T r e − β 2 HÔ (t)e − β 2 HÔ (0) where Z (β) is the thermal partition function, the dimension of the Hilbert space is finite and exponentially large in the entropy and we omit to denote dependence ofÔ on spacial coordinates. At early times G β (t) decays exponentially in time as the effect of thermalization, which in the bulk corresponds to the black hole quasi-normal modes relaxation. However, the finite sum in (1) cannot not go to zero at large times, instead it keeps oscillating with an amplitude exponentially small in the entropy [8][9][10]. On the other hand, the above requirement is not satisfied by a semiclassical bulk theory. The reflection coefficient for an incoming wavepacket scattered from the black hole classical horizon diminishes with boundary time, because of the increasing blueshift of the scattered particle, with an increasing penetrating power beyond the horizon of its wave function. From this intuitive argument, one expects the two point function to go to zero in the separation time infinite limit, which is indeed the case. This process is described by black hole quasi normal modes [11]. On the other hand, one expects a departure from classical dynamics for the black hole horizon at times exponentially large in the entropy, when discreteness of the spectrum of the black hole microstates becomes relevant. This effect should conspire to reproduce the expected erratic fluctuations in the two point function. This is a non perturbative quantum gravity effect that is hard to be reproduced. It is interesting to study this problem in the SYK model, as its collective mode describes the dynamics of a boundary graviton in a nearly two dimensional anti de Sitter spacetime (n Ad S 2 ), which accounts for the dynamics of the horizon of a low temperature nearly extremal black hole in four dimensions. The above formulation of the black hole information paradox in an holographic setup can be rephrased in terms of a somehow simpler and more universal correlation function then (1), which does not contain matrix elements of specific boundary operators | n|Ô|m | 2 but still exhibits the same phenomenon [12]. This is the so called spectral form factor While at short times |Z (β, t)| 2 is of the order of the square of the thermal partition function Z (β) 2 , at very long times, of the order of the inverse of the mean energy level spacing, the spectral form factor reaches a limiting value Z (2β) usually called the plateau, due to a cancellation between the off diagonal contributions in the sum (2). A direct computation at large times of the spectral form factor in theories with a well defined gravity dual such as super Yang Mills is currently impossible. However, the Sachdev-Ye-Kitaev model (SYK) [13][14][15][16][17][18][19] offers both a numerical and an analytic handle for studying the spectral form factor. An accurate numerical analysis of the SYK spectral form factor was done in [20], (see also [21] for earlier related work), while a subsequent work [1] explains part of the observed behavior in terms of a large N collective fields description.
SYK is a statistical mechanics model over an ensemble of quantum mechanical many body systems of Majorana fermions, with random couplings of even order q ≥ 4 all to all interactions. A statistical average over an ensemble of quantum systems is not equivalent to a single quantum mechanical model, a fact that may cast doubts upon using SYK for discussing delicate issues like unitarity in a black hole holographic description. Yet, there are qualitative and quantitative features that survive the ensemble disorder average that make SYK an interesting playground for discussing certain quantum gravity issues related to black holes. 1 1 There are also the so called colored tensor models, quantum mechanical models without disorder that exhibit the same diagrammatic as SYK [22]. SYK has a collective fields description that in the large N limit at low temperatures/strong coupling exhibits an interesting quasi conformal behavior, with aspects of a gravity dual. In particular, as an holographic model for nearly extremal black holes it provides an interesting arena for testing various ideas and proposals and sharpen open problems. One of its attractive features is a dominant soft mode dynamics at low temperatures [13,14,23,24], that from a bulk perspective encodes the full black hole gravitational backreaction [25][26][27][28], and, under certain circumstances, allows for a quantum mechanical description of a black hole interior [28][29][30][31][32][33]. One of the exciting recent results, that is understood also from the SYK perspective, is a gravitational description of the quantum teleportation protocol [28,34,35]. This is based on recent observations on how certain double trace deformations that violate the average null energy condition make a wormhole temporarily traversable [36]. In the teleportation protocol those double trace deformations implement the transmission of classical instruction for the quantum protocol between the two boundaries through external space-time and make the wormhole connecting them temporarily traversable.
Another interesting feature is that SYK exhibits quantum chaotic behavior both at short and long times scales. At short times scales, out of time ordered (OTO) four point thermal correlators involving a perturbation on a typical operator saturate [13,14,37] a universal chaos bound [38] for their Liapunov exponent. In AdS/CFT, saturation of the chaos bound in the boundary is understood [39] by a corresponding bulk near horizon scattering dynamics that involves Dray and 't Hooft shock waves [40]. The requirement for the boundary theory to be a fast scrambler [41] avoids black hole quantum cloning in certain gedanken experiments [5,42]. The fast scrambler property is satisfied by boundary theories that are q-local quantum many body systems with all to all interactions [41]. In the SYK case there are no shock waves in the two dimensional bulk and the only bulk gravitational degrees of freedom are boundary modes that undergo classical chaotic dynamics by moving in hyperbolic space [25,28,43,44]. On the other hand, SYK exhibits also chaotic behavior at large time scales, 2 related to the fine details of correlations among its energy eighenvalues. The connected part of the two point correlation function for the SYK spectral density controls the large time behavior of the spectral form factor, a quantity of interest to diagnose the existence of quantum information loss. Numerics on the SYK spectral form factor exhibit the emergence at long time scales 3 of typical behaviour of quantum chaotic systems, in particular random matrix theory (RMT) universality. We refer to [20] for relevant numerical plots and typical values of the SYK spectral form factor, here we would like to mention the Fig. 1 A log-log plot [20] of the SYK spectral form factor for q = 4, N = 34, β J = 5. In red the time plot for a single realization of the disorder, wild erratic oscillations are quite visible in the ramp and the plateau parts of the plot. Oscillations are washed out by disorder average over a large enough number of samples without tough throwing away the main trend. Statistical average (in black) still exhibits an interested lack of information loss in the form of a random matrix universality. The black plot refers to numerics from an ensemble of 90 samples [20] main qualitative features of the time behavior of the spectral form factor. The plot in Fig. 1, starts at early times with a decaying phase from its initial value Z (β) 2 with a characteristic decreasing power law up to a minimum value. This decaying part of the plot is called the slope and the minimum value is called the dip, which occurs at t di p ∼ e N /2 . After the dip, the SFF starts a linear rising behavior called the ramp. The ramp occurs up to a saturation time t plateu ∼ e N where |Z (β, T )| 2 reaches the limiting value Z (2β). The ramp and the plateau are typical of quantum chaotic systems in particular in random matrix theory(RMT). The emergence of RMT universality in itself is not surprising from a Hilbert space perspective, the interest, motivated by the black hole information problem, is to understand the behavior of the spectral form factor by using a large N collective fields description. This problem was discussed in [1] for the SYK model and by the same authors in [45] for Jackiw Teitelboim (JT) gravity [46,47], where a non perturbative completion of JT gravity in terms of RMT is proposed. In [48] the relations between JT gravity and RMT are extended to the case where the boundary theory has time-reversal symmetry and have fermions with or without supersymmetry.
An interesting point that emerges from the numerics on the SYK spectral form factor is that for a single realization of the disorder, the time plot exhibits erratic oscillations in its ramp and plateau regions, (Fig. 1). Oscillations are washed out either by averaging SYK over a large enough ensemble or by taking a suitable time average of the spectral form factor. Although part of the information that corresponds to the erratic wild oscillations is washed out by the disorder average, the main trend that does not exhibit information loss is still present after ensemble average. This provides an interest for a quantitative study of the SYK spectral form factor in terms of large N collective fields, in relation to the black hole information paradox.
It is easy to check that the connected part of the disorder averaged two point function for the analytically continued thermal partition function exhibits a linear ramp for a contribution of the form 2 in the connected component of the two point correlator for the spectral density of states ρ(E). In the above expression brackets in the l.h.s. denote disorder average, while brackets in the r.h.s. denotes statistical correlation. It is indeed the disorder average over a statistical ensemble that provides a non zero connected component Z (β + i T )Z (β − i T ) c to the analytically continued partition function. In a single quantum mechanical model the partition function is just a number and there would be no two points connected contributions whatsoever. The above relation indicates that in order to study the connected component of the spectral form factor, in terms of a functional integral in collective fields, one should look to a two replica system. In fact, disorder average creates interactions between distinct replica, which leads to a connected component for the spectral form factor. Another indication that the non decaying contributions for the SYK spectral form factor may be understood by a two replica systems may come from the ER = EPR conjecture [7,49,50]. Indeed, it turns out [1] that the family of off replica diagonal saddle points responsible for the SYK ramp are obtained with good approximation from a sum over images of SYK correlators on the thermofield double state. On the other hand, from the JT gravity side, the ramp is reproduced by a JT double trumpet instanton [45], a connected Euclidean baby universe that connects two identical Euclidean black holes at the same temperature.
In this paper we derive a time reparametrization soft modes effective action that corrects the SYK critical (replica non-diagonal) saddles of the spectral form factor (3). This effective action governs the late times dynamics for the spectral form factor. We follow a method inspired by [26] to study enhanced responses to the conformal breaking kinetic operator in SYK. This two sides Schwarzian effective action that we obtain, appear also in [1]. In that work, the form of the effective action was guessed by indirect reasoning, supported by numerical evidence [51]. On the other hand, we carry on a constructive analysis by outlining the dominant quantum effects that correct the spectral form factor correlation func-tions at late times. Our results as a byproduct allow to obtain an explicit derivation of the late times effective action. A point of interest of the analysis carried here is a direct observation of when and how non Schwarzian contributions to the effective action become ineffective for the spectral form factor. This kind of analysis has also some interesting points of contact with the recently appeared work [52].
The organization of the paper is the following, in Sect. 2 we go through the construction of a two replica functional integral for the SYK spectral form factor, by going into details in the derivation of the replica non-diagonal conformal saddles. In Sect. 3 we compute a large N approximation for the spectral form factor by steepest descent method, around conformal replica non-diagonal saddles. We then analyse the effects of a smeared version of the two replica kinetic operator and single out that a specific projection of the kinetic operator dominates the late dynamics of the spectral form factor. In section 4 we compute the soft modes effective action for the spectral form factor that governs the late times dynamics. In appendix, we review the method presented in [26] for the case of regular SYK.

Two replica path integral representation for the connected part of the spectral form factor
The Sachdev-Ye-Kitaev model (SYK) [13][14][15][16][17][18][19] is a statistical mechanics model over an ensemble of many body systems given by N (even) Majorana fermions on a complete hypergraph of even order q ≥ 4. Each system of the ensemble has Hamiltonian where the couplings J i 1 ,...,i q are random variables, independently taken from a Gaussian probability distribution with zero mean J i 1 ,...,i q = 0 and variance J 2 The parameter J fixes the characteristic energy scales of the model. In the large N limit, at high temperatures/weak coupling β J << 1, the model is asymptotically free, while at low temperatures/strong coupling β J >> 1, SYK develops an interesting quasi-conformal dynamics. Conformal symmetry is slightly explicitly broken by a dominant O β J N correction, with leading contribution coming from a local Schwarzian effective action. The same pattern of symmetry breaking and a Schwarzian effective action occur in a certain limit of two dimensional Jackiw-Teitelboim (JT) gravity. Moreover, JT gravity describes the classical dynamics of the near horizon region of a low temperature nearly extremal black hole. In the following we mainly focus on the q = 4 case, extensions of the results for generic even q are usually straightforward. 4 We consider the following connected contribution to the SYK disorder averaged analytically continued thermal partition function two point correlator where brackets denote disorder average. As already remarked in the introduction, in order to have a non vanishing connected component we shall look for a system of two replicas on which to perform disorder average. In fact, disorder average creates an interaction among non interacting replica. By following [1] we construct a functional integral representation for Z (β − i T )Z (β + i T ) c in terms of two copies or replica of SYK, SY K L and SY K R . We consider the following two replica representation for the spectral form factor for one particular realization of the disorder where and The parametrization employed in Eq. (7) requires Euclidean time with fermions antiperiodic by period T .
By disorder averaging (7) one finds the following path integral representation where replica indexes I, J = L , R are summed up. The above connected component has to be contrasted with the disconnected part of the two point function Z SY K (β − i T ) Z SY K (β +i T ) , obtained by the standard SYK disorder averaged thermal partition function Z SY K (β) by analytic continuation β → β ± i T . Z SY K (β) can be computed at full quantum level in the Schwarzian approximation, since the Schwarzian path integral over soft modes is one loop exact [54]. This same result for Z SY K (β) can also be obtained by solving a quantum mechanical problem for a particle scattered by a Liouville potential [55,56]. Accurate methods for computing Schwarzian amplitudes at full quantum level are developed in [57,58]. The full quantum answer is given by which gives the following disconnected contribution to the spectral form factor This disconnected contribution reproduces accurately the decaying slope of the spectral form factor in Fig. 1. It is not surprising that this contribution manifests information loss, since it is given by a product of analytic continuations of the SYK thermal partition function. From the two replica system perspective, it is the result of a replica diagonal saddle plus the one loop determinant from fluctuations, that together give the full quantum answer. In contrast, the connected contribution from off diagonal replica saddle is somehow related to a purification of the thermal density matrix, obtained by doubling the system. Indeed it turns out that the connected saddles can be written in terms of an antisymmetrized version of the thermofield double correlators [1], (see Eq. (35) and related discussion).
Concerning the functional integral in (10), before switching from the representation in terms of Majorana fermions to a more convenient description in terms of O(N ) singlets collective fields, let us notice that in any regime where the effects of the bilocal kinetic operator in (10) can be neglected, the eight fermions interaction vertex in (10) is invariant under the following transformation for two independent time reparametrization diffeomorphisms, (aka left and right soft modes), f I (τ I ), I = L , R. In fact, it can be checked easily that this transformation corresponds to a change of integration variables in the double integral interaction term in (10). Therefore in any regime wherê σ I J (τ, τ ) (13) can be neglected, the system develops the time reparametrization symmetry (14), where fermions are primary fields of weight Δ = 1 4 . Let us notice also that fermions appear in the interaction term in (10) as the follow- As a consequence of (14), in any regime where the kinetic operatorσ I J (τ, τ ) (13) can be ignored, the action is invariant under the following reparametrization of the bilocal field In order to study Z (β + i T )Z (β − i T ) c in the large N limit, it is convenient to recast the functional integral (10) in terms of the collective field G I J (τ, τ ) (15) and a corresponding Lagrangian multiplier Σ I J (τ, τ ), by integrating out the Majorana fermion fields. This is achieved by inserting in the path integral (10) the identity where integration over Σ I J is performed along an imaginary direction in field space. By integrating out fermions one finds where Notice that for notational convenience from now on we use lowercase indexes to denote left and right replicas entries, by switching our previous notation I = L , R to i = L , R. This should not be source of confusion, since fermions have been integrated out. We also use the notation J i j = J i J i where which follows from the definitions (9).
The two replica action (19) gives the following saddle point equations. A variation with respect to Σ i j gives where * is the convolution product While a variation with respect to G i j gives We now show that in the large time T regime, for large enough time separation intervals Δτ , an almost conformal regime emerges. Let us consider the Fourier transform of the saddle equation (22) We look at the regime where the iωδ i j term can be neglected, and the equation simplifies to The following ansatz for the self energy does the job in (26) for small enough ω, because of the square root. The square root ansatz (28), solves also (25) for a suitable constant invertible matrix Φ i j , if one take also By plugging (28) and (29) into (25), one has the following equation which is solved by The conformal solution (28), (29), with Φ i j given by (31) is valid in the regime which corresponds to a negligible iωδ i j term in Eq. (26). In order to have a quasi conformal regime (32) it is therefore necessary that Whenever T >> β, the above condition simplifies further to J T >> 1.
Let us summarize what we have just discussed: for large time argument T , the SYK spectral form factor saddle is well approximated by a conformal replica-non-diagonal complex saddle. This approximation breaks down for short enough time separation intervals |τ − τ |, according to (32).
By taking the conformal complex saddle we just calculated, make it periodical by T , and analytically continue it to real time, one finds Let us pause for a moment, in order to comment on how Eq. (34) relates to Ref. [1]. Here we have computed the conformal complex saddle directly from the two replica (conformal) action, (Eq. (19) without the derivative term). On the other hand Ref. [1] makes a broader analysis on the SYK spectral form factor and, for his purpose, employs an auxiliary SYK thermal system at inverse temperature β aux . The introduction of the auxiliary system follows from the observation that in the infinite temperature limit (β → 0) and large time T , the Keldysh contour in the path integral for the SYK spectral form factor is similar to the one for an auxiliary thermal system, the difference being the way the two contours are closed. Saad et al. [1] claims that the spectral form factor and the auxiliary thermal system have similar saddle point equations, since saddle point equations are insensitive to boundary conditions. On the other hand, since the spectral form factor does not depend on β aux , for consistency the SFF action needs to vanish on the auxiliary system saddle, a fact that [1] checks explicitly.
For non vanishing β [1] provides the following approximated saddle point for the spectral form factor where the above correlators need to be made periodic by period T , by summing over images. In the β > 0 case, however the spectral form factor action is no longer vanishing on the saddle (35) and a runaway in the β aux parameter appears [1]. One can check that Eq. (35) is recovered by taking the saddle solution obtained here (28) and (29) and make it periodic by the following rescaled auxiliary thermal system inverse temperaturẽ For illustrative purposes, let us pause for a moment in order to recall the meaning of the parameters β aux and Δ. According to [1], β aux and Δ come from phase space considerations. By using the correspondence with JT gravity in the bulk, [1] expects the phase space for the spectral form factor backgrounds to be two dimensional. Two possible parameters to describe the phase space, are the auxiliary β aux > 0 and the compact coordinate Δ ∈ [−T, T ]. Δ is a relative time shift between the two holographic boundaries clocks. On the other hand, β aux is a uv cutoff for the boundary theory, arising by holographic renormalization.
Let us notice that the saddle (35) corresponds to a complex version of SYK correlators on the auxiliary thermal system at inverse temperatureβ aux (36), computed on the double field thermal state From a bulk gravity perspective, the auxiliary thermal system inverse temperature β aux fixes the uv energy cutoff of the two holographic boundaries in n Ad S 2 . In the n Ad S 2 wormhole, described holographically by the double field thermal state (37), on each of the two thermal Rindler pathces, the metric and the dilaton have the following form In the bulk solution argued to be relevant for the spectral form factor [1], there is a Rindler time periodic identificationt R ∼t R +T andt L ∼t L −T . This is compatible with the asymmetric analytic continuation in the spectral form factor Z (β − i T )Z (β + i T ) c . This Rindler time identifications give rise to a Lorentzian manifold with the topology of a double cone with closed time curves [1]. β aux arises in relating Rindler timet to boundary time t in the following way. By using an holographic renormalization parameter one can relate the boundary proper time t to the bulk Rindler timet from the bulk metric (38) at a large fixed ρ = ±ρ c This gives the following relation between SYK boundary time t and bulk Rindler timet Given, that β aux provides a uv cutoff for the boundary theory, it plays a somehow similar role to the time argument T in the spectral form factor. In fact, T fixes the uv coarse graining energy scale at which correlation between pairs of energy eigenvalues of SYK are probed Saad et al. [1] points out in their analysis that for non vanishing β, there is an instability in the parameter β aux , which appears as a runaway direction. They suggest that such instability can be cured by going at fixed total energy, in a microcanonical description.
In Ref. [1], the form of the correlarators in (35) was also checked numerically, by computing geodesics distances from boundaries points in a Ad S 2 wormhole geometry. The idea is that in the large N limit the bulk theory classicizes and boundary CFT two points correlators reduce to simple functions of the length of the geodesic connecting the pair of points on which the two point function is evaluated.
As we already mentioned, besides β aux , there is at least a second phase space parameter, since phase space is always even dimensional. In the case of pure JT gravity the dimension of phase space is two, (see for example [59] for a detailed account). This missing phase space parameter is the compact parameter Δ ∈ [−T, T ] that appears in the off diagonals saddle correlator G (35) and it is responsible for the linear ramp behavior in the SYK spectral form factor [1]. Δ corresponds to a relative shift between the time coordinate origins on the right and left boundary. It is not a surprise that it appears only on the left-right diagonal correlator in (35), since only there a relative shift on the origins of times coordinates is relevant. The fact that Δ is responsible for the ramp goes as follows [1]. It turns out that the conformal saddles (35) have zero action. In order to compute the functional integral in the large N limit, one has to still integrate over the phase space parameters β aux and Δ. Integration over β aux gives just a constant overall constant to the spectral form factor. On the other hand, Δ integration on [−T, T ] gives a linear T factor which reproduces the ramp. The above discussion was at the conformal level, for the spectral form factor family of conformal saddle points of the critical action. When effects of the two replica bilocal kinetic operatorσ i j (τ, τ ) (13) are taken into account, at lowest order in the perturbative expansion an effective action for the left and right time reparametrization soft modes f L (τ ), f R (τ ) occurs. This effective action at the saddle point level turns out to be independent on T , consistently with the existence of a linear ramp. However, the value of the action in the saddle point is non zero and a runaway direction for the functional integral representation of the spectral form factor arises, toward large β aux . The β aux instability can be cured in the microcanonical ensemble [1], this is a consequence of β aux being actually related to the energy of the Schwarzian modes. On the other hand, β aux might be stabilized in the canonical ensemble, at the price of slightly changing the model, by adding a small non local coupling between left and right soft modes f L (τ ) and f R (τ ) [1]. This is an interesting possibility that deserves further studies. The two boundaries non local coupling term that stabilizes β aux has some formal analogy to non-local double trace interaction terms that make a wormhole temporarily traversable [28,36]. Analysis related to this discussion from the perspective of two boundaries JT gravity is found in [60]. Beyond the Schwarzian saddle point approximation, the one loop determinant receives contributions both from the two time reparametrization soft modes f L (τ ), f R (τ ), and from fluctuations of the phase space parameters δβ aux (t), and δΔ(t) [1]. These latter fluctuations are controlled by an hydrodynamic action and computation of the one loop determinant consistently reproduces the ramp [1].

A large N, large T , approximation for the spectral form factor path integral
The breaking of conformal invariance is due to the presence of the kinetic operatorσ i j = δ i j δ(t − t )∂ t in the action (19). The action can be recast in the form (42) below, where in the last termσ i j (t, t ) appears as a source coupled to G i j (t, t ). This term breaks conformal invariance explicitly, while the previous three terms in (42) are conformal invariant. As it was discussed in the previous section in the part leading to Eq. (34), for large enough T , the spectral form factor saddle is well approximated by the conformal one given in (34), for sufficiently large time separations |Δt| = |t − t |. The breaking of conformal invariance induced by the source term in the action (42) is therefore a uv effect. In the large N limit and large T regime a steepest descent approximation for the spectral form factor path integral can be taken through the conformal saddle. We take a regularized version of the sourcê σ i j = δ i j δ(t − t )∂ t , by replacing the Dirac delta by a smearing function that cuts off time intervals below the scale 1/J . On the other hand, we replace the linear derivative operator ∂ τ by a more generic combination of terms with different scalings, in order to account for non linear effects. By steepest descent, the approximated action depends on the regularized source, in particular, through a term where the four point function two replica ladder kernel emerges, (see Eq. (50) below). By looking at Eq. (50) one expects an enhanced effect along the four point kernel invariant eigenfunction in functional space. Due to the vanishing of the denominator it seems that a divergence occurs. However, the divergence is cured by the regularization of the source. We provide an argument that shows how for large enough T , all eigenfunctions terms in the expansion for the source does not contribute, except that the invariant eigenfunction of the ladder kernel. This argument provides an explicit form for the regularized source in the large T regime. By plugging this form in the last source term in (42), we are finally able to compute the large T effective action for the spectral form factor in the f L (t), f R (t) time reparametrization modes. We start by performing a translation in field space

The action (19) turns into
As it was remarked in the previous section, (see the discussion that leads to Eq. (16)), the first three terms of the above action are invariant under the simultaneous time transforms as a primary bilocal field with weights 1 4 and Σ i j (t 1 , t 2 ) transforms with weights 3 4 . The spectral form factor conformal saddle (34) breaks spontaneously the twofold conformal invariance down to the diagonal SL diag (2, R) subgroup of SL L (2) × SL R (2). On the other hand, the source term in (42) breaks explicitly conformal invariance. In order to study the effects of the kinetic operatorσ i j (t, t ), we will regularize it, by smearing the singular Dirac delta kernel on the account that times shorter than the time scale 1/J cannot be resolved.
In the large N limit, at strong coupling, the path integral for the spectral form factor can be approximated by using the steepest descent method through the conformal saddle (34) where the normalization in front of the fluctuations is chosen for convenience.
In the large N limit, the path integral is well approximated by evaluating the action in the saddle points in the complex plane and by Gaussian integrating fluctuations along the steepest descent direction through the saddle point.
An expansion of the action up to quadratic order in the fluctuations gives where I c is the action evaluated in the conformal saddle, replica indexes i, j = R, L are summed over and f |g denotes the following scalar product in the space of bilocal functions The four indexes integral kernel appearing in (44) is the two replica version of the four point function symmetrized ladder kernel in regular SYK [23]. The latter for generic even coupling q ≥ 4 has the form We now introduce a convenient alternative indexing for the replica entries. We define α = 1, 2, 3, 4 such that α = 1 = L L, α = 2 = R R, α = 3 = L R, α = 4 = RL. The action (44) in this new indexing reads Since fluctuations are of order 1/ √ N , in the large N limit, one can extend the integration contour on the complex plane along the steepest direction to make a full Gaussian integral. By δΣ α integration on finds By then integrating along the steepest descent contour in the δG α complex plane through the saddle G c α , one finds the following approximation at large N for the action of the spectral form factor where s α (t, t) is a rescaled version of the regularized source. So far, by steepest descent method, we obtained the following large N approximation for the spectral form factor Equation (52) allows to compute various correlators by functional deriving w.r.t. regularized source σ α . In particular, the two point function is given by We now discuss the form of the regularized source. Let us start by considering the eigenvalues problem for the four points conformal ladder Kernel Where, h labels the eigenvalues of the Carimir C diag of SL diag (2), the diagonal subgroup of SL L (2) × S R (2) which is the symmetry group that survives the spontaneous breaking induced by the conformal replica non-diagonal saddle (34). The h dependence on the eigenvaluesk SF F (h) is due to the fact thatK c αβ (t 1 , t 2 ; t, t ) commutes with C diag .
Let us expand the source in terms of the ladder kernel eigenfunctions where U (ξ ) is a smooth smearing function of ξ α = log(ζ α (|t − t |)) that regularizes the singular Dirac kernel δ(t − t ) in the two replica kinetic operatorσ i j = δ i j δ(t − t )∂ t . U cuts off the short times intervals region, below the scale related to the SYK energy scale 1/J . It works as follows, U (ξ ) is demanded to vanish for ξ < 0, and for |t − t | > 2T . ζ α (|t − t |) is defined to be the denominator of the α component of the weight h eigenfunction of the two replica four point kernel Since ψ h α (t, t ) is also an eighenfunction of the SL diag (2) Casimir C diag , ζ α (|t − t |) is typically an hyperbolic sine or cosine. One has 0 < ζ α (|t − t |) ≤ 1 for |t − t | < 1 J , which makes ξ α < 0 because of the logarithm and the smearing function U vanishes correspondingly. This implements the impossibility to resolve time intervals below 1/J . On the other hand, for large enough |t − t |, ζ α (|t − t |) grows exponentially in the separation time. This makes the support of the smearing function U (log(ζ α (|t − t |))) to be of order O(T ).
The smearing function is normilized as folllow Let us write the smearing function as a Fourier integral Since U (ξ ) has support Δξ ∼ O(T ),Ũ (η) has support Δη ∼ O 1 T . Therefore, in the large T limit, the Fourier transformŨ (η) becomes very narrow.
SinceŨ (η) has support Δη = O 1 T , in the large T limit, one can expandk SF F (h + iη) up to the first order in etak SF F (h +iη) ∼k SF F (h)+iηk SF F (h). By residues theorem on the η integral, one finds with Since the support ofŨ (η) is Δη ∼ O 1 T , for large enough T , the only non vanishing term in (61) is the η h = 0 one. We haveŨ (0) = 1 and by Eq.
To summarize, we have shown that in the large N limit, for T large enough, the effects of breaking of conformal invariance in the spectral form factor are fully determined by the regularized source along the invariant direction under the four point ladder kernel. In particular, the large T limit of the two point function is given by It follows that at large T , the regularized source that encodes the explicit breaking of conformal invariance has the following short separation times |t − t | → 0 behavior where P is the projector over the two diagonal directions α = 1 = L L and α = 2 = R R in replica indexes. This projection is required, since the source is diagonal in replica indexes. As already remarked, the breaking of conformal invariance is a |t − t | → 0 effect, that's why we focus on the short time separation uv behavior of the quantities of interest. Since s 1 (t, t ) = s L L (t, t ) is expected to have the same form as s 2 (t, t ) = s R R (t, t ), it follows that the source that dominates the large T dynamics needs to be invariant under the operator P −1K c αβ P =K c 11 . By using (34) in (46) whereK c is the regular SYK, q = 4, four point function symmetrized ladder kernel (47) Therefore, the regularized source whose effects dominate the non conformal dynamics for the spectral form factor at large T can be obtained by looking for the invariant function under the four point ladder conformal kernel in regular SYK The spectrum of eigenvalues and eigenfunctions ofK c is knowñ in the following we recall some of their properties [23].
Since the SL(2) Casimir operator C commutes with the four points ladder conformal kernelK c , C eigenfunctions are alsoK c eigenfunctions. C has a spectrum with both a continuum and a discrete component [23]. The continuum component is given by the points on the critical line h s = 1 2 + is, s ∈ R, while the discrete component is given by h n = 2n, n ∈ N >0 . The Casimir C eigenfunctions have the following form [23] ψ h (t 0 , t 1 , Notice that in the t 0 large limit In the eigenvalue equatioñ by SL(2) symmetry, one can fix t 0 = ∞, t 1 = 1 and t 2 = 0 and compute explicitlyk c (h) by the a double integral For the four order coupling q = 4, the result is The discrete component of the spectrum Spec(C) thus is giveñ We are looking for invariant functions under the action of K c . From the above results on the spectrum, we have that k c (2) = 1.
On the other hand, the continuous component of the spectrum does not provide any further invariant eigenfunction, besides.
To summarize, we have found that in the large T regime, the source that induces the breaking of conformal invariance for the spectral form factor has the following form where a 0 is a coefficient to be fitted numerically,

Effective action for the time reparametrization soft modes
By using the |t − t | → 0 behavior of the source (78), we are finally able to compute the effective action for the breaking of conformal invariance in the large T regime, by the last term in (42) The |t −t | → 0 behavior of the T → ∞ dominant source is given by On the other hand, σ R R (t, t ) has the same form (80) but with J R R replacing J L L .
The components of the Green function by the two times reparametrization t i → f i (t i ) transform as where As already remarked in various occasions, the breaking of conformal invariance in the SYK model is due to a uv short-time effect, as the conformal forms for the correlators describe accurately the low energy ir dynamics. We find a local effective local action, by going from (t, t ) to (t + , t − ), slow average time t + = t+t 2 , and t − = t − t . We expand the fields up to lowest order in t − and then we integrate out t − . This gives a local effective action in t + in the two time reparametrization soft modes f i (t + i ), i = L , R.
A Taylor expansion in t − gives On the other hand, the diagonal components (34) of spectral form factor Green functions have the short time separation form By inserting in Eq. (79) the short time t − expansions (83), the rescaled source along the enhanced direction (80) and the conformal Green function (84) one finds In the above expression α S is found by numerical fitting.

Conclusions
In this work we obtain a large N , late times T , steepest descent approximation for the SYK spectral form factor. We worked out a regularization for the replica diagonal kinetic operator and find explicitly its dominant form for the large T regime. By those results, we derive the effective action in the two time reparametrization soft modes, for the large time T regime. Our work puts on a stronger ground [1], where the form of the two Schwarzian effective action was heuristically guessed by indirect arguments, supported by numerical evidences [51]. the perturbation. In (87) the smearing function U , replaces the singular kernel δ(t − t ) of the kinetic SYK operator. U (ξ ) has a smooth compact support and vanishes for ξ = log J |τ − τ | ≤ 0. This implements the uv cutoff for time resolutions less then 1/J . In the finite temperature case, the infrared cutoff is given by the inverse temperature β, and the smearing function U needs to have period β. Moreover, the following normalization condition is imposed β 0 dξ U (ξ ) = 1.