1. Introduction
At present, we are living in a cosmological era where, on the one hand, we have several cosmological data that include the scalar spectral index, the tensor-to-scalar ratio describing the early stage of the universe, as well as the late-time equation of state parameter, the Om(z) parameter, in regard to the dark energy era of the universe. However, on the other hand, modern cosmology is still riddled with questions on whether the universe started its expansion from a Big Bang singularity or from a non-singular bouncing scenario. Inflation is one of the cosmological scenarios that successfully describes the early stage of the universe; in particular, it solves the horizon and flatness problems, and most importantly, it predicts an almost scale-invariant curvature perturbation power spectrum that is rather consistent with the recent Planck data [
1,
2,
3,
4,
5]. However, extrapolating backwards in time, the inflationary scenario implies an initial singularity of the universe known as the Big Bang singularity, where, due to geodesic incompleteness, the spacetime curvature diverges at the point of singularity. Bouncing cosmology is one of the alternatives of inflation that can generate a scale-invariant curvature power spectrum. In addition, the bounce scenario leads to a singular-free evolution of the universe [
6,
7,
8,
9,
10,
11,
12,
13,
14,
15,
16,
17,
18,
19,
20,
21,
22,
23,
24,
25,
26,
27,
28,
29,
30,
31,
32,
33,
34,
35,
36,
37,
38,
39,
40,
41,
42,
43]. In the present work, we are interested in bounce cosmology that appears with certain features, which we will explain.
Among the various bounce models proposed so far, the matter bounce scenario (MBS) has gained a lot of popularity due to the fact that it produces a scale invariant primordial power spectrum and also leads to a matter-dominated universe at the late expanding phase [
11,
18,
20,
44,
45,
46]. However, the MBS is plagued by some problems such as: (1) the tensor-to-scalar ratio in the MBS generally becomes too large to be consistent with the observational data [
47]; (2) the evolution of the universe during the contracting stage becomes unstable due to the growth of anisotropies, which leads to the BKL instability in the MBS [
48]; (3) the scale factor describing the MBS clearly depicts the universe undergoing a late-time deceleration phase, which is not consistent with the dark energy observation [
49,
50,
51]. Here, it is worth mentioning that such problems can be well-resolved by suitable modifications to the matter bounce scenario, and the success of the bounce cosmology becomes quite illuminating. In particular, the authors of [
52,
53] proposed an extended bounce scenario in the Lagrange multiplier F(R) theory or in the ghost-free Gauss–Bonnet theory, where the first problem is resolved; however, the second and third problems persist. The article in [
21,
54] proposed an ekpyrotic bounce scenario which is free from the BKL instability and predicts an almost scale-invariant curvature power spectrum, but the models are unable to explain the dark energy issue. Recently, a smooth unified scenario from an ekpyrotic bounce to the dark energy era has been proposed in [
55], which stands to be a viable bounce in respect to the Planck data and concomitantly resolves the BKL instability.
In addition to the Big Bang or the initial singularity, there are other types of finite-time singularities in cosmology, which were classified and studied for the first time in [
56]. The most severe is the Big Rip or the Type-I singularity. However, there are three more of such finite-time singularities—the Type-II, Type-III, and Type-IV singularities. For studies on these types of singularities, see [
57,
58,
59,
60,
61,
62,
63,
64,
65,
66,
67,
68,
69,
70]. In the case of the Type-I, Type-II, or Type-III singularities, the effective energy density and/or the effective pressure of the universe diverge at the point of singularity. However, on the contrary, in the Type-IV singularity, the effective energy density and the effective pressure of the universe remain finite when the singularity occurs. Thus, the Type-IV singularity is the mildest one among the finite-time singularities. In particular, unlike the Big Rip singularity, geodesic incompleteness does not appear in the Type-IV singularity. As a result, the universe can smoothly pass through a Type-IV singularity (if any), and moreover, the Type-IV does not lead to catastrophic events at observable quantities.Therefore, it is possible that the universe faced a Type-IV singularity during its evolution in the past, through which it smoothly passed without any geodesic incompleteness. However, the presence of a Type-IV singularity may have had a significant influence on the evolution of the universe as well as on the generation era of the primordial perturbation modes, as the singularity could globally affect the Hubble parameter. Therefore, in the realm of bouncing cosmology, the important question that immediately arises is the following:
We will address these questions in the present work. For the gravity theory, we will consider the well-formulated
theory, which turns out to be ghost-free with the presence of the Lagrange multiplier in the gravitational action, as developed in [
71] (
R is the Ricci scalar, and
is the Gauss–Bonnet term). The cosmology of the
gravity from various perspectives has been discussed in [
72,
73,
74,
75,
76,
77,
78,
79,
80]. The GB coupling function in the present context will be considered in such a way that it will satisfy a constraint equation such as
(where
is the GB coupling function, and
is the Hubble parameter of the universe). Owing to such conditions on
, the speed of the gravitational wave approaches unity, and the model becomes compatible with the event GW170817. In such a scenario, our main aim will be to examine how the presence of a Type-IV singularity (a mild singularity) can influence the dynamics of a bouncing universe. At this point, we would like to mention that some of our authors have studied bounce cosmology with a Type-IV singularity in the context of F(R) gravity in [
64], particularly the authors of [
64] showed that the presence of a Type-IV singularity destroys the viability of a bounce scenario with respect to the Planck data. However, in the present work, we will consider a ghost-free Gauss–Bonnet theory of gravity, and more importantly, we will address how a bounce scenario that appears with a Type-IV singularity can be made viable with the recent observational data. These make the present work essentially different from [
64]. We will consider two different bouncing scenarios in the current work: (1) In the first scenario, the Type-IV singularity “globally” affects the spacetime (where the term “global” means that although the Type-IV singularity occurs at a finite time
, it can have an effect even on the asymptotic evolution of the Hubble parameter in the distant past as well as in the distant future). (2) In the second scenario, the Type-IV singularity affects the spacetime “locally”, at around the time when it occurs. These two scenarios are qualitatively different, particularly in the former scenario where the Type-IV singularity globally affects the spacetime, the comoving Hubble radius asymptotically goes to zero on both sides of the bounce, and thus the perturbation modes generate near the bounce; in the second bounce scenario, where the Type-IV singularity locally affects the spacetime, the comoving Hubble radius asymptotically diverges to infinity in the distant past, and consequently, the primordial perturbation modes generated far away from the bounce in the deep contracting phase. We will perform the scalar and tensor perturbations for these two scenarios, which are quite different when compared to one another. The possible implications will be discussed accordingly.
This paper is organized as follows: in
Section 2, we briefly discuss the ghost-free
gravity compatible with GW170817. Then, in
Section 3 and
Section 4, we will describe two different bounce scenarios, as mentioned above. Consequently, we will examine the possible effects of the Type-IV singularity on the bouncing dynamics. The paper ends with some conclusions in
Section 5. We would like to clarify the notations and conventions that we will adopt in the subsequent sections. We will work with natural units, and the metric signature will be mostly positive, i.e.,
. A suffix “b” with a particular quantity will refer to the quantity at the instant of bounce, and
will denote the time when the Type-IV singularity occurs. An over prime with some arguments will indicate the derivative with respect to the argument; otherwise, an over prime will represent
. Moreover, an over dot will denote the derivative with respect to the cosmic time.
2. Brief Review of Ghost-Free Gravity Compatible with the GW170817 Event
In this section, we shall recall the essential features of the ghost-free
gravity theory developed in [
71]. We consider
, which, owing to the presence of
, contains ghosts with respect to perturbations of the spacetime metric. However, the ghost modes may be eliminated by introducing a Lagrange multiplier
in the standard
gravity action [
71], which leads to a ghost-free action, as follows:
where
is a constant with mass dimension [+1]. By varying the action with respect to the Lagrange multiplier
, one obtains the following constraint equation:
The kinetic term of
seems to be a constant, and thus it can be absorbed in the scalar potential as:
and as a result, the action of Equation (
1) can be equivalently expressed as follows:
The equations of motion for
and
from the action in (
4) take the following forms:
where
is given by the following:
while having
in mind. The trace of Equation (
6) (i.e., multiplied with
) becomes the following equation:
and solving Equation (
7) with respect to
yields the equation below:
The spatially flat Friedmann–Robertson–Walker (FRW) metric ansatz will fulfill our purpose in the present context, and hence we consider the line element as follows:
Assuming that the functions
and
are homogeneous in cosmic time, and also that
, then Equation (
2) immediately leads to the following simple solution:
where
is a constant and will be identified with the cosmic instance of bounce. Consequently, the
and
components of Equation (
6) result to the following equations:
and furthermore, from Equation (
5), we obtain the following:
It may be worth mentioning that the three equations above are not independent; in fact, by combining Equations (
13) and (
14), one can obtain Equation (
12). It is evident that Equation (
11) is an algebraic equation with respect to
, and that in solving this equation, we obtain the following:
Similarly, the scalar potential
can be obtained by solving Equation (
12) and is given by:
Recall that
, and thus
is equal to the following:
The functional form of the Lagrange multiplier from Equation (
14) becomes:
Therefore, a form of the Hubble parameter and the Gauss–Bonnet coupling function in turn fix the scalar potential and the Lagrange multiplier. We will consider the Hubble parameter in such a way that it leads to a bounce scenario with a Type-IV singularity. The resulting theory with the Lagrangian (
4) is a form of the scalar Einstein–Gauss–Bonnet gravity combined with a Lagrange multiplier, through which it is well-known that the speed of the gravitational waves can be expressed as follows [
81,
82,
83,
84,
85]:
with
being the Hubble parameter. Clearly, the
is different from unity, and the deviation of the
from unity is controlled by the Gauss–Bonnet coupling function. However, in regard to the GW170817 event (which validates the fact that the gravitational waves have the same propagation speed as the electromagnetic waves, which is unity in natural units), we consider the Gauss–Bonnet coupling function in the present context in such a way that it results in
. For this purpose, we need to consider the coupling function in the following fashion [
86]:
Therefore, in order to be compatible with GW170817, we will consider such Gauss–Bonnet coupling functions which satisfy Equation (
19). At this stage, it is worth mentioning that the constraint above on
indeed fits the original equations of motion. This is clear from the fact that there are two independent equations, namely the
component of the gravitational equation and the equation for
; however, there are two unknown functions (
(t),
) to yet to be determined. Finally, owing to the condition of Equation (
19), the scalar potential and the Lagrange multiplier can be simplified from Equations (
16) and (
17), respectively, and are given by the following:
Thus, as a whole, Equations (
19)–(
21) are the main equations which, with a suitable form of
, lead to the corresponding forms of
,
, and
.
3. Realization of a Bounce with a Type-IV Singularity
As mentioned in the introductory section, we are interested in examining the possible effects of a mild singularity, particularly those of a Type-IV singularity, in a bounce scenario. Therefore, in this paper, we intend to realize a bounce scenario in presence of a Type-IV singularity, and depending on the time of the singularity, three different cases may arise: (1) when the singularity occurs before the bounce happens, (2) when the singularity occurs after the bounce, and (3) the case where the Type-IV singularity occurs at the instant of the bounce. The scale factor we consider is given by the following equation:
The scale factor is taken as a product of two factors—
and
—where
ensures the occurrence of a finite-time singularity.
is sufficient for leading to a non-singular bouncing universe, where the bounce occurs at
. However, for
, the bounce scenario becomes free from any finite-time singularity. Thus, due to our particular interest in the present context, i.e., to examine the possible effects of a finite-time singularity in an otherwise non-singular bouncing dynamics, we consider the scale factor as presented in Equation (
22), where
is multiplied by
. The presence of
results in a finite-time singularity at
in the bouncing dynamics, which we will discuss after Equation (
24). We will show that the presence of
does not harm the bouncing character of the universe; however, it slightly shifts the bouncing time from
to a negative or a positive time, depending on whether
or
, respectively. Moreover, the scale factor of Equation (
22) leads to an asymmetric bounce scenario (as
). The parameter
present in
is considered to have the form similar to
(with
p and
q being positive integers), so that the term
acquires positive values during the entire cosmic range (we take the root:
, so that the scale factor and the corresponding Hubble parameter are real-valued functions). Otherwise, for
, the term
becomes negative during
, and consequently,
in the distant past, which is not healthy for a non-singular cosmological evolution of the universe. Thus, we take
in the subsequent calculation. This is also important for obtaining a bounce in the present context, as we will demonstrate below with Equation (
27).
Equation (
22) immediately leads to the Hubble parameter and its first derivative (with respect to the cosmic time) as:
and
respectively. The above expression of
refers to different types of finite-time singularity depending on the values of
. In particular:
For , a Type-I singularity appears at , i.e., the scale factor, the Hubble parameter, and the derivative(s) of the Hubble parameter simultaneously diverge at . The divergence of the Hubble parameter and its first derivative indicates the divergence of the effective energy density and the effective pressure, respectively.
For , a Type-III singularity occurs at , i.e., the scale factor tends to a finite value, while the Hubble parameter and its derivative(s) diverge at .
For , a Type-II singularity appears at . In this case, the scale factor and tend to a finite value, while (and also the higher derivatives) diverges at .
For and a non-integer, a Type-IV singularity appears at , in which case the scale factors and tend to a finite value . However, the higher derivatives of the Hubble parameter diverge at the singularity point. Clearly, in the case of the Type-IV singularity, the effective energy density and the effective pressure are finite at the time of the singularity.
Therefore, among these finite-time singularities, the Type-IV singularity is the mildest singularity, and in the present context, we want to investigate the possible effects of such a mild singularity on an otherwise non-singular bounce cosmology. For this purpose, we take
. Here, it is worth mentioning that the parameter
satisfies the following conditions—(a)
should be greater than unity in order to have a Type-IV singularity at a finite time, (b)
should be of a form similar to
, where
p and
q are positive integers, and (c) we consider the real root of
(note that
has two complex branches and one real negative), in particular,
. The last two conditions ensure the requirement that both the scale factor and the Hubble parameter (see Equations (
22) and (
23)) are real-valued functions during the entire cosmic range in the present context. In particular, owing to such form of
, the term
present in the scale factor becomes a real-positive during
and is a real-negative during
. Therefore, complex values for the scale factor are avoided with the aforementioned choice of
. In this regard, we would like to mention that a similar kind of scale factor has been considered in the earlier literature (see [
55,
64,
65]).
We would now like to mention that the cosmological evolution predicted from the scale factor of Equation (
22) can indeed be realized in the present context of the Gauss–Bonnet (GB) theory of gravity with suitable forms of the scalar field potential and the GB coupling function. By integrating both sides of Equation (
19), we get the following equation:
where
is a constant having a mass dimension
. By using the above expression of
, we evaluate the scalar field potential (
) from Equation (
20) and the Lagrange multiplier (
) from Equation (
21):
where
,
, and
are given above in Equations (
22)–(
24), respectively. Since we are dealing with
in order to have a Type-IV singularity at a finite time, the Hubble parameter and its first derivative are regular even at the time of singularity, which in turn ensures the regular behavior of both the scalar field potential and the Lagrange multiplier during the entire cosmic evolution of the universe. It may be noted that both the scalar field potential and the Lagrange multiplier contain
, and thus their derivatives diverge at the Type-IV singularity at
.
In a bouncing universe, the universe initially contracts where the Hubble parameter is negative, and after it bounces off, the universe enters an expanding phase when the Hubble parameter becomes positive. Therefore, at the bouncing point, the Hubble parameter satisfies the following conditions:
and
, respectively. Next, we will examine whether the scale factor of Equation (
22) leads to a bouncing universe, and for this purpose, we separately consider
,
, and
.
For : Here, we take
. Therefore, in the cosmic regime
, both terms of
in Equation (
23) are negative. Thus, there is no possibility to have
(or equivalently, a bounce) in this regime. However, during
, the first term of
is negative, while the second term containing
becomes positive. Thus, there is a possibility for
during
, which may lead to a bouncing universe. In particular, during the negative cosmic time, i.e., for
, the Hubble parameter from Equation (
23) can be written as follows:
where
and
are the first and the second terms, respectively, sitting on the right-hand side of the expression above. The evolutions of
and
and their comparisons are shown in the following
Table 1.
Therefore,
starts from zero in the distant past, and having a maximum within
, it again reaches zero at
. On the other hand,
seems to be an increasing function during
, in particular,
in the distant past, and experiences a zero crossing at
. These clearly argue that there exists a time (say,
) within
when
, i.e., from Equation (
27):
Thus, as a whole:
, or equivalently , at ;
, or equivalently , at ;
, or equivalently , at .
This implies that
(
) is the bouncing time that can be determined from Equation (
28). A closed solution of
may not be possible, however, we numerically obtain the solution for
from Equation (
28) for a suitable set of parameter values. In particular, for
,
,
,
,
, and
(where ’By’ stands for Billion years), we obtain
. In general, Equation (
28) clearly indicates that for
, the Type-IV singularity occurs before the instant of the bounce, i.e.,
or
.
For : Performing the same procedure we followed in the previous case, we argue that for
, the scale factor of Equation (
22) leads to a bouncing universe, and that the bounce happens within
. If the instant of bounce is denoted by
, then from Equation (
23), we obtain the following equation:
Once again, one may numerically solve
from the above algebraic equation with a suitable set of parameter values. In particular, for
,
,
,
,
, and
, the bounce time is at
. From Equation (
29), it is clear that for
, the Type-IV singularity occurs after the bounce, i.e.,
.
For : The expression of
in Equation (
23) reveals that for
, the universe experiences a bounce at
. Thus, the condition
leads to the fact that the Type-IV singularity occurs at the instant of the bounce.
Therefore, for all the three cases, the scale factor we considered in Equation (
22) leads to a bouncing universe along with a Type-IV singularity for
. It is worth mentioning that the presence of the Type-IV singularity affects the time of the bounce. In particular, (1) if the singularity occurs at
, then the bounce happens at some negative instant of time, and moreover, it suggests that the singularity occurs before the bounce happens. (2) The case of
leads to the bounce occurring at some positive time, a nd consequently, the singularity takes place after the bounce. Finally, (3) for
, the bounce and the Type-IV singularity appear at the same instant of time
.
The comoving Hubble radius is defined by
, which, by definition, diverges at the bounce. By using Equations (
22) and (
23), we give the plots for
vs.
t for the above three cases (see
Figure 1 and
Figure 2).
The left plot of
Figure 1 corresponds to the case
, which clearly depicts that
diverges or equivalently, that the bounce happens at a
. Moreover, the Type-IV singularity is found to occur before the bounce (see the discussion in the caption of the figure). The other plots of
are also in accordance with the arguments mentioned above. One important point to consider is that irrespective of whether
or
or
, the comoving Hubble radius asymptotically goes to zero on both sides of the bounce. Currently,
in the distant past or in the distant future comes in the following form:
which is independent of
and clearly demonstrates that
asymptotically vanishes at
. Moreover, this argument holds for any value of the parameter
n. In a bouncing universe, the primordial perturbation modes generate either near the bounce or long before the bounce, depending on the asymptotic evolution of the comoving Hubble radius. In the present context, we find that
at
, which suggests that the perturbation modes generate near the bounce when all the modes lie within the sub-Hubble regime. Here, it may be mentioned that for
(i.e., in absence of the Type-IV singularity), the scale factor is described by
, and consequently,
asymptotically diverges to infinity (for
). Therefore, it turns out that for
, the primordial perturbation modes generate far away from the bounce in the deep contracting phase, unlike the case where
. In fact, for
, the presence of the Type-IV singularity globally affects the evolution of the Hubble parameter (compared to the case where the singularity is absent), which in turn leads to the comoving Hubble radius tending to zero at
and makes the generation era of the perturbation modes occur near the bounce.
3.1. Cosmological Perturbation
Since the perturbation modes generate and cross the horizon near the bounce, we are interested in examining the evolution of scalar and tensor perturbations near the bouncing phase. For this purpose, the useful quantities are the scale factor and the Hubble parameter near
(recall that
is the time of the bounce), and they are given by the following:
where we use the Taylor series expansion around
. Hence, by using Equations (
22) and (
23), we get:
The Gauss–Bonnet coupling function in the present context obeys
from the requirement of compatibility with the GW170817 event. By integrating
with respect to the cosmic time, one gets the following equation:
where the integration constant is taken as
from the dimensional analysis. The above forms of
and
along with the solution of
lead to the scalar field potential and the Lagrange multiplier from Equation (
21), which are expressed as:
respectively. In the context of the Gauss–Bonnet theory of gravity, we may introduce an effective potential for the scalar field as follows:
which, due to Equation (
31) along with
, can be expressed around
as:
We will use the above expressions to examine the evolution of scalar and tensor perturbations, and consequently in the determination of various observable quantities such as the scalar spectral index and the tensor-to-scalar ratio.
3.1.1. Scalar Perturbation
The scalar perturbation over the FRW metric in the longitudinal gauge is as follows:
where
is the scalar perturbation variable, and
is the conformal time coordinate. Here, it is worth noting that the background evolution has no anisotropic stress, so we work with one scalar perturbation variable, namely
.
The scalar field is perturbed as follows:
where
is the background evolution of the scalar field, and given by
in terms of cosmic time, with
being the instant of bounce. The scalar perturbation equations up to the leading order in the longitudinal gauge are given by [
87]:
where
is the Hubble parameter in the
coordinate, and
denotes the variation of the energy–momentum tensor due to the perturbation of the spacetime metric and the scalar field mentioned above. The variation of the matter energy–momentum tensor in the present context comes with the following expressions:
where we use Equation (
38), and
is obtained earlier in Equation (
36). Plugging back the above expressions of
into Equation (
40) yields the following set of equations:
The second equality of Equation (
42) helps to extract
in terms of
and
, which when used with the other two equalities, one gets the following equation for
:
Clearly,
depends on the background evolution through
and
, which were determined in the previous section in terms of the cosmic time (t). Therefore, it will be more useful if we transform the above equation into cosmic time, for which we need the following relations:
Accordingly, Equation (
42) is given by:
where
is the Hubble parameter in cosmic time. Here, we may argue how the dynamics of the scalar perturbation depends on various terms present in Equation (
43); in particular, the second term leads to an oscillation character of
, the third term produces a friction-like force, and the fourth term acts similarly to a restoring force. As we mentioned earlier, owing to the presence of the Type-IV singularity, the comoving Hubble radius asymptotically goes to zero on both sides of the bounce. Thus, the perturbation modes generate and cross the horizon near the bounce when all the modes lie within the Hubble radius. Therefore, we intend to solve Equation (
43) near the bounce, and we use the near-bounce expressions of
and
from Equations (
31) and (
36), respectively. Consequently, we get the following equation:
where we retain up to the leading order in
. With the above expression, Equation (
43) can be written as follows:
where
p and
q are given by:
respectively. The Fourier transformation of Equation (
45) is as follows:
where
is the Fourier transformed variable of
. Equation (
47) has the following solution for
:
where
is the Hermite polynomial with the order
s. The order of the Hermite polynomial in the above expression is given by the following equation.
Moreover,
is the integration constant which can be determined from the Bunch–Davies condition given by
, where
is the canonical scalar Mukhanov–Sasaki variable. The Bunch–Davies condition is well-justified by the fact that the perturbation modes near the bounce lie within the Hubble radius. The Bunch–Davies vacuum condition on the Mukhanov–Sasaki variable can be transformed into the corresponding condition on
through the following relation [
87]:
where we use
. Due to
from Equation (
35), the above equation can be equivalently written as follows:
Consequently, the integration constant
gets the following form:
Accordingly, the solution of the scalar perturbation from Equation (
48) results in the following equation:
with
p and
q being shown in Equation (
46). Consequently, the scalar power spectrum for the
k-th mode is determined as follows:
We are interested in determining the observable quantities such as the scalar spectral index (
) and the tensor-to-scalar ratio (
r), and we will examine the possible effects of the Type-IV singularity on such observable indices. The horizon-crossing condition for the
k-th mode is
, where
, i.e., we intend to calculate
and
r over the large-scale modes. The occurrence of the Type-IV singularity in the present context leads to the large-scale modes crossing the horizon near the bounce, and thus, by using the near-bounce expression of
(see Equation (
31)), the horizon-crossing condition can be written as follows:
where
symbolizes the horizon-crossing instant of the
k-th mode. Therefore, the scalar power spectrum at the horizon crossing is given by:
which clearly depends on
k through the term containing
as well as through the factor
s (the order of the Hermite polynomial, see Equation (
49)). In particular, the spectral tilt of the scalar power spectrum is defined by the following equation:
However, before estimating the , let us perform the tensor perturbation that is useful for the observable quantity, namely the tensor-to-scalar ratio.
3.1.2. Tensor Perturbation
The tensor perturbation variable satisfies the following equation:
where
is the tensor perturbation variable, and
, in the context of the Lagrange multiplier in the Gauss–Bonnet gravity, is given by [
81,
82,
83]:
Using Equations (
31) and (
33), we determine
as follows:
The Fourier-transformed tensor perturbation variable is defined as
, where
and
represent two polarization modes. Therefore, the above form of
along with Equation (
57) leads to the tensor-perturbed equation in terms of the Fourier-transformed variable, expressed as follows:
where we retain the terms up to the leading order in
and recall that
(see Equation (
46)). Here, it is worth noting that both tensor polarizations (× and + modes) in the present context obey the same differential Equation (
60), due to which we do not put any polarization index in the tensor perturbation variable. However, we will multiply the final expression of the tensor power spectrum by a factor of “2” due to their equal contribution to the spectrum. Solving Equation (
60), we then get the following equation:
where
is the order of the Hermite polynomial and is given by:
Moreover, the integration constant
can be determined from the Bunch–Davies vacuum state near the bounce when the relevant perturbation modes lie within the sub-Hubble regime. In particular, the Bunch–Davies vacuum state is defined by
. Due to
from Equation (
59), the Bunch–Davies condition results in the following equation:
Accordingly, the final solution of
is as follows:
Consequently, the tensor power spectrum for the
k-th mode is given by:
Here, we consider the contribution from both polarization modes of the tensor perturbation. Using Equation (
54), the tensor power spectrum at the horizon crossing comes in the following form:
Therefore, the tensor power spectrum is not scale invariant due to the term containing
as well as due to
(the order of the Hermite polynomial, see Equation (
62)).
We now calculate the scalar spectral tilt and the tensor-to-scalar ratio; the scalar tilt is defined in Equation (
56), while the tensor-to-scalar ratio is given by the following equation:
where the suffix “h” denotes the horizon-crossing instant. Clearly, the
and
r depend on the parameter
, i.e., the instant when the Type-IV singularity occurs. As we mentioned earlier, depending on whether
or
or
, the Type-IV singularity appears before the bounce or after the bounce or at the bounce, respectively. Therefore, in the following, we will estimate
and
r separately for these three cases.
For : In this case, we consider
(for other positive values of
, the main arguments will not change). The theoretical predictions for
and
r, with respect to the parameter
n, are given in
Table 2, which clearly demonstrates that the scalar power spectrum is highly red-tilted, and the tensor-to-scalar ratio gets a large value in respect to the Planck results.
For : Here, the Type-IV singularity occurs after the bounce happens, and we safely consider
. As a result, the
and
r are predicted, and they are shown in
Table 3.
For : Here, the theoretical estimations for
and
r (with respect to
n) are shown in
Table 4.
Therefore, in all three cases, the scalar power spectrum is found to be highly red-tilted. Moreover, the model predicts a large value for the tensor-to-scalar ratio that lies outside of the Planck data, which indicates that the model is not viable with the observational data. However, as we observe in [
53], the
bounce without any finite-time singularity, where the scalar factor is described by
, indeed leads to the simultaneous compatibility of
and
r with the Planck data. Therefore, we may argue that the occurrence of the Type-IV singularity considerably affects the bouncing dynamics in the present context, which in turn results in the non-viability of the model. At present, the appearance of the Type-IV singularity “globally” affects the dynamics of the universe as compared to the case when the singularity is absent. The term “global” means that although the singularity occurs at a finite time
, it controls the asymptotic evolution of the comoving Hubble radius; in particular, the comoving Hubble radius asymptotically goes to zero due to the presence of the singularity. Consequently, the perturbation modes generate near the bounce, unlike the scenario when the singularity is absent and the perturbation modes generate long before the bounce in the deep contracting phase. Such a generation era for the perturbation modes near the bounce is the main reason why the scalar power spectrum gets red-tilted and the tensor-to-scalar ratio has a large value in the present bounce scenario.
Thus, as a whole, the following arguments can be made for the bounce that appears with a Type-IV singularity: (1) If the singularity appears at , then depending on whether or or , the singularity appears before the bounce or after the bounce or at the instant of the bounce, respectively. (2) In all these three cases, the scalar power spectrum becomes red-tilted, and the tensor-to-scalar ratio is too large to be consistent with the Planck data. As we have just mentioned, such inconsistency of the observable quantities is due to the occurrence of the Type-IV singularity, particularly due to the “global” effects of the singularity on the evolution of the universe.
4. Realization of a Bounce with a Type-IV Singularity That Locally Affects the Spacetime
In the previous section, we demonstrated that in the case when the Type-IV singularity “globally” affects the spacetime, the perturbation modes are generated near the bounce. As a result, the observable quantities do not lie within the Planck constraints. Based on these findings, it becomes important to examine a bouncing scenario where the Type-IV singularity “locally” affects the universe’s evolution at around the time when it occurs. This is the subject of the present section. To induce the local effects of the singularity, we introduce a regulating Gaussian factor within the expression of
in Equation (
22). Such regulating factor actually controls when the singularity becomes effective. In particular, the scale factor we consider is given by:
where
acts as the regulating factor, which peaks at around
, i.e., at the time when the singularity occurs. The above expression is similar to the previous form of the scale factor (see Equation (
22)), except the presence of the regulating factor. Once again, the scale factor is written as a product of
and
, where
triggers a bounce scenario and
ensures the occurrence of a finite-time singularity at
. Despite the presence of
, the whole scale factor, i.e.,
, predicts a bouncing universe near
. Therefore, the presence of
results in a finite-time singularity without jeopardizing the bouncing behavior of the universe. Moreover, the term
sitting in the expression of
clearly indicates that
becomes effective only at around
; otherwise,
away from
, and the universe’s evolution is controlled entirely by
. As a result, we may argue that the finite-time singularity locally affects the spacetime at around the time when it occurs. This emphasizes the importance of the regulating factor in producing a
effect for the finite-time singularity on the bouncing dynamics, which in turn reflects the significance of the scale factor considered in Equation (
68) for our present interest. This will be cleared further by the expression of the Hubble parameter defined by
. Equation (
68) immediately leads to the Hubble parameter as:
The above expression of
refers to a Type-IV singularity for
. The appearance of
in the expression of Equation (
69) clearly indicates that the second term in
, which is actually responsible for the singularity, affects the evolution of the Hubble parameter only at around
, i.e., the Type-IV singularity locally affects the spacetime at around the time when it occurs. Therefore, the Hubble parameter of Equation (
69) predicts a bounce at
. Moreover, depending on whether
or
, the Type-IV singularity occurs before the bounce or after the bounce, respectively. Using Equation (
69), we provide the plot for
vs.
t in
Figure 3, where the left and right plots correspond to
and
, respectively.
Figure 3 demonstrates the following points about the Hubble parameter: (1)
becomes zero and a increasing function with the cosmic time at
, which indicates a bounce of the universe at
. (2)
is not symmetric with respect to the bounce point; this is due to the occurrence of the Type-IV singularity in the present cosmological scenario. In fact, the term containing
in the expression of
, which is responsible for the singularity, yields the asymmetric nature of the Hubble parameter. (3) The symmetric nature of
seems to deviate only around
. This, however, is expected because, as we mentioned earlier, the singularity affects the Hubble parameter locally at around
. From Equation (
68), the factor
away from
, and the universe’s evolution is controlled by the symmetric scale factor
.
The comoving Hubble radius is defined by
, where
and
are shown above. Therefore, in the distant past, the comoving Hubble radius turns out to be:
Depending on whether or , the asymptotic nature of becomes different, which in turn reveals the generation era of primordial perturbation modes. In particular, for , the comoving Hubble radius in the present context diverges to infinity in the distant past, and as a consequence, the primordial perturbation modes are generated far away from the bounce in the deep contracting phase. This ensures the Bunch–Davies vacuum state of the perturbation at , and as a result, the horizon problem is resolved. However, for , the comoving Hubble radius asymptotically goes to zero on both sides of the bounce, and hence the perturbation modes generate near the bounce when the Hubble radius is infinite in size, enough to contain all the modes within it. In such a case, the perturbation modes in the distant past lie outside of the Hubble radius, and thus the horizon problem persists for . Based on the above arguments, we will consider so that the perturbation modes are within the sub-Hubble regime in the distant past and the horizon problem can be resolved.
4.1. Cosmological Perturbation and Phenomenology of the Bounce
As mentioned in the previous section, we take
, which leads to the generation era of the perturbation modes long before the bounce in the sub-Hubble regime. Therefore, the useful quantities are the background scale factor, the Hubble parameter as well as its derivative (with respect to the cosmic time), and the Ricci scalar, during the contracting era. They are given by the following:
Using the expression of
, we can write the scale factor as well as the Hubble parameter and its derivative in terms of the Ricci scalar as follows:
where in the last equality, we write
from Equation (
19), and
is a constant that has a mass dimension of [1 + 2n]. Moreover,
and
in the above expressions. Using
and
from Equation (
71), we determine the functions
in the context of the ghost-free Gauss–Bonnet theory of gravity [
81,
82,
83] as:
respectively, and we will use these expressions frequently in the perturbation calculations. Recall that
has a mass dimension of [1 + 2n], and thus from dimensional analysis, we can consider the following:
where
is a positive constant. The parameters
and
are considered to be related by the equation above, and
can be regarded as being the replacement of both of them. Consequently, Equation (
21) along with Equation (
71) immediately lead to the Lagrange multiplier function:
We will use the these expressions in addressing the evolution of scalar and tensor perturbations in the following two subsections.
4.1.1. Scalar Perturbation
As the perturbation modes generate in the distant past, it will be useful to work in the comoving gauge, in which case the second order of perturbed action for the curvature perturbation (symbolized by
) is given by [
81,
82,
83]:
Here, in the context of the Lagrange multiplier
gravity,
and
have the following forms [
81]:
and
respectively, with the function
being defined earlier and
. From Equation (
75), it is clear that the kinetic term of the curvature perturbation comes with a positive sign under the condition
. Later, we will show that this condition, i.e.,
, indeed holds in the present context, which in turn ensures the stability of the curvature perturbation. By using Equations (
71) and (
72), we determine the various terms present in the expression of
as follows:
and
respectively. Consequently, the form of
from Equation (
76) becomes:
where
and
have the following forms:
and
respectively. Note the dependence of
on the parameter
, which actually arises due to the occurrence of the Type-IV singularity during the universe’s evolution. Our intention is to examine how the observable quantities (such as the scalar spectral index and the tensor-to-scalar ratio) depend on
when the Type-IV singularity affects the spacetime locally at around
. As demonstrated earlier, the perturbation modes are generated during the late contracting phase, when the Ricci scalar satisfies the condition such as
as
at
(the numerical estimation of
is given after Equation (
91)). As a consequence, Equation (
78) leads to
, which makes the curvature perturbation stable.
It will be more useful if we transform the time coordinate to the conformal time defined by
. Due to Equation (
68), the scale factor in the distant past behaves similarly to
, and hence the corresponding conformal time is expressed as follows:
Recall that
in order to resolve the horizon issue, due to which
appears to be a monotonic increasing function of the cosmic time. With Equation (
81), we obtain the Ricci scalar in terms of
as:
Using the expression of
above in Equation (
78), we can obtain
as follows:
with
and
. Consequently, we determine the factor
(which is essential for solving the scalar Mukhanov–Sasaki equation):
where
, and we use
. Furthermore, with the aforementioned expressions of
(see Equation (
72)), the speed of the scalar perturbation from Equation (
77) turns out to be the following:
At this stage, we now introduce the scalar Mukhanov–Sasaki (MS) equation for the variable
(also known as the MS variable):
where
is the Fourier mode for
. Clearly, the evolution of
depends on the background evolution through the factors
and
(the over prime denotes the derivative with respect to
). Due to the condition
(as depicted after Equation (
80)),
and
can be expressed as:
respectively, where we retain the terms up to the order
. Furthermore,
along with
(in order to generate the perturbation modes long before the bounce and consequently, to resolve the horizon issue) clearly depict that the term
within the parenthesis can be safely considered as small during the contracting era. As a result,
becomes proportional to
, i.e.,
, with:
which is approximately a constant during the era when the perturbation modes generate deep inside the Hubble radius. Accordingly, along with
, we solve
from Equation (
86), and it is given by the following equation:
with
, and moreover,
and
are the Hermite functions (having an order of
) of the first and the second kind, respectively. Here,
,
are integration constants which can be determined from the initial condition of the MS variable. The Bunch–Davies vacuum state is considered to be the initial state for
, in particular,
. The Bunch–Davies initial condition is ensured by the fact that the perturbation modes generate in the deep sub-Hubble regime (when all the perturbation modes lie within the Hubble radius). The Bunch–Davies condition immediately leads to
and
, respectively. Consequently, the scalar power spectrum for the
k-th mode is as follows:
where in the second equality, we use the solution of
. The
k-th mode satisfies the relation
at the instant of horizon crossing, which, by using Equation (
71), is obtained as follows:
where the suffix “
h” represents the horizon-crossing instant. Equation (
91) estimates the horizon-crossing time for large-scale modes, particularly for
(around which we will determine the observable quantities), as
. This is, however, expected because the large-scale modes re-enter the horizon at around the present epoch, when the age of the universe is nearly
. Since the universe is almost symmetrical with respect to the bounce point (except at around
, when the Type-IV singularity occurs), one can already guess that the large-scale modes cross the horizon during the contracting phase, nearly at
, which is also reflected by Equation (
91). Consequently, we estimate the Ricci scalar at the horizon crossing of the large-scale modes; in particular, we get
(where we consider
,
, and
. We will show that such considerations of
n,
, and
are indeed consistent with the viability of the observable quantities in respect to the Planck data). This justifies the condition
, which we considered earlier in the expression of
.
Equation (
91) depicts the sub-Hubble and super-Hubble regime of the
k-th mode as:
As a result, the scalar power spectrum (from Equation (
90)) in the super-Hubble regime can be expressed as follows:
with recall that
. By using Equation (
93), we can determine the spectral tilt of the primordial curvature perturbations (symbolized by
). Clearly,
depends on
, and thus the occurrence of the Type-IV singularity affects the scalar power spectrum as well as the corresponding spectral tilt. However, before proceeding to calculate
, we will determine the tensor power spectrum, which is necessary for the prediction of the tensor-to-scalar ratio.
4.1.2. Tensor Perturbation
The tensor perturbation over the FRW metric is expressed as follows:
where
is the tensor perturbation variable, and the corresponding tensor perturbed action (up to the quadratic order) is given by [
81,
82,
83,
84,
85]:
In the Lagrange multiplier Gauss–Bonnet gravity theory, the function
is [
81]:
where
, and
is given in Equation (
72). Equation (
95) indicates that the speed of the tensor perturbation (or equivalently, the gravitational waves) is equal to unity—this is due to the fact that the Gauss–Bonnet coupling in the present context satisfies
, which in turn makes
. Consequently, the model becomes compatible with the GW170817 event. The scale factor from Equation (
71) immediately leads to the following expression of
:
Therefore,
is positive, which results in the stability of the tensor perturbation. By using Equation (
82), i.e.,
, we determine
in terms of the conformal time as:
where
is given by the following equation:
Accordingly, we calculate
:
where we recall that
, and we use
in order to arrive at the above expression. The above expression will be useful for solving the tensor Mukhanov–Sasaki equation. Due to the condition
(as demonstrated earlier) along with
, the term containing
can be safely regarded as small during the contracting phase. As a result,
becomes proportional to
, i.e.,
, with:
In effect, the tensor Mukhanov–Sasaki (MS) equation becomes:
where
is the Fourier-transformed quantity of the tensor MS variable, which is defined by
. Considering the Bunch–Davies initial condition for
, i.e.,
, we solve Equation (
102) as,
where
, and
represents the Hermite function of the second kind and having order of
. Consequently, the tensor power spectrum for the
k-th mode in the superhorizon scale (when the relevant modes are outside of the Hubble radius and satisfy
from Equation (
92)) comes with the following expression:
where we consider the contributions from both polarization modes.
Having set the stage, we now calculate the observable quantities such as the scalar spectral index (
) and the tensor-to-scalar ratio (
r), respectively. They are defined by the following equation:
where the suffix “
h” represents the horizon-crossing instant of the large-scale modes (
), around which we will estimate the observable indices. According to the recent Planck data,
and
r are constrained by [
47]:
respectively. Due to Equations (
93) and (
104), we determine the final forms of
and
r in the present context as follows:
where all the quantities are evaluated at the horizon crossing of large-scale modes, in particular:
Here, we would like to mention that the dependence of
and
r on the parameter
actually decodes the possible effects of the Type-IV singularity on the observable quantities. Clearly, the above expressions contain
, which is the Ricci scalar at the horizon crossing of the large-scale modes. Hence, from Equation (
82), one may write the following equation:
where
is shown in Equation (
91), in particular,
Here, we use
, which crosses the horizon during the contracting phase, nearly around
. Plugging back the above expression of
into Equation (
109), we get
in terms of
n and
:
Thus, as a whole, the theoretical expressions of
and
r depend on the parameters
n,
, and
. We would like to mention that the scalar tilt as well as the tensor-to-scalar ratio do not depend on the parameter
(the time when the Type-IV singularity occurs). This, however, is expected because the singularity locally affects the spacetime at around the finite time
, and the perturbation modes are generated in the deep contracting phase where the singularity provides almost no effects on the universe’s evolution. This is unlike the previous scenario where the Type-IV singularity globally affects the spacetime, and as a result, the observable quantities are found to depend on
(see the discussion after
Table 4).
It turns out that the theoretical predictions of the scalar spectral index and the tensor-to-scalar ratio in the present case simultaneously become compatible with the Planck 2018 data for a small range of the parameters, given by:
,
and
. Therefore, the viable range of
n seems to be less than that of the matter bounce scenario, where
; this result is in agreement with [
53]. The parametric plot
vs.
r is depicted in
Figure 4.
Therefore, in the context of the ghost-free Gauss–Bonnet theory of gravity, the bouncing scenario in presence of a Type-IV singularity—where the Type-IV singularity affects the spacetime at around the time when the singularity occurs—proves to be viable in respect to the Planck observations.
5. Conclusions
We examined how the presence of a Type-IV singularity can influence the dynamics of a bouncing universe; that is, we worked with bounce cosmology that appears with a Type-IV singularity at a finite time. In the case of a Type-IV singularity, the scale factor, the Hubble parameter, and its first derivative are finite; however, the higher derivatives of the Hubble parameter diverge at the time of the singularity. Therefore, the Type-IV singularity is not a crushing type, unlike the initial or the Big Rip singularity, and the universe can smoothly pass through a Type-IV singularity. However, the presence of such a Type-IV singularity can severely influences the dynamics of the universe, as we found here for an otherwise non-singular bounce scenario. The theory of gravity considered was the well-known ghost-free Gauss–Bonnet (GB) gravity, where the ghost-free nature is ensured by the presence of a Lagrange multiplier in the gravitational action, as developed in [
71]. Moreover, we chose a class of Gauss–Bonnet coupling function (
) that satisfies a constraint equation such as
(where
H is the Hubble parameter of the universe), which in turn leads to the speed of the gravitational wave as unity. Consequently, the model became compatible with the event GW170817. By using the reconstruction technique, we obtained the explicit form of the scalar field potential as well as the GB coupling function, which trigger a bouncing scenario with a Type-IV singularity at a finite time
.
We found three different cases depending on whether or or : (1) For , the bounce seems to happen at some negative time (, where the suffix stands for bounce), and the singularity occurs before the bounce. (2) For , the bounce shows at a positive time (); moreover, the singularity occurs after the bounce. (3) Finally, for , the bounce and the Type-IV singularity occur at the same instant of time, i.e., . This was the first indication that the Type-IV singularity indeed affects the dynamics of the bouncing scenario. Consequently, we analyzed the evolutions of the scalar and tensor perturbations in this context, and determined various observable quantities such as the scalar spectral index () and the tensor-to-scalar ratio (r)—both of which were found to depend on . Therefore, the perturbation power spectra become different depending on the cases, whether or or . However, in all of these cases, the scalar power spectrum proves to be highly red-tilted, and the tensor-to-scalar ratio becomes too large to be consistent with the Planck data. Through a rigorous analysis of the scalar and the tensor perturbation, we were able to show that it is difficult to obtain a matching and r up to the observed CMB spectra, in the case when the Type-IV singularity globally affects the spacetime. Such inconsistency of the observable quantities arises due to the appearance of the Type-IV singularity, which “globally” affects the evolution of the Hubble parameter compared to the bouncing scenario where such a Type-IV singularity is absent. By the term “global”, we mean that although the singularity occurs at a finite time , it significantly affects the asymptotic evolution of the universe in the distant past as well as in the distant future. In particular, the presence of the Type-IV singularity results in the comoving Hubble radius going to zero asymptotically on both sides of the bounce. As a result, the perturbation modes generate near the bounce when all the relevant modes lie within the sub-Hubble regime. This is the reason why the scalar power spectrum shows a red-tilted behavior, and the tensor-to-scalar ratio becomes too large in respect to the observational constraints.
Based on the findings above, we investigated a different bouncing scenario, which appears with a Type-IV singularity. However, the Type-IV singularity “locally” affects the spacetime at around the time when the singularity occurs. As a result, and unlike the previous scenario, the comoving Hubble radius in this bounce scenario diverges to infinity in the distant past. Thus, the primordial perturbation modes generate far away from the bounce in the deep contracting phase. We also calculated the scalar spectral index and the tensor-to-scalar in this context, which are found to be simultaneously compatible with the recent Planck data for the suitable regime of parameter values. This ensures the viability of the bounce model where the Type-IV singularity shows local effects on the spacetime at around the time of the singularity. It is important to mention here that the scalar tilt as well as the tensor-to-scalar ratio were found to be independent of the parameter . This, however, is expected because the singularity affects the spacetime locally at around the time , and the perturbation modes generate in the deep contracting phase where the singularity provides almost no effects on the universe’s evolution.
Thus, as a whole, this work clearly reveals that the presence of a Type-IV singularity has significant effects on an otherwise non-singular bounce scenario. We have shown that a bounce model that appears with a Type-IV singularity is viable if the singularity “locally” affects the spacetime at around the time when it occurs; otherwise, the observable quantities are found to be problematic (with respect to the Planck constraints) when the Type-IV singularity “globally” affects the spacetime. Therefore, in the realm of bouncing cosmology, if the universe had faced a Type-IV singularity during its evolution in the past, then the singularity should have “locally” affected the spacetime.