We consider FRW cosmology in $f(R)= R+ \gamma R^2+\delta R^3$ modified framework. The Palatini approach reduces its dynamics to the simple generalization of Friedmann equation. Thus we study the dynamics in two-dimensional phase space with some details. After reformulation of the model in the Einstein frame, it reduces to the FRW cosmological model with a homogeneous scalar field and vanishing kinetic energy term. This potential determines the running cosmological constant term as a function of the Ricci scalar. As a result we obtain the emergent dark energy parametrization from the covariant theory. We study also singularities of the model and demonstrate that in the Einstein frame some undesirable singularities disappear.

We classify singularities in FRW cosmologies, which dynamics can be reduced to the dynamical system of the Newtonian type. This classification is performed in terms of geometry of a potential function if it has poles. At the sewn singularity, which is of a finite scale factor type, the singularity in the past meets the singularity in the future. We show, that such singularities appear in the Starobinsky model in $f(\hat{R})=\hat{R}+\gamma \hat{R}^2$ in the Palatini formalism, when dynamics is determined by the corresponding piece-wise smooth dynamical system. As an effect we obtain a degenerated singularity. Analytical calculations are given for the cosmological model with matter and the cosmological constant. The dynamics of model is also studied using dynamical system methods. From the phase portraits we find generic evolutionary scenarios of the evolution of the Universe. For this model, the best fit value of $\Omega_\gamma=3\gamma H_0^2$ is equal $9.70\times 10^{-11}$. We consider model in both Jordan and Einstein frames. We show that after transition to the Einstein frame we obtain both form of the potential of the scalar field and the decaying Lambda term.

We investigate further (cf. arXiv:1512.01199, JCAP01 (2016) 040) Starobinsky cosmological model $R+\gamma R^2$ in the Palatini formalism with Chaplygin gas and baryonic matter as a source. For this aim we use dynamical system theory. The dynamics is reduced to the 2D sewn dynamical system of a Newtonian type (a piecewise-smooth dynamical system). We classify all evolutional paths in the model as well as trajectories in the phase space. We demonstrate the presence of a degenerate freeze singularity (glued freeze type singularities) for the positive $\gamma$. In this case it is a generic feature of early evolution of the universe. We point out that a degenerate type III of singularity can be considered as an endogenous model of inflation between the matter dominating epoch and the dark energy phase. We also investigate cosmological models with negative $\gamma$. It is demonstrated that $\gamma$ equal zero is a bifurcation parameter and dynamics qualitatively changes in comparison to positive $\gamma$. Instead of the big bang the sudden bounce singularity of a finite scale factor appears and there is a generic class of bouncing solutions sewn along the line $a=a_{\text{sing}}$. And we argue that the presence of sudden singularities in an evolutional scenario of the Universe falsifies the negative $\gamma$ in the Palatini cosmology. Only very small values of $\Omega_{\gamma}$ parameter are admissible if we requires that agreements physics with the $\Lambda$CDM model. From the statistical analysis of astronomical observations, we deduce that the case of negative values of $\Omega_\gamma$ can be rejected even if it may fit better to the data.

We present a simple generalisation of the $\Lambda$CDM model which on the one hand reaches very good agreement with the present day experimental data and provides an internal inflationary mechanism on the other hand. It is based on Palatini modified gravity with quadratic Starobinsky term and generalized Chaplygin gas as a matter source providing, besides a current accelerated expansion, the epoch of endogenous inflation driven by type III freeze singularity. It follows from our statistical analysis that astronomical data favors negative value of the parameter coupling quadratic term into Einstein-Hilbert Lagrangian and as a consequence the bounce instead of initial Big-Bang singularity is preferred.

Symmetries play a crucial role in physics and, in particular, the Noether symmetries are a useful tool both to select models motivated at a fundamental level, and to find exact solutions for specific Lagrangians. In this work, we consider the application of point symmetries in the recently proposed metric-Palatini Hybrid Gravity in order to select the $f({\cal R})$ functional form and to find analytical solutions for the field equations and for the related Wheeler-DeWitt (WDW) equation. We show that, in order to find out integrable $f({\cal R})$ models, conformal transformations in the Lagrangians are extremely useful. In this context, we explore two conformal transformations of the forms $d\tau=N(a) dt$ and $d\tau=N(\phi) dt$. For the former conformal transformation, we found two cases of $f({\cal R})$ functions where the field equations admit Noether symmetries. In the second case, the Lagrangian reduces to a Brans-Dicke-like theory with a general coupling function. For each case, it is possible to transform the field equations by using normal coordinates to simplify the dynamical system and to obtain exact solutions. Furthermore, we perform quantization and derive the WDW equation for the minisuperspace model. The Lie point symmetries for the WDW equation are determined and used to find invariant solutions.

We study new FRW type cosmological models of modified gravity treated on the background of Palatini approach. These models are generalization of Einstein gravity by the presence of a scalar field non-minimally coupled to the curvature. The models employ Starobinsky's term in the Lagrangian and dust matter. Therefore, as a by-product, an exhausted cosmological analysis of general relativity amended by quadratic term is presented. We investigate dynamics of our models, confront them with the currently available astrophysical data as well as against LCDM model. We have used the dynamical system methods in order to investigate dynamics of the models. It reveals the presence of a final sudden singularity. Fitting free parameters we have demonstrated by statistical analysis that this class of models is in a very good agreement with the data (including CMB measurements) as well as with the standard LCDM model predictions. One has to use statefinder diagnostic in order to discriminate among them. Therefore Bayesian methods of model selection have been employed in order to indicate preferred model. Only in the light of CMB data the concordance model remains invincible.

We explain the effect of dark matter (flat rotation curve) using modified gravitational dynamics. We investigate in this context a low energy limit of generalized general relativity with a nonlinear Lagrangian ${\cal L}\propto R^n$, where $R$ is the (generalized) Ricci scalar and $n$ is parameter estimated from SNIa data. We estimate parameter $\beta$ in modified gravitational potential $V(r) \propto -\frac{1}{r}(1+(\frac{r}{r_c})^{\beta})$. Then we compare value of $\beta$ obtained from SNIa data with $\beta$ parameter evaluated from the best fitted rotation curve. We find $\beta \simeq 0.7$ which becomes in good agreement with an observation of spiral galaxies rotation curve. We also find preferred value of $\Omega_{m,0}$ from the combined analysis of supernovae data and baryon oscillation peak. We argue that although amount of "dark energy" (of non-substantial origin) is consistent with SNIa data and flat curves of spiral galaxies are reproduces in the framework of modified Einstein's equation we still need substantial dark matter. For comparison predictions of the model with predictions of the $\Lambda$CDM concordance model we apply the Akaike and Bayesian information criteria of model selection.

Recent supernovae of type Ia measurements and other astronomical observations suggest that our universe is in accelerating phase of evolution at the present epoch. While a dark energy of unknown form is usually proposed as the most feasible mechanism for the acceleration, there are appears some alternative conception that some effects arising from generalization of Einstein equation can mimic dark energy through a modified Friedmann equation. In this work we investigate some observational constraints on modified Friedmann equation obtained from generalized Lagrangian ${\cal L} \propto R^n$ in minimal coupling with matter in Palatini formalism. We mainly concentrate on the constraints of model parameters from distant supernovae but other constraint from baryon oscillation prior is also considered. We obtain the confidence levels on two additional model parameter ($n$,$\Omega_{m,0}$). We conclude that the FRW model of First-Order Non-linear gravity survives several observational test like SNIa observation and baryon oscillation peaks. We find preferred value of $\Omega_{m,0} \simeq 0.3$ from combined analysis of supernovae data and baryon oscillation peak. For deeper statistical analysis we apply Akaike and Bayesian information criteria of model selection for comparison prediction of the model with prediction of concordance $\Lambda$CDM model.

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

The evidence of the acceleration of universe at present time has lead to investigate modified theories of gravity and alternative theories of gravity, which are able to explain acceleration from a theoretical viewpoint without the need of introducing dark energy. In this paper we study alternative gravitational theories defined by Lagrangians which depend on general functions of the Ricci scalar invariant in minimal interaction with matter, in view of their possible cosmological applications. Structural equations for the spacetimes described by such theories are solved and the corresponding field equations are investigated in the Palatini formalism, which prevents instability problems. Particular examples of these theories are also shown to provide, under suitable hypotheses, a coherent theoretical explanation of earlier results concerning the present acceleration of the universe and cosmological inflation. We suggest moreover a new possible Lagrangian, depending on the inverse of sinh(R), which gives an explanation to the present acceleration of the universe.