arxiv: v1 [gr-qc] 3 Jan 2008

Transcription

1 Exact solutions for Big Bounce in loop quantum cosmology Jakub Mielczarek Astronomical Observatory, Jagiellonian University, - Kraków, Orla 7, Poland and The Niels Bohr Institute, Copenhagen University, Blegdamsvej 7, DK- Copenhagen, Denmark Tomasz Stachowiak Astronomical Observatory, Jagiellonian University, - Kraków, Orla 7, Poland Marek Szyd lowski Department of Theoretical Physics, Catholic University of Lublin, Al. Rac lawickie, -95 Lublin, Poland and Marc Kac Complex Systems Research Centre, Jagiellonian University, Reymonta, -59 Kraków, Poland arxiv:8.5v gr-qc] Jan 8 In this paper we study the cosmological FRW model k = with holonomy corrections of Loop Quantum Gravity. The considered universe contains a massless scalar field and cosmological constant Λ. We find analytical solutions for this model in different configurations and investigate its dynamical behaviour in the whole phase space. We show the explicit influence of Λ on the qualitative and quantitative character of solutions. For the case of positive Λ the oscillating solutions without the initial and final singularity appear as a generic case for some quantisation schemes. I. INTRODUCTION Loop Quantum Gravity (LQG) has drastically changed our classical view of the early Universe. Firstly, the serious problem of initial singularity has been solved by quantum effects on the Planck scale which lead to bouncing solutions. Secondly, the new phase of super-inflation was introduced. In this paper we concentrate on the former effect. We consider the flat FRW model with a free scalar field and with the cosmological constant. Loop Quantum Gravity effects are introduced in terms of quantum correction to classical theory. Generally one consider two types of of quantum correction: correction from inverse volume and holonomy corrections. The leading effect of the volume corrections is the appearance of the super-inflationary phase. The effect of holonomy corrections is the appearance of a bounce resulting in singularity avoidance. The aim of this paper is to investigate analytically these effects in a flat FRW model. It means that we neglect corrections from inverse volume, these effects however has been extensively studied elsewhere. On the other hand, both types of corrections are not equally important in the same regimes. The inverse volume corrections are mainly important for small values of the scale factor, while holonomy corrections are mainly important for large values of the Hubble parameter. In other words, when the minimal scale factor (during the bounce) is large enough, the effects of inverse volume corrections can be neglected. The flat FRW model in Loop Quantum Gravity has been derived in the papers, ] and then improved and studied in ]. This is the simplest cosmological model and has been studied as the firs one in LQG. Later on, other cosmological models were studied, like the FRW k ± models in, 5, 6] and Bianchi I in 7, 8]. We must mention that cosmological models in LQG were considered earlier in the pioneer works of Bojowald. However holonomy correction, in the form which we consider here, were derived later in the cited works of Ashtekar,Paw lowski and Singh. In these models with holonomy corrections, there occurs an unambiguity in the choice of the elementary area. In this paper we consider two kind of approaches to this problem called µ scheme and µ scheme (for a more detailed description see Appendix A). We find analytical solutions for considered models in these two schemes. The Hamiltonian of the considered model is given by H eff = ] sin ( µc) p 8πGγ + p φ µ p + Λ / p / 8πG. () In Appendix A we show the derivation of this Hamiltonian from Loop Quantum Gravity. The canonical variables for gravitational field are (c, p) and for the scalar field (φ, p φ ). The equations of motions can now be derived with the use of the Hamilton equation f = {f, H eff } () where the Poisson bracket is defined as follows {f, g} = 8πGγ f g c p f ] g p c f g + f ] g. () φ p φ p φ φ Electronic address: jakubm@poczta.onet.pl Electronic address: toms@oa.uj.edu.pl Electronic address: uoszydlo@cyf-kr.edu.pl From this definition we can immediately retrieve the elementary brackets {c, p} = 8πGγ and {φ, p φ } =. ()

2 With use of the Hamiltonian () and equation () we obtain equations of motion for the canonical variables ṗ = p γ µ sin( µc)cos( µc), { ċ = ] } p sin ( µc) γ p µ sgn(p) κγ φ = p / p φ, p φ + sgn(p)λγ 5/ p, p p φ = (5) where κ = 8πG. The Hamiltonian constraint H eff = implies ] sin ( µc) γ = κ p φ p µ p + Λ. (6) The variable µ corresponds to the length of the edge of the elementary area and can be written in the general form µ = ξ p n (7) where / n and ξ is a constant ξ >. The choice of n and ξ depends on the particular scheme in the holonomy corrections. Description of these ambiguities can be found in Appendix A. Combining equations (6), (7) and the first one from the set of equations (5) we obtain ( ) dp = p + I p II p n dt V p n Ω V p n+ (8) where new parameters are defined as follow = κp φ, (9) I = Λ, () II = 9 κ γ ξ p φ, () V = 9 κλp φ γ ξ, () Ω V = 9 Λ γ ξ. () The equation (8) is, in fact, a modified Friedmann equation H = 8πG eff ρ + Λ eff () where effective constants are expressed as follow ] G eff = G ρρc, (5) Λ eff = Λ ρ Λ ] (6) ρ c κρ c and p φ ρ = p, (7) ρ c = κγ µ p. (8) In this paper we study solutions of the equation (8) for both µ scheme and µ scheme. The organisation of the text is the following. In section II we consider models with Λ =. We find solutions of the equations (8) for both µ scheme and µ scheme. Next, in section III we add to our considerations non vanishing cosmological constant Λ. We carry out an analysis similar to the case without lambda. We find analytical solutions of the equation(8) for µ scheme. Then, we study the behaviour of this case in µ scheme. In section V we summarise the results. II. MODELS WITH Λ = In this section we begin our considerations with the model without Λ. Then, equation (8) simplifies to the form ( ) dp = p II p n. (9) dt We solve this equation for both µ scheme and µ scheme in the present section. A. µ scheme (n = /) In the µ scheme, as it is explained in the Appendix A, the µ is expressed as µ = () p where πγl Pl is the area gap. So ξ = and n = /. To solve (9) in the considered scheme we introduce a new dependent variables u in the form p(t) = u / (t). () With use of the variable u the equation (9) takes the form ( ) du = 9 u(t) 9II. () dt and has a solution in the form of second order polynomial u(t) = II + 9 t 8 C t + 9 C ()

3 where C is a constant of integration. We can choose now C =, so that the minimum of u(t) occurs for t =. Going back to the canonical variable p we obtain a bouncing solution ΩIII p(t) = sgn(p) + 9 ] / t. () The main property of this solutions is that p(t) never reaches zero value for non vanishing p φ. The minimal value of p(t) is given by p min = γl Pl 8 π ( pφ l Pl ) ] /. (5) We show the solution () in the Fig.. The dynamical 5 Applying this new parametrisation to the equation (9) leads to the equation in the form which has a solution ( ) dy = y(u) II (9) du y(u) = II + u C u + C. () We can now choose normalisation taking C =, then the minimal value of y(u) is for u =. Now we can go back to the initial parameters p and t, then ΩIII p(u) = sgn(p) + u Ω () I u ] / t(u) = du ΩIII + u () Introducing the variable p x = ΩIII u () we can rewrite integral () to the simplest form t FIG. : Typical solution with Λ = in the µ scheme. behaviour in this model is simple. For negative times we have contracting pre Big Bang Universe. For t > we have a Big Bang evolution from minimal p min (8). It is, however, convenient to call this kind of stage Big Bounce rather than Big Bang because of initial singularity avoidance. B. µ Scheme (n = ) In the µ Scheme the µ is expressed as µ = µ =. (6) So ξ = µ and n =. To solve (9) in this scheme, we change the time dt = p(u) du (7) and introduce a new variable y as follow y(u) = p(u). (8) t(u) = ( ΩIII ΩIII ) / ΩIII u and the solution of such integral is given by x + y ] / dy = x( + x ) / dx + x ] / () + x F,, ; x ] where F is the hypergeometric function, defined as pf q a,..., a p ; b,..., b q ; x] = k= (a ) k...(a p ) k x k (b ) k...(b q ) k k! (5) (6) where (a) k is the Pochhammer symbol defined as follow (a) k Γ(a + k) Γ(a) = a(a + )...(a + k ). (7) This solution is very similar to the on in µ scheme. However, the time parametrisation is expressed in a more complex way. We show this solution in the Fig.. In this case minimal value of p(t) is expressed as p min = π γl Pl ( pφ l Pl ). (8)

4 5 and introduce a new variable p as follow y(u) = p(u). (6) This change of variables leads to the equation in the form p ( ) dy = 9γ I y (y y )(y y ) (7) du t FIG. : Typical solution with Λ = in the µ scheme. III. MODELS WITH Λ where y = κp φ γ ξ 6 Λ ] (8) y = κp φ Λ γ ξ (9) We have three general types of solutions corresponding to the values of the parameters (γ I, y, y ). We summarise these possibilities in the table below. In this section, we investigate the general model with non vanishing cosmological constant. It will be useful to write equation (8) in the form ( ) dp = p + I p ] dt II p n Ω ] V p n+ (9) I We see that when we perform multiplication in this equation and define V = III + Ω V I () we recover equation (8). In this and the next section we us equivalently Λ = ξ γ and α = α 6 p φ to simplify notation. A. µ scheme (n = /) In this case equation (9) can be rewritten to the form ( ) dp = γ I p + γ II p γ III p () dt where the parameters are defined as γ I = I Ω V = Λ Λ ] γ ξ, () γ II = V = κp φ Λ ] γ ξ, () γ III = II = 9 κ γ ξ p φ. () To solve equation () we re-parametrise the time variable dt = p(u) du (5) Λ < > > Λ γ ξ > > < γ I < > < y > > < y > < < It is important to notice that the product Υ γ I y y = α Λ is negative in all cases. This property will be useful to solve equation (7). In the Fig. we show values of the roots y and y in the function of the Λ. We y y y y 5 5 Λ FIG. : Parameters y and y as functions of Λ. Region corresponds to oscillatory solution, region to bouncing solution. There is no physical solutions in region. see that there are two values of cosmological constants where the signs of the roots change, namely Λ = and Λ, µ = m Pl πγ =.m Pl. (5)

5 5 Now we perform a change of variables in the equation (7) in the form and w(u) = y(u) We also introduce the parameters (5) χ = Υ u. (5) w = y, w = y (5) and then equation (7) takes the form ( ) dw = (w w)(w w). (5) dχ The equation (5) is the equation of shifted harmonic oscillator and its solution is w(χ) = Λ, µ Λ ] + cos(χ). (55) 6α Λ, µ where we set the integration constant to zero. It is now possible to go back to the initial parameters p and t which are expressed as follows p(u) = sgn(p) t(u) = u 6α Λ, µ Λ Λ, µ + cos ( ) (56) Υ u du p (u ) (57) The integral (57) can be easily solved but solutions depend on the parameters in (56). We define now an integral I(x, a) = x dx a + cosx (58) so the expression (57) can now be written as ( t(u) = Υ I Υ u, Λ ) Λ, µ (59) Solutions of integral (58) depend on the value of parameter a and read ( (a )tan x )] I(x, a > ) = a arctan a ( x I(x, a = ) = tan ) ( (a )tan x )] I(x, < a < ) = arctanh a a ( ) I(x, a = ) = ln cosx tanx In Fig. we show the solution for Λ < or equivalently with parameter a >. In Fig. 5 we show the solution for < Λ < Λ, µ which corresponds to < a < u FIG. : Parametric solution with Λ <. Oscillating curve (green) is p(u) and the increasing curve (red) is t(u)..5 p t FIG. 5: Solution p(t) with Λ < u FIG. 6: Parametric solution with < Λ < Λ, µ. Top curve (green) shows p(u) and bottom curve (red) is t(u).

6 6 p For Λ < equation (6) have bouncing type solutions. This can be seen when we redefine equation (6) to the point particle form. This approach is useful for qualitative analysis and will be used in the next section. However in this section we will use it to select types of solutions. Introducing new time variable dt = p / + I p du, (69) t FIG. 7: Solution p(t) with < Λ < Λ, µ. B. µ Scheme (n = ) In this section we study the last case in which n = and ξ = /. Equation (9) can then be written in the form ( ) dp ( = p + Ω ) II p W( p ) (6) dt where polynomial W( p ) = p a + p a + p a + a (6) and coefficients are expressed as a = 6 κp φγ ξ = α, Λ,µ (6) a =, (6) a =, (6) a = Λ γ ξ = Λ. Λ,µ (65) The determinant of polynomial (6) is defined as = a a a a a a + 8a a a a 7a a a.(66) Inserting the values of parameters {a, a, a, a } listed above we obtain ( ) ( ) ] Λ α α Λ = 7. (67) Λ,µ Λ,µ For >, so for Λ,µ Λ,µ ( ) ( ) α Λ 8 > (68) Λ,µ Λ,µ polynomial (6) have three real roots. When relation (68) is fulfilled and Λ > oscillatory solutions occur. we can rewrite equation (6) to the form H = p + V (p) = E = α (7) Λ,µ where the potential function V (p) = p ( Λ Λ,µ p ). (7) Equation (7) have a form of Hamiltonian constraint for point particle in potential hole. We see that for Λ < potential (7) have only one extrema (for p = ) and only bouncing solution is possible. For Λ < potential (7) have minimum for p min = Λ,µ Λ. In this case physical solutions correspond to condition E > V ( p min ) and the energy of imagined particle in the potential hole is greater than minimum of potential hole. Imagined particle oscillate then between boundaries of the potential hole. The condition E > V ( p min ) is equivalent with relation (68) calculated from determinant. Now we can introducing variables and g := a a a, (7) g := 8 a a a 6 a 6 a a, (7) p = a v a a (7) then equation (6) takes on the form of the Weierstrass equation ( ) dv = v g v g. (75) du Detailed analysis and plots of solutions of this equation for different values of parameters can by find in appendix to the work 9]. Solution of equation (75) is the Weierstrass -function where v(u) = (u u ; g, g ). (76) g =, (77) g = 6 + ( ) Λ α. (78) 6 Λ,µ Λ,µ

7 7 So the parametric solution for the parameter p is p(u) = sgn(p) Λ ],µ Λ (u u ; g, g ) (79) Time variable can be expressed as the integral t(u) = u ( ) / du a v(u ) a a ( ). (8) + I a v(u ) a a In the Fig. 8 and 9 we show an exemplary parametric bounce solution with Λ < and a possible oscillatory solution with Λ >. p p 6 8 u FIG. 8: Parametric solution p(u) with Λ < u FIG. 9: Parametric solution p(u) with Λ > and relation (68) satisfied. In general the solution is expressible as an explicit function of time, by means of the so called Abelian functions. However, we chose not to employ them here as the appropriate formulae are not as clear, and the commonly used numerical packages do not allow for their direct plotting. The fact of existence of such solutions allows us to assume that the above integral is well defined, and so is the solution itself. IV. QUALITATIVE METHODS OF DIFFERENTIAL EQUATIONS IN STUDY EVOLUTIONAL PATHS The main advantage of using qualitative methods of differential equations (dynamical systems methods) is the investigation of all solutions for all admissible initial conditions. We demonstrate that dynamics of the model can be reduced to the form of two dimensional autonomous dynamical system. In our case the phase space is D (p, ṗ). First we can find the solutions corresponding to vanishing of the right hand sides of the system which are called critical points. The information about their stability and character is contained in the linearisation matrix around a given critical point. In the considered case, the dynamical system is of the Newtonian type. For such a system, the characteristic equation which determines the eigenvalues of the linearisation matrix at the critical point is of the form λ + V (p)/ p where V (p) is a potential function and p = a, a being the scale factor. As it is well known for dynamical systems of the Newtonian type, only two types of critical points are admissible. If the diagram of the potential function is upper convex then eigenvalues are real and of opposite signs, and the corresponding critical point is of the saddle type. In the opposite case, if V (p)/ p > then the eigenvalues are purely imaginary and conjugate. The corresponding critical point is of the centre type. ṗ = p n (α + Λ p )( p n+ γ ξ (α + Λ p )) where α = κ 6 p φ p n+ (8) α + Λ ṗ = p = p ( p n γ ξ Λ p p ) γ ξ α (8) where we have made the following time reparametrisation t u d du = p n+/ d α + Λ dt p (8) Now we are able to write the Hamiltonian constraint in the form analogous to the particle of the unit mass moving in the one dimensional potential well H = p + V (p) = E = γ ξ α where the potential function (8) V (p) = p ( p n γ ξ Λ p ). (85) As we can see, the constant γ ξ α/ plays the role of the total energy of the fictitious particle. The domain admissible for motion in the configuration space is determined by the condition V (p) + γ ξ α/ <.

8 8 A dynamical system of the Hamiltonian type has the following form p = y, y = V p (86) where the prime denotes differentiation with respect to a new re-parametrised time variable which is a monotonous function of the original, cosmological time. The structure of the phase plane is organised by the number and location of critical points. In our case, critical points in the finite domain are located only on the line p = y = and the second coordinate is determined form the equation V/ p = which is p ( ( n) p n γ ξ Λ p ) = (87) The number of critical points depends on the value of Λ and n. We can distinguish two cases: for / < n : Λ < : ( p, p ) = (, ); Λ > : ( p, p ) = (, ) and for n = /: ( p, p ) = ( n+ ( n) γ ξ Λ, ); Λ γ ξ : ( p, p ) = (, ); Λ = γ ξ : degenerate V (p) =. The full analysis of the behaviour of trajectories requires investigation also at the infinity. To this aim we introduce radial coordinates on the phase plane for compactification of the plane by adjoining the circle at infinity p = r r cosθ, y = r r sin θ. The phase portraits for both cases are shown at Figs. and. For the general case of Λ > we can write the parametric equation of the boundary α = of the physically admissible region in the phase space, namely ṗ = Λ p ( γ ξ Λ p n+ ) (88) where the dot denotes differentiation with respect to cosmological time t. This equation greatly simplifies for the special case n = /, and we receive the value of the Hubble parameter at the boundary p = H = Λ ( γ ξ Λ ) (89) ṗ In this special case α = and n = we can integrate eq.() for Λ > and choosing the integration constant equal to zero p(t) = sgn(p) Λ,µ Λ cosh ( Λ t ), (9) which gives the maximal value o parameter p p max = Λ Λ. (9) In the same case with Λ < and choosing integration constant equal to zero we have p(t) = sgn(p) Λ,µ Λ sinh ( Λ t ). (9) For α = and Λ < Λ, µ in the case n = / we obtain the de Sitter solution { p(t) = sgn(p) p exp Λ } Λ Λ, µ t. (9) These solutions represent the lines on the boundaries of the physically admissible domains. V. SUMMARY In this paper we have studied dynamics of the flat Friedmann-Robertson-Walker cosmological model with free scalar field and the cosmological constants, modified by the holonomy corrections of Loop Quantum Gravity. We have found analytical solutions for models with both cosmological constant Λ and the free scalar field. We performed calculations in two setups called µ scheme and µ scheme explained in the appendix. We investigate the evolutional paths of the model, from the point of view of qualitative methods of dynamical systems of differential equations. We found that in the special case of n = / boundary trajectory (α = ) approaching the de Sitter state; and demonstrate that in the case of positive cosmological constant there are two types of dynamical behaviours in the finite domain. For the case of n = there appear oscillating solutions without the initial and final singularities, and that they change into bouncing for > n /. The results of this paper can give helpful background dynamics to study variety of physical phenomenas during the bounce in loop quantum cosmology. For example, the interesting question of the fluctuations like gravitational waves,, ] or scalar perturbations ] during this period. The exact solutions for bouncing universe with dust and cosmological constant has been previously studied in ]. The dynamics of the flat FRW model with scalar field has been numerically investigated in 5]. Acknowledgments Authors are greatefull to Orest Hrycyna for discussion. This work was supported in part by the Marie Curie Actions Transfer of Knowledge project COCOS (contract MTKD-CT--5786).

9 9 a) b) FIG. : The phase space diagram for the case n = and a) Λ > and b) Λ <. The physical domain admissible for motion is shaded. Note that for the case (a) the boundary of admissible for motion is bounded by homoclinic orbit and all solutions in this area are oscillating without initial and final singularities. All trajectories situated in physical region posses the minimum value of scale factor. APPENDIX A: FLAT FRW MODEL AND HOLONOMY CORRECTIONS IN LOOP QUANTUM GRAVITY In this appendix, we derive the form of the Hamiltonian () considered in the paper. The FRW k = spacetime metric can be written as ds = N (x)dt + q ab dx a dx b (A) where N(x) is the lapse function and the spatial part of the metric is expressed as q ab = δ ij ω i a ωj b = a (t) o q ab = a (t)δ ij o ω i a o ω j b. (A) In this expression o q ab is fiducial metric and o ωa i are cotriads dual to the triads o e a i, o ω i ( o e j ) = δj i where o ω i = o ωadx i a and o e i = o e a i a. From these triads we construct the Ashtekar variables A i a Γi a + γki a = c o ω i a, (A) E a i detq e a i = p o q o e a i (A) where p = a, c = γȧ. (A5) (A6) Note that the Gaussian constraint implies that p p leads to the same physical results. The factor γ is called Barbero-Immirzi parameter, γ = ln /(π ). In the definition (A) the spin connection is defined as Γ i a = ǫijk e b j ( ae k b] + ec k el a ce l b] ) and the extrinsic curvature is defined as K ab = ] qab D (a N b) N (A7) (A8) what corresponds to K i a := K abe b i. The scalar constraint, in Ashtekar variables, has the

10 a) b) FIG. : The phase space diagram for the case n = / and a) Λ = γ ξ < γ ξ and b) Λ = 5 γ ξ > γ ξ. The physical domain admissible for motion is shaded. In the second case there is no physically allowed region for which α >. This case is distinguished by behaviour at infinity when Hubble function is finite like for the de Sitter solution (see formula (89)). This state is the global attractor in the future. All solutions are of the bouncing type. form H G = d xn(x) Ea i Eb j ε ij k Fab k 6πG Σ dete ] ( + γ )Ka i Kj b] where field strength is expressed as F k ab = a A k b b A k a + ǫ k ija i aa j b. (A9) (A) With use of (A),(A) and (A) the Hamiltonian (A9) assumes the form H G = V 8πGγ p c (A) where we have assumed a gauge of N(x) =. The constant V is the volume of fiducial cell. It is convenient to absorb the factor V by redefinition p = pv /, c = cv / (A) The holonomy along a curve α Σ is defined as follows h α = P exp τ i A i a dxa α (A) where iτ i = σ i and σ i are the Pauli matrices. From this definition we can calculate holonomy in the direction o e a i a and the length µv / h (µ) i / µv = exp = Icos ( µc τ i cv / o ωadx i a = exp τ i µc ) + τ i sin ( µc ) (A) where we used definition of the Ashtekar variable A (A). From such a particular holonomies we can construct holonomy along the closed curve α = ij. This is schematically presented on the diagram below.

11 o e a i a k i = i o e a j a h (µ) j h (µ) i h (µ) i h (µ) j o e a j a j The area i should in the limit i which corresponds to µ µ. Now, we must connect the area i with the length µ. We can choose that area i to correspond to physical area a µ or to the area µ. So in the case i = a µ we have the limit o e a i a µ µ = p. (A8) This holonomy can be written as h (µ) ij = h (µ) = e = exp i h (µ) j h (µ) i h (µ) j / µv A o a e a i e µv / A o a e a j e µv / A o a e a i e µv / A o a e a j µ V / A l a Am b o e a i o e b j τ l, τ m ] + O(µ ) = I + µ V / F k ab τ k o e a i o e b j + O(µ ) (A5) where we have used Baker-Campbell-Hausdorff formula and the fact that for flat FRW, field strength (A) simplifies to the form Fab k = ǫk ij Ai a Aj b. Now, equation (A5) can be simply rewritten to the form ( )] tr τ k h (µ) Fab k ij I = lim o ω µ µ V / a i o ω j b. (A6) The trace in this equation can be calculated with use of definition (A) ( )] tr τ k h (µ) ij I = ǫ kij sin (µc). (A7) In Loop Quantum Gravity the limit µ in the formula (A6) does not exist because of existence of the area gap. The area gap corresponds to the minimal quanta of area = πγlpl. So instead of the limit in equation (A6) we should stop shrinking the loop at the appropriate minimal area. This area corresponds to the area intersected by the loop. For the holonomy in the direction o e a i a the area is i, as explained in the diagram below. ] This approach we call µ scheme. In the case when i = µ we obtain µ µ = µ =, (A9) Which we call µ scheme. In the quantum version, we can combine equations (A6) and (A7), and write Fab k = sin ( µc) ǫ o µ V / kij ωa i o ω j b, whereas in the classical case we have F k ab = c V / ǫ kij o ω i a o ω j b. (A) (A) So, from these two equations we see that quantum effects can be introduced by a replacement c sin ( µc) µ (A) in the classical expressions. The Hamiltonian (A) with holonomy correction takes the form H eff = ] sin ( µc) p 8πGγ. (A) µ ] A. Ashtekar, T. Pawlowski and P. Singh, Phys. Rev. Lett. 96 (6) arxiv:gr-qc/686]. ] A. Ashtekar, T. Pawlowski and P. Singh, Phys. Rev. D 7 (6) 8 arxiv:gr-qc/6]. ] A. Ashtekar, T. Pawlowski and P. Singh, Phys. Rev. D 7 (6) 8 arxiv:gr-qc/679]. ] A. Ashtekar, T. Pawlowski, P. Singh and K. Vandersloot, Phys. Rev. D 75 (7) 5 arxiv:gr-qc/6]. 5] L. Szulc, W. Kaminski and J. Lewandowski, Class. Quant. Grav. (7) 6 arxiv:gr-qc/6].

12 6] L. Szulc, arxiv:77.86 gr-qc]. 7] D. W. Chiou, arxiv:7.6 gr-qc]. 8] D. W. Chiou, arxiv:gr-qc/7. 9] M. P. Dabrowski and T. Stachowiak, Annals Phys. (6) 77 arxiv:hep-th/99]. ] M. Bojowald and G. M. Hossain, arxiv:79.65 gr-qc]. ] J. Mielczarek and M. Szydlowski, Phys. Lett. B 657 (7) arxiv:75.9 gr-qc]]. ] J. Mielczarek and M. Szydlowski, arxiv:7.7 grqc]. ] M. Bojowald, H. H. Hernandez, M. Kagan, P. Singh and A. Skirzewski, Phys. Rev. D 7 (6) 5 arxiv:gr-qc/6957]. ] T. Stachowiak and M. Szydlowski, Phys. Lett. B 66 (7) 9 arxiv:gr-qc/6]. 5] P. Singh, K. Vandersloot and G. V. Vereshchagin, Phys. Rev. D 7 (6) 5 arxiv:gr-qc/66].