Cosmological Singularities in Conformal Weyl Gravity

In this work, we study the issue of the past and future cosmological singularities in the context of the fourth-order conformal Weyl gravity. In particular, we investigate the emergent universe scenario proposed by Ellis {\it et al}, and find the stability conditions of the corresponding Einstein static state using the fixed point approach. We show that depending on the values of the parameters of the conformal Weyl gravity theory, there are possibilities for having initially stable emergent states for an FRW universe with both the positive and negative spatial curvatures. This represents that the conformal Weyl gravity can be free of the initial singularity problem. Then, following Barrow {\it et al}, we address the possible types of the finite-time future cosmological singularities. We discuss how these singularities also can be avoided in the context of this theory.


Introduction
The standard model of cosmology based on the Einstein general theory of relativity (GR) has a good agreement with the recent high resolution observations and provides an accurate overview of the cosmic history of our universe. However, it suffers from some fundamental problems such the initial big bang singularity problem, the flatness problem, the horizon problem, the magnetic monopoles and etc. Although the inflation scenario provides solutions for some of these problems, the initial big bang singularity problem still remains unsolved. Indeed, if the inflation scenario is realized by the dynamics of scalar fields coupled to Einstein gravity, then the Hawking-Penrose singularity theorem [1,2] proves that an inflationary universe is geodesically incomplete in its past. Therefore, there exists a singularity before the onset of inflation and consequently the inflationary scenario is not capable to yield the complete history of the very early universe. To deal with this shortcoming, a number of attempts have been made in construction of cosmological models which are initially non-singular. There are two main proposals for this aim named as the "bouncing universe" and "emergent universe" scenarios. The bouncing universe scenario is based on a smooth transition of the universe at a finite radius from a contracting phase to an expanding phase. Then, by going backward in time, the universe collapses up to a finite value for the scale factor, before which it starts to expand again [3]. Although the bouncing scenario avoids the big bang singularity, the existence of the non-singular bounce requires non-standard matter fields violating energy conditions [4]. In the emergent universe scenario, the universe emerges from a static state, namely as the "Einstein static universe (ESU)", characterized by a non-zero spatial volume with a positive curvature [5]. This scenario avoids the initial big bang singularity while preserving the standard energy conditions. However, due to the existence of varieties of perturbations, such as the quantum fluctuations and the variations of the average energy density of the system [6], the emergent universe model is unstable and suffers from a fine-tuning problem [7]. This unstability problem can be amended by the modifications of the cosmological field equations of Einstein's GR. Indeed, there are two conditions which any gravitational theory should satisfy for admitting an emergent universe paradigm: i) the existence and stability of the ESU, and ii) the possibility of joining the standard cosmological history by a graceful exit mechanism from the ESU. By this motivation, the emergent universe scenario and the stability issue of its Einstein static state has been explored in the context of various modified theories of gravity such as the loop quantum gravity [8], f (R) gravity [9], f (T ) gravity [10], Einstein-Cartan theory [11], Braneworlds [12] and Massive gravity [13] as well as in the presence of vacuum energy corresponding to conformally invariant fields [14] among the others. Most of the mentioned references just focused on the first condition without addressing the later one.
On the other hand, after the discovery of the accelerating expansion of the universe, deeper studies of the features of the cosmic fluid responsible for this accelerating expansion, the so called "dark energy", showed the plethora of new types of singularities different than the initial big bang singularity. These new types of singularities are future singularities for a universe and they are characterized by the violation of some of the energy conditions. This violation in turn results in the divergencies in some of the physical quantities such as the scale factor, energy density and pressure profiles. As instances for these type of singularities, one may refer to Big Rip singularity [15], Sudden singularity [16] and Big Freeze singularity [17]. There are also other classifications for the possible future singularities, see for example [18].
The theory of conformal Weyl gravity has been proposed as an alternative theory to Einstein's GR for the high energy limit and is based on the assumption of the local conformal invariance in the geometry of spacetime [19]. There are many attempts for obtaining the Einstein's GR as the low energy limit of the conformal Weyl theory, by dynamically breaking the conformal invariance, with varying degrees of success and difficulty. In this direction, the main problem is that by breaking the conformal symmetry, the quantum fluctuations bring back the Einstein-Hilbert action with a repulsive rather than attractive gravity. However, one may refer to the recent work of Maldacena [20] where it is shown that in the four dimension, the conformal gravity with a Neumann boundary condition is classically equivalent to the Einstein gravity with a cosmological constant. It has been also suggested by Kazanas and Mannheim that the conformal Weyl gravity can be considered an independent theory on its own right instead of looking for a low energy GR limit. By this motivation, the exact vacuum solutions as well as the cosmological solutions of this theory have been found [21,22,23,24]. It is shown that in the cosmological setup, the conformal Weyl gravity can solve naturally both the cosmological constant and the flatness problems in GR. Beside these successes of the conformal Weyl gravity, it is argued that this theory fails to fulfill simultaneously the observational constraints on the present cosmological parameters and on the primordial light element abundances [25]. For the static setup of the theory, it is shown that the de Sitter space can be found as the vacuum solution of this theory and there is no need to an ad hoc introduction of the cosmological constant to the gravitational action and its corresponding field equations [21]. Moreover, in the Newtonian limit of the theory, there are also linear modification terms to the exterior Schwarzschild-de Sitter solutions that could explain the rotation curves of galaxies without any need to the mysterious dark matter [21,26,27]. Another argument about this theory is given in [28]. It was argued that the conformal Weyl gravity does not agree with the predictions of general relativity in the limit of weak fields and non-relativistic velocities. Nevertheless, it was counter-argued in [29], that in the presence of macroscopic long range scalar fields, the standard Schwarzschild phenomenology of GR is still recovered. In order to check the viability of the conformal Weyl gravity, the deflection of light and time delay in the exterior of a static spherically symmetric source were investigated in [30], and it was found that the parameters of the theory fit the experimental constraints. Also, it is shown in [31] that one can find a class of wormhole geometries satisfying the energy conditions in the neighborhood of the throat of wormhole. This is in a clear contrast to the solutions in GR where exotic matter fields are needed to support a wormhole geometry.
In this work, regarding the above mentioned interesting features of the conformal Weyl gravity, we are motivated to study the issue of the past and future cosmological singularities in the context of this theory. Such a study can be also well motivated even one does not consider the Weyl gravity theory as an independent theory on its own right, but as the high energy limit alternative to GR where the being of singularities are important and inevitable issues. The organization of this paper is as follows. In Section 2, we review the conformal Weyl gravity theory and its cosmological field equations. In Section 3, we explore the existence and stability conditions for an ESU within the emergent universe scenario. In Section 4, we address the possible future finite-time cosmological singularities. Finally, in Section 5, we give our concluding remarks.

Conformal Weyl Gravity Theory
The action of the conformal Weyl gravity is given by where C µνρλ is the Weyl tensor and L M is the Lagrangian of the matter fields. Here, α is a dimensionless gravitational coupling constant. The conformal Weyl tensor and matter Lagrangian read respectively as where S(x) is a conformal scalar field that yields symmetry breaking and makes the particles to be massive, and ψ(x) is a fermion field representing the normal matter fields. Also, Γ(x), R µ µ , h and λ are the fermion spin connection, scalar curvature and two dimensionless coupling constants, respectively. The gravitational field equations of this theory can be obtained by the variation of the action (1) with respect to the metric g µν as [21] 4αW µν = T µν , where T µν is the energy-momentum tensor corresponding to the Lagrangian of the matter fields and W µν is given by where W (1) µν and W (2) µν are Due to the conformal invariance of the theory both of the T µν and W µν are kinematically traceless. Considering the FRW metric which is a conformally flat metric, one finds the Friedmann equations [23] where ρ(t) and p(t) are the density and pressure profiles of the perfect fluid T µν and S(x) = S is a spacetime independent constant due to the conformal invariance of the theory. The covariant conservation of T µν leads to [23] where H =ȧ(t)/a(t) is the Hubble parameter.
In the next section, using the cosmological field equations (8), (9) and (10), we study the problem of initial big bang singularity which exists in the standard model of cosmology based on Einstein's GR. We show that the conformal Weyl gravity admits the emergent universe scenario and then it can be free of the initial singularity problem.

Emergent Universe Paradigm: The Initially Non-singular State and its Stability Analysis
As mentioned in the introduction, the emergent universe paradigm in GR is unstable against prevailing perturbations in very early universe and is expected to decay rapidly [7]. Therefore this scenario was developed in modified gravity theories with the hope to improve the stability at the high energy regimes. Thus, the conformal Weyl gravity as a successful theory in solving some of the fundamental problems in standard model of cosmology also deserves for studying the singularity problems in detail, which is the aim of this and the next sections. To study the emergent universe paradigm and its stability in the context of the conformal Weyl gravity, here we consider that the perfect fluid has the barotropic equation of state p = ωρ. Then, we can combine the Friedmann equations (8) and (9) as the Raychadhuri equation The ESU is characterized by the conditionä(t) = 0 =ȧ(t) [7]. Then, to begin with, one has to obtain the existence condition for an ESU solution in the conformal Weyl theory. The corresponding matter density and scale factor for ESU can be obtained from (8) and (11) as Using (13), the existence condition for an ESU solution reduces to the reality of the corresponding scale factor a ES . Then, one can obtain the existence conditions as they are classified in Table 1. Then, the following points can be realized regarding (12), (13) and Table 1.
• One notes that since the scale factor of the initial ESU a ES has very small size, the corresponding density in (12) will be extremely high unless the scalar field S is very small too. In turn, this represents that the matter density ρ ES can be arbitrary small too as the consequence of the smallness of S. This is one of the unique and interesting properties of the fourth order conformal Weyl gravity. The interpretation is that the Universe can born from nothing and any cosmological inhomogeneities such as those due to the large scale structures represent a perturbation around the vacuum state T µν = 0, which has mentioned also by Mannheim [23]. Indeed, this is a direct result of the vanishing rank four Weyl tensor C µναβ and consequently the rank two gravitational tensor W µν in the field equation of conformal Weyl gravity (4) for an FRW metric. There is no such an interesting property for an FRW universe in the context of the second order Einstein's GR theory.
• In the context of Weyl gravity there exists the possibility of having ESU for the ordinary matter fields possessing the equation of state ω ≥ 0, and there is no need for exotic matter fields.
• The case of ω = −1, representing the cosmological constant field, does not correspond to an small size initial emergent ESU for the conformal Weyl gravity theory. Indeed, as the equation of state parameter ω approaches to −1, we have a ES → ∞. This fact has two meanings: i) there is an extremely large static de Sitter limit for the theory, and ii) by approaching the equation of state parameter ω to −1, possibly due to the perturbation in T µν in very early universe [23], the initial small size non-singular Einstein static state enlarges and enters to an inflationary phase which makes hint to a natural graceful exit mechanism. At the end of this section, we will discuss this case again.
• There is no non-singular spatially flat ESU for the conformal Weyl theory, i.e., for k = 0.
• The case of ω > −1/3, λ > 0 and k < 0 has a particular importance in the sense that it can provide a solution for the flatness problem in the standard model of cosmology, see [23].
• The case ω = −1/3 represents a singular state, i.e., a ES = 0. Then, any perturbation around T µν = 0 toward the equation of state ω = −1/3 for matter fields yields an initially singular universe for the conformal Weyl gravity. Now, to study the stability of ESU, we consider the Raychadhuri equation (11) as the following two dimensional dynamical system with the phase-space variables x 1 = a(t) and x 2 =ȧ(t) Then, for this dynamical system, ESU is corresponding to the fixed point (x 1 , x 2 ) = (a ES , 0). The stability of this equilibrium point can be realized by determining eigenvalues σ 2 of the Jacobian matrix J ij = ∂ẋi ∂xj | x2=0 . The case σ 2 > 0 represents a hyperbolic fixed point, and hence it will be unstable. For this case, the trajectories in the neighborhood of this point diverge from it. In contrast, σ 2 < 0 represents a stable center equilibrium point in which a small departure form this point results in oscillations about it. Indeed this case corresponds to the initial stable Einstein static state in which the seed universe oscillates around it indefinitely. The Jacobian matrix associated to our dynamical system (14) reads as where its eigenvalues can be obtained as where we have used (13) in (16). In Figure 1, we give the two dimensional contour plot representing σ 2 (λ, ω). The presence of σ 2 (λ, ω) < 0 regions imply the stability of the solutions. Then, one finds that there are possible stable solutions satisfying the following stability conditions Stability conditions: Then, regarding the conditions in Table 1 and in (17), one can conclude the required conditions guaranteing both the existence and stability of an ESU in the conformal Weyl gravity as are given in Table 2.
In Figure 2, we have plotted the phase-space (ȧ(t), a(t)) diagram and the evolution of the scale factor for the corresponding closed trajectories in the phase-space. The upper plot represents some closed trajectories in the phase-space for the given typical values of the parameters satisfying the existence and stability conditions in Table 2. In general, a closed trajectory in the phase space, as a limit cycle, represents an oscillatory system  1 . The lower plot represents the evolution of the corresponding oscillatory scale factors, versus the time, to each of the closed trajectories in the phase-space. The corresponding cases are denoted by the same colors and the same set of parameters. Until here, we have proved that the conformal Weyl gravity admits an stable ESU for the emergent universe scenario. As we addressed before, the matter density supporting this ESU can be arbitrary small, representing the possibility for being born from nothing [23]. Beside the possibility of being born from nothing for the universe, the ESU scenario in the conformal Weyl gravity has one another important theoretical advantage in comparison to the ESU in Einstein gravity as well as to some of the other modified gravity theories. This advantage is related to the fact that the de Sitter solution is a vacuum solution of the conformal Weyl theory. This vacuum solution is capable of triggering the inflation such that under the prevailing perturbations the universe can transit from the oscillations around its small size ESU to its subsequent inflationary phase and then produces whole of the cosmic history [21]. This advantage is supported by the fact that for both the FRW solution and the de Sitter solution (represented by R µν = Λg µν ), both of the Weyl tensor C µναβ and the gravitational tensor W µν vanish. Then, as mentioned by Mannheim [23], the FRW cosmology is also the case that the energy momentum tensor is everywhere zero, and then it is another vacuum solution of the theory. Thus, any cosmological inhomogeneities such as those in the large-scale structure of the universe is a consequence of a perturbation around T µν = 0. In the presence of normal matter, such a property can be achieved nontrivially by an interplay between the various fields of the theory. In the conformal Weyl gravity, the perturbations around T µν = 0 are responsible for both the generation the matter fields (ρ(t) and p(t)) associated to the inhomogeneities in the large scale structures and the transition from one non-singular vacuum state to its subsequent another vacuum state deriving the inflation. This unique feature of having these two vacuum solutions in the conformal Weyl gravity may provides a natural mechanism for the graceful exit from the ESU to its subsequent inflationary phase and then to the radiation and matter dominated eras. We aim to study and elaborate on this important advantage of the conformal Weyl theory and report it elsewhere [32].
Before concluding this section, we consider an specific solution: a non-singular solution for ω = 1/3 and ρ(t) = ρ 0 /a 4 (t) which represents a radiation filled universe. For this case, the Friedmann equation (8) can 1 The differential equation (14) possesses the Jacobi elliptic function x 1 = A cn(ωt + δ) as its solution which satisfies the equations where k 1 is an arbitrary constant. In the case of k 1 = 0, the solution reduces to the form x 1 = A cos(ωt + δ). be integrated which gives the following solution for the scale factor This solution represents the oscillation of the scale factor of the radiation universe between two finite minimum and maximum sizes. Then, one can interpret this solution in two different ways as i) An ever oscillating radiation dominated universe model with the finite minimum and maximum sizes.
ii) An eternally oscillating small size emergent universe standing as the seed of a large scale universe. In this case, the large scale universe can born and exit from this finite size state and then experience an inflationary phase by the prevailing perturbations around T µν = 0.
Believing in the history of the universe having subsequent phase transitions from the radiation and dust to the dark energy dominated phases seems to be more in the favor of the above second interpretation. This is also supported by the current observations about the accelerating expansion of the universe in the sense that the repulsive dark energy profile makes unlikely the being of a contracting phase. Here, regarding the square roots in (18), one sees that to have the oscillatory modes around the scale factor a 0 = − k 4λS 2 , the conditions k < 0 and λ > 0 are required. These conditions are in agreement with our results for an emergent ESU given in Table 2. Also it is seen that ρ should satisfy the condition ρ 0 < k 2 /16λ which represents that the initial density of ESU state is finite for a finite and non-zero values of k and λ parameters, respectively.

The Possible Future Singularities
In this section, following Barrow et al [33,34], we study the possible types of finite-time singularities may happen in the future of an FRW universe in the context of the conformal Weyl gravity. To do this, based on our study in the past section for the possibility of having an initially non-singular cosmology in the context of this theory, we construct initially non-singular scale factors and then investigate the behavior of H(t),ä(t), ρ(t) and p(t) versus the cosmic time.
In Table 3, we give some types of possible future singularities for an FRW universe regarding the behavior of the scale factor, matter density and pressure in a finite cosmological time [35]. Types I, II, III and IV singularities are known as "Big Rip or Cosmic Doomsday" [36], "Sudden" [33], "Big Freeze" [37] and "Big Brake or Big Separation" [38] singularities, respectively. In type IV singularity, although a(t), ρ(t) and p(t) remain finite but the higher derivatives of H diverge. Following [33,34], and to keep the generality of the study in this section, we do not impose any specific equation of state on the pressure and density profiles in the field equations (8) and (9). We consider the solution for the scale factor a(t) introduced by Barrow [33,34] where B, q, C and n are positive free constants to be determined, and t s represents the occurring moment of the singularity. For an initially singular universe, i.e possessing the big bang singularity, one can fix the zero of time by the condition a(0) = 0 which leads to C = −1/t n s and the following form for the scale factor where q > 0 and n > 0. In the present work, regarding the previous section showing the possibility of having initially non-singular universe in the context of the conformal Weyl gravity, we consider the following generalization for the above scale factor where we supposed 0 < a ES ≪ 1 < a s and it describes a non-singular universe evolving from a(0) = a ES , where a ES represents the initially non-singular Einstein static state, toward a(t s ) = a s . One can simply recover the initially singular universe scenario by setting a ES = 0. Then, using (21), we obtaiṅ On the other hand, using the field equations (8) and (9), we find the density and pressure as Here, one notes to the difference between the field equations of the fourth-order conformal Weyl gravity and the second-order Einstein gravity. The sign differences in ρ(t) and p(t) is a natural result of the difference between the actions of these two theories. Using this difference in (24), it was shown by Mannheim [23] that the conformal Weyl gravity is free of the flatness problem. Another interesting feature of the conformal Weyl gravity is the coupling way of the scalar field S to the dynamical quantities a(t) andȧ(t) which can play an important role in the future singularity avoidance as we will discuss later.
Considering that the coupling λ and the scalar field S are supposed to be constant in the conformal Weyl gravity [23], we can classify the possible finite time future singularities in the context of this theory with respect to n and q values as in the following cases.
• The case of t → t s with 0 < n < 1 and 0 < q ≤ 1.
For this case we have This case represents the type III singularity as in the table (3), or Big Freeze singularity. For this case, the evolution of the scale factor a(t), density ρ(t) and pressure p(t) are plotted in Figure 3 for two typical values of q.
• The case of t → t s with n = 1 and 0 < q ≤ 1. For this case we find where all the a s ,ȧ s , H s ,ä s , ρ s and p s are finite. Then, there is no finite-time future singularity for this case. Here, the zero acceleration case happens for q = 1. Similarly, the evolution of the scale factor a(t), density ρ(t) and pressure p(t) are plotted in Figure 4 for two typical values of q.
• The case of t → t s with 1 < n < 2 and 0 < q ≤ 1.
For this case we have  This case represents a type II singularity, i.e. a cosmological sudden singularity. For this case, the evolution of the scale factor a(t), density ρ(t) and pressure p(t) are plotted in Figure 5 for typical values of q.
• The case of t → t s with n ≥ 2 and 0 < q ≤ 1.
For this case we have where all a s ,ȧ s , H s ,ä s , ρ s and p s are finite. Then, for this case all the physical quantities remain finite. The evolution of the scale factor a(t), density ρ(t) and pressure p(t) are plotted in Figure 6 for two typical values of q. Then, we see that for an initially non-singular universe which evolves according to (21), in general there are type II and III cosmological singularities in the context of conformal Weyl gravity. The type IV singularity, i.e t → t s , a → a s , ρ → 0 and |p| → 0, in the context of this theory with the scale factor (21) can be achieved for the very small values of the scalar field S values (S → 0) such that Sȧ s → 0 and S 2ä s → 0. Even in the cases of divergentȧ(t) andä(t), the limiting values of Sȧ s → 0 and S 2ä s → 0 can be obtained as the result of 0 × ∞ type limits which requireä (t) a 2 (t) and ... a (t) a 2 (t) to be finite, respectively. Finally, for the scale factor (21), for 0 ≤ t ≤ t s we have a ES ≤ a(t) ≤ a s , then the big rip singularity defined as t → t s , a → ∞, ρ → ∞ and |p| → ∞ can not be achieved for a finite time. However, one may consider the scale factor for a universe evolving from an initial non-singular state a ES where 0 ≤ t ≤ t s and n ≥ 1. Then, through the field equations (24) and (25), for t → t s , we find a → ∞, ρ → −∞ and |p| → ∞ representing a Big Rip type singularity. We conclude our analysis for the possible future cosmological singularities in the context of the conformal Weyl gravity in the following two points. i) In general, there are possibilities for having the types II, III and I cosmological singularities for an initially non-singular universe evolving according to the scale factors (21) and (34), respectively. The type IV singularity requires that the scalar field S takes very small values andä (t) a 2 (t) and ... a (t) a 2 (t) to be finite. Then, the final fate of the cosmos in the conformal Weyl gravity can be a sudden, big freeze, big rip or big brake singularity depending on the evolution of the scale factor a(t) and scalar field S value.
ii) The conformal Weyl gravity possesses one interesting property: the unique type of coupling the scalar field S to the dynamical quantities a(t), H(t), andä(t). Indeed, these dynamical quantities are the origins of various type of future cosmological singularities. For the singularities that the divergencies in ρ(t) and p(t) are due to the divergencies in H andä(t), one may consider the following finite-time finetunings: S H(t s ) → f inite and Sä(t s ) → f inite. Then, by these fine-tunings, the scalar field S can control and regularize the corresponding singularities and makes free the conformal theory from any singularity problem. The smallness of S as a result of being the FRW solution as the vacuum solution for the conformal Weyl theory guarantees the avoidance of finite time future singularities. There is no such a possibility in Einstein's GR for a scalar field as the matter source of the field equations. Then, the conformal Weyl gravity can provide a cosmological model which is free of both the past and future singularities.

Conclusion
In the present work, we studied the issue of the past and future cosmological singularities in the context of the fourth-order conformal Weyl gravity. For the past singularity problem, we investigated the emergent universe scenario proposed by Ellis et al. We obtained the stability conditions for the corresponding ESU using the fixed point approach. We showed that depending on the values of the parameters of the conformal Weyl gravity theory, there are possibilities for having initially stable ESUs for an FRW universe with both the positive and negative spatial curvatures. In particular, it is found that there exists the possibility of having ESU for the ordinary matter fields possessing the equation of state ω ≥ 0, but the conformal Weyl theory is not capable to have a non-singular spatially flat ESU. In contrast to GR, here we found that the matter density of the initial ESU, i.e ρ ES , can be arbitrary small which represents the possibility of being born form nothing. Moreover, the fact that both the FRW solution and de Sitter solution are vacuum solutions of the conformal Weyl theory provides a unique feature for this theory in the sense that the initially non-singular emergent universe can gracefully exit from one vacuum to the other vacuum solution. We discussed that by approaching the equation of state parameter ω to −1, possibly due to the perturbation in T µν = 0 in the early universe [23], the initial small size non-singular ESU enlarges and enters to an inflationary phase. In the conformal Weyl theory, the perturbations around T µν = 0 are responsible for both the generation of the matter fields (ρ(t) and p(t)) associated to the inhomogeneities in the large scale structures and for the transition from one non-singular vacuum state to its subsequent another vacuum state deriving the inflation. This unique feature of having these two vacuum solutions in the conformal Weyl gravity may provides a natural mechanism for a graceful exit from the ESU to inflationary phase and then to its subsequent radiation and matter dominated eras. We aim to study and elaborate on this point and report such a possibility elsewhere [32]. Then, following Barrow et al, we addressed the possible types of the finite-time future cosmological singularities such as Big Rip, Sudden, Big Freeze and Big Brake singularities. We discussed that these future cosmological singularities can be the fate of the universe in the context of this theory in a general setup. However, the conformal Weyl gravity possesses one interesting property: the unique type of coupling of the scalar field S to the dynamical quantities a(t), H(t), andä(t). Indeed, these dynamical quantities are the origins of various type of future cosmological singularities. For the singularities that the divergencies in ρ(t) and p(t) are due to the divergencies in H(t) andä(t), by considering the finite-time fine-tunings S H(t s ) → f inite and Sä(t s ) → f inite, the scalar field S can control and regularize the corresponding singularities and make free the conformal Weyl theory from any singularity problem. The smallness of S as a result of being the FRW solution as the vacuum solution for the conformal Weyl theory guarantees the avoidance of finite time future singularities. There is no such a possibility in Einstein's GR for a scalar field as the matter source of the field equations. Then, the conformal Weyl gravity can provide a cosmological model which is free of both the past and future singularities.