Numerical analysis of the self-energy in covariant LQG

Pietropaolo Frisoni pfrisoni@uwo.ca Dept. of Physics & Astronomy, Western University, London, ON N6A 3K7, Canada

Francesco Gozzini gozzini@cpt.univ-mrs.fr Aix Marseille University, Université de Toulon, CNRS, CPT, 13288 Marseille, France

Francesca Vidotto fvidotto@uwo.ca Dept. of Physics & Astronomy, Western University, London, ON N6A 3K7, Canada Dept. of Philosophy and Rotman Institute, Western University, London, ON N6A 3K7, Canada

Abstract

We study numerically the first order radiative corrections to the self-energy, in covariant loop quantum gravity. We employ the recently developed sl2cfoam-next spinfoam amplitudes library, and some original numerical methods. We analyze the scaling of the divergence with the infrared cutoff, for which previous analytical estimates provided widely different lower and upper bounds. Our findings suggest that the divergence is approximately linear in the cutoff. We also investigate the role of the Barbero-Immirzi parameter in the asymptotic behavior, the dependence of the scaling on some boundary data and the expectation values of boundary operators.

I Introduction

The spinfoam, or covariant, formulation of loop quantum gravity Perez:2012wv , is based on the definition of truncated amplitudes that define the background independent dynamics of quantum spacetime Oeckl:2003vu ; Rovelli2004 in a way similar to the Feynman graphs’ definition of the dynamics of quantum field theory. In a Feynman graph, radiative corrections appear as loops (set of propagators forming a circle) and may give rise to ultraviolet divergences. In a spinfoam, radiative corrections appear as bubbles (set of faces forming a sphere) and may give rise to infrared divergences when the cosmological constant is taken to be zero¹. Here we investigate the first order correction to a single edge of the spinfoam, which is analog to the first order contribution to the self-energy in quantum field theory, hence are referred to as self-energy radiative corrections. We use the version of the amplitude developed in Pereira:2007nh ; Livine:2007vk ; Livine:2007ya ; Freidel:2007py ; Engle2007 ; Engle:2007uq ; Engle:2007wy ; Kaminski:2009fm ; Ding:2010fw generally denoted as the EPRL amplitude. For reviews and notation, see Rovelli:2014ssa ; Rovelli2011c ; Dona:2010cr .

The general definition of the spinfoam amplitudes includes a non-zero cosmological constant Fairbairn:2011aa ; Han:2011nx ; Haggard:2015yda that plays the role of an infrared cutoff and makes the amplitudes finite.

The spinfoam infrared divergences have been studied analytically from a number of different perspectives dittrich_2016_continuum_limit ; Dittrich_decorated_tensor_network ; Bahr_Investigation_of_the_spinfoam ; Bonzom_random_tensor_models ; Benedetti_phase_transition_in_dually ; Carrozza_renormalization_of ; Geloun_functional_renormalization ; Bonzom:2010zh ; Bonzom:2011br ; Baratin_melonic_Phase_Transition . EPRL divergences have been studied in the Euclidean Perini:2008pd ; Krajewski:2010yq and Lorentzian theories Riello:2013bzw ; Dona:2018infrared . In particular the last two references by Riello and Doná provide lower and upper bounds to the degree of divergence of the self-energy, resulting in a wide window. In Riello:2013bzw , Riello gives a detailed analysis of the critical point that contributes to the amplitude in the large-spin scaling limit, based on the techniques introduced in article_Barrett_etal_2010_lorentzian_spinfoam_amplitude , in the case of non-degenerate geometric configurations. He finds a logarithmic divergence in the cutoff K. He then estimates in a final appendix that the contribution of the degenerate sector might be linear in K, hence dominating the total divergence, postponing a detailed analysis of this sector to future studies. On the other hand, in Dona:2018infrared , Doná disregards the interference between the different terms in the sum, and this allows him to estimate an upper bound on the divergence that scales as K⁹. He also proves that the upper bounds derived with his algorithm provide excellent estimates of different kinds of infrared divergences restricted to the 3D model, suggesting that this might be also true in 4D, underlining the need for a complete numerical analysis for future studies.

These two results leave a window of possibilities which spans more than nine powers of the cutoff. In this work, we address the issue numerically. What we find is consistent with Riello’s rough estimate: for small values of the Barbero-Immirzi parameter, the numerical analysis is compatible with a linear divergence in K.

The numerical calculation is possible thanks to the recently developed sl2cfoam-next library Francesco_draft_new_code . Our computations were performed on the clusters Cedar and Graham of Compute Canada (www.computecanada.ca) and on the Mesocentre of Aix-Marseille University (mesocentre.univ-amu.fr). The analysis and visualization of the computed data have been done with Julia bezanson2017julia and Mathematica Mathematica .

In section II we describe the structure of the self-energy amplitude. In section III we recall the calculation of the degree of divergence in topological SU(2) BF model, which can be found exactly. This allows us to compare the numerical results with the analytical ones in the literature. In section IV we define the EPRL amplitude. In section V we describe the results of the numerical analysis of the EPRL self-energy amplitude. Finally, in the appendices we report some further numerical studies and consistency checks.

II Self-energy

Consider two four-simplices joined by four tetrahedra. The boundary of this cellular complex is formed by two tetrahedra. In quantum gravity, this cellular complex can be associated to the transition from a single quantum of space to a single quantum of space, via a splitting into four virtual quanta, which then recombine Oeckl:2003vu ; Rovelli2004 . Notice the similarity with the self-energy Feynman graph in a quantum field theory. This however is only an analogy: there is no notion of energy or radiation involved here. This graph should be understood as a higher order term in the spinfoam expansion.

Figure 1: The two-complex C of the self-energy. All vertices and edges, but only one internal and two external faces are depicted.

The dual of this cellular complex is a two-complex C. It is formed by two pentavalent vertices v_± joined by four internal edges eₐ, with a=1,2,3,4. Each of the two vertices has one other edge e_±, ending on the boundary. The two-complex has four external faces fₐ that end on the boundary, and six internal faces fₐb, dual to the surfaces that separate the four internal tetrahedra. See Figure 1. The six faces fₐb, taken together, form the bubble we are interested in.

The boundary of the two-complex C is the graph Γ formed by the two 4-valent nodes n_± where the edges e_± end (dual to the initial and final quanta of space), joined by four links lₐ that bound the faces fₐ (see Figure 2). The kinematical LQG boundary Hilbert space is therefore H_Γ=L₂ [SU(2)⁴/SU(2)² ]_Γ. The graph Γ is represented in Figure 2.

Figure 2: The boundary graph.

A basis in H_Γ is given by the spinnetwork states |jₐ,i_±⟩, where jₐ are the spins attached to the links lₐ. In the two intertwiner spaces we choose the recoupling basis i_± defined by pairing the links j₁ and j₂ at the node n_±, i.e. the basis that diagonalizes the modulussquare of the sum of the SU(2) generators in the SU(2) representation j₁ and j₂ Bianchi:2011hp ; Bianchi:2011ub ; Freidel:2010aq . The label i indicates the spin of the virtual link in this recoupling basis.

We focus on the subspaces Hⱼ ⊂H_Γ where the four spins are equal jₐ=j, and denote the basis states as

$$|j,i_{\pm}\rangle=|j,i_{+}\rangle\otimes|j,i_{-}\rangle .$$

(1)

We are interested in the EPRL amplitude Rovelli:2014ssa

$$W(j,i_{\pm})=\langle W_{\mathcal{C}}|j,i_{\pm}\rangle$$

(2)

relative to the two complex C. To compute this amplitude numerically, we need an infrared cutoff to bound the sum over the spins associated to the internal faces. We denote this cutoff K, and we define the cutoff amplitude as

$$W(j,i_{\pm},K)=\langle W_{\mathcal{C},K}|j,i_{\pm}\rangle.$$

(3)

We are interested in studying the dependence of this amplitude on K.

III BF amplitude

As a first step, we consider the amplitude

$$W_{BF}(j,i_{\pm})=\langle W^{BF}_{\mathcal{C}}|j,i_{\pm}\rangle$$

(4)

in the topological SU(2) BF theory, where the amplitude can be computed analytically Dona:2018infrared . The BF self-energy amplitude can be written as

$$W_{BF}$$

(5)

where $d_{\lbrace j_{ab}\rbrace }=\prod\limits_{(a,b)}(2j_{ab}+1)$ and $d_{\lbrace i_{a}\rbrace }=\prod\limits_{a}(2i_{a}+1)$, with a,b=1,...4, a≠b. The BF vertex amplitude {15j} _j_AB,i_A is the 15-j Wigner symbol defined contracting the five intertwiners i_A with the ten spins J_AB, with A,B=1,...5, A≠B. See the Appendix C for notation and conventions. A graphical representation of this amplitude that will be useful later on is the following:

$$W_{BF}\qquad {\includegraphics[width=199.16928pt]{correctedbf}}.$$

(6)

Here for completeness we should include in the expression for the amplitude the dimensional factor of the boundary faces and intertwiners, but this is a constant overall factor that can be ignored in the analysis since it does not affect the scaling of the divergence, that is what interests us. In (6), the six internal faces jₐb are highlighted in red. If we choose a recoupling basis for the internal intertwiners as follows

$$W_{BF}\qquad {\includegraphics[width=199.16928pt]{correctedBF2}},$$

(7)

then this picture gives directly the expression of the two BF vertex amplitudes in terms of Wigner 3-j symbols: each intersection of three lines is one such symbols and the lines specify the contractions’ scheme.

Figure 3: Log-log plot of the BF self-energy amplitude as a function of the cutoff K on the internal spins jₐb. The continuous curve represents the function ${}4.8\cdotK^{9}$.

The multiple sum over the spins jₐb and the virtual spins iₐ of the intertwiners are constrained by the Mandelstam identities that need to be satisfied for the Wigner symbols not to vanish. But it is clear from the last picture that the fixed values of the external spins jₐ do not prevent the internal spins jₐb for growing arbitrarily large.

This is the origin of the divergence. To study it numerically, we introduce a cutoff jₘₐₓ=K on the six internal faces jₐb, therefore bounding all the sums by K. That is, we define

$$W_{BF}$$

We have computed this quantity numerically for a fixed value of the jₐ and i_±, for different values of K. In the BF case it is known analytically that the degree of divergence does not depend on jₐ and i_±. The result is shown in Figure 3, where the continuous curve is W_BF=4.8· K⁹.

This result is consistent with the analytical result Dona:2018infrared , which gives

$$W_{BF}\left(j,i_{\pm};K\right)\propto K^{9}.$$

(8)

We briefly recall how to derive this result. Consider the simple case where all external spins are zero. The BF amplitude can be expressed in Fourier transform as a multiple integral over SU(2), given by one integration per link and one delta function per face. In this case we have

$$W_{BF}=\int\mathrm{d}g_{a}\prod_{a\gt b}\delta\left(g_{a}g_{b}^{-1}\right).$$

It is easy to see that one integration is redundant, and the other three can be used to eliminate three delta functions, leaving

$$W_{BF}=\delta\left(\mathbb{1}\right)^{3}.$$

(9)

To regularize this expression, one can write the delta function as a sum over characters χʲ of irreducible representations of SU(2) and introduce an infrared cutoff K, as follows:

$$\delta\left(U\right)=\lim_{K}\delta_{K}\left(U\right)=\lim_{K}\sum_{j}^{K}(2j+1)\chi^{j}\left(U\right) .$$

Then, easily, δ_K(1)≈ K³ and (8) follows from (9).

The numerical analysis in Figure 3 shows that the asymptotic scaling of the self-energy BF divergence (8) is reached for very low values of the infrared cutoff K. The computation involved in producing Figure 3 is nowadays trivial and can be easily performed on a laptop thanks to very efficient routines for computing Wigner (3,6,9)-j symbols johanssonFastAccurateEvaluation2015 ; Efficient_Storage_Scheme . However, we note that the computational complexity grows for a large power of the cutoff index K, approximately equal to the number of internal faces of the foam. We plot in Figure 4 the total number of combinations of internal spins that must be computed for each value of the cutoff K, till K=20. We can see that the triangular inequalities reduce the complexity from K⁶ to effectively ∼ K⁵.⁶. This still gives a large number of configurations to consider, and provides a considerable numerical challenge in the EPRL case, whose vertex amplitude is much more computationally complex than the BF one, as we show in the next section.

Figure 4: Log-log plot of the total number of different configurations of internal spins jₐb to be summed over for increasing values of the cutoff K. The power-law fit is ∼ 1.7· K⁵.⁶.

IV EPRL

Let us now come to the quantum gravity theory. The EPRL vertex amplitude is defined using a restriction of the unitary irreducible representations in the principal series of SL(2,ℂ) Engle:2007wy ; Engle:2007uq . In order to evaluate it numerically, a more tractable formulation has been derived in Speziale2016 and discussed in Dona2018 ; article:Dona_etal_2019_numerical_study_lorentzian_EPRL_spinfoam_amplitude . In this formulation, the vertex is expressed as a sum over SU(2) 15-j symbols weighted by one “booster function” B4 per edge (see the Appendix C.2 for explicit formulas). The amplitude of a vertex bounded by spins j_f and intertwiners iₑ can thus be written as a sum

$$V\left(j_{f},\,i_{e},\,\gamma\right)=\sum_{\substack{l_{f}\geq j_{f}\quad k_{e}}}d_{\lbrace k_{e}\rbrace }\left(\prod_{e}B_{4}(j_{f},l_{f};i_{e},k_{e};\gamma)\right)\lbrace 15j\rbrace _{l_{f},k_{e}}$$

(10)

with virtual intertwiners kₑ and virtual spins l_f≥j_f, e=2...5. One booster function must be replaced by the identity, corresponding to the suppression of a redundant SL(2,C) integration². Notice that we keep the dependence on the Barbero-Immirzi parameter explicit. In this expression the infinite sum over l_f, called the “shell expansion”, can be shown analytically to converge, but in view of the numerical calculation it is convenient to truncate it at a finite value defining

The 15-j symbol in (10) depends also on the 4 spins which label the gauge fixed edge, even if in the analytic notation it is not emphasized

$$V\left(j_{f},\,i_{e},\gamma\right)=\lim_{\Delta l\to\infty}V\left(j_{f},\,i_{e},\gamma;\Delta l\right)$$

(11)

where

V\left(j_{f},\,i_{e},\gamma;\Delta l\right)=\qquad {}=\sum_{\substack{j_{f}\leq l_{f}\leq j_{f}+\Delta\quad k_{e}}}d_{\lbrace k_{e}\rbrace }\left(\prod_{e}B_{4}(j_{f},l_{f};i_{e},k_{e};\gamma)\right)\lbrace 15j\rbrace _{l_{f},k_{e}}\qquad {}=\sum_{\substack{j_{f}\leq l_{f}\leq j_{f}+\Delta\quad k_{e}}}d_{\lbrace k_{e}\rbrace }\raisebox{-56.9055pt}{\includegraphics[width=142.26378pt]{vertexAmplitude}}.

(12)

As before, we have neglected dimensional factors attached to the boundary intertwiners iₑ and to the boundary spins j_f. In the following, we will refer to the cutoff Δ l as the number of shells. Even if the convergence of the amplitudes a function of Δ l is assured Dona2018 ; article:Dona_etal_2019_numerical_study_lorentzian_EPRL_spinfoam_amplitude , it is not possible to have a unique prescription to set the optimal Δ l to get an acceptable convergence, since it depends on the details of data such as the face spins j_f and the Barbero-Immirzi parameter. In addition, the convergence depends strongly on the structure of the 2-complex, and there is no known general procedure to estimate the error made in truncating the sum over the auxiliary spins. In the next Section we describe a technique to extrapolate the limit Δ l→∞ for the self-energy.

With the form of the vertex described above, we can write the self-energy amplitude explicitly, in a form suitable for the numerical analysis. This gives

$$W\left(j,i_{\pm},\gamma;K,\Delta l\right)=\sum\limits_{\substack{j_{ab}\leq K\quad j_{ab}\leq l_{ab}\leq j_{ab}+\Delta l\quad i_{a},k_{q\pm}}}d_{\lbrace j_{ab}\rbrace }d_{\lbrace i_{a}\rbrace }d_{\lbrace k_{q\pm}\rbrace }\times\qquad {} \includegraphics[width=227.62204pt]{EPRLintertwiners},$$

(13)

with q=1...4. Compared to the diagrammatic expressions used for the EPRL vertex (12), we have rearranged some links to emphasize the structure of the contractions between the vertices, which defines the spinfoam associated with the triangulation described in Section II.

The EPRL self-energy amplitude (IV) is the main object of the analysis performed in this paper. The interested reader can derive the complete analytical expressions by comparing the graphical counterparts with the definitions given in Appendix C.

When the two boundary intertwiners i₊ and i₋ have the same value i, the EPRL self-energy amplitude simplifies as

$$W\left(j,i,\gamma;K,\Delta l\right)= \sum_{\substack{j_{ab}\leq K\quad i_{a}}}d_{\lbrace j_{ab}\rbrace }\,d_{\lbrace i_{a}\rbrace }V(j,j_{ab};i,i_{a};\gamma;\Delta l)^{2}$$

(14)

since the two vertex amplitudes entering the sum over the internal faces jₐb are identical. Except for the coherent state (Appendix A.1) and the boundary observable analysis (Appendix A.2), we shall always consider the two boundary intertwinwers i_± to have the same value. In the following we use i=(0,0). Result with different values are substantially identical, differing only by a slightly overall shift of the points, as we show in Appendix A.1.

Numerically, the vertex amplitude V(j,jₐb;i,iₐ;γ;Δ l) in (14) can be efficiently computed using the recently developed library sl2cfoam-next Francesco_draft_new_code given the list of boundary spins (jᵢ=j,jₐb) and the number of shells Δ l. The squared sum over the internal intertwiner indices iₐ can be implemented by contracting the resulting vertex tensor with itself. All the configurations of internal spins are automatically distributed by the library across the available cluster nodes, and each node parallelizes the computation of the various shells over the available local cpus. The particular code employed for this work is available at the repository self_energy_repo .

From (12) it is clear that the computation of a single EPRL vertex amplitude is much more complex than the corresponding BF vertex amplitude with the same boundary spins and intertwiner indices. In fact, a single EPRL amplitude requires the computation of roughly (Δ l+1)⁶ different BF amplitudes and (Δ l+1)³ booster functions, for each different set of intertwiners iₐ that bound the vertex (whose number is ∼(2j+1)⁵ if the boundary spins are of order ∼ j). In the present case, a different EPRL vertex tensor (i.e. a bundle of amplitudes at fixed boundary spins with running boundary intertwiner indices) must be computed for each configuration of the internal spins jₐb, whose number as a function of the cutoff K has been estimated in the previous Section (see Figure 4). From these simple estimates, we see that the numerical complexity of computing (14), using our methods, grows quickly when any of the parameters K, Δ l or j increases, roughly with a power law between 5 to 6 powers of the increasing parameter. We provide in the next Section some data for the resources employed in a computation of (14) up to K=10 and Δ l=20.

V Divergence Analysis

We have studied the amplitude (14) numerically and considered the following questions:

1. What is the asymptotic scaling of the EPRL divergence?
2. What is the dependence on the Barbero-Immirzi parameter?
3. At fixed boundary spins and Barbero-Immirzi parameter, does the scaling change using boundary spinnetwork states or intrinsic coherent states?
4. What are the configurations which contribute most to the divergence itself?

The first two questions are addressed in the main text, the third in Appendix A.1 and the fourth in Appendix B. The main difficulty towards answering these questions is the limit Δ l→∞ that cannot be taken numerically. To address this problem, we have relied on a property of the dependence of the amplitude on Δ l at fixed K, that has emerged from the numerical analysis itself.

Figure 5: The EPRL divergence (14) as a function of the cutoff parameter K, for a number of shells Δ l=0 to 10. Darker colors correspond to larger Δ l. A portion of the plot around K=5 is highlighted to show the convergence properties of the shells approximation. Left: γ=0.1. Right: γ=1.0.

Extrapolation scheme

We show in Figure 5 the amplitude (14) as a function of the cutoff K, for various values of the number of shells Δ l and for γ=0.1 and γ=1.0. From a simple qualitative analysis we infer that the convergence of the shell approximation is faster for a reduced bulk spins cutoff, while it becomes slower as K increases. The highest curve is a lower bound to the amplitude as a function of K, which is recovered in the limit Δ l→∞. Convergence in the number of shells appears to depend on γ, with the case γ=1.0 converging faster.

The convergence of (14) in the parameter Δ l can be extrapolated from the data reported in Figure 5. To this aim, we have computed the ratios of the differences between adjacent curves of Figure 5, at fixed K, as a function of Δ l. That is, we have studied the function:

$$f(\gamma,K,\Delta l)\equiv\frac{W(\gamma,K,\Delta l+2)-W(\gamma,K,\Delta l+1)}{W(\gamma,K,\Delta l+1)-W(\gamma,K,\Delta l)} ,$$

(15)

for each fixed value of j and i. The result is shown in Figure 6 for K≥5 and γ=1.0.

Figure 6: Function (15) for γ=1.0. The curves appear to tend to a constant value c_[K,γ].

It appears from Figure 6 that the function f(γ,K,Δ l) tends to a constant for larger enough values of K, as a function of Δ l (however, see also the more complete analysis for the case γ=0.1, that is studied in a following section). Since (15) is the ratio between successive terms of a series, this means that the sum over the parameter Δ l can be well approximated by a geometric series when K is large. Defining

$$c_{K,\gamma}\equiv f\left(\gamma,K,N-2\right)$$

(16)

where N is the largest Δ l parameter that has been reached numerically, our assumption implies for Δ l≥N

$$W(\gamma,K,\Delta l+1)-W(\gamma,K,\Delta l)\approx\left(c_{K,\gamma}\right)^{\Delta l} .$$

(17)

We thus approximate the divergence in the limit Δ l→∞ as

$$W(\gamma,K) \equiv\lim_{\Delta l\to\infty}W(\gamma,K,\Delta l)\qquad {} \approx W(\gamma,K,N-1)+\frac{W(\gamma,K,N)-W(\gamma,K,N-1)}{1-c_{K,\gamma}},$$

(18)

where we used the elementary limit of the geometric series. We note that the approximation improves when (i) more shells can be computed numerically, i.e. the index N increases, and (ii) the ratios (15) are closer to true constants in Δ l. We applied this extrapolation scheme to the amplitudes computed numerically for different values of γ. In particular, in the case γ=0.1, where the convergence in Δ l is the slowest, hence the approximation is the least accurate, we pushed the Δ l parameter to very large values in order to study the accuracy of our extrapolation scheme and obtain very precise results for a physically relevant case.

Case γ=0.1

Figure 7: The complete analysis for the case γ=0.1 till N=Δ lₘₐₓ=20. Left: All the curves at various Δ l with the extrapolation from Δ l=20. Right: The extrapolations from Δ l=10 and Δ l=20 with the corresponding linear fits. The slopes of the fits are 8.1· 10⁻⁶ and 8.6· 10⁻⁶, respectively. The standard errors of the fits, obtained with Julia’s LsqFit package, are 2.3· 10⁻⁸ and 4.6· 10⁻⁹, respectively.

Figure 8: Function (15) for γ=0.1 plotted till N=Δ lₘₐₓ=20. It is evident that the approximation f≈const for Δ l>N provides a lower bound to the exact limit Δ l→∞.

In Figure 7 we show the amplitude (14) for the case γ=0.1, with the computed data for Δ l≤20 and the extrapolated curves from N=10 and from N=20. From these plots we can verify that the approximation gives results that are already accurate starting with low values of N, since the difference between the two extrapolations is small, around 5% for K=10. Plotting the ratios (15) till very large values of Δ l, as we did in Figure 8, we can see that for the case γ=0.1 the approximation f(γ,K,Δ l)≈const is more accurate for larger values of K and less accurate for smaller cutoffs. This is not an issue since convergence in Δ l is reached sooner for small cutoffs (since the ratios are smaller, the series converges faster). Also, Figure 8 shows that the extrapolation is actually a lower bound to the true amplitude, although very close to exact result.

Importantly, the qualitative behavior of the divergence, either extrapolated from N=10 or N=20, is clearly linear. In Figure 7 we show the fit of the extrapolated curves with straight lines. It is clear that the extrapolated points fall exactly on a straight line for K≳ 5 in both cases, with the better approximation at N=20 being also closer to linear behavior than the one at N=10. We conclude in this case that the overall divergence is very close to linear in the cutoff parameter K.

A few comments are in order. First, our results are in accordance with the rough analytical estimate of Riello Riello:2013bzw . It is important to notice that his main analysis focuses on what Riello calls the non-degenerate sector of the divergence, and disregards the contributions from degenerate configurations when the 4-simplices dual to the two spinfoam vertices v_± have non-maximal dimension. He found that this sector contributes with a leading factor of log K to the divergence. However, in Appendix C of Riello:2013bzw the remaining degenerate sector is estimated to contribute with a dominant factor of K to the overall divergence, which is what our numerical analysis shows. Our result also clearly shows that the interference effects neglected in Dona:2018infrared result in a strong suppression of the divergence, of approximately 8 powers of K.

Second, our numerical analysis does not completely excludes the possibility that the behavior for very large values of K is different from what is inferred in the regime K≲ 10. There are however multiple arguments supporting our conclusion: (i) our findings are compatible with past analytical estimates Riello:2013bzw ; (ii) from Figure 7 we see that improving the approximation also improves the matching with the linear fit, which appears perfect for as much as about 12 data points; (iii) we can compare the EPRL divergence with the BF case, where the asymptotic power-law divergence is already manifest for spins of order ∼ 10 (see Figure 3) and analytically it is seen to be independent from the details of the chosen boundary data; (iv) recent numerical investigations Francesco_draft_new_code ; asanteEffectiveSpinFoam2020 strongly support the hypothesis that the semiclassical (i.e. asymptotic) regime of theory is reached for relatively small spins if the Barbero-Immirzi parameter γ is small (e.g. around 0.1). This appears to be related to the slower frequency of asymptotic oscillations of the amplitude whenever γ is small asanteEffectiveSpinFoam2020 . Hence, we can assume that we are looking at the asymptotic divergence for K≳ 5, in this case. This is also supported by our analysis for larger values of γ, which do not appear to converge to a fixed (γ-independent) asymptotic behavior in the considered range of the cutoff, as we show in the next paragraph.

Figure 9: Fit of the amplitude $W(\gamma,K)$ extrapolated via equation (18) with N=10 for decreasing values of the Barbero-Immirzi parameter γ≠0.1 and with N=20 for γ=0.1. The fit is made cutting the values of $W(\gamma,K)$ for K≤4 and superimposing the resulting curve with all the extrapolated points.

As a final comment for the γ=0.1 case, we report some technical details about the hardware resources that were employed. The computation of the divergence at K≤10 and Δ l=20 ran on 640 cpus for ∼ 60 hours, for a total of ∼ 40000 cpu-hours. The whole computation at γ=0.1 from Δ l=0 to 20 took about 125000 cpu-hours. Recall that the computational complexity of the EPRL model scales as (Δ l+1)⁶, thus the plots at largest number of shells took most of the allocated cpu-time. Symmetries in the internal spins jₐb peculiar to this case (j=1/2,i_±=0) have been implemented in the code to reduce the computational time of about a factor ∼ 6.

Case γ≳ 0.1

The extrapolation method outlined above allowed us to investigate also larger values of the Barbero-Immirzi parameter γ without replicating the already substantial computational effort required for the case γ=0.1. We have studied the scaling of the asymptotic divergence for 9 values of γ, ranging between 0.1 and 10. The result is shown in Figure 9, for which we used N=10 in equation (18). Notice that the extrapolation scheme is more effective in the cases γ>0.1 than in the case γ=0.1, as the curves of the ratios given by the function (15) become closer to constants as γ increases (compare Figure 6 with Figure 8). The data of Figure 9 are fitted with a function:

$$W(\gamma,K)=a+bK^{c}$$

(19)

where a,b,c are real coefficients. The best fit for the values of these parameters are shown in Table 1.

Table 1: Fit coefficients table

Table 1 shows that the amplitude is suppressed for large γ. The divergence is increasingly well fitted by an approximately linear scaling as the Barbero-Immirzi parameter decreases, while for γ≳ 1 the curve is approximated by a quadratic function in the range K∈[0,10]. We interpret this apparent dependence of the scaling on γ as a sign that the asymptotic or semiclassical regime is not reached for K≤10 when γ≳ 1. A possible educated guess is that the non-linear behavior for large γ could correspond to the initial non-linear part of the curve seen also at small γ, namely that the dependence could become linear for any γ at sufficiently high K. This intriguing hypothesis could have interesting consequences for a renormalization procedure. We leave the accurate testing of the large-γ sector for future works. As the computational complexity dramatically increases when K increases, it is likely that different numerical methods or approximations would be needed in order to test this sector, which is interesting from a theoretical point of view but however does not appear interesting physically, given the current status of the theory.

VI Conclusions

In this work we applied new computational techniques to the study of the infrared divergence of the spinfoam self-energy graph in the EPRL model. We have computed the divergent amplitudes with a running cutoff under an approximation (the “shells approximation”) which we then lifted by introducing an extrapolation scheme suited for the case considered. We tested our assumptions by refining the approximation up to a large number of shells, using considerable computational resources, for the most interesting case with the Barbero-Immirzi parameter γ=0.1. We also explored the divergence at other values of γ and investigated the dependence and reliability of our results by using different boundary data and by computing geometrical boundary observables.

Our findings considerably refine the upper and lower bounds given in the previous literature Riello:2013bzw ; Dona:2018infrared . We have shown that the numerical evidence clearly points to a linear divergence in the infrared cutoff parameter. On one hand, our result shows that the divergence scaling is much less severe than what expected by approximate power-counting arguments of BF and Euclidean theories Perini:2008pd ; Krajewski:2010yq or upper-bound estimates to the Lorentzian theory Dona:2018infrared . On the other hand, our result also shows that the total divergence is of higher order than the analytical estimate of the so-called “non-degenerate” sector Riello:2013bzw and thus provides the first strong evidence that virtual “degenerate” (i.e. of non-maximal dimension) configurations effectively dominate the first order quantum corrections to the bare spinfoam propagator. This is a point that was hinted at in previous works (see especially Appendix C of Riello:2013bzw ) but left unverified.

As in ordinary quantum field theories, understanding quantum corrections to the bare spinfoam propagator is tightly related to the issue of renormalization, which is still an active area of research in the context of spinfoam and group field theories. Our work fills one of the very first step towards this program by providing a precise estimate of the degree of divergence. Our results about the divergence at various values of γ also suggest that a possible renormalization scheme might involve the running of the Barbero-Immirzi parameter. The numerical testing of this hypothesis is left for future work, as the possible generalization of our extrapolation scheme to different divergent graphs.

This is one of the first works in the field of covariant loop quantum gravity to specifically exploit high-performance codes on a computer cluster. The computations accomplished in this paper would not have been possible using a single or a few machines. We believe that complex computational projects will prove to be more and more fundamental to the advancement of the field in the coming years and we hope that our work and our codes will provide a useful and encouraging ground for progressing in this direction. In this regard, the application of MCMC methods in covariant LQG to study the quantum regime in spinfoams composed of multiple vertex amplitudes glued together will appear soon Cosm_project .

Acknowledgments

We thank Pietro Donà for numerous contributions, especially during the early stages of this work. We are grateful to Carlo Rovelli for discussing at length this manuscript and providing many suggestions that improved it. The Centre de Calcul Intensif d’Aix-Marseille and the Shared Hierarchical Academic Research Computing Network (SHARCNET) are acknowledged for granting access to their high-performance computing resources. We thank in particular the Compute/Calcul Canada staff for the constant support provided with the Cedar and Graham clusters. This work was supported by the Natural Science and Engineering Council of Canada (NSERC) through the Discovery Grant ‘Loop Quantum Gravity: from Computation to Phenomenology’. We acknowledge support also from the QISS JFT grant 61466. FV’s research is supported by the Canada Research Chairs Program. We acknowledge the Anishinaabek, Haudenosaunee, Lūnaapéewak and Attawandaron peoples, on whose traditional lands Western University is located.

***

Appendix A Consistency checks

A.1 Divergence with coherent boundary states

In this Section we test whether the results presented in Section V depend on the chosen boundary data. The basis chosen diagonalizes one of the two dihedral angles of the boundary tetrahedra, thus leaving the other maximally spread. Would the results change using more semiclassical boundary state? To answer this question, we have repeated the computation using intrinsic coherent intertwiners Rovelli:2014ssa as boundary states.

The definition of the intrinsic semiclassical intertwiners is based on the relation between polyhedra and SU(2) invariants Barrett:2009gg ; article:Dona_etal_2018_SU2_graph_invariants and coherent intertwiners Livine:2007vk . A 4-valent coherent intertwiner is defined by the group averaging:

$$||j_{a},\vec{n}_{a}\rangle\equiv\int dg\bigotimes_{a=1}^{4}g\ket{j_{a},\vec{n}_{a}}=\sum\limits_{i}c_{i}(\vec{n}_{a})|j_{a},i\rangle,$$

(20)

where $\ket{j_{a},\vec{n}_{a}}$ are SU(2) coherent states, g∈ SU(2) The vectors n⃗ₐ can be parametrized as n⃗≡(sinΘcosΦ, sinΘ sinΦ,cosΘ). The decomposition of a 4-valent coherent intertwiner, in the recoupling channel i with outgoing links, reads

c_{i}(\vec{n}_{a}) \equiv\bra{j_{a},i}j_{a},\vec{n}_{a}\rangle\qquad {}=\sum_{m_{a}}\left(\begin{array}{cccc}j_{1} j_{2} j_{3} j_{4}\qquad m_{1} m_{2} m_{3} m_{4}\end{array}\right)^{(i)}\prod_{a=1}^{4}D^{j_{a}}_{m_{a},j_{a}}(\vec{n}_{a})\qquad {}=\raisebox{-28.45274pt}{\includegraphics[width=113.81102pt]{SingleCoherent}} ,

(21) (24)

where Dʲₘ,ⱼ(n⃗)=Dʲₘ,ⱼ(ϕ,Θ,-ϕ) is the SU(2) Wigner matrix parametrized as:

$$D^{j}_{m,j}(\phi,\Theta,-\phi)=e^{-im\phi}d^{j}_{m,j}(\Theta)e^{ij\phi} ,$$

with dʲₘ,ⱼ as a small Wigner matrix book:varshalovic .

In order to be interpreted as quantum polyhedra, the boundary tetrahedra of the spinfoam must satisfy the closure condition:

$$\sum_{f}\Lambda_{f} \vec{n}_{f}=0$$

(25)

where Λ_f is the eigenvalue of the area operator associated with each face of the tetrahedron.

An example of the numerical analysis of the amplitude obtained by contracting the amplitude (IV) with coherent states (20), for increasing boundary spin values, is shown³ in Figure 10, compared with boundary spinnetwork states. The scaling of the self-energy amplitude turns out to be unchanged whether using boundary spinnetwork states |j,i_±⟩ or coherent states ||j,n⃗_± f⟩, apart for a global factor that multiplies all curves. Also, using coherent states, the scaling does not seem to depend on the orientation of the normals n⃗_± f.

normals n⃗ₐ satisfy condition (25) and we considered the absolute value in order to obtain real numbers

Figure 10: Self-energy amplitude (IV) with fixed intertwiners (first two panels) versus coherent (third panel) boundary data (20). Using coherent intertwiners, we find no difference in scaling by changing the orientation of the normals n⃗ₐ.

A.2 Boundary observables

Another check of consistency of the geometrical picture and the numerical calculation can be obtained by looking at the geometry of the two boundary tetrahedra. Their geometrical properties are jointly fixed by choosing the boundary data and by the dynamics. Here we show the resulting expectation values of some geometrical boundary operators, that is:

$$\langle O\rangle=\frac{\langle W|O|\Psi\rangle}{\langle W|\Psi\rangle}$$

(26)

where the bra W contains the propagator, namely the dynamics, while the ket Ψ turns out to be the tensor product of the in and out states of the LQG Hilbert space. With the term ‘propagator’ we refer to the square matrices (they are such since the self-energy triangulation has two boundary tetrahedra) in which the element a,b corresponds to the amplitude (IV) with a=i₊,b=i₋. Computing the expectation values (26) of geometric operators describing boundary tetrahedra is a significant numerical investigation, especially considering that computations carried out with the EPRL model are still in their primordial stages Dona2019 ; article:Dona_etal_2019_numerical_study_lorentzian_EPRL_spinfoam_amplitude ; Gozzini_primordial ; Spinfoam_on_a_Lefschetz_thimble ; han2021complex ; frisoni2021studying , and, as far as we know, at present time there are no such numerical computations with spinfoams containing a bubble.

Since in the calculation of (26) the dimensional factor of the boundary intertwiners is relevant, in this section it is convenient to define a ‘normalized’ amplitude as:

$$W_{N}(j,i_{\pm},K)\equiv W(j,i_{\pm},K)\sqrt{d_{\lbrace i_{\pm}\rbrace }} ,$$

(27)

where $d_{\lbrace i_{\pm}\rbrace }=\prod\limits_{\pm}(2i_{\pm}+1)=(2i_{+}+1)(2i_{-}+1)$.

A.2.1 Angles

The shape of the tetrahedra in twisted geometry is measured by the angle operator:

$$A_{ab}|j,i_{\pm}\rangle=\cos(\theta_{ab})_{i_{\pm}}|j,i_{\pm}\rangle$$

(28)

which is interpreted as the cosine of the external dihedral angle between the faces a and b of the tetrahedron defined on the node associated with the intertwiner i_±. The spin-network basis states diagonalize the dihedral angle θₐb between faces a and b. The equation for the dihedral angle cos(θₐb) in terms of intertwiner spin i_± reads:

$$.\cos(\theta_{ab})_{i_{\pm}}=\frac{i_{\pm}(i_{\pm}+1)-j_{a}(j_{a}+1)-j_{b}(j_{b}+1)}{2\sqrt{j_{a}(j_{a}+1)j_{b}(j_{b}+1)}}.$$

(29)

We consider the expectation value (26) of the angle operator (28) in any of the two (equal) boundary regular tetrahedra of the triangulation using the spinnetwork state. According to the recoupling basis (j₁,j₂), we focus on the angle between faces 1 and 2. The expectation value on node + can be computed as:

$$\frac{\langle W|A_{12}|W\rangle}{\langle W|W\rangle}=\frac{\sum\limits_{i_{\pm}}\left[W_{N}(j,i_{\pm},K)\right]^{2}\cos(\theta_{12})_{i_{+}}}{\sum\limits_{i_{\pm}}\left[W_{N}(j,i_{\pm},K)\right]^{2}},$$

(30)

where we considered the case in which all boundary spins are equal to j. Carrying out the numerical computation of (30), we obtain a value that is in agreement with the geometric value of the external angle of a regular tetrahedron up to the tenth significant digit, as shown in Figure 11. We repeated the calculation for increasing values of boundary spins j up to j=5, also varying Δ l, finding exactly the same value and the unaltered precision.

Figure 11: Dynamic expectation value of the cosine of the external dihedral angle operator (30). We find an excellent agreement with the value -0.3.

This is consistent with the fact that we are looking at the only angle which is completely sharp, while the others turn out to be spread.

A.2.2 Volumes

The volume operator of a tetrahedron reads:

$$V_{abc}=\frac{\sqrt{2}}{3}\left(8\pi G\hbar\gamma\right)^{\frac{3}{2}}\sqrt{\left|\vec{J}_{a}\cdot\vec{J}_{b}\times\vec{J}_{c}\right|} ,$$

(31)

where J⃗ₐ is the vector of SU(2) generators on link a. The operator (31) is not diagonal in the spinnetwork states. The dynamic expectation value (see Rovelli:2014ssa for the derivation) on node + turns out to be:

$$\frac{\langle W|V|W\rangle}{\langle W|W\rangle}=\qquad {}=\frac{\sum\limits_{i_{\pm},i^{\prime}_{+}}\sqrt{d_{\lbrace i_{+}\rbrace }d_{\lbrace i_{+}^{\prime}\rbrace }}W_{N}(j,i_{\pm},K)W_{N}(j,i_{+}^{\prime},K)v\left(i_{+},i_{+}^{\prime},j\right)}{\sum\limits_{i_{\pm}}\left[W_{N}(j,i_{\pm},K)\right]^{2}}$$

(32)

where v(i₊,i₊′,j) is the matrix element:

$$v\left(i_{+},i_{+}^{\prime},j\right)=\frac{\sqrt{2}}{3}\left(8\pi G\hbar\gamma\right)^{\frac{3}{2}}\sum\limits_{v}\langle j,i_{+}^{\prime}|v\rangle\sqrt{|q_{v}|}\langle v|j,i_{+}\rangle,$$

(33)

in which qᵥ are the eigenvalues of the volume, while |v⟩ corresponds to the eigenvector relative to the eigenvalue qᵥ. In Figure 12 we report the numerical values obtained from (32). We see that the scaling of the expectation value of the boundary volume as a function of the spin j corresponds to that actually existing between the volume of a regular tetrahedron and the area of one of its faces, that is V∝ j³⁄². As for the angle operator, we repeated the calculation by varying the number of shells, finding the same trend.

Finally, both for the angle operator and for the volume, we repeated the calculation with boundary coherent states (20) verifying that the expectation value is dominated by the classical contribution.

The picture that emerges is that the boundary observables have an expectation value that coincides with the classical one, and are not affected by the presence of divergence in the bulk of the spinfoam. Namely, even if the amplitude (IV) is divergent, the expectation values (30), (32) are finite and fully consistent with the geometry of a regular tetrahedron.

Figure 12: Dynamic expectation value of the volume operator (32). We neglected all the constant factors in (32) and (33). The volumes scale as j³⁄².

Figure 13: The largest individual contributions to the divergent amplitude (14) come from “semiclassical” bulk spins configurations with j̅∼ j and low σ. Triangular inequalities result in a triangular shape in the plane (j̅,σ).

Appendix B Internal contributions

We add here a tentative analysis of the kind of configurations that contribute to the divergence. For this aim, it is convenient to rewrite equation (14) as

$$W\left(j,i,\gamma;K,\Delta l\right)=\sum\limits_{j_{ab}\leq K}\omega(j_{ab};j,i,\gamma,\Delta l).$$

(34)

Now, consider the average and the standard deviation of the sextuplet of internal spins {jₐb}:

$$\bar{j}(j_{ab})=\frac{1}{6}\sum_{a\neq b}j_{ab}, \sigma(j_{ab})=\sqrt{\frac{1}{6}\sum_{a\neq b}(j_{ab}-\bar{j}_{ab})^{2}}$$

(35)

and define

$$w(\bar{j},\sigma;j,i,\gamma,K,\Delta l)=\sum\limits_{\substack{\bar{j}(j_{ab})=\bar{j}\quad \sigma(j_{ab})=\sigma}}\omega(j_{ab};j,i,\gamma,\Delta l).$$

(36)

That is, we write the contributions to the divergence as a function of the average and standard deviation of the internal spins, adding the contributions with different sextuplets that have the same average and standard deviation. The function 100-|logw(j̅,σ;j,i,γ,K,Δ l|) is plotted for fixed j,i,γ,Δ l and various values of K in Figure 13, where the factor 100 is only for visual purposes. Increasing Δ l has the effect of increasing the value of each point, while increasing K also implies adding more configurations, as shown in Figure 4. Figure 13 indicates that the largest individual contributions to the amplitude (14) come from “semiclassical” bulk spins configurations with j̅∼ j and low σ, as one might expect. While the divergence comes from the proliferation of smaller contribution with arbitrary high spin and deviations.

We also note that the divergence of the self-energy, manifestly shown in Figure 9, occurs in the limit Δ l→∞ only. In fact, in Dona:2018infrared it is shown that the sum over the internal spins of the amplitude (14) is convergent in the simplified model Δ l=0. Even if an analytical proof is not yet available, it seems likely that the amplitude has a finite limit K→∞ for each finite Δ l, since the internal spins configurations in the limit become more and more suppressed. Since also the limit Δ l→∞ at fixed K is finite, we can schematize the role of the two limits with the following four node diagram:

$${W(K,\Delta l)}\quad {\text{divergence}.}\quad \scriptstyle{K\rightarrow\infty}\quad \scriptstyle{\Delta l\rightarrow\infty}\quad \scriptstyle{\Delta l\rightarrow\infty}\quad \scriptstyle{K\rightarrow\infty}$$

Appendix C Conventions and basic formulas

We collect here the basic definitions needed for the calculation in this paper.

C.1 SU(2) Wigner j-symbols

We use the definition of the Wigner’s (3jm) symbol given in book:varshalovic with the orthogonality property

$$\sum_{m_{1},m_{2}}\left(\begin{array}{ccc}j_{1}&j_{2}&j_{3}\\ m_{1}&m_{2}&m_{3}\end{array}\right)\left(\begin{array}{ccc}j_{1}&j_{2}&j_{3}\\ m_{1}&m_{2}&n_{3}\end{array}\right)=\frac{\delta_{j_{3}l_{3}}\delta_{m_{3}n_{3}}}{2j_{3}+1}.$$

(37)

These are normalized and vanish if triangular inequalities are not satisfied. We define (4jm) symbol as the contraction of two (3jm) symbol using an intertwiner k:

i^{(k)}_{m_{1}m_{2}m_{3}m_{4}}\equiv\left(\begin{array}{cccc}j_{1} j_{2} j_{3} j_{4}\qquad m_{1} m_{2} m_{3} m_{4}\end{array}\right)^{(k)}\qquad {} \equiv\sum_{m_{i}}(-1)^{k-m_{i}}\left(\begin{array}{ccc}j_{1} j_{2} k\qquad m_{1} m_{2} m_{i}\end{array}\right)\left(\begin{array}{ccc}k j_{3} j_{4} {}-m_{i} m_{3} m_{4}\end{array}\right) ,

(40) (45)

with the orthogonality relations:

$$\sum_{m_{1},m_{2},m_{3}}\left(\begin{array}{cccc}j_{1} j_{2} j_{3} j_{4}\qquad m_{1} m_{2} m_{3} m_{4}\end{array}\right)^{(k_{1})}\left(\begin{array}{cccc}j_{1} j_{2} j_{3} j_{4}\qquad m_{1} m_{2} m_{3} m_{4}\end{array}\right)^{(k_{2})}\qquad {}=\frac{\delta_{k_{1}k_{2}}}{2k_{1}+1}\frac{\delta_{j_{4}l_{4}}\delta_{m_{4}n_{4}}}{2j_{4}+1} ,$$

(50) (51)

where (δ_k₁ k₂)/(2k₁+1) is the normalization factor. A synthetic notation for (4jm) symbols is also:

$$\left(\begin{array}{c}j_{f}\\ m_{f}\end{array}\right)^{(k)}\equiv\left(\begin{array}{cccc}j_{1}&j_{2}&j_{3}&j_{4}\\ m_{1}&m_{2}&m_{3}&m_{4}\end{array}\right)^{(k)}.$$

| | (52) |

With these definitions and the conventional definitions of Wigner’s matrices for SU(2), D⁽ʲ⁾ₘₙ(g), we have:

$${} \int dg D^{j_{1}}_{m_{1}n_{1}}\left(g\right)D^{j_{2}}_{m_{2}n_{2}}\left(g\right)D^{j_{3}}_{m_{3}n_{3}}\left(g\right)D^{j_{4}}_{m_{4}n_{4}}\left(g\right)\qquad {}=\sum_{k}d_{k}\left(\begin{array}{cccc}j_{1} j_{2} j_{3} j_{4}\qquad m_{1} m_{2} m_{3} m_{4}\end{array}\right)^{(k)}\left(\begin{array}{cccc}j_{1} j_{2} j_{3} j_{4}\qquad n_{1} n_{2} n_{3} n_{4}\end{array}\right)^{(k)}.$$

(53) (58)

We denote the 4-valent intertwiner with the ket |k⟩, defined as:

$$|k\rangle\equiv i^{(k)}_{m_{1}m_{2}m_{3}m_{4}}|j_{1},m_{1}\rangle|j_{2},m_{2}\rangle|j_{3},m_{3}\rangle|j_{4},m_{4}\rangle$$

(59)

where the (4jm) symbol is defined in (40). The states (59) constitute a basis in the intertwiner space of a 4-valent node:

$$\mathcal{I}_{4}\equiv Inv\left[V^{j_{1}}\otimes V^{j_{2}}\otimes V^{j_{3}}\otimes V^{j_{4}}\right]$$

(60)

where Vʲᵢ is the irreducible representation of spin jᵢ. The resolution of the identity in the intertwiner space I₄ and the orthogonality condition can be written as:

$$\mathbf{1}^{(\mathcal{I}_{4})} =\sum_{k=0}^{2j}(2k+1)|k\rangle\langle k|,\qquad \langle k|k^{\prime}\rangle =\frac{\delta_{k,k^{\prime}}}{(2k+1)}.$$

(61) (62)

The 15-j symbol we use in this work is the irreducible symbol of first type, following the convention of GraphMethods . Its graphical representation and its definition in terms of Wigner’s 6-j symbols is

\raisebox{-42.67912pt}{ \includegraphics[width=113.81102pt]{Single15j}} = \left\lbrace \begin{array}{ccccc}j_{1} j_{2} j_{3} j_{4} j_{5}\qquad l_{1} l_{2} l_{3} l_{4} l_{5}\qquad k_{1} k_{2} k_{3} k_{4} k_{5}\end{array}\right\rbrace =\qquad {}= (-1)^{\sum_{i=1}^{5}j_{i}+l_{i}+k_{i}}\sum_{s}d_{s}\left\lbrace \begin{array}{ccc}j_{1} k_{1} s k_{2} j_{2} l_{1}\end{array}\right\rbrace \left\lbrace \begin{array}{ccc}j_{2} k_{2} s\qquad k_{3} j_{3} l_{2}\end{array}\right\rbrace \left\lbrace \begin{array}{ccc}j_{3} k_{3} s\qquad k_{4} j_{4} l_{3}\end{array}\right\rbrace \left\lbrace \begin{array}{ccc}j_{4} k_{4} s\qquad k_{5} j_{5} l_{4}\end{array}\right\rbrace \left\lbrace \begin{array}{ccc}j_{5} k_{5} s j_{1} k_{1} l_{5}\end{array}\right\rbrace .

(64) (75)

C.2 Booster functions

The booster functions Dona2018 , article:Dona_etal_2019_numerical_study_lorentzian_EPRL_spinfoam_amplitude , also known as B4 functions Speziale2016 , are the non compact remnants of the SL(2,ℂ) integrals. These functions turn out to encode all the details of the EPRL model, such as the Y-map. We define them as:

B_{4}\left(j_{f},l_{f};i,k\right)\equiv\frac{1}{4\pi}\sum_{p_{f}}\left(\begin{array}{c}j_{f}\\ p_{f}\end{array}\right)^{(i)}\left(\int_{0}^{\infty}\mathrm{d}r\sinh^{2}r\,\prod_{f=1}^{4}d^{(\gamma j_{f},j_{f})}_{j_{f}l_{f}p_{f}}(r)\right)\left(\begin{array}{c}l_{f}\\ p_{f}\end{array}\right)^{(k)} =\raisebox{-35.56593pt}{ \includegraphics[width=68.28644pt]{SingleBooster}},

(76)

where d^(ρ,k) are the boost matrix elements for γ-simple irreducible representations of SL(2,ℂ) in the principal series, γ is the Barbero-Immirzi parameter and the (4jm) symbols are defined in section C.1.

The booster function B₄(j_f,l_f;i,k) is interpreted as a quantum tetrahedron being boosted among adjacent frames: the two sets j_f and l_f describe the four areas of the tetrahedron in the two frames connected by a boost, and the two intertwiners i and k describe the quantum intrinsic shape of the tetrahedron Speziale2016 . For a precise interpretation of the booster functions and their semiclassical limit we refer to dona2020 . The explicit form of the boost matrix elements can be found in the literature. For the general form, see Ruhl:1970fk ; Speziale2016 . We just report here the case of simple irreducible representations:

$$d^{(\gamma j,j)}_{jlp}(r)=\qquad \times e^{-(j-i\gamma j+p+1)r}\sum_{s}\frac{(-1)^{s}\,e^{-2sr}}{s!(l-j-s)!}\,{}_{2}F_{1}[l+1-i\gamma j,j+p+1+s,j+l+2,1-e^{-2r}] .$$

(77)

References

• (1)

  A. Perez, “The Spin-Foam Approach to Quantum Gravity,” Living Reviews in Relativity 16 (2013) .

• (2)

  R. Oeckl, “A ’general boundary’ formulation for quantum mechanics and quantum gravity,” Phys. Lett. B575 (2003) 318–324, arXiv:0306025.

• (3)

  C. Rovelli, D. Colosi, L. Doplicher, W. Fairbairn, L. Modesto, and K. Noui, “Background independence in a nutshell,” Class. Quant. Grav. 22 (2005), 2971-2990. arXiv:0408079.

• (4)

  W. J. Fairbairn and C. Meusburger, “q-Deformation of Lorentzian spin foam models,” PoS QGQGS2011 (2011) 17. arXiv:1112.2511.

• (5)

  M. Han, “Cosmological Constant in LQG Vertex Amplitude,” Phys. Rev. D 84 (2011) 64010. arXiv:1105.2212

• (6)

  H. M. Haggard, M. Han, W. Kamiński, and A. Riello, “Four-dimensional Quantum Gravity with a Cosmological Constant from Three-dimensional Holomorphic Blocks,” Phys. Lett. B752 (2016) 258–262. arXiv:1509.00458.

• (7)

  R. Pereira, “Lorentzian LQG vertex amplitude,” Class. Quant. Grav. 25 (2008) 85013. arXiv:0710.5043

• (8)

  E. R. Livine and S. Speziale, “A new spinfoam vertex for quantum gravity,” Phys. Rev. D76 (2007) 84028. arXiv:0705.0674

• (9)

  E. R. Livine and S. Speziale, “Consistently Solving the Simplicity Constraints for Spinfoam Quantum Gravity,” Europhys. Lett. 81 (2008) 50004. arXiv:0708.1915.

• (10)

  L. Freidel and K. Krasnov, “A New Spin Foam Model for 4d Gravity,” Class. Quant. Grav. 25 (2008) 125018, arXiv:0708.1595.

• (11)

  J. Engle, R. Pereira, and C. Rovelli, “Flipped spinfoam vertex and loop gravity,” Nucl. Phys. B798 (2008) 251-290. arXiv:0708.1236

• (12)

  J. Engle, R. Pereira, and C. Rovelli, “The loop-quantum-gravity vertex-amplitude,” Phys. Rev. Lett. 99 (2007) 161301. arXiv:0705.2388

• (13)

  J. Engle, E. Livine, R. Pereira, and C. Rovelli, “LQG vertex with finite Immirzi parameter,” Nucl. Phys. B799 (2008) 136-149. arXiv:0711.0146

• (14)

  W. Kaminski, M. Kisielowski, and J. Lewandowski, “Spin-Foams for All Loop Quantum Gravity,” Class. Quant. Grav. 27 (2010) 95006. arXiv:0909.0939

• (15)

  Y. Ding, M. Han, and C. Rovelli, “Generalized Spinfoams,” Phys. Rev. D83 (2011) 124020. arXiv:1011.2149.

• (16)

  C. Rovelli and F. Vidotto, Covariant Loop Quantum Gravity. Cambridge University Press, 2014.

• (17)

  C. Rovelli, “Zakopane lectures on loop gravity,” PoS QGQGS2011 (2011) 003. arXiv:1102.3660

• (18)

  P. Donà and S. Speziale, “Introductory lectures to loop quantum gravity,” in TVC 79. Gravitation : théorie et expérience, A. Bounames and A. Makhlouf, eds. Hermann, 2013.

• (19)

  B. Dittrich, “The continuum limit of loop quantum gravity - a framework for solving the theory,” in 100 Years of General Relativity: Volume 4, Loop Quantum Gravity: The First 30 Years, A. Bounames and A. Makhlouf, eds. World Scientific, 2017. arXiv:1409.1450.

• (20)

  B. Dittrich, S. Mizera, and S. Steinhaus, “Decorated tensor network renormalization for lattice gauge theories and spin foam models,” New Journal of Physics 18 no. 5, (May, 2016) 053009. arXiv:1409.2407

• (21)

  B. Bahr and S. Steinhaus, “Investigation of the spinfoam path integral with quantum cuboid intertwiners,” Phys. Rev. D93 no. 10, (May, 2016) . arXiv:1508.07961

• (22)

  V. Bonzom, R. Gurau, and V. Rivasseau, “Random tensor models in the largeNlimit: Uncoloring the colored tensor models,” Phys. Rev. D85 no. 8, (Apr, 2012) . arXiv:1202.3637

• (23)

  D. Benedetti and R. Gurau, “Phase transition in dually weighted colored tensor models,” Nucl. Phys. B855 no. 2, (Feb, 2012) 40-437. arXiv:1108.5389

• (24)

  S. Carrozza, D. Oriti, and V. Rivasseau, “Renormalization of a SU(2) Tensorial Group Field Theory in Three Dimensions,” Communications in Mathematical Physics 330 no. 2, (Mar, 2014) 581-637. arXiv:1303.6772

• (25)

  J. B. Geloun, T. A. Koslowski, D. Oriti, and A. D. Pereira, “Functional renormalization group analysis of rank-3 tensorial group field theory: The full quartic invariant truncation,” Phys. Rev. D97 no. 12, (Jun, 2018). arXiv:1805.01619

• (26)

  V. Bonzom and M. Smerlak, “Bubble divergences from twisted cohomology,” Comm. Math. Phys.312 (2), 399-426 (2012). arXiv:1008.1476

• (27)

  V. Bonzom and M. Smerlak, “Bubble divergences: sorting out topology from cell structure,” Annales Henri Poincaré 13 (2012) 185-208. arXiv:1103.3961

• (28)

  A. Baratin, S. Carrozza, D. Oriti, J. Ryan, and M. Smerlak, “Melonic Phase Transition in Group Field Theory,” Letters in Mathematical Physics 104 no. 8, (May, 2014) 1003-1017. arXiv:1307.5026

• (29)

  C. Perini, C. Rovelli, and S. Speziale, “Self-energy and vertex radiative corrections in LQG,” Phys. Lett. B682 (2009) 78–84. arXiv:0810.1714

• (30)

  T. Krajewski, J. Magnen, V. Rivasseau, A. Tanasa, and P. Vitale, “Quantum Corrections in the Group Field Theory Formulation of the EPRL/FK Models,” Phys. Rev. D82 (2010) 124069. arXiv:1007.3150

• (31)

  A. Riello, “Self-Energy of the Lorentzian EPRL-FK Spin Foam Model of Quantum Gravity,” Phys. Rev. D88 (2013) 24011. arXiv:1302.1781

• (32)

  P. Donà, “Infrared divergences in the EPRL-FK Spin Foam model,” Class. Quant. Grav. 35 no. 17, (2018) 175019. arXiv:1803.00835

• (33)

  J. W. Barrett, R. J. Dowdall, W. J. Fairbairn, F. Hellmann, and R. Pereira, “Lorentzian spin foam amplitudes: graphical calculus and asymptotics,” Classical and Quantum Gravity 27 (2010). arXiv:0907.2440

• (34)

  F. Gozzini, “A high-performance code for EPRL spin foam amplitudes,” Classical and Quantum Gravity 38 no. 22, (Oct, 2021) 225010. arXiv:2107.13952

• (35)

  J. Bezanson, A. Edelman, S. Karpinski, and V. B. Shah, “Julia: A fresh approach to numerical computing,” SIAM review 59 no. 1, (2017) 65–98.

• (36)

  W. R. Inc., “Mathematica, Version 12.2.”. Champaign, IL, 2020.

• (37)

  E. Bianchi and Y. Ding, “Lorentzian spinfoam propagator,” Phys.Rev. D86 (2012) 104040. arXiv:1109.6538

• (38)

  E. Bianchi and H. M. Haggard, “Discreteness of the volume of space from Bohr-Sommerfeld quantization,” Phys. Rev. Lett. 107 (2011) 011301. arXiv:1102.5439

• (39)

  L. Freidel and S. Speziale, “Twisted geometries: A geometric parametrisation of SU(2) phase space,” Phys.Rev. D82 (2010) 84040. arXiv:1001.2748

• (40)

  H. T. Johansson and C. Forssén, “Fast and Accurate Evaluation of Wigner 3j, 6j, and 9j Symbols Using Prime Factorisation and Multi-Word Integer Arithmetic,” SIAM Journal on Scientific Computing 38 no. 1, (Apr., 2015) A376–A384.

• (41)

  J. Rasch and A. C. H. Yu, “Efficient Storage Scheme for Precalculated Wigner 3j, 6j and Gaunt Coefficients,” SIAM Journal on Scientific Computing 25 no. 4, (2004) 1416–1428, arXiv:1807.03066

• (42)

  S. Speziale, “Boosting Wigner’s nj-symbols,” Journal of Mathematical Physics 58 no. 3, (Sep, 2017) . arXiv:1609.01632

• (43)

  P. Donà and G. Sarno, “Numerical methods for EPRL spin foam transition amplitudes and Lorentzian recoupling theory,” General Relativity and Gravitation 50 no. 10, (2018). arXiv:1807.03066

• (44)

  P. Donà, M. Fanizza, G. Sarno, and S. Speziale, “Numerical study of the Lorentzian Engle-Pereira-Rovelli-Livine spin foam amplitude,” Physical Review D 100 (2019). arXiv:1903.12624.

• (45)

  github.com/PietropaoloFrisoni/Self-energy.

• (46)

  S. K. Asante, B. Dittrich, and H. M. Haggard, “Effective Spin Foam Models for Four-Dimensional Quantum Gravity,” Phys. Rev. Lett. 125 no. 23, (Dec., 2020) 231301. arXiv:2004.07013

• (47)

  F. V. Pietropaolo Frisoni, Francesco Gozzini, “Markov Chain Monte-Carlo methods for graph refinement in covariant Loop Quantum Gravity,” in preparation (2022) .

• (48)

  J. W. Barrett, R. J. Dowdall, W. J. Fairbairn, H. Gomes, and F. Hellmann, “Asymptotic analysis of the EPRL four-simplex amplitude,” J. Math. Phys. 50 (2009) 112504. arXiv:0902.1170

• (49)

  P. Donà, M. Fanizza, G. Sarno, and S. Speziale, “SU(2) graph invariants, Regge actions and polytopes,” Classical and Quantum Gravity 35 (2018). arXiv:1708.01727.

• (50)

  V. D. Aleksandroviic, A. N. Moskalev, and K. V. Kel’manoviic, Quantum theory of angular momentum: irreducible tensors, spherical harmonics, vector coupling coefficients, 3nj symbols. World scientific, 1988.

• (51)

  P. Donà, M. Fanizza, G. Sarno, and S. Speziale, “Numerical study of the Lorentzian Engle-Pereira-Rovelli-Livine spin foam amplitude,” Physical Review D 100 no. 10, (Mar, 2019). arXiv:1903.12624

• (52)

  F. Gozzini and F. Vidotto, “Primordial Fluctuations From Quantum Gravity,” Frontiers in Astronomy and Space Sciences 7 (Feb, 2021) . arXiv:1906.02211

• (53)

  M. Han, Z. Huang, H. Liu, D. Qu, and Y. Wan, “Spinfoam on a Lefschetz thimble: Markov chain Monte Carlo computation of a Lorentzian spinfoam propagator,” Physical Review D 103 no. 8, (Apr, 2021) . arXiv:2012.11515

• (54)

  M. Han, Z. Huang, H. Liu, and D. Qu, “Complex critical points and curved geometries in four-dimensional Lorentzian spinfoam quantum gravity,” 2021. arXiv:2110.10670.

• (55)

  P. Frisoni, “Studying the EPRL spinfoam self-energy,” Proceedings of The Sixteenth Marcel Grossmann Meeting, 2021. arXiv:2112.08528

• (56)

  A. P. Yutsin, I. B. Levinson, and V. V. Vanagas, Mathematical Apparatus of the Theory of Angular Momentum. Israel Program for Scientific Translation, Jerusalem, Israel, 1962.

• (57)

  P. Donà and S. Speziale, “Asymptotics of lowest unitary SL(2,C) invariants on graphs,” 2020. arXiv:2007.09089

• (58)

  W. Ruhl, The Lorentz group and harmonic analysis. W.A. Benjamin, Inc, New York, 1970.