Harold Erbin
email: harold.erbin@gmail.com

Vincent Lahoche
email: vincent.lahoche@cea.fr

Mohamed Tamaazousti
email: mohamed.tamaazousti@cea.fr

Constructive expansion for vector field theories I. Quartic models in low dimensions

HAL is

I. INTRODUCTION

Constructive quantum field theory (CQFT) is a branch of mathematical physics that aims to investigate whether quantum or statistical field theory models are well-defined in a rigorous mathematical sense. Despite their experimental successes, quantum field theories (QFTs) may be badly defined mathematically. This defect was pointed out since their origin almost a century ago. Historically, a vast amount of QFT axioms have been designed since 1950 [START_REF] Bardakci | Nature of the axioms of relativistic quantum field theory I[END_REF][START_REF] Summers | A perspective on constructive quantum field theory[END_REF] , the most popular and minimal being the Wightman axioms [START_REF] Streater | PCT, Spin and Statistics, and All That[END_REF][START_REF] Flato | Wightman formulation for the quantization of the gravitational field[END_REF][START_REF] Glimm | The Wightman Axioms and the Mass Gap for the P (Phi) in Two-Dimensions Quantum Field Theory[END_REF] . The proof of Wightman axioms for QFT models is reputed to be a difficult task, the difficulty increasing generally with the realism of the considered field theory. All these pioneer works took place in the operator formalism. The development of Feynman's functional integral provided an important change at the beginning of 1970 [START_REF] Glaser | On the equivalence of the Euclidean and Wightman formulation of field theory[END_REF][START_REF] Fröhlich | Schwinger functions and their generating functionals I[END_REF] . In Euclidean time, the connection between Wightman and Schwinger functions, a consequence (for instance) of the Osterwalder-Schrader (OS) analytic continuation procedure [START_REF] Glaser | On the equivalence of the Euclidean and Wightman formulation of field theory[END_REF][START_REF] Osterwalder | Axioms for Euclidean Green's Functions[END_REF][START_REF] Benfatto | Functional Integral Construction of the Thirring model: Axioms verification and massless limit[END_REF][START_REF] Glimm | A note on eflection positivity[END_REF][START_REF] Glimm | A Tutorial Course in Constructive Field Theory[END_REF] , leads to the well-known OS axiomatization. These axioms provide a set of minimal conditions (OS positivity, clustering, ergodicity, etc.) ensuring the possibility to analytically continue from Euclidean Schwinger functions to Wightman functions satisfying all the Wightman conditions. Nowadays, a large literature exists on this subject, especially on the connections between the Euclidean and Minkowskian theories.

At this stage of history, a new type of question emerged in CQFT: "Can we obtain a rigorous definition of quantities like Schwinger functions for interacting models?" From a perturbative point of view, the expansion usually fails to be well-defined, and even for models without divergences or renormalizable [START_REF] Rivasseau | From Perturbative to Constructive Renormalization[END_REF] . For instance, the free energy F may be formally computed in terms of amplitudes A G indexed by Feynman diagrams G:

F = G connected A G . (I.1)
However, in general, this sum fails to converge G |A G | = ∞ because the number of connected Feynman diagrams grows very fast with the number of vertices [START_REF] Dyson | Divergence of Perturbation Theory in Quantum Electrodynamics[END_REF] -a bad consequence coming from the illegal permutation between sums and integrals in the perturbative expansion. The central success of constructive theory on this question was to realize that, using sophisticated methods such as Mayer expansion and lattice clustering [START_REF] Brydges | Mayer Expansions and the Hamilton-Jacobi Equation[END_REF][START_REF] Sokal | An Improvement of Watson's Theorem on Borel Summability[END_REF] , the existence of a non-empty radius of convergence and Borel summability of the perturbative series can gen-erally be proven. The loop vertex expansion (LVE) is a recent constructive technique [START_REF] Rivasseau | Constructive Matrix Theory[END_REF][START_REF] Magnen | Constructive φ 4 field theory without tears[END_REF][START_REF] Gurau | Tree Quantum Field Theory[END_REF][START_REF] Rivasseau | How to Resum Feynman Graphs[END_REF][START_REF] Gurau | Renormalization: An Advanced Overview[END_REF] :

starting with a Hubbard-Stratonovich transformation [START_REF] Stratonovich | On a Method of Calculating Quantum Distribution Functions[END_REF][START_REF] Hubbard | Calculation of Partition Functions[END_REF] , it replaces the perturbative expansion in terms of Feynman graphs by an expansion indexed by trees, whose number increases very slowly with the number of vertices in comparison [START_REF] Rivasseau | How to Resum Feynman Graphs[END_REF] . The strategy may be summarized as follows: for any pair made out of a connected Feynman diagram G and a spanning tree T ⊂ G, there exists a universal normalized non-trivial weight (G, T ) such that:

T ⊂G (G, T ) = 1 , (I.2)
which allows to rewrite the original sum (I.1) as a formal sum over trees

T A T := G T ⊂G (G, T )A G , (I.3)
and, in the best case, to prove that T |A T | < ∞ inside a non-empty domain. In more details, the general strategy for proving Borel summability with LVE uses the Nevanlinna-Sokal theorem, which requires essentially two conditions:

• the existence of a finite radius of convergence;

• a bound for the remainder in the Taylor expansion.

Since its beginning, the LVE has been successfully applied to a vast amount of problems, in particular, to scalar and matrix models [START_REF] Rivasseau | Constructive Matrix Theory[END_REF][START_REF] Magnen | Constructive φ 4 field theory without tears[END_REF][START_REF] Gurau | Tree Quantum Field Theory[END_REF][START_REF] Rivasseau | How to Resum Feynman Graphs[END_REF][START_REF] Gurau | Renormalization: An Advanced Overview[END_REF][START_REF] Rivasseau | Loop Vertex Expansion for Φ 2k Theory in Zero Dimension[END_REF][START_REF] Lionni | Note on the Intermediate Field Representation of Φ 2k Theory in Zero Dimension[END_REF][START_REF] Gurau | Analyticity Results for the Cumulants in a Random Matrix Model[END_REF] , non-commutative field theories [START_REF] Wang | Constructive Renormalization of the 2-Dimensional Grosse-Wulkenhaar Model[END_REF][START_REF] Rivasseau | Corrected Loop Vertex Expansion for Φ 4 2 Theory[END_REF] , and, more recently, tensor models and tensor field theories, and discrete gravity models [START_REF] Gurau | The 1/N Expansion of Tensor Models Beyond Perturbation Theory[END_REF][START_REF] Delepouve | Constructive Tensor Field Theory: The T 4 3 Model[END_REF][START_REF] Delepouve | Universality and Borel Summability of Arbitrary Quartic Tensor Models[END_REF][START_REF] Lahoche | Constructive Tensorial Group Field Theory I: The U (1)-T 4 3 Model[END_REF][START_REF] Lahoche | Constructive Tensorial Group Field Theory II: The U (1)-T 4 4 Model[END_REF][START_REF] Rivasseau | Constructive Tensor Field Theory[END_REF][START_REF] Geloun | A Renormalizable SYK-Type Tensor Field Theory[END_REF] .

In particular, a non-trivial improvement over the orthodox LVE, called multi-scale loop vertex expansion (MLVE), has been introduced [START_REF] Gurau | The Multiscale Loop Vertex Expansion[END_REF][START_REF] Abdesselam | Trees, Forests and Jungles: A Botanical Garden for Cluster Expansions[END_REF][START_REF] Abdesselam | Explicit Fermionic Tree Expansions[END_REF] to deal with models having a finite set of divergent graphs (see specifically [START_REF] Delepouve | Constructive Tensor Field Theory: The T 4 3 Model[END_REF][START_REF] Lahoche | Constructive Tensorial Group Field Theory II: The U (1)-T 4 4 Model[END_REF][START_REF] Rivasseau | Constructive tensor field theory: The T 4 4 model[END_REF] ). Finally, the LVE and MLVE do not seem to be the end of the history: a recent paper [START_REF] Rivasseau | Loop Vertex Expansion for Higher Order Interactions[END_REF] proposed a generalization called loop vertex representation (LVR) which allows to deal with higher-order interactions and which has been successfully applied to random matrix theories [START_REF] Krajewski | Constructive Matrix Theory for Higher Order Interaction[END_REF][START_REF] Krajewski | Constructive Matrix Theory for Higher Order Interaction II: Hermitian and Real Symmetric Cases[END_REF][START_REF] Krajewski | Constructive Matrix Theory for Higher Order Interaction II: Hermitian and Real Symmetric Cases[END_REF] .

This paper is the first of a series aiming to prove Borel summability theorems for a large set of Euclidean vector models using specifically the LVE and its generalizations (MLVE, LVR, and so on). These models and the results we obtain are relevant for machine learning [START_REF] Lin | Why Does Deep and Cheap Learning Work so Well?[END_REF][START_REF] Bradde | PCA Meets RG[END_REF] , condensed matter [START_REF] Weinberg | The Quantum Theory of Fields[END_REF] and SYK models [START_REF] Kitaev | A Simple Model of Quantum Holography[END_REF][START_REF] Kitaev | A Simple Model of Quantum Holography (Part II)[END_REF][START_REF] Maldacena | Comments on the Sachdev-Ye-Kitaev Model[END_REF][START_REF] Polchinski | The Spectrum in the Sachdev-Ye-Kitaev Model[END_REF] (see [START_REF] Rosenhaus | An Introduction to the SYK Model[END_REF][START_REF] Sárosi | AdS 2 holography and the SYK model[END_REF] for recent reviews), (toy models of) string field theory [START_REF] Lacroix | Closed Superstring Field Theory and Its Applications[END_REF] , and multi-Higgs models in particle physics phenomenology [START_REF] Martin | A Supersymmetry Primer[END_REF][START_REF] Binétruy | Supersymmetry[END_REF][START_REF] Bento | Multi-Higgs Doublet Models: Physical Parametrization, Sum Rules and Unitarity Bounds[END_REF] . For this first work, we focus on the LVE method and general quartic just-renormalizable models with Bose and Fermi statistics. The LVE provides a unified framework to study models having Schrödinger (non-relativistic) and Klein-Gordon / Dirac (relativistic) types of propagators, and arbitrary interactions such that our analysis exhaust all the possible cases of renormalizable models with a purely quartic interaction (up to constraints arising from the stability of the theory). We adopt a presentation where the difficulty is increased step by step with the spatial dimension, the UV divergences becoming more and more difficult to control. For this paper, we focus on models having a finite number of UV divergences:

we will study the analyticity and Borel summability of super-renormalizable bosonic and fermionic models in d = 0, 1, and relativistic bosonic models in d = 2.

Note that this work is not the first one on the subject, the constructive aspects of several bosonic and fermionic models having been investigated in the past. Relativistic bosons, for instance, have been studied in [START_REF] Loeffel | Pade Approximants and the Anharmonic Oscillator[END_REF][START_REF] Simon | Coupling Constant Analyticity for the Anharmonic Oscillator[END_REF][START_REF] Graffi | Borel Summability: Application to the Anharmonic Oscillator[END_REF] for d = 0; in [START_REF] Billionnet | Analyticity and Borel summability of the phi 4 models. 1. The Dimensions of d = 1[END_REF] for d = 1, in [START_REF] Rivasseau | Corrected Loop Vertex Expansion for Φ 4 2 Theory[END_REF][START_REF] Simon | Borel Summability of the Ground-State Energy in Spatially Cutoff φ 4 2[END_REF][START_REF] Eckmann | Decay properties and Borel summability for the Schwinger functions in P (φ) 2 theories[END_REF][START_REF] Eckmann | Borel summability of the mass and the S-matrix in φ 4 models[END_REF] for d = 2, and in [START_REF] Magnen | Constructive φ 4 field theory without tears[END_REF][START_REF] Magnen | Phase space cell expansion and borel summability for the Euclidean φ 4 3 theory[END_REF][START_REF] De Calan | Local existence of the Borel transform in Euclidean φ 4 4[END_REF] for d = 3, 4. A general but less rigorous argument for all super-renormalizable relativistic scalar fields with generic polynomial interactions has been given in [START_REF] Serone | The Power of Perturbation Theory[END_REF][START_REF] Serone | λφ 4 Theory I: The Symmetric Phase Beyond NNNNNNNNLO[END_REF] . The Borel summability in the large N limit is discussed in [START_REF] Billionnet | Analytic interpolation and Borel summability of the λ/N |φ N | :4 ) 2 models[END_REF][START_REF] Fröhlich | Borel summability of the 1/N expansion for the N -vector [O(N ) non-linear σ] models[END_REF] , and the non-relativistic case in [START_REF] Ginibre | The Classical Field Limit of Nonrelativistic Bosons. I. Borel Summability for Bounded Potentials[END_REF] . Fermions at finite temperature were investigated for instance in [START_REF] Gawedzki | Massless Lattice Phi**4 in Four-Dimensions Theory: A Nonperturbative Control of a Renormalizable Model[END_REF][START_REF] Gawedzki | Gross-Neveu model through convergent perturbation expansions[END_REF][START_REF] Gawedzki | Exact renormalization for the Gross-Neveu model of quantum fields[END_REF][START_REF] Benfatto | Ward identities and chiral anomaly in the Luttinger liquid[END_REF][START_REF] Lesniewski | Effective Action for the Yukawa(2) Quantum Field Theory[END_REF] . Other standard references and textbooks on this subject are [START_REF] Glimm | A Tutorial Course in Constructive Field Theory[END_REF][START_REF] Rivasseau | From Perturbative to Constructive Renormalization[END_REF][START_REF] Mastropietro | Non-perturbative renormalization[END_REF][START_REF] Salmhofer | Renormalization: An introduction[END_REF] The main addition of our paper is to offer a completely unified perspective using the modern method of LVE.

In this paper, we prove the existence of a constructive expansion for general superrenormalizable fermionic and bosonic quartic models in low dimensions. Then, we extend the results to the quenched Sachdev-Ye-Kitaev (SYK) model [START_REF] Kitaev | A Simple Model of Quantum Holography[END_REF][START_REF] Kitaev | A Simple Model of Quantum Holography (Part II)[END_REF][START_REF] Maldacena | Comments on the Sachdev-Ye-Kitaev Model[END_REF][START_REF] Polchinski | The Spectrum in the Sachdev-Ye-Kitaev Model[END_REF] . This is the main new result of this paper since the existence of a constructive expansion for SYK models has not been addressed previously. While the cases of fermionic and bosonic fields in low dimensions have been largely investigated in the past, we stress again that our paper is the first to deal with them using the LVE and to provide a unified perspective with general quartic interactions.

As explained, this paper is the first of a series, and some questions remain open. Among them, the case of higher-order interactions, which we expect to be treated using improved versions of the LVR (which do not exist today); and the case of just renormalizable interactions. The last case is a difficult task and is also an open problem in tensorial field theories [START_REF] Delepouve | Constructive Tensor Field Theory: The T 4 3 Model[END_REF][START_REF] Lahoche | Constructive Tensorial Group Field Theory II: The U (1)-T 4 4 Model[END_REF][START_REF] Rivasseau | Constructive tensor field theory: The T 4 4 model[END_REF] where the LVE has been used successfully. [START_REF]Because the φ 4 model in d = 4 is not asymptotically free[END_REF] The details on the models that we consider and the main statements are given in Section II, the proofs being relegated in Section III. Appendix A describes the BKAR forest formula, and Appendix B Nevanlinna theorem.

II. THE MODELS, BASICS AND MAIN STATEMENTS

In this section, we define the class of models considered in this paper and fix the notations.

For convenience, the proofs are given in Section III. We recall some useful properties of ddimensional bosonic and fermionic vector models, including in particular a description of the propagators and the respective power-counting theorems. Finally, we outline the main statements of this paper.

A. The models and power counting

This paper is devoted to a family of statistical models for random N -component vector

fields {χ i } in d dimensions, χ i : R d → R, C, i = 1, • • • , N , collectively denoted as Ξ := {χ i }. ?
The value at the point x ∈ R d is denoted as χ i (x). The field components χ i can be either ordinary commutative numbers (χ i χ j = χ j χ i ) or non-commutative (Grassmann valued) numbers (χ i χ j = -χ j χ i ), depending if we describe bosonic or fermionic random degrees of freedom. For the rest of this paper and to clarify the notations, we will denote as φ i the bosonic variables and as ψ i the fermionic ones. Φ = {φ i } then denotes the collective bosonic vector, and Ψ := {ψ i } the analogous fermionic vector.

We focus on statistical models, that is to say, quantum field theory in imaginary (Euclidean) time. Without loss of generality, the probability law p(Ξ) for a configuration Ξ is fixed by the choice of a classical action S(Ξ):

p(Ξ) = 1 Z e -S(Ξ) , (II.1)
the partition function Z ensuring the normalization of the probability distribution,

Z := N i=1 dχ i e -S(Ξ) , (II.2)
where dχ i denotes the generalized Lebesgue measure.

In this paper, we limit our investigations to quartic models, the first non-Gaussian contributions compatible with Z 2 -symmetry. For real fields χ i : R d → R, we choose a classical action of the form:

S(Ξ) = dx   1 2 i,j χ i (x)C -1 ij χ j (x) + u 4! ijkl W ijkl χ i (x)χ j (x)χ k (x)χ l (x)   . (II.3)
In this expression, C ij denotes the covariance matrix, u the coupling constant and W ijkl the coupling tensor. Note that the interaction is local in the usual meaning in field theory.

Moreover, in these notations, C ij is understood as a differential operator with respect to the space variable x:

C -1 ij := M ij + a (1) ij,α ∂ ∂x α + a (2) ij,αβ ∂ 2 ∂x α ∂x β + • • • , (II.4)
where the Greek indices α, β run from 1 to d and the Einstein convention for repeated indices is assumed. Note that the mass matrix M ij as well as the weight matrices a

(n) ij do not depend on the positions x. For bosonic fields, χ i χ j and χ i χ j χ k χ l are completely symmetric tensors, whereas they are completely antisymmetric tensors for fermionic fields such that:

• bosonic fields: C ij = C ji and W ijkl is completely symmetric.

• fermionic fields: C ij = -C ji and W ijkl is completely anti-symmetric.

For complex fields χ i : R d → C, we choose the classical action to be of the form:

S(Ξ) = dx   i,j χi (x)C -1 ij χ j (x) + u 4! ijkl W ijkl χ i (x) χj (x)χ k (x) χl (x)   . (II.5)
Once again, the permutation symmetries of the covariance and coupling tensors depend on the nature of the fields. Requiring the classical action to be a real functional of the fields, we get the conditions:

• bosonic fields: C ij = C † ij and:

W ijkl = W klij = W ilkj = W kjil = W * jilk . (II.6)
• fermionic fields: C ij = -C † ij and:

W ijkl = W klij = -W ilkj = -W kjil = W * jilk . (II.7)
Without much loss of generalities, we assume that the propagator takes the form of a direct product between a purely spatial part and a purely internal flavor space:

C -1 ij := K -1 ij C-1 , (II.8)
choosing the spatial part of the propagator C among three kinds of models relevant for physics and compatible with standard axioms of field theory:

• Schrödinger propagator: non-relativistic bosons and fermions in (d -1)-dimensional Euclidean space:

C-1 = ∂ ∂t -κ * ∆ + m (II.9)
where ∆ is the usual Laplacian over R d-1 : ∆

:= d-1 i=1 ∂ 2 ∂x 2 i
, and κ * := 2 2m * is a positive constant involving the effective mass m * .

• Klein-Gordon propagator: relativistic bosons in d-dimensional Euclidean space-time:

C-1 = - ∂ 2 ∂t 2 -∆ + m 2 .
(II.10)

• Dirac propagator: relativistic fermions in d-dimensional Euclidean space-time:

C-1 = γ α ∂ α + m. (II.11)
The have to be understood as 2 n -anti-commutative vectors:

ψ T i = (ψ i,1 , • • • , ψ i,2 n ).
In all these cases, the mass squared is positive, m 2 > 0 (unbroken phase). 78 The matrix K ij , on the other hand, is a real invertible N × N matrix. Note that, in addition to the external indices i and j, the matrix K may depend on some internal indices like spin, as it will be the case for Majorana fermions or Dirac fermions. The choice of this matrix depends on the specificity of the model that we consider, however, we can expect it to be bounded:

|K ij | ≤ κ , ∀i, j , κ ∈ R . (II.12)
For large N , additional conditions can be necessary to make the sums over loops not too large. With this respect, it may be reasonable to impose a condition such that K ij goes to zero for i, j 1. The sums are expected to converge for exponential decay, but it has to be checked for slower decays, especially for power decays:

|K ij | i,j 1 min 1 i ε , 1 j ε , > 1.
(II. 13) In this paper, to simplify the proofs, we will assume that the field components are identically distributed K ij = δ ij , but the conclusions can be easily extended to the generic case K ij = δ ij .

We will be interested as well in the case where the coupling tensor has its own degrees of freedom. In particular, we will consider the large N limit of the SYK model, introduced initially in [START_REF] Kitaev | A Simple Model of Quantum Holography[END_REF][START_REF] Kitaev | A Simple Model of Quantum Holography (Part II)[END_REF][START_REF] Maldacena | Comments on the Sachdev-Ye-Kitaev Model[END_REF][START_REF] Polchinski | The Spectrum in the Sachdev-Ye-Kitaev Model[END_REF] (see [START_REF] Rosenhaus | An Introduction to the SYK Model[END_REF][START_REF] Sárosi | AdS 2 holography and the SYK model[END_REF]

S = dt 1 2 ψ i d dt ψ i + 1 4! J ijkl ψ i ψ j ψ k ψ l . (II.14)
where the index i runs from 1 to N , with random couplings (quenching), which is achieved by integrating over J ijkl with a Gaussian measure. One considers the quenched partition function:

Z(u) = dJ e -S(J ) dψ e -S(J ,ψ) =: dJ e -S(J ) Z(J ) , (II. 15) where we call Z(J ) the partial partition function, integrating over the fermionic degrees of freedom only. Once again, note that J is assumed to be real, with strictly positive eigenvalues ξ 2 I . The measure for J , dJ e -S(J ) is Gaussian:

S(J ) := α(N ) 2ḡ 2 ijkl J ijkl J ijkl , (II.16)
where α(N ) is a certain power of N . Note that the generalized SYK model involves an interaction of degree q, J i 1 •••iq ψ i 1 • • • ψ iq : our conclusion about the Borel summability holds only for the original q = 4 SYK model. The function α(N ) is fixed to N q-1 usually, and we will set α(N ) = N 3 for our considerations.

The classical action for the fermions is chosen to be:

S(J , ψ) := β/2 -β/2
L(Ψ, J )dt (II.17)

with the Lagrangian density:

L(Ψ, J ) := 1 2 i ψ i d dt ψ i - u 4! ijkl J ijkl ψ i ψ j ψ k ψ l . (II.18)
Note the minus sign in front of the interaction term, which arises from -1 = i q/2 for q = 4, will be crucial for the convergence of the constructive bound.

In this paper, we focus on UV finite and super-renormalizable field theories in low dimensions. Power-counting provides a criterion to decide if a theory is UV safe or not (we do not consider IR problems in this paper). The multi-scale analysis is a rigorous method to count the divergences of the theory [START_REF] Rivasseau | From Perturbative to Constructive Renormalization[END_REF] . It introduces a slicing in the interior of the Feynman amplitudes, decomposing each propagator edge into slices. More concretely, regularizing all the propagators using the Schwinger parametrization, we introduce a slice index i ≥ 1 and a thickness parameter M , in such a way that the full propagator C in Fourier space is given by the sum:

C = imax i=1 C(i) + C(0) , (II.19)
the integer i max being fixed. The sliced propagator C(i) for i ≥ 1 is given by:

C(i) (ω, p 2 ) := (-iω + E c ) M -2(i-1) M -2i dα e -α(ω 2 +E 2 c ) , (II.20)
for non-relativistic bosons and fermions, the kinetic energy being E c := κ p 2 + m;

C(i) (ω 2 , p 2 ) := M -2(i-1) M -2i
dα e -α(ω 2 + p 2 +m 2 ) , (II.21)

for relativistic bosons, and

C(i) (ω, p ) := (-γ µ ∂ µ + m) M -2(i-1) M -2i
dα e -α(ω 2 + p 2 +m 2 ) , (II. 22) for relativistic fermions.

In the amplitude decomposition, M imax has to be interpreted as the UV cut-off, ensuring the finiteness of the loop integrals for large momenta. The remainder C(0) ij includes integration over Schwinger parameter α from 1 to ∞, and is completely regular in the UV. From the slice decomposition of the amplitude, it is not hard to show the following statement: Theorem 1 Let A G be the amplitude of the connected Feynman graph G with V vertices and L internal propagator edges. Introducing a slicing for each propagator edge, we denote by µ = {i 1 , • • • , i L } a scale assignment for edges labelled from 1 to L. The amplitude can then be decomposed as a sum over scale assignments such that

A G = µ A G, µ .
(II.23)

For fixed µ, let G i = k(i)
k=1 G k i ⊂ G be the subgraph of G made of edges with scale assignments higher or equal to i, and G k i its kth connected component, and where k(i) is the total number of connected components. We have the following statement in the UV sector, that is, for scale assignments higher or equal to 1:

|A G, µ | ≤ K L(G) i,k M Ω(G k i ) , (II.24)
the bound being uniform for some constant K. For a quartic model we have:

• non-relativistic bosons and fermions:

Ω(G k i ) = 1 2 (d -3)V (G k i ) - d -1 4 N (G k i ) + 1 2 (d + 1) , (II.25)
• relativistic bosons:

Ω(G k i ) = (d -4)V (G k i ) -N (G k i ) d 2 -1 + d , (II.26)
• relativistic fermions:

Ω(G k i ) = (d -2)V (G k i ) - d -1 2 N (G k i ) + d , (II.27)
In these expressions, N (G) and V (G) denote respectively the numbers of external edges and vertices of the diagram G.

The proof of this theorem is classic, and the reader may consult [START_REF] Rivasseau | From Perturbative to Constructive Renormalization[END_REF] for details. We have to specify the ranges of the momenta ω and p, which may be different. For non-relativistic models at finite temperature β, the time coordinate is periodic t ∼ t + β and a fundamental

interval is t ∈ [-β/2, β/2].
Then, the conjugate momentum ω takes discrete values:

ω = 2πn β (bosons), ω = (2n + 1)π β (fermions) , (II.28)
where n ∈ N. These two constraints for bosons and fermions come respectively from the usual periodicity and anti-periodicity imposed on Fourier basis functions. The zerotemperature limit β → ∞ corresponds to non-compact Euclidean time t ∈ R (without periodic conditions). The spatial momentum p, however, is not restricted, except if we impose to the system to leave into a box of finite size. For our purpose, we consider p ∈ R d-1 .

We deduce the following renormalizability criteria for each case studied in this paper:

• Non-relativistic bosons and fermions are super-renormalizable for d < 3, just-renormalizable for d = 3 and non-renormalizable otherwise.

• Relativistic bosons are super-renormalizable for d < 4, just renormalizable for d = 4

and non-renormalizable otherwise.

• Relativistic fermions are super-renormalizable for d < 2, just renormalizable for d = 2, and non-renormalizable otherwise.

B. Main statements

We now outline our main results. We call tree expansion the expansion of the free energy obtained from BKAR forest formula (Appendix A). We have:

1. First, for fixed and positive (see definition 1) coupling tensor W,

|u| ≤ O (1) 1 N λ 2 0 cos 2 ϕ 2 , (II.29)
and corresponds to the Borel sum of the perturbative expansion of the free energy F (u).

2. Secondly, for disordered coupling with standard deviation J 2 = ḡ/N 3 and Majorana fermions (SYK model), we have:

Proposition 2 In the large N limit, the tree expansion of the free energy F (u) of the SYK model, for u = ρe iϕ is analytic in the interior of the cardioid domain:

ρ ≤ O (1) 1 ḡ cos 2 ϕ 2 , (II.30)
and corresponds to the Borel sum of the perturbative expansion in u.

The rest of this paper is devoted to the proof of these statements.

Note that we do not consider the non-relativistic bosonic models in dimension d = 2, even though it is power counting super-renormalizable. The reason is that we expect such divergences to be unphysical because the non-relativistic approximation is valid only with a cut-off. In dimensions d ≤ 1, however, no physical divergence is expected. Indeed, the power counting Ω = 1 -V provides a logarithmic divergence for V = 1, coming from formal sums as ω 1/(iω + m). These divergences however can be removed by remembering that the path integral construction requires time discretization with finite time step , providing additional oscillations e iω . These oscillations allows to use the standard trick to compute sums like ω g(iω) by replacing them by a contour integration into the complex plane. For β = 2π, the net result is:

n∈Z 1 iω n + y ≡ ∓ 2π 1 e -2πy ∓ 1 . (II.31)

III. PROOFS

In this section, we give the proofs of the two statements 1 and 2 using the LVE. Our A. Bosons in dimensions d ≤ 1

Existence of a finite analyticity domain

We start by proving that the first requirement of Nevanlinna's theorem holds, that is, the existence of a finite analyticity domain. The zero-dimensional model is simpler and captures the dependence on the internal indices. Since there are no UV divergences, the case d = 1 follows easily.

Let Z(u) be the partition function of the bosonic model in d = 0:

Z(u) = dµ C (φ) e -u 4! ijkl W ijkl φ i φ j φ k φ l , (III.1)
where dµ C (φ) is a shorthand notation for the normalized Gaussian measure with covariance

C -1 , dµ C (φ) := N i=1 dφ i e -1 2 φ i C -1 ij φ j . For d = 0, C ij ≡ δ ij for i.i.d variables.
The interaction tensor can be viewed as a matrix W IJ , where the big indices I, J ≡ (ij), (kl) run from 1 to N (N + 1)/2. As a symmetric matrix with real coefficients, W IJ can be diagonalized with eigenvalues {Λ I }:

W IJ = L Λ L O IL O T LJ , (III.2)
where O is an orthogonal matrix. Defining the new field Ψ L := I=(ij) O IL φ i φ j , the partition function becomes:

Z(u) = dµ C (φ) L e -u 4! Λ L [Ψ L ] 2 . (III.3)
If we suppose u to be positive, a simple way to ensure stability is to assume that the eigenvalues Λ L are all positive. We adopt the following definition:

Definition 1 A coupling tensor W is said to be positive definite if all the eigenvalues Λ L are positive.

In fact, we will remove this condition formally, venturing into the complex plane u = ρe iϕ , and show that summability makes sense when u becomes complex. Now, we introduce the matrix-like intermediate fields σ (L) , indexed with a pair of indices L, with the normalized Gaussian measure:

dν(σ) := L e -1 2 (σ (L) ) 2 dσ (L) L e -1 2 (σ (L) ) 2 dσ (L) , (III.4)
in order to break the quartic interaction [Ψ L ] 2 as a three-body interaction:

Z(u) = dµ C (φ)dν(σ) L e √ -u 12 λ L Ψ L σ (L) , (III.5)
where λ 2 L := Λ L . Note that λ L Ψ L := A L is a symmetric matrix in the original little indices L ≡ (kl). Therefore, only the symmetric part of the matrix σ (kl) is relevant, and the antisymmetric part of the intermediate field matrix can be integrated out. ? The integration over the original field variable is Gaussian with effective covariance C -1 I --u 3 O I (σ). Performing this integration leads to the factor:

det   1 - -u 3 CO(σ)   -1/2 = e -1 2 Tr ln 1- √ -u 3 CO(σ) , (III.6) with O(σ) := L O IL λ L σ (L)
. Then, we expand the partition function in powers of the ef-

fective vertex V(O(σ)) := 1 2 Tr ln 1 --u 3 CO(σ)
. Next, we introduce replicas, replacing the intermediate fields σ (L) by n vector fields σ (L) p , with p running from 1 to n and with inverse covariance matrix 1, with all entries equals to 1. Denoting as dν 1 (σ) the corresponding normalized Gaussian measure, and exchanging sums and integrations, the partition function can be rewritten as:

Z(u) = n dν 1 (σ) (-1) n n! n p=1 V(O(σ) p ) . (III.7)
From the Brydges-Kennedy-Abdesselam-Rivasseau (BKAR) forest formula (see Appendix A), and the elementary properties of Gaussian integration, one arrives at the following expression for the free energy: 82

F (u) = n 1 n! Tn ∈Tn 1 0 dx e 1 2 pq ∂ ∂χ (I) p W IJ Xpq(x ) ∂ ∂χ (J) q ∈Tn 1 2 ∂ ∂χ (I) i( ) W IJ ∂ ∂χ (J) j( ) n p=1 V(χ p ) χ (L) p =0 , (III.8)
where we defined χ (L) p := λ L σ (L) p , and performed the change of variable χ by χ → χ := Oχ. The term with n = 0 will be treated separately, and we define

F (u) = ∞ n=1 1 n! F n (u).
Computing the derivatives, we get:

F n (u) := (-1) n 2 n-1
Tn ∈Tn

1 0 dx dµ ⊗x (χ (I) p ) v∈Tn (c(v) -1)! - u 3 c(v) ∈Tn I t( ) I s( ) v∈Tn V I v,1 ,••• ,I v,c(v) (χ p ) χ (L) p =0
, (III.9)

where the s( ) (resp. t( )) are sources (resp. targets) of the edge , which can be written as

couples s( ) = (v( ), i) 1 ≤ i ≤ c(v) (and similarly t( ) = (v( ), i)). The couple (v, i) with 1 ≤ i ≤ c(v)
denotes the boundaries of each of the c(v) edges hooked to v. The vertices

V I v,1 ,••• ,I v,c(v)
are defined as:

V I v,1 ,••• ,I v,c(v) :=   |v∈∂ ∂ ∂χ (J) j( )   V(χ v ). (III.10)
Explicitly, if we denote as (i v,n , j v,n ) ≡ I v,n the pair of indices labelled with I v,n , we have:

V I v,1 ,••• ,I v,c(v) = 1 2   c(v) n=1 δ jv,ni v,n+1   [R c(v)-1 (χ v )] i v,1 j v,c(v) , (III.11)
where we introduced the resolvent

R(χ) := 1 1 --u 3 χ . (III.12)
The following classical Lemma allows to bound the resolvent, and to extend the domain of convergence in the complex plane:

Lemma 1 Let E ≡ C N ⊗ L 2 (S 1 , R),
Let u = |u|e iϕ , ϕ ∈ (-π, π]
and R be the resolvent, defined as

R -1 = id + i √ uH , (III.14)
where id denote the identity operator on E. We have the following uniform bound for the standard induced operator norm on E:

R ≤ cos -1 (ϕ/2) . (III.15)
Our aim is to use this Lemma to bound each term in the tree expansion of F . We define:

B Tn := ∈Tn I t( ) I s( ) v∈Tn V I v,1 ,••• ,I v,c(v) (χ p ) . (III.16)
Let us consider a tree T and an edge ∈ T with boundary ∂ = {v, w}. Let us denote as

V (v) I and V (w)
I the two connected components connected together with the bridge , so that the amplitude A T reads:

B T = I,J V (v) I W IJ V (w) J = L O(V ) (v) L (λ L ) 2 O( V ) (w)
L , (III.17)

where

O(V ) (v) L := I V (v)
I O IL . Denoting as λ 2 0 ≤ I W II the highest eigenvalue λ L , we get:

B T ≤ (λ 0 ) 2 L O(V ) (v) L O( V ) (w) L = (λ 0 ) 2 L V (v) L V (w) L . (III.18)
Finally, using the standard Cauchy-Schwarz inequality:

L V (v) L V (w) L ≤ I |V (v) I | J |V (w) J |
recursively for all edges of the tree, we get the first bound:

B Tn ≤ (λ 0 ) 2(n-1) v∈Tn 1 2 Tr |R c(v) (χ p )| . (III.19) From Lemma 1, we have Tr |R c(v) (χ p )| ≤ cos ϕ 2 -c(v)
Tr I, where I is the N × N identity matrix. We thus get the second bound:

B Tn ≤ (λ 0 ) 2(n-1) 2 n 1 cos 2 ϕ 2 n-1 (Tr I) n = (λ 0 ) 2(n-1) 2 n 1 cos 2 ϕ 2 n-1 N n (III.20)
where we used the topological relation v c(v) = 2(n -1). From Cayley formula, we easily deduce that the number of labelled trees with n vertices for a fixed configuration for

coordination numbers {c(v)} is bounded by Ω(n, {c(v)}) ≤ n! v (c(v)-1
)! , and the remaining sum over coordination numbers is trivially bounded with the area of the (n -1)-dimensional sphere of radius 2(n -1), leading to:

1 n! |F n | ≤ 2 π 24e 2 | cos(ϕ/2)| 2 u λ 2 0 N λ 2 0 cos 2 ϕ 2 u 6 π e n , (III.21)
where we used of the Stirling formula. Finally, for F 0 , a simple integration by parts leads to:

dµ ⊗x (χ (I) ) Tr ln   1 - -u 3 χ   ≤ (λ 0 ) 2 dµ ⊗x (χ (I) ) 1 0 dt ρ 3 (1 --u 3 χt) 2 , (III.22)
which is trivially bounded. Thus, we explicitly showed the existence of a finite analyticity domain, the interior of the cardioid ρ ≤ 6

N λ 2 0 e π cos 2 ϕ 2 .
Remark 1 For the derivation of the bound, we assumed the i.i.d Gaussian measure K ij = δ ij . This condition can be easily relaxed. First of all, Lemma 1 still holds and the resolvent R(χ) remains bounded as cos -1 (ϕ/2). Each effective vertex involves one K -1 insertion per derivative with respect to the intermediate field, replacing Tr I → Tr K -c(v) . Now, we move on to the relativistic model in d = 1. Without lost of generality, we consider a complex model φ i (x) ∈ C at finite temperature β, such that the field can be expanded in Fourier modes as

φ i = 1 √ β n∈Z ϕ i (n) exp i 2πn β t . The propagator of the theory is then C(ω, ω ) = C(ω)δ(ω -ω ), with: C(ω) = 1 ω 2 + m 2 .
(III.23)

As in the previous case, we introduce Ψ L := I=ij O IL φ i φj , and break the quartic interaction dt|Ψ L | 2 introducing a time-dependent intermediate field σ (L) with normalized Gaussian measure dν(σ):

Z(u) = dµ C (φ, φ)dν(σ) L e √ -u 2 λ L dtΨ L σ (L) . (III.24)
Integrating over the complex field Φ, we get an effective matrix-field theory:

Z(u) = dν(σ) e -Tr ln K[σ] , (III.25)
where the big trace Tr sums over internal indices and frequencies and K is the operator with entries (W

1 2 ijkl := O ij,mn λ mn O T mn,kl ): K ij (ω, ω ) := δ ij δ ωω - -u 2β 1 ω 2 + m 2 kl W 1 2 ijkl σ kl (ω -ω ) .
(III.26)

The only changes come from the dependence on ω of the intermediate fields, as well as the additional bound coming from ω-integrations. Introducing the fields χ(ω) like in (III.12), the new resolvent is:

R[χ](ω, ω ) := K -1 [χ](ω, ω ) . (III.27)
Note that the operator CΣ, with Σ having entries Σ ij,ωω :

= kl W 1 2
ijkl σ kl (ω -ω ), is not Hermitian. However, this operator appearing in the trace, we can replace it by C 1/2 ΣC 1/2 , which is obviously Hermitian. The resolvent is thus bounded by cos -1 (ϕ/2) thanks to the Lemma 1. Defining

F n = (-1) n 2 n-1
Tn A Tn , where A Tn denotes the amplitude of the tree T n with n vertices, we have the following Lemma:

Lemma 2

The amplitude A Tn for the amplitude associated to the tree T n with n vertices admits the following bound:

|A Tn | ≤ u λ 2 0 4β n-1 (Tr I) n cos 2n-2 (ϕ/2) β 2m 2n-2 coth 2n-2 βm 2 ω 1 . (III.28)
Proof. This statement can be easily proved recursively in the number of vertex. Let us consider the case n = 2. The amplitude reads explicitly:

A T 2 = 1 2 ω 1 0 dx dν X [χ] ∂ ∂χ (I) 1 (ω) W IJ ∂ ∂χ (J) 2 (-ω) 2 p=1 V(χ p ) χ (L) p =0 , (III.29)
where the Gaussian measure is for the propagator

X IJ ij (ω, ω ) = W IJ x ij δ ω,-ω .
Computing the derivatives with respect to χ

V(χ 1 ) = - ω ω -u 2β 1 ω 2 + m 2 δ(ω -ω + ω )V I [χ 1 ](ω , ω ) .
Formally, the last term has the structure of a trace over ω, --u 2β Tr BV I [χ 1 ], involving the strictly positive matrix

B ω ω (ω) := 1 ω 2 + m 2 δ(ω -ω + ω ).
(III.30)

For A bounded and B strictly positive | Tr(AB)| ≤ A Tr B. However, B is not diagonalizable. In fact, only the kth diagonal of the matrix is non-vanishing (counting from the principal diagonal): B ij = b i δ i,j-k with b i > 0 and k = 0. For our purpose, A ≡ R, and the resolvent is unitary diagonalizable as a function of the Hermitian operator H. Denoting as U the unitary operator diagonalizing A and as a i its eigenvalues, we get:

| Tr(AB)| ≤ | l,i a l U i-k,l U † l,i b i | ≤ A l,i |U i-k,l U † l,i b i | . (III.31)
Using Cauchy-Schwarz inequality, the last product can be bounded as:

i,l |U i-k,l U † l,i b i | = i b i l |U i-k,l ||U i,l | ≤ i b i l |U i-k,l | 2 l |U i,l | 2 , (III.32)
that is to say, because U U † = id: | Tr(AB)| ≤ A ( i b i ), and we get:

83 I ∂ ∂χ (I) 1 (ω) V(χ 1 ) ≤ 1 cos(ϕ/2) u 2β 1/2 (Tr I) β 2m coth βm 2 , (III.33)
where we used Lemma 1. Similarly to the d = 0 case, from the Cauchy-Schwarz inequality and taking into account the normalization of the integrals, we get:

|A T 2 | ≤ u λ 2 0 16 (Tr I) 2 cos 2 (ϕ/2) β m 2 coth 2 βm 2 ω 1 . (III.34)
Now, we assume that the statement holds for trees having n -1 vertices. Any tree with n vertices can be obtained from a tree with n -1 vertices by adding a leaf. Up to Gaussian and x integrations, we get the explicit expression for the amplitude A Tn :

A Tn ∝ 1 2 n u 2β n ω R K,ω W KL U L; ω , (III.35)
R K,ω denoting the rest of the diagram, up to the added leaf, the Cauchy-Schwartz inequality leads to the bound:

| ω R K,ω W KL U L; ω | ≤ λ 2 0 R K,ω U L; ω (Tr I) ω 1 . (III.36)
We introduced the notation:

U I 1 ,•,In; ω 1 ,••• , ωn :=   3β -u   n n i=1 ∂ ∂χ I i ( i ω i ) V(χ p ) , (III. 37
)
where i = ±. Explicitly, the structure of this vertex is the following:

U I 1 ,•,In; ω 1 ,••• , ωn = R j 1 i 2 ; ω1 ω 1 (Cχ) ω 1 ω2 ;I 1 ω 1 R j 2 i 3 ; ω2 ω 2 • • • R jni 1 ; ωn ω n (Cχ) ω n ω1 ;Inωn (III.38)
where we integrated over repeated ω-indices, and explicitly:

(Cχ) ω i ωi+1 ;I i ω i = 1 ω 2 i + m 2 δ( i ω i -ω i + ωi+1 ) .
(III.39)

The effective vertex can be easily bounded using successive Cauchy-Schwartz inequalities as:

I 1 ,••• ,I k |U I 1 ,••• ,I k ; ω 1 ,••• ,ω k | ≤ Tr I cos k (ϕ/2) δ k i=1 i ω i β 2m k coth k βm 2 . (III.40)
Therefore, from the recursion hypothesis, the positivity of the propagator and the momentum conservation along the loops ensured by the Dirac delta, we must have for the one-point amplitude R K,ω :

|R K,ω | ≤ u λ 2 0 4β n-2 (Tr I) n-1 cos 2n-4 (ϕ/2) β 2m 2n-4 coth 2n-4 βm 2 .
(III.41)

The statement following from inequality (III.36).

Thus, the series is absolutely convergent in the interior of the cardioid domain:

ρ ≤ 16 e/π m 2 N λ 2 0 β coth 2 βm 2 cos 2 ϕ 2 .
(III.42)

Note that the massless limit of (III.42) is ill-defined because the RHS goes to zero. Moreover, in the continuum limit β → ∞ (zero-temperature limit), the bound still exists. Indeed, in the low temperature limit, t ∈ R, and the sums over ω become integrations:

ω f (ω) → β 2π dω f (ω) .
(III. 43) In that limit, β completely disappears, and the Fourier series become continuous Fourier transformations:

φ i = 1 √ 2π +∞ -∞ dω ϕ i (ω)e iωt , (III.44)
and from the integral +∞ -∞ dω ω 2 +m 2 = π m , the rule to pass from finite to infinite β is to replace β by 2π everywhere, and set coth(βm/2) = 1. Thus the cardioid domain becomes:

ρ ≤ 8 e/π 3 m 2 N λ 2 0 cos 2 ϕ 2 .
(III.45)

Borel summability

The last piece to prove Borel summability is to check the second requirement of Nevanlinna's theorem. The proof is classic and follows immediately from the previous results, and we provide only the main steps in d = 0, the extension to d = 1 being straightforward.

The remainder of the expansion of the free energy reads (g = u/4!):

R r F (g) := g r+1 1 0

(1 -t) r r! (F ) (r+1) (tg)dt , (III.46) with F := ∞ n=1 F n = F -F 0 .
The remainder can be expanded as a sum over labelled trees thanks to the BKAR formula, and we get:

R r F (g) := - ∞ n=1 1 n! Tn ∈Tn 1 0 dx dµ ⊗x (χ (I) p )R r [Y Tn ] , (III.47)
where

Y Tn := g n-1 v∈Tn (c(v) -1)! ∈Tn I t( ) I s( ) v∈Tn V I v,1 ,••• ,I v,c(v) (χ p ) χ (L) p =0,xpq=1 . (III.48) First, note that when n -2 ≥ r, R r [Y Tn ] = Y Tn .
Each of these terms admits a bound of the form K n |λ| n , and the sum converges at least in the interior of the cardioid. When n -2 < r, however, the remainder has to be computed following a Taylor expansion of the resolvents in Y T . This expansion dresses the trees with loops. Concretely, we extract the global factor g n-1 in front of Y T , and we write:

Y Tn = g n-1
ȲTn , (III.49)

such that: R r [Y Tn ] = g n-1 R r-n+1 [ ȲTn ] . (III.50)
Defining z := -u/3, we have:

d dz Tr R(χ) = n! Tr[χ n R n+1 (χ)] . (III.51)
As announced, the tree T n becomes a spanning tree for a graph with loops, generated from the new Gaussian integrations of the fields χ v coming from the Taylor expansion of the resolvent. Each power of the resolvent can be bounded using Lemma 1, and we deduce:

R r-n+1 [ ȲTn ] ≤ z 2(r-n) 2 n 1 0 dt (1 -t) 2r-2n+2 (2r -2n + 2)! {k l }| k l =2r-2n+3 (2r -2n + 3)! n l=1 k l ! × n v=1 (c(v) + k v -1)! × Tr   I v∈Tn χ kv v   χ (L) p =0
.

(III.52)

We can now plug this expression in equation (III.47). First, note that because

C p n ≤ 2 n , the factor c(m)(c(m) + k m )!/[k m !c(m)!]
is bounded by e ln(c(m)) 2 c(m)+km ≤ (2e) c(m) 2 km , and the product over m gives a factor (2e) 2n-2 2 2r-2n+3 . Secondly, we can perform the Gaussian integration. Because there are 2r-2n fields, the number of Wick contractions is (2(r-n)!! = 2 r-n (r -n)! ≤ 2 r r!. Finally, the remaining integration over t gives (2r -2n + 3) -1 , which together with the denominator (2r -2n + 2)! exactly compensates the combinatorial factor (2r-2n+3)!. As in the previous section, using Cayley's theorem for the number of trees with n vertices and Stirling's formula, we find a bound of the form: AB n 1 B r 2 r! for some constants A, B 1 and B 2 . Because n -2 < r, summing over n, we find the final bound: A |g| r B r r! for the contributions in (III.47) for which n -2 < r. As explained before, the contributions for n -2 ≥ r are all bounded by bounds of the form: |g| n K n , and the sum behaves as:

A |g| r K r . Ultimately, because, for positive constants k 1 and k 2 , k 1 r! + k 2 ≤ (k 1 + k 2 )r!, we find that |R r F (g)| ≤ A (B ) r |g| r r!.

B. Bosons in dimension d = 2

In this section, we prove the existence of a finite analytic domain with a cardioid shape.

It requires improving the standard LVE to deal with UV divergences, which is achieved by the so-called MLVE [START_REF] Gurau | Renormalization: An Advanced Overview[END_REF][START_REF] Gurau | The Multiscale Loop Vertex Expansion[END_REF] . It allows us to subtract scale by scale the divergent subgraphs, from an expansion requiring "forests into forests", and known as jungle expansion.

Renormalized model and jungle expansion

From the power-counting theorem (II.26), the degree of divergence for two-dimensional relativistic bosons becomes independent of the number of external edges N , and show that divergences arises only from diagrams with a single vertex, that is for (V, N ) = (1, 2) and (V, N ) = (1, 0), and are both logarithmic. We only consider the case of complex bosons.

From an elementary perturbative calculation using Feynman rules for the interaction:

S int := u 4! ijkl +β/2 -β/2 dt +∞ -∞ dx W ijkl φ i (x, t) φj (x, t)φ k (x, t) φl (x, t) , (III.53)
we get the explicit expressions of the two divergent diagrams:

1. The vacuum diagram (V, N ) = (1, 0): = - u 2πβ ω +∞ -∞ dp ω 2 + p 2 + m 2 2 ij W iijj , (III.54) 2. The 2-point diagram (V, N ) = (1, 2): i j = - u πβ ω +∞ -∞ dp ω 2 + p 2 + m 2 k W ijkk .
(III.55)

Both of these diagrams scale as ln(Λ), with Λ designating some arbitrary UV cut-off in the sums and integrals. A renormalized (i.e. UV finite) action can be obtained including two counter-terms in order to systematically subtract the divergent diagrams. The two counterterms corresponding respectively to 2-point and vacuum divergences are denoted by δV 1 and δV 2 and are respectively given by:

δV 1 ( Φ, Φ) := ω +∞ -∞ dp φi (ω, p)δm 2 ij φ j (ω, p) , (III.56)
with δm 2 ij =: ∆ 2 k W jikk , and:

∆ 2 := - u πβ ω +∞ -∞ dp ω 2 + p 2 + m 2 , (III.57)
and

δV 2 := - u 2πβ ω +∞ -∞ dp ω 2 + p 2 + m 2 2 ij W iijj .
(III.58)

The mass divergence arising from the tadpole can be subtracted by adding δV 1 to the original action. However, we have to be careful, because this counter-term generates a new vacuum divergence as well. Nonetheless, it is easy to check that this new divergence can be subtracted by a counter-term equals to -2δV 2 . Therefore, the cancellation of all vacuum divergences requires the counter-term δV ac = δV 2 -2δV 2 = -δV 2 . Thus, the completely renormalized model reads:

Z R (J, J) = e δV 2 dΦd Φe -S PR ( Φ,Φ)+ i ji φ i + i φi j i , (III.59) with:

S int,PR ( Φ, Φ) = 1 2 u 2πβ W ijkl φ i φj φ k φl + ∆ 2 k W ijkk φ i φj .
Completing the square, we get:

S int ( Φ, Φ) = 1 2 u 2πβ I M I + 2πβ∆ 2 u X I 2 - πβ 2u (∆ 2 ) 2 I X 2 I , (III.60)
where we defined M kl := i,j W 1 2

ijkl φ i φj and

X ij := k W 1 2
ijkk , and the last term is explicitly:

πβ 2u (∆ 2 ) 2 I X 2 I = u 2πβ ω +∞ -∞ dp ω 2 + p 2 + m 2 2 ij W iijj = -δV 2 .
(III.61)

Next, introducing intermediate fields to break the square, integrating over the fields Φ and Φ, and rewriting the determinant as a Tr-log effective interaction, we finally get the effective model for the intermediate fields σ = {σ I }:

Z(u) = dν(σ)e -Tr ln K[σ]-2πβ -u ∆ 2 I X I σ I , (III.62)
where K is the matrix with elements: We refer to "momentum space" as the set of p into the disk of radius Λ. In order to use the multi-scale expansion, we need to introduce a slicing into the momentum space. To this end, we consider a pair of real and integer numbers (M, j max ) and assume that Λ = M jmax .

K ij (p, p ) := δ ij δ(p, p ) - -u 2πβ 1 p 2 + m 2 W
Intermediate scales are then defined as M j for j ≤ j max , and we introduce the slice function χ j (p) such that:

χ j (p) := χ ≤j (p) -χ ≤j-1 (p) j ≥ 2 , (III.67)
where χ ≤j (p) denotes the step function χ ≤j (p) := ϕ(M 2jp 2 ). We define

V ≤j := Tr ln 2 (id + χ ≤j Hχ ≤j ) , (III.68)
and

V j := V ≤j -V ≤j-1 , (III.69)
such that j V j = V ≡ Tr ln 2 (id + H). In terms of this slice decomposition, the partition function (III.65) reads:

Z(u) = dν(σ) jmax j=0 e -V j [σ] . (III.70)
The first MLVE trick is to view the product as a determinant, and to rewrite it as a Grassmann integral. More precisely, defining W j := e V j -1, we may rewrite (III.70) as:

Z(u) = dν(σ) jmax j=0
dµ 1 ( χj , χ j )e -jmax j=0 χj W j χ j , (III.71)

with the Gaussian Grassmann measure with identity kernel: dµ 1 ( χj , χ j ) := d χj dχ j exp (-χj χ j ) . (III.72)

We are now in the position to perform the MLVE of the model. We only reproduce the main step of the general method, referring to the standard reference [START_REF] Gurau | The Multiscale Loop Vertex Expansion[END_REF] for technical details. Let us define S = [0, j max ] to be the set of scales and I S the S × S identity matrix, allowing to rewrite the partition function Z(u) in the compact form

Z(u) = dν S e -W " = " ∞ n=0 1 n! dν s (-W ) n , (III.73)
where dν s := dν(σ)dµ Is ({ χj , χ j }). The quotes around the equality symbol refer to the illegal permutation of sum and Gaussian integrations. In fact, the aim of our procedure is to give a sense for this equality. As for the d = 0 and d = 1 cases, we introduce replicas for the vertices in the set

V := {1, • • • , n}: Z(u) = ∞ n=0 1 n! dν S,V n a=1 (-W a ) , (III.74)
in such a way that the vertex W a express in terms of bosonic intermediate fields σ a , and the Gaussian measure dν S,V becomes dν S,V = dν 1 V (σ)dµ Is ({ χj , χ j }), the covariance becoming the V × V matrix 1 V with all entries equals to 1. The MLVE trick for factorizing the integral is to use forest formula two times. The first one, over the bosonic intermediate field,

follows in the same way as for the previous cases. We introduce the covariance x ab = x ba , with x aa = 1 between the bosonic replicas, and use the forest formula (2) to write, in the derivative representation (holding as well for Grassmann fields):

Z(u) = ∞ n=0 1 n! F 1 0   ∈F dx   e 1 2 n a,b X ab (x ) ∂ ∂σ a ∂ ∂σ b + jmax j=0 ∂ ∂ χj ∂ ∂χ j × ∈F ∂ ∂σ a( ) ∂ ∂σ b( ) n a=1   - jmax j=0 χj W j (σ a )χ j   σ=χ= χ = 0 . (III.75)
where a( ) and b( ) denote the end points of the edge . The forest F defines a natural partition of the set V into blocks building of its connected spanning trees. We call such a block as B, and denote by V /F the reduced set building of bosonic spanning trees B. Next, we use the same replica trick by assigning the fermionic fields to these blocks, such that χ j → χ B j with covariance y BB = y B B and y BB = 1. Following standard notations, we call L F a generic fermionic edge between two blocks B and B . Using the forest formula a second time, we obtain an expansion indexed by two-level jungles J rather than forests [START_REF] Abdesselam | Trees, Forests and Jungles: A Botanical Garden for Cluster Expansions[END_REF] :

Z(u) = ∞ n=0 1 n! J j 1 • • • jn 1 0 dz J dν J ∂ J B a∈B χ B ja W ja (σ a )χ B ja , (III.76)
where:

• The two-level jungles J = (F B , F F ) are ordered pairs of bosonic and fermionic disjoint forests of the set V , denoted as F B and F F . We denote as B and L F the bosonic and fermionic edges of the two components of the jungle J . The notation being kept to denote generic edges of J .

• dz J means integration from 0 to 1 over parameters z = x for ∈ F B and z = y for ∈ F F , for each ∈ J .

• ∂ J is a compact notation for:

∂ J := B =(a,b)∈F B ∂ ∂σ a ∂ ∂σ b × L F =(p,q)∈F F δ jpjq   ∂ ∂ χB(a) ja ∂ ∂χ B(b) j b + ∂ ∂ χB(b) j b ∂ ∂χ B(a) ja   .
(III.77)

• The Gaussian measure dν J has covariance X(x ) ⊗ 1 S for bosons and Y (x ) ⊗ I S for fermions, that is to say, in derivative representation:

dν J ≡   e 1 2 n a,b X ab (x ) ∂ ∂σ a ∂ ∂σ b + BB Y BB (x ) a∈B,b∈B δ jaj b ∂ ∂ χB ja ∂ ∂χ B j b   σ=χ= χ=0
.

(III.78)

• X ab (x ) is the infimum of the parameters x for the bosonic component of the jungle, in the unique path from a to b. It is set to be equal to zero if the path does not exist, and to 1 if a = b.

• Y BB (y ) is the infimum of the fermionic parameters along the fermionic edges L F , between the blocks B and B .

As for the cases discussed in the previous section, our aim is now to bound the amplitudes indexed with jungles. More precisely, we will prove the existence of a finite analyticity domain for the free energy F (u) = ln Z(u), expanding in terms of connected jungles:

F (u) = ∞ n=0 1 n! J trees j 1 • • • jn 1 0 dz J dν J ∂ J B a∈B χ B ja W ja (σ a ) χB ja .
(III.79)

Bounds and convergence

In this section, we will bound the bosonic and fermionic contributions of the decomposition (III.79). The technical parts of the bounds are essentially the same as for the general treatment given in [START_REF] Gurau | The Multiscale Loop Vertex Expansion[END_REF] , and we refer to this paper for the technical subtleties, only indicating the main steps of the proof and focusing on the specificity of the model.

Fermionic integrals

The fermionic part of the expansion (III.79) is exactly the same as in [START_REF] Gurau | The Multiscale Loop Vertex Expansion[END_REF] . Due to the standard properties of Grassmann integration, the Gaussian integration over these variables can be written as: Taking into account the hard core constraint inside each block, meaning that the integral (III.80) vanishes as soon as two vertices belong to the same bosonic block B with the same scale attribution, one can rewrite the equation (III.80) as: This lemma can be easily checked, and the proof may be found in [START_REF] Gurau | The Multiscale Loop Vertex Expansion[END_REF] .

B a∈B d χB ja dχ B ja e -n a,b χB(a) ja Y ab χ B(b) ja L F ∈F F δ j a(L F ) j b(L F ) (III.80) × χB(a(L F )) j a(L F ) χ B(b(L F )) j b(L F ) + χB(b(L F )) j b(L F ) χ B(a(L F )) j a(L F )
B a,b∈B a =b (1 -δ jaj b ) L F ∈F F δ j a(L F ) j b(L F ) Y p 1 ,...,p k a 1 ,...,a k + Y m 1 ,...,p k p 1 ,...,m k + • • • + Y m 1 ,...,m k p 1 ,...,

Bosonic integrals

We now move on to the bosonic integrals. From the two-level trees decomposition of the free energy, it is clear that bosonic integrations factorizes over each block B. As a result, we can only consider a single block B. It involves the Gaussian integration:

dν B F B (σ) = e 1 2 a,b∈B X ab (w l ) ∂ ∂σa ∂ ∂σ b F B (σ) σ=0 (III.84)
with F B (σ) defined as:

F B (σ) = B ∈B dp B ∂ 2 ∂σ a( B ) (-p B )∂σ b( B ) (p B ) a∈B W ja (σ a ) , (III.85)
where we introduced the explicit momentum dependence and the integrations over momenta along each bosonic edge. The derivatives ∂/∂σ can be evaluated using the famous Faà di Bruno formula, which extends the standard derivation rule for composed functions. Computing the first derivative of the potential V j , we get, for k > 0:

k l=1 ∂ σ a(l) (-V j ) = π Tr (∂ σ a(π(1)) H)χ ≤j R ≤j • • • (∂ σ a(π(k)) H)χ ≤j R ≤j k times -(∂ σ a(π(1)) H)χ ≤j-1 R ≤j-1 • • • (∂ σ a(π(k)) H)χ ≤j-1 R ≤j-1 k times , (III.86)
where:

R ≤j := (1 + χ ≤j Hχ ≤j ) -1 , (III.87)
and the sum over π runs over permutation of k elements, up to cyclic permutations. Then, the k-th derivative of W j can be deduced from the Faà di Bruno formula (III.86). For k > 0 we get:

k l=1 ∂ σ a(l) (-W j ) = e -V j {m l } l≥1 lm l =k k! l≥1 m l !(l!) m l l≥1 [∂ l σ (-V j )] m l , (III.88)
where we used the compact notation l≥1 [∂ l σ (-V j )] m l to represent a block of derivative of size l. In (III.85), we can rewrite the product as a product over the arcs of the vertices:

F B ( σ) = v∈B   c(v) k=1 ∂ ∂σ a(v) ( k p k )   W jv (σ a(v) ), (III.89)
where once again c(v) denotes the coordination number of the vertex v, equal to the number of half lines of the intermediate-fields hooked to this vertex, and k = ±1 is a sign depending of the orientation of the momentum. As a result, the bosonic integral (III.84) becomes:

dν B m∈B e -V j (λ) c(m) {x (m) l } l≥1 lx (m) l =c(m) c(m)! l≥1 x (m) l !l x (m) l l≥1 [∂ l σ (-V j )] m l , (III.90)
where we used again the compact notations for the derivatives, and the integrations over momenta are implicit. Note that the bar over the product means that we extracted the factor (λ) l×x

dν B F B ( τ ) ≤ λ 2 cos 2 (φ/2) |B|-1 dν B m∈B e -V jm ( τm) × {x (m) l } l≥1 lx (m) l =c(m) c(m)! l≥1 x (m) l !l x (m) l l≥1 |∂ l σ (-V j )| x l (m) R=id , (III.91)
where in the right hand side we set R = id. Defining:

G B := m∈B {x (m) l } l≥1 lx (m) l =c(m) c(m)! l≥1 x (m) l !l x (m) l l≥1 |∂ l σ (-V j )| x l (m) R=id , (III.92)
and since the Gaussian measure dν B is positive, we can use the Cauchy-Schwarz inequality to get:

dν B m∈B e -V jm ( τm) G B ≤ dν B m∈B e -2V jm ( τm) 1/2 dν B |G B | 2 1/2 . (III.93)
We called the first term non-perturbative factor, and the second term perturbative factor, and shall treat separately each of them.

Final bounds

a. Non-perturbative bound. We begin with the first term of the bosonic integral, the non-perturbative contribution (following the notations of 32 ): ( τa) , (III.94) dν B being the restriction of the measure dν J to the block B. We aim at proving the following lemma:

B 1 := dν B a∈B e -2V ja
Lemma 4 For |u| small enough, B 1 defined in (III.94) satisfies the following uniform bound:

B 1 ≤ e O(1) |u|N 2πβ cos -1 (ϕ/2) . (III.95)
Proof. First of all, note that: e -2V ja ( τa) is uniformly bounded by ≤ e 2|V ja ( τa)| . Secondly, because of the identity:

ln 2 (1 -x) = 1 0 dt tx 2 1 -tx , (III.96)
and using the bound (III.66), we get the inequality:

|V j | ≤ 1 cos(ϕ/2) Tr(H ≤j ) 2 -Tr(H ≤j-1 ) 2 , (III.97)
with the definition H ≤j := χ ≤j Hχ ≤j . The cyclicity of the trace gives the relation:

Tr(H ≤j ) 2 -Tr(H ≤j-1 ) 2 = Tr H(χ ≤j Hχ ≤j -χ ≤j-1 Hχ ≤j-1 ) . (III.98)
Moreover, because of the definition χ j := χ ≤j -χ ≤j-1 , the difference may be rewritten as:

χ ≤j Hχ ≤j -χ ≤j-1 Hχ ≤j-1 = χ j Hχ ≤j-1 + χ ≤j-1 Hχ j (III.99)
and we get, using again the fact that the trace is cyclic:

Tr(H ≤j ) 2 -Tr(H ≤j-1 ) 2 = 2 Tr (Hχ j Hχ ≤j-1 ) .

(III.100)

Because χ 2 ≤j = χ ≤j and χ 2 j = χ j , and defining C j := χ j Cχ j and C ≤j := χ ≤j Cχ ≤j , the trace may be explicitly computed as:

Tr (Hχ j Hχ ≤j-1 ) = -uN 2πβ dp σ(-p ) dp C ≤j (p)C j (p -p ) σ(p ) , (III.101)
the factor N arising from the trace trace over internal indices Tr I = N . The covariance for bosonic variables is then modified as

X -1 B → X -1 B -8 cos -1 (ϕ/2)I B ⊗ L, L being the bounded operator of L 2 (S 1 , R): L(p ) := O(1) |u|N 2πβ dp C ≤j (p)C j (p -p ) , L ≤ O(1) |u|N 2πβ M -2j , (III .102)
O(1) denoting some numerical constant of order 1. Computing the Gaussian integration gives the determinant:

det id - 8 cos(ϕ/2) |u|N 2πβ X B ⊗ L -1/2 = e -1 2 Tr ln(id- 8 cos(ϕ/2) |u|N 2πβ X B ⊗ L) .
(III.103)

Given a positive bounded operator Â, we have the following bound:

-Tr ln(1 -Â) = ∞ n=1 Tr( Â) n n ≤ Tr( Â) × ∞ n=2 Â n n ≤ Tr( Â) × ∞ n=2 Â n . (III.104)
Setting  := 8|u|N 2πβ cos -1 (ϕ/2)I B ⊗ L, it is easy to check that:

Tr( Â) ≤ O(1) 8|u|N 2πβ cos -1 (ϕ/2)|B| , Â ≤ O(1) 8|u|N 2πβ cos -1 (ϕ/2) , (III.105)
leading to:

| Tr ln(1 -Â)| ≤ |B| ∞ n=1 O(1) 8|u|N 2πβ cos -1 (ϕ/2) n .
(III.106) b. Perturbative bound. We now move on to the perturbative bound:

B 2 := dν B |G B | 2 1/2 . (III.107)
and we have to prove the lemma:

Lemma 5 For |λ| small enough, the perturbative contribution B 2 satisfies the following bound, up to a global conservation factor:

|B 2 | ≤ (4|B| -4)!! × m∈B (O(1)π) c(m) c(m)!M -2jm . (III.108)
Proof. For l > 1 we have:

dp C l j (p) ≤ dp C 2 j (p) ≤ 1 M 4(j-1) dp χ j (p) .
(III.109)

The last integral is nothing but the volume between the disks of radius p 2 = M 2(j-1) and

p 2 = M 2j , as a result: dp C l j (p) ≤ O(1) π M 2j-4 .
(III.110)

The result may be generalized for integrals of the type:

I p 1 ,••• ,p l := dp (C ≤j (p)C ≤j (p -p 1 ) • • • C ≤j (p -p l ) l times -(j → j -1) . (III.111)
The domain of the integral is the intersection of disks centered in 0, p 1 , • • • , p l . In the worst case, p 1 = p 2 = • • • = 0, we get:

I p 1 ,••• ,p l ≤ M -2(l+1)(j-1) dp χ j (p) ≤ O(1) π M -2lj+2 ≤ O(1) π M 2j-4 .
(III.112)

To compute the Gaussian integral, we have to retain only the leaves of the trees, each of them involving a single intermediate field. In more details, because of the inequality (III.112), the Gaussian block has to be uniformly bounded by a Gaussian integral of the form:

H B := dν B a∈B O(1) π M 2ja-4 ka × (σ a ) ka .
(III.113)

Because we bounded the integrals over each loop with the worst bound, the remaining Gaussian integrations are momentum-independent. Then, the variables σ a do not depend on the momentum, and the Gaussian measure dν B has covariance

X B . Because ||X B || ≤ 1,
it follows that the Gaussian integration (III.113) is bounded by:

|H B | ≤ a∈B O(1) π M 2ja-4 ka × a∈B k a !! .
(III.114)

From equation (III.107), the double factorial is trivially bounded by (4|B| -4)!!. As a result, assuming j m > 2, we get the following bound 84 (up to a global conservation factor): [START_REF] Flato | Wightman formulation for the quantization of the gravitational field[END_REF] .

|B 2 | ≤ (4|B| -4)!! × m∈B {x (m) l } l≥1 lx (m) l =c(m) c(m)! l≥1 x (m) l !l x (m) l |λ| x (m) 1 2 (O(1)π) c(m) M -(2jm-
Assuming |λ| ≤ 1, {x (m) l } l≥1 lx (m) l =c(m) 1 l≥1 x (m) l !l x (m) l (III.115) is nothing but the coefficient of x c(m) in the Taylor expansion of k e x k /k = 1/(1 -x).
As a result, taking into account Lemmas 4 and 5, we deduce the final bound:

dν B m∈B e -V jm ( τm) G B ≤ e O(1) λ cos(φ/2) |B| (4|B| -4)!! × m∈B c(m)! (O(1)π) c(m) M -2jm , (III.116)
and the following proposition for the bosonic integration: Proposition 3 For |λ| small enough, the bosonic integration admits the bound:

dν B F B ( τ ) ≤ e O(1) λ cos(ϕ/2) |B| λ cos 2 (ϕ/2) |B|-1 N |B| × (4|B| -4)!! × m∈B c(m)! (O(1)π) c(m) M -2jm .
(III.117)

Finally, collecting the results for bosonic and Grassmann bounds, we get the uniform bound for the free energy ln(Z) (to simplify the expression we forget the exponential factor

| ln Z[J, J, λ]| ≤ ∞ n=1 1 n! J tree n k=1 jmax j k =0 2 L(F F ) f ∈F F δ j s( f ) j t( f ) × B m,m ∈B m =m (1 -δ jmj m )N |B| |λ| cos 2 (ϕ/2) |B|-1 × (4|B| -4)!! m∈B c(m)! (O(1)π) c(m) M -2jm , (III.118)
the factor 2 L(F F ) involving the number of fermionic edges L(F F ). Because of Cayley's Theorem, the number of trees with n labelled vertices and coordination numbers c i for each vertex i = 1, . . . , n is (n -2)!/ i (c i -1)!, and the sum involved in (III.118) obeys to

c(m) m c(m)=2|B|-2 m∈B c(m) = (3|B| -3)! (|B| -2)!(2|B| -1)! . (III.119)
Collecting all the factorials and using Stirling's formula, we get:

2 (4|B| -4)!! (3|B| -3)! (2|B| -1)! ≤ (|B| -1)! 3 |3B| e -|B| |B| |B| . (III.120)
We now move on to sum over scale attributions, taking into account the hard core constraint.

As explained in full details in [START_REF] Gurau | The Multiscale Loop Vertex Expansion[END_REF] , the hard core constraint imposes that the scale assignments of vertices in a same block are all different, implying that:

m∈B (j m -2) ≥ 1 2 m∈B j m + j min -2 2 |B| + |B|(|B| -5) 4 , (III.121)
where we have introduced explicitly the minimal scale j min > 2. Then:

{jm} m,m ∈B m =m (1 -δ jmj m ) m∈B M -jm ≤ jmax j=j min M -j/2 |B| 1 M j min -2 2 |B|+ |B|(|B|-5) 4
(III.122) which, for j min > 2 and M > 4, is uniformly bounded by M -|B| 2 /4 . Note that the upper bound j max can now be sent to infinity without any divergence, ensuring that the theory is well-defined in the ultraviolet limit.

The final step is to sum over the fermionic forest. Such a forest can be partitioned into components of cardinal b k , associated to connected blocks of size k, having k sub-vertices.

The number of fermionic edges is then k b k -1, and for each component with k sub-vertices, there are n b k ways to hook a fermionic edge. Moreover, Cayley's theorem (without constraint on the coordinate number) states that the number of trees with v labelled vertices is v v-2 , leading to a contribution n k b k -2 . Finally, because of the constraint k kb k = n, when the number of (sub) vertices is fixed to n, and from Stirling formula:

n ( k b k )-2 ≤ ( k b k )!e n ,
we find:

| ln Z[J, J, λ]| ≤ N n 1 n! {b k } k kb k =n n! k b k !(k!) b k 2 k b k -1 k ( k b k )-2 k n b k × k O(1) π 2 |λ|N cos 2 (ϕ/2) k-1 (4k -4)!! (3k -3)! (2k -1)! M -k 2 /4 b k . (III.123)
Taking into account the bound (III.120) and performing the sum over the {b i }, we get finally:

| ln Z[J, J, λ]| ≤ b≥0 n≥1 π 2 |λ|N cos 2 (ϕ/2) n-1 3 3n n n M -n 2 /4 b .
(III.124)

The power of M ensures that, for M sufficiently large, this factor compensates the bad divergence associated to n n , ensuring analyticity in the interior of a cardioid domain |λ| ≤ ρ N cos(ϕ/2) for small enough ρ. The final check to prove Borel summability follows straightforwardly the Section III A 2.

C. Fermions in dimensions d ≤ 1

In this section, we investigate the Borel summability of quartic fermionic models. Note that the partition function of fermions in d = 0 at finite N is polynomial (and thus analytic) in the coupling constant since the fermions are Grassmann valued. As a consequence, the partition function Z is analytic (but not polynomial): thus, the main non-trivial content of our result lies in the large N limit. We briefly discuss the case d = 0, as an instructive step to extend the proof for d = 1. Moreover, we focus on Majorana fields to anticipate the proof of the constructive bound for the SYK model in the next section, but the extension to Dirac fields is straightforward. To get a well i.i.d real field model, we have to use a physical feature of fermionic fields. Indeed, physically relevant fermionic fields share some internal discrete degrees of freedom, like spin. Assuming that such a spin exists for our fields, we can introduce a new index σ = ± such that ψ i becomes a bi-vector: ψ T i ≡ (ψ i,+ , ψ i,-); and we will denote the components as ψ i,σ . Thus, the total number of degrees of freedom equals 2N in this simplest model; with the following kinetic action:

S kin (Ψ) := i,σ,j,σ ψ i,σ K -1 ij τ σ,σ ψ j,σ , (III.125)
where the matrix τ = τ -1 couples the up and down spins, and we must set τ = σ 1 , the first Pauli matrix. In zero dimension, real fields correspond to Majorana spinors, and τ may be interpreted as a charge conjugation matrix. We can then use the same strategy as for the bosonic case. We can diagonalize the matrix W IJ for the blocks indices I = (ij) and J = (kl), denoting the eigenvalues as ξ 2 L :

J IJ = L ξ 2 L O IL O T LJ .
(III.126)

For d = 0 and K -1 ij = δ ij , the strategy used for the bosonic bound holds. Introducing several intermediate fields, the partition function becomes:

Z(u) = dµ C (ψ)dν(σ) L e √ -u 12 ξ L Ξ L σ (L) , ( III
C(ω) = 1 iω . (III.132)
Note that the model does not require a mass term to be UV regularized because of the boundary condition ψ i = -ψ i (t + β), which excludes the zero-frequency mode ω = 0. The spectrum for fermions is given by formula (II.28):

ω 1 iω = β 2iπ n 1 n + 1/2 = β 2 . (III.133)
As for bosons, bounding the amplitude requires to bound sums like:

f (ω) := ω ,ω ,i 1 ω δ(ω -ω + ω )V ii [χ 1 ](ω , ω ) := Tr B(ω)R[χ 1 ] , (III.134) with B ω ω (ω) := 1 ω δ(ω -ω + ω ) . (III.135)
There is an additional subtlety coming from the fact that the propagator is not real. Thus, redefining the resolvent as:

R -1 ωω = δ ω,ω - u 3β C 1/2 (ω)Σ ωω C 1/2 (ω ) = δ ω,ω + i u 3β 1 √ ω Σ ωω 1 √ ω .
The operator

H := 1 √ ω Σ ωω 1 √ ω , (III.136)
is Hermitian due to the reality condition Σ * ωω = Σ -ω,-ω . As a result, R admits the operator bound: R ≤ cos -1 ϕ 2 , (III.137) with ϕ = arg(u). Note that, due to the sign in front of u in the definition (II.18), the cardioid domain is rotated by π. Because B(ω) is not positive defined, we do not use the same trick as in the previous section to bound f (ω). However, we can use the properties of the operator bound. In particular, let U be the unitary transformation acting on the Hilbert space C N ⊗ L 2 (S 1 , R) which diagonalizes V I , and let r iω be its eigenvalues. f (ω) may be rewritten as:

f (ω) = ω ,i r iω Bi;ω (ω) (III.138)
where Bi;ω (ω) := (U † B(ω)U ) i,i;ω ω . Obviously, Bi;ω ∈ C N ⊗ L 2 (S 1 , R). Indeed:

B(ω) 2 = i,ω B † i;ω Bi;ω ≤   i,ω | Bi;ω |   2 ≤ β 2 2 . (III.139)
As a result, because R is bounded and from definition of the operator norm, we get:

|f (ω)| ≤ R × B(ω) 2 = β 2 cos -1 ϕ 2 . (III.140)
Then, the bound deduced in the previous section for 1d bosons holds, up to the replacements β coth(βm/2)/2m → β/2 and u/2β → u/3β; and taking into account the 1/2 per vertex coming from the substitution e -Tr ln → e 1 2 Tr ln . We thus proved the existence of a finite analytic domain with shape of a cardioid.

D. Disordered Majorana fermions in d = 1: the SYK model

We proved the existence of a finite analyticity domain for the fermionic part of the SYK model with a fixed coupling tensor. In this section, we will complete the proof of the Proposition 2, taking into account the disordered nature of the coupling tensor J ijkl .

In particular, we show that the Gaussian integration preserves the existence of a uniform analyticity domain in the large N limit.

Let us return on the building of the first bound, see equation (III.17). Let T be the tree pictured on the Figure 1, with n vertices, and let O(V [START_REF] Bardakci | Nature of the axioms of relativistic quantum field theory I[END_REF] ) I , O(V 2 ) IJ and O(V [START_REF] Streater | PCT, Spin and Statistics, and All That[END_REF] ) J three connected components such that the amplitude A T reads (we are only interested by the internal index structure of the amplitude, and we forget the sums over frequencies ω):

A T ∝ IJ O(V (1) ) I ξ 2 I O(V 2 ) IJ ξ 2 J O(V (3) ) J . (III.141)
Due to Wick theorem, the computation of the Gaussian integral over J generates "pairings" between ξ 2 I variables, and then adds loops on the tree. There are two sources for ξ 2 I variables: the resolvent and the links, generating three types of pairing: 3. Wick contractions between resolvent and edge variables, as on Figure 2 (c).

O(V (1) )I O(V (2) )IJ O(V (3) )J
Wick contractions of type ( 2) and ( 3) are effective, in the sense that the resolvent is a function of ξ 2 I . To clarify the graphical notations, we will picture the ordinary Wick contractions with dashed-dotted lines and the effective ones with dashed lines on Figure 2. Ultimately, we are interested to deal with the limit N → ∞, and some of the allowed contractions will be discarded in this limit. For instance, the contractions of type (1) between edge variables introduce additional Kronecker delta, which reduces the ranges of the sums over internal indices, and ultimately the final dependence on N of the amplitude A T . The same thing holds for contractions between edge variable and resolvent which are not hooked on the boundary of the corresponding edge, as well as between different resolvents. The contractions of the edge variables which optimize the N dependence of A T are then necessarily between one of the two resolvents hooked on the boundary of the corresponding edge; and the remaining contractions have to be "self loops" between resolvents. As a result, for each edge, there are two allowed contractions: + , (III.142) a typical tree contributing to A T J -the averaged amplitude is pictured on Figure 3. 85 A T J is such a sum over all permutations corresponding to allowed contractions describing in equation (III.142). There are two type of contractions per edges, and n -1 edges.

We denote as T J the corresponding 2 n-1 graphs obtained from T . Therefore, A Along each of the c(v) edges, we can bound the expression using Cauchy-Schwarz inequality.

For k = 1, we define

V (n) I 1 = {I k },k =1,I V I;I 1 ,•••I c(v) c(v) k=2 ξ 2 I k V (m(k)) I k ξ 2 I RI . (III.148)
Then, assuming that we averaged on J only on the connected component V (m(1)) I 1

, we have two configurations, to the two allowed Wick contractions along the edge 1.

However, even in the case when the contraction is on the component V , so that in any case we have the bound:

I 1 V (n) I 1 J ξ 2 I 1 V (m(1)) I 1 J ≤ I 1 |V (n) I 1 | J I 1 |ξ 2 I 1 V (m(1)) I 1 | 2 J .
Recursively, we get: Finally, because of our recursion hypothesis, all the C k have to be of order 1 in N , implying

I 1 V (n) I 1 J ξ 2 I 1 V (m( 1
| A T n+1 J | ∼ √ N .
Therefore, the only changes from the bound for the fermionic part coming from the Gaussian integration are 1) The replacement N ξ 2 0 → ḡ and 2) the combinatorial factor 2 n-1 coming from the Wick contractions. As a result, the averaged amplitude A Tn J admits the for some constant K, with (ρ ) -1 := 24 e/π/βḡ. Once again, the final check to prove Borel summability follows the Section III A 2.

Let [1, . . . , n] be the finite set of points considered above. An edge l between two elements i, j ∈ [1, . . . , n] is a pair (i, j) for 1 ≤ i < j ≤ n, and the set of such edges can be identified with the set of lines of K n , the complete graph with n vertices. Consider the vector space S n of n × n symmetric matrices, whose dimension is n(n + 1)/2 and the compact and convex subset P S n of positive symmetric matrices whose diagonal coefficients are all equal to 1, and off-diagonal elements are between 0 and 1. Any X ∈ P S n can be parametrized by n(n-1)/2 elements X l , where l run over the edges of the complete graph K n . Let us consider a smooth function f defined in the interior of P S n with continuous extensions to P S n itself. The BKAR forest formula states that:

Theorem 2 (The BKAR forest formula) We have

f (1) = F dw F ∂ F f [X F (w F )] (A.1)
where 1 is the matrix with all entries equal to 1, and:

• The sum is over the forests F over n labelled vertices, including the empty forest.

• The integration over dw F means integration from 0 to 1 over one parameter for each edge of the forest. Note that there is no integration for the empty forest since by convention an empty product is 1.

• ∂ F := l∈F ∂ l means a product of partial derivatives with respect to the variables X l associated to the edge l of F.

• The matrix X F (w F ) ∈ P S n is such that X F ii (w F ) = 1 ∀i, and for i = j X F ij (w F ) is the infimum of the w l variables for l in the unique path from i to j in F. If no such path exists, by definition X F ij (w F ) = 0.

• f (λ) admits a Taylor expansion at the origin:

f (λ) = r-1 k=0 a k λ k + R r f (λ), |R r f (λ)| ≤ Kσ r r!|λ| r , (B.1)
for some constants K and σ independent of r.

If f (λ) is Borel summable in λ, then:

B = ∞ n=0 1 n! a n t n (B.

  Dirac matrices γ α with α = 1, • • • , d satisfy the Dirac algebra {γ α , γ β } = 2δ αβ and the condition tr(γ α ) = 0 for d > 0. Dirac matrices are 2 n × 2 n matrices where d = 2n (resp. d = 2n + 1) for even (resp. odd) d meaning that relativistic fermions ψ i

Proposition 1

 1 Let F (u) := ln(Z(u)) be the free energy of a quartic model in d dimensions with generic but fixed tensor coupling and K ij = δ ij . Then, for the relativistic fermionic models in d ≤ 1, non-relativistic bosonic models in dimensions d ≤ 1 and the relativistic bosonic models in d ≤ 2, the tree expansion of F (u) is analytic in the interior of the cardioid domain:

  general strategy for the proofs is to show that the two requirements of Nevanlinna's theorem (recalled in appendix B) hold 81 . We start with bosons and with dimensions d ≤ 1. The case d = 2 is technically more difficult asince it requires an improved version of the LVE, called multi-scale loop vertex expansion (MLVE).

1 2 ijkl

 12 σ kl (pp ) , (III.63) with δ(p, p ) := δ ωω δ(p -p ) and Tr in formula (III.62) means Tr K := i dp K ii (p, p) , (III.64) with the shorthand notation dp := ω dp. Note that the linear terms in σ in the equation (III.62) exactly compensate the first term in the Taylor expansion of the logarithm, in such a way that the effective partition function may be rewritten as: Z(u) = dν(σ) e -Tr ln 2 K[σ] . (III.65) We recall the standard definition ln 2 (1 -x) = ln(1 -x) + x. The notation ln 2 in this expression makes sense because K has the form K =: Id/H, the explicit expression for the operator H being given on equation (III.63). Defining the resolvent R = K -1 , we have once again the resolvent bound: R ≤ cos -1 ϕ 2 (III.66) where ϕ = arg(u), |u| ≤ π. This bound follows from the Lemma 1 and from the same arguments used in the case d = 1 to transform CΣ as C 1/2 ΣC 1/2 . Now, we move on to the jungle expansion. First, we fix our different conventions and notations. A natural choice, taking into account the SO(2) invariance of the Laplacian is to choose a disk: 0 ≤ p 2 ≤ Λ 2 , with the shorthand notations p := (ω, p) and p 2 := ω 2 + p 2 .

,

  with Y ab := Y B(a)B(b) δ jaj b . We thus define: Y p 1 ,...,p k m 1 ,...,m k := of the matrix Y, having the lines p 1 • • • p k and the columns m 1 • • • m k deleted.

Lemma 3

 3 p k (III.82) where the sum runs over the 2 k ways to exchange the upper and lower indices, and k := |F F | is the cardinal of the fermionic forest, and the first product over factors (1-δ jaj b ) implements the hard core constraint. To bound the fermionic contribution, we have the important lemma: Due to the positivity of the covariance Y, for any {m i } and {p i } the minor Y p 1 ,...,p k m 1 ,...,m k defined in (III.81) satisfies: |Y p 1 ,...,p k m 1 ,...,m k | ≤ 1. (III.83)

, with λ 2 :

 2 = -u/2πβ. Using the constraint: m c(m) = 2(|B| -1), with |B| the number of vertices of B, and because of the bound R ≤ cos -1 (ϕ/2), the equation (III.90) satisfies the inequality:

I 2 .

 2 .127) where Ξ L := I=ij O IL ψ i ψ j and dµ C (ψ), dν(σ) are normalized Gaussian measures. Using the standard result on Berezin integrals:2n i dψ i exp Ψ T AΨ = 2 n √ det A , (III.128)we get the following effective bosonic theory:Z(u, J) = dν(χ)e V(χ) . (III.129)Except for the difference of sign in front of the effective interaction V, it is exactly what we obtained for the bosonic fields. In particular, our conclusions on the bound of absolute values for terms involved in the tree expansion of the free energy hold. The only difference comes from the presence of the matrix τ . From Lemma 1, it follows that the norm of the resolvent R(χ) remains bounded as cos -1 (ϕ/2). However, the trace over external indices,"Tr", involved in the definition of the vertex function V(χ) becomes a trace over the whole index structure, including spin indices. Each derivative involved in the definition of the effective vertex introduces a τ insertion. Let v and w be two effective vertices and a link between them. Denoting by tr the trace over spin indices, the trace tr[R (v) τ R (w) τ ] corresponds to the edge . The external index structure with respect to the indices i, j remaining unchanged, the bound given for bosons still holds. Fixing the external indices, and because the operator norm of τ is τ = 1, the traces over spin indices can be bounded as:| tr[R (v) τ R (w) τ ]| ≤ tr |R (v) R (w) | ≤ tr |R (v) | ×tr |R (w) | . (III.130) As a consequence, the trace Tr |R c(v) | involved in the bosonic bounds becomes Tr(tr |R|) c(v) ≤ R tr I 2 × Tr I = R × 2N ; where we introduced a subscript 2 for the 2 × 2 identity matrix Taking into account this new factor 2 coming from traces over internal indices, this implies the existence of a non-vanishing analyticity domain: can be easily extended to Dirac fermions, or for an internal space with dimension larger than 2. Let us move to the d = 1 case. The Fourier representation of the propagator C(ω) of the theory is:

Figure 1 .

 1 Figure 1. A typical tree contributing to the constructive expansion of the fermionic integral.

1 .

 1 Wick contractions between two edge variables represented as the dashed-dotted line on Figure 2 (a).

2 .

 2 Wick contractions between two resolvents, as on Figure 2 (b).

Figure 2 .

 2 Figure 2. Allowed Wick-contractions on the tree T . Contraction between two edge variables (a), between two corners (b) and between edge and corner (c).

Figure 3 .

 3 Figure 3. A typical tree contributing to the averaged amplitude A T J .

T J is bounded by 2 n- 1 ×Lemma 6 IRI J ξ 2 I 2 N

 1622 sup{|A G |, G ∈ T J }, where A G denotes the amplitude for the graph G. The final step to bound A T J is then to bound the amplitudes A G , and we have the following statement: The averaged amplitudes A G , for G ∈ T J must be of order √ N : | A Tn J | = O( √ N ) . (III.143) Proof. To prove this, we proceed recursively on n. Let us consider the simpler tree T 2 with two vertices, on the right of equation (III.142). Writing Rij = R ij δ ij , the amplitude becomes proportional to: | A T 2 J | ∝ RI normalization of the Gaussian integration, the first factor I | RI | 2 J is uniformly bounded as cos -2 ϕ 2 Tr I. For the second factor, the Cauchy-Schwarz inequality leads to: -3 N 2 × N . (III.145) As a result | A T 2 J | = O( √ N ).Now, assuming that the property holds for any tree of size n, we have to prove that it survives after adding a leaf. The amplitude A T n+1 reads: of the tree with n vertices. Denoting as v the vertex to which the leaf is added and as c(v) + 1 its coordination number, V

I 1 ,

 1 we can replace the contracted resolvent with the resolvent of the component V (m(1)) I 1

)) I 1 J 2 J 2 I

 122 ≤ {I k } |V I;I 1 ,•••I c(v). Finally, returning on the sum (III.146), we get: k }|V I;I 1 ,•••I c(v) | J I |ξ RI | 2 J . (III.150)From equation (III.145), the second term into the square-root is of order 1. Moreover, the first term is nothing but a big trace, like for (III.19), as a result:I,{I k } |V I;I 1 ,•••I c(v) | J 1 2 cos -c(v)-1 ϕ 2 × N .(III.151)

  bound:| A Tn J | ≤ K√ N × |u| n (ρ ) n cos -2n ϕ 2 , (III.152)meaning that the expansion in analytic in the interior of the cardioid domain ρ ≤ ρ cos 2 (ϕ/2)

2 )

 2 is an analytic function for |t| < σ -1 which admits an analytic continuation in the strip{z | |Im(z)| < σ -1 } such that |B| ≤ Be t/Rfor some constant B and f (λ) is represented by the absolutely convergent integral:f (λ) = 1 λ +∞ 0 dt Be -t/λ . (B.3)

ACKNOWLEDGEMENTS

We are grateful to Vincent Rivasseau for discussions and comments on the last version of the draft. V.L. is also indebted to Vincent Rivasseau for having introduced him to constructive field theory -and, in particular, to the forest formula -and for his continuous support.

The work of H.E. is conducted under a Carl Friedrich von Siemens Research Fellowship of the Alexander von Humboldt Foundation for postdoctoral researchers.

DATA AVAILABILITY

Data sharing is not applicable to this article as no new data were created or analyzed in this study.

Appendix A: The Brydges-Kennedy-Abdesselam-Rivasseau forest formula

A forest formula expands a quantity defined on n points in terms of forests built on these points. There are many forest formulas, but the BKAR formula seems to be the only one which is both symmetric under permutation of the n points and positive. In this appendix, we only recall the theorem without proof, much additional material may be found in [START_REF] Gurau | The Multiscale Loop Vertex Expansion[END_REF][START_REF] Abdesselam | Trees, Forests and Jungles: A Botanical Garden for Cluster Expansions[END_REF] .

Appendix B: Nevanlinna's theorem

Theorem 3 A series n an n! λ n is Borel summable to a function f (λ) if the following conditions are met:

Borel summable when integrating along with the flow (Landau pole). This reflects the fact that the theory does not exist by itself, but it can be coupled with other fields. For this reason, the structure of the interaction must be slightly modified to make the beta function negative, which can be understood as writing an effective theory. See 87 for a comparison of the vector, matrix, and tensor cases. 78 The method employed in this paper relies heavily on this condition. Indeed, for m 2 < 0, the propagator is not positive definite which does not allow to bound the resolvent as in Lemma 1 or (III.66). Very recently, the Borel summability for scalar theories with m 2 < 0 have been investigated in [START_REF] Serone | The Power of Perturbation Theory[END_REF][START_REF] Serone | λφ 4 Theory I: The Symmetric Phase Beyond NNNNNNNNLO[END_REF] .

21, 261-263 (1980).

82 Indeed, decomposing σ in symmetric and anti-symmetric parts as σ = σ S + σ AS , the kinetic action for the matrix σ writes as

, so that the 48 anti-symmetric part factorizes exactly in front of the path integral, and is canceled by the normalization.

83 Note that the free energy, which is the logarithm of the partition function, is expanded as a sum of trees rather than forests.

84 We used the well known sum n∈Z 1 n 2 +y 2 = π y coth(πy). 85 Note that this restriction is unnecessary, the first slices can be bounded with a simple loop vertex expansion. 86 The notation X J means Gaussian averaging with respect to the tensor coupling J , normalized to 1: 1 J = 1. 87 V. Rivasseau, "Why Are Tensor Field Theories Asymptotically Free?" EPL (Europhysics Letters) 111, 60011 (2015), 1507.04190.