We investigate the compatibility of minimally coupled scalar, fermion and gauge fields with asymptotically safe quantum gravity, using nonperturbative functional renormalization group methods. We study $d=4,5$ and 6 dimensions and within certain approximations find that for a given number of gauge fields there is a maximal number of scalar and fermion degrees of freedom compatible with an interacting fixed point at positive Newton coupling. The bounds impose severe constraints on grand unification with fundamental Higgs scalars. Supersymmetry and universal extra dimensions are also generally disfavored. The standard model and its extensions accommodating right-handed neutrinos, the axion and dark-matter models with a single scalar are compatible with a fixed point.

Intertwiners are the building blocks of spin-network states. The space of intertwiners is the quantization of a classical symplectic manifold introduced by Kapovich and Millson. Here we show that a theorem by Minkowski allows us to interpret generic configurations in this space as bounded convex polyhedra in R3: A polyhedron is uniquely described by the areas and normals to its faces. We provide a reconstruction of the geometry of the polyhedron: We give formulas for the edge lengths, the volume, and the adjacency of its faces. At the quantum level, this correspondence allows us to identify an intertwiner with the state of a quantum polyhedron, thus generalizing the notion of the quantum tetrahedron familiar in the loop quantum gravity literature. Moreover, coherent intertwiners result to be peaked on the classical geometry of polyhedra. We discuss the relevance of this result for loop quantum gravity. In particular, coherent spin-network states with nodes of arbitrary valence represent a collection of semiclassical polyhedra. Furthermore, we introduce an operator that measures the volume of a quantum polyhedron and examine its relation with the standard volume operator of loop quantum gravity. We also comment on the semiclassical limit of spin foams with nonsimplicial graphs.

In a quantum system in a pure state, a subsystem generally has a nonzero entropy because of entanglement with the rest of the system. Is the average entanglement entropy of pure states also the typical entropy of the subsystem? We present a method to compute the exact formula of the momenta of the probability $P(S_A)\mathrm{d}{S}_A$ that a subsystem has entanglement entropy $S_A$. The method applies to subsystems defined by a subalgebra of observables with a center. In the case of a trivial center, we reobtain the well-known result for the average entropy and the formula for the variance. In the presence of a nontrivial center, the Hilbert space does not have a tensor product structure and the well-known formula does not apply. We present the exact formula for the average entanglement entropy and its variance in the presence of a center. We show that for large systems the variance is small, $\mathrm{\Delta }{S}_A/⟨S_A⟩\ll 1$, and therefore the average entanglement entropy is typical. We compare exact and numerical results for the probability distribution and comment on the relation to previous results on concentration of measure bounds. We discuss the application to physical systems where a center arises. In particular, for a system of noninteracting spins in a magnetic field and for a free quantum field, we show how the thermal entropy arises as the typical entanglement entropy of energy eigenstates.

We investigate some aspects of the renormalization group flow of gravity in the presence of fermions, which have remained somewhat puzzling so far. The first is the sign of the fermionic contribution to the running of Newton's constant, which depends on details of the cutoff. We argue that only one of the previously used schemes correctly implements the cutoff on eigenvalues of the Dirac operator, and it acts in the sense of screening Newton's constant. We also show that K\"ahler fermions give the same contribution to the running of the cosmological and Newton constant as four Dirac spinors. We then calculate the graviton contributions to the beta functions by imposing the cutoffs on the irreducible spin components of the tetrad. In this way we can probe the gauge dependence of the off-shell flow. The results resemble closely those of the metric formalism, except for an increased scheme\char22{}and (off shell) gauge dependence.

We explicitly compute the tree-level on-shell four-graviton amplitudes in four, five and six dimensions for local and weakly nonlocal gravitational theories that are quadratic in both, the Ricci and scalar curvature with form factors of the d’Alembertian operator inserted between. More specifically we are interested in renormalizable, super-renormalizable or finite theories. The scattering amplitudes for these theories turn out to be the same as the ones of Einstein gravity regardless of the explicit form of the form factors. As a special case the four-graviton scattering amplitudes in Weyl conformal gravity are identically zero. Using a field redefinition, we prove that the outcome is correct for any number of external gravitons (on-shell n−point functions) and in any dimension for a large class of theories. However, when an operator quadratic in the Riemann tensor is added in any dimension (with the exception of the Gauss-Bonnet term in four dimensions) the result is completely altered, and the scattering amplitudes depend on all the form factors introduced in the action.

The Lorentzian Engle-Pereira-Rovelli-Livine (EPRL) spin foam amplitude for loop quantum gravity is a multidimensional noncompact 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.