Cosmic voids and filaments from quantum gravity

Using computer simulations, we study the geometry of a typical quantum universe, i.e., the geometry one might expect before a possible period of inflation. We display it using coordinates defined by means of four classical scalar fields satisfying the Laplace equation with nontrivial boundary conditions. They are a close analogue of the harmonic coordinate condition used in the context of GR (Kuchar and Torre in Phys Rev D 43:419–441, 1991). It is highly nontrivial that these ideas can be applied to understand the structures which appear in very irregular and fluctuating geometries. The field configurations reveal cosmic web structures surprisingly similar to the ones observed in the present-day universe.

Introduction -One major unsolved problem in theoretical physics is how to unite the theory of general relativity and quantum mechanics.It is hoped that such a unification will allow us to understand physics at the Planck scale, where the assumed quantum nature of gravity most likely plays a dominant role.Furthermore, the idea of an inflationary period in the history of our Universe has taught us that these quantum fluctuations at or near the Planck scale can, owing to an exponential growth of the size of the Universe, freeze and be expanded into seeds for macroscopic large-scale structures.Results presented in this article suggest that the scenario of cosmic filaments and voids that we observe in the Universe today might be much more generic in quantum theories of geometry.
Lattice Quantum Universes -In order to discuss universes of the size of a few Planck lengths and their fluctuating quantum geometry one needs a non-perturbative model of quantum gravity.We will discuss here a particular model called Causal Dynamical Triangulations (abbreviated to CDT; see [1] for a comprehensive introduction and an explanation of its somewhat technical name), but we believe that our results are generic and will be present in any reasonable quantum model of gravity.In CDT, spacetime is a triangulation built by joining together fixed-size four-dimensional simplices in a way that satisfies certain topological requirements.The edge length of the four-simplices acts as an ultraviolet (UV) cutoff.Its choice also fixes the geometry of a triangulation.A natural geometric way of calculating the classical Einstein-Hilbert action on such piecewise linear manifolds leads to the so-called Regge action.The lattice regularized path integral of quantum gravity is then given by (1) * ambjorn@nbi.dk.† zbigniew.drogosz@doctoral.uj.edu.pl‡ jakub.gizbert-studnicki@uj.edu.pl§ andrzej.goerlich@uj.edu.pl¶ jerzy.jurkiewicz@uj.edu.pl* * nemeth.daniel.1992@gmail.comwhere M H is a globally defined hyperbolic Lorentzian manifold and D M H [g L ] denotes the integration over equivalence classes [g L ] of Lorentzian metrics on M H . T E is a suitable set of Wick-rotated Euclidean triangulations.The Regge action S R [T E ] for a triangulation T E ∈ T E contains the bare couplings related to the cosmological and Newton constants.In principle, we want to adjust the bare coupling constants such that we can take the UV cutoff to zero while keeping physics unchanged (see [2] for a recent review).In accordance with the imposed global hyperbolicity, CDT introduces a time foliation of the four-dimensional manifolds into threedimensional leaves, which are three-dimensional spatial sub-manifolds with a global time t and a fixed topology.The explicit CDT construction permits a Wick rotation of the time coordinate t (still denoted t), whereby CDT becomes a statistical model that can be studied using Monte-Carlo simulations.This allows us, inter alia, to measure the time dependence of some global quantities, such as the spatial volume.It was shown [3][4][5] that, for suitable choices of the bare coupling constants, both the average spatial volume and its fluctuations can, with a large degree of accuracy, be described by the Hartle-Hawking minisuperspace model, which assumes isotropy and homogeneity of the Universe such that the only dynamical variable is the scale factor a(t).It should be emphasized that the isotropy and homogeneity is not put in by hand in CDT but follows from integrating out all degrees of freedom other than the scale factor, and the typical geometries encountered in the quantum path integral are not at all close to the classical homogeneous solution of GR.The approximate agreement with the classical minisuperspace solution is obtained from an average of an ensemble of highly fluctuating geometric states and is caused by a non-trivial interplay between the physical action and the entropy of configurations.
The CDT model is background independent and, in the spatial directions, coordinate free.There is no background geometry in the definition of the path integral.However, a good choice of coordinates can be very useful for a description of geometry [6][7][8], and below we will introduce a coordinate system that is suitable for describing the highly fluctuating geometries we encounter in the path integral.
Boundary Conditions -We will now consider a version of CDT where the piecewise linear manifolds are periodic both in time and space directions.Such a toroidal topology can be viewed as a periodically repeated fourdimensional elementary cell, bounded by a set of four independent non-contractible three-dimensional boundaries.These boundaries are not physical entities, are not unique and can be locally deformed as long as they still form an elementary cell, and yet they can serve as a reference frame for a coordinate system on a given triangulation.The non-trivial fractal structure of the encountered geometries makes it difficult to introduce spatial coordinates in a chosen elementary cell in a constructive geometric way.However, below we will show how to use four massless classical scalar fields, which satisfy the Laplace equation with non-trivial boundary conditions, to parametrize the fractal geometry.The fields depend on the geometry, but they do not modify it.One can say that they act as a microscope which uncovers the complicated four-dimensional structure of density fluctuations.Even for a configuration with a very irregular geometry, such fields allow us to define periodic pseudo-continuous coordinates and provide a generalized foliation in all space-time directions.Consequently, it becomes possible to visualize and measure multi-dimensional correlations in all directions.As will be reported below, what we see is a quantum universe which seems surprisingly similar to our present day macroscopic Universe.
Scalar fields as coordinates with values on S 1 -We want to find nontrivial harmonic maps between two Riemannian manifolds M(g µν ) → N (h αβ ), where g µν is an arbitrary metric and h αβ is a flat one.If N has the topology of T 4 , then it can be defined by four scalar fields φ α , α = 1, 2, 3, 4, where φ α (x) is a map M → S 1 , such that the following action is minimized: Because we have chosen the trivial metric h ρσ on N , eq. ( 2) splits up into four independent equations for the four scalar fields φ σ .Minimizing eq. ( 2) leads to the set of Laplace equations: Thus φ σ becomes a harmonic map M → (S 1 ) 4 .Let us consider a trivial one-dimensional example.In this case let M be S 1 with a unit circumference and a positive and strictly periodic density g(x).We want x → φ(x) to be a non-trivial map S 1 → S 1 such that φ can serve as a coordinate instead of x.One way to implement this is to find a solution satisfying which maps the circle with a unit circumference to a circle with a circumference δ.The solution to the Laplace equation in this case satisfies By rescaling the field, we can always enforce δ = 1.The solution φ(x) is fixed by picking x 0 where φ(x 0 ) = 0.The map x → φ(x) becomes a monotonically increasing invertible map in the whole domain R.If we parametrize the one-dimensional manifold M in terms of φ instead of x, we will find the volume density in the range φ to φ + dφ to be proportional to g(x) dx, so that effectively g(φ) = 1.We can also consider a function ψ(x) = mod(φ(x) − φ(x 1 ), δ).This function satisfies the Laplace equation in the range between x 1 (where ψ(x) = 0) and x 1 + 1 (where ψ(x) = δ).The equation satisfied by ψ(x) becomes a Poisson equation with the extra inhomogeneous local term, producing jumps at boundary points x = x 1 and x = x 1 + 1.It can still be considered to be a Laplace equation with a non-trivial boundary "jump" condition.Generalizing this to M with the topology of T 4 , we want a solution to the Laplace eq.
(3) that wraps around S 1 in a particular direction once, and, in addition, we want the points x in M that satisfy φ σ (x) = c to form hypersurfaces H σ (c) whose union for c varying in a range of length 1 covers the whole M.
Classical scalar fields with a jump -In CDT, the fourdimensional manifolds are represented by regular fourdimensional triangulations constructed by gluing together four-simplices so that each face is shared by exactly two simplices.Denote the number of four-simplices in the triangulation by N 4 .Each triangulation we consider is generated by a Monte Carlo simulation, using the CDT partition function.We call such a generated triangulation a configuration.For each configuration we keep information about the position of the four boundaries of the elementary cell.This information is non-dynamical: it does not influence the Monte Carlo process.Each boundary is a connected set of three-dimensional faces, each of which separates two simplices, for instance i and j, belonging to the two neighboring elementary cells.The connection i → j can have either positive or negative orientation, depending on the direction in which we cross the boundary.The boundary between the neighboring elementary cells in a direction σ can be parametrized by the N 4 × N 4 anti-symmetric matrix B σ ij = −B σ ji with the elements The number of directed boundary faces of a simplex i is given by b σ i = j B σ ij , with the obvious constraints −5 < b σ i < 5 and i b σ i = 0.For any simplex i adjacent to a boundary, the values B σ ij are all positive or zero (on one side of the boundary), or all negative or zero (on the other side).We consider four scalar fields φ σ i located in the centers of simplices and solve the minimization problem for the following discrete version of the continuous action in eq. ( 2), for each field φ σ i : In (7) the sum is over all pairs of neighboring foursimplices in the triangulation T E representing the manifold M(g µν ) in eq. ( 2).The parameter δ plays the same role as in the one-dimensional example considered previously, and here too, by rescaling the field, we can always set δ = 1.The action (7) has two important symmetries.The first one is the invariance under a constant shift of the scalar field (the Laplacian zero mode).The second is a local invariance under a modification of the boundary B σ ij and a shift by ±1 (depending on the side of the boundary) of the field value in a simplex i adjacent to the boundary.This is equivalent to moving the simplex to the other side of the boundary and compensating for the change of the field in its center.After such a move, the number of faces belonging to the boundary will in general be changed, but the action (7) will remain constant.
The classical field, henceforth denoted as φ σ i , minimizes the action (7), and thus has to satisfy the nonhomogeneous Poisson-like equation where L = 51 − A is the N 4 × N 4 Laplacian matrix, and A ij is the adjacency matrix with entries of value 1 if simplices i and j are neighbors and 0 otherwise.The Laplacian matrix L has a constant zero mode but can be inverted if we fix a value of the field φ σ i0 = 0 for an arbitrary simplex i 0 .Although L is a sparse matrix, inverting it is a major numerical challenge for a system typically of size N 4 ≈ 10 6 .After we nevertheless obtain the classical solution φ σ , the multi-dimensional analogue of the one-dimensional function ψ(x) is given by ψ σ i : A new boundary is defined by bσ = Lψ σ .This allows us to reconstruct a new three-dimensional hypersurface H, separating the elementary cell from its copies in the direction σ and characterized by the fact that the field jumps from 0 to 1 when crossing H.This hypersurface can be moved to another position if we consider a family of hypersurfaces H(α σ ) obtained from Changing 0 ≤ α σ < 1, we shift the position of the hypersurface and cover the whole elementary cell defined by the boundary (10), and in this way we obtain a foliation in the direction σ.We may now use ψ σ i = ψ σ i (0) as a coordinate in the σ-direction.The same construction can be repeated in all directions σ ∈ {x, y, z, t} for any configuration obtained in the numerical simulations, and in this way every simplex i will be assigned a unique set of coordinates {ψ x i , ψ y i , ψ z i , ψ t i }, all in the range between 0 and 1.A solution to the Laplace equation has the property that coordinates of each simplex are equal to the mean value of coordinates of its neighbors (up to the shift of the field at the boundary), i.e., it preserves the triangulation structure.This is the required map from our configuration with a topology of T 4 to (S 1 ) 4 (which of course also has the topology T 4 ).Note that the coordinate ψ t i is not the same as the one coming from the original foliation of the CDT model.The parametrization defined above permits to analyze the distribution of the four-volume (the number of simplices) contained in hypercubic blocks with sizes {∆ψ x i , ∆ψ y i , ∆ψ z i , ∆ψ t i }, which is equivalent to measuring the integrated g(ψ): We can measure the full four-dimensional distribution N (ψ).In Figs. 1 and 2 we show projections of the volume density distribution of a typical configuration on two-dimensional parameter subspaces, the xy-plane and the tx-plane respectively, integrating over the remaining two directions.One observes a remarkable pattern of voids and filaments, which qualitatively looks quite similar to pictures of voids and filaments observed in our real Universe (see e.g.[9]; the plots can be found on the web [10,11]).Using the new coordinates, we observe a pattern of volume concentrations in the spatial directions.The higher-density domains tend to attract each other, forming denser clouds, which survive in time evolution (see Fig. 2).There seems to appear a sequence of scales, characterizing a gradual condensation of gravitationally interacting "objects", but one should remember that there is no matter in this system, only pure geome-FIG.2. The projection of four-volume on the tx-plane for a CDT configuration.There is a strong correlation between the original t-foliation (color) and new time coordinate ψ t (horizontal axis).
try, which behaves as if quantum fluctuations could produce massive interacting objects, somewhat analogous to dark matter.Of course, we are talking about quantum objects of Planckian size, but if a more extended model exhibited inflation, one could imagine that aspects of these objects would be frozen when entering the horizon, like the standard Gaussian fluctuations in simple inflationary models, and then would re-enter the horizon at a later stage, after reheating, as classical densities.
Discussion -CDT presents us with a model of what we believe are generic fluctuations of geometry at the Planck scale.Hopefully, a detailed analysis of these fluctuations can be used in inflationary models, and such a study is now possible thanks to the coordinates we introduced here in spacetimes with toroidal topology.Apart from then being able to discuss the nature of geometric fluctuations before inflation we hope that measurements of such correlations will allow us to determine "experimentally" (i.e., using Monte Carlo data) the effective continuum action that governs our lattice model.The construction of such an effective action will help us to understand if CDT is an UV-complete quantum field theory of gravity, as imagined in the so-called asymptotic safety scenario, or only an effective quantum theory of geometries.

FIG. 1 .
FIG. 1.The projection of four-volume on the xy-plane, as defined by (11) for a CDT configuration.Different colors correspond to different times t of the original t-foliation.