Eddington-inspired-Born-Infeld tensorial instabilities neutralized in a quantum approach

The recent direct detection of gravitational waves has highlighted the huge importance of the tensorial modes in any extended gravitational theory. One of the most appealing approaches to extend gravity beyond general relativity is the Eddington-inspired-Born-Infeld gravity which is formulated within the Palatini approach. This theory can avoid the Big Bang singularity in the physical metric although a Big Bang intrinsic to the affine connection is still there, which in addition couples to the tensorial sector and might jeopardize the viability of the model. In this paper, we suggest that a quantum treatment of the affine connection, or equivalently of its compatible metric, is able to rescue the model. We carry out such an analysis and conclude that from a quantum point of view such a Big Bang is unharmful. We expect therefore that the induced tensorial instability, caused by the Big Bang intrinsic to the affine connection, can be neutralized at the quantum level.

(Dated: January 16,2020) The recent direct detection of gravitational waves has highlighted the huge importance of the tensorial modes in any extended gravitational theory. One of the most appealing approaches to extend gravity beyond general relativity is the Eddington-inspired-Born-Infeld gravity which is formulated within the Palatini approach. This theory can avoid the big bang singularity in the physical metric although a big bang intrinsic to the affine connection is still there, which in addition couples to the tensorial sector and might jeopardize the viability of the model. In this paper, we suggest that a quantum treatment of the affine connection, or equivalently of its compatible metric, is able to rescue the model. We carry out such an analysis and conclude that from a quantum point of view such a big bang is unharmful. We expect therefore that the induced tensorial instability, caused by the big bang intrinsic to the affine connection, can be neutralized at the quantum level. a imanol@ubi.pt b mariam.bouhmadi@ehu.eus c b97202056@gmail.com d pisinchen@phys.ntu.edu.tw

I. INTRODUCTION
It is commonly recognized that Einstein's general relativity (GR), though very successful in describing our universe, nonetheless suffers from several fundamental puzzles. On the very early stage of the universe where the energy scale and the curvature scale were huge, say, close to the Planck scale, a purely classical description of gravity based on GR would not be sufficient. Actually, it is expected that a fundamental quantum theory of gravity is necessary such that some pathologies of GR at high energy scales can be resolved, such as the non-renormalizability of the theory and the issue regarding spacetime singularities. Whereas, it is still not clear so far how such a fundamental quantum gravity theory should be built in a self-consistent way. The development of a complete quantum theory of gravity is still an open question and it is certainly one of the most active research directions in modern theoretical physics.
From a more conservative point of view, to escape from the aforementioned theoretical swamp, one may resort to other modified theories beyond GR and regard them as effective theories of the unknown quantum theory of gravity [1]. It is likely that such extended theories of gravity, even presumably not complete, are already able to ameliorate the UV problems in GR. Among the plethora of extended theories of gravity, the Eddington-inspired-Born-Infeld gravity (EiBI) proposed in [2] is appealing in several theoretical aspects. First, it reduces exactly to GR in vacuum and deviates from it when matter fields are included. Second, due to the square root structure in the gravitational action, the curvature scale and the energy scale seem to be bounded from above and the big bang singularity is naturally avoided in the EiBI gravity. Third, the theory is simple in the sense that it only contains one free additional parameter, the Born-Infeld constant κ compared with GR. Fourth, it is free of ghost instabilities because the theory is constructed through the Palatini rather than the metric variational principle. Actually, the idea of including the Born-Infeld structure into the gravitational theory was proposed in Ref. [3]. However, the theory is built with the metric variational principle and it has ghost because of the higher order derivative terms in the field equations. The EiBI theory, on the contrary, is formulated via the Palatini variational principle. The field equations only contain second order derivatives and therefore no ghost is present in the theory. The applications and several properties of the EiBI gravity have been studied widely in the literature . Some attempts to quantize the EiBI gravity have been proposed in Refs. [40][41][42][43][44][45]. See also [46] for a nice review on the EiBI gravity and other interesting Born-Infeld inspired theories of gravity. A further motivation to consider the EiBI theory is that the Born-Infeld type of theories, of which EiBI gravity is a subclass, have intrinsic Noether symmetries as shown in [47]. This is not surprising as the same happens in other modified theories of gravity [48]. For an interesting review on Noether symmetry, please see [49].
In this paper, we will highlight an important issue regarding the viability of the EiBI theory. In the EiBI gravity, the big bang singularity in the physical metric is avoided by hiding the divergences of quantities in the second spacetime structure defined by the affine connection. The physical metric g µν is non-singular while the other metric, which we will call the auxiliary metric q µν later, turns out to be singular. Since the matter field is assumed to be coupled only to the physical metric, the singularity in the auxiliary metric seems to be unharmful for a physical observer. However, it has been proven in Refs. [8,16] that the metric perturbations, especially the tensor perturbations, are actually unstable for the non-singular solutions in the EiBI gravity, jeopardizing the validity of the theory. A more careful analysis in Ref. [30] reveals that the propagation of gravitational waves does see the structure of the auxiliary metric and it is the singularity in the auxiliary metric that gives rise to the linear instabilities of the theory. It should be noted that the problem of tensor instabilities mentioned above can be ameliorated for a positive Born-Infeld coupling constant if a time-dependent equation of state parameter is considered [9].
In order to resolve this problem, we will suggest a quantum treatment to the EiBI gravity in the framework of quantum geometrodynamics. In this approach, the construction of the Wheeler DeWitt (WDW) equation is crucial since the WDW equation describes the quantum evolution of the universe as a whole [50]. The derivation of the WDW equation stems from a self-consistent classical Hamiltonian, which and all the phase space functions are then promoted to quantum operators. The Hamiltonian, being a first class constraint of the system, turns out to be a restriction on the Hilbert space which is exactly the WDW equation. The strategy is to see whether the wave function would vanish near the configuration of the singularity in the auxiliary metric, satisfying the DeWitt (DW) boundary condition [51]. If the answer is yes, it is then expected that the singularity can be avoided in the quantum world and the linear instabilities, which result from this singularity, can be naturally resolved. For the sake of completeness, we will consider two kinds of matter descriptions, one is the perfect fluid description and the other is the scalar field description. For the perfect fluid description, the matter field is governed by a perfect fluid with a constant equation of state. The system contains only one degree of freedom corresponding to the scale factor of the metric. As for the scalar field description, a scalar field degree of freedom is included into the system and the WDW equation turns out to be a partial differential equation with two independent variables. For each description and each non-singular solution, we will solve the corresponding WDW equation and we will exhibit that the DW condition can always be satisfied, indicating the resolution of the singularity in the auxiliary metric as well as the instabilities via quantum effects.  TABLE I. This table summarizes how the EiBI theory of gravity cures the big bang singularity in a radiation dominated universe. If κ > 0, the big bang singularity in the physical metric is replaced with a loitering stage in which the universe gets its minimum size in the infinite past. If κ < 0, the physical metric is described by a bouncing scenario in the past. However, there is still a big bang singularity in the auxiliary metric for both cases.
This paper is outlined as follows. In section II, we briefly review the non-singular cosmological solutions in the EiBI gravity, depending on different signs of the Born-Infeld parameter κ. We will also exhibit how the tensor instabilities are related to the singularity of the auxiliary metric. In section III, we use the perfect fluid description and derive the WDW equations of the non-singular solutions with regard the physical metric for each sign of κ. The WDW equations within the scalar field description are obtained in section IV. After deriving the WDW equation, we will obtain the wave function and see whether the solution satisfies the DW boundary condition near the singularity of the auxiliary metric. For the perfect fluid description and the scalar field description, the WDW equations will be solved, respectively, in sections V and VI. We finally conclude in section VII.

II. THE CLASSICAL UNIVERSE: BIG BANG IN THE AUXILIARY METRIC
The EiBI gravity is formulated by the following action [2]: where S M is the matter Lagrangian of the matter field Ψ and it is assumed to be coupled only to the physical metric g µν . There is only the symmetric part of the Ricci curvature tensor R (µν) appearing in the action and the curvature is constructed by the affine connection Γ, which is in principle independent of g µν . Based on these assumptions, the theory respects the projective symmetry and one does not have to consider torsion fields since they can be removed by imposing a gauge fixing condition [46]. Moreover, in the action (2.1), |g µν + κR (µν) | is the absolute value of the determinant of the tensor g µν + κR (µν) , where κ characterizes the theory and has inverse dimensions to that of a cosmological constant. 1 Finally, the dimensionless constant λ qnantifies the effective cosmological constant at the low curvature limit. Since the theory is formulated within the Palatini variational principle, one has to vary the action with respect to the physical metric and the affine connection separately. It turns out that one can define an auxiliary metric satisfying λq µν = g µν + κR (µν) such that q µν is compatible with the affine connection. One of the field equations relates algebraically the matter field to the two metrics, and the other equation is a second order differential equation of q µν . It can be seen that when the curvature vanishes, the two metrics are identical up to a constant conformal rescaling, rendering the equivalence of the EiBI theory and GR in the zero curvature regime.
It is well known that the EiBI gravity reduces to Einstein GR when matter fields are absent. However, the theory could have significant differences from GR when, say, the curvature and the energy density of the matter field take large values. Essentially, that is how the big bang singularity is removed in the EiBI gravity. The existence of the affine structure and its corresponding auxiliary metric actually plays a crucial role in the avoidance of singularities. The divergences of the physical metric at the singularity are transferred to the auxiliary spacetime, leaving the physical metric g µν non-singular. Since the matter field only sees the spacetime structure of the physical metric, the hidden singularity in the auxiliary metric seems unharmful for physical observers. Depending on the sign of the parameter κ, the big bang singularity can be replaced with a bouncing solution in the physical metric when κ < 0, or can be replaced with a loitering stage in which the universe acquires its minimum size in the asymptotic past when κ > 0. Table I briefly summarizes how the EiBI gravity cures the big bang singularity in the physical metric, and also points out the singularity appearing in the auxiliary metric.
In the following subsections, we will briefly review how the big bang singularity is avoided in the EiBI gravity with different signs of κ, and we shall point out the fact that the big bang singularity actually migrates to the auxiliary spacetime, i.e., the curvature invariants defined by the auxiliary metric diverge and the scale factor of the auxiliary metric is zero. We will illustrate it by considering a homogeneous and isotropic universe filled with a perfect fluid with a constant and positive equation of state w > 0. Then, we will mention how the instability issues in the physical metric arise alongside the auxiliary singularity, which motivates us for the quantum analysis in this paper.
A. The big bang singularity in the auxiliary metric with κ < 0 We first consider a homogeneous and isotropic universe which can be described by the following metric ansatz: where N (t) and a(t) are the lapse function and the scale factor of the physical metric g µν , while M (t) and b(t) represent the lapse function and the scale factor of the auxiliary metric q µν . These functions are functions of the cosmic time t and they can be determined by the Euler-Lagrange equations of motion. For later convenience, we will define two new variables If the universe is dominated by a perfect fluid with energy density ρ and pressure p = wρ, the equations of motion obtained from varying the action with respect to g µν relate algebraically the two metrics to the matter sector as follows From the above equations (2.4), one can see that when κ < 0 the energy density of the perfect fluid ρ is bounded from above by Since the energy density ρ ∝ a −3(1+w) is bounded from above, the physical scale factor has a minimum value a m1 satisfying ρ (a m1 ) = λ/|κ|. The Hubble rate defined by the physical metric H ≡ȧ/a, where the dot denotes the derivative with respect to t, reads [12] when a → a m1 . By assuming a constant lapse function N , it can be proven that the big bang singularity in the physical metric is replaced with a bouncing solution in the sense that Eq. (2.6) can be integrated to get a − a m1 ∝ t 2 .
To study the behavior of the auxiliary metric when a → a m1 , we rewrite Eqs. (2.4) as follows When ρ → λ/|κ| and a → a m1 , Eqs. (2.7) can be written as Therefore, at the bounce where a = a m1 , the auxiliary scale factor b vanishes. On the other hand, the scalar curvature defined by the auxiliary metric is given by When a = a m1 , it can be seen that R[q] diverges because Y → ∞. Also, by suitably choosing the lapse functions, it can be shown that this divergence happens at a finite time t. Therefore, there is a big bang singularity in the auxiliary metric when b = 0.
B. The big bang singularity in the auxiliary metric with κ > 0 If κ is positive, the big bang singularity in the physical metric is again avoided in the EiBI gravity but in a different manner. In this case, according to Eqs. (2.4), the energy density of the perfect fluid ρ is bounded from above by Therefore, the physical scale factor has a minimum value a m2 satisfying ρ (a m2 ) = λ/ (wκ). The Hubble rate defined by the physical metric can be approximated as [12] when a → a m2 . By assuming a constant lapse function N , Eq. (2.11) can be integrated to get a − a m2 ∝ e t . It can be seen that the scale factor a takes its minimum value when t → −∞. The big bang singularity in the physical metric is thus avoided.
Regarding the asymptotic behavior of the auxiliary metric when a → a m2 , we use Eqs. (2.7) and consider the limit Therefore, the auxiliary scale factor b vanishes and it can be shown that the auxiliary curvature diverges when b → 0. Also, by suitably choosing the lapse functions, it can be proven that this divergence happens at a finite time t. Therefore, there is a big bang singularity in the auxiliary metric.

C. The instability of linear perturbations
In the EiBI gravity, the big bang singularity in the physical metric can be avoided in the sense that the matter field is minimally coupled with the physical metric, hence the physical observers can only see the geometry of that metric, which is free of the big bang singularity. However, non-singular behaviors of the physical metric in the EiBI gravity are still problematic because of the tensor instabilities. Actually, it has been proven in Ref. [30] that these tensor instabilities are highly related to the singular behaviors of the auxiliary metric. In other words, the propagation of gravitational waves would be affected by the geometry of the auxiliary metric. The tensor instabilities in the EiBI gravity were firstly found in Ref. [8]. In addition, the instabilities of scalar mode and vector mode perturbations have been discovered in Ref. [16]. In this subsection, we will briefly review the tensor instabilities of the non-singular solutions in the EiBI gravity and it will become clear that these instabilities are indeed related to the singularity in the auxiliary metric.
Considering the tensor perturbations of the metrics such that the perturbed metrics are δg ij = a 2 h ij and δq ij = b 2 γ ij , it has been proven in Ref. [52] that in the absence of any anisotropic stress, the transverse-traceless tensor perturbations of the two metrics are equivalent, that is, h ij = γ ij . The evolution of the tensor perturbation is described by the following equation [8,16,30,32]: It can be seen that the propagation of the tensor mode is able to see the geometry of the auxiliary metric. When κ < 0, by using Eqs. (2.6) and (2.8), it can be proven that the coefficient of the friction term in Eq. (2.13), that is, the The coefficient of the last term, i.e., M 2 /b 2 , is proportional to 1/b 4 in the same limit. Therefore, these two terms diverge and the tensor perturbation is severely unstable at the physical bounce. On the other hand, for a positive κ, the coefficients of both the friction term and the last term are proportional to b 4 when b → 0. In this regard, the tensor mode behaves linearly in time. Therefore, when approaching the loitering stage at t → −∞, the tensor mode grows linearly backward in time and it diverges in the asymptotic past.
As can be seen above and according to the results in Ref. [30], the instabilities of the tensor perturbations in the EiBI gravity indeed result from the divergence of the auxiliary metric. Even though the physical observers can only see the non-singular metric, the linear instabilities still jeopardize the validity of the theory. This motivates us to study whether the hidden singularity in the auxiliary metric can be resolved by including some sorts of quantum effects and we will address this issue in the following sections.

III. THE WDW EQUATION: PERFECT FLUID
As mentioned in the previous section, the instability of linear perturbations in the EiBI gravity is highly associated with the divergences appearing in the auxiliary metric. Therefore, as long as such a singularity can be ameliorated by quantum effects, the instability problems can be naturally resolved. To address this issue, we shall consider a quantum geometricodynamical approach in which the WDW equation plays a central role. To derive the WDW equation, one is supposed to start with the correct classical Hamiltonian H T , which gives the classical equations of motion, and then promote the Hamiltonian to a quantum operator: H T →Ĥ T . In this regard, it can be proven that the Hamiltonian stands for a first class constraint, indicating that it corresponds to a restriction on the Hilbert space, more precisely, In Refs. [6,46], it was shown that the EiBI action (2.1) can be transformed to its Einstein frame via a Legendre transformation. After such a transformation, the new action reads On the above action (3.1), the fundamental variables are g µν and the auxiliary metric q µν . It can be proven that the equations of motion derived from the original action (2.1) can be obtained unambiguously by varying the action (3.1) with respect to g µν and q µν . In our previous papers [40,42,43,45], we have used this alternative action (3.1) to deduce the classical Hamiltonian and the corresponding WDW equation in the EiBI gravity. It turns out that the construction of the WDW equation is much more straightforward because of the absence of the square root structure of the curvature present in the original action (2.1).
Using the alternative action (3.1) and assuming that the matter sector is described by a perfect fluid with energy density ρ and pressure p, the reduced Lagrangian associated with the action (3.1) can be written as Note that the energy density is expressed as a function of the physical scale factor a and the relation a = bY has been imposed. According to the definition of conjugate momenta, we have three primary constraints: where ∼ denotes the weak equality, i.e., the equality on the constraint surface. The total Hamiltonian is defined as follows where p b is the conjugate momentum of the phase space variable b. In the last few terms, λ X , λ Y , and λ M are Lagrange multipliers associated with each primary constraint. In Refs. [42,43], a thorough constraint analysis of this system has been carried out. In Ref. [45], an improved investigation has been done in which the matter sector is assumed to be a scalar field rather than a perfect fluid. As expected, the Hamiltonian itself is a first class constraint and at the quantum level, it would be treated as a restriction on the Hilbert space, giving rise to the WDW equation. In addition, the equations of motion (2.4), which relate algebraically the metrics and the energy-momentum tensor, are exactly the secondary constraints of the system and they are second class constraints. In the presence of second class constraints, one has to resort to the Dirac brackets to promote the phase space functions to quantum operators [53]. By doing so, the second class constraints can be directly regarded as zero operators and the WDW equation can be significantly simplified.
As a result, we can find a basis b| to write down the WDW equation b|Ĥ T |ψ = 0 as follows Note that the variables X and Y can be expressed as functions of b by imposing the second class constraints given by Eqs. (2.4) in this model (see also Eqs. (3.9) and (3.10) in Ref. [43]). Therefore, the potential V (b) can be expressed as follows

5)
A. The WDW equation for κ < 0 In the perfect fluid description, Eq. (3.4) and the potential (3.5) stands for a general expression of the WDW equation of the EiBI gravity. The explicit expression of X(b) and Y (b) are given by the constraints (2.4) and they depend on the cosmological solutions under consideration. In this subsection, we focus on the approximated cosmological solutions of the bouncing scenario for a negative κ, which has been discussed in subsection II A. In this case, we insert the approximated behaviors (2.8) to the potential V (b) and the potential can be approximated as when b → 0. It can be also seen that when b → 0, the last term on the right hand side dominates. We choose the following factor orderingp and introduce a new variable The WDW equation when b → 0 (y → 0) can be written as In this subsection, we shall derive the approximated WDW equation when b → 0 in the EiBI gravity with a positive κ. In this case, we insert the asymptotic equations (2.12) into the potential (3.5), and the potential can be approximated as It can be seen that when b → 0, the second term on the right hand side dominates.
To proceed, we choose the same factor ordering as Eq. (3.7), and introduce a new variable The WDW equation when b → 0 (z → 0) can be written as where c 2 is a positive constant and it is defined as Note that only the second term on the right hand side of Eq. (3.10) is considered since it dominates the potential when b → 0.

IV. THE WDW EQUATION: SCALAR FIELD
In the previous section, we have derived the WDW equations near the singularity of the auxiliary metric by using a perfect fluid description. The quantum system has only one degree of freedom in the sense that the WDW equation turns out to be an ordinary differential equation of a single variable, the scale factor b. However, the assumption of the perfect fluid description is just for convenience and may not be complete to describe the quantum evolution of the universe in a satisfactory manner. For the sake of completeness, in this section we will introduce an additional degree of freedom, the scalar field φ minimally coupled to the EiBI gravity, to describe the matter sector of the gravitational system. In the classical regime, it is well-known that the properties of a perfect fluid, including its equation of state and evolution, can be described by a scalar field when a corresponding potential V (φ) is chosen. In the quantum regime, on the other hand, the two degrees of freedom from the geometrical sector and from the matter sector do not necessarily relate to each other as in the classical regime. In this regard, the WDW equation becomes a partial differential equation of two variables b and φ. To have a more complete picture of the quantum behavior of the universe near the singularity, we will solve the wave function and investigate how the wave function evolves in the two dimensional (b, φ) space. We shall mention that in Ref. [45], we have studied the quantum avoidance of the big rip singularity in the EiBI gravity by solving the WDW equation with two degrees of freedom, one from the geometrical sector and the other is the phantom scalar field from the matter sector.
Considering a homogeneous and isotropic metric and assuming a scalar field minimally coupled to the gravity sector, the reduced Lagrangian can be written as follows [45]: After the Legendre transformation, the total Hamiltonian can be obtained as follows where p φ is the conjugate momentum of the variable φ. Note that λ M , λ X , and λ Y are the Lagrange multipliers associated with the primary constraints as in the case of a perfect fluid description. The complete constraint analysis of this system has been carried out in Ref. [45]. It turns out that the Hamiltonian is again a first class constraint as expected and it becomes a restriction in the Hilbert space at the quantum level. p M is another first class constraint and it corresponds to a gauge degree of freedom which can be fixed by assuming the lapse function M to be a constant. It should be stressed that in the EiBI theory, there are two additional second class constraints which correspond to two algebraic relations in the theory. In the perfect fluid description, these second class constraints are given by Eqs. (2.4). In the scalar field description, on the other hand, these constraints are also given by Eqs. (2.4) but one has to substitute the energy density and the pressure by their scalar field counterparts, ρ φ and p φ , respectively. The explicit expressions of these second class constraints are given in Eqs. (3.9) and (3.10) in Ref. [45]. Essentially, once we introduce the Dirac brackets to promote the phase space functions to quantum operators, the second class constraints can be regarded as zero operators [53]. As a result, the total Hamiltonian can be significantly simplified [45] Note that at the quantum level, the last two terms can be omitted since p X and p Y are also second class constraints. We shall emphasize that Eq. (4.3) is still not the final expression of the WDW equation that we are going to study. The final expression of the WDW equation is expected to be a partial differential equation of b and φ. After inserting p X ∼ 0 and p Y ∼ 0, there remain two variables X and Y in Eq. (4.3). Technically, we have to use again the second class constraints to relate these two variables to the phase space variables b, p b , φ, and p φ . It can be expected that these relations would depend on the cosmological models under consideration and also on the scalar field potential that we choose in the model. Different choices of the cosmological solutions and potentials certainly change the expressions of the second class constraints, hence change the quantization of the system and also the expression of the WDW equation. In the following two subsections, we will first consider the cosmological solution near the singularity of the auxiliary metric for a negative κ and rewrite the WDW equation as a partial differential equation from which the wave function can be solved. A similar study for a positive κ will be presented in the subsection IV B.
A. The WDW equation for κ < 0 The energy density and its pressure can be described by a scalar field φ and its potential V (φ) via the following relations Note that we have used a lower index φ to highlight that the energy density and the pressure are described through the dynamics of a scalar field. For the EiBI gravity with a negative κ, the physical Hubble rate near the bounce a → a m1 (b → 0) is approximated as in Eq. (2.6). On the other hand, the conservation of the energy-momentum tensor implies that the energy density and the pressure, if expressed as a function of the scale factor a, read ρ φ ≈ ρ 0 a −3(1+w) and p φ ≈ wρ 0 a −3(1+w) , respectively, where ρ 0 is an integration constant and w stands for the equation of state defined by w ≡ p φ /ρ φ . Using the approximated Hubble rate Eq. (2.6) and Eqs. (4.4) and (4.5), we obtain the asymptotic expression of the scalar field as a function of the scale factor a near the bounce: where δa ≡ a − a m1 and On the above equations, φ 0 is an integration constant and it is the value of the scalar field when δa = 0. Furthermore, using again Eqs. (4.4) and (4.5), the scalar field potential can be expressed with the equation of state w and it approaches a constant (1 − w) , (4.8) when δa → 0. Now, we shall rewrite the WDW equation (4.3) in such a way that it only contains the phase space variables b, φ, and their conjugate momenta. According to Eqs. (2.8), we get the following approximated equations: when b → 0. Using Eq. (4.9), the Hamiltonian can be written as We would like to stress again that if we use the Dirac brackets to promote the phase space functions to quantum operators, the second class constraints can be treated as zero operators. That is why we can use Eqs. (4.9) to simplify the WDW equation. Actually, the second class constraints that correspond to the above substitutions are secondary constraints, which are also the equations of motion of the theory, i.e., Eqs. (2.7). More strictly speaking, one should first use the Dirac brackets to promote the phase space functions to quantum operators, obtaining the Hamiltonian operatorĤ T . Then one regards the second class constraints as zero quantum operators and make the substitutions mentioned above. In the aforementioned procedures, it looks like we are doing it the other way around, that is, making substitutions at the classical level then promoting the phase space functions to operators. In either case, the final expression of the WDW equation is the same. Therefore, we will use the Hamiltonian (4.10) to construct the WDW equation. After promoting the phase space functions to quantum operators and simplifying the Hamiltonian operator with the second class constraints, we construct the WDW equation as follows and choose the following factor ordering (4.12) Note that we have defined a new variable to label the scale factor of the auxiliary metric. Finally, the WDW equation becomes Near the singularity (x → 0) of the auxiliary metric, the WDW equation (4.14) can be further approximated as In the next section, we will solve this equation (4.15) to get the asymptotic behavior of the wave function ψ(x, φ) near the singularity (x → 0, φ → φ 0 ) with a negative κ.
B. The WDW equation for κ > 0 As shown explicitly in section II, the EiBI theory of gravity with a positive κ resolves the big bang singularity quite differently as compared with the situation for a negative κ. When κ > 0, the asymptotic behavior of the physical Hubble function is given by Eq. (2.11) when the physical scale factor approaches its minimum value a → a m2 (b → 0). Note that this happens at t → −∞ for a constant lapse function N . If the density and pressure of the fluid are related via a constant equation of state effectively, that is, p φ = wρ φ , their relations with the physical scale factor can be obtained from the conservation equation: ρ φ ≈ ρ 0 a −3(1+w) and p φ ≈ wρ 0 a −3(1+w) . Combining Eqs. (2.11), (4.4) and (4.5), we obtain the asymptotic expression of the scalar field as a function of δa ≡ a − a m2 as follows It can be seen that φ → −∞ when δa and b vanish. The scalar field potential approaches a constant when a → a m2 and its value depends on the equation of state w: (1 − w) . Similar to what we have done in the previous subsection, we have to express X/Y 3 and Y 2 in Eq. (4.3) in terms of b, φ, and their conjugate momenta. To do this, we use Eqs. (2.12) to get Substituting (4.19) into the Hamiltonian, we obtain To proceed, we use the following factor ordering: (4.23) Note that we have defined a new variable and it can be seen that z → −∞ when b → 0. Finally, the WDW equation can be expressed as

V. THE WAVE FUNCTIONS IN THE PERFECT FLUID DESCRIPTION
In the previous sections, we have obtained the asymptotic expressions of the WDW equation near the singularity (b → 0) of the auxiliary metric or equivalently the physical connection. For a negative κ, we have derived the WDW equations (3.9) and (4.15), by assuming that the matter field is governed by a perfect fluid and a scalar field, respectively. On the other hand, for a positive κ, the corresponding WDW equations with a perfect fluid and a scalar field have been obtained in Eqs. (3.12) and (4.25), respectively. We will solve the wave functions for all these WDW equations and see whether the wave functions satisfy the DW boundary condition, i.e., the wave functions vanish, near the configuration of the singularity of the auxiliary metric. Let us first consider the cases in which the matter field is described by a perfect fluid and solve the WDW equations (3.9) and (3.12). A.
The κ < 0 case When the matter content is governed by a perfect fluid, the WDW equation for a negative κ that we will take into account is given by Eq. (3.9). The general solution can be written as a linear combination of two independent solutions as follows ] approaches a non-zero constant. Therefore, the wave function (5.1) satisfies the DW condition near the singularity as long as one assumes C 2 = 0.

B.
The κ > 0 case On the other hand, the WDW equation for a positive κ within the perfect fluid description is given by Eq. (3.12). Depending on the value of the parameter c 2 , which is positive (see Eq. (3.13)), the general solution can be categorized as follows where D 1 , D 2 ,D 1 ,D 2 ,D 1 , andD 2 are integration constants. Next, we will consider different values of c 2 and investigate whether the wave function is able to satisfy the DW boundary condition near the singularity where z → 0.
• If 0 < c 2 < 1/4, the general solution of the wave function is given by Eq. (5.3). It can be shown that for each independent solution, the variable z has a positive power. Therefore, the general solution vanishes when z → 0, satisfying the DW condition at the singularity.
• If c 2 = 1/4, the general solution is given by Eq. (5.4) and the DW condition at z → 0 is unambiguously satisfied due to the factor √ z.
• If c 2 > 1/4, the general solution is given by Eq. (5.5) and the power of z is complex. In consequence, the wave function acquires an oscillating behavior described by the imaginary part of the power of z. However, the modulus of the wave function behaves as |ψ| ≈ √ z. Therefore, when z → 0, the modulus of the wave function vanishes and the DW condition is fulfilled.
Consequently, in a perfect fluid description and when κ > 0, the wave function always satisfies the DW condition at the singularity of the auxiliary metric.

VI. THE WAVE FUNCTIONS IN THE SCALAR FIELD DESCRIPTION
For the scalar field description, the asymptotic expressions of the WDW equations near the singularity are given by Eqs. (4.15) and (4.25), corresponding to a negative and a positive value of κ, respectively. The WDW equations are partial differential equations with two independent variables. We will prove that even in these general cases in which one more degree of freedom is included into the system, the DW boundary condition can still be satisfied near the singularity of the auxiliary metric. A.
The κ < 0 case The WDW equation in the scalar field description with a negative κ is given by the partial differential Eq. (4.15). The general solution to Eq. (4.15) can be obtained by using the separation of variables such that the total wave function can be decomposed as a series of products of two single variable functions where C k (x) and ϕ k (φ) are the solutions to the following two ordinary differential equations and k corresponds to the decoupling constant. The above ordinary differential equations (6.2) can be solved to get the solution of the gravitational part C k (x) 3) and the solution of the matter part ϕ k (φ) It should be stressed that near the singularity of the auxiliary metric, the scalar field approaches a constant and the auxiliary scale factor vanishes, i.e., φ → φ 0 and x → 0. Therefore, the matter part of the wave function ϕ k (φ) is welldefined and the gravitational part tends to a constant value E 1,k + E 2,k . Since the total wave function is constructed by the product of C k and ϕ k , the DW boundary condition can be satisfied as long as one requires E 1,k = −E 2,k .

B.
The κ > 0 case For a positive value of κ, the asymptotic expression of the WDW equation in the scalar field description near the singularity of the auxiliary metric is given by Eq. (4.25). Again, the partial differential equation can be solved by using the separation of variables. The total wave function can be decomposed as a series of products of the solutions corresponding to the gravitational part and matter part: where Q m (z) and ξ m (φ) are, respectively, the gravitational and matter part of the wave function. In this regard, the WDW equation (4.25) can be decoupled into two ordinary differential equations as follows where m is the value of the decoupling constant. The general solution to the gravitational part can be written in terms of the modified Bessel functions I µ [g (z)] and K µ [g (z)] as follows [54]: where G 1,m and G 2,m are integration constants. The order µ and the argument g (z) of the modified Bessel functions can be explicitly expressed as 2 κ e 2z . (6.9) On the other hand, the solution of the matter part can be solved as follows where H 1,m and H 2,m are integration constants. Note that the scalar field φ → −∞ near the singularity of the auxiliary metric.
To further proceed, we assume that the decoupling constant m is a real number. This assumption is fully physical as m has dimension of energy. Under this assumption, the order µ acquires either a non-negative real value when m ≥ 0, or a purely imaginary value when m < 0. Depending on the value of m, the asymptotic expressions of the modified Bessel functions at small arguments (z → −∞ and g(z) → 0) are given as follows [54] I µ [g (z)] ≈ 1 Γ (µ + 1) g (z) 2 µ when µ = −1, −2, −3, ... , (6.11) when µ is real and positive , (6.12) when µ = 0 , (6.13) K µ [g (z)] ≈ − π ν sinh (πν) 1 2 sin ν ln g(z) 2 when µ = iν is purely imaginary , (6.14) where Γ (α) stands for the Gamma function. In the following, we will investigate the behaviors of the total wave function for different values of the decoupling constant m.
• If m < 0, the matter part of the wave function (6.10) turns out to be a plane wave solution whose oscillating amplitude is constant. As for the gravitational part, the order µ becomes imaginary, and therefore according to Eqs. (6.11) and (6.14), the modified Bessel functions I µ [g (z)] and K µ [g (z)] are both rapidly oscillating functions with a non-zero constant modulus. In consequence, for m < 0 the total wave function does not vanish at the singularity and the DW condition cannot be satisfied.
• If m > 0, the matter part of the wave function turns out to be exponential functions. If we assume H 2,m = 0, the growing part of the solution when φ → −∞ is removed. On the other hand, the order µ of the modified Bessel functions in the gravitational part is a positive and real number. In this case, it can be seen from Eqs. (6.11) and (6.12) that the modified Bessel function I µ [g (z)] vanishes when g(z) → 0, while K µ [g (z)] diverges. Consequently, one has to further choose G 2,m = 0 in order to ensure the DW condition near the singularity of the auxiliary metric.
In summary, we have found that if m ≤ 0, it is impossible to obtain a wave function satisfying the DW boundary condition at the singularity of the auxiliary metric. In fact, one is supposed to impose an additional condition on the decoupling constant, i.e., m > 0, such that the DW condition is able to be satisfied.

VII. CONCLUSIONS
In the context of the EiBI gravity, it has been shown that the propagation of gravitational waves would be affected by the geometry of the auxiliary metric, which is compatible with the affine connection of the theory. Therefore, even though the big bang singularity can be resolved, the singularity is present in the auxiliary metric and it has an important consequence on the behavior of the linear perturbations. The linear perturbations, including the tensor modes, turn out to be unstable in the non-singular solutions within the EiBI theory. In this paper, we consider the quantum geometrodynamical approach in the context of the EiBI gravity. Note that the Born-Infeld type of theories seem to have intrinsic Noether symmetries as shown recently in [47]. This also supports the choice of the EiBI action in this paper. Our motivation is to see whether or not the singularity in the auxiliary metric can be ameliorated by the quantum effects. It turns out the answer is yes and therefore, the linear instabilities of the physical metric, which are associated with the singular behavior of the auxiliary metric, would be resolved by the same token.
For the sake of completeness, we have considered two descriptions regarding the matter sector of the theory. In the perfect fluid description, the matter field is governed by a perfect fluid with a positive and constant equation of state. In the homogeneous and isotropic universe, the system is characterized by a single variable b, the scale factor of the auxiliary metric. In the second description, that is, the scalar field description, we introduce a scalar field degree of freedom to incorporate the matter sector which, in the classical level, describes the evolution of the corresponding perfect fluid in the perfect fluid description. In this setup, the system contains two canonical degrees of freedom, the scale factor b and the scalar field φ, spanning a two dimensional configuration space.
In the framework of quantum geometrodynamical approach, the building block is the WDW equation describing the quantum evolution of the universe as a whole. Essentially, we start with the alternative EiBI action in the Einstein frame and derive the classical Hamiltonian for both descriptions mentioned above. The Hamiltonian constraint, which is a first class constraint, is regarded as a restriction on the Hilbert space and the WDW equation is derived by promoting all phase space functions to quantum operators. The commutation relations are constructed by using the Dirac brackets which are necessary for a system containing second class constraints. We have derived the asymptotic expressions of the WDW equations for the two matter descriptions, and for positive and negative values of the Born-Infeld parameter κ. For a negative value of κ, the physical metric bounces in the past at the classical level. The asymptotic expressions of the WDW equations near the bounce are given by Eqs. (3.9) and (4.15), for the perfect fluid and the scalar field descriptions, respectively. For a positive value of κ, the physical metric acquires its minimum scale factor in the asymptotic past (the loitering effect). The approximated WDW equations are given by Eqs. (3.12) and (4.25), for the perfect fluid and the scalar field descriptions, respectively.
After deriving the WDW equations, we have studied the quantum behavior of the universe by solving the wave function as a solution to the WDW equations. We have found that for each WDW equation under consideration, wave functions which satisfy the DW boundary conditions at the singularity of the auxiliary metric can always be obtained. Therefore, the hidden singularity in the auxiliary metric is expected to be avoided at the quantum level and the linear instabilities are not harmful in the quantum world.