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

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

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.

We analyze gravitational waves propagating in an isotropic cosmic fluid endowed with a bulk viscosity $\zeta$ and a shear viscosity $\eta$, assuming these coefficients to vary with fluid density $\rho$ as $\rho^\lambda$, with $\lambda=1/2$ favored by experimental evidence. We give the general governing equation for the gravitational waves, and focus thereafter on two examples. The first concerns waves in the very late universe, close to the Big Rip, where the fate of the comic fluid is dependent highly on the values of the parameters. Our second example considers the very early universe, the lepton era; the motivation for this choice being that the microscopical bulk viscosity as calculated from statistical mechanics is then at maximum. We find that the gravitational waves on such an underlying medium are damped, having a decay constant equal to the inverse of the conformal Hubble parameter. Our results turn out to be in good agreement with other viscosity-based approaches.

Recently we have proposed models of topological field theory including gravity in Mod. Phys. Lett. A 31 (2016) no.37, 1650213 and Phys. Rev. D 96 (2017) no.2, 024009, in order to solve the problem of the cosmological constant. The Lagrangian densities of the models are BRS (Becchi-Rouet-Stora) exact and therefore the models can be regarded as topological theories. In the models, the coupling constants, including the cosmological constant, look as if they run with the scale of the universe and its behavior is very similar to the renormalization group. Motivated by these models, we propose new models with an the infrared fixed point, which may correspond to the late time universe, and an ultraviolet fixed point, which may correspond to the early universe. Especially we construct a model with the solutions corresponding to the de Sitter space-time both in the ultraviolet and the infrared fixed points.

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.

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.

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.

Future singularities arising in a family of models for the expanding Universe, characterized by sharing a convenient parametrization of the energy budget in terms of the deceleration parameter, are classified. Finite-time future singularities are known to appear in many cosmological scenarios, in particular, in the presence of viscosity or non-gravitational interactions, the last being known to be able to suppress or just change in some cases the type of the cosmological singularity. Here, a family of models with a parametrization of the energy budget in terms of the deceleration parameter are studied in the light of Gaussian processes using reconstructed data from $40$-value $H(z)$ datasets. Eventually, the form of the possible non-gravitational interaction between dark energy and dark matter is constructed from these smoothed $H(z)$ data. Using phase space analysis, it is shown that a non-interacting model with dark energy $\omega_\mathrm{de} = \omega_{0} + \omega_{1}q$ ($q$ being the deceleration parameter) may evolve, after starting from a matter dominated unstable state, into a de Sitter Universe (the solution being in fact a stable node). Moreover, for a model with interaction term $Q = 3 H b \rho_\mathrm{dm}$ ($b$ is a parameter and $H$ the Hubble constant) three stable critical points are obtained, what may have important astrophysical implications. In addition, part of the paper is devoted to a general discussion of the finite-time future singularities obtained from direct numerical integration of the field equations, since they appear in many cosmological scenarios and could be useful for future extended studies of the models here introduced. Numerical solutions for the new models, produce finite-time future singularities of Type I or Type III, or an $\omega$-singularity, provided general relativity describes the background dynamics.

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.

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.

In arXiv:1601.02203 and arXiv:1702.07063, we have proposed a topological model with a simple Lagrangian density and have tried to solve one of the cosmological constant problems. The Lagrangian density is the BRS exact and therefore the model can be regarded as a topological theory. In this model, the divergence of the vacuum energy coming from the quantum corrections from matters can be absorbed into the redefinition of the scalar field. In this paper, we consider the extension of the model in order to apply the mechanism to other kinds of divergences coming from the quantum correction and consider the cosmology in an extended model.

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.

We make a theoretical prediction for the ratio of the dark energy to other components in the Universe based on the scenario of the sequestering mechanism which was recently proposed as one possible way to solve the cosmological constant problem. In order to evaluate the value of dark energy and the others, we assume a specific scale factor which describes the Big-Crunch scenario in the scalar-tensor theory. We specify the parameter region where the one can explain the observed dark energy-matter ratio.

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.

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 arXiv:1601.02203, a simple model has been proposed in order to solve one of the problems related with the cosmological constant. The model is given by a topological field theory and the model has an infinite numbers of the BRS symmetries. The BRS symmetries are, however, spontaneously broken in general. In this paper, we investigate the BRS symmetry in more details and show that there is one and only one BRS symmetry which is not broken and the unitarity can be guaranteed. In the model, the quantum problem of the vacuum energy, which may be identified with the cosmological constant, reduces to the classical problem of the initial condition. In this paper, we investigate the cosmology given by the model and specify the region of the initial conditions which could be consistent with the evolution of the universe. We also show that there is a stable solution describing the de Sitter space-time, which may explain the accelerating expansion in the current universe.

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.

In the previous works (arXiv:1202.5375 and 1402.1346), the dynamical domain wall, where the four dimensional FRW universe is embedded in the five imensional space-time, has been realized by using two scalar fields. In this paper, we consider the localization of vector field in three formulations. The first formulation was investigated in the previous paper (arXiv:1510.01099) for the $U(1)$ gauge field. In the second formulation, we investigate the Dvali-Shifman mechanism (hep-th/9612128), where the non-abelian gauge field is confined in the bulk but the gauge symmetry is spontaneously broken on the domain wall. In the third formulation, we investigate the Kaluza-Klein modes coming from the five dimensional graviton. In the Randall-Sundrum model, the graviton was localized on the brane. We show that the $(5,\mu)$ components $\left(\mu=0,1,2,3\right)$ of the graviton are also localized on the domain wall and can be regarded as the vector field on the domain wall. There are, however, some corrections coming from the bulk extra dimension if the domain wall universe is expanding.

We propose several covariant models which may solve one of the problems in the cosmological constant. One of the model can be regarded as an extension of sequestering model. Other models could be regarded as extensions of the covariant formulation of the unimodular gravity. The contributions to the vacuum energy from the quantum corrections from the matters are absorbed into a redefinition of a scalar field and the quantum corrections become irrelevant to the dynamics. In a class of the extended unimodular gravity models, we also consider models which are regarded as topological field theories. The models can be extended and not only the vacuum energy but any quantum corrections to the gravitational action could become irrelevant for the dynamics. We find, however, that the BRS symmetry in the topological field theories is broken spontaneously and therefore the models might not be consistent.

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.

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 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 consider the Palatini formalism of the Born-Infeld gravity. In the Palatini formalism, the propagating mode is only graviton, whose situation is different from that in the metric formalism. We discuss about the FRW cosmology by using an effective potential. Especially we consider the condition that the bouncing could occur. We also give some speculations about the black hole formation

Based on the previous works (arXiv:1202.5375 and 1402.1346), we investigate the localization of the fields on the dynamical domain wall, where the four dimensional FRW universe is realized on the domain wall in the five dimensional space-time. Especially we show that the chiral spinor can localize on the domain wall, which has not been succeeded in the past works as the seminal work in arXiv:0810.3746.

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 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 Phys.\ Rev.\ D 90 043006 (2014), we proposed a new ghost-free massive spin-two model in flat spacetime. Furthermore, as some extension, we couple the new model with a non-dynamical curved background in Phys.\ Rev.\ D 90 123013 (2014) and constructed new interaction terms without appearance of an extra mode. The characteristic property of the new model is the existence of nonlinear potential terms which give the nontrivial vacua. The presence of the nontrivial vacua, however, does not mean that the particle can be defined around all vacua. Therefore, in this paper, we discuss the condition for the new model to have stable vacua in flat spacetime and curved spacetime. Then, we couple this spin-two theory with a dynamical background and obtain the solutions. Moreover, we investigate the effect of this new spin-two model to the Einstein gravity by calculating the black hole entropy since the gravity coupled massive spin-two theory admits a black hole solution in addition to (anti-) de Sitter space solution.

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

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

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.

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

We have proposed a new ghost-free model with interactions of massive spin two particles in Phys.\ Rev.\ D {\bf 90} (2014) 043006 [arXiv:1402.5737 [hep-th]]. Although the model is ghost-free on the Minkowski space-time, it is not obvious whether or not this desirable property is preserved on curved space-time. In fact, Buchbinder et al. already pointed out that the Fierz-Pauli theory is not ghost-free on curved space-time without non-minimal coupling terms. In this paper, we construct a new theory of massive spin two particles with non-minimal coupling on curved space-time and show that the model can be ghost-free. Furthermore, we propose new non-minimal coupling terms.

We propose a new ghost-free model describing massive spin two field. This model consists of a kinetic term and interaction terms without derivative. We report on the properties of this model, especially we consider what could happen when this model couples with gravity. Although the model does not generate any ghost on the Minkowski space-time, it is not so clear whether or not this property is preserved even on curved space-time. In fact, Buchbinder et al. have found that the ghost appears even in the Fierz-Pauli theory on curved space-time if we do not include non-minimal coupling terms. We report on the model with interactions on curved space-time and show that we can construct a model without ghost by including non-minimal coupling terms.

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

In this paper, we propose a toy model of the renormalizable theory describing massive spin two field. Although the model is renormalizable, we show that the model contains ghost. The coupling of the theory with gravity can be regarded as a new kind of bimetric gravity or bigravity. We show that the massive spin two field plays the role of the cosmological constant.

The bigravity models coupled with two scalar fields are constructed. We show that a wide class of the expansion history of the universe, especially corresponding to dark energy and/or inflation, can be described by a solution of the bigravity model. We discuss the stability of the solution and give the conditions for the stability. We also explicitly construct a model which gives a stable solution. By using the stable model, for an arbitrary evolution of the universe expansion, we construct the Brans-Dicke like model which reproduces the evolution.

In this paper we investigate the case of non-minimal coupling in the (extended) nonlinear massive gravity theories. We first consider massive gravity in the Brans-Dicke background such that the graviton mass is replaced by $A^2(\sigma)m$ where $\sigma$ is the Brans-Dicke field and $A(\sigma)$ is conformal coupling and show that there is no viable thermal history of the universe in this case. We then invoke a cubic galileon term as nonlinear completion of the $\sigma$ Lagrangian and show that there is a stable de Sitter solution in this case. However, the de Sitter is blocked by the matter phase which is also a simultaneous attractor of the dynamics. The de Sitter phase can, however, be realized by invoking unnatural fine tunings. We next investigate cosmology of quasi-dilaton nonlinear massive gravity with non-minimal coupling. As a generic feature of the non-minimal coupling, we show that the model exhibits a transient phantom phase which is otherwise impossible. While performing the observational data analysis on the models, we find that a small value of coupling constant is allowed for quasi-dilaton nonlinear massive gravity. For both the cases under consideration, it is observed that we have an effective pressure of matter which comes from the constraint equation. For mass-varying nonlinear massive gravity in the Brans-Dicke background, the effective pressure of matter is non zero which affects the evolution of the Hubble parameter thereby spoiling consistency of the model with data. As for, quasi-dilaton nonlinear massive gravity, the effective pressure of matter can be kept around zero by controlling the coupling constant, the model is shown to be fit well with observations.

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 a four-dimensional domain wall universe by using the Brans-Dicke type gravity with two scalar field. We give a formulation where for arbitrarily given warp factor and scale factor, we construct an action which reproduces both of the warp and scale factors as a solution of the Einstein equation and the field equations given by the action. This formulation could be called as reconstruction. We show that the model does not contain ghost with negative kinetic term and there occur the localization of gravity as in the Randall-Sundrum model. It should be noted that in the equation of the graviton, there appear extra terms related with the extra dimension, which might affect the tensor mode in the fluctuations

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.

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 explore the possibility of further gravitational wave modes in $F(T)$ gravity, where $T$ is the torsion scalar in teleparallelism. It is explicitly demonstrated that gravitational wave modes in $F(T)$ gravity are equivalent to those in General Relativity. This result is achieved by calculating the Minkowskian limit for a class of analytic function of $F(T)$. This consequence is also confirmed by the preservative analysis around the flat background in the weak field limit with the scalar-tensor representation of $F(T)$ gravity.

The Nine-Year WMAP results combined with other cosmological data seem to indicate an enhanced favor for the phantom regime, comparing to previous analyses. This behavior, unless reversed by future observational data, suggests to consider the phantom regime more thoroughly. In this work we provide three modified gravitational scenarios in which we obtain the phantom realization without the appearance of ghosts degrees of freedom, which plague the naive approaches on the subject, namely the Brans-Dicke type gravity, the scalar-Einstein-Gauss-Bonnet gravity, and the $F(R)$ gravity, which are moreover free of perturbative instabilities. The phantom regime seems to favor the gravitational modification instead of the universe-content alteration.

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 symmetron model in a generalized background with a hope to make it compatible with dark energy. We observe a "no go" theorem at least in case of a conformal coupling. Being convinced of symmetron incapability to be dark energy, we try to retain its role for spontaneous symmetry breaking and assign the role of dark energy either to standard quintessence or $F(R)$ theory which are switched on by symmetron field in the symmetry broken phase. The scenario reduces to standard Einstein gravity in the high density region. After the phase transition generated by symmetron field, either the $F(R)$ gravity or the standard quintessence are induced in the low density region. we demonstrate that local gravity constraints and other requirements are satisfied although the model could generate the late-time acceleration of Universe.

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.

We develop a formulation where for an arbitrarily given warp factor, we construct a scalar field action which reproduces the warp factor as a solution of the Einstein equation and field equation corresponding to the action. This formulation could be called as the reconstruction. By using the formulation of the reconstruction, we construct models which have an exact solution describing the domain wall. The shape of the domain wall can be flat, de Sitter space-time, or anti-de Sitter space-time. In the constructed domain wall solutions, there often appears ghost with negative kinetic energy. We give, however, an example of the de Sitter domain wall solution without ghost, which could be a toy model of inflation. We also investigate the localization of gravity as in the Randall-Sundrum model. It is demonstrated that the four dimensional Newton law could be reproduced even in the de Sitter space-time domain wall solution. We show that we can construct a space-time, where the domain wall is the general Friedmann-Robertson-Walker (FRW) universe and the warp factor can be arbitrary. For such a construction, we use two scalar fields. It is also illustrated that the scalar field equations are equivalent to the Bianchi identities: $\nabla^\mu (R_{\mu\nu} - \frac{1}{2}R g_{\mu\nu})=0$.

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.

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.

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.

Explicit formulae of the equations in the generalized Galileon models are given. We also develop the formulation of the reconstruction. By using the formulation, we can explicitly construct an action which reproduces an arbitrary development of the expansion of the universe. The conditions how the reconstructed solution becomes stable and therefore it becomes an attractor solution are also given. Working in the static and spherically symmetric space-time, we investigate how the Vainshtein mechanism works in the generalized Galileon model and the correction to the Newton law becomes small. It is also shown that any spherically symmetric and static geometry can be realized by properly choosing the form of the action, which may tell that the solution could have fourth hair corresponding to the scalar field.

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.

Subhorizon approximation is often used in cosmological perturbation theory. In this paper, however, it is shown that the subhorizon approximation is not always a good approximation at least in case of $k$-essence model. We also show that the sound speed given by $k$-essence model exerts a huge influence on the time evolution of the matter density perturbation, and the future observations could clarify the differences between the $\Lambda$CDM model and $k$-essence model.

We construct cosmological models with two scalar fields, which has the structure as in the ghost condensation model or k-essence model. The models can describe the stable phantom crossing, which should be contrasted with one scalar tensor models, where the infinite instability occurs at the crossing the phantom divide. We give a general formulation of the reconstruction in terms of the e-foldings N by including the matter although in the previous two scalar models, which are extensions of the scalar tensor model, it was difficult to give a formulation of the reconstruction when we include matters. In the formulation of the reconstruction, we start with a model with some arbitrary functions, and find the functions which generates the history in the expansion of the universe. We also give general arguments for the stabilities of the models and the reconstructed solution. The viability of a model is also investigated by comparing the observational data.

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.

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.

We review on the models of gravity with a constraint by the Lagrange multiplier field. The constraint breaks general covariance or Lorentz symmetry in the ultraviolet region. We report on the $F(R)$ gravity model with the constraint and the proposal of the covariant (power-counting) renormalized gravity model by using the constraint and scalar projectors. We will show that the model admits flat space solution, its gauge-fixing formulation is fully developed, and the only propagating mode is (higher derivative) graviton, while scalar and vector modes do not propagate. The preliminary study of FRW cosmology indicates to the possibility of inflationary universe solution is also given.

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.

By using the formulation of the reconstruction, we explicitly construct models of $k$-essence, which unify the inflation in the early universe and the late accelerating expansion of the present universe by a single scalar field. Due to the higher derivative terms, the solution describing the unification can be stable in the space of solutions, which makes the restriction for the initial condition relaxed. The higher derivative terms also eliminate tachyon. Therefore we can construct a model describing the time development, which cannot be realized by a usual inflaton or quintessence models of the canonical scalar field due to the instability or the existence of tachyon. We also propose a mechanism of the reheating by the quantum effects coming from the variation of the energy density of the scalar field.

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.

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.

Based on the formulation of the reconstruction for the k-essence model, which was recently proposed in arXiv:1001.0220 [hep-th], we explicitly construct cosmological model to unifying the late-time acceleration and the inflation in the early universe.

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.

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 review how we can construct the gravity models which reproduces the arbitrary development of the universe. We consider the reconstruction in the Einstein gravity coupled with generalized perfect fluid, scalar-Einstein gravity, scalar-Einstein-Gauss-Bonnet gravity, Einstein-$F(G)$-gravity, and $F(R)$-gravity. Very explicit formulas are given to reconstruct the models, which could be used when we find the detailed data of the development of the universe by future observations. Especially we find the formulas using e-foldings, which has a direct relation with observed redshift. As long as we observe the time development of the Hubble rate $H$, there exists a variety of models describing the arbitrary development of universe.

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.

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.

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 review cosmology in non-minimal Yang-Mills/Maxwell theory, in which the Yang-Mills/electromagnetic field couples to a function of the scalar curvature. We show that power-law inflation can be realized due to the non-minimal gravitational coupling of the Yang-Mills field which may be caused by quantum corrections. Moreover, we study non-minimal vector model in the framework of modified gravity and demonstrate that both inflation and the late-time accelerated expansion of the universe can be realized. We also discuss the cosmological reconstruction of the Yang-Mills theory. Furthermore, we investigate late-time cosmology in non-minimal Maxwell-Einstein theory. We explore the forms of the non-minimal gravitational coupling which generate the finite-time future singularities and the general conditions for this coupling in order that the finite-time future singularities cannot appear.

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 study the $f(R)$ theory of gravity using metric approach. In particular we investigate the recently proposed model by Hu-Sawicki, Appleby $-$ Battye and Starobinsky. In this model, the cosmological constant is zero in flat space time. The model passes both the Solar system and the laboratory tests. But the model parameters need to be fine tuned to avoid the finite time singularity recently pointed in the literature. We check the concordance of this model with the $H(z)$ and baryon acoustic oscillation data. We find that the model resembles the $\Lambda$CDM at high redshift. However, for some parameter values there are variations in the expansion history of the universe at low redshift.

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 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 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 review the several models of the dark energy, which may generate the accelerated expansion of the present universe. We also discuss the the Big Rip singularity, which may occur when the equation of the state parameter w is less than -1. We show that the quantum correction would be very important near the singularity.

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.