Quantum and classical cosmology in the Brans-Dicke theory

In this paper we discuss classical and quantum aspects of cosmological models in Brans-Dicke theory. First, we review cosmological bounce solution in Brans-Dicke theory that obeys energy conditions (without ghost) for a universe filled with radiative fluid. Then we quantize this classical model in a canonical way, establishing the corresponding Wheeler-DeWitt equation in the minisuperspace, and analyze the quantum solutions. When the energy conditions are violated, corresponding to the case $\omega<-3/2$, the energy is bounded from below and singularity-free solutions are found. However, in the case $\omega>-3/2$ we cannot compute the evolution of the scale factor by evaluating the expectation values because the wave function is not finite (energy spectrum is not bounded from below). But we can analyze this case using Bohmian mechanics and the de Broglie-Bohm interpretation of quantum mechanics. Using this approach, the classical and quantum results can be compared for any value of $\omega$.


Introduction
General Relativity (GR) theory is a highly successful theory that describes gravitational interactions of the universe. It has successfully survived observational tests in the solar system, and its predictions about the emission of gravitational waves by binary systems were confirmed. However, there are strong indications that the theory is incomplete. For example, it cannot account for the observed gravitational anomalies, mainly in the cosmological context, and thus it is necessary to introduce the hypothesis of dark matter and dark energy to explain cosmological observational data within GR's framework. Another reason is the initial singularity at the beginning of the universe and at the final stage of some class of stars, which are generally predicted by the GR theory. These drawbacks led to many possible extensions of the GR that can address those problems. The prototype of alternative theory of gravity is Brans-Dicke theory (BD). Historically, it is one of the most important modifications to the standard GR theory, which was introduced by Brans and Dicke [1] as a possible implementation of Mach's principle in a relativistic theory (eventually, it did not work out as we will discuss later).
It is expected that a quantum formulation of GR might solve some of its problems, especially those related to the existence of singularities. However, there are many obstacles to this quantization [2,3]. Several attempts have been made to overcome the difficulties that appear when combining the principles of the GR theory and Quantum Mechanics (QM), be it via canonical methods or other procedures, like loop quantization or string theory. A simplified approach, like the quantization of the Einstein-Hilbert action in the minisuperspace in presence of matter fields, shows that it is possible to obtain cosmological models without singularities. The construction of a quantum cosmological model encounters many problems, even when the minisuperspace restriction is used. The first one is the absence of an explicit time coordinate due to the invariance by time reparametrizations in the classical theory [4,5]. There are different ways to solve this problem. One of them is to allow the matter fields to play the role of time, which can be achieved, for example, through Schutz's description of a fluid [6]: its corresponding canonical formulation results in a Schrödinger-type equation since the conjugate momentum associated with the matter variables appears linearly in the Hamiltonian. We will employ this approach in the analysis to be exposed in the present text.
Another point of discussion is a choice of the suitable formalism to interpret the quantum theory and thus obtain specific predictions. This is a very sensible question. The usual Copenhagen interpretation is based on a probabilistic formalism, using concepts such as decoherence and a measurement mechanism through the spectral theorem. It is not ideal for a system consisting of a unique realization as it is the universe. Nevertheless, many adaptations of the Copenhagen interpretation are possible, like the Many World [7] or the Consistent Histories [8]. One alternative is the de Broglie-Bohm (dBB) interpretation of quantum mechanics [9,10], one we will address in this work. The dBB approach keeps the concept of trajectories of a given system, and a probabilistic analysis is not fundamental in this scheme [11].
In this work, we will investigate the BD theory. The reason to analyze this theory, which has been thoroughly studied in many contexts, is that it reserves some interesting and even unexpected features which deserve to be discussed in more details. Classically, cosmological models constructed from the BD theory may lead to non-singular scenarios if the energy conditions are violated, as it happens in the GR theory. Such violations of the energy conditions occur when the Brans-Dicke parameter varies in the domain of ω < −3/2. Surprisingly, it is also possible to obtain a non-singular model in the BD theory even if the energy conditions are satisfied. We will discuss this possibility using radiation as the matter content. The situation becomes more complex when we turn to quantum models: a consistent quantum model, from the point of view of the Copenhagen interpretation, is possible when ω < −3/2, since only in this case the energy is bounded from below; if the energy conditions are satisfied, the spectrum of energy is not bounded from below. However, in both cases, the analysis becomes possible even if we use the dBB interpretation of quantum mechanics. Moreover, the use of the dBB formalism allows a comparison between the classical and quantum models in the BD theory, at least to some simple configurations. This possibility will be explored in the present work, revealing many peculiarities at the classical and quantum levels.
The paper is organized as follows. In section 2 we review the classical BD theory. In section 3 we present a non-singular model with radiative fluid in BD theory that obeys the energy conditions and does not contain ghosts. Section 4 presents the quantization of the model with a Lagrangian with a non-minimally coupled scalar field and matter fluid, establishing the corresponding Wheeler-DeWitt equation in the minisuperspace. Finally, we analyze quantum solution via dBB interpretation in section 5. Our conclusions are exposed in the Section 6.

A short review of the Brans-Dicke theory
The Brans-Dicke theory was proposed as a modification of the relativistic theory of gravity to include two new features: the possibility of a dynamical coupling to gravity and Mach's principle. Dirac had suggested earlier that the gravitational coupling in Einstein's theory might not be a constant, implying it could vary with time, at least in the cosmological context [12]. This idea was further developed by Jordan [13], but the rigorous implementation was made by Brans and Dicke in their seminal 1961 article [1]. They replaced the gravitational coupling by the inverse of a scalar field φ, such that, A simple but elegant solution to implement this idea in a relativistic context is to consider in the action a non-minimal coupling between the new scalar field φ and the geometry represented by the Ricci scalar. The introduction of a kinetic term coupled by an arbitrary parameter ω complemented the theory, preserving the minimal coupling between gravity and matter and assuring the invariance by the full diffeomorphism group and the consequent conservation of the energy-momentum tensor. The presence of a long-range scalar field was connected initially with the intention to implement the Mach's principle. The most common of many formulations of Mach's principle states that the inertial property of a given body results from its interaction with all matter present in the universe. It is not easy to implement such an appealing idea. General Relativity was considered to be a Machian theory since it relates the geometry of space-time to the matter distribution, but it has non-Machian features, such as the locality of the relation between space and matter and initial conditions for the matter fields. Perhaps, these issues could be circumvented by introducing a long-range scalar field as it occurs in the BD theory. In a Machian relativistic theory, the only solution in the absence of matter should be the Minkowski one. This is not the case for GR or for BD theory. In other words, both theories failed as a proposal to implement the Mach principle. However, introducing the scalar field φ nonminimally coupled to gravity led to many phenomenological applications and intriguing results: it opens new possibilities in comparison with the GR original framework.
The BD action reads as, where ω is a coupling constant and L m is the matter term [1] and Ψ indicating generically the matter fields.
It is commonly understood that in the ω → ∞ limit BD theory coincides with GR [14,15]. Althought it is true in most situations, this statement is not valid in general. When ω 1 the field equations show that φ = O 1 ω , so we have where G N is Newton's gravitational constant and G µν is the Einstein tensor. But there are some examples [16]- [18] where exact solutions cannot be continuously deformed into the corresponding GR solutions by taking the ω → ∞ limit. In this case, the solutions decay as Moreover, there is a particle solution that admits the appropriate asymptotic behavior given by Eq. (3) but no GR limit [19]. From the action (2), the whole theory is derived. The field equations are obtained through the variation of the action (2) with respect to the metric and scalar field, The coupling with the curvature in the action action (2) can be avoided if we perform a conformal transformation on the metric g µν such that, With this, we change our frame of reference. With the reference frameg µν , one could interpret the scalar field as a matter field, and thus recover general relativity. We, however, want to remain closer to Brans and Dicke's original idea of the gravitation not being purely geometrical, and thus the scalar field does not account for a matter content. The original frame is known as Jordan's frame, and, after the conformal transformation (8), it is called Einstein's frame.
In this latter, the action becomes, where we have defined, The resulting field equations are (ignoring the tildes in the redefined geometric terms), Remark that this set of equations implies a non-usual expression for the conservation laws, which is due to the non-minimal coupling between the scalar field and matter. Thus, in the Einstein's frame, the case = 1 (ω > − 3 2 ) corresponds to an ordinary scalar field with positive energy density, while for = −1 (ω < − 3 2 ) the kinetic term of the scalar field changes sign, and it becomes a phantom field with negative energy density. In the special case ω = −3/2 with a redefinition of the matter fields, the GR theory is recovered. When the matter field is given by radiation, the non-minimal coupling between the scalar field and the matter component is also broken since T = 0 for a radiative fluid. We will use this fact later.
Currently, the limits on the parameter ω are very stringent [20]. Cosmological constraints lead to values for ω of the order of some hundreds. Binary pulsars push these bounds to dozen of thousands. Hence, the BD theory essentially becomes GR in most of the cases, as mentioned before. Despite this, BD theory continues to be relevant for several reasons. For example, for a particular value of ω = −1, the action (9) coincides with the effective string action for the dilatonic sector. In this case, the action (9) acquires some duality properties that have been explored in the Pre-Big Bang scenarios, in which there is a contracting phase in the evolution of the universe before the actual expanding phase. In fact, for ω = −1, the action (2) in the cosmological context is invariant by the transformation, where a is the scale factor. Therefore, an expanding universe can be mapped into a contracting universe. For an excellent recent review of Pre-Big Bang scenarios, we refer to [21]. Even if there is a clear discrepancy with the observational constraint on the value of ω, it is possible to take into account a possible dependence of ω on energy scales that may reconciles observations with theoretical considerations. This may be achieved, for example, by supposing that ω is a function of the scalar field φ, ω = ω(φ). Perhaps, one of the most delicate aspects concerning the pre-big bang scenario deals with the evolution of the perturbations in the primordial universe, which becomes highly discrepant with GR and, at the same time, strictly constrained by observations. In general, predictions for the evolution of perturbations depend on the transition from the contracting phase to an expanding phase. A possible solution to this problem may be found by studying perturbations behavior, taking into account the string configuration underlining the pre-big bang model since the transition may occur at the large curvature regime. Such analysis is not an easy task, technically and conceptually. Similar proposals, like the ekpyrotic scenario based on brane-world structures coming from string theories, have been developed in the literature [22].
In what concerns the primordial universe, we may also mention the extended inflationary model using the original proposal based on a cosmological constant [23] but implemented in the context of the BD theory. It is possible to obtain the transition to the radiative phase, but only if the value of the parameter ω is relatively small, in contradiction with the constraints mentioned above. In such a case, it is necessary to find a mechanism to change the value of ω during the universe's evolution since the energy scale may depend on this parameter, as already evoked above.
Multidimensional theories, when compactified to four dimensions leads to action similar to the BD, with In this case, gauge fields can appear, depending on the original configuration, with non-trivial coupling [24,25]. The proposal of modified gravity theories spiked a revival in the interest in the BD theory. Modified gravity theories were conceived mainly to address the problem of the acceleration of the universe without introducing a new exotic component in the universe called the dark energy. There are many types of the modified gravity theories. For example, in the Horndesky theories, the most general Lagrangian includes the scalar field in a non-trivial way and leads to second-order equations of motion. In fact, the BD theory can be considered the first of the modified theories, and it was formulated long before the Horndesky classification. Another class of modified gravity theories is the f (R) theories. They are based on a non-linear generalization of the Einstein-Hilbert action.
Interestingly, the f (R) theories can be recast as the BD theory, with ω = 0 and a potential term that depends on the form of the f (R) function. In general, all modified gravity theories must introduce a screening mechanism to reconcile the large-scale and small-scale constraints. According to these screening mechanisms, the supplementary degree of freedom associated with the scalar field does not propagate in a dense region (compared to the averaged cosmological density). The BD theory, when reformulated in the Einstein's frame, gives the prototype of the chameleon mechanism, one of the most important screening mechanisms: the conformal transformation used to reformulate the BD theory in the Einstein's frame introduces a non-minimal coupling between matter and scalar field, making the mass of the scalar field depend on the density of the medium, as it can be verified introducing a potential term in the set of equations (13,14).
Intriguing results on the quantization of the BD theory adds to the relevance of the theory. Analysis on quantum cosmological scenarios using the BD theory was made in Refs. [26,27]. In Ref. [26], the authors considered the Einstein's frame and used the WKB formalism to recover the notion of time. The result reveals a curious behavior: the quantum effects become relevant in a late phase of the expansion of the universe, contrary to expectation. Due to this property, the initial singularity existing in the FLRW models can not be avoided. In [27], the study was extended to the Jordan's frame through the de Broglie-Bohm formalism to compute the quantum evolution of the universe. In this case, the opposite scenario was found: the quantum effects become important mainly in the early universe, and the initial singularity can be avoided.
We must also mention the study of the resulting quantum Schrödinger-like equation from the point of view of its self-adjoint properties [28], which are essential to employ the spectral theorem. The results indicated that for the case ω < −3/2, which corresponds to a phantom scalar field, it is possible to recover the self-adjointness of the quantum operator. We will discuss this in more detail in the forthcoming sections.

Bouncing solutions and the energy conditions
The most common solution to avoiding a singularity in cosmological models is the introduction of exotic types of matter fields, for example, a scalar field with negative energy density. On the other hand, to obtain a bouncing solution in classical General Relativity, violation of the energy conditions is required. In this section, after briefly review the bouncing scenarios, we present a non-singular model with radiative fluid in BD theory that obeys the energy conditions and does not contain ghosts [29]. To do so, we first briefly analyze the solutions determined by Gurevich et al [30] for the cosmological isotropic and homogeneous flat universe with a perfect fluid with an equation of state p = nρ, where the parameter n is given by 0 ≤ n ≤ 1. Then we focus on the analysis of these solutions in the case of radiative fluid.
The study of bouncing models is motivated by the search for cosmological solutions without singularities. Bouncing models have been widely studied to solve the initial-singularity problem and as an addition or alternative to inflation to describe the primordial universe because they can explain, in their own way, the horizon and flatness problems and justify the power spectrum of primordial cosmological perturbations inferred by observations [31]- [33].
In such models, an initial singularity is replaced with a bounce-a smooth transition from contraction to expansion. To obtain a bounce, one needs to change the value of the Hubble parameter H ≡ȧ (t) a(t) , which appears to be negative during the contracting phase, to a positive value for the following expanding phase. One of the options to switch the sign of the Hubble function lies within GR. It usually requires the violation of the null energy condition (NEC) [34] ρ + p ≥ 0, where ρ is energy density, and p is pressure as usual. In most cases, the energy conditions reflect the nature of matter fields, "ordinary" (attractive effects) or "exotic" (repulsive effects). However, as we will see later, the bounce can be achieved through some non-standard coupling between the fields existing in a given theory. At the early stages of the expansion of the universe, the curvature is not essential since the matter components are more relevant for the dynamic than the curvature term due to their dependence on the scale factor, and we can restrict our further analysis to the quasi-euclidian variant of the isotropic model. In this case, for a flat FLRW metric with N 2 = 1, the classical field equations (6) and (7) reduce to ρ + 3ȧ a (ρ + p) = 0 .
The general solutions in [30] are obtained for ω < − 3 2 and ω > − 3 2 . In the first case, there is violation of the energy conditions for the scalar field in Einstein's frame. In the latter case the energy conditions for the scalar field are satisfied as we will discuss in this section.
The solutions for the scale factor and scalar field for ω < − 3 2 are where σ = 1 + ω(1 − n), 2A = (1 − 3n) + 3σ(1 − n), a 0 is an arbitrary constant, and θ − > θ + are integration constants. The time coordinate θ is connected with the cosmic time t by definition, When θ → ∞, the scale factor a does not vanish. The infinite contraction has a minimum a min , and it is followed by the expansion. Thus, this model admits a cosmological bounce.
Henceforth, we shall focus on the universe filled with radiative fluid, represented by the equation of state p = 1 3 ρ. In this case, the solutions for ω > − 3 2 are given by the following expressions: where η is the conformal time and η ± are constants such that η + > η − . In Figure  1, we plot the scale factor and scalar field for the lower sign. The solutions for In the case of the lower sign in the equations (28)-(29), we obtain bounce solutions for − 3 2 < ω < 0. Nevertheless, there is a curvature singularity at η = η + for − 4 3 < ω < 0, even if the scale factor diverges at this point. But for − 3 2 < ω ≤ − 4 3 bounce solutions are always regular without curvature singularity. In this case, there are two possible scenarios of the evolution of the universe due to the time reversal invariance. In the first one the universe begins at η = η + , with a → ∞ and an infinite value for the gravitational coupling (φ = 0). It evolves to the other asymptotic limit with a → ∞, although with φ constant and finite. The second option is the reversal behaviour of the first one for −∞ < η < −η + . In both cases, the cosmic times varies as −∞ < t < ∞. The dual solution in the Einstein's frame for − 3 2 < ω ≤ − 4 3 is given by with b = φ 1/2 a and contains an initial singularity. This can be considered as a specific case of "conformal continuation" in the scalar-tensor gravity proposed in [35]. Let us show that the energy conditions for the scalar field are satisfied for ω > − 3 2 . In general, in order to have a bounce solution, violation of the energy conditions is required. Using the Friedmann and Raychaudhuri equations we represent the strong and null energy conditions in GR as We reformulate the BD theory in the Einstein's frame so that it would be possible to use the energy condition in this form. One can see that both energy conditions are satisfied as long as ω > − 3 2 . This is in agreement with the fact that in the Einstein's frame the cosmological models are singular unless ω < − 3 2 . However, in the original Jordan's frame there are non-singular models for − 3 2 < ω < − 4 3 . But, in this range, the scalar field and the matter component obey the energy condition. The effects that lead to the absence of the singularity come from the non-minimal coupling. In the Figure 2, we present the effective energy condition, defined in the left-hand side of Equations (33)- (34), considering the effects of the non-minimal coupling. If we analyze only the left-hand side in the Equations (33)- (34), the effects of the interaction due to the non-minimal coupling are included, and the energy conditions can be violated even if the matter terms do not violate them. For more details, see Ref. [29]. It is important to observe that only in the case of radiative fluid it is possible to obtain a model without singularity preserving the energy conditions, at least in the BD theory. It is true also for the model with flat spatial sections. For a non-flat universe, one can obtain a singularity-free scenario even in General Relativity if the strong energy condition-but not necessarily the null energy condition-is violated.

Canonical quantization of the BD theory
In the early Universe, when the cosmo is compressed in a Planck scale, quantum effects become relevant, and the quantization of the classical models may be a necessity. Quantum cosmology stands on this principle, considering a wave description of the universe that satisfy the Hamiltonian constraint of the theory. In this section, we will investigate the canonical quantization of the BD theory in the minisuperspace, considering the FLRW metric (18). We shall quantize the Hamiltonian constraint H tot ≈ 0 ( where "≈" means weak equality. See, for example, [36]) to obtain the Wheeler-DeWitt EquationĤ tot Ψ = 0, with Ψ being the wave function of the universe, and using the canonical operators ( = 1), The total Hamiltonian H = H tot is formed by the gravitational part and the Hamiltonian of the matter content, H = H G + H M . To find H M , let us consider a model of early-universe filled with a radiative fluid. Using Schutz formalism [6], the super-Hamiltonian of the fluid is given by where T is directly related to the entropy of the fluid [37]. With this, the canonical quantization of the total Hamiltonian in the original Jordan's frame is, where p, q are ordering factors for the quantization of the momenta squared. For more details on the computation of (37), see Ref. [37]. Notice, however, that here we use a different sign convention. The calculation is similar in Einstein's frame, considering the transformation (8). In this case, the canonical quantization results in Here, again, we have thatp andq are ordering factors. Notice that Equations (37) and (38) are Schrödinger-like, that is, if we consider the matter field playing the role of time. We can, therefore, treat it as quantum system to analyze this cosmological scenario in the early-universe. The first step is to verify the conditions for the effective Hamiltonian operators of these models to be self-adjoint. It is known [37] that the Hamiltonian operator of this quantized BD model with a radiative fluid can be self-adjoint only if q =q = 1. On the other hand, we have the quantum equivalence between Jordan's and Einstein's frames [38] if, and only if, p =p = 1, and thus we can use Einstein's frame. Now, in Einstein's frame, we choose the coordinate instead of φ. With this, the relation between the scale factor in Jordan's and Einstein's frames is, In terms of the new variable σ, the Schrödinger equation (39) is, The measure of the Hilbert space is such that The regular solution for the equation (42) is, where k is a separation constant. There is also another solution written in terms of the Bessel function J −ν (x), which is not regular at the origin, at least when = −1. For this reason, we will disregard it at the moment. We will choose = −1, for this case the Hamiltonian operator is bounded from bellow and it is essentially self-adjoint for −1 ≤p ≤ 3 (see [37]), which is our case. Thus, ν = |k|. Let us choose the coefficient A(k, E) such that the wave packet becomes, where N is a normalization factor, and with γ a positive real parameter. The integration in x gives us [39], To find the normalization factor N, we should remember that, and thus, 4α * e −ik σ dk bdb dσ Integrating over σ, k and k, and defining u = b/2|α|, we obtain where Φ(x) is the error function, and Since −1 < Φ(x) < 1 and ∞ 0 u 2 e −γu 2 is strictly positive, we have Therefore, N 2 is positive, as expected.
We can now calculate the expected values of the scale factor b and scalar field σ. For the scale factor, Similarly to the calculation of the normalization factor, integrating in σ and k , we have, where, Since |α| = γ 2 + T 2 , and defining which is strictly positive, the expected value of b is By definition, γ > 0. Therefore, b > 0, that is, there is no singularity at T = 0. Moreover, for T γ, at late times, b → T .
The integral over k is zero because the function is odd in respect with this parameter. Thus, This does not mean, however, that the scalar field φ is a constant, since the expectation value allows fluctuations. Notice that we have a symmetrical bouncing in both frames because of Equation (41).

Analysis of the solution via the de Broglie-Bohm approach
In the previous section, we introduced a quantum model of the BD theory, identifying the universe with a wave function that obeys a Schrödinger-like equation. Thus, we proceeded with the usual methods of quantum mechanics to analyze the behavior of the scale factor, which is directly connected with the volume of the universe. However, in the case of quantum cosmology, we must give up on the usual Copenhagen interpretation of QM, since it requires an external observer to cause the wave function to collapse into one state. An alternative is the many-world interpretation, which considers every state of the wave function as real, existing in parallel with each other. Our universe, therefore, is one of many. This interpretation does not require the collapse of the wave function, and so we can investigate different states separately from the wave packet. The de Broglie-Bohm (dBB) interpretation, on the other hand, does not consider a wave description of the universe at all. Instead, it is a dynamical theory in which real trajectories can be obtained in the configuration space of the quantum system. Those trajectories are observer-independent, and, therefore, the de Broglie-Bohm interpretation does not rely on any collapse mechanism. This interpretation was formulated as a causal alternative to the probabilistic Copenhagen interpretation still in the early years of the quantum era [40,41], and it has replicated most of classic results of quantum mechanics. In the dBB interpretation, observers are also described by quantum operators to be applied on a wave function that satisfies the Schrödinger equation. The wave function dictates the equations of motion to find the trajectories. In a Bohmian mechanics, Ψ(x i , t) is decomposed as, and the probability density of the trajectory is given by ρ(x i , t) = |Ψ(x i , t)| = R 2 (x i , t). The Bohmian trajectories are realized through the momenta defined as Since it is not probabilistic, the observers in the dBB do not necessarily need to be represented by self-adjoint operators. This ontological interpretation of the quantum theory has its critics, who argue that other interpretations may be better-suited [42], or take issue with the so-called hidden variables which determine the behavior of the trajectories in a many-body configuration [43]. Still, because of its distinct characteristic, the dBB interpretation is a suitable candidate to be applied in the quantization of cosmological scenarios [44,45], where the wave function becomes a guide to the possible observer-independent evolution of the universe. Notice that, contrary to the Copenhagen interpretation, the dBB interpretation of quantum mechanics does not require necessarily to work in the Hilbert space.
Let us compare these two approaches in particular cases. In the previous section we have already made the computation using the expectation values finding, quite generally, that the expectation value for the scalar field is zero and the expectation value for the scale factor is the same as found in the corresponding case in GR. We will now analyze the predictions for the evolution of the universe using the dBB formulation. We will work initially in the Einstein's frame. Remember the usual solution to the Eq. (42) given in Eq. (44): The construction of the wave packet exposed in the previous section is not convenient here. The reason is that the integration is made in the interval −∞ < k < ∞, but the terms coming from the order of the Bessel functions depend on |k|. Hence, the integration on the whole interval does not allow to obtain a tractable form for the wave packet in view of using the de Broglie-Bohm approach. Let us consider the general wave packet given by, with the definition x = √ E and ν = |k|. The integration in x leads to, With this, we investigate the behaviors of the scale factor and the scalar field in some special cases, and compare them with the results of previous section.

The scalar field is absent
Let us consider the case where the factor A is given by a delta function. If we choose A(k) = δ(k), which means the world where ν = 0, the contribution of the scalar field vanishes and the wave function is such that, Notice that this does not mean we recover General Relativity. This would be one world of the many in a BD theory of gravitation. Considering Eq. (48), it is straightforward to verify that the norm of the wave function (61) is finite and time-independent. Similarly to the calculation in last section, we can compute the expected value of the scale factor.
where N here is the normalization factor of the wave function (61), naturally. Therefore, in this world, a bounce occurs, and for late times, that is, when For the Bohmian trajectories (57) we use the phase of the wavefunction (61) which is, in this case, Remembering that the conjugate momentum is p b =ḃ/2, we obtain for the Bohmian trajectories,ḃ The solution of this differential equation is b = b 0 γ 2 + T 2 , which is the same result as before using expectation values. In a Bohmian analysis, this universe also has a bounce. The results exposed in the previous section show quite generically that wave packets with constant finite norm lead to a zero expectation value for σ. From now on we explore possibilities where we circumvent this restriction but at the price to have wave packets that are not finite, what is not an obstacle when dBB formulation is used, see Ref. [11].

A single scalar mode
Let us consider now the superposition function such that A(k) = δ(k − k 0 ), where k 0 is a positive constant. This implies to consider a single scalar mode behaving like a plane wave. The wave function reads as, with ν 0 = √ − k 0 . We will analyze the Bohmian scenario for both = 1 and = −1.
Let us start with = −1. In this case, the phase of the wave function becomes The presence of the term k 0 σ changes the previous analysis, since now the variable σ is featured in the wave function. The guidance equations become: and their solutions are, Notice the solutions (75) and (76) are non-singular only if σ 0 > π 2k 0 and σ 0 > 0. The classical solution is recovered asymptotically. This fact reveals again the very peculiar features when the energy conditions are satisfied.
We can compute the quantum potential that results from the modified Hamilton-Jacobi equation in the dBB formulation of quantum mechanics [11], which is given by, The Laplacian operator is defined in the minisuperspace with variables b and σ. For the cases studied in this subsection, the result is the expected: the quantum potential reaches its maximum value at the bounce, and decreases to zero asymptotically where the classical solution are recovered.
In performing the study above, we have used the Bessel function given in (44). It can be explicitly verified that if we have used the Bessel function of negative order, which is not regular at the origin, the results would be the same with a reversal of the time, T → −T .

Multiple scalar modes
We can combine different scalar modes. One example can be achieved by the combination, with η = ±1. With a similar computation, we obtain, for the positive sign, cos(k 0 σ) instead of e ik0 and for negative sign sin(k 0 σ). In both of these cases, the scalar field is not present in the phase of the wave function, and we recover the same solutions already given in the case the scalar field is absent.

Conclusions
The Brans-Dicke theory is one of the oldest proposals of a modification to the theory of general relativity. Although it has been studied for over sixty years, the Brans-Dicke theory continues to reveal intriguing aspects. Some of those were discussed in the present work concerning classical and quantum scenarios for the early universe. First of all, it is possible to obtain a singularity-free cosmological solution if the Brans-Dicke parameter ω varies as −3/2 < ω < −4/3. In this range, the energy conditions are satisfied in the Einstein's frame: the avoidance of the singularity is driven by the non-minimal coupling. When we turn to the quantum scenario in the minisuperspace, other curious features appear. The energy is bounded only for ω < −3/2, that is, when the energy condition is violated in the Einstein's frame. On the other hand, in the case of ω > −3/2 the energy conditions are satisfied, but the energy is not bounded from below, so it becomes problematic to employ the usual interpretation scheme based on the Copenhagen formulation of quantum mechanics: the wave function is not finite anymore.
This issue motivated us to consider the de Broglie-Bohm interpretation of quantum mechanics. The previous result are found again if ω < −3/2. However, when ω > −3/2, we have either a singular or non-singular solution. This implies that in the interval − 3 2 < ω < − 4 3 the classical model displays singularity-free scenarios, while the quantum models may display either singular or non-singular solution. This result raises the following question: can such results occur in more general modified gravity theories belonging to the Horndesky class?
We add to the previous discussion some additional remarks. First of all, it must be verified to what extent the results reported here depend on the wave packet construction and the choice of the time variable. On the other hand, the problems of convergence of the wave function, when the energy conditions are obeyed, is somehow general in the presence of scalar fields since it leads to a hyperbolic signature in the Hamiltonian and, consequently, in the Schrödingertype equation. However, such a well-known fact acquires new features that were briefly described above. The analysis displayed in the present work shows that the well-studied Brans-Dicke cosmological models present interesting properties at classical and quantum levels as well, no matter which interpretation is chosen.