Spin Foam Amplitudes Research Articles

Geometry transition in spinfoams

Abstract We show how the fixed-spin asymptotics of the EPRL model can be used to perform the spin-sum for spin foam amplitudes defined on fixed two-complexes without interior faces and contracted with coherent spin-network states peaked on a discrete simplicial geometry with macroscopic areas. We work in the representation given in Ref 1. We first rederive the latter in a different way suitable for our purposes. We then extend this representation to 2-complexes with a boundary and derive its relation to the coherent state representation. We give the measure providing the resolution of the identity for Thiemann's state in the twisted geometry parametrization. We then piece together these with other results in the literature and show how the spin sum can be performed analytically using the model asymptotics. These results are relevant to analytic investigations regarding the transition of a black hole to a white hole geometry. In particular, this work was the basis of the calculation presented in Ref 2.

A new 2+1 coherent spin-foam vertex for quantum gravity

Abstract We propose an explicit spin-foam amplitude for Lorentzian gravity in three dimensions, allowing for both space- and time-like boundaries. The model is based on two main requirements: that it should be structurally similar to its well-known Euclidean analog, and that geometricity should be recovered in the semiclassical regime. To this end we introduce new coherent states for space-like boundary edges, derived from the continuous series of unitary SU ( 1 , 1 ) representations. We show that the relevant objects in the amplitude can be written in terms of the defining representation of the group, just as so happens in the Euclidean case. We derive an expression for the semiclassical amplitude at large spins, showing it relates to the Lorentzian Regge action.

Summing bulk quantum numbers with Monte Carlo in spin foam theories

We introduce a strategy to compute EPRL spin foam amplitudes with many internal faces numerically. We work with sl2cfoam-next, the state-of-the-art framework to numerically evaluate spin foam transition amplitudes. We find that uniform sampling Monte Carlo is exceptionally effective in approximating the sum over internal quantum numbers of a spin foam amplitude, considerably reducing the computational resources necessary. We apply it to compute large volume divergences of the theory and find surprising numerical evidence that the EPRL vertex renormalization amplitude is instead finite.

Markov chain Monte Carlo methods for graph refinement in spinfoam cosmology

We study the behavior of the Lorentzian Engle-Pereira-Rovelli-Livine spinfoam amplitude with homogeneous boundary data, under a graph refinement going from five to twenty boundary tetrahedra. This can be interpreted as a wave function of the Universe, for which we compute boundary geometrical operators, correlation functions, and entanglement entropy. The numerical calculation is made possible by adapting the Metropolis-Hastings algorithm, along with recently developed computational methods appropriate for the deep quantum regime. We confirm that the transition amplitudes are stable against such refinement. We find that the average boundary geometry does not change, but the new degrees of freedom correct the quantum fluctuations of the boundary and the correlations between spatial patches. The expectation values are compatible with their geometrical interpretation and the correlations between neighboring patches decay when computed across different spinfoam vertices.

Experimental simulation of loop quantum gravity on a photonic chip

The unification of general relativity and quantum theory is one of the fascinating problems of modern physics. One leading solution is Loop Quantum Gravity (LQG). Simulating LQG may be important for providing predictions which can then be tested experimentally. However, such complex quantum simulations cannot run efficiently on classical computers, and quantum computers or simulators are needed. Here, we experimentally demonstrate quantum simulations of spinfoam amplitudes of LQG on an integrated photonics quantum processor. We simulate a basic transition of LQG and show that the derived spinfoam vertex amplitude falls within 4% error with respect to the theoretical prediction, despite experimental imperfections. We also discuss how to generalize the simulation for more complex transitions, in realistic experimental conditions, which will eventually lead to a quantum advantage demonstration as well as expand the toolbox to investigate LQG.

Geometry from local flatness in Lorentzian spin foam theories

Local flatness is a property shared by all the spin foam models. It ensures that the theory's fundamental building blocks are flat by requiring locally trivial parallel transport. In the context of simplicial Lorentzian spin foam theory, we show that local flatness is the main responsible for the emergence of geometry independently of the details of the spin foam model. We discuss the asymptotic analysis of the EPRL spin foam amplitudes in the large quantum number regime, highlighting the interplay with local flatness.

Complex critical points and curved geometries in four-dimensional Lorentzian spinfoam quantum gravity

This paper focuses on the semiclassical behavior of the spinfoam quantum gravity in four dimensions. There has been long-standing confusion, known as the flatness problem, about whether the curved geometry exists in the semiclassical regime of the spinfoam amplitude. The confusion is resolved by the present work. By numerical computations, we explicitly find curved Regge geometries that contribute dominantly to the large-$j$ Lorentzian Engle-Pereira-Rovelli-Livine (EPRL) spinfoam amplitudes on triangulations. These curved geometries are with small deficit angles and relate to the complex critical points of the amplitude. The dominant contribution from the curved geometry to the spinfoam amplitude is proportional to ${e}^{i\mathcal{I}}$, where $\mathcal{I}$ is the Regge action of the geometry plus corrections of higher order in curvature. As a result, in the semiclassical regime, the spinfoam amplitude reduces to an integral over Regge geometries weighted by ${e}^{i\mathcal{I}}$, where $\mathcal{I}$ is the Regge action plus corrections of higher order in curvature. As a by-product, our result also provides a mechanism to relax the cosine problem in the spinfoam model. Our results provide important evidence supporting the semiclassical consistency of the spinfoam quantum gravity.

Numerical analysis of the self-energy in covariant loop quantum gravity

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.

Finiteness of spinfoam vertex amplitude with timelike polyhedra and the regularization of full amplitude

This work focuses on Conrady-Hnybida's 4-dimensional extended spinfoam model with timelike polyhedra, while we restrict all faces to be spacelike. Firstly, we prove the absolute convergence of the vertex amplitude with timelike polyhedra, when SU(1,1) boundary states are coherent states or the canonical basis, or their finite linear combinations. Secondly, based on the finite vertex amplitude and a proper prescription of the SU(1,1) intertwiner space, we construct the extended spinfoam amplitude on arbitrary cellular complex, taking into account the sum over SU(1,1) intertwiners of internal timelike polyhedra. We observe that the sum over SU(1,1) intertwiners is infinite for the internal timelike polyhedron that has at least 2 future-pointing and 2 past-pointing face-normals. In order to regularize the possible divergence from summing over SU(1,1) intertwiners, we develop a quantum cut-off scheme based on the eigenvalue of the ``shadow operator''. The spinfoam amplitude with timelike internal polyhedra (and spacelike faces) is finite, when 2 types of cut-offs are imposed: one is imposed on $j$ the eigenvalue of area operator, the other is imposed on the eigenvalue of shadow operator for every internal timelike polyhedron that has at least 2 future-pointing and 2 past-pointing face-normals.

Analytic continuation of spinfoam models

The Lorentzian Engle-Pereira-Rovelli-Livine/Freidel-Krasnov (EPRL/FK) spinfoam model and the Conrady-Hnybida (CH) timelike-surface extension can be expressed in the integral form $\int e^S$. This work studies the analytic continuation of the spinfoam action $S$ to the complexification of the integration domain. Our work extends our knowledge from the real critical points well-studied in the spinfoam large-$j$ asymptotics to general complex critical points of $S$ analytic continued to the complexified domain. The complex critical points satisfying critical equations of the analytic continued $S$. In the large-$j$ regime, the complex critical points give subdominant contributions to the spinfoam amplitude when the real critical points are present. But the contributions from the complex critical points can become dominant when the real critical point are absent. Moreover the contributions from the complex critical points cannot be neglected when the spins $j$ are not large. In this paper, we classify the complex critical points of the spinfoam amplitude, and find a subclass of complex critical points that can be interpreted as 4-dimensional simplicial geometries. In particular, we identify the complex critical points corresponding to the Riemannian simplicial geometries although we start with the Lorentzian spinfoam model. The contribution from these complex critical points of Riemannian geometry to the spinfoam amplitude give $e^{-S_{Regge}}$ in analogy with the Euclidean path integral, where $S_{Regge}$ is the Riemannian Regge action on simplicial complex.

Four-dimensional spinfoam quantum gravity with a cosmological constant: Finiteness and semiclassical limit

We present an improved formulation of 4-dimensional Lorentzian spinfoam quantum gravity with cosmological constant. The construction of spinfoam amplitudes uses the state-integral model of PSL(2,$\mathbb{C}$) Chern-Simons theory and the implementation of simplicity constraint. The formulation has 2 key features: (1) spinfoam amplitudes are all finite, and (2) With suitable boundary data, the semiclassical asymptotics of the vertex amplitude has two oscillatory terms, with phase plus or minus the 4-dimensional Lorentzian Regge action with cosmological constant for the constant curvature 4-simplex.

A high-performance code for EPRL spin foam amplitudes

We present sl2cfoam-next, a high-performance software library for computing Lorentzian EPRL spin foam amplitudes. The library improves on previous codes by many orders of magnitude in single-core performance, can be parallelized on a large number of CPUs and on the GPU, and can be used interactively. We describe the techniques used in the code and provide many usage examples. As first applications, we use sl2cfoam-next to complete the numerical test of the Lorentzian single-vertex asymptotics and to confirm the presence of the ‘flatness problem’ of spin foam models in the BF and EPRL cases.

Efficient Simulation of Loop Quantum Gravity: A Scalable Linear-Optical Approach

The problem of simulating complex quantum processes on classical computers gave rise to the field of quantum simulations. Quantum simulators solve problems, such as boson sampling, where classical counterparts fail. In another field of physics, the unification of general relativity and quantum theory is one of the greatest challenges of our time. One leading approach is loop quantum gravity (LQG). Here, we connect these two fields and design a linear-optical simulator such that the evolution of the optical quantum gates simulates the spin-foam amplitudes of LQG. It has been shown that computing transition amplitudes in simple quantum field theories falls into the bounded-error quantum polynomial time class, which strongly suggests that computing transition amplitudes of LQG are classically intractable. Therefore, these amplitudes are efficiently computable with universal quantum computers, which are, alas, possibly decades away. We propose here an alternative special-purpose linear-optical quantum computer that can be implemented using current technologies. This machine is capable of efficiently computing these quantities. This work opens a new way to relate quantum gravity to quantum information and will expand our understanding of the theory.

Numerical computations of next-to-leading order corrections in spinfoam large- j asymptotics

We numerically study the next-to-leading order corrections of the Lorentzian Engle-Pereira-Rovelli-Livine (EPRL) 4-simplex amplitude in the large-$j$ expansions. We perform large-$j$ expansions of Lorentzian EPRL 4-simplex amplitudes with two different types of boundary states, the coherent intertwiners and the coherent spin-network, and numerically compute the leading-order and the next-to-leading order $O(1/j)$ contributions of these amplitudes. We also study the dependences of these $O(1/j)$ corrections on the Barbero-Immirzi parameter $\gamma$. We show that they, as functions of $\gamma$, stabilize to finite real constants as $\gamma\to\infty$. Lastly, we obtain the quantum corrections to the Regge action because of the $O(1/j)$ contribution to the spinfoam amplitude.

Numerical analysis of spin foam dynamics and the flatness problem

In this paper, we apply a recently proposed numerical algorithm for finding stationary phase points in spin foam amplitudes. We study a spin foam amplitude with three vertices and a bulk face in 4d BF theory. We fix the boundary coherent states to three possible triangulations, one with zero deficit angle on the bulk face and two with non-zero deficit angle. We compute the amplitude numerically and we find a stationary phase point already at low spins in all the three cases. We comment on how this result contrasts with the claims of flatness problem in spin foam theories. We point out where these arguments may be misleading and we propose further computations to definitively answer the question.

Searching for classical geometries in spin foam amplitudes: a numerical method

We develop a numerical method to investigate the semiclassical limit of spin foam amplitudes with many vertices. We test it using the Ponzano–Regge model, a spin foam model for three-dimensional euclidean gravity, and a transition amplitude with three vertices. We study the summation over bulk spins, and we identify the stationary phase points that dominate it and that correspond to classical geometries. We complement with the numerical analysis of a four vertex transition amplitude and with a modification of the model that includes local curvature. We discuss the generalization of our results to the four-dimensional EPRL spin foam model, and we provide suggestions for new computations.

Numerical study of the Lorentzian Engle-Pereira-Rovelli-Livine spin foam amplitude

The Lorentzian EPRL spin foam amplitude for loop quantum gravity is a multi-dimensional non-compact integral of highly oscillating functions. Using a method based on the decomposition of Clebsch-Gordan coefficients for the unitary infinite-dimensional representations of SL(2,C) in terms of those of SU(2), we are able to provide for the first time numerical evaluations of the vertex amplitude. The values obtained support the asymptotic formula obtained by Barrett and collaborators with a saddle point approximation, showing, in particular, a power-law decay and oscillations related to the Regge action. The comparison offers a test of the efficiency of the method. Truncating the decomposition to the first few terms provides a qualitative matching of the power-law decay and oscillations. For vector and Euclidean Regge boundary data, a qualitative matching is obtained with just the first term, which corresponds to the simplified EPRL model. We comment on future developments for the numerics and extension to higher vertices. We complete our work with some analytic results: We provide an algorithm and explicit configurations for the different geometries that can arise as boundary data, and explain the geometric consequences of the decomposition used.

Semiclassical behavior of spinfoam amplitude with small spins and entanglement entropy

Spinfoam amplitudes with small spins can have interesting semiclassical behavior and relate to semiclassical gravity and geometry in four dimensions. We study the generalized spinfoam model [spinfoams for all loop quantum gravity (LQG) [Kaminski et al., Spin-foams for all loop quantum gravity, Classical Quantum Gravity 27, 095006 (2010); Erratum, Classical Quantum Gravity 29, 049502(E) (2012), Ding et al., Generalized spinfoams, Phys. Rev. D 83, 124020 (2011)] with small spins $j$ but a large number of spin degrees of freedom (d.o.f.), and find that it relates to the simplicial Engle-Pereira-Rovelli-Livine-Freidel-Krasnov model with large spins and Regge calculus by coarse-graining spin d.o.f. $\mathrm{Small}\text{\ensuremath{-}}j$ generalized spinfoam amplitudes can be employed to define semiclassical states in the LQG kinematical Hilbert space. Each of these semiclassical states is determined by a four-dimensional Regge geometry. We compute the entanglement R\'enyi entropies of these semiclassical states. The entanglement entropy interestingly coarse grains spin d.o.f. in the generalized spinfoam model, and satisfies an analog of the thermodynamical first law. This result possibly relates to the quantum black hole thermodynamics in [Ghosh and Perez, Black Hole Entropy and Isolated Horizons Thermodynamics, Phys. Rev. Lett. 107, 241301 (2011); 108, 169901(E) (2012)].

Asymptotic analysis of spin foam amplitude with timelike triangles

The large $j$ asymptotic behavior of $4$-dimensional spin foam amplitude is investigated for the extended spin foam model (Conrady-Hnybida extension) on a simplicial complex. We study the most general situation in which timelike tetrahedra with timelike triangles are taken into account. The large $j$ asymptotic behavior is determined by critical configurations of the amplitude. We identify the critical configurations that correspond to the Lorentzian simplicial geometries with timelike tetrahedra and triangles. Their contributions to the amplitude are phases asymptotically, whose exponents equal to Regge action of gravity. The amplitude may also contains critical configurations corresponding to non-degenerate split signature $4$-simplices and degenerate vector geometries. But for vertex amplitudes containing at least one timelike tetrahedron and one spacelike tetrahedron, critical configurations only give Lorentzian $4$-simplices, while the split signature and degenerate $4$-simplices do not appear.

Infrared divergences in the EPRL-FK spin foam model

We provide an algorithm to estimate the divergence degree of the Lorentzian EPRL-FK spin foam amplitudes for arbitrary 2-complexes. We focus on the ‘self-energy’ and ‘vertex renormalization’ diagrams and find an upper bound estimate. We argue that our upper bound must be close to the actual value, and explain what numerical improvements are needed to verify this numerically. For the self-energy, this turns out to be significantly more divergent than the lower bound estimate present in the literature. We support the validity of our algorithm using 3-stranded versions of the amplitudes (corresponding to a toy 3d model) for which our estimates are confirmed numerically. We also apply our methods to the simplified EPRLs model, finding an utterly convergent behavior, and to BF theory, independently recovering the divergent estimates present in the literature.