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.

The knowledge of the universe future is of fundamental importance for any advanced civilization. We study the future of singular dark universe where thermal effects due to the Hawking radiation on the apparent horizon of the FRW universe are taken in consideration. It is shown that dark universe which ends up at finite-time Type I and Type III singularity or infinite-time Little Rip singularity transits to finite-time Type II singularity due to account of thermal effects. On the same time, the Type II and IV singular universe does not change its qualitative behavior. The combined account of quantum and thermal effects shows that depending on specific features of the universe only one of effects is dominant. When (conformal matter) quantum effects are dominant, the future singularity is usually removed while for dominant thermal effects the universe final state is the Type II singularity.

In this paper, two inhomogeneous single fluid models for the Universe, which are able to naturally solve the $H_{0}$ tension problem, are discussed. The analysis is based on a Bayesian Machine Learning approach that uses a generative process. The adopted method allows to constrain the free parameters of each model by using the model itself, only. The observable is taken to be the Hubble parameter, obtained from the generative process. Using the full advantages of our method, the models are constrained for two redshift ranges. Namely, first this is done with mock $H(z)$ data over $z\in [0,2.5]$, thus covering known $H(z)$ observational data, which are most helpful to validate the fit results. Then, aiming to extend to redshift ranges to be covered by the most recent ongoing and future planned missions, the models are constrained for the range $z\in[0,5]$, too. Full validation of the results for this extended redshift range will have to wait for the near future, when higher redshift $H(z)$ data become available. This makes our models fully falsifiable. Finally, our second model here is able to explain the BOSS reported value for $H(z)$ at $z=2.34$.

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 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 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.

The present paper represents an attempt for a very generic string inspired theory of gravitation, based on a stringy action in the teleparallel gravity which includes a specific functional which depends on the scalar field and its kinetic energy, as well as the torsion and boundary terms, embedding also possible effects from the teleparallel Gauss--Bonnet invariants. After we deduce the field equations for the associated generic theory of gravitation, we focus on string inspired couplings which are studied by considering different analytical techniques. The first analytical technique is based on the linear stability theory, by introducing proper dimensionless variables which enables us to study the structure of the phase space and the associated physical effects. In this case we have obtained different cosmological solutions which corresponds to matter and dark energy dominated solutions, achieving a possible transition between matter and dark energy dominated epochs. For each type of cosmological solutions we have discussed the corresponding physical features, attaining viable constraints for the coupling constants due to dynamical effects. The dynamical study of the physical features included also a numerical analysis by fine--tuning the initial conditions deep into the matter era, obtaining possible trajectories for the effective equation of state for specific coupling functions.

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 investigated realistic neutron stars in axion $R^{2}$ gravity. The coupling between curvature and axion field $\phi$ is assumed in the simple form $\sim R^2\phi$. For the axion mass in the range $m_{a}\sim 10^{-11}-10^{-10}$ eV the solitonic core within neutron star and corresponding halo with size $\sim 100$ km can exist. Therefore the effective contribution of $R^2$ term grows inside the star and it leads to change of star parameters (namely, mass and radius). We obtained the increase of star mass independent from central density for wide range of masses. Therefore, maximal possible mass for given equation of state grows. At the same time, the star radius increases not so considerably in comparison with GR. {Hence, our model may predict possible existence of supermassive compact stars with masses $M\sim 2.2-2.3M_\odot$ and radii $R_{s}\sim 11$ km for realistic equation of state (we considered APR equation of state). In General Relativity one can obtain neutron stars with such characteristics only for unrealistic, extremely stiff equations of state.} Note that this increase of mass occurs due to change of solution for scalar curvature outside the star. In GR curvature drops to zero on star surface where $\rho=p=0$. In the model under consideration the scalar curvature dumps more slowly in comparison with vacuum $R^2$ gravity due to axion "galo" around the star.

Some cosmological scenarios with bulk viscosity for the dark energy fluid are considered. Based on some considerations related to hydrodynamics, two different equations of state for dark energy are assumed, leading to power-law and logarithmic effective corrections to the pressure. The models are tested with the latest astronomical data from Type Ia supernovae (Pantheon sample), measurements of the Hubble parameter $H(z)$, Baryon Acoustic Oscillations and Cosmic Microwave Background radiation. In comparison with $\Lambda$CDM model, some different results are obtained and their viability is discussed. The power-law model shows some modest results, achieved under negative values of bulk viscosity, while the logarithmic scenario provide good fits in comparison to $\Lambda$CDM model.

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.

We show that there is a correspondence between cosmology from non-extensive thermodynamics and cosmology with fluids of redefined and generalized equation of state. We first establish the correspondence in the case of basic non-extensive thermodynamics, and then we proceed by investigating the more consistent case, from the quantum field theoretical point of view, of varying exponent, namely depending on the scale. The obtained duality provides a way of explaining the complicated phenomenological forms of the effective fluid equation-of-state parameters that are being broadly used in the literature, since their microphysical origin may indeed lie in the non-extensive thermodynamics of spacetime. Finally, concerning the cosmological behavior, we show that at late times the effective fluid may drive the universe acceleration even in the absence of an explicit cosmological constant, and even if the initial fluid is the standard dust matter one. Similarly, at early times we obtain an effective cosmological constant which is enhanced through screening, and hence it can drive a a successful inflation without spoiling the correct late-time acceleration.

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.

We investigate the bounce realization arising from the application of the holographic principle in the early universe, inspired by its well-studied late-time application. We first consider as Infrared cutoffs the particle and future event horizons, and we show that the decrease of the horizons at early times naturally increases holographic energy density at bouncing scales, while we additionally obtain the necessary null energy condition violation. Furthermore, adding a simple correction to the horizons due to the Ultraviolet cutoff we analytically obtain improved nonsingular bouncing solutions, in which the value of the minimum scale factor is controlled by the UV correction. Finally, we construct generalized scenarios, arisen from the use of extended Infrared cutoffs, and as specific examples we consider cutoffs that can reproduce $F(R)$ gravity, and the bounce realization within it.

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 $.

We apply the holographic principle at the early universe, obtaining an inflation realization of holographic origin. Such a consideration has equal footing with its well-studied late-time application, and moreover the decrease of the horizons at early times naturally increases holographic energy density at inflationary scales. Taking as Infrared cutoff the particle or future event horizons, and adding a simple correction due to the Ultraviolet cutoff, whose role is non-negligible at the high energy scales of inflation, we result in a holographic inflation scenario that is very efficient in incorporating inflationary requirements and predictions. We first extract analytically the solution of the Hubble function in an implicit form, which gives a scale factor evolution of the desired e-foldings. Furthermore, we analytically calculate the Hubble slow-roll parameters and then the inflation-related observables, such as the scalar spectral index and its running, the tensor-to-scalar ratio, and the tensor spectral index. Confronting the predictions with Planck 2018 observations we show that the agreement is perfect and in particular deep inside the 1$\sigma$ region.

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}$.

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.

We present a modified cosmological scenario that arises from the application of non-extensive thermodynamics with varying exponent. We extract the modified Friedmann equations, which contain new terms quantified by the non-extensive exponent, possessing standard $\Lambda$CDM cosmology as a subcase. Concerning the universe evolution at late times we obtain an effective dark energy sector, and we show that we can acquire the usual thermal history, with the successive sequence of matter and dark-energy epochs, with the effective dark-energy equation-of-state parameter being in the quintessence or in the phantom regime. The interesting feature of the scenario is that the above behaviors can be obtained even if the explicit cosmological constant is set to zero, namely they arise purely from the extra terms. Additionally, we confront the model with Supernovae type Ia and Hubble parameter observational data, and we show that the agreement is very good. Concerning the early-time universe we obtain inflationary de Sitter solutions, which are driven by an effective cosmological constant that includes the new terms of non-extensive thermodynamics. This effective screening can provide a description of both inflation and late-time acceleration with the same parameter choices, which is a significant advantage.

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.

The intriguing question, why the present scale of the universe is free from any perceptible footprints of rank-2 antisymmetric tensor fields? (generally known as Kalb-Ramond fields) is addressed. A quite natural explanation of this issue is given from the angle of higher-curvature gravity, both in four- and in five-dimensional spacetime. The results here obtained reveal that the amplitude of the Kalb-Ramond field may be actually large and play a significant role during the early universe, while the presence of higher-order gravity suppresses this field during the cosmological evolution, so that it eventually becomes negligible in the current universe. Besides the suppression of the Kalb-Ramond field, the extra degree of freedom in $F(R)$ gravity, usually known as scalaron, also turns out to be responsible for inflation. Such F(R) gravity with Kalb-Ramond fields may govern the early universe to undergo an inflationary stage at early times (with the subsequent graceful exit) for wider range of F(R) gravity than without antisymmetric fields.. Furthermore, the models---in four- and five-dimensional spacetimes---are linked to observational constraints, with the conclusion that the corresponding values of the spectral index and tensor-to-scalar ratio closely match the values provided by the Planck survey 2018 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.

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.

This paper is devoted to the analysis of a class of $F(R)$ gravity, where additional logarithmic corrections are assumed. The gravitational action includes an exponential term and a $R^2$ inflationary term, both with logarithmic corrections. This model can unify an early time inflationary era and also the late time acceleration of the universe expansion. This model is deeply analysed, confronting with recent observational data coming from the largest Pantheon Type Ia supernovae sample, the latest measurements of the Hubble parameter $H(z)$, manifestations of Baryon Acoustic Oscillations and Cosmic Microwave Background radiation. The viability of the model is studied and the corresponding constraints on the free parameters are obtained, leading to an statistical analysis in comparison to $\Lambda$CDM model. The inflationary era is also analysed within this model and its compatibility with the latest observational data for the spectral index of primordial curvature perturbations and the scalar-to-tensor ratio. Finally, possible corrections on the Newton's law and constraints due to primordial nucleosynthesis are analysed.

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.

The cosmological dynamics of a non-locally corrected gravity theory, involving a power of the inverse d'Alembertian, is investigated. Casting the dynamical equations into local form, the fixed points of the models are derived, as well as corresponding de Sitter and power-law solutions. Necessary and sufficient conditions on the model parameters for the existence of de Sitter solutions are obtained. The possible existence of power-law solutions is investigated, and it is proven that models with de Sitter solutions have no power-law solutions. A model is found, which allows to describe the matter-dominated phase of the Universe evolution.

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.

We revisit the concept of turnaround radius in cosmology, in the context of modified gravity. While preliminary analyses were limited to scalar-tensor/$F(R)$ gravity, we extend the definition and the study of this quantity to a much broader class of theories including also quantum $R^2$ gravity. The turnaround radius is computed in terms of the parameters of the theory and it is shown that a deviation not larger than 10\% of this quantity from its value in Einstein's theory could constrain the model parameters and even rule out some current theories.

The propagation of gravitational waves is explored in the cosmological context. It is explicitly demonstrated that the propagation of gravitational waves could be influenced by the medium. It is shown that in the thermal radiation, the propagation of gravitational waves in general relativity is different from that in the scalar-tensor theory. The propagation of gravitational waves is investigated in the uniform magnetic field. As a result, it is found that cosmic magnetic fields could influence on the propagation of gravitational waves to non-negligible extent. The corresponding estimation for the spiral galaxy NGC 6946 effect is made.

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.

Energy conditions can play an important role in defining the cosmological evolution. Specifically acceleration/deceleration of cosmic fluid, as well as the emergence of Big Rip singularities, can be related to the constraints imposed by the energy conditions. Here we discuss this issue for $f(R)$ gravity considering also conformal transformations. Cosmological solutions and equations of state can be classified according to energy conditions. The qualitative change of some energy conditions for transformations from the Jordan frame to the Einstein frame is also observed.

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.

We consider the evolution of cosmological gravitational waves in scalar-tensor theory and $F(R)$ gravity theory as typical models of the modified gravity. Although the propagation speed is not changed from the speed of light, the propagation phase changes when we compare the propagation in these modified gravity theories with the propagation in the $\Lambda$CDM model. The phase change might be detected in future observations.

Here we analysed a particular type of $F(R)$ gravity, the so-called exponential gravity which includes an exponential function of the Ricci scalar in the action. Such term represents a correction to the usual Hilbert-Einstein action. By using Supernovae Ia, Barionic Acoustic Oscillations, Cosmic Microwave Background and $H(z)$ data, the free parameters of the model are well constrained. The results show that such corrections to General Relativity become important at cosmological scales and at late-times, providing an alternative to the dark energy problem. In addition, the fits do not determine any significant difference statistically with respect to the $\Lambda$CDM model. Finally, such model is extended to include the inflationary epoch in the same gravitational Lagrangian. As shown in the paper, the additional terms can reproduce the inflationary epoch and satisfy the constraints from Planck data.

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.

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 investigated the regular multi-horizon black holes in the Einstein gravity, $F(R)$ gravity and the 5 dimensional Gauss-Bonnet gravity, all of them coupled with non-linear electrodynamics. We presented several explicit examples of the actions which admit the solutions describing regular black hole space-time with multi-horizons. Thermodynamics of the obtained black hole solutions is studied. The explicit expressions of the temperature, the entropy, the thermodynamical energy and the free energy are obtained. Although the temperature vanishes in the extremal limit where the radii of the two horizons coincide with each other as in the standard multi-horizon black hole like the Reissner-Nordstr\" om black hole or the Kerr black hole, the larger temperature corresponds to the larger horizon radius. This is different from the standard black holes, where the larger temperature corresponds to the smaller horizon radius. We also found that the specific heat becomes positive for the large temperature, which is also different from the standard black holes, where the specific heat is negative. It should be also noted that the thermodynamical energy is not identical with the ADM mass. Furthermore in case of the Gauss-Bonnet gravity, it is demonstrated that the entropy can become negative.

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 the context of f(R)=R + alpha R^2 gravity, we study the existence of neutron and quark stars with no intermediate approximations in the generalised system of Tolman-Oppenheimer-Volkov equations. Analysis shows that for positive alpha's the scalar curvature does not drop to zero at the star surface (as in General Relativity) but exponentially decreases with distance. Also the stellar mass bounded by star surface decreases when the value alpha increases. Nonetheless distant observers would observe a gravitational mass due to appearance of a so-called gravitational sphere around the star. The non-zero curvature contribution to the gravitational mass eventually is shown to compensate the stellar mass decrease for growing alpha's. We perform our analysis for several equations of state including purely hadronic configurations as well as hyperons and quark stars. In all cases, we assess that the relation between the parameter $\alpha$ and the gravitational mass weakly depend upon the chosen equation of state. Another interesting feature is the increase of the star radius in comparison to General Relativity for stars with masses close to maximal, whereas for intermediate masses around 1.4-1.6 solar masses, the radius of star depends upon alpha very weakly. Also the decrease in the mass bounded by star surface may cause the surface redshift to decrease in R^2-gravity when compared to Einsteinian predictions. This effect is shown to hardly depend upon the observed gravitational mass. Finally, for negative values of alpha our analysis shows that outside the star the scalar curvature has damped oscillations but the contribution of the gravitational sphere into the gravitational mass increases indefinitely with radial distance putting into question the very existence of such relativistic stars.

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 proposed the generalized holographic dark energy model where infrared cutoff is identified with the combination of the FRW universe parameters: the Hubble rate, particle and future horizons, cosmological constant, the universe life-time (if finite) and their derivatives. It is demonstrated that with the corresponding choice of the cutoff one can map such holographic dark energy to modified gravity or gravity with general fluid. Explicitly, F(R) gravity and general perfect fluid are worked out in detail and corresponding infrared cutoff is found. Using this correspondence, we get realistic inflation or viable dark energy or unified inflationary-dark energy universe in terms of covariant holographic dark energy.

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.

We have investigated the attractor structure for the CMB fluctuations in composite inflation scenario within the gauged Nambu-Jona-Lasinio (NJL) model. Such composite inflation represents an attractor which can not be found in a fundamental scalar model. As is known, the number of inflationary models contains the attractor classified by the $\alpha$-attractor model. It is found that the attractor inflation in the gauged NJL model corresponds to the $\alpha = 2$ case.

We investigate the inflationary realization in the context of unimodular $F(T)$ gravity, which is based on the $F(T)$ modification of teleparallel gravity, in which one imposes the unimodular condition through the use of Lagrange multipliers. We develop the general reconstruction procedure of the $F(T)$ form that can give rise to a given scale-factor evolution, and then we apply it in the inflationary regime. We extract the Hubble slow-roll parameters that allow us to calculate various inflation-related observables, such as the scalar spectral index and its running, the tensor-to-scalar ratio, and the tensor spectral index. Then, we examine the particular cases of de Sitter and power-law inflation, of Starobinsky inflation, as well as inflation in a specific model of unimodular $F(T)$ gravity. As we show, in all cases the predictions of our scenarios are in a very good agreement with Planck observational data. Finally, inflation in unimodular $F(T)$ gravity has the additional advantage that it always allows for a graceful exit for specific regions of the model parameters.

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.

We discuss the production and evolution of cosmological gravitons showing how the cosmological background affects their dynamics. Besides, the detection of cosmological gravitons could constitute an extremely important signature to discriminate among different cosmological models. Here we consider the cases of scalar-tensor gravity and $f(R)$ gravity where it is demonstrated the amplification of graviton amplitude changes if compared with General Relativity. Possible observational constraints are discussed.

The possibility to construct an inflationary universe scenario for the finite-scale gauged Nambu-Jona-Lasinio model is investigated. This model can be described by the Higgs-Yukawa type interaction model with the corresponding compositeness scale. Therefore, the one-loop Higgs-Yukawa effective potential is used with the compositeness condition for the study of inflationary dynamics. We evaluate the fluctuations in the cosmic microwave background for the model with a finite compositeness scale in the slow-roll approximation. We find the remarkable dependence on the gauge group and the number of fermion flavors. It is also proved that the model has similar behavior with the $\phi^{4n}$ chaotic inflation and the Starobinsky model at the flat and steep limits, respectively. It is demonstrated that realistic inflation consistent with Planck data is possible for a range of theory parameters.

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.

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 study the power spectra of f(R) inflation using a new technique in which the norm-squared of the mode functions is evolved. Our technique results in excellent analytic approximations for how the spectra depend upon the function $f(R)$. Although the spectra are numerically the same in the Jordan and Einstein frames for the same wave number $k$, they depend upon the geometries of these frames in quite different ways. For example, the power spectra in the two frames are different functions of the number of e-foldings until end of inflation. We discuss how future data on reheating can be used to distinguish f(R) inflation from scalar-driven inflation.

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 discuss the possibility that a Born-Infeld condensate coupled to neutrinos can generate both neutrino masses and an effective cosmological constant. In particular, an effective field theory is provided capable of dynamically realizing the neutrino superfluid phase firstly suggested by Ginzburg and Zharkov. In such a case, neutrinos acquire a mass gap inside the Born-Infeld ether forming a long-range Cooper pair. Phenomenological implications of the approach are also discussed.

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.

We study relativistic stars in the simplest model of the de Rham-Gabadadze-Tolley massive gravity which describes the massive graviton without ghost propagating mode. We consider the hydrostatic equilibrium, and obtain the modified Tolman-Oppenheimer-Volkoff equation and the constraint equation coming from the potential terms in the gravitational action. We give analytical and numerical results for quark and neutron stars and discuss the deviations compared with General Relativity and $F(R)$ gravity. It is shown that theory under investigation leads to small deviation from the General Relativity in terms of density profiles and mass-radius relation. Nevertheless, such deviation may be observable in future astrophysical probes.

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.

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.

The possibility to construct an inflationary scenario for renormalization-group improved potentials corresponding to the Higgs sector of finite gauge models is investigated. Taking into account quantum corrections to the renormalization-group potential which sums all leading logs of perturbation theory is essential for a successful realization of the inflationary scenario, with very reasonable parameter values. The inflationary models thus obtained are seen to be in good agreement with the most recent and accurate observational data. More specifically, the values of the relevant inflationary parameters, $n_s$ and $r$, are close to the corresponding ones in the $R^2$ and Higgs-driven inflation scenarios. It is shown that the model here constructed and Higgs-driven inflation belong to the same class of cosmological attractors.

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 explore a fluid description of the inflationary universe. In particular, we investigate a fluid model in which the equation of state (EoS) for a fluid includes bulk viscosity. We find that the three observables of inflationary cosmology: the spectral index of the curvature perturbations, the tensor-to-scalar ratio of the density perturbations, and the running of the spectral index, can be consistent with the recent Planck results. We also reconstruct the explicit EoS for a fluid from the spectral index of the curvature perturbations compatible with the Planck analysis. In the reconstructed models of a fluid, the tensor-to-scalar ratio of the density perturbations can satisfy the constraints obtained from the Planck satellite. The running of the spectral index can explain the Planck data. In addition, it is demonstrated that in the reconstructed models of a fluid, the graceful exit from inflation can be realized. Moreover, we show that the singular inflation can occur in a fluid model. Furthermore, we show that a fluid description of inflation can be equivalent to the description of inflation in terms of scalar field theories.

We discuss the Mass -Radius diagram for static neutron star models obtained by the numerical solution of modified Tolman-Oppenheimer-Volkoff equations in $f(R)$ gravity where the Lagrangians $f(R)=R+\alpha R^2 (1+\gamma R)$ and $f(R)=R^{1+\epsilon }$ are adopted. Unlike the case of the perturbative approach previously reported, the solutions are constrained by the presence of an extra degree of freedom, coming from the trace of the field equations. In particular, the stiffness of the equation of state determines an upper limit on the central density $\rho_c$ above which the the positivity condition of energy-matter tensor trace $T^{\rm m}=\rho - 3 p$ holds. In the case of quadratic f(R)-gravity, we find higher masses and radii at lower central densities with an inversion of the behavior around a pivoting $\rho_c$ which depends on the choice of the equation of state. When considering the cubic corrections, we find solutions converging to the required asymptotic behavior of flat metric only for $\gamma < 0$. A similar analysis is performed for $f(R)=R^{1+\epsilon }$ considering $\epsilon$ as the leading parameter. We work strictly in the Jordan frame in order to consider matter minimally coupled with respect to geometry. This fact allows us to avoid ambiguities that could emerge in adopting the Einstein frame.

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 investigate the gauged Nambu-Jona-Lasinio model in curved spacetime at the large $N_c$ limit and in slow-roll approximation. The model can be described by the renormalization group corrected gauge-Higgs-Yukawa theory with the corresponding compositeness conditions. Evaluating the renormalization group (RG) improved effective action, we show that such model can produce CMB fluctuations and find inflationary parameters: spectral index, tensor-to-scalar-ratio and running of the spectral index. We demonstrate that the model can naturally satisfy the Planck 2015 data and maybe considered as an alternative candidate for Higgs inflation.

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.

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 explore inflationary cosmology in a theory where there are two scalar fields which non-minimally couple to the Ricci scalar and an additional $R^2$ term, which breaks the conformal invariance. Particularly, we investigate the slow-roll inflation in the case of one dynamical scalar field and that of two dynamical scalar fields. It is explicitly demonstrated that the spectral index of scalar mode of the density perturbations and the tensor-to-scalar ratio can be consistent with the observations acquired by the recent Planck satellite. The graceful exit from the inflationary stage is achieved as in convenient $R^2$ gravity. We also propose the generalization of the model under discussion with three scalar fields.

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 investigate the chaotic inflationary model using the two-loop effective potential of a self-interacting scalar field theory in curved spacetime. We use the potential which contains a non-minimal scalar curvature coupling and a quartic scalar self-interaction. We analyze the Lyapunov stability of de Sitter solution and show the stability bound. Calculating the inflationary parameters, we systematically explore the spectral index $n_s$ and the tensor-to-scalar ratio $r$, with varying the four parameters, the scalar-curvature coupling $\xi_0$, the scalar quartic coupling $\lambda_0$, the renormalization scale $\mu$ and the e-folding number $N$. It is found that the two-loop correction on $n_s$ is much larger than the leading-log correction, which has previously been studied. We show that the model is consistent with the observation by Planck with WMAP and a recent joint analysis of BICEP2.

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.

We consider higher-derivative quantum gravity where renormalization group improved effective action beyond one-loop approximation is derived. Using this effective action, the quantum-corrected FRW equations are analyzed. De Sitter universe solution is found. It is demonstrated that such de Sitter inflationary universe is instable. The slow-roll inflationary parameters are calculated. The contribution of renormalization group improved Gauss-Bonnet term to quantum-corrected FRW equations as well as to instability of de Sitter universe is estimated. It is demonstrated that in this case the spectral index and tensor-to-scalar ratio are consistent with Planck data.

We explore cosmology with a bounce in Gauss-Bonnet gravity where the Gauss-Bonnet invariant couples to a dynamical scalar field. In particular, the potential and and Gauss-Bonnet coupling function of the scalar field are reconstructed so that the cosmological bounce can be realized in the case that the scale factor has hyperbolic and exponential forms. Furthermore, we examine the relation between the bounce in the string (Jordan) and Einstein frames by using the conformal transformation between these conformal frames. It is shown that in general, the property of the bounce point in the string frame changes after the frame is moved to the Einstein frame. Moreover, it is found that at the point in the Einstein frame corresponding to the point of the cosmological bounce in the string frame, the second derivative of the scale factor has an extreme value. In addition, it is demonstrated that at the time of the cosmological bounce in the Einstein frame, there is the Gauss-Bonnet coupling function of the scalar field, although it does not exist in the string frame.

We investigate the cosmological evolution of mimetic matter model with arbitrary scalar potential. The cosmological reconstruction, which is the way to construct a model for arbitrary evolutions of the scale factor, is explicitly done for different choices of potential. The cases that mimetic matter model shows the evolution as Cold Dark Matter(CDM), wCDM model, dark matter and dark energy with dynamical $Om(z)$, where $Om(z) \equiv [ (H(z)/H_0)^2-1 ]/[(1+z)^3 -1],$ or phantom dark energy with phantom-non-phantom crossing are presented in detail. The cosmological perturbations for such evolution are studied in mimetic matter model. For instance, the evolution behavior of the matter density contrast which is different from usual one, i.e. $\ddot \delta + 2 H \dot \delta - \kappa ^2 \rho \delta /2 = 0$ is investigated. The possibility of peculiar evolution of $\delta$ in the model under consideration is shown. Special attention is paid to the behavior of matter density contrast near to future singularity where decay of perturbations may occur much earlier the singularity.

We review inflationary cosmology in modified gravity such as $R^2$ gravity with its extensions in order to generalize the Starobinsky inflation model. In particular, we explore inflation realized by three kinds of effects: modification of gravity, the quantum anomaly, and the $R^2$ term in loop quantum cosmology. It is explicitly demonstrated that in these inflationary models, the spectral index of scalar modes of the density perturbations and the tensor-to-scalar ratio can be consistent with the Planck results. Bounce cosmology in $F(R)$ gravity is also explained.

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.

We investigate inflation in frames of two classes of $F(R)$ gravity and check its consistency with Planck data. It is shown that $F(R)$ inflation without self-reproduction may be constructed in close analogy with the corresponding scalar example proposed by Mukhanov for the resolution the problems of multiverse, predictability and initial conditions.

Quark star models with realistic equation of state in nonperturbative $f(R)$ gravity are considered. The mass-radius relation for $f(R)=R+\alpha R^2$ model is obtained. Considering scalar curvature $R$ as an independent function, one can find out, for each value of central density, the unique value of central curvature for which one has solutions with the required asymptotic $R\rightarrow 0$ for $r\rightarrow\infty$. In another words, one needs a fine-tuning for $R$ to achieve quark stars in $f(R)$ gravity. We consider also the analogue description in corresponding scalar-tensor gravity. The fine-tuning on $R$ is equivalent to the fine-tuning on the scalar field $\phi$ in this description. For distant observers, the gravitational mass of the star increases with increasing $\alpha$ ($\alpha>0$) but the interpretation of this fact depends on frame where we work. Considering directly $f(R)$ gravity, one can say that increasing of mass occurs by the "gravitational sphere" outside the star with some "effective mass". On the other hand, in scalar-tensor theory, we also have a dilaton sphere (or "disphere") outside the star but its contribution to gravitational mass for distant observer is negligible. We show that it is possible to discriminate modified theories of gravity from General Relativity due to the gravitational redshift of the thermal spectrum emerging from the surface of the star.

We propose the mimetic $F(R)$ theory and investigate the early-time and late-time acceleration in such theory. It is demonstrated that inflation consistent with observable data may be realized in such theory. The reconstruction of realistic $\Lambda$CDM era is also possible as well as the unification of early-time inflation with late-time acceleration or bounce universe. It is stressed that specific universe evolution is governed by mimetic $F(R)$ theory which is different from convenient $F(R)$ gravity. The corresponding examples are presented. Mimetic $F(R)$ gravity is generalized by the addition of the scalar potential in the formulation of convenient $F(R)$ gravity with specific Lagrange multiplier constraint. It is demonstrated that such theory may admit the arbitrary universe evolution via the corresponding choice of the scalar potential and/or function $F(R)$.

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 discuss neutron stars with strong magnetic mean fields in the framework of Extended Theories of Gravity. In particular, we take into account models derived from $f(R)$ and $f(\cal G)$ extensions of General Relativity where functions of the Ricci curvature invariant $R$ and the Gauss-Bonnet invariant ${\cal G}$ are respectively considered. Dense matter in magnetic mean field, generated by magnetic properties of particles, is described by assuming a model with three meson fields and baryons octet. As result, the considerable increasing of maximal mass of neutron stars can be achieved by cubic corrections in $f(R)$ gravity. In principle, massive stars with $M> 4 M_{\odot}$ can be obtained. On the other hand, stable stars with high strangeness fraction (with central densities $\rho_{c}\sim 1.5-2.0$ GeV/fm$^{3}$) are possible considering quadratic corrections of $f(\cal {G})$ gravity. The magnetic field strength in the star center is of order $6-8\times 10^{18}$ G. In general, we can say that other branches of massive neutron stars are possible considering the extra pressure contributions coming from gravity extensions. Such a feature can constitute both a probe for alternative theories and a way out to address anomalous self-gravitating compact systems.

We consider two types of modifications of Born-Infeld gravity in the Palatini formulation and explore their dynamics in the early universe. One of these families considers $f(R)$ corrections to the Born-Infeld Lagrangian, which can be seen as modifications of the dynamics produced by the quantum effects of matter, while the other consists on different powers of the elementary building block of the Born-Infeld Lagrangian, which we denote by $|\hat\Omega|$. We find that the two types of nonsingular solutions that arise in the original Born-Infeld theory are also present in these extensions, being bouncing solutions a stable and robust branch. Singular solutions with a period of approximate de Sitter inflation are found even in universes dominated by radiation.

We investigate the descriptions for the observables of inflationary models, in particular, the spectral index of curvature perturbations, the tensor-to-scalar ratio, and the running of the spectral index, in the framework of perfect fluid models and $F(R)$ gravity theories through the reconstruction methods. Furthermore, the perfect fluid and $F(R)$ gravity descriptions of inflation are compared with the recent cosmological observations such as the Planck satellite and BICEP2 experiment. It is demonstrated with explicit examples that perfect fluid may lead to the inflationary universe consistent with the Planck data. It is also shown that several $F(R)$ gravity models, especially, a power-law model gives the best fit values compatible with the spectral index and tensor-to-scalar ratio within the allowed ranges suggested by the Planck and BICEP2 results.

We study the instabilities and related anti-evaporation of the extremal Reissner-Nordstr\"om (RN) black hole in $F(R)$ gravity. It is remarkable that the effective electric charge can be generated for some solutions of $F(R)$ gravity without electromagnetic field. The anti-evaporation effect occurs but it emerges only in the strong coupling limit of the effective gravitational coupling. The instabilities of RN black hole are also investigated when the electromagnetic sector is added to the action of $F(R)$ gravity. We show the anti-evaporation occurs in the Maxwell-$F(R)$ gravity with the arbitrary gravitational coupling constant although it does not occur in the Maxwell-Einstein gravity. Furthermore, general spherically-symmetric solution of $F(R)$ gravity in the Einstein frame is obtained.

We study inflation induced by (power-low) scalar curvature corrections to General Relativity. The class of inflationary scalar potentials $V(\sigma)\sim\exp[n\,\sigma]$, $n$ general parameter, is investigated in the Einsein frame and the corresponding actions in the Jordan frame are derived. We found the conditions for which these potentials are able to reproduce viable inflation according with the last cosmological data and lead to large scalar curvature corrections which emerge only at a mass scale larger than the Planck mass. Cosmological constant may appear or be set equal to zero in the Jordan frame action without changing the behaviour of the model during inflation. Moreover, polynomial corrections to General Relativity are analyzed in detail. When de Sitter space-time emerges as an exact solution of the models, it is necessary to use perturbative equations in the Jordan framework to study their dynamics during the inflation. In this case, we demonstrate that the Ricci scalar decreases after a correct amount of inflation, making the models consistent with the observable evolution of the universe.

The possibility to construct inflationary models for the renormalization-group improved potentials corresponding to scalar electrodynamics and to $SU(2)$ and $SU(5)$ models is investigated. In all cases, the tree-level potential, which corresponds to the cosmological constant in the Einstein frame, is seen to be non-suitable for inflation. Rather than adding the Hilbert-Einstein term to the action, quantum corrections to the potential, coming from to the RG-equation, are included. The inflationary scenario is analyzed with unstable de Sitter solutions which correspond to positive values of the coupling function, only. We show that, for the finite $SU(2)$ model and $SU(2)$ gauge model, there are no de Sitter solutions suitable for inflation, unless exit from it occurs according to some weird, non-standard scenarios. Inflation is realized both for scalar electrodynamics and for $SU(5)$ RG-improved potentials, and the corresponding values of the coupling function are seen to be positive. It is shown that, for quite reasonable values of the parameters, the inflationary models obtained both from scalar electrodynamics and from the $SU(5)$ RG-improved potentials, are in good agreement with the most recent observational data coming from the Planck2013 and BICEP2 collaborations.

We reconstruct scalar field theories to realize inflation compatible with the BICEP2 result as well as the Planck. In particular, we examine the chaotic inflation model, natural (or axion) inflation model, and an inflationary model with a hyperbolic inflaton potential. We perform an explicit approach to find out a scalar field model of inflation in which any observations can be explained in principle.

Inflation models can be examined by the cosmological observations, WMAP, Planck, BICEP2 and so on. These observations directly constrain the spectral index, $n_s$, and the tensor-to-scalar ratio, $r$. Besides, from a theoretical point of view, it has been shown that any inflation models asymptote a universal attractor in $(n_s,r)$ plane for a larger scalar-gravity coupling. In this work we consider a simple chaotic inflation model with a scalar quartic and a scalar-curvature interactions. The quantum corrections are introduced for these interactions through the renormalization group. The inflationary parameters, $n_s$ and $r$, are numerically calculated with shifting the bare scalar-gravity coupling $\xi_0$, the quartic scalar bare coupling $\lambda_0$, the renormalization scale $\mu$ and the e-folding number $N$. The Planck data is consistent with the $\phi^4$ theory with a finite scalar-curvature coupling. It is found that the RG running induces a non-universal contribution for $n_s$ and $r$. It can increase the tensor-to-scalar ratio, $r$, which is consistent with Planck data so that it may approach the BICEP2 data for a small renormalization scale.

We discuss the Noether Symmetry Approach in the framework of Gauss-Bonnet cosmology showing that the functional form of the $F(R, {\cal G})$ function, where $R$ is the Ricci scalar and ${\cal G}$ is the Gauss-Bonnet topological invariant, can be determined by the presence of symmetries. Besides, the method allows to find out exact solutions due to the reduction of cosmological dynamical system and the presence of conserved quantities. Some specific cosmological models are worked out

We explore conformal-anomaly driven inflation in $F(R)$ gravity without invoking the scalar-tensor representation. We derive the stress-energy tensor of the quantum anomaly in the flat homogeneous and isotropic universe. We investigate a suitable toy model of exponential gravity plus the quantum contribution due to the conformal anomaly, which leads to the de Sitter solution. It is shown that in $F(R)$ gravity model, the curvature perturbations with its enough amplitude consistent with the observations are generated during inflation. We also evaluate the number of $e$-folds at the inflationary stage and the spectral index $n_\mathrm{s}$ of scalar modes of the curvature perturbations by analogy with scalar tensor theories, and compare them with the observational data. As a result, it is found that the Ricci scalar decreases during inflation and the standard evolution history of the universe is recovered at the small curvature regime. Furthermore, it is demonstrated that in our model, the tensor-to-scalar ratio of the curvature perturbations can be a finite value within the $68\%\,\mathrm{CL}$ error of the very recent result found by the BICEP2 experiment.

Motivated by the properties of matter quantum fields in curved space-times, we work out a gravity theory that combines the Born-Infeld gravity Lagrangian with an $f(R)$ piece. To avoid ghost-like instabilities, the theory is formulated within the Palatini approach. This construction provides more freedom to address a number of important questions such as the dynamics of the early universe and the cosmic accelerated expansion, among others. In particular, we consider the effect that adding an $f(R)=a R^2$ term has on the early-time cosmology. We find that bouncing solutions are robust against these modifications of the Lagrangian whereas the solutions with {\it loitering} behavior of the original Born-Infeld theory are very sensitive to the $R^2$ term. In fact, these solutions are modified in such a way that a plateau in the $H^2$ function may arise yielding a period of (approximately) de Sitter inflationary expansion. This inflationary behavior may be found even in a radiation dominated universe.

Neutron stars with strong magnetic fields are considered in the framework of f(R) gravity. In order to describe dense matter in magnetic field, the model with baryon octet interacting through $\sigma$$\rho$$\omega$-fields is used. The hyperonization process results in softening the equation of state (EoS) and in decreasing the maximal mass. We investigate the effect of strong magnetic field in models involving quadratic and cubic corrections in the Ricci scalar $R$ to the Hilbert-Einstein action. For large fields, the Mass-Radius relation differs considerably from that of General Relativity only for stars with masses close to the maximal one. Another interesting feature is the possible existence of more compact stable stars with extremely large magnetic fields ($\sim 6\times 10^{18}$ G instead of $\sim 4\times 10^{18}$ G as in General Relativity) in the central regions of the stars. Due to cubic terms, a significant increasing of the maximal mass is possible.

Motivated by issues on inflation, a generalized modified gravity model is investigated, where the model Lagrangian is described by a smooth function $f(R, K, \phi)$ of the Ricci scalar $R$, the kinetic term $K$ of a scalar field $\phi$. In particular, the one-loop effective action in the de Sitter background is examined on-shell as well as off-shell in the Landau gauge. In addition, the on-shell quantum equivalence of $f(R)$ gravity in the Jordan and Einstein frames is explicitly demonstrated. Furthermore, we present applications related to the stability of the de Sitter solutions and the one-loop quantum correction to inflation in quantum-corrected $R^2$ gravity. It is shown that for a certain range of parameters, the spectral index of the curvature perturbations can be consistent with the Planck analysis, but the tensor-to-scalar ratio is smaller than the minimum value within the 1 $\sigma$ error range of the BICEP2 result.

We explore bounce cosmology in $F(\mathcal{G})$ gravity with the Gauss-Bonnet invariant $\mathcal{G}$. We reconstruct $F(\mathcal{G})$ gravity theory to realize the bouncing behavior in the early universe and examine the stability conditions for its cosmological solutions. It is demonstrated that the bouncing behavior with an exponential as well as a power-law scale factor naturally occurs in modified Gauss-Bonnet gravity. We also derive the $F(\mathcal{G})$ gravity model to produce the ekpyrotic scenario. Furthermore, we construct the bounce with the scale factor composed of a sum of two exponential functions and show that not only the early-time bounce but also the late-time cosmic acceleration can occur in the corresponding modified Gauss-Bonnet gravity. Also, the bounce and late-time solutions in this unified model is explicitly analyzed.

We review recent progress on cosmological issues and theoretical properties of modified gravity theories. In particular, we explicitly explore the conformal transformation, the Starobinsky inflation, and a unified scenario of inflation and late time acceleration in $F(R)$ gravity and $F(T)$ gravity (extended teleparallel gravity). Furthermore, we examine neutron stars and the hyperon problem in $F(R)$ gravity. Moreover, for loop quantum cosmology (LQC), the natures of finite-time future singularities in $F(T)$ gravity are presented. In addition, we investigate $F(T)$ gravity theories from the Kaluza-Klein (KK) and Randall-Sundrum (RS) theories.

We explore trace-anomaly driven inflation in modified gravity. It is explicitly shown that in $T^2$ teleparallel gravity, the de Sitter inflation can occur, although quasi de Sitter inflation happens in $R^2$ gravity. Furthermore, we investigate the influence of the trace anomaly on inflation. It is found that in $f(T)$ gravity, the de Sitter inflation can end because it becomes unstable due to the trace anomaly, whereas also in higher derivative gravity, the de Sitter inflation can be realized and it will be over thanks to the trace anomaly for smaller parameter regions in comparison with those in teleparallelism. The instability of the de Sitter inflation in $T^2$ gravity and $R^2$ gravity (both with taking account of the trace anomaly) is examined. In addition, we study trace-anomaly driven inflation in minimal massive bigravity, where the contribution from the massive graviton acts as negative cosmological constant. It is demonstrated that the de Sitter inflation can occur and continue for long enough duration.

We construct an $F(R)$ gravity theory corresponding to the Weyl invariant two scalar field theory. We investigate whether such $F(R)$ gravity can have the antigravity regions where the Weyl curvature invariant does not diverge at the Big Bang and Big Crunch singularities. It is revealed that the divergence cannot be evaded completely but can be much milder than that in the original Weyl invariant two scalar field theory.

We reconstruct $F(R)$ gravity models with exponential and power-law forms of the scale factor in which bounce cosmology can be realized. We explore the stability of the reconstructed models with analyzing the perturbations from the background solutions. Furthermore, we study an $F(R)$ gravity model with a sum of exponentials form of the scale factor, where the bounce in the early universe as well as the late-time cosmic acceleration can be realized in a unified manner. As a result, we build a second order polynomial type model in terms of $R$ and show that it could be stable. Moreover, when the scale factor is expressed by an exponential form, we derive $F(R)$ gravity models of a polynomial type in case of the non-zero spatial curvature and that of a generic type in that of the zero spatial curvature. In addition, for an exponential form of the scale factor, an $F(R)$ bigravity model realizing the bouncing behavior is reconstructed. It is found that in both the physical and reference metrics the bouncing phenomenon can occur, although in general the contraction and expansion rates are different each other.

Neutron star models in perturbative $f(R)$ gravity are considered with realistic equations of state. In particular, we consider the FPS, SLy and other equations of state and a case of piecewise equation of state for stars with quark cores. The mass-radius relations for $f(R)=R+R(e^{-R/R_{0}}-1)$ model and for $R^2$ models with logarithmic and cubic corrections are obtained. In the case of $R^2$ gravity with cubic corrections, we obtain that at high central densities ($\rho>10\rho_{ns}$, where $\rho_{ns}=2.7\times 10^{14}$ g/cm$^{3}$ is the nuclear saturation density), stable star configurations exist. The minimal radius of such stars is close to $9$ km with maximal mass $\sim 1.9 M_{\odot}$ (SLy equation). A similar situation takes place for AP4 and BSK20 EoS. Such an effect can give rise to more compact stars than in General Relativity. If observationally identified, such objects could constitute a formidable signature for modified gravity at astrophysical level. Another interesting result can be achieved in modified gravity with only a cubic correction. For some EoS, the upper limit of neutron star mass increases and therefore these EoS can describe realistic star configurations (although, in General Relativity, these EoS are excluded by observational constraints).

The modified Lane-Emden equation for stellar hydrostatic equilibrium in f(R)-gravity is numerically solved by using an iterative procedure. Such an integro-differential equation can be obtained in the weak field limit approximation of f(R)-gravity by considering a suitable polytropic equation of state in the modified Poisson equation. The approach allows, in principle, to deal with still unexplored self-graviting systems that could account for exotic stellar structures that escape standard stellar theory.

We propose new version of massive $F(R)$ gravity which is natural generalization of convenient massive ghost-free gravity. Its Hamiltonian formulation in scalar-tensor frame is developed. We show that such $F(R)$ theory is ghost-free. The cosmological evolution of such theory is investigated. Despite the strong Bianchi identity constraint the possibility of cosmic acceleration (especially, in the presence of cold dark matter) is established. Ghost-free massive $F(R,T)$ gravity is also proposed.

We investigate the future evolution of Nariai black hole which is extremal limit of Schwarzschild-de Sitter one in modified gravity. The perturbations equations around Nariai black hole are derived in static and cosmological patches for general $F(R)$-gravity. The analytical and numerical study of several realistic $F(R)$-models shows the occurence of rich variety of scenarios: instabilities, celebrated Hawking evaporation and anti-evaporation of black hole. The realization of specific scenario depends on the model under consideration. It is remarkable that the presence of such primordial black holes at current universe may indicate towards the modified gravity which supports the anti-evaporation as preferrable model. As some generalization we extend the study of Nariai black hole evolution to modified Gauss-Bonnet gravity. The corresponding perturbations equations turn out to be much more complicated than in the case of $F(R)$ gravity. For specific example of modified Gauss-Bonnet gravity we demonstrate that Nariai solution maybe stable.

We discuss conformal issues of pure and extended teleparallel gravity. In particular, we present formulations of conformal transformation in teleparallel gravity. Furthermore, we propose conformal scalar and gauge field theories in teleparallel gravity and study conformal torsion gravity. We explicitly demonstrate that a power-law acceleration (including the $\Lambda$CDM universe) as well as the de Sitter expansion of the universe can be realized in extended teleparallel gravity with a conformal scalar field, and that pure teleparallel gravity with a conformal scalar field can lead to the $\Lambda$CDM model. It is also shown that there can exist the de Sitter solution in conformal torsion gravity.

We propose general f(R,T,R_{\mu\nu}T^{\mu\nu}) theory as generalization of covariant Horava-like gravity with dynamical Lorentz symmetry breaking. FRLW cosmological dynamics for several versions of such theory is considered. The reconstruction of the above action is explicitly done, including the numerical reconstruction for the occurrence of \Lambda CDM universe. De Sitter universe solutions in the presence of non-constant fluid are also presented. The problem of matter instability in f(R,T,R_{\mu\nu}T^{\mu\nu}) gravity is discussed.

We explore the four-dimensional effective $F(T)$ gravity with $T$ the torsion scalar in teleparallelism originating from higher-dimensional space-time theories, in particular the Kaluza-Klein (KK) and Randall-Sundrum (RS) theories. First, through the KK dimensional reduction from the five-dimensional space-time, we obtain the four-dimensional effective theory of $F(T)$ gravity with its coupling to a scalar field. Second, taking the RS type-II model in which there exist the five-dimensional Anti-de Sitter (AdS) space-time with four-dimensional Friedmann-Lema\^{i}tre-Robertson-Walker (FLRW) brane, we find that there will appear the contribution of $F(T)$ gravity on the four-dimensional FLRW brane. It is demonstrated that inflation and the dark energy dominated stage can be realized in the KK and RS models, respectively, due to the effect of only the torsion in teleparallelism without that of the curvature.

We consider modified gravity which may describe the early-time inflation and/or late-time cosmic acceleration of the universe. In particular, we discuss the properties of $F(R)$, $F(G)$, string-inspired and scalar-Einstein-Gauss-Bonnet gravities, including their FRW equations and fluid or scalar-tensor description. Simplest accelerating cosmologies are investigated and possibility of unified description of the inflation with dark energy is described. The cosmological reconstruction program which permits to get the requested universe evolution from modified gravity is developed. As some extension, massive $F(R)$ bigravity which is ghost-free theory is presented. Its scalar-tensor form turns out to be the easiest formulation. The cosmological reconstruction method for such bigravity is presented. The unified description of inflation with dark energy in $F(R)$ bigravity turns out to be possible.

We study the accelerating cosmology in massive $F(R)$ bigravity via the reconstruction scheme. The consistent solution of the FRW equations is presented: it includes Big and Little Rip, quintessence, de Sitter and decelerating universes described by the physical $g$ metric while the corresponding solution of the universe described by the reference $f$ metric is also found. It is demonstrated that in general the cosmological singularities of $g$ metric are not always manifested as cosmological singularities of the reference $f$ metric. We give two consistent ways to describe the Big and Little Rip, quintessence, de Sitter and decelerating universes. In one of the consistent solutions, the two metrics $g$ and $f$ coincide with each other, which may indicate the connection with the convenient single metric background formulation. For another solution, where two metrics $g$ and $f$ do not coincide with each other, here always appears super-luminal mode.

We studied the anti-evaporation of degenerate Schwarzschild-de Sitter black hole (so-called Nariai space-time) in modified $F(R)$ gravity. The analysis of perturbations of the Nariai black hole is done with the conclusion that anti-evaporation may occur in such a theory already on classical level. For several power-law $F(R)$ gravities which may describe the inflation and/or dark energy eras we presented the theory parameters bounds for occurrence of anti-evaporation and conjectured creation of infinite number of horizons.

A string-inspired effective theory of gravity, containing Gauss-Bonnet invariant interacting with a scalar field, is considered in view of obtaining cosmological dark energy solutions. A Lagrange multiplier is inserted into the action in order to achieve the cosmological reconstruction by selecting suitable forms of couplings and potentials. Several cosmological exact solutions (including dark energy of quintessence, phantom or Little Rip type) are derived in presence and in absence of the Lagrange multiplier showing the difference in the two dynamical approaches. In the models that we consider, the Lagrange multiplier behaves as a sort of dust fluid that realizes the transitions between matter dominated and dark energy epochs. The relation between Lagrange multipliers and Noether symmetries is discussed.

We review the accelerating (mainly, dark energy) cosmologies in modified gravity. Special attention is paid to cosmologies leading to finite-time future singularities in $F(R)$, $F(G)$ and $\mathcal{F}(R,G)$ modified gravities. The removal of the finite-time future singularities via addition of $R^2$-term which simultaneously unifies the early-time inflation with late-time acceleration is also briefly mentioned. Accelerating cosmology including the scenario unifying inflation with dark energy is considered in $F(R)$ gravity with Lagrange multipliers. In addition, we examine domain wall solutions in $F(R)$ gravity. Furthermore, covariant higher derivative gravity with scalar projectors is explored.

The fermionic condensate and vacuum expectation value of the energy-momentum tensor, for a massive fermionic field on the background of anti-de Sitter spacetime, in the geometry of two parallel boundaries with bag boundary conditions, are investigated. Vacuum expectation values, expressed as series involving the eigenvalues of the radial quantum number, are neatly decomposed into boundary-free, single-boundary-induced, and second-boundary-induced parts, with the help of the generalized Abel-Plana summation formula. In this way, the renormalization procedure is very conveniently reduced to the one corresponding to boundary-free AdS spacetime. The boundary-induced contributions to the fermionic condensate and to the vacuum expectation value of the energy density are proven to be everywhere negative. The vacuum expectation values are exponentially suppressed at distances from the boundaries much larger than the curvature radius of the AdS space. Near the boundaries, effects related with the curvature of the background spacetime are shown to be subdominant and, to leading order, all known results for boundaries in the Minkowski bulk are recovered. Zeta function techniques are successfully used for the evaluation of the total vacuum energy in the region between the boundaries. It is proven that the resulting interaction forces between them are attractive and that, for large separations, they also decay exponentially. Finally, our results are extended and explicitly translated to fermionic Casimir densities in braneworld scenarios of Randall-Sundrum type.

Several dark energy models on the brane are investigated. They are compared with corresponding theories in the frame of 4d Friedmann-Robertson-Walker cosmology. To constrain the parameters of the models considered, recent observational data, including SNIa apparent magnitude measurements, baryon acoustic oscillation results, Hubble parameter evolution data and matter density perturbations are used. Explicit formulas of the so-called {\it state-finder} parameters in teleparallel theories are obtained that could be useful to test these models and to establish a link between Loop Quantum Cosmology and Brane Cosmology. It is concluded that a joint analysis as the one developed here allows to estimate, in a very convenient way, possible deviation of the real universe cosmology from the standard Friedmann-Robertson-Walker one.

A generic feature of viable exponential $F(R)$-gravity is investigated. An additional modification to stabilize the effective dark energy oscillations during matter era is proposed and applied to two viable models. An analysis on the future evolution of the universe is performed. Furthermore, a unified model for early and late-time acceleration is proposed and studied.

We demonstrate how holonomy corrections in loop quantum cosmology (LQC) prevent the Big Rip singularity by introducing a quadratic modification in terms of the energy density $\rho$ in the Friedmann equation in the Friedmann-Lema\^{i}tre-Robertson-Walker (FLRW) space-time in a consistent and useful way. In addition, we investigate whether other kind of singularities like Type II,III and IV singularities survive or are avoided in LQC when the universe is filled by a barotropic fluid with the state equation $P=-\rho-f(\rho)$, where $P$ is the pressure and $f(\rho)$ a function of $\rho$. It is shown that the Little Rip cosmology does not happen in LQC. Nevertheless, the occurrence of the Pseudo-Rip cosmology, in which the phantom universe approaches the de Sitter one asymptotically, is established, and the corresponding example is presented. It is interesting that the disintegration of bound structures in the Pseudo-Rip cosmology in LQC always takes more time than that in Einstein cosmology. Our investigation on future singularities is generalized to that in modified teleparallel gravity, where LQC and Brane Cosmology in the Randall-Sundrum scenario are the best examples. It is remarkable that $F(T)$ gravity may lead to all the kinds of future singularities including Little Rip.

We confront dark energy models which are currently similar to $\Lambda$CDM theory with observational data which include the SNe data, matter density perturbations and baryon acoustic oscillations data. DE cosmology under consideration may evolve to Big Rip, type II or type III future singularity, or to Little Rip or Pseudo-Rip universe. It is shown that matter perturbations data define more precisely the possible deviation from $\Lambda$CDM model than consideration of SNe data only. The combined data analysis proves that DE models under consideration are as consistent as $\Lambda$CDM model. We demonstrate that growth of matter density perturbations may occur at sufficiently small background density but still before the possible disintegration of bound objects (like clusters of galaxies, galaxies, etc) in Big Rip, type III singularity, Little Rip or Pseudo-Rip universe. This new effect may bring the future universe to chaotic state well before disintegration or Rip.

Brane dark energy cosmologies, leading to various possible evolutions of our universe, are investigated. The discussion shows that while all these models can be made arbitrarily close to the standard $\Lambda$CDM cosmology at present, their future evolutions can be very different, even diverge with time in a number of ways. This includes asymptotic de-Sitter evolution, Little Rip with dissolution of bound structures, and various possible singularities, as the Big Rip, a sudden future singularity (Type II), and Type III and Type IV cases. Specifically, some interesting effects coming from the brane tension are investigated. It is shown, in particular, that the Little Rip occurs faster on the brane model than in usual FRW cosmology. And in the asymptotic de-Sitter regime the influence of the brane tension leads to a deviation of the value of the effective cosmological constant from that corresponding to ordinary dark energy. As a consequence, the value of the inertial force from the accelerating expansion can greatly exceed the corresponding inertial force in ordinary cosmological models.

We develop the general reconstruction scheme in two scalar model. The quintom-like theory which may describe (different) non-singular Little Rip or de Sitter cosmology is reconstructed. (In)stability of such dark energy cosmologies as well as the flow to fixed points is studied. The stability of Little Rip universe which leads to dissolution of bound objects sometime in future indicates that no classical transition to de Sitter space occurs.

A generic feature of viable $F(R)$ gravity is investigated: It is demonstrated that during the matter dominated era the large frequency oscillations of the effective dark energy may influence the behavior of higher derivatives of the Hubble parameter with the risk to produce some singular unphysical solutions at high redshift. This behavior is explicitly analyzed for realistic $F(R)$ models, in particular, exponential gravity and a power form model. To stabilize such oscillations, we consider the additional modification of the models via a correction term which does not destroy the viability properties. A detailed analysis on the future evolution of the universe and the evolution history of the growth index of the matter density perturbations are performed. Furthermore, we explore two applications of exponential gravity to the inflationary scenario. We show how it is possible to obtain different numbers of $e$-folds during the early-time acceleration by making different choices of the model parameters in the presence of ultrarelativistic matter, which destabilizes inflation and eventually leads to the exit from the inflationary stage. We execute the numerical analysis of inflation in two viable exponential gravity models. It is proved that at the end of the inflation, the effective energy density and curvature of the universe decrease and thus a unified description between inflation and the $\Lambda$CDM-like dark energy dominated era can be realized.

We propose a bigravity analogue of the $F(R)$ gravity. Our construction is based on recent ghost-free massive bigravity where additional scalar fields are added and the corresponding conformal transformation is implemented. It turns out that $F(R)$ bigravity is easier to formulate in terms of the auxiliary scalars as the explicit presentation in terms of $F(R)$ is quite cumbersome. The consistent cosmological reconstruction scheme of $F(R)$ bigravity is developed in detail, showing the possibility to realize nearly arbitrary physical universe evolution with consistent solution for second metric. The examples of accelerating universe which includes phantom, quintessence and $\Lambda$CDM acceleration are worked out in detail and their physical properties are briefly discussed.

We review different dark energy cosmologies. In particular, we present the $\Lambda$CDM cosmology, Little Rip and Pseudo-Rip universes, the phantom and quintessence cosmologies with Type I, II, III and IV finite-time future singularities and non-singular dark energy universes. In the first part, we explain the $\Lambda$CDM model and well-established observational tests which constrain the current cosmic acceleration. After that, we investigate the dark fluid universe where a fluid has quite general equation of state (EoS) [including inhomogeneous or imperfect EoS]. All the above dark energy cosmologies for different fluids are explicitly realized, and their properties are also explored. It is shown that all the above dark energy universes may mimic the $\Lambda$CDM model currently, consistent with the recent observational data. Furthermore, special attention is paid to the equivalence of different dark energy models. We consider single and multiple scalar field theories, tachyon scalar theory and holographic dark energy as models for current acceleration with the features of quintessence/phantom cosmology, and demonstrate their equivalence to the corresponding fluid descriptions. In the second part, we study another equivalent class of dark energy models which includes $F(R)$ gravity as well as $F(R)$ Ho\v{r}ava-Lifshitz gravity and the teleparallel $f(T)$ gravity. The cosmology of such models representing the $\Lambda$CDM-like universe or the accelerating expansion with the quintessence/phantom nature is described. Finally, we approach the problem of testing dark energy and alternative gravity models to general relativity by cosmography. We show that degeneration among parameters can be removed by accurate data analysis of large data samples and also present the examples.

We discuss the Hamiltonian dynamics for cosmologies coming from Extended Theories of Gravity. In particular, minisuperspace models are taken into account searching for Noether symmetries. The existence of conserved quantities gives selection rule to recover classical behaviors in cosmic evolution according to the so called Hartle criterion, that allows to select correlated regions in the configuration space of dynamical variables. We show that such a statement works for general classes of Extended Theories of Gravity and is conformally preserved. Furthermore, the presence of Noether symmetries allows a straightforward classification of singularities that represent the points where the symmetry is broken. Examples of nonminimally coupled and higher-order models are discussed.

Multiple $\Lambda$CDM cosmology is studied in a way that is formally a classical analog of the Casimir effect. Such cosmology corresponds to a time-dependent dark fluid model or, alternatively, to its scalar field presentation, and it motivated by the string landscape picture. The future evolution of the several dark energy models constructed within the scheme is carefully investigated. It turns out to be almost always possible to choose the parameters in the models so that they match the most recent and accurate astronomical values. To this end, several universes are presented which mimick (multiple) $\Lambda$CDM cosmology but exhibit Little Rip, asymptotically de Sitter, or Type I, II, III, and IV finite-time singularity behavior in the far future, with disintegration of all bound objects in the cases of Big Rip, Little Rip and Pseudo-Rip cosmologies.

Dark energy models with various scenarios of evolution are considered from the viewpoint of the formalism for the equation of state. It is shown that these models are compatible with current astronomical data. Some of the models presented here evolve arbitrarily close to $\Lambda$CDM up to the present, but diverge in the future into a number of different possible asymptotic states, including asymptotic de-Sitter (pseudo-rip) evolution, little rips with disintegration of bound structures, and various forms of finite-time future singularities. Therefore it is impossible from observational data to determine whether the universe will end in a future singularity or not. We demonstrate that the models under consideration are stable for a long period of time (billions of years) before entering a Little Rip/Pseudo-Rip induced dissolution of bound structures or before entering a soft finite-time future singularity. Finally, the physical consequences of Little Rip, Type II, III and Big Crush singularities are briefly compared.

We propose the most general modified first-order Ho\v{r}ava-Lifshitz (HL) gravity, whose action does not contain time derivatives higher than the second order. The Hamiltonian structure of this theory is studied in all the details in the case of the spatially-flat FRW space-time, demonstrating many of the features of the general theory. It is shown that, with some plausible assumptions, including the projectability of the lapse function, this model is consistent. As a large class of such theories, the modified HL $F(R)$ gravity is introduced. The study of its ultraviolet properties shows that its $z=3$ version seems to be renormalizable in the same way as the original HL proposal. The Hamiltonian analysis of the modified HL $F(R)$ gravity shows that it is in general a consistent theory. The $F(R)$ gravity action is also studied in the fixed-gauge form, where the appearance of a scalar field is particularly illustrative. Then the spatially-flat FRW cosmology for this $F(R)$ gravity is investigated. It is shown that a special choice of parameters for this theory leads to the same equations of motion as in the case of traditional $F(R)$ gravity. Nevertheless, the cosmological structure of the modified HL $F(R)$ gravity turns out to be much richer than for its traditional counterpart. The emergence of multiple de Sitter solutions indicates to the possibility of unification of early-time inflation with late-time acceleration within the same model. Power-law $F(R)$ theories are investigated in detail. It is analytically shown that they have a quite rich cosmological structure: early/late-time cosmic acceleration of quintessence, as well as of phantom types. Also it is demonstrated that all the four known types of finite-time future singularities may occur in the power-law HL $F(R)$ gravity. Finally, a covariant proposal for (renormalizable) $F(R)$ gravity within the HL spirit is presented.

We propose a general approach for the construction of modified gravity which is invariant under foliation-preserving diffeomorphisms. Special attention is paid to the formulation of modified $F(R)$ Ho\v{r}ava-Lifshitz gravity (FRHL), whose Hamiltonian structure is studied. It is demonstrated that the spatially-flat FRW equations of FRHL are consistent with the constraint equations. The analysis of de Sitter solutions for several versions of FRHL indicates that the unification of the early-time inflation with the late-time acceleration is possible. It is shown that a special choice of parameters for FRHL leads to the same spatially-flat FRW equations as in the case of traditional $F(R)$-gravity. Finally, an essentially most general modified Ho\v{r}ava-Lifshitz gravity is proposed, motivated by its fully diffeomorphism-invariant counterpart, with the restriction that the action does not contain derivatives higher than the second order with respect to the time coordinate.

We propose the reconstruction of $F(R)$ gravity in such a way that corresponding theory admits cyclic and ekpyrotic universe solutions. The number of explicit examples of such $F(R)$ model is found. The comparison with the reconstructed scalar-tensor theory is made. We also present $F(R)$ gravity which provides the little rip evolution and gives the realistic gravitational alternative for $\Lambda$CDM cosmology. The time for little rip dissolution of bound structures in such theory is estimated. We demonstrate that transformed little rip $F(R)$ solution becomes qualitatively different cosmological solution with Big Bang type singularity in Einstein frame.

Cosmological models with an SU(2) Yang-Mills field are studied. For a specific model with a minimally coupled Yang-Mills Lagrangian, which includes an arbitrary function of the second-order term and a fourth-order term, a corresponding reconstruction program is proposed. It is shown that the model with minimal coupling has no de Sitter solutions, for any nontrivial function of the second-order term. To get de Sitter solutions, a gravitational model with nonminimally coupled Yang-Mills fields is then investigated. It is shown that the model with non-minimal coupling has in fact a de Sitter solution, even in absence of the cosmological constant term.

We investigate de Sitter solutions in non-local gravity as well as in non-local gravity with Lagrange constraint multiplier. We examine a condition to avoid a ghost and discuss a screening scenario for a cosmological constant in de Sitter solutions. Furthermore, we explicitly demonstrate that three types of the finite-time future singularities can occur in non-local gravity and explore their properties. In addition, we evaluate the effective equation of state for the universe and show that the late-time accelerating universe may be effectively the quintessence, cosmological constant or phantom-like phases. In particular, it is found that there is a case in which a crossing of the phantom divide from the non-phantom (quintessence) phase to the phantom one can be realized when a finite-time future singularity occurs. Moreover, it is demonstrated that the addition of an $R^2$ term can cure the finite-time future singularities in non-local gravity. It is also suggested that in the framework of non-local gravity, adding an $R^2$ term leads to possible unification of the early-time inflation with the late-time cosmic acceleration.

We construct phantom energy models with the equation-of-state parameter $w$ such that $w<-1$, but finite-time future singularity does not occur. Such models can be divided into two classes: (i) energy density increases with time ("phantom energy" without "Big Rip" singularity) and (ii) energy density tends to constant value with time ("cosmological constant" with asymptotically de Sitter evolution). The disintegration of bound structure is confirmed in Little Rip cosmology. Surprisingly, we find that such disintegration (on example of Sun-Earth system) may occur even in asymptotically de Sitter phantom universe consistent with observational data. We also demonstrate that non-singular phantom models admit wormhole solutions as well as possibility of big trip via wormholes.

Dynamics and collapse of collisionless self-gravitating systems is described by the coupled collisionless Boltzmann and Poisson equations derived from $f(R)$-gravity in the weak field approximation. Specifically, we describe a system at equilibrium by a time-independent distribution function $f_0(x,v)$ and two potentials $\Phi_0(x)$ and $\Psi_0(x)$ solutions of the modified Poisson and collisionless Boltzmann equations. Considering a small perturbation from the equilibrium and linearizing the field equations, it can be obtained a dispersion relation. A dispersion equation is achieved for neutral dust-particle systems where a generalized Jeans wave-number is obtained. This analysis gives rise to unstable modes not present in the standard Jeans analysis (derived assuming Newtonian gravity as weak filed limit of $f(R)=R$). In this perspective, we discuss several self-gravitating astrophysical systems whose dynamics could be fully addressed in the framework of $f(R)$-gravity.

We examine in more detail specific models which yield a little rip cosmology, i.e., a universe in which the dark energy density increases without bound but the universe never reaches a finite-time singularity. We derive the conditions for the little rip in terms of the inertial force in the expanding universe and present two representative models to illustrate in more detail the difference between little rip models and those which are asymptotically de Sitter. We derive conditions on the equation of state parameter of the dark energy to distinguish between the two types of models. We show that coupling between dark matter and dark energy with a little rip equation of state can alter the evolution, changing the little rip into an asymptotic de Sitter expansion. We give conditions on minimally-coupled phantom scalar field models and on scalar-tensor models that indicate whether or not they correspond to a little rip expansion. We show that, counterintuitively, despite local instability, a little-rip cosmology has an infinite lifetime.

We construct a domain wall solution in $F(R)$ gravity. We reconstruct a static domain wall solution in a scalar field theory. We also reconstruct an explicit $F(R)$ gravity model in which a static domain wall solution can be realized. Moreover, we show that there could exist an effective (gravitational) domain wall in the framework of $F(R)$ gravity. In addition, it is demonstrated that a logarithmic non-minimal gravitational coupling of the electromagnetic theory in $F(R)$ gravity may produce time-variation of the fine structure constant which may increase with decrease of the curvature, and that this model would be ruled out by the constraints on the time variation of the fine structure constant from quasar absorption lines. We also present cosmological consequences of the coupling of the electromagnetic field to a scalar field as well as the scalar curvature and discuss the relation between variation of the fine structure constant and the breaking of the conformal invariance of the electromagnetic field.

Oscillations of the $F(R)$ dark energy around the phantom divide line, $\omega_{DE}=-1$, both during the matter era and also in the de Sitter epoch are investigated. The analysis during the de Sitter epoch is revisited by expanding the modified equations of motion around the de Sitter solution. Then, during the matter epoch, the time dependence of the dark energy perturbations is discussed by using two different local expansions. For high values of the red shift, the matter epoch is a stable point of the theory, giving the possibility to expand the $F(R)$-functions in terms of the dark energy perturbations. In the late-time matter era, the realistic case is considered where dark energy tends to a constant. The results obtained are confirmed by precise numerical computation on a specific model of exponential gravity. A novel and very detailed discussion is provided on the critical points in the matter era and on the relation of the oscillations with possible singularities.

Dark energy of phantom or quintessence nature with an equation of state parameter $w$ almost equal to -1 often leads the universe evolution to a finite-time future singularity. An elegant solution to this problem has been recently proposed \cite{frampton11} under the form of the so-called Little Rip cosmology which appears to be a realistic alternative to the $\Lambda$CDM model. A viscous Little Rip cosmology is here proposed. Whereas generically bulk viscosity tends to promote the Big Rip, we find that there are a number of situations where this is not the case and where the formalism nicely adjusts itself to the Little Rip scenario. We prove, in particular, that a viscous fluid (or, equivalently, one with an inhomogeneous (imperfect) equation of state) is perfectly able to produce a Little Rip cosmology as a purely viscosity effect. The possibility of its induction as a combined result of viscosity and a general (power-like) equation of state is also investigated in detail. To finish, a physical, inertial force interpretation of the dissolution of bound structures in the Little Rip cosmology is presented.

We consider f(R,T) modified theories of gravity, where the gravitational Lagrangian is given by an arbitrary function of the Ricci scalar R and of the trace of the stress-energy tensor T. We obtain the gravitational field equations in the metric formalism, as well as the equations of motion for test particles, which follow from the covariant divergence of the stress-energy tensor. Generally, the gravitational field equations depend on the nature of the matter source. The field equations of several particular models, corresponding to some explicit forms of the function f(R,T), are also presented. An important case, which is analyzed in detail, is represented by scalar field models. We write down the action and briefly consider the cosmological implications of the $f(R,T^{\phi})$ models, where $T^{\phi}$ is the trace of the stress-energy tensor of a self-interacting scalar field. The equations of motion of the test particles are also obtained from a variational principle. The motion of massive test particles is non-geodesic, and takes place in the presence of an extra force orthogonal to the four-velocity. The Newtonian limit of the equation of motion is further analyzed. Finally, we provide a constraint on the magnitude of the extra-acceleration by analyzing the perihelion precession of the planet Mercury in the framework of the present model.

Classical generalization of general relativity is considered as gravitational alternative for unified description of the early-time inflation with late-time cosmic acceleration. The structure and cosmological properties of number of modified theories, including traditional $F(R)$ and Ho\v{r}ava-Lifshitz $F(R)$ gravity, scalar-tensor theory, string-inspired and Gauss-Bonnet theory, non-local gravity, non-minimally coupled models, and power-counting renormalizable covariant gravity are discussed. Different representations and relations between such theories are investigated. It is shown that some versions of above theories may be consistent with local tests and may provide qualitatively reasonable unified description of inflation with dark energy epoch. The cosmological reconstruction of different modified gravities is made in great detail. It is demonstrated that eventually any given universe evolution may be reconstructed for the theories under consideration: the explicit reconstruction is applied to accelerating spatially-flat FRW universe. Special attention is paid to Lagrange multiplier constrained and conventional $F(R)$ gravities, for last theory the effective $\Lambda$CDM era and phantom-divide crossing acceleration are obtained. The occurrence of Big Rip and other finite-time future singularities in modified gravity is reviewed as well as its curing via the addition of higher-derivative gravitational invariants.

We formulate higher derivative gravity with Lagrange multiplier constraint and scalar projectors. Its gauge-fixed formulation as well as vector fields formulation is developed and corresponding spontaneous Lorentz symmetry breaking is investigated. We show that the only propagating mode is higher derivative graviton while scalar and vector modes do not propagate. Despite to higher derivatives structure of the action, its first FRW equation is the first order differential equation which admits the inflationary universe solution.

As first part of this work, experimental information about the decay of isotropic turbulence in ordinary hydrodynamics, u^2(t) proportional to t^{-6/5}, is used as input in FRW equations in order to investigate how an initial fraction f of turbulent kinetic energy in the cosmic fluid influences the cosmological development in the late, quintessence/phantom, universe. First order perturbative theory to the first order in f is employed. It turns out that both in the Hubble factor, and in the energy density, the influence from the turbulence fades away at late times. The divergences in these quantities near the Big Rip behave essentially as in a non-turbulent fluid. However, for the scale factor, the turbulence modification turns out to diverge logarithmically. As second part of our work, we consider the full FRW equation in which the turbulent part of the dark energy is accounted for by a separate term. It is demonstrated that turbulence occurrence may change the future universe evolution due to dissipation of dark energy. For instance, phantom-dominated universe becomes asymptotically a de Sitter one in the future, thus avoiding the Big Rip singularity.

A theory of exponential modified gravity which explains both early-time inflation and late-time acceleration, in a unified way, is proposed. The theory successfully passes the local tests and fulfills the cosmological bounds and, remarkably, the corresponding inflationary era is proven to be unstable. Numerical investigation of its late-time evolution leads to the conclusion that the corresponding dark energy epoch is not distinguishable from the one for the $\Lambda$CDM model. Several versions of this exponential gravity, sharing similar properties, are formulated. It is also shown that this theory is non-singular, being protected against the formation of finite-time future singularities. As a result, the corresponding future universe evolution asymptotically tends, in a smooth way, to de Sitter space, which turns out to be the final attractor of the system.

We investigate a curvature singularity appearing in the star collapse process in $F(R)$ gravity. In particular, we propose an understanding of the mechanism to produce the curvature singularity. Moreover, we explicitly demonstrate that $R^\alpha$ ($1 < \alpha \leq 2$) term addition could cure the curvature singularity and viable $F(R)$ gravity models could become free of such a singularity. Furthermore, we discuss the realization process of the curvature singularity and estimate the time scale of its appearance. For exponential gravity, it is shown that in case of the star collapse, the time scale is much shorter than the age of the universe, whereas in cosmological circumstances, it is as long as the cosmological time.

We investigate the hydrostatic equilibrium of stellar structure by taking into account the modi- fied La\'e-Emden equation coming out from f(R)-gravity. Such an equation is obtained in metric approach by considering the Newtonian limit of f(R)-gravity, which gives rise to a modified Poisson equation, and then introducing a relation between pressure and density with polytropic index n. The modified equation results an integro-differential equation, which, in the limit f(R) \rightarrow R, becomes the standard La\'e-Emden equation. We find the radial profiles of gravitational potential by solving for some values of n. The comparison of solutions with those coming from General Relativity shows that they are compatible and physically relevant.

This paper is devoted to the study of various aspects of projectable F(R) Horava-Lifshitz (HL) gravity. We show that some versions of F(R) HL gravity may have stable de Sitter solution and instable flat space solution. In this case, the problem of scalar graviton does not appear because flat space is not vacuum state. Generalizing the U(1) HL theory proposed in arXiv:1007.2410, we formulate U(1) extension of scalar theory and of F(R) Horava-Lifshitz gravity. The Hamiltonian approach for such the theory is developed in full detail. It is demonstrated that its Hamiltonian structure is the same as for the non-relativistic covariant HL gravity. The spectrum analysis performed around flat background indicates towards the consistency of the theory because it contains graviton with only transverse polarization. Finally, we analyze the spatially-flat FRW equations for U(1) invariant F(R) Horava-Lifshitz gravity.

We discuss a possible mechanism to screen a cosmological constant in non-local gravity. We find that in a simple model of of non-local gravity with the Lagrangian of the form, $R+f(\Box^{-1}R)-2\Lambda$ where $f(X)$ is a quadratic function of $X$, there is a flat spacetime solution despite the presence of the cosmological constant $\Lambda$. Unfortunately, however, we also find that this solution contains a ghost in general. Then we discuss the condition to avoid a ghost and find that one can avoid it only for a finite range of `time'. Nevertheless our result suggests the possibility of solving the cosmological constant problem in the context of non-local gravity.

The review of basic cosmological properties of four-dimensional F(R)-gravity, including FRW equations of motion and its accelerating solutions, generalized fluid and scalar-tensor representation of the theory is done. Cosmological reconstruction equation is written and conditions for stability of cosmological solution are discussed. The overview of realistic F(R)-models unifying inflation with dark energy epoch is made. The avoidance of finite-time future singularities in such theories via the introduction of R^2-term is studied. New realistic non-singular F(R)-gravity unifying early-time inflation with late-time acceleration is presented. The exit from inflationary era in such model may be caused by the gravitational scenario. It is demonstrated that five-dimensional F(R)-gravity considered as non-perturbative stringy effective action leads to universal relation for viscous bound ratio.

We derive a Cardy-Verlinde-like formula which relates the entropy of the closed FRW universe to its energy, and Casimir energy, for a multicomponent coupled fluid. The generalized fluid obeys an inhomogeneous equation of state. A viscous dark fluid is included, and we include also modified gravity using its fluid representation. It is demonstrated how such an expression reduces to the standard Cardy-Verlinde formula corresponding to the 2d CFT entropy only in some special cases. The dynamical entropy bound for a closed FRW universe with dark components is obtained. The universality of the dynamical entropy bound near a future singularity (of all four types), as well as near the Big Bang singularity, is investigated. It is demonstrated that except from some special cases of Type II and Type IV singularities the dynamical entropy bound is violated near the singularity even if quantum effects are taken into account. The dynamical entropy bound seems to be universal for the case of a regular universe, including the asymptotic de Sitter one.

We study FRW cosmology for a non-linear modified F(R) Horava-Lifshitz gravity which has a viable convenient counterpart. A unified description of early-time inflation and late-time acceleration is possible in this theory, but the cosmological dynamic details are generically different from the ones of the convenient viable F(R) model. Remarkably, for some specific choice of parameters they do coincide. The emergence of finite-time future singularities is investigated in detail. It is shown that these singularities can be cured by adding an extra, higher-derivative term, which turns out to be qualitatively different when compared with the corresponding one of the convenient F(R) theory.

We study scalar-tensor theory, k-essence and modified gravity with Lagrange multiplier constraint which role is to reduce the number of degrees of freedom. Dark Energy cosmology of different types ($\Lambda$CDM, unified inflation with DE, smooth non-phantom/phantom transition epoch) is reconstructed in such models. It is shown that mathematical equivalence between scalar theory and $F(R)$ gravity is broken due to presence of constraint. The cosmological dynamics of $F(R)$ gravity is modified by the second $F_2(R)$ function dictated by the constraint. Dark Energy cosmology is defined by this function while standard $F_1(R)$ function is relevant for local tests (modification of newton regime). A general discussion on the role of Lagrange multipliers to make higher-derivative gravity canonical is developed.

We continue the study of covariant power-counting renormalizable gravity constrained by scalar Lagrange multiplier. Lorentz symmetry breaking is investigated in such a theory in comparison with the one in ghost condensation model. Covariant power-counting renormalizable vector gravity which is invariant under $U(1)$ gauge symmetry is proposed. Several forms of vector Lagrange multiplier in this theory are discussed. It is shown that covariant scalar/vector gravity under consideration may have power-law or de Sitter accelerating cosmological solution corresponding to inflationary era. Simplest black hole solution is obtained and dispersion relations for graviton are presented.

Several different explicit reconstructions of f(R) gravity are obtained from the background FRW expansion history. It is shown that the only theory whose Lagrangian is a simple function of the Ricci scalar R, that admits an exact LCDM expansion history is standard General Relativity with a positive cosmological constant and the only way to obtain this behaviour of the scale factor for more general functions of R is to add additional degrees of freedom to the matter sector.

The class of covariant gravity theories which have nice ultraviolet behavior and seem to be (super)-renormalizable is proposed. The apparent breaking of Lorentz invariance occurs due to the coupling with the effective fluid which is induced by Lagrange multiplier constrained scalar field. Spatially-flat FRW cosmology for such covariant field gravity may have accelerating solutions. Renormalizable versions of more complicated modified gravity which depends on Riemann and Ricci tensor squared may be constructed in the same way.

We discuss the cosmological reconstruction in modified Gauss-Bonnet and F(R) gravities. Two alternative representations of the action (with and without auxiliary scalar) are considered. The approximate description of deceleration-acceleration transition cosmologies is reconstructed. It is shown that cosmological solution containing Big Bang and Big Rip singularities may be reconstructed only using the representation with the auxiliary field. The analytical description of the deceleration-acceleration transition cosmology in modified Gauss-Bonnet gravity is demonstrated to be impossible at sufficiently general conditions.

We study all four types of finite-time future singularities emerging in late-time accelerating (effective quintessence/phantom) era from $\mathcal{F}(R,G)$-gravity, where $R$ and $G$ are the Ricci scalar and the Gauss-Bonnet invariant, respectively. As an explicit example of $\mathcal{F}(R,G)$-gravity, we also investigate modified Gauss-Bonnet gravity, so-called $F(G)$-gravity. In particular, we reconstruct the $F(G)$-gravity and $\mathcal{F}(R,G)$-gravity models where accelerating cosmologies realizing the finite-time future singularities emerge. Furthermore, we discuss a possible way to cure the finite-time future singularities in $F(G)$-gravity and $\mathcal{F}(R,G)$-gravity by taking into account higher-order curvature corrections. The example of non-singular realistic modified Gauss-Bonnet gravity is presented. It turns out that adding such non-singular modified gravity to singular Dark Energy makes the combined theory to be non-singular one as well.

An explicit model of $F(R)$ gravity with realizing a crossing of the phantom divide is reconstructed. In particular, it is shown that the Big Rip singularity may appear in the reconstructed model of $F(R)$ gravity. Such a Big Rip singularity could be avoided by adding $R^2$ term or non-singular viable $F(R)$ theory to the model because phantom behavior becomes transient.

We consider the diffeomorphism invariant gravity coupled with the ideal fluid in the non-standard way. The Lorentz-invariance of the graviton propagator in such a theory considered as perturbation over flat background turns out to be broken due to non-standard coupling with the ideal fluid. As a result the behavior of the propagator in the ultraviolet/infrared region indicates that some versions of such theory are (super-)renormalizable ones (with appearance of only physical transverse modes). The FRW cosmology in same cases may be different from the one in General Relativity with the possible quintessence/phantom-like inflationary stage.

We explicitly show that the equations of motion for modified gravity theories of $F(R)$-gravity, the scalar-Gauss-Bonnet gravity, $F(\mathcal{G})$-gravity and the non-local gravity are equivalent to the Clausius relation in thermodynamics. In addition, we discuss the relation between the expression of the entropy and the contribution from the modified gravity as well as the matter to the definition of the energy flux (heat).

The fundamental problem of the occurrence/removal of finite-time future singularity in the universe evolution for coupled dark energy (DE) is addressed. It is demonstrated the existence of the (instable or local minimum) de Sitter space solution which may cure the Type II or Type IV future singularity for DE coupled with DM as the result of tuning the initial conditions. In case of phantom DE, the corresponding coupling may help to resolve the coincidence problem but not the Big Rip (Type I) singularity issue. We show that modified gravity of special form or inhomogeneous DE fluid may offer the universal scenario to cure the Type I,II,III or IV future singularity of coupled (fluid or scalar) DE evolution.

The 1+3 covariant approach and the covariant gauge-invariant approach to perturbations are used to analyze in depth conformal transformations in cosmology. Such techniques allow us to obtain very interesting insights on the physical content of these transformations, when applied to non-standard gravity. The results obtained lead to a number of general conclusions on the change of some key quantities describing any two conformally related cosmological models. In particular, it is shown that the physics in the Einstein frame has characteristics which are completely different from those in the Jordan frame. Even if some of the geometrical properties of the cosmology are preserved (homogeneous and isotropic Universes are mapped into homogeneous and isotropic universes), it can happen that decelerating cosmologies are mapped into accelerated ones. Differences become even more pronounced when first-order perturbations are considered: from the 1+3 equations it is seen that first-order vector and tensor perturbations are left unchanged in their structure by the conformal transformation, but this cannot be said of the scalar perturbations, which include the matter density fluctuations. Behavior in the two frames of the growth rate, as well as other evolutionary features, like the presence or absence of oscillations, etc., appear to be different too. The results obtained are then explicitly interpreted and verified with the help of some clarifying examples based on $f(R)$-gravity cosmologies.

Using the fluid representation, we formulate the conditions for the appearance of all four types finite-time future singularity in modified gravity in accelerating FRW universe. It stressed that number of standard quintessence/phantom DE theories (including scalar, fluid, DBI ones, etc) brings the accelerating cosmology to future singularity precisely in the same way as singular modified gravity DE. The viable non-singular modified gravity unifying the early-time inflation with late-time acceleration is considered. It is shown that adding such non-singular theory to another realistic singular modified gravity which has the accelerating solution with future singularity may cure the singularity of resulting combined model. This universal scenario may be naturally applied to standard singular DE models as well as to inflationary theories with initial singularity. This suggests the additional fundamental reasoning for modification of General Relativity.

The one-loop effective action for non-minimal scalar modified gravity on de Sitter background with a constant scalar field is found. The corresponding induced cosmological constant is evaluated. It is shown that quantum effects in non-minimal modified gravity may induce an early-time de Sitter universe even in the situation when such solution does not occur on the classical level. Classical reconstruction of the theory is presented in such a way that the resulting theory has a cosmological solution unifying early-time inflation with late-time acceleration.

The cosmological reconstruction scheme for modified $F(R)$ gravity is developed in terms of e-folding (or, redshift). It is demonstrated how any FRW cosmology may emerge from specific $F(R)$ theory. The specific examples of well-known cosmological evolution are reconstructed, including $\Lambda$CDM cosmology, deceleration with transition to phantom superacceleration era which may develop singularity or be transient. The application of this scheme to viable $F(R)$ gravities unifying inflation with dark energy era is proposed. The additional reconstruction of such models leads to non-leading gravitational correction mainly relevant at the early/late universe and helping to pass the cosmological bounds (if necessary). It is also shown how cosmological reconstruction scheme may be generalized in the presence of scalar field.

The future finite-time singularities emerging in alternative gravity dark energy models are classified and studied in Jordan and Einstein frames. It is shown that such singularity may occur even in flat spacetime for the specific choice of the effective potential. The conditions for the avoidance of finite-time singularities are presented and discussed. The problem is reduced to the study of a scalar field evolving on an effective potential by using the conformal transformations. Some viable modified gravity models are analyzed in detail and the way to cure singularity is considered by introducing the higher-order curvature corrections. These results maybe relevant for the resolution of the conjectured problem in the relativistic star formation in such modified gravity where finite-time singularity is also manifested.

We reconstruct an explicit model of modified gravity in which a crossing of the phantom divide can be realized. It is shown that the (finite-time) Big Rip singularity appears in the model of modified gravity (i.e., in the so-called Jordan frame), whereas that in the corresponding scalar field theory obtained through the conformal transformation (i.e., in the so-called Einstein frame) the singularity becomes the infinite-time one. Furthermore, we investigate the relations between the scalar field theories with realizing a crossing of the phantom divide and the corresponding modified gravitational theories by using the inverse conformal transformation. It is demonstrated that the scalar field theories describing the non-phantom phase (phantom one with the Big Rip) can be represented as the theories of real (complex) $F(R)$ gravity through the inverse (complex) conformal transformation. We also study a viable model of modified gravity in which the transition from the de Sitter universe to the phantom phase can occur. In addition, we explore the stability for the obtained solutions of the crossing of the phantom divide under a quantum correction coming from conformal anomaly.

We consider ideal fluid and equivalent scalar field dark energy universes where all four known types of finite-time, future singularities occur at some parameter values. It is demonstrated that pressure/energy density of such quintessence/phantom dark energy diverges in spherically-symmetric spacetime at finite radius or at the center. This may cause the instability of the relativistic star or black hole in such universe. The resolution of the problem via the extra modification of the equation of state is briefly discussed.

Viable models of modified gravity which satisfy both local as well as cosmological tests are investigated. It is demonstrated that some versions of such highly non-linear models exhibit multiply de Sitter universe solutions, which often appear in pairs, being one of them stable and the other unstable. It is explicitly shown that, for some values of the parameters, it is possible to find several de Sitter spaces (as a rule, numerically); one of them may serve for the inflationary stage, while the other can be used for the description of the dark energy epoch. The numerical evolution of the effective equation of state parameter is also presented, showing that these models can be considered as natural candidates for the unification of early-time inflation with late-time acceleration through dS critical points. Moreover, based on the de Sitter solutions, multiply SdS universes are constructed which might also appear at the (pre-)inflationary stage. Their thermodynamics are studied and free energies are compared.

We study accelerating cosmological solutions of a general class of non-linear gravities which depend on Gauss-Bonnet and other higher derivative invariants. To achieve this goal a local formulation with auxiliary scalars for arbitrary higher-derivative non-local gravity is developed. It is demonstrated that non-local Gauss-Bonnet gravity can be reduced, in the local formulation, to a model of string-inspired scalar-Gauss-Bonnet gravity. A natural unification, in the theory here developed, of the early-time inflation epoch with a late-time acceleration stage can also be realized.

We discuss different kinds of Killing horizons possible in static, spherically symmetric configurations and recently classified as "usual", "naked" and "truly naked" ones depending on the near-horizon behavior of transverse tidal forces acting on an extended body. We obtain necessary conditions for the metric to be extensible beyond a horizon in terms of an arbitrary radial coordinate and show that all truly naked horizons, as well as many of those previously characterized as naked and even usual ones, do not admit an extension and therefore must be considered as singularities. Some examples are given, showing which kinds of matter are able to create specific space-times with different kinds of horizons, including truly naked ones. Among them are fluids with negative pressure and scalar fields with a particular behavior of the potential. We also discuss horizons and singularities in Kantowski-Sachs spherically symmetric cosmologies and present horizon regularity conditions in terms of an arbitrary time coordinate and proper (synchronous) time. It turns out that horizons of orders 2 and higher occur in infinite proper times in the past or future, but one-way communication with regions beyond such horizons is still possible.

We investigate the future evolution of the dark energy universe in modified gravities including $F(R)$ gravity, string-inspired scalar-Gauss-Bonnet and modified Gauss-Bonnet ones, and ideal fluid with the inhomogeneous equation of state (EoS). Modified Friedmann-Robertson-Walker (FRW) dynamics for all these theories may be presented in universal form by using the effective ideal fluid with an inhomogeneous EoS without specifying its explicit form. We construct several examples of the modified gravity which produces accelerating cosmologies ending at the finite-time future singularity of all four known types by applying the reconstruction program. Some scenarios to resolve the finite-time future singularity are presented. Among these scenarios, the most natural one is related with additional modification of the gravitational action in the early universe. In addition, late-time cosmology in the non-minimal Maxwell-Einstein theory is considered. We investigate the forms of the non-minimal gravitational coupling which generates the finite-time future singularities and the general conditions for this coupling in order that the finite-time future singularities cannot emerge. Furthermore, it is shown that the non-minimal gravitational coupling can remove the finite-time future singularities or make the singularity stronger (or weaker) in modified gravity.

We review the unification of early-time inflation with late-time acceleration in several local modified gravity models which pass Solar System and cosmological tests. It is also demonstrated that account of non-local gravitational corrections to the action does not destroy the possibility of such unification. Dark matter effect is caused by composite graviton degree of freedom in such models.

We analyze seven different viable $f(R)$-gravities towards the Solar System tests and stochastic gravitational waves background. The aim is to achieve experimental bounds for the theory at local and cosmological scales in order to select models capable of addressing the accelerating cosmological expansion without cosmological constant but evading the weak field constraints. Beside large scale structure and galactic dynamics, these bounds can be considered complimentary in order to select self-consistent theories of gravity working at the infrared limit. It is demonstrated that seven viable $f(R)$-gravities under consideration not only satisfy the local tests, but additionally, pass the above PPN-and stochastic gravitational waves bounds for large classes of parameters.

We review modified $F(R)$ gravity as realistic candidate to describe the observable universe expansion history. We show that recent cosmic acceleration, radiation/matter-dominated epoch and inflation could be realized in the framework of $F(R)$-gravity in the unified way. For some viable classes of $F(R)$-gravity, the Newton law is respected and there is no so-called matter instability (the very heavy positive mass for additional scalar degree of freedom is generated). The reconstruction program in modified gravity is also reviewed and it is demonstrated that {\it any} time-evolution of the universe expansion could be realized in $F(R)$-gravity. These models remain to be realistic also in the presence of non-minimal gravitational coupling with usual matter. It is shown that same model which passes local tests and predicts the unification of inflation with cosmic acceleration also describes dark matter thanks to presence of additional scalar degree of freedom and chameleon mechanism.

We study the future evolution of quintessence/phantom dominated epoch in modified $F(R)$-gravity which unifies the early-time inflation with late-time acceleration and which is consistent with observational tests. Using the reconstruction technique it is demonstrated that there are models where any known (Big Rip, II, III or IV Type) singularity may classically occur. From another side, in Einstein frame (scalar-tensor description) only IV Type singularity occurs. Near the singularity the classical description breaks up, it is demonstrated that quantum effects act against the singularity and may prevent its appearance. The realistic $F(R)$-gravity which is future singularity free is proposed. We point out that additional modification of any $F(R)$-gravity by the terms relevant at the early universe is possible, in such a way that future singularity does not occur even classically.

We study inflationary cosmology and the late-time accelerated expansion of the universe in non-minimal Yang-Mills (YM) theory, in which the YM field couples to a function of the scalar curvature. It is shown that power-law inflation can be realized due to the non-minimal YM field-gravitational coupling which maybe caused by quantum corrections. Moreover, it is demonstrated that both inflation and the late-time accelerated expansion of the universe can be realized in a modified YM-$F(R)$ gravity which is consistent with solar system tests. Furthermore, it is shown that this result can be realized also in a non-minimal vector-$F(R)$ gravity. In addition, we consider the duality of the non-minimal electromagnetic theory and that of the non-minimal YM theory, and also discuss the cosmological reconstruction of the YM theory.

In this paper we have investigated the cosmological dynamics of non-locally corrected gravity involving a function of the inverse d'Alembertian of the Ricci scalar, $f(\Box^{-1} R))$. Casting the dynamical equations into local form, we derive the fixed points of the dynamics and demonstrate the existence and stability of a one parameter family of dark energy solutions for a simple choice, $f(\Box^{-1} R)\sim \exp(\alpha \Box^{-1} R)$. The effective EoS parameter is given by, $w_{\rm eff}=({\alpha-1})/({3\alpha-1})$ and the stability of the solutions is guaranteed provided that $1/3<\alpha<2/3$. For $1/3<\alpha<1/2$ and $1/2<\alpha<2/3$, the underlying system exhibits phantom and non-phantom behavior respectively; the de Sitter solution corresponds to $\alpha=1/2$. For a wide range of initial conditions, the system mimics dust like behavior before reaching the stable fixed point. The late time phantom phase is achieved without involving negative kinetic energy fields. A brief discussion on the entropy of de Sitter space in non-local model is included.

We study inflation and late-time acceleration in the expansion of the universe in non-minimal electromagnetism, in which the electromagnetic field couples to the scalar curvature function. It is shown that power-law inflation can be realized due to the non-minimal gravitational coupling of the electromagnetic field, and that large-scale magnetic fields can be generated due to the breaking of the conformal invariance of the electromagnetic field through its non-minimal gravitational coupling. Furthermore, it is demonstrated that both inflation and the late-time acceleration of the universe can be realized in a modified Maxwell-$F(R)$ gravity which is consistent with solar system tests and cosmological bounds and free of instabilities. At small curvature typical for current universe the standard Maxwell theory is recovered. We also consider classically equivalent form of non-minimal Maxwell-$F(R)$ gravity, and propose the origin of the non-minimal gravitational coupling function based on renormalization-group considerations.

We review on the viability of $F(R)$-gravity. We show that recent cosmic acceleration, radiation/matter-dominated epoch and inflation could be realized in the framework of $F(R)$-gravity in the unified way. For some classes of $F(R)$-gravity, the correction to the Newton law is extremely small and there is no so-called matter instability (the very heavy positive mass for additional scalar degree of freedom is generated). The reconstruction program in modified gravity is also reviewed and it is demonstrated that {\it any} time-evolution of the universe expansion could be realized in $F(R)$-gravity. Special attention is paid to modified gravity which unifies inflation with cosmic acceleration and passes local tests. It turns out that such a theory may describe also dark matter.

A general approach to viable modified $f(R)$ gravity is developed in both the Jordan and the Einstein frames. A class of exponential, realistic modified gravities is introduced and investigated with care. Special focus is made on step-class models, most promising from the phenomenological viewpoint and which provide a natural way to classify all viable modified gravities. One- and two-steps models are explicitly considered, but the analysis is extensible to $N$-step models. Both inflation in the early universe and the onset of recent accelerated expansion arise in these models in a natural, unified way. Moreover, it is demonstrated that models in this category easily pass all local tests, including stability of spherical body solution, non-violation of Newton's law, and generation of a very heavy positive mass for the additional scalar degree of freedom.

We consider modified $f(R)$ gravity which may unify $R^m$ early-time inflation with late-time $\Lambda$CDM epoch. It is shown that such model passes the local tests (Newton law, stability of Earth-like gravitational solution, very heavy mass for additional scalar degree of freedom) and suggests the realistic alternative for General Relativity. Various scenarios for future evolution of $f(R)$ $\Lambda$CDM era are discussed.

We consider FRW cosmology in non-local modified gravity. Its local scalar-tensor formulation is developed. It is explicitly demonstrated that such theory may lead to the unification of early-time inflation with late-time cosmic acceleration. The quintessence or phantom era may emerge for specific form of the action. The coupled non-local-F(R) gravity is also investigated. It is shown that such theory being consistent with Solar System tests may lead to the known universe history sequence: inflation, radiation/matter dominance and dark epoch.

We consider the modified gravity non-minimally coupled with matter Lagrangian for the description of early-time and late-time universe. Such $F(R)$ ($F(G)$) gravity in the absence of non-minimal coupling is viable theory which passes the local tests and reproduces the $\Lambda$CDM era. For qualitatively similar choice of non-minimal gravitational coupling function it is shown that the unified description of early-time inflation and late-time cosmic acceleration is possible. It is interesting that matter (scalar) which supports the inflationary era is gravitationally screened at late times. Hence, it may be effectively invisible at current universe.

We investigate the possibility of a dark energy universe emerging from an action with higher-order string loop corrections to Einstein gravity in the presence of a massless dilaton. These curvature corrections (up to $R^4$ order) are different depending upon the type of (super)string model which is considered. We find in fact that Type II, heterotic, and bosonic strings respond differently to dark energy. A dark energy solution is shown to exist in the case of the bosonic string, while the other two theories do not lead to realistic dark energy universes. Detailed analysis of the dynamical stability of the de-Sitter solution is presented for the case of a bosonic string. A general prescription for the construction of a de-Sitter solution for the low-energy (super)string effective action is also indicated. Beyond the low-energy (super)string effective action, when the higher-curvature correction coefficients depend on the dilaton, the reconstruction of the theory from the universe expansion history is done with a corresponding prescription for the scalar potentials.

We suggest two realistic f(R) and one F(G) modified gravities which are consistent with local tests and cosmological bounds. The typical property of such theories is the presence of the effective cosmological constant epochs in such a way that early-time inflation and late-time cosmic acceleration are naturally unified within single model. It is shown that classical instability does not appear here and Newton law is respected. Some discussion of possible anti-gravity regime appearence and related modification of the theory is done.

We consider class of modified $f(R)$ gravities with the effective cosmological constant epoch at the early and late universe. Such models pass most of solar system tests as well they satisfy to cosmological bounds. Despite their very attractive properties, it is shown that one realistic class of such models may lead to significant Newton law corrections which become small at the future universe only. The model with acceptable Newton law regime shows the matter instability. This suggests that more complicated version of such theory (or extended parameters space) should be investigated.

The modified F(R)-scalar-Gauss-Bonnet gravity is proposed as dark energy model. The reconstruction program for such theory is developed. It is explicitly demonstrated that the known classical universe expansion history (deceleration epoch, transition to acceleration and effective quintessence, phantom or cosmological constant era) may naturally occur in such unified theory for some (reconstructed) classes of scalar potentials. Gauss-Bonnet assisted dark energy is also proposed. The possibility of cosmic acceleration is studied there.

We consider scalar-Gauss-Bonnet and modified Gauss-Bonnet gravities and reconstruct these theories from the universe expansion history. In particular, we are able to construct versions of those theories (with and without ordinary matter), in which the matter dominated era makes a transition to the cosmic acceleration epoch. It is remarkable that, in several of the cases under consideration, matter dominance and the deceleration-acceleration transition occur in the presence of matter only. The late-time acceleration epoch is described asymptotically by de Sitter space but may also correspond to an exact $\Lambda$CDM cosmology, having in both cases an effective equation of state parameter $w$ close to -1. The one-loop effective action of modified Gauss-Bonnet gravity on the de Sitter background is evaluated and it is used to derive stability criteria for the ensuing de Sitter universe.

We discuss FRW universe with time-dependent EoS dark fluid which leads to multiple de Sitter space. This model (as well as its scalar-tensor version) may be considered as some classical analog of cosmological landscape. The universe expansion history may look as transitions between different deSitter eras what suggests the interesting solution for cosmological constant problem. For specific time-dependent EoS dark fluid the possibility of transitions between universe regions with positive and negative cosmological constant is also established.

It is shown that phantom scalar models can be mapped into a mathematically equivalent, modified $F(R)$ gravity, which turns out to be complex, in general. Only for even scalar potentials is the ensuing modified gravity real. It is also demonstrated that, even in this case, modified gravity becomes complex at the region where the original phantom dark energy theory develops a Big Rip singularity. A number of explicit examples are presented which show that these two theories are not completely equivalent, from the physical viewpoint. This basically owes to the fact that the physical metric in both theories differ in a time-dependent conformal factor. As a result, an FRW accelerating solution, or FRW instanton, in the scalar-tensor theory may look as a deccelerating FRW solution, or a non-instantonic one, in the corresponding modified gravity theory.

We develop the reconstruction program for the number of modified gravities: scalar-tensor theory, $f(R)$, $F(G)$ and string-inspired, scalar-Gauss-Bonnet gravity. The known (classical) universe expansion history is used for the explicit and successful reconstruction of some versions (of special form or with specific potentials) from all above modified gravities. It is demonstrated that cosmological sequence of matter dominance, decceleration-acceleration transition and acceleration era may always emerge as cosmological solutions of such theory. Moreover, the late-time dark energy FRW universe may have the approximate or exact $\Lambda$CDM form consistent with three years WMAP data. The principal possibility to extend this reconstruction scheme to include the radiation dominated era and inflation is briefly mentioned. Finally, it is indicated how even modified gravity which does not describe the matter-dominated epoch may have such a solution before acceleration era at the price of the introduction of compensating dark energy.

We develop the general program of the unification of matter-dominated era with acceleration epoch for scalar-tensor theory or dark fluid. The general reconstruction of single scalar-tensor theory is fulfilled. The explicit form of scalar potential for which the theory admits matter-dominated era, transition to acceleration and (asymptotically deSitter) acceleration epoch consistent with WMAP data is found. The interrelation of the epochs of deceleration-acceleration transition and matter dominance-dark energy transition for dark fluids with general EOS is investigated. We give several examples of such models with explicit EOS (using redshift parametrization) where matter-dark energy domination transition may precede the deceleration-acceleration transition. As some by-product, the reconstruction scheme is applied to scalar-tensor theory to define the scalar potentials which may produce the dark matter effect. The obtained modification of Newton potential may explain the rotation curves of galaxies.

We develop the general scheme for modified $f(R)$ gravity reconstruction from any realistic FRW cosmology. We formulate several versions of modified gravity compatible with Solar System tests where the following sequence of cosmological epochs occurs: a. matter dominated phase (with or without usual matter), transition from decceleration to acceleration, accelerating epoch consistent with recent WMAP data b. $\Lambda$CDM cosmology without cosmological constant. As a rule, such modified gravities are expressed implicitly (in terms of special functions) with late-time asymptotics of known type (for instance, the model with negative and positive powers of curvature). In the alternative approach, it is demonstrated that even simple versions of modified gravity may lead to the unification of matter dominated and accelerated phases at the price of the introduction of compensating dark energy.

In this paper we investigate the cosmological effects of modified gravity with string curvature corrections added to Einstein-Hilbert action in the presence of a dynamically evolving scalar field coupled to Riemann invariants. The scenario exhibits several features of cosmological interest for late universe. It is shown that higher order stringy corrections can lead to a class of dark energy models consistent with recent observations. The model can give rise to quintessence, deSitter or phantom dark energy, in last case without recourse to negative kinetic energy field. The detailed treatment of reconstruction program for general scalar-Gauss-Bonnet gravity is presented for any given cosmology. The explicit examples of reconstructed scalar potentials are given for accelerated (quintessence, cosmological constant or phantom) universe. Finally, the relation with modified $F(G)$ gravity is established on classical level and is extended to include third order terms on curvature.

We suggest to generalize the dark energy equation of state (EoS) by introduction the relaxation equation for pressure which is equivalent to consideration of the inhomogeneous EoS cosmic fluid which often appears as the effective model from strings/brane-worlds. As another, more wide generalization we discuss the inhomogeneous EoS which contains derivatives of pressure. For several explicit examples motivated by the analogy with classical mechanics the accelerating FRW cosmology is constructed. It turns out to be the asymptotically de Sitter or oscillating universe with possible transition from deceleration to acceleration phase. The coupling of dark energy with matter in accelerating FRW universe is considered, it is shown to be consistent with constrained (or inhomogeneous) EoS.

We show that f(R)-gravity can, in general, give rise to cosmological viable models compatible with a matter-dominated epoch evolving into a late accelerated phase. We discuss the various representations of f(R)-gravity as an ideal fluid or a scalar-tensor gravity theory, taking into account conformal transformations. We point out that mathematical equivalence does not correspond, in several cases, to the physical equivalence of Jordan frame and Einstein frame. Finally, we show that wide classes of f(R)-gravity models, including matter and accelerated phases, can be phenomenologically reconstructed by means of observational data. In principle, any popular quintessence models could be "reframed" as an f(R)-gravity model.

We consider the oscillating dark energy with periodic equation of state in two equivalent formulations: ideal fluid or scalar-tensor theory. It is shown that such dark energy suggests the natural way for the unification of early-time inflation with late-time acceleration. We demonstrate how it describes the transition from deceleration to acceleration or from non-phantom to phantom era and how it solves the coincidence problem. The occurence of finite-time future singularity for the oscillating (phantom) universe is also investigated.

The unifying approach to early-time and late-time universe based on phantom cosmology is proposed. We consider gravity-scalar system which contains usual potential and scalar coupling function in front of kinetic term. As a result, the possibility of phantom-non-phantom transition appears in such a way that universe could have effectively phantom equation of state at early time as well as at late time. In fact, the oscillating universe may have several phantom and non-phantom phases. As a second model we suggest generalized holographic dark energy where infrared cutoff is identified with combination of FRW parameters: Hubble constant, particle and future horizons, cosmological constant and universe life-time (if finite). Depending on the specific choice of the model the number of interesting effects occur: the possibility to solve the coincidence problem, crossing of phantom divide and unification of early-time inflationary and late-time accelerating phantom universe. The bound for holographic entropy which decreases in phantom era is also discussed.

We review various modified gravities considered as gravitational alternative for dark energy. Specifically, we consider the versions of $f(R)$, $f(G)$ or $f(R,G)$ gravity, model with non-linear gravitational coupling or string-inspired model with Gauss-Bonnet-dilaton coupling in the late universe where they lead to cosmic speed-up. It is shown that some of such theories may pass the Solar System tests. On the same time, it is demonstrated that they have quite rich cosmological structure: they may naturally describe the effective (cosmological constant, quintessence or phantom) late-time era with a possible transition from decceleration to acceleration thanks to gravitational terms which increase with scalar curvature decrease. The possibility to explain the coincidence problem as the manifestation of the universe expansion in such models is mentioned. The late (phantom or quintessence) universe filled with dark fluid with inhomogeneous equation of state (where inhomogeneous terms are originated from the modified gravity) is also described.

Dark energy cosmology is considered in a modified Gauss-Bonnet (GB) model of gravity where an arbitrary function of the GB invariant, $f(G)$, is added to the General Relativity action. We show that such theory is endowed with a quite rich cosmological structure: it may naturally lead to an effective cosmological constant, quintessence or phantom cosmic acceleration, with a possible transition from deceleration to acceleration. It is demonstrated in the paper that this theory is perfectly viable, since it is compliant with Solar System constraints. Specific properties of $f(G)$ gravity in a de Sitter universe, such as dS and SdS solutions, their entropy and its explicit one-loop quantization are studied. The issue of a possible solution of the hierarchy problem in modified gravities is addressed too.

We investigate the effects of viscosity terms depending on the Hubble parameter and its derivatives in the dark energy equation of state. Such terms are possible if dark energy is a fictitious fluid originating from corrections to the Einstein general relativity as is the case for some braneworld inspired models or fourth order gravity. We consider two classes of models whose singularities in the early and late time universe have been studied by testing the models against the dimensionless coordinate distance to Type Ia Supernovae and radio-galaxies also including priors on the shift and the acoustic peak parameters. It turns out that both models are able to explain the observed cosmic speed up without the need of phantom (w < -1) dark energy.

Dark energy dynamics of the universe can be achieved by equivalent mathematical descriptions taking into account generalized fluid equations of state in General Relativity, scalar-tensor theories or modified F(R) gravity in Einstein or Jordan frames. The corresponding technique transforming equation of state description to scalar-tensor or modified gravity is explicitly presented. We show that such equivalent pictures can be discriminated by matching solutions with data capable of selecting the true physical frame.

The solution of dark energy problem in the models without scalars is presented. It is shown that late-time accelerating cosmology may be generated by the ideal fluid with some implicit equation of state. The universe evolution within modified Gauss-Bonnet gravity is considered. It is demonstrated that such gravitational approach may predict the (quintessential, cosmological constant or transient phantom) acceleration of the late-time universe with natural transiton from deceleration to acceleration (or from non-phantom to phantom era in the last case).

Phantom cosmology allows to account for dynamics and matter content of the universe tracing back the evolution to the inflationary epoch, considering the transition to the non-phantom standard cosmology (radiation/matter dominated eras) and recovering the today observed dark energy epoch. We develop the unified phantom cosmology where the same scalar plays the role of early time (phantom) inflaton and late-time Dark Energy. The recent transition from decelerating to accelerating phase is described too by the same scalar field. The (dark) matter may be embedded in this scheme, giving the natural solution of the coincidence problem. It is explained how the proposed unified phantom cosmology can be fitted against the observations which opens the way to define all the important parameters of the model.

We suggest the modified gravity where some arbitrary function of Gauss-Bonnet (GB) term is added to Einstein action as gravitational dark energy. It is shown that such theory may pass solar system tests. It is demonstrated that modified GB gravity may describe the most interesting features of late-time cosmology: the transition from deceleration to acceleration, crossing the phantom divide, current acceleration with effective (cosmological constant, quintessence or phantom) equation of state of the universe.

The current accelerated universe could be produced by modified gravitational dynamics as it can be seen in particular in its Palatini formulation. We analyze here a specific non-linear gravity-scalar system in the first order Palatini formalism which leads to a FRW cosmology different from the purely metric one. It is shown that the emerging FRW cosmology may lead either to an effective quintessence phase (cosmic speed-up) or to an effective phantom phase. Moreover, the already known gravity assisted dark energy dominance occurs also in the first order formalism. Finally, it is shown that a dynamical theory able to resolve the cosmological constant problem exists also in this formalism, in close parallel with the standard metric formulation.

The dark energy universe equation of state (EOS) with inhomogeneous,Hubble parameter dependent term is considered. The motivation to introduce such a term comes from time-dependent viscosity considerations and modifications of general relativity. For several explicit examples of such EOS it is demonstrated how the type of future singularity changes, how the phantom epoch emerges and how crossing of phantom barrier occurs. Similar cosmological regimes are considered for the universe with two interacting fluids and for universe with implicit EOS. For instance, the crossing of phantom barrier is realized in easier way, thanks to the presence of inhomogeneous term. The thermodynamical dark energy model is presented where the universe entropy may be positive even at phantom era as a result of crossing of w=-1 barrier.

The functional formulation and one-loop effective action for scalar self-interacting theory non-linearly coupled with some power of the curvature are studied. After the explicit one-loop renormalization at weak curvature, we investigated numerically the phase structure for such unusual phi^4 theory. It is demonstrated the possibility of curvature-induced phase transitions for positive values of the curvature power, while for negative values the radiative symmetry breaking does not take place. The dynamical mechanism for the explanation of the current smallness of the cosmological constant is presented for several models from the class of theories under consideration.

We propose the Gauss-Bonnet dark energy model inspired by string/M-theory where standard gravity with scalar contains additional scalar-dependent coupling with Gauss-Bonnet invariant. It is demonstrated that effective phantom (or quintessence) phase of late universe may occur in the presence of such term when the scalar is phantom or for non-zero potential (for canonical scalar). However, with the increase of the curvature the GB term may become dominant so that phantom phase is transient and $w=-1$ barrier may be passed. Hence, the current acceleration of the universe may be caused by mixture of scalar phantom and (or) potential/stringy effects. It is remarkable that scalar-Gauss-Bonnet coupling acts against the Big Rip occurence in phantom cosmology.

We aim at the construction of dark energy models without exotic matter but with a phantom-like equation of state (an effective phantom phase). The first model we consider is decaying vacuum cosmology where the fluctuations of the vacuum are taken into account. In this case, the phantom cosmology (with an effective, observational $\omega$ being less than -1) emerges even for the case of a real dark energy with a physical equation of state parameter $\omega$ larger than -1. The second proposal is a generalized holographic model, which is produced by the presence of an infrared cutoff. It also leads to an effective phantom phase, which is not a transient one as in the first model. However, we show that quantum effects are able to prevent its evolution towards a Big Rip singularity.

The properties of future singularities are investigated in the universe dominated by dark energy including the phantom-type fluid. We classify the finite-time singularities into four classes and explicitly present the models which give rise to these singularities by assuming the form of the equation of state of dark energy. We show the existence of a stable fixed point with an equation of state $w<-1$ and numerically confirm that this is actually a late-time attractor in the phantom-dominated universe. We also construct a phantom dark energy scenario coupled to dark matter that reproduces singular behaviors of the Big Rip type for the energy density and the curvature of the universe. The effect of quantum corrections coming from conformal anomaly can be important when the curvature grows large, which typically moderates the finite-time singularities.

Motivated by the dark energy issue, the one-loop quantization approach for a family of relativistic cosmological theories is discussed in some detail. Specifically, general $f(R)$ gravity at the one-loop level in a de Sitter universe is investigated, extending a similar program developed for the case of pure Einstein gravity. Using generalized zeta regularization, the one-loop effective action is explicitly obtained off-shell, what allows to study in detail the possibility of (de)stabilization of the de Sitter background by quantum effects. The one-loop effective action maybe useful also for the study of constant curvature black hole nucleation rate and it provides the plausible way of resolving the cosmological constant problem.

We discuss the modified gravity which includes negative and positive powers of the curvature and which provides the gravitational dark energy. It is shown that in GR plus the term containing negative power of the curvature the cosmic speed-up may be achieved, while the effective phantom phase (with $w$ less than -1) follows when such term contains the fractional positive power of the curvature. The minimal coupling with matter makes the situation more interesting: even 1/R theory coupled with the usual ideal fliud may describe the (effective phantom) dark energy. The account of $R^2$ term (consistent modified gravity) may help to escape of cosmic doomsday.

We review some modified gravity models which describe the gravitational dark energy and the possibility of cosmic speed-up. The new consistent version of such theory which contains inverse and HD curvature terms as well as new type of coupling with matter is proposed. The accelerating cosmologies are discussed there. The structure of finite-time (sudden) singularities is investigated.

As it follows from the classical analysis, the typical final state of the dark energy universe where dominant energy condition is violated is finite time, sudden future singularity (Big Rip). For a number of dark energy universes (including scalar phantom and effective phantom theories as well as specific quintessence model) we demonstrate that quantum effects play the dominant role near Big Rip, driving the universe out of future singularity (or, at least, making it milder). As a consequence, the entropy bounds with quantum corrections become well-defined near Big Rip. Similarly, black holes mass loss due to phantom accretion is not so dramatic as it was expected: masses do not vanish to zero due to transient character of phantom evolution stage. Some examples of cosmological evolution for negative, time-dependent equation of state are also considered with the same conclusions. The application of negative entropy (or negative temparature) occurence in the phantom thermodynamics is briefly discussed.

The creation of the inflationary brane universe in 5d bulk Einstein and Einstein-Gauss-Bonnet gravity is considered. We demonstrate that emerging universe is ambigious due to arbitrary function dependence of the junction conditions (or freedom in the choice of boundary terms). We argue that some fundamental physical principle (which may be related with AdS/CFT correspondence) is necessary in order to fix the 4d geometry in unique way.

It is proposed that dark energy may become dominant over standard matter due to universe expansion (curvature decrease). Two models: non-linear gravity-matter system and modified gravity may provide the effective phantom or effective quintessence dark energy which complies with the conjecture. It is interesting that future of such universe is not necessary finite time singularity (Big Rip). The effective quintessence naturally describes current cosmic speed-up.

We consider finite-time, future (sudden or Big Rip type) singularities which may occur even when strong energy condition is not violated but equation of state parameter is time-dependent. Recently, example of such singularity has been presented by Barrow, we found another example of it. Taking into account back reaction of conformal quantum fields near singularity, it is shown explicitly that quantum effects may delay (or make milder) the singularity. It is argued that if the evolution to singularity is realistic, due to quantum effects the universe may end up in deSitter phase before scale factor blows up. This picture is generalized for braneworld where sudden singularity may occur on the brane with qualitatively similar conclusions.

We consider late-time cosmology in a (phantom) scalar-tensor theory with an exponential potential, as a dark energy model with equation of state parameter close to -1 (a bit above or below this value). Scalar (and also other kinds of) matter can be easily taken into account. An exact spatially-flat FRW cosmology is constructed for such theory, which admits (eternal or transient) acceleration phases for the current universe, in correspondence with observational results. Some remarks on the possible origin of the phantom, starting from a more fundamental theory, are also made. It is shown that quantum gravity effects may prevent (or, at least, delay or soften) the cosmic doomsday catastrophe associated with the phantom, i.e. the otherwise unavoidable finite-time future singularity (Big Rip). A novel dark energy model (higher-derivative scalar-tensor theory) is introduced and it is shown to admit an effective phantom/quintessence description with a transient acceleration phase. In this case, gravity favors that an initially insignificant portion of dark energy becomes dominant over the standard matter/radiation components in the evolution process.

We study the entropy of a FRW universe filled with dark energy (cosmological constant, quintessence or phantom). For general or time-dependent equation of state $p=w\rho$ the entropy is expressed in terms of energy, Casimir energy, and $w$. The correspondent expression reminds one about 2d CFT entropy only for conformal matter. At the same time, the cosmological Cardy-Verlinde formula relating three typical FRW universe entropies remains to be universal for any type of matter. The same conclusions hold in modified gravity which represents gravitational alternative for dark energy and which contains terms growing at low curvature. It is interesting that BHs in modified gravity are more entropic than in Einstein gravity. Finally, some hydrodynamical examples testing new shear viscosity bound, which is expected to be the consequence of the holographic entropy bound, are presented for the early universe in the plasma era and for the Kasner metric. It seems that the Kasner metric provides a counterexample to the new shear viscosity bound.

We propose a multi-graviton theory with non-nearest-neighbor couplings in the theory space. The resulting four-dimensional discrete mass spectrum reflects the structure of a latticized extra dimension. For a plausible mass spectrum motivated by the discretized Randall-Sundrum brane-world, the induced cosmological constant turns out to be positive and may serve as a quite simple model for the dark energy of our accelerating universe.

The modified gravity with $\ln R$ or $R^{-n} (\ln R)^m$ terms which grow at small curvature is discussed. It is shown that such a model which has well-defined newtonian limit may eliminate the need for dark energy and may provide the current cosmic acceleration. It is demonstrated that $R^2$ terms are important not only for early time inflation but also to avoid the instabilities and the linear growth of the gravitational force. It is very interesting that the condition of no linear growth for gravitational force coincides with the one for scalar mass in the equivalent scalar-tensor theory to be very large. Thus, modified gravity with $R^2$ term seems to be viable classical theory.

One possibility to explain the current accelerated expansion of the universe may be related with the presence of cosmologically evolving scalar whose mass depends on the local matter density (chameleon cosmology). We point out that matter quantum effects in such scalar-tensor theory produce the chameleon scalar field dependent conformal anomaly. Such conformal anomaly adds higher derivative terms to chameleon field equation of motion. As a result, the principal possibility for instabilities appears. These instabilities seem to be irrelevant at small curvature but may become dangerous in the regions where gravitational field is strong.

We discuss the modified gravity which may produce the current cosmic aceleration of the universe and eliminates the need for dark energy. It is shown that such models where the action quickly grows with the decrease of the curvature define the FRW universe with the minimal curvature. It is required the infinite time to reach the minimal curvature during the universe evolution. It is demonstrated that quantum effects of conformal fields may strongly suppress the instabilities discovered in modified gravity. We also briefly speculate on the modification of gravity combined with the presence of the cosmological constant dark energy.

It is suggested \cite{CDTT} that current cosmic acceleration arises due to modification of General Relativity by the terms with negative powers of curvature. We show that time-dependent (hyperbolic) compactifications of string/M-theory lead to the effective 4d gravity which naturally contains such terms. The same may be achieved in braneworld by the proper choice of the boundary action. Hence, such a model which seems to eliminate the need for dark energy may have the origin in M-theory.

The modified gravity, which eliminates the need for dark energy and which seems to be stable, is considered. The terms with positive powers of the curvature support the inflationary epoch while the terms with negative powers of the curvature serve as effective dark energy, supporting current cosmic acceleration. The equivalent scalar-tensor gravity may be compatible with the simplest solar system experiments.

Some cosmological implications of ultraviolet quantum effects leading to a condensation of Born-Infeld matter are considered. It is shown that under very general conditions the quantum condensate can not act as phantom matter if its energy density is positive. On the other hand, it behaves as an effective cosmological constant in the limit where quantum induced contributions to the energy-momentum tensor dominate over the classical effects.

We discuss the model consisting of tachyon which may have the negative kinetic energy plus scalar phantom and plus conformal quantum matter. It is demonstrated that such a model naturally admits two deSitter phases where the early universe inflation is produced by quantum effects and the late time accelerating universe is caused by phantom/tachyon. The energy conditions bounds for such cosmology are derived. It is interesting that effective equation of state may change its sign which depends from the proper choice of the combination of phantom/tachyon and quantum effects.

We consider deSitter universe and Nariai universe induced by quantum CFT with classical phantom matter and perfect fluid. The model represents the combination of trace-anomaly driven inflation and phantom driven deSitter universe. The similarity of phantom matter with quantum CFT indicates that phantom scalar may be the effective description for some quantum field theory. It is demonstrated that it is easier to achieve the acceleration of the scale factor preserving the energy conditions in such unified model. Some properties of unified theory (anti-gravitating solutions, negative ADM mass Nariai solution, relation with steady state) are briefly mentioned.

Five-dimensional braneworld cosmology with deSitter (inflationary) brane universe induced by classical and quantum matter is discussed. It is shown that negative energy phantom field with quantum CFT supports the creation of deSitter universe. On the same time, pure phantom or dust with quantum effects, or Chaplygin gas with quantum effects may naturally lead to the occurence of Anti-deSitter brane universe but not deSitter one. It is also interesting that unlike to four-dimensional gravity, for phantom with (or without) quantum contribution the standard cosmological energy conditions may be effectively satisfied.

Thermodynamics of 5d SdS black hole is considered. Thermal fluctuations define the (sub-dominant) logarithmic corrections to black hole entropy and then to Cardy-Verlinde formula and to FRW brane cosmology. We demonstrate that logarithmic terms (which play the role of effective cosmological constant) change the behavior of 4d spherical brane in dS, SdS or Nariai bulk. In particularly, bounce Universe occurs or 4d dS brane expands to its maximum and then shrinks. The entropy bounds are also modified by next-to-leading terms. Out of braneworld context the logarithmic terms may suggest slight modification of standard FRW cosmology.

The AdS/CFT correspondence is employed to derive logarithmic corrections to the Cardy-Verlinde formula when thermal fluctuations in the Anti-de Sitter black hole are accounted for. The qualitative effect of these corrections on the braneworld cosmology is investigated. The role of such terms in enabling a contracting universe to undergo a bounce is demonstrated. Their influence on the stability of black holes in AdS space and the Hawking-Page-Witten phase transitions is also discussed.