Unitarity and information in quantum gravity: a simple example

In approaches to quantum gravity, where smooth spacetime is an emergent approximation of a discrete Planckian fundamental structure, any effective smooth field theoretical description would miss part of the fundamental degrees of freedom and thus break unitarity. This is applicable also to trivial gravitational field (low energy) idealizations realized by the use of the Minkowski background geometry which, as any other spacetime geometry, corresponds, in the fundamental description, to infinitely many different and closely degenerate discrete microstates. The existence of such microstates provides a large reservoir for information to be coded at the end of black hole evaporation and thus opens the way to a natural resolution of the black hole evaporation information puzzle. In this paper we show that these expectations can be made precise in a simple quantum gravity model for cosmology motivated by loop quantum gravity. Concretely, even when the model is fundamentally unitary, when microscopic degrees of freedom irrelevant to low-energy cosmological observers are suitably ignored, pure states in the effective description evolve into mixed states due to decoherence with the Planckian microscopic structure. Moreover, in the relevant physical regime these hidden degrees freedom do not carry any `energy' and thus realize in a fully quantum gravitational context the idea (emphasized before by Unruh and Wald) that decoherence can take place without dissipation, now in a concrete gravitational model strongly motivated by quantum gravity. All this strengthen the perspective of a quite conservative and natural resolution of the black hole evaporation puzzle where information is not destroyed but simply degraded (made unavailable to low energy observers) into correlations with the microscopic structure of the quantum geometry at the Planck scale.


INTRODUCTION
The mathematical models that so far define our successful physical theories are all reversible in the sense that they can predict the future value of the variables they use from their initial values, while conversely the past can be uniquely reconstructed from the values of these variables in the future. The memory of the initial conditions is not lost in the dynamics and their information content remains. This is true for classical mechanics and field theory and it is also true for quantum mechanics and quantum field theory as long as we do not invoke the postulate of the collapse of the wave function (i.e., as long as we do not intervene from the outside via a measurement 1 ). In the quantum mechanical setting, this property boils down to the fact that evolution to the future is given by a unitary operator which can always be undone via its adjoint transformation.
This property of our fundamental models has always troubled naive intuition when faced with situations that appear to be irreversible. For example, what would happen to these words if the computer collapsed at this very moment? What if, after being printed, this paper is burned? Common sense would answer that the information in these pages (if of any relevance) would be lost. However, the physicist, trained to firmly believe in the statement of the previous paragraph, would say that the information in these words is not lost but simply hidden (to the point of becoming unrecoverable) in the humongous number of microscopic variables that would describe the whole system. In the case of burning the paper, these words remain "written" (it would be claimed) in the multiple correlations between the degrees of freedom of the molecules in the gas of the combustion diffusing in the atmosphere while transferring the information to even larger and yet pristine portions of the very large phase space of an unbounded universe. In the case where the computer collapses, a similar story can be told involving the dissipation of the bits into the environment. Of course the physicist cannot prove this; however, it is a consistent story in view of the strongly cherished principle of unitarity.
Such effective irreversibility is clearly captured in the second law of thermodynamics stating that (for suitable situations involving a large number of degrees of freedom) entropy can only increase. At the classical level this clashes at first sight with Liouville's theorem stating that the phase space volume of the support of a probability distribution is preserved by dynamical evolution. However, nothing restricts the shape of this volume to evolve into highly intricate forms that a macroscopic observer might be unable to resolve. More precisely, suitable initial conditions that the observer agent regards as special (for instance the macroscopic configurations of ink particles defining words in this paper before the fire reached them) come with an uncertainty in accordance to the observers limited measurement capabilities. This is idealized by a distribution in phase space occupying an initial phase space volume of a regular shape (this ensemble of points represents the system in what follows). Now as time goes by the apparent phase space volume (not the real volume which remains constant) would seem to grow to the agent just because of its intrinsic inability to separate the points in phase space that the systems occupies from the close neighboring ones where the system is not. The arrow of time (characterizing large systems) is only emergent macroscopically due to the special initial conditions, and the intrinsic coarse graining introduced by a macroscopic observer with its limitations. We will argue that the general lines of this story remain the same when black hole evaporation is considered.
General relativity combined with quantum field theory, in a regime where both are expected to be good approximations, imply that large isolated black holes behave like thermodynamical systems in equilibrium. They are objects close to equilibrium at the Hawking temperature that lose energy extremely slowly via Hawking radiation. When perturbed they come back to equilibrium to a new state and the process satisfies the first law of thermodynamics with an entropy equal to 1/4 of the area A of the black hole horizon in Planck units. Under such perturbation (which in particular can also be associated with their slow evaporation), the total entropy of the universe can only increase namely where δS matter represents the entropy of whatever is outside the black hole (including, for instance, the emitted radiation). This quasi equilibrium phase-which is extremely long lasting for macroscopic black holes but, at the same time, is only an intermediate situation before complete evaporation-is predicted by general relativity as the result of gravitational collapse taking place for suitable initial conditions. Indeed in order to make a black hole (BH), the past must be special (low entropy) too. Thus, the irreversibility captured by (1) can once more be associated with the same ingredients present in our previous example: the special nature of the initial conditions (low curvature and low densities in the past), high curvature, and huge new phase space regions available in the future; more precisely, in the Planckian regime that the singularity theorems of general relativity predict to develop inside the black hole horizon (the would-be-singularity from now on).
Such a perspective resonates with the one emphasized by Penrose (see for instance Penrose, 2005), among others: in full generality (now including gravity) the arrow of time comes from the fact that the universe is special in the past with a spacetime that was well-approximated by a homogeneous geometry and matter distribution (gas and dust) with tiny perturbations that would eventually grow and form galaxies and stars that one day can collapse to form black holes 2 . Before the formation of a black hole, the story of our system exploring larger and larger portions of the available phase space is the usual standard involving molecules, atoms, and fundamental particles. The perspective we put forward here is that the story continues to be the same when black holes form and a new huge portion of phase space is opened by the physics of gravitational collapse. This new channel for entropy growth is opened by the appearance of the internal would-be-singularity of the classical description beyond the event horizon. The gravitational collapse ignites interactions with the Planckian regime inside the black hole horizon (see footnote 3), like the lighter setting the paper on fire and thus degrading the ink in the words when burning the paper in our previous discussion. The singularity predicted by general relativity brings the system in contact with the quantum gravity scale, and (as in 2 To these two specialty conditions one might also have to add one concerning the state of the hypothetical microscopic constituents at the Planck scale if the view we are advocating here and in Perez (2015Perez ( , 2017 is correct. the burning paper) this provides the key for resolving the puzzle of information in black hole evaporation.
This perspective, first advocated in Perez (2015Perez ( , 2017, can be described with the help of Figure 1 (which applies to the general black hole formation and evaporation background scenario of reference; Ashtekar and Bojowald, 2005). The first assumption in the diagram is that there is evolution across the wouldbe-singularity (predicted by the classical dynamics) inside the black hole. This assumption is intrinsic in the representation in the figure; however, the scenario still makes sense if instead a baby universe is formed, i.e., if the would-be-singularity remains causally disconnected from the outside after evaporation. In that case, the correlations established with the baby universe remain hidden forever to outside observers. The virtue of the present scenario in such a case would be to give an identity to the degrees of freedom involved. The idea that the spacetime representation of the situation resembles the one in Figure 1 comes from the various results in symmetry-reduced models for both cosmology (Bojowald, 2001;Ashtekar et al., 2006) as well as for black holes (Modesto, 2004(Modesto, , 2006Bojowald, 2005;Bojowald et al., 2005;Gambini and Pullin, 2013;Corichi and Singh, 2016;Ashtekar et al., 2018), and was first pictured in Ashtekar and Bojowald (2005). In such a context, a "scattering theory" representation (where an in-state evolves into an out-state) is possible even though the result (as we will argue) cannot be translated into the language of effective quantum field theory.
But what do we mean by a black hole in this evaporating context? In the asymptotically flat idealization, the black hole region is defined in classical general relativity as the portion of the spacetime M that is not part of the past of I + . Such a definition needs to be modified in quantum gravity. In order to do that, we introduce the notion of the semiclassical past J − C (I + ) of I + as the collection of events in the spacetime that can be connected to I + by causal curves along which the Kretschmann scalar K ≡ R abcd R abcd ≤ Cℓ −4 p for some constant C of order unity. The black hole region can then be defined as Its dependence on the constant C is not an important limitation in the discussion about information. Different C would lead to BH regions that coincide with Planckian corrections. The most clear physical picture emerges from the analysis of the Penrose diagram on the left panel of Figure 1, from the point of view of observers at future null infinity I + . These observers are assumed to be at the center of the mass Bondi frame of the BH formed via gravitational collapse. We also assume that the Bondi mass of the BH is initially M ≫ m p at some delayed time u on I + representing the time where the BH has achieved its quasi-equilibrium state and starts evaporating slowly via Hawking radiation, i.e., the BH is initially macroscopic. Under such conditions the evaporation is very slow and we can trust the semiclassical description that suggests that the Bondi mass M(u) will slowly decrease with u from this initial value M until time u = u 0 (see figure) with M(u 0 ) m p in a time of the order of M 3 in Planck units. From this time on, the details depend on a full quantum gravity calculation because the curvature around the BH horizon has become Planckian. Nevertheless, independently of such details we can safely say that the spacetime and the matter degrees of freedom encoded on I + for u > u 0 must be in a superposition of states, all of which are very close to flat spacetime, as far as their geometry is concerned, with matter excitations very close to the vacuum because there is only at best an energy of the order of E late ≈ m p to substantiate both. In addition these excitations must be correlated with the early Hawking radiation with energy E early ≈ M − m p if unitarity is to hold. The late degrees of freedom are often referred to as purifying degrees of freedom.
One possibility is to assume that such purifying degrees of freedom are particle excitations coming from what is left of the BH (a remnant). Now, due to the fact that these particles must be extremely low-energy particles as only a total energy of E late ≈ m p is available for purification, a simple estimate of the time (denoted as τ p ) that the process would have to last if this is the main channel for purification yields τ p ≈ (M/m p ) 4 . This is the scenario of an extremely long lasting point-particle-like remnant with a huge internal degeneracy which is claimed to be problematic from the point of view of effective quantum field theory (Banks et al., 1993).
Instead we propose a different alternative: if smooth spacetime and matter fields are emergent notions from underlying discrete microscopic physics, then the coarse low-energy notion of classical geometry with smooth fields living on it would correspond in the fundamental Hilbert space to a very large set with (possibly) an infinite number of microscopic states. For instance the Minkowski vacuum unicity in standard quantum field theory would fail in the sense that the requirement that states look flat for (coarse-grained) low-energy observers-which are those for which an effective quantum field theory description in terms of smooth fields living on a smooth geometry is a suitable approximation-would still admit a highly degenerate ensemble of microscopic states (all states with total mass indistinguishable from zero by these macroscopic observers). Now, such lowenergy modes cannot be identified with effective field theoretical excitations as the infrared excitations of fields mentioned in the previous paragraph (say low-energy photons). Like the molecular structure that escapes the smooth characterization of the Navier-Stokes effective theory of fluids, the degrees of freedom of interest here correspond to defects in the Planckian fabric of quantum gravity bound to be missed by coarse low-energy agents and their effective field theory mathematical models based on smoothness.
Why should information be hidden in the UV but not in the IR modes as in the remnant scenario mentioned above? It is often believed that because the volume inside the black hole actually becomes very large (according to suitable definitions; Christodoulou and De Lorenzo, 2016) then modes that are correlated with the Hawking radiation are redshifted and become highly IR inside. Although this is true for spherically symmetric Hawking quanta in the spherically symmetric Schwarzshild background-where such modes are indeed infinitely redshifted as detected by regular observers when they approach the singularity at r = 0-this conclusion fails when one considers no-spherical modes no matter how small the deviation from spherical symmetry is 3 . Therefore, generically all modes become UV close to the singularity [this is the central weakness of the bag-of-gold scenarios and the perspective proposed in Ashtekar (2020)].
We can draw a formal analogy with the Unruh effect as follows. The Unruh effect arises from the structure of the vacuum state |0 of a quantum field on the Minkowski spacetime when written in terms of the modes corresponding to Rindler accelerated observers with their intrinsic positive frequency notion. The vacuum takes the form where |n, k L and |n, k R define the particle modes-as viewed by an accelerated observer with uniform acceleration a-on the left and the right of the Rindler wedge (Wald, 1995). Here we see from the form of the previous expansion that even when we are dealing with a pure state (if we define the density matrix |0 0|), the reduced density obtained by tracing over one of the two wedges would produce a thermal state with T = a/(2π). The statement in the perspective we propose on the purification of information in BH evaporation can be schematically represented (the following is certainly not a precise equation) by where an initial state of a flat quantum geometry |flat, 0 tensor product with a state representing initially diluted matter fields |φ evolves unitarily via U into the formation of a BH and the subsequent evaporation (Figure 1) which after complete evaporation is written as a superposition of flat quantum geometry states |flat, n -which are all indistinguishable from |flat, 0 to low-energy agents and differ among them by quantum numbers n corresponding to quantities that are only measurable if one probes the state down to its Planckian structuretensor product with normal n-particle excitations of matter 3 In the Schwarzshild background, the frequency measured by a radially freely falling observer normal to the r =constant hypersurfaces goes like where E = −k · ξ and ℓ = k · ψ are the conserved quantities associated with the massless particle with wave vector k a , and ξ a and ψ a are the stationarity and rotation killing fields of the background. The qualitative behavior approaching r = 0 would be the same for any other observer measuring ω (the divergence of ω is observer-independent). Only exactly spherically symmetric modes with ℓ = 0 would become IR at the singularity. However, this conclusion is no longer true if the BH rotates or if we consider that at the fundamental level, states with exact spherical symmetry inside the BH are of measure zero. Notice that such UV divergence in the non-spherically symmetric Hawking partners implies large deviations from spherical symmetry near the singularity (if their back reaction would be taken into account). This should be kept in mind when modeling the situation with spherically symmetric mini superspace quantum gravity models. fields representing Hawking radiation. As mentioned above, the previous equation is only a schematic. Its main inappropriateness is the fact that the reduced density matrix obtained by tracing over the quantum geometry hidden degrees of freedom would give a thermal state at a fixed temperature T. This is at odds with the expectation that the Hawking radiation should contain a superposition of the thermal radiation emitted at different Hawking temperatures during the long history of the evaporation of the BH. But the point that this equation and the discussion of the previous paragraph should make clear is that the purification mechanism proposed here has nothing to do with the pointlike remnant scenario with all its problems associated with a long lasting particle-like remnant. In the present scenario, to the future of the would-be-singularity in Figure 1, we simply have a quantum superposition of different quantum geometry states that all look flat to low-energy observers. There are no localized remnants hiding in the huge internal degeneracy; there is only a large superposition of states that are inequivalent in the fundamental quantum gravity Hilbert space but seem identical when tested with low-energy probes. Such degrees of freedom cannot be captured by any effective description in terms of smooth fields (EQFT) for the simple reason that they are discrete in their fundamental nature.
Notice that the degrees of freedom, where information would be coded after BH evaporation, do not satisfy the usual Einstein-Planck relationship E =hω or equivalently p = h/λ (for some "wavelength" λ or "frequency" ω), and this might deceive intuition 4 . These are Planckian defects yet they do not carry Planckian energy. The point is that such a relationship only applies under suitable conditions which happen to be met in many cases but do not need to always be valid. One case is the one of degrees of freedom that can be thought of as waves moving on a preexistent spacetime. This is the case for particle excitations in the Fock space of quantum field theory or effective quantum field theories (both of which are defined in terms of a preexistent spacetime geometry). There is no clear meaning to the above intuitions in the full quantum gravity realm where the present discussion is framed. Even when such relations (linked to the usual uncertainty principle of quantum mechanics) should hold in a suitable sense-if the structure suggested by canonical quantization survives in quantum gravity as it should to a certain degree-they would apply to phase space variables encoding to 4 A nice counter example of this intuition is given by the case of a non-relativistic charged particle in a two-dimensional infinite perfect conductor in a uniform magnetic field normal to the conducting plane. The energy eigenvalues are given by the Landau levels E n =hω B (n + 1/2) where ω B = qB/(mc) is the Bohr magneton frequency, but they are infinitely degenerate. There are canonically conjugated variables (P, Q) associated with the particle that are cyclic, i.e., do not appear in the Hamiltonian. In this case, one can produce wave packets that are as "localized" in the variable Q without changing the energy of the system. Interestingly, this is a perfect example of a system where one could have an apparent loss of information of the type we are proposing here (for a more realistic analog gravity model discussing the information paradox along the lines of the present scenario; see Liberati et al., 2019). If one scatters a second particle interacting softly with the charged particle on the plate so that the interaction does not produce a jump between different Landau levels, then correlations with the cyclic variables would be established without changing the energy of the system. This is the perfect model to illustrate the possibility of decoherence without dissipation.  (Perez, 2015). The shaded region represents the would-be-singularity where high fluctuations in geometry and fields are present and where the low-energy degrees of freedom of the Hawking pairs are forced to interact with the fundamental Planckian degrees of freedom. (Right) Same situation as a scattering process from an initial to a final Cauchy hypersurface. This figure contains basically the same information as the Penrose diagram. Its additional merit, if any, is the more intuitive representation of the shrinking black hole as well as the time scales involved (collapse time is very short with respect to the evaporation time). Both these features are absent in the conformal representation on the left. There are limitations to this kind of spacetime representation of a process that is fundamentally quantum and hence only understandable in terms of superpositions of different spacetime geometries.
FIGURE 2 | Diagram (effectively) illustrating the natural scenario, suggested by the fundamental features of LQG, for the resolution of the information puzzle in black hole evaporation (Perez, 2015). As in Figure 1, one should keep in mind the limitations of this kind of spacetime representation of a process that is fundamentally quantum and hence only understandable in terms of superpositions of different spacetime geometries.
the true degrees of freedom of gravity that we expect (from general covariance) to be completely independent of a preexistent background geometry. We will see that such degrees of freedom with such a peculiar nature actually arise naturally in the toy model of quantum gravity that we analyze in this article.
It is presently hard to prove that such a scenario is viable in a quantum theory of gravity simply because there is no such theoretical framework developed enough for tackling BH formation and evaporation in detail. However, the application of loop quantum gravity to quantum cosmology leads to a model with similar features, where evolution across the classical singularity is well defined (Bojowald, 2001). The results have been reported in Amadei and Perez (2019). In this article we present the main features of this model in more detail and show that the conclusions of Amadei and Perez (2019), some of which are drawn from some simplified models of matter coupling, are generic and remain true in more physically realistic models.
The results can be described briefly by making reference to Figure 2 which should be compared with the right panel in Figure 1. We will show that the evolution in loop quantum cosmology from a universe in an initially contracting state in the past of the would-be-singularity to an expanding universe in its future is perfectly unitary in its fundamental description. Nevertheless, states in the Hilbert space of loop quantum cosmology contain quantum degrees of freedom which are hidden to low-energy coarse-grained observers. If these degrees of freedom are traced out in the initial density matrix then we will see that pure states (in the sense of the reduced density matrix) generically evolve into mixed states across the wouldbe-singularity. Information is lost in correlations with degrees of freedom that are Planckian and thus inaccessible to macroscopic observers. These correlations are established in an inevitable way during the strong curvature phase of evolution across the big bang (just as expected in the BH scenario described above). As energy is conserved (energy is a delicate notion in cosmology but happens to be well-defined in our model as we will see), the defects that purify the final state do not enter into the energy balance which realizes another crucial necessary ingredient of the general scenario (decoherence happens with negligible dissipation; Unruh, 2012).
The paper is organized as follows. In the first part (section 2), we show that the scenario we have described in general terms so far is realized in unimodular quantum cosmology following the standard quantization prescription of loop quantum cosmology. Aside from a different choice of time variable, the model of this section is exactly equivalent to other models studied in the standard literature (Ashtekar and Singh, 2011). In the second part of the paper, we observe that there is natural extension of loop quantum cosmology based on the regularization ambiguity associated with the quantization of the Hamiltonian. This extension provides another interesting realization of our mechanism. Although this second option is not necessary to illustrate our point (already realized in the standard theory in the first part), it gives a different identity to the defects which could lead to independent and thus useful insights. We have included a series of appendices where some calculations are shown. Appendix D is especially important because some of the over-simplified model (analytic) calculations in the body of the paper are completed (numerically) using a more physically realistic case of the coupling of gravity with a massless scalar field.

UNIMODULAR GRAVITY: FOUR-VOLUME AS TIME
In this section we introduce the basic structure of quantum unimodular gravity in the framework of quantum cosmology. This theory will provide us with a toy model to study the unitary evolution of the state of the universe across the big-bang singularity. We will see that the Hilbert space of loop quantum cosmology contains the type of microscopic degrees of freedom evoked in the general discussion of the introduction. This is a minimalistic model where our scenario can be explicitly realized.
Unimodular gravity as a concept is nearly as old as general relativity itself, it was introduced by Einstein in 1919 (Einstein, 1919) as an attempt to describe nuclear structure geometrically. In his work, Einstein also identified an appealing feature of the theory, which is the fact that the cosmological constant arises as a dynamical constant of motion that needs to be added to the initial values of the theory. In unimodular gravity, the cosmological constant is a constant of integration and not a universal or fundamental constant of nature. Interest in the theory was regained in the late 80's with the observation of Weinberg (1989) that, for the above reason, semiclassical unimodular gravity provides a trivial resolution of the cosmological constant problem as vacuum energy simply does not gravitate. Unimodular gravity is the natural low-energy description that emerges from the formal thermodynamical ideas of Jacobson (1995), and represents the expected low energy regime of the causal set approach (Bombelli et al., 1987).
Another property of unimodular gravity (specially important for us here) is that it completely resolves the problem of time (Unruh, 1989) in the cosmological FLRW context. More precisely, the theory comes with a preferred time evolution and a preferred Hamiltonian (the energy of the universe is well defined and directly linked with the value of the cosmological constant). The quantum theory is described by a Schrodinger-like equation where states of the universe are evolved by a unitary evolution operator. Therefore, unlike the general situation in quantum gravity, the notion of unitarity is unambiguously defined in unimodular quantum cosmology. This is the main reason why unimodular gravity provides the perfect framework for the discussion of the central point of this work.
Here we specialize in homogeneous and isotropic cosmologies that are spatially flat (k = 0), i.e., the spatial manifold is topologically R 3 . What follows is the standard construction. For a detailed account of the Hamiltonian analysis in the cosmological framework (see Chiou and Geiller, 2010). The FLRW metric is where q ab denotes the fiducial spacial metric. Since is non-compact, some integrals are infrared divergent and are regularized by restricting them to a fixed fiducial cell V of fiducial volume V 0 with respect to the fiducial spacial metric where e i a denotes a fiducial triad and the physical metric is given by q ab = a 2 (t)q ab . The action of unimodular gravity in the FLRW mini superspace setup is given by where λ is a Lagrange multiplier imposing the unimodular constraint N = a −3 (i.e., |g| = 1), and the first term is the Einstein-Hilbert action restricted to the FLRW geometries 5 . In order to use loop quantum cosmology techniques (for a discussion of the quantization in the full loop quantum gravity context; see Smolin, 2009Smolin, , 2011, one introduces the new canonical variables c and p via the basic Ashtekar-Barbero connection variables A i a and E a i , namely where ω i a is a fiducial reference connection. These variables are related to those in (8) via the equations The action becomes and c and p are canonically conjugated in the sense that The unimodular condition N = a −3 fixes the lapse to N = V 0 /|p| 3/2 and the unimodular Hamiltonian becomes The proportionality of the Hamiltonian with V 0 , and the fact that the four-volume bounded by V 0 at two different times is given by v (4) = V 0 t, implies that time evolution can be parameterized in terms of the four-volume elapsed from some reference initial slice. The associated and corresponds to the cosmological constant.

Quantization
The loop quantum cosmology quantization uses a non standard representation of the canonical variables where the variable c does not exist as a quantum operator, and the definition of the Hamiltonian requires a special regularization procedure known as theμ-scheme (Ashtekar and Singh, 2011). The quantization prescription is greatly simplified by the introduction of new canonically conjugated dynamical variables b and ν defined as (Ashtekar et al., 2008) b ≡ c 5 There is an overall minus sign in the definition action with respect to standard treatments. This is done so that the pure-geometry Hamiltonian is positive definite.
with Poisson brackets 6 The variable ν corresponds to the physical volume of the fiducial cell divided by ℓ 2 p ; it has units of distance. The variable b is simply its conjugate momentum. In terms of these variables the gravitational (unimodular) Hamiltonian (13) integrated in a fiducial cell V becomes Note the extreme simplicity of the previous expression: the unimodular Hamiltonian is just the analog of that of a free particle in one dimension with a mass parameter of m = 4πγ 2 /(3V 0 ) and momentum b. In the absence of matter, the Hamiltonian can be quantized in the Wheeler-DeWitt representation where the evolution in unimodular time is unitary and there is no singularity (the classical solutions correspond to De-Sitter geometries with arbitrary but positive cosmological constants). The singularity in the classical theory becomes real when matter is introduced.
In the loop quantum cosmology polymer representation, just as for c, there is no operator corresponding to b but only the operators corresponding to finite ν translations (Ashtekar et al., 2006); from here on referred to as shift operators For k = q √ ℓ p and q ∈ N, states that diagonalize the previous shift operators, denoted |b 0 ; Ŵ ǫ , are labeled by a real value b 0 and by a graph Ŵ ǫ . The graph is a 1D lattice of points in the real line of the form ν = 4n √ ℓ p + ǫ with ǫ ∈ [0, 4 √ ℓ p ) and n ∈ N. The corresponding wave function is given by where the symbol δ Ŵ ǫ means that the wavefunction vanishes when ν / ∈ Ŵ ǫ . It follows from (18) that The states |b; Ŵ ǫ are eigenstates of the shift operators that preserve the lattice Ŵ ǫ . Notice that the eigenvalues are independent of the parameter ǫ, i.e., they are infinitely degenerate and span a non separable subspace of the quantum cosmology Hilbert space H lqc . A scaleμ is needed in order to define a regularization of (17) representing the Hamiltonian in H lqc . The reason is that there are no operators associated with b but only approximants constructed via the shift operators (18). The so-calledμ-scheme (Ashtekar and Singh, 2011) introduces a dynamical length scalē µ defined as 6 The factorh −1 appears on the right hand side of the Poisson brackets due to the introduction ofh (via ℓ 2 p ) in the definition of the new variable ν. This is done to match standard definitions (Ashtekar and Singh, 2011). where represents the so-called "area-gap" which plays the role of a UV regulator. It is normally associated with the smallest non-vanishing area quantum in the full theory of loop quantum gravity. For the moment (as in the standard treatment), this is just a fixed parameter 7 . When translated into the variables (15), µ corresponds to considering approximants to b constructed out of shift operators (18) with fixed k ≡ √ ℓ p . In terms of these, one obtains the following regularization of the Hamiltonian (17) which is a well-defined self-adjoint operator 8 acting on H lqc which coincides with (17) leading to (zero) order in ℓ 2 p . From (14), we obtain an operator associated to the (here dynamical) cosmological constant, namely In the pure gravity case, the cosmological constant is positive definite and bounded from above by the maximum value λ max = 1/(γ 2 ℓ 2 p ). Negative cosmological constant solutions are possible when matter is added (see Appendix B). The states (19) with energy eigenvalues States |b 0 ; Ŵ ǫ are also eigenstates of the cosmological constant with eigenvalue λ (b 0 ) = (8πG)E (b 0 )/V 0 . Notice that the energy eigenvalues do not depend on ǫ ∈ [0, 4 √ ℓ p ). Thus, the energy levels are infinitely degenerate with energy eigenspaces that are non-separable. This is not something peculiar of our model but a general property of the non-standard representation of the canonical commutation relations used in loop quantum cosmology.

On the Interpretation of the ǫ Sectors
It is customary in the loop quantum cosmology literature to restrict to a fixed value of ǫ in concrete cosmological models, 7 In section 3, we will turn this quantity into a quantum operator acting on the microscopic sector of the Hilbert space that will be introduced. 8 The HamiltonianĤ 0 (21) is symmetric, that is 1 ,Ĥ 0 2 = Ĥ 0 1 , 2 , with respect to the inner product 1 , 2 = ν 1 (ν) 2 (ν). The action of the Hamiltonian on (ν) is given by: The key property is 1 (ν), 2 (ν + 2λ) = 1 (ν − 2λ), 2 (ν) where ν is in the support of both 1 (ν) and 2 (ν). This is the statement of the unitarity of the shift operators e −i2λb 1 , 2 = 1 , e i2λb 2 . The symmetric nature of the shift operators appearing in H 0 implies the result. as the dynamical evolution does not mix different ǫ sectors. The terminology "superselected sectors" is used in a loose way in discussions. However, these sectors are not superselected in the strict sense of the term because they are not preserved by the action of all the possible observables in the model, i.e., there are non trivial Dirac observables mapping states from one sector to another. The explicit construction of such observables might be very involved in general (as is the usual case with Dirac observables); nevertheless, it is possible to exhibit them directly at least in one simple situation: the pure gravity case. In that case, the shift operator (18) with shift parameter δ commutes with the pure gravity Hamiltonian (the Hamiltonian constraint if we were using standard loop quantum cosmology) and maps the ǫ sector to the ǫ − 4δ sector. The analogous Dirac observables in a generic matter model can be formally described with techniques of the type used for the definition of evolving constants of motion (Rovelli, 1991). No matter how complicated this might be in practice, the point is well-illustrated by our explicit example in the matter free case 9 .
Thus, different ǫ sectors are not superselected and therefore the infinite degeneracy of the energy eigenvalues of the Hamiltonian (which again we exhibit explicitly in the previous discussion only in the vacuum case) must be taken at face value. How can we understand this large degeneracy from the fact that there would be only a two-fold degeneracy (associated with a contracting or expanding universe) if we had quantized the model using the standard Schrodinger representation or, in other words, the standard Wheeler-DeWitt quantization?
The answer is to be found, we claim, in the notion of coarse graining: low-energy observers only distinguish a two-fold degeneracy for energy (or cosmological constant) eigenvalues: one the universe has a given cosmological constant, and two it is expanding or contracting. These are the quantum numbers in the Wheeler-DeWitt quantization which play a role in our context of the low-energy effective quantum field theory formulation. Such coarse observers are declared to be insensitive to the huge additional degeneracy of energy eigenstates encoded in the quantum number ǫ. All these infinitely many states in the quantum cosmology representation must be considered as equivalent to the two-fold degeneracy mentioned above.
In what follows, and for concreteness, we will consider combinations of states with two different values for ǫ only, i.e., on two different lattices. The idea of the previous paragraph will naturally produce a notion of coarsegraining entropy associated with the intrinsic statistical uncertainty due to the inability for a low-energy agent to distinguish these microscopically orthogonal states. Arbitrary superposition with N different ǫ sectors would lead to similar results [the entropy capacity growing with the usual log(N)]. The N = 2 case treated here renders some explicit calculations straightforward.

Matter Couplings and a Model Capturing Its Essential Features
Here we discuss two simple matter models in order to isolate the generic features of the influence of matter. At the end of the section, we will define a simple and trivially solvable model capturing these features.
Perhaps the simplest matter model that would serve our purposes is minimal and isotropic coupling to a Dirac fermion defined in de Berredo-Peixoto et al. (2012). After symmetry reduction, the action for matter is from which we read the fermionic contribution to the Hamiltonian where ( , ) are the fermionic canonical variables ≡ (V 0 a) 3/2 ψ † and ≡ (V 0 a) 3/2 ψ (Thiemann, 2001), and in the second line, we use the unimodular condition N = V −1 0 a −3 . In the quantum theory, the non trivial anti-commutator is {η, p η } = 1 with the rest equal to zero. This is achieved by writing η = s a s u s e −imt + b † s v s e imt with non trivial anti-commutation relations for the creation and annihilation of operators {a r , a † s } = δ rs = {b r , b † s }, with u s e −imt and v s e imt as a complete basis of solutions of the Dirac equation for positive and negative frequency, respectively (Peskin and Schroeder, 1995). In our model, we can have either the vacuum state, or one or two fermions which saturate the Pauli exclusion principle. If we assume normal ordering, the contribution to the unimodular energy is where n = 0, 1, 2 is the occupation number for the fermion. If instead of the condition N = V −1 0 a −3 , we had used N = 1 (where time is comoving time) then the energy contribution would have been just m for which we have a clear physical intuition: a single fermion homogeneously distributed in the universe contributes to the Hamiltonian with its total mass. The factor 1/a 3 in the previous expression comes from the unimodular condition.
In the case of Wheeler-DeWitt quantization, the contribution of the fermion becomes singular at the big bang a = 0. In loop quantum cosmology, such a quantity remains bounded above due to loop quantum gravity discreteness. Indeed, using the inverse volume quantization given in reference (Ashtekar and Singh, 2011), one haŝ where h F (ν; λ) ≡ 1 4λ 2 |ν + 2λ| We notice that h F (ν; √ ℓ p ) < 1 and decays like 1/ν for ν → ∞ 10 . One could, in principle, add this term to the free Hamiltonian and solve the unimodular time-independent Schrodinger equation Solutions can be interpreted in the sense of scattering theory starting with free wave packets for large ν picked around a value of the cosmological constant (22) or energy (24). The case of coupling with a scalar field is formally very similar, especially in the simplified case where we assume it to be massless. Following Ashtekar and Singh (2011), and using the unimodular condition N = a −3 , we get This leads tô where The momentum p φ commutes with the Hamiltonian and thus it is a constant of motion. If we consider an eigenstate of p φ then the problem reduces again to a scattering problem with a potential decaying like 1/ν 2 when we consider solving the time-independent Schrodinger equation Therefore, both the fermion as well as the scalar field models (which are closer to a possibly realistic scenario) seem tractable with a slight generalization of the standard scattering theory to the discrete loop quantum cosmology setting. However, the main objective in this section is to illustrate an idea in terms of a concrete and simple toy model. With this idea in mind, we will modify the structure suggested by the fermion coupling and the scalar field coupling and simply add an interaction term where the "long-distance interaction" term represented by the function F(ν; λ) is replaced by a short-range analog F(ν; λ) ∝ δ ν,0 . The qualitative properties of the scattering will be the same and the model becomes sufficiently trivial for straightforward analytic computations. The results for the more realistic free scalar field model have been dealt with numerically and are presented in Appendix D.
For that we consider an interaction that begins at ν = 0: 10 There is a great degree of ambiguity in writing the inverse volume operators. Perhaps the simplest is the one introduced in WilsonEwing (2012) that we will actually use in the concrete computations of Appendix D. For more discussion on this see Singh and Wilson-Ewing (2014) and references therein.
FIGURE 3 | The function h φ (ν; λ) evaluated on an ǫ sector containing ν = 0, for p φ = 10 in natural units and λ = 1/2 is plotted using two different ranges. On the left, we see that the function is finite near ν = 0. On the right, we can see that it behaves like −ν −2 for large values of ν. This function can be seen as the effective potential where an asymptotically free state of the universe (pure gravity with cosmological constant state or asymptotically de Sitter state) scatters. If the cosmological constant is negative there are bound states whose superposition can be used to define semiclassical universes oscillating in an endless series of big bangs and big crunches (see section B).
where µ is a dimensionless coupling,Ĥ 0 is given in (58), andĤ int isĤ We have added by hand an interaction Hamiltonian that switches on only when the universe evolves through the would-besingularity at the zero volume state. The key feature of theĤ int is that-as its more realistic relatives, matter Hamiltonians (28) and (32)-it breaks translational invariance and thus, it leads to different dynamical evolution for different ǫ sectors.

Solutions as a Scattering Problem
The scattering problem is very similar to the standard one in onedimensional quantum mechanics; however, one needs to take into account the existence of the peculiar degeneracy of energy eigenvalues contained in the ǫ sectors; see sections 2.1 and 2.2. We will consider, for simplicity, the superposition of only two states supported on two lattices respectively: the lattice Ŵ ǫ with ǫ = 0 for the first one and the one with ǫ = 2 √ ℓ p for the second one. The degenerate eigenstates of the shift operators (19) with eigenvalues exp(i2kb) will be denoted as |b, 1 ≡ |b; Ŵ 0 , and |b, 2 ≡ |b; respectively, while we will denote Ŵ 1 and Ŵ 2 as the corresponding underlying lattices. The immediate observation is that states supported on Ŵ 2 (superpositions of |b, 2 ) will propagate freely because they are supported on a lattice that does not contain the point ν = 0 where the interaction is non trivial. On the other hand, states supported on Ŵ 1 (superpositions of |b, 1 ) will be affected by the interaction at the big bang. Before and after the big bang, the universe's evolution of the second type of states is well described by the eigenstates of the Hamiltonian (58) described in section 2.1. Such asymmetry of the interaction on different ǫ sectors is not an artifact of the simplicity of the interaction Hamiltonian. This is just a consequence of the necessary breaking of the shift invariance for any realistic matter interaction as we argued in the previous section and we show explicitly in Appendix D (see Figure 3). Therefore, the non trivial scattering problem concerns only states on the lattice Ŵ 1 = {ν = 4n √ ℓ p | n ∈ Z} that is preserved by the Hamiltonian and contains the point ν = 0. In order to solve the scattering problem, we consider an in-state of the form where ν ∈ Ŵ 1 , and A(k) and B(k) are coefficients depending on k. For suitable coefficients, such states are eigenstates of the full Hamiltonian (35). Arbitrary solutions (wave packets) can then be constructed in terms of appropriate superpositions of these "plane-wave" states. We can compute the scattering coefficients A(k) and B(k) from the discrete (finite difference) time-independent Schrodinger equation which amounts to the following finite difference equation in the ν basis: Frontiers in Astronomy and Space Sciences | www.frontiersin.org The matching conditions on ν = 0 are given by: where the first equation comes from continuity at ν = 0, the second equation from the time-independent Schrodinger equation. The solution of the previous equations is We consider an in-state of the form (valid for early times) where ψ(b; b 0 , ν 0 ) is a wave function picked at a b = b 0 value and ν = ν 0 . Notice that we are superimposing two wave packets supported on lattices Ŵ 1 and Ŵ 2 , respectively. We can now write the pure in-density matrix which scatters into the out-density matrix Let us assume that ψ(b) is highly picked at a b 0 so that we can substitute the integration variables b and b ′ by b 0 and have a finite dimensional representation of the reduced density matrix after the scattering (this step is rather formal, it involves an approximation but it helps when visualizing the result). In the relevant 4 × 4 sector (with basis elements ordered as |1, b 0 , |1, −b 0 , |2, b 0 , |2, −b 0 ), we get the matrix representation

Matter Coupling Produces a Coarse-Graining Entropy Jump at the Big Bang
A reduced density matrix encoding the notion of coarse graining associated with the low-energy equivalence of the ǫ sectors is defined by tracing over the discrete degrees of freedom labeling the component of the state in either the Ŵ 1 or the Ŵ 2 lattices. In other words, tracing over the two (macroscopically indistinguishable) ǫ sectors, namely In other words, the subspace of the Hilbert space we are working with is the one supported on two different ǫ sectors H(Ŵ 1 ) ⊕ H(Ŵ 2 ) ⊂ H lqc which, as the two terms are isomorphic H(Ŵ 1 ) ≈ H(Ŵ 2 ) ≈ H 0 , H(Ŵ 1 ) ⊕ H(Ŵ 2 ) ⊂ H lqc can be written as The coarse graining is defined by tracing over the C 2 factor. This implies that from the previous 4 × 4 matrix, we obtain 2 × 2 reduced density matrices. The reduced density matrix ρ R in remains pure, explicitly Nevertheless, the reduced density matrix ρ R out is now mixed, namely We can now compute the entanglement entropy. To first order the cosmological constant, the result is The behavior as a function of b is shown in Figure 4. In Appendix C1, we discuss an alternative definition of coarse graining with the same qualitative implications.  (93) is apparent. The entropy is periodic for bℓ p ∈ [0, π] as expected from (86). The dotted line represents the maximum possible entropy which is log[2] in our model.

QUANTUM COSMOLOGY ON A SUPERPOSITION OF BACKGROUNDS
In the first part of this paper, we have seen how the fact that the Hilbert space of loop quantum cosmology is vastly larger than the standard Schrodinger representation implies (via coarse graining) that the coarse-graining entropy would rise generically through the evolution across the big bang would-besingularity. In this section, we explore another closely related feature that leads to an apparent non-unitary evolution when dynamics are probed by a low-energy agent. The quantum dynamics in loop quantum cosmology depend on a UV regulator known as the area-gap (section 2.1). We will see here that the loop quantum cosmology model can be extended naturally to admit superpositions of dynamics with different regulators. Such an extension generally leads to the dynamical development of correlations between the macroscopic and the microscopic degrees of freedom. If the microscopic degrees of freedom are assumed to remain hidden to low-energy observers, then such correlations lead to an apparent violation of unitarity in the low-energy description where pure states evolve into mixed states. In other words, the UV data needed to define the quantum dynamics open an independent microscopic channel for information to be degraded.

The UV Input in Quantum Cosmology: Revisiting the µ Scheme
Theμ scheme was designed to avoid an inconsistency in an early model of loop quantum cosmology with the low-energy limit (or large universe limit) of loop quantum cosmology (Green and Unruh, 2004). The problem arises from the effective compactification of the connection variable c due to the polymer regularization of the Hamiltonian with a fixed fiducial scale µ which implies that c and c + 4π/µ are dynamically identified. This leads to anomalous deviations from classical behavior in situations where the variable c is classically expected to be unbounded for large universes. This can be seen clearly in the present situation where the unimodular Hamiltonian (17) is given, in (c, p) variables, by For non-vanishing energies (or equivalently non-vanishing cosmological constants), the conservation of the Hamiltonian implies that c grows as |p| ∝ a 2 , i.e., c grows without limits as the universe expands so that no matter how small µ is, anomalous effects due to the compactification of the c become relevant at macroscopic scales (Noui et al., 2005). As no quantum gravity effects seem acceptable in the large universe regime for a model with finitely many degrees of freedom, this anomaly is seen as an inconsistency of the model. Theμ scheme solves this inconsistency by "renormalizing" the regulating scale µ as the universe grows (Equation 20). The interesting thing is that such a renormalization is justified by quantum geometry arguments that link the mini superspace model of loop quantum cosmology to the geometry of a microscopic background state in the full theory. The argument explicitly uses the idea that the low-energy degrees of freedom (dynamical variable of loop quantum cosmology) arise from the coarse graining of the fundamental ones in loop quantum gravity.
Here we review the construction of theμ scheme as described in Ashtekar and Singh (2011). Consider a fundamental quantum geometry state |s in the Hilbert space of loop quantum gravity, representing a microscopic state on top of which the quantum cosmological coarse-grained dynamics will eventually be defined. Such an underlying fundamental state will have to be approximately homogeneous and isotropic up to a scale of L > ℓ p with respect to the preferred foliation defining the comoving FLRW observers at low energies. If that is the case then such space slices can be divided into (approximately) cubic 3-cells of physical side length L which all have approximately equivalent quantum geometries. The area of a face of such cubic cells in Planck units will be denoted s so that L 2 = ℓ 2 p s . Note that s is a property of the underlying microstate: an area eigenstate if the microstate is an eigenstate, or an area expectation value if the state is sufficiently peaked on a quantum geometry and has small fluctuations around it. A simple realization is the one where s is an area eigenvalue, and the important assumption is that s is the same for all cells (this encodes the homogeneity of the microscopic state). Consider the area of a large twodimensional surface (the face of a fiducial cell V), whose area is measured by the low-energy (coarse-grained) quantity p used as the configuration variable in loop quantum cosmology. We naturally would expect that |p| ≫ ℓ 2 p or alternatively that where N denotes the number of microscopic cells contained in the coarse-grained surface (a face of V), and N ≫ 1. The fiducial cell has a fiducial coordinate volume V 0 and hence fiducial side coordinate length V 1/3 0 . Therefore, the fiducial coordinate length i.e., the fiducial scaleμ is dynamical: as the universe grows (and |p| becomes larger), the underlying fiducial length scale decreases. The fiducial regularization scale (53) depends on the fundamental state |s via the quantity s , hence we denote it asμ s . When such a dynamical scale is used in the regularization of the quantum cosmology Hamiltonian, the effective compactification scale for c grows like |p| and the inconsistency previously discussed is avoided. This is transparent in terms of the new canonical pair (b, v). From Equation (15), we can see that b = cμ s /( √ s ℓ p ), in contrast with c (see 17), remains constant (see 51) in the de Sitter universe. The quantization of the Hamiltonian presented in section 2.1 introduces an effective compactification of the variable b whose dynamical effect is now only relevant when the cosmological constant approaches one in Planck units. This can be seen from (21). The cosmological constant is bounded from above by its natural value in Planck units due to the underlying quantum geometry structure while the anomalous IR behavior is avoided (the problems exhibited in the model studied in Green and Unruh, 2004 are also resolved).
The above is the standard account of the motivation of thē µ scheme of Ashtekar et al. (2006) with the little twist (which is very important for us here) that s need not be the lowest area eigenvalue of loop quantum gravity. In the usual argument, the microscopic state is thought to be built from a special homogeneous spin network (geometry eigenstate) with all spins equal to the fundamental representation. This implies that, in the above construction, s = 1/2 ≡ 2πγ √ 3. The observation here is that s can take different values according to the microscopic properties of the underlying quantum geometry state. One could take for instance all spins equal to the vector representation and then have s = 1 ≡ 4πγ √ 2 instead, or take j as arbitrary and use s = j . It is important to point out that such a possibility can arise naturally in quantum cosmology models obtained in the group field theory framework (Gielen et al., 2013;Oriti et al., 2016Oriti et al., , 2017. As we have seen in section 3.1, the field strength regularization, and hence the Hamiltonian, depend on the value s of the background (approximately homogeneous) spin network state |s through the dynamical scaleμ s . In this way, the dynamics of loop quantum cosmology establish correlations with a microscopic degree of freedom in the underlying loop quantum gravity fundamental state. As such a degree of freedom (the area eigenvalue s of the minimal homogeneity cells) is quantum, it is natural to model the system by a tensor product Hilbert space H ≡ H m ⊗ H lqc , where H m is the Hilbert space representing the microscopic degree of freedom encoded in the minimal homogeneous cell operator (whose eigenvalues we denote as s ), and H lqc the standard kinematical Hilbert space of loop quantum cosmology.
General states in H can be expressed as linear combinations of product states |s ⊗ ψ in the respective factor Hilbert spaces. The quantum Hamiltonian has a natural definition on such states and therefore on the whole of H, namelŷ whereĤ s is the usual loop quantum cosmology Hamiltonian in theμ s scheme, which in our particular case is defined in Equation (21) with regulator = s . Notice that the previous extension of the standard loop quantum cosmology framework to the larger Hilbert space H is also natural from the perspective of the full theory. Indeed, the generally accepted regularization procedure of the Hamiltonian constraint in loop quantum gravity (first introduced by Thiemann, 1998 and further developed in recent analysissee Ashtekar and Pullin, 2017 and references therein) is state dependent in that the loops defining the regulated curvature of the connection are added on specific nodes of the state where the Hamiltonian is acting upon. This feature finds its analog in the action (55) where the regulating scale s depends on the state |s ∈ H m .
In order to simplify the following discussion we will restrict states in H m even further and consider a subspace h = C 2 ⊂ H m , i.e., we will model the situation where the underlying microscopic state is an arbitrary superposition of only two fixed microscopic homogeneous spin-network states. For example we take where |± ∈ H m are two suitable orthogonal background states (these two states will be conveniently picked below). From the infinitely dimensional Hilbert space H m we are now selecting a single q-bit subspace C 2 . The Hilbert space of our model is The factor h represents additional microscopic (hidden to lowenergy observers) UV degrees of freedom, while H lqc encodes the data that under suitable circumstances (e.g., when the universe is large) represent the low-energy cosmological degrees of freedom.
In this way we see that in addition to the intrinsic degeneracy of energy eigenvalues analyzed in the first part of this paper, there is another candidate for the microscopic degree of freedom associated with the regularization of the Hamiltonian action via theμ-scheme. Both mechanisms are correct for the present loop quantum cosmology toy model but reflect generic properties of the full theory of loop quantum gravity. More generally, we expect similar features to be present in any quantum gravity approach where smooth geometry is only emergent from a discrete fundamental theory. From now on we adopt the convenient notation |s with s = ± for such preferred basis elements of h. With this notation, and using (21), the Hamiltonian (55) becomeŝ where (ν) ≡ ν|ψ . Now, the only special feature of the basis |± is that it is preferred from the perspective of the regularization of the effective (unimodular) loop quantum cosmology Hamiltonian. Consequently, a natural question for the quantum theory is how the dynamics would look if the initial state is arbitrary in factor h? More precisely, what if we consider the linear combination of two background spin networks 1 √ 2 (|+ + |− ) ∈ h times a loop quantum cosmology wave function as depicted in Figure 5?
To answer this question we consider a special initial state where correlations between the low-energy and the UV degrees of freedom are not present. This will lead to a reduced density matrix-tracing out the microscopic space h in (57)-that is pure initially, the form of such a state is illustrated in Figure 5. At the same time we want to be able to diagonalize the Hamiltonian with such uncorrelated initial states; more precisely this boils down to diagonalizing both H + and H − in H lqc . This implies that the factor ψ(ν) ∈ H lqc , in Figure 5, must be supported on a lattice Ŵ ǫ that is left invariant by the action of both H + and H − (left invariant in the sense that the shift operators in the definition of the Hamiltonian only relate to points of Ŵ ǫ and never map points out). This can be achieved by assuming that √ + = m √ − for some natural number m. For simplicity we will take m = 2 from now on 11 . The parameter ǫ will be taken so that the lattice Ŵ ǫ contains the point ν = 0. This is a standard choice. With all this, the invariant lattice, denoted as Ŵ − , is 11 One might be worried that this is hard to achieve if one sticks to the form of the area spectrum of loop quantum gravity. This is however simply a model and the link with the full theory (remember) must be taken at the heuristic level. Nevertheless, solutions do exist for instance m = 4 for j + = 3 and j − = 1/2.
Note that in the notation described below (37), we have that The choices made above are not mandatory. One could have chosen a different initial state. The previous choice is particularly interesting here because it would lead to a reduced initial density matrix that is pure and hence initially had vanishing entanglement entropy. Other states would involve correlations and would therefore carry a non vanishing entropy load from the beginning. For the discussion that interests us here and for the analogy with black hole evaporation, it is more transparent to set the entropy to zero initially.
An arbitrary (unimodular) loop quantum cosmology state associated with such a choice of background state can be expressed as: where ψ(k; b 0 , ν 0 ) is a properly normalized function peaked at k = b 0 and ν = ν 0 . The initial state in the momentum representation is given by: Where, in the first line, we used the natural extension of the notation introduced in (37) where |b, 1 ∪ 2 means an eigenstate of the corresponding shift operators (19) supported on the lattice Ŵ − = Ŵ 1 ∪ Ŵ 2 . Notice that we can also write |b, 1 ∪ 2 = |b, 1 + |b, 2 , keeping in mind that terms on the r.h.s. are individually eigenstates of the shift operators with twice the lattice spacing of Ŵ − . We also used We can then write where is the Haar measure on the circle of circumference π/ √ 2 − ℓ p . We notice from (64) that even when our initial state contains no correlations between the low energy degrees of freedom represented by b and the microscopic degrees of freedom encoded in |s at t = 0, quantum correlations between the two will develop with time due to the non trivial dependence of the energy spectrum with s. Even when this is quite clear from (64), one can state this fact in an equivalent way by analysing the (pure) density matrix ρ in (t) ≡ | in (t) in (t)|, whose matrix elements in the b basis are: As coarse-grained observers are assumed to be insensitive to the microscopic structure encoded here in the "spin" quantum number s, low-energy physical information is encoded in the reduced density matrix which can be simply be written as Notice that (69) is only pure at t = 0 and becomes mixed due to the correlations evoked above as time passes. One can compute the entanglement entropy S(t) ≡ −Tr ρ inR (t) log(ρ inR (t)) which turns out to be given by the simple analytic expression (see Appendix A) where For generic wave packets ψ s (b), the entanglement entropy is a monotonic growing function of time which grows asymptotically to the maximally mixed situation S max = log(2) (see an example in Figure 6).
FIGURE 6 | Here we plot S(t) as a function of time for a Gaussian wave packet centered at b = 2.5 10 −2 , b = 5 10 −2 , b = 7 10 −2 , and b = 10 −1 with width σ = b, respectively. Numerical integration plus the approximation (73) was used with the assumption that 2( + − − )/γ 2 = 1, all in Planck units. As b grows the scalar curvature (the cosmological constant) grows and the rate at which entropy increases grows as well. For b ≪ 1, an effective unitary evolution is recovered.
A more intuitive picture can be obtained from a suitable expansion of the energy eigenvalues (24) in powers of the label bℓ p Such an expansion makes sense in that it allows for the identification of the low energy effective Hamiltonian (the one that one would define in a purely Wheeler-DeWitt quantization) plus corrections that involve interactions with the underlying discrete structure of LQG here represented by the spin s degree of freedom. Namely, we can read from the previous expansion that whereĤ 0 eff (b) ≡ 6 γ 2b 2 is the Wheeler-DeWitt Hamiltonian and the additional term is an interaction with the environment represented by the underlying discrete structure represented by the dependence on s (a hidden degree of freedom from the low-energy continuum perspective). Of course, the hats in the previous equation denote operators in a different representation (the continuum Schrodinger representation) that is not unitarily equivalent to the "fundamental" polymer representation introduced in section 2.1 and used in the LQC setup (recall for instance that the operatorb does not even exist in the polymer representation). The lack of purity for t > 0 of the reduced density matrix (69) is due to correlations that develop between the low-energy degree of freedom b and the hidden microscopic degree of freedom s via this non trivial interaction Hamiltonian. This means that generically [i.e., for arbitrary initial states ψ s (b)], the fundamental evolution would seem to violate unitarity, from the perspective of low-energy observers, due to the decoherence with the microscopic quantum geometric structure. Notice however that for states ψ s (b) picked at sufficiently smallk, i.e.,k √ s ℓ p ≪1, we have from (22) and the density matrix (69) is pure at all times. More precisely, we can translate the criterion for the absence of decoherence with the underlying microscopic discrete structure in terms of the value of the cosmological constant of the given state. For an eigenstate of the Hamiltonian the relation is given by ≡ E s (b). Therefore the criterion for the absence of decoherence in terms of the cosmological constant is Interestingly, for states with low values of the cosmological constant in natural units-equivalent semi-classically to the scalar curvature R in our matter-free model-define a decoherence free subspace. When the cosmological constant does not satisfy condition (76), decoherence with the microscopic structure is turned on and maximized for of order one in Planck units: notice incidentally that due to the polymer quantization the cosmological constant is bounded by For low values of , unitarity is recovered in the effective description that ignores the microscopic structure. Decoherence takes place here due to an interaction between the low-energy coarse degrees of freedom and the microscopic discreteness in the underlying quantum geometry background, but in a way (in our simple model) that the energy and hence the cosmological constant is conserved. However, the presence of decoherence suggests the possibility for a natural deviation of this idealized absence of dissipation: generally decoherence and dissipation often come together. Therefore, a surprising and unexpected consequence of our analysis is the suggestion of a natural channel for the relaxation of a large cosmological constant due to the possibility of dissipative effects associated with the decoherence pointed out here.
Incidentally, all this shows that only in the limit of low values of E (small cosmological constant), the coarse graining that leads from the full theory of loop quantum gravity to the mini superspace description of loop quantum cosmology is well-defined. This is not surprising and only confirms the usual intuition that drives the construction of models of loop quantum cosmology. However, it opens the door for a qualitative understanding of the necessity of decoherence effects in more general situations. For instance, the standardμ s construction suggests that coarse graining is weaker at the big bang where the Hamiltonian evolution (58) takes the universe through the ν = 0 states. During this high (spacial) curvature phase it is natural to expect that the higher corrections in (73) (describing the interaction with the microscopic Planckian structure) can no longer be neglected.
Interestingly, there is another way to make decoherence disappear. This is due to the asymptotic behavior of the separation of area eigenvalues in loop quantum gravity which imply that for large s there are states such that s − s ′ ≈ s exp(−π √ 2 s /3) (Fernando Barbero et al., 2018). Therefore, in the continuum limit s − s ′ ≪ 1 the dynamical entanglement growth of our model can be made as small as wanted.

Matter Coupling Produces an Entanglement Entropy Jump at the Big Bang
In the pure gravity case, we can make decoherence as small as needed by choosing states with a cosmological constant that is sufficiently small. Here we show that this is no longer possible once matter is added and that there is a generic development of correlations with the UV degrees of freedom in the evolution across the would-be-singularity: an initially pure state (reduced low-energy density matrix) evolves generically into a mixed state (reduced low-energy density matrix) after the big bang.
In order to see this in more detail, we just need to write out the matter Hamiltonians acting in the Hilbert space (57). One needs the natural generalization of the expressions written in section 2.3 for the present context. For instance, for scalar field coupling, Equation (32) becomeŝ where The momentum p φ commutes with the Hamiltonian and thus is a constant of motion. As before, if we consider an eigenstate of p φ then the problem reduces again to a scattering problem with a potential decaying like 1/ν 2 when solving the time-independent Schrodinger equation From the discussion in section 2.3, we can capture the basic qualitative effect of matter interaction by considering a simple solvable model where the matter contribution is concentrated at a single event at the big bang. None of the qualitative conclusions that follow depend on this simplification, and a more realistic free scalar field model can be dealt with (some results are shown in Appendix D). With some extra effort one could actually analyze a more realistic model [say the one defined by (78)] but the conclusion will remain the same. Therefore, we consider where µ is a dimensionless coupling,Ĥ 0 is given in (58), andĤ int is the generalization of (36) whereÔ is a self adjoint operator in h = C 2 . A natural and simple model for this operator is to choosê This choice is formulated in the notation introduced below (58) and inspired by the analogy with a spin system. We have added by hand an interaction Hamiltonian that switches on only when the universe evolves through the would-be-singularity at the zero volume state. This encodes the idea of the intrinsic uncertainty of the peculiar construction of the mini superspace model of loop quantum cosmology that we discussed in section 3.1. The discrete local degrees of freedom must be important close to the big bang and symmetry reduction must fail in some way that can only be correctly described if a full quantum gravity theory is available.
Here we model such unknown dynamics in the simplest fashion available to us here, which consists of including the possibility for the background state |s (representing in spirit the underlying quantum geometry) to be modified by the dynamics viaĤ int .
Here we proceed as in section 2.4 while keeping in mind that, in the present case, there are two distinct cases at hand given by the two possible values ± . Let us consider an in-state of the form where A s (k) and B s (k) are coefficients depending on k and (in contrast with the case in section 2.4) now also on s = ±1 (with |± the eigenstates ofŜ z ). For suitable coefficients, such states are eigenstates of the Hamiltonian H 0 as well as the full Hamiltonian (35). Arbitrary solutions (wave packets) can then be constructed in terms of appropriate superpositions of these "plane-wave" states. where One can superimpose the previous eigenstates to produce wave packets (semiclassical states) for the wave function of the universe that are picked at value ν 0 of the rescaled volume (see footnote 6). Wave packets will evolve in time according to the Schrodinger equation which in our case is just a discrete analog of the one corresponding to a free particle in quantum mechanics with an interaction term at the "origin" ν = 0. If we start with a state that is sufficiently picked around ν 0 for ν ≫ ℓ p initially, then the state can be described in terms of the superposition (64) where the explicit values of the coefficients A s (b) and B s (b) do not appear. Equation (42) is generalized to The coefficients (85) enter the expression of the scattered wave packet at a later time which becomes Note that the solution of the scattering problem for the E + (b) eigenvalues is asymmetric with respect to the components of the in-state supported on Ŵ 1 and Ŵ 2 . Indeed the states |b, 2 are eigenstates of the Hamiltonian directly because they are not supported on ν = 0 and hence they do not "see" the interaction: this is captured by trivial scattering coefficients for this component.

Entropy Associated With the Entanglement With the UV Degrees of Freedom
From the previous initial state we can calculate [by tracing over the factor h, see (57)] the initial reduced density matrix The reduced density matrix after the big bang is where α + = 1/4 and α − = 1 and δ s+ is unity when s = + and vanishes when s = −. Then the non vanishing entries of the reduced density matrix are = ρ R 11 out (−b 0 , b 0 ) The matrix ρ R out is positive definite, Tr[ρ R out ] = 1 and ρ R out = ρ R † out . In the case b 0 ℓ p ≪ 1 we have We can now compute the entanglement entropy jump δS to the first leading order in b 0 ℓ p /µ. The result [expressed in terms of the cosmological constant in this regime, namely (75)] is δS = δ 0 S − 3 − ℓ 2 p log(3) 128π 2 γ 2 µ 2 + O( 2 ℓ 4 p ),  (93) is apparent. The entropy is periodic for bℓ p ∈ [0, π] as expected from (86).
where δ 0 S = 2 log(2) − 3 4 log(3). The previous equation shows that the entropy jump is non trivial at crossing the big bang would-be-singularity, even in the low cosmological (low-energy) limit where (according to the analysis of the previous section) decoherence with the microscopic Planckian structure can be neglected during the time the universe is large. Information is unavoidably degraded (it seems lost for low-energy observers) during the singularity crossing.
The general entropy jump for arbitrary (not necessarily small) can be computed explicitly. Its value is bounded by log(2) in our model. Finally, the energy is conserved through the big bang and during all the dynamical evolutions for the arbitrary values of b 0 . The decoherence and entanglement which can be interpreted as an information loss happens without energy spending as required by the scheme put forward in Perez (2015).

DISCUSSION
We have seen that one can precisely realize the scenario put forward in Perez (2015) for the resolution of Hawking's information loss paradox in quantum gravity in the context of loop quantum cosmology. The key feature making this possible is the existence of additional degrees of freedom with no macroscopic interpretation which unavoidably entangle with the macroscopic degrees of freedom during the dynamical evolution and lead to a reduced density matrix whose entropy grows. The fundamental description is unitary but the effective descriptionthat does not take the microscopic degrees of freedom into account and hence is analogous to the QFT description of BH evaporation-evolves pure states into mixed states. The microscopic degrees of freedom in the toy model are not introduced by hand, their existence is intimately related to the peculiar choice of representation of the fundamental phase space variables that leads to singularity resolution (Bojowald, 2001). Moreover, such a representation mimics the one used in the full theory of loop quantum gravity (Lewandowski et al., 2006) where also one expects such extra residual and microscopic degrees of freedom to exist and remain hidden to low-energy coarse-grained observers describing physics in terms of an effective QFT.
From a more general perspective, we expect this scenario to transcend the framework of loop quantum gravity: in any approach to quantum gravity, where spacetime geometry is emergent 12 from more fundamental discrete degrees of freedom, the effect (precisely illustrated here by our toy model) would generically occur.
These results extrapolated to the context of black hole formation and evaporation suggest a simple resolution of the information paradox that avoids the pathological features of other proposals. For instance, the possible development of firewalls (Almheiri et al., 2013;Braunstein et al., 2013) or the risks of information cloning that the holographic type of scenarios must deal with (Marolf, 2017) are completely absent here. As decoherence in our model takes place without diffusion Unruh (2012), the usual difficulties (Banks et al., 1984) with energy conservation in the purification process are avoided along the lines of Unruh and Wald (1995), Unruh (2012) in a concrete quantum gravity framework (hence without the problems faced by the QFT approach Hotta et al., 2015;Wald, 2019).
We notice that the possibility of decoherence illustrated in the present model also suggests the possibility of diffusion into the underlying Planckian structure, such diffusion might have, in suitable situations, important consequences at large scales as argued in a series of recent papers (Josset et al., 2017;Perez et al., 2018;Perez and Sudarsky, 2019). The present model is very simplistic and realizes an example where such diffusion is not possible due to (unimodular) energy conservation and the fact that the microscopic degrees of freedom do not contribute independently to the Hamiltonian. Nevertheless, one could generalize these models easily in order to include diffusion. 12 For instance in the causal sets approach (Bombelli et al., 1987), or in the context of Jacobson's ideas about emergence (Jacobson, 1995) (where, incidentally, in both cases unimodular gravity is the natural emergent structure), in causal dynamical triangulations (Ambjorn et al., 2004), and group field theory (Oriti, 2011), etc. This possibility is under current investigation and we plan to report the results elsewhere.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

AUTHOR'S NOTE
This is a toy model of quantum cosmology that illustrates a novel mechanism for the resolution of Hawking's information puzzle for black hole evaporation.