1 Introduction

Die Grenzen meiner Sprache bedeuten die Grenzen meiner Welt.

(The limits of my language mean the limits of my world.)

Ludwig Wittgenstein

Tractatus logico-philosophicus

While general relativity is very successful in describing the gravitational interaction and the structure of space and time on large scales [205], quantum gravity is needed for the small-scale behavior. This is usually relevant when curvature, or in physical terms energy densities and tidal forces, becomes large. In cosmology this is the case close to the Big Bang, and also in the interior of black holes. We are thus able to learn about gravity on small scales by looking at the early history of the universe.

Starting with general relativity on large scales and evolving backward in time, the universe becomes smaller and smaller and quantum effects eventually become important. That the classical theory by itself cannot be sufficient to describe the history in a well-defined way is illustrated by singularity theorems [123] which also apply in this case: After a finite time of backward evolution the classical universe will collapse into a single point and energy densities diverge. At this point, the theory breaks down and cannot be used to determine what is happening there. Quantum gravity, with its different dynamics on small scales, is expected to solve this problem.

The quantum description does not only present a modified dynamical behavior on small scales but also a new conceptual setting. Rather than dealing with a classical space-time manifold, we now have evolution equations for the wave function of a universe. This opens a vast number of problems on various levels from mathematical physics to cosmological observations, and even philosophy. This review is intended to give an overview and summary of the current status of those problems, in particular in the new framework of loop quantum cosmology.

2 The Viewpoint of Loop Quantum Cosmology

Loop quantum cosmology is based on quantum Riemannian geometry, or loop quantum gravity [172, 22, 195, 174], which is an attempt at a non-perturbative and background independent quantization of general relativity. This means that no assumptions of small fields or the presence of a classical background metric are made, both of which is expected to be essential close to classical singularities where the gravitational field would diverge and space degenerates. In contrast to other approaches to quantum cosmology there is a direct link between cosmological models and the full theory [38, 66], as we will describe later in Section 6. With cosmological applications we are thus able to test several possible constructions and draw conclusions for open issues in the full theory. At the same time, of course, we can learn about physical effects which have to be expected from properties of the quantization and can potentially lead to observable predictions. Since the full theory is not completed yet, however, an important issue in this context is the robustness of those applications to choices in the full theory and quantization ambiguities.

The full theory itself is, understandably, extremely complex and thus requires approximation schemes for direct applications. Loop quantum cosmology is based on symmetry reduction, in the simplest case to isotropic geometries [46]. This poses the mathematical problem as to how the quantum representation of a model and its composite operators can be derived from that of the full theory, and in which sense this can be regarded as an approximation with suitable correction terms. Research in this direction currently proceeds by studying symmetric models with less symmetries and the relations between them. This allows to see what role anisotropies and inhomogeneities play in the full theory.

While this work is still in progress, one can obtain full quantizations of models by using basic features as they can already be derived from the full theory together with constructions of more complicated operators in a way analogous to what one does in the full theory (see Section 5). For those complicated operators, the prime example being the Hamiltonian constraint which dictates the dynamics of the theory, the link between model and the full theory is not always clear-cut. Nevertheless, one can try different versions in the model in explicit ways and see what implications this has, so again the robustness issue arises. This has already been applied to issues such as the semiclassical limit and general properties of quantum dynamics. Thus, general ideas which are required for this new, background independent quantization scheme, can be tried in a rather simple context in explicit ways to see how those constructions work in practice.

At the same time, there are possible phenomenological consequences in the physical systems being studied, which is the subject of Section 4. In fact it turned out, rather surprisingly, that already very basic effects such as the discreteness of quantum geometry and other features briefly reviewed in Section 3, for which a reliable derivation from the full theory is available, have very specific implications in early universe cosmology. While quantitative aspects depend on quantization ambiguities, there is a rich source of qualitative effects which work together in a well-defined and viable picture of the early universe. In such a way, as illustrated later, a partial view of the full theory and its properties emerges also from a physical, not just mathematical perspective.

With this wide range of problems being investigated we can keep our eyes open to input from all sides. There are mathematical consistency conditions in the full theory, some of which are identically satisfied in the simplest models (such as the isotropic model which has only one Hamiltonian constraint and thus a trivial constraint algebra). They are being studied in different, more complicated models and also in the full theory directly. Since the conditions are not easy to satisfy, they put stringent bounds on possible ambiguities. From physical applications, on the other hand, we obtain conceptual and phenomenological constraints which can be complementary to those obtained from consistency checks. All this contributes to a test and better understanding of the background independent framework and its implications.

Other reviews of loop quantum cosmology at different levels can be found in [56, 55, 199, 50, 69, 51, 96]. For complementary applications of loop quantum gravity to cosmology see [140, 141, 2, 114, 152, 1].

3 Loop Quantum Gravity

Since many reviews of full loop quantum gravity [172, 195, 22, 174, 161] as well as shorter accounts [9, 10, 173, 190, 167, 198] are already available, we describe here only those properties which will be essential later on. Nevertheless, this review is mostly self-contained; our notation is closest to that in [22]. A recent bibliography can be found in [93].

3.1 Geometry

General relativity in its canonical formulation [6] describes the geometry of space-time in terms of fields on spatial slices. Geometry on such a spatial slice Σ is encoded in the spatial metric qab, which presents the configuration variables. Canonical momenta are given in terms of extrinsic curvature Kab which is the derivative of the spatial metric under changing the spatial slice. Those fields are not arbitrary since they are obtained from a solution of Einstein’s equations by choosing a time coordinate defining the spatial slices, and space-time geometry is generally covariant. In the canonical formalism this is expressed by the presence of constraints on the fields, the diffeomorphism constraint and the Hamiltonian constraint. The diffeomorphism constraint generates deformations of a spatial slice or coordinate changes, and when it is satisfied spatial geometry does not depend on which coordinates we choose on space. General covariance of space-time geometry also for the time coordinate is then completed by imposing the Hamiltonian constraint. This constraint, furthermore, is important for the dynamics of the theory: Since there is no absolute time, there is no Hamiltonian generating evolution, but only the Hamiltonian constraint. When it is satisfied, it encodes correlations between the physical fields of gravity and matter such that evolution in this framework is relational. The reproduction of a space-time metric in a coordinate dependent way then requires to choose a gauge and to compute the transformation in gauge parameters (including the coordinates) generated by the constraints.

It is often useful to describe spatial geometry not by the spatial metric but by a triad \(e_i^a\) which defines three vector fields which are orthogonal to each other and normalized in each point. This yields all information about spatial geometry, and indeed the inverse metric is obtained from the triad by \({q^{ab}} = e_i^ae_i^b\) where we sum over the index i counting the triad vector fields. There are differences, however, between metric and triad formulations. First, the set of triad vectors can be rotated without changing the metric, which implies an additional gauge freedom with group SO(3) acting on the index i. Invariance of the theory under those rotations is then guaranteed by a Gauss constraint in addition to the diffeomorphism and Hamiltonian constraints.

The second difference will turn out to be more important later on: We can not only rotate the triad vectors but also reflect them, i.e., change the orientation of the triad given by sgn det \(e_i^a\). This does not change the metric either, and so could be included in the gauge group as O(3). However, reflections are not connected to the unit element of O(3) and thus are not generated by a constraint. It then has to be seen whether or not the theory allows to impose invariance under reflections, i.e., if its solutions are reflection symmetric. This is not usually an issue in the classical theory since positive and negative orientations on the space of triads are separated by degenerate configurations where the determinant of the metric vanishes. Points on the boundary are usually singularities where the classical evolution breaks down such that we will never connect between both sides. However, since there are expectations that quantum gravity may resolve classical singularities, which indeed are confirmed in loop quantum cosmology, we will have to keep this issue in mind and not restrict to only one orientation from the outset.

3.2 Ashtekar variables

To quantize a constrained canonical theory one can use Dirac’s prescription [105] and first represent the classical Poisson algebra of a suitable complete set of basic variables on phase space as an operator algebra on a Hilbert space, called kinematical. This ignores the constraints, which can be written as operators on the same Hilbert space. At the quantum level the constraints are then solved by determining their kernel, to be equipped with an inner product so as to define the physical Hilbert space. If zero is in the discrete part of the spectrum of a constraint, as e.g., for the Gauss constraint when the structure group is compact, the kernel is a subspace of the kinematical Hilbert space to which the kinematical inner product can be restricted. If, on the other hand, zero lies in the continuous part of the spectrum, there are no normalizable eigenstates and one has to construct a new physical Hilbert space from distributions. This is the case for the diffeomorphism and Hamiltonian constraints.

To perform the first step we need a Hilbert space of functionals ψ[q] of spatial metrics. Unfortunately, the space of metrics, or alternatively extrinsic curvature tensors, is mathematically poorly understood and not much is known about suitable inner products. At this point, a new set of variables introduced by Ashtekar [7, 8, 30] becomes essential. This is a triad formulation, but uses the triad in a densitized form (i.e., it is multiplied with an additional factor of a Jacobian under coordinate transformations). The densitized triad \(E_i^a\) is then related to the triad by \(E_i^a = \vert \det e_j^b{\vert ^{- 1}}e_i^a\) but has the same properties concerning gauge rotations and its orientation (note the absolute value which is often omitted). The densitized triad is conjugate to extrinsic curvature coefficients \(K_a^i: = {K_{ab}}e_i^b\):

$$\{K_a^i(x),E_j^b(y)\} = 8\pi G\delta _a^b\delta _j^i\delta (x,y)$$
(1)

with the gravitational constant G. Extrinsic curvature is then replaced by the Ashtekar connection

$$A_a^i = \Gamma _a^i + \gamma K_a^i$$
(2)

with a positive value for γ, the Barbero-Immirzi parameter [30, 133]. Classically, this number can be changed by a canonical transformation of the fields, but it will play a more important and fundamental role upon quantization. The Ashtekar connection is defined in such a way that it is conjugate to the triad,

$$\{A_a^i(x),E_j^b(y)\} = 8\pi \gamma G\delta _a^b\delta _j^i\delta (x,y)$$
(3)

and obtains its transformation properties as a connection from the spin connection

$$\Gamma _a^i = - {\epsilon ^{ijk}}e_j^b({\partial _{\left[ a \right.}}e_{\left. e \right]}^k + {1 \over 2}e_k^ce_a^l{\partial _{\left[ c \right.}}e_{\left. b \right]}^l).$$
(4)

Spatial geometry is then obtained directly from the densitized triad, which is related to the spatial metric by

$$E_i^aE_i^b = {q^{ab}}\det \,q.$$

There is more freedom in a triad since it can be rotated without changing the metric. The theory is independent of such rotations provided the Gauss constraint

$$G[\Lambda ] = {1 \over {8\pi \gamma G}}\int\nolimits_\Sigma {{{\rm{d}}^3}x{\Lambda ^i}{D_a}E_i^a} = {1 \over {8\pi \gamma G}}\int\nolimits_\Sigma {{{\rm{d}}^3}x{\Lambda ^i}({\partial _a}E_i^a + {\epsilon _{ijk}}A_a^jE_k^a) \approx 0}$$
(5)

is satisfied. Independence from any spatial coordinate system or background is implemented by the diffeomorphism constraint (modulo Gauss constraint)

$$D[{N^a}] = {1 \over {8\pi \gamma G}}\int\nolimits_\Sigma {{{\rm{d}}^3}x{N^a}F_{ab}^iE_i^b} \approx 0,$$
(6)

with the curvature \(F_{ab}^i\) of the Ashtekar connection. In this setting, one can then discuss spatial geometry and its quantization.

Space-time geometry, however, is more complicated to deduce since it requires a good knowledge of the dynamics. In a canonical setting, dynamics is implemented by the Hamiltonian constraint

$$H[N] = {1 \over {16\pi \gamma G}}\int\nolimits_\Sigma {{{\rm{d}}^3}xN\vert \det E{\vert ^{- 1/2}}\left({{\epsilon _{ijk}}F_{ab}^iE_j^aE_k^b - 2(1 + {\gamma ^2})K_{[a}^iK_{b]}^jE_i^aE_j^b} \right)} \approx 0,$$
(7)

where extrinsic curvature components have to be understood as functions of the Ashtekar connection and the densitized triad through the spin connection.

3.3 Representation

The key new aspect is now that we can choose the space of Ashtekar connections as our configuration space whose structure is much better understood than that of a space of metrics. Moreover, the formulation lends itself easily to a background independent quantization. To see this we need to remember that quantizing field theories requires one to smear fields, i.e., to integrate them over regions in order to obtain a well-defined algebra without δ-functions as in Equation (3). Usually this is done by integrating both configuration and momentum variables over three-dimensional regions, which requires an integration measure. This is no problem in ordinary field theories, which are formulated on a background such as Minkowski or a curved space. However, doing this here for gravity in terms of Ashtekar variables would immediately spoil any possible background independence since a background would already occur at this very basic step.

There is now a different smearing available that does not require a background metric. Instead of using three-dimensional regions we integrate the connection along one-dimensional curves e and exponentiate in a path-ordered manner, resulting in holonomies

$${h_e}(A) ={\mathcal P}\exp {\int\nolimits_e {{\tau _i}A_a^i\dot{e}} ^a}{\rm{d}}t$$
(8)

with tangent vector ėa to the curve e and \({\tau _j} = - {1 \over 2}i{\sigma _j}\) in terms of Pauli matrices. The path ordered exponentiation needs to be done in order to obtain a covariant object from the non-Abelian connection. The prevalence of holonomies or, in their most simple gauge invariant form as Wilson loops trhe(A) for closed e, is the origin of loop quantum gravity and its name [175]. Similarly, densitized vector fields can naturally be integrated over 2-dimensional surfaces, resulting in fluxes

$${F_S}(E) = \int\nolimits_S {{\tau ^i}E_{i}^{a}} {n_a}{{\rm{d}}^2}y$$
(9)

with the co-normal na to the surface.

The Poisson algebra of holonomies and fluxes is now well-defined and one can look for representations on a Hilbert space. We also require diffeomorphism invariance, i.e., there must be a unitary action of the diffeomorphism group on the representation by moving edges and surfaces in space. This is required since the diffeomorphism constraint has to be imposed later. Under this condition, there is even a unique representation that defines the kinematical Hilbert space [179, 180, 164, 183, 113, 146].

We can construct the Hilbert space in the representation where states are functionals of connections. This can easily be done by using holonomies as “creation operators” starting with a “ground state” which does not depend on connections at all. Multiplying with holonomies then generates states that do depend on connections, but only along the edges used in the process. These edges can be collected in a graph appearing as a label of the state. An independent set of states is given by spin network states [178] associated with graphs whose edges are labeled by irreducible representations of the gauge group SU(2), in which to evaluate the edge holonomies, and whose vertices are labeled by matrices specifying how holonomies leaving or entering the vertex are multiplied together. The inner product on this state space is such that these states, with an appropriate definition of independent contraction matrices in vertices, are orthonormal.

Spatial geometry can be obtained from fluxes representing the densitized triad. Since these are now momenta, they are represented by derivative operators with respect to values of connections on the flux surface. States as constructed above depend on the connection only along edges of graphs such that the flux operator is non-zero only if there are intersection points between its surface and the graph in the state it acts on [145]. Moreover, the contribution from each intersection point can be seen to be analogous to an angular momentum operator in quantum mechanics which has a discrete spectrum [20]. Thus, when acting on a given state we obtain a finite sum of discrete contributions and thus a discrete spectrum of flux operators. The spectrum depends on the value of the Barbero-Immirzi parameter, which can accordingly be fixed using implications of the spectrum such as black hole entropy, which gives a value of the order of but smaller than one [11, 12, 108, 155]. Moreover, since angular momentum operators do not commute, flux operators do not commute in general [17]. There is thus no triad representation, which is another reason why using a metric formulation and trying to build its quantization with functionals on a metric space is difficult.

There are important basic properties of this representation, which we will use later on. First, as already noted, flux operators have discrete spectra and, secondly, holonomies of connections are well-defined operators. It is, however, not possible to obtain operators for connection components or their integrations directly but only in the exponentiated form. These are direct consequences of the background independent quantization and translate to particular properties of more complicated operators.

3.4 Function spaces

A connection 1-form \(A_a^i\) can be reconstructed uniquely if all its holonomies are known [118]. It is thus sufficient to parameterize the configuration space by matrix elements of he for all edges in space. This defines an algebra of functions on the infinite dimensional space of connections \({\mathcal A}\), which are multiplied as \(\mathbb C\)-valued functions. Moreover, there is a duality operation by complex conjugation, and if the structure group G is compact a supremum norm exists since matrix elements of holonomies are then bounded. Thus, matrix elements form an Abelian C*-algebra with unit as a subalgebra of all continuous functions on \({\mathcal A}\).

Any Abelian C*-algebra with unit can be represented as the algebra of all continuous functions on a compact space \({\bar {\mathcal A}}\). The intuitive idea is that the original space \({\mathcal A}\), which has many more continuous functions, is enlarged by adding new points to it. This increases the number of continuity conditions and thus shrinks the set of continuous functions. This is done until only matrix elements of holonomies survive when continuity is imposed, and it follows from general results that the enlarged space must be compact for an Abelian unital C*-algebra. We thus obtain a compactification A, the space of generalized connections [23], which densely contains the space \({\mathcal A}\).

There is a natural diffeomorphism invariant measure dμAL on \({\bar {\mathcal A}}\), the Ashtekar-Lewandowski measure [19], which defines the Hilbert space \({\mathcal H} = {L^2}(\bar {\mathcal A},{\rm{d}}{\mu _{{\rm{AL}}}})\) of square integrable functions on the space of generalized connections. A dense subset Cyl of functions is given by cylindrical functions \(f({h_{{e_1}}}, \ldots, {h_{{e_n}}})\), which depend on the connection through a finite but arbitrary number of holonomies. They are associated with graphs g formed by the edges e1, …, en. For functions cylindrical with respect to two identical graphs the inner product can be written as

$$\langle f\vert g\rangle = \int\nolimits_{\bar{\mathcal A}} {{\rm{d}}{\mu _{AL}}(A)f(A)\ast g(A)} = \int\nolimits_{{\rm{SU}}{{(2)}^n}} {\prod\limits_{i = 1}^n {{\rm{d}}{\mu _{\rm{H}}}({h_i})f({h_1}, \ldots, {h_n})\ast g}} ({h_1}, \ldots, {h_n})$$
(10)

with the Haar measure on G. The importance of generalized connections can be seen from the fact that the space \({\mathcal A}\) of smooth connections is a subset of measure zero in \({\bar {\mathcal A}}\) [154].

With the dense subset Cyl of \({\mathcal H}\) we obtain the Gel’fand triple

$${\rm{Cyl}} \subset {\mathcal H} \subset {\rm{Cyl}}\ast$$
(11)

with the dual Cyl* of linear functionals from Cyl to the set of complex numbers. Elements of Cyl* are distributions, and there is no inner product on the full space. However, one can define inner products on certain subspaces defined by the physical context. Often, those subspaces appear when constraints with continuous spectra are solved following the Dirac procedure. Other examples include the definition of semiclassical or, as we will use in Section 6, symmetric states.

3.5 Composite operators

From the basic operators we can construct more complicated ones which, with growing degree of complexity, will be more and more ambiguous for instance from factor ordering choices. Quite simple expressions exist for the area and volume operator [177, 20, 21], which are constructed solely from fluxes. Thus, they are less ambiguous since no factor ordering issues with holonomies arise. This is true because the area of a surface and volume of a region can be written classically as functionals of the densitized triad alone, \({A_S} = \int\nolimits_S {\sqrt {E_i^a{n_a}E_i^b{n_b}} {{\rm{d}}^2}y}\) and \({V_R} = \int\nolimits_R {\sqrt {\vert\det \,E_i^a\vert}}{{\rm{d}}^3}x\). At the quantum level, this implies that, just as fluxes, also area and volume have discrete spectra showing that spatial quantum geometry is discrete. (For discrete approaches to quantum gravity in general see [150].) All area eigenvalues are known explicitly, but this is not possible even in principle for the volume operator. Nevertheless, some closed formulas and numerical techniques exist [l49, 103, 102, 83].

The length of a curve, on the other hand, requires the co-triad which is an inverse of the densitized triad and is more problematic. Since fluxes have discrete spectra containing zero, they do not have densely defined inverse operators. As we will describe below, it is possible to quantize those expressions but requires one to use holonomies. Thus, here we encounter more ambiguities from factor ordering. Still, one can show that also length operators have discrete spectra [192].

Inverse densitized triad components also arise when we try to quantize matter Hamiltonians such as

$${H_\phi} = \int {{{\rm{d}}^3}x\left({{1 \over 2}{{p_\phi ^2 + E_i^aE_i^b{\partial _a}\phi {\partial _b}\phi} \over {\sqrt {\vert \det E_j^c\vert}}} + \sqrt {\vert \det E_j^c\vert} V(\phi)} \right)}$$
(12)

for a scalar field Φ with momentum pΦ and potential V(Φ) (not to be confused with volume). The inverse determinant again cannot be quantized directly by using, e.g., an inverse of the volume operator which does not exist. This seems, at first, to be a severe problem not unlike the situation in quantum field theory on a background where matter Hamiltonians are divergent. Yet, it turns out that quantum geometry allows one to quantize these expressions in a well-defined manner [193].

To do this, we notice that the Poisson bracket of the volume with connection components,

$$\{A_a^i,\int {\sqrt {\vert \det E\vert} {d^3}x\} = 2\pi \gamma G{\epsilon}^{ijk}{\epsilon _{abc}}{{E_j^bE_k^c} \over {\sqrt {\vert \det E\vert}}},}$$
(13)

amounts to an inverse of densitized triad components and does allow a well-defined quantization: we can express the connection component through holonomies, use the volume operator and turn the Poisson bracket into a commutator. Since all operators involved have a dense intersection of their domains of definition, the resulting operator is densely defined and amounts to a quantization of inverse powers of the densitized triad.

This also shows that connection components or holonomies are required in this process, and thus ambiguities can arise even if initially one starts with an expression such as \({\sqrt {\vert\det E\vert} ^{- 1}}\), which only depends on the triad. There are also many different ways to rewrite expressions as above, which all are equivalent classically but result in different quantizations. In classical regimes this would not be relevant, but can have sizeable effects at small scales. In fact, this particular aspect, which as a general mechanism is a direct consequence of the background independent quantization with its discrete fluxes, implies characteristic modifications of the classical expressions on small scales. We will discuss this and more detailed examples in the cosmological context in Section 4.

3.6 Hamiltonian constraint

Similarly to matter Hamiltonians one can also quantize the Hamiltonian constraint in a well-defined manner [194]. Again, this requires to rewrite triad components and to make other regularization choices. Thus, there is not just one quantization but a class of different possibilities.

It is more direct to quantize the first part of the constraint containing only the Ashtekar curvature. (This part agrees with the constraint in Euclidean signature and Barbero-Immirzi parameter γ = 1, and so is sometimes called Euclidean part of the constraint.) Triad components and their inverse determinant are again expressed as a Poisson bracket using the identity (13), and curvature components are obtained through a holonomy around a small loop a of coordinate size Δ and with tangent vectors \(s_1^a\) and \(s_2^a\) at its base point [176]:

$$s_1^{a}s_2^{b}F_{ab}^i{\tau _i} = {\Delta}^{-1}({h_\alpha} - 1) - O(\Delta).$$
(14)

Putting this together, an expression for the Euclidean part HE[N] can then be constructed in the schematic form

$${H^E}[N] \propto \sum\limits_{v} N (v){\epsilon ^{IJK}}{\rm tr}({h_{\alpha IJ}}{h_{sK}}\{h_{{s_K}}^{- 1},V\}) + O(\Delta),$$
(15)

where one sums over all vertices of a triangulation of space whose tetrahedra are used to define closed curves αIJ and transversal edges sK.

An important property of this construction is that coordinate functions such as A disappear from the leading term, such that the coordinate size of the discretization is irrelevant. Nevertheless, there are several choices to be made, such as how a discretization is chosen in relation to a graph the constructed operator is supposed to act on, which in later steps will have to be constrained by studying properties of the quantization. Of particular interest is the holonomy hα since it creates new edges to a graph, or at least new spin on existing ones. Its precise behavior is expected to have a strong influence on the resulting dynamics [189]. In addition, there are factor ordering choices, i.e., whether triad components appear to the right or left of curvature components. It turns out that the expression above leads to a well-defined operator only in the first case, which in particular requires an operator non-symmetric in the kinematical inner product. Nevertheless, one can always take that operator and add its adjoint (which in this full setting does not simply amount to reversing the order of the curvature and triad expressions) to obtain a symmetric version, such that the choice still exists. Another choice is the representation chosen to take the trace, which for the construction is not required to be the fundamental one [116].

The second part of the constraint is more complicated since one has to use the function Γ(E) in \(K_a^i\). As also developed in [194], extrinsic curvature can be obtained through the already constructed Euclidean part via K ∼ {HE, V}. The result, however, is rather complicated, and in models one often uses a more direct way exploiting the fact that Γ has a more special form. In this way, additional commutators in the general construction can be avoided, which usually does not have strong effects. Sometimes, however, these additional commutators can be relevant, which can always be decided by a direct comparison of different constructions (see, e.g., [125]).

3.7 Open issues

For an anomaly-free quantization the constraint operators have to satisfy an algebra mimicking the classical one. There are arguments that this is the case for the quantization as described above when each loop a contains exactly one vertex of a given graph [191], but the issue is still open. Moreover, the operators are quite complicated and it is not easy to see if they have the correct expectation values in appropriately defined semiclassical states.

Even if one regards the quantization and semiclassical issues as satisfactory, one has to face several hurdles in evaluating the theory. There are interpretational issues of the wave function obtained as a solution to the constraints, and also the problem of time or observables emerges [143]. There is a wild mixture of conceptual and technical problems at different levels, not at least because the operators are quite complicated. For instance, as seen in the rewriting procedure above, the volume operator plays an important role even if one is not necessarily interested in the volume of regions. Since this operator is complicated, without an explicitly known spectrum, it translates to complicated matrix elements of the constraints and matter Hamiltonians. Loop quantum gravity should thus be considered as a framework rather than a uniquely defined theory, which however has important rigid aspects. This includes the basic representation of the holonomy-flux algebra and its general consequences.

All this should not come as a surprise since even classical gravity, at this level of generality, is complicated enough. Most solutions and results in general relativity are obtained with approximations or assumptions, one of the most widely used being symmetry reduction. In fact, this allows access to the most interesting gravitational phenomena such as cosmological expansion, black holes and gravitational waves. Similarly, symmetry reduction is expected to simplify many problems of full quantum gravity by resulting in simpler operators and by isolating conceptual problems such that not all of them need to be considered at once.

4 Loop Cosmology

Je abstrakter die Wahrheit ist, die du lehren willst, um so mehr mußt du noch die Sinne zu ihr verfuhren.

(The more abstract the truth you want to teach is, the more you have to seduce to it the senses.)

Friedrich Nietzsche

Beyond Good and Evil

The gravitational field equations, for instance in the case of cosmology where one can assume homogeneity and isotropy, involve components of curvature as well as the inverse metric. (Computational methods to derive information from these equations are described in [5].) Since singularities occur, these components will become large in certain regimes, but the equations have been tested only in small curvature regimes. On small length scales such as close to the Big Bang, modifications to the classical equations are not ruled out by observations and can be expected from candidates of quantum gravity. Quantum cosmology describes the evolution of a universe by a constraint equation for a wave function, but some effects can be included already at the level of effective classical equations. In loop quantum gravity, the main modification happens through inverse metric components which, e.g., appear in the kinematic term of matter Hamiltonians. This one modification is mainly responsible for all the diverse effects of loop cosmology.

4.1 Isotropy

Isotropy reduces the phase space of general relativity to be 2-dimensional since, up to SU(2)-gauge freedom, there is only one independent component in an isotropic connection and triad, respectively, which is not already determined by the symmetry. This is analogous to metric variables, where the scale factor a is the only free component in the spatial part of an isotropic metric

$${\rm{d}}{s^2} = - N{(t)^2}{\rm{d}}{t^2} + a{(t)^2}({(1 - k{r^2})^{- 1}}{\rm{d}}{r^2} + {r^2}{\rm{d}}{\Omega ^2}).$$
(16)

The lapse function N(t) does not play a dynamical role and correspondingly does not appear in the Friedmann equation

$${\left({{{\dot{a}} \over a}} \right)^2} + {k \over {{a^2}}} = {{8\pi G} \over 3}{a^{- 3}}{H_{{\rm{matter}}}}(a)$$
(17)

with the matter Hamiltonian Hmatter and the gravitational constant G, and the parameter k taking the discrete values zero or ±1 depending on the symmetry group or intrinsic spatial curvature.

Indeed, N(t) can simply be absorbed into the time coordinate by defining proper time r through dτ = N(t)dt. This is not possible for the scale factor since it depends on time but multiplies space differentials in the line element. The scale factor can only be rescaled by an arbitrary constant, which can be normalized at least in the closed model where k = 1.

One can understand these different roles of metric components also from a Hamiltonian analysis of the Einstein-Hilbert action

$${S_{{\rm{EH}}}} = {1 \over {16\pi G}}\int {{\rm{d}}t{{\rm{d}}^3}xR[g]}$$

specialized to isotropic metrics (16) whose Ricci scalar is

$$R = 6\left({{{\ddot{a}} \over {{N^2}a}} + {{{{\dot{a}}^2}} \over {{N^2}{a^2}}} + {k \over {{a^2}}} - {{\dot{a}} \over a}{{\dot{N}} \over {{N^3}}}} \right).$$

The action then becomes

$$S = {{{V_0}} \over {16\pi G}}\int {{\rm{d}}t\;N{a^3}R = {{3{V_0}} \over {8\pi G}}} \int {{\rm{d}}t\;N} \left({- {{a{{\dot a}^2}} \over {{N^2}}}ka} \right)$$

(with the spatial coordinate volume \({V_0} = \int\nolimits_\Sigma{{\rm{d}}^3}x\)) after integrating by parts, from which one derives the momenta

$${p_a} = {{\partial L} \over {\partial \dot a}} = - {{3{V_0}} \over {4\pi G}}{{a\dot a} \over N},\quad {p_N} = {{\partial L} \over {\partial \dot N}} = 0$$

illustrating the different roles of a and N. Since pN must vanish, N is not a degree of freedom but a Lagrange multiplier. It appears in the canonical action S = (16πG)−1 ∫ dtpaNH)) only as a factor of

$$H = - {{2\pi G} \over 3}{{p_a^2} \over {{V_0}a}} - {3 \over {8\pi G}}{V_0}ak$$

such that variation with respect to N forces H, the Hamiltonian constraint, to be zero. In the presence of matter, H also contains the matter Hamiltonian, and its vanishing is equivalent to the Friedmann equation.

4.2 Isotropy: Connection variables

Isotropic connections and triads, as discussed in Appendix B.2, are analogously described by single components \({\tilde c}\) and \({\tilde p}\), respectively, related to the scale factor by

$$\vert \tilde{p}\vert = {\tilde{a}^2} = {{{a^2}} \over 4}$$
(18)

for the densitized triad component \({\tilde p}\) and

$$\tilde{c} = \tilde\Gamma + \gamma \dot{\tilde{a}} = {1 \over 2}(k + \gamma \dot a)$$
(19)

for the connection component \({\tilde c}\). Both components are canonically conjugate:

$$\{\tilde{c},\tilde{p}\} = {{8\pi \gamma G} \over 3}{V_0}.$$
(20)

It is convenient to absorb factors of V0 into the basic variables, which is also suggested by the integrations in holonomies and fluxes on which background independent quantizations are built [15]. We thus define

$$p = V_0^{2/3}\tilde{p},\quad c = V_0^{1/3}\tilde{c}$$
(21)

together with \(\Gamma = V_0^{1/3}\tilde \Gamma\). The symplectic structure is then independent of V0 and so are integrated densities such as total Hamiltonians. For the Hamiltonian constraint in isotropic Ashtekar variables we have

$$H = - {3 \over {8\pi G}}({\gamma ^{- 2}}{(c - \Gamma)^2} + {\Gamma ^2})\sqrt {\vert p\vert} + {H_{{\rm{matter}}}}(p) = 0,$$
(22)

which is exactly the Friedmann equation. (In most earlier papers on loop quantum cosmology some factors in the basic variables and classical equations are incorrect due, in part, to the existence of different and often confusing notations in the loop quantum gravity literature.Footnote 1)

The part of phase space where we have p = 0 and thus a = 0 plays a special role since this is where isotropic classical singularities are located. On this subset the evolution equation (17) with standard matter choices is singular in the sense that Hmatter, e.g.,

$${H_\phi}(a,\phi, {p_\phi}) = {1 \over 2}\vert p{\vert ^{- 3/2}}p_\phi ^2 + \vert p{\vert ^{3/2}}V(\phi)$$
(23)

for a scalar Φ with momentum and potential V(Φ), diverges and the differential equation does not pose a well-defined initial value problem there. Thus, once such a point is reached the further evolution is no longer determined by the theory. Since, according to singularity theorems [123, 80], any classical trajectory must intersect the subset a = 0 for the matter we need in our universe, the classical theory is incomplete.

This situation, certainly, is not changed by introducing triad variables instead of metric variables. However, the situation is already different since p = 0 is a submanifold in the classical phase space of triad variables where p can have both signs (the sign determining whether the triad is left or right handed, i.e., the orientation). This is in contrast to metric variables where a = 0 is a boundary of the classical phase space. There are no implications in the classical theory since trajectories end there nonetheless, but it will have important ramifications in the quantum theory (see the sections following Section 5.13).

4.3 Isotropy: Implications of a loop quantization

We are now dealing with a simple system with finitely many degrees of freedom, subject to a constraint. It is well known how to quantize such a system from quantum mechanics, which has been applied to cosmology starting with DeWitt [104]. Here, one chooses a metric representation for wave functions, i.e., ψ(a), on which the scale factor acts as multiplication operator and its conjugate pa, related to ȧ, as a derivative operator. These basic operators are then used to form the Wheeler-DeWitt operator quantizing the constraint (17) once a factor ordering is chosen.

This prescription is rooted in quantum mechanics which, despite its formal similarity, is physically very different from cosmology. The procedure looks innocent, but one should realize that there are already basic choices involved. Choosing the factor ordering is harmless, even though results can depend on it [142]. More importantly, one has chosen the Schrödinger representation of the classical Poisson algebra, which immediately implies the familiar properties of operators such as the scale factor with a continuous spectrum. There are inequivalent representations with different properties, and it is not clear that this representation, which works well in quantum mechanics, is also correct for quantum cosmology. In fact, quantum mechanics is not very sensitive to the representation chosen [18] and one can use the most convenient one. This is the case because energies and thus oscillation lengths of wave functions described usually by quantum mechanics span only a limited range. Results can then be reproduced to arbitrary accuracy in any representation. Quantum cosmology, in contrast, has to deal with potentially infinitely high matter energies, leading to small oscillation lengths of wave functions, such that the issue of quantum representations becomes essential.

That the Wheeler-DeWitt representation may not be the right choice is also indicated by the fact that its scale factor operator has a continuous spectrum, while quantum geometry which is a well-defined quantization of the full theory, implies discrete volume spectra. Indeed, the Wheeler-DeWitt quantization of full gravity exists only formally, and its application to quantum cosmology simply quantizes the classically reduced isotropic system. This is much easier, and also more ambiguous, and leaves open many consistency considerations. It would be more reliable to start with the full quantization and introduce the symmetries there, or at least follow the same constructions of the full theory in a reduced model. If this is done, it turns out that indeed we obtain a quantum representation inequivalent to the Wheeler-DeWitt representation, with strong implications in high energy regimes. In particular, just as the full theory such a quantization has a volume or operator with a discrete spectrum, as derived in Section 5.2.

4.4 Isotropy: Effective densities and equations

The isotropic model is thus quantized in such a way that the operator \({\hat p}\) has a discrete spectrum containing zero. This immediately leads to a problem since we need a quantization of |p|−3/2 in order to quantize a matter Hamiltonian such as (23) where not only the matter fields but also geometry are quantized. However, an operator with zero in the discrete part of its spectrum does not have a densely defined inverse and does not allow a direct quantization of |p|−3/2.

This leads us to the first main effect of the loop quantization: It turns out that despite the non-existence of an inverse operator of \({\hat p}\) one can quantize the classical |p|−3/2 to a well-defined operator. This is not just possible in the model but also in the full theory where it even has been defined first [193]. Classically, one can always write expressions in many equivalent ways, which usually result in different quantizations. In the case of |p|−3/2, as discussed in Section 5.3, there is a general class of ways to rewrite it in a quantizable manner [41] which differ in details but have all the same important properties. This can be parameterized by a function d(p)j,l [47, 50] which replaces the classical |p|3/2 and strongly deviates from it for small p while being very close at large p. The parameters \(j \in {1 \over 2}\mathbb N\) and 0 < l < 1 specify quantization ambiguities resulting from different ways of rewriting. With the function

$$\begin{array}{*{20}c}{{p_l}(q) = {3 \over {2l}}{q^{1 - l}}{{\left({{1 \over {l + 2}}({{(q + 1)}^{l + 2}} - \vert q - 1\vert} \right.}^{l + 1}})} \\ {\left. {\quad \quad \quad - {1 \over {l + 2}}q({{(q + 1)}^{l + 1}} - {\mathop{\rm sgn}} (q - 1)\vert q - 1{\vert ^{l + 1}}} \right)} \\ \end{array}$$
(24)

we have

$$d{(p)_{j,l}}: = \vert p{\vert ^{- 3/2}}pl{(3\vert p\vert/\gamma j\ell _{\rm P}^2)^{3/(2 - 2l)}},$$
(25)

which indeed fulfills d(p)j,l ∼ |p|3/2 for \(\vert p\vert \gg {p_ast}: = {1 \over 3}j\gamma \ell _{\rm{P}}^2\), but is finite with a peak around p* and approaches zero at p = 0 in a manner

$$d{(p)_{j,l}} \sim {3^{3(3 - l)/(2 - 2l)}}{(l + 1)^{ - 3/(2 - 2l)}}{(\gamma j)^{ - 3(2 - l)/(2 - 2l)}}\ell _P^{ - 3(2 - l)/(1 - l)}|p{|^{3/(2 - 2l)}}$$
(26)

as it follows from pl(q) ∼ 3q2−l/(1 + l). Some examples displaying characteristic properties are shown in Figure 9 in Section 5.3

The matter Hamiltonian obtained in this manner will thus behave differently at small p. At those scales also other quantum effects such as fluctuations can be important, but it is possible to isolate the effect implied by the modified density (25). We just need to choose a rather large value for the ambiguity parameter j such that modifications become noticeable already in semiclassical regimes. This is mainly a technical tool to study the behavior of equations, but can also be used to find constraints on the allowed values of ambiguity parameters.

We can thus use classical equations of motion, which are corrected for quantum effects by using the effective matter Hamiltonian

$$H_\phi ^{({\rm{eff}})}(p,\phi, {p_\phi}):={1 \over 2}d{(p)_{j,l}}p_\phi ^2 + \vert p{\vert ^{3/2}}V(\phi)$$
(27)

(see Section 5.5 for details on effective equations). This matter Hamiltonian changes the classical constraint such that now

$$H = - {3 \over {8\pi G}}({\gamma ^{- 2}}{(c - \Gamma)^2} + {\Gamma ^2})\sqrt {\vert p\vert} + H_\phi ^{({\rm{eff}})}(p,\phi, {p_\phi}) = 0.$$
(28)

Since the constraint determines all equations of motion, they also change: We obtain the effective Friedmann equation from H = 0,

$${\left({{{\dot a} \over a}} \right)^2} + {k \over {{a^2}}} = {{8\pi G} \over 3}\left({{1 \over 2}\vert p{\vert ^{- 3/2}}d{{(p)}_{j,l}}p_\phi ^2 + V(\phi)} \right)$$
(29)

and the effective Raychaudhuri equation from ċ = {c, H},

$${{\ddot a} \over a} = - {{4\pi G} \over {3\vert p{\vert ^{3/2}}}}\left({{H_{{\rm{matter}}}}(p,\phi, {p_\phi}) - 2p{{\partial {H_{{\rm{matter}}}}(p,\phi, {p_\phi})} \over {\partial p}}} \right)$$
(30)
$$= - {{8\pi G} \over 3}\left({\vert p{\vert ^{- 3/2}}d(p)_{j,l}^{- 1}{{\dot \phi}^2}\left({1 - {1 \over 4}a{{d\log (\vert p{\vert ^{3/2}}d{{(p)}_{j,l}})} \over {da}}} \right) - V(\phi)} \right).$$
(31)

Matter equations of motion follow similarly as

$$\begin{array}{*{20}c}{\dot \phi = \{\phi, H\} = d{{(p)}_{j,l}}{p_\phi}} \\ {{{\dot p}_\phi} = \{{p_\phi},H\} = - \vert p{\vert ^{3/2}}V^{\prime}(\phi),} \\\end{array}$$

which can be combined to the effective Klein-Gordon equation

$$\ddot{\phi} = \dot{\phi} \dot{a}{{{\rm{d}}\log d{{(p)}_{j,l}}} \over {{\rm{d}}a}} - \vert p{\vert ^{3/2}}d{(p)_{j,l}}V^{\prime}(\phi).$$
(32)

Further discussion for different forms of matter can be found in [186].

4.5 Isotropy: Properties and intuitive meaning

As a consequence of the function d(p))j,l, the effective equations have different qualitative behavior at small versus large scales p. In the effective Friedmann equation (29), this is most easily seen by comparing it with a mechanics problem with a standard Hamiltonian, or energy, of the form

$$E = {1 \over 2}{\dot a^2} - {{2\pi G} \over {3{V_0}}}{a^{- 1}}d{(p)_{j,l}}p_\phi ^2 - {{4\pi G} \over 3}{a^2}V(\phi) = 0$$

restricted to be zero. If we assume a constant scalar potential V(Φ), there is no Φ-dependence and the scalar equations of motion show that is constant. Thus, the potential for the motion of a is essentially determined by the function d(p)j,l.

In the classical case, d(p) = |p|3/2 and the potential is negative and increasing, with a divergence at p = 0. The scale factor a is thus driven toward a = 0, which it will always reach in finite time where the system breaks down. With the effective density d(p)j,l, however, the potential is bounded from below, and is decreasing from zero for a = 0 to the minimum around p*. Thus, the scale factor is now slowed down before it reaches a = 0, which depending on the matter content could avoid the classical singularity altogether.

The behavior of matter is also different as shown by the effective Klein-Gordon equation (32). Most importantly, the derivative in the \({\dot \phi}\)-term changes sign at small a since the effective density is increasing there. Thus, the qualitative behavior of all the equations changes at small scales, which as we will see gives rise to many characteristic effects. Nevertheless, for the analysis of the equations as well as conceptual considerations it is interesting that solutions at small and large scales are connected by a duality transformation [147], which even exists between effective solutions for loop cosmology and braneworld cosmology [90].

We have seen that the equations of motion following from an effective Hamiltonian are expected to display qualitatively different behavior at small scales. Before discussing specific models in detail, it is helpful to observe what physical meaning the resulting modifications have.

Classical gravity is always attractive, which implies that there is nothing to prevent collapse in black holes or the whole universe. In the Friedmann equation this is expressed by the fact that the potential as used before is always decreasing toward a = 0 where it diverges. With the effective density, on the other hand, we have seen that the decrease stops and instead the potential starts to increase at a certain scale before it reaches zero at a = 0. This means that at small scales, where quantum gravity becomes important, the gravitational attraction turns into repulsion. In contrast to classical gravity, thus, quantum gravity has a repulsive component which can potentially prevent collapse. So far, this has only been demonstrated in homogeneous models, but it relies on a general mechanism which is also present in the full theory.

Not only the attractive nature of gravity changes at small scales, but also the behavior of matter in a gravitational background. Classically, matter fields in an expanding universe are slowed down by a friction term in the Klein-Gordon equation (32) where ȧd log a3/da = −3ȧ/a is negative. Conversely, in a contracting universe matter fields are excited and even diverge when the classical singularity is reached. This behavior turns around at small scales where the derivative dlogd(a)j,l/da becomes positive. Friction in an expanding universe then turns into antifriction such that matter fields are driven away from their potential minima before classical behavior sets in. In a contracting universe, on the other hand, matter fields are not excited by antifriction but freeze once the universe becomes small enough.

These effects do not only have implications for the avoidance of singularities at a = 0 but also for the behavior at small but non-zero scales. Gravitational repulsion can not only prevent collapse of a contracting universe [187] but also, in an expanding universe, enhance its expansion. The universe then accelerates in an inflationary manner from quantum gravity effects alone [45]. Similarly, the modified behavior of matter fields has implications for inflationary models [77].

4.6 Isotropy: Applications

There is now one characteristic modification in the matter Hamiltonian, coming directly from a loop quantization. Its implications can be interpreted as repulsive behavior on small scales and the exchange of friction and antifriction for matter, and it leads to many further consequences.

4.6.1 Collapsing phase

When the universe has collapsed to a sufficiently small size, repulsion becomes noticeable and bouncing solutions become possible as illustrated in Figure 1. Requirements for a bounce are that the conditions ȧ = 0 and ä > 0 can be fulfilled at the same time, where the first one can be evaluated with the Friedmann equation, and the second one with the Raychaudhuri equation. The first condition can only be fulfilled if there is a negative contribution to the matter energy, which can come from a positive curvature term k = 1 or a negative matter potential V(Φ) < 0. In those cases, there are classical solutions with ȧ = 0, but they generically have ä < 0 corresponding to a recollapse. This can easily be seen in the flat case with a negative potential where (30) is strictly negative with dloga3d(a)j,l/da ≈ 0 at large scales.

Figure 1
figure 1

Examples for bouncing solutions with positive curvature (left) or a negative potential (right, negative cosmological constant). The solid lines show solutions of effective equations with a bounce, while the dashed lines show classical solutions running into the singularity at a = 0 where Φ diverges.

The repulsive nature at small scales now implies a second point where ȧ = 0 from (29) at smaller a since the matter energy now decreases also for a → 0. Moreover, the modified Raychaudhuri equation (30) has an additional positive term at small scales such that ä > 0 becomes possible.

Matter also behaves differently through the modified Klein-Gordon equation (32). Classically, with ȧ < 0 the scalar experiences antifriction and Φ diverges close to the classical singularity. With the modification, antifriction turns into friction at small scales, damping the motion of Φ such that it remains finite. In the case of a negative potential [68] this allows the kinetic term to cancel the potential term in the Friedmann equation. With a positive potential and positive curvature, on the other hand, the scalar is frozen and the potential is canceled by the curvature term. Since the scalar is almost constant, the behavior around the turning point is similar to a de Sitter bounce [187, 203]. Further, more generic possibilities for bounces arise from other correction terms [100, 97].

4.6.2 Expansion

Repulsion can not only prevent collapse but also accelerates an expanding phase. Indeed, using the behavior (26) at small scales in the effective Raychaudhuri equation (30) shows that ä is generically positive since the inner bracket is smaller than −1/2 for the allowed values 0 < l < 1. Thus, as illustrated by the numerical solution in the upper left panel of Figure 2, inflation is realized by quantum gravity effects for any matter field irrespective of its form, potential or initial values [45]. The kind of expansion at early stages is generically super-inflationary, i.e., with equation of state parameter w < −1. For free massless matter fields, w usually starts very small, depending on the value of l, but with a non-zero potential such as a mass term for matter inflation w is generically close to exponential: weff ≈ −1. This can be shown by a simple and elegant argument independently of the precise matter dynamics [101]: The equation of state parameter is defined as w = P/ρ where P = −∂E/∂V is the pressure, i.e., the negative change of energy with respect to volume, and ρ = E/V energy density. Using the matter Hamiltonian for E and V = |p|3/2, we obtain

$${P_{{\rm{eff}}}} = - {1 \over 3}\vert p{\vert ^{- 1/2}}d^{\prime}(p)p_\phi ^2 - V(\phi)$$

and thus in the classical case

$$w = {{{1 \over 2}\vert p{\vert ^{- 3}}p_\phi ^2 - V(\phi)} \over {{1 \over 2}\vert p{\vert ^{- 3}}p_\phi ^2 + V(\phi)}}$$

as usually. In the modified case, however, we have

$${w_{{\rm{eff}}}} = {{{1 \over 3}\vert p\vert ^{{- 1/2}}d^{\prime}(p)p_\phi ^2 + V(\phi)} \over {{1 \over 2}\vert p\vert ^{{- 3/2}}d(p)p_\phi ^2 + V(\phi)}}.$$
Figure 2
figure 2

Example for a solution of a(t) and Φ(t), showing early loop inflation and later slow-roll inflation driven by a scalar that is pushed up its potential by loop effects. The left hand side is stretched in time so as to show all details. An idea of the duration of different phases can be obtained from Figure 3 .

In general, we need to know the matter behavior to know w and weff. But we can get generic qualitative information by treating pΦ and V(Φ) as unknowns determined by w and weff. In the generic case, there is no unique solution for \(p_\phi ^2\) and V(Φ) since, after all, pΦ and Φ change with t. They are now subject to two linear equations in terms of w and weff, whose determinant must be zero resulting in

$${w_{{\rm{eff}}}} = - 1 + {{\vert p{\vert ^{3/2}}(w + 1)(d(p) - {2 \over 3}\vert p\vert d^{\prime}(p))} \over {1 - w + (w + 1)\vert p{\vert ^{3/2}}d(p)}}.$$

Since for small p the numerator in the fraction approaches zero faster than the second part of the denominator, weff approaches minus one at small volume except for the special case w = 1, which is realized for V(Φ) = 0. Note that the argument does not apply to the case of vanishing potential since then \(p_\phi ^2 = {\rm{const}}\) and ) = 0 presents a unique solution to the linear equations for w and weff. In fact, this case leads in general to a much smaller \({w_{{\rm{eff}}}} = - {2 \over 3}\vert p\vert d(p)^\prime/d(p) \approx - 1/(1 - l) < - 1\) [45].

One can also see from the above formula that weff, though close to minus one, is a little smaller than minus one generically. This is in contrast to single field inflaton models where the equation of state parameter is a little larger than minus one. As we will discuss in Section 4.15, this opens the door to characteristic signatures distinguishing different models.

Again, also the matter behavior changes, now with classical friction being replaced by antifriction [77]. Matter fields thus move away from their minima and become excited even if they start close to a minimum (Figure 2). Since this does not only apply to the homogeneous mode, it can provide a mechanism of structure formation as discussed in Section 4.15. But also in combination with chaotic inflation as the mechanism to generate structure does the modified matter behavior lead to improvements: If we now view the scalar Φ as an inflaton field, it will be driven to large values in order to start a second phase of slow-roll inflation which is long enough. This is satisfied for a large range of the ambiguity parameters j and l [67] and can even leave signatures [197] in the cosmic microwave spectrum [134]: The earliest moments when the inflaton starts to roll down its potential are not slow roll, as can also be seen in Figures 2 and 3 where the initial decrease is steeper. Provided the resulting structure can be seen today, i.e., there are not too many e-foldings from the second phase, this can lead to visible effects such as a suppression of power. Whether or not those effects are to be expected, i.e., which magnitude of the inflaton is generically reached by the mechanism generating initial conditions, is currently being investigated at the basic level of loop quantum cosmology [27]. They should be regarded as first suggestions, indicating the potential of quantum cosmological phenomenology, which have to be substantiated by detailed calculations including inhomogeneities or at least anisotropic geometries. In particular the suppression of power can be obtained by a multitude of other mechanisms.

Figure 3
figure 3

mpg-Movie (1.12 MB) Still from a Movie showing the initial push of a scalar Φ up its potential and the ensuing slow-roll phase together with the corresponding inflationary phase of a. (For video see appendix)

4.6.3 Model building

It is already clear that there are different inflationary scenarios using effects from loop cosmology. A scenario without inflaton is more attractive since it requires less choices and provides a fundamental explanation of inflation directly from quantum gravity. However, it is also more difficult to analyze structure formation in this context while there are already well-developed techniques in slow role scenarios.

In these cases where one couples loop cosmology to an inflaton model one still requires the same conditions for the potential, but generically gets the required large initial values for the scalar by antifriction. On the other hand, finer details of the results now depend on the ambiguity parameters, which describe aspects of the quantization that also arise in the full theory.

It is also possible to combine collapsing and expanding phases in cyclic or oscillatory models [148]. One then has a history of many cycles separated by bounces, whose duration depends on details of the model such as the potential. There can then be many brief cycles until eventually, if the potential is right, one obtains an inflationary phase if the scalar has grown high enough. In this way, one can develop ideas for the pre-history of our universe before the Big Bang. There are also possibilities to use a bounce to describe the structure in the universe. So far, this has only been described in effective models [137] using brane scenarios [151] where the classical singularity has been assumed to be absent by yet to be determined quantum effects. As it turns out, the explicit mechanism removing singularities in loop cosmology is not compatible with the assumptions made in those effective pictures. In particular, the scalar was supposed to turn around during the bounce, which is impossible in loop scenarios unless it encounters a range of positive potential during its evolution [68]. Then, however, generically an inflationary phase commences as in [148], which is then the relevant regime for structure formation. This shows how model building in loop cosmology can distinguish scenarios that are more likely to occur from quantum gravity effects.

Cyclic models can be argued to shift the initial moment of a universe in the infinite past, but they do not explain how the universe started. An attempt to explain this is the emergent universe model [110, 112] where one starts close to a static solution. This is difficult to achieve classically, however, since the available fixed points of the equations of motion are not stable and thus a universe departs too rapidly. Loop cosmology, on the other hand, implies an additional fixed point of the effective equations which is stable and allows to start the universe in an initial phase of oscillations before an inflationary phase is entered [160, 53]. This presents a natural realization of the scenario where the initial scale factor at the fixed point is automatically small so as to start the universe close to the Planck phase.

4.6.4 Stability

Cosmological equations displaying super-inflation or antifriction are often unstable in the sense that matter can propagate faster than light. This has been voiced as a potential danger for loop cosmology, too [94, 95]. An analysis requires inhomogeneous techniques at least at an effective level, such as those described in Section 4.12. It has been shown that loop cosmology is free of this problem, because the modified behavior for the homogeneous mode of the metric and matter is not relevant for matter propagation [129]. The whole cosmological picture that follows from the effective equations is thus consistent.

4.7 Anisotropies

Anisotropic models provide a first generalization of isotropic ones to more realistic situations. They thus can be used to study the robustness of effects analyzed in isotropic situations and, at the same time, provide a large class of interesting applications. An analysis in particular of the singularity issue is important since the classical approach to a singularity can be very different from the isotropic one. On the other hand, the anisotropic approach is deemed to be characteristic even for general inhomogeneous singularities if the BKL scenario [31] is correct.

A general homogeneous but anisotropic metric is of the form

$${\rm{d}}{s^2} = N{(t)^2}{\rm{d}}{t^2} + \sum\limits_{I,J = 1}^3 {{q_{IJ}}(t){\omega ^I} \otimes \;} {\omega ^J}$$

with left-invariant 1-forms ωI on space Σ, which, thanks to homogeneity, can be identified with the simply transitive symmetry group S as a manifold. The left-invariant 1-forms satisfy the Maurer-Cartan relations

$${\rm{d}}{\omega ^I} = - {1 \over 2}C_{JK}^I{\omega ^J} \wedge {\omega ^K}$$

with the structure constants \(C_{JK}^I\) of the symmetry group. In a matrix parameterization of the symmetry group, one can derive explicit expressions for ωI from the Maurer-Cartan form ωITI = θMC = g1dg with generators TI of S.

The simplest case of a symmetry group is an Abelian one with \(C_{JK}^I=0\), corresponding to the Bianchi I model. In this case, S is given by ℝ3 or a torus, and left-invariant 1-forms are simply ωI = dxI in Cartesian coordinates. Other groups must be restricted to class A models in this context, satisfying \(C_{JI}^I=0\) since otherwise there is no Hamiltonian formulation. The structure constants can then be parameterized as \(C_{JK}^I = \epsilon _{JK}^I{n^{(I)}}\).

A common simplification is to assume the metric to be diagonal at all times, which corresponds to a reduction technically similar to a symmetry reduction. This amounts to \({q_{IJ}} = a_{(I)}^2{\delta _{IJ}}\) as well as KIJ = K(I)δIJ for the extrinsic curvature with KI = ȧI. Depending on the structure constants, there is also non-zero intrinsic curvature quantified by the spin connection components

$${\Gamma _I} = {1 \over 2}\left({{{{a_J}} \over {{a_K}}}{n^J} + {{{a_K}} \over {{a_J}}}{n^K} - {{a_I^2} \over {{a_J}{a_K}}}{n^I}} \right)\quad {\rm{for}}\quad {\epsilon _{IJK}} = 1.$$
(33)

This influences the evolution as follows from the Hamiltonian constraint

$$\begin{array}{*{20}c}{- {1 \over {8\pi G}}({a_1}{{\dot{a}}_2}{{\dot{a}}_3} + {a_2}{{\dot{a}}_1}{{\dot{a}}_2} + {a_3}{{\dot{a}}_1}{{\dot{a}}_2} - ({\Gamma _2}{\Gamma _3} - {n^1}{\Gamma _1}){a_1} - ({\Gamma _1}{\Gamma _3} - {n^2}{\Gamma _2}){a_2}} \\ {- ({\Gamma _1}{\Gamma _2} - {n^3}{\Gamma _3}){a_3}) + {H_{{\rm{matter}}}}({a_I}) = 0.\quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad} \\\end{array}$$
(34)

In the vacuum Bianchi I case the resulting equations are easy to solve by \({a_I} \propto {t^{\alpha I}}\) with \({\sum\nolimits_I}{\alpha _I} = {\sum\nolimits_I}\alpha _I^2 = 1\) [135]. The volume a1a2a3t vanishes for t = 0 where the classical singularity appears. Since one of the exponents αI must be negative, however, only two of the aI vanish at the classical singularity while the third one diverges. This already demonstrates how different the behavior can be from the isotropic one and that anisotropic models provide a crucial test of any mechanism for singularity resolution.

4.8 Anisotropy: Connection variables

A densitized triad corresponding to a diagonal homogeneous metric has real components pI with |pI| = aJaK if ϵIJK = 1 [48]. Connection components are cI = ΓI + γKI = ΓI + γȧI and are conjugate to the \({p_I},\{{c_I},{p^J}\} = 8\pi \gamma G\delta _I^J\). In terms of triad variables we now have spin connection components

$${\Gamma _I} = {1 \over 2}\left({{{{p^K}} \over {{p^J}}}{n^J} + {{{p^J}} \over {{p^K}}}{n^K} - {{{p^J}{p^K}} \over {{{({p^I})}^2}}}{n^I}} \right)$$
(35)

and the Hamiltonian constraint (in the absence of matter)

$$\begin{array}{*{20}c}{H = {1 \over {8\pi G}}\left\{{[({c_2}{\Gamma _3} + {c_3}{\Gamma _2} - {\Gamma _2}{\Gamma _3})(1 + {\gamma ^{- 2}}) - {n^1}{c_1} - {\gamma ^{- 2}}{c_2}{c_3}]\sqrt {\left\vert {{{{p^2}{p^3}} \over {{p^1}}}} \right\vert}} \right.} \\ {\quad \quad + [({c_1}{\Gamma _3} + {c_3}{\Gamma _1} - {\Gamma _1}{\Gamma _3})(1 + {\gamma ^{- 2}}) - {n^2}{c_2} - {\gamma ^{- 2}}{c_1}{c_3}]\sqrt {\left\vert {{{{p^1}{p^3}} \over {{p^2}}}} \right\vert}} \\ {\left. {\quad \quad + [({c_1}{\Gamma _2} + {c_2}{\Gamma _1} - {\Gamma _1}{\Gamma _2})(1 + {\gamma ^{- 2}}) - {n^3}{c_3} - {\gamma ^{- 2}}{c_1}{c_2}]\sqrt {\left\vert {{{{p^1}{p^2}} \over {{p^3}}}} \right\vert}} \right\}.} \\ \end{array}$$
(36)

Unlike in isotropic models, we now have inverse powers of pI even in the vacuum case through the spin connection, unless we are in the Bianchi I model. This is a consequence of the fact that not just extrinsic curvature, which in the isotropic case is related to the matter Hamiltonian through the Friedmann equation, leads to divergences but also intrinsic curvature. These divergences are cut off by quantum geometry effects as before such that also the dynamical behavior changes. This can again be dealt with by effective equations where inverse powers of triad components are replaced by bounded functions [62]. However, even with those modifications, expressions for curvature are not necessarily bounded unlike in the isotropic case. This comes from the presence of different classical scales, aI, such that more complicated expressions as in ΓI are possible, while in the isotropic model there is only one scale and curvature can only be an inverse power of p, which is then regulated by effective expressions like d(p).

4.9 Anisotropy: Applications

4.9.1 Isotropization

Matter fields are not the only contributions to the Hamiltonian in cosmology, but also the effect of anisotropies can be included in this way to an isotropic model. The late time behavior of this contribution can be shown to behave as a−6 in the shear energy density [156], which falls off faster than any other matter component. Thus, toward later times the universe becomes more and more isotropic.

In the backward direction, on the other hand, this means that the shear term diverges most strongly, which suggests that this term should be most relevant for the singularity issue. Even if matter densities are cut off as discussed before, the presence of bounces would depend on the fate of the anisotropy term. This simple reasoning is not true, however, since the behavior of shear is only effective and uses assumptions about the behavior of matter. It can thus not simply be extrapolated to early times. Anisotropies are independent degrees of freedom which affect the evolution of the scale factor. But only in certain regimes can this contribution be modeled simply by a function of the scale factor alone; in general one has to use the coupled system of equations for the scale factor, anisotropies and possible matter fields.

4.9.2 Bianchi IX

Modifications to classical behavior are most drastic in the Bianchi IX model with symmetry group S ≅ SU(2) such that nI = 1. The classical evolution can be described by a 3-dimensional mechanics system with a potential obtained from (34) such that the kinetic term is quadratic in derivatives of aI with respect to a time coordinate τ defined by dt = a1a2a3dτ. This potential

$$\begin{array}{*{20}c}{W({p^I}) = ({\Gamma _2}{\Gamma _3} - {n^1}{\Gamma _1}){p^2}{p^3} + ({\Gamma _1}{\Gamma _3} - {n^2}{\Gamma _2}){p^1}{p^3} + ({\Gamma _1}{\Gamma _2} - {n^3}{\Gamma _3}){p^1}{p^2}} \\ {\quad \quad \quad = {1 \over 4}\left({{{\left({{{{p^2}{p^3}} \over {{p^1}}}} \right)}^2} + {{\left({{{{p^1}{p^3}} \over {{p^2}}}} \right)}^2} + {{\left({{{{p^1}{p^2}} \over {{p^3}}}} \right)}^2} - 2{{({p^1})}^2} - 2{{({p^2})}^2} - 2{{({p^3})}^2}} \right)} \\ \end{array}$$
(37)

diverges at small pI, in particular (in a direction dependent manner) at the classical singularity where all pI = 0. Figure 4 illustrates the walls of the potential, which with decreasing volume push the universe toward the classical singularity.

Figure 4
figure 4

mpg-Movie (1.66 MB) Still from a Movie illustrating the Bianchi IX potential ( 37 ) and the movement of its walls, rising toward zero p1 and p2 and along the diagonal direction, toward the classical singularity with decreasing volume \(V = \sqrt {\vert {p^1}{p^2}{p^3}\vert}\). The contours are plotted for the function W(p1, p2, V2/(p1p2)). (For video see appendix)

As before in isotropic models, effective equations where the behavior of eigenvalues of the spin connection components is used do not have this divergent potential. Instead, if two pj are held fixed and the third approaches zero, the effective quantum potential is cut off and goes back to zero at small values, which changes the approach to the classical singularity. Yet, the effective potential is unbounded if one pI diverges while another one goes to zero and the situation is qualitatively different from the isotropic case. Since the effective potential corresponds to spatial intrinsic curvature, curvature is not bounded in anisotropic effective models. However, this is a statement only about curvature expressions on minisuperspace, and the more relevant question is what happens to curvature along trajectories obtained by solving equations of motion. This demonstrates that dynamical equations must always be considered to draw conclusions for the singularity issue.

The approach to the classical singularity is best analyzed in Misner variables [157] consisting of the scale factor \(\Omega := - {1 \over 3}\log V\) and two anisotropy parameters β± defined such that

$${a_1} = {e^{- \Omega + {\beta _ +} + \sqrt 3 {\beta _ -}}},\quad {a_2} = {e^{- \Omega + {\beta _ +} - \sqrt 3 {\beta _ -}}},\quad {a_3} = {e^{- \Omega - 2{\beta _ +}}}.$$

The classical potential then takes the form

$$W(\Omega, {\beta _ \pm}) = {1 \over 2}{e^{- 4\Omega}}\left({{e^{- 8{\beta _ +}}} - 4{e^{- 2{\beta _ +}}}\cosh (2\sqrt 3 {\beta _ -}) + 2{e^{4\beta}} + (\cosh (4\sqrt 3 {\beta _ -}) - 1)} \right),$$

which at fixed Ω has three exponential walls rising from the isotropy point β± =0 and enclosing a triangular region (Figure 5).

Figure 5
figure 5

mpg-Movie (1.65 MB) Still from a Movie illustrating the Bianchi IX potential in the anisotropy plane and its exponentially rising walls. Positive values of the potential are drawn logarithmically with solid contour lines and negative values with dashed contour lines. (For video see appendix)

A cross section of a wall can be obtained by taking β =0 and β+ to be negative, in which case the potential becomes \(W(\Omega, {\beta _ +},0) \approx {1 \over 2}{e^{- 4\Omega - 8{\beta _ +}}}\). One thus obtains the picture of a point moving almost freely until it is reflected at a wall. In between reflections, the behavior is approximately given by the Kasner solution described before. This behavior with infinitely many reflections before the classical singularity is reached, can be shown to be chaotic [32], which suggests a complicated approach to classical singularities in general.

With the effective modification, however, the potential for fixed \(\Omega\) does not diverge and the walls, as shown in Figure 6, break down already at a small but non-zero volume [60]. As a function of densitized triad components the effective potential is illustrated in Figure 7, and as a function on the anisotropy plane in Figure 8. In this scenario, there are only finitely many reflections, which does not lead to chaotic behavior but instead results in asymptotic Kasner behavior [61].

Figure 6
figure 6

Approximate effective wall of finite height [60] as a function of x = −β+, compared to the classical exponential wall (upper dashed curve). Also shown is the exact wall W(p1, p1, (V/p1)2) (lower dashed curve), which for x smaller than the peak value coincides well with the approximation up to a small, nearly constant shift.

Figure 7
figure 7

mpg-Movie (1.59 MB) Still from a Movie illustrating the effective Bianchi IX potential and the movement and breakdown of its walls. The contours are plotted as in Figure 4 . (For video see appendix)

Figure 8
figure 8

mpg-Movie (1.66 MB) Still from a Movie illustrating the effective Bianchi IX potential in the anisotropy plane and its walls of finite height, which disappear at finite volume. Positive values of the potential are drawn logarithmically with solid contour lines and negative values with dashed contour lines. (For video see appendix)

Comparing Figure 5 with Figure 8 shows that in their center they are very close to each other, while strong deviations occur for large anisotropies. This demonstrates that most of the classical evolution, which mostly happens in the inner triangular region, is not strongly modified by the effective potential. Quantum effects are important only when anisotropies become too large, for instance when the system moves deep into one of the three valleys, or the total volume becomes small. In those regimes the quantum evolution will take over and describe the further behavior of the system.

4.9.3 Isotropic curvature suppression

If we use the potential for time coordinate t rather than τ, it is replaced by W/(p1p2p3), which in the isotropic reduction \({p^1} = {p^2} = {p^3} = {1 \over 4}{a^2}\) gives the curvature term ka−2. Although the anisotropic effective curvature potential is not bounded it is, unlike the classical curvature, bounded from above at any fixed volume. Moreover, it is bounded along the isotropy line and decays when a approaches zero. Thus, there is a suppression of the divergence in ka−2 when the closed isotropic model is viewed as embedded in a Bianchi IX model. Similarly to matter Hamiltonians, intrinsic curvature then approaches zero at zero scale factor.

This is a further illustration for the special nature of isotropic models compared to anisotropic ones. In the classical reduction, the pI in the anisotropic spin connection cancel such that the spin connection is a constant and no special steps are needed for its quantization. By viewing isotropic models within anisotropic ones, one can consistently realize the model and see a suppression of intrinsic curvature terms. Anisotropic models, on the other hand, do not have, and do not need, complete suppression since curvature functions can still be unbounded.

4.10 Anisotropy: Implications for inhomogeneities

Even without implementing inhomogeneous models the previous discussion allows some tentative conclusions as to the structure of general singularities. This is based on the BKL picture [31] whose basic idea is to study Einstein’s field equations close to a singularity. One can then argue that spatial derivatives become subdominant compared to time-like derivatives such that the approach should locally be described by homogeneous models, in particular the Bianchi IX model since it has the most freedom in its general solution.

Since spatial derivatives are present, though, they lead to small corrections and couple the geometries in different spatial points. One can visualize this by starting with an initial slice which is approximated by a collection of homogeneous patches. For some time, each patch evolves independently of the others, but this is not precisely true since coupling effects have been ignored. Moreover, each patch geometry evolves in a chaotic manner, which means that two initially nearby geometries depart rapidly from each other. The approximation can thus be maintained only if the patches are subdivided during the evolution, which goes on without limits in the approach to the singularity. There is, thus, more and more inhomogeneous structure being generated on arbitrarily small scales, which leads to a complicated picture of a general singularity.

This picture can be taken over to the effective behavior of the Bianchi IX model. Here, the patches do not evolve chaotically even though at larger volume they follow the classical behavior. The subdivision thus has to be done also for the initial effective evolution. At some point, however, when reflections on the potential walls stop, the evolution simplifies and subdivisions are no longer necessary. There is thus a lower bound to the scale of structure whose precise value depends on the initial geometries. Nevertheless, from the scale at which the potential walls break down one can show that structure formation stops at the latest when the discreteness scale of quantum geometry is reached [60]. This can be seen as a consistency test of the theory since structure below the discreteness could not be supported by quantum geometry.

We have thus a glimpse on the inhomogeneous situation with a complicated but consistent approach to a general classical singularity. The methods involved, however, are not very robust since the BKL scenario, which even classically is still at the level of a conjecture for the general case [32, 168], would need to be available as an approximation to quantum geometry. For more reliable results the methods need to be refined to take into account inhomogeneities properly.

4.11 Inhomogeneities

Allowing for inhomogeneities inevitably means to take a big step from finitely many degrees of freedom to infinitely many ones. There is no straightforward way to cut down the number of degrees of freedom to finitely many ones while being more general than in the homogeneous context. One possibility would be to introduce a small-scale cut-off such that only finitely many wave modes arise (e.g., through a lattice as is indeed done in some coherent state constructions [184]). This is in fact expected to happen in a discrete framework such as quantum geometry, but would at this stage of defining a model simply be introduced by hand.

For the analysis of inhomogeneous situations there are several different approximation schemes:

  • Use only isotropic quantum geometry and in particular its effective description, but couple to inhomogeneous matter fields. Problems in this approach are that back-reaction effects are ignored (which is also the case in most classical treatments) and that there is no direct way how to check modifications used in particular for gradient terms of the matter Hamiltonian. So far, this approach has led to a few indications of possible effects.

  • Start with the full constraint operator, write it as the homogeneous one plus correction terms from inhomogeneities, and derive effective classical equations. This approach is more ambitious since contact to the full theory is realized. So far, there are not many results since a suitable perturbation scheme has to be developed.

  • There are inhomogeneous symmetric models, such as Einstein-Rosen waves or the spherically symmetric one, which have infinitely many kinematical degrees of freedom but can be treated explicitly. Also here, contact to the full theory is present through the symmetry reduction procedure of Section 6. This procedure itself can be tested by studying those models between homogeneous ones and the full theory, but results can also be used for physical applications involving inhomogeneities. Many issues that are of importance in the full theory, such as the anomaly problem, also arise here and can thus be studied more explicitly.

4.12 Inhomogeneous matter with isotropic quantum geometry

Inhomogeneous matter fields cannot be introduced directly to isotropic quantum geometry since after the symmetry reduction there is no space manifold left for the fields to live on. There are then two different routes to proceed: One can simply take the classical field Hamiltonian and introduce effective modifications modeled on what happens to the isotropic Hamiltonian, or perform a mode decomposition of the matter fields and just work with the space-independent amplitudes. The latter is possible since the homogeneous geometry provides a background for the mode decomposition.

The basic question, for the example of a scalar field, then is how the metric coefficient in the gradient term of Equation (12), \(E_i^aE_i^b/\sqrt \vert{\det \,E}\vert\), would be replaced effectively. For the other terms, one can simply use the isotropic modification, which is taken directly from the quantization. For the gradient term, however, one does not have a quantum expression in this context and a modification can only be guessed. The problem arises since the inhomogeneous term involves inverse powers of E, while in the isotropic context the coefficient just reduces to \(\sqrt {\vert p\vert}\), which would not be modified at all. There is thus no obvious and unique way to find a suitable replacement.

A possible route would be to read off the modification from the full quantum Hamiltonian, or at least from an inhomogeneous model, which requires a better knowledge of the reduction procedure. Alternatively, one can take a more phenomenological point of view and study the effects of possible replacements. If the robustness of these effects to changes in the replacements is known, one can get a good picture of possible implications. So far, only initial steps have been taken and there is no complete programme in this direction.

Another approximation of the inhomogeneous situation has been developed in [70] by patching isotropic quantum geometries together to support an inhomogeneous matter field. This can be used to study modified dispersion relations to the extent that the result agrees with preliminary calculations performed in the full theory [115, 3, 4, 181, 182] even at a quantitative level. There is thus further evidence that symmetric models and their approximations can provide reliable insights into the full theory.

4.13 Inhomogeneity: Perturbations

With a symmetric background, a mode decomposition is not only possible for matter fields but also for geometry. The homogeneous modes can then be quantized as before, while higher modes are coupled as perturbations implementing inhomogeneities [120]. As with matter Hamiltonians before, one can then also deal with the gravitational part of the Hamiltonian constraint. In particular, there are terms with inverse powers of the homogeneous fields which receive modifications upon quantization. As with gradient terms in matter Hamiltonians, there are several options for those modifications which can only be restricted by relating them to the full Hamiltonian. This would require introducing the mode decomposition, analogously to symmetry conditions, at the quantum level and writing the full constraint operator as the homogeneous one plus correction terms.

An additional complication compared to matter fields is that one is now dealing with infinitely many coupled constraint equations since the lapse function N(x) is inhomogeneous, too. This function can itself be decomposed into modes \({\sum\nolimits_n}{N_n}{Y_n}(x)\), with harmonics \({Y_n}(x)\) according to the symmetry, and each amplitude Nn is varied independently giving rise to a separate constraint. The main constraint arises from the homogeneous mode, which describes how inhomogeneities affect the evolution of the homogeneous scale factors.

4.14 Inhomogeneous models

The full theory is complicated at several different levels of both conceptual and technical nature. For instance, one has to deal with infinitely many degrees of freedom, most operators have complicated actions, and interpreting solutions to all constraints in a geometrical manner can be difficult. Most of these complications are avoided in homogeneous models, in particular when effective classical equations are employed. These equations use approximations of expectation values of quantum geometrical operators which need to be known rather explicitly. The question then arises whether one can still work at this level while relaxing the symmetry conditions and bringing in more complications of the full theory.

Explicit calculations at a level similar to homogeneous models, at least for matrix elements of individual operators, are possible in inhomogeneous models, too. In particular the spherically symmetric model and cylindrically symmetric Einstein-Rosen waves are of this class, where the symmetry or other conditions are strong enough to result in a simple volume operator. In the spherically symmetric model, this simplification comes from the remaining isotropy subgroup isomorphic to U(1) in generic points, while the Einstein-Rosen model is simplified by polarization conditions that play a role analogous to the diagonalization of homogeneous models. With these models one obtains access to applications for black holes and gravitational waves, but also to inhomogeneities in cosmology.

In spherical coordinates x, ϑ, ϕ a spherically symmetric spatial metric takes the form

$${\rm{d}}{s^2} = {q_{xx}}(x,t){\rm{d}}{x^2} + {q_{\varphi \varphi}}(x,t){\rm{d}}{\Omega ^2}$$

with dΩ2 = dϑ2 + sin2ϑdϕ2. This is related to densitized triad components by [196, 136]

$$\vert {E^x}\vert = {q_{\varphi \varphi}},\quad {({E^\varphi})^2} = {q_{xx}}{q_{\varphi \varphi}},$$

which are conjugate to the other basic variables given by the Ashtekar connection component Ax and the extrinsic curvature component Kϕ:

$$\{{A_x}(x),{E^x}(y)\} = 8\pi G\gamma \delta (x,y),\quad \{\gamma {K_\varphi}(x),{E^\varphi}(y)\} = 16\pi G\gamma \delta (x,y).$$

Note that we use the Ashtekar connection for the inhomogeneous direction x but extrinsic curvature for the homogeneous direction along symmetry orbits [75]. Connection and extrinsic curvature components for the ϕ-direction are related by \(A_\varphi ^2 = \Gamma _\varphi ^2 + {\gamma ^2}K_\varphi ^2\) with the spin connection component

$${\Gamma _\varphi} = - {{{E^{x^{\prime}}}} \over {2{E^\varphi}}}.$$
(38)

Unlike in the full theory or homogeneous models, Aϕ is not conjugate to a triad component but to [52]

$${P^\varphi} = \sqrt {4{{({E^\varphi})}^2} - A_\varphi ^{- 2}{{({P^\beta})}^2}}$$

with the momentum Pβ conjugate to a U(1)-gauge angle β. This is a rather complicated function of both triad and connection variables such that the volume \(V = 4\pi \int \sqrt {\vert {E^x}\vert} {E^\varphi}{\rm{d}}x\) would have a rather complicated quantization. It would still be possible to compute the full volume spectrum, but with the disadvantage that volume eigenstates would not be given by triad eigenstates such that computations of many operators would be complicated [74]. This can be avoided by using extrinsic curvature which is conjugate to the triad component [75]. Moreover, this is also in accordance with a general scheme to construct Hamiltonian constraint operators for the full theory as well as symmetric models [194, 42, 58].

The constraint operator in spherical symmetry is given by

$$H(N) = - {(2G)^{- 1}}\int\nolimits_B {{\rm{d}}x} N(x)\vert {E^x}{\vert ^{- 1/2}}((K_\varphi ^2{E^\varphi} + 2{K_\varphi}{K_x}{E^x}) + (1 - \Gamma _\varphi ^2){E^\varphi} + 2\Gamma _\varphi ^{\prime}{E^x})$$
(39)

accompanied by the diffeomorphism constraint

$$D[{N^x}] = {(2G)^{- 1}}\int\nolimits_B {{N^x}} (x)(2{E^\varphi}K_\varphi ^{\prime} - {K_x}{E^{x^{\prime}}}).$$
(40)

We have expressed this in terms of Kx for simplicity, keeping in mind that as the basic variable for quantization we will later use the connection component Ax.

Since the Hamiltonian constraint contains the spin connection component Γϕ given by (38), which contains inverse powers of densitized triad components, one can expect effective classical equations with modifications similar to the Bianchi IX model. However, the situation is now much more complicated since we have a system constrained by many constraints with a non-Abelian algebra. Simply replacing the inverse of with a bounded function as before will change the constraint algebra and thus most likely lead to anomalies. It is currently open if a more refined replacement can be done where not only the spin connection but also the extrinsic curvature terms are modified. This issue has the potential to shed light on many questions related to the anomaly issue. It is one of the cases where models between homogeneous ones, where the anomaly problem trivializes, and the full theory are most helpful.

4.15 Inhomogeneity: Results

There are some results obtained for inhomogeneous systems. We have already discussed glimpses from the BKL picture, which used loop results only for anisotropic models. Methods described in this section have led to some preliminary insights into possible cosmological scenarios.

4.15.1 Matter gradient terms and small-a effects

When an inhomogeneous matter Hamiltonian is available it is possible to study its implications on the cosmic microwave background with standard techniques. With modifications of densities there are then different regimes since the part of the inflationary era responsible for the formation of currently visible structure can be in the small-a or large-a region of the effective density.

The small-a regime below the peak of effective densities has more dramatic effects since inflation can here be provided by quantum geometry effects alone and the matter behavior changes to be anti-frictional [45, 77]. Mode evolution in this regime has been investigated for a particular choice of gradient term and using a power-law approximation for the effective density at small a, with the result that there are characteristic signatures [130]. As in standard inflation models the spectrum is nearly scale invariant, but its spectral index is slightly larger than one (blue tilt) as compared to slightly smaller than one (red tilt) for single-field inflaton models. Since small scale factors at early stages of inflation generate structure which today appears on the largest scales, this implies that low multipoles of the power spectrum should have a blue tilt. The running of the spectral index in this regime can also be computed but depends only weakly on ambiguity parameters.

The main parameter then is the duration of loop inflation. In the simplest scenario, one can assume only one inflationary phase, which would require huge values for the ambiguity parameter j. This is unnatural and would imply that the spectrum is blue on almost all scales, which is in conflict with present observations. Thus, not only conceptual arguments but also cosmological observations point to smaller values for j, which is quite remarkable.

In order to have sufficient inflation to make the universe big enough one then needs additional stages provided by the behavior of matter fields. One still does not need an inflaton since now the details of the expansion after the structure generating phase are less important. Any matter field being driven away from its potential minimum during loop inflation and rolling down its potential thereafter suffices. Depending on the complexity of the model there can be several such phases.

4.15.2 Matter gradient terms and large-a effects

At larger scale factors above the peak of effective densities there are only perturbative corrections from loop effects. This has been investigated with the aim of finding trans-Planckian corrections to the microwave background, also here with a particular gradient term. In this model, cancellations have been observed that imply that corrections appear only at higher orders of the perturbation series and are too weak to be observable [126].

A common problem of both analyses is that the robustness of the observed effects has not yet been studied. This is in particular a pressing problem since one modification of the gradient term has been chosen without further motivation. Moreover, the modifications in both examples were different. Without a more direct derivation of the modifications from inhomogeneous models or the full theory one can only rely on a robustness analysis to show that the effects can be trusted. In particular the cancellation in the second example must be shown to be realized for a larger class of modifications.

4.15.3 Non-inflationary structure formation

Given a modification of the gradient term, one obtains effective equations for the matter field, which for a scalar results in a modified Klein-Gordon equation. After a mode decomposition, one can then easily see that all the modes behave differently at small scales with the classical friction replaced by anti-friction as in Section 4.5. Thus, not only the average value of the field is driven away from its potential minimum but also higher modes are being excited. The coupled dynamics of all the modes thus provides a scenario for structure formation, which does not rely on inflation but on the anti-friction effect of loop cosmology.

Even though all modes experience this effect, they do not all see it in the same way. The gradient term implies an additive contribution to the potential proportional to k2 for a mode of wave number k, which also depends on the metric in a way determined by the gradient term modification. For larger scales, the additional term is not essential and their amplitudes will be pushed to similar magnitudes, suggesting scale invariance for them. The potential relevant for higher modes, however, becomes steeper and steeper such that they are less excited by antifriction and retain a small initial amplitude. In this way, the structure formation scenario provides a dynamical mechanism for a small-scale cut-off, possibly realizing older expectations [165, 166].

4.15.4 Stability

As already noted, inhomogeneous matter Hamiltonians can be used to study the stability of cosmological equations in the sense that matter does not propagate faster than light. The modified behavior of homogeneous modes has led to the suspicion that loop cosmology is not stable [94, 95] since other cosmological models displaying super-inflation have this problem. A detailed analysis of the loop equations, however, shows that the equations as they arise from modifications are automatically stable. While the homogeneous modes display super-inflationary and anti-frictional behavior, they are not relevant for matter propagation. Modes relevant for propagation, on the other hand, are modified differently in such a manner that the total behavior is stable [129]. Most importantly, this is an example where an inhomogeneous matter Hamiltonian with its modifications must be used and the qualitative result of stability can be shown to be robust under possible changes of the effective modification. This shows that reliable conclusions can be drawn for important issues without a precise definition of the effective inhomogeneous behavior.

4.16 Summary

Loop cosmology is an effective description of quantum effects in cosmology, obtained in a framework of a background independent and non-perturbative quantization. There is mainly one change compared to classical equations coming from modified densities in matter Hamiltonians or also anisotropy potentials. These modifications are non-perturbative as they contain inverse powers of the Planck length and thus the gravitational constant, but also perturbative corrections arise from curvature terms, which are now being studied.

The non-perturbative modification alone is responsible for a surprising variety of phenomena, which all improve the behavior in classical cosmology. Nevertheless, the modification had not been motivated by phenomenology but derived through the background independent quantization. Details of its derivation in cosmological models and its technical origin will now be reviewed in Section 5, before we come to a discussion of the link to the full theory in Section 6.

5 Loop quantization of symmetric models

Analogies prove nothing, but they can make one feel more at home.

Sigmund Freud

Introductory Lectures on Psychoanalysis

In full loop quantum gravity, the quantum representation is crucial for the foundation of the theory. The guiding theme there is background independence, which requires one to smear the basic fields in a particular manner to holonomies and fluxes. In this section, we will see what implications this has for composite operators and the physical effects they entail. We will base this analysis on symmetric models in order to be able to perform explicit calculations.

Symmetries are usually introduced in order to simplify calculations or make them possible in the first place. However, symmetries can sometimes also lead to complications in conceptual questions if the additional structure they provide is not fully taken into account. In the present context, it is important to realize that the action of a symmetry group on a space manifold provides a partial background such that the situation is always slightly different from the full theory. If the symmetry is strong, such as in homogeneous models, other representations such as the Wheeler-DeWitt representation can be possible even though the fact that a background has been used may not be obvious. While large scale physics is not very sensitive to the representation used, it becomes very important on the smallest scales, which we have to take into account when the singularity issue is considered.

Instead of looking only at one symmetric model, where one may have different possibilities to choose the basic representation, one should thus keep the full view on different models as well as the full theory. In fact, in loop quantum gravity it is possible to relate models and the full theory such that symmetric states and basic operators, and thus the representation, can be derived from the unique background independent representation of the full theory. We will describe this in detail in Section 6, after having discussed the construction of quantum models in the present section. Without making use of the relation to the full theory, one can construct models by analogy. This means that quantization steps are modeled on those which are known to be crucial in the full theory, which starts with the basic representation and continues to the Hamiltonian constraint operator. One can then disentangle places where additional input as compared to the full theory is needed and which implications it has.

5.1 Symmetries and backgrounds

It is impossible to introduce symmetries in a completely background independent manner. The mathematical action of a symmetry group is defined by a mapping between abstract points, which do not exist in a diffeomorphism invariant setting (if one, for instance, considers only equivalence classes up to arbitrary diffeomorphisms).

More precisely, while the full theory has as background only a differentiable or analytic manifold Σ, a symmetric model has as background a symmetric manifold (Σ, S) consisting of a differentiable or analytic manifold Σ together with an action of a symmetry group S: Σ → Σ. How strong the additional structure is depends on the symmetry used. The strongest symmetry in gravitational models is realized with spatial isotropy, which implies a unique spatial metric up to a scale factor. The background is thus equivalent to a conformal space.

All constructions in a given model must take its symmetry into account since otherwise its particular dynamics, for instance, could not be captured. The structure of models thus depends on the different types of background realized for different symmetry groups. This can not only lead to simplifications but also to conceptual differences, and it is always instructive to keep the complete view on different models as well as the full theory. Since the loop formalism is general enough to encompass all relevant models, there are many ways to compare and relate different systems. It is thus possible to observe characteristic features of (metric) background independence even in cases where more structure is available.

5.2 Isotropy

Isotropic models are described purely in terms of the scale factor a(t) such that there is only a single kinematical degree of freedom. In connection variables, this is parameterized by the triad component p conjugate to the connection component c.

If we restrict ourselves to invariant connections of a given form, it suffices to probe them with only special holonomies. For an isotropic connection \(A_a^i = \tilde c\Lambda _I^iw_a^I\) (see Appendix B.2) we can choose holonomies along one integral curve of a symmetry generator XI. They are of the form

$${h_I} = \exp \int {A_a^iX_I^a{\tau _i} = \cos {1 \over 2}\mu c} + 2\Lambda _I^i{\tau _i}\sin {1 \over 2}\mu c$$
(41)

where μ depends on the parameter length of the curve and can be any real number (thanks to homogeneity, path ordering is not necessary). Since knowing the values \(\cos {1 \over 2}\mu\) and \(\sin {1 \over 2}\mu c\) for all μ uniquely determines the value of c, which is the only gauge invariant information contained in the connection, these holonomies describe the configuration space of connections completely.

This illustrates how symmetric configurations allow one to simplify the constructions behind the full theory. But it also shows which effects the presence of a partial background can have on the formalism [15]. In the present case, the background enters through the left-invariant 1-forms ωI defined on the spatial manifold whose influence is contained in the parameter μ. All information about the edge used to compute the holonomy is contained in this single parameter, which leads to degeneracies compared to the full theory. Most importantly, one cannot distinguish between the parameter length and the spin label of an edge: Taking a power of the holonomy in a non-fundamental representation simply rescales μ, which could just as well come from a longer parameter length. That this is related to the presence of a background can be seen by looking at the roles of edges and spin labels in the full theory. There, both concepts are independent and appear very differently. While the embedding of an edge, including its parameter length, is removed by diffeomorphism invariance, the spin label remains well-defined and is important for ambiguities of operators. In the model, however, the full diffeomorphism invariance is not available such that some information about edges remains in the theory and merges with the spin label. Issues like that have to be taken into account when constructing operators in a model and comparing with the full theory.

The functions appearing in holonomies for isotropic connections define the algebra of functions on the classical configuration space which, together with fluxes, is to be represented on a Hilbert space. This algebra does not contain arbitrary continuous functions of c but only almost periodic ones of the form [15]

$$f(c) = \sum\limits_\mu {{f_\mu}\exp (i\mu c/2)}$$
(42)

where the sum is over a countable subset of ℝ. This is analogous to the full situation, reviewed in Section 3.4, where matrix elements of holonomies define a special algebra of continuous functions of connections. As in this case, the algebra can be represented as the set of all continuous functions on a compact space, called its spectrum. This compactification can be imagined as being obtained from enlarging the classical configuration space ℝ by adding points, and thus more continuity conditions, until only functions of the given algebra survive as continuous ones. A well-known example is the one point compactification, which is the spectrum of the algebra of continuous functions f for which \({\lim\nolimits_{x \rightarrow - \infty}}f(x) = {\lim\nolimits_{x \rightarrow - \infty}}f(x)\) exists. In this case, one just needs to add a single point at infinity.

In the present case, the procedure is more complicated and leads to the Bohr compactification \({\overline{\mathbb R}_{{\rm{Bohr}}}}\), which contains ℝ densely. It is very different from the one point compactification, as can be seen from the fact that the only function which is continuous on both spaces is the zero function. In contrast to the one point compactification, the Bohr compactification is an Abelian group, just like ℝ itself. Moreover, there is a one-to-one correspondence between irreducible representations of ℝ and irreducible representations of \({\overline{\mathbb R}_{{\rm{Bohr}}}}\), which can also be used as the definition of the Bohr compactification. Representations of \({\overline{\mathbb R}_{{\rm{Bohr}}}}\) are thus labeled by real numbers and given by \({\rho _\mu}:{\overline{\mathbb R}_{{\rm{Bohr}}}} \rightarrow {\mathbb C},c \mapsto {e^{i\mu c}}\). As with any compact group, there is a unique normalized Haar measure dμ(c) given by

$$\int\nolimits_{{{\overline{\mathbb R}}_{{\rm{Bohr}}}}} {f(c)} \;{\rm{d}}\mu (c) = \mathop {\lim}\limits_{T \to \infty} {1 \over {2T}}\int\nolimits_{- T}^T {f(c){\rm{d}}c,}$$
(43)

where on the right hand side the Lebesgue measure on ℝ is used.

The Haar measure defines the inner product for the Hilbert space \({L^2}({\overline{\mathbb R}_{{\rm{Bohr}}}},{\rm{d}}\mu (c))\) of square integrable functions on the quantum configuration space. As one can easily check, exponentials of the form ⟨c|μ⟩ = eiμc/2 are normalized and orthogonal to each other for different μ,

$$\langle {\mu _1}\left\vert {{\mu _2}} \right.\rangle = {\delta _{{\mu _1},{\mu _2}}},$$
(44)

which demonstrates that the Hilbert space is not separable.

Similarly to holonomies, one needs to consider fluxes only for special surfaces, and all information is contained in the single number p. Since it is conjugate to c, it is quantized to a derivative operator

$$\hat{p}= - {1 \over 3}i\gamma \ell _{\rm{P}}^2{{\rm{d}} \over {{\rm{d}}c}},$$
(45)

whose action

$$\hat{p}\left\vert {\mu \rangle} \right. = {1 \over 6}\gamma \ell _{\rm{P}}^2\mu \left\vert {\mu \rangle = :} \right.{p_\mu}\left\vert {\mu \rangle} \right.$$
(46)

on basis states |μ⟩ can easily be determined. In fact, the basis states are eigenstates of the flux operator, which demonstrates that the flux spectrum is discrete (all eigenstates are normalizable).

This property is analogous to the full theory with its discrete flux spectra, and similarly it implies discrete quantum geometry. We thus see that the discreteness survives the symmetry reduction in this framework [37]. Similarly, the fact that only holonomies are represented in the full theory but not connection components is realized in the model, too. In fact, we have so far represented only exponentials of c, and one can see that these operators are not continuous in the parameter μ. Thus, an operator quantizing c directly does not exist on the Hilbert space. These properties are analogous to the full theory, but very different from the Wheeler-DeWitt quantization. In fact, the resulting representations in isotropic models are inequivalent. While the representation is not of crucial importance when only small energies or large scales are involved [18], it becomes essential at small scales which are in particular realized in cosmology.

5.3 Isotropy: Matter Hamiltonian

We now know how the basic quantities p and c are quantized, and can use the operators to construct more complicated ones. Of particular importance, also for cosmology, are matter Hamiltonians where now not only the matter field but also geometry is quantized. For an isotropic geometry and a scalar, this requires us to quantize |p|−3/2 for the kinetic term and |p|3/2 for the potential term. The latter can be defined readily as \(\vert\hat p\vert^{3/2}\), but for the former we need an inverse power of p. Since \({\hat p}\) has a discrete spectrum containing zero, a densely defined inverse does not exist.

At this point, one has to find an alternative route to the quantization of d(p) = |p|−3/2, or else one could only conclude that there is no well-defined quantization of matter Hamiltonians as a manifestation of the classical divergence. In the case of loop quantum cosmology it turns out, following a general scheme of the full theory [193], that one can reformulate the classical expression in an equivalent way such that quantization becomes possible. One possibility is to write, similarly to (13)

$$d(p) = {\left({{1 \over {3\pi \gamma G}}\sum\limits_{I = 1}^3 {{\rm{tr}}} \left({{\tau _I}{h_I}\{h_I^{- 1},\sqrt V \}} \right)} \right)^6}$$

where we use holonomies of isotropic connections and the volume V = |p|3/2. In this expression we can insert holonomies as multiplication operators and the volume operator, and turn the Poisson bracket into a commutator. The result

$$\widehat {d(p)} = {\left({8i{\gamma ^{- 1}}\ell _{\rm{P}}^{- 2}(\sin {1 \over 2}c\sqrt {\hat{V}} \cos {1 \over 2}c - \cos {1 \over 2}c\sqrt {\hat{V}} \sin {1 \over 2}c)} \right)^6}$$
(47)

is not only a densely defined operator but even bounded, which one can easily read off from the eigenvalues [41]

$$\widehat {d(p)}\left\vert {\mu \rangle} \right. = {\left({4{\gamma ^{- 1}}\ell _{\rm{P}}^{- 2}(\sqrt {{V_{\mu + 1}}} - \sqrt {{V_{\mu - 1}}})} \right)^6}\left\vert {\mu \rangle} \right.$$
(48)

with Vμ = |pμ|3/2 and pμ from (46).

Rewriting a classical expression in such a manner can always be done in many equivalent ways, which in general all lead to different operators. In the case of |p| 3/2, we highlight the choice of the representation in which to take the trace (understood as the fundamental representation above) and the power of |p| in the Poisson bracket (\(\sqrt V = \vert p\vert^{3/4}\) above). This freedom can be parameterized by two ambiguity parameters \(j \in {1 \over 2}{\mathbb N}\) for the representation and 0 < l < 1 for the power such that

$$d(p) = {\left({{3 \over {8\pi \gamma Glj(j + 1)(2j + 1)}}\sum\limits_{I = 1}^3 {{\rm{t}}{{\rm{r}}_j}} ({\tau _I}{h_I}\{h_I^{- 1},\vert p{\vert ^l}\})} \right)^{3/(2 - 2l)}}.$$

Following the same procedure as above, we obtain eigenvalues [47, 50]

$$\widehat {d(p)}_{j,l}^{(\mu)} = {\left({{9 \over {\gamma \ell _{\rm{P}}^2lj(j + 1)(2j + 1)}}\sum\limits_{k = - j}^j {k\vert {p_{\mu + 2k}}\vert}^{l}} \right)^{3/(2 - 2l)}},$$

which, for larger j, can be approximated by (25), see also Figure 9. This provides the basis for loop cosmology as described in Section 4.

Figure 9
figure 9

Discrete subset of eigenvalues of \(\widehat {d(p)}\) (left) for two choices of j (and \(l = {3 \over 4}\)), together with the approximation d(p)j,l from Equation ( 25 ) and small-p power laws. The classical divergence at small p, where the behavior is strongly modified, is cut off. The right panel shows the dependence of the initial increase on l.

Notice that operators for the scale factor, volume or their inverse powers do not refer to observable quantities. It can thus be dangerous, though suggestive, to view their properties as possible bounds on curvature. The importance of operators for inverse volume comes from the fact that this appears in matter Hamiltonians, and thus the Hamiltonian constraint of gravity. Properties of those operators such as their boundedness or unboundedness can then determine the dynamical behavior (see, e.g., [54]).

5.4 Isotropy: Hamiltonian constraint

Dynamics is controlled by the Hamiltonian constraint, which classically gives the Friedmann equation. Since the classical expression (28) contains the connection component c, we have to use holonomy operators. In the quantum algebra we only have almost periodic functions at our disposal, which does not include polynomials such as c2. Quantum expressions can therefore only coincide with the classical one in appropriate limits, which in isotropic cosmology is realized for small extrinsic curvature, i.e., small c in the flat case. We thus need an almost periodic function of c, which for small approaches c2. This can easily be found, e.g., the function sin2c. Again, the procedure is not unique since there are many such possibilities, e.g., δ−2 sin2δc, and more quantization ambiguities ensue. In contrast to the density | p| −3/2, where we also used holonomies in the reformulation, the expressions are not equivalent to each other classically but only in the small curvature regime. As we will discuss shortly, the resulting new terms have the interpretation of higher order corrections to the classical Hamiltonian.

One can restrict the ambiguities to some degree by modeling the expression on that of the full theory. This means that one does not simply replace c2 by an almost periodic function, but uses holonomies tracing out closed loops formed by symmetry generators [42, 46]. Moreover, the procedure can be embedded in a general scheme that encompasses different models and the full theory [194, 42, 58], further reducing ambiguities. In particular models with non-zero intrinsic curvature on their symmetry orbits, such as the closed isotropic model, can then be included in the construction. One issue to keep in mind is the fact that “holonomies” are treated differently in models and the full theory. In the latter case, they are ordinary holonomies along edges, which can be shrunk and then approximate connection components. In models, on the other hand, one sometimes uses direct exponentials of connection components without integration. In such a case, connection components are approximated only when they are small; if they are not, the corresponding objects such as the Hamiltonian constraint receive infinitely many correction terms of higher powers in curvature (similarly to effective actions). The difference between both ways of dealing with holonomies can be understood in inhomogeneous models, where they are both realized for different connection components.

In the flat case the construction is easiest, related to the Abelian nature of the symmetry group. One can directly use the exponentials hI in (41), viewed as 3-dimensional holonomies along integral curves, and mimic the full constraint where one follows a loop to get curvature components of the connection \(A_a^i\). Respecting the symmetry, this can be done in the model with a square loop in two independent directions I and J. This yields the product \({h_I}{h_J}h_I^{- 1}h_J^{- 1}\), which appears in a trace, as in (15), together with a commutator \({h_K}[h_K^{- 1},\hat V]\) using the remaining direction K. The latter, following the general scheme of the full theory reviewed in Section 3.6, quantizes the contribution \(\sqrt {\vert p\vert}\) to the constraint, instead of directly using the simpler \(\sqrt {\vert \hat p\vert}\).

Taking the trace one obtains a diagonal operator

$$\sin \left({{1 \over 2}\delta c} \right)\hat V\cos \left({{1 \over 2}\delta c} \right) - \cos \left({{1 \over 2}\delta c} \right)\hat V\sin \left({{1 \over 2}\delta c} \right)$$

in terms of the volume operator, as well as the multiplication operator

$${\sin ^2}({1 \over 2}\delta c){\cos ^2}({1 \over 2}\delta c) = {\sin ^2}(\delta c).$$

In the triad representation where instead of working with functions ⟨c|ψ⟩ = ψ(c) one works with the coefficients ψμ in an expansion \(\vert \psi \rangle = {\sum\nolimits_\mu}{\psi _\mu}\vert \mu \rangle\), this operator is the square of a difference operator. The constraint equation thus takes the form of a difference equation [46, 77, 15]

$$\begin{array}{*{20}c}{\;\;({V_{\mu + 5\delta}} - {V_{\mu + 3\delta}}){e^{ik}}{\psi _{\mu + 4\delta}}(\phi) - (2 + {k^2}{\gamma ^2}{\delta ^2})({V_{\mu + \delta}} - {V_{\mu - \delta}}){\psi _\mu}(\phi)} \\{+ ({V_{\mu - 3\delta}} - {V_{\mu - 5\delta}}){e^{- ik}}{\psi _{\mu - 4\delta}}(\phi) = - {{16\pi} \over 3}G{\gamma ^3}{\delta ^3}\ell _{\rm{P}}^2{{\hat{H}}_{{\rm{matter}}}}(\mu){\psi _\mu}(\phi)} \\ \end{array}$$
(49)

for the wave function ψμ, which can be viewed as an evolution equation in internal time μ. (Note that this equation is not valid for k = −1 since the derivation via a Hamiltonian formulation is not available in this case.) Thus, discrete spatial geometry implies a discrete internal time [43]. The equation above results in the most direct way from a non-symmetric constraint operator with gravitational part acting as

$$\hat{H}\left\vert {\mu \rangle} \right. = {3 \over {16\pi G{\gamma ^3}{\delta ^3}\ell _{\rm{P}}^2}}({V_{\mu + \delta}} - {V_{\mu - \delta}})({e^{- ik}}\vert \mu + 4\delta \rangle - (2 + {k^2}{\gamma ^2}{\delta ^2})\left\vert {\mu \rangle} \right. + {e^{ik}}\vert \mu - 4\delta \rangle).$$

One can symmetrize this operator and obtain a difference equation with different coefficients, which we do here after multiplying the operator with \(\widehat{{\mathop{\rm sgn}} p}\) for reasons that will be discussed in the context of singularities in Section 5.15. The resulting difference equation is

$$\begin{array}{*{20}c}{\quad \quad (\vert {\Delta _\delta}V\vert (\mu + 4\delta) + \vert {\Delta _\delta}V(\mu)){e^{ik}}{\psi _{\mu + 4\delta}}(\phi) - (2 + {k^2}{\gamma ^2}{\delta ^2})\vert {\Delta _\delta}V\vert (\mu){\psi _\mu}(\phi)} \\{+ (\vert {\Delta _\delta}V\vert (\mu - 4\delta) + \vert {\Delta _\delta}V\vert (\mu)){e^{- ik}}{\psi _{\mu - 4\delta}}(\phi) = - {{32\pi} \over 3}G{\gamma ^3}{\delta ^3}\ell _{\rm{P}}^2{{\hat H}_{{\rm{matter}}}}(\mu){\psi _\mu}(\phi)} \\\end{array}$$
(50)

where |ΔδV:= sgn(μ)(Vμ+δVμδ) = |Vμ+δVμδ|.

Since \(\sin c\vert \mu \rangle = - {1 \over 2}i(\vert \mu + 2\rangle - \vert \mu - 2\rangle)\), the difference equation is of higher order, even formulated on an uncountable set, and thus has many independent solutions. Most of them, however, oscillate on small scales, i.e., between μ and μ + with small integer m. Others oscillate only on larger scales and can be viewed as approximating continuum solutions. The behavior of all the solutions leads to possibilities for selection criteria of different versions of the constraint since there are quantization choices. Most importantly, one chooses the routing of edges to construct the square holonomy, again the spin of a representation to take the trace [116, 201], and factor ordering choices between quantizations of c2 and \(\sqrt {\vert p\vert}\). All these choices also appear in the full theory such that one can draw conclusions for preferred cases there.

When the symmetry group is not Abelian and there is non-zero intrinsic curvature, the construction is more complicated. For non-Abelian symmetry groups integral curves as before do not form a closed loop and one needs a correction term related to intrinsic curvature components [46, 62]. Moreover, the classical regime is not as straightforward to specify since connection components are not necessarily small when there is intrinsic curvature. A general scheme encompassing intrinsic curvature, other symmetric models and the full theory will be discussed in Section 5.14.

5.5 Semiclassical limit and correction terms

When replacing c2 by holonomies we have modified the constraint as a function on the classical phase space. This is necessary since otherwise the function cannot be quantized, but is different from the quantization of densities because now the replacements are not equivalent to the original constraint classically. Also the limit \({\lim\nolimits_{\delta \rightarrow 0}}{\delta ^{- 2}}{\sin ^2}\delta c\), which would give the classical result, does not exist at the operator level.

This situation is different from the full theory, again related to the presence of a partial background [15]. There, the parameter length of edges used to construct appropriate loops is irrelevant and thus can shrink to zero. In the model, however, changing the edge length with respect to the background does change the operator and the limit does not exist. Intuitively, this can be understood as follows: The full constraint operator (15) is a vertex sum obtained after introducing a discretization of space used to choose loops αIJ. This classical regularization sums over all tetra-hedra in the discretization, whose number diverges in the limit where the discretization size shrinks to zero. In the quantization, however, almost all these contributions vanish since a tetrahedron must contain a vertex of a state in order to contribute non-trivially. The result is independent of the discretization size once it is fine enough, and the limit can thus be taken trivially.

In a homogeneous model, on the other hand, contributions from different tetrahedra of the triangulation must be identical owing to homogeneity. The coordinate size of tetrahedra drops out of the construction in the full background independent quantization, as emphasized in Section 3.6, which is part of the reason for the discretization independence. In a homogeneous configuration the number of contributions thus increases in the limit, but their size does not change. This results in an ill-defined limit as we have already seen within the model itself.

The difference between models and the full theory is thus only a consequence of the symmetry and not of different approaches. This will also become clear later in inhomogeneous models where one obtains a mixture between the two situations. Moreover, in the full theory one has a situation similar to symmetric models if one does not only look at the operator limit when the regularization is removed but also checks the classical limit on semiclassical states. In homogeneous models, the expression in terms of holonomies implies corrections to the classical constraint when curvature becomes larger. This is in analogy to other quantum field theories where effective actions generally have higher curvature terms. In the full theory, those correction terms can be seen when one computes expectation values of the Hamiltonian constraint in semiclassical states peaked at classical configurations for the connection and triad. When this classical configuration becomes small in volume or large in curvature, correction terms to the classical constraint arise. In this case, the semiclassical state provides the background with respect to which these corrections appear. In a homogeneous model, the symmetry already provides a partial background such that correction terms can be noticed already for the constraint operator itself.

5.5.1 WKB approximation

There are different procedures to make contact between the difference equation and classical constraints. The most straightforward way is to expand the difference operators in a Taylor series, assuming that the wave function is sufficiently smooth. On large scales, this indeed results in the Wheeler-DeWitt equation as a continuum limit in a particular ordering [44]. From then on, one can use the WKB approximation or Wigner functions as usually.

That this is possible may be surprising because as just discussed the continuum limit δ → 0 does not exist for the constraint operator. And indeed, the limit of the constraint equation, i.e., the operator applied to a wave function, does not exist in general. Even for a wave function the limit δ→0 does not exist in general since some solutions are sensitive to the discreteness and do not have a continuum limit at all. When performing the Taylor expansion we already assumed certain properties of the wave function such that the continuum limit does exist. This then reduces the number of independent wave functions to that present in the Wheeler-DeWitt framework, subject to the Wheeler-DeWitt equation. That this is possible demonstrates that the constraint in terms of holonomies does not have problems with the classical limit.

The Wheeler-DeWitt equation results at leading order, and in addition higher order terms arise in an expansion of difference operators in terms of δ or γ. Similarly, after the WKB or other semiclassical approximation there are correction terms to the classical constraint in terms of γ as well as [99].

This procedure is intuitive, but it is not suitable for inhomogeneous models where the Wheeler-DeWitt representation becomes ill-defined. One can evade this by performing the continuum and semiclassical limit together. This again leads to corrections in terms of γ as well as , which are mainly of the following form [29]: Matter Hamiltonians receive corrections through the modified density d(p), and there are similar terms in the gravitational part containing \(\sqrt {\vert p\vert}\). These are purely from triad coefficients, and similarly connection components lead to higher order corrections as well as additional contributions summarized in a quantum geometry potential. A possible interpretation of this potential in analogy to the Casimir effect has been put forward in [128]. A related procedure to extract semiclassical properties from the difference operator, based on the Bohmian interpretation of quantum mechanics, has been discussed in [185].

5.5.2 Effective formulation

In general, one does not only expect higher order corrections for a gravitational action but also higher derivative terms. The situation is then qualitatively different since not only correction terms to a given equation arise, but also new degrees of freedom coming from higher derivatives being independent of lower ones. In a WKB approximation, this could be introduced by parameterizing the amplitude of the wave function in a suitable way, but it has not been worked out yet. An alternative approach makes use of a geometrical formulation of quantum mechanics [26], which not only provides a geometrical picture of the classical limit but also a clear-cut procedure for computing effective Hamiltonians in analogy to effective actions [73].

Instead of using linear operators on a Hilbert space, one can formulate quantum mechanics on an infinite-dimensional phase space. This space is directly obtained from the Hilbert space where the inner product defines a metric as well as a symplectic form on its linear vector space (which in this way even becomes Kähler). This formulation brings quantum mechanics conceptually much closer to classical physics, which also facilitates a comparison in a semiclassical analysis.

We thus obtain a quantum phase space with infinitely many degrees of freedom, together with a flow defined by the Schrödinger equation. Operators become functions on this phase space through expectation values. Coordinates can be chosen by suitable parameterizations of a general wave function, in particular using the expectation values \(q = \langle \hat q\rangle\) and \(p = \langle \hat p\rangle\) together with uncertainties and higher moments. The projection \(\pi :{\mathcal H} \rightarrow {\mathbb R}^2,\psi \mapsto (\langle \psi \vert \hat q\vert \psi \rangle, \langle \psi \vert \hat p\vert \psi \rangle)\) defines the quantum phase space as a fiber bundle over the classical phase space with infinite-dimensional fibers. Sections of this bundle can be defined by embedding the classical phase space into the quantum phase space by means of suitable semiclassical states.

For a harmonic oscillator this embedding can be done by coherent states which are preserved by the quantum evolution. This means that the quantum flow is tangential to the embedding of the classical phase space such that it agrees with the classical flow. The harmonic oscillator thus does not receive quantum corrections as is well known from effective actions for free field theories. Other systems, however, behave in a more complicated manner where in general states spread. This means that additional coordinates of the quantum phase space are dynamical and may become excited. If this is the case, the quantum flow differs from the classical flow and an effective Hamiltonian arises with correction terms that can be computed systematically. This effective Hamiltonian is given by the expectation value ⟨Ĥ⟩ in approximate coherent states [16, 206, 188]. In these calculations, one can include higher degrees of freedom along the fibers, which, through the effective equations of motion, can be related to higher derivatives or higher curvature in the case of gravity.

For a constrained system, such as gravity, one has to compute the expectation value of the Hamiltonian constraint, i.e., first go to the classical picture and then solve equations of motion. Otherwise, there would simply be no effective equations left after the constraints would already have been solved. This is the same procedure as in standard effective actions, which one can also formulate in a constrained manner if one chooses to parameterize time. Indeed, also for non-constrained systems agreement between the geometrical way to derive effective equations and standard path integral methods has been shown for perturbations around a harmonic oscillator [73].

5.6 Homogeneity

A Hamiltonian formulation is available for all homogeneous models of Bianchi class A [111], which have structure constants \(C_{JK}^I\) fulfilling \(C_{JI}^I = 0\). The structure constants also determine left-invariant 1-forms ωI in terms of which one can write a homogeneous connection as \(A_a^i = \tilde \phi _I^iw_a^I\) (see Appendix B.1) where all freedom is contained in the x-independent \(\tilde \phi _I^i\). A homogeneous densitized triad can be written in a dual form with coefficients \(\tilde p_i^I\) conjugate to \(\tilde \phi _I^i\). As in isotropic models, one absorbs powers of the coordinate volume to obtain variables \(\phi_I^i\) and \(p_i^I\).

The kinematics is the same for all class A models, except possibly for slight differences in the diffeomorphism constraint [25, 36]. Connection components define a distinguished triple of su(2) elements \(\tilde \phi _I^i{\tau _i}\), one for each independent direction of space. Holonomies in those directions are then obtained as \(h_I^{({\mu _I})} = \exp ({\mu _I}\phi _I^i{\tau _i}) \in {\rm{SU(2)}}\) with parameters μI for the edge lengths. Cylindrical functions depend on those holonomies, i.e., are countable superpositions of terms \(f(h_1^{({\mu _1})},h_2^{({\mu _2})},h_3^{({\mu _3})})\). A basis can be written down as spin network states

$$f(h_1^{({\mu _1})},h_2^{({\mu _2})},h_3^{(\mu 3)}) = {\rho _{j1}}(h_1^{({\mu _1})})_{{B_1}}^{{A_1}}{\rho _{j2}}(h_2^{({\mu _2})})_{{B_2}}^{{A_2}}{\rho _{j3}}(h_3^{({\mu _3})})_{{B_3}}^{{A_3}}K_{{A_1}{A_2}{A_3}}^{{B_1}{B_2}{B_3}}$$

where the matrix K specifies how the representation matrices are contracted to a gauge invariant function of \(\phi _I^i\). There are uncountably many such states for different μI and thus the Hilbert space is non-separable. In contrast to isotropic models, the general homogeneous theory is genuinely SU(2) and therefore not much simpler than the full theory for individual calculations.

As a consequence of homogeneity we observe the same degeneracy as in isotropic models where both spin and edge length appear similarly as parameters. Spins are important to specify the contraction K and thus appear, e.g., in the volume spectrum. For this one needs to know the spins, and it is not sufficient to consider only products jIδI. On the other hand, there is still a degeneracy of spin and edge length and keeping both jI and δI independent leaves too many parameters. It is therefore more difficult to determine what the analog of the Bohr compactification is in this case.

5.7 Diagonalization

The situation simplifies if one considers diagonal models, which is usually also done in classical considerations since it does not lead to much loss of information. In a metric formulation, one requires the metric and its time derivative to be diagonal, which is equivalent to a homogeneous densitized triad \(p_i^I = {p^{(I)}}\Lambda _i^I\) and connection \(\phi _I^i = {c_{(I)}}\Lambda _I^i\) with real numbers cI and pI (where coordinate volume has been absorbed as described in Appendix B.1) which are conjugate to each other, \(\{{c_I},{p^J}\} = 8\pi \gamma G\delta _I^J\), and internal directions \(\Lambda _I^i\) as in isotropic models [48]. In fact, the kinematics becomes similar to isotropic models, except that there are now three independent copies. The reason for the simplification is that we are able to separate off the gauge degrees of freedom in \(\Lambda _I^i\) from gauge invariant variables cI and pI (except for remaining discrete gauge transformations changing the signs of two of the pI and cI together). In a general homogenous connection, gauge-dependent and gauge-invariant parameters are mixed together in \(\phi _I^i\), which both react differently to a change in μI. This makes it more difficult to discuss the structure of relevant function spaces without assuming diagonalization.

As mentioned, the variables pI and cI are not completely gauge invariant since a gauge transformation can flip the sign of two components pI and cI while keeping the third fixed. There is thus a discrete gauge group left, and only the total sign sgn(p1p2p3) is gauge invariant in addition to the absolute values.

Quantization can now proceed simply by using as Hilbert space the triple product of the isotropic Hilbert space, given by square integrable functions on the Bohr compactification of the real line. This results in states \(\vert\psi \rangle = \sum\nolimits_{{\mu _1},{\mu _2},{\mu _3}} {{\psi _{{\mu _1},{\mu _2},{\mu _3}}}\vert{\mu _1},{\mu _2},{\mu _3}\rangle}\) expanded in an orthonormal basis

$$\langle {c_1},{c_2},{c_3}\vert {\mu _1},{\mu _2},{\mu _3}\rangle = {e^{i({\mu _1}{c_1} + {\mu _2}{c_2} + {\mu _3}{c_3})/2}}.$$

Gauge invariance under discrete gauge transformations requires \({\psi _{{\mu _1},{\mu _2},{\mu _3}}}\) to be symmetric under a flip of two signs in μI. Without loss of generality one can thus assume that ψ is defined for all real μ3 but only non-negative μ1 and μ2.

Densitized triad components are quantized by

$${\hat{p}^I}\vert {\mu _1},{\mu _2},{\mu _3}\rangle = {1 \over 2}{\mu _{I\gamma}}\ell _{\rm{P}}^2\vert {\mu _1},{\mu _2},{\mu _3}\rangle,$$

which directly give the volume operator \(\hat V = \sqrt {\vert {{\hat p}^1}{{\hat p}^2}{{\hat p}^3}\vert}\) with spectrum

$${V_{{\mu _1},{\mu _2},{\mu _3}}} = ({1 \over 2}\gamma \ell _{\rm{P}}^2)\sqrt {\vert {\mu _1},{\mu _2},{\mu _3}\vert},$$

Moreover, after dividing out the remaining discrete gauge freedom the only independent sign in triad components is given by the orientation \({\mathop{\rm sgn}} ({{\hat p}^1}{{\hat p}^2}{{\hat p}^3})\), which again leads to a doubling of the metric minisuperspace with a degenerate subset in the interior, where one of the pI vanishes.

5.8 Homogeneity: Dynamics

The Hamiltonian constraint can be constructed in the standard manner and its matrix elements can be computed explicitly thanks to the simple volume spectrum. There are holonomy operators for all three directions, and so in the triad representation the constraint equation becomes a partial difference equation for \({\psi _{{\mu _1},{\mu _2},{\mu _3}}}\) in three independent variables. Its (lengthy) form can be found in [48] for the Bianchi I model and in [62] for all other class A models.

Simpler cases arise in so-called locally rotationally symmetric (LRS) models, where a non-trivial isotropy subgroup is assumed. Here, only two independent parameters μ and ν remain, where only one, e.g., ν can take both signs if discrete gauge freedom is fixed, and the vacuum difference equation is, e.g., for Bianchi I,

$$\begin{array}{*{20}c}{\quad 2\delta \sqrt {\vert \nu + 2\delta \vert} ({\psi _{\mu + 2\delta, \nu + 2\delta}} - {\psi _{\mu - 2\delta, \nu + 2\delta}})} \\ {\quad + {1 \over 2}(\sqrt {\vert \nu + \delta \vert} - \sqrt {\vert \nu - \delta \vert})\;((\mu + 4\delta){\psi _{\mu + 4\delta, \nu}} - 2\mu {\psi _{\mu, \nu}} + (\mu - 4\delta){\psi _{\mu - 4\delta, \nu}})} \\ {\quad - 2\delta \sqrt {\vert \nu - 2\delta \vert} ({\psi _{\mu + 2\delta, \nu - 2\delta}} - {\psi _{\mu - 2\delta, \nu - 2\delta}})} \\ {= 0} \\ \end{array}$$
(51)

from the non-symmetric constraint and

$$\begin{array}{*{20}c}{\quad 2\delta (\sqrt {\vert \nu + 2\delta \vert} + \sqrt {\vert \nu \vert})({\psi _{\mu + 2\delta, \nu + 2\delta}} - {\psi _{\mu - 2\delta, \nu + 2\delta}})} \\ {\quad + (\sqrt {\vert \nu + \delta \vert} - \sqrt {\vert \nu - \delta \vert})\;((\mu + 2\delta){\psi _{\mu + 2\delta, \nu}} - \mu {\psi _{\mu, \nu}} + (\mu - 2\delta){\psi _{\mu - 2\delta, \nu}})} \\ {\quad - 2\delta (\sqrt {\vert \nu - 2\delta \vert} + \sqrt {\vert \nu \vert})({\psi _{\mu + 2\delta, \nu - 2\delta}} - {\psi _{\mu - 2\delta, \nu - 2\delta}})} \\ {= 0} \\ \end{array}$$
(52)

from the symmetric version (see also [14]). This leads to a reduction between fully anisotropic and isotropic models with only two independent variables, and provides a class of interesting systems to analyze effects of anisotropies.

5.9 Inhomogeneous models

Homogeneous models provide a rich generalization of isotropic ones, but inhomogeneities lead to stronger qualitative differences. To start with, at least at the kinematical level one has infinitely many degrees of freedom and is thus always dealing with field theories. Studying field theoretical implications does not require going immediately to the full theory since there are many inhomogeneous models of physical interest.

We will describe some 1-dimensional models with one inhomogeneous coordinate x and two others parameterizing symmetry orbits. A general connection is then of the form (with coordinate differentials ωy and ωz depending on the symmetry)

$$A = {A_x}(x){\Lambda _x}(x){\rm{d}}x + {A_y}(x){\Lambda _y}(x){\omega _y} + {A_z}(x){\Lambda _z}(x){\omega _z} + {\rm{field}}\;{\rm{independent}}\;{\rm{terme}}$$
(53)

with three real functions AI(x) and three internal directions AI(x) normalized to \({\rm{tr(}}\Lambda _1^2{\rm{) = -}}{1 \over 2}\), which in general are independent of each other. The situation in a given point x is thus similar to general homogeneous models with nine free parameters. Correspondingly, there are not many simplifications from this general form, and one needs analogs of the diagonalization employed for homogeneous models. What is required mathematically for simplifications to occur is a connection with internally perpendicular components, i.e., \({\rm{tr(}}{\Lambda _I}{\Lambda _J}{\rm{) = -}}{1 \over 2}{\delta _{IJ}}\) in each point. This arises in different physical situations.

5.10 Einstein-Rosen waves

One class of 1-dimensional models is given by cylindrically symmetric gravitational waves, with connections and triads

$$A = {A_x}(x){\tau _3}{\rm{d}}x + ({A_1}(x){\tau _1} + {A_2}(x){\tau _2}){\rm{d}}z + ({A_3}(x){\tau _1} + {A_4}(x){\tau _2}){\rm{d}}\varphi,$$
(54)
$$E = {E^x}(x){\tau _3}{\partial \over {\partial x}} + ({E^1}(x){\tau _1} + {E^2}(x){\tau _2}){\partial \over {\partial z}} + ({E^3}(x){\tau _1} + {E^4}(x){\tau _2}){\partial \over {\partial \varphi}}$$
(55)

in cylindrical coordinates. This form is more special than (53), but still not simple enough for arbitrary A1, A2, A3 and A4. Einstein-Rosen waves [109, 34] are a special example of cylindrical waves subject to the polarization condition A2A4 + A1A3 = 0, and analogously for triad components. This is just what is needed to restrict the model to internally perpendicular connection components and is thus analogous to diagonalization in a homogeneous model.

5.10.1 Canonical variables

A difference to homogeneous models, however, is that the internal directions of a connection and a triad do not need to be identical, which in homogeneous models with internal directions \(\Lambda _I^i\) is the case as a consequence of the Gauss constraint \({\epsilon ^{ijk}}\phi _I^jp_k^I = 0\). With inhomogeneous fields, now, the Gauss constraint reads

$${E^{x^{\prime}}} + {A_1}{E^2} - {A_2}{E^1} + {A_3}{E^4} - {A_4}{E^3} = 0$$
(56)

or, after splitting off norms and internal directions

$${A_z}: = \sqrt {A_1^2 + A_2^2}, \quad {A_\varphi}: = \sqrt {A_3^2 + A_4^2}$$
(57)
$$\Lambda _z^A: = {{{A_1}{\tau _1} + {A_2}{\tau _2}} \over {{A_z}}},\quad \Lambda _\varphi ^A: = {{{A_3}{\tau _1} + {A_4}{\tau _2}} \over {{A_\varphi}}}$$
(58)

and analogously \({E^z},{E^\varphi},\Lambda _E^z\) and \(\Lambda _E^\varphi\),

$${E^{x^{\prime}}} + ({A_z}{E^z} + {A_\varphi}{E^\varphi})\sin \alpha = 0$$
(59)

with \(\sin \alpha : = - 2{\rm{tr}}(\Lambda _z^A\Lambda _E^z{\tau _3})\). If Ex is not constant, α cannot be zero and thus connections and triads have different internal directions.

As a consequence, Ez is not conjugate to Az, anymore, and instead the momentum of Az is Ez cos α [52]. This seems to make a quantization more complicated since the momenta will be quantized to simple flux operators, but do not directly determine the geometry such as the volume \(V = 4\pi \int {\rm{d}}x\sqrt {\vert {E^x}{E^z}{E^\varphi}\vert}\). For this, one would need to know the angle α which depends on both connections and triads. Moreover, it would not be obvious how to obtain a discrete volume spectrum since then volume does not depend only on fluxes.

It turns out that there is a simple canonical transformation, which allows one to work with canonical variables Ez and Eϕ playing the role of momenta of Azcosα and Aϕ cosα[75]. This seems to be undesirable, too, since now the connection variables are modified which play an important role for holonomies. That these canonical variables are very natural, however, follows after one considers the structure of spin connections and extrinsic curvature tensors in this model. The new canonical variables are then simply given by Azcosα = γKz, Aϕcosα = γKϕ, i.e., proportional to extrinsic curvature components. Thus, in the inhomogeneous model we simply replace connection components with extrinsic curvature in homogeneous directions (note that Ax remains unchanged) while momenta remain elementary triad components. This is part of a broader scheme which is also important for the Hamiltonian constraint operator (Section 5.14).

5.10.2 Representation

With the polarization condition the kinematics of the quantum theory simplifies. Relevant holonomies are given by \({h_e}(A) = \exp \left({{1 \over 2}i\int\nolimits_e {{A_x}} (x){\rm{d}}x} \right)\) along edges in the 1-dimensional manifold and

$$h_v^z(A) = \exp (i\gamma {\nu _v}{K_z}(v)),\quad h_v^\varphi (A) = \exp (i\gamma {\mu _v}{K_\varphi}(v))$$

in vertices v with real μv, νv ≥ 0. Cylindrical functions depend on finitely many of those holonomies, whose edges and vertices form a graph in the 1-dimensional manifold. Flux operators, i.e., quantized triad components, act simply by

$${\hat{E}^x}(x){T_{g,k,u}} = {{\gamma \ell _{\rm{P}}^2} \over {8\pi}}{{{k_{{e^ +}(x)}} + {k_{{e^ -}(x)}}} \over 2}{T_{g,k,u}}$$
(60)
$$\int\nolimits_{\mathcal I} {{{\hat{E}}^z}} {T_{g,k,u}} = {{\gamma \ell _{\rm{P}}^2} \over {4\pi}}\sum\limits_{v \in {\mathcal I}} {{\nu _v}} {T_{g,k,u}}$$
(61)
$$\int\nolimits_{\mathcal I} {{{\hat E}^\varphi}} {T_{g,k,u}} = {{\gamma \ell _{\rm{P}}^2} \over {4\pi}}\sum\limits_{v \in {\mathcal I}} {{\mu _v}} {T_{g,k,u}}$$
(62)

on a spin network state

$$\begin{array}{*{20}c}{{T_{g,k,u}}(A) = \prod\limits_{e \in g} {{\rho _{{k_e}}}} ({h_e})\prod\limits_{v \in V(g)} {{\rho _{{\mu _v}}}} (\gamma {K_\varphi}(v)){\rho _{{\nu _v}}}(\gamma {K_z}(v)){\rho _{{k_v}}}(\beta (v))} \\ {\quad \quad \quad \;\; = \prod\limits_{e \in g} {\exp \left({{1 \over 2}i{k_e}\begin{array}{*{20}c}\int \\ e \\ \end{array} {A_x}(x){\rm{d}}x} \right)} \prod\limits_{v \in V(g)} {{e^{i\gamma {\mu _v}{K_\varphi}(v)}}{e^{i\gamma {\mu _v}{K_z}(v)}}{e^{i{k_v}\beta (v)}}},} \\ \end{array}$$
(63)

which also depend on the gauge angle β determining the internal direction of \(\Lambda _z^E\). If we solve the Gauss constraint at the quantum level, the labels kv will be such that a gauge invariant spin network only depends on the gauge invariant combination Ax + β′.

Since triad components have simple quantizations, one can directly combine them to get the volume operator and its spectrum

$${V_{k,\mu, \nu}} = {{{\gamma ^{3/2}}\ell _{\rm{P}}^3} \over {4\sqrt \pi}}\sum\limits_v {\sqrt {{\mu _v}{\nu _v}\vert {k_{{e^ +}(v)}} + {k_{{e^ +}(v)}}\vert}.}$$
(64)

The labels μv and νv are always non-negative, and the local orientation is given through the sign of edge labels ke.

Commutators between holonomies and the volume operator will technically be similar to homogeneous models, except that there are more possibilities to combine different edges. Accordingly, one can easily compute all matrix elements of composite operators such as the Hamiltonian constraint. The result is only more cumbersome because there are more terms to keep track of. Again as in diagonal homogeneous cases, the triad representation exists and one can formulate the constraint equation there. Now, however, one has infinitely many coupled difference equations for the wave function since the lapse function is inhomogeneous (one obtains one difference equation for each vertex).

There are obvious differences to cases considered previously owing to inhomogeneity. For instance, each edge label can take positive or negative values, or go through zero during evolution corresponding to the fact that a spatial slice does not need to intersect the classical singularity everywhere. Also the structure of coefficients of the difference equations, though qualitatively similar to homogeneous models, is changed crucially in inhomogeneous models, mainly due to the volume eigenvalues (64). Now, \({k_{e +}}\), say, and thus Ex can be zero without volume eigenvalues in neighboring vertices having zero volume.

5.11 Spherical symmetry

For spherically symmetric models, a connection has the form (Appendix B.3)

$$A = {A_x}(x){\tau _3}{\rm{d}}r + ({A_1}(x){\tau _2}){\rm{d}}\vartheta + ({A_1}(x){\tau _2} - {A_2}(x){\tau _1})\sin \vartheta {\rm{d}}\varphi + {\tau _3}\cos \vartheta {\rm{d}}\varphi$$
(65)

whose field-dependent terms automatically have perpendicular internal directions. In this case, it is not diagonalization as in the polarization condition for Einstein-Rosen waves but a non-trivial isotropy subgroup which leads to this property. The kinematical quantization is then simplified as discussed before, with the only difference being that there is only one type of vertex holonomy

$${h_v}(A) = \exp (i\gamma {\mu _v}{K_\varphi}(v))$$

as a consequence of a non-trivial isotropy subgroup. The Hamiltonian constraint can again be computed explicitly [75].

Spherically symmetric models are usually used for applications to non-rotating black holes, but they can also be useful for cosmological purposes. They are particularly interesting as models for the evolution of inhomogeneities as perturbations, which can be applied to gravitational collapse but also cosmology. In such a context one often reduces the spherically symmetric configuration even further by requiring a spatial metric

$${\rm{d}}{s^2} = {q_{xx}}(x,t){\rm{d}}{x^2} + {q_{\varphi \varphi}}(x,t){\rm{d}}{\Omega ^2},$$

where qxx is related to qϕϕ by \({\partial _x}\sqrt {{q_{xx}}} = \sqrt {{q_{\varphi \varphi}}}\). One example for such a metric is the spatial part of a flat Friedmann-Robertson-Walker space-time, where qϕϕ(x, t) = x2a(t)2. This allows one to study perturbations around a homogeneous space-time, which can also be done at the quantum level.

5.12 Loop inspired quantum cosmology

The constructions described so far in this section follow all the steps in the full theory as closely as possible. Most importantly, one obtains quantum representations inequivalent to those used in a Wheeler-DeWitt quantization, which results in many further implications. This has inspired investigations where not all the steps of loop quantum gravity are followed, but only the same type of representation, i.e., the Bohr Hilbert space in an isotropic model, is used. Other constructions, based on ADM rather than Ashtekar variables, are then done in the most straightforward way rather than a way suggested by the full theory [131].

In isotropic models the results are similar, but already here one can see conceptual differences. Since the model is based on ADM variables, in particular using the metric and not triads, it is not clear what the additional sign factor sgn(μ), which is then introduced by hand, means geometrically. In loop quantum cosmology it arose naturally as orientation of triads, even before its role in removing the classical singularity, to be discussed in Section 5.15, had been noticed. (The necessity of having both signs available is also reinforced independently by kinematical consistency considerations in the full theory [117].) In homogeneous models the situation is even more complicated since sign factors are still introduced by hand, but not all of them are removed by discrete gauge transformations as in Section 5.7 (see [158] as opposed to [14]). Those models are useful to illuminate possible effects, but they also demonstrate how new ambiguities, even with conceptual implications, arise if guidance from a full theory is lost.

In particular the internal time dynamics is more ambiguous in those models and thus not usually considered. There are then only arguments that the singularity could be avoided through boundedness of relevant operators, but those statements are not generic in anisotropic models [62] or even the full theory [85]. Moreover, even if all curvature quantities could be shown to be bounded, the evolution could still stop (as happens classically where not any singularity is also a curvature singularity).

5.13 Dynamics

Because irrational numbers are always the result of calculations, never the result of direct measurement, might it not be possible in physics to abandon irrational numbers altogether and work only with the rational numbers? That is certainly possible, but it would be a revolutionary change.…

At some future time, when much more is known about space and time and the other magnitudes of physics, we may find that all of them are discrete.

Rudolf Carnap

An Introduction to the Philosophy of Science

So far we have mainly described the kinematical construction of symmetric models in loop quantum gravity up to the point where the Hamiltonian constraint appears. Since many dynamical issues in different models appear in a similar fashion, we discuss them in this section with a common background. The main feature is that dynamics is formulated by a difference equation that by itself, compared to the usual appearance of differential equations, implies new properties of evolution. Depending on the model there are different classes, which even within a given model are subject to quantization choices. Yet, since there is a common construction procedure many characteristic features are very general.

Classically, curvature encodes the dynamics of geometry and does so in quantum gravity, too. On the other hand, quantum geometry is most intuitively understood in eigenstates of geometry, e.g., a triad representation if it exists, in which curvature is unsharp. Anyway, only solutions to the Hamiltonian constraint are relevant, which in general are peaked neither on spatial geometry nor on extrinsic curvature. The role of curvature thus has a different, less direct meaning in quantum gravity. Still, it is instructive to quantize classical expressions for curvature in special situations, such as a−2 in isotropy. Since the resulting operator is bounded, it has played an influential role on the development of statements regarding the fate of classical singularities.

However, one has to keep in mind that isotropy is a very special case, as emphasized before, and already anisotropic models shed quite a different light on curvature quantities. Isotropy is special because there is only one classical spatial length scale given by the scale factor a, such that intrinsic curvature can only be a negative power such as a−2 just for dimensional reasons. That the modification is not obvious by quantization in the model is illustrated by comparing the intrinsic curvature term ka−2, which remains unmodified and thus unbounded in the purely isotropic quantization, with the term coming from a matter Hamiltonian where the classical divergence of a−3 is cut off.

In an anisotropic model we do have different classical scales and thus dimensionally also terms like \({a_1}a_2^{- 3}\) are possible. It is then not automatic that the quantization is bounded even if \(a_2^{- 3}\) were to be bounded. As an example for such quantities consider the spatial curvature scalar given by W(p1, p2, p3)/(p1p2p3) with W in (37) through the spin connection components. When quantized and then reduced to isotropy, one does obtain a cut-off to the intrinsic curvature term ka−2 as mentioned in Section 4.9, but the anisotropic expression remains unbounded on minisuperspace. The limit to vanishing triad components is direction dependent and the isotropic case picks out a vanishing limit. However, in general this is not the limit taken by the dynamical trajectories. Similarly, in the full theory one can show that inverse volume operators are not bounded even, in contrast to anisotropic models, on states where the volume eigenvalue vanishes [85]. However, this is difficult to interpret since nothing is known about its relevance for dynamics, and even the geometrical role of spin labels, and thus of the configurations considered, is unclear.

It is then quantum dynamics that is necessary to see what properties are relevant and how degenerate configurations are approached. This should allow one to check if the classical boundary a finite distance away is removed in quantum gravity. This can only happen if quantum gravity provides candidates for a region beyond the classical singularity, and means to probe how to evolve there. The most crucial aim is to prevent incompleteness of space-time solutions or their quantum replacements. Even if curvature would be finite, by itself it would not be enough since one could not tell if the singularity persists as incompleteness. Only a demonstration of continuing evolution can ultimately show that singularities are absent.

5.14 Dynamics: General construction

Not all steps in the construction of the full constraint can be taken over immediately to a model since symmetry requirements have to be respected. It is thus important to have a more general construction scheme that shows how generic different steps are, and whether or not crucial input in a given symmetric situation is needed.

We have already observed one such issue, which is the appearance of holonomies but also simple exponentials of connection components without integration. This is a consequence of different transformation properties of different connection components in a reduced context. Components along remaining inhomogeneous directions, such as Ax for Einstein-Rosen waves, play the role of connection components in the model, giving rise to ordinary holonomies. Other components, such as Az and Aϕ in Einstein-Rosen waves or all components in homogeneous models, transform as scalars and thus only appear in exponentials without integration. In the overall picture, we have the full theory with only holonomies, homogeneous models with only exponentials, and inhomogeneous models in between where both holonomies and exponentials appear.

Another crucial issue is that of intrinsic curvature encoded in the spin connection. In the full theory, the spin connection does not have any covariant meaning and in fact can locally be made to vanish. In symmetric models, however, some spin connection components can become covariantly well-defined since not all coordinate transformations are allowed within a model. In isotropic models, for instance, the spin connection is simply given by a constant proportional to the curvature parameter. Of particular importance is the spin connection when one considers semiclassical regimes because intrinsic curvature does not need to become small there in contrast to extrinsic curvature. Since the Ashtekar connection mixes the spin connection and extrinsic curvature, its semiclassical properties can be rather complicated in symmetric models.

The full constraint is based on holonomies around closed loops in order to approximate Ashtekar curvature components when the loop becomes small in a continuum limit. For homogeneous directions, however, one cannot shrink the loop and instead works with exponentials of the components. One thus approximates the classical components only when arguments of the exponential are small. If these arguments were always connection components, one would not obtain the right semiclassical properties because those components can remain large. In models one thus has to base the construction for homogeneous directions on extrinsic curvature components, i.e., subtract off the spin connection from the Ashtekar connection. For inhomogeneous directions, on the other hand, this is not possible since one needs a connection in order to define a holonomy.

At first sight this procedure seems rather ad hoc and even goes half a step back to ADM variables since extrinsic curvature components are used. However, there are several places where this procedure turns out to be necessary for a variety of independent reasons. We have already seen in Section 5.10 that inhomogeneous models can lead to a complicated volume operator when one insists on using all Ashtekar connection components. When one allows for extrinsic curvature components in the way just described, on the other hand, the volume operator becomes straightforward. This appeared after performing a canonical transformation, which rests non-trivially on the form of inhomogeneous spin connections and extrinsic curvature tensors.

Moreover, in addition to the semiclassical limit used above as justification one also has to discuss local stability of the resulting evolution equation [59]: Since higher order difference equations have additional solutions, one must ensure that they do not become dominant in order not to spoil the continuum limit. This is satisfied with the above construction, while it is generically violated if one were to use only connection components.

There is thus a common construction scheme available based on holonomies and exponentials. As already discussed, this is responsible for correction terms in a continuum limit, but also gives rise to the constraint equation being a difference equation in a triad representation, whenever it exists. In homogeneous models the structure of the resulting difference equation is clear, but there are different open possibilities in inhomogeneous models. This is intimately related to the issue of anomalies, which also appears only in inhomogeneous models.

With a fixed choice, one has to solve a set of coupled difference equations for a wave function on superspace. The basic question then always is what kind of initial or boundary value problem has to be used in order to ensure the existence of solutions with suitable properties, e.g., in a semiclassical regime. Once this is specified one can already discuss the singularity problem since one needs to find out if initial conditions in one semiclassical regime together with boundary conditions away from classical singularities suffice for a unique solution on all of superspace. A secondary question is how this equation can be interpreted as evolution equation for the wave function in an internal time. This is not strictly necessary and can be complicated owing to the problem of time in general. Nevertheless, when available, an evolution interpretation can be helpful for interpretations.

5.15 Singularities

Il n’est rien de plus précieux que le temps, puisque c’est le prix de l’éternité.

(There is nothing more precious than time, for it is the price of eternity.)

Louis Bourdaloue

Sermon sur la perte de temps

In the classical situation, we always have trajectories on superspace running into singular submanifolds where some or all densitized triad components vanish. In semiclassical regimes one can think of physical solutions as wave packets following these trajectories in internal time, but at smaller triad components spreading and deformations from a Gaussian become stronger. Moreover, discreteness becomes essential and properties of difference equations need to be taken into account in order to see what is happening at the singular submanifolds.

The simplest situation is given by isotropic models where superspace is one dimensional with coordinate p. Minisuperspace is thus disconnected classically with two sides separated by the classical singularity at p = 0. At this point, classical energy densities diverge and there is no well-defined initial value problem to evolve further. (Sometimes, formal extensions of solutions beyond a classical singularity exist [122], but they are never unique and unrelated to the solution preceding the singularity. This shows that a resolution of singularities has not only to provide a new region, but also an evolution there uniquely from initial values at one side.) A Wheeler-DeWitt quantization would similarly lead to diverging matter Hamiltonian operators and the initial value problem for the wave function generically breaks down. In isotropic loop quantum cosmology we have already seen that the matter Hamiltonian does not have diverging contributions from inverse metric components even at the classical singularity. Nevertheless, the evolution could break down if highest order coefficients in the difference equation become zero. This indeed happens with the non-symmetric constraint (49) or (51), but in these cases can be seen not to lead to any problems: some coefficients can become zero such that the wave function at μ = 0 remains undetermined by initial conditions, but the wave function at the other side of the classical singularity is still determined uniquely. There is no breakdown of evolution, and thus no singularity [39]. As one can see, this relies on crucial properties of the loop representation with well-defined inverse metric components and a difference rather than differential equation [43].

Also the structure of difference equations is important, depending on some choices. Most important is the factor ordering or symmetrization chosen. As just discussed, the ordering used earlier leads to non-singular evolution but with the wave function at the classical singularity itself remaining undetermined. In anisotropic models one can symmetrize the constraint and obtain a difference equation, such as (52), whose leading order coefficients never vanish. Evolution then never stops and even the wave function at the classical singularity is determined. In the isotropic case, direct symmetrization would lead to a break-down of evolution, which thus provides an example for singular quantum evolution and demonstrates the non-triviality of continuing evolution: The leading order coefficient would then be Vμ−3δVμ−5δ + Vμ+δVμ−δ, which vanishes if and only if μ = 2δ. Thus, in the backward evolution ψ−2δ remains undetermined, just as ψ0 is undetermined in the non-symmetric ordering. However, now ψ−δ would be needed to evolve further. Since it is not determined by initial data, one would need to prescribe this value, or else the evolution stops. There is thus a new region at negative μ, but evolution does not continue uniquely between the two sides. In such a case, even though curvature is bounded, the quantum system would be singular. Similar behavior happens in other orderings such as when triads are ordered to the left. Note that also in the full theory one cannot order triads to the left since otherwise the constraint would not be densely defined [194].

The breakdown of the symmetric ordering in isotropic models is special and related to the fact that all directions degenerate. The breakdown does not happen for a symmetric ordering in anisotropic or even inhomogeneous systems. One can avert it in isotropic cases by multiplying the constraint with sgn p before symmetrizing, so that the additional factor of sgn μ leads to non-zero coefficients as in (50).

This is the general scheme, which also applies in more complicated cases. The prime example for the general homogeneous behavior is given by the Kasner evolution of the Bianchi I model. Here, the approach to the singularity is not isotropic but given in such a way that two of the three diagonal metric components become zero while the third one diverges. This would lead to a different picture than described before since the classical singularity then lies at the infinite boundary of metric or co-triad minisuperspace. Also unlike in the isotropic case, densities or curvature potentials are not necessarily bounded in general as functions on minisuperspace, and the classical dynamical approach is important. In densitized triad variables, however, we have a situation as before since here all components approach zero, although at different rates. Now the classical singularity is in the interior of minisuperspace and one can study the evolution as before, again right through the classical singularity. Note that densitized triad variables were required for a background independent quantization, and now independently for non-singular evolution.

Other homogeneous models are more complicated since for them Kasner motion takes place with a potential given by curvature components. Approximate Kasner epochs arise when the potential is negligible, intermitted by reflections at the potential walls where the direction of Kasner motion in the anisotropy plane changes. Still, since in each Kasner epoch the densitized triad components decrease, the classical singularity remains in the interior and is penetrated by the discrete quantum evolution.

One can use this for indications as to the general inhomogeneous behavior by making use of the BKL scenario. If this can be justified, in each spatial point the evolution of geometry is given by a homogeneous model. For the quantum formulation this indicates that also here classical singularities are removed. However, it is by no means clear whether the BKL scenario applies at the quantum level since even classically it is not generally established. If the scenario is not realized (or if some matter systems can change the local behavior), diverging p are possible and the behavior would qualitatively be very different. One thus has to study the inhomogeneous quantum evolution directly as done before for homogeneous cases.

In the 1-dimensional models described here classical singularities arise when Ex becomes zero. Since this is now a field, it depends on the point x on the spatial manifold where the slice hits the classical singularity. At each such place, midisuperspace opens up to a new region not reached by the classical evolution, where the sign of Ex(x) changes and thus the local orientation of the triad. Again, the kinematics automatically provides us with these new regions just as needed, and quantum evolution continues. Also, the scheme is realized much more non-trivially, and now even the non-symmetric ordering is not allowed. This is a consequence of the fact that ke = 0 for a single edge label does not imply that neighboring volume eigenvalues vanish. There is thus no obvious decoupling in a non-singular manner, and it shows how less symmetric situations put more stringent restrictions on the allowed dynamics. Still, the availability of other possibilities, maybe with leading coefficients which can vanish and result in decoupling, needs to be analyzed. Most importantly, the symmetric version still leads to non-singular evolution even in those inhomogeneous cases which have local gravitational degrees of freedom [57].

There is thus a general scheme for the removal of singularities: In the classical situation, one has singular boundaries of superspace which cannot be penetrated. Densitized triad variables then lead to new regions, given by a change in the orientation factor sgn det E which, however, does not help classically since singularities remain as interior boundaries. For the quantum situation one has to look at the constraint equation and see whether or not it uniquely allows to continue a wave function to the other side (which does not require time parameters even though they may be helpful if available). This usually depends on factor ordering and other choices that arise in the construction of constraint operators and play a role also for the anomaly issue. One can thus fix ambiguities by selecting a non-singular constraint if possible. However, the existence of nonsingular versions, as realized in a natural fashion in homogeneous models, is a highly non-trivial and by no means automatic property of the theory showing its overall consistency.

In inhomogeneous models the issue is more complicated. We thus have a situation where the theory, which so far is well-defined, can be tested by trying to extend results to more general cases. It should also be noted that different models should not require different quantization choices unless symmetry itself is clearly responsible (as happens with the orientation factor in the symmetric ordering for an isotropic model, or when non-zero spin connection components receive covariant meaning in models), but that there should rather be a common scheme leading to nonsingular behavior. This puts further strong conditions on the construction, and is possible only if one knows how models and the full theory are related.

5.16 Initial/boundary value problems

In isotropic models the gravitational part of the constraint corresponds to an ordinary difference operator which can be interpreted as generating evolution in internal time. One thus needs to specify only initial conditions to solve the equation. The number of conditions is large since, first, the procedure to construct the constraint operator usually results in higher order equations and, second, this equation relates values of a wave function ψμ defined on an uncountable set. In general, one thus has to choose a function on a real interval unless further conditions are used.

This can be achieved, for instance, by using observables that can reduce the kinematical framework back to wave functions defined on a countable discrete lattice [202]. Similar restrictions can come from semiclassical properties or the physical inner product [162], all of which has not yet been studied in generality.

The situation in homogeneous models is similar, but now one has several gravitational degrees of freedom only one of which is interpreted as internal time. One has a partial difference equation for a wave function on a minisuperspace with boundary, and initial as well as boundary conditions are required [48]. Boundary conditions are imposed only at non-singular parts of minisuperspace such as μ = 0 in LRS models (51). They must not be imposed at places of classical singularities, of course, where instead the evolution must continue just as at any regular part.

In inhomogeneous models, then, there are not only many independent kinematical variables but also many difference equations for only one wave function on midisuperspace. These difference equations are of a similar type as in homogeneous models, but they are coupled in complicated ways. Since one has several choices in the general construction of the constraint, there are different possibilities for the way how difference equations arise and are coupled. Not all of them are expected to be consistent, i.e., in many cases some of the difference equations will not be compatible such that there would be no non-zero solution at all. This is related to the anomaly issue since the commutation behavior of difference operators is important for properties and the existence of common solutions.

So far, the evolution operator in inhomogeneous models has not been studied in detail, and solutions in this case remain poorly understood. The difficulty of this issue can be illustrated by the expectations in spherical symmetry where there is only one classical physical degree of freedom. If this is to be reproduced for semiclassical solutions of the quantum constraint, there must be a subtle elimination of infinitely many kinematical degrees of freedom such that in the end only one physical degree of freedom remains. Thus, from the many parameters needed in general to specify a solution to a set of difference equations, only one can remain when compatibility relations between the coupled difference equations and semiclassicality conditions are taken into account.

How much this cancellation depends on semiclassicality and asymptotic infinity conditions remains to be seen. Some influence is to be expected since classical behavior should have a bearing on the correct reproduction of classical degrees of freedom. However, it may also turn out that the number of solutions to the quantum constraint is more sensitive to quantum effects. It is already known from isotropic models that the constraint equation can imply additional conditions for solutions beyond the higher order difference equation, as we will discuss in Section 5.18. This usually arises at the place of classical singularities where the order of the difference equation can change. Since the quantum behavior at classical singularities is important here, the number of solutions can be different from the classically expected freedom, even when combined with possible semiclassical requirements far away from the singularity. We will now first discuss these requirements in semiclassical regimes, followed by more information on possibly arising additional conditions for solutions.

5.17 Pre-classicality and boundedness

The high order of difference equations implies that there are in general many independent solutions, most of which are oscillating on small scales, i.e., when the labels change only slightly. One possibility to restrict the number of solutions then is to require suppressed or even absent oscillations on small scales [40]. Intuitively, this seems to be a pre-requisite for semiclassical behavior and has thus been called pre-classicality. It can be motivated by the fact that a semiclassical solution should not be sensitive to small changes of, e.g., the volume by amounts of Planck size. However, even though the criterion sounds intuitively reasonable, there is so far no justification through more physical arguments involving observables or measurement processes to extract information from wave functions. The status of pre-classicality as a selection criterion is thus not final.

Moreover, pre-classicality is not always consistent in all disjoint classical regimes or with other conditions. For instance, as discussed in the following section, there can be additional conditions on wave functions arising from the constraint equation at the classical singularity. Such conditions do not arise in classical regimes, but they nevertheless have implications for the behavior of wave functions there through the evolution equation [87, 86]. Pre-classicality also may not be possible to impose in all disconnected classical regimes. If the evolution equation is locally stable — which is a basic criterion for constructing the constraint — choosing initial values in classical regimes, which do not have small-scale oscillations, guarantees that oscillations do not build up through evolution in a classical regime [59]. However, when the solution is extended through the quantum regime around a classical singularity, oscillations do arise and do not in general decay after a new supposedly classical regime beyond the singularity is entered. It is thus not obvious that indeed a new semiclassical region forms even if the quantum evolution for the wave function is nonsingular. On the other hand, evolution does continue to large volume and macroscopic regions, which is different from other scenarios such as [124] where inhomogeneities have been quantized on a background.

A similar issue is the boundedness of solutions, which also is motivated intuitively by referring to the common probability interpretation of quantum mechanics [119] but must be supported by an analysis of physical inner products. The issue arises in particular in classically forbidden regions where one expects exponentially growing and decaying solutions. If a classically forbidden region extends to infinite volume, as happens for models of recollapsing universes, the probability interpretation would require that only the exponentially decaying solution is realized. As before, such a condition at large volume is in general not consistent in all asymptotic regions or with other conditions arising in quantum regimes.

Both issues, pre-classicality and boundedness, seem to be reasonable, but their physical significance has to be founded on properties of the physical inner product. They are rather straightforward to analyze in isotropic models without matter fields, where one is dealing with ordinary difference equations. However, other cases can be much more complicated such that conclusions drawn from isotropic models alone can be misleading. Moreover, numerical investigations have to be taken with care since in particular for boundedness an exponentially increasing contribution can easily arise from numerical errors and dominate the exact, potentially bounded solution.

One thus needs analytical or at least semi-analytical techniques to deal with these issues. For pre-classicality one can advantageously use generating function techniques [87] if the difference equation is of a suitable form, e.g., has only coefficients with integer powers of the discrete parameter. The generating function \(G(x): = {\sum\nolimits_n}{\psi _n}{x^n}\) for a solution ψn on an equidistant lattice then solves a differential equation equivalent to the difference equation for ψn. If G(x) is known, one can use its pole structure to get hints for the degree of oscillations in ψn. In particular the behavior around x = −1 is of interest to rule out alternating behavior where ψn is of the form ψn = (−1)ξn with ξn > 0 for all n (or at least all n larger than a certain value). At x = −1 we then have \(G(- 1) = {\sum\nolimits_n}{\xi _n}\), which is less convergent than the value for a non-alternating solution ξn =ξn resulting in \(G(- 1) = {\sum\nolimits_n}{(- 1)^n}{\xi _n}\). One can similarly find conditions for the pole structure to guarantee boundedness of ψn, but the power of the method depends on the form of the difference equation. More general techniques are available for the boundedness issue, and also for alternating behavior, by mapping the difference equation to a continued fraction which can be evaluated analytically or numerically [71]. One can then systematically find initial values for solutions that are guaranteed to be bounded.

5.18 Dynamical initial conditions

I am Aton when I am alone in the Nun, but I am Re when he appears, in the moment when he starts to govern what he has created.

Book of the Dead

The traditional subject of quantum cosmology is the imposition of initial conditions for the wave function of a universe in order to guarantee its uniqueness. In the Wheeler-DeWitt framework this is done at the singularity a = 0, sometimes combined with final conditions in the classical regime. One usually uses intuitive pictures as guidance, akin to Lemaitre’s primitive atom whose decay is supposed to have created the world, Tryon’s and Vilenkin’s tunneling event from nothing, or the closure of space-time into a Euclidean domain by Hartle and Hawking.

In the latter approaches, which have been formulated as initial conditions for solutions of the Wheeler-DeWitt equation [204, 121], the singularity is still present at a = 0, but re-interpreted as a meaningful physical event through the conditions. In particular, the wave function is still supported at the classical singularity, i.e., ψ(0) ≠ 0, in contrast to DeWitt’s original idea of requiring ψ(0) = 0 as a means to argue for the absence of singularities in quantum gravity [104]. DeWitt’s initial condition is in fact, though most appealing conceptually, not feasible in general since it does not lead to a well-posed initial value formulation in more complicated models: the only solution would then vanish identically. Zeh tried to circumvent this problem, for instance by proposing an ad hoc Planck potential, which is noticeable only at the Planck scale and makes the initial problem well-defined [89]. However, the problem remains that in general there is no satisfying origin of initial values.

In all these ideas, the usual picture in physics has been taken that there are dynamical laws describing the general behavior of a physical system, and independently initial or boundary conditions to select a particular situation. This is reasonable since usually one can prepare a system, corresponding to choosing initial and boundary values, and then study its behavior as determined by the dynamical laws. For cosmology, however, this is not appropriate since there is no way to prepare the universe.

At this point, there is a new possibility opened up by loop quantum cosmology where the dynamical law and initial conditions can be part of the same entity [40, 78, 49]. This is a specialty of difference equations whose order can change locally, in contrast to differential equations. Mathematically, such a difference equation would be called singular since its leading order coefficients can become zero. However, physically we have already seen that the behavior is non-singular since the evolution does not break down.

The difference equation follows from the constraint equation, which is the primary object in canonical quantum gravity. As discussed before, it is usually of high order in classical regimes, where the number of solutions can be restricted, e.g., by pre-classicality. But this, at most, brings us to the number of solutions that we have for the Wheeler-DeWitt equation such that one needs additional conditions as in this approach. The new aspect now is that this can follow from the constraint equation itself: since the order of the difference equation can become smaller at the classical singularity, there are less solutions than expected from the semiclassical behavior. In the simplest models, this is just enough to result in a unique solution up to norm, as appropriate for a wave function describing a universe. In those cases, the dynamical initial conditions are comparable to DeWitt’s initial condition, albeit in a manner that is well-posed even in some cases where DeWitt’s condition is not [64].

In general, the issue is not clear but should be seen as a new option presented by the discrete formulation of loop quantum cosmology. Since there can be many conditions to be imposed on wave functions in different regimes, one has to see in each model whether or not suitable nonzero solutions remain at all. In fact, some first investigations indicate that different requirements taken together can be very restrictive [86], which seems to relate well with the non-separability of the kinematical Hilbert space [98]. So far, only homogeneous models have been investigated in detail, but the mechanism of decoupling is known not to be realized in an identical manner in inhomogeneous models.

Inhomogeneous models can add qualitatively new ingredients also to the issue of initial conditions due to the fact that there are many coupled difference equations. There can then be consistency conditions for solutions to the combined system which can strongly restrict the number of independent solutions. This may be welcome, e.g., in spherical symmetry where a single physical parameter remains, but the restriction can easily become too strong even below the classically expected one. Since the consistency between difference equations is related to the anomaly issue, there may be an important role played by quantum anomalies. While classically anomalies should be absent, the quantum situation can be different since it also takes the behavior at the classical singularity into account and is supposed to describe the whole universe. Anomalies can then be precisely what one needs in order to have a unique wave function of a universe even in in-homogeneous cases where initially there is much more freedom. This does not mean that anomalies are simply ignored or taken lightly since it is difficult to arrange having the right balance between many solutions and no non-zero solution at all. However, quantum cosmology suggests that it is worthwhile to have a less restricted, unconventional view on the anomaly issue.

5.19 Summary

There is a general construction of a loop representation in the full theory and its models, which is characterized by compactified connection spaces and discrete triad operators. Strong simplifications of some technical and conceptual steps occur in diverse models. Such a general construction allows a view not only on the simplest case, isotropy, but on essentially all representative systems for gravity.

Most important is the dynamics, which in the models discussed here can be formulated by a difference equation on superspace. A general scheme for a unique extension of wave functions through classical singularities is realized, such that the quantum theory is non-singular. This general argument, which has been verified in many models, is quite powerful since it does not require detailed knowledge of or assumptions about matter. It is independent of the availability of a global internal time, and so the problem of time does not present an obstacle. Moreover, a complicated discussion of quantum observables can be avoided since once it is known that a wave function can be continued uniquely one can extract relational information at both sides of the classical singularity. (If observables would distinguish both sides with their opposite orientations, they would strongly break parity even on large scales in contradiction with classical gravity.) Similarly, information on the physical inner product is not required since there is a general statement for all solutions of the constraint equation. The uniqueness of an extension through the classical singularity thus remains even if some solutions have to be excluded for the physical Hilbert space or factored out if they have zero norm.

This is far from saying that observables or the physical inner product are irrelevant for an understanding of dynamical processes. Such constructions can, fortunately, be avoided for a general statement of non-singular evolution in a wide class of models. For details of the transition and to get information of the precise form of space-time at the other side of classical singularities, however, all those objects are necessary and conceptual problems in their context have to be understood.

So far, the transition has often been visualized by intuitive pictures such as a collapsing universe turning its inside out when orientation is reversed. An hourglass presents a picture for the importance of discrete quantum geometry close to the classical singularity and the emergence of continuous geometry on large scales: Away from the bottleneck of the hourglass, its sand seems to be sinking down almost continuously. Directly at the bottleneck with its small circumference, however, one can see that time measured by the hourglass proceeds in discrete steps — one grain at a time.

The main remaining issue for the mechanism to remove singularities then is the question how the models, where it has been demonstrated, are related to the full theory and to what extent they are characteristic for full quantum geometry.

6 Models within the Full Theory

If he uses a model at all, he is always aware that it pictures only certain aspects of the situation and leaves out other aspects. The total system of physics is no longer required to be such that all parts of its structure can be clearly visualized. …

A physicist must always guard against taking a visual model as more than a pedagogical device or makeshift help. At the same time, he must also be alert to the possibility that a visual model can, and sometimes does, turn out to be literally accurate. Nature sometimes springs such surprises.

Rudolf Carnap

An Introduction to the Philosophy of Science

In the preceding section, the link between models and the full theory was given through the same basic variables and the same kind of representation used, as well as a general construction scheme for the Hamiltonian constraint operator. The desired simplifications were realized thanks to the symmetry conditions, but not too strongly since basic features of the full theory are still recognizable in models. For instance, even though possible in many ways and often made use of, we did not employ special gauges or coordinate or field dependent transformations obscuring the relation. The models are thus as close to the full theory as possible while making full use of simplifications in order to have explicit applications.

Still, there are always some differences not all of which are easy to disentangle. For instance, we have discussed possible degeneracies between spin labels and edge lengths of holonomies, which can arise in the presence of a partial background and lead to new ambiguity parameters not present in the full theory. The question thus arises what the precise relation between models and the full theory is, or even how and to what extent a model for a given symmetry type can be derived from the full theory.

This is possible for the basic representation: The symmetry and the partial background it provides can be used to define natural subalgebras of the full holonomy/flux algebra by using holonomies and fluxes along symmetry generators and averaging in a suitable manner. Since the full representation is unique and cyclic, it induces uniquely a representation of models that is taken directly from the full theory. This will now be described independently for states and basic operators to provide the idea and to demonstrate the role of the extra structure involved. See also [57] and [63] for illustrations in the context of spherical symmetry and anisotropy, respectively.

6.1 Symmetric states

One can imagine to construct states that are invariant under a given action of a symmetry group on space by starting with a general state and naively summing over all its possible translates by elements of the symmetry group. For instance on spin network states, the symmetry group acts by moving the graph underlying the spin network, keeping the labels fixed. Since states with different graphs are orthogonal to each other, the sum over uncountably many different translates cannot be normalizable. In simple cases, such as for graphs with a single edge along a symmetry generator, one can easily make sense of the sum as a distribution. But this is not clear for arbitrary states, in particular for states whose graphs have vertices, which on the other hand would be needed for sufficient generality. A further problem is that any such action of a symmetry group is a subgroup of the diffeomorphism group. At least on compact space manifolds where there are no asymptotic conditions for diffeomorphisms in the gauge group, it then seems that any group averaged diffeomorphism invariant state would already be symmetric with respect to arbitrary symmetries, which is obviously not sensible.

In fact, symmetries and (gauge) diffeomorphisms are conceptually very different, even though mathematically they are both expressed by group actions on a space manifold. Gauge diffeomorphisms are generated by first class constraints of the theory, which in canonical quantum gravity are imposed in the Dirac manner [105] or following refined algebraic quantization [28], conveniently done by group averaging [153]. Symmetries, however, are additional conditions imposed on a given theory to extract a particular sector of special interest. They can also be formulated as constraints added to the theory, but these constraints must be second class for a well defined framework: One obtains a consistent reduced theory, e.g., with a non-degenerate symplectic structure, only if configuration and momentum variables are required to be symmetric in the same (or dual) way.

In the case of gravity in Ashtekar variables, the symmetry type determines, along the lines of Appendix A the form of invariant connections and densitized triads defining the phase space of the reduced model. At the quantum level, however, one cannot keep connections and triads on the same footing since a polarization is required. One usually uses the connection representation in loop quantum gravity such that states are functionals on the space of connections. In a minisuperspace quantization of the classically reduced model states would then be functionals only of invariant connections for the given symmetry type. This suggests to define symmetric states in the full theory to be those states whose support contains invariant connections as a dense subset [66, 38] (one requires only a dense subset because possible generalized connections must be allowed for). As such, they must necessarily be distributional, as already expected from the naive attempt at construction. Symmetric states thus form a subset of the distributional space Cyl*. In this manner, only the reduced degrees of freedom are relevant, i.e., the reduction is complete, and all of them are indeed realized, i.e., the reduction is not too strong. Moreover, an “averaging” map from a non-symmetric state to a symmetric one can easily be defined by restricting the non-symmetric state to the space of invariant connections and requiring it to vanish everywhere else.

This procedure defines states as functionals, but since there is no inner product on the full Cyl* this does not automatically result in a Hilbert space. Appropriately defined subspaces of Cyl*, nevertheless, often carry natural inner products, which is also the case here. In fact, since the reduced space of invariant connections can be treated by the same mathematical techniques as the full space, it carries an analog of the full Ashtekar-Lewandowski measure and this is indeed induced from the unique representation of the full theory. The only difference is that in general an invariant connection is not only determined by a reduced connection but also by scalar fields (see Appendix A). As in the full theory, this space \({{\mathcal A}_{{\rm{inv}}}}\) of reduced connections and scalars is compactified to the space \({{\bar {\mathcal A}}_{{\rm{inv}}}}\) of generalized invariant connections on which the reduced Hilbert space is defined. One thus arrives at the same Hilbert space for the subset of symmetric states in Cyl* as used before for reduced models, e.g., using the Bohr compactification in isotropic models. The new ingredient now is that these states have meaning in the full theory as distributions, whose evaluation on normalizable states depends on the symmetry type and partial background structure used.

That the symmetric Hilbert space obtained in this manner is identical to the reduced loop quantization of Section 5 does not happen by definition but is a result of the procedure. The support of a distribution is by definition a closed subset of the configuration space, and would thus be larger than just the set of generalized invariant connections if \({{\bar {\mathcal A}}_{{\rm{inv}}}}\) would not be a closed subset in \({\bar {\mathcal A}}\). In such a case, the reduction at the quantum level would give rise to more degrees of freedom than a loop quantization of the classically reduced model. As shown in [66], however, the set of invariant connections is a closed subset of the full space of connections such that loop quantum cosmology can be interpreted as a minisuperspace quantization.

6.2 Basic operators

In the classical reduction, symmetry conditions are imposed on both connections and triads, but so far at the level of states only connections have been taken into account. Configuration and momentum variables play different roles in any quantum theory since a polarization is necessary. As we based the construction on the connection representation, symmetric triads have to be implemented at the operator level. (There cannot be additional reduction steps at the state level since, as we already observed, states just implement the right number of reduced degrees of freedom.)

Classically, the reduction of phase space functions is simply done by pull back to the reduced phase space. The flow generated by the reduced functions then necessarily stays in the reduced phase space and defines canonical transformations for the model. An analog statement in the corresponding quantum theory would mean that the reduced state space would be fixed by full operators such that their action (or dual action on distributions) could directly be used in the model without further work. This, however, is not the case with the reduction performed so far. We have considered only connections in the reduction of states, and also classically a reduction to a subspace \({{\mathcal A}_{{\rm{inv}}}} \times \varepsilon\), where connections are invariant but not triads, would be incomplete. First, this would not define a phase space of its own with a non-degenerate symplectic structure. More important in this context is the fact that this subspace would not be preserved by the flow of reduced functions.

As an example (see also [52] for a different discussion in the spherically symmetric model) we consider a diagonal homogeneous model, such as Bianchi I for simplicity, with connections of the form \(A_a^i{\rm{d}}{x^a} = {{\tilde c}_{(I)}}\Lambda _I^i{w^I}\) and look at the flow generated by the full volume \(V = \int {{\rm{d}}^3}x\sqrt {\vert \det \,E\vert}\). It is straightforward to evaluate the Poisson bracket

$$\{A_a^i(x),V\} = 2\pi \gamma G{\epsilon _{abc}}{\epsilon ^{ijk}}E_j^bE_k^c/\sqrt {\vert \det E\vert}$$

already used in (13). A point on \({{\mathcal A}_{{\rm{inv}}}} \times \varepsilon\) characterized by \({{\tilde c}_{(I)}}\Lambda _I^i\) and an arbitrary triad thus changes infinitesimally by

$$\delta ({\tilde{c}_{(I)}}\Lambda _I^I) = 2\pi \gamma G{\epsilon _{Ibc}}{\epsilon ^{ijk}}E_j^bE_k^c/\sqrt {\vert \det E\vert},$$

which does not preserve the invariant form: First, on the right hand side we have arbitrary fields E such that \(\delta ({{\tilde c}_{(I)}}\Lambda _I^i)\) is not homogeneous. Second, even if we would restrict ourselves to homogeneous \(E,\delta ({{\tilde c}_{(I)}}\Lambda _I^i)\) would not be of the original diagonal form. This is the case only if \(\delta ({{\tilde c}_{(I)}}\Lambda _I^i) = \Lambda _I^i\delta ({{\tilde c}_{(I)}})\) since only the \({\tilde c}_{I}\) are canonical variables. The latter condition is satisfied only if

$${\epsilon ^{ijk}}\Lambda _I^j\delta ({\tilde{c}_{(I)}}\Lambda _I^i) = 4\pi \gamma G{\epsilon _{Ibc}}\Lambda _{(I)}^jE_i^bE_j^c/\sqrt {\vert \det E\vert}$$

vanishes, which is not the case in general. This condition is true only if \(E_i^a \propto \Lambda _i^a\), i.e., if we restrict the triads to be of diagonal homogeneous form just as the connections.

A reduction of only one part of the canonical variables is thus incomplete and leads to a situation where most phase space functions generate a flow that does not stay in the reduced space. Analogously, the dual action of full operators on symmetric distributional states does not in general map this space to itself. Thus, an arbitrary full operator maps a symmetric state to a non-symmetric one and cannot be used to define the reduced operator. In general, one needs a second reduction step that implements invariant triads at the level of operators by an appropriate projection of its action back to the symmetric space. This can be quite complicated, and fortunately there are special full operators adapted to the symmetry for which this step is not necessary.

From the above example, it is clear that those operators must be linear in the momenta \(E_i^a\), for otherwise one would have a triad remaining after evaluating the Poisson bracket, which on \({{\mathcal A}_{{\rm{inv}}}} \times \varepsilon\) would not be symmetric everywhere. Fluxes are linear in the momenta, so we can try \({p^K}({z_0}): = {\int\nolimits_{{S_{{z_0}}}}}{{\rm{d}}^2}y\Lambda _{(K)}^kE_k^aw_a^K\) where Sz0 is a surface in the IJ-plane at position z = z0 in the indirection. By choosing a surface along symmetry generators XI and XJ this expression is adapted to the symmetry, even though it is not fully symmetric yet since the position z0 has to be chosen. Again, we compute the Poisson bracket

$$\{A_a^i(x),{p^K}({z_0})\} = 8\pi \gamma G\Lambda _{(K)}^i\int\nolimits_{{S_{{z_0}}}} {\delta (x,y)\omega _a^K(y){{\rm{d}}^{\rm{2}}}y}$$

resulting in

$$\delta ({\tilde{c}_{(I)}}\Lambda _I^i) = 8\pi \gamma G\Lambda _I^i\delta (z,{z_0}).$$

Also here the right hand side is not homogeneous, but we have \({\epsilon ^{ijk}}\Lambda _I^j\delta ({{\tilde c}_{(I)}}\Lambda _I^k) = 0\) such that the diagonal form is preserved. The violation of homogeneity is expected since the flux is not homogeneous. This can easily be remedied by “averaging” the flux in the K-direction to

$${p^K}: = \mathop {\lim}\limits_{N \to \infty} {N^{- 1}}\sum\limits_{\alpha = 1}^N {{p^K}(\alpha {N^{- 1}}{L_0}),}$$

where L0 is the coordinate length of the K-direction if it is compact. For any finite N the expression is well-defined and can directly be quantized, and the limit can be performed in a well-defined manner at the quantum level of the full theory.

Most importantly, the resulting operator preserves the form of symmetric states for the diagonal homogeneous model in its dual action, corresponding to the flux operator of the reduced model as used before. In averaging the full operator the partial background provided by the group action has been used, which is responsible for the degeneracy between edge length and spin in one reduced flux label. Similarly, one can obtain holonomy operators along the I-direction that preserve the form of symmetric states after averaging them along the J and K directions (in such a way that the edge length is variable in the averaging limit). Thus, the dual action of full operators is sufficient to derive all basic operators of the model from the full theory. The representation of states and basic operators, which was seen to be responsible for most effects in loop quantum cosmology, is thus directly linked to the full theory. This, then, defines the cosmological sector of loop quantum gravity.

6.3 Quantization before reduction

When quantizing a model after a classical reduction, there is much freedom even in choosing the basic representation. For instance, in homogeneous models one can use the Wheeler-DeWitt formulation based on the Schrödinger representation of quantum mechanics. In other models one could choose different smearings, e.g., treating triad components by holonomies and connection components by fluxes, since transformation properties can change from the reduced point of view (see, e.g., [52]). There is thus no analog of the uniqueness theorem of the full theory, and models constructed in this manner would have much inherent freedom even at a basic level. With the link to the full theory, however, properties of the unique representation there are transferred directly to models, resulting in analogous properties such as discrete fluxes and an action only of exponentiated connection components. This is sufficient for a construction by analogy of composite operators, such as the Hamiltonian constraint according to the general scheme.

If the basic representation is taken from the full quantization, one makes sure that many consistency conditions of quantum gravity are already observed. This can never be guaranteed when classically reduced models are quantized since then many consistency conditions trivialize as a consequence of simplifications in the model. In particular, background independence requires special properties, as emphasized before. A symmetric model, however, always incorporates a partial background and within a model alone one cannot determine which structures are required for background independence. In loop quantum cosmology, on the other hand, this is realized thanks to the link to the full theory. Even though a model in loop quantum cosmology can also be seen as obtained by a particular minisuperspace quantization, it is distinguished by the fact that its representation is derived by quantizing before performing the reduction.

In general, symmetry conditions take the form of second class constraints since they are imposed for both connections and triads. It is often said that second class constraints always have to be solved classically before the quantization because of quantum uncertainty relations. This seems to make impossible the above statement that symmetry conditions can be imposed after quantizing. It is certainly true that there is no state in a quantum system satisfying all second class constraints of a given reduction. Also using distributional states, as required for first class constraints with zero in the continuous spectrum, does not help. The reduction described above thus does not simply proceed in this way by finding states, normalizable or distributional, in the full quantization. Instead, the reduction is done at the operator algebra level, or alternatively the selection of symmetric states is accompanied by a reduction of operators which, at least for basic ones, can be performed explicitly. In general terms, one does not look for a sub-representation of the full quantum representation, but for a representation of a suitable subalgebra of operators related to the symmetry. This gives a well-defined map from the full basic representation to a new basic representation for the model. In this map, non-symmetric degrees of freedom are removed irrespective of the uncertainty relations from the full point of view.

Since the basic representations of the full theory and the model are related, it is clear that similar ambiguities arise in the construction of composite operators. Some of them are inherited directly, such as the representation label j one can choose when connection components are represented through holonomies [116]. Other ambiguities are reduced in models since many choices can result in the same form or are restricted by adaptations to the symmetry. This is for instance the case for positions of new vertices created by the Hamiltonian constraint. However, also new ambiguities can arise from degeneracies such as that between spin labels and edge lengths resulting in the parameter δ in Section 5.4. Also factor ordering can appear more ambiguously in a model and lead to less unique operators than in the full theory. As a simple example we can consider a system with two degrees of freedom (q1,p1; q2,p2) constrained to be equal to each other: C1 = q1q2, C2 = p1p2. In the unconstrained plane (q1, q2), angular momentum is given by J = q1p2q2p1 with an unambiguous quantization. Classically, J vanishes on the constraint surface C1 = 0 = C2, but in the quantum system ambiguities arise: q1 and p2 commute before but not after reduction. There is thus a factor ordering ambiguity in the reduction which is absent in the unconstrained system. Since angular momentum operators formally appear in the volume operator of loop quantum gravity, it is not surprising that models have additional factor ordering ambiguities in their volume operators. Fortunately, they are harmless and result, e.g., in differences as an isotropic volume spectrum |μ,|3/2 compared to \(\sqrt {(\vert \mu \vert - 1)\vert \mu \vert (\vert \mu \vert + 1)}\), where the second form [37] is closer to SU(2) as compared to U(1) expressions.

6.4 Minisuperspace approximation

Most physical applications in quantum gravity are obtained in mini- or midisuperspace approximations by focusing only on degrees of freedom relevant for a given situation of interest. Other degrees of freedom and their interactions with the remaining ones are ignored so as to simplify the complicated full dynamics. Their role in particular for the evolution, however, is not always clear, and so one should check what happens if they are gradually tuned in.

There are examples, in the spirit of [144], where minisuperspace results are markedly different from less symmetric ones. In those analyses, however, already the classical reduction is unstable, or back reaction is important, and thus solutions that start almost symmetric move away rapidly from the symmetric submanifold of the full phase space. The failure of a minisuperspace quantization in those cases can thus already be decided classically and is not a quantum gravity issue. Even a violation of uncertainty relations, which occurs in any reduction at the quantum level, is not automatically dangerous but only if corresponding classical models are unstable.

As for the general approach to a classical singularity, the anisotropic behavior and not so much inhomogeneities is considered to be essential. Isotropy can indeed be misleading, but the anisotropic behavior is more characteristic. In fact, relevant features of full calculations on a single vertex [85] agree with the anisotropic [48, 62], but not the isotropic behavior [46]. Also patching of homogeneous models to form an inhomogeneous space reproduces some full results even at a quantitative level [70]. The main differences and simplifications of models can be traced back to an effective Abelianization of the full SU(2)-gauge transformations, which is not introduced by hand in this case but a consequence of symmetries. It is also one of the reasons why geometrical configurations in models are usually easier to interpret than in the full theory. Most importantly, it implies strong conceptual simplifications since it allows a triad representation in which the dynamics can be understood more intuitively than in a connection representation. Explicit results in models have thus been facilitated by this property of basic variables, and therefore a comparison with analogous situations in the full theory is most interesting in this context, and most important as a test of models.

If one is using a quantization of a classically reduced system, it can only be considered a model for full quantum gravity. Relations between different models and the full theory are important in order to specify to what degree such models approximate the full situation, and where additional correction terms by the ignored degrees of freedom have to be taken into account. This is under systematic investigation in loop quantum cosmology.

6.5 Quantum geometry from models to the full theory

By now, many models are available explicitly and can be compared with each other and the full theory. Original investigations were done in isotropic models, which in many respects are special, but important aspects of the loop quantization are now known to be realized in all models and sometimes the full theory without contradictions so far. There is thus a consistent picture of singularity-free dynamical behavior together with candidates for characteristic phenomenology.

There are certainly differences between models, which can be observed already for geometrical spectra such as area or volume. Akin to level splitting in atoms or molecules, spectra become more complicated when symmetry is reduced [66, 74, 65]. Also the behavior of densities or curvatures on arbitrary geometrical configurations can be different in different models. In isotropic models, densities are bounded, which is a kinematical statement but in this case important for a singularity free evolution. It is important here since minisuperspace is just one-dimensional and so dynamical trajectories could not pass regions of unbounded curvature should they exist. Anisotropic models are more characteristic for the approach to classical singularities, and here curvature expressions in general remain unbounded if all of minisuperspace is considered. Again, this is only kinematical, and here the dynamics tells us that evolution does not proceed along directions of unbounded curvature. This is similar in inhomogeneous models studied so far.

In the full theory the situation becomes again more complicated since here densities can be unbounded even on degenerate configurations of vanishing volume eigenvalue [85]. In this case, however, it is not known what the significance for evolution is, or even the geometrical meaning of the degenerate configurations.

As an analogy one can, as before, take spectroscopy of atoms and level splitting. Essential properties, such as the stability of the hydrogen atom in quantum mechanics as opposed to the classical theory, are unchanged if complicated interactions are taken into account. In this context, it is important to take into account that stability can and does change if arbitrary interactions would be considered, rather than realistic ones, which one already fixed from other observations. Hydrogen then remains stable under those realistic interactions, but its properties would change drastically if any possible interaction term would be considered. Similarly, it is not helpful to consider the behavior of densities on arbitrary geometries unless it is known which configurations are important for dynamics or at least their geometrical role is clear. Dynamics in the canonical picture is encoded in the Hamiltonian constraint, and including it (or suitable observables) in the analysis is analogous, in the picture of atomic spectra, to making use of realistic gravitational interaction terms. In the full theory, such an analysis is currently beyond reach, but it has been extensively studied in loop quantum cosmology. Since the non-singular behavior of models, whether or not curvature is bounded, is a consequence of basic effects and the representation derived from the full theory, it can be taken as reliable information on the behavior in quantum geometry.

7 Philosophical Ramifications

In the context of loop quantum cosmology or loop quantum gravity in general, some wider issues arise that have already been touched briefly. This has to be seen in the general context of what one should expect from quantum theories of gravity for which there are several quite different approaches. These issues deal with questions about the uniqueness of theories or solutions and what information is accessible in one universe. Also the role of time plays a more general role, and the related question of unitarity or determinism.

7.1 Unique theories, unique solutions

It is often the case that, before quantitative concepts can be introduced into a field of science, they are preceded by comparative concepts that are much more effective tools for describing, predicting, and explaining than the cruder classificatory concepts.

Rudolf Carnap

An Introduction to the Philosophy of Science

The rise of loop quantum gravity presents an unprecedented situation in physics where full gravity is tackled in a background independent and non-perturbative manner. Not surprisingly, the result is often viewed skeptically since it is very different from other well-studied quantum field theories. Usually, intuition in quantum field theory comes either from models which are so special that they are completely integrable, or from perturbative expansions around free field theories. Since no relevant ambiguities arise in this context, ambiguities in other frameworks are usually viewed with suspicion. A similar treatment is not possible for gravity because a complete formulation as a perturbation series around a free theory is unavailable and would anyway not be suitable in important situations of high curvature. In fact, reformulations as free theories exist only in special, non-dynamical backgrounds such as Minkowski space or planar waves which, if used, immediately introduce a background.

If this is to be avoided in a background independent formulation, it is necessary to deal with the full non-linear theory. This leads to complicated expressions with factor ordering and other ambiguities, which are usually avoided in quantum field theory but not unfamiliar from quantum theory in general. Sometimes it is said that such a theory looses its predictive power or even suggested to stop working on applications of the theory until all ambiguities are eliminated. This view, of course, demonstrates a misunderstanding of the scientific process where general effects play important roles even if they can be quantified only at later stages. What is important is to show that qualitative effects are robust enough such that their implications do not crucially depend on one choice among many.

So far, applications of loop quantum gravity and cosmology are in comparative stages where reliable effects can be derived from basic properties and remaining ambiguities preclude sharp quantitative predictions in general (notable exceptions are fundamental properties, such as the computation of γ through black hole entropy [11, 12, 108, 155]). These ambiguities have to be constrained by further theoretical investigations of the overall consistency, or by possible observations.

Ambiguities certainly mean that a theory cannot be formulated uniquely, and uniqueness often plays a role in discussions of quantum gravity. In the many approaches different kinds of uniqueness have been advertised, most importantly the uniqueness of the whole theory, or the uniqueness of a solution appropriate for the one universe we can observe. Both expectations seem reasonable, though immodest. But they are conceptually very different and even, maybe surprisingly, inconsistent with each other as physical properties: For let us assume that we have a theory from which we know that it has one and only one solution. Provided that there is sufficient computational access to that theory, it is falsifiable by comparing properties of the solution with observations in the universe. Now, our observational access to the universe will always be limited and so, even if the one solution of our theory does agree with observations, we can always find ways to change the theory without being in observational conflict. The theory thus cannot be unique. Changing it in the described situation may only violate other, external conditions that are not observable.

The converse, that a unique theory cannot have a unique solution, follows by logically reversing the above argument. However, one has to be careful about different notions of uniqueness of a theory. It is clear from the above argument that uniqueness of a theory can be realized only under external, such as mathematical, conditions, which always are a matter of taste and depend on existing knowledge. Nevertheless, the statement seems to be supported by current realizations of quantum gravity. String theory is one example where the supposed uniqueness of the theory is far outweighed by the non-uniqueness of its solutions. It should also be noted that the uniqueness of a theory is not falsifiable, and therefore not a scientific claim, unless its solutions are sufficiently restricted within the theory. Otherwise, one can always find new solutions if one comes in conflict with observations. A theory itself, however, is falsifiable if it implies characteristic effects for its solutions even though it may otherwise be ambiguous.

7.2 The role of time

Dies alles dauerte eine lange Zeit, oder eine kurze Zeit: denn, recht gesprochen, gibt es für dergleichen Dinge auf Erden keine Zeit.

(All this took a long time, or a short time: for, strictly speaking, for such things no time on earth exists.)

Friedrich Nietzsche

Thus Spoke Zarathustra

Often, time is intuitively viewed as coordinate time, i.e., one direction of space-time. However, this does not have invariant physical meaning in general relativity, and conceptually an internal time is more appropriate. Evolution is then measured in a relational manner of some degrees of freedom with respect to others [33, 170, 106]. In quantum cosmology, as we have seen, this concept is even more general since internal time keeps making sense at the quantum level also around singularities where the classical space-time dissolves.

The wave function thus extends to a new branch beyond the classical singularity, i.e., to a classically disconnected region. Intuitively this leads to a picture of a collapsing universe preceding the Big Bang, but one has to keep in mind that this is the picture obtained from internal time where other time concepts are not available. In such a situation it is not clear, intuitive pictures notwithstanding, how this transition would be perceived by observers were they able to withstand the extreme conditions. It can be said reliably that the wave function is defined at both sides, “before” and “after”, and every computation of physical predictions, e.g., using observables, we can do at “our” side can also be done at the other side. In this sense, quantum gravity is free of singularities and provides a transition between the two branches. The more complicated question is what this means for evolution in a literal sense of our usual concept of time (see also [200]).

Effective equations displaying bounces in coordinate time evolution indicate that indeed classical singularities are replaced by a bouncing behavior. However, this does not occur completely generally and does not say anything about the orientation reversal which is characteristic for the quantum transition. In fact, effective equations describe the motion of semiclassical wave packets, which becomes less reliable at very small volume. And even if the effective bounce happens far away from the classical singularity will there in general be a part of the wave function splitting off and traversing to the other orientation as can be seen in the example of Figure 10.

Figure 10
figure 10

mpg-Movie (866 KB) Still from a Movie showing the coordinate time evolution [72] of a wave packet starting at the bottom and moving toward the classical singularity (vertical dotted line) for different values of an ambiguity parameter. Some part of the wave packet bounces back (and deforms) according to the effective classical solution (dashed), but other parts penetrate to negative μ. The farther away from a = 0 the effective bounce happens, depending on the ambiguity parameter, the smaller the part penetrating to negative μ is. The coordinate time evolution represents a physical state obtained after integrating over t [72]. (For video see appendix)

It is not clear in general that a wave function penetrating a classical singularity enters a new classical regime even if the volume becomes large again. For instance, there can be oscillations on small scales, i.e., violations of pre-classicality, picked up by the wave function when it travels through the classical singularity. As discussed in Section 5.17, the question of what conditions on a wave function to require for a classical regime is still open, but even if one can confidently say that there is such a new classical region does the question arise if time continues during the transition through the pure quantum regime. At least in the special model of a free massless scalar in isotropic cosmology the answer to both questions is affirmative, based on the availability of a physical inner product and quantum observables in this model [24].

Also related to this context is the question of unitary evolution. Even if one uses a selfadjoint constraint operator, unitary evolution is not guaranteed. First, the constraint splits into a time generator part containing derivatives or difference operators with respect to internal time and a source part containing, for instance, the matter Hamiltonian. It is then not guaranteed that the time generator will lead to unitary evolution. Secondly, it is not obvious in what inner product to measure unitarity since the constraint is formulated in the kinematical Hilbert space but the physical inner product is relevant for its solutions. This shows that the usual expectation of unitary evolution, commonly motivated by preservation of probability or normalization of a wave function in an absolute time parameter, is not reliable in quantum cosmology. It must be replaced by suitable conditions on relational probabilities computed from physical wave functions.

7.3 Determinism

Hat die Zeit nicht Zeit? (Does time not have time?)

Friedrich Nietzsche

Beyond Good and Evil

Loosely related to unitarity, but more general, is the concept of determinism. This is usually weakened in quantum mechanics anyway since in general one makes only probabilistic statements. Nevertheless, the wave function is determined at all times by its initial values, which is sometimes seen as the appropriate substitute for deterministic behavior. In loop quantum cosmology the situation again changes slightly since, as discussed in Section 5.18, the wave function may not be determined by the evolution equation everywhere, i.e., not at points of classical singularities, and instead acquire new conditions on its initial values. This could be seen as a form of indeterministic behavior, even though the values of a wave function at classical singularities would not have any effect on the behavior for non-degenerate configurations.Footnote 2 (If they had such an effect, the evolution would be singular.) In this situation one deals with determinism in a background independent context, which requires a new view.

In fact, rather than interpreting the freedom of choosing values at classical singularities as indeterministic behavior, it seems more appropriate to see this as an example for deterministic behavior in a background independent theory. The internal time label μ first appears as a kinematical object through the eigenvalues of the triad operator (46). It then plays a role in the constraint equation (49) when formulated in the triad representation. Choosing internal time is just made for convenience, and it is the constraint equation that must be used to see if this choice makes sense in order to formulate evolution. This is indeed the case at non-zero μ where we obtain a difference operator in the evolution parameter. At zero μ, however, the operator changes and does not allow us to determine the wave function there from previous values. Now, we can interpret this simply as a consequence of the constraint equation rejecting the internal time value μ = 0. The background independent evolution selects the values of internal time it needs to propagate a wave function uniquely. As it turns out, μ = 0 is not always necessary for this and thus simply decouples. In hindsight, one could already have split off |0⟩ from the kinematical Hilbert space, thereby removing the classical singularity by hand. Since we did not do this, it is the evolution equation that tells us that this is happening anyway. Recall, however, that this is only one possible scenario obtained from a non-symmetric constraint. For the evolution (50) following from the symmetric constraint, no decoupling happens and μ = 0 is just like any other internal time value.

8 Research Lines

Currently, the development of loop quantum cosmology proceeds along different lines, at all levels discussed before. We present here a list of the main ones, ordered by topics rather than importance or difficulty.

8.1 Conceptual issues

The list of conceptual issues is not much different from but equally pressing as in quantum gravity in general. Here, mainly the issue of time (its interpretation, different roles and explicit implementation into physics), the interpretation of the wave function in quantum theory, and technical as well as conceptual questions related to the physical inner product need to be addressed.

8.2 Mathematical development of models

The main open issue, requiring new insights at all levels, is that of inhomogeneities. While inho-mogeneous models have been formulated and partly analyzed, the following tasks are still to be completed:

  • Exact models: In particular the dynamics of inhomogeneous models is much more complicated to analyze than in homogeneous ones. Understanding may be improved by an interesting cross-relation with black holes. This allows one to see if the different ingredients and effects of a loop quantization fit together in a complete picture, which so far seems to be the case [14, 159, 13, 76, 57]. Moreover, the dynamics can possibly be simplified and understood better around slowly evolving horizons [79, 76]. Other horizon conditions are also being studied in related approaches [132, 107].

  • Consistency: Not directly related to physical applications but equally important is the issue of consistency of the constraints. The constraint algebra trivializes in homogeneous models, but is much more restrictive with inhomogeneities. Here, the feasibility of formulating a consistent theory of quantum gravity can be tested in a treatable situation. Related to consistency of the algebra, at least at a technical level, is the question of whether or not quantum gravity can predict initial conditions for a universe, or at least restrict its set of solutions.

  • Relation between models and the full theory: By strengthening the relation between models and the full theory, ideally providing a complete derivation of models, physical applications will be put on a much firmer footing. This is also necessary to understand better effects of reductions such as degeneracies between different concepts or partial backgrounds. One aspect not realized in models so far is the large amount of non-Abelian effects in the full theory, which can be significant also in models [54].

  • Numerical quantum gravity: Most systems of difference equations arising in loop quantum gravity are too complicated to solve exactly or even to analyze. Special techniques, such as those in [72, 87, 71, 24, 88] have to be developed so as to apply to more general systems. In particular for including inhomogeneities, both for solving equations and interpreting solutions, a new area of numerical quantum gravity has to be developed.

  • Perturbations: If the relation between different models is known, as presently realized for isotropic within homogeneous models [63], one can formulate the less symmetric model perturbatively around the more symmetric one. This then provides a simpler formulation of the more complicated system, easing the analysis and uncovering new effects. In this context, also alternative methods to introduce approximate symmetries, based on coherent states as e.g., advocated in [84], exist.

  • Effective equations: Finding effective equations that capture the quantum behavior of basic difference equations, at least in some regimes, will be most helpful for a general analysis. However, their derivation is much more complicated for inhomogeneous systems owing to the consistency issue. On the other hand, trying to derive them will provide important tests for the framework, in addition to giving rise to new applications.

8.3 Applications

Once available, equations for inhomogeneous systems have the prospect of applications such as the following.

  • Structure formation: There are diverse scenarios for the early universe with a potential for viable structure formation, which can only be checked with a reliable handle on inhomo-geneities. This applies to inflaton models with loop effects, inflation models without inflaton, and the generation of structure before and subsequent propagation through a bounce.

  • Robustness: All results obtained so far have to be regarded as preliminary and their validity in the presence of perturbative inhomogeneities has to be established. A detailed analysis of their robustness to quantization freedom such as ambiguities or choosing matter fields is still to be undertaken.

  • New effects: Some cosmological issues that have not been addressed so far from loop quantum gravity, and which most likely require inhomogeneities, are: the initial state of the inflaton (Gaussianity) or the present acceleration of cosmic expansion. The latter could be a result of small, local quantum corrections adding up to a sizeable effect on the whole universe. From a technical point of view, contact to quantum gravity phenomenology in a particle physics context can be made (as initially in [70]).

  • Ansätze: For the time being, those questions can be addressed preliminaryly by choosing suitable forms of inhomogeneous equations motivated by operators in full loop quantum gravity.

8.4 Homogeneous models

There are still several open areas in homogeneous models, which later can be extended to inhomogeneous ones:

  • Conceptual issues: This has already been mentioned above. Isotropic models provide simpler settings to analyze, e.g., the physical inner product [127, 162, 24], observables, different interpretations of quantum aspects or the emergence of a classical world.

  • Effective equations: Even in isotropic models effective equations have not yet been derived completely. A general scheme exists, shown to be analogous to standard effective action techniques [73], but it remains to be applied completely to quantum cosmology. This will then lead to a complete set of correction terms and their ranges of validity and importance. Also the question of whether an effective action for quantum cosmology exists and what its form is can be addressed.

  • Matter systems: Matter systems provide a rich source of diverse scenarios, but a full analysis is yet to be done. This includes adding different kinds of fluids [163], fermions or anisotropy parameters (shear term).

8.5 Outlook

All these developments will certainly also aid and suggest developments in the full theory, and reciprocally be assisted by new ideas realized there. At the other side, guidance as well as means for testing can be expected from future observations.