{'fraction_non_alphanumeric': 0.03186213115575372, 'fraction_numerical': 0.008366314312910175, 'mean_word_length': 4.465160848733744, 'pattern_counts': {'":': 0, '<': 7, '<?xml version=': 0, '>': 3, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 7, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Classical and quantum massive cosmology for the open FRW universe

In an open Friedmann-Robertson-Walker (FRW) space background, we study the classical and quantum cosmological models in the framework of the recently proposed nonlinear massive gravity theory. Although the constraints which are present in this theory prevent it from admitting the flat and closed FRW models as its cosmological solutions, for the open FRW universe, it is not the case. We have shown that, either in the absence of matter or in the presence of a perfect fluid, the classical field equations of such a theory adopt physical solutions for the open FRW model, in which the mass term shows itself as a cosmological constant. These classical solutions consist of two distinguishable branches: One is a contacting universe which tends to a future singularity with zero size, while another is an expanding universe having a past singularity from which it begins its evolution. A classically forbidden region separates these two branches from each other. We then employ the familiar canonical quantization procedure in the given cosmological setting to find the cosmological wave functions. We use the resulting wave function to investigate the possibility of the avoidance of classical singularities due to quantum effects. It is shown that the quantum expectation values of the scale factor, although they have either contracting or expanding phases like their classical counterparts, are not disconnected from each other. Indeed, the classically forbidden region may be replaced by a bouncing period in which the scale factor bounces from the contraction to its expansion eras. Using the Bohmian approach of quantum mechanics, we also compute the Bohmian trajectory and the quantum potential related to the system, which their analysis shows are the direct effects of the mass term on the dynamics of the universe.


Introduction
General Relativity (GR) introduced by Einstein began a renaissance in scientific thought which changed our viewpoint on the concept of space-time geometry and gravity. The interpretation of gravitational force as a modification of geometrical structure of space-time made and makes this force distinguishable from other fundamental interactions, although there are arguments which support the idea that the other interactions may also have geometrical origin. Because of the unknown behavior of gravitational interaction at short distances, this distinction may have some roots in the heart of problems with quantum gravity. Therefore, any hope of dealing with such concepts would be in vain unless a reliable quantum theory of gravity can be constructed. In the absence of a full theory of quantum gravity, it would be then useful to describe its quantum aspects within the context of modified theories of gravity. From a field theory point of view, the gravitational force in GR can be represented as a field theory in which the space-time metric plays the role of the fields and the particle that is responsible to propagate gravity is named graviton. Then, naturally in comparison to the other field theories, one may ask about the different properties of such a particle. The answer to this question is deduced by the linearized form of GR and expansion of the space-time metric g µν , around a fixed background geometry η µν , as g µν = η µν + h µν , where h µν is the field representation of the graviton. Eventually, it is possible to show that the graviton is a massless spin-2 particle.
Then, since our knowledge about the behavior of gravity at very long distances is also incomplete, a question arises: Is it possible to consider a small nonvanishing mass for the graviton, i.e., a massive spin-2 particle? In the first attempts to deal with this question, it seemed that adding a mass term to the action may be sufficient. This was done by Fierz and Pauli [1]. However, it was shown that by considering the number of degrees of freedom, this model suffers from the existence of a ghost field, the so-called Boulware-Deser ghost [2], after studying the non-linear terms. This fact made massive gravity an abandoned theory for a while. Recently, de Rahm and Gabadadze proposed a new scenario in which they have shown that it is possible to have a ghost-free massive gravity even at the non-linear level [3]. That was a positive signal in this area, and the early results in this subject have been followed by a number of works that address different aspects of massive gravity [4]. As in the case of the other modified theories of gravity, it is important to seek cosmological solutions in the newly proposed massive theory of gravity. This is done by the authors of Ref. [5], who show that the existence of some constraints prevent the theory from having the nontrivial homogeneous and isotropic cosmological solutions. Indeed, what is shown in Ref. [5] is that, beginning with the flat FRW ansatz in the context of massive gravity, the corresponding field equations result in nothing but the Minkowski metric. However, by re-examination of the conditions, the authors of Refs. [6,7] have shown that, for the open FRW model, this is not the case, and the nonlinear massive gravity admits the open FRW as a compatible solution for its field equations. Another progresses to find the massive cosmologies lie in the field of the bi-metric theories of gravity; see, for instance, Ref. [8] based on the works of Hassan and Rosen [9], in which they show that a bi-metric representation for massive gravity exists.
Our purpose in the present paper is to continue the works of the authors of Ref. [6,7] in greater detail, based on the Hamiltonian formalism of the open FRW cosmology in the framework of massive gravity. We obtain the solutions to the vacuum and perfect fluid classical field equations and investigate their different aspects, such as the roll of the graviton's mass as a cosmological constant, the appearance of singularities, and the late time expansion. We then consider the problem at hand in the context of canonical quantum cosmology to see how the classical picture will be modified. Our final results show that the singular behavior of the classical cosmology will be replaced by a bouncing one when quantum mechanical considerations are taken into account. This means that the quantization of the model suggests the existence of a minimal size for the corresponding universe. We shall also study the quantum model by the Bohmian approach of quantum mechanics to show how the mass term exhibits its direct effects on the evolution of the system. The structure of the paper is as follows. In section 2, we briefly present the basic elements of the issue of massive gravity and its canonical Hamiltonian for a given open FRW universe. In section 3, classical cosmological dynamics is introduced for the vacuum and perfect fluid. Quantization of the model is the subject of section 4, and in section 5, the Bohmian approach of quantum mechanics is applied to the model. Finally, the conclusions are summarized in section 6.

Preliminary set-up
In this section, we start by briefly studying the nonlinear massive gravity action presented in Refs. [6,7] for the open FRW model, where the metric is given by with N (t) and a(t) being the lapse function and the scale factor, respectively, and K = −1 denoting the curvature index. Here we work in units where c =h = 16πG = 1. In the massive gravity scenario one considers a metric perturbation as [5,10] where η ab = diag(−1, 1, 1, 1) and φ a (x) are four scalar fields known as Stückelberg scalars and are introduced to keep the principle of general covariance also in massive general relativity [11]. It is clear that the first term in (2) is a representation of the Minkowski space-time in terms of the coordinate system (φ 0 , φ i ) and thus the tensor H µν is responsible for describing the propagation of gravity in this space. The action of the model consists of the gravitational part S g and the matter action S m as The matter part of the action is independent of the massive corrections to the gravity part. Also, the gravity part can be expressed in terms of the usual Einstein-Hilbert, with an additional correction term coming from the massive graviton; that is [5] in which all of the modifications due to the mass and also the interactions between the tensor fields H µν and g µν are summarized in the potential U (g, H). By using ghost-free conditions for the theory in Ref. [11], we propose the following form for the potential term: [10] where in which the tensor K µν is defined as and the notations < K >= g µν K µν , < K 2 >= g αβ g µν K αµ K βν ,... are used for the corresponding traces. Now, equations (4)-(7) describe the gravitational part of the action for a massive gravity theory. Since its explicit form directly depends on the choice of scalar fields φ a (x), it is appropriate to concentrate on this point first.
Interesting forms for such fields should involve terms which would describe a suitable coordinate transformation on the Minkowski space-time. In a flat FRW background, for instance, one may select φ 0 = f (t) and φ i = x i , as is used in [5]. Here, for the open FRW metric (1), we use the following ansatz proposed in [6] φ Upon substitution of these scalar fields and also the definition of the Ricci scalar into the relations (4)-(7), we are led to a point-like form for the gravitational Lagrangian in the minisuperspace {N, a, f } as where in which an overdot represents differentiation with respect to the time parameter t. It is seen that this Lagrangian does not involveṄ , which means that the momentum conjugate to this variable vanishes.
In the usual canonical formalism of general relativity, we know this issue as the primary constraint in the sense that the variable N is not a dynamical variable but a Lagrange multiplier in the Hamiltonian formalism. On the other hand, Lagrangian (9) seems to show an additional constraint related to the Stückelberg scalars whose dynamics are encoded in the function f (t). We see that in spite of the common Lagrangians in which the first derivative of the configuration variables are of second order, f appears linearly in the Lagrangian (9). Therefore, by computing the momentum conjugate to f ; that is, P f = ∂Lg ∂ḟ , we obtain Now, it is clear that this relation is not invertible to obtainḟ (f, P f ). In such a case, the Lagrangian is said to be singular and the relations like (11), which hinder the inversion, are known as primary constraints. One may use the method of Lagrange multipliers to analyze the dynamics of the system by adding to the Lagrangian all of the primary constraints multiplied by arbitrary functions of time. However, to deal with our constrained system, we act differently and proceed as follows. We vary the Lagrangian (9) with respect to f to obtain The solutionȧ = |K|N of this equation is nothing but what we obtain from the variation of the usual Einstein-Hilbert Lagrangian with respect to N . Since its counterpart in massive gravity iṡ we cannot accept the relationȧ = |K|N as a physical solution. Therefore, the constraint corresponding to the dynamic of f (t) shows itself in the equation where using the same notation as in [6], its solutions can be written as As is argued in [6], in the limit where α 3 and α 4 are of the order of a small quantity ǫ, the expression of X + goes to infinity while X − → 3/2. Because of this limiting behavior, we use the subscript − in the following for numerical values of constants with subscript ±. Now we may insert the constraints (15) into the relations (10) to reduce the degrees of freedom of the system and obtain a minimal number of dynamical variables. If we do so, we obtain in terms of which the Lagrangian (9) takes its reduced form with only one physical degree of freedom a. The momentum conjugate to a is Noting that one gets Now, the Hamiltonian of the model can be obtained from its standard definition H =ȧP a − L, with result in which we have defined We see that the lapse function enters in the Hamiltonian as a Lagrange multiplier as expected. Thus, when we vary the Hamiltonian with respect to N , we get H g = 0, which is called the Hamiltonian constraint. On a classical level this constraint is equivalent to the Friedmann equation, wherein our problem at hand can be easily checked by comparing it with the equation of motion (4.5) in [6]. On a quantum level, on the other hand, the operator version of this constraint annihilates the wave function of the corresponding universe, leading to the so-called Wheeler-DeWitt equation. Now, let us deal with the matter field with which the action of the model is augmented. As we have mentioned, the matter part of the action is independent of modifications due to the mass terms. Therefore, the matter may come into play in a common way and the total Hamiltonian can be made by adding the matter Hamiltonian to the gravitational part of (20). To do this, we consider a perfect fluid whose pressure p is linked to its energy density ρ by the equation of state where −1 ≤ ω ≤ 1 is the equation of the stated parameter. According to Schutz's representation for the perfect fluid [12], its Hamiltonian can be viewed as (see [13] for details) where T is a dynamical variable related to the thermodynamical parameters of the perfect fluid and P T is its conjugate momentum. Finally, we are in a position in which can write the total Hamiltonian H = H g + H m as The setup for constructing the phase space and writing the Lagrangian and Hamiltonian of the model is now complete. In the following section, we shall deal with classical and quantum cosmologies which can be extracted from a theory with the previously mentioned Hamiltonian.

Cosmological dynamics: classical point of view
The classical dynamics are governed by the Hamiltonian equations. To achieve this purpose, we divide this section into two parts. We first consider the case in which the matter is absent, i.e., the vacuum, and then include the matter.

The vacuum classical cosmology
In this case, we can construct the equations of motion by the Hamiltonian equations with use of the Hamiltonian (20). Equivalently, one may directly write the Friedmann equation from the Hamiltonian constraint H = 0 which, as we mentioned previously, reflects the fact that the corresponding gravitational theory is a parameterized theory in the sense that its action is invariant under time reparameterization. Noting from (19) thaṫ equation (20) gives in which we have chosen the gauge N = 1, so that the time parameter t becomes the cosmic time τ . As is indicated in [6], this equation looks like the Friedmann equation for the open FRW universe with an effective cosmological constant Λ ± = c ± m 2 and admits the following solutions where τ * is an integration constant and we have taken Λ = c − m 2 . For a positive τ * , the condition a(τ ) ≥ 0 implies that the expressions of a + (τ ) and a − (τ ) are valid for τ ≥ τ * and τ ≤ −τ * respectively, such that a ± (τ * ) = 0. It is seen that the evolution of the corresponding universe with the scale factor a + (τ ) begins with a big-bang-like singularity at τ = τ * and then follows an exponential law expansion at late time of cosmic evolution in which the mass term shows itself as a cosmological constant. For a universe with the scale factor a − (τ ), on the other hand, the behavior is opposite. The universe decreases its size from large values of scale factor at τ = −∞ and ends its evolution at τ = −τ * with a zero size. In figure 1 we have plotted theses scale factors for typical values of the parameters. As this figure shows, although the behavior of a + (τ ) (a − (τ )) is like a de Sitter (a(τ ) ∼ e √ Λ/3τ ) universe at τ → ∞ (τ → −∞), in spite of the de Sitter, it begins (ends) its evolution with a singularity.
In summary, what we have shown previously is that in the framework of an open FRW background geometry, the vacuum solutions of the massive theory are equivalent to the solutions of the usual GR with a cosmological constant. Accordingly, the zero-size singularity of both theories has the same nature. In this sense we would like to emphasize that the metric (1) with the scale factor (27) is indeed a section of the de Sitter hyperboloid embedded in a 5-dimensional Minkowski space To see this, one may parameterize the hyperboloid in terms of the spherical coordinates (r, θ, φ) as which, upon substitution into the metric (29), yields the open FRW metric with the scale factor a(τ ) = sinh τ . This means that the point a = 0 can be viewed as a coordinate singularity. However, we have to note that in the presence of any kind of matter field the point a(τ * ) = 0 represents a true singularity. Thus, our following analysis to quantize the model is based on the minisuperspace coordinate system in terms of which the dynamical representation of the metric, i.e. (1), is written.
In the next section, we shall see how the previous picture may be modified when one takes into account quantum mechanical considerations.

Perfect fluid classical cosmology
Now, we assume that a perfect fluid in its Schutz's representation is coupled with gravity. In this case the Hamiltonian (24) describes the dynamics of the system. The equations of motion for T and P T read asṪ A glance at the above equations shows that with choosing the gauge N = a 3ω , we shall have which means that variable T may play the role of time in the model. Therefore, the Friedmann equation H = 0 can be written in the gauge N = a 3ω as follows where we take P T = P 0 = const. from the second equation of (31). Since it is not possible to find the analytical solutions of the above differential equation for any arbitrary ω, we present its solutions only in some special cases.
• ω = − 1 3 : cosmic string. In this case we obtain where t 0 is an integration constant. We see that the evolution of the universe based on (34) has big-bang-like singularities at t = t 0 ± t * where t * = √ 3(P 0 +3|K|) Λ . Indeed, the condition a 2 (t) ≥ 0 separates two sets of solutions a (I) (t) and a (II) (t), each of which is valid for t ≤ t 0 − t * and t ≥ t 0 + t * , respectively. For the former, we have a contracting universe which decreases its size according to a power law relation and ends its evolution in a singularity at t = t 0 − t * , while for the latter, the evolution of the universe begins with a big-bang singularity at t = t 0 + t * and then follows the power law expansion at late time of cosmic evolution.
One may translate these results in terms of the cosmic time τ . Using its relationship with the time parameter t in this case, that is, dτ = a −1 (t)dt, we are led to where τ * = 1 2 3 Λ ln 3(P 0 + 3|K|). Again, it is seen that there is a classically forbidden region τ 0 − τ * < τ < τ 0 + τ * , for which we have no valid classical solutions. For τ ≤ τ 0 − τ * , the universe has a exponential decreasing behavior which ends its evolution in a singular point with zero size at τ = τ 0 − τ * , while in the region τ ≥ τ 0 + τ * it begins with the big-bang singularity at τ = τ 0 + τ * and then grows exponentially forever.
• ω = −1: cosmological constant. Performing the integration, we get the following implicit relation between t and a(t): In terms of the cosmic time τ , it is easy to see that this solution returns to (27), in which the cosmological term is replaced by Λ → Λ + cons. This is expected because the solutions (27) were equivalent to an open FRW universe with a cosmological constant. Therefore, adding a new cosmological term (a perfect fluid with ω = −1) only makes a shift in the corresponding cosmological constant.

Cosmological dynamics: quantum point of view
In this section we look for the quantization of the model presented above via the method of canonical quantization. As is well known, this procedure is based on the Wheeler-DeWitt equationĤΨ = 0, whereĤ is the operator version of the Hamiltonian constraint and Ψ is the wave function of the universe, a function of the 3-geometries and the matter fields. As in the case of the classical cosmology, we consider the matter of free and perfect fluid quantum cosmology separately.
Before going to the subject, a remark is in order related to the Hamiltonians (20) and (24). The term in the round bracket in these Hamiltonians is like the Hamiltonian of a charged particle moving in an electromagnetic field. From this analogy, one may define the transformation to simplify the form of the classical Hamiltonian. It is clear that this is a canonical transformation both classically and quantum mechanically [15]. Since going back from a new set of variables to the old ones in a classical canonical transformation can be made without any ambiguity, applying this transformation may not be important for the classical dynamics presented in the previous section.
In the context of quantum mechanics, on the other hand, the subject is of little difference. The transition to the quantum version of the theory is achieved by promoting observables to operators which are not necessarily commuting. Thus, by replacing the canonical variables (a, P a ) by their operator counterparts (â,P a = −id/da), we obtain the quantum Hamiltonian where ... denotes the terms out of the round bracket in expressions (20) or (24). When calculating the square, it should be noted that the operatorsâ andP a do not commute. Although the order of these operators does not matter in the classical analysis, quantum mechanically this issue is quite crucial. Indeed, this is the operator ordering problem and, unfortunately, there is no well defined principle which specifies the order of operators in the passage from classical to quantum theory. There are, however, some simple rules which one uses conventionally. If, for instance, we order the products of a andP a inΠ a 2 such that the momentum stands to the right of the scale factor, we obtain in which we have used the commutation relation [â,P a ] = i. With this expression at hand, there is still another factor ordering ambiguity in the termsâ −1P a 2 andâ 2P a to construct the quantum Hamiltonian (38). As Hawking and Page have shown [16], the choice of different factor ordering will not affect semiclassical calculations in quantum cosmology, so for convenience one usually chooses a special place for it in the special models. However, in general, the behavior of the wave function depends on the chosen factor ordering [17]. In what follows, as one usually does in the minisuperspace approximation to the cosmological models, we work in the framework of a special factor ordering in which, in addition to the expression (39) forΠ a 2 , we also use the orderingsâ −1P a and a 2P a =âP aâ to make the Hamiltonian hermitian 1 .

The vacuum quantum cosmology
In this case, with the help of the Hamiltonian (20) and use of the abovementioned choice of ordering, the Wheeler-DeWitt equation reads 1 With the canonical transformation (37) at hand, one may uses the transformed Hamiltonian H = − 1 12a Π 2 a + ..., to quantize the system, where again ... denotes the terms out of the round bracket in expressions (20) or (24). Using this Hamiltonian in the hermitian form a −1 Π 2 a = Πaa −1 Πa and also representing Πa by −i∂a, this is equivalent to our above treatment in which the last term in (39) is absent. Therefore, one may have some doubts on the validity of the main following results due to the effects of the chosen factor ordering. To overcome this problem, we have made some calculations based on the above mentioned transformed Hamiltonian and have verified that the general patterns of the resulting wave functions follow the behavior shown in following sections.
This equation does not seem to have analytical solutions. However, we can get some properties of its solutions in special regions where there is interest in classical and quantum regimes. First of all, let us rewrite this equation in the form in which the quantum potential is defined as It is well-known that the quantum effects are important for small values of the scale factor and in the limit of the large scale factor can be neglected. Therefore, in the semiclassical approximation region we can omit the Q term in (42) and obtain dS da = 3Λa 3 ± a 36 + 12Λa 2 .
In the WKB method, the correlation between classical and quantum solutions is given by the relation P a = ∂S ∂a . Thus, using the definition of P a in (19), the equation for the classical trajectories becomeṡ from which one finds which shows that the late time behavior of the classical cosmology (27) is exactly recovered. The meaning of this result is that for large values of the scale factor, the effective action corresponding to the expanding and contracting universes is very large and the universe can be described classically.
On the other hand, for small values of the scale factor we cannot neglect the quantum effects, and the classical description breaks down. Since the WKB approximation is no longer valid in this regime, one should go beyond the semiclassical approximation. In the quantum regime, if we neglect the term Λ 2 a 6 in (41), the two linearly independent solutions to this equation can be expressed in terms of the Hermite H ν (x) and hypergeometric F 1 1 (a, b; z) functions, leading to the following general solution: At this step we take a quick glance at the question of the boundary conditions on the solutions to the Wheeler-DeWitt equation. Note that the minisuperspace of the above model has only one degree of freedom denoted by the scale factor a in the range 0 < a < ∞. According to [18], its nonsingular boundary is the line a = 0, while at the singular boundary this variable is infinite.
Since the minisuperspace variable is restricted to the abovementioned domain, the minisuperspace quantization deals only with wave functions defined on this region. Therefore, to construct the quantum version of the model, one should take into account this issue. This is because in such cases, one usually has to impose boundary conditions on the allowed wave functions; otherwise the relevant operators, especially the Hamiltonian, will not be self-adjoint. The condition for the Hamiltonian operatorĤ associated with the classical Hamiltonian function (20) and (24) to be self-adjoint is (ψ 1 ,Ĥψ 2 ) = (Ĥψ 1 , ψ 2 ) or Following the calculations in [19] and dealing only with square integrable wave functions, this condition yields a vanishing wave function at the nonsingular boundary of the minisuperspace. Hence, we impose the boundary condition on the solutions (46) such that at the nonsingular boundary (at a = 0), the wave function vanishes. This makes the Hamiltonian hermitian and self-adjoint and can avoid the singularities of the classical theory, i.e. there is zero probability for observing a singularity corresponding to a = 0. 2 Therefore, we require Note that equation (41) is a Schrödinger-like equation for a fictitious particle with zero energy moving in the field of the superpotential with the real part U (a) = −(36a 2 + 12Λa 4 ). Usually, in the presence of such a potential the minisuperspace can be divided into two regions, U > 0 and U < 0, which could be termed the classically forbidden and classically allowed regions, respectively. In the classically forbidden region the behavior of the wave function is exponential, while in the classically allowed region the wave function behaves oscillatorily. In the quantum tunneling approach [18], the wave function is so constructed as to create a universe emerging from nothing by a tunneling procedure through a potential barrier in the sense of usual quantum mechanics. Now, in our model, the superpotential is always negative, which means that there is no possibility of tunneling anymore, since a zero energy system is always above the superpotential. In such a case, tunneling is no longer required as classical evolution is possible. As a consequence the wave function always exhibits oscillatory behavior. In figure 2, we have plotted the square of the wave functions for typical values of the parameters. It is seen from this figure that the wave function has a well-defined behavior near a = 0 and describes a universe emerging out of nothing without any tunneling. (See [22], in which such wave functions also appeared in the case study of the probability of quantum creation of compact, flat, and open de Sitter universes.) On the other hand, the emergence of several peaks in the wave function may be interpreted as a representation of different quantum states that may communicate with each other through tunneling. This means that there are different possible universes (states) from which the present universe could have evolved and tunneled in the past, from one universe (state) to another.

Perfect fluid quantum cosmology
In this case, the Wheeler-DeWitt equation can be constructed by means of the Hamiltonian (24). With the same approximations as we used in the previous subsection, we obtain We separate the variables in this equation as  Figure 2: The square of the wave function for the quantum universe. We take the numerical value Λ = 1.5.
leading to The solutions of the above differential equation may be written in the form for ω = −1/3 and for ω = −1. Now the eigenfunctions of the Wheeler-DeWitt equation can be written as We may now write the general solution to the Wheeler-DeWitt equation as a superposition of its eigenfunctions; that is, where A(E) is a suitable weight function to construct the wave packets. The above relations seem to be too complicated to extract an analytical expression for the wave function. Therefore, in the following (for the case ω = −1/3), we present an approximate analytic method which is valid for very small values of scale factor, i.e., in the range that we expect the quantum effects to be important. In this regime if we keep only the a −1 and a 2 terms in the second and third terms of (51), the solutions to this equation can be viewed as a superposition of the functions sin √ 36+E−6iΛ 2 a 2 and cos √ 36+E−6iΛ 2 a 2 . If we impose the boundary condition ψ(a = 0) = 0 on these solutions, we are led to the following eigenfunctions: Now, by using the equality we can evaluate the integral over E in (55), and the simple analytical expression for this integral is found if we choose the function A(E) to be a quasi-Gaussian weight factor A(E) = e −γE (γ is an arbitrary positive constant and E = 36 + E − 6iΛ), which results in Using the relation (57) yields the following expression for the wave function where N is a numerical factor. Now, having this expression for the wave function of the universe, we are going to obtain the predictions for the behavior of the dynamical variables in the corresponding cosmological model. To do this, one may calculate the time dependence of the expectation value of a dynamical variable q as Following this approach, we may write the expectation value for the scale factor as which yields This relation may be interpreted as the quantum counterpart of the classical solutions (34). However, in spite of the classical solutions, for the wave function (59), the expectation value (62) of a never vanishes, showing that these states are nonsingular. Indeed, in (62) T varies from −∞ to +∞, and any T 0 is just a specific moment without any particular physical meaning like big-bang singularity. The above result may be written in terms of the cosmic time τ . By the definition dτ = a −1 (T )dT , we obtain the quantum version of the relations (35) as In figure 3, we have plotted the classical scale factors (34) and (35) and their quantum counterparts (62) and (63). As is clear from this figure, for a perfect fluid with ω = −1/3, the corresponding classical cosmology admits two separate solutions which are disconnected from each other by a classically forbidden region. One of these solutions represents a contracting universe ending in a singularity while another describes an expanding universe which begins its evolution with a big-bang singularity.
On the other hand, the evolution of the scale factor based on the quantum-mechanical considerations shows a bouncing behavior in which the universe bounces from a contraction epoch to a reexpansion era. Indeed, the classically forbidden region is where the quantum bounce has occurred. We see that in the late time of cosmic evolution in which the quantum effects are negligible, these two behaviors coincide with each other. This means that the quantum structure which we have constructed has a good correlation with its classical counterpart.

Bohmian trajectories
In the previous sections, we saw how the classical singular behavior of the universe was replaced with a bouncing one in a quantum picture. Now, a natural question may arise: Why will the bounce occur? Clearly, it is due to the quantum mechanical effects which show themselves when the size of the universe tends to very small values. However, we would like to know whether the massive correction to the underlying gravity theory has any contribution to this phenomenon. To deal with this question, let us return to the wave function (59) and write it in the polar form Ψ(a, T ) = Ω(a, T )e iS(a,T ) , where Ω(a, T ) and S(a, T ) are real functions, which simple algebra gives as .
According to the Bohm-de Broglie interpretation of quantum mechanics [23] and also its usage in quantum cosmology [24], upon using this form of the wave function in the corresponding wave equation, we arrive at the modified Hamilton-Jacobi equation as where P i are the momentum conjugate to the dynamical variables q i and Q is the quantum potential. With beginning of the wave equation (49), for which we have used the same approximations as in the previous section, the above mentioned procedure gives the quantum potential as On the other hand, the Bohmian equations of motion can be obtained by P a = ∂S ∂a , where by means of the relation (19) reads .
The solution to this equation denotes the Bohmian representation of the scale factor; that is where c is an integration constant and Ei(b; z) is the exponential integral function defined by The bouncing behavior of the scale factor is again its main property near the classical singularities as we have shown in figure 4. To achieve an expression for the quantum potential in terms of the scale factor, we note that all of our above calculations are in the vicinity of T ∼ 0, where the scale factor is small. In this regime, a numerical analysis shows that the Bohmian scale factor (69) behaves as a(T ) ∼ (γ 2 + T 2 ) 1/2 , in agreement with the expectation value (62). Thus, substituting in (64), we get the quantum potential from (67) as In figure 4 we also have plotted the qualitative behavior of the quantum potential versus the scale factor. As this figure shows, this potential goes to zero for the large values of the scale factor. This behavior is expected, since in this regime the quantum effects can be neglected and the universe evolves classically. On the other hand, for the small values of the scale factor the potential takes a large magnitude and the quantum mechanical considerations come into the scenario. This is where the quantum potential can produce a huge repulsive force which may be interpreted as being responsible of the avoidance of singularity. In figure 4 the horizontal line represents a constant energy level which in intersecting with the potential curves gives the turning points at which the bounce will occur. The solid curve in this figure is plotted in the case of Λ = 0; i.e., for the massive theory, while the dashed curve is for Λ = 0; i.e., for when the massive corrections are absent. It is seen that, although the mass term Λ is not the only reason for the bouncing behavior in the vicinity of the classical singularity, it may shift the bouncing point into the smaller values of the scale factor. This means that if we consider the bouncing point as the minimum size of the universe (which is suggested by quantum cosmology), then the massive version of the underlying gravity theory predicts a smaller value for this minimal size in comparison with the usual Einstein-Hilbert model. These facts and also other considerable possibilities such as quantum tunneling between different classically allowed regimes (as can be seen from figure 4) through the potential barrier support the idea that the massive corrections to the classical cosmology are some signals from quantum gravity.

Conclusions
In this paper we have applied the recently proposed nonlinear massive theory of gravity to an open FRW cosmological setting. Although the absence of homogeneous and isotropic solutions is one of the main challenges related to this kind of gravitational theory, we moved along the lines of [6,7], in which the existence of open FRW cosmologies is investigated. By using the constraint corresponding to the Stückelberg scalars, we reduced the number of degrees of freedom, according to which the total Hamiltonian of the model is deduced. We then presented in detail, the classical cosmological solutions either for the empty universe or in the case where the universe is filled by a perfect fluid (in its Schutz representation) with the equation of state parameter ω = −1/3, −1. We saw that in both of these cases, the solutions consist of a contraction universe which finalizes its evolution in a singular point and an expanding universe which begins its dynamic with a big-bang singularity. These two branches of solutions are disconnected from each other by a classically forbidden region. Also, the common feature of the vacuum and matter classical solutions is that the mass term plays a role which resembles the role of cosmological constant in the usual de Sitter universe. In this sense we may relate the massive corrections of GR to the problem of dark energy.
In another part of the paper, we dealt with the quantization of the model described above via the method of canonical quantization. For an empty universe, we have shown that by applying the WKB approximation on the Wheeler-DeWitt equation, one can recover the late time behavior of the classical solutions. For the early universe, we obtained oscillatory quantum states free of classical singularities by which two branches of classical solutions may communicate with each other. In the presence of matter, we focused our attention on the approximate analytical solutions to the Wheeler-DeWitt equation in the domain of small scale factor, i.e. in the region which the quantum cosmology is expected to be dominant. Using Schutzs representation for the perfect fluid, under a particular gauge choice, we led to the identification of a time parameter which allowed us to study the time evolution of the resulting wave function. Investigation of the expectation value of the scale factor shows a bouncing behavior near the classical singularity. In addition to singularity avoidance, the appearance of bounce in the quantum model is also interesting in its nature due to prediction of a minimal size for the corresponding universe. We know the idea of existence of a minimal length in nature is supported by almost all candidates of quantum gravity. Finally, we repeated the quantum calculations by means of the Bohmian approach to quantum mechanics. The analysis of the quantum potential shows the importance of the mass term in the action of the model. Indeed, we have shown that in the presence of the massive graviton, the quantum potential changes its behavior from an infinite barrier to a finite one, and hence the minimal size of the universe, from which the bounce occurs, will be shifted to the smaller values. Also, the massive theory of quantum cosmology exhibits some other possibilities; for example, tunneling between different classically allowed regions, for cosmic evolution in the early universe epoch.