Loop Quantum Cosmology with a noncommutative quantum deformed photon gas

The noncommutativity of the space-time had important implications for the very early Universe, when its size was of the order of the Planck length. An important implication of this effect is the deformation of the standard dispersion relation of special relativity. Moreover, in the Planck regime gravity itself must be described by a quantum theory. We consider the implications of the modified dispersion relations for a photon gas, filling the early Universe, in the framework of Loop Quantum Cosmology, a theoretical approach to quantum gravity. We consider three types of deformations of the dispersion relations of the photon gas, from which we obtain the Planck scale corrections to the energy density and pressure. The cosmological implications of the modified equations of state are explored in detail for all radiation models in the framework of the modified Friedmann equation of Loop Quantum Cosmology. By numerically integrating the evolution equations we investigate the evolution of the basic cosmological parameters (scale factor, Hubble function, radiation temperature, and deceleration parameter) for a deformed photon gas filled Universe. In all models the evolution of the Universe shows the presence of a (nonsingular) bounce, corresponding to the transition from a contracting to an expanding phase.

Cosmology is undergoing presently a very rapid development. Very precise observational data, obtained from the study of both Cosmic Microwave Background (CMB) [1][2][3], and the matter distribution in the Universe [4,5], did offer the possibility of testing sophisticated cosmological models. Moreover, the high level of accuracy of cosmological data allows the testing of some of the predictions of elementary particle physics, and string theory.
The inflationary cosmological paradigm, introduced in 1981 by Alan Guth [6] is one of the cornerstones of modern cosmology. It attempts to give a solution to some of the major, and yet unsolved, problems of Standard Big Bang Cosmology (horizon, flatness, magnetic monopole, entropy etc.) by postulating the existence, during the very early period of the existence of the Universe, of an epoch of de Sitter type exponential expansion. Inflation also provides a theoretical explanation for the origin of inhomogeneities in the Universe [7][8][9][10][11], and it predicts that the spectrum of the primordial perturbations is very close to scale invariance, a result which is consistent with present day observations [3]. Hence the primordial power spectrum predicted by inflation can thus explain the late formation of large scale structures in the Universe, as well as the small inhomogeneities present in the Cosmic Microwave Background. For reviews of the inflationary paradigm and of the corresponding cosmological models of the early Universe see [12][13][14][15][16].
However, the inflationary paradigm also faces a number of fundamental theoretical problems, like, for example, the initial singularity problem, and the trans-Planckian problem, which originated from the fact that the origin and early evolution of the cosmological fluctuations must have taken place in the trans-Planckian regime, where both General Relativity as well as quantum field theory break down [17]. For a review of the fundamental conceptual problems that confront inflationary cosmology see [18,19].
In the standard cosmological models based on the theory of general relativity a fundamental assumption is that (ordinary) matter satisfies the standard energy conditions. In this case the Penrose-Hawking singularity theorem [20,21] is valid, and therefore every expanding cosmological solution (isotropic or anisotropic, homogeneous or inhomogeneous) of general relativity must had began with a singularity. However, from a physical point of view, the existence of the strong curvature singularity may be just the result of the extension of general relativity, an essentially classical theory, beyond its domain of applicability. According to the basic rules of quantum theory, general relativity, as a classical theory, must break down at length and time scales smaller than the Planck length and time, respectively. Hence, assuming the correctness of the quantum theory, it follows that general relativity cannot correctly describe the physics near the singular state.
To solve some of these theoretical problems raised by the description of the early Universe, an alternative scenario, called Bouncing Cosmology, was developed within the framework of the classical theory of gravity [22][23][24][25] (for recent reviews see [17,[26][27][28][29]). There are two basic approaches to Bouncing Cosmology, the introduction of a nontrivial stress-energy tensor in the matter sector of the Einstein field equations, or the modification of the gravitational Lagrangian in the Einstein-Hilbert action. The classical bouncing scenarios solve some of the puzzles of the standard model of cosmology [23], and they also can alleviate the singularity problem by postulating that the Universe begins its evolution with a contracting phase, and then it bounces back to the expanding phase before attaining the Planck length. However, these model also have serious drawbacks, like, for example, the instability problem [29]. But, most importantly, since the Universe attains extremely high densities very early, it follows that quantum gravity effects would necessarily become extremely important at such small length scales, and high energies. Therefore, a quantum gravitational description must be necessarily used to explore the high density regimes. Hence, if one would to solve the singularity problem, the use of quantum gravity and/or quantum mechanics becomes inevitable.
One of the promising candidates for a quantum gravity theory is Loop Quantum Gravity (LQG), whose starting point is the Ashtekar formalism of general relativity and the introduction of loop presentation [31]. LQG starts from the fundamental assumption that even at the quantum level the gravitational force is an expression of the geometry of the spacetime. From this postulate LQG advances further and develops a systematic and mathematically rigorous theory of quantum geometry, with a specific emphasis on the Riemannian one. For an in depth introduction to LQG see [32]. Generally, LQG represents a nonperturbative background approach to quantize gravity [33][34][35]. The theory fully takes into account the possible existence of fundamental discrete structures at the Planck scale, where all the geometric observables, including physical surfaces and volumes of bounded regions, are essentially quantum in their nature [36][37][38].
The application of the LQG quantization methods to spacetimes with reduced symmetry, and in particular to homogenous and isotropic spacetimes, led to the development of the field of Loop Quantum Cosmology (LQC) theory [39]. One of the major successes of LQC is the possibility of solving the singularity problem. This is due to the important result according to which the macroscopic physical observables, like, for example, the energy density, and the curvature, which in general relativity do diverge at the initial moment of the Big Bang, have a finite upper bound in LQC. The finite nature of the macroscopic quantities follows from the fact that they depend on the fundamental microscopic parameter of the LQC theory, the fundamental area gap, whose smallest eigen value is always nonzero [40][41][42]. Consequently, since a maximum value of the energy density can be obtained in LQC, a contracting Friedmann-Lemaitre-Robertson-Walker) (FLRW) Universe will bounce back to an expanding Universe, thus preventing the presence of a physical singularity. It is important to note that this result is derived without introducing any nontrivial or exotic forms of matter into the cosmological models. This represents an important difference as compared to the case of the classical bouncing scenarios. The quantum bounce is generated by the quantum geometrical effects that act as a novel repulsive force, a result that follows immediately from the quantum corrected Friedmann and Raychaudhuri equations. Also, it is important to point out that in the framework of LQC in all classes of spacetimes with different sets of symmetries, including Bianchi and Gowdy models, the singularity problem is successfully solved in most cases [43]- [64]. The resolution of sudden singularities, corresponding to the case when the pressure diverges, but the energy density approaches a finite value depends on the ratio of the latter to a critical energy density of the order of the Planck density [43]. On the other hand if the value of this ratio is greater than unity, the Universe avoids the sudden future singularity, and enters into an oscillatory phase.
On the other hand one can extend the Loop Quantum Cosmology matter bounce scenario by including a dark energy era, which ends abruptly at a Big Rip singularity, where both the scale factor and the Hubble function do diverge [65]. The scalar tensor reconstruction method in LQC as well as some scenarios for which the Hubble function becomes unbounded at finite time, were discussed in [66]. This corresponds to a case in which a Big Rip occurs. For reviews of the solution of the singularity problem in LQC see [67] and [68].
In addition to the singularity problem, cosmological models face an other important challenge. In order for the inflationary Universe models to be consistent with the present day observations, the Universe must have expanded at least 50 e-folds. On the other hand, if, for example, the Universe had expanded more than 70 efolds, (a value that is required by a large class of in-flationary cosmological models [69]), it follows that the wavelengths of all fluctuation modes, which currently are inside the Hubble radius, must have been smaller than the Planck length at the beginning of the inflationary era. This represents the so-called trans-Planckian problem [70]. Hence, the assumption that spacetime can still be treated classically, but that matter fields are essentially quantum in nature cannot be maintained. Therefore a quantum treatment of geometry must be also necessary. Moreover, the dynamics of the Universe before the beginning of the inflation cannot be neglected anymore, even when the modes were well inside the Hubble radius [70][71][72][73].
Hence all these problems requires to look for a quantum theory of gravity and cosmology. In this framework LQC deals successfully with one of the major problems in cosmology, by replacing the singularity with a quantum bounce, followed by a slow roll inflation [74]. Presently, it is already possible to look for experimental tests of LQC, and to search for observational signatures of the quantum bounce and of the pre-inflationary dynamics in current or forthcoming observations [75][76][77]. Calculations of cosmological perturbations for dressed metrics and for deformed algebras [67,78] have been performed and investigated numerically in [79,80]. As a result it was found, for example, that the approach based on deformed algebras is inconsistent with some of the present day observations [81].
The assumption of the quantum nature of the spacetime at the Big Bang immediately raises the problem of its non-commutativity. Theoretical models derived within the framework of space-time non-commutativity are of particular theoretical importance. In 1947, spacetime non-commutativity was already suggested by Snyder [82]. Snyder put forward a model to quantize spacetime, by defining space and time coordinate as operators that keep invariant a five-dimensional quadratic form when operating on it. The five-dimensional quadratic form considered by Snyder can be regarded as a real four-dimensional space with constant curvature, which is a De Sitter space [83]. The theory suggested by Snyder is invariant under Lorentz transformations, however, it is not invariant under space-time translations. Later, C. N. Yang [84] proposed another theory of quantization of space-time, by defining L i / , in which L i is the angular momentum, M i / , in which M i is the invariance quantity corresponding to Lorentz boost, Rp µ / , in which p µ is the momentum 4-vector and x µ /a, in which x µ is the coordinate 4-vector, as the 15 operators that do not change a six-dimensional quadratic form when operating on it, thus maintaining both Lorentz invariance and translation invariance when implementing discreteness to space-time. Recently, it has been found that space-time non-commutativity is a very promising way of solving some of the fundamental problems of current physics [85]. It plays an important role in the description of space-time at scales comparable to the Planck length [85]. Space-time noncommutativity determines the exis-tence of a minimum length scale, with important implications in quantum-mechanics, quantum-electrodynamics, thermodynamics, black-hole physics and cosmology [86][87][88][89] etc. Its connection with varying speed of light theories was investigated in [90]. Space-time noncommutativity also leads to modified energy-momentum dispersion relations, with important implications in physics, cosmology, and astrophysics [91][92][93][94][95]. For example, the thermodynamics of the radiation gas with modified dispersion relation E 2 −k 2 = m 2 (1 − E/κ) 2 , where E is the photon energy, k is its momentum and κ is a constant parameter was studied in [92], while the implications on the thermodynamics of the modified dispersion relation n ], in which α and E p are constant parameters were studied in [93]. Modified dispersion relation at energy scales near the Planck energy scale have been used to study black hole radiation in [94]. Newtonian white dwarfs composed of degenerate relativistic electron gas describing by modified dispersion relation is studied in, leading to the deformation of mass-radius relation [95].
An improvement of the usual LQC computations of the primordial power spectra in the inflationary framework was realized in [96], by taking into account the transplanckian effects by using some modified dispersion relations. In particular, it was found that the unphysical behavior of the deformed algebra spectrum with a quadratic potential and without backreaction can be fixed by an Unruh-like modified dispersion relation. On the other hand if an exponentially decreasing modified dispersion relation is considered, the nature of the dressed metric spectrum can become unphysical. The consideration of the modified dispersion relations can bring about significant changes in the power spectrum, which, depending on the details of the considered models, may become either compatible or incompatible with the observational data.
The Loop Quantum Cosmology of a flat Friedmann-Lematre-Robertson-Walker space-time with a Maxwell type photon field was studied in [97]. It turns out that many of the qualitative properties derived for the case of a massless scalar field are still valid for the case of radiation, with the big-bang singularity replaced by a quantum bounce. Moreover, the operator corresponding to the matter energy density is bounded above by the same critical energy density. Evolution equations that very closely approximate the full quantum dynamics of sharply peaked states at all times were also derived, and the analytical and numerical methods used to study other forms of matter fields were improved.
It was suggested in [98] that in the limit of high energies loop quantum gravity may lead to the deformation of the local Poincaré algebra. The deformation is a direct consequence of the quantum modification of the effective off-shell hypersurface deformation algebra. Interestingly enough, the form of the deformation indicates that at large curvatures the signature of the space-time may change from Lorentzian to Euclidean. Particular realizations of the loop-deformed Poincaré algebra were also constructed, and it turns out that they can be related to the curved momentum space. By using the corresponding energy dependence of the group velocity one can obtained some constraints on loop quantum gravity effects by investigating the time lags of high energy photons detected from distant astrophysical sources, like, for example, the gamma ray bursts.
It is the main goal of the present paper to investigate the evolution of the very early Universe in the framework of the Loop Quantum Cosmology by assuming as a matter source a deformed photon gas. When the energy scale approaches the Planck energy the quantum effects on the very structure of the space-time must also be taken into consideration. One of the important physical implications of the space-time noncommutativity is the modification of the dispersion relation of particles, which has important implications on the statistical mechanics of particles.
In the following we will investigate the statistical physics and LQC implications of three different models of non-commutative space-times, which lead to modified dispersion relations. The first one has the space-time commutators described by [x i , t] = iλx i and [x i , x j ] = 0, whose astrophysics properties were studied in [99]. The second and the third model are described by Modi- in which β and α are defined by β 0 /M 2 p c 2 and α 0 /M p c [100]. The modified Heisenberg relation is related to the non-commutative space via the Jacobi identity. The dimensionless parameter β 0 and α 0 can be constrained by many condensed matter experiments, astrophysical experiments and high energy physics experiments [100][101][102][103] etc. For example, one of the astrophysical experiment that determines the constraints on α 0 and β 0 is the gravitational wave event GW150914 [100]. As a next step in our study we consider the statistical mechanics of the deformed photon gas satisfying modified energy-momentum dispersion relations. We begin by writing down the general expression of the partition function Z for a thermodynamic system, which allows us to obtain the general expressions of the energy density and pressure of the photon gas satisfying modified dispersion relations. The general formalism allows us to obtain the equations of state of the radiation in the cosmologically relevant limit of high temperatures and densities, or, more exactly, near the Planck scale regimes. The presence of the noncommutativity of the space-time leads to significant modifications in the temperature dependence of the photon gas energy density and pressure, as compared to the standard quantum case of the black body radiation.
At the very early stages of evolution of the Universe, the general relativistic Friedmann equations are not valid since the quantum effect of gravity must be taken into consideration. Loop Quantum Cosmology, as an application of Loop-Quantum-Gravity, has been widely used to investigate the Big-Bang singularity. By investigating different LQC models, it was found that the Big-Bang singularity is replaced by a Big Bounce [44]- [61], which offers a solution to one of the most important problems in theoretical cosmology. In the present paper we investigate the LQC-modified Friedmann equations, representing effective cosmological evolution equation derived from LQG, together with the assumption that the Universe was initially filled with a photon gas described the deformed statistical mechanics induced by the noncommutativity of the space-time. We consider explicitly four cases of deformed dispersion relations, and in each case the energy density and pressure of the photon gas are obtained in the high (Planck) temperature limit. Then the modified Friedmann equations are investigated numerically by using the deformed radiation fluid equations of state. The time behaviors of the scale factor, Hubble function, temperature (energy density), and deceleration parameter are considered in detail for different numerical values of the model parameters. The dynamical evolution is studied beginning with an initial contracting phase, in which the scale factor is a monotonically decreasing function of time. After the scale factor reaches a minimum finite value, the Universe enters in an expanding phase. The maximum finite temperature (energy density) is reached at the moment of the bounce, and in the expanding phase the energy density of the photon gas is monotonically decreasing. The behavior of the decelerating parameters indicates that in the contracting phase the Universe did experience a super-accelerated contraction, with the deceleration parameter having high negative values, while during the expanding phase the deceleration parameter reaches positive values, indicating a decelerating behavior.
The present paper is organized as follows. In Section II, we briefly review the modified dispersion relations to be considered, as well as the general formulation of the deformed quantum statistical mechanics. We introduce the basic ideas of LQG and the LQC-modified Friedmann equations in Section III, and formulate them in a dimensionless form. In Section IV, the cosmological evolution of the deformed radiation filled very early Universe is investigated for three different photon dispersion relations by using the LQC-modified Friedmann equations. In each case the cosmological evolution equations are integrated numerically, and the time behaviors of the scale factor, of the Hubble function, of the temperature of the Universe, and of the deceleration parameter are studied in detail. We discuss and conclude our results in Section V. A closed form representation of the equations of state of the deformed photon gas with the dispersion relation E = kc 1 − 2βk 2 is presented Appendix A. The invariant phase space volume for the photon dispersion relation E = kc(1 − αk) is calculated in Appendix B.

II. MODIFIED DISPERSION RELATION FROM NON-COMMUTATIVE SPACE-TIME, AND DEFORMED STATISTICAL MECHANICS
The relation between non-commutative nature of the space-time and the modified dispersion relations can be understood in various way. Quantum group methods can be used to find the modified dispersion relation with given commutator of space-time [104]. Another approach to understand this relation is via the minimum length implied by the non-commutativity. The four-dimensional momentum operator should be modified since the spacetime transformation is only allowed for the multiples of the minimum length. Imposing Lorentz invariance with modified four-dimensional momentum leads to deformed dispersion relations. Another approach to investigate the relation between non-commutative space-time and modified dispersion relation is using non-commutative geometry [105,106].

A. Modified dispersion relations
In this work, three models of non-commutative spacetime will be discussed, which lead to three deformed dispersion relations.

The first non-commutative space-time model
The first non-commutative space-time model is described by the commutators [99] in which λ is a constant parameter. The dispersion relation of above commutators can be studied by κdeformed Poincaré group, with space-time commutators [x 0 ,x i ] = i κx i and [x l ,x j ] = 0, which was studied in [90]. Using the invariant wave operator in the κ-deformed Minkowski space [90], in which for massless particles, the dispersion relation is obtained in [90] as Considering the varying speed of light generalization E 2 − p 2 c 2 f 2 = 0 of the above dispersion relation [90], taking f (E) to be f = 1 + λE [107], the dispersion relation corresponding to the first non-commutative spacetime model considered in this paper can be written as [90,99,107] 2. The second non-commutative space-time model The second non-commutative space-time model is described by the commutators [100,108] The operators for position x i and momentum p i are defined as [100] x i = x 0i , in which x 0i and p 0i satisfy the canonical commutation relation. Using the relation and retaining the terms up to O(β 2 ) only, the above equation can be written as [100] p a p a = −(p 0 ) 2 + p 2 + 2βp 2 p 2 = −m 2 c 2 + 2βp 2 p 2 . (9) Equivalently, the equation can be reformulated as [100] ( Thus we obtain the dispersion relation [100] 3. The third non-commutative space-time model The last non-commutative space-time model investigated in this paper is described by the commutators [100,109,110] (12) By using the similar approach as that we used for the second non-commutative space-time model, the deformed dispersion relation for this model can be obtained as [100,111]

B. Deformed quantum statistical mechanics
In the present work we will investigate the ways the non-commutative space-time influences the cosmological evolution in Loop Quantum Cosmology. In our analysis we will concentrate mostly on the dynamics of a quantum Universe filled with a radiation fluid. As a first step in our study we need to consider the formulation of the quantum statistical mechanics with modified energymomentum dispersion relations.
According to the basic rules of statistical physics, the grand canonical partition function Z can be written as [112] where n j , j = 1, 2, ...N are the particle numbers, E j are the corresponding energies of the quantum states, µ is the chemical potential, while T and k B denote the temperature, and Boltzmann's constant, respectively. For radiation, µ = 0, and therefore we immediately obtain and respectively. When the space-time quantum commutators are deformed, the volume of the invariant phase space is also deformed. We assume that the deformed volume of the invariant phase space can be written as where k denotes the particle momentum, and α is a small parameter in the modified dispersion relation. The function f k, α must satisfy the condition lim α→9 f k, α = 1, a condition ensuring that the invariant phase space volume goes back to d 3 kd 3 x in the limiting case when the dispersion relation takes its standard form E = kc. The energy degeneracy of each state is obtained as in which the factor 2 is due to two polarization states of the photon. In the limit of a continuous system, consisting of identical particles, ln Z can be written in an integral form as The energy density u of the photon gas can be computed as The pressure of the gas can be obtained in the form Since the last term in the above equation is rather small as compared to the first one, in the following we will neglect it in the expression of the pressure of the deformed radiation fluid. Thus, the basic thermodynamic parameters of a photon gas in the deformed statistical mechanics are obtained as The value of the maximum momentum k max , the volume of the invariant phase space d 3 xd 3 k/f k, α , as well as the expression of the dispersion relation E = E(k) are different for the various models of the non-commutative space-time. In the following, by using the expressions (22) and (23) of the thermodynamic parameters of the photon gas in the deformed statistical mechanics, we will investigate the influence of the non-commutative structure of the space-time on the cosmological evolution of Loop Quantum Cosmological models.

III. LOOP QUANTUM COSMOLOGY OF THE DEFORMED RADIATION GAS-THE MODIFIED FRIEDMANN EQUATIONS
In a radiation dominated flat FLRW universe, the basic cosmological equations describing the early evolu-tion of the Universe during the quantum gravity dominated phase are the Loop Quantum Cosmology modified Friedmann equations. Loop Quantum Cosmology, which solves the Big-Bang singularity problem by predicting a quantum gravity "bounce", connecting the expanding epoch with the contracting epoch, is an application and an effective theory of Loop Quantum Gravity. The LQC modified Friedmann equations are only qualitatively valid in the early Universe, and they can be applied when the energy scale is comparable to the Planck scale, since in this effective theory the degrees of freedom with large quantum fluctuations are artificially suppressed.

A. Radiation fields in cosmology
Before going into Loop Quantum Cosmology, let's first have a look at the classical formalism of cosmology. In the following we assume that the geometry of the Universe can be described by the flat Friedmann-Lemaitre-Robertson-Walker (FLRW) metric, given by where a is the dimensionless scale factor. We also introduce the Hubble function, defined as H =ȧ/a. In such a k = 0, FLRW universe, the scale factor ν and its canonical variable b suffice for describing this symmetryreduced system [97]. The definitions of ν and b are [97] where α = 2πγl 2 P l √ ∆, and l P l is the Planck length, defined as l P l = G /c 3 . γ is the Barbero-Immirzi parameter, satisfying the constraint ln 2/π ≤ γ ≤ ln 3/π [113,114]. However, in the following for explicit applications we fix the value of γ as γ = 0.2375. ∆ is the smallest nonzero eigenvalue of the LQG area operator, with the order of magnitude ∆ ≈ l 2 P l [115]. The two variables describing the gravitational field in this system satisfy the relations {ν, b} = −2/ [97]. The Hamiltonian constraint density of the gravitational part is given by − 3πG 2 2α |ν|b 2 [97]. In a radiation dominated homogeneous and isotropic Universe, the photon field can be considered as the superposition of three U (1) vector fields, with field strengths given by [97] The Hamiltonian constraint density of matter part is given by The total Hamiltonian constraint is given by the combination of the Hamiltonian constraint of gravitational part and matter part, which is [97] The Hilbert space of quantized gravitational space-time consists of all those states that can be annihilated by the quantum Hamiltonian constraint.

B. Loop Quantum Cosmology and modified Friedmann equations
Next we discuss the Loop Quantum Cosmology of the spatially flat, homogeneous and isotropic Universe filled with a radiation fluid. In the flat FLRW universe dominated by radiation, the kinematic Hilbert space is the product of gravitational and matter Hilbert space [97] . Taking ν as the elementary variable for quantization, the gravitational Hilbert space can be formed by the eigenfunctions ofν, which satisfy the condition < ν|ν ′ >= δ ν,ν ′ [97] . For the matter Hilbert space, the elementary variable for quantization is the field A γ , giving the Hilbert space formed by the eigenfunctions of operatorÂ γ , represented as |A γ > [97] . Working in the above representation, the Hamiltonian constraint in Loop Quantum Cosmology can be written as [97] where The physical states are those annihilated by the quantum Hamiltonian constraint. Thus, these states are given by It is possible to derive an effective theory of the above LQC in the k = 0, FLRW radiation dominated universe. The effective Hamiltonian constraint can be written as [97] From the effective Hamiltonian constraint we can obtain the Friedmann equation [97] dν which is equivalent to with (34) Another equation that must be considered during the cosmological evolution is the energy conservation equation, given byρ Taking the time derivative of the Friedmann equation (33), and using the energy conservation equation (35), it can be easily shown that the time evolution of the Hubble parameter is given bẏ

Dimensionless form of the LQC evolution equation
In order to obtain a dimensionless form of the Friedmann and energy conservation equations, we introduce a set of dimensionless variables (r, P, h, τ ), defined as or, equivalently, where the Planck time t P l is defined as t P l = l P l /c = 5.391 × 10 −44 s. Then the cosmological evolution equations (33), (35) and (36) take the following dimensionless form, An important cosmological quantity, the deceleration parameter q, is defined according to In this work, Eqs. (40), (41) and (42) will be used to investigate the cosmological evolution of the radiationdominated flat FLRW Universe filled with a deformed radiation fluid, in its early epochs described by the Loop Quantum Cosmology approach, and in the presence of a bounce.

IV. COSMOLOGICAL EVOLUTION OF THE LOOP QUANTUM COSMOLOGY MODELS WITH DEFORMED RADIATION GAS
In the present Section we will consider the cosmological evolution of some classes of Loop Quantum Cosmological models, whose dynamics is described by the modified Friedmann equations, under the assumption that the matter content of the early Universe consisted of radiation only. Due to the quantum nature of the very early Universe, we will also take into account the quantum effects due to the noncommutative nature of the space-time, which induces a modification of the energymomentum dispersion relation for photons, and, consequently, of the quantum statistical mechanics of the radiation gas. In our analysis we will investigate the cosmological implications of three such modifications of the statistical properties of the photon gas, which are closely connected with the modifications of the dispersion relations due to the quantum nature of the space-time.
A. LQC implications of the dispersion relation E = kc (1 + λE) As a first example of a cosmological model in LQC with a modified photon gas we consider the case in which the deformed dispersion relation for photons is given by where λ is a constant. This dispersion relation is particular case of the more general relation [90] for the case of the massless particles, with m = 0.

Thermodynamics of the deformed photon gas
The energy density of the radiation fluid assumed to fill the Universe when space-time noncommutativity effects are important is given by [99] By introducing a new variable x = λE, Eq. (46) can be rewritten as The pressure of the deformed radiation gas is given by [99] or, equivalently, We will consider the above deformed photon energy density and pressures relations in two limiting cases.
a. The case λk B T >> 1. First, we assume that the condition λk B T >> 1 holds. Then, by expanding the exponential in a power series, and keeping only the first order approximation, gives For the pressure of the gas, under the assumption λk B T >> 1, we obtain where in order to avoid the divergence in the integration we have extended the limit to the maximum value of x = λE only. By adopting for x max a value of the order of x max = λk B T >> 1, we finally obtain the approximate equation of state of the deformed radiation gas as given by b. The case λk B T << 1. In the opposite temperature limit, with λk B T << 1, we expand f (x) = x 3 / (1 + x) 4 in Taylor series around zero. Far away from the origin, the other terms will be exponentially suppressed by the term e x/λkB T − 1 −1 . By taking into account the formula [116] ∞ 0 where Γ (ν) and ζ (µ) are the gamma and zeta functions, respectively, we obtain For the pressure of the radiation fluid, by taking into account again that λk B T << 1, we easily find where g(x) = x 3 / (1 + x) 3 .

Cosmological dynamics in the case λkBT >> 1
We proceed now to the investigation of the cosmological evolution of Loop Quantum Gravity models with deformed radiation satisfying the dispersion relation (44). We will consider first the case λk B T >> 1. In this case the energy density and the pressure of the radiation can be expressed in the dimensionless form where we have denoted Then Eq. (41) can be written as The energy conservation equation (42) takes the form By using Eqs. (43) and (90), the expression of the deceleration parameter q can be obtained as To solve the cosmological evolution equations numerically, we need first to estimate the value of the dimensionless parameter A 1 . In the expression of A 1 , as given by Eq. (57), ρ max is the maximum density of the Universe, which in the framework of LQC can be obtained as [97] where the Planck length l P l is defined as l P l = G /c 3 . Substituting the expression of the Planck length, and putting the value of the speed of light back into the expression of ρ max , we obtain for ρ max the expression, As for A 1 and θ, they can be written as and GeV is the Planck energy, with M P l denoting the Planck mass M P l = c/G. The condition λk B T >> 1 can be written as and it can be satisfied in two cases: if T ≈ T P l , λE P l >> 1, or if T >> T P l , λE P l ≈ 1, respectively. Taking γ = 0.2375 [114] and λ = 10 −19 GeV −1 [117], the numerical value of A 1 can be obtained as The cosmological evolution of a Universe filled with a deformed radiation fluid satisfying the generalized dispersion relation (44) and in the limit λk B T >> 1 is represented in Figs. 1 and 2, respectively. To obtain the cosmological parameters of interest (scale factor, Hubble function, temperature and deceleration parameter) we have numerically integrated Eqs. (90) and (91), for different numerical values of the dimensionless parameter A 1 . In all plots we have fixed the initial value of the temperature θ 0 , and obtained the initial value of the dimensionless Hubble function as h 0 = − A 1 θ 0 (1 − A 1 θ 0 ).
As one can see from the left panel of Fig. 1, the evolution of the Universe begins in a contracting phase, with the scale factor a monotonically decreasing in time. At the time τ = τ 0 (the bouncing time) the scale factor reaches its minimum finite value, and the Universe bounces into an expanding phase, with the scale factor becoming a monotonically expanding function of time. The cosmological evolution is strongly dependent on the numerical values of the parameter A 1 , which combines the LQC parameter γ and the deformation parameter λ. The Hubble function, depicted in the right panel of Fig. 1, shows a complex behavior. Generally, it decreases at the initial stages of the contraction, when it takes negative values. After reaching a local minimum, the Hubble function begins to increase, and becomes zero at the bouncing time, when the scale factor takes its minimum value. After the bounce the Hubble function still continues to increase from the zero to a maximum value, corresponding to a moment τ = τ max . After that the Hubble function becomes a monotonically decreasing function of time. During the increasing phase, the cosmological evolution is strongly dependent on the numerical values of A 1 , but in the decreasing phase there is no significant influence of the deformation parameter on the time variation of h. The temperature of the Loop Quantum Cosmological Universe, shown in the left panel of Fig. 2, increasing during the contracting phase, and reaches its maximum at the bounce. After that the temperature of the Universe decreases monotonically. The deceleration parameter q, plotted in the left panel of Fig. 2, indicates a complicated dynamical evolution. During the contracting phase q takes very high negative values, of the order of q ≈ −10 7 , indicating the existence of a superaccelerating phase. However, it increases rapidly, and around the bounce time τ 0 it reaches the marginally ac-celerating value q ≈ 0, after which q enters in the region of positive values, indicating a transition to the decelerating phase. After reaching a maximum in the positive domain, the deceleration parameter begins again to decrease. There is a strong influence of the numerical values of the parameter A 1 on q in both the contracting and expanding phases.

The cosmological evolution for λkBT << 1
In the case λk B T << 1, the dimensionless equation of state of the deformed radiation fluid is obtained from Eq. (55) as where we have kept only the first two terms in the general series expansion. The dimensionless energy density can be written as By introducing the new variable θ = λk B T , and by denoting Eqs. (67) and (66) can be written as and respectively. Substituting Eqs. (69) and (70) into Eqs. (41) and (42), the equation describing the evolution of h and θ can be obtained as and respectively. The deceleration parameter, defined by Eq. (43), is obtained as Under the above thermodynamics relations, it can be easily shown that only when A 2 ≥ 1.03 × 10 6 does the equation r − 1 = 0 have real roots in the real and positive domain. This means that only when A 2 ≥ 1.03 × 10 6 , the Big Bounce can happen. Equivalently, it means that the upper bound of λ is at λ ≈ 10 −21 GeV −1 . This gives for θ the numerical expression Substituting the expression of ρ max as given by Eq. (62) back to Eq. (68), and taking γ = 0.2375 and λ = 10 −21 GeV −1 , we obtain for the numerical value of A 2 the estimation A 2 ≈ 2 × 10 7 .
To solve the cosmological evolution equations numerically, the initial value of θ 0 is also needed. The case being considered now is λk B T << 1, or a condition which is satisfied for two choices of the parameters: T ≈ T P l , λE P l << 1, or T << T P l , λE P l ≈ 1. For A 2 ≈ 2 × 10 7 the condition is equivalent to T << 1.16 × 10 34 K. Also, when the numerical value of the parameter A 2 is of the order A 2 ≈ 10 7 , the maximum of θ occurs for values of θ ≈ 10 −2 . Considering that the Planck temperature is T P l = 1.42 × 10 32 K, the bounce occurs at the temperature scale of Planck temperature when the parameter A 2 is of the order of A 2 ≈ 10 7 . From the above considerations, the initial value of θ 0 in the contracting epoch is taken to be 0.009. The time evolution of the scale factor a, of the dimensionless Hubble function h, of the temperature of the Universe θ, and of the deceleration parameter q are presented in Figs. 3 and 4, respectively. They are obtained by numerically integrating Eqs. (71) and (72) for different values of the model parameter A 2 . In order to integrate the evolution equations we fix the value of the initial temperature θ (τ 0 ) = θ 0 , and we obtain the initial value h (τ 0 ) = h 0 of the dimensionless Hubble function in the contracting phase as As one can see from Figs. 3 and 4, the time evolution of the deformed radiation filled Universe with λk B T << 1 is qualitatively similar to the opposite temperature limit, considered previously. The Universe begins its evolution in a contracting phase, with the scale factor, shown in the left panel of Fig. 3 decreasing rapidly in time. At a moment τ = τ 0 , the scale factor reaches a minimum value, and the Universe bounces to an expanding phase. The  Hubble function h, presented in the right panel of Fig. 3, decreases during the contracting phase, and reaches a local minimum at the bounce time τ 0 , after which the decreasing phase is followed by an increasing one, in which h enters the positive value range, and it reaches a local maximum at a finite moment τ max . After reaching its maximum value, the Hubble function monotonically decreases in time. The temperature of the Universe, depicted in the left panel of Fig. 4, decreases rapidly during the contracting phase, and reaches its maximum at the bounce time, and then it decreases monotonically. The deceleration parameter, represented in the right panel of Fig. 4, takes very large negative values during the contracting phase, and reaches the region of positive values approximately at the moment when the Hubble function reaches its maximum value. During the expanding phase that follows after h has reached its maximum, the deceleration parameter is approximately constant, with values of the order of q ≈ 1, indicating a decelerating expansion of the Universe. Generally, in the early contracting and immediately post contracting phases, the cosmological evolution is strongly dependent on the numerical values of A 2 . However, for large times, the dynamics of h and q is independent on A 2 , but the scale factor and temperature evolutions still depend on A 2 .
B. LQC implications of the dispersion relation E = kc 1 − 2βk 2 As the second example of a cosmological model in LQC with a modified photon gas we consider the case in which the deformed dispersion relation for photons is given by where β is a constant. This dispersion relation is particular case of the more general relation [100]

Thermodynamics of the photon gas
For the case of the massless particles, with m = 0, the energy density of the photon gas is given by There are two things one should be mentioned about this integral expression: first, the classical degree of degeneracy of each phase unit d 3 xd 3 k 3 is modified to [118], due to the modification of the invariant phase space; secondly, the upper limit of integration is taken to be 1/2 √ β. When k = 1/(2 √ β), E = E max ; when k = 1/ √ 2β, E = 0. Which should be the maximum value of the momentum, k = 1/(2 √ β) or k = 1/ √ 2β ? The correct upper bound for the photon momentum is If the upper bound of the momentum is taken as k = 1/ √ 2β, the pressure will become negative, which is not physical.
By introducing a new variable x = kc/k B T , and de-notingT the energy density can be rewritten as where we have denoted The pressure of the deformed radiation gas is given by ) or, equivalently, by where we have denoted The variations of the functionT 2 I T and K T are represented in Fig. 5.
The numerical results can be fitted with a high precision by two simple functions, proportional to T . Hence we can represent the energy density and the pressure of the photon gas as and respectively. Exact closed form analytical representations of the energy density and pressure of the photon gas satisfying the modified dispersion relation E = kc 1 − 2βk 2 are presented in Appendix A.

Loop Quantum Cosmological evolution of the deformed radiation filled Universe with
We proceed now to the investigation of the cosmological evolution of Loop Quantum Gravity models with deformed radiation described by the modified dispersion relation E = kc 1 − 2βk 2 . In this case the energy density and the pressure of the radiation can be expressed in the dimensionless form where we have denoted Then Eq. (41) can be written as The energy conservation equation (42) takes the form dθ dτ + 3.5679hθ = 0.
By using Eqs. (43) and (90), the expression of the deceleration parameter q can be obtained as To solve the cosmological evolution equations (90) and (91) numerically, we need to obtain first the value of A 3 . Substituting the expression of ρ max as given by Eq. (62) back into A 3 , we find which is a dimensionless quantity. In the expression of A 3 , β can be written as β 0 /M 2 P l c 2 , where M P l is the Planck mass [100]. The Planck energy M P l c 2 is given by 1.2 × 10 28 eV. For the dimensionless quantity θ we find For γ = 0.2375 and β 0 = 10 21 (the upper bound of β 0 is obtained from electron tunnelling experiments) [101,102], the numerical value of A 3 can be obtained approximately as To investigate the cosmological evolution, the initial value of θ is needed. In the contracting epoch, the initial temperature is taken to be 0.1T P l , corresponding to θ (τ 0 ) = 3.22 × 10 9 . In the expanding epoch, the initial temperature (the temperature of the bounce) is taken to be the Planck temperature T P l = 1.42 × 10 32 K, corresponding to θ (τ 0 ) = 3.22 × 10 10 . Hence we investigate the evolution of the Universe within the Loop Quantum Cosmological framework with a modified radiation dispersion relation in the range of temperatures of 0.1T P l < T < T P l . The initial value of the Hubble function h (τ 0 ) = h 0 is determined by the initial temperature θ (τ 0 ) = θ 0 of the Universe as The time evolution of a, h, θ and of the deceleration parameter q are investigated numerically, for different values of A 3 . The results of the numerical analysis are presented in Figs. 6 and 7, respectively.
The cosmological evolution is qualitatively similar to the previously considered LQC models. The Universe begins its evolution in a contracting phase, with the scale factor, depicted in the left panel of Fig. 6, reaching a finite minimum value at some time τ min , after which the Universe bounces, and enters in an expansionary phase, in which the evolution of a is strongly dependent on the numerical values of the parameter A 3 . The Hubble function, shown in the right panel of Fig. 6, decreases during the contracting phase, and after reaching a local minimum it increases to a maximum value, followed by a monotonically decreasing phase. The temperature of the Universe, plotted in the left panel of Fig. 7, increases during the contracting phase, with the Universe reaching its maximum (Planck) temperature at the bounce. After that the temperature as well as the radiation energy density become some decreasing functions of the cosmic time. The time evolution of the deceleration parameter, presented in the right panel of Fig. 7, shows that in the contracting phase the Universe is strongly accelerating, with the deceleration parameter taking very high (negative values). However, after the bounce, the deceleration parameter increases, and enters in the positive values region, with values of the order of q ≈ 0.5, indicating a decelerating behavior. In the large time limit the variation of q is independent of the numerical values of the parameter A 3 .

C. LQC implications of the dispersion relation
As the third example of a cosmological model in LQC with a modified photon gas we consider the case in which the deformed dispersion relation for photons is given by  where α is a constant. This dispersion relation is particular case of the more general relation [100,111]

Energy density and pressure of the radiation fluid
For the case of the massless particles, with m = 0, the energy density of the photon gas is given by dk.
(99) There are two things one should be mentioned about this integral expression: first, due to the modification of the invariant phase space, the classical degree of degeneracy of each phase unit d 3 xd 3 k 3 is modified to The derivation of this relation is presented in detail in Appendix B. Secondly, the upper limit of integration is taken to be 1/2α. When k = 1/2α, E = E max ; when k = 1/α, E = 0. The correct upper bound for the photon momentum is thus k = 1/2α. If the upper bound of the momentum is taken as k = 1/α, the pressure will again become negative, which is not physical. By introducing a new variable x = kc/k B T , and defining the energy density of the photon gas can be written as (102) where we have denoted The pressure of the deformed radiation gas is given by or, equivalently, where we have introduced the function N θ), defined as dx.
The numerical results can be fitted with a high precision by two linear functions, proportional to θ. Hence we can represent the energy density and the pressure of the photon gas with the modified dispersion relation E = kc(1 − αk) as u = 0.04943 α 4 π 2 3 αk B T, p = 0.006317 π 2 3 α 4 αk B T.
2. Evolution equations for the radiation filled Universe with photon dispersion relation E = kc(1 − αk) a. The dimensionless equations of state. We investigate now the cosmological evolution of Loop Quantum Cosmology models with deformed radiation described by the modified dispersion relation E = kc(1 − αk). In this case the energy density and the pressure of the radiation can be expressed in the dimensionless form where we have denoted b. The generalized Friedmann equations. Then Eq. (41) can be written as The energy conservation equation (42) takes the form The expression of the deceleration parameter q can be obtained as To solve the cosmological evolution of the radiation fluid in the case E = kc(1 − αk) numerically, we need to obtain the value of A 4 . Substituting the expression of ρ max as given in (62) back into A 4 , the expression of A 3 can be obtained as which is a dimensionless quantity. In the expression of A 4 , α = α 0 /M p c [100]. By adopting for γ and α the values γ = 0.2375 and α 0 = 10 14 , respectively, with the upper bound of α 0 obtained from 87 Rb cold-atom-recoil experiments [103], the numerical value of A 4 can be estimated as For the dimensionless variable θ we obtain

The cosmological dynamics
The time evolutions of a, h, θ and of the deceleration parameter q, obtained by numerically integrating the generalized Friedmann equations (110) and (111) for different values of the parameter A 4 are presented in Figs. 9 and 10, respectively.
The cosmological evolution is qualitatively very similar to the previously considered LQC models, so we will not discuss it in detail. There is a typical bounce at a finite time, with the scale factor and the Hubble function reaching a local minimum, followed by a transition to an expanding state. However, despite the qualitative similarities, there are significant quantitative differences between the different models we have considered. In the present case, for a Universe filled with a radiation fluid satisfying the dispersion relation E = kc(1−αk), the Universe ends in a marginally inflating phase, with q ≈ 0. The post-bounce behavior of the deceleration parameter is practically independent on the numerical values of the model parameter A 4 . The photon gas reaches at the bounce the Planck temperature, which rapidly decreases in the post-bounce phase.

V. DISCUSSIONS AND FINAL REMARKS
In the present work we have investigated the very early stages of the cosmological evolution of the Universe. Since in this phase the energy scale of the matter reaches the Planck scale, the quantum effects of the gravity must be taken into consideration. Since a complete theory of quantum gravity is still missing, to describe the cosmological evolution at the Planck sale we have adopted the effective LQC approach that results in a set of modified Friedmann equations. Another class of quantum effects on the space time geometry,the non-commutative nature of the space-time at the Planck scale are also taken into consideration by using, for obtaining the thermodynamic parameters of the photon gas, the statistical mechanics of the deformed dispersion relations. Non-commutative space-time changes the statistical mechanics by modifying the photon dispersion relation, by restricting the range of the momentum by introducing an upper bound, and by modifying the invariant phase space, respectively.
From a cosmological perspective we have adopted the fundamental assumption that the early Universe consisted of a photon gas only. This physical model has been much less investigated in Loop Quantum Cosmology as compared to the alternative model in which the early composition of the Universe is modeled by a scalar field only [40][41][42]. However, it is well known that at some moment the Universe was dominated by radiation [119], and we can see the remnants of this radiation in the form of the CMB. In standard cosmology all the matter content of the Universe was produced during the postinflationary reheating phase by the decay of the inflaton field [120]. However, the idea of a Universe that may have been dominated from the very early beginning by quantum deformed radiation may represent an attractive alternative to the standard cosmological scenarios implying the presence of a primordial scalar field.
The first drastic effect of the noncommutativity of the space-time on the photon gas at the Planck scale is the modification of the equation of state, and especially of the temperature dependence of the fundamental thermodynamic parameters, the energy density and pressure. In standard quantum physics the radiation gas with a black body spectrum is described by the equation of state p = u/3, and with u ∝ T 4 . In the first model we have discussed, the dispersion relation is given by E = kc(1+λE), and we have discussed its thermodynamic properties in the limits λk B T >> 1 and λk B T << 1. In the high temperature limit the quantum properties of the space time modify the equation of state of the photons to u ∝ T , and ∝ T (ln T + A), where A is a constant. This equation of state is valid for the photon gas at the Planck temperature. In the opposite limit of "small" temperatures, T ≈ 0.1T P l , corrections to the standard photon equation of state do also appear, so that u ∝ π 2 /15 − AT T 4 , and p ∝ u (1/3 + BT ), where A and B are constants. The second model we have considered is described by the dispersion relation E = kc 1 − 2βk 2 , with the up-per bound of momentum given by 1/2 √ β, while the invariant phase space is given by d 3 xd 3 k (1+βk 2 ) 3 . In this case near the Planck temperature the thermodynamic parameters of the photon gas take the simple form u ∝ T , p ∝ T , and p ∝ u, with the proportionality constant in the equation of state different from 1/3. The last model is described by the dispersion relation E = kc(1−αk), with the upper bound of the momentum given by 1/2α, and the invariant phase space obtained as The behavior of the thermodynamic parameters is similar to the second model, with both radiation energy density and pressure proportional to the temperature, but with different proportionality coefficients. Even that the equations of state of the radiation fluid can be represented in a closed analytical form in all three cases, we have followed an alternative numerical approach, which is based on the high precision fit of the exact numerical results. For all the above described equations of state we have studied numerically the modified Friedmann equations of the Loop Quantum Cosmology. In all cases the behavior of the Universe shows qualitatively similar features. We began the investigation of the evolution of the Universe by fixing the initial conditions, which implies giving the values of the initial temperature and Hubble function of the Universe. These initial values fully determine the dynamics of the Universe. By assuming a negative value for H (τ 0 ) (contracting initial state), the bouncing evolution of the radiation gas naturally follows from the LQC modified Friedmann equations. In all the considered models the presence of a singular (geometric or physical) state is completely avoided, and the contracting phase is always followed by an expansionary one. While the scale factor has only one local minimum (the bounce), the Hubble function has both a local minimum and a local maximum. The moment of the maximum of the Hubble parameter occurs later as compared to the moment of the bounce of the scale factor. The energy density of the radiation, which also gives the temperature of the Universe, reaches its maximum value at the bounce. Generally, the behavior of all geometrical and physical quantities show a strong dependence on the numerical values of the basic dimensionless model parameter that combines the noncommutativity and the LQC parameters.
Even that qualitatively similar, the analyzed LQC cosmological models show important quantitative differences. These differences can be seen very clearly by analyzing the behavior of the deceleration parameter. As a common feature of all models, the contracting phase is a super-accelerating one, with q having very high negative values, of the order of q ≈ −10 8 . However, the Universe decelerates quickly, and in the post-bounce phase it is still in an accelerating, de Sitter type phase, with q of the order of q ≈ −1. The deceleration process continues, and in the large time limit the Universe enters in the positive deceleration parameters phase. This general picture is strongly dependent on the type of the noncommutativity, as well as of the model parameters. In the present analysis, far from the bounce, the decelera-tion parameter can take values in the range q ∈ [0, 3], and by slightly modifying the model parameters other ranges are also possible. Hence, once the behavior and the numerical values of the deceleration parameter in the very early Universe are known, the comparison of the observational data and the theoretical predictions could lead to the testing of quantum gravity and noncommutative space-time theories at the Planck scale. In the framework of a simple physical picture we may interpret the highly accelerating phase before the bounce as kind of pre-inflationary era that is necessary to "push" the Universe in the post-bounce inflationary phase. In the present model at the bounce and immediately after the Universe is still accelerating at a de Sitter rate, and this acceleration is the consequence of the high acceleration acquired in the previous phase, which can be described as a kind of "inertial" behavior when experiencing the transition from one phase to another.
One of the important tests of cosmological models is the study of their linear perturbations. The quantitative predictions of primordial curvature perturbations in LQC depend on the adopted cosmological scenario, and on the initial conditions. Generally, the evolution of the background can be divided into three different phases, bouncing, transition and slow-roll inflation. In the present analysis we have imposed some initial conditions for the cosmological evolution in the contracting phase, before the bounce. With this choice all the modes are inside the Hubble horizon and the Bunch-Davies vacuum for perturbations can be imposed. There are some other possibilities for fixing the initial conditions, like, for example, at the bounce, corresponding to the appearance of the Big Bang "singularity". However, the Bunch-Davies vacuum cannot be imposed here, as some modes are outside of the Hubble horizon. Another possibility is to adopt as the initial moment the so-called silent point, where ρ = ρ max /2. At this point the perturbations in the algebraic approach change signature (the hyperbolic equation becomes elliptic for ρ < ρ max /2) [121]. Moreover, it was found that in the dressed metric approach of [80] imposing the initial conditions in the contracting phase and at the bounce lead to the same results [121].
The evolution of the matter density and curvature perturbations in inflation and in the matter bounce scenarios are quite different [122]. Loop Quantum Cosmology models naturally induce an inflationary phase, in which the primeval seed perturbations are created [123]. The LQC evolution generates excited states at the onset of the slow-roll inflation [124]. The power spectrum of density perturbation generated during the induced inflationary phase is broadly scale-invariant [123]. There are two scales that determine the form of the LQC-corrected power spectrum. The first one is associated with the spacetime curvature at the bounce (k I ), while the second one is associated with the spacetime curvature at the onset of inflation (k LQC ). The amplitude of the power spectrum is amplified with respect to the standard predictions of slow-roll inflation for modes k I < k < k LQC [125]. The initial value of the inflaton field (which is a massive scalar field) ψ B , and the mass of the massive scalar field m φ control whether the modes affected by the pre-inflationary dynamics of LQC fall within the window of modes that are observable today. The region that shows non-negligible LQC modifications is a subset which satisfied the condition ψ B ≈ 3 × 10 −6 /m φ [125]. LQC corrections to the power spectra tend to make the tensor spectral index n t more negative, produce a positive running α s of the scalar spectral index, reduce the tensorto-scalar ratio r, and modify the inflationary consistency relation to r/n t = −8 [125]. The particular choice of the initial data for the quantum scalar and tensor perturbations has very little impact on the above conclusions, at least for the reasonable choices of initial vacuum state that is considered above [125].
If we only consider the suitable initial vacuum state that satisfy the requirement that the initial conditions for each mode minimize the time variation of the Mukhanov-Sasaki variable from the bounce to the beginning of the inflation, LQC can generate a power spectrum of the curvature perturbations in which large oscillations with the averaged enhanced power are not present. Instead, in LQC one can obtain for that range of scales a behaviour compatible with the one obtained from the slow-roll inflation formula. Moreover, LQC also predicts a strong power suppression for large scales, a result that provides a better fitting of the current observations [122,126]. The same power suppression effect is predicted to occur in the power spectrum of tensor modes as well [122].
The non-Gaussianity of the primordial curvature perturbations in Loop Quantum Cosmology have been also studied, and it was found that the corrections due to quantum effects are of the same order of magnitude as the inflationary slow-roll parameters in the observable scales. Thus they are well within current observational constraints [124]. The non-Gaussianity in the squeezed limit can be enhanced at superhorizon scales. These effects may yield a large statistical anisotropy on the power spectrum [124].
In the matter bounce scenario, the fluctuation of the curvature and the tensor perturbations become scaleinvariant in a contracting FLRW space-time, where the matter content has vanishing pressure [122]. Then, if a bounce can be generated to provide a non-singular transition from contraction to expansion, these scale-invariant perturbations provide suitable initial conditions for the expanding Universe that can explain the scale-invariance observed in the CMB, under the assumption that the bounce does not modify the power spectrum [122]. LQC can generate the bounce that is required for this scenario to be viable [122]. During the LQC predicted bounce, the tensor-to-scalar ratio will be suppressed. The precise numerical factor of the tensor-to-scalar ratio suppression depends on the dominant matter fields during the bounce. The possible observational dependence of the suppression factor may provide important information about the physics of the bounce [122].
In the present approach of a radiation filled Loop Quantum Cosmological Universe, one of the theoretical possibilities is that the primordial curvature perturbations are generated during the induced inflation period. Hence we may assume that the generation of the primordial curvature perturbation takes place in the postbounce phase. However, the presence of the matter bounce makes also possible the generation of the primordial curvature perturbations in the early contracting phase. During the contraction towards the bounce, the vacuum fluctuations and the curvature perturbations become scale-invariant. Then the existence of the nonsingular bounce predicted by LQC leads to the transition from contraction to expansion, a process that does not modify the scale-invariance of the power spectrum. The scale-invariant power spectrum, maintained after the bounce, provides a suitable initial condition for the expanding Universe, which matches our observations based on the Cosmic Microwave Background Radiation.
If the primordial curvature perturbation are explained by the post-bounce inflationary scenario, the accelerating state before the bounce can be considered as a preinflationary era, whose main role is to generate excited states at the onset of inflation. Hence, in this scenario, the bouncing state corresponds to the starting point of the pre-inflationary phase. The boundary of the preinflation era is characterized by the quantities k LQC and k L , which also describe the power spectrum of density perturbations generated during the inflation. If the generation of primordial curvature perturbation is explained during the contracting phase, the accelerating era before the bounce can be considered as the period when the vacuum fluctuations in curvature and tensor perturbations become scale-invariant. In this scenario, the bouncing point, representing the non-singular transition point from the contracting to the expanding era maintains the scale-invariance of the power spectra of the density perturbations. On the other hand, since the generation of the matter density/curvature perturbations during the inflationary era and during the matter pre-bounce phase predict different power spectra, and different tensor-toscalar ratios, it is possible to test the viability of these scenarios via observations of CMB, as described in [122].
In the present paper we have introduced a theoretical model that combines two major approaches that may lead to a better description of the quantum properties of the gravitational force, Loop Quantum Cosmology and noncommutative space-time theories, respectively. These results may open the possibility of a better understanding of the theoretical phenomena at the Planck scale, and of the in depth comparison of the theoretical predictions with the observational data.