Lattice Quantum Gravity and Asymptotic Safety

Transcription

1 Lattice Quantum Gravity and Asymptotic Safety Jack Laiho (Scott Bassler, Simon Catterall, Raghav Jha, Judah Unmuth-Yockey) Syracuse University June 18, 2018

2 Asymptotic Safety Weinberg proposed idea that gravity might be Asymptotically Safe in 1976 [Erice Subnucl. Phys. 1976:1]. This scenario would entail: Gravity is effectively renormalizable when formulated non-perturbatively. Problem lies with perturbation theory, not general relativity. In a Euclidean lattice formulation the fixed point would show up as a continuous phase transition point, the approach to which would define a continuum limit.

3 Lattice gravity Euclidean dynamical triangulations (EDT) is a lattice formulation that was introduced in the 90 s. [Ambjorn, Carfora, and Marzuoli, The geometry of dynamical triangulations, Springer, Berlin, 1997] Lattice geometries are approximated by triangles with fixed edge lengths. The dynamics is contained in the connectivity of the triangles, which can be added or deleted. In lattice gravity, the lattice itself is a dynamical entity, which evolves in Monte Carlo time. The dimension of the building blocks can be fixed, but the effective fractal dimension must be calculated from simulations. EDT works in 2d, where it reproduces the results of non-critical bosonic string theory. The EDT formulation in 4d was shown to have two phases, a collapsed" phase with infinite Hausdorff dimension and a branched polymer phase, with Hausdorff dimension 2. The critical point separating them was shown to be first order, so that new continuum physics is not expected. [Bialas et al, Nucl. Phys. B472, 293 (1996), hep-lat/ ; de Bakker, Phys. Lett. B389, 238 (1996), hep-lat/ ]

4 Einstein Hilbert Action Continuum Euclidean path-integral: Z = Dg e S[g], (1) S[g µν ] = k 2 d d x det g(r 2Λ), (2) where k = 1/(8πG N ).

5 Discrete action Discrete Euclidean (Regge) action is S E = k 2V 2 δ λ V 4, (3) where δ = 2π θ is the deficit angle around a triangular face, V i is the volume of an i-simplex, and λ = kλ. Can show that ( S E = 2 πkn 2 + N 4 2 k arccos 1 ) λ (4) where N i is the total number of i-simplices in the lattice. Conveniently written as S E = κ 2 N 2 + κ 4 N 4. (5)

6 Measure term Continuum calculations suggest a form for the measure Z = Going to the discretized theory, we have β Dg det g e S[g], (6) x N β 2 det g x j=1 O(t j ) β, (7) where O(t j ) is the order of triangle t j, i.e. the number of 4-simplices to which a triangle belongs. Can incorporate this term in the action by taking exponential of the log. β is a free parameter in simulations. Can interpret as an ultra-local measure term, since it looks like a product over local 4-volumes.

7 New Idea Revisiting the EDT approach because other formulations (renormalization group and other lattice approaches) suggest that gravity is asymptotically safe. New work done in collaboration with students (past and present) and postdoc: JL, S. Bassler, D. Coumbe, Daping Du, J. Neelakanta, (arxiv: ). Key new idea that inspired this study is that a fine-tuning of bare parameters in EDT is necessary to recover the correct continuum limit. This is in analogy to using Wilson fermions in lattice gauge theory to study quantum chromodynamics (QCD) with light or massless quarks. Striking similarities are seen. Coincidence? Previous work did not implement this fine-tuning, leading to negative results.

8 Main problems to overcome Must show recovery of semiclassical physics in 4 dimensions. Must show existence of continuum limit at continuous phase transition.

9 Phase diagram EDT vs. QCD with Wilson fermions β A C Branched Polymer Phase κ 2 Collapsed Phase Crinkled Region B κ D EDT β QCD

10 Three volume distribution k 8k 16k n 4 (ρ) ρ

11 Three volume distribution β=1.5 β=0 β=-0.6 β = -0.8 n 4,r (ρ r ) ρ r

12 Spectral Dimension χ 2 /dof=1.25, p-value=17% D S ( ) = ± 0.041, D S (0) = ± D S (σ) k σ

13 Infinite volume, continuum extrapolation χ 2 /dof=0.52, p-value=59% D S ( ) = 3.94 ± 0.16 D S (V,a) β=0 β=0.8 β= /V

14 Relative lattice spacing β=-0.8 β=-0.6 β=0 β=0.8 β= β= 0.8 β=-0.6 β=0 β=0.8 β=1.5 P(σ) 0.03 P(σ r ) σ σ r Return probability left and rescaled return probability right.

15 Causal dynamical triangulations Euclidean de Sitter space solution from arxiv: (Ambjorn et al.)

16 Semiclassical fluctuations Looking at quantum fluctuations about de Sitter space allow one to fix M Planck. A simple minisuperspace model fits the CDT data well. Ambjorn et al. (arxiv: ) look at the correlator C(i,j) = 1 K k (N (k) 3 (i) N 3(i))(N (k) 3 (j) N 3(j)), (8) where one can show that C(i,j) G N. The size of these quantum fluctuations compared to the width of the de Sitter universe can be used to fix the lattice spacing.

17 Causal dynamical triangulations Semiclassical fluctuations about de Sitter, arxiv: (Ambjorn et al.)

18 Semiclassical fluctuations from EDT 4k, β = (C N4 (i,i)/n 4 ) 1/ x

19 Semiclassical fluctuations from EDT 8k, β = (C N4 (i,i)/n 4 ) 1/ x

20 Finite volume effects 2.5 β=0 2 (C N4 (i,i)/n 4 ) 1/ /V

21 Finite volume effects β=0 β=1.5 β=0.8 β=-0.6 β=-0.8 (C N4 (i,i)/n 4 ) 1/ /V

22 Relative lattice spacing 1.2 P(σ) fluctuations 1 a/a fid κ 2

23 Relative lattice spacing The agreement between these two approaches for obtaining the relative lattice spacing bolsters the interpretation of the theory as quantum gravity, and strengthens the link between EDT and CDT. Could they be in the same universality class?

24 Adding matter We have looked at adding both scalar and fermion matter fields. Adding scalars and U(1) gauge fields dynamically (unquenched) was done already over twenty years ago. For now we are revisiting things in the quenched approximation, where matter loops are neglected, since this allows us to reuse existing lattice ensembles.

25 Scalar field In the continuum we can add to the Einstein Hilbert action the action for the scalar field: S[g,φ] = d 4 x ( 1 g 2 gµν µ φ ν φ + 1 ) 2 m2 0φ 2, (9) Lattice discretization is straightforward. Looking at scalar propagators on quenched configurations. This is being done by Judah Unmuth-Yockey and Raghav Jha.

26 Scalar propagator k, β =0, m 0 = origins 5 origins G(r) r G(r) = ( + m 2 0) 1 0r (10)

27 Mass dependence k, β = m r m 0 The shift symmetry of the action ensures that m r 0 as m 0 0.

28 Z m /Z φ factor 2.0 β =0 1.5 m r /m k 8k 16k m 0 m r /m 0 = (Z φ /Z m ) 1 2, (11) γ m = 1 2 d ln(z m /Z φ ) d lna (12)

29 Fermions Will leave this to Judah Unmuth-Yockey to discuss during his talk. (Work done in collaboration with Judah and Simon Catterall.) We have an approach to fermions that has not been looked at in the context of dynamical triangulations. Numerically, the formulation passes some nontrivial tests.

30 Conclusions Important to test the picture presented here against other approaches, renormalization group and other lattice formulations. If this holds, lattice provides a nonperturbative definition of a renormalizable quantum field theory of gravity. Beginning to probe the matter sector (in the quenched approximation) by looking at effects of the geometry on scalar and fermion fields.

31 Back-up Slides

32 Simulations Methods for doing these simulations were introduced in the 90 s. We wrote new code from scratch. The Metropolis Algorithm is implemented using a set of local update moves. We introduce a new algorithm for parallelizing the code, which we call parallel rejection. Exploits the low acceptance of the model, and partially compensates for it. Checked that it reproduces the scalar code configuration-by-configuration. Buys us a factor of 10.

33 Visualization of geometries Coarser to finer, left to right, top to bottom.

34 Diffusion process and the spectral dimension Spectral dimension is defined by a diffusion process d log P(σ) D S (σ) = 2 d logσ, (13) where σ is the diffusion time step on the lattice, and P(σ) is the return probability, i.e. the probability of being back where you started in a random walk after σ steps.

35 Infinite volume, continuum extrapolation χ 2 /dof=0.17, p-value=84% D S (0) = 1.44 ± β=1.5 β=0.8 β=0 2 D S (V,a) /V

36 What does it mean? Interesting results that suggest that the correct classical result might be restored in the continuum, large-volume limit. Analogy with Wilson fermions that inspired this study may tell us more. We have to perform a fine-tuning, and long distance physics gets messed up by discretization effects. These things happen when the regulator breaks a symmetry of the quantum theory. In this case, natural to identify the symmetry as continuum diffeomorphism invariance. If true, then β would not be a relevant parameter in a continuum formulation with diffeo invariance unbroken. (Would still need a measure term if the regulator broke scale invariance, but β would be fixed.) Interesting because if true, number of relevant couplings in continuum theory could be less than 3.