Arxiv:1909.08412V4

arXiv:1909.08412v4 [gr-qc] 20 Feb 2020 hrmsknky 5 ocw 129 Russia 117259, Moscow, 25, Cheremushkinskaya ∗ aumeeg nqatmfil hoymyb on in the found of be Calculation may theory field [18]. quantum value in energy observed magnitude vacuum the of its of orders than value many the larger is at (of arrive that cutoff we constant mass) ultraviolet cosmological Plank the the Taking of of order value vacuum. the finite of di- the energy account ultraviolet cosmo- into the the called to give - field contributions fluctuations so quantum vergent vacuum the the in the of Namely, content theory the problem. is constant value logical appearance a The such appears. of constant cosmological of value Ein- [15–17]. been of in has entanglement emergence considered quantum has possible the gravity from The [14] equations In stein entropy. to 13]. thermody- [12, linked from in been discussed equations in been Einstein example, has of namics for Emergence discussed, been [11]. has concept gravity the emergent to of correspondence AdS/CFT the of Relation 10]). [9, also (see su- model [8] to the gravity recently called quantum used were be the they do may and deformations systems crystals, Such elastic perplastic the energy. if require [7] not appear gravitational transfor- would the coordinate mations of general the form under the Invariance and have theory. disclinations not of the does theory in dislocations However, dynamical for the [2–6]. background materials excitations the real various as of serves motion which Riemann the emergent geometry, to (discli- Cartan it rise lattices earlier - give crystal Even dislocations) the of and [1]. defects nations Sakharov the that A.D. gravity known emergent by was on written at papers first been dynamically the has not appear of One is to energies. is theory low therefore, a transformations. invariance, Such coordinate This general theory. the approximation under energy microscopic invariant low certain the a as to appears gravity of ory nlaeo bec rmNC”ucao nttt”-IE,B. ITEP, - Institute” ”Kurchatov NRC from absence of leave On yial ntemdl iheegn rvt h huge the gravity emergent with models the in Typically gravity. emergent of scenarios other many exist There the- quantum the gravity emergent of scenario the In mretgaiyi uepatccytl n cosmologica and crystals superplastic in gravity Emergent o hne ftesse smdfidsotl.Tecosmologi The numbers: PACS smoothly. modified is coord well. v on system as invariant whole depends the topological the weakly if over integrated metric changed being background not that the way, that a this provided In such gravity. in to defined coupled be fermions (massive) gapped the of edsuseegn rvt ntesprlsi rsas W crystals. superplastic the in gravity emergent discuss We .INTRODUCTION I. 1 hsc eatet re nvriy re 00,Israel 40700, Ariel University, Ariel Department, Physics Dtd eray2,2020) 21, February (Dated: ..Zubkov M.A. oooia nains nSc.I eedwt h con- the with end we IV are clusions. trace- Sect. tensor In the energy invariants. stress that fermionic topological prove - the Wigner to of the us components describe less allows super- we which the III calculus, in Section Weyl gravity In emergent crystals. of plastic description the theory.recall field quantum [28–31] the formalism Weyl to - [32–36] Wigner in the techni- adapted use the will we On we energy, side, modi- vacuum cal constant. system. smooth vanishing cosmological the a with vanishing world is obtain of the world variation of observed smooth fication the the that to Assuming, there- robust constant, is cosmological The fore, system the a smoothly. when is modified changed, it not is volume is whole It the invariant. such over topological integrated in defined being that be way may a tensor energy that coordinates, stress on demonstrate fermionic depends We the weakly metric tensor. background component energy the trace if - the the stress cosmological of as vacuum The form appears of pattern the relativity. this receive in general may constant of gravitational it the action energies to Einstein low rise discli- At gives the and action. dynamics dislocations Their the and nations. by torsion emergent caused at The excita- are look fermionic crystal. curvature to the superplastic us a at in allows as tions fermions scenario Model This super- Standard the 8]. the of [7, model the crystal consider in We plastic scenario problem. gravity constant emergent the cosmological the of tion found. therein). been not references has solution and complete the [20–27] the far in example, so However, proposed for been cos- (see, the have past solve problem possibly, constant may, that mological scenarios Several [19]. ie rsa ycoordinates by crystal given ie eetetre( three the here sider h ae sognzda olw.I eto Iwe II Section In follows. as organized is paper The solu- alternative an propose we paper present the In e sprmtrz h netre atc fthe of lattice unperturbed the parameterize us Let ∗ .Eatcdfrain n testensor stress and deformations Elastic A. 1 I LSI EOMTOSAND DEFORMATIONS ELASTIC II. prahtesrs nrytno may tensor energy - stress the approach etitorevsb h consideration the by ourselves restrict e nts neulbimisvleis value its equilibrium In inates. lm ti oooia invariant topological a is it olume a osati hspteni a is pattern this in constant cal MRETGRAVITY EMERGENT D osatproblem constant l )dmninlsse.The system. dimensional 3) = x k k , 1 = . D .., , ewl con- will We . 2 displacement of each point has three components ua(x), time in the units of ∆t. The elasticity tetrads are gradi- where a =1, 2, 3. The resulting coordinates of the crystal ents of the phase function lattice sites are a a Eµ (y)= ∂µX (y), a =1, ..., D + 1 (6)

yk(x) = xk + uk(x), k =1, 2, 3 (1) a µ The inverse matrix to Eµ is denoted by ea . For a = In the absence of the displacements, when u = 0, the 1, .., D those quantities have units of crystal momentum, a a µ [Eµ ]=1/[l], [ea ] = [l] while for a = D + 1 the quantity crystal is unperturbed. Metric of elasticity theory is a given by E has the dimension 1/[t]. In the absence of dislocations, the tetrads Ea = a µ 1 Eµ dy are exact differentials. They can be expressed gik = δik +2uik, uik = ∂iuk + ∂kui + ∂iua∂kua (2), 2 in general form in terms of a system of deformed crys- a =1, 2, 3, i,k =1, ..., D (3) tallographic coordinate planes, i.e. surfaces of constant phase Xa(y)=2πna, na Z with a =1, ..., D. ∈ a The stress tensor T ij is defined as the response of the In the absence of dislocations, tensor Eµ (y) satisfies thermodynamical potential to elastic deformation. the integrability condition:

a 1 a a µ ν ij δ dE = (∂ E (y) ∂ E (y))dx dy =0 . (7) T (x)= log Z 2 µ ν − ν µ ∧ −δuij(x) µ a Matrix ea inverse to Eµ represents the inverse vielbein. It determines the distance along the curve that connects B. Elasticity tetrads the two points in space - time correspondingC to the two values of X (X1,X2) as follows: The alternative way to describe elastic deformations is given in terms of the elasticity tetrads E a(x), which ρ (X ,X )= eµeν δ dXadXb µ C 1 2 Z a b µν represent the hydrodynamical variables of the elasticity C q theory [4, 9, 10]. where δµν = diag (1, 1, 1, 1). In dimensionless coordinates The three crystallographic planes may be chosen in the X metric is defined as three dimensional (D = 3) crystals1. The coordinates µ ν Xa (a =1, ..., D) are introduced in such a way, that the µ ν ∂y ∂y gab(X)= e e δµν = δµν planes given by equations Xa(x)=2πna, na Z with a b ∂Xa ∂Xb a = 1, ..., D represent the mentioned crystallographic∈ It is useful to introduce the time displacement uD+1 0. planes. In the presence of elastic deformations the chosen This allows to express metric as ≡ crystallographic planes become curved. The corresponding displacement vectors satisfy the following equation ∂uµ ∂uν g (X)= [e(0)]µ + [e(0)]ν + δ ab a ∂Xa b ∂Xb µν Xa(x + u)=2πna where The intersections of the D surfaces µ µ (0) µ ∂x (0) −1 [e ]a = a = [E ] X1(r,t)=2πn1 , ... , XD(r,t)=2πnD , (4) ∂X a and coincide with the positions of certain crystal lattice sites a (not necessarily all of them may be found in this way). (0) a ∂X [E ]µ = Those intersection points constitute the rectangular lat- ∂xµ tice Metric in coordinates X is related to metric in coordi- D a D nates x as follows: L = r = R(n1, ..., nD) r R ,n Z . (5) { | ∈ ∈ } ∂Xa ∂Xb g (X)= g (x)= δ + u In order to describe dynamics we define all quantities as ∂xµ ∂xν ab µν µν µν dependent on imaginary time t. (We assume from the very beginning that the Wick rotation to space - time of that is Euclidean signature is performed.) The extra dimension- a b (0) (0) less coordinate XD+1 = t/∆t is defined as the imaginary E E gab(X)= gµν (x)= δµν + uµν h iµh iν Recall also, that in coordinates y metric is given simply by 1 Those crystallographic planes are not necessarily orthogonal for (0) the unperturbed crystals. gµν (y)= δµν 3

ij ..ν Similar to the stress tensor T the response of the ther- vκ = Tκν modynamical potential to the elasticity tetrads gives rise i to the tensor Θa: ..λ ..λ λ ..ρ λρ tµν = Tµν + δ[µTν]ρ g Tρ[µν] δ − Θi (x)= log Z a −δEa(x) For certain particular values of α,β,γ we will obtain the i so called teleparallel equivalent of general relativity. In the similar way We consider the quantum vacuum as the spacetime superplastic crystal, where Xa(y)=2πna, na Z with δ a a = 1, 2, 3, 4 is the system of the four deformed∈ crystal- Θi (x)= i log Z −δea(x) lographic coordinate hyperplanes in 3 + 1 D. One may Those tensors are related as follows: assume that in the superplastic vacuum, the eﬀective ac- a tion for Eµ is the gravitational one given by Eq. (9). ij 1 {j i} ab 1 a k Metric g , which originates from the elasticity tetrads, T = eb Θa δ ,Tij = gk{j Θi}ea µν 2 2 is:

g = δ E aE b . (10) C. Teleparallel gravity µν ab µ ν

In the presence of the single dislocation (placed at the line ξ) the coordinates X become ambiguous. The viel- a D. Riemann - Cartan geometry bein Eµ being integrated along the closed path sur- rounding ξ gives the Burgers vector specific for theC given lattice and the given type of the dislocation 2: Above we discussed the case, when the torsion is present while the curvature is absent. This occurs in the dyµEa = Ba (8) presence of dislocations. The disclinations give rise to the Z µ C curvature. This occurs as follows. The disclination corre- In the differential form this equation is given by sponds to the removal of a sector ended at the given line. The remaining parts of the lattice are glued together. a a a ρ (3) Therefore, the path around the given line captures the ∂µE (y) ∂ν E (y)= B ǫµνρ dξ (τ)δ (y ξ(τ)),D =3 ν − µ Z − angle different from 2π. This means that the curvature is nonzero. In the unperturbed lattice each path corre- The dislocations distributed continuously give the viel- sponds to vanishing parallel transporter. In the presence bein Ea that (being averaged) corresponds to torsion µ of a disclination the parallel transporter along the closed ..a a a Tµν (y)= ∂µEν (y) ∂ν Eµ(y) path, which encloses the disclination, is equal to the ma- − trix of rotation. The distribution of disclinations may a We may speak of the field Tµν as of the continuous func- be described by continuous SO(3) SO(4) connection. tion of coordinates. This connection modifies derivative⊂ as usual: it becomes a The curved geometry given by the vielbein Eµ is, actu- covariant. In the absence of dislocations we deal, there- ally, the geometry of the teleparallel gravity. This is the fore, with Riemannian geometry. version of Riemann - Cartan gravity with vanishing cur- If both the disclinations and the dislocations are vature (and spin connection), and nonzero torsion. The present, there is the emergent Riemann - Cartan geome- corresponding space is parallelizable Weitzenbock space. a try. The basic variables are the vielbein Eµ, and the spin The possible terms in the action of the dynamical the- connection ωµ so(D+1), where D = 3. The mentioned ory of this type of gravity may be found in [37]. If we above parallel transporter∈ around the disclination results restrict ourselves by the terms with the derivatives up to in the rotation matrix: the second, and if there is no selected direction in D +1 µ φlΣˆ dimensional space - time, while the P and T symmetries P exp ωµdy = e (11) hold, then we are left with ZC Here φ is the angle specific for the given disclination, S = dD+1y det E λ + αv2 + βa2 + γt2 (9) Z while l is the unity vector along the axis of rotation, while Σˆ is the generator of rotation. As usual, curvature Rµν Here a,v,t are the irreducible components of torsion: is defined as 1 aκ = ǫκαβγT µ 1 µ ν 4 det E αβγ lim|A|→0P exp ωµdy = exp Rµν dy dy Z∂A 2 ZA ∧ The covariant derivative is given by

2 Its value is quantized according to the lattice type Dµ = ∂µ + ωµ 4 while the definition of torsion is modified as follows: elastic deformations, which are not caused by the appearance of the disclinations and the dislocations, do not ..a a a Tµν (y)= DµEν (y) Dν Eµ(y) require energy. As a result, the low energy effective the- − ory appears to be invariant under the diffeomorphisms. There are a lot of possible local terms in the action that This is the case of the superplastic crystal. contain derivatives up to the second. Those terms may Thus, we may suppose that the Standard Model with be composed in different way of torsion, curvature, and gravity appears as the low energy description of the su- various covariant derivatives of the vielbein. These terms perplastic crystal with the fermionic excitations inside were classified in [37]. it. The real gravity appears as the theory of the elastic vielbein. The presence of the dislocations results in torsion while the presence of disclinations results in cur- E. The fields experienced by the fermions in the presence of dislocations and disclinations vature. Notice, that the mentioned here tight - binding model may be thought of as an effective model, while the true microscopic theory remains continuous. Below we Let us consider the tight - binding model with Hamil- will concentrate on the formalism that allows us to deal tonian with such a microscopic theory. In this pattern of space - ¯ A AB B time the time still remains distinguished from space. The H = ΨX1 tX1X2 ΨX2 (12) X curvature and the torsion do not have nontrivial time X1X2 components. The description of the space - time gravity Hopping parameters tAB correspond to the jumps of and the Standard Model fermions as an effective theory X1X2 may be thought of as the theory written in the reference the fermions between the sites X1 and X2 of the lattice. D+1 D+1 1 The simplest dependence of the hopping parameters on frame with Eµ = δµ ∆t . This is the theory with the elastic deformations is the partially fixed gauge. (The symmetry with respect to the space - time diffeomorphisms is reduced down to AB AB AB µ ν a a b b the diffeomorphisms of space.) tX1X2 = t + t β ea e δµν (X1 X2 )(X1 X2) q b − − There is an alternative way to look at the problem (0) µ (0) ν a a b b of the distinguished (imaginary) time direction. Notice, [e ]a [e ]b δµν (X1 X2 )(X1 X2) (13) a −q − − that Eµ enters all expressions in such a way, that there is no essential difference between the components with Here β is the so - called Gruneisen parameter. Suppose a,µ =1, 2, 3 and a (or µ) equal to 4. We may extend the that the model admits the topologically protected Fermi notion of the elastic deformations (and the corresponding points. In their vicinities the emergent Weyl fermions ap- displacement vectors u ) from the 3 - component vectors pear (this is the case of the so - called Weyl semimetals). µ defined in 3 - dimensional space to the 4 - component In the absence of multiple dislocations and disclinations vectors defined in 4 - dimensional space - time. We may the above mentioned modification of the hopping param- actually think that space - time is the four - dimensional eters results in the appearance of emergent gauge field superplastic crystal. This restores the complete symme- and emergent teleparallel gravity. The emergent gauge try between space and (imaginary) time components of field appears in the first approximation to the complete all vectors and tensors. The dislocations will appear as theory, while the emergent gravity appears in the next the two - dimensional defects, such that for the closed approximation. The emergent gravity is characterized by contours that enclose them, the integrals of Eq. (8) are torsion while the emergent gauge field is characterized by nontrivial. The disclinations also appear as the two - di- its field strength. The presence of multiple dislocations mensional defects. The integrals over the contours that and disclinations modifies the emergent gravity experi- surround them given by Eq. (11) are nontrivial. Σ rep- enced by the fermions. Namely, it has been shown that resents the generators of the SO(4) group (recall that we the dislocations carry the emergent singular torsion flux perform the Wick rotation to the Euclidean signature of and the emergent singular U(1) flux. In the similar way space - time). one may demonstrate, that the disclinations carry the singularity of curvature. The especially interesting case is when the emergent torsion carried by the dislocations is much larger than the III. CONTINUUM MICROSCOPIC THEORY ON THE LANGUAGE OF WIGNER - WEYL emergent torsion caused by the above mentioned modifi- FORMALISM cation of the hopping parameters. This occurs when the Gruneisen parameter β is small β 1 3. In this situation we deal with the Weyl fermions≪ that experience the A. Wigner - Weyl formalism described above gravity with torsion. In this case those In this section we consider the Wigner - Weyl formalism to be used for the description of the underlined microscopic theory. For more details of this formalism we 3 Notice that in graphene, for example, β ∼ 2. refer to [32–35]. Let us start from the discussion of the 5 model with the noninteracting fermions in D + 1 space As it was mentioned above, we made here rotation to - time dimensions. The corresponding partition function space - time of Euclidean signature, i.e. we consider the has the form theory in imaginary time. Here QW depends on parameter η, and we introduced the Wigner transformed Green dD+1y(Ψ¯ QˆΨ) Z = DΨ¯ DΨeR (14) function G as the Weyl symbol of Gˆ, where Z W where operator Qˆ contains an information about the microscopic physics including the periodic potential of the fermions in the crystal. It is useful to introduce the 1 Tr ˆ ˆ G(x, y)= x Gˆ y = DΨ¯ DΨe W QΨ(¯ x)Ψ(y) (Grassmann - valued) distribution function: h | | i Z Z W (x, y)= x Wˆ y = Ψ(¯ x)Ψ(y) h | | i we also denote the matrix elements of operator Qˆ by As usual in field theory, in order to calculate this inte- Q(x, y)= x Qˆ y h | | i gral in practise we need to put the system into the large, As a result, the partition function may be rewritten as but finite spatial volume. Correspondingly, the imagi- follows nary time extent (the inverse temperature) also should be large, but finite. Then, strictly speaking, the sums Tr ˆ ˆ Z = DΨ¯ DΨe W Q (15) over the discrete momenta appear instead of the integral. Z Though, we are able to calculate the integrals instead Here in the exponent the functional trace appears. The of the sums because the overall volume is large. Corre- next step is the introduction of the Weyl symbol of op- spondingly, the boundary conditions are to be imposed. erator Aˆ: The experience with the conventional quantum field theory prompts that for the large volume (and the large −iyp imaginary time extent) the type of the boundary con- AW (p, x) dye A(x + y/2, x y/2) (16) ≡ Z − dition is irrelevant for the physics in the bulk. For the ˆ ˆ definiteness we assume the periodical spacial boundary We introduce the Weyl symbols of operators Q and W , conditions, and the anti - periodic boundary conditions and come to for the imaginary time direction. dD+1p Z = DΨ¯ DΨexp dD+1x Z Z Z (2π)D+1 The Wigner transformed Green function obeys the Groenewold equation: TrWW [Ψ, Ψ](¯ p, x) ⋆QW (p, x) (17) We denote here by the star (Moyal) product the following expression

i ←− −→ ←−−→ QW ⋆ GW =1 = e 2 ∂ x∂p−∂p ∂ x ∗

The properties of the star product guarantee the com- B. Derivative of the partition function as a mutativity topological invariant

Let us suppose that operator Qˆ introduced above depends on a parameter η. We consider derivative of the dD+1p logarithm of the partition function: dD+1x TrA (p, x) ⋆B (p, x) Z Z (2π)D+1 W W D+1 D+1 d 1 D+1 d p d p log Z = DΨ¯ DΨ d x = dD+1x TrB (p, x) ⋆ A (p, x) dη Z Z Z Z (2π)D+1 Z Z (2π)D+1 W W D+1 D+1 d p = TrAˆBˆ (19) TrWW [ψ, ψ¯](p, x) ⋆QW (p, x) d x Z Z (2π)D+1 d TrWW [ψ, ψ¯](p, x) ⋆ QW (p, x) dη D+1 Tr D+1 d p d Here is the functional trace of an operator that obeys = d x TrGW (p, x) ⋆ QW (p, x)(18) ˆ ˆ ˆ ˆ Z Z (2π)D+1 dη TrAB = TrBA. Therefore, since δ(QW ⋆ GW ) = 0, the 6

d variation of dη log Z is equal to zero: for the given class of variations δQW (p, x). Notice, that the theories containing gapless fermions in the presence d dD+1p δ log Z = dD+1x of interactions suffer from the infrared divergencies, we dη Z Z (2π)D+1 do not consider them here. d d An example, when the derivative of partition function Tr GW (p, x) ⋆ δQW (p, x)+ δGW (p, x) ⋆ QW (p, x) dη dη with respect to a parameter is topological invariant is given by the electric current density integrated over the dD+1p D+1 whole system volume. In this case the role of parameter η = d x D+1 Z Z (2π) is played by constant external electromagnetic potential d Ak along the k - th axis. Then in Eq. (21) we substitute Tr GW (p, x) ⋆ δQW (p, x) d ∂ dη by , and this condition is satisfied if we chose ap- dη − ∂pk d propriate boundary conditions in momentum space. The GW (p, x) ⋆δQW (p, x) ⋆ GW (p, x) ⋆ QW (p, x) − dη topological invariance of the total electric current derived in this way represents a field theoretical generalization of dD+1p = dD+1x the quantum - mechanical Bloch theorem [39]. Z Z (2π)D+1 d Tr GW (p, x) ⋆ δQW (p, x) dη C. Contribution of fermions to the stress - energy d tensor GW (p, x) ⋆δQW (p, x) ⋆ GW (p, x) ⋆ QW (p, x) − dη Here we will consider example of the topological in- dD+1p d D+1 variant given by the traceless components of stress - en- = d x D+1 Tr GW (p, x) ⋆ δQW (p, x) Z Z (2π) dη ergy tensor integrated over the whole volume. In the d absence of nontrivial background metric this stress ten- + GW (p, x) ⋆δQW (p, x) dη sor is given by response of partition function to the vari- D+1 ations of elastic deformations δuij (x). In the presence D+1 d p d = d x Tr GW (p, x) ⋆δQW (p, x) (20) of the background metric originated from distribution of Z Z (2π)D+1 dη disclinations and dislocations we consider the response We come to an interesting conclusion that the deriva- of thermodynamical potential to the variation of effective of (the logarithm of) the partition function is topo- tive metric given by logical invariant provided that the last integral in Eq. ij k ij ik j jk i (20) vanishes. δg (x)= δu ∂kg (x)+g (x)∂kδu (x)+g (x)∂kδu (x) In the above proof it has been implied that the inte- − grals are convergent. The absence of the infrared diver- and effective vielbein gencies means that neither Gˆ nor Qˆ have poles at finite δei (x)= δuk∂ ei (x)+ ej (x)∂ δui(x) values of momenta. As for the ultraviolet divergencies, a − k a a j we assume that the discussed theory does not contain ij them. In particular, we may imagine, that this theory In the considered theory the response to variation of g represents stress - energy tensor Tij (x). Response to the has the form of a lattice model (the imaginary time is i discretized as well). The formalism described in Sect. variation of ea gives tetrad components of the stress - a II, which leads to the Groenewold equation, formally has energy tensor Θµ. In our consideration it is more con- been derived for the continuous theory. It remains valid venient to deal with the latter. Relations between those for the lattice models as well if the external field (in our two tensors are given by case this is the elastic deformation tensor) varies slowly, 1 so that its variation at the distance of the order of the Θj =Θaej ,T = (Θkg +Θkg ) i i a ij 2 i kj j ki lattice spacing may be neglected (for the details see [32– 34]). This condition is assumed in the theory of elasticity. The total stress - energy tensor is given by an integral of j a j a However, we do not rely on the assumption that the con- eaΘi (x) over the whole volume averaged in time eaΘi : sidered theory is defined on the lattice. This is only one h i of the possible options. For our purposes it is enough to Θj = T dD+1x ej Θa(x) E(x) (22) assume that the theory is free from the ultraviolet diver- h i i Z h a i | |i gencies. =2T dD+1z E(z) dD+1x Thus Eq. (18) at zero temperature is topological in- Z | | Z variant for the gapped systems that do not contain ul- dD+1p δ traviolet divergencies provided that Tr G (p, x) ⋆ej (z) Q (p, x) Z (2π)D+1 W a δei (z) W dD+1p d a dD+1x Tr G (p, x) ⋆δQ (p, x)=0 Z Z (2π)D+1 dη W W Here T is temperature, which is assumed to be small. If (21) QW and GW do not depend on (imaginary) time, then 7

T dxD+1 becomes equal to unity. In Eq. (22) we ig- For the variation of stress - energy tensor we get noreR the possible contribution to stress - energy tensor that may appear due to the dependence of measure over dD+1p δΘj = T dD+1x E(x) Tr fermionic fields on E(x) . This contribution will be dis- h i i Z | | Z (2π)D+1 cussed later. | | ∂ j Using technique developed in Appendix A it is shown pi δi GW (p, x) ⋆δQW (p, x) in Appendix C that we are able to substitute here ∂pj − D+1 j δ 1 ∂ j i d p e a by p δ δ(z x) provided that e de- D+1 a δe (z) 2 j ∂pi i a T d x E(x) Tr i − − Z Z D+1 pend on coordinates only slightly. In Appendix B we − | | (2π) show, in addition, that the covariant derivative of stress a µ ∂ j Eµδea GW (p, x) ⋆ pi δi QW (p, x) - energy tensor is a topological invariant. ∂pj − The effective low energy theory for the model with a Fermi point gives rise to the following expression for QW One can see, that under a suitable ultraviolet regularization the traceless components of stress - energy ten- E(y) µ a i QW (p,y)= | | e (y) γ , (pµ Aµ(y) ωµ(y)) (23) sor Θ (those with i = j) are topological invariants 2 a { − − } j provided that E(x) remains6 constant during the given provided that the Gruneisen parameter is much smaller variations. This| may| always be achieved via a suitable than unity. This situation may be realized in the ideal reparametrization that accompanies the modification of superplastic crystals, where the elastic deformations do the system. not require energy. Here ω is the emergent spin connec- A few comments are in order. First of all, we still did tion resulted from the disclinations, while the elasticity not take into account the possible dependence of inte- vielbein leads to the presence of torsion. In this case the gration measure on E . We suppose that there is no space components of the emergent stress - energy tensor dependence on the other| | components of ei . Such a de- appear to be identical to the stress tensor of elasticity. a pendence may influence only the trace components of the As it was mentioned above, we assume that the appear- stress - energy tensor. It will be taken into account in ance of the Fermi point is broken by a small mass term. section III E following the requirement that the whole As a result the fermions become gapped, and the sys- stress - energy tensor is a topological invariant. Next, tem already does not suffer from the infrared singularities the term in the brackets in Eq. (26) proportional to the (which appear for any gapless fermions in the presence product δj Q (p, x) does not appear if we use the naive of interactions). Thus we come to i W expression for the stress - energy tensor (it is propor- E(y) µ a tional to the averaged fermion lagrangian, it vanishes on Q (p,y)= | |e (y) γ , (p A (y) ω (y)) m E(y) W 2 a { µ− µ − µ }− | | the classical solutions). The whole expression Eq. (26) (24) suffers from all possible types of divergencies - both ul- It was noticed above that Eq. (24) represents the low en- traviolet and infrared. The infrared divergency is to be ergy approximation for the theory with the Fermi point regularized through the finite overall volume (and finite (broken weakly by a small mass term). For momenta far temperature). The naive ultraviolet regularization (via from the Fermi point the form of Qˆ is essentially differ- the ultraviolet cutoff in the integral over momentum) ent. The whole theory may be actually considered as the will then give the ultraviolet divergent expression, which regularization of the relativistic theory with fermions and does not have a form of topological invariant. As it was gravity. mentioned above, we need the ultraviolet regularization Taking into account all mentioned above, we may try that obeys the requirement that in the integrals over mo- to consider the following expression as the definition of mentum p, its shift by a constant a may be performed: the stress - energy tensor (again, up to the possible con- p p+a. The superplastic crystal by itself represents an tribution from the integration measure): ultraviolet→ regularization that satisfies this requirement. It is worth mentioning, that metric is a quantity, which Θj = T dD+1x Θj (x) E(x) (25) h i i Z h i | |i is well defined in continuum theory. Correspondingly the D+1 stress energy tensor being the response of the system to D+1 d p = T d x E(x) Tr GW (p, x) ⋆ the variation of metric, is also well - defined in continuum Z | | Z (2π)D+1 theory only. ∂ j pi δi QW (p, x) ∂pj − Applying this expression formally to Eq. (24) we would D. Lattice regularization obtain: dD+1p On the lattice we need an appropriate definition, which Θj = T dD+1x E(x) Tr G (p, x) ⋆ h i i Z | | Z (2π)D+1 W tends to a continuum definition in continuum limit. One may try to take as the definition of stress - energy tensor j j a δi QW (p, x)+ E eaγ (pi Ai(x)) (26) on rectangular lattice (up to the terms coming from the − | | − 8 integration measure) the following expression: E. Cosmological constant as the topological invariant dD+1p Θj = T E(x) Tr G (p, x) ⋆ (27) i Z D+1 W It has been noticed above that the traceless compo- h i X | | (2π) j j j x nents θ =Θ 1 δ Θk of stress - energy tensor of Eq. i i − D+1 i k ∂ j (27) (and of its lattice regularization of Eq. (29)) are sin pi δi QW (p, x) ∂pj − robust to the smooth modifications of the system accom- panied by suitable reparametrizations aimed at keeping In this expression the lattice spacing is equal to unity E(x) unchanged. However, we did not take into account while momentum space is a torus (0; 2π]D+1. We assume the| contribution| to the stress energy tensor of the integra- here that all fields vary slowly, so that their variation at tion measure over the fermionic fields. This integration the distance of lattice spacing may be neglected. Un- measure is assumed to depend on the determinant E(x) a | | der these conditions sin pm commutes with ⋆. Under the and not on the other components of Eµ. Therefore, the same conditions we may substitute the sum over x by corresponding contribution to the stress - energy tensor an integral. We obtain the following expression for the cannot modify the traceless components θj . We propose j i variation of Θ (as above we assume that E(x) is not the hypothesis that the topological invariance is an es- i | | changed during this variation): sential property of the source of the gravitational excitations. And therefore, it is natural to require that not dD+1p only the traceless components, but also the whole stress - δΘj = T dD+1x Tr (28) h i i Z Z (2π)D+1 energy tensor is a topological invariant. The term coming from the integration measure are supposed to be respon- ∂ sin pi GW (p, x) ⋆δQW (p, x) sible for this. This hypothesis leads us to the redefinition ∂pj of the whole stress - energy tensor in such a way, that its trace is the topological invariant as well. This gives the One can see, that this expression is not a topological following expression invariant even for i = j. Therefore, we redefine the lattice 6 expression for the stress - energy tensor as follows dD+1p Θj = T dD+1x E(x) Tr G (p, x)(31)⋆ h i i Z | | Z (2π)D+1 W D+1 j d p ∂ Θ = T E(x) Tr GW (p, x) ⋆ (29) j h i i | | Z (2π)D+1 pi δi QW (p, x) Xx ∂pj − ∂ j j +T Tr δ LogQˆ (32) sin pi δi QW (p, x) i ∂pj − j +T Tr(cosp î 1)δ LogQˆ LogQˆ0 and its lattice regularized form − i − D+1 j d p ˆ Θ = T E Tr GW (p, x) ⋆ (33) Here Q0 is built with a certain particular fixed vielbein h i i | | Z (2π)D+1 i(0) Xx ea . By Log Qˆ we understand the functional Logarithm Tr ∂ j and by we understand the functional trace. Their sin pi δi QW (p, x) definitions and properties are discussed in Appendix D. ∂pj − To avoid the uncertainties in the definition we define both Tr j ˆ ˆ +T cos piδi LogQ LogQ0 in the gauge E(x) = 1. Their values for arbitrary E(x) − | | | | are to be calculated for the configuration of the vielbein By Log and Tr we understand here again the functional with E(x) = 1 connected by reparametrization to the | | logarithm and functional trace considered in Appendix given one. If the fields vary slowly at the distance of the D. With this new definition the fermionic stress - en- order of lattice spacing, and operation ⋆ commutes with ergy tensor integrated over all space is the topological cos p, then the variation of this expression is invariant, provided that the fermions are massive while background metric depends on coordinates weakly. dD+1p δΘj = T dD+1x E Tr (30) We apply this consideration to the fermions of the h i i Z | | Z (2π)D+1 Standard Model. Then the vacuum average Θj h i,vacuumi ∂ j appears to be the topological invariant if the background sin pi δi metric does not depend on coordinates. Phenomenolog- ∂pj − ical gravitational action, which appears at low energies G (p, x) ⋆δQ (p, x) W W as the effective description of the disclinations and the dislocations, is assumed to have the standard form One can see, that being defined in this way the traceless components of stress energy tensor are topological S = m2 R√ gd4x 2λm2 √ gd4x invariants. g − P Z − − P Z − 9

It contains the cosmological constant λ. Variation of the are responsible for the topological invariance of the trace corresponding term in the action with respect to gµν gives components of the stress - energy tensor. The appear- the following ”vacuum” stress - energy tensor: ance of those terms reflects difficulties with the definition of measure over the fermionic fields in the presence Θµ = λm2 δµ ν,vacuum − P ν of nontrivial background metric. In fact, measure over the fermionic fields in the superplastic crystal is trivial One can see, that λ is the topological invariant. Smooth for the trivial metric only. It is then defined as the in- modification of the system does not change value of the tegration over the Grassmann variables attached to the cosmological constant. Assuming that the smooth modi- lattice sites. The presence of nontrivial metric means fication connects our world with the world with vanishing that there are dislocations and disclinations. In the pres- vacuum energy, we come to the conclusion that the cos- ence of these defects it is not sufficiently clear ab initio mological constant λ vanishes. Its small deviation from how to define the integration measure. Thus, our above zero results, therefore, from the contributions that were definition assumes, that the measure is defined in such a not taken into account in our consideration. way, that its variation with respect to variation of metric gives the mentioned above extra terms in Eqs. (31) and IV. CONCLUSIONS (33), so that not only the traceless components of Θ, but also its trace components are the topological invariants. In the present paper we discuss the scenario, in which In the proposed scheme the bosonic fields of the Stan- the Standard Model of elementary particles and the dard Model were not taken into account (as well as in quantum gravity appear in the effective low energy de- [8]). There should be the reason, why this contribution scription of a condensed matter system. The latter sys- vanishes or appears to be small. In this sense the solu- tem describes the ideal superplastic crystal, in which tion of the cosmological constant problem proposed here smooth elastic deformations do not require energy. As is incomplete, and requires the further elaboration. How- a result, the emergent theory is invariant under the gen- ever, this remains out of the scope of the present paper eral coordinate transformations of space. In this theory 4. the disclinations result in curvature, while the disloca- The author is grateful for the discussions to tions result in torsion. The action, which describes the G.E.Volovik, who actually proposed him to look for the dynamics of dislocations and disclinations, may have at topological nature of the cosmological constant using an low energies the form of the Einstein action of general analogy with the superplastic crystals. relativity. Alternatively, the superplastic crystals may be considered as the ultraviolet regularization of quantum gravity coupled to fermions. The advantage of this regularization is that the integration measure over fermions does not depend on the gauge fields and on the gravitational fields. Appendix A. Variation of Weyl symbol of an Instead, the dependence of operator QW (x, p) on mo- operator with respect to reparametrizations menta is responsible for the quantum anomalies. There- fore, the variational derivative of the integration measure does not contribute the total stress energy tensor. Let us consider variation of Wigner transformed ma- We show that in the discussed theory the fermionic trix elements of an operator Aˆ caused by the following stress energy tensor defined properly is a topological in- transformation: variant, provided that the fermions are massive and the background metric depends on coordinates weakly. With our definition of the stress - energy tensor the cosmolog- A(y1,y2) A(y1 δu(y1),y2 δu(y2)) ical constant becomes topological invariant. Finally, if → − − our World appears to be a smooth modification of the system with vanishing vacuum energy, we come to vanishing value of the cosmological constant. The extremely 4 We only mention one of the possible ways to understand the small observable value of the cosmological constant may reason, why the gauge field contribution to the vacuum stress follow from the deviations of the system from this state energy tensor vanishes. One may accept the scheme, in which the gauge fields appear dynamically. Namely, the variation of related to dependence of metric on both time and space the position of the Fermi point describes the emergent gauge coordinates. One of such contributions may come from field [38]. The fluctuations of the Fermi point occur due to the the scale anomaly, i.e. the contribution due to curva- interactions between the fermions. Those interactions may be ture to the trace of the stress - energy tensor (for the caused by the emergent gravity. In principle, all gauge fields of discussion of such contributions see [40] and references the Standard Model (or of its ultraviolet completion) may appear in this way. Those fields appear as the composite excitations therein). existing at the relatively low energies (the energies much smaller The extra terms proportional to the functional trace than the Plank mass). At the same time, at high energies, we of the logarithm of Qˆ added to Eq. (31) and Eq. (33) deal with the system discussed in the present paper. 10 with small elastic deformation displacement δuk(x): variational derivative as

AW (p, x)+ δAW (p, x) δAW (p, x) dk 2i(z−x)k k = 2 4 e dy dy e−i(y1−y2)pA(y δu(y ),y δu(y )) δu (z) − Z (2π) ≡ Z 1 2 1 − 1 2 − 2 (i(kk + pk)+ ∂xk /2)AW (k + p, x) δ(x (y1 + y2)/2) (34) dk − 2 e2i(x−z)k − Z (2π)4 This gives ( i(k + p )+ ∂ k /2)A (k + p, x) − k k x W dk 2i(z−x)k+kj ∂pj −i(y1−y2)p = 2 e δAW (p, x)= dy1dy2e )δ(x (y1 + y2)/2) − Z (2π)4 Z − (i(p )+ ∂ k /2)A (p, x) k k k x W δu (y1)∂ k δu (y2)∂ k y1 y2 dk − − 2i(x−z)k+kj ∂pj 2 4 e A(y1,y2) − Z (2π)

dq −i(y1−y2)p ( i(pk)+ ∂xk /2)AW (p, x) = dy1dy2 e δ(x (y1 + y2)/2) 4 − ←− −→ Z (2π) − dk i ∂ ∂ /2 = 2 e2i(z−x)ke xj pj k k − Z (2π)4 δu (y1)∂ k δu (y2)∂ k − y1 − y2 (i(pk)+ ∂xk /2)AW (p, x) i(y1−y2)q e AW (q, (y1 + y2)/2) ←− −→ dk 2i(x−z)k −i ∂ j ∂ p /2 dq 2 e e x j = dy dy e−i(y1−y2)pδ(x (y + y )/2) − Z (2π)4 Z 1 2 (2π)4 − 1 2 ( i(pk)+ ∂xk /2)AW (p, x) k k − ←− −→ δu (y1)(iqk + ∂xk /2) δu (y2)( iqk + ∂xk /2) i ∂ j ∂ p /2 − − − = δ(z x)e x j (i(pk) − − i(y1−y2)q e AW (q, x) +∂xk /2)AW (p, x) ←− −→ dq i(y −y )(q−p) −i ∂ xj ∂ pj /2 = dy dy e 1 2 δ(x (y + y )/2) δ(z x)e ( i(pk) Z 1 2 (2π)4 − 1 2 − − − +∂xk /2)AW (p, x) k k δu (y1)(iqk + ∂xk /2) δu (y2)( iqk + ∂xk /2) = δ(z x) ⋆ (i(p ) − − − − − k AW (q, x) +∂xk /2)AW (p, x) (35) δ(z x)( i(pk) − − − +∂ k /2)A (p, x) ⋆δ(z x) x W − Therefore, Coming back to δA we may write it in the following δA (p, x) dq W W = dy ei(z−y2)(q−p)δ(x (z + y )/2) form: δuk(z) Z 2 (2π)4 − 2

( iqk ∂xk /2)AW (q, x) − − k dq δAW (p, x) = i δu (x) ⋆ pkAW (p, x) i(y1−z)(q−p) − + dy2 4 e δ(x (z + y1)/2) h i Z (2π) − k pkAW (p, x) ⋆δu (x) (iqk ∂xk /2)AW (q, x) −h i − dq 1 k = 2 e2i(z−x)(q−p) δu (x) ⋆ ∂xk AW (p, x) Z (2π)4 −2 k ( iq ∂ k /2)A (q, x) +∂xk AW (p, x) ⋆δu (x) (37) − k − x W dq +2 e2i(x−z)(q−p) Z (2π)4

(iq ∂ k /2)A (q, x) k − x W (36) Appendix B. Covariant derivative of stress - energy tensor in terms of Wigner transformed Green We assume that boundary conditions in momentum function space are organized in such a way, that we may shift the integration variable q q p = k. This may always be achieved via the appropriate→ − ultraviolet regularization. We know that with respect to reparametriza- As a result we are able to represent the above written tions y x(y) = y δu the action S = → − 11

D+1 D+1 ij k ij d y1d y2Ψ(¯ y1) y1 Qˆ[g] y2 Ψ(y2) is invariant: metric equal to δg = δu ∂kg . Response of log Z to R h | | i the variation δuk gives the− covariant derivative of stress - D+1 D+1 S = d x1d x2Ψ(¯ x1) x1 Qˆ[e] x2 Ψ(x2) energy tensor. According to Eq. (20) the corresponding j Z h | | i response of QW (p, x) to δu = const is = dD+1x(y )dD+1x(y )(Ψ(¯ x(y )) (0) k 1 2 1 δ Q (p, x) = δu ∂ k Q (p, x) (42) Z W − x W x(y1) Qˆ[g] x(y2) Ψ(x(y2)) Then h | | i D+1 D+1 ¯ ∂x = d y1d y2Ψ(x(y1)) (y1) D T = T dD+1x D T (x) g(x) (43) Z ∂y i ij Z i ij h i h ip| | ˆ ∂x dD+1p ∂ x(y1) Q[g] x(y2) (y2) Ψ(x(y2)) = T dD+1x Tr G (p, x) ⋆ Q (p, x) h | | i ∂y D+1 W j W − Z Z (2π) ∂δu D+1 D+1 ¯ D+1 = d y1d y2Ψ(x(y1)) D+1 d p Z = T d x Tr G (p, x) ⋆ ∂ j Q (p, x)(44) Z Z (2π)D+1 W x W x y1 Qˆ[e ] y2 Ψ(x(y2)) (38) h | | i According to Eq. (20) variation of the covariant deriva- Here ex is the corresponding transformation of vielbein. tive of total stress - energy tensor is given by Thus the reparametrizations inside Qˆ D+1 ˆ ˆ δ DiTij = T d xδ DiTij (x) g(x) (45) y1 Q y2 y1 δu(y1) Q y2 δu(y2) Z h | | i → h − | | − i h i h ip| | ∂δu ∂δu dD+1p ∂ 1 (y1) 1 (y2) (39) = T dD+1x Tr G (p, x) ⋆δQ (p, x) − ∂y − ∂y D+1 j W W − Z Z (2π) ∂δu D+1 compensate transformation Ψ(y) Ψ(y δu) and cor- D+1 d p → − = T d x ∂ j Tr G (p, x) ⋆δQ (p, x)(46) respond to the transformation of vielbein. Recall that Z Z (2π)D+1 x W W spinors are subject to SU(2) SU(2) transformations that are not related to the reparametrizations⊗ of coor- The last integral vanishes if we take appropriate bound- dinates. Spin connection is transformed as vector with ary conditions in coordinate space. Therefore, the covari- respect to the transformations of coordinates. However, ant derivative of stress - energy tensor integrated over the in Riemannian geometry the spin connection is expressed whole volume is robust to the smooth modiﬁcation of the through the vielbein, and therefore, its transformation is system. This means that the global Einstein anomaly is caused by the appropriate transformation of the latter. a topological invariant. (Recall that at D = 3 the ﬁeld Therefore, indeed the transformation of Eq. (39) is re- theory in conventional regularization predicts no Einstein i (Weyl) anomaly.) duced to the transformation of ea. This transformation gives:

i i QW (1 ∂iδu ) ⋆ (QW + δQW ) ⋆ (1 ∂iδu ) (40) Appendix C. Stress - energy tensor in terms of → − − Wigner - Weyl calculus where δQW has been calculated in Appendix A. Let us denote the resulting overall variation (that takes into ac- Equipped with Eq. (41) we derive here response of ˆ count both the transformation of matrix elements of Q partition function to the variation of an induced vielbein. ∂δu and the extra factors 1 ) by δ2. We get Recall, that the latter variation has the form − ∂y i k i k i i i δe (x)= δu ∂ke (x)+ e (x)∂kδu (x) δ2QW = δQW ∂iδu ⋆QW QW ⋆ ∂iδu a − a a − − i in the presence of background vielbein ea. First we cal- and j culate the response of δ2QW (p, x) to ∂kδu . This cal- k δ2QW = i δu (x) ⋆ pkQW (p, x) culation does not give an immediate expression for the − h i response to the vielbein variation. But it allows to calcu- k pkQW (p, x) ⋆δu (x) late it indirectly. In order to obtain this response we use −h i derivative expansion in Eq. (37). The term proportional 1 k δu (x) ⋆ ∂xk QW (p, x) to the ﬁrst order in derivatives gives: −2 k (1) k j +∂xk QW (p, x) ⋆δu (x) δ2 QW (p, x)= ∂pj pkQW ∂j δu (x) 2QW ∂j δu (x) − i i a l k a l j ∂iδu ⋆QW QW ⋆ ∂iδu (41) = ∂pj pkQW Ej ea∂lδu (x) 2QW Ej ea∂lδu (x) − − − First we demonstrate the power of the proposed tech- a k m k = ∂pj pkQW (p, x) Ej (δea(x)+ δu (x)∂mea(x)) nique and calculate response of the thermodynamical po- h i tential to constant δuk. Such a δuk causes variation of 2Q (p, x)Eaδ(ej (x)+ δum(x)∂ ej (x)) (47) − W j a m a 12

Let us suppose, that vielbein is given by a slight variation exp Tr Log Qˆ. One can also define this logarithm using a i i (0)i (0)i suitably regularized representation: of the constant one ea(x) = ha(x)+ ea , where ea = (0) const(x). We have for δ QW (p, x): ∞ (0) k ds D+1 δ Q (p, x) = δu ∂ k Q (p, x) Tr O (Log Qˆ Log Qˆ ) = lim d x 2 W − x W − 0 − ǫ→0 Z s Z δ ǫ k i ˆ ˆ δu dz∂ k h (z) Q (p, x) isQ isQ0 ≈ − Z x a δei (z) W x e e x (49) a Here O is a certain operator.h | Convergence− | i of the last k ˆ We know that δu may enter the expression for δQW only expression is guaranteed if all eigenvalues of Q are real, i and we add to each of them the small imaginary part iδ, through the variation δea. Combining the two leading terms in derivative expansion we get up to the terms which gives the mentioned above regularization. Thus it linear in h and its first derivatives the following equation: is assumed, that after the Wick rotation Qˆ is Hermitian. Under such a regularization the integral over s becomes a convergent at s . δ2QW (p, x) = ∂pj pkQW (p, x) Ej → ∞ h i k m k One can check that Eq. (49) is manifestly not invariant (δea(x)+ δu (x)∂mea(x)) under the reparametrizations. One can see this expand- a j m j 2QW (p, x)Ej δ(ea(x)+ δu (x)∂mea(x)) ing the exponent as a series in powers of Qˆ. An example − ˆ k i δ of Q given by Eq. (24) shows how the invariance under δu dz∂ k h (z) Q (p, x) − Z x a δei (z) W the reparametrizations is broken due to the product of de- a terminants E(x) . In order to have the expression invari- | | Solution of this equation gives ant under the reparametrizations we, therefore, have to fix the gauge E(x) = 1. Then, Eq. (49) gives the proper | | δ2QW (p, x) a a definition of the functional logarithm in this gauge. It has i = Ek pi∂pk Ek QW (p, x)δ(x z)(48) δea(z) h − i − the same value for the configurations related to the given one by reparametrizations, where already E(x) = 1. The total stress - energy tensor receives the form | | 6 Thus dealing with Eq. (49) we assume the above mentioned gauge fixing. If O commutes with Qˆ and δQˆ, one ej (x)Θa(x) = T dD+1xdD+1z E(z) a i Z can easily calculate (in the gauge E(x) = 1) h i | | | | D+1 d p j δ D+1 Tr GW (p, x) ⋆ea i QW (p, x) ∞ Z (2π) δea(z) ˆ δTr Oˆ Log Qˆ = i dD+1x x dsOeisQδQˆ x dD+1p − Z h | Z | i = T dD+1x E Tr G (p, x) 0 Z | | Z (2π)D+1 W ∞ = i ds¯ dD+1xeis¯ x OˆQˆ−1δQˆ x j − Z0 Z h | | i ⋆ pi∂pj δi QW (p, x) h − i = dD+1x x OˆQˆ−1δQˆ x This expression being understood naively is divergent Z h | | i both in ultraviolet and in infrared. Infrared regulariza- dD+1p = dD+1x O ⋆ G ⋆δQ (50) tion goes through the consideration of the system in large Z (2π)D+1 W W W but finite volume. The superplastic crystal by itself represents an ultraviolet regularization. It is worth mentioning, that metric is a quantity, which is well defined For OW = O(p) that is the function of p only the last in continuum theory. Correspondingly the stress energy expression may be easily generalized to its covariant form tensor being the response of the system to the variation adding the extra factor E(x) : of metric, is also well - defined in continuum theory only. | |

dD+1p δTr Oˆ Log Qˆ = dD+1x E(x) O(p) Appendix D. Functional logarithm Z | |(2π)D+1

GW ⋆δQW (51) In the main text the functional Logarithm is used. It enters expression for the partition function Z DetQˆ = ≡

[1] A. D. Sakharov, ”Vacuum quantum ﬂuctuations Dokl.Akad.Nauk SSSR 177,70-71 (1967) [Sov.Phys.Dokl. in curved space and the theory of gravitation,” 12,1040-1041 (1968). AlsoS14, pp. 167-169] Usp. Fiz. 13

Nauk 161,64-66 (1991) Progress in Physics 80, no. 2 (2017) 022001 [2] B.A. Bilby and E. Smith, Continuous distributions of [21] H. Martel, P.R. Shapiro, S. Weinberg, ”Likely Val- dislocations: A new application of the methods of non- ues of the Cosmological Constant”. The Astrophysical Riemannian geometry, Proc. Roy. Soc. Sect. A 231, 263– Journal. 492 (1) (1998) 2940, arXiv:astro-ph/9701099 273 (1955); Continuous distributions of dislocations. III, doi:10.1086/305016 Proc. Roy. Soc. Sect. A 236, 481–505 (1956). [22] P. Bull, Yashar Akrami, Julian Adamek, Tessa Baker, [3] E. Kröner, Allgemeine Kontinuumstheorie der Versetzon- Emilio Bellini, Jose Beltrn Jimenez, Eloisa Bentivegna gen and Ligenspannunge, Arch. Rational Mech. Anal. 4, et al, ”Beyond ΛCDM: Problems, solutions, and the road 18–334 (1960). ahead,” Physics of the Dark Universe 12 (2016): 56-99. [4] I.E. Dzyaloshinskii, and G.E. Volovick, Poisson brackets [23] Q. Wang, Z. Zhu, W. Unruh, ”How the huge en- in condensed matter, Ann. Phys. 125 67–97 (1980). ergy of quantum vacuum gravitates to drive the [5] H. Kleinert and J. Zaanen, World nematic crystal model slow accelerating expansion of the Universe”. Physi- of gravity explaining the absence of torsion, Phys. Lett. cal Review D. 95 (10) (2017) 103504 arXiv:1703.00543 A 324, 361–365 (2004). doi:10.1103/PhysRevD.95.103504 [6] Petr Jizba, Hagen Kleinert and Fabio Scardigli, Un- [24] G. Ellis, ”The trace-free Einstein equations and infla- certainty relation on world crystal and its applications tion”, General Relativity and Gravitation 46 (2014) 1619 to micro black holes, Phys. Rev. D 81, 084030 (2010), arXiv:1306.3021 doi:10.1007/s10714-013-1619-5 arXiv:0912.2253 [25] R. Peracci, ”Unimodular quantum gravity and the cos- [7] H. Kleinert, ”Gravity as Theory of Defects in mological constant,” Foundations of Physics. 48 (10) a Crystal with Only Second-Gradient Elastic- (2018) 13641379 arXiv:1712.09903 doi:10.1007/s10701- ity”. Annalen der Physik. 44 (2) (1987) 117, 018-0189-5 doi:10.1002/andp.19874990206. [26] S. J. Brodsky, C. D. Roberts, R. Shrock and P. C. Tandy, [8] F.R. Klinkhamer and G.E. Volovik, Tetrads and q-theory, ”Essence of the vacuum quark condensate”, Phys.Rev. arXiv:1812.07046. C82 (2010) 022201 [arXiv:1005.4610] [9] J. Nissinen and G.E. Volovik, Tetrads in solids: from S. J. Brodsky, C. D. Roberts, R. Shrock and P. C. elasticity theory to topological quantum Hall systems Tandy, ”Confinement contains condensates”, Phys.Rev. and Weyl fermions, ZhETF 154, 1051–1056 (2018), C85 (2012) 065202 [arXiv:1202.2376] arXiv:1803.09234. [27] O. Luongo, M. Muccino, ”Speeding up the Universe using [10] J. Nissinen and G.E. Volovik, 3+1d QHE, elasticity dust with pressure,” Physical Review D. 98 (10) (2018) tetrads and mixed axial-gravitational anomalies, Phys. 23 arXiv:1807.00180 doi:10.1103/physrevd.98.103520 Rev. Research 1, 023007 (2019), arXiv:1812.03175. [28] H. J. Groenewold, Physica, 12, 405 (1946) . [11] M. Van Raamsdonk, ”Building up spacetime with [29] J. E. Moyal, Proceedings of the Cambridge Philosophical quantum entanglement”, General Relativity and Grav- Society, 45, 99 (1949). itation. 42 (10) (2010): 23232329, arXiv:1005.3035, [30] F.A. Berezin and M.A. Shubin, in: Colloquia Mathemat- doi:10.1007/s10714-010-1034-0 ica Societatis Janos Bolyai (North-Holland, Amsterdam) [12] T. Jacobson, ”Thermodynamics of Spacetime: The p. 21, (1972). Einstein Equation of State”. Physical Review Let- [31] T. L. Curtright and C. K. Zachos, Asia Pacific ters. 75 (7) (1995) 12601263, arXiv:gr-qc/9504004, Physics Newsletter, issue 01, pages 37 − 46 (2012), doi:10.1103/PhysRevLett.75.1260 arXiv:1104.5269. [13] T. Padmanabhan, ”Thermodynamical Aspects of Grav- [32] M. A. Zubkov, Annals Phys. 373, 298 (2016). ity: New insights”, Reports on Progress in Physics. 73 [arXiv:1603.03665 [cond-mat.mes-hall]]. (4) (2010) 046901, arXiv:0911.5004, doi:10.1088/0034- [33] I. V. Fialkovsky, M. A. Zubkov, arXiv:1905.11097 4885/73/4/046901 [34] M. Suleymanov and M. A. Zubkov, “Wigner – Weyl [14] E. Verlinde, ”On the origin of gravity and the laws of formalism and the propagator of Wilson fermions Newton”. Journal of High Energy Physics (4)(2011) 29, in the presence of varying external electromagnetic arXiv:1001.0785, doi:10.1007/jhep04(2011)029 field,” Nucl. Phys. B 938 , 171 (2019) (Corrigen- [15] Jae-Weon Lee, Hyeong-Chan Kim, Jungjai Lee, ”Gravity dum: https://doi.org/10.1016/j.nuclphysb.2019.114674 ) from quantum information”, Journal of the Korean Phys- [arXiv:1811.08233 [hep-lat]]. ical Society. 63 (5) (2013) 10941098, arXiv:1001.5445, [35] M. A. Zubkov and Z. V. Khaidukov, “Topology of the doi:10.3938/jkps.63.1094 momentum space, Wigner transformations, and a chiral [16] B. Swingle, M. Van Raamsdonk, ”Universality of Gravity anomaly in lattice models,” JETP Lett. 106, 172 (2017) from Entanglement”, arXiv:1405.2933 [hep-th]. [Pisma Zh. Eksp. Teor. Fiz. 106 (2017) no.3, 166]. [17] E. Oh, I.Y. Park, S.J. Sin,”Complete Einstein equations [36] Z. V. Khaidukov and M. A. Zubkov, “Chiral from the generalized First Law of Entanglement”, Phys- torsional effect,” JETP Lett. 108 (2018) no.10, ical Review D. 98 (2) (2018) 026020, arXiv:1709.05752, 670 doi:10.1134/S0021364018220046 [arXiv:1812.00970 doi:10.1103/PhysRevD.98.026020 [cond-mat.mes-hall]]. [18] J. Martin, ”Everything you always wanted to know about [37] D. Diakonov, A. G. Tumanov and A. A. Vladimirov, the cosmological constant problem (but were afraid to “Low-energy General Relativity with torsion: A ask),” Comptes Rendus Physique 13.6-7 (2012) 566-665 Systematic derivative expansion,” Phys. Rev. D [19] Y.B. Zeldovich, ”Cosmological Constant and Elementary 84 (2011) 124042 doi:10.1103/PhysRevD.84.124042 Particles’, JETP letters 6 (1967), 316-317; ”The Cosmo- [arXiv:1104.2432 [hep-th]]. logical Constant and the Theory of Elementary Parti- [38] G.E. Volovik, The Universe in a Helium Droplet, Claren- cles”, Soviet Physics Uspekhi 11 (1968), 381-393. don Press, Oxford (2003) [20] A. Linde, ”A brief history of the multiverse,” Reports on [39] C. X. Zhang and M. A. Zubkov, “Note on the Bloch 14

theorem,” Phys. Rev. D 100 (2019) no.11, 116021 the QCD trace anomaly?,” Phys. Rev. Lett. 89 doi:10.1103/PhysRevD.100.116021 [arXiv:1909.12128 (2002) 081302 doi:10.1103/PhysRevLett.89.081302 [cond-mat.mes-hall]]. [gr-qc/0204018]. [40] R. Schutzhold, “Small cosmological constant from