Finite entanglement entropy and spectral dimension in quantum gravity

What are the conditions on a field theoretic model leading to a finite entanglement entropy density? We prove two very general results: 1) Ultraviolet finiteness of a theory does not guarantee finiteness of the entropy density; 2) If the spectral dimension of the spatial boundary across which the entropy is calculated is non-negative at all scales, then the entanglement entropy cannot be finite. These conclusions, which we verify in several examples, negatively affect all quantum-gravity models, since their spectral dimension is always positive. Possible ways out are considered, including abandoning the definition of the entanglement entropy in terms of the boundary return probability or admitting an analytic continuation (not a regularization) of the usual definition. In the second case, one can get a finite entanglement entropy density in multi-fractional theories and causal dynamical triangulations.


Introduction: entanglement entropy and return probability
Quantum correlations in the vacuum state of a quantum field give rise to entanglement entropy once we trace over the degrees of freedom associated, for example, with a region of space. If Σ is the boundary of such region (a codimension-2 hypersurface), customary derivations using path-integral techniques in Euclidean time and the replica trick [1][2][3][4] lead to an entanglement entropy which is proportional to the "area" of the boundary times a divergent integral where ǫ is an ultraviolet (UV) cut-off for small times or coarse resolutions which regularizes the integral and σ is the square of a fictitious diffusion-time parameter, or the inverse of the resolution ∝ 1/ √ σ in a measurement. The function P d (σ) in the integrand is the return probability for a diffusion process on the (d = D − 2)dimensional boundary hypersurface, where D is the topological dimension of spacetime. Equivalently, P d (σ) is given by the trace of the heat kernel 1 associated with the propagator of a field theory living on the boundary [6]. In particular, (2) represents the effective action of the field restricted to the boundary which, for a free theory, is proportional to the vacuum energy.
Notice that the integral can be interpreted as the regularized version of the entanglement entropy density The short-distance correlations responsible for the entanglement entropy are thus intimately related to the properties of a diffusion process on the boundary. Is there a way to formalize this relation? In [4], Solodukhin showed that there is indeed a relation between the short-distance scaling of the two-point function and the UV divergence of the entanglement entropy, but their power-law behaviour coincides only for a standard quadratic energy-momentum dispersion relation. The relevance of entanglement entropy across a boundary in quantum gravity is well exemplified in Jacobson's derivation of the Einstein equations as an equation of state [7]. The key assumption of his result is that one can associate a finite (universal) entropy density to any local Rindler horizon and that such entropy is due to the correlations present in the vacuum state of a quantum field. This scenario has led to suggest (see, e.g., [8]) that quantum entanglement across the Rindler horizon can be obtained from correlations of the quantum geometry itself and, due to the fundamental discreteness of the geometry, the associated entropy density is finite.
In view of this, it is interesting to investigate whether there exist unconventional quantum field theories which incorporate putative quantum-gravity effects and naturally predict a finite entropy density. Field theories with deformed dispersion relations which either break Lorentz invariance or preserve it (being the deformed dispersion relation a non-linear function of the ordinary Lorentz-invariant mass-shell relation) have already been ruled out in [3] as possible candidates. In that work, the authors showed that such theories can make the short-distance behaviour of the twopoint function finite but, no matter how regular this behaviour is, the return probability associated with a deformed dispersion relation always diverges in the limit σ → 0, hence the need of a UV cut-off ǫ in (2).
In this paper, we turn our attention to more general classes of field-theoretical models and explore whether they lead to a finite entanglement entropy. 2 As a starting point, in Sect. 2 we discuss a theory with compact momentum space, linked to a well-studied model of non-commutative field theory, with all the good UV properties one might desire (a well-behaved Euclidean propagator, calculated here for the first time, and a finite return probability at σ = 0) which, however, fails to have a finite entanglement entropy density. This counter-example is important because it is a top-down theory that does not fulfill the expectation of a finite entanglement entropy density, despite realizing precisely the features of the ad hoc toy model analyzed in [3] and bypassing the no-go result for modified dispersion relations presented there. In Sects. 3 and 4, we present the main result of this paper, a theorem (proven almost in full) stating a background-independent necessary condition under which a finite entanglement entropy density ρ d can be obtained: the spectral dimension d b S of the spatial boundary of spacetime must be negative. We split the proof in two parts, a technically involved 2 Quantum entanglement is a hot topic on which new papers appear almost on a daily basis. After the submission of this manuscript, we became aware of work on the entanglement entropy with cosmological applications: see [9] and references therein for more details.
one on the special value d b S = 0 (Sect. 3) and a much easier one on the case d b S > 0 with strict equality (Sect. 4).
Retrospectively, from the results of Sect. 3 it becomes easy to understand the failure of the compactmomentum-space example of Sect. 2 as due to a violation of this condition. Interestingly, many kinematical states of the class of theories based on a discrete combinatorial structure (group field theory (GFT) and its realizations as loop quantum gravity and spin foams) have a vanishing spectral dimension at the scale of the combinatorial structure [10] and, according to our conjecture, their entanglement entropy density is either infinite or zero.
An interesting byproduct of our analysis regards the treatment of the spectral dimension d S of spacetime. While d S is traditionally described, in the quantumgravity literature, in terms of a diffusion process on the spacetime geometry, here we employ the related but essentially alternative view of d S as the real pole(s) in the spectral zeta function of the theory. Technically, this entails a transformation of the heat kernel from the usual position or momentum space to the Laplace-Mellin momentum space. This tool is very well known in the spectral analysis of fractal geometries [5,[14][15][16] but it has not received much attention by quantum-gravity researchers. We will show its usefulness in several applications. This is not the end of the story. An ever more powerful result establishes the impossibility to get a finite entanglement entropy density if the spectral dimension is positive definite (Sect. 4). This excludes essentially all quantum-gravity and string-related models, where d b S ≥ 0, but it also clashes with some positive results obtained with independent techniques, which we will discuss later. To understand such a contradiction, we should look for ways out of this no-go theorem. The most obvious one is to abandon the definition (3) and concentrate on other formulations of the entanglement entropy (such as those employed, precisely, in string theory). However, here we propose a more interesting alternative. Without recurring to the regularization (2), one can calculate the integral (3) via an analytic continuation in the parameter space of the theory. Once integration is performed, one can revert back to the allowed parameter range and find a finite, positive ρ d . We will show this for the multi-fractional theory with q-derivatives (T q in short) [11] as well as for others theories such as causal dynamical triangulations (CDT), and comment on the viability of this solution.
Our top-down multi-scale examples (noncommutative and multi-fractional) show that, contrary to common wisdom, there is no clear relation between (a) good UV properties and (b) a finite entropy density. The compact-momentum-space case realizes (a) but not (b), while T q realizes (b) but not (a), at least not in an obvious way. The necessary condition found here does not have the full status of a theorem because its proof entails the properties of the spectral zeta function, which can be studied analytically when the spectral dimension is constant but that become slightly more complicated in a multi-scale geometry (varying spectral dimension, also known as dimensional flow). We will employ an extension to multi-scale geometries that seems very robust. None of this depends on a specific background or dynamics: we use only intrinsic properties of dimensional flow, in the same spirit of recent developments in quantum gravities and multi-fractional theories [11,12].
Before starting, it may be useful to comment on the application of our results to other quantum-gravity models. First, elements such as back-reaction and the renormalization of Newton's constant usually play a role in full calculations of the entanglement entropy at the semi-classical level, while, as we just noted, our approach is more general in the sense that it solely relies on the spectral dimension. The spectral dimension is always well defined in the semi-classical limit, but also in other regimes with heavy non-perturbative quantum effects. This is shown by many examples in the literature (most being fully non-perturbative), among which we mention causal dynamical triangulations, asymptotic safety, loop quantum gravity and spin foams, and string field theory. Therefore, any method employing the spectral dimension may be valid even beyond the semi-classical level, depending on the specific theory.
Having said that, it is important to specify in which way the method captures quantum-gravity features. We start from the entanglement entropy density (3) defined for a generic quantum field theory and we examine under what conditions this quantity is rendered finite. To model quantum-gravity effects, these field theories, representing generic degrees of freedom in two neighboring regions separated by a surface or horizon, have deformed dispersion relations and live on a spacetime with fractal-like features. This abstract framework does not aim to represent full blown quantum gravity with all its details and caveats. However, it does encode dimensional flow, the change of spacetime dimension with the probed scale induced by quantum-gravity effects in all known theories. Dimensional flow is a direct, global manifestation of physics beyond Einstein's gravity and, in particular, is tightly related to microscopic quantum uncertainties in spacetime texture [13]. Although, as is obvious, dimensional flow does not capture all features of quantum gravity, it is enough to extract interest-ing phenomenology. In the present case, it is sufficient to draw a consequence regarding the entanglement entropy. In fact, the abstract quantum field theory on which one applies (3) can be regarded as an effective description of some degrees of freedom (matter or gravitational) living in a spacetime with varying spectral dimension. The underlying fundamental theory, whatever it is, induces effective dispersion relations or kinetic terms in the effective field theory and, most of all, governs the way the spectral dimension changes. Since no quantum gravity to date predicts a negative spectral dimension, the result of the theorem follows.
By this, we do not mean that the entanglement entropy of the full theory is certainly ill defined. We will consider ways out of this restrictive theorem via possible discreteness or quantum mechanisms in specific theories (Sect. 5). Our model-independent findings are enough to conclude that dimensional flow per se is not sufficient to get a finite results, while we leave the role of discreteness as an open question to be explored in future work.
2 Infinite entropy density in a UV-finite example: compact momentum space Before starting our discussion of UV-modified field theories in Euclideanized spacetime, let us recall the relationship between the return probability and the Euclidean two-point function G(x, y) of a scalar theory in ordinary four-dimensional Euclidean space. It is well known that G(x, y) can be written in terms of the heatkernel K(x, y; σ) as with where C(p) is the Fourier transform of the wave operator of the theory which, for a massive scalar field, takes the familiar form C(p) = p 2 + m 2 . The return probability is proportional to trace of the heat kernel: and thus it is related with the short-distance behaviour of the two-point function, Such limit usually diverges as the inverse power of the Euclidean distance between x and y. Let us look, for example, at the two-point function of a massive scalar field: The integral over momentum space can be performed by transforming to spherical coordinates, in particular putting y = 0. The two-point function can be expressed as where I 1 (ipx) is a modified Bessel function of the first kind, and the final result reads which, for a massless field, simplifies to Equations (10) and (11) are divergent in the UV limit x → 0, as expected. Intuitively, one might think that a theory with a regular two-point function in the UV, reflecting a finite behaviour of quantum correlations at short distances, should also exhibit a finite entanglement entropy density. Probably, the easiest way to modify a field theory in order to get a finite UV behaviour is to introduce "by hand" a constant length scale to the distance dependence in the denominator of the twopoint function [17]. However, it can be shown that such modification amounts to a choice of deformed energymomentum dispersion relation [3] falling into the class of modified field theories with a diverging return probability and entanglement entropy density. As a further example of models with an intrinsic UV scale, we hereby consider field theories with a non-trivial geometry of momentum space. This structure naturally introduces a fundamental mass scale in the two-point function without spoiling the Lorentz invariance of the theory. A famous example of such models is a field theory defined on a de Sitter momentum space [18][19][20][21][22][23], the momentum-space counterpart of a field theory on κ-Minkowski non-commutative spacetime, associated with κ-deformations of relativistic symmetries [23][24][25][26][27][28][29]. For our purposes, we need to consider the Euclidean version of this model, which we assume to be characterized by a compact momentum space given by a 4-sphere of radius κ 3 . The radius κ provides a UV energy scale, the Euclidean analogue of the de Sitter momentum-space "cosmological constant". 4 In order to determine whether the geometry of momentum space reflects in a finite UV behaviour of the Euclidean two-point function, we explicitly evaluate the Green function from its integral representation. Focusing for simplicity on the massless case, the two-point function reads where the factor 1 − p 2 /κ 2 is obtained from the integration measure over the 4-sphere. Upon evaluation, the integral gives which remains finite in the UV limit: while for large distances x → ∞ the correlations vanish as expected. This shows that the compact geometry of Euclidean momentum space naturally renders the UV correlations of the theory finite. One might expect that this finite behaviour in the UV is related to a finite return probability when σ vanishes, which would suggest the possibility of a finite entanglement entropy. In general, in theories with compact momentum space or a non-trivial momentumspace measure the technical reason why this can happen is rather simple. For such models, the return probability is a momentum integral with a certain measure w(p) (non-trivial both in theories with compact momentum space and, e.g., in the multi-scale spacetimes described by multi-fractional theories [11]) and a modified dispersion relation C(p) = p 2 K(p) (where K(p) = 1 in the ordinary case): where the support of w(p) may limit the integration range to a semi-interval or an interval. After the change of variables k 2 := σp 2 , one gets Now, in all models related to quantum gravity, C(∞) = ∞ and C(0) = 0. When C(p) = p 2 , the integrand is σ-independent apart from the measure weight, while in the general case the exponential goes to zero when σ → 0. On the other hand, in all known instances either w(0) = 1 or, as in compact-momentum-space theories, it has a compact support where the upper limit ∞ is replaced by √ σκ, which collapses to zero when σ → 0. Therefore, in the limit σ → 0 the integral is either σ-independent or tends to zero. In the first case, the return probability diverges as P D (σ) ∼ σ −D/2 , while in the second case one has a 0/0 expression that could give a finite result.
Unfortunately, having a finite P D (0) is not sufficient to have a finite entanglement entropy. In this section, we will show how a known three-dimensional example [32] and a four-dimensional model with compact momentum space do have a finite P D (0) but cannot generate a finite entanglement entropy. These counter-examples bar the possibility to cure the divergence in the entanglement entropy by a general mechanism involving, among other possibilities, compact momentum spaces.
In [32], the diffusion process for an SU (2) momentum space, the (Euclideanized) momentum space of a particle coupled to three-dimensional gravity, was studied to determine the associated spectral dimensional flow. In particular, it was found that, in the limit σ → 0, the associated spectral dimension goes to zero. These results were based on the return probability determined by the heat kernel in SU (2) momentum space, where κ is a momentum UV cut-off related to the radius of the SU (2) group manifold. The integral can be carried out analytically and the result is given in terms of modified Bessel functions of the first kind I n (x): It is immediate to see that, unlike the return probability in ordinary flat momentum space, P κ 3 (σ) does not diverge in the limit σ → 0 and, indeed, we have The normalization in (17) has been chosen so that the heat kernel is normalized to 1 when integrated over all position space and, consistently, a finite κ corresponds to a finite P κ 3 (0). This should be compared with the The same holds for the return probability for a diffusion process on a one-dimensional spacetime boundary, where the momentum space is assumed to be a circular one-dimensional sub-manifold of SU (2): for which This seems very promising, but the real question is whether this finite behaviour of the return probability suffices to ensure a finite entropy density. We can immediately see that this is not the case. Indeed, the integral is logarithmically divergent and a UV regulator ǫ is needed in order to obtain a finite result. Analogously, in the example of a momentum space given by a 4-sphere discussed above, we can consider a two-dimensional boundary diffusion process with momentum space given by a 2-sphere: where F is the Dawson F -function and we left out a σ-independent overall factor. This return probability is finite in the short-diffusion-time limit, since However, as in the SU (2) momentum-space example above, the entropy density obtained from the Mellin transform of the return probability, is divergent and a UV regulator ǫ is required to render it finite. As in the example of the circle, the divergence goes as ln ǫ.
To summarize, contrary to naive expectations, having a finite return probability in the limit of diffusion scale going to zero does not guarantee a finite result for the entanglement entropy density. It may be possible that this first approach to the calculation of the entropy density in a deformed field theory is too simple-minded because it relies on two assumptions: (a) the calculation of the entanglement entropy leading to the relation (3) between entropy and return probability carries through in a field theory on a curved momentum space in the same way as in an ordinary field theory, and (b) the momentum space for the boundary diffusion process is a compact section of the larger manifold momentum space under consideration. However, both assumptions are quite natural. In the first case, the path integral for a free field theory on the non-trivial momentum spaces considered poses no technical problem [23] and should go through as in the ordinary case. The main step that, in principle, should require some care would be to make sure that the boundary conditions for the field can be imposed unambiguously. Regarding the second point, we do not really see any particular concern or a plausible alternative approach.
It thus seems that demanding a finite entropy density is quite restrictive and it would be very interesting to explore more in depth what kind of non-standard features a field-theoretical model should exhibit in order to comply with such requirement. For this purpose, we now turn to a different perspective.

Entanglement entropy and spectral dimension
Since the finiteness of the entanglement entropy is related to the analytic properties of the return probability, and since the latter also determines the value of the spectral dimension d S of the space or spacetime on which the fictitious diffusion process is considered, the study of d S may offer a new insight in the problem of the entanglement entropy.

Spectral dimension from the zeta function
Let us recall that the spectral dimension of spacetime is defined as the scaling of the return probability, A constant d S corresponds to a power law P(σ) ∼ σ −dS/2 ; in ordinary Euclidean(ized) space(time), d S coincides with the topological dimension D. We will indicate as d b S the spectral dimension of the (D − 2)dimensional spatial boundary determining the entan-glement entropy. Integrating (26) and plugging the resulting expression (with D → d) into (3), one gets The intuitive and simple case where d S is constant illustrates the role of the sign of the spectral dimension. If A whole literature has been spent on the dependence of the spectral dimension on the resolution 1/ √ σ as an indicator of the geometric properties of exotic spacetimes. Usually in quantum gravity, Eq. (26) is the only expression used for the purpose. However, here we take a different and much less explored view, which is more familiar to the fields of fractal geometry and statistical mechanics [5]. Given a function f (σ) living in Ê, let us consider its associated zeta function Applying this formula to P D (σ), we have the zeta function of the return probability: This is just a special way to write the Mellin transform of P D . Using the identity and the definition of the return probability, it is easy to convince oneself that another way to cast the zeta function of the heat kernel is where λ n are the eigenvalues of the Laplacian. In the continuum case, with the return probability given by (15), and for a non-compact spacetime, λ n = p 2 K(p) and where Ω D is the angular integral. Now, the zeta function can be analytically continued to the complex plane, where it has poles s m ∈ . In scale-invariant geometries, the largest real part in this set is nothing but half the spectral dimension d S in a small-σ expansion of the heat kernel. 5 To understand 5 When d S is constant, it can be defined by d S = −2 lim σ→0 ln P D (σ)/ ln σ, which explains the small-σ expansion mentioned in the text. When d S is not constant, its correct definition is (26).
where this result comes from, we recall that the inverse Mellin transform of the zeta function (28) is Evaluating it by the residue theorem, we immediately recognize the power-law behavior of the return probability P D (σ) ∼ σ −dS/2 at small σ.
We give two examples of calculation of constant d S via (28) instead of using (26). In the standard case where I is a divergent integral that can be made finite by introducing a regulator p uv . The only pole s = D/2 is real, so that d S = D. Another instance is w = 1 and K = |p| 2(γ−1) , corresponding to a massless dispersion relation |p| 2γ = 0. This type of dispersion relation, where γ is some real parameter, is typical of the UV regime of many quantum gravities and phenomenological models (the full dispersion relation would be of the form p 2 +β|p| 2γ = 0, where β is a constant). Clearly, the standard case corresponds to γ = 1. The zeta function of the return probability is where p uv is a UV cut-off introduced to regularize the integral when D − 2sγ > 0. The pole now is at s = D/(2γ) = d S /2, i.e., correctly reproducing the spectral dimension of spacetimes with this modified dispersion relation [33,34]. To obtain Eq. (32), we used Eq. (30) and integrated in momenta p after performing the integration in σ of (28). Switching integration order (first on p and then on σ), one finds P D (σ) ∝ σ −D/(2γ) and ζ D (s) ∝ δ[D/(2γ)−s]/Γ (s), thus getting again Eq. (33). We will come back to this point later.
In no-scale geometries (sets or spaces with constant dimensionality), the poles of the spectral zeta function are known to describe the underlying geometric properties of the set on which ζ D is calculated. In particular, the pole with largest real part Re(s) = d S /2 determines the spectral dimension d S in a small-σ expansion [5,[14][15][16]. In the most general case, ζ D (s) has many complex poles and the imaginary part of those such that Re(s n ) = d S /2 defines a complex dimension d CS . This situation is typical of fractals [5,16,35] and also of multi-fractional spacetimes [11,36]. Here we are going to see that the description of the dimension of spacetime in terms of spectral geometry is valid only at plateaux in dimensional flow, i.e., only in those scale ranges where d S is approximately constant. We argue that in multi-scale geometries (variable d S ), one has to consider all the poles of the zeta function, not just the one with largest real value. The examples provided here (the compact-momentum-space model of the previous section and the multi-fractional theory of the next section) strongly indicate a principle which, in the most general terms, can be stated as follows. Meaning by "poles" the representative n, m, . . . In other words, whenever the spectral dimension of spacetime varies continuously along a resolution (probing) scale σ, flat regions in dimensional flow appear (d S (σ) almost flat and constant in some regions in ln σ). This claim is a conjecture only because we do not have a formal mathematical proof at hand. However, its validity is strongly supported by the observation that the multi-scale-zeta principle is just an application, to all plateau regions in dimensional flow, of the well-known rigorous case where d S is exactly constant. There are no compelling reasons why the spectral zeta function should not encode all the main values of d S in a multiscale geometry. The multi-scale-zeta conjecture agrees with the general notion that the plateaux in dimensional flow carry the main information on the underlying geometry and that transient regions between different plateaux are strongly model-dependent and associated with non-physical features such as regularization schemes [10,34].
It is not difficult to recognize an immediate consequence of the multi-scale-zeta principle: Zero-dimension lemma. If d S = 0 at some point or plateau in dimensional flow, then ζ D (0) is finite.
This lemma deals with the only case not explicitly contemplated by the multi-scale-zeta principle. It states that a zero spectral dimension does not correspond to a pole of ζ D (s). To prove that ζ D (0) remains finite, one would have to show that, in the limit s → 0, the zeta function always has a divergence ∼ Γ (s) in the numerator, cancelling exactly the factor Γ (s) in the denominator in (28). Spectral theory supports this result. A wellknown characteristic of the zeta function is that its pole structure is usually given by a combination of the Riemann zeta (in turn related to Euler gamma function) and of gamma functions. Therefore, it is no wonder that a pole in the numerator would scale as Γ (s). The zerodimension lemma is nothing but a consequence of this cancellation due to the definition (28).
A multi-scale example satisfying the zero-dimension lemma is the case of the SU (2) momentum space. Taking the Mellin transform of Eq. (18), we get which, again, is finite at s = 0. In this case, the poles are s n = 3/2 + n, where n ∈ AE, and the smallest 6 pole gives the spectral dimension of space d S = 3; this is correct only in the IR. The full profile d S (σ) runs from 0 in the UV up to some peak > 3 at a mesoscopic scale, then dropping down to 3 in the IR [32].
Another check of the principle and of the lemma will be given in Sect. 4.

Boundary zeta function and entanglement entropy
Let us go back to the entropy density problem. On the spatial boundary of spacetime, all the above formulae hold with D replaced by d = D − 2. The spectral zeta function on the boundary reads Up to a divergent normalization factor 1/Γ (0), we see that ζ d (0) coincides with the entanglement entropy density (3): The constants λ n in (29) are the eigenvalues of the Laplacian restricted to the spatial boundary and the real part of the poles of ζ d is the spectral dimension d b S at different plateaux, where "b" stands for boundary. For non-compact spaces, the boundary can be regarded as an abstract (D − 2)-surface. It becomes clear now why the examples with compact momentum space did not give rise to a finite ρ d . The one-dimensional momentum space S 1 has a zeta function Crucially, the factor 1/Γ (s) in the definition of the zeta function compensates a factor Γ (s) coming from integration of the return probability. This is the reason why 6 Here no small-σ expansion is made.
Eq. (37) is finite at s = 0 and Eq. (22), which lacks the 1/Γ (0) prefactor, is not. Incidentally, d b S = 1 for this boundary space, but one realizes that this is true only in the IR. The full profile d b S (σ) obtained from the restriction of (26) to the boundary of spacetime runs monotonically from 0 in the UV to 1 in the IR.
Turning our attention to the two-dimensional boundary and transforming Eq. (23), one has up to an s-independent finite normalization constant. This expression diverges at the family of poles s n = 1 + n, and n = 0 (the representative of the family) corresponds to d b S = 2. Again, ζ 2 (0) is finite while Eq. (25) is ill defined.
From this exercise, we learn two things: (a) the entanglement entropy density ρ d of our compact-space examples always diverges as Γ (s) in the limit of Laplace-Mellin momentum s → 0; (b) the zeta function ζ d of the heat kernel on the boundary does not have any pole at s = 0 and thus misses the UV behavior of the spatial geometry thereon (d b S = 0 in the UV in these examples). The crucial factor here is the presence of a vanishing spectral dimension, i.e., the existence of a pole in s = 0 of the zeta function. The combination Γ (s)ζ d (s) in (36) diverges both at the poles of ζ d (s) and at s = 0, provided the latter is not a zero of the zeta function. However if ζ d (0) ∝ 1/Γ (0) = 0, then Eq. (36) is finite. Conversely, if ζ d (0) is finite or diverges, then the entropy density diverges.
What is the physical meaning of this? When applied to D → d and Eq. (36), the results obtained so far (multi-scale-zeta principle and zero-dimension lemma) yield a novel and important consequence, a Necessary condition for a finite entropy density. If the entanglement entropy density (36) is finite, then the spectral dimension d b S of the spatial boundary never vanishes at any scale. This is the physical result we were looking for. It is impossible to have a finite entanglement entropy density when the spectral dimension of the space boundary is zero at some scale.
At this point, we can understand the failure of compact momentum space models in getting a finite entanglement entropy: it is because their spectral dimension vanishes in the UV. Technically, this happens because the zeta function does not scale as ∼ 1/Γ (s) for small s. The fact that the meromorphic function ζ d is analytic at s = 0 leads to a divergent ρ d .
To summarize how we reached the formulation of the above necessary condition, we traded the usual definition (3) of the entanglement entropy density for the combination of (35) and (36). This is completely equivalent to (3) or (27), the only difference being that in (3) or (27) the limit s → 0 is taken before integrating in σ, while in Eq. (36) it is taken afterwards. This operation commutes in all the cases considered here (compare Eqs. (22), (25) and (47) with those that can be found in the limit σ → 0 from, respectively, (37), (38) and (48)) and it does not entail any apparent problem such as regularization artifacts; in particular, we do not use any regularization such as that in (2). If ρ d is finite and non-zero, it will be exactly the same when computed from (3), (27), or (36). The main result is then obtained from the fact that the spectral dimension is, by definition, the (real part of the) pole of the boundary zeta function ζ d (s). In particular, if the boundary spectral dimension vanishes, then ζ d (0) diverges. But if ζ d (0) diverges, then the right-hand side of Eq. (36) cannot be finite.
Let us stress that the non-vanishing d b S condition is necessary but not sufficient. Even when the spectral dimension does not vanish, it may be possible that the entanglement entropy density is ill defined for other reasons. The case with constant spectral dimension (33) is a good example of this problem: in (32) (integration first in σ and then in p) we introduced a UV cut-off p uv to render the zeta function finite and, indeed, the entanglement entropy density is infinite: ρ d = lim s→0 Γ (s)/(d − 2sγ) = +∞. Integrating first in p and then in σ, for γ finite we get zero instead of infinity, ρ d = lim s→0 δ[D/(2γ) − s] = 0.

No-go result and ways out
After finding, through general arguments and via explicit examples, that UV finiteness is not sufficient and that a non-vanishing spectral dimension is necessary, we discuss a third and more serious obstacle towards a finite entropy density, the following necessary condition: the boundary spectral dimension must be negative definite at short scales. We first prove it and then consider its consequences.
From the definition (26), In essentially all quantum gravity and string related models, d S ≥ 0 at all scales. Since σ ≥ 0 and the return probability is positive definite, then the right-hand side is negative semi-definite, which implies that P D (σ) increases as σ decreases. On the other hand, the integral (3) is finite at σ = 0 provided P d (σ) vanishes faster than σ at short scales. This contradicts the previous condition when applied to the boundary, unless d b S (σ) < 0 at σ ∼ 0. This configuration makes (3) finite also at σ = +∞ provided P d (σ) does not grow faster than σ at large scales. We do not need to review all the calculations of the spectral dimension in quantum gravities to convince the reader that the condition d b S (σ) < 0 is unphysical and never realized. However, we cannot accept this result uncritically because we do have at least two examples of theories with finite entanglement entropy density.
-A class of theories of quantum gravity is based on non-local operators of exponential form exp [39][40][41][42][43][44]. Kinetic terms in field actions have the typical appearance The same type of operators appear also in the gravitational sector and render it non-local. Due to the strong suppression ∼ e −p 2 /M 2 /p 2 of the free propagator, the spectral dimension of non-local quantum gravity flows to 0 in the UV [45], in flagrant violation of the first necessary condition found here. Therefore, we predict a divergent entropy density in models with this type of free propagator. At the interacting level, the situation is not so clear. The dressed propagator, which is a modification of the free one by interactions, gives rise to a different dispersion relation, hence to a different effective return probability, spectral dimension and entanglement entropy. In [46], it is argued that in UV-finite interacting theories characterized by non-local terms of the type (40), the absence of divergent contributions to the renormalized Newton's constant indeed makes the entropy finite. The key elements here are the interaction terms and UV-finiteness rather than non-locality. -The low-energy limit of string field theory has exactly the same type of non-locality (40) but the spectral dimension in the UV is bounded from below by the worldsheet [45]. Therefore, at short scales and the spectral dimension never vanishes. From this, according to our spectral analysis the entanglement entropy in string theory may be finite. This statement, as imprecise as it is, agrees with explicit calculations of the entropy regarded as an effect of quantum entanglement [47,48] and with the statistical-mechanics derivation of the blackhole entropy-area law from the counting of microscopic degrees of freedom (BPZ states) [49]. Only explicit calculations in certain regimes of the theory, such as those where the AdS/CFT correspondence applies, could verify the above conclusion. Thus, either interactions or different definitions of the entanglement entropy may avoid the no-go theorem. However, instead of giving up the definition (3) we propose a third alternative: to continue it analytically whenever possible. This solution is less radical, but perhaps subtler, than inserting an ad hoc UV regularization as in expression (2). For both these reasons, it deserves our attention.

Finite entropy density with analytic continuation
In this section, we illustrate the validity of the multiscale-zeta principle, of the zero-dimension lemma and of the analytically continued finite-ρ d necessary condition with a fully analytic example: the multi-fractional theory T q with q-derivatives. This framework, reviewed in [11], is a field theory defined on a multi-fractal geometry with varying Hausdorff and spectral dimension. Dimensional flow follows a multi-parametric universal profile determined only by very general (backgroundand dynamics-independent) properties of the spacetime dimension [12]. Spacetime is characterized by a hierarchy of length scales ℓ * ≡ ℓ 1 > ℓ 2 > . . .; above ℓ * , ordinary geometry is recovered, while at scales ℓ * the spacetime dimensionality drops below D. At ultra-microscopic scales, a discrete structure naturally emerges, encoded in logarithmic oscillations of the geometry. Physical observables are affected by this scale dependence. The generality of the measure, together with the fact that all quantum gravities have dimensional flow, justifies the interest in these theories as an efficient framework wherein to explore all the main physical consequences of dimensional flow and to constrain them with experiments and observations ranging from particle physics to cosmology. There is also the hope to improve the perturbative quantization properties of gravity thanks to dimensional flow itself, as it may happen in other theories.
Depending on the symmetries of the field Lagrangian, there are different versions of multi-fractional dynamics. Here we concentrate on one of the best studied cases, the theory T q . Here the return probability can be calculated exactly and it reads [37] where the profile τ (σ) depends on the scale hierarchy in the same way as the spacetime measure. In general, one can take the polynomial form τ (σ) = N n=1 g n (σ/σ n ) αn F n (σ), where g n are constants, 0 < α n ≤ 1 are called fractional exponents and F n is a modulation factor affecting all scales and related to the discrete structure in the UV. Without any loss of information, one can consider a binomial profile with N = 2 and only one fundamental scale σ * = σ 1 ∝ ℓ 2 * and one modulation factor F ω dependent on a frequency parameter ω. Overall, where 0 < α < 1. In this paper, we coarse grain over log oscillations, which is equivalent to set F ω = 1. Calculating (26) with (42) and (43), one gets [37] which is positive definite because α ≥ 0. The calculation of the entanglement entropy density can be done via three completely equivalent routes. We develop all the details of the first to highlight a delicate analytic continuation one should be careful about; 7 the other two routes pass through the same steps. Plugging Eqs. (42) and (43) into (3), where the absolute value accounts for a change of sign of the integral depending on whether α ≷ 1. This integral coincides with the one in formula 3.194.4 of [63], where β = 1, ν = d/2, µ = −dα/[2(1 − α)] and B is given in formula 8.384.1. For the integral (46) to converge to the right-hand side, it must be |argβ| < π and Reν > Reµ > 0. The first condition is trivial, while the second requires 1 > −α/(1 − α) > 0, i.e., α < 0.
Assuming that, one finds Before commenting this expression, let us derive it also from the spectral dimension and from the zeta function.
Considering the boundary dimension d b S and plugging (44) (with the spacetime dimension D replaced by the boundary dimension d; obviously, this is not a regularization) into (27), one obtains (47). Setting F ω = 1 allows us to find ζ d analytically. 8 Alternatively, combining Eqs. (28), (42) and (43), and continuing analytically the integral to the regions of the parameter space of the theory, we get From (36), one gets again (47). In none of these routes leading to (47) did we use any regularization, but we did make an analytic continuation to the region α < 0.
Negative values of α correspond to a negative spectral dimension, which is consistent with the no-go theorem. At this point, however, we make an experiment: we analytically continue (47) and (48) to values 0 < α < 1 and see whether we find a sensible result. There are several consequences we can draw from (48).
-Let α = 0. The poles of ζ D are s n = D/2 + n(1 − α) and s m = Dα/2 − m(1 − α), where n, m ∈ AE. For n = 0 = m, we obtain the IR and UV spectral dimension, respectively: in agreement with (44). Thus, we have shown that the spectral zeta function brings about information as regards the spectral dimension of both the IR and the UV plateau, when neither of them is zero. Generalizing to many fractional coefficients α n , one will find as many poles as the number of non-zero plateaux d S ≃ Dα n . Although we were unable to find an analytic expression of ζ D beyond the binomial case, there is really no conceptual obstacle against this expectation. The multi-scale-zeta principle is thus confirmed. -ζ D (s) diverges in the limit α → 1, in agreement with the regularized divergence in (32) with γ = 1. -When α = 0, the UV spectral dimension obtained from Eqs. (42) and (43) is d UV S = 0, in agreement with Eq. (49). However, this result can no longer be obtained from the zeta function (48), which reduces to Here we can appreciate the cancellation between the Γ (s) factor in the denominator and the factor Γ [(s− Dα/2)/(1 − α)] → Γ (s) in the numerator. At s = 0, ζ D,0 (0) = (4πσ * ) −D/2 is finite, in agreement with the zero-dimension lemma. Consistently, in the limit σ * → 0 we get the power-law divergence of (32) with γ = 1. -For D → d (zeta function of the boundary diffusion process), the entanglement entropy density is exactly (47). For a generic α, this expression is finite, consistently with the fact that the residue of (31) at s = 0 is zero. The d-dimensional entropy density diverges at α = 1 + d/(2n) and at α where n, m ∈ AE + . The first set of poles is excluded because α < 1. The second set corresponds to geometries with infinite entropy density. Therefore, in D = 4 the exponents α = 1/2, 2/3, 3/4, . . . are ruled out if we insist to have a finite ρ 2 . This leads to two major results.
-The value α = 1/2, suggested by quantumgravity arguments and realizing the famous twodimensional UV limit d UV S = 2 so popular in quantum gravity, does not give a finite entanglement entropy. In a sense to be better specified below, this may indicate that it is difficult to get a finite ρ d in quantum gravity. The compactmomentum-space model and some results mentioned in the introduction [3] support this view.
-The excluded values are a countable subset of 1/2 ≤ α < 1. This range is special in another multi-fractional theory called T γ=α , characterized by multi-scale fractional derivatives in the action and of which T q is an approximation [11]. The spectral dimension in T γ=α is basically the same as in T q and we may transpose all the present discussion to that case. In T γ=α , values in the interval 0 < α < 1/2 are excluded if we want to have a normed space at all scales. The limit case α = 1/2 is acceptable but it corresponds to a taxicab geometry where the shortest path between two points is not unique [38]. On the other hand, when log oscillations are turned on and the full discrete structure of spacetime is considered, the norm condition breaks down in the UV for any α and one may simply give it up altogether. In fact, space is normed in the IR regardless of what happens in the UV. In turn, log oscillations generate a rather peculiar behaviour of the UV propagator that can have important consequences for the renormalization of gravity not explorable with naive power-counting arguments [11]. The bottom line of all this is that, at least in multi-fractional theories but probably 0.05 0.50 5 50 Fig. 1 The typical kink profile of the boundary spectral dimension in many quantum gravities (solid curve). This is reproduced by the multi-fractional theory with q-derivatives, whose zeta function (48)  also in quantum gravity at large, a finite entanglement entropy density (range 0 < α < 1/2) may be tightly related to geometries highly nontrivial in the UV and where gravity is well behaved.
Applied to the spatial boundary (D → d), Eq. (49) and the general profile (typical of many theories of quantum gravity) derived from Eq. (48) is depicted in Fig. 1.
Let us now comment on the analytic continuation 0 < α < 1 → α < 0 → 0 < α < 1. Figure 2 shows the behaviour of the entanglement entropy density (47) as a function of α. As expected, ρ d is positive definite in the healthy but unphysical case α < 0. However, ρ d is finite and positive also in the intervals Although this might be just an artifact of the analytic continuation, it is intriguing to speculate that such analytic continuation might be simply part of the definition of the entanglement entropy density when calculated via the boundary return probability. If this were true, it would be a genuine intrinsic property of the theory, not a UV regularization as (2). To check the absence of inconsistencies, one should calculate ρ d from other definitions, for instance using microstates or thermodynamics, and verify whether the entanglement entropy is really finite in the ranges (52) or if, on the contrary, it diverges for the allowed parameter space. As this would go beyond the scope of this work, we will stop here, but not before noting that this solution would rescue not only the multi-fractional theory with q-derivatives, but also other quantum gravity models with a dimensional flow similar to (44). In fact, the profile of the spectral dimension of multi-fractional theories with binomial measure, where d S (σ) runs between two non-zero values, is typical also in many other models of quantum gravity. Quantum gravities such as CDT [50][51][52][53] and asymptotic safety [54,55] respect our weak necessary condition (their scale-dependent spectral dimension never vanishes), but not the no-go theorem because d S > 0 in these cases. However, for the case of CDT we can say something more. The CDT profile of the spectral dimension has a form very similar to (44), d CDT S (σ) ≃ (b + 4σ)/(c + σ) > 0 [50,52,53], where b, c > 0 are constants found numerically. Calculating the entropy density for this profile on a d- This expression is finite only if b/(2c) = 1, which is avoided if the theoretical values of b and c are irrational. Again, one goes through an analytic continuation of the parameters.

Conclusions: role of dimensional flow and discreteness
In this paper, we have determined a necessary condition to obtain a finite entanglement entropy density ρ d : the spectral dimension of the spatial boundary of the spacetime geometry must not be positive at any scale. If it does, then ρ d diverges. We have seen the negative example of theories with compact momentum space (necessary condition on non-vanishing spectral dimension violated) 9 and the example of multi-fractional spacetimes, successful only after augmenting the definition of ρ d with an analytic continuation. We also considered the case (33) of continuous scale-invariant geometries (constant spectral dimension), where ρ d diverges.
When d b S < 0, P d (σ) ∼ σ −d b S /2 grows with σ and the spectral dimension makes little sense both as the scaling of the return probability of a diffusing process (which should decrease with the diffusion time) and as an indicator of spacetime dimensionality (which should be positive semi-definite). A physical interpretation of why a vanishing spectral dimension is problematic could be that, when d S → 0, the boundary effectively collapses to a point. Trying to match the quantum entanglement and statistical-mechanics interpretations of the entropy density, one might identify this type of degenerate configuration as the responsible for an ill-defined counting of degrees of freedom on the boundary, eventually resulting in an infinite entropy density. However, the boundary configuration may be much more complicated than a degenerate point. If quantum geometry hinders diffusion of a particle too much, for instance when its structure at some given scale is discrete or too coarse, then the spectral dimension of spacetime vanishes at that scale and the entropy density blows up. The condition d b S = 0 implies a constant return probability P d (σ) ∼ σ −d b S /2 = const at the diffusion scales of the zero plateau, meaning that a test particle does not diffuse on the boundary at those scales. This may happen if the boundary geometry is too irregular or disconnected. Therefore, a finite entropy can arise only in geometries sufficiently smooth, non-degenerate, or not pathologically rough. Note that one cannot define the concept of "irregularity" as a synonym of "nowhere differentiable," since nowhere-differentiable geometries such as those in the multi-fractional theory T γ with fractional derivatives correspond to multi-fractal dimensional profiles where d S = 0 and, as the above calculation for its approximation T q shows, the entanglement entropy should be finite.
We can start to draw some conclusions from the above results. We begin to recognize dimensional flow as playing an important role in the mechanism responsible for a finite entropy density. However, as the compact momentum space examples indicates, dimensional flow per se is not sufficient to get ρ d < ∞. These geometries are multi-scale but ρ d is infinite. The same conclusion holds for non-local quantum gravity, as we have seen in the previous section, and unless one considers interactions. On the other hand, string theory avoids the vanishing-d S problem, while multi-fractional theories and CDT work only if we allow for an analytic continuation. Therefore, there exist theories (all multiscale, i.e., with varying d S ) based on either a nowheredifferentiable or an ordinary differentiable structure that may realize a finite entanglement entropy density. Differentiability is not the main requirement here, the only common element being that the spacetime described by these theories is multi-scale. All these examples are on a continuum and one may wonder whether discreteness renders the entropy density finite.
-A graphical instance of the virtues of discreteness is provided by deterministic fractals. These sets are characterized by a discrete scale invariance and have a constant spectral dimension. Their spectral zeta function, however, does not vanish at s = 0 [5,16] and, thus, their entanglement entropy density has a divergent term [35], as in the case (33) of continuous scale-invariant geometries. Thus, theories of quantum gravity based on a multi-fractal geometry and solving the entropy-density problem might work but not because of an irregular spacetime texture. -Multi-fractional spacetimes have a measure with precisely such properties: both T q and T γ are field theories on a multi-fractal spacetime with a UV discrete scale invariance. However, as we commented in the previous section, discreteness is not necessary to get ρ d < ∞, since a finite ρ d can be obtained in the continuum regime of the theory T q (coarse-grained log oscillations). At the same time, a finite ρ d is also curiously related to a region of the parameter space of the multi-fractional theory T γ corresponding to geometries where a continuum description is flawed (0 < α < 1/2), but this does not strictly imply that discreteness must enter the game. -Discreteness plays an important role also in quantum gravities based upon pre-geometric combinatorial structures, such as loop quantum gravity, spin foams and their mother theory GFT. These theories exhibit a finite number of degrees of freedom on boundaries, which reproduce the entropy-area law for black holes [56][57][58][59]. This statistical-mechanics interpretation of the entropy density (in terms of state counting) should agree, at some point, with the interpretation of the same quantity as due to the quantum entanglement between states inside and outside the boundary or horizon. The results of [8] are encouraging in this respect, but we are now in the position to make an observation that jeopardizes this hope. Numerical investigations of the properties of a wide class of kinematical states of quantum ge-ometry appearing in GFT-spin foams-loop quantum gravity found that the UV is dominated by effects of the underlying discrete combinatorial structure, so that d S → 0 at short scales [10]. Therefore, according to our necessary condition, these quantum geometries do not have a finite entanglement entropy density, in contradiction with the suggestions of [8] but in syntony with independent arguments against a resolution of the information-loss problem (presumably related to the entanglement-entropy problem) in loop quantum gravity [60]. Before reaching a final verdict, however, we would like to emphasize that: (a) keeping in mind our cautionary remarks made at the end of Sect. 2, we should recall that our conclusion is based on the assumption that the relation (3) between entanglement entropy and the return probability of a boundary diffusion process remains valid in the particular quantum-geometry setting under consideration; (b) a resolution of this apparent contradiction might lie in the fact that, for example, the states analyzed in [10] are kinematical; states in the dynamical Hilbert space, so far left alone due to technical difficulties, might fare better. The limit d S → 0 is caused by a discreteness artifact of the theory, the combinatorial structure of complexes. In a way not really clear to us at this stage, it may be possible that such an artifact is smoothed out in fully dynamical geometric states. In other words, if we accepted that no physical measurement can take place at the discreteness scale of underlying simplicial complexes, then the regime d UV S = 0 would be unreachable for all purposes, thus avoiding our no-go result. In the case of black holes, this condition would amount to forbid objects with a Planckian event horizon. Future studies in loop quantum gravity should attempt to bypass the necessary condition found here and to address this issue.
-Field theories with a hard UV cut-off implemented as a finite spatial "bandwidth" based on Shannon sampling theory [61] share various features with discrete lattice models, including a minimal wavelength for the field oscillators. Recent studies show that also in these theories the presence of a UV cut-off does not lead to a finite entanglement entropy even though the field degrees of freedom in such model seem to occupy an incompressible volume [62]. From our perspective, it would be of interest to study the information density carried by the field degrees of freedom in the various models analyzed in this work along the lines of [62], to see which features are responsible for a finite entanglement entropy density.
To conclude, the main agent responsible for a finite entropy density in quantum gravity may be dimensional flow. This factor is not sufficient by itself, as shown by compact-momentum-space and non-local examples. The role of discreteness is less clear, but it might turn out to be a liability rather than an asset. String theory and the multi-fractional example in the absence of log oscillations show that discreteness is not necessary, while conflicting evidence in loop quantum gravity and GFT is questioning discreteness as an efficient agent. To better understand the relationships in this love triangle among entanglement entropy, varying geometry and discrete or irregular structures, we will have to wait for further advances. The most crucial element of disturbance, the positive sign of the spectral dimension forbidding a finite ρ d , kills virtually all quantum-gravity and string-related models. Changing the definition of the entanglement entropy, as we did minimalistically by an analytic continuation, might be advisable but it should be better understood. The use of the spectral zeta function promoted in this work is a promising starting point from where to refocus the problem.