# Astrophysical and cosmological signatures of Loop Quantum Gravity

Aurélien Barrau (2017), Scholarpedia, 12(10):33321. | doi:10.4249/scholarpedia.33321 | revision #183481 [link to/cite this article] |

The **loop quantum gravity** theory is a non-perturbative and background-independent quantization of general relativity. It has been developed both in the canonical and covariant formalisms. As for all tentative theories of quantum gravity, one of the key question is to produce clear experimental or observational predictions. As the discrete nature of spacetime at the Planck scale is extremely difficult to probe directly, phenomenological attempts do focus mostly on the cosmological sector, on black holes and on possible high energy astrophysics effects.

The main prediction of loop quantum gravity in the cosmological sector is that the Big Bang is replaced by a Big Bounce: a contracting phase should have taken place before the expanding phase we are now living in. In the black hole sector, the Hawking evaporation spectrum should receive substantial corrections and the very existence of an event horizon can be questioned, potentially leading to black holes bouncing into white holes. Finally, on the astroparticle physics side, the quantum discreteness of spacetime could lead to interesting effects for high energy cosmic rays.

## Contents |

## Brief introduction to loop quantum gravity

This article uses Planck units except in some formulae where the constants where reinserted to help the understanding.

### Canonical approach

In the canonical formulation, classical general relativity can be expressed in terms of a spatial metric \(q_{ab}\), describing a slice \(\Sigma\), and the extrinsic curvature \(K_{ab}\) which is the derivative of this metric with respect to a change of spatial slice. The Hamiltonian and diffeomorphism constraints ensure that those fields are indeed solutions to the Einstein’s equations, by generating deformations of the spatial slice and ensuring general covariance. In this framework, it is natural to use instead of the spatial metric the triad field \(e^a_i\) such that \(q^{ab}=\delta^{ij}e^a_ie^b_i\). This brings an additional SO(3) gauge freedom which is accounted for by the Gauss constraint.

To quantize such a constrained system, one can follow the usual Dirac prescription. To perform this task, a Hilbert space of functionals \(\psi[q]\) of spatial metrics is needed. However, the space of metrics - or extrinsic curvature tensors - is mathematically poorly understood and defining a suitable inner products (the key ingredient of the Hilbert space) is difficult. This is why new variables, introduced by Ashtekar, play an important role (an introduction can be found in Bojowald 2006). Their definition requires to use the densitized triad \(E^a_i\) which is related to the triad by \(E^a_i=\left|\det e^b_j\right|^{-1} e^a_i\) and its conjugate, the extrinsic curvature coefficients \(K_a^i:=K_{ab}e^b_i\) with \( \{K_a^i(x),E^b_j(y)\}= 8\pi G\delta^b_a\delta^i_j\delta(x-y) \). The Ashtekar connection is then introduced through \(A_a^i=\Gamma_a^i+\gamma K_a^i\), where \(\gamma\) is the Barbero-Immirzi parameter which has no classical relevance but is important for quantization, and \[ \Gamma_a^i = -\epsilon^{ijk}e^b_j (\partial_{[a}e_{b]}^k+ {\textstyle\frac{1}{2}} e_k^ce_a^l\partial_{[c}e_{b]}^l), \tag{1}\] is the spin connection of the triad. The Ashtekar connection is conjugate to the densitized triad: \[ \{A_a^i(x),E^b_j(y)\}=8\pi\gamma G\delta_a^b\delta_j^i\delta(x,y). \tag{2}\] The dynamics is implemented though the Hamiltonian constraint (\(N\) being the laps): \[ H[N]=\frac{1}{16\pi\gamma G} \int_{\Sigma} \mathrm{d}^3x N\left|\det E\right|^{-1/2} \left(\epsilon_{ijk}F_{ab}^iE^a_jE^b_k -2(1+\gamma^2) K^i_{[a}K^j_{b]}E^a_iE^b_j\right)\approx 0, \tag{3} \] where \(F_{ab}^i\) is the curvature of the Ashtekar connection.

The Gauss constraint reads \[ G[\Lambda]=\frac{1}{8\pi\gamma G}\int_{\Sigma}\mathrm{d}^3x \Lambda^iD_aE^a_i= \frac{1}{8\pi\gamma G} \int_{\Sigma}\mathrm{d}^3x\Lambda^i(\partial_aE^a_i+ \epsilon_{ijk}A^j_aE^a_k)\approx 0, \tag{4} \] and the diffeomorphism constraint is given by \[ D[N^a]= \frac{1}{8\pi\gamma G}\int_{\Sigma}\mathrm{d}^3x N^aF_{ab}^iE_i^b \approx 0. \tag{5} \]

The Algebra of constraints if closed and constitues a first class Dirac Algebra.

When quantizing a field theory, one needs to smear fields, that is to integrate them in order to get rid of the \(\delta\)-functions. In this framework, it should be done using Ashtekar variables but without spoiling the background invariance, i.e. without relying on a given a priori metric. This can be achieved by integrating the connection along a one-dimensional curve $e$ and then exponentiate in a path-ordered manner (in order to obtain a covariant object from the non-Abelian connection). This defines the holonomy: \[ h_e(A)={\cal P}\exp\int_e\tau_i A_a^i\dot{e}^a\mathrm{d}t\tag{6}, \] with tangent vector \(\dot{e}^a\) to the curve \(e\) and \(\tau_j=-\frac{1}{2}i\sigma_j\) where \(\sigma_j\)’s are Pauli matrices. Following the same logic, the densitized triads can be integrated out over a 2-surface \(S\), resulting in the flux \[ F_S(E)=\int_S \tau^i E^a_in_a\mathrm{d}^2y\tag{7}, \] where $n_a$ is the co-normal to the considered surface.

The Poisson algebra of holonomies and fluxes is well-defined and can be represented on an appropriate Hilbert space. When diffeomorphism invariance is imposed there is a unique representation which defines the kinematical Hilbert space. The Hilbert space is basically the space of an \(SU(2)\) lattice Yang-Mills theory, which can be interpreted as a description of quantized geometries. The dynamics is given by the quantum Hamiltonian constraint (see for example Ashtekar 2012). These are the key ingredient of loop quantum gravity (LQG).

### Covariant approach

There is also path-integral formulation of LQG, explicitly covariant, relying on spin foams (see Rovelli 2011). This approach builds on the formal "sum over 4-geometries" \[ Z \sim \int Dg \ \ e^{\frac{i}{\hbar}\! \int \! \!R\sqrt{g}\,d^4\!x}. \tag{8} \] This expression is however hard to use from the practical viewpoint.

In four-dimensional Lorentzian LQG, the partition function is defined by \[ Z_{\cal C}=\sum_{j_f,\mathtt{v}_e} \ \prod_f (2j_f+1)\ \prod_v A_v(j_f,\mathtt{v}_e), \tag{9} \] where \(\cal C\) is a 2-complex with faces \(f\), edges \(e\) and vertices \(v\), whereas \(\mathtt{v}_e\) are the intertwines, and \[ A_v(j_f,\mathtt{v}_e)=Tr\left[\otimes_{e\in v} (f_\gamma \mathtt{v}_{e})\right], \tag{10} \] where \(f_\gamma\) sends \(SU(2)\) spin networks into \(SL(2,\C)\) spin networks. This basically determines the covariant dynamics of LQG. It is worth noticing that the theory is determined by the imbedding of \(SU(2)\) functions into \(SL(2,\C)\) functions.

A key LQG result for phenomenology - that will be used in the following - is that areas are now quantized with a spectrum \[ A_j=8\pi \gamma l_{Pl}^2\sqrt{j(j+1)}, \tag{11} \] where \(j\) is a half integer.

## Loop quantum cosmology

### Basic ideas

Loop quantum cosmology (LQC) is the application of loop quantum gravity ideas to the Universe itself (see Ashtekar & Singh 2004). Let us begin by considering isotropic and homogeneous loop quantum cosmology. The classical FLRW (Friedmann-Lemaître-Robertson-Walker) universe is described by the scale factor \(a(t)\), which is a solution to the Einstein's equations. In Ashtekar variables, the scale factor and its time derivative are expressed by a densitized triad component \(\tilde{p}\) with \(|\tilde{p}|=a^2\), determined by the isotropic densitized triad \(E^a_i=\tilde{p}\delta^a_i\), and a connection variable \(\tilde{c}=\gamma\dot{a}\), determined by the isotropic connection \(A_a^i=\tilde{c}\delta_a^i\). The spatial volume of a given region is \(V=\int_Rd^3x\sqrt{|\det E|}\) which, for an isotropic triad, reduces to \(V=V_0|\tilde{p}|^{3/2}=|p|^{3/2}\) where the coordinate volume \(V_0\) is absorbed in \(p\). The re-scaled variable \(p\), and similarly \(c=V_0^{1/3}\tilde{c}\), is then independent of the fiducial cell. The Poisson bracket is given by: \[ \{c,p\}=\frac{8\pi\gamma G}{3}\,. \tag{12} \]

Following LQG ideas, these fundamental objects are represented on a Hilbert space with an orthonormal basis \(\{|\mu\rangle\}_{\mu\in\R}\) of states given as functions of the connection component by \(\langle c|\mu\rangle = e^{i\mu c/2}\). Terms like \(e^{i\mu'c/2}\) are analogous to holonomies of the full theory and define a shift operator: \[ \widehat{e^{i\mu'c/2}}|\mu\rangle = |\mu+\mu’\rangle, \tag{13} \] while the triad component $p$ acts by multiplication: \[ \hat{p}|\mu\rangle = \frac{1}{6}\gamma l_{Pl}^2\mu|\mu\rangle\,. \tag{14} \] An important property of the mother theory is preserved: only exponentials of \(c\) are represented as it is not possible to obtain an operator for the component \(c\) directly, and the triad operator \(\hat{p}\) has a discrete spectrum since its eigenstates are normalizable. From \(\hat{p}\) one directly obtains the volume operator \(\hat{V}=|\hat{p}|^{3/2}\).

In LQC, the basic equation is a difference equation and not a differential equation as in the Wheeler-deWitt theory (see Bojowald 2012). Written in terms of the wave function of the Universe (filled with a scalar field \(\phi\)) it can be shown to be (Bojowald 2003,Bojowald 2004) in the isotropic and flat case:

\[
(V_{\mu+5}-V_{\mu+3})e^{ik}\psi_{\mu+4}(\phi)- 2(V_{\mu+1}-V_{\mu-1})\psi_{\mu}(\phi)
+ (V_{\mu-3}-V_{\mu-5})e^{-ik}\psi_{\mu-4}(\phi)
= -\frac{4}{3}\pi \gamma^3G\ell_{\rm P}^2\hat{H}_{\rm matter}(\mu)\psi_{\mu}(\phi),
\tag{15}
\]
where \(V_{\mu}\) are the volume eigenvalues and \(\hat{H}_{\rm matter}(\mu)\) is the matter Hamiltonian.

This equation is non-singular (Bojowald, 2001). Starting from a large volume, it is possible to evolve the wave function of the Universe backward through the classical singularity. The evolution does not stop, and one obtains a collapsing branch preceding the classical singularity (see details in Bojowald 2006). This is the Big Bounce scenario. The existence of this bounce, replacing the Big Bang, has been shown in LQC to be true for exactly solvable models (a flat, isotropic universe filled with a massless scalar field), in the presence of spatial curvature, for Bianchi models, with both a positive and a negative cosmological constant, for \(\phi^2\) inflationary potentials, for barotropic matter with different equations of state, and for wide states corresponding to highly quantum universes. In addition, numerical simulations (see Singh 2012, Diener et al. 2014) have shown that the quantum evolution is quite generically well described by effective semi-classical equations.

The main basic LQC result is that the effective Friedmann equation now reads: \[ H^2=\frac{8\pi G\rho}{3}\left(1-\frac{\rho}{\rho_c}\right), \tag{16} \] where \(\rho\) is the density and \(\rho_c\) is the critical density (expected to be of the order of the Planck density and more precisely equal to \(0.4\rho_{Pl}\) under quite natural hypotheses). The usual Friedmann equation receives a quadratic correction with the correct sign so that the Universe now bounces.

At this stage, loop quantum cosmology is a framework that import techniques and ideas from LQG, but it is not the full quantum theory. There could be higher order effects missing in the effective cosmological equations. However, the effective equations seem to capture the most relevant quantum effects, and the conclusions are supported by results obtained using techniques in the full theory. This is the case of the singularity resolution in full covariant LQG (see Bianchi et al. 2010), or the techniques of gauge fixing (referred as "reduced LQG", see Alesci et al. 2014) that indeed give an Hamiltonian yielding the same modification of the Friedman equations as in standard LQC.

### Predictions for the background

The clear background dynamics description is one of the key successes of LQC (see Ashtekar et al. 2006). Interestingly, this result has also been recovered in the framework of group field theory (see Oriti et al. 2016). The large scale cosmological dynamics is then described by the hydrodynamics of condensate states in the Gross-Pitaevskii approximation and the same modified equations is derived (under some specific assumptions).

However the model becomes also interesting because it does more than predicting a bounce (see Fig. 1 for the dynamics of the field). Under some hypotheses, it might also predict the duration of inflation. This is something quite specific (see Ashtekar & Sloan 2011). The fact that inflation occurs generically in LQC once the appropriate content - e.g. a massive scalar field - is assumed should not come as a surprise: inflation is an attractor and starting with an energy density close to the Planck density and the appropriate content yields nearly inevitably to inflation. The interest of LQC is here more subtle.

One can define the fractions of potential and kinetic energy, normalized to the maximum energy density, \[ x := \frac{m\phi}{\sqrt{2\rho_{\text{c}}}} \quad \text{and} \quad y :=\frac{\dot{\phi}}{\sqrt{2\rho_{\text{c}}}}, \tag{17} \] so that \[ \rho=\rho_c\left(x^2+y^2\right). \tag{18} \] Interestingly, in the contracting phase of the Universe, far away (before) from the Big Bounce, the density can be shown to be well approximated by \[ \rho=\rho_0\left(1-\frac{1}{2}\sqrt{3\kappa\rho_0}\left( t+\frac{1}{2m}\sin(2mt+2\delta)\right)\right)^{-2}.\tag{19} \] The \(\delta\) parameter is just a phase and can be assigned a flat a priori probability distribution function (PDF). In addition of being the expected distribution for any random oscillatory process of this kind, a flat PDF for \(\delta\) will be preserved over time within the prebounce oscillation phase, making it a very natural choice for initial conditions. The important point is that starting from a flat PDF for \(\delta\) and a small enough initial density \(\rho_0\), the distribution of the fraction of potential energy at the bounce \(x_B\) can be calculated and is far from being flat between 0 and 1. It is sharply peaked around \(3.5\times 10^{-6}\) (this value scales with \(m\), the mass of the scalar field, as \(m\log\left(\frac{1}{m}\right)\) and is here given for the usual value around \(10^{-6}\) in Planck units).

Remarkably, this distribution for \(x_B\) can be translated into a PDF for the number \(N\) of e-folds of inflation (defined as the logarithm of the ratio of the scale factor evaluated just after and just before inflation). The number of e-folds could a priori be any number up to a huge value bounded by the total energy available. If one assumes a flat PDF for \(x_B\), as sometimes done in cosmology, the most probable value of \(N\) is around \(10^{12}\). But if the causal evolution is taken seriously and if the PDF for \(x_B\) expected from LQC is taken into account, the prediction changes drastically and leads to a value of \(N\) peaked around 140 (see Linsefors & Barrau 2013a).

This result assumes isotropy. However, in bouncing cosmologies, the issue of anisotropies is crucial for a simple reason: the shear term scales as \(1/a^6\) where \(a\) is the mean scale factor of the Universe. Therefore, when the Universe is in its contraction phase, it is expected that the shear term eventually dominates and drives the dynamics. The LQC background dynamics has also been developed in anisotropic settings (see Ashtekar & Wilson-Ewing 2009, Linsefors & Barrau 2013b). The key result is that the bounce prediction remains correct even if the shear drives the dynamics (see Linsefors & Barrau 2014). In a Bianchi-I universe, the LQC-modified effective Friedmann equation reads as \[ H^2=\sigma_Q+\frac{\kappa}{3}\rho-\lambda^2\gamma^2\left(\frac{3}{2}\sigma_Q+\frac{\kappa}{3}\rho\right)^2, \tag{20} \] where \(\sigma_Q\) defines the "quantum shear" as: \[ \sigma_Q:=\frac{1}{3\lambda^2\gamma^2}\left(1-\frac{1}{3}\Big[\cos(\bar{\mu}_1c_1-\bar{\mu}_2c_2)+\cos(\bar{\mu}_2c_2-\bar{\mu}_3c_3)+\cos(\bar{\mu}_3c_3-\bar{\mu}_1c_1) \Big]\right), \tag{21} \] with \[ \bar{\mu}_1 = \lambda\sqrt{\frac{p_1}{p_2p_3}} \qquad \text{and cyclic expressions,} \tag{22} \] where \(\lambda\) is the square root of the minimum area eigenvalue of the LQG area operator (\(\lambda=\sqrt{\Delta}\)) and the \(p_i\)’s are conjugate to the \(c_i\)’s (both are diagonal elements of the Ashtekar variables) so that: \[ \{c_i,p_j\}=\kappa\gamma\delta_{ij}\quad ,\quad \{\phi_n,\pi_m\}=\delta_{mn}. \tag{23} \] In this framework, it is possible to evaluate the PDF for the number of e-folds of inflation. As expected it is lowered when compared to the isotropic setting. This is promising as this might make the actual number of e-folds close to the experimental lower bound which corresponds to cases with a rich phenomenology. (If the number of e-folds is much greater than 70, the interesting quantum gravity effects are beyond the horizon.)

### Predictions for the perturbations

The situation for perturbations is less clear than for the background (see Barrau et al. 2013). Different approaches to treat cosmological perturbations have been developed in LQC.

The dressed metric approach (see Agullo et al. 2012, Agullo et al. 2013) relies on a minisuperspace strategy in which the homogeneous and isotropic degrees of freedom, together with the inhomogoneous degrees of freedom are quantized. The former is obtained by the loop quantization and the latter is obtained through a Fock-like quantization on a quantum background (see Ashtekar et al. 2009). The inhomogeneous degrees of freedom considered as small perturbations are given by the Mukhanov-Sasaki variables derived from the linearized classical constraints. This field is essential in the theory of cosmological perturbations because it is gauge-invariant. The second order Hamiltonian is promoted to be an operator and the quantization is carried out relying on specific techniques suitable for the quantization of a test field evolving in a quantum background. The Hilbert space is the tensor product \(\Psi(\nu,v_\mathrm{S(T)},\varphi)=\Psi_\mathrm{FLRW}(\nu,\bar\varphi)\otimes\Psi_\mathrm{pert}(v_\mathrm{S},v_\mathrm{T},\bar\varphi)\) with \(\nu\) the homogeneous and isotropic degrees of freedom and \(v_\mathrm{S(T)}\) the degrees of freedom associated with perturbations. As long as the backreaction of the perturbations on \(\Psi_{\rm FLRW}\) is negligible, the Schrödinger equation for the perturbations was shown to be the same as the Schrödinger equation for the quantized perturbations evolving in a classical background but using a dressed metric. It reads for tensor modes: \[ i\hbar\partial_{\bar\varphi}\Psi_\mathrm{pert}=\frac{1}{2}\displaystyle\int \frac{d^3k}{(2\pi)^3}\left\{\frac{32\pi G}{\tilde{p}_{\varphi}}\left|\hat\pi_{\mathrm{T},\vec{k}}\right|^2\Psi_\mathrm{pert}+\frac{k^2}{32\pi G}\frac{\tilde{a}^4({\bar\varphi})}{\tilde{p}_{\varphi}}\left|\hat{v}_{\mathrm{T},\vec{k}}\right|^2\Psi_\mathrm{pert}\right\}, \tag{24} \] with \[ (\tilde{p}_{\varphi})^{-1}=\left<\hat{H}^{-1}_\mathrm{FLRW}\right> ~\mathrm{and}~\tilde{a}^4=\frac{\left<\hat{H}^{-1/2}_\mathrm{FLRW}\hat{a}^4({\bar\varphi})\hat{H}^{-1/2}_\mathrm{FLRW}\right>}{\left<\hat{H}^{-1}_\mathrm{FLRW}\right>}, \tag{25} \] where \((\hat{v}_{\mathrm{T},\vec{k}},\hat\pi_{\mathrm{T},\vec{k}})\) are the configuration and momentum operators of the perturbations while \(\hat{H}_\mathrm{FLRW}\) is the Hamiltonian operator for the background. A Fock quantization is finally performed and the modes are solutions of \[ Q''_k+2\left(\frac{\tilde{a}'}{\tilde{a}}\right)Q'_k+\left(k^2+\tilde{U}\right)Q_k=0, \\ h''_k+2\left(\frac{\tilde{a}'}{\tilde{a}}\right)h'_k+k^2h_k=0. \tag{26} \] The gauge-invariant variable \(Q_k\) is related with the Mukhanov-Sasaki variables for scalar modes via \(Q_k=v_{\mathrm{S},k}/a\), and \(\tilde{U}\) is a dressed potential-like term given by \[ \tilde{U}({\bar\varphi})=\frac{\left<\hat{H}^{-1/2}_\mathrm{FLRW}\hat{a}^2({\bar\varphi})\hat{U}({\bar\varphi})\hat{a}^2({\bar\varphi}) \hat{H}^{-1/2}_\mathrm{FLRW}\right>}{\left<\hat{H}^{-1/2}_\mathrm{FLRW}\hat{a}^4({\bar\varphi})\hat{H}^{-1/2}_\mathrm{FLRW}\right>}, \tag{27} \] the quantum counterpart of \[ U({\bar\varphi})=a^2\left(fV({\bar\varphi})-2\sqrt{f}\partial_{\bar\varphi} V+\partial^2_{\bar\varphi} V\right), \tag{28} \] with \(f=24\pi G (\dot{\bar\varphi}^2/\rho)\) the fraction of kinetic energy in the scalar field.

This approach leads to spectra that are compatible with the current observations (see Fig. 2). In addition, the large scale modes (appearing in the IR part of the power spectrum) can be excited and not be anymore in the Bunch-Davies vacuum at the bounce time. This opens the exciting possibility to a better fit the data in the low-l (large scale) part of the CMB spectrum. Furthermore, there are specific predictions for future missions that distinguish this LQG based mechanism from non-primordial mechanisms for power suppression.

A second approach (see Bojowald 2006 and Barrau et al. 2014), usually referred to as the "deformed algebra" one, has also been considered in LQC. It is probably "less quantum" by construction and less deeply rooted in the profound principles of loop quantum gravity. But it is more focused on the important issue of the consistency of the resulting effective theory.
When quantum corrections are taken into account at the effective level, it is hard to know whether the delicate consistency conditions summarized in the first-class nature of the classical constraint algebra remain intact.
Especially, in background-independent frameworks, it is questionnable to rely on standard covariance arguments because the very notion of space-time should emerge from solutions to the fundamental equations. The consistency of the equations must in principle be ensured before they can be solved. The preserved symmetry implies consistent equations based on first-class constraints (see Bojowald et al. 2008).
The "deformed algebra" approach puts a specific emphasis on such gauge issues. Gauge fixing before quantization was shown to be often harmless but the case of gravity is quite different from Yang-Mills theories (dynamics is part of the gauge system).
In the "deformed algebra" approach, one considers LQG-inspired corrections to the constraints. Holonomy corrections (see, e.g., Wilson-Ewing 2011) are basically implemented through the replacement
\[
\bar{k} \rightarrow \mathbb{K}[n] := \frac{\sin(n\bar{\mu} \gamma \bar{k})}{n\bar{\mu}\gamma},
\tag{29}
\]
where \(n\) is some unknown integer, \(\bar{k}\) is still the mean Ashtekar connection, and \(\bar{\mu}\) is the coordinate size of a loop.
The quantum-corrected constraints are noted \(\mathcal{C}^Q_I\). If the previous replacement is performed naively, the algebra reads
\[
\{ \mathcal{C}^Q_I, \mathcal{C}^Q_J \} = {f^K}_{IJ}(A^j_b,E^a_i) \mathcal{C}^Q_K+
\mathcal{A}_{IJ},
\tag{30}
\]
where \(\mathcal{A}_{IJ}\) stand for anomaly terms. The consistency condition (algebra closure) imposes \(\mathcal{A}_{IJ}=0\). In turn, this condition imposes restrictions on the shape of the quantum corrections.
The quantum corrected constraints are explicitly written for the perturbations. Then the Poisson brackets are calculated and anomalies are evaluated. Counterterms, which are required to vanish at the classical limit, are added to ensure the anomaly freedom. The resulting theory is not only consistent but is also quite uniquely defined. The resulting final algebra is simple and depends on a unique structure function \(\Omega=1-2\rho/\rho_c\). Finally, the unknown integers \(n_i\) entering the correction functions \(\mathbb{K}[n_i]\) can all be determined uniquely. The solution closes the algebra non-pertubatively. This strategy has been followed for vector and scalar perturbations. It has been shown (see Cailleteau et al. 2011) that a single algebra structure can be consistently written for all perturbations. It basically reads as:
\[
\left\{D[M^a],D [N^a]\right\} = D[M^b\partial_b N^a-N^b\partial_b M^a], \\
\left\{D[M^a],S^Q[N]\right\} = S^Q[M^a\partial_b N-N\partial_a M^a], \\
\left\{S^Q[M],S^Q[N]\right\} = \Omega
D\left[q^{ab}(M\partial_bN-N\partial_bM)\right], \tag{31}
\]
where \(\Omega\) is the "deformation factor" (equal to one in the classical theory with Lorentzian signature). A similar structure was found in lattice LQC (see Wilson-Ewing 2012).
This algebra leads, in particular, to an effective signature change at high density reminiscent of the Hartle-Hawking proposal. Perturbations in this framework have been analyzed according to two different philosophical positions.
In any case, the equation of motion is now more complex. It reads for tensor modes in conformal time
\[
v''_k(\eta)+\left(\Omega k^2-\frac{z_T''}{z_T}\right)v_k(\eta)=0, \tag{32}
\]
where
\(z_T\equiv \frac{a}{\sqrt{\Omega}}\). The mode functions are related to the amplitude of the tensor modes of the metric perturbation, \(h_k\), via \(v_k=z_T h_k/\sqrt{32\pi G}\). The dynamics of the modes if not anymore driven only by the hierarchy between \(k\) and \(a'’/a\) (and, hence, by the ratio between the considered length scale and the curvature radius) and quite a lot of new phenomena can happen.
The first approach considered is to use the "silent surface" (that is the surface corresponding to \(\Omega=0\) or \(\rho=\rho_c/2\)) as the natural place where to impose initial conditions. Space points are here decoupled (as anticipated by the BKL conjecture which is here recovered by the quantum corrections) and fluctuations are naturally described by a white noise spectrum. The resulting cosmological power spectrum can accommodate with the data. It is also possible to take advantage of the mathematical knowledge of the so-called Tricomi problem (a change of regime from hyperbolic to elliptic, see Bojowald & Mielczarek 2015) which leads to an interesting balance between singular Big Bang models and deterministic cyclic models.
The second approach to solving equations in this deformed algebra framework is to set initial conditions in the remote past of the contracting branch and to evolve the Fourier transform of the modes through the bounce. The Bunch Davies vacuum initial conditions can be set in the usual manner. Although this is more tricky this has also been solved for scalar modes. With \(v_S:=\sqrt{\bar{p}}\,\left(\delta\phi + \frac{{\bar{\phi}}'}{\mathscr{H}}\Phi\right)\) and \(v_S :=\sqrt{\bar{p}}\, \frac{{\bar{\phi}}'}{\mathscr{H}}\), \(\mathscr{H}\) being the conformal Hubble parameter, the Mukhanov equation of motion reads
\[
\ddot{v}_k + H \dot{v}_k + f_{k}^{\scriptscriptstyle{(v)}}(t) v_k = 0,
\tag{33}
\]
with \(z= a \frac{\dot{\bar{\phi}}}{H}\) and
\[
f_{k}^{\scriptscriptstyle{(v)}}(t) := {\bf \Omega}\frac{k^2}{a^2} - \frac{\dot{z}}{z} H - \frac{\ddot{z}}{z},
\tag{34}
\]
the effective frequency term. In terms of the scalar curvature perturbation \(\mathcal{R} := v/z\), this leads to (Schander et al. 2016)
\[
\ddot{\mathcal{R}}_k - \left( 3H + 2 m^2 \frac{\bar{\phi}}{\dot{\bar{\phi}}} + 2 \frac{\dot{H}}{H} \right) \dot{\mathcal{R}}_k + {\bf \Omega}\frac{k^2}{a^2} \mathcal{R}_k = 0.
\tag{35}
\]
The solutions to this equation are regular and were shown to be quite insensitive to the detailed choice of the initial state (e.g. instantaneous or WKB vacua). The resulting power spectrum is in disagreement with CMB data. This does not disprove LQG or LQC but underlines that a naive propagation through the euclidean phase (where there is strictly speaking no time anymore) is just incorrect. It could also be that the well known transplanckian problem (for a typical duration of inflation all the relevant scales for cosmology were smaller than the Planck length before inflation) needs to be seriously addressed in this framework, together with anisotropies. Preliminary results do suggest that the modified dispersion relation allows to recover the correct spectrum.

Important results have also been recently derived in hybrid LQC (Gomar et al. 2015, Navascues et al. 2014).

## Black holes in loop quantum gravity

In LQG a black hole is basically described as an isolated horizon (that is a horizon defined using local spacetime structures only, see Ashtekar 2004) punctured by the edges of a spin network - a graph with edges labeled with \(SU(2)\) representations and nodes characterized by intertwiners. A surface punctured by \(N\) edges has a spectrum given by \[ A_j=8 \pi \gamma l_{Pl}^2\sum_{n=1}^N\sqrt{j_n(j_n+1)}, \tag{36} \] where the sum is carried out over all intersections of the edges with the surface. Each state with spin \(j\) has a degeneracy \((2j+1)\). Many works have been devoted to derivations of the Bekenstein-Hawking entropy from the LQG formulation (see Rovelli 1996). Recent results from holographic models are very encouraging (see Ghosh et al. 2014).

It is in principle difficult to find astrophysical effects for large macroscopic black holes (although it has recently be claimed, using quite general arguments, that quantum gravity effects might be found outside the horizon, at a distance of the order of \(\frac{7}{6}R_S\), where \(R_S\) is the Schwarzschild radius, see Haggard & Rovelli 2016). Most of the phenomenology therefore focuses on small black holes.

### Evaporating black holes

One approach consists in studying how the usual Hawking evaporation spectrum will be modified by the LQG structure of the area spectrum (see Ashtekar & Bojowald 2005). Monte-Carlo simulations were performed, including the time-consuming projection constraint \[ \tag{37} \sum_{p=1}^N{m_p=0}, \] where the \(m_p\)’s carry information about the curvature of a surface whose area is given by the spin \(j\). The simulation was started at 200 \(A_{Pl}\) and the greybody factors encoding the backscattering probability within the gravitation potential were taken into account. A Kolmogorov-Smirnov test was carried out to see if it is possible to distinguish, from the measured spectrum, between the LQG predictions and the standard Hawking spectrum. The conclusion of the study is that with an experimental resolution of 5%, around \(2\times10^5\) black holes are required for a 5\(\sigma\) detection whereas \(1.2\times10^6\) are required with a 20% resolution.

A more recent study was devoted to the discrimination between different approaches to compute the BH entropy within the theory. The first class of models are those where the entropy is calculated using only quantum geometry excitations, leading to \[ \tag{38} S=\frac{\gamma_0}{\gamma} \frac{A}{4} + {o}(A), \] where \(\gamma_0\) is a numerical factor of order one depending on the details of the state counting. The quantum corrections \(o(A)\) are usually logarithmic at the leading order. The compatibility with the first law of thermodynamics imposes \(\gamma=\gamma_0\). The second class of models uses the qualitative behavior of matter degeneracy suggested by standard quantum field theory (QFT) with a cut-off at the vicinity of the horizon (i.e. an exponential growth of vacuum entanglement in terms of the BH area). The entropy then becomes \[ \tag{39} S=\frac{A}{4}+ \sqrt{\frac{\pi A}{6\gamma}} + o(\sqrt{A}). \] The Immirzi parameter \(\gamma\) affects only quantum corrections. The numerical simulations show that the spectra have two distinct parts: a continuous background which corresponds to the semi-classical stages of the evaporation and series of discrete peaks which constitute a signature of the deep quantum structure of the black hole (see Barrau et al.2015). The Immirzi parameter has an effect on both parts. The number of black holes and the instrumental resolution required to experimentally distinguish between the considered models were also studied and it was shown that the different models can be discriminated one from the other.

Finally, a new idea has also emerged, reviving the hope for possible low energy quantum gravity effects in the evaporating black hole spectrum. The density of energy levels in LQG reads as \[ \rho(M)\sim \rm{exp}(M\sqrt{4\pi G/3}), \tag{40} \] which means that the spectral lines are virtually dense in frequency. This leads to a quasi-continuum for large-mass black holes. However it might be that the evolution of black holes at each emission of a quanta is actually associated with the change of quantum state of only one puncture in the spin network (see Barrau 2016). In that case, somehow in agreement with the old Bekenstein-Mukhanov hypothesis, one should expect clear spectral lines, even for black hole masses arbitrarily far away from the Planck mass. This effect is neither washed out by the dynamics of the process, nor by the existence of a mass spectrum of primordial black holes up to a given width, nor by the secondary component induced by the decay of neutral pions emitted during the time-integrated evaporation.

### Bouncing black holes

When matter or radiation reaches the Planck density, quantum gravity might generate a sufficient pressure to counterbalance the classically attractive gravitational force. In a black hole, matter's collapse could stop before the central singularity is formed (there are quite general singularity resolution theorems in LQG, see Rovelli & Vidotto 2013). The standard event horizon of the black hole would then be replaced by an apparent horizon which is locally equivalent to an event horizon, but from which matter can eventually bounce out. In this scenario, the contracting classical black hole solution is glued to the expanding white hole solution. The important point is that such a simple minimal evolution was shown to be possible: the white hole horizon can be in the future of the black hole horizon, bounding the same external Schwarzschild region with nothing dramatic happening in the surrounding universe. This unexpected possibility is obtained by carving out the relevant solution from a double covering of the Kruskal metric (see Haggard & Rovelli 2014).

One can use coordinates \((u,v,\theta,\phi)\) with \(u\) and \(v\) null coordinates in the \(r\)-\(t\) plane and the metric is determined by two functions: \[ ds^2=-F(u,v) du dv + r^2(u,v)(d\theta^2+\sin^2\theta d\phi^2). \tag{41} \] In these Kruskal-Szekeres coordinates, the Kruskal metric can be obtained by taking \[ F(u,v)=\frac{32M^3}{r}e^{\frac{r}{2m}}, \tag{42} \] with \(r\) the function of \((u,v)\) defined by \[ \left(1-\frac{r}{2m}\right)e^{\frac{r}{2M}}=uv. \tag{43} \] The region of interest is bounded by a constant \(v={v}_o\) null line. This constant is a fundamental parameter of the metric under consideration. The key result is that, by gluing the different parts of the effective metric, and computing the minimal time for quantum gravitational effects to pile up in the region outside the horizon, one obtains an estimate for the bounce duration: \[ \tau=-8m \ln v_o> \tau_{q} = 4p\, M^2, \tag{44} \] where \(p\) was estimated to be of the order of 0.05. The bounce time is therefore proportional to \(M^2\), whereas the duration of the Hawking evaporation is proportional to \(M^3\). Let us write the actual bounce time as \(x\tau_{q}\) with \(x>1\). This leads to \[ \tau=4k M^2% \frac{M^2}{l_{Pl}}, \tag{45} \] with \(k\equiv xp>0.05\). At this stage \(k\) can be taken as a free parameter whose upper bound is such that the bounce time remains smaller than the Page time (which is smaller than the Hawking time).

Two different components are expected from bouncing black holes. The first one, referred to as the low-energy component, is just determined by dimensional analysis. The Schwarzschild radius of the black hole is the only length scale in the problem and the emitted radiation can be expected to have a mean wavelength of this order. The second one, referred to as the high-energy component, is based on the fact that the model is by construction time-symmetric: what goes out is exactly what went in (the gravitational redshift is compensated by - nearly - exactly the same amount of gravitational blueshift). The energy of the radiation emitted is therefore the same as the one of the null shell assumed to have formed the black holes in the early Universe (as only primordial black holes are phenomenologically interesting in this scenario). It has been suggested that the low energy component could explain the mysterious fast radio bursts observed in the Universe (see Barrau at al. 2014).

The detailed electromagnetic component from both those signals was studied taking into account the emission of hadrons eventually decaying into neutral pions that will subsequently produce gamma-rays (see Barrau et al. 2015). The maximum distance at which a single bouncing black hole can be detected was evaluated accounting for the absorption and the instrumental sensitivity. The result does depend on the value of \(k\). For the lowest authorized value, the maximum distance at which a bouncing black hole can be detected is larger than the Hubble radius for the low energy channel and of the order of one tenth on the Hubble scale for the high energy channel. It crosses the galactic scale for the low energy channel around \(k\sim 10^{14}\), and for the high energy channel around \(k\sim 10^{8}\). The integrated signal from a possible distribution of bouncing black holes contributing to the cold dark matter was also evaluated. The result of this study is that the measured signal should be quite close to the individual signal emitted by a single bouncing black hole, but with a slight low-energy distorsion due to redshift effects (see Barrau et al. 2016).

Finally, it has been suggested that the gamma-ray excess from the galactic center, measured by the Fermi satellite, could be accounted for by bouncing black holes. This requires a specific tuning of the free parameter of the theory and of the width of the emitted signal. This prediction is interesting because a specific redshift dependance for this signal can be estimated. The measure wavelength reads as \[ \lambda_{obs}^{BH}\sim \frac{2GM}{c^2} (1+z) \times \sqrt{\frac{H_0^{-1}}{6\,k\Omega_\Lambda^{\,1/2}}Argsinh\!\!\left[ \left(\frac{\Omega_\Lambda}{\Omega_M}\right)^{\!1/2} (z + 1)^{-3/2}\right]}, \tag{46} \] where we have reinserted the Newton constant \(G\) and the speed of light \(c\) to make things easier to read. The important consequence of this formula is that the redshift dependance of this specify model is much weaker that for other origins (either from astrophysical origin or due to decaying dark matter): the curve becomes quite flat for z>2.

Recently, more rigorous calculations for the transition amplitude \(W(m,T)\), defined such that its modulus squared determines the probability density for the bounce to happen in a given time \(T\) for a given mass \(m\), were carried out in loop quantum gravity (see Christodoulou et al. 2016). They are based on the fact that, in LQG, the amplitude associated to a two complex \(C\) is given by \[ W_{\mathcal C}(h_\ell)=N_{\mathcal C}\int_{{}_{SU(2)}} dh_{fv}\ \prod_f \delta\Big(\!\prod_{v\in f} h_{fv}\!\Big) \ \prod_v A_v(h_{fv}), \tag{47} \] where \(f\) and \(v\) denote the faces and the vertices, \(h_\ell\in SU(2)\) for every link \(\ell\), and (with \(h_\ell=h_{vf}\)), \[ A_v (h_{\ell})= \int_{{}_{SL(2,\mathbb{C})}}\!\!\!\!\! d{g_{e}}' \, \prod_{\ell} \sum_{j}\ d_{j} \ D_{jn\, jm}^{(\gamma j, j)} (g_{e}g_{e' }^{-1}) \ D_{mn}^{(j)}(h_{\ell}), \tag{48} \] the integration being over one \(g_e\) for each node (edge of \(v\)), except one. The product is over 10 faces \(f\) per each vertex, and \(D^{(j)}\) and \(D^{(p\, k)}\) are matrix elements of the \(SU(2)\) and \(SL(2,\mathbb{C})\) representations. Calculating those amplitudes for the considered case leads to clear evidences that the bouncing time is indeed of the order of \(M^2\). On should notice that these equations are generic for a spin foam amplitude. The calculation of this phenomenon cannot be done over a background because it is a quantum tunnelling and it is therefore doable in covariant LQG.

## Lorentz invariance violation

It has been argued that the discreetness of space, understood as a genuine property of quantum geometry, could be tested by searching for a Lorentz invariance violation in vacuum. Counter arguments were also given, saving Lorentz invariance, in the sense that a boosted frame might lead to a different mean expectation value for the length operator but not to a different spectrum. The situation is therefore unclear. Several studies (see Girelli et al. 2012) were however devoted to possible Lorentz violations in LQG. The phenomenological studies are based on an energy-momentum relation modified by \(E\simeq p+m^2/2p\pm\xi(E^2/M_{QG})^n\) where \(M_{QG}\) is the energy scale at which quantum gravity begins to play an important role, \(\xi>0\), and \(n\) is an integer.

This idea can be implemented in several frameworks, from non-commutative space-times to effective field theories and non-linear Poincaré symmetries. In LQG, there is an additional helicity dependance than can be added to the previously given relation. It is possible to carry out calculations by analyzing the Hamitonian of the electromagnetic field in a semi-classical state which is a discrete approximation to the flat geometry. As the densitized triad operator enters the Hamiltonian, its expectation value on this state will receive LQG corrections. The associated modified dispersion relation for photons becomes \(\omega^2_{\pm}=k^2\mp4\chi k^3/M_{Pl}\) with \(\chi\sim1\). Photons would be submitted to birefringence in vacuum.

It might be possible to detecting chiral gravity with the pure pseudospectrum reconstruction of the cosmic microwave background polarized anisotropies (see Ferté & Grain 2014).

Other studies are based on a different approach. Classically, the action \(S[A]=\int_\Sigma\mathcal{S}[A]\) can be used to define a slicing. Quantum mechanically, this slicing fluctuates. The calculation was performed using a Born-Oppenheimer state \(\Psi_0[A]\chi[A,\phi]\) where \(\Psi_0\) is a semiclassical state. The densitized triad on this state, evaluated around the classical trajectory, is deformed according to \(E^{(0)~a}_{~~~i}(x,t,\omega)=E^{(0)~a}_{~~~i}(x,t)(1-\alpha L_{Pl}\omega)\). As the triad now depends on \(\omega\), this induces an \(\omega\)-dependent metric and thus a modified dispersion relation: \(m^2=\omega^2-k^2/(1-\alpha L_{Pl}\omega)\). Quantum gravity fluctuations might then lead to an effective preferred frame. Quantum fluctuations would lead to a frame \(\tilde{e}_a^\mu\) which is non-linearly related to \(e_a^\mu\) through a \(\pi_\mu\) dependance, \(\tilde{e}_a^\mu=F(e_a^\mu,\pi_\mu)\). As the momenta are measured by \(\tilde{p}_a=\tilde{e}_a^\mu\pi_\mu\), the transformation law would not be given anymore by usual the Lorentz matrices. This leads to a deformation of the Poincaré algebra.

There are also indications that different choices of the Immirzi parameter can, in some cases, lead to different outcomes for the modifications of the dispersion relations, depending on the quantization scheme chosen (see Brahma et al. 2016). This might allowsone to differentiate between the quantization schemes via testable phenomenological predictions.

## References

- Rovelli, Carlo (1996). "Black hole entropy from loop quantum gravity".
*Phys. Rev. Lett.*77: 3288-3291. doi:"10.1103/PhysRevLett.77.3288.

- Bojowald, Martin (2001). "Absence of singularity in loop quantum cosmology".
*Phys. Rev. Lett.*86: 5227-5230. doi:10.1103/PhysRevLett.86.5227.

- Bojowald, Martin (2003).
*"Quantum gravity and the big bang"*.

- Ashtekar, Abhay and Krishnan, Badri (2004). "Isolated and dynamical horizons and their applications".
*Living Rev. Rel.*7: 10. doi:10.12942/lrr-2004-10.

- Ashtekar, Abhay and Lewandowski, Jerzy (2004). "Background independent quantum gravity: A Status report".
*Class. Quant. Grav.*21: R53. doi:10.1088/0264-9381/21/15/R01.

- Bojowald, Martin (2004). "Loop quantum cosmology: Recent progress".
*Pramana*63: 765-776. doi:10.1007/BF02705198.

- Bojowald, Martin (2005). "Loop quantum cosmology".
*Living Rev. Rel.*8: 11. doi:10.12942/lrr-2005-11.

- Ashtekar, Abhay and Bojowald, Martin (2005). "Black hole evaporation: A Paradigm".
*Class. Quant. Grav.*22: 3349-3362. doi:10.1088/0264-9381/22/16/014.

- Bojowald, Martin (2006). "Loop quantum cosmology and inhomogeneities".
*Gen. Rel. Grav.*38: 1771-1795. doi:10.1007/s10714-006-0348-4.

- Ashtekar, Abhay and Pawlowski, Tomasz and Singh, Parampreet (2006). "Quantum Nature of the Big Bang: Improved dynamics".
*Phys. Rev.*D74: 084003. doi:10.1103/PhysRevD.74.084003.

- Bojowald, Martin and Hossain, Golam Mortuza and Kagan, Mikhail and Shankaranarayanan, S. (2008). "Anomaly freedom in perturbative loop quantum gravity".
*Phys. Rev.*D78: 063547. doi:10.1103/PhysRevD.78.063547.

- Ashtekar, Abhay and Wilson-Ewing, Edward (2009). "Loop quantum cosmology of Bianchi I models".
*Phys. Rev.*D79: 083535. doi:10.1103/PhysRevD.79.083535.

- Ashtekar, Abhay and Kaminski, Wojciech and Lewandowski, Jerzy (2009). "Quantum field theory on a cosmological, quantum space-time".
*Phys. Rev.*D79: 064030. doi:10.1103/PhysRevD.79.064030.

- Bianchi, Eugenio and Rovelli, Carlo and Vidotto, Francesca (2010). "Towards Spinfoam Cosmology".
*Phys. Rev.*D82: 084035. doi:10.1103/PhysRevD.82.084035.

- Mielczarek, Jakub and Cailleteau, Thomas and Grain, Julien and Barrau, Aurelien (2010). "Inflation in loop quantum cosmology: dynamics and spectrum of gravitational waves".
*Phys. Rev.*D81: 104049. doi:10.1103/PhysRevD.81.104049.

- Rovelli, Carlo and Vidotto, Francesca (2013). "Evidence for Maximal Acceleration and Singularity Resolution in Covariant Loop Quantum Gravity".
*Phys. Rev. Lett.*111: 091303. doi:10.1103/PhysRevLett.111.091303.

- Ashtekar, Abhay and Singh, Parampreet (2011). "Loop Quantum Cosmology: A Status Report".
*Class. Quant. Grav.*28: 213001. doi:10.1088/0264-9381/28/21/213001.

- Ashtekar, Abhay and Sloan, David (2011). "Probability of Inflation in Loop Quantum Cosmology".
*Gen. Rel. Grav.*43: 3619-3655. doi:10.1007/s10714-011-1246-y.

- Rovelli, Carlo (2011). "Zakopane lectures on loop gravity".
*PoS*QGQGS2011: 003.

- Agullo, Ivan and Ashtekar, Abhay and Nelson, William (2012). "A Quantum Gravity Extension of the Inflationary Scenario".
*Phys. Rev. Lett.*109: 251301. doi:10.1103/PhysRevLett.109.251301.

- Singh, Parampreet (2012). "Numerical loop quantum cosmology: an overview".
*Class. Quant. Grav.*29: 24400.

- Wilson-Ewing, Edward (2012). "Holonomy Corrections in the Effective Equations for Scalar Mode Perturbations in Loop Quantum Cosmology".
*Class. Quant. Grav.*29: 085005. doi:10.1088/0264-9381/29/8/085005.

- Wilson-Ewing, Edward (2012). "Lattice loop quantum cosmology: scalar perturbations".
*Class. Quant. Grav.*29: 215013. doi:10.1088/0264-9381/29/21/215013.

- Cailleteau, Thomas and Mielczarek, Jakub and Barrau, Aurelien and Grain, Julien (2012). "Anomaly-free scalar perturbations with holonomy corrections in loop quantum cosmology".
*Class. Quant. Grav.*29: 095010. doi:10.1088/0264-9381/29/9/095010.

- Girelli, Florian and Hinterleitner, Franz and Major, Seth (2012). "Loop Quantum Gravity Phenomenology: Linking Loops to Observational Physics".
*SIGMA*8: 098. doi:10.3842/SIGMA.2012.098.

- Agullo, Ivan and Ashtekar, Abhay and Nelson, William (2013). "The pre-inflationary dynamics of loop quantum cosmology: Confronting quantum gravity with observations".
*Class. Quant. Grav.*30: 085014. doi:10.1088/0264-9381/30/8/085014.

- Agullo, Ivan and Ashtekar, Abhay and Nelson, William (2013). "Extension of the quantum theory of cosmological perturbations to the Planck era".
*Phys. Rev.*D87 (4): 043507. doi:10.1103/PhysRevD.87.043507.

- Ashtekar, Abhay (2013). "Introduction to loop quantum gravity and cosmology".
*Lect. Notes Phys.*863: 31-56. doi:10.1007/978-3-642-33036-0_2.

- Bojowald, Martin (2013). "Mathematical structure of loop quantum cosmology: Homogeneous models".
*SIGMA*9: 082. doi:10.3842/SIGMA.2013.082.

- Linsefors, Linda and Barrau, Aurelien (2013). "Duration of inflation and conditions at the bounce as a prediction of effective isotropic loop quantum cosmology".
*Phys. Rev.*D87 (12): 123509. doi:10.1103/PhysRevD.87.123509.

- Barrau, A. and Cailleteau, T. and Grain, J. and Mielczarek, J. (2014). "Observational issues in loop quantum cosmology".
*Class. Quant. Grav.*31: 053001. doi:10.1088/0264-9381/31/5/053001.

- Barrau, Aurélien and Rovelli, Carlo and Vidotto, Francesca (2014). "Fast Radio Bursts and White Hole Signals".
*Phys. Rev.*D90 (12): 127503. doi:10.1103/PhysRevD.90.127503.

- Diener, Peter and Gupt, Brajesh and Singh, Parampreet (2014). "{Numerical simulations of a loop quantum cosmos: robustness of the quantum bounce and the validity of effective dynamics".
*Class. Quant. Grav.*31: 105015. doi:10.1088/0264-9381/31/10/105015.

- Ferté, A. and Grain, J. (2014). "Detecting chiral gravity with the pure pseudospectrum reconstruction of the cosmic microwave background polarized anisotropies".
*Phys. Rev.*D89 (10): 103516. doi:10.1103/PhysRevD.89.103516.

- Ghosh, Amit and Noui, Karim and Perez, Alejandro (2014). "Statistics, holography, and black hole entropy in loop quantum gravity".
*Phys. Rev.*D89 (8): 084069. doi:10.1103/PhysRevD.89.084069.

- Linsefors, Linda and Barrau, Aurelien (2014). "Modified Friedmann equation and survey of solutions in effective Bianchi-I loop quantum cosmology".
*Class. Quant. Grav.*31: 015018. doi:10.1088/0264-9381/31/1/015018.

- Alesci, Emanuele and Cianfrani, Francesco (2015). "Loop quantum cosmology from quantum reduced loop gravity".
*Europhys. Lett.*111 (4): 40002. doi:10.1209/0295-5075/111/40002.

- Ashtekar, Abhay and Barrau, Aurelien (2015). "Loop quantum cosmology: From pre-inflationary dynamics to observations".
*Class. Quant. Grav.*32 (23): 234001. doi:10.1088/0264-9381/32/23/234001.

- Barrau, Aurelien and Cao, Xiangyu and Noui, Karim and Perez, Alejandro (2015). "Black hole spectroscopy from Loop Quantum Gravity models".
*Phys. Rev.*D92 (12): 124046. doi:10.1103/PhysRevD.92.124046.

- Linsefors, Linda and Barrau, Aurelien (2015). "Exhaustive investigation of the duration of inflation in effective anisotropic loop quantum cosmology".
*Class. Quant. Grav.*32 (3): 035010. doi:10.1088/0264-9381/32/3/035010.

- Gomar, Laura Castelló and Martín-Benito, Mercedes and Marugán, Guillermo A. Mena (2015). "Gauge-Invariant Perturbations in Hybrid Quantum Cosmology".
*JCAP*1506 (06): 045. doi:10.1088/1475-7516/2015/06/045.

- Haggard, Hal M. and Rovelli, Carlo (2015). "Quantum-gravity effects outside the horizon spark black to white hole tunneling".
*Phys. Rev.*D92 (10): 104020. doi:10.1103/PhysRevD.92.104020.

- "Elizaga Navascués, Beatriz and Martín-Benito, Mercedes and Mena Marugán, Guillermo A. (2015). "Modeling effective FRW cosmologies with perfect fluids from states of the hybrid quantum Gowdy model".
*Phys. Rev.*D91 (2): 024028. doi:10.1103/PhysRevD.91.024028.

- Barrau, Aurelien and Bojowald, Martin and Calcagni, Gianluca and Grain, Julien and Kagan, Mikhail (2015). "Anomaly-free cosmological perturbations in effective canonical quantum gravity".
*JCAP*1505 (05): 051. doi:10.1088/1475-7516/2015/05/051.

- Barrau, Aurélien and Bolliet, Boris and Vidotto, Francesca and Weimer, Celine (2016). "Phenomenology of bouncing black holes in quantum gravity: a closer look".
*JCAP*1602 (02): 022. doi:10.1088/1475-7516/2016/02/022.

- Bojowald, Martin and Mielczarek, Jakub (2015). "Some implications of signature-change in cosmological models of loop quantum gravity".
*JCAP*1508 (08): 052. doi:10.1088/1475-7516/2015/08/052.

- Barrau, Aurélien (2016). "Evaporation Spectrum of Black Holes from a Local Quantum Gravity Perspective".
*Phys. Rev. Lett.*117 (27): 271301. doi:10.1103/PhysRevLett.117.271301.

- Bolliet, Boris and Barrau, Aurélien and Grain, Julien and Schander, Susanne (2016). "Observational exclusion of a consistent loop quantum cosmology scenario".
*Phys. Rev.*D93: 124011. doi:10.1103/PhysRevD.93.124011.

- Christodoulou, Marios and Rovelli, Carlo and Speciale, Simone and Vilensky, Ilya (2016). "Planck star tunneling time: An astrophysical relevant observable from background-free quantum gravity".
*Phys. Rev.*D94 (8): 084035. doi:10.1103/PhysRevD.94.084035.

- Haggard, Hal M. and Rovelli, Carlo (2016).
*"Quantum Gravity Effects around Sagittarius A*"*.**D25**. pp. 1644021. doi:10.1142/S0218271816440211.

- Schander, Susanne and Barrau, Aurélien and Jolliet, Boris and Linsefors, Linda and Mielczarek, Jakub and Grain, Julien (2016). "Primordial scalar power spectrum from the Euclidean Big Bounce".
*Phys. Rev.*D93 (2): 023531. doi:10.1103/PhysRevD.93.023531.

- Brahma, Suddhasattwa and Ronco, Michele and Amelino-Camelia, Giovanni and Marciano, Antonino (2017). "Linking loop quantum gravity quantization ambiguities with phenomenology".
*Phys. Rev.*D95 (4): 044005. doi:10.1103/PhysRevD.95.044005.

- Oriti, Daniele and Sindoni, Lorenzo and Wilson-Ewing, Edward (2017). "Bouncing cosmologies from quantum gravity condensates".
*Class. Quant. Grav."*34 (4): 04LT0.