We study the inflationary phenomenology of a non-minimally coupled Einstein Gauss-Bonnet gravity theory, in the presence of a scalar potential, under the condition that the gravitational wave speed of the primordial gravitational waves is equal to unity, that is $c_T^2=1$, in natural units. The equations of motion, which are derived directly from the gravitational action, form a system of differential equations with respect to Hubble's parameter and the inflaton field which are very complicated and cannot be solved analytically, even in the minimal coupling case. In this paper, we present a variety of different approximations which could be used, along with the constraint $c_T^2=1$, in order to produce an inflationary phenomenology compatible with recent observations. All the different approaches are able to lead to viable results if the model coupling functions obey simple relations, however, different approaches contain different approximations which must be obeyed during the first horizon crossing, in order for the model to be rendered correct. Models which may lead to a non-viable phenomenology are presented as well in order to understand better the inner framework of this theory. Furthermore, since the velocity of the gravitational waves is set equal to $c_T^2=1$, as stated by the striking event of GW170817 recently, the non-minimal coupling function, the Gauss-Bonnet scalar coupling and the scalar potential are related to each other. Here, we shall assume no particular form of the scalar potential and we choose freely the scalar functions coupled to the Ricci scalar and the Gauss-Bonnet invariant. Certain models are also studied in order to assess the phenomenological validity of the theory, but we need to note that all approximations must hold true in order for a particular model to be valid.

After the recent GW170817 event of the two neutron stars merging, many string corrected cosmological theories confronted the non-viability peril. This was due to the fact that most of these theories produce massive gravitons primordially. Among these theories were the ones containing a non-minimal kinetic coupling correction term in the Lagrangian, which belong to a subclass of Horndeski theories. In this work we demonstrate how these theories may be revived and we show how these theories can produce primordial gravitational waves with speed $c_T^2=1$ in natural units, thus complying with the GW170817 event. As we show, if the gravitational action of an Einstein-Gauss-Bonnet theory also contains a kinetic coupling of the form $\sim \xi(\phi) G^{\mu\nu}\partial_\mu\phi\partial_\nu\phi$, the condition of having primordial massless gravitons, or equivalently primordial gravitational waves with speed $c_T^2=1$ in natural units, results to certain conditions on the scalar field dependent coupling function of the Gauss-Bonnet term, which is also the non-minimal coupling of the kinetic coupling. We extensively study the phenomenological implications of such a theory focusing on the inflationary era, by only assuming slow-roll dynamics for the scalar field. Accordingly, we briefly study the case that the scalar field evolves in a constant-roll way. By using some illustrative examples, we demonstrate that the viability of the theoretical framework at hand may easily be achieved. Also, theories containing terms of the form $\sim \xi(\phi)\Box\phi g^{\mu\nu}\partial_\mu\phi\partial_\nu\phi$ and $\sim \xi(\phi)\left(g^{\mu\nu}\partial_\mu\phi\partial_\nu\phi\right)^2$ also lead to the same gravitational wave speed as the theory we shall study in this paper, so this covers a larger class of Horndeski theories.

In this work we shall present models of $F(R)$ gravity which realize in a unified way the inflationary era along with a post-inflationary early dark energy era, with the late-time dark energy era. We shall use two approach methods in order to realize the unified cosmological eras, firstly we specify a Hubble rate which may describe the three distinct acceleration eras, and then by using well known $F(R)$ gravity reconstruction techniques, we shall find the differential equation which may yield the $F(R)$ gravity that realizes the cosmologies. In our second approach, we shall present in a qualitative way, several $F(R)$ gravities which unify the inflationary era with the early and late-time dark energy eras, and we discuss several qualitative issues related to the terms that realize the post-inflationary early dark energy era. We quantify our analysis by numerically solving the Friedmann equation, using the redshift as the main variable, and expressing all the physical quantities as functions of the statefinder $y_\mathrm{H}(z)$, which depends on the redshift and the Hubble rate $H(z)$. For the model studied numerically, we present the behavior of some statefinder quantities, like the deceleration parameter, and we calculate the dark energy density parameter and the dark energy equation of state parameter at present time. After demonstrating that the dark energy era is viable, we investigate when the early dark energy term does not affect the late-time era, and this restricts the free parameters of the model.

In this letter we shall demonstrate that the viable $F(R)$ gravities can be classified mainly into two classes of inflationary attractors, either the $R^2$ attractors or the $\alpha$-attractors. To show this, we shall derive the most general relation between the tensor-to-scalar ratio $r$ and the spectral index of primordial curvature perturbations $n_s$, namely the $r-n_s$ relation, by assuming that the slow-roll condition constrains the values of the slow-roll indices. As we show, the relation between the tensor-to-scalar ratio and the spectral index of the primordial curvature perturbations has the form $r\simeq \frac{48 (1-n_s)^2}{(4-x)^2}$, where the dimensionless parameter $x$ contains higher derivatives of the $F(R)$ gravity function with respect to the Ricci scalar, and it is a function of the $e$-foldings number $N$ and may also be a function of the free parameters of the various $F(R)$ gravity models. For $F(R)$ gravities which have a spectral index compatible with the observational data and also yield $x\ll 1$, these belong to the $R^2$-type of attractors, with $r\sim 3 (1-n_s)^2$, and these are viable theories. Moreover, in the case that $x$ takes larger values in specific ranges and is constant for a given $F(R)$ gravity, the resulting $r-n_s$ relation has the form $r\sim 3 \alpha (1-n_s)^2$, where $\alpha$ is a constant. Thus we conclude that the viable $F(R)$ gravities may be classified into two limiting types of $r-n_s$ relations, one identical to the $R^2$ model at leading order in $x$, and one similar to the $\alpha$-attractors $r-n_s$ relation, for the $F(R)$ gravity models that yield $x$ constant. Finally, we also discuss the case that $x$ is not constant.

In this letter we shall provide a formal derivation of the inflationary slow-roll indices for $F(R)$ gravity from the power spectrum without the condition $\dot{\epsilon}_1=0$, frequently used in the literature, where $\epsilon_1=-\frac{\dot{H}}{H^2}$ is the first slow-roll index. We shall employ Karamata's theorem for regularly varying functions, and as we show, we shall derive the same expressions for the scalar and tensor spectral index of $F(R)$ gravity, also appearing in the current literature, without the misleading and rather strong condition $\dot{\epsilon}_1=0$. The only conditions that are needed are the slow-roll assumption, and the smallness of the slow-roll indices and of their derivatives during inflation.

In this work we shall study the late-time behavior of $k$-Essence $f(R)$ gravity without scalar potential, in the presence of matter and radiation perfect fluids. We quantify the late-time study by using the statefinder function $Y_H(z)=\frac{\rho_{DE}}{\rho_m^{(0)}}$, which is a function of the redshift and of the Hubble rate. By appropriately rewriting the Friedmann equation in terms of the redshift and of the function $Y_H(z)$, we numerically solve it using appropriate initial conditions, and we critically examine the effects of the $k$-Essence higher order kinetic terms. As we demonstrate, the effect of the higher order scalar field kinetic terms on the late-time dynamics is radical, since the dark energy oscillations are absent, and in addition, the cosmological physical quantities are compatible with the latest Planck data and also the model is almost indistinguishable from the $\Lambda$ Cold Dark Matter model. This is in contrast to the standard $f(R)$ gravity case, where the oscillations are present. Furthermore, by choosing a different set of values of two of the free parameters of the model, and specifically the coefficient of the higher order kinetic term and of the exponent of $R^{\delta}$ appearing in the $f(R)$ action, we demonstrate that it is possible to obtain $\rho_{DE}<0$ for redshifts $z\sim 2-3.8$, which complies phenomenologically with, and seems to explain, the observational data for the same redshifts, and also to obtain a viable cosmological evolution at $z\sim 0$, at least when the dark energy equation of state parameter and the dark energy density parameters are considered.

In this work we introduce a new theoretical framework for Einstein-Gauss-Bonnet theories of gravity, which results to particularly elegant, functionally simple and transparent gravitational equations of motion, slow-roll indices and the corresponding observational indices. The main requirement is that the Einstein-Gauss-Bonnet theory has to be compatible with the GW170817 event, so the gravitational wave speed $c_T^2$ is required to be $c_T^2\simeq 1$ in natural units. This assumption was also made in a previous work of ours, but in this work we express all the related quantities as functions of the scalar field. The constraint $c_T^2\simeq 1$ restricts the functional form of the scalar Gauss-Bonnet coupling function $\xi(\phi)$ and of the scalar potential $V(\phi)$, which must satisfy a differential equation. However, by also assuming that the slow-roll conditions hold true, the resulting equations of motion and the slow-roll indices acquire particularly simple forms, and also the relation that yields the $e$-foldings number is $N=\int_{\phi_i}^{\phi_f}\xi''/\xi'd \phi$, a fact that enables us to perform particularly simple calculations in order to study the inflationary phenomenological implications of several models. As it proves, the models we presented are compatible with the observational data, and also satisfy all the assumptions made during the process of extracting the gravitational equations of motion. More interestingly, we also investigated the phenomenological implications of an additional condition $\xi'/\xi''\ll 1$, which is motivated by the slow-roll conditions that are imposed on the scalar field evolution and on the Hubble rate, in which case the study is easier. Our approach opens a new window in viable Einstein-Gauss-Bonnet theories of gravity.

We provide a compact review on recent developments on axion $F(R)$ gravity. The axion field is a string theory originating theoretical particle that is a perfect candidate for low-mass particle dark matter. In this review we present how a viable inflationary phenomenology and a viable late-time evolution can be described by an axion $F(R)$ gravity theory, in which the $F(R)$ gravity part can drive in a geometric way the inflationary and the late-time era, and the axion field behaves as dark matter, with its energy density $\rho_a$ behaving as a function of the scale factor as $\rho_a\sim a^{-3}$. We also briefly discuss the effect of a non-trivial axion Chern-Simons coupling on the inflationary phenomenology of the $R^2$ model. Finally, we briefly discuss the effects of a non-minimal coupling of the axion field with the curvature on neutron stars, and also the propagation of gravity waves in Chern-Simons axion gravity.

In the context of a ghost free $f(R,\mathcal{G})$ model, an extended matter bounce scenario is studied where the form of the scale factor is given by $a(t) = (a_0t^2 + 1)^n$. The ghost free character of the model is ensured by the presence of a Lagrange multiplier, as developed in \cite{Nojiri:2018ouv}. The conditions under which, in this model, the speed of gravitational waves becomes equal to the speed of light (equal to one, in natural units), thus becoming compatible with the striking event GW170817, is investigated. It is shown that this happens for a class of Gauss-Bonnet (GB) coupling functions ($h(t)$) which satisfies a constraint equation of the form $\ddot{h} = \dot{h}H$, with $H$ the Hubble parameter. This constraint is then imposed on the ghost free $f(R,\mathcal{G})$ gravity theory to be consistent with the GW170817 event, subsequently, the corresponding non-singular bouncing cosmology with the aforementioned scale factor is extensively studied. The forms of the coupling function and Lagrange multiplier in the "low curvature limit" of the theory are reconstructed, yielding a viable approximation for $n < 1/2$. Correspondingly, by solving the cosmological perturbation equation, the main observable quantities, namely the spectral index, tensor to scalar ratio, and the running index are determined and confronted with the latest Planck 2018 data. Consistency with the data is proven for those parametric regimes that which correspond to $n < 1/2$. This makes the low curvature approximation a viable one for calculating the scalar and tensor power spectra.

We present a model of $F(R)$ gravity in the presence of a string theory motivated misalignment axion like particle materialized in terms of a canonical scalar field minimally coupled with gravity, and we study the cosmological phenomenology of the model, emphasizing mainly on the late-time era. The main result of the paper is that inflation and the dark energy era may be realized in a geometric way by an $F(R)$ gravity, while the axion is the dark matter constituent of the Universe. The $F(R)$ gravity model consists of an $R^2$ term, which as we show dominates the evolution during the early time, thus producing a viable inflationary phenomenology, and a power law term $\sim R^{\delta}$ with $\delta\ll 1 $ and positive, which eventually controls the late-time era. The axion field remains frozen during the inflationary era, which is an effect known for misalignment axions, but as the Universe expands, the axion starts to oscillate, and its energy density scales eventually as we show, as $\rho_a\sim a^{-3}$. After appropriately rewriting the gravitational equations in terms of the redshift $z$, we study in detail the late-time phenomenology of the model, and we compare the results with the $\Lambda$CDM model and the latest Planck 2018 data. As we show, the model for small redshifts $0<z<5$ is phenomenologically similar to the $\Lambda$CDM model, however at large redshifts and deeply in the matter domination era, the results are different from those of the $\Lambda$CDM model due to the dark energy oscillations. For the late-time study we investigate the behavior of several well-known statefinder quantities, like the deceleration parameter, the jerk and $Om(z)$, and we demonstrate that the statefinders which contain lower derivatives of the Hubble rate have similar behavior for both the $\Lambda$CDM and the axion $F(R)$ gravity model.

In this work we shall investigate the phenomenological implications of the constant-roll condition on a $k$-Inflation theory of gravity. The latter theories are particularly promising, since these remained robust to the results of GW170817, since these have a gravitational wave speed $c_T=1$ in natural units. We shall mainly focus on the phenomenology of the $k$-Inflation models, with the only assumption being the slow-roll condition imposed on the first and fourth slow-roll parameters, and the constant-roll condition for the evolution of the scalar field. We present in detail the final form of the gravitational equations of motion that control the dynamics of the cosmological system, with the constant-roll condition imposed, and by using a conveniently, from the perspective of analytical manipulations, chosen potential, we express the slow-roll indices and the resulting observational indices of the theory as functions of the $e$-foldings number. The results of our analysis indicate that the constant-roll $k$-Inflation theory can be compatible with the Planck 2018 data, for a wide range of the free parameters. Also we examine in a quantitative way the effects of the constant-roll condition on the parameter $f_{NL}^{equil}$ on which the bispectrum is proportional, in the equilateral momentum approximation, and we demonstrate that the effect of the constant-roll condition is non-trivial. In effect, non-Gaussianities in the theory may be enhanced, a phenomenon which is known to be produced by constant-roll scalar theories of gravity in general.

In this work we shall study ghost-free non-local $F(R)$ gravity models. Firstly we shall demonstrate how the ghost degrees of freedom may occur in the non-local $F(R)$ gravity models, and accordingly we shall modify appropriately the gravitational action of non-local $F(R)$ gravity models in order to eliminate the ghosts. Also we shall investigate how the (anti-)de Sitter and the Minkowski spacetime cosmological solutions may arise in the theory, and we investigate when these solutions are stable. Moreover, we shall examine the inflationary phenomenology of the Jordan frame ghost-free non-local $F(R)$ gravity. We shall study two $F(R)$ gravity models, the power law $F(R)$ gravity model $\sim R^n$ with $1<n<2$, $n\neq 2$ and the $R^2$ model, assuming that the slow-roll condition holds true for the Hubble rate during the inflationary era $\dot{H}\ll H^2$ and that the general constant-roll condition $\ddot{\phi}=3\beta H\dot{\phi}$ holds true for the evolution of the scalar field, which includes the slow-roll case for $\beta=0$. As we shall demonstrate, the power-law non-local $F(R)$ gravity case can produce a viable inflationary era, compatible with observations, in the constant-roll case, unlike for the $R^2$ model. It is conceivable that the results are model dependent, as in the ordinary vacuum $F(R)$ gravity.

In this work we study non-singular{ bounce cosmology} in the context of the Lagrange multiplier generalized $F(R)$ gravity theory of gravity. We specify our study by using a specific variant form of the well known matter{ bounce cosmology}, with scale factor $a(t) = \left(a_0t^2 + 1 \right)^n$, and we demonstrate that for $n < 1/2$, the primordial curvature perturbations are generated deeply in the contraction era. Particularly, we show explicitly that the perturbation modes exit the horizon at a large negative time during the contraction era, which in turn makes the ``low-curvature'' regime, the era for which the calculations of observational indices related to the primordial power spectrum can be considered reliable. Using the reconstruction techniques for the Lagrange multiplier $F(R)$ gravity, we construct the form of effective $F(R)$ gravity that can realize such a cosmological evolution, and we determine the power spectrum of the primordial curvature perturbations. Accordingly, we calculate the spectral index of the primordial curvature perturbations and the tensor-to-scalar ratio, and we confront these with the latest observational data. We also address the issue of stability of the primordial metric perturbations, and to this end, we determine the form of $F(R)$ which realizes the non-singular cosmology for the whole range of cosmic time $-\infty < t < \infty$, by solving the Friedmann equations without the ``low-curvature'' approximation. This study is performed numerically though, due to the high complexity of the resulting differential equations. We also investigate the energy conditions in the present context. The phenomenology of the non-singular bounce is also studied in the context of a standard $F(R)$ gravity.

In this paper we study a class of quintessential Einstein Gauss-Bonnet models, focusing on their early and late-time phenomenology. With regard to the early-time phenomenology, we formalize the slow-roll evolution of these models and we calculate in detail the spectral index of the primordial curvature perturbations and the tensor-to-scalar ratio. As we demonstrate, the resulting observational indices can be compatible with both the Planck and the BICEP2/Keck-Array observational constraints on inflation. With regard to the late-time behavior, by performing a numerical analysis we demonstrate that the class of models for which the coupling function $\xi (\phi)$ to the Gauss-Bonnet scalar satisfies $\xi (\phi)\sim \frac{1}{V(\phi )}$, produce a similar pattern of evolution, which at late-times is characterized by a decelerating era until some critical redshift, at which point the Universe super-decelerates and subsequently accelerates until present time, with a decreasing rate though. The critical redshift crucially depends on the initial conditions chosen for the scalar field and for all the quintessential Einstein Gauss-Bonnet models studied, the late-time era is realized for large values of the scalar field.

This paper focuses on the realization of an inflationary model from a canonical scalar field theory with a Woods-Saxon potential, in the slow-roll approximation. Our analysis indicates that the observable quantities derived theoretically from our model, namely the spectral index of the primordial scalar curvature perturbations and the tensor-to-scalar ratio, are compatible with the latest Planck collaboration data. We also discuss the qualitative features of the potential, and we show that the value of the scalar field for which the graceful exit occurs, coincides with the inflection point of the scalar potential. We also attempt to study the post-inflation reheating phase of the model, in order to further examine the viability of the Woods-Saxon scalar field model, and as we demonstrate the results indicate viability of the model for this era too, however the instantaneous reheating is not allowed for the model at hand.

In this work we shall study Einstein Gauss-Bonnet theories and we investigate when these can have their gravitational wave speed equal to the speed of light, which is unity in natural units, thus becoming compatible with the striking event GW170817. We demonstrate how this is possible and we show that if the scalar coupling to the Gauss-Bonnet invariant is constrained to satisfy a differential equation, the gravitational wave speed becomes equal to one. Accordingly, we investigate the inflationary phenomenology of the resulting restricted Einstein Gauss-Bonnet model, by assuming that the slow-roll conditions hold true. As we demonstrate, the compatibility with the observational data coming from the Planck 2018 collaboration, can be achieved, even for a power-law potential. We restricted ourselves to the study of the power-law potential, due to the lack of analyticity, however more realistic potentials can be used, in this case though the calculations are not easy to be performed analytically. We also pointed out that a string-corrected extension of the Einstein Gauss-Bonnet model we studied, containing terms of the form $\sim \xi(\phi) G^{ab}\partial_a\phi \partial_b \phi $ can also provide a theory with gravity waves speed $c_T^2=1$ in natural units, if the function $\xi(\phi)$ is appropriately constrained, however in the absence of the Gauss-Bonnet term $\sim \xi(\phi) \mathcal{G}$ the gravity waves speed can never be $c_T^2=1$. Finally, we discuss which extensions of the above models can provide interesting cosmologies, since any combination of $f(R,X,\phi)$ gravities with the above string-corrected Einstein Gauss-Bonnet models can yield $c_T^2=1$, with $X=\frac{1}{2}\partial_{\mu}\phi\partial^{\mu}\phi $.

In this work we investigate several theoretical and phenomenological implications of a scalar -$F(R)$ gravity containing a non-minimal coupling to the scalar curvature. This kind of model is a generalization of axion-$F(R)$ gravity models, so we shall examine several implications of the latter theory. Firstly we study in detail the Einstein frame picture of the model, and also we discuss the dynamics of the cosmological system. By appropriately using the equations of motion, we demonstrate that an arbitrary cosmological evolution can be realized. Also we study the gravitational waves of the theory, and we demonstrate that the speed of their propagation is the same as in $F(R)$ gravity, but there is the possibility of enhancement or dissipation of the gravitational waves, an effect quite similar to the propagation of gravity waves in a viscous fluid. Finally, we examine the energy momentum tensor and we investigate which quantities related to it are conserved. We also present the constraints imposed by the radiation domination era on the non-minimal coupling of the axion scalar field to the scalar curvature.

In this work we investigate the inflationary phenomenological implications of a recently developed ghost-free Gauss-Bonnet theory of gravity. The resulting theory can be viewed as a scalar Einstein-Gauss-Bonnet theory of gravity, so by employing the formalism for cosmological perturbations for the latter theory, we calculate the slow-roll indices and the observational indices, and we compare these with the latest observational data. Due to the presence of a freely chosen function in the model, in principle any cosmological evolution can be realized, so we specify the Hubble rate and the freely chosen function and we examine the phenomenology of the model. Specifically we focus on de Sitter, quasi-de Sitter and a cosmological evolution in which the Hubble rate evolves exponentially, with the last two being more realistic choices for describing inflation. As we demonstrate, the ghost-free model can produce inflationary phenomenology compatible with the observational data. We also briefly address the stability of first order scalar and tensor cosmological perturbations, for the exponential Hubble rate, and as we demonstrate, stability is achieved for the same range of values of the free parameters that guarantee the phenomenological viability of the models.

In this work we introduce an effective model of $f(R)$ gravity containing a non-minimal coupling to the axion scalar field. The axion field is described by the misalignment model, in which the primordial $U(1)$ Peccei-Quinn symmetry is broken during inflation and the $f(R)$ gravity is described by the $R^2$ model, and in addition, the non-minimal coupling has the form $\sim h(\phi)R^{\gamma}$, with $0<\gamma<0.75$. By appropriately constraining the non-minimal coupling at early times, the axion field remains frozen in its primordial vacuum expectation value, and the $R^2$ gravity dominates the inflationary era. As the Universe expands, when $H$ equals the axion mass $m_a$ and for cosmic times for which $m_a\gg H$, the axion field oscillates. By assuming a slowly varying evolution of the axion field, the axion energy density scales as $\rho_a\sim a^{-3}$, where $a$ is the scale factor, regardless of the background Hubble rate, thus behaving as cold dark matter. At late times, the axion still evolves as $\rho_a\sim a^{-3}$, however the Hubble rate of the expansion and thus the dynamical evolution of the Universe is controlled by terms containing the higher derivatives of $\sim R^{\gamma}$, which are related to the non-minimal coupling, and as we demonstrate, the resulting solution of the Friedman equation at late times is an approximate de Sitter evolution. The late-time de Sitter Hubble rate scales as $H\sim \Lambda^{1/2}$, where $\Lambda$ is an integration constant of the theory, which has its allowed values very close to the current value of the cosmological constant. Finally, the theory has a prediction for the existence of a pre-inflationary primordial stiff era, in which the energy density of the axion scales as $\rho_a\sim a^{-6}$.

Among the remaining viable theories that can successfully describe the late-time era is the $k$-Essence theory and in this work we study in detail the phase space of $k$-Essence $f(R)$ gravity in vacuum. This theory can describe in a viable way the inflationary era too, so we shall study the phase space in detail, since this investigation may reveal general properties regarding the inflationary attractors. By appropriately choosing the dimensionless variables corresponding to the cosmological system, we shall construct an autonomous dynamical system, and we find the fixed points of the system. We focus on quasi-de Sitter attractors, but also to radiation and matter domination attractors, and study their stability. As we demonstrate, the phase space is mathematically rich since it contains stable manifold and unstable manifold. With regard to the inflationary attractors, these exist and become asymptotically unstable, a feature which we interpret as a strong hint that the theory has an inherent mechanism for graceful exit from inflation. We describe in full detail the underlying mathematical structures that control the instability of the inflationary attractors, and also we also address the same problem for radiation and matter domination attractors. The whole study is performed for both canonical and phantom scalar fields, and as we demonstrate, the canonical scalar $k$-Essence theory is structurally more appealing in comparison to the phantom theory, a result also demonstrated in the related literature on $k$-Essence $f(R)$ gravity.

In this paper we study in detail the phase space of a cosmological system consisting of two coupled fluids, namely a dark energy fluid coupled with a superfluid dark matter fluid. The dark matter fluid is assumed to have a superfluid equation of state, hence it is not pressureless and our aim is to find the impact of this non-trivial equation of state on the phase space of the coupled system. We shall use two theoretical contexts, namely that of classical cosmology and that of loop quantum cosmology. In the classical case, we investigated the existence and stability of fixed points, and as we will show, no de Sitter fixed points occur, however matter and radiation domination fixed points occur, which are hyperbolic and unstable. We also show that there exist limited sets of initial conditions for which singular solutions occur in the phase space. With regard to the loop quantum cosmology case, we demonstrate that stable de Sitter fixed points exist, for some values of the free parameters of the theory, and interestingly enough, for the same values, singular solutions corresponding to general sets of initial conditions occur. To our knowledge this feature does not occur so frequently in loop quantum cosmological frameworks. We also demonstrate that non-singular solutions corresponding to a general set of initial conditions occur, however these occur when the dark matter superfluid has negative pressure, so it is a rather physically unappealing situation.

In this work we shall investigate the singularity structure of the phase space corresponding to an exponential quintessence dark energy model recently related to swampland models. The dynamical system corresponding to the cosmological system is an autonomous polynomial dynamical system, and by using a mathematical theorem we shall investigate whether finite-time singularities can occur in the dynamical system variables. As we demonstrate, the solutions of the dynamical system are non-singular for all cosmic times and this result is general, meaning that the initial conditions corresponding to the regular solutions, belong to a general set of initial conditions and not to a limited set of initial conditions. As we explain, a dynamical system singularity is not directly related to a physical finite-time singularity. Then, by assuming that the Hubble rate with functional form $H(t)=f_1(t)+f_2(t)(t-t_s)^{\alpha}$, is a solution of the dynamical system, we investigate the implications of the absence of finite-time singularities in the dynamical system variables. As we demonstrate, Big Rip and a Type IV singularities can always occur if $\alpha<-1$ and $\alpha>2$ respectively. However, Type II and Type III singularities cannot occur in the cosmological system, if the Hubble rate we quoted is considered a solution of the cosmological system.

It is quite well known for some time that string inspired axionic terms of the form $\nu (\phi)\tilde{R}R$, known also as Chern-Simons terms, do not affect the scalar perturbations and the background evolution for a flat Friedman-Robertson-Walker Universe. In this paper we study and quantify the implications of the presence of the above term in the context of vacuum $f(R)$. Particularly, we assume that axionic dark matter is present during inflation, and we examine in a quantitative way the effects of axionic Chern-Simons terms on the tensor perturbations. The axion field is quantified in terms of a canonical scalar field, with broken Peccei-Quinn symmetry. The model perfectly describing axions as potential dark matter candidates is based on the so-called misalignment mechanism, in which case the axion is frozen near its non-zero vacuum expectation value during early times in which $H\gg m_a$. In effect, the inflationary era is mainly controlled by the $f(R)$ gravity and the Chern-Simons term. As we demonstrate, the Chern-Simons term may achieve to make a non-viable $f(R)$ gravity theory to be phenomenologically viable, due to the fact that the tensor-to-scalar ratio is significantly reduced, and the same applies to the spectral index of the tensor perturbations $n_T$. Also by studying the Starobinsky model in the presence of the Chern-Simons term, we demonstrate that it is possible to further reduce the amount of primordial gravitational radiation. The issues of having parity violating gravitational waves, also the graceful exit from inflation due to axion oscillations and finally the unification of dark energy-inflation and axion dark matter in the same $f(R)$ gravity-axion dark matter model, are also briefly discussed.

In this work we develop a theoretical framework for Gauss-Bonnet modified gravity theories, in which ghost modes can be eliminated at the equations of motion level. Particularly, after we present how the ghosts can occur at the level of equations of motion, we employ the Lagrange multipliers technique, and by means of constraints we are able to eliminate the ghost modes from Gauss-Bonnet theories of the form $f(\mathcal{G})$ and $F(R,\mathcal{G})$ types. Some cosmological realizations in the context of the ghost free $f(\mathcal{G})$ gravity are presented, by using the reconstruction technique we developed. Finally, we explore the modifications to the Newton law of gravity generated by the ghost-free $f(\mathcal{G})$ theory.

In this work we study a modified version of $f(R)$ gravity in which higher order kinetic terms of a scalar field are added in the action of vacuum $f(R)$ gravity. This type of theory is a type of $k$-essence $f(R)$ gravity, and it belongs to the general class of $f(R,\phi ,X)$ theories of gravity, where $\phi$ is a scalar field and $X=\frac{1}{2}\partial^{\mu}\phi\partial_{\mu}\phi$. We focus on the inflationary phenomenology of the model, in the slow-roll approximation, and we investigate whether viable inflationary evolutions can be realized in the context of this theory. We use two approaches, firstly by imposing the slow-roll conditions and by using a non-viable vacuum $f(R)$ gravity. As we demonstrate, the spectral index of the primordial scalar perturbations and the tensor-to-scalar ratio of the resulting theory can be compatible with the latest observational data. In the second approach, we fix the functional form of the Hubble rate as a function of the $e$-foldings number, and we modify well-known vacuum $f(R)$ gravity reconstruction techniques, in order to find the $k$-essence $f(R)$ gravity which can realize the given Hubble rate. Accordingly, we calculate the slow-roll indices and the corresponding observational indices, and we also provide general formulas of these quantities in the slow-roll approximation. As we demonstrate, viability can be obtained in this case too, however the result is strongly model dependent. In addition, we discuss when ghosts can occur in the theory, and we investigate under which conditions ghosts can be avoided by using a particular class of models. Finally, we qualitatively discuss the existence of inflationary attractors for the non-slow-roll theory, and we provide hints towards finding general de Sitter attractors for the theory at hand.

In this paper we shall study the inflationary aspects of a logarithmic corrected $R^2$ Starobinsky inflation model, in the presence of a Kalb-Ramond field in the gravitational action of $F(R)$ gravity. Our main interest is to pin down the effect of this rank two antisymmetric tensor field on the inflationary phenomenology of the $F(R)$ gravity theory at hand. The effects of the Kalb-Ramond field are expected to be strong during the inflationary era, however as the Universe expands, the energy density of the Kalb-Ramond field scales as $\sim a^{-6}$ so dark matter and radiation dominate over the Kalb-Ramond field effects. In general, antisymmetric fields constitute the field content of superstring theories, and thus their effect at the low-energy limit of the theory is expected to be significant. As we will show, for a flat Friedmann-Robertson-Walker metric, the Kalb-Ramond field actually reduces to a scalar field, so it is feasible to calculate the observational indices of inflation. We shall calculate the spectral index and the tensor-to-scalar ratio for the model at hand, by assuming two conditions for the resulting Kalb-Ramond scalar field, the slow-roll and the constant-roll condition. As we shall demonstrate, in both the slow-roll and constant-roll cases, compatibility with the latest observational data can be achieved. Also the effect of the Kalb-Ramond field on the inflationary phenomenology is to increase the amount of the predicted primordial gravitational radiation, in comparison to the corresponding $f(R)$ gravities, however the results are still compatible with the observational data.

In this paper we study the effects of spatial curvature of the metric on the phase space of vacuum $f(R)$ gravity. Particularly, we appropriately choose the variables of the dynamical system, in order for this to be autonomous, and we study the phase space of the resulting theory, focusing on de Sitter, matter and radiation domination fixed points. Our analysis indicates that the effect of spatial curvature on the phase space is radical, since it destabilizes all the stable de Sitter vacua of the flat spacetime vacuum $f(R)$ gravity phase space, making the phase space having non-trivial unstable submanifolds. This instability occurs regardless if the spacetime has elliptic or hyperbolic spatial sections, and it is also robust towards the choice of initial conditions. We investigate the source of the instability in the system, and also we discuss the stability of the matter and radiation domination vacua, which, as we demonstrate, are also highly unstable. Our results for de Sitter attractors indicate that the stable de Sitter attractors of the vacuum $f(R)$ gravity theory for a flat Universe, are destabilized by the presence of curvature, and this shows that inflation for vacuum $f(R)$ gravity in non-flat spacetime is problematic, at least at the phase space level. This result holds true for both elliptic and hyperbolic spacetimes.

We systematically review some standard issues and also the latest developments of modified gravity in cosmology, emphasizing on inflation, bouncing cosmology and late-time acceleration era. Particularly, we present the formalism of standard modified gravity theory representatives, like $F(R)$, $F(\mathcal{G})$ and $F(T)$ gravity theories, but also several alternative theoretical proposals which appeared in the literature during the last decade. We emphasize on the formalism developed for these theories and we explain how these theories can be considered as viable descriptions for our Universe. Using these theories, we present how a viable inflationary era can be produced in the context of these theories, with the viability being justified if compatibility with the latest observational data is achieved. Also we demonstrate how bouncing cosmologies can actually be described by these theories. Moreover, we systematically discuss several qualitative features of the dark energy era by using the modified gravity formalism, and also we critically discuss how a unified description of inflation with dark energy era can be described by solely using the modified gravity framework. Finally, we also discuss some astrophysical solutions in the context of modified gravity, and several qualitative features of these solutions. The aim of this review is to gather the different modified gravity techniques and form a virtual modified gravity "toolbox", which will contain all the necessary information on inflation, dark energy and bouncing cosmologies in the context of the various forms of modified gravity.

In this work we study a modified version of vacuum $f(R)$ gravity with a kinetic term which consists of the first derivatives of the Ricci scalar. We develop the general formalism of this kinetic Ricci modified $f(R)$ gravity and we emphasize on cosmological applications for a spatially flat cosmological background. By using the formalism of this theory, we investigate how it is possible to realize various cosmological scenarios. Also we demonstrate that this theoretical framework can be treated as a reconstruction method, in the context of which it is possible to realize various exotic cosmologies for ordinary Einstein-Hilbert action. Finally, we derive the scalar-tensor counterpart theory of this kinetic Ricci modified $f(R)$ gravity, and we show the mathematical equivalence of the two theories.

We investigate the cosmological applications of fluids having an equation of state which is the analog to the one related to the isotropic deformation of crystalline solids, that is containing logarithmic terms of the energy density, allowing additionally for a bulk viscosity. We consider two classes of scenarios and we show that they are both capable of triggering the transition from deceleration to acceleration at late times. Furthermore, we confront the scenarios with data from Supernovae type Ia (SN Ia) and Hubble function observations, showing that the agreement is excellent. Moreover, we perform a dynamical system analysis and we show that there exist asymptotic accelerating attractors, arisen from the logarithmic terms as well as from the viscosity, which in most cases correspond to a phantom late-time evolution. Finally, for some parameter regions we obtain a nearly de Sitter late-time attractor, which is a significant capability of the scenario since the dark energy, although dynamical, stabilizes at the cosmological constant value.

In this paper the focus is on inflationary dynamics in the context of Einstein Gauss-Bonnet gravitational theories. We investigate the implications of the slow-roll condition on the slow-roll indices and we investigate how the inflationary dynamical evolution is affected by the presence of the Gauss-Bonnet coupling to the scalar field. For exemplification of our analysis we investigate how the dynamics of inflationary cubic, quartic order and also exponential scalar potentials are affected by the non-trivial Gauss-Bonnet coupling to the scalar field. As we demonstrate it is possible to obtain a viable phenomenology compatible with the observational data, although the canonical scalar field theory with cubic and quartic order potentials does not yield phenomenologically acceptable results. In addition, with regard to the exponential potential example, the Einstein Gauss-Bonnet extension of the single canonical scalar field model has an inherent mechanism that can trigger the graceful exit from inflation. Furthermore we introduce a bottom-up reconstruction technique, in the context of which by fixing the tensor-to-scalar ratio and the Hubble rate as a function of the $e$-foldings number, one is capable of reproducing the Einstein Gauss-Bonnet theory which generates the aforementioned quantities. We illustrate how the method works by using some relatively simple examples.

In this letter we will study the cosmological dynamical system of an $f(R)$ gravity in the presence of a canonical scalar field $\phi$ with an exponential potential, by constructing the dynamical system in a way that it is render autonomous. This feature is controlled by a single variable $m$, which when it is constant, the dynamical system is autonomous. We focus on the $m=0$ case which, as we demonstrate by using a numerical analysis approach, leads to an unstable de Sitter attractor, which occurs after $N\sim 60$ $e$-foldings. This instability can be viewed as a graceful exit from inflation, which is inherent to the dynamics of de Sitter attractors.

We investigate in detail the implications of the constant-roll condition on the inflationary era of a scalar field coupled to a teleparallel $f(T)$ gravity. The resulting cosmological equations constitute a reconstruction technique which enables us to find either the $f(T)$ gravity which corresponds to a given cosmological evolution, or the Hubble rate of the cosmological evolution generated by a fixed $f(T)$ gravity. We also analyze in some detail the phase space of the constant-roll teleparallel gravity and we discuss the physical significance of the resulting fixed points and trajectories. Also we calculate the observational indices of a theory with given $f(T)$ gravity, and we discuss all the implications of the constant-roll condition on these. As we demonstrate, the resulting theory can be compatible with the current observational data, for a wide range of values of the free parameters of the theory.

In this work we study certain classes of $F(R,{\cal G})$ gravity which have appealing phenomenological features, with respect to the successful realization of the dark energy and of the inflationary era. Particularly, we discuss the general formalism and we demonstrate how several inflationary and dark energy evolutions can be described in the context of $F(R,{\cal G})$ gravity. Also we propose a unified model, in the context of which the early and late-time dynamics are controlled by the $F(R,{\cal G})$ gravity, thus producing inflation and the dark energy era, while the intermediate era is approximately identical with standard Einstein-Hilbert gravity. Also we calculate the power spectrum of the primordial curvature perturbations corresponding to the unified $F(R,{\cal G})$ gravity model we propose, which as we demonstrate is nearly scale invariant and compatible with the latest observational data constraints.

In this work we shall investigate the occurrence of future cosmological finite-time singularities in the dynamical system corresponding to two cosmological theories, namely that of vacuum $f(R)$ gravity and that of three fluids. The vacuum $f(R)$ gravity is an example for which the variables we will choose to quantify the phase space dynamics, do not necessarily blow-up near a cosmological singularity. After appropriately choosing the variables, we shall investigate the behavior of the corresponding dynamical system near some types of cosmological finite-time singularities, for some limiting cases in which we can produce analytic solutions for the dynamical variables. The most interesting case from both a mathematical and physical point of view, is the Big Rip case, and particularly in the limiting case of a very strong singularity. The physically appealing outcome is that the resulting non-autonomous dynamical system is attracted asymptotically to an accelerating attractor solution, with equation of state parameter $w_{eff}=-1$. Our analytic results, show that an extremely strong Big Rip singularity in vacuum $f(R)$ gravity theories is always related to an accelerating solution, or tends to acceleration. The converse statement though may not be true. The second cosmology we shall study is a multifluid cosmology, consisting of three fluids, the interacting dark matter and dark energy fluids, and the baryonic fluid. By appropriately choosing the variables, we will show that the dynamical system can become an autonomous polynomial dynamical system, in which case, by using a dominant balance analysis, we shall investigate the occurrence of finite-time singularities in this system.

In this work we study the occurrence of finite-time cosmological singularities in a cosmological system comprising from three fluids. Particularly, the system contains two dark fluids, namely that of dark energy and dark matter, which are interacting, and of a non-interacting baryonic fluid. For the study we adopt the phase space approach by constructing the cosmological dynamical system in such a way so that it rendered to be an autonomous polynomial dynamical system, and in order to achieve this, we appropriately choose the variables of the dynamical system. By employing a rigid mathematical framework, that of dominant balances analysis, we demonstrate that there exist non-singular solutions of the dynamical system, which correspond to a general set of initial conditions, which proves that no Big Rip or Type III finite-time singularities occur in this LQC multifluid dynamical system. Thus the new feature of this work is that we are able to do this using an analytic technique instead of adopting a numerical approach. In addition, we perform a fixed point analysis of the cosmological dynamical system, and we examine the behavior of the total effective equation of state parameter, at the fixed points, as a function of the free parameters of the system. Finally, we investigate the phenomenological implications of the dark energy equation of state which we assumed that it governs the dark energy fluid.

In this work we construct a bottom-up reconstruction technique for Loop Quantum Cosmology scalar-tensor theories, from the observational indices. Particularly, the reconstruction technique is based on fixing the functional form of the scalar-to-tensor ratio as a function of the $e$-foldings number. The aim of the technique is to realize viable inflationary scenarios, and the only assumption that must hold true in order for the reconstruction technique to work is that the dynamical evolution of the scalar field obeys the slow-roll conditions. We shall use two functional forms for the scalar-to-tensor ratio, one of which corresponds to a popular inflationary class of models, the $\alpha$-attractors. For the latter, we shall calculate the leading order behavior of the spectral index and we shall demonstrate that the resulting inflationary theory is viable and compatible with the latest Planck and BICEP2/Keck-Array data. In addition, we shall find the classical limit of the theory, and as we demonstrate, the Loop Quantum Cosmology corrected theory and the classical theory are identical at leading order in the perturbative expansion quantified by the parameter $\rho_c$, which is the critical density of the quantum theory. Finally, by using the formalism of slow-roll scalar-tensor Loop Quantum Cosmology, we shall investigate how several inflationary potentials can be realized by the quantum theory, and we shall calculate directly the slow-roll indices and the corresponding observational indices. In addition, the $f(R)$ gravity frame picture is presented.

In this paper we investigate the inflationary dynamics of an $f(R)$ gravity in the presence of a canonical scalar field. We specifically choose the cosmological evolution to be a quasi-de Sitter evolution and also the $f(R)$ gravity is assumed to be a modified version of the $R^2$ gravity. We investigate which scalar field potential can produce the quasi-de Sitter evolution for the choice of the $f(R)$ gravity we made, and also we study in detail the inflationary dynamics of the resulting theory. As we demonstrate, the spectral index is identical to the one corresponding to the ordinary $R^2$ gravity, while the scalar-to-tensor ratio is found to be smaller than the $R^2$ inflation one consequently, compatibility with both the Planck 2015 and BICEP2/Keck-Array data is achieved.

In this paper we shall consider an exponential inflationary model in the context of vacuum $F(R)$ gravity. By using well-known reconstruction techniques, we shall investigate which $F(R)$ gravity can realize the exponential inflation scenario at leading order in terms of the scalar curvature, and we shall calculate the slow-roll indices and the corresponding observational indices, in the context of slow-roll inflation. We also provide some general formulas of the slow-roll and the corresponding observational indices in terms of the $e$-foldings number. In addition, for the calculation of the slow-roll and of the observational indices, we shall consider quite general formulas for which the assumption that all the slow-roll indices are much smaller than unity, is not necessary to hold true. Finally, we investigate the phenomenological viability of the model by comparing it with the latest Planck and BICEP2/Keck-Array observational data. As we demonstrate, the model is compatible with the current observational data for a wide range of the free parameters of the model.

By using a bottom-up reconstruction technique for non-minimally coupled scalar-tensor theories, we realize the Einstein frame attractor cosmologies in the $\Omega (\phi)$-Jordan frame. For our approach, what is needed for the reconstruction method to work is the functional form of the non-minimal coupling $\Omega(\phi)$ and of the scalar-to-tensor ratio, and also the assumption of the slow-roll inflation in the $\Omega (\phi)$-Jordan frame. By appropriately choosing the scalar-to-tensor ratio, we demonstrate that the observational indices of the attractor cosmologies can be realized directly in the $\Omega (\phi)$-Jordan frame. We investigate the special conditions that are required to hold true in for this realization to occur, and we provide the analytic form of the potential in the $\Omega (\phi)$-Jordan frame. Also, by performing a conformal transformation, we find the corresponding Einstein frame canonical scalar-tensor theory, and we calculate in detail the corresponding observational indices. The result indicates that although the spectral index of the primordial curvature perturbations is the same in the Jordan and Einstein frames, at leading order in the $e$-foldings number, the scalar-to-tensor ratio differs. We discuss the possible reasons behind this discrepancy, and we argue that the difference is due to some approximation we performed to the functional form of the potential in the Einstein frame, in order to obtain analytical results, and also due to the difference in the definition of the $e$-foldings number in the two frames, which is also pointed out in the related literature. Finally, we find the $F(R)$ gravity corresponding to the Einstein frame canonical scalar-tensor theory.

In this work, we extend the bottom-up reconstruction framework of $F(R)$ gravity to other modified gravities, and in particular for $f(\phi)R$ and mimetic $F(R)$ gravities. We investigate which are the important conditions in order for the method to work, and we study several viable cosmological evolutions, focusing on the inflationary era. Particularly, for the $f(\phi)R$ theory case, we specify the functional form of the Hubble rate and of the scalar-to-tensor ratio as a function of the $e$-foldings number and accordingly, the rest of the physical quantities and also the slow-roll and the corresponding observational indices can be calculated. The same method is applied in the mimetic $F(R)$ gravity case, and in both cases we thoroughly analyze the resulting free parameter space, in order to show that the viability of the models presented is guaranteed and secondly that there is a wide range of values of the free parameters for which the viability of the models occurs. In addition, the reconstruction method is also studied in the context of mimetic $F(R)=R$ gravity. As we demonstrate, the resulting theory is viable, and also in this case, only the scalar-to-tensor ratio needs to be specified, since the rest follow from this condition. Finally, we discuss in brief how the reconstruction method could function for other modified gravities.

In this work we use the Loop Quantum Cosmology modified scalar-tensor reconstruction techniques in order to investigate how bouncing and inflationary cosmologies can be realized. With regard to the inflationary cosmologies, we shall be interested in realizing the intermediate inflation and the Type IV singular inflation, while with regard to bouncing cosmologies, we shall realize the superbounce and the symmetric bounce. In all the cases, we shall find the kinetic term of the LQC holonomy corrected scalar-tensor theory and the corresponding scalar potential. In addition, we shall include a study of the effective equation of state, emphasizing at the early and late time eras. As we demonstrate, in some cases it is possible to have a nearly de Sitter equation of state at the late-time era, a result that could be interpreted as the description of a late-time acceleration era. Also, in all cases we shall examine the dynamical stability of the LQC holonomy corrected scalar-tensor theory, and we shall confront the results with those coming from the corresponding classical dynamical stability theory. The most appealing cosmological scenario is that of a Type IV singular inflationary scenario, in which the singularity may occur at the late-time era. As we demonstrate, for this model, during the dark energy era, a transition from non-phantom to a phantom dark energy era occurs.

In this work we investigate which Loop Quantum Cosmology corrected Gauss-Bonnet $F(\mathcal{G})$ gravity can realize two singular cosmological scenarios, the intermediate inflation and the singular bounce scenarios. The intermediate inflation scenario has a Type III sudden singularity at $t=0$, while the singular bounce has a soft Type IV singularity. By using perturbative techniques, we find the holonomy corrected $F(\mathcal{G})$ gravities that generate at leading order the aforementioned cosmologies and we also argue that the effect of the holonomy corrections is minor to the power spectrum of the primordial curvature perturbations of the classical theory.

In this work, we shall incorporate a stiff era in the Universe's evolution in the context of $F(R)$ gravity. After deriving the vacuum $F(R)$ gravity, which may realize a stiff evolution, we combine the stiff $F(R)$ gravity with an $R^2$ model, and we construct a qualitative model for the inflationary and stiff era, with the latter commencing after the end of the inflationary era. We assume that the baryogenesis occurs during the stiff era, and we calculate the baryon to entropy ratio, which effectively constraints the functional form of the stiff $F(R)$ gravity. Further constraints on the stiff $F(R)$ gravity may come from the primordial gravitational waves, and particularly their scalar mode, which is characteristic of the $F(R)$ gravity theory. The stiff era presence does not contradict the standard cosmology era, namely, inflation, and the radiation-matter domination eras. Furthermore, we investigate which $F(R)$ gravity may realize a dust and stiff matter dominated Einstein-Hilbert evolution.

In this paper we shall analyze the $f(\mathcal{G})$ gravity phase space, in the case that the corresponding dynamical system is autonomous. In order to make the dynamical system autonomous, we shall appropriately choose the independent variables, and we shall analyze the evolution of the variables numerically, emphasizing on the inflationary attractors. As we demonstrate, the dynamical system has only one de Sitter fixed point, which is unstable, with the instability being traced in one of the independent variables. This result holds true both in the presence and in the absence of matter and radiation perfect fluids. We argue that this instability could loosely be viewed as an indication of graceful exit in the $f(\mathcal{G})$ theory of gravity.

We present a class of inflationary potentials which are invariant under a special symmetry, which depends on the parameters of the models. As we show, in certain limiting cases, the inverse symmetric potentials are qualitatively similar to the $\alpha$-attractors models, since the resulting observational indices are identical. However, there are some quantitative differences which we discuss in some detail. As we show, some inverse symmetric models always yield results compatible with observations, but this strongly depends on the asymptotic form of the potential at large $e$-folding numbers. In fact when the limiting functional form is identical to the one corresponding to the $\alpha$-attractors models, the compatibility with the observations is guaranteed. Also we find the relation of the inverse symmetric models with the Starobinsky model and we highlight the differences. In addition, an alternative inverse symmetric model is studied and as we show, not all the inverse symmetric models are viable. Moreover, we study the corresponding $F(R)$ gravity theory and we show that the Jordan frame theory belongs to the $R^2$ attractor class of models. Finally we discuss a non-minimally coupled theory and we show that the attractor behavior occurs in this case too.

In this work we shall investigate the cosmological dynamical system of $f(R)$ gravity, by constructing it in such a way so that it is rendered autonomous. We shall study the vacuum $f(R)$ gravity case, but also the case that matter and radiation perfect fluids are present along with the $f(R)$ gravity. The dynamical system is constructed in such a way so that the time-dependence of the system is contained in a single parameter which depends on the Hubble rate and it's second derivative. The autonomous structure of the dynamical system is achieved when this parameter is constant, therefore we focus on these cases. For the vacuum $f(R)$ case, we investigate two cases with the first leading to a stable de Sitter attractor fixed point but also to an unstable de Sitter fixed point, and the second is related to a matter dominated era. The stable de Sitter attractor is also found for the $f(R)$ gravity in the presence of matter and radiation perfect fluids. In all the cases we performed a detailed numerical analysis of the dynamical system and we also investigate in detail the stability of the resulting fixed points. Also, we present an exceptional feature of the $R^2$ gravity model, in the absence of perfect fluids. Finally, we investigate what is the approximate form of the $f(R)$ gravities near the stable and the unstable de Sitter fixed points.

We propose two new versions of ghost-free generalized $F(R)$ gravity with Lagrange multiplier constraint. The first version of such theory for a particular degenerate choice of the Lagrange multiplier, corresponds to mimetic $F(R)$ gravity. The second version of such theory is just the Jordan frame description of mimetic gravity with potential. As we demonstrate, it is possible to realize several cosmological scenarios in such theory. In particulary, de Sitter solutions may also be found.

We propose the study of constant-roll inflation in $F(R)$ gravity. We use two different approaches, one that relates an $F(R)$ gravity to well known scalar models of constant-roll and a second that examines directly the constant-roll condition in $F(R)$ gravity. With regards to the first approach, by using well known techniques, we find the $F(R)$ gravity which realizes a given constant-roll evolution in the scalar-tensor theory. We also perform a conformal transformation in the resulting $F(R)$ gravity and we find the Einstein frame counterpart theory. As we demonstrate, the resulting scalar potential is different in comparison to the original scalar constant-roll case, and the same applies for the corresponding observational indices. Moreover, we discuss how cosmological evolutions that can realize constant-roll to constant-roll eras transitions in the scalar-tensor description, can be realized by vacuum $F(R)$ gravity. With regards to the second approach, we examine directly the effects of the constant-roll condition on the inflationary dynamics of vacuum $F(R)$ gravity. We present in detail the formalism of constant-roll $F(R)$ gravity inflationary dynamics and we discuss how the inflationary indices become in this case. We use two well known $F(R)$ gravities in order to illustrate our findings, the $R^2$ model and a power-law $F(R)$ gravity in vacuum. As we demonstrate, in both cases the parameter space is enlarged in comparison to the slow-roll counterparts of the models, and in effect, the models can also be compatible with the observational data. Finally, we briefly address the graceful exit issue.

In this work, we introduce a bottom-up $F(R)$ gravity reconstruction technique, in which we fix the observational indices and we seek for the $F(R)$ gravity which may realize them. Particularly, as an exemplification of our method, we shall assume that the scalar-to-tensor ratio has a specific form, and from it we shall reconstruct the $F(R)$ gravity that may realize it, focusing on special values of the parameters in order to obtain analytical results. The observational indices we study are compatible with the latest observational data, and we discuss how the functional form of the observational indices may affect the viability of the model.

In this work, we investigate how a smooth transition from a constant-roll to a slow-roll inflationary era may be realized, in the context of a canonical scalar field theory. We study in some detail the dynamical evolution of the cosmological system, and we investigate whether a stable attractor exists, both numerically and analytically. We also investigate the slow-roll era and as we demonstrate, partially compatibility of the resulting scalar theory, may be achieved, with the potential of the latter belonging to a class of modular inflationary potentials. The novel features of the constant-roll to slow-roll transition which we achieved, are firstly that the it is not compelling for the slow-roll era to last for $N\sim 50-60$ $e$-foldings, but it may last for a smaller number of $e$-foldings, since some $e$-foldings may occur during the constant-roll era. Secondly, when the slow-roll era occurs after the constant-roll era, the graceful exit from inflation may occur, a feature absent in the constant-roll scenario, due to the stability properties of the final attractor in the constant-roll case.

In this paper we investigate how to describe in a unified way a constant-roll inflationary era with a dark energy era, by using the theoretical framework of $F(R)$ gravity. To this end, we introduce some classes of appropriately chosen $F(R)$ gravity models, and we examine in detail how the unification of early and late-time acceleration eras can be achieved. We study in detail the inflationary era, and as we demonstrate it is possible to achieve a viable inflationary era, for which the spectral index of primordial curvature perturbations and the scalar-to-tensor ratio can be compatible with the latest observational data. In addition, the graceful exit issue is briefly discussed for a class of models. Finally, we discuss the dark energy oscillations issue, and we investigate which model from one of the classes we introduced, can produce oscillations with the smallest amplitude.

We study in detail the phase space of a Friedmann-Robertson-Walker Universe filled with various cosmological fluids which may or may not interact. We use various expressions for the equation of state, and we analyze the physical significance of the resulting fixed points. In addition we discuss the effects of the stability or an instability of some fixed points. Moreover we study an interesting phenomenological scenario for which there is an oscillating interaction between the dark energy and dark matter fluid. As we demonstrate, in the context of the model we use, at early times the interaction is negligible and it starts to grow as the cosmic time approaches the late-time era. Also the cosmological dynamical system is split into two distinct dynamical systems which have two distinct de Sitter fixed points, with the early-time de Sitter point being unstable. This framework gives an explicit example of the unification of the early-time with late-time acceleration. Finally, we discuss in some detail the physical interpretation of the various models we present in this work.

We study the gravitational baryogenesis mechanism for generating baryon asymmetry in the context of running vacuum models. Regardless if these models can produce a viable cosmological evolution, we demonstrate that they produce a non-zero baryon-to-entropy ratio even if the universe is filled with conformal matter. This is a sound difference between the running vacuum gravitational baryogenesis and the Einstein-Hilbert one, since in the latter case, the predicted baryon-to-entropy ratio is zero. We consider two well known and most used running vacuum models and show that the resulting baryon-to-entropy ratio is compatible with the observational data. Moreover, we also show that the mechanism of gravitational baryogenesis may constrain the running vacuum models.

In this paper we study canonical scalar field models with a varying second slow-roll parameter, that allow transitions between constant-roll eras. In the models with two constant-roll eras it is possible to avoid fine-tunings in the initial conditions of the scalar field. We mainly focus on the stability of the resulting solutions and we also investigate if these solutions are attractors of the cosmological system. We shall calculate the resulting scalar potential and by using a numerical approach, we examine the stability and attractor properties of the solutions. As we show, the first constant-roll era is dynamically unstable towards linear perturbations and the cosmological system is driven by the attractor solution to the final constant-roll era. As we demonstrate, it is possible to have a nearly-scale invariant power spectrum of primordial curvature perturbations in some cases, however this is strongly model dependent and depends on the rate of the final constant-roll era. Finally, we present in brief the essential features of a model that allows oscillations between constant-roll eras.

In this work we address the reheating issue in the context of $F(R)$ gravity, for theories that the inflationary era does not obey the slow-roll condition but the constant-roll condition is assumed. As it is known, the reheating era takes place after the end of the inflationary era, so we investigate the implications of a constant-roll inflation era on the reheating era. We quantify our considerations by calculating the reheating temperature for the constant-roll $R^2$ model and we compare to the standard reheating temperature in the context of $F(R)$ gravity. As we demonstrate, the new reheating temperature may differ from the standard one, and in addition we show how the reheating era may restrict the constant-roll era by constraining the constant-roll parameter.

In this paper we investigate the implications of having a varying second slow-roll index on the canonical scalar field inflationary dynamics. We shall be interested in cases that the second slow-roll can take small values and correspondingly large values, for limiting cases of the function that quantifies the variation of the second slow-roll index. As we demonstrate, this can naturally introduce a smooth transition between slow-roll and constant-roll eras. We discuss the theoretical implications of the mechanism we introduce and we use various illustrative examples in order to better understand the new features that the varying second slow-roll index introduces. In the examples we will present, the second slow-roll index has exponential dependence on the scalar field, and in one of these cases, the slow-roll era corresponds to a type of $\alpha$-attractor inflation. Finally, we briefly discuss how the combination of slow-roll and constant-roll may lead to non-Gaussianities in the primordial perturbations.

An alternative to the Big Bang cosmologies is obtained by the Big Bounce cosmologies. In this paper, we study a bounce cosmology with a Type IV singularity occurring at the bouncing point, in the context of $F(R)$ modified gravity. We investigate the evolution of the Hubble radius and we examine the issue of primordial cosmological perturbations in detail. As we demonstrate, for the singular bounce, the primordial perturbations originating from the cosmological era near the bounce, do not produce a scale invariant spectrum and also the short wavelength modes, after these exit the horizon, do not freeze, but grow linearly with time. After presenting the cosmological perturbations study, we discuss the viability of the singular bounce model, and our results indicate that the singular bounce must be combined with another cosmological scenario, or should be modified appropriately, in order that it leads to a viable cosmology. The study of the slow-roll parameters leads to the same result, indicating the singular bounce theory is unstable at the singularity point, for certain values of the parameters. We also conformally transform the Jordan frame singular bounce, and as we demonstrate, the Einstein frame metric leads to a Big Rip singularity. Therefore, the Type IV singularity in the Jordan frame, becomes a Big Rip singularity in the Einstein frame. Finally, we briefly study a generalized singular cosmological model, which contains two Type IV singularities, with quite appealing features.

In this paper we study the occurrence of accelerating universe versus decelerating universe between the F(R) gravity frame (Jordan frame) and non-minimally coupled scalar field theory frame, and the minimally coupled scalar field theory frame (Einstein frame) for various models. As we show, if acceleration is imposed in one frame, it will not necessarily correspond to an accelerating metric when transformed in another frame. As we will demonstrate, this issue is model and frame-dependent but it seems there is no general scheme which permits to classify such cases.

We investigate how baryogenesis can occur by the presence of an $f(T)$-related gravitational term. We study various cases of $f(T)$ gravity and we discuss in detail the effect of the novel terms on the baryon-to-entropy ratio. Additionally, we study the constraints imposed by the observational values of the baryon-to-entropy ratio and we discuss how more generalized cosmologies can contribute successfully, in a viable and consistent way, in the gravitational baryogenesis mechanism.

In this paper we study some classes of $\alpha$-attractors models in the Jordan frame and we find the corresponding $F(R)$ gravity theory. We study analytically the problem at leading order and we investigate whether the attractor picture persists in the $F(R)$ gravity equivalent theory. As we show, if the slow-roll conditions are assumed in the Jordan frame, the spectral index of primordial curvature perturbations and the scalar-to-tensor ratio are identical to the corresponding observational indices of the $R^2$ model, a result which indicates that the attractor property is also found in the corresponding $F(R)$ gravity theories of the $\alpha$-attractors models. Moreover, implicit and approximate forms of the $F(R)$ gravity inflationary attractors are found.

In this paper, we demonstrate that a unified description of early and late-time acceleration is possible in the context of mimetic $F(R)$ gravity. We study the inflationary era in detail and demonstrate that it can be realized even in mimetic $F(R)$ gravity where traditional $F(R)$ gravity fails to describe the inflation. By using standard methods we calculated the spectral index of primordial curvature perturbations and the scalar-to-tensor ratio. We use two $F(R)$ gravity models and as it turns out, for both the models under study the observational indices are compatible with both the latest Planck and the BICEP2/Keck array data. Finally, the graceful exit from inflation is guaranteed by the existence of growing curvature perturbations when the slow-roll era ends.

Loop Quantum Cosmology is an appealing quantum completion of classical cosmology, which brings along various theoretical features which in many cases offer remedy or modify various classical cosmology aspects. In this paper we address the gravitational baryogenesis mechanism in the context of Loop Quantum Cosmology. As we demonstrate, when Loop Quantum Cosmology effects are taken into account in the resulting Friedmann equations for a flat Friedmann-Robertson-Walker Universe, then even for a radiation dominated Universe, the predicted baryon-to-entropy ratio from the gravitational baryogenesis mechanism is non-zero, in contrast to the Einstein-Hilbert case, in which case the baryon-to-entropy ratio is zero. We also discuss various other cases apart from the radiation domination case, and we discuss how the baryon-to-entropy ratio is affected from the parameters of the quantum theory. In addition, we use illustrative exact solutions of Loop Quantum Cosmology and we investigate under which circumstances the baryon-to-entropy ratio can be compatible with the observational constraints.

We extend the Loop Quantum Cosmology matter bounce scenario in order to include a dark energy era, which ends abruptly at a Rip singularity where the scale factor and the Hubble rate diverge. In the "deformed matter bounce scenario", the Universe is contracting from an initial non-causal matter dominated era until it reaches a minimal radius. After that it expands in a decelerating way, until at late times, where it expands in an accelerating way, thus the model is described by a dark energy era that follows the matter dominated era. Depending on the choice of the free parameters of the model, the dark energy era is quintessential like which follows the matter domination era, and eventually it crosses the phantom divide line and becomes phantom. At the end of the dark energy era, a Rip singularity exists, where the scale factor and Hubble rate diverge, however the physical system cannot reach the singularity, since the effective energy density and pressure become complex. This indicates two things, firstly that the ordinary Loop Quantum Cosmology matter bounce evolution stops, thus ending the infinite repetition of the ordinary matter bounce scenario. Secondly, the fact that both the pressure and the density become complex indicates probably that the description of the cosmic evolution within the theoretical context of Loop Quantum Cosmology, ceases to describe the physics of the system and possibly a more fundamental theory of quantum gravity is needed near the would be Rip singularity. We describe the qualitative features of the model and we also investigate how this cosmology could be realized by a viscous fluid and scalar field in the context of LQC.

In this paper we address the problem of dark energy oscillations in the context of mimetic $F(R)$ gravity with potential. The issue of dark energy oscillations can be a problem in some models of ordinary $F(R)$ gravity and a remedy that can make the oscillations milder is to introduce additional modifications in the functional form of the $F(R)$ gravity. As we demonstrate the power-law modifications are not necessary in the mimetic $F(R)$ case, and by appropriately choosing the mimetic potential and the Lagrange multiplier, it is possible to make the oscillations almost to vanish at the end of the matter domination era and during the late-time acceleration era. We examine the behavior of the dark energy equation of state parameter and of the total effective equation of state parameter as functions of the redshift and we compare the resulting picture with the ordinary $F(R)$ gravity case. As we also show, the present day values of the dark energy equation of state parameter and of the total effective equation of state parameter are in better agreement with the observational data, in comparison to the ordinary $F(R)$ gravity case. Finally, we study the evolution of the growth factor as a function of the redshift for all the mimetic models we shall use.

We extend the scalar-tensor reconstruction techniques for classical cosmology frameworks, in the context of loop quantum cosmology. After presenting in some detail how the equations are generalized in the loop quantum cosmology case, we discuss which new features and limitations does the quantum framework brings along, and we use various illustrative examples in order to demonstrate how the method works. As we show the energy density has two different classes of solutions, and one of these yields the correct classical limit while the second captures the quantum phenomena. We study in detail the scalar tensor reconstruction method for both these solutions. Also we discuss some scenarios for which the Hubble rate becomes unbounded at finite time, which corresponds for example in a case that a Big Rip occurs. As we show this issue is non-trivial and we discuss how this case should be treated in a consistent way. Finally, we investigate how the classical stability conditions for the scalar-tensor solutions are generalized in the loop quantum framework.

In this letter we study some variant forms of gravitational baryogenesis by using higher order terms containing the partial derivative of the Gauss-Bonnet scalar coupled to the baryonic current. This scenario extends the well known theory that uses a similar coupling between the Ricci scalar and the baryonic current. One appealing feature of the scenario we study is that the predicted baryon asymmetry during a radiation domination era is non-zero. We calculate the baryon to entropy ratio for the Gauss-Bonnet term and by using the observational constraints we investigate which are the allowed forms of the $R+F(\mathcal{G})$ gravity controlling the evolution. Also we briefly discuss some alternative higher order terms that can generate a non-zero baryon asymmetry, even in the conformal invariance limit.

In this paper we demonstrate that with vacuum $F(G)$ gravity it is possible to describe the unification of late and early-time acceleration eras with the radiation and matter domination era. The Hubble rate of the unified evolution contains two mild singularities, so called Type IV singularities, and the evolution itself has some appealing features, such as the existence of a deceleration-acceleration transition at late times. We also address quantitatively a fundamental question related to modified gravity models description of cosmological evolution: Is it possible for all modified gravity descriptions of our Universe evolution, to produce a nearly scale invariant spectrum of primordial curvature perturbations? As we demonstrate, the answer for the $F(G)$ description is no, since the resulting power spectrum is not scale invariant, in contrast to the $F(R)$ description studied in the literature. Therefore, although the cosmological evolution can be realized in the context of vacuum $F(G)$ gravity, the evolution is not compatible with the observational data, in contrast to the $F(R)$ gravity description of the same cosmological evolution.

We present some cosmological models which unify the late and early-time acceleration eras with the radiation and the matter domination era, and we realize the cosmological models by using the theoretical framework of $F(R)$ gravity. Particularly, the first model unifies the late and early-time acceleration with the matter domination era, and the second model unifies all the evolution eras of our Universe. The two models are described in the same way at early and late times, and only the intermediate stages of the evolution have some differences. Each cosmological model contains two Type IV singularities which are chosen to occur one at the end of the inflationary era and one at the end of the matter domination era. The cosmological models at early times are approximately identical to the $R^2$ inflation model, so these describe a slow-roll inflationary era which ends when the slow-roll parameters become of order one. The inflationary era is followed by the radiation era and after that the matter domination era follows, which lasts until the second Type IV singularity, and then the late-time acceleration era follows. The models produce a nearly scale invariant power spectrum of primordial curvature perturbations and a scalar-to-tensor ratio which are compatible with the most recent observational data and it seems that the deceleration acceleration transition is crucially affected by the presence of the second Type IV singularity which occurs at the end of the matter domination era. We perform an analysis of the Hubble horizon evolution, and in addition, we investigate which $F(R)$ gravity can successfully realize each of the four cosmological epochs.

We combine the unimodular gravity and mimetic gravity theories into a unified theoretical framework, which is proposed to provide a suggestive proposal for a framework that may assist in the discussion and solution search of the cosmological constant problem and of the dark matter issue. After providing the formulation of the unimodular mimetic gravity and investigating all the new features that the vacuum unimodular gravity implies, by using the underlying reconstruction method, we realize some well known cosmological evolutions, with some of these being exotic for the ordinary Einstein-Hilbert gravity. Specifically we provide the vacuum unimodular mimetic gravity description of the de Sitter cosmology and of the perfect fluid with constant equation of state cosmology. As we demonstrate, these cosmologies can be realized by vacuum mimetic unimodular gravity, without the existence of any matter fluid source. Moreover, we investigate how cosmologically viable cosmologies, which are compatible with the recent observational data, can be realized by the vacuum unimodular mimetic gravity. Since in some cases, the graceful exit from inflation problem might exist, we provide a qualitative description of the mechanism that can potentially generate the graceful exit from inflation in these theories, by searching for the unstable de Sitter solutions in the context of unimodular mimetic theories of gravity.

We propose a covariant ghost-free unimodular $F(R)$ gravity theory, which contains a three-form field and study its structure using the analogy of the proposed theory with a quantum system which describes a charged particle in uniform magnetic field. Newton's law in non-covariant unimodular $F(R)$ gravity as well as in unimodular Einstein gravity is derived and it is shown to be just the same as in General Relativity. The derivation of Newton's law in covariant unimodular $F(R)$ gravity shows that it is modified precisely in the same way as in the ordinary $F(R)$ theory. We also demonstrate that the cosmology of a Friedmann-Robertson-Walker background, is equivalent in the non-covariant and covariant formulations of unimodular $F(R)$ theory.

In this paper we investigate how to realize various quite well known cosmological bouncing models in the context of the recently developed unimodular $F(R)$ gravity. Particularly, we shall study the matter bounce scenario, the singular bounce, the superbounce and a symmetric bounce scenario. We present the behavior of the Hubble radius for each of the bouncing models we shall take into account and we investigate which era of the bouncing model is responsible for the cosmological perturbations. As we shall demonstrate, the various bouncing models do not behave in the same way, so the cosmological perturbations for each model may correspond to a different era, in comparison to other models. Also we present which unimodular $F(R)$ gravity realizes each model. We also show that Newton's law is not modified in the unimodular $F(R)$ gravity, which also is proven to be a ghost-free theory, and in addition we discuss the matter stability issue. Finally, we demonstrate how it is possible to solve a cosmological constant problem in the context of unimodular $F(R)$ gravity.

In this brief note we investigate the conditions under which a Schwarzschild de Sitter black hole spacetime is a solution of the mimetic $F(R)$ gravity with Lagrange multiplier and potential. As we demonstrate, the resulting mimetic $F(R)$ gravity is a slight modification of the ordinary $F(R)$ gravity case, however the resulting perturbation equations are not in all cases identical to the ordinary $F(R)$ gravity case. In the latter case, the perturbation equations are identical to the ones corresponding to the Reissner-Nordstr\"{o}m anti-de Sitter black hole.

We extend the formalism of the Einstein-Hilbert unimodular gravity in the context of modified $F(R)$ gravity. After appropriately modifying the Friedmann-Robertson-Walker metric in a way that it becomes compatible to the unimodular condition of having a constant metric determinant, we derive the equations of motion of the unimodular $F(R)$ gravity by using the metric formalism of modified gravity with Lagrange multiplier constraint. The resulting equations are studied in frames of reconstruction method, which enables us to realize various cosmological scenarios, which was impossible to realize in the standard Einstein-Hilbert unimodular gravity. Several unimodular $F(R)$ inflationary scenarios are presented, and in some cases, concordance with Planck and BICEP2 observational data can be achieved.

We study the finite time singularity correspondence between the Jordan and Einstein frames for various $F(R)$ gravity theories. Particularly we investigate the ordinary pure $F(R)$ gravity case and the unimodular $F(R)$ gravity cases, in the absence of any matter fluids. In the ordinary $F(R)$ gravity cases, by using specific illustrative examples, we show that it is possible to have various correspondences of finite time singularities, and in some cases it is possible a singular cosmology in one frame might be non-singular in the other frame. In the unimodular $F(R)$ gravity case, the unimodular constraint is affected from the conformal transformation, so this has an effect on the metric we choose. Moreover, we study the Einstein frame counterpart theory of the unimodular $F(R)$ gravity case, and we investigate the correspondences of the singularities in the two theories by considering specific illustrative examples. Finally, a brief dynamical system analysis is performed for the vacuum unimodular $F(R)$ gravity and we demonstrate how the dynamical system behaves near the future Big Rip singularity.

We propose the unimodular-mimetic $F(R)$ gravity theory, to resolve cosmological constant problem and dark matter problem in a unified geometric manner. We demonstrate that such a theory naturally admits accelerating universe evolution. Furthermore, we construct unimodular-mimetic $F(R)$ inflationary cosmological scenarios compatible with the Planck and BICEP2/Keck-Array observational data. We also address the graceful exit issue, which is guaranteed by the existence of unstable de Sitter vacua.

In this paper we demonstrate that in the context of mimetic $F(R)$ gravity with Lagrange multiplier, it is possible to realize cosmologies which are compatible with the recent BICEP2/Keck Array data. We provide some characteristic examples for which the predicted scalar to tensor ratio can be quite smaller in comparison to the upper limit imposed by the BICEP2/Keck Array observations.

We study mimetic $F(R)$ gravity with potential and Lagrange multiplier constraint. In the context of these theories, we introduce a reconstruction technique which enables us to realize arbitrary cosmologies, given the Hubble rate and an arbitrarily chosen $F(R)$ gravity. We exemplify our method by realizing cosmologies that are in concordance with current observations (Planck data) and also well known bouncing cosmologies. The attribute of our method is that the $F(R)$ gravity can be arbitrarily chosen, so we can have the appealing features of the mimetic approach combined with the known features of some $F(R)$ gravities, which unify early-time with late-time acceleration. Moreover, we study the existence and the stability of de Sitter points in the context of mimetic $F(R)$ gravity. In the case of unstable de Sitter points, it is demonstrated that graceful exit from inflation occurs. We also study the Einstein frame counterpart theory of the Jordan frame mimetic $F(R)$ gravity, we discuss the general properties of the theory and exemplify our analysis by studying a quite interesting from a phenomenological point of view, model with two scalar fields. We also calculate the observational indices of the two scalar field model, by using the two scalar field formalism. Furthermore, we extensively study the dynamical system that corresponds to the mimetic $F(R)$ gravity, by finding the fixed points and studying their stability. Finally, we modify our reconstruction method to function in the inverse way and thus yielding which $F(R)$ gravity can realize a specific cosmological evolution, given the mimetic potential and the Lagrange multiplier.

We investigate how the gravitational baryogenesis mechanism can potentially constrain the form of a Type IV singularity. Specifically, we study two different models with interesting phenomenology, that realize two distinct Type IV singularities, one occurring at the end of inflation and one during the radiation domination era or during the matter domination era. As we demonstrate, the Type IV singularities occurring at the matter domination era or during the radiation domination era, are constrained by the gravitational baryogenesis, in such a way so that these do not render the baryon to entropy ratio singular. Both the cosmological models we study cannot be realized in the context of ordinary Einstein-Hilbert gravity, and hence our work can only be realized in the context of $F(R)$ gravity and more generally in the context of modified gravity only.

The possibility to have singular accelerated evolution in the context of $F(R)$ bimetric gravity is investigated. Particularly, we study two singular models of cosmological evolution, one of which is a singular modified version of the Starobinsky $R^2$ inflation model. As we demonstrate, for both models in some cases, the slow-roll parameters become singular at the Type IV singularity, a fact that we interpret as a dynamical instability of the theory under study. This dynamically instability may be an indicator of graceful exit from inflation and we thoroughly discuss this scenario and the interpretation of the singular slow-roll parameters. Furthermore, it is demonstrated that for some versions of $F(R)$ bigravity, singular inflation is realized in consistent way so that inflationary indices are compatible with Planck data. Moreover, we study the late-time behavior of the two singular models and we show that the unified description of early and late-time acceleration can be achieved in the context of bimetric $F(R)$ gravity.

We study how a cosmological bounce with a Type IV singularity at the bouncing point, can be generated by a classical vacuum $F(G)$ gravity. We focus our investigation on the behavior of the vacuum $F(G)$ theory near the Type IV singular bouncing point and also we address the stability of the resulting solution, by treating the equations of motion as a dynamical system. In addition, we investigate how the scalar perturbations of the background metric evolve, emphasizing to cosmological times near the Type IV singular bouncing point. Finally, we also investigate which mimetic vacuum $F(G)$ gravity can describe the singular bounce cosmology.

Using Mukhanov-Chamseddine mimetic approach, we study $F(R)$ gravity with scalar potential and Lagrange multiplier constraint. As we demonstrate, for a given $F(R)$ gravity and for suitably chosen mimetic potential, it is possible to realize inflationary cosmology consistent with Planck observations. We also investigate the de Sitter solutions of the mimetic $F(R)$ theory and study the stability of the solutions, when these exist, towards linear perturbations, with the unstable solutions, which can provide a mechanism for graceful exit from inflation. Finally, we describe a reconstruction method which can yield the $F(R)$ gravity that can generate realistic inflationary cosmological evolution, given the mimetic potential and the Hubble rate.

Unlike crushing singularities, the so-called Type IV finite-time singularity offers the possibility that the Universe passes smoothly through it, without any catastrophic effects. Then the question is if the effects of a Type IV singularity can be detected in the process of cosmic evolution. In this paper we address this question in the context of $F(R)$ gravity. As we demonstrate, the effects of a Type IV singularity appear in the Hubble flow parameters, which determine the dynamical evolution of the cosmological system. So we study various inflation models incorporating a Type IV singularity, with the singularity occurring at the end of inflation. Particularly we study a toy model and a singular version of the $R^2$ gravity Hubble rate. As we evince, some of the Hubble flow parameters become singular at the singularity, an effect which indicates that at that point a dynamical instability occurs. This dynamical instability eventually indicates the graceful exit from inflation. We demonstrate that the toy model has an unstable de Sitter point at the singularity, so indeed graceful exit could be triggered. In the case of the singular inflation model, graceful exit proceeds in the standard way. In the case of the singular inflation model, we found various scenarios for singular evolution, most of which are compatible with observations, and only one leads to severe instabilities. We also compare the ordinary Starobinsky with the singular inflation model, and we point out the qualitative and quantitative differences. Finally, we study the late-time dynamics of the toy model and of the singular inflation model and we demonstrate that the unification of early and late-time acceleration can be achieved. We also show that it is possible to achieve late-time acceleration similar to the $\Lambda$-Cold Dark Matter model.

Using the reconstruction technique with an auxiliary field, we investigate which $F(R)$ gravities can produce the matter bounce cosmological solutions. Owing to the specific functional form of the matter bounce Hubble parameter, the reconstruction technique leads, after some simplifications, to the same Hubble parameter as in the matter bounce scenario. Focusing the study to the large and small cosmic time $t$ limits, we were able to find which $F(R)$ gravities can generate the matter bounce Hubble parameter. In the case of small cosmic time limit, which corresponds to large curvature values, the $F(R)$ gravity is $F(R)\sim R+\alpha R^2$, which is an inflation generating gravity, and at small curvature, or equivalently, large cosmic time, the $F(R)$ gravity generating the corresponding limit of the matter bounce Hubble parameter, is $F(R)\sim \frac{1}{R}$, a gravity known to produce late-time acceleration. Thus we have the physically appealing picture in which a Jordan frame $F(R)$ gravity that imitates the matter bounce solution at large and small curvatures, can generate Starobinsky inflation and late-time acceleration. Moreover, the scale factor corresponding to the reconstruction technique coincides almost completely to the matter bounce scenario scale factor, when considered in the aforementioned limiting curvature cases. This is scrutinized in detail, in order to examine the validity of the reconstruction method in these limiting cases, and according to our analysis, exact agreement is achieved.

We investigate in which cases a singular evolution with a singularity of Type IV, can be consistently incorporated in deformations of the $R^2$ inflationary potential. After demonstrating the difficulties that the single scalar field description is confronted with, we use a general two scalar fields model without other matter fluids, to describe the Type IV singular evolution, with one of the two scalar fields being canonical. By appropriately choosing the non-canonical scalar field, we show that the canonical scalar field corresponds to a potential that is nearly the $R^2$ inflation potential. If the Type IV singularity occurs at the end of inflation, the Universe's dynamical evolution near inflation is determined effectively by the canonical scalar field and at late-time the evolution is effectively determined by the non-canonical scalar. We also discuss the evolution of the Universe in terms of the effective equation of state and we show that the Type IV singularity, that occurs at the end of inflation, drives late-time acceleration. If however the singularity occurs at late-time, this might affect the inflationary era. We also investigate which Jordan frame pure $F(R)$ gravity corresponds to the nearly $R^2$ inflation scalar potentials we found. The stability of the solutions in the two scalar fields case is also studied and also we investigate how Type IV singularities can be incorporated in certain limiting cases of $R+R^p$ gravity in the Einstein frame. Finally, we briefly discuss a physical appealing scenario triggered by instabilities in the dynamical system that describes the evolution of the scalar fields.

We investigate the realization of two bouncing paradigms, namely of the superbounce and the loop quantum cosmological ekpyrosis, in the framework of various modified gravities. In particular, we focus on the $F(R)$, $F(G)$ and $F(T)$ gravities, and we reconstruct their specific subclasses which lead to such universe evolutions. These subclasses constitute from power laws, polynomials, or hypergeometric ansatzes, which can be approximated by power laws. The qualitative similarity of different effective gravities which realize the above two bouncing cosmologies, indicates to some universality lying behind such a bounce. Finally, performing a linear perturbation analysis, we show that the obtained solutions are conditionally or fully stable.

We demonstrate that finite time singularities of Type IV can be consistently incorporated in the Universe's cosmological evolution, either appearing in the inflationary era, or in the late-time regime. While using only one scalar field instabilities can in principle occur at the time of the phantom-divide crossing, when two fields are involved we are able to avoid such instabilities. Additionally, the two-field scalar-tensor theories prove to be able to offer a plethora of possible viable cosmological scenarios, at which various types of cosmological singularities can be realized. Amongst others, it is possible to describe inflation with the appearance of a Type IV singularity, and phantom late-time acceleration which ends in a Big Rip. Finally, for completeness, we also present the Type IV realization in the context of suitably reconstructed $F(R)$ gravity.

In this paper we propose and extensively study mimetic $f({\cal G})$ modified gravity models, with various scenarios of cosmological evolution, with or without extra matter fluids. The easiest formulation is based on the use of Lagrange multiplier constraint. In certain versions of this theory, it is possible to realize accelerated expansion of the Universe or even unified evolution which includes inflation with dark energy, and at the same time in the same theoretical framework, dark matter is described by the theory. This is achieved by the re-parametrization of the metric tensor, which introduces a new degree of freedom in the cosmological equations and leads to appearance of the mimetic "dark matter" component. In the context of mimetic $f({\cal G})$ theory, we also provide some quite general reconstruction schemes, which enable us to find which $f({\cal G})$ gravity generates a specific cosmological evolution. In addition, we also provide the general reconstruction technique for Lagrange multiplier $f({\cal G})$ gravity. All our results are accompanied by illustrative examples, with special emphasis on bouncing cosmologies.

We investigate how a Type IV future singularity can be included in the cosmological evolution of a well-known exponential model of inflation. In order to achieve this we use a two scalar field model, in the context of which the incorporation of the Type IV singularity can be consistently done. In the context of the exponential model we study, when a Type IV singularity is included in the evolution, an instability occurs in the slow-roll parameters, and in particular on the second slow-roll parameter. Particularly, if we abandon the slow-roll condition for both the scalars we shall use, then the most consistent description of the dynamics of the inflationary era is provided by the Hubble slow-roll parameters $\epsilon_H$ and $\eta_H$. Then, the second Hubble slow-roll parameter $\eta_H$, which measures the duration of the inflationary era, becomes singular at the point where the Type IV singularity is chosen to occur, while the Hubble slow-roll parameter $\epsilon_H$ is regular there. Therefore, this infinite singularity indicates that the occurrence of the finite-time singularity is responsible for the instability in the scalar field model we study. This sort of instability has it's imprint on the dynamical system that can be constructed from the cosmological equations, with the dynamical system being unstable. In addition, the instability due to the singularity mechanism we propose, is discussed in the context of other inflationary scalar potentials. Finally, we discuss the implications of such a singularity in the Hubble slow-roll parameters and we also critically discuss qualitatively, what implications could this effect have on the graceful exit problem of the exponential model.

We investigate which Jordan frame $F(R)$ gravity can describe a Type IV singular bouncing cosmological evolution, with special emphasis given near the point at which the Type IV singularity occurs. The cosmological bounce is chosen in such a way so that the bouncing point coincides exactly with Type IV singularity point. The stability of the resulting $F(R)$ gravity is examined and in addition, we study the Einstein frame scalar-tensor theory counterpart of the resulting Jordan frame $F(R)$ gravity. Also, by assuming that the Jordan frame metric is chosen in such a way so that, when conformally transformed in the Einstein frame, it yields a quasi de Sitter or de Sitter Friedmann-Robertson-Walker metric, we study the observational indexes which turn out to be consistent with Planck 2015 data in the case of the Einstein frame scalar theory. Finally, we study the behavior of the effective equation of state corresponding to the Type IV singular bounce and after we compare the resulting picture with other bouncing cosmologies, we critically discuss the implications of our analysis.

We study models with a generalized inhomogeneous equation of state fluids, in the context of singular inflation, focusing to so-called Type IV singular evolution. In the simplest case, this cosmological fluid is described by an equation of state with constant $w$, and therefore a direct modification of this constant $w$ fluid, is achieved by using a generalized form of an equation of state. We investigate from which models with generalized phenomenological equation of state, a Type IV singular inflation can be generated and what the phenomenological implications of this singularity would be. We support our results with illustrative examples and we also study the impact of the Type IV singularities on the slow-roll parameters and on the observational inflationary indices, showing the consistency with Planck mission results. The unification of singular inflation with singular dark energy era for specific generalized fluids is also proposed.

We provide a detailed quantitative description of singular inflation. Its close analogy with finite-time future singularity which is associated to dark energy era is described. Calling and classifying the singularities of such inflation as finite-time cosmological singularities we investigate their occurrence, with special emphasis on the Type IV singularity. The study is performed in the context of a general non-canonical scalar-tensor theory. In addition, the impact of finite time singularities on the slow-roll parameters is also investigated. Particularly, we study three cases, in which the singularity occurs during the inflationary era, at the end, and also we study the case that the singularity occurs much more later than inflation ends. Using the obtained slow-roll parameters, for each case, we calculate explicitly the spectral index of primordial curvature perturbations $n_s$, the associated running of the spectral index $a_s$ and of the scalar-to-tensor ratio $r$ and compare the resulting values to the Planck and BICEP2 data. As we demonstrate, in some cases corresponding to the Type IV singularity, there might be the possibility of agreement with the observational data, when the singularity occurs at the end, or after inflation. However, absolute concordance of all observational indices is not achieved. On the contrary, if the singularity occurs during the inflationary era, this is catastrophic for the theory, since the observational indices become divergent. We also show how a Type IV singularity may be consistently accommodated in the Universe's late time evolution.

As is known, in modified cosmological theories of gravity many of the cosmologies which could not be generated by standard Einstein gravity, can be consistently described by $F(R)$ theories. Using known reconstruction techniques, we investigate which $F(R)$ theories can lead to a Hubble parameter describing two types of cosmological bounces, the superbounce model, related to supergravity and non-supersymmetric models of contracting ekpyrosis and also the Loop Quantum Cosmology modified ekpyrotic model. Since our method is an approximate method, we investigate the problem at large and small curvatures. As we evince, both models yield power law reconstructed $F(R)$ gravities, with the most interesting new feature being that both lead to accelerating cosmologies, in the large curvature approximation. The mathematical properties of the some Friedmann-Robertson-Walker spacetimes $M$, that describe superbounce-like cosmologies are also pointed out, with regards to the group of curvature collineations $CC(M)$.

We provide an $F(R)$ gravity description of a $\Lambda$CDM bouncing model, without the need for matter fluids or for cosmological constant. As we explicitly demonstrate, the two cosmological eras that constitute the $\Lambda$CDM bouncing model, can be generated by $F(R)$ gravity which can lead to accelerating cosmologies. The resulting F(R) gravity has Einstein frame inflationary properties that have concordance to the latest Planck observational data. Both the $F(R)$ gravity stability properties are thoroughly investigated and also, the gravitational particle production, a feature necessary for the viability of the $\Lambda$CDM bounce scenario, is also addressed. As we will show, the $\Lambda$CDM bounce model can be successfully described by pure $F(R)$ gravity, with appealing phenomenological attributes, which we extensively discuss.

Using the reconstruction method, we investigate which $F(R)$ theories, with or without the presence of matter fluids, can produce the matter bounce scenario of holonomy corrected Loop Quantum Cosmology. We focus our study in two limits of the cosmic time, the large cosmic time limit and the small cosmic time limit. For the former, we found that, in the presence of non-interacting and non-relativistic matter, the $F(R)$ gravity that reproduces the late time limit of the matter bounce solution is actually the Einstein-Hilbert gravity plus a power law term. In the early time limit, since it corresponds to large spacetime curvatures, assuming that the Jordan frame is described by a general metric that when it is conformally transformed to the Einstein frame, produces an accelerating Friedmann-Robertson-Walker metric, we found explicitly the scalar field dependence on time. After demonstrating that the solution in the Einstein frame is indeed accelerating, we calculate the spectral index derived from the Einstein frame scalar-tensor counterpart theory of the $F(R)$ theory and compare it with the Planck experiment data. In order to implement the resulting picture, we embed the $F(R)$ gravity explicitly in a Loop Quantum Cosmology framework by introducing holonomy corrections to the $F(R)$ gravity. In this way, the resulting inflation picture corresponding to the $F(R)$ gravity can be corrected in order it coincides to some extent with the current experimental data.

We study some aspects of cosmological evolution in a universe described by a viable curvature corrected exponential $F(R)$ gravity model, in the presence of matter fluids consisting of collisional matter and radiation. Particularly, we express the Friedmann-Robertson-Walker equations of motion in terms of parameters that are appropriate for describing the dark energy oscillations and compare the dark energy density and the dark energy equation of state parameter corresponding to collisional and non-collisional matter. In addition to these, and owing to the fact that the cosmological evolution of collisional and non-collisional matter universes, when quantified in terms of the Hubble parameter and the effective equation of states parameters, is very much alike, we further scrutinize the cosmological evolution study by extending the analysis to the study of matter perturbations in the matter domination era. We quantify this analysis in terms of the growth factor of matter perturbations, in which case the resulting picture of the cosmological evolution is clear, since collisional and non-collisional universes can be clearly distinguished. Interestingly enough, since it is known that the oscillations of the effective equation of state parameter around the phantom divide are undesirable and unwanted in $F(R)$ gravities, when these are considered for redshifts near the matter domination era and before, in the curvature corrected exponential model with collisional matter which we study here there exist oscillations that never cross the phantom divide. Therefore, this rather unwanted feature of the effective equation of state parameter is also absent in the collisional matter filled universe.

We study the late time cosmological evolution of $f(R)$ theories of modified gravity, with the matter content of the universe being that of collisional self interacting matter. We assume that the universe is described by a flat Friedmann-Lemaitre-Robertson-Walker metric and that it is matter and dark energy dominated. The results of our numerical analysis for a collisional matter $f(R)$ theory are compared with those resulting from pressure-less matter $f(R)$ theory and from the $\mathrm{\Lambda}\mathrm{CDM}$ model. As we shall demonstrate, the resulting picture can vary from model to model, indicating that the effect of collisional matter in $f(R)$ theories is strongly model dependent. Particularly, in a few cases, may give better fit to the $\mathrm{\Lambda}\mathrm{CDM}$ model. In all studied cases, the effective equation of state parameter does not cross the phantom divide, both in the collisional matter and pressure-less matter $f(R)$ theories. Finally, we thoroughly study the effects of collisional matter on one of the $f(R)$ models that is known to provide a unified description of early time inflation and late time acceleration. The overall picture of the evolution of the universe is not drastically affected, apart from the matter era which is further enhanced with an additional matter energy density term, which is of leading order. However, a fully consistent description of the universe's evolution requires the introduction of a dark energy compensate in the total energy density, a concept very well known from the literature.

The event rates for the direct detection of dark matter for various types of WIMPs are presented. In addition to the neutralino of SUSY models, we considered other candidates (exotic scalars as well as particles in Kaluza-Klein and technicolour theories) with masses in the TeV region. Then one finds reasonable branching ratios to excited states. Thus the detection of the WIMP can be made not only by recoil measurements, but by measuring the de-excitation $\gamma$-rays as well.