Gravitational waves with generalized holonomy corrections

The cosmological tensor perturbation equation with generalized holonomy corrections is derived in the framework of effective loop quantum gravity. This results in a generalized dispersion relation for gravitational waves, encompassing holonomy corrections. Furthermore, we conduct an examination of the constraint algebra concerning vector modes with generalized holonomy corrections. The requirement of anomaly cancellation for vector modes imposes constraints on the possible functional forms of the generalized holonomy corrections. What's more, we estimate the theoretical value of the effective graviton mass and discuss the potential detectability of this effective mass in future observations.

Cosmological perturbation theory serves an important role in modern cosmology as the nexus between theory and observation, providing a fundamental framework and tool for investigating the large-scale structures of the universe, along with the origin and evolution of the anisotropy of the cosmic microwave background (CMB) radiation.Within the framework of cosmological perturbation theory, the universe's inhomogeneity and anisotropy are treated as fluctuations or perturbations of the spatially homogeneous and isotropic cosmic background.These equations of motion can be systematically solved order by order.Observations of the CMB radiation demonstrate that the universe only slightly deviates from this background on large scales with a magnitude of approximately 10 −5 .When studying such small fluctuations, the perturbation equations up to the linear order can already offer a sufficiently precise depiction of the early universe.The decomposition theorem [1][2][3] states that the linear gravity perturbations can be decomposed into distinct modes, namely scalar, vector and tensor modes.This decomposition is based on the different transformation properties of these modes under spatial rotations.These modes evolve autonomously and remain uncoupled, allowing for independent analysis of distinct perturbation modes.In this paper, we will study the tensor and vector modes within the framework of effective loop quantum gravity (LQG) with generalized holonomy corrections.
In cosmology, delving into the study of gravitational waves (GWs) that emerged from the intricate physical processes of the early universe offers an unparalleled opportunity to peer into the universe's primordial stages.These GWs, often referred to as tensor mode perturbations, carry with them an imprint of the conditions and events that prevailed during the universe's formative stages.Currently, considerable scientific endeavors are dedicated to detecting the subtle signatures left by these tensor mode perturbations on the structure of spacetime.This is achieved by measuring the polarization in the CMB.
Given the potential significance of quantum gravity effects in the universe's early stage, there exists a compelling motivation to delve into the exploration of quantum gravity's possible impact on the propagation of GWs throughout these epochs.The study of such quantum gravity effects on the propagation of GWs holds the promise of shedding light on the fundamental nature of spacetime itself during these early stages.
Among various prospective quantum gravity theories, LQG, a nonperturbative and background-independent quantum gravity theory [4][5][6][7], has recently emerged as a prominent contender and attracted significant attention.The technique of LQG has been effectively employed to achieve quantization within the framework of symmetry-reduced cosmological spacetime, which is commonly known as loop quantum cosmology (LQC) [8][9][10][11][12][13][14].The effective theory of LQC can be formulated by incorporating two types of primary quantum gravitational effects: the inverse volume corrections and the holonomy corrections.These effects can be introduced through both the canonical approach [15][16][17][18][19][20] and the path integral perspectives [21][22][23][24][25][26][27].Significant advancements have been achieved in the applications of LQC within the early universe in recent years.Notably, LQC has successfully replaced the classical big bang singularity from general relativity (GR) with a non-singular big bounce at the effective level, as demonstrated by key works such as Refs.[8][9][10].It can also provide potential origins for the inflationary model [28] and set appropriate initial conditions [29].
Furthermore, quantum gravitational effects may also leave detectable imprints, which can be observed in the CMB [29,30].In summary, LQC is the symmetry-reduced model of LQG, thus inheriting its fundamental characteristics.Within the realm of LQG research, LQC plays a twofold role: firstly, it explains the physical properties of the universe, and secondly, it severs as a natural laboratory to test LQG.
In this paper, we delve into the investigation of LQG effects concerning tensor modes.
Specifically, our focus lies in the exploration of the generalized holonomy corrections on the dynamics of GWs.An exemplary effective holonomy correction involves substituting the classical background connection variable, denoted as k, with the holonomy function sin(γ μk ) , where μ signifies the polymerisation scale and γ is the Barbero-Immirzi parameter [31,32].
In the classical limit where μ → 0, the effective LQG system recovers the classical counterpart.It's noteworthy that this substitution takes place due to our adoption of the minimal spin number j = 1/2 in the SU(2) representation.Owing to the inherent quantization ambiguities, the explicit formulations of holonomy corrections are far from unique.For instance, diverse representations with different j values have been investigated extensively in prior research works [33][34][35][36].Additionally, alternative formulations, incorporating quantization ambiguities or higher powers of μ, have been explored, as evidenced by references such as [17,37].
Recently, the generalized holonomy corrections utilizing the substitution k → g k, p , where p is the conjugate momentum of k and g k, p is a generalized undetermined function, have been put forth in the work of Ref. [38].Through the adoption of the anomaly-free constraint algebra approach or the deformed algebra approach proposed in Refs.[39,40], and by strategically incorporating specific counter terms, they successfully achieve an anomalyfree constraint algebra.Notably, they are able to determine and address all the necessary counter terms within the constraints, thereby elucidating the restriction on this undetermined function g k, p .And then they work out the gauge invariant perturbation equations together with the corresponding inflationary power spectra.In a more recent study [41], the authors employ a comprehensive investigation into the background dynamics of the inflationary phase, considering generalized holonomy corrections.Their research also delves into the impact of these corrections on the primordial scalar power spectrum.It is noteworthy, however, that they employ the classical scalar perturbation equation at the perturbative level, which may lead to a loss of information regarding the generalized holonomy corrections at the level of perturbations.Notably, in a subsequent work [42], the authors explore the primordial power spectra incorporating generalized holonomy corrections based on perturbation equations derived within the framework of deformed algebra presented in Ref. [38].This approach retains information about the generalized holonomy corrections at the perturbative level.
However, the inclusion of the counter terms within the constraints has the potential to introduce quantum gravity corrections that extend beyond the scope of holonomy corrections.It is undoubtedly important to focus solely the effects of holonomy corrections on GW dynamics.Therefore, we will exclude counter terms and consider only holonomy corrections, which may otherwise bring unsuspected corrections.Since we considered a generalized holonomy corrected function, which is undertermined, it is possible to achieve anomaly freedom in perturbative LQG, even though excluding the incorporation of counter terms.As a first step, we also study the anomaly-free constraint algebra of the vector modes in this work.
The organization of the paper is as follows.In Sec.II, we recall the canonical formulation of GR and deduce its classical dynamics.In Sec.III, the quantum-corrected tensor mode equation with generalized holonomy corrections along with the corresponding dispersion relation of GW propagation is derived.In Sec.IV, we consider anomaly cancellation in the constraint algebra for vector modes and obtain the restriction on the function g k, p .Besides, a specific example of the correction function g k, p is introduced to provide a more intuitive and clear interpretation of the quantum dynamics and anomaly cancellation condition obtained earlier.In Sec.V, we estimate the theoretical value of the effective graviton mass and explore the potential detectability of this mass in future observations.In Sec.VI, we summarize our results and conclusions, and then present our outlook.Throughout this paper, we adopt the Planck unit system, where c = G = ℏ = 1, with c denoting the speed of light, G representing the gravitational constant, and ℏ signifying the reduced Planck constant.

II. CANONICAL FORMULATION AND CLASSICAL DYNAMICS
In this section, following the idea in Ref. [43], we will review the canonical formulation of GR and deduce its classical dynamics, which will help us derive the quantum dynamics with generalized holonomy corrections in the section that follows.

A. Spacetime metric in canonical formulation
A spatially flat Friedmann-Robertson-Walker (FRW) background spacetime can be described by the following metric: where t denotes cosmic time, i.e. the proper time of isotropic observers, η represents conformal time, and a (t) or a (η) is the scale factor.The indices a, b, . . .denote spatial coordinates.
When only taking into account the perturbations of the tensor modes, the metric can be expressed as follows: In this equation, the symmetric metric perturbation field h ab adheres to the conditions of being transverse and traceless, i.e. ∂ a h ab = δ ab h ab = 0, which effectively removes any vectorial or scalar contributions.In a canonical formulation based on the 3+1 decomposition of spacetime, the metric can be expressed in terms of the spatial metric q ab , the lapse function N , and the shift vector N a as follows: This canonical formulation serves as a powerful framework for effectively depicting the dynamics of spacetime.
In this formulation, the spatial metric q ab is replaced by the densitized triad E a i , as defined by: Here, e a i and e i a represent the triad and its dual, respectively, satisfying q ab = e i a e i b .The indices i, j, . . .correspond to internal indices.
Another important basic variable is the Ashtekar-Barbero connection: Here, K i a represents the extrinsic curvature, and the Barbero-Immirzi parameter γ is a positive value, i.e., γ > 0 [31,32].The spin connection Γ i a , which corresponds to the covariant derivative D a satisfying D a e a i = 0 by definition, is determined as follows: .
The Ashtekar-Barbero connection A i a is the canonically conjugate variable of the densitized triad E a i and characterizes spatial curvature.After perturbing the spacetime, for any physical variable V , the background part and the perturbation part are respectively represented by V and δV , that is, V = V + δV .By matching Eqs. ( 1) and ( 3) and utilizing relevant definitions, we can derive the background variables as follows: where p = a 2 .
By comparing Eqs. ( 2) and (3), and leveraging the relevant definitions and background variables in Eq. ( 7), the perturbative variables are computed at the linear level to yield: Regarding the extrinsic curvature K i a , determining the explicit expressions for both the background component k and the perturbation component δK i a necessitates solving the Hamiltonian canonical equations, a process that will be elaborated on later.
Within the context of linear perturbed theory, the independent phase space variables are k, p and (δK i a , δE a i ), whereas the symplectic structure is given by the following non-trivial Poisson brackets [39,40,43,[45][46][47][48]: Here V 0 := Σ d 3 x is a fiducial volume introduced to ensure a finite symplectic structure for the background variables.This is accomplished by confining the integration of the action to a finite cell, rather than extending it across the entirety of space R 3 .The restriction to a cell doesn't result in information loss due to the background's homogeneity.It's important to note that this fiducial quantity only appears in the basic variables and their symplectic structure, without influencing the ultimate physical outcomes.Consequently, there are separate canonical structures for the background and the perturbations, yet these variables will be dynamically coupled.The dynamical homogeneous background, in particular, would receive back-reaction effects beyond the linear order.

C. Classical dynamics
Since GR is a fully constrained theory, its Hamiltonian is formulated as a summation of smeared constraints.In the context of a canonical triad formulation, these constraints encompass three categories: the Hamiltonian constraint, the diffeomorphism constraint, and the Gauss constraint.However, when one focuses on linear perturbations involving solely the tensor mode, only the Hamiltonian constraint remains [43].
In the framework of the connection dynamical formalism, the Hamiltonian constraint is expressed as: By utilizing the expressions of the perturbed basic variables and the curvature , one can perturb the Hamiltonian constraint up to the quadratic terms in perturbations, yielding: It should be mentioned here that the terms dependent on γ are eliminated from the Hamiltonian constraint when employing the spin connection and taking into account the symmetry of the densitized triad and extrinsic curvature for the tensor mode.
In the standard covariant formulation, the equations of motion can be derived by taking the variation of the action with respect to the metric.In a canonical formulation, however, the equations of motion governing any phase space function f can be elegantly expressed as In this context, the dot signifies the derivative with respect to the coordinate time x 0 , which depends on the chosen background lapse function N .
Subsequently, we will compute the Hamiltonian canonical equations for both the background variables and the perturbed variables, respectively.Using the expression of the Hamiltonian constraint (11), we yield the following equations of motion for the background variables: Since the matter Hamiltonian constraint H M [N ] is independent of the extrinsic curvature, Eq. ( 13) remains unaffected by the matter Hamiltonian constraint, in contrast to Eq. ( 12).
When choosing N = √ p, where the coordinate time x 0 corresponds to the conformal time η, according to Eq. ( 13) and the relation p = a 2 , one can obtain the background extrinsic curvature k = ṗ 2p = ȧ a =: H , which is nothing more than the conformal Hubble parameter.We proceed to calculate the Hamiltonian canonical equation for the perturbed densitized triad δE a i , which is derived as: While deriving the aforementioned equation, we have employed the expression of δE a i as provided in Eq. (8).By combining Eq. ( 14) with the time derivative of the expression of δE a i in Eq. ( 8) and using Eq. ( 13), the expression for the perturbed extrinsic curvature can be formulated as follows: Then, using the δE a i expression from Eq. ( 8) along with Eq. ( 15), we derive the Hamiltonian canonical equation for the perturbed extrinsic curvature δK i a as follows: Moreover, by amalgamating Eq. ( 16) with the time derivative of Eq. ( 15) and incorporating Eqs. ( 12) and ( 13), one can deduce the equation governing the motion of tensor mode perturbations: where The quantity Π i a characterizes linear source terms that are both transverse and traceless, and these terms can be associated with the transverse and traceless component of the perturbed stress-energy tensor through the relation Π i a = pδT (TT) i a .For an in-depth discussion, please refer to Ref. [43].
Furthermore, the choice of N = √ p results in Eq. ( 17) transforming into: Alternatively, setting N = 1 to correspond with cosmic time t yields Eq. ( 17) taking the form of: 1 2 When the source terms Π i a are absent, the aforementioned tensor mode equation demonstrates propagating wave solutions recognized as GWs in the given cosmological background.
For the investigation of wave propagation, it is often advantageous to deduce the associated dispersion relation, illustrating the interrelation between the frequency ω and the proper wave vector ⃗ k/a.Notice that here the frequency ω corresponds to cosmic time t, and accordingly, the tensor perturbation equation of the version in Eq. ( 20) is to be applied.As our main interest is on local propagation taking no account of cosmic scales, we ignore the friction term, which is the first-order time derivative term originating from cosmological expansion and proportional to k and, consequently, to the Hubble parameter.
By employing a plane wave ansatz h i a ∝ exp iωt − i ⃗ k • ⃗ x , one arrives at: which is the well-known classical dispersion relation.This relation gives rise to the group velocity of GWs as follows: As expected, the group velocity of classical GWs is exactly equal to the speed of light.

III. QUANTUM DYNAMICS WITH GENERALIZED HOLONOMY CORREC-TIONS A. Corrected tensor mode equation
So far, the explicit form of quantum corrections in LQC lacks uniqueness and subjects to quantization ambiguities.Building upon the idea described in Ref. [38], in order to introduce holonomy corrections, we introduce a generalized undefined function g k, p .In the classical limit, g k, p is required to satisfy g k, p → k, thereby ensuring that the standard GR is recovered.After replacing k → g k, p , the classical Hamiltonian constraint (11) becomes: We would like to emphasize that, in our approach, the background component k is substituted by the same quantum correction function g k, p , regardless of whether it pertains to background or perturbation terms.This is the most reasonable consideration.However, when we consider the conventional holonomy corrections as k → sin(γ μk ) γ μ , we encounter issues with achieving an anomaly-free constraint algebra, even when applying it exclusively to vector modes, if the parameter m assumes identical values in both background and perturbation terms [43,46].In the subsequent section, we will demonstrate that by employing the aforementioned generalized holonomy corrections, we can achieve an anomaly-free constraint algebra in the context of vector modes.
Following the same procedure as in the classical case discussed in subsection II C, we will derive the tensor mode equation with generalized holonomy corrections based on the quantum-corrected Hamiltonian constraint (23).Firstly, we calculate the Hamiltonian canonical equations of k, p and δE a i , which are given as follows: It is necessary to mention that the expression of δE a i in Eq. ( 8) still remains applicable in the present scenario with quantum corrections.By combining the results from the aforementioned Hamiltonian canonical equations and the time derivative of the expression of δE a i in Eq. ( 8), we arrive at the modified formulation for perturbed extrinsic curvature: Then the Hamiltonian canonical equation of δK i a can be computed as: By performing a time derivative on Eq. ( 27) and then comparing it with Eq. ( 28), we can derive the tensor mode equation incorporating generalized holonomy corrections: where and As expected, both the friction term and the source terms arising from the matter Hamiltonian receive the generalized holonomy corrections.Moreover, an additional term emerges, proportionally related to the field perturbation h i a .It's worth highlighting that the nature of these corrections is intricately linked to the specific functional form of the correction function g k, p .
In conclusion, a generalized correction function g k, p has the potential to bridge the theoretical framework with empirical observations.Furthermore, the emergence of the extra term, linked to the field perturbation h i a , adds a layer of complexity to the corrected equations.This term's significance and implications warrant a comprehensive analysis.Its existence might entail novel physical interpretations or provide insights into previously uncharted aspects of the phenomena being studied.

B. Modified dispersion relation
In this subsection, we proceed to study the modified dispersion relation of GWs resulting from the inclusion of generalized holonomy corrections.Only the additional term proportionally related to the field perturbation h i a receives the holonomy corrections once the source terms are removed and the friction term is ignored.By performing the computations of the classical scenario but employing the holonomy corrected wave equation (Eq.( 29)), we arrive at the dispersion relation with generalized holonomy corrections: where The modified dispersion relation (Eq.( 32)) reveals that under quantum corrections, generalized holonomy corrections introduce a novel additive term m G 2 within the dispersion relation, as compared to the classical scenario described by equation ( 21).This signifies that GWs subject to these generalized holonomy corrections embody an effective mass m G .
Remarkably, this modified dispersion relation (32), regardless of the specific form of correction function g k, p , is formally consistent with that with the conventional holonomy corrections in Ref. [43].However, the entire impact arising from distinct specific choices of the correction function g k, p becomes encapsulated within the effective mass squared term For visualization purposes, we present the modified dispersion relation with different m G 2 in Fig. 1.
FIG. 1: Dispersion relation for GWs with generalized holonomy corrections.
The group velocity of GWs corresponding to the modified dispersion relation (32) becomes: This indicates that different modes of GWs propagate with distinct group velocities that remain below unity.As a result, causality is maintained, provided that the condition of positivity holds for the effective mass squared term m G 2 .
Obviously, the dependence of the group velocities v Q GW on the effective mass m G of GWs is evident, contingent upon the specific form of the generalized holonomy-correction function g k, p .When considering a fixed proper wave number k/a, the behavior of v Q GW becomes apparent: it diminishes as m G increases.This behavior aligns with expectations, as greater effective mass m G corresponds to more pronounced quantum corrections received by the GWs.Upon closer examination of GWs with identical effective mass m G , a notable observation arises: distinct k/a modes of GWs exhibit varying group velocities v Q GW , consistently remaining below unity.Notably, the group velocities of lower k/a modes deviate significantly from the classical case, while those with higher k/a values gradually approach unity-the speed of light.This finding indicates that, for modes featuring large proper wave numbers k/a, the impact of holonomy corrections becomes negligible.

IV. ANOMALY FREEDOM OF THE VECTOR MODES
We've successfully derived the GW equation with generalized holonomy corrections and discussed its associated dispersion relation.Moving forward, our focus will shift to addressing the constraint algebra for vector modes, along with investigating a particular instance of the generalized holonomy-correction function g k, p .Specifically, we will delve into the anomaly cancellation condition of the constraint algebra for vector modes and the parameter restrictions on this specific form of g k, p .

A. Anomaly cancellation in the constraint algebra for vector modes
In the classical framework of GR, the Poisson brackets among the smeared constraints weakly vanish on the constraint surfaces, thereby forming a first-class constraint algebra.However, when quantum corrections are factored in, to ensure a consistent theory, it becomes imperative to maintain closedness within this algebraic structure.In the context of background dynamics alone, it is relatively straightforward to observe that the quantumcorrected constraints continue to yield a closed algebra.Nevertheless, the scenario evolves as we include perturbations into the constraints.In this case, the resultant constraint algebra often deviates from closure, featuring extra terms referred to as anomaly terms.Achieving a well-defined, closed algebra necessitates the vanish of these anomaly terms.
Regarding the tensor mode with generalized holonomy corrections under consideration, as mentioned previously, only the Hamiltonian constraint remains.Its effective form H Q G [N ] has been presented in Eq. ( 23).It's evident that the Poisson bracket between two effective ] is identically zero due to the absence of lapse perturbations δN , aligning with classical results.Thus, the constraint algebra for the tensor mode are automatically anomaly-free.
Let us now shift our focus to examining the constraint algebra for vector modes.The gravitational contribution of the perturbed Hamiltonian constraint, considering only vector modes and accounting for generalized holonomy corrections, is expressed as follows: On the other hand, the perturbed gravitational diffeomorphism constraint, considered exempt from quantum corrections due to its direct quantization through phase space transformations [49], is formulated as follows: Once again, it's straightforward to observe that the Poisson bracket ] is trivial.However, the non-trivial Poisson bracket between H Q vG [N ] and D vG [N a ] turns out to be: where The anomaly cancellation condition A c j = 0 imposes constraints on the functional form of g k, p , leading to a first-order quadratic partial differential equation: Because of being a nonlinear partial differential equation, obtaining an explicit analytical general solution for Eq. ( 39) remains challenging.Subsequently, our focus will shift towards examining a specific instance of the generalized holonomy-correction function g k, p .

B. Parameter restrictions for a specific correction function
In Ref. [41], the authors propose a particular expression for g k, p , given by: where μ (p) = ∆/p following the standard μ-scheme.Here, ∆ = 4 √ 3πγl P 2 , where l P represents the Planck length-the fundamental length unit in the Planck unit system, in which it equals unity.This quantity serves as the LQG area gap, signifying the minimum non-zero eigenvalue of the quantum area operator.The function f takes the form: Here, both α and C are undetermined parameters, with the specific note that C is a nonnegative integer.The independent variable of this function f is only considered to range over [0, π] as indicated in Ref. [41].
The provided function form g k, p satisfies the prerequisite that g k, p → k in the classical limit (μ → 0).Furthermore, it's evident that by setting the parameter C to zero, Eq. ( 40) restores the standard μ-scheme of the holonomy corrections, i.e., g k, p = sin(γ μk ) γ μ .
Additional details and discussions regarding this specific form of g k, p can be found in Ref. [41].This function form of g k, p can serve as a simplified toy model to provide greater clarity and intuition for the quantum-corrected dynamics explored in Sec.III, as well as for comprehending the constraints on the function form of g k, p arising from the anomaly cancellation of vector modes, as discussed in the previous subsection.Moving forward, we will delve into the analysis of constraints on the free parameters α and C in Eq. ( 41) based on the previously obtained results.
We would like to underscore that, despite adopting the specific form of the function f as outlined in Eq. ( 41), obtaining an explicit analytical general solution for the anomaly cancellation condition (39) remains challenging.Fortunately, our main focus is on the propagation of GWs within the low-curvature phase, which is far away from the high-curvature stage of the very early universe.In this context, the background extrinsic curvature k (and therefore the Hubble parameter) is considered to be much lower than O (1).This means that the argument γ μk of the function f in Eq. ( 40 Substituting Eq. ( 40) and Eq. ( 41) into Eq.( 39) resulting from the anomaly cancellation of vector modes, gives rise to the following equation: This implies that the anomaly cancellation, up to the order of k3 , enforces a stringent constraint on the parameter α, namely, α = 1.The choice of α = 1 corresponds precisely to the situation maintaining the asymptotic behavior of LQC, as described in Ref. [41].
Furthermore, the modified dispersion relation (Eq.( 32)) incorporates an effective mass squared term m G 2 defined in Eq. ( 33) with Eq. ( 30).Notably, this m G 2 term, governed by the yet undetermined function g k, p , lacks a guarantee of strict positivity.The stipulation for m G 2 to be positive indeed imparts an additional constraint on the form of g k, p .Upon substituting Eq. ( 40) and Eq. ( 41), we arrive at the expression: where the second equality is obtained by imposing the condition α = 1.This indicates that m G 2 is bound to be positive on the premise of setting α as unity.In other words, the requirement of a positive m G 2 term does not introduce any supplementary constraints on the free parameters, thereby leaving the parameter C unconstrained within our considerations.
Moreover, it is evident that m G 2 exhibits a linear increase with the parameter C, up to the order of k4 .

V. CONSTRAINTS ON THE EFFECTIVE MASS OF THE GRAVITON
In this section, we will further delve into a detailed exploration of the effective mass of the graviton when considering the specific form (40) with Eq. ( 41) of the generalized holonomy corrections.Our emphasis will be on scrutinizing the constraints imposed on the effective mass through astronomical observations.
As a preliminary, we first derive the modified Friedman equation with the generalized holonomy corrections.To achieve this, we express the effective total Hamiltonian constraint for the homogeneous and isotropic background in the following manner: where the function f given in Eq. ( 41) now has been applied with the restriction condition α = 1 from anomaly cancellation and becomes: Then varying this Hamiltonian constraint with respect to the background lapse N gives rise to the Hamiltonian constraint equation, i.e. the effective Friedmann equation: where ρ is the background energy density defined as ρ := . The above effective Friedmann equation can be further reduced to: In the given equation, we introduce a constant ρ c := 3 8πγ 2 ∆ , representing the critical energy density at which the bounce happens in conventional holonomy-correction LQC [9][10][11].
From Eq. ( 48), it becomes evident that the energy density ρ exhibits a new finite upper critical value of ρ c / (C + 1), decreasing as C increases.This contrasts with the conventional holonomy-correction scenario, where the upper bound of the energy density ρ remains the definite constant ρ c .For more discussions on this effective Friedmann equation with generalized holonomy corrections, readers can also refer to Ref. [41].Now, let's delve into the discussion of the effective mass m G associated with tensor mode perturbations.By combining Eq. ( 43) with Eq. ( 47), we derive the following expression: The equation above indicates that the effective mass m G depends on the background energy density ρ, which evolves as the universe expands.Consequently, m G evolves with time.
Nevertheless, this temporal evolution can be ignored because we only focus on the lowcurvature phase.We would like to emphasize once more that m G 2 exhibits nearly linear growth with respect to the parameter C. We would like to emphasize that, as we all known, there is currently a lack of comprehensive analysis on the quantum gravity effects in scalar-induced gravitational waves, even only incorporating the usual holonomy corrections.In addition, a complete analysis of the angular power spectra is also absent.We aspire to address this gap and complete the puzzle.This effort represents a unique opportunity to bridge the gap between theoretical quantum gravity frameworks and observable cosmological phenomena.Work on these endeavors is currently in progress.

Contents I. Introduction 2 II. Canonical formulation and classical dynamics 6 A. Spacetime metric in canonical formulation 6 B. Canonical variables in Ashtekar's formulation 7 C. Classical dynamics 8 III.
Quantum dynamics with generalized holonomy corrections A. Corrected tensor mode equation B. Modified dispersion relation IV.Anomaly freedom of the vector modes A. Anomaly cancellation in the constraint algebra for vector modes B. Parameter restrictions for a specific correction function V. Constraints on the effective mass of the graviton VI.Summary and outlook I. INTRODUCTION wherein H [N ] is the total Hamiltonian constraint, encompassing both the gravitational component H G [N ] and the matter component H M [N ].When f represents the fundamental variables, its equations of motion are just the Hamiltonian canonical equations.
Then, we can estimate the value of the effective mass m G of the graviton at the present epoch according to Eq.(49).Given the orders of magnitude of γ ∼ O (1), ∆ ∼ O (1),ρ c ∼ O (1)and the energy density of the present universe ρ(t 0 ) ∼ O (10 −120 ), one obtains that m G (t 0 ) ∼ O √ C O (10 −120 ), which corresponds to m G (t 0 ) ∼ O √ C O (10 −92 ) eV/c 2 after filling in the physical constants c, G and ℏ.Latest observational data of pulsar timing arrays (PTAs) and their analysis have established the state-of-the-art upper bound on the graviton mass of m G < O (10 −24 ) eV/c 2 [50, 51].And this upper bound could continue to be lowered in future gravitational wave measurements.Consequently, the parameter C obtains a conservative constraint of C < O (10 136 ).This constraint appears to be too loose, and we believe it is mainly limited by the accuracy of the observational upper bound of the graviton mass.Unfortunately, due to the extremely large order of magnitude of the upper bound on C obtained at present, it is unlikely to obtain an adequately tight restriction by observations in the near future.On the other hand, provided that the parameter C ∼ O (1), our theoretically estimated effective graviton mass m G (t 0 ) will be far below the observational upper bound currently.Nonetheless, considering the dependence of the effective graviton mass m G on the background energy density ρ, this effective mass might play a significant role in some physical phenomena in the early universe, such as inflation.