Comments on the Random Thirring Model

The Thirring model with random couplings is a translationally invariant generalisation of the SYK model to 1+1 dimensions, which is tractable in the large N limit. We compute its two point function, at large distances, for any strength of the random coupling. For a given realisation, the couplings contain both irrelevant and relevant marginal operators, but statistically, in the large N limit, the random couplings are overall always marginally irrelevant, in sharp distinction to the usual Thirring model. We show the leading term to the $\beta$ function in conformal perturbation theory, which is quadratic in the couplings, vanishes, while its usually subleading cubic term matches our RG flow.

In this paper, we consider a 1+1 dimensional generalisation of the SYK model, by considering the Thirring model with random couplings. i.e.
with random couplings J ij;kl drawn from a gaussian ensemble with standard deviation J (we also discuss generalisations of this model to higher powers of ν andν). One can use much of the technology developed in the SYK context, though unlike previous generalisations [17][18][19], the couplings are not taken to be spatially disordered. Some further work on SYK can be seen in .
The model is also an example of a conformal field theory -in this case free fermions -which is perturbed by a large set of marginal operators with random coefficients, some of them being relevant and some irrelevant. In the usual lore, couplings along marginally relevant directions increase, and couplings along marginally irrelevant directions decrease, as one flows towards the IR. In this case, however, RG mixes the two and the anomalous dimension of each operator depends on all of the other random couplings. If we fix the scale at which we define the theory, one expects it will eventually flow in the direction of some relevant operators, but the details and rates at which this happens are unclear.
It is also interesting to ask whether the theory is renormalizable, since we are necessarily turning on marginally irrelevant perturbations -but on the other hand they are in the same statistical ensemble as the marginally relevant ones. For the model in (1.2), the two point function cannot be solved completely, but we compute it at large distances (in the sense that we make precise below) for all values of J, without exhibiting a singularity, or the generation of a mass gap, in the physical domain 1 . From it we can extract an effective β function β(J) = 4π 2 J 3 (1. 3) and study the above questions. The β function indicates that the interaction is marginally irrelevant so the theory is not conformal. It will also turn out not to be conformal in the limit of J → ∞. Rather, the theory requires regularization for any value of J, and this regularization breaks conformal invariance already at the level of the 2-pt function. Useful reference points are the Thirring or Gross-Neveu models [42][43][44] (which are equivalent in our case), where the interaction term is non-random and diagonal between the left and right moving fermions. The β function of these models is well understood -in 2D, the leading term is quadratic in the 4-fermion coupling, and the interaction is marginally relevant for one sign of it. The theory then develops a mass gap (which is seen in [44] in perturbation theory via the appearance of a tachyonic pole at very low energies).
At finite N we can expect a similar behavior in the extreme IR for each realisation, since each contains both relevant and irrelevant operators, and the former will eventually grow to be large, most likely triggering a mass gap. This, however, is the equivalent of the statement that in the 0+1 SYK the very late time behavior, or the behavior of states close enough to the ground state, differ from realisation to realisation [29] The large N ensemble average does not capture this expectation, but it captures the flow prior to chiral SO(N ) L × SO(N ) R symmetry breaking, where it is modified relative to the Thirring and Gross-Neveu models -in our models we see no pole in the physical regime. Furthermore, the leading quadratic term in the β function, which drove the RG in the Thirring and Gross-Neveu models, vanishes for the ensemble (in the large N limit). The β function (1.3) is positive for any value of J, and it is driven by the cubic terms in the random couplings of the β function. The latter is non-universal in the Thirring or Gross-Neveu models, but it is universal in our case for schemes preserving the statistical chiral symmetry. In fact, we reproduce the result (1.3) from suitable average considerations of the perturbative conformal field theory results reported in [46,47].
The 2-pt function therefore does not show signs of the chiral phase transition, at leading order in 1/N , for any value of the coupling. It is known that in low dimensions quenched disorder can smooth out thermal phase transitions [48]. Our results generalize that statement to translationally invariant theories, with the chiral phase transition being smoothed out at this order. To our knowledge this is the first demonstration of smoothing out of a phase transition by translationally invariant ensemble average in a field theory.
One final note about the application of this model to black hole physics. The obvious analogy is to consider this model as related to a black hole in AdS 3 . However, we do not want to consider this as the full dual to the entire AdS 3 space because the β function is positive and the theory will not be asymptote to empty AdS 3 . Perhaps a better way to think about this model is as describing the effective interactions of degrees of freedom inside a BTZ black hole [19]. One then needs to 1) "prevent" these degrees of freedom from appearing at high momenta, and 2) couple additional "probe" fields to the model, which will mimic the degrees of freedom outside the black hole 3 . This paper is organised as follows. In section 2, we define the Random Thirring model. In section 3, we compute the fermion 2-pt function by solving the Schwinger-Dyson (SD) equations in the IR, for all values of J, and show how the 2-pt function flows. We also discuss the reparameterization invariance of the model. In section 4, we compare our model to the flow of the usual Thirring model (or the Gross-Neveu model in our case). We compute the β function and show which terms are washed out due to the randomness of the coupling, and which remain active. Section 5 contains some concluding remarks.
During the completion of this work, the preprint [49] appeared which contains a 1+1 generalisation of the SYK model with the same interaction term. Although some of the intermediate steps are similar, the UV starting point of our model is different than the one described there, as is our solution and its implications for the RG flow of the model.

The Random Thirring, or SYK, model in 1+1 dimensions
Consider a 2d non-chiral theory with N left moving ν i and N right movingν i (i = 1, . . . , N ). Its free theory action and partition function in Euclidean space are given by The free propagators are where x denotes the pair z,z. Their Fourier transforms equal For more details on our conventions, see appendix A. The free theory is SO(N ) L × SO(N ) R invariant, with conserved currents defined for a pair [ij] with i < j. Using these currents, we can add an interaction of the form The model we discuss in this work consists in treating the coupling constants J ij;kl as random variables drawn independently from a gaussian distribution defined by where denotes ensemble average over J ij;kl . The resulting model can either be viewed as a translationally invariant continuum generalisation of the SYK model in 1+1 dimensions, or as a random cousin of the massless SO(N ) Thirring model, where by the latter we mean the theory whose interaction Lagrangian involves a fixed diagonal coupling J ij;kl ∼ g δ ki δ lj We will refer to the model in equation (2.1) and (2.5) as the random Thirring model. Though much of our discussion will be about (2.5), there are two natural generalisations to consider. The first involves a higher order interaction term, as already discussed in [49]. The second involves the inclusion of a low-pass filter multiplying each fermion in Fourier space L q,int = dp 1 ..dp q du 1 ..du q δ( where the filter F (k) decays at large momenta. The use of such filters, introduced in this context in [19], might be useful for holography in order to force a distinction between UV and IR degrees of freedom, but clearly it will be very restricted from a local field theory point of view. Our model (2.5) has two types of symmetries : exact and statistical. The first ones hold for any realisation of the couplings, whereas the second ones only arise after carrying the coupling average. Poincaré symmetry and a Z 2 symmetry ν → −ν andν →ν for even q, or its composition for odd q, are exact symmetries of the Lagrangian. Once the ensemble average is carried out, and assuming self-averaging, i.e. that any realisation behaves similarly to the ensemble average, there is an emergent SO(N ) L × SO(N ) R symmetry, rotating left and right Majorana fermions independently, together with a parity exchanging ν i ↔ν i .
As in the case of the SYK model, we will sum over different realisations of the couplings J ij;kl . If the model is self-averaging then this should reproduce, up to small corrections, the results of any specific realisation. However, for a specific realisation the random Thirring model and its q > 2 generalisations are ordinary field theories. For q > 2 the interaction term is irrelevant and we don't expect it to generate new interesting dynamics in the IR (although it can generate one in some intermediate scale). For q = 2 the situation is more interesting, as we perturb the theory by a large number of operators, some of them relevant and some others irrelevant. In fact all the couplings mix together under RG flow but the qualitative picture that some are (marginally) relevant and some are (marginally) irrelevant is expected to hold. Further discussion on these points will appear in section 4.

Two Point Function : Schwinger-Dyson Equations
To derive the Schwinger-Dyson (SD) equations for our models, we use the standard replica method [11] to perform the disorder average. If we further assume the replica symmetry is unbroken, the single replica action acquires the following bi-local interaction term (3.1) Next, we introduce two Lagrange multipliers Σ ν ( x 1 , x 2 ), Σν( x 1 , x 2 ), whose equations of motion impose the constraints This is achieved by adding to the action the extra term The resulting action is quadratic in the fermions. Performing their gaussian integral, we are left with the action: The equations of motion for G ν ( x 1 , x 2 ) and Gν( x 1 , x 2 ) are whereas the equations of motion for Σ ν ( x 1 , x 2 ) and Σν( x 1 , x 2 ) are then: The same procedure can be carried out for the 2q fermion interaction model (2.8). The corresponding SD equations are (3.7)

A comment on reparametrization invariance
Consider the J → ∞ limit of the collective field action (3.4) for the 2q fermion model (2.8). This results in dropping the dependence on the free propagator in the equations of motion (3.7). It is convenient to rewrite the first two equations in real space to discuss their symmetries. As pointed out in [49], these equations of motion appear to be invariant under reparametrization.
For the case q = 2, which we will argue below is the only physically interesting one for us, we can also keep the free kinetic term and reintroduce it into equations (3.8). In this case the equations are invariant under conformal transformation z → f (z) andz →f (z). To see this, consider the 2-pt function transformations The last two equations in (3.7) determine the transformation for the self-energies Σ ν and Σν, and the invariance of equations (3.8) (with the free kinetic term reinstated) requires that (∆ L , ∆ R ) = ( 1 2 , 0) and (∆ L ,∆ R ) = (0, 1 2 ) (for a general q, dropping the kinetic term, invariance under conformal transformations implies only that ∆ L + ∆ R =∆ L +∆ R = 1 q ). The conclusion from this analysis would suggest to solve the SD equations with a scale invariant ansatz for any J for q = 2 (and for large J for q > 2). Below, we will argue this is not necessarily the case. This quantum field theory requires regularization, which breaks conformal invariance. A-priori this symmetry may or may not be restored in the IR (or may be restored upon fine tuning some operators away). Actually, here it will not be restored. We will see this explicitly for the q = 2 theory -the solution would not be scale invariant (although in a mild sense). It is important to emphasize that the regulator needs to be included already in the action, and the breaking of scale invariance is explicit (all the way to the IR) and not spontaneous.
As an aside we would also like to comment that the conformal symmetry of the action is also broken in the IR, beyond the effects of the regulator. The reason is that the action (3.4) should be appended by boundary conditions, or, in the language of the equations of motions, the solutions are legitimate only when the functions decay fast enough at infinity. This is nothing but the usual argument why ordinary 2D CFTs can have only an SL(2)×SL(2) symmetry on the plane -all other generators of the Virasoro algebra introduce singularities at infinity, which change an n-point function into an n-point function in the presence of an operator at infinity. This situation is different from reparametrization invariance in the SYK model in 0+1 dimensions. In that case, at the level of the 2-pt function, the action was consistently reparametrization invariant and a scale invariant solution could be found. Reparametrization was then broken spontaneously. This pattern of symmetry breaking was used heavily in computing the 4-pt function, although in order to make sense of it one had to eventually re-introduce the explicit breaking of reparametrization. Here the cut-off breaks reparametrization explicitly and the solution is not scale invariance. The net result of an explicit breaking of scale invariant (let alone reparametrization) is similar here and in the 0+1 SYK model, but the differences in the intermediate stage suggests that computing the 4-pt function would be considerably more involved.

Solutions of the SD equations -general q
In this subsection we will look for solutions to the SD equations (3.7) for various values of q. We will show that for q > 2 the interaction term is subleading in the IR, and for q = 2 it introduces an unusual form of a logarithmic running into the 2-pt function. We will assume translation invariance and, to evaluate the importance of various terms, we will carry out a scaling argument in the IR. We will also assume that Lorentz symmetry (Euclidean rotation) is unbroken throughout the flow and hence the UV Lorentz numbers of the fermions remain fixed. We will also assume that parity is unbroken. Under these assumptions and one can similarly define Σ(z,z). Then the SD equations (3.7) collapse to Consider a scaling solution with 2(∆ L − ∆ R ) = 1, as dictated by invariance under Euclidean rotations. These solve the SD equations if We treat the q > 2 and q = 2 cases separately.
For q > 2 we obtain ∆ R < 0, which means the solutions are not physical. If the model is self-averaging (which we assumed before) then correlators should be close to being correlators in a unitary field theory. This necessitates that both ∆ L and ∆ R are non-negative. Or, related to it, the free propagator G (0) (p z , pz) ∼ 1 p already scales asλ and is leading at low momenta. Basically this means the theory is free in the IR 4 . This is what one expects from field theory considerations as well. Deforming a free fermion theory by operators of the form ν qνq , which are operators of dimension (q/2, q/2), the perturbation is irrelevant 5 . Thus, we will not consider the q > 2 case further.

Solutions of the SD equations for q=2
The scaling argument (3.12) suggests that the solutions to the SD equation (3.11) for the q = 2 model are for some constant C. The difficulty with this ansatz originates in the Fourier transform of Σ(z,z). The latter is not well defined and requires the introduction of a scale. Following the conventions in appendix A, it is given by One might still be able to work in intermediate energy regimes, where the above approximation might be correct, but it will break down at sufficiently low energies. 5 This is also potentially where we differ from [49] which add the interaction term to an altogether different theory.
We would therefore like to see what the solution to the SD equations is. Our approach will be to modify our ansatz to include some dependence on logarithms, which will turn out to be essential to encode the RG flow of the model. We will present the large J solution in the IR before proceeding to the all-J solution (still at large distances).

The Large J IR solution
We first consider the large J limit, where the free propagator terms in the SD equations can be neglected. If we interpret the presence of the logarithms above as inducing some soft RG flow on top of some leading power law term, as dictated by the scaling argument (3.12), it may appear natural to consider the ansatz This depends on two parameters A and α, to be determined by solving the SD equations, and a cut-off scale Λ, whose role as a UV cut-off will become more apparent below. In this ansatz, we solve the "interaction" SD equation exactly, and then we analyse whether there is any regime of low momentum where our ansatz is self-consistent, i.e., the "propagator" SD equation is solved (at least in some approximation).
Working in momentum space, we observe that the Fourier transforms G(p,p) and Σ(p,p) of our ansatz (3.16) satisfy (3.17) Hence, the integral we need to evaluate is where we used the integral representation of the Bessel function J 0 (x) to perform the integration over the angular variable and |p| 2 = 4pp = p 2 according to our conventions in appendix A.
To evaluate this integral, we start from the convergent formula 6
valid for 0 < < 3 2 and where ψ(x) is the Digamma function. Notice we can compute (3.20) To proceed, we can write G (|p|/Λ) for small as where γ is the Euler's constan. Plugging this leading logarithm into (3.20), and dropping the divergent momentum independent −2 term and the subleading momentum dependent logarithms, we get Removing the 1/ 2 divergences in the procedure above resembles dimensional regularization, but since we do not understand the renormalization of the model well enough we provide an explicit argument. In appendix B, we provide a specific regularisation of our original integral (3.18) leading to the same conclusion in the same IR regime. The gist of the argument is that since Λ is our cut-off, we should not do the integral below a distance Λ −1 . Since the integral is convergent for > 0, we can write Changing variables to y = Λr, working in the limit of small |p|Λ and expanding J 0 (since the range of integration is finite in y), we obtain in the second integral a divergent 1/ 2 piece which cancels the divergence from the first integral, justifying the derivation of equation (3.22). We would like to extend our conclusion for integer α to a general one by the following qualitative argument (which becomes more precise for α's which are negative enough). Working with the dimensionless variable y = Λ r, When |p|/Λ is small, we can split the integral (3.25) into [1, C Λ/|p|] and [C Λ/|p|, ∞), where C is finite and smaller than 1. In the first interval, the argument of the Bessel function is generically small and the function can be approximated by one. The resulting integral immediately reproduces (3.24). We are left to argue the contribution from the second interval is subleading. This is so for two reasons. First, when the argument z of the Bessel function is large, the latter can be approximated by z −1/2 cos(z − π/4) (up to coefficients of order one). The integral involves a non-oscillating function multiplied by this oscillating decaying function. Second, for α negative and large enough in absolute value, the logarithm is small in this range. In any case, the integral can also be evaluated numerically for the specific values of α that interest us in the same regime |p|/Λ 1. These confirm the expression (3.24). In the next subsection we present additional arguments for the validity of this large J solution.
Using (3.24) determines the two Fourier transforms . (3.26) The SD equation G(p,p) Σ(p,p) = −1 requires As a quick consistency check that this is a sensible result, note that A is positive. This is a necessary condition because it has to correspond to a 2-pt function (where ν is an Hermitian operator) which is positive.

Solution for any J
Another regime to consider in coupling space is small J, by carrying out perturbation theory in J. We can actually do better and find an all-J solution, by remaining in the IR, i.e. large Λ 2 zz or Λ 2 /|p| 2 . We will assume an ansatz of the form with a 0 = 1 to match the free theory propagator. Note that this is not the most general expansion in J and in logarithms that we could consider. In fact if we solve (3.5)-(3.6) perturbatively in J, the coefficient of J n in (3.29) will have corrections which are lower powers of log(Λ 2 zz) than indicated. However, at lengths much larger than the cutoff scale, these terms are subleading. We will show that the above ansatz is self-consistent, and also solves the SD equations with the right limits J → 0 and large J, i.e. it re-sums the perturbative series allowing to interpolate between both expansions, as long as the IR criteria above is satisfied.
Instead of plugging our expansions in the SD equations to determine the coefficients a n , b n recursively, we notice that the relevant Fourier transforms, summarized in eq. (A.8-A.9) (whose operation is denoted by F below) can be recast as Here F G , F Σ should be thought of as a Taylor expansion in its arguments. This allows us to write our SD equations (3.5) and (3.6) as where the integration constant is determined by the boundary condition at J = 0, i.e. the absence of self energy (Σ = 0) when the interaction is turned off (J = 0). Taking derivatives, we derive the ODE equation The ODE has a unique solution, giving rise to the 2-pt function G(z,z) = 1 z

(3.33)
The latter interpolates from small to large J, and has the correct limit at J = 0 and large J.

Beta Function
Our analysis confirms the existence of an RG flow. We can extract the beta function β(J) for the coupling J that controls our random Thirring model from the Callan-Symanzik equation satisfied by the fermion propagator G(z,z) whose solution is given by A positive β function tells us the effective coupling becomes smaller in the IR. We are perturbing the action by a large number of marginal operators, which have non trivial 3-pt functions between them. This means that some of the coefficients of these operators will grow and some will decrease as we flow to the IR. In fact, which ones grow and which one decrease change as the couplings themselves evolve. If we think about J as measuring the root mean square of the couplings, then some of them clearly decrease, generating a positive β function. Still, the result is a bit peculiar if we want to think about J → ∞, at least at finite N , which is not the case here. If we have a theory at finite cut-off and we turn on some relevant and irrelevant operators, then eventually some relevant directions will grow (and the irrelevant ones will decrease). This means that we expect that eventually the β function will turn negative -we see no traces of this, presumably because of the large N scaling taken.
In the next section, we will gain a better understanding for the positivity of the beta function (3.35) by using conformal perturbation theory around the free field fixed point.

The β function in conformal perturbation theory
To obtain a better understanding of the RG flow, we compare our theory to the Gross-Neveu [44] and Thirring models (for reviews, see [42,43]). These have been studied extensively, also in the context of the large N limit. Furthermore, since our model also deforms a free theory by 4-fermion marginal operators, conformal perturbation theory remains a very useful tool for any specific realisation of the couplings J ij,kl , when these are small. The main difference is that our couplings are as random as possible (although they have a large statistical symmetry) and the key issue is the interplay between randomness and conformal perturbation theory.
In this section, we show that even though the β function is quadratic in the couplings for a given realisation of the couplings J ij,kl , this flow is suppressed at large N for the ensemble average in our models. We will argue that contrary to what happens to most known models, the next order contribution, i.e the cubic coupling term in the β function, is universal 7 , drives the RG flow discovered in the previous section and in fact agrees with it for all values of J, and not just at weak coupling, in the |p|/Λ 1 regime.

RG in the Thirring and Gross-Neveu models
As discussed in section 2, our model is the "disordered" version of the Thirring model, given in equation (2.7). The latter is usually treated by introducing a gauge field A a µ and rewriting the interaction as where we switched to a vector notation for the original currents and suppressed gauge indices [42,43]. The A 2 µ coupling is driven g → ∞ in the IR. For the case of a single fermion in the fundamental representation, the theory is driven to a massive theory. If the starting UV field theory is a higher level WZW model, then it is expected 8 to flow in the IR to a coset theory where we quotient by SO(N ).
In our case, it is actually more useful to compare to the Gross-Neveu model [44]. For Majorana fermions, the Thirring model can be written as Gross-Neveu model by rewriting Following the notation of Dirac fermions in Minkowski space used in [44], this model is It is convenient to introduce an auxiliary bosonic field σ because at large N the fermion propagator is not renormalized and the boson propagator becomes .

(4.5)
One can use the Callan-Symanzik equation to extract the β function at large N to be Hence, if g > 0, the coupling is marginally relevant and the theory is asymptotically free in the UV. If g < 0, the coupling is marginally irrelevant and the theory becomes free in the IR. We expect a similar structure to exist in the random Thirring model -some of the couplings will be relevant and some irrelevant, their nature changing due to the nonlinearity of the β functions (as we will see below). In this situation it is more convenient to think about the theory as defined at some cut-off scale which is held fixed, in which case, at least for finite N , we expect the theory to emerge along some relevant directions and eventually develop a mass gap. However, this may not be accessible in the large N , in the phase where the statistical SO(N ) L × SO(N ) R is preserved.
For completeness, in the Gross-Neveu model the dynamical scale is given by where Λ is a UV cutoff used to regularize the theory. In the IR, the theory breaks the chiral symmetry and flows to a massive theory. In fact, the theory is integrable and the complete spectrum of particles, labelled by n (appearing in full multiplets of the symmetry) can be computed exactly and it is given by .
β function from conformal perturbation theory. The RG flows for the Thirring and Gross-Neveu models can also be understood using conformal perturbation theory. Consider a 2d CFT (the free fermion theory in our case) deformed by some marginal operators O α with conformal dimensions (h,h) = (1, 1) Conformal perturbation theory determines the first contribution to the β function of these couplings to be 9 [51] β when we normalise the 2-pt functions canonically. In our discussion, the deformation is of the form J a andJā are currents labelled by some index in the adjoint representation of SO(N ) (see appendix C for more details on these conventions). (the same letter J is used to denote both the currents and the coupling constants, but they will be distinguished by their index structure). The operators that we are perturbing by are O aā (z,z) = J a (z)Jā(z). Their three point functions are of the form C aā,bb,cc ∝ f abc fābc . For the Gross-Neveu model, the couplings equal J a;ā = δ aā g √ N and f abc are the structure constants of SO(N ) given in appendix C. Using the identity (C.4), the beta function (4.13) reduces to β(g) ∼ −g 2 , (4.14) reproducing the behaviour in (4.6).

β function in the random Thirring model
The standard use of conformal perturbation theory is for specific realisations of the set of couplings J a;ā . In the following, we explore its extension to 2d CFT with disordered marginal deformations.
Flow of ensemble : The 0+1 SYK model, the random Thirring model, or any of their cousins, have an infinite number of couplings in a given realisation. In principle we should track all of them. Once we declare that couplings are drawn from some class of distributions f ({J a;ā }), e.g. Gaussian in our case, the number of parameters is greatly reduced. However, the distribution may flow, and its functional form may change with scale, i.e the ensemble may become a function of scale f ({J a;ā }, µ), in effect introducing more parameters. There are two equivalent ways of thinking of the RG flow in disordered theories. The first is to think of the flow of ensemble of couplings. In this, the ensemble average is done with a scale dependent distribution denoted by µ . For example n α=1 The second is to consider the flow of couplings for a specific realisation drawn from a fixed ensemble. This we denote by the usual . For example Given the beta function of couplings in this way, one can deduce the corresponding flow of ensemble by demanding that (4.15) be the same as (4.16) which is equivalent to where µ is the RG scale. We can characterize these distributions by their moments, and follow their RG flows by the infinite set of β functions We will refer to them as the distribution β functions. Working with the distribution β functions might be a bit counter-intuitive at times. For example, the simplest case is when n = 1. Using the leading β function, quadratic in J a;ā , in (4.13) and performing the ensemble average with J b;b J c;c ∼ J 2 N 3 δ bc δ b c , the average β function equals However, this does not mean the theory does not flow.
In general, if we choose some functional form for the distribution, the RG flow can take us out of this subspace of functions. For example, working at finite N in our model, the distribution does not remain Gaussian, or even symmetric under J a;ā → −J a;ā , along the RG flow. To see this, consider (4.18) for n = 3 (4.20) Hence, β 3 = 0 at finite N , which is not compatible with a gaussian distribution. We observe, though, that this departure from gaussianity is subleading in N , compared to the nominal scaling in which J a,ā ∼ J/N 3/2 .
Flow of the couplings: The various β n are statistical averages of various couplings in the theory. In principle, all of them appear in different Callan-Symanzik equations and one needs to track the RG flow for all of them. At large N , and as long as the statistical symmetry is unbroken, some simplifications occur. Here, we will focus on some additional aspects of the flow generated by the leading term in the β function given in (4.13). We will see this term gives vanishing contributions (in the large N limit) to some of the interesting β n 's.
• First, the ensemble average of the β function was already computed in (4.19) Note this result also follows from symmetry -the β function transforms in the ad(SO(N ) L ) × ad(SO(N ) R ) representation and SO(N ) L × SO(N ) R is restored after the ensemble average. Hence the averaged β function must vanish.
• Although we find that the averaged β function vanishes, it could be that there are large fluctuations around this mean value. To gain some intuition on what the physics may be for a given realisation, we compute the standard deviation of this averaged beta function. This is given by β 2 a,a (not summed over a, a ) bb cc bb ĉĉ The associated nominal scaling is and hence to leading order in the large N expansion (with J held fixed), the leading term in the β function should be taken to be zero.
• If the ensemble remains gaussian to leading order in N , the only parameter characterising the ensemble is J (and in any case, the standard deviation is the leading statistical moment). Hence it is interesting to study its flow by studying At this order of the β function, i.e. quadratic order in the couplings, β 2 also vanishes i.e., These are actually different measures why the leading quadratic term in the β function (4.13) does not contribute to the flow, due to cancellation between the different couplings. These arguments are correct in perturbation theory, or as long as the SO(N ) L × SO(N ) R is unbroken. They may be limited in teaching us about global issues of the RG flow because the latter is non-linear in the parameters. For example, if we fix the couplings at some high scale Λ, then since relevant couplings are turned on, then eventually an IR dynamical scale would be generated (at least at finite N). Computing this dynamical scale for a specific realisation, as a function of J a;ā , and then averaging over the J a;ā is different from first averaging the β function and then trying to deduce the flow's end.

β function from conformal perturbation theory
Conformal perturbation theory provides an expansion of the β function in higher powers of the couplings J a;ā . We saw the quadratic terms are effectively zero at large N in our models. Yet the distribution parameter J 2 still flows as we see explicitly in our computations (3.35).
In this section, we argue this RG flow originates from the cubic terms in the β function. Before starting this discussion, let us comment on the status of these cubic terms, which are usually scheme dependent and hence non-universal. For example, if we redefine the couplings appearing in (4.10) bỹ the higher order terms in the β function expansion will change. However, we can argue that such terms are universal in schemes preserving the SO(N ) L ×SO(N ) R statistical symmetry. To see this, notice that the redefinition (4.27) would change the cubic term in the beta function as follows 10 where C αβγ are given in (4.12). If the redefinition preserves the SO(N ) L × SO(N ) R , the matrix A has to intertwine SO(N ) The only option is to take A αβγ ∝ C αβγ but then the cubic term in the β function does not change.
In the following, we work with such a set of schemes, dictated by symmetry and having a universal cubic term contribution to the β function. We use the results in [46,47]. The formula for the β function evaluated to this order is given in equations (5.12) and (5.13) of [47] which when adapted to our notation gives 11 10 Hereλ ≡ µ∂µλ. 11 Our convention for the β function has an extra sign compared to that of [47]. In what follows, we will raise and lower adjoint indices a, b etc by δ ab , δ ab respectively.
Flow of ensemble : Finally, it is worth noting that this cubic term preserves the Gaussianity of the distribution of the couplings at large N . To prove this we can show that the RG flow of any statistical average of 2n J's is given by taking into account, at leading order in N, by the RG flow of each pair. i.e., Gaussianity implies that where the permutation σ ⊂ S 2n has n 2-cycles -basically all possible Wick contractions. To see this, we evaluate the left hand side first The evaluation of the right hand side has two kinds of terms. One in which two of the J b i ;b i contract among themselves. Notice these are exactly the terms captured by (4.36). Their nominal scaling is N −3n . The second kind involves no contractions among the J b i ;b i . Consider one such term below (where say J b i ;b i contracts with J a i ;ā i for i = 1, 2, 3 and k > 3). We show it has subleading nominal scaling. Hence, it is negligible.
Since r fr a kā1 frā 2ā3 does not scale with N , we conclude the nominal scaling of this term is N −3n−3 , which is subleading compared to the N −3n captured by the gaussianity preserving term in (4.36).

Discussion and Future Directions
We computed the 2-pt function for a set of N 1+1 dimensional Majorana fermions with a random Thirring interaction. We obtained that the β function is positive, which means that the theory is not renormalizable. The running is, however, weaker than the usual logarithmic running of marginal operators.
As an effective theory below some scale, it is perfectly valid, and its RG exhibits some new features. This is the context in which the theory might be useful for black holes in AdS -the model above might serve as a heuristic model for the degrees of freedom inside a black hole at some finite energy, and if we want to understand its surrounding weakly curved AdS outside the horizon, we need to couple it to additional "probe" degrees of freedom [19]. The energy scale is therefore set by the temperature and the degrees of freedom of the black hole "never make it" to high energies.
There are some natural future directions 1. A better understanding of the solution to the Schwinger-Dyson equation: We treated the Schwinger-Dyson equations in an approximation in which the we solve the spatial "interaction" equation exactly, and the momentum "propagator" equation only approximately. This was convenient since the "interaction" equation is more nonlinear. The fact that we were able to approximately solve the equations consistently for all values of J for sufficiently low energies lends support to this approximation scheme. However, it would be worthwhile to see if the approximation can be justified, or tested on other examples.
2. Thermal partition function: One natural extension is to compute the thermal partition function. We expect that we will be able to do it at low temperature, where we can use a long distance approximations of the type that we used above.
3. 4-pt function: In the 0+1 dimensional SYK there was a truncation of the action, in terms of 2-pt functions, to a reparametrization invariant action, and there was no need to regulate the action in order to make sense of the 2-pt function. This enabled partial computation of the 4-pt function, although the final expression required a scale (or the re-introduction of the RG flow). Here we need to include the RG flow already at the level of the 2-pt function, which means the computation of the 4-pt function might be considerably more difficult. It is, however, the main indication whether the theory can be related to a bulk theory of any sort.
4. Renormalizability and stability: We would like to see whether the theory is renormalizable, and at what scales it develops a mass gap in the IR. The former might be easier in the bi-local description of the mode. To check for a mass gap one can look for poles in futher correlation functions. In the Gross-Neveu model, for example, one indication for the mass scale is a tachyonic pole in the 2-pt function for σ, which feeds into the a pole in the 4-pt function of the fermions (there is no pole in the 2-pt function of the fermion which is not renormalized when writing the theory with the σ field). To see this effect in our model, we again need to compute the 4-pt function of the fermions in the model.

The
Thirring model is dual to Sine-Gordon models via bosonization dualities. This suggests that one can obtain an SYK like models with bosonic fields too, atleast in 2 dimensions. It is interesting to then explore if there are other interesting bosonic SYK models.
6. Finally, it would be interesting to couple this model, at finite temperature, to additional "probe" degrees of freedom and examine whether it resembles a black hole in AdS 3 .
We hope to return to these issues in future work.

A 2d Euclidean CFT conventions
Given the Euclidean action S[ν,ν] for the Majorana fermions ν,ν, the euclidean partition function is given by Complex coordinates are defined as and similarly for momenta One can check that zz = | x| 2 and |p| 2 ≡ 4pp = | p| 2 p. x = pz +pz , δ 2 ( x) = 2δ(z)δ(z) , δ 2 ( p) = δ(p)δ(p) 2 (A.5) The definition of Fourier transform in complex coordinates translates to  Since the integral diverges at small r, we need to cut it off. It is natural to do so at r = 1/Λ (or larger values so that the logarithm does not become negative). This is another indication that the scale Λ behaves like a UV cut-off. The resulting integral reduces to F α (|p|/Λ) = π 2 α+1 where we introduced the dimensionless variable y = Λ r. Let us consider the α = n ∈ Z to get some intuition. We are interested in studying the momentum dependence of this integral. For this purpose, define τ ≡ |p|/Λ and t = τ y, so that (B.2) equals Its derivative with respect to τ satisfies Notice the contribution from the lower limit of integration does not contribute since the integrand vanishes at t = τ due to the logarithm. To proceed, we use an integral appearing in [52] I 0 (τ ) = This result is valid for any τ . In the following, we will explore the physical regime corresponding to the IR, τ ≡ |p|/Λ 1, where the dominant contribution to I 0 (τ ) is captured by the logarithm and all the analytic terms are dropped.
It is easy to solve the recurrence relation (B.3) when focusing on this dominant logarithmic contribution. Indeed, by induction, we can prove that the ansatz solves (B.3). Notice that the integer nature of α was used in the induction step, so that I 0 (τ ) belongs to our series. The above derivation for α ∈ Z also determines the logarithmic nature of the subleading contributions F α (|p|/Λ) = π α + 1 log(Λ 2 /|p| 2 ) α+1 (1 + O(1/ log(|p|/Λ))) , (B. 6) in agreement with the result (3.22) in the main text.

C Normalisation of currents
When matching the free fermionic notation with the SO(N ) at level one, the currents are written in terms of adjoint indices a J a (z) = 1 2 i,j ν i t a ij ν j (z) , (C.1) with a similar expression forJ b (z) in terms ofν k , where a stands for a pair of integers (kl) satisfying 1 ≤ k < l ≤ N and the matrices t a ij satisfy 12 t a ij ≡ i δ k i δ l j − δ k j δ l i Tr t a t b = 2δ ab ≡ 2 δ ki δ lj − δ kj δ li a t a ij t a kl = −δ ik δ jl + δ jk δ il , where in the second line a = (ij) and b = (kl). Notice a δ aa = N 2 2 (1 − N −1 ) due to the range of the indices (kl). These matrices have a commutator t a , t b = c if abc t c whose structure constants are explicitly given by It is easy to show these structure constants satisfy the identity Finally, the OPE between currents is with an analogous expression for the opposite chirality. These determine the 3-pt functions (4.12).