The cosmological observations suggest that the presently accelerating universe should be filled by an exotic form of matter, violating the strong energy condition, of unknown nature and origin. We propose the viscous dark matter of a source of acceleration in the form of Chaplygin gas which is characterized by equation of state in the phenomenological form $p=-\frac{A}{\rho^{\alpha}}$, where $p$ and $\rho$ are pressure and energy density respectively ($A$ and $\alpha$ are constants). Chaplygin gas is interpreted in terms of viscous matter and without the cosmological constant. The acceleration effect is caused only by viscosity in this class of cosmological models. We show that bulk viscosity effects introduced to the standard FRW cosmology give rise to the natural unification of both dark matter and dark energy. We show that dust viscous cosmological models are structurally stable if $m < 1/2$ ($1+\alpha=1/2-m$).

In this paper, we apply a method of reducing the dynamics of FRW cosmological models with the barotropic form of the equation of state to the dynamical system of the Newtonian type to detect the finite scale factor singularities and the finite-time singularities. In this approach all information concerning the dynamics of the system is contained in a diagram of the potential function $V(a)$ of the scale factor. Singularities of the finite scale factor manifest by poles of the potential function. In our approach the different types of singularities are represented by critical exponents in the power-law approximation of the potential. The classification can be given in terms of these exponents. We have found that the pole singularity can mimick an inflation epoch. We demonstrate that the cosmological singularities can be investigated in terms of the critical exponents of the potential function of the cosmological dynamical systems. We assume the general form of the model contains matter and some kind of dark energy which is parameterized by the potential. We distinguish singularities (by ansatz about the Lagrangian) of the pole type with the inflation and demonstrate that such a singularity can appear in the past.

We consider a cosmology with decaying metastable dark energy and assume that a decay process of this metastable dark energy is a quantum decay process. Such an assumption implies among others that the evolution of the Universe is irreversible and violates the time reversal symmetry. We show that if to replace the cosmological time $t$ appearing in the equation describing the evolution of the Universe by the Hubble cosmological scale time, then we obtain time dependent $\Lambda (t)$ in the form of the series of even powers of the Hubble parameter $H$: $\Lambda (t) = \Lambda (H)$. Out special attention is focused on radioactive like exponential form of the decay process of the dark energy and on the consequences of this type decay.

We investigate cosmological dynamics based on $f(R)$ gravity in the Palatini formulation. In this study we use the dynamical system methods. We show that the evolution of the Friedmann equation reduces to the form of the piece-wise smooth dynamical system. This system is is reduced to a 2D dynamical system of the Newtonian type. We demonstrate how the trajectories can be sewn to guarantee $C^0$ extendibility of the metric similarly as `Milne-like' FLRW spacetimes are $C^0$-extendible. We point out that importance of dynamical system of Newtonian type with non-smooth right-hand sides in the context of Palatini cosmology. In this framework we can investigate singularities which appear in the past and future of the cosmic evolution. We consider cosmological systems in both Einstein and Jordan frames. We show that at each frame the topological structures of phase space are different.

In the framework of polynomial Palatini cosmology, we investigate a simple cosmological homogeneous and isotropic model with matter in the Einstein frame. We show that in this model during cosmic evolution, it appears the early inflation and the accelerating phase of the expansion for the late times. In this frame we obtain the Friedmann equation with matter and dark energy in the form of a scalar field with the potential whose form is determined in a covariant way by the Ricci scalar of the FRW metric. The energy density of matter and dark energy are also parametrized through the Ricci scalar. The early inflation is obtained only for an infinitesimally small fraction of energy density of matter. Between the matter and dark energy, there exists interaction because the dark energy is decaying. For characterization of inflation we calculate the slow roll parameters and the constant roll parameter in terms of the Ricci scalar. We have found a characteristic behaviour of the time dependence of density of dark energy on the cosmic time following the logistic-like curve which interpolates two almost constant value phases. From the required numbers of $N$-folds we have found a bound on model parameter.

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 study the cosmology with the running dark energy. The parametrization of dark energy with the respect to the redshift is derived from the first principles of quantum mechanics. Energy density of dark energy is obtained from the quantum process of transition from the false vacuum state to the true vacuum state. This is the class of the extended interacting $\Lambda$CDM models. We consider the energy density of dark energy parametrization $\rho_\text{de}(t)$, which follows from the Breit-Wigner energy distribution function which is used to model the quantum unstable systems. The idea that properties of the process of the quantum mechanical decay of unstable states can help to understand the properties of the observed universe was formulated by Krauss and Dent and this idea was used in our considerations. In the cosmological model with the mentioned parametrization there is an energy transfer between the dark matter and dark energy. In such a evolutional scenario the universe is starting from the false vacuum state and going to the true vacuum state of the present day universe. We find that the intermediate regime during the passage from false to true vacuum states takes place. The intensity of the analyzed process is measured by a parameter $\alpha$. For the small value of $\alpha$ ($0<\alpha <0.4$) this intermediate (quantum) regime is characterized by an oscillatory behavior of the density of dark energy while the for $\alpha > 0.4$ the density of the dark energy simply jumps down. In both cases (independent from the parameter $\alpha$) the today value of density of dark energy is reached at the value of $0.7$. We estimate the cosmological parameters for this model with visible and dark matter. This model becomes in good agreement with the astronomical data and is practically indistinguishable from $\Lambda$CDM model.

We study dynamics of cosmological models with diffusion effects modeling dark matter and dark energy interactions. We show the simple model with diffusion between the cosmological constant sector and dark matter, where the canonical scaling law of dark matter $(\rho_{dm,0}a^{-3}(t))$ is modified by an additive $\epsilon(t)=\gamma t a^{-3}(t)$ to the form $\rho_{dm}=\rho_{dm,0}a^{-3}(t)+\epsilon(t)$. We reduced this model to the autonomous dynamical system and investigate it using dynamical system methods. This system possesses a two-dimensional invariant submanifold on which the DM-DE interaction can be analyzed on the phase plane. The state variables are density parameter for matter (dark and visible) and parameter $\delta$ characterizing the rate of growth of energy transfer between the dark sectors. A corresponding dynamical system belongs to a general class of jungle type of cosmologies represented by coupled cosmological models in a Lotka-Volterra framework. We demonstrate that the de Sitter solution is a global attractor for all trajectories in the phase space and there are two repellers: the Einstein-de Sitter universe and the de Sitter universe state dominating by the diffusion effects. We distinguish in the phase space trajectories, which become in good agreement with the data. They should intersect a rectangle with sides of $\Omega_{m,0}\in [0.2724, 0.3624]$, $\delta \in [0.0000, 0.0364]$ at the 95\% CL. Our model could solve some of the puzzles of the $\Lambda$CDM model, such as the coincidence and fine-tuning problems. In the context of the coincidence problem, our model can explain the present ratio of $\rho_{m}$ to $\rho_{de}$, which is equal $0.4576^{+0.1109}_{-0.0831}$ at a 2$\sigma$ confidence level.

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.

We discussed the dynamics of cosmological models in which the cosmological constant term is a time dependent function through the scale factor $a(t)$, Hubble function $H(t)$, Ricci scalar $R(t)$ and scalar field $\phi(t)$. We considered five classes of models; two non-covariant parametrization of $\Lambda$: 1) $\Lambda(H)$CDM cosmologies where $H(t)$ is the Hubble parameter, 2) $\Lambda(a)$CDM cosmologies where $a(t)$ is the scale factor, and three covariant parametrization of $\Lambda$: 3) $\Lambda(R)$CDM cosmologies, where $R(t)$ is the Ricci scalar, 4) $\Lambda(\phi)$-cosmologies with diffusion, 5) $\Lambda(X)$-cosmologies, where $X=\frac{1}{2}g^{\alpha\beta}\nabla_{\alpha}\nabla_{\beta}\phi$ is a kinetic part of density of the scalar field. We also considered the case of an emergent $\Lambda(a)$ relation obtained from the behavior of trajectories in a neighborhood of an invariant submanifold. In study of dynamics we use dynamical system methods for investigating how a evolutional scenario can depend on the choice of special initial conditions. We showed that methods of dynamical systems offer the possibility of investigation all admissible solutions of a running $\Lambda$ cosmology for all initial conditions, their stability, asymptotic states as well as a nature of the evolution in the early universe (singularity or bounce) and a long term behavior at the large times. We also formulated an idea of the emergent cosmological term derived directly from an approximation of exact dynamics. We show that some non-covariant parametrizations of Lambda term like $\Lambda(a)$, $\Lambda(H)$ give rise to pathological and nonphysical behaviour of trajectories in the phase space. This behaviour disappears if the term $\Lambda(a)$ is emergent from the covariant parametrization.

Dynamical systems methods are used to investigate global behavior of the spatially flat Friedmann-Robertson-Walker cosmological model in gravitational theory with a non-minimally coupled scalar field and a constant potential function. We show that the system can be reduced to an autonomous three-dimensional dynamical system and additionally is equipped with an invariant manifold corresponding to an accelerated expansion of the universe. Using this invariant manifold we find an exact solution of the reduced dynamics. We investigate all solutions for all admissible initial conditions using theory of dynamical systems to obtain a classification of all evolutional paths. The right-hand sides of the dynamical system depend crucially on the value of the non-minimal coupling constant therefore we study bifurcation values of this parameter under which the structure of the phase space changes qualitatively. We found a special bifurcation value of the non-minimal coupling constant which is distinguished by dynamics of the model and may suggest some additional symmetry in matter sector of the theory.

We study dynamics of $\Lambda(t)$ cosmological models which are a natural generalization of the standard cosmological model (the $\Lambda$CDM model). We consider a class of models: the ones with a prescribed form of $\Lambda(t)=\Lambda_{\text{bare}}+\frac{\alpha^2}{t^2}$. This type of a $\Lambda(t)$ parametrization is motivated by different cosmological approaches. We interpret the model with running Lambda ($\Lambda(t)$) as a special model of an interacting cosmology with the interaction term $-d\Lambda(t)/dt$ in which energy transfer is between dark matter and dark energy sectors. For the $\Lambda(t)$ cosmology with a prescribed form of $\Lambda(t)$ we have found the exact solution in the form of Bessel functions. Our model shows that fractional density of dark energy $\Omega_e$ is constant and close to zero during the early evolution of the universe. We have also constrained the model parameters for this class of models using the astronomical data such as SNIa data, BAO, CMB, measurements of $H(z)$ and the Alcock-Paczy{\'n}ski test. In this context we formulate a simple criterion of variability of $\Lambda$ with respect to $t$ in terms of variability of the jerk or sign of estimator $(1-\Omega_{\text{m},0}-\Omega_{\Lambda,0})$. The case study of our model enable us to find an upper limit $\alpha^2 < 0.012$ ($2\sigma$ C.L.) describing the variation from the cosmological constant while the LCDM model seems to be consistent with various data.

We investigate the dynamics of the generalized $\Lambda$CDM model, which the $\Lambda$ term is running with the cosmological time. The $\Lambda(t)$ term emerges from the covariant theory of the scalar field $\phi$ with the self-interacting potential $V(\phi)$. On the example of the model $\Lambda(t)=\Lambda_{\text{bare}}+\frac{\alpha^2}{t^2}$ we show the existence of a mechanism of the modification of the scaling law for energy density of dark matter: $\rho_{\text{dm}}\propto a^{-3+\lambda(t)}$. We discuss the evolution of $\Lambda(t)$ term and pointed out that during the cosmic evolution there is a long phase in which this term is approximately constant. This effect justifies Alcaniz and Lima's approach to $\Lambda(H)$ cosmologies. We also present the statistical analysis of both the $\Lambda(t)$CDM model with dark energy and decaying dark matter and the $\Lambda$CDM standard cosmological model. We divide the observational data into two groups: low $z$ data (SNIa, BAO, $H(z)$ and AP test) and high $z$ data (Planck, WP and lensing). While for the former we find the best fit value of the parameter $\lambda$ is positive ($\lambda=0.0338$, energy transfer is from the dark energy to dark matter sector), for the latter we find that $\lambda$ is $-0.0199$ which is an evidence that the energy transfer is from decaying dark matter. This disagreement correlates with estimated values of $H_0$ ($67.77$ km/(s Mpc) and $65.62$ km/(s Mpc) respectively). The decaying dark matter causes to lowering a mass of dark matter particles which are lighter than CDM particles and remain relativistic. The rate of the process of decaying matter is estimated. We show that in the models of decaying dark matter, the cosmological constant problem disappears naturally. The model with decaying dark matter possesses one parameter more but in light of the AIC it is better than the $\Lambda$CDM standard cosmological model.

Within the quantum mechanical treatment of the decay problem one finds that at late times $t$ the survival probability of an unstable state cannot have the form of an exponentially decreasing function of time $t$ but it has an inverse power-like form. This is a general property of unstable states following from basic principles of quantum theory. The consequence of this property is that in the case of false vacuum states the cosmological constant becomes dependent on time: $\Lambda - \Lambda_{\text{bare}}\equiv \Lambda(t) -\Lambda_{\text{bare}} \sim 1/t^{2}$. We construct the cosmological model with decaying vacuum energy density and matter for solving the cosmological constant problem and the coincidence problem. We show the equivalence of the proposed decaying false vacuum cosmology with the $\Lambda(t)$ cosmologies (the $\Lambda(t)$CDM models). The cosmological implications of the model of decaying vacuum energy (dark energy) are discussed. We constrain the parameters of the model with decaying vacuum using astronomical data. For this aim we use the observation of distant supernovae of type Ia, measurements of $H(z)$, BAO, CMB and others. The model analyzed is in good agreement with observation data and explain a small value of the cosmological constant today.

We construct the cosmological model to explain the cosmological constant problem. We built the extension of the standard cosmological model $\Lambda$CDM by consideration of decaying vacuum energy represented by the running cosmological term. From the principles of quantum mechanics one can find that in the long term behavior survival probability of unstable states is a decreasing function of the cosmological time and has the inverse power-like form. This implies that cosmological constant $\rho_{\text{vac}} = \Lambda(t) = \Lambda_{\text{bare}} + \frac{\alpha}{t^2}$ where $\Lambda_{\text{bare}}$ and $\alpha$ are constants. We investigate the dynamics of this model using dynamical system methods due to a link to the $\Lambda(H)$ cosmologies. We have found the exact solution for the scale factor as well as the indicators of its variability like the deceleration parameter and the jerk. From the calculation of the jerk we obtain a simple test of the decaying vacuum in the FRW universe. Using astronomical data (SNIa, $H(z)$, CMB, BAO) we have estimated the model parameters and compared this model with the $\Lambda$CDM model. Our statistical results indicate that the decaying vacuum model is a little worse than the $\Lambda$CDM model. But the decaying vacuum cosmological model explains the small value of the cosmological constant today.

Recent astronomical observations have indicated that the Universe is in the phase of accelerated expansion. While there are many cosmological models which try to explain this phenomenon, we focus on the interacting $\Lambda$CDM model where the interaction between the dark energy and dark matter sectors takes place. This model is compared to its simpler alternative---the $\Lambda$CDM model. To choose between these models the likelihood ratio test was applied as well as the model comparison methods (employing Occam's principle): the Akaike information criterion (AIC), the Bayesian information criterion (BIC) and the Bayesian evidence. Using the current astronomical data: SNIa (Union2.1), $h(z)$, BAO, Alcock--Paczynski test and CMB we evaluated both models. The analyses based on the AIC indicated that there is less support for the interacting $\Lambda$CDM model when compared to the $\Lambda$CDM model, while those based on the BIC indicated that there is the strong evidence against it in favor the $\Lambda$CDM model. Given the weak or almost none support for the interacting $\Lambda$CDM model and bearing in mind Occam's razor we are inclined to reject this model.

We investigate observational constraints on the Brans-Dicke cosmological model using observational data coming from distant supernovae type Ia, the Hubble function $H(z)$ measurements, information coming from the Alcock-Paczy{\'n}ski test, and baryon acoustic oscillations. Our analysis is based on the modified Friedmann function resulting form dynamical investigations of Brans-Dicke cosmology in the vicinity of a de Sitter state. The qualitative theory of dynamical systems enables us to obtain three different behaviors in the vicinity of this state. We find for a linear approach to the de Sitter state $\omega_{\textrm{BD}}=-0.8606^{+0.8281}_{-0.1341}$, for an oscillatory approach to the de Sitter state $\omega_{\textrm{BD}}=-1.1103^{+0.1872}_{-0.1729}$, and for the transient de Sitter state represented by a saddle-type critical point $\omega_{\textrm{BD}}=-2.3837^{+0.4588}_{-4.5459}$. We obtain the mass of the Brans-Dicke scalar field at the present epoch as $m_{\phi}\sim H_{0}$. The Bayesian methods of model comparison are used to discriminate between obtained models. We show that observational data point toward vales of the $\omega_{\textrm{BD}}$ parameter close to the value suggested by the low-energy limit of the bosonic string theory.

The dynamics of the Brans-Dicke theory with a quadratic scalar field potential function and barotropic matter is investigated. The dynamical system methods are used to reveal complexity of dynamical evolution in homogeneous and isotropic cosmological models. The structure of phase space crucially depends on the parameter of the theory $\omega_{\textrm{BD}}$ as well as barotropic matter index $w_{m}$. In our analysis these parameters are treated as bifurcation parameters. We found sets of values of these parameters which lead to generic evolutional scenarios. We show that in isotropic and homogeneous models in the Brans-Dicke theory with a quadratic potential function the de Sitter state appears naturally. Stability conditions of this state are fully investigated. It is shown that these models can explain accelerated expansion of the Universe without the assumption of the substantial form of dark matter and dark energy. The Poincare construction of compactified phase space with a circle at infinity is used to show that phase space trajectories in a physical region can be equipped with a structure of a vector field on nontrivial topological closed space. For $\omega_{\textrm{BD}}<-3/2$ we show new types of early and late time evolution leading from the anti-de Sitter to the de Sitter state through an asymmetric bounce. In the theory without a ghost we find bouncing solutions and the coexistence of the bounces and the singularity. Following the Peixoto theorem some conclusions about structural stability are drawn.

The dynamics of the Brans-Dicke theory with a scalar field potential function is investigated. We show that the system with a barotropic matter content can be reduced to an autonomous three-dimensional dynamical system. For an arbitrary potential function we found the values of the Brans-Dicke parameter for which a global attractor in the phase space representing de Sitter state exists. Using linearized solutions in the vicinity of this critical point we show that the evolution of the Universe mimics the $\Lambda$CDM model. From the recent Planck satellite data, we obtain constraints on the variability of the effective gravitational coupling constant as well as the lower limit of the mass of the Brans-Dicke scalar field at the de Sitter state.

We study the Scalar Field Cosmology (SFC) using the geometric language of the phase space. We define and study an ensemble of dynamical systems as a Banach space with a Sobolev metric. The metric in the ensemble is used to measure a distance between different models. We point out the advantages of visualisation of dynamics in the phase space. It is investigated the genericity of some class of models in the context of fine tuning of the form of the potential function in the ensemble of SFC. We also study the symmetries of dynamical systems of SFC by searching for their exact solutions. In this context we stressed the importance of scaling solutions. It is demonstrated that scaling solutions in the phase space are represented by unstable separatrices of the saddle points. Only critical point itself located on two dimensional stable submanifold can be identified as scaling solution. We have also found a class of potentials of the scalar fields forced by the symmetry of differential equation describing the evolution of the universe. A class of potentials forced by scaling (homology) symmetries was given. We point out the role of the notion of a structural stability in the context of the problem of indetermination of the potential form of the Scalar Field Cosmology. We characterise also the class of potentials which reproduces the \Lambda CDM model, which is known to be structurally stable. We show that the structural stability issue can be effectively used is selection of the scalar field potential function. This enables us to characterise a structurally stable and therefore a generic class of SFC models. We have found a nonempty and dense subset of structurally stable models. We show that these models possess symmetry of homology.

Properties of unstable false vacuum states are analyzed from the point of view of the quantum theory of unstable states. Some of false vacuum states survive up to times when their survival probability has a non-exponential form. At times much latter than the transition time, when contributions to the survival probability of its exponential and non-exponential parts are comparable, the survival probability as a function of time t has an inverse power-like form. We show that at this time region the instantaneous energy of the false vacuum states tends to the energy of the true vacuum state as $1/t^{2}$ for $t \to \infty$.

Dynamical systems methods are used to study the evolution of the Bianchi I model with a scalar field. We show that inclusion of the non-minimal coupling term between the scalar field and the curvature changes evolution of the model compared with the minimally coupled case. In the model with the non-minimally coupled scalar field there is a new type of singularity dominated by the non-minimal coupling term. We examine the impact of the non-minimal coupling on the anisotropy evolution and demonstrate the existence of its minimal value in the generic case.

We study evolution of cosmological models filled with the scalar field and barotropic matter. We consider the scalar field minimally and non-minimally coupled to gravity. We demonstrated the growth of degree of complexity of evolutional scenario through the description of matter content in terms of the scalar field. In study of all evolutional paths for all initial conditions methods of dynamical systems are used. Using linearized solutions we present simple method of derivation corresponding form of the Hubble function of the scale factor $H(a)$.

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 investigate dynamics of a flat FRW cosmological model with a barotropic matter and a non-minimally coupled scalar field (both canonical and phantom). In our approach we do not assume any specific form of a potential function for the scalar field and we are looking for generic scenarios of evolution. We show that dynamics of universe can be reduced to a 3-dimensional dynamical system. We have found the set of fixed points and established their character. These critical points represent all important epochs in evolution of the universe : (a) a finite scale factor singularity, (b) an inflation (rapid-roll and slow-roll), (c) a radiation domination, (d) a matter domination and (e) a quintessence era. We have shown that the inflation, the radiation and matter domination epochs are transient ones and last for a finite amount of time. The existence of the radiation domination epoch is purely the effect of a non-minimal coupling constant. We show the existence of a twister type solution wandering between all these critical points.

We study the dynamical evolution of cosmological models with the Robertson-Walker symmetry with a scalar field non-minimally coupled to gravity and barotropic matter. For this aim we use dynamical system methods. We have found a type of evolutional path which links between all important events during the evolution, the cosmological singularity of finite time, inflation, radiation and matter dominating epoch and the accelerated phase expansion of the universe. We point out importance of finding the new generic solution called a twister solution for a deeper description of the evolution of the Universe. We demonstrate that including the non-minimal coupling leads to a new, richer evolutional cosmological scenario in comparison to the case of minimal coupling.

We study generic solutions in a non-minimally coupled to gravity scalar field cosmology. It is shown that dynamics for both canonical and phantoms scalar fields with the potential can be reduced to the dynamical system from which the exact forms for an equation of the state parameter can be derived. We have found the stationary solutions of the system and discussed their stability. Within the large class of admissible solutions we have found a non-degenerate critical points and we pointed out multiple attractor type of trajectory travelling in neighborhood of three critical points at which we have the radiation dominating universe, the barotropic matter dominating state and finally the de Sitter attractor. We have demonstrated the stability of this trajectory which we call the twister solution. Discovered evolutional path is only realized if there exist the non-minimal coupling constant. We have found simple duality relations between twister solutions in phantom and canonical scalar fields in the radiation domination phase. For the twister trajectory we have found an oscillating regime of approaching the de Sitter attractor.

We study the possibility of the vortices formation in axion condensate on the galactic scale. Such vortices can occur as a result of global rotation of the early universe. We study analytical models of vortices and calculate exemplary galaxy rotation curves. Depending on the setup it is possible to obtain a variety of shapes which give a good qualitative agreement with observational results. However, as we show, the extremely low velocity dispersions of the axions velocity are required to form the single vortex on the galactic scales. We find that the required velocity dispersion is of the order of $\sigma \approx 10^{-12} \frac{m}{s}$. This is much smaller that predicted within the present understanding of the axion physics. The vortices in axion condensate can be however formed on the much smaller scales and give seeds to the galaxy formation and to their angular momenta. On the other hand, the vortices can be formed on the galactic scales, but only if the mass of the axion-like particles is of the order of $10^{-30}$eV. In this case, the particle de Broglie wavelength is comparable with the galactic diameter. This condition must be fulfilled in order to keep the coherence of the quantum condensate on the galactic scales.

We develop here a relatively simple description of dark energy based on the dynamics of non-minimally coupled to gravity phantom scalar field which, in limit, corresponds to cosmological constant. The dark energy equation of state, obtained directly from the dynamics of the model, turns out to be an oscillatory function of the scale factor. This parameterisation is compared to other possible dark energy parameterisations, among them, the most popular one, linear in the scale factor. We use the Bayesian framework for model selection and make a comparison in the light of SN Ia, CMB shift parameter, BAO A parameter, observational H(z) and growth rate function data. We find that there is evidence to favour a parameterisation with oscillations over {\it a priori} assumed linear one.

In this paper we study possible observational consequences of the bouncing cosmology. We consider a model where a phase of inflation is preceded by a cosmic bounce. While we consider in this paper only that the bounce is due to loop quantum gravity, most of the results presented here can be applied for different bouncing cosmologies. We concentrate on the scenario where the scalar field, as the result of contraction of the universe, is driven from the bottom of the potential well. The field is amplified, and finally the phase of the standard slow-roll inflation is realized. Such an evolution modifies the standard inflationary spectrum of perturbations by the additional oscillations and damping on the large scales. We extract the parameters of the model from the observations of the cosmic microwave background radiation. In particular, the value of inflaton mass is equal to $m=(2.6 \pm 0.6) \cdot 10^{13}$ GeV. In our considerations we base on the seven years of observations made by the WMAP satellite. We propose the new observational consistency check for the phase of slow-roll inflation. We investigate the conditions which have to be fulfilled to make the observations of the Big Bounce effects possible. We translate them to the requirements on the parameters of the model and then put the observational constraints on the model. Based on assumption usually made in loop quantum cosmology, the Barbero-Immirzi parameter was shown to be constrained by $\gamma<1100$ from the cosmological observations. We have compared the Big Bounce model with the standard Big Bang scenario and showed that the present observational data is not informative enough to distinguish these models.

We present the Bayesian analysis of four different types of backreation models, which are based on the Buchert equations. In this approach, one considers a solution to the Einstein equations for a general matter distribution and then an average of various observable quantities is taken. Such an approach became of considerable interest when it was shown that it could lead to agreement with observations without resorting to dark energy. In this paper we compare the LambdaCDM model and the backreation models with SNIa, BAO, and CMB data, and find that the former is favoured. However, the tested models were based on some particular assumptions about the relation between the average spatial curvature and the backreaction, as well as the relation between the curvature and curvature index. In this paper we modified the latter assumption, leaving the former unchanged. We find that, by varying the relation between the curvature and curvature index, we can obtain a better fit. Therefore, some further work is still needed -- in particular the relation between the backreaction and the curvature should be revisited in order to fully determine the feasibility of the backreaction models to mimic dark energy.

We study the dynamics of a non-minimally coupled scalar field cosmology with a potential function. We use the framework of dynamical systems theory to investigate all evolutional paths admissible for all initial conditions. Additionally, we assume the presence of barotropic matter and show that the dynamics can be formulated in terms of an autonomous dynamical system. We have found fixed points corresponding to three main stages of the evolution of the universe, namely, radiation, matter and quintessence domination epochs. Using the linearization of the dynamical systems in the vicinity of the critical points we explicitly obtain formulas determining the effective equation of state parameter for the universe at different epochs. In our approach the form of $w(z)$ parametrisation is derived directly from the dynamical equations rather than postulated {\it a priori}.

The main aim of this paper is to extend the early approach to quantum cosmogenesis provided by Fomin. His approach was developed independently to the well-known Tryon description of the creation of the closed universe as a process of quantum fluctuation of vacuum. We apply the Fomin concept to derive the cosmological observables. We argue that Fomin's idea from his 1973 work, in contrast to Tryon's one has impact on the current Universe models and the proposed extension of his theory now can be tested by distant supernovae SNIa. Fomin's idea of the creation of the Universe is based on the intersection of two fundamental theories: general relativity and quantum field theory with the contemporary cosmological models with dark energy. As a result of comparison with contemporary approaches concerning dark energy, we found out that Fomin's idea appears in the context of the present acceleration of the Universe explanation: cosmological models with decaying vacuum. Contemporary it appears in the form of Ricci scalar dark energy connected with the holographic principle. We show also that the Fomin model admits the bounce instead of the initial singularity. We demonstrate that the Fomin model of cosmogenesis can be falsified and using SNIa data the values of model parameters is in agreement with observations.

In this paper we study dynamics of the closed FRW model with holonomy corrections coming from loop quantum cosmology. We consider models with a scalar field and cosmological constant. In case of the models with cosmological constant and free scalar field, dynamics reduce to 2D system and analysis of solutions simplify. If only free scalar field is included then universe undergoes non-singular oscillations. For the model with cosmological constant, different behaviours are obtained depending on the value of $\Lambda$. If the value of $\Lambda$ is sufficiently small, bouncing solutions with asymptotic de Sitter stages are obtained. However if the value of $\Lambda$ exceeds critical value $\Lambda_{\text{c}} =\frac{\sqrt{3}m^2_{\text{Pl}}}{2\pi\gamma^3} \simeq 21 m^2_{\text{Pl}}$ then solutions become oscillatory. Subsequently we study models with a massive scalar field. We find that this model possess generic inflationary attractors. In particular field, initially situated in the bottom of the potential, is driven up during the phase of quantum bounce. This subsequently leads to the phase of inflation. Finally we find that, comparing with the flat case, effects of curvature do not change qualitatively dynamics close to the phase of bounce. Possible effects of inverse volume corrections are also briefly discussed.

In this publication we investigate dynamics of a flat FRW cosmological model with a non-minimally coupled scalar field with the coupling term $\xi R \psi^{2}$ in the scalar field action. The quadratic potential function $V(\psi)\propto \psi^{2}$ is assumed. All the evolutional paths are visualized and classified in the phase plane, at which the parameter of non-minimal coupling $\xi$ plays the role of a control parameter. The fragility of global dynamics with respect to changes of the coupling constant is studied in details. We find that the future big rip singularity appearing in the phantom scalar field cosmological models can be avoided due to non-minimal coupling constant effects. We have shown the existence of a finite scale factor singular point (future or past) where the Hubble function as well as its first cosmological time derivative diverges.

In this letter we apply dynamical system methods to study all evolutional paths admissible for all initial conditions of the FRW cosmological model with a non-minimally coupled to gravity scalar field and a barotropic fluid. We choose "energy variables" as phase variables. We reduce dynamics to a 3-dimensional dynamical system for an arbitrary potential of the scalar field in the phase space variables. After postulating the potential parameter $\Gamma$ as a function of $\lambda$ (defined as $-V'/V$) we reduce whole dynamics to a 3-dimensional dynamical system and study evolutional paths leading to current accelerating expansion. If we restrict the form of the potential then we will obtain a 2-dimensional dynamical system. We use the dynamical system approach to find a new generic quintessence scenario of approaching to the de Sitter attractor which appears only for the case of non-vanishing coupling constant.

In the paper we investigate observational constraints on coupling to gravity constant parameter $\xi$ using distant supernovae SNIa data, baryon oscillation peak (BOP), the cosmic microwave background radiation (CMBR) shift parameter, and $H(z)$ data set. We estimate the value of this parameter to constrain the extended quintessence models with nonminimally coupled to gravity phantom scalar field. The combined analysis of observational data favors a value of $\xi$ which lies in close neighborhood of the conformal coupling. While our estimations are model dependent they give rise to a indirect bound on the Equivalence Principle.

Recent astronomical observations of SNIa, CMB, as well as BAO in the Sloan Digital Sky Survey, suggest that the current Universe has entered a stage of an accelerated expansion with the transition redshift at $z \simeq 0.5$. While the simplest candidates to explain this fact is cosmological constant/vacuum energy, there exists a serious problem of coincidence. In theoretical cosmology we can find many possible approaches alleviating this problem by applying new physics or other conception of dark energy. We consider state of art candidates for the description of accelerating Universe in the framework of the Bayesian model selection. We point out advantages as well as troubles of this approach. We find that the combination of four data bases gives a stringent posterior probability of the $\Lambda$CDM model which is 74%. This fact is a quantitative exemplification of a turmoil in modern cosmology over the $\Lambda$ problem.

The contemporary cosmic expansion is considered in the context of Modified Friedmann Dynamics (MOFD). We discuss some relativistic model exploring analogy to MOND modification of Newtonian dynamics. We argue that MOFD cosmologies can explain fraction of dark matter in the accelerating Universe. We discuss some observational constraints on possible evolutional MOFD scenarios of cosmological models coming from SN Ia distant supernovae. We show that Modified Newtonian Dynamics can be obtained as a Newtonian limit of more general relativistic models with polytropic component of Equation of State. They constitute a special subclass of generalized Cardassian models basing on generalization of the Raychaudhuri equation rather than on generalization of the Friedmann first integral. We demonstrate that MOND cosmologies are compatible with observed accelerated phase of expansion of current universe only for high value of cosmological constant. The Bayesian framework of model selection favored this model over $\Lambda$CDM model if $\Omega_{m,o}$ is fixed but this evidence is not significant. Moreover obtained from statistical analysis value of the MOND characteristic $\beta$ parameter is far from value required for explanation of the flat rotation curves of spiral galaxies.

We investigate evolutional paths of an extended quintessence with a non-minimally coupled phantom scalar field $\psi$ to the Ricci curvature. The dynamical system methods are used to investigate typical regimes of dynamics at the late time. We demonstrate that there are two generic types of evolutional scenarios which approach the attractor (a focus or a node type critical point) in the phase space: the quasi-oscillatory and monotonic trajectories approach to the attractor which represents the FRW model with the cosmological constant. We demonstrate that dynamical system admits invariant two-dimensional submanifold and discussion that which cosmological scenario is realized depends on behavior of the system on the phase plane $(\psi, \psi')$. We formulate simple conditions on the value of coupling constant $\xi$ for which trajectories tend to the focus in the phase plane and hence damping oscillations around the mysterious value $w=-1$. We describe this condition in terms of slow-roll parameters calculated at the critical point. We discover that the generic trajectories in the focus-attractor scenario come from the unstable node. It is also investigated the exact form of the parametrization of the equation of state parameter $w(z)$ (directly determined from dynamics) which assumes a different form for both scenarios.

In our paper we constrain Friedman Robertson Walker (FRW) model with ``radiation-like'' contribution to the Friedmann equation against the astronomical data. We analyze the observational limitations on a $(1+z)^4$ term from SNIa data, FRIIb radio galaxy (RG) data, baryon oscillation peak and CMBR observations. We argue that it is not possible to determine the energy densities of individual components scaling like radiation from a kinematic astronomical test. The bounds for density parameter for total radiation-like term can be obtained. We find different interpretations of the presence of scaling like radiation term: FRW universe filled with a massless scalar field in a quantum regime (Casimir effect), Friedmann-Robertson-Walker (FRW) model in a semi-classical approximation of loop quantum gravity, FRW model in the Randall Sundrum scenario with dark radiation or cosmological model with global rotation. In this paper we mainly concentrate on Casimir effect arising from quantum effects of the scalar field. This contribution can describe decaying part of cosmological constant. We discuss the back reaction of gravity on Casimir-type force which is a manifestation of the vacuum fluctuations of the quantum scalar field at low temperature. It is shown that while the Casimir energy gives rise to accelerating Universe, the cosmological constant term is still required. We argue that a small negative contribution of a radiation-like term can reconcile the tension between the observed primordial ${}^4He$ and $D$ abundance. Moreover the presence of such contribution can also remove the disagreement between Hubble parameter $H_0$ values obtained from both SNIa and Wilkinson Microwave Anisotropy Probe (WMAP) satellite data.

All evolutional paths for all admissible initial conditions of FRW cosmological models with dissipative dust fluid (described by dark matter, baryonic matter and dark energy) are analyzed using dynamical system approach. With that approach, one is able to see how generic the class of solutions leading to the desired property -- acceleration -- is. The theory of dynamical systems also offers a possibility of investigating all possible solutions and their stability with tools of Newtonian mechanics of a particle moving in a 1-dimensional potential which is parameterized by the cosmological scale factor. We demonstrate that flat cosmology with bulk viscosity can be treated as a conservative system with a potential function of the Chaplygin gas type. We also confront viscous models with SNIa observations. The best fitted models are obtained by minimizing the $\chi^{2}$ function which is illustrated by residuals and $\chi^{2}$ levels in the space of model independent parameters. The general conclusion is that SNIa data supports the viscous model without the cosmological constant. The obtained values of $\chi^{2}$ statistic are comparable for both the viscous model and LCDM model. The Bayesian information criteria are used to compare the models with different power law parameterization of viscous effects. Our result of this analysis shows that SNIa data supports viscous cosmology more than the LCDM model if the coefficient in viscosity parameterization is fixed. The Bayes factor is also used to obtain the posterior probability of the model.

The oscillating scenario of route to Lambda was recently proposed by us arXiv:0704.1651 as an alternative to a cosmological constant in a explanation of the current accelerating universe. In this scenario phantom scalar field conformally coupled to gravity drives the accelerating phase of the universe. In our model $\Lambda$CDM appears as a global attractor in the phase space. In this paper we investigate observational constraints on this scenario from recent measurements of distant supernovae type Ia, CMB R shift, BAO and $H(z)$ observational data. The Bayesian methods of model selection are used in comparison the model with concordance $\Lambda$CDM one as well as with model with dynamical dark energy parametrised by linear form. We conclude that $\Lambda$CDM is favoured over FRW model with dynamical oscillating dark energy. Our analysis also demonstrate that FRW model with oscillating dark energy is favoured over FRW model with decaying dark energy parametrised in linear way.

We develop a simple method of dark energy reconstruction using a geometrical form of the luminosity-distance relation. In this method the FRW dynamical system with dark energy is reconstructed instead of the equation of state parameter. We give several examples which illustrate the usefulness of our method in fitting the redshift transition from the decelerating to accelerating phase as the value of the Hubble function at the transition.

In this letter we investigate acceleration in the flat cosmological model with a conformally coupled phantom field and we show that acceleration is its generic feature. We reduce the dynamics of the model to a 3-dimensional dynamical system and analyze it on a invariant 2-dimensional submanifold. Then the concordance FRW model with the cosmological constant $\Lambda$ is a global attractor situated on a 2-dimensional invariant space. We also study the behaviour near this attractor, which can be approximated by the dynamics of the linearized part of the system. We demonstrate that trajectories of the conformally coupled phantom scalar field with a simple quadratic potential crosses the cosmological constant barrier infinitely many times in the phase space. The universal behaviour of the scalar field and its potential is also calculated. We conclude that the phantom scalar field conformally coupled to gravity gives a natural dynamical mechanism of concentration of the equation of state coefficient around the magical value $w_{\text{eff}}=-1$. We demonstrate route to Lambda through the infinite times crossing the $w_{\text{eff}}=-1$ phantom divide.

We study the general chaotic features of dynamics of the phantom field modelled in terms of a single scalar field conformally coupled to gravity. We demonstrate that the dynamics of the FRW model with dark energy in the form of phantom field can be regarded as a scattering process of two types: multiple chaotic and classical non-chaotic. It depends whether the spontaneously symmetry breaking takes place. In the first class of models with the spontaneously symmetry breaking the dynamics is similar to the Yang-Mills theory. We find the evidence of a fractal structure in the phase space of initial conditions. We observe similarities to the phenomenon of a multiple scattering process around the origin. In turn the class of models without the spontaneously symmetry breaking can be described as the classical non-chaotic scattering process and the methods of symbolic dynamic are also used in this case. We show that the phantom cosmology can be treated as a simple model with scattering of trajectories which character depends crucially on a sign of a square of mass. We demonstrate that there is a possibility of chaotic behavior in the flat Universe with a conformally coupled phantom field in the system considered on non-zero energy level. We obtain that the acceleration is a generic feature in the considered model without the spontaneously symmetry breaking. We observe that the effective EOS coefficient oscillates and then approach to $w=-1$.

We discuss certain issues related to the limitations of density parameters for a ``radiation''-like contribution to the Friedmann equation using kinematical or geometrical measurements. We analyse the observational constraint of a negative $(1+z)^4$-type contribution in cosmological models. We argue that it is not possible to determine the energy densities of individual components of matter scaling like radiation from astronomical observations. We find three different interpretations of the presence of the radiation term: 1) the FRW universe filled with a massless scalar field in a quantum regime (the Casimir effect), 2) the Friedmann-Robertson-Walker (FRW) model in the Randall Sundrum scenario with dark radiation, 3) the cosmological model with global rotation. From supernovae type Ia (SNIa) data, Fanaroff-Riley type IIb (FRIIb) radio galaxy (RG) data, baryon oscillation peak and cosmic microwave background radiation (CMBR) observations we obtain bounds for the negative radiation-like term. A small negative contribution of dark radiation can reconcile the tension in nucleosynthesis and remove also the disagreement between $H_0$ values obtained from both SNIa and the Wilkinson Microwave Anisotropy Probe (WMAP) satellite data.

ecent observations of type Ia supernovae indicate that the Universe is in an accelerating phase of expansion. The fundamental quest in theoretical cosmology is to identify the origin of this phenomenon. In principle there are two possibilities: 1) the presence of matter which violates the strong energy condition (a substantial form of dark energy), 2) modified Friedmann equations (Cardassian models -- a non-substantial form of dark matter). We classify all these models in terms of 2-dimensional dynamical systems of the Newtonian type. We search for generic properties of the models. It is achieved with the help of Peixoto's theorem for dynamical system on the Poincar{\'e} sphere. We find that the notion of structural stability can be useful to distinguish the generic cases of evolutional paths with acceleration. We find that, while the $\Lambda$CDM models and phantom models are typical accelerating models, the cosmological models with bouncing phase are non-generic in the space of all planar dynamical systems. We derive the universal shape of potential function which gives rise to presently accelerating models. Our results show explicitly the advantages of using a potential function (instead of the equation of state) to probe the origin of the present acceleration. We argue that simplicity and genericity are the best guide in understanding our Universe and its acceleration.

We discuss observational constrains coming from supernovae imposed on the behaviour of the Randall-Sundrum models. The data indicates that, under the assumption that we admit zero pressure dust matter on the brane, the cosmological constant is still needed to explain current observations. The observations from supernovae give a large value of the density parameter for brane matter $\Omega_{\lambda,0} \simeq 0.01$ as the best fit. For high redshifts $z > 1.2$, the difference between the brane model and the $\Lambda$CDM (Perlmutter) model becomes detectable observationally. From the maximum likelihood method we obtained the favored value of $\Omega_{\lambda,0} = 0.004 \pm 0.016$ for $\Omega_{k,0}=0$ and $\Omega_{\mathrm{m},0}=0.3$. This gives the limit $\Omega_{\lambda,0} < 0.02$ at $1 \sigma$ level. While the model with brane effects is preferred by the supernovae type Ia data, the model without brane fluid is still statistically admissible. For comparison the limit on brane effects which comes from CMB anisotropies and BBN is also obtained. The uncertainty in the location of the first peak gives a stronger limit $\Omega_{\lambda,0} < 1.0 \cdot 10^{-12}$, whereas from BBN we obtain that $\Omega_{\lambda,0} < 1.0 \cdot 10^{-27}$. However, both very strict limits are obtained with the assumption that brane effects do not change the physics in the pre-recombination era, while the SN Ia limit is model independent.We demonstrate that the fit to supernovae data can also be obtained if we admit the phantom matter $p = - (4/3) \varrho$ on the brane, where this matter mimics the influence of the cosmological constant. Finally, we propose to check for dark radiation and brane tension by the application of the angular diameter of galaxies minimum value test.

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.

Recent astronomical observations indicate that the Universe is presently almost flat and undergoing a period of accelerated expansion. Basing on Einstein's general relativity all these observations can be explained by the hypothesis of a dark energy component in addition to cold dark matter (CDM). Because the nature of this dark energy is unknown, it was proposed some alternative scenario to explain the current accelerating Universe. The key point of this scenario is to modify the standard FRW equation instead of mysterious dark energy component. The standard approach to constrain model parameters, based on the likelihood method, gives a best-fit model and confidence ranges for those parameters. We always arbitrary choose the set of parameters which define a model which we compare with observational data. Because in the generic case, the introducing of new parameters improves a fit to the data set, there appears the problem of elimination of model parameters which can play an insufficient role. The Bayesian information criteria of model selection (the AIC and BIC) are dedicated to promotion a set of parameters which should be incorporated to the model. We divide class of all accelerating cosmological models into two groups according to the two types of explanation acceleration of the Universe. Then the Bayesian framework of model selection is used to determine the set of parameters which gives preferred fit to the SNIa data. We find a few of flat cosmological models which can be recommend by both the Bayes factor and Akaike information criterion.

In this paper we proposed to use the group of analysis of symmetries of the dynamical system to describe the evolution of the Universe. This methods is used in searching for the unknown equation of state. It is shown that group of symmetries enforce the form of the equation of state for noninteracting scaling multifluids. We showed that symmetries give rise the equation of state in the form $p=-\Lambda+w_{1}\rho(a)+w_{2}a^{\beta}+0$ and energy density $\rho=\Lambda+\rho_{01}a^{-3(1+w)}+\rho_{02}a^{\beta}+\rho_{03}a^{-3}$, which is commonly used in cosmology. The FRW model filled with scaling fluid (called homological) is confronted with the observations of distant type Ia supernovae. We found the class of model parameters admissible by the statistical analysis of SNIa data. We showed that the model with scaling fluid fits well to supernovae data. We found that $\Omega_{\text{m},0} \simeq 0.4$ and $n \simeq -1$ ($\beta = -3n$), which can correspond to (hyper) phantom fluid, and to a high density universe. However if we assume prior that $\Omega_{\text{m},0}=0.3$ then the favoured model is close to concordance $\Lambda$CDM model. Our results predict that in the considered model with scaling fluids distant type Ia supernovae should be brighter than in $\Lambda$CDM model, while intermediate distant SNIa should be fainter than in $\Lambda$CDM model. We also investigate whether the model with scaling fluid is actually preferred by data over $\Lambda$CDM model. As a result we find from the Akaike model selection criterion prefers the model with noninteracting scaling fluid.

The recent astronomical measurements of distant supernovae as well as other observations indicate that our universe is presently accelerating. There are different proposals for the explanation of this acceleration, such as the cosmological constant $\Lambda$, decaying vacuum energy, an evolving scalar field (quintessence), phantom energy, etc. Most of these proposals require the existence of exotic matter with negative pressure violating the strong energy condition. On the other hand, there have appeared many models which offer dramatically different mechanisms for the current acceleration, in which dark energy emerges from the gravity sector rather than from the matter sector. In this paper, we compare the concordance $\Lambda$CDM model with the Sahni-Shtanov brane-world models of dark energy by using the Akaike and Bayesian information criteria. We show that new parameters in the brane model are not statistically significant in terms of the information criteria, although the best fit method gives an improved fit to the SNIa data, because of the additional parameters. This is because the information criteria of model selection compensate for this advantage by penalizing models having more free parameters. We conclude that only new future observational data are accurate enough to give an advantage to dark-energy models of the brane origin, i.e., a very high-significance detection is required to justify the presence of new parameters. In our statistical analysis both Riess et al.'s and Astier et al.'s SNIa samples are used. For stringent constraining parameters of the models the baryon oscillation peak (BOP) test is used.

It is described dynamics of a large class of accelerating cosmological models in terms of dynamical systems of the Newtonian type. The evolution of the models is reduced to the motion of a particle in a potential well parameterized by the scale factor. This potential function can be reconstructed from distant supernovae type Ia data and many cosmological models represented in terms of the potential becomes in a good agreement with current observational data. It is proposed to use the information criteria to overcome this degeneracy within a class of A) dark energy models and B) simple models basing on modification of the FRW equation. Two class of models can be recommended by the Akaike (AIC) and Schwarz (BIC) information criteria: the phantom and $\Lambda$CDM models.

We investigate the possibility of testing cosmological models with interaction between matter and energy sector. We assume the standard FRW model while the so called energy conservation condition is interpreted locally in terms of energy transfer. We analyze two forms of dark energy sectors: the cosmological constant and phantom field. We find a simple exact solution of the models in which energy transfer is described by a Cardassian like term in the relation of $H^{2}(z)$, where $H$ is Hubble's function and $z$ is redshift. The considered models have two additional parameters $(\Omega_{\text{int}},n)$ (apart the parameters of the $\Lambda$CDM model) which can be tested using SNIa data. In the estimation of the model parameters Riess et al.'s sample is used. We also confront the quality of statistical fits for both the $\Lambda$CDM model and the interacting models with the help of the Akaike and Bayesian informative criteria. Our conclusion from standard best fit method is that the interacting models explains the acceleration of the Universe better but they give rise to a universe with high matter density. However, using the tools of information criteria we find that the two new parameters play an insufficient role in improving the fit to SNIa data and the standard $\Lambda$CDM model is still preferred. We conclude that high precision detection of high redshift supernovae could supply data capable of justifying adoption of new parameters.

We study the Friedmann-Robertson-Walker model with phantom fields modelled in terms of scalar fields. We apply the Ziglin theory of integrability and find that the flat model is nonintegrable. Then we cannot expect to determine simple analytical solutions of the Einstein equations. We demonstrate that there is only a discrete set of parameters where this model is integrable. For comparison we describe the phantoms fields in terms of the barotropic equation of state. It is shown that in the contrast to the phantoms modelled as scalar fields, the dynamics is always integrable and phase portraits are constructed. In this case we find the duality relation.

Recent measurements of distant type Ia supernovae (SNIa) as well as other observations indicate that our universe is in accelerating phase of expansion. In principle there are two alternative explanation for such an acceleration. While in the first approach an unknown form of energy violating the strong energy condition is postulated, in second one some modification of FRW dynamics is postulated. The both approaches are in well agreement with present day observations which is the manifestation of the degeneracy problem appearing in observational cosmology. We use the Akaike (AIC) and Bayesian (BIC) information criteria of model selection to overcome this degeneracy and to determine a model with such a set of parameters which gives the most preferred fit to the SNIa data. We consider five representative evolutional scenarios in each of groups. Among dark energy proposal the $\Lambda$CDM model, CDM model with phantom field, CDM model with topological defect, model with Chaplygin gas, and the model with the linear dynamical equation of state parameter. As an alternative prototype scenarios we consider: brane world Dvali Gabadadze Porrati scenario, brane models in Randall-Sundrum scenario, Cardassian models with dust matter and radiation, bouncing model with the cosmological constant and metric-affine gravity (MAG) inspired cosmological models. Applying the model selection criteria we show that both AIC and BIC indicates that additional contribution arises from nonstandard FRW dynamics are not necessary to explain SNIa. Adopting the model selection information criteria we show that the AIC indicates the flat phantom model while BIC indicates both flat phantom and flat $\Lambda$CDM models.

The observations of SNIa suggest that we live in the acceleration epoch when the densities of the cosmological constant term and matter are almost equal. This leads to the cosmic coincidence conundrum. As the explanation for this problem we propose the FRW model with dark matter and dark energy which interact each other exchanging energy. We show that the cubic correction to the Hubble law, measured by distant supernovae type Ia, probes this interaction. We demonstrate that influences between nonrelativistic matter and vacuum sectors are controlled by third and higher derivatives of the scale factor. As an example we consider flat decaying $\Lambda(t)$ FRW cosmologies. We point out the possibility of measure of the energy transfer by the cubic and higher corrections to Hubble's law. The statistical analysis of SNIa data is used as an evidence of energy transfer. We find that there were the transfer from the dark energy sector to the dark matter one without any assumption about physics governing this process. We confront this hypothesis about the transfer with SNIa observations and find that the transfer the phantom and matter sector is admissible for $\Omega_{\text{m},0}=0.27$. We also demonstrate that it is possible to differentiate between the energy transfer model and the variable coefficient equation of state model.

We study metric-affine gravity (MAG) inspired cosmological models. Those models were statistically estimated using the SNIa data. We also use the cosmic microwave background observations and the big-bang nucleosynthesis analysis to constrain the density parameter $\Omega_{\psi,0}$ which is related to the non-Riemannian structure of the underlying spacetime. We argue that while the models are statistically admissible from the SNIa analysis, complementary stricter limits obtained from the CMB and BBN indicate that the models with density parameters with a $a^{-6}$ scaling behaviour are virtually ruled out. If we assume the validity of the particular MAG based cosmological model throughout all stages of the Universe, the parameter estimates from the CMB and BBN put a stronger limit, in comparison to the SNIa data, on the presence of non-Riemannian structures at low redshifts.

In the present paper we raise the question whether initial cosmological singularity can be proved from the cosmological tests. The classical general relativity predict the existence of singularity in the past if only some energy conditions are satisfied. On the other hand the latest quantum gravity applications to cosmology suggest of possibility of avoiding the singularity and replace it with the bounce. The distant type Ia supernovae data are used to constraints on bouncing evolutional scenario where square of the Hubble function $H^2$ is given by formulae $H^2=H^2_0[\Omega_{m,0}(1+z)^{m}-\Omega_{n,0}(1+z)^{n}]$, where $\Omega_{m,0}, \Omega_{n,0}>0$ are density parameters and $n>m>0$. We show that the on the base of the SNIa data standard bouncing models can be ruled out on the $4\sigma$ confidence level. If we add the cosmological constant to the standard bouncing model then we obtain as the best-fit that the parameter $\Omega_{n,0}$ is equal zero which means that the SNIa data do not support the bouncing term in the model. The bounce term is statistically insignificant the present epoch. We also demonstrate that BBN offer the possibility of obtaining stringent constraints of the extra term $\Omega_{n,0}$. The other observational test methods like CMB and the age of oldest objects in the Universe are used. We also use the Akaike informative criterion to select a model according to the goodness of fit and we conclude that this term should be ruled out by Occam's razor, which makes that the big bang is favored rather then bouncing scenario.

We discuss the relation between angular momenta and masses of galaxy structures base on the Li model of the universe with global rotation. In our previous paper (God{\l}owski et al 2002) it was shown that the model predicts the presence of a minimum in this relation. In the present paper we discuss observational evidence allowing us to verify this relation. We find null angular momentum J=0 for the masses corresponding to mass of galaxy grups and non-vanishing angular momenta for other galactic structures. We check these theoretical predictions analysing Tully's galaxy grups. The existing data comparing alignment in different galactic structure are consistent with obtained theoretical relation $J(M)$ if we interpret the groving alignment as the galactic increasing angular momenta in the galactic structure.

The so called Generalized Chaplygin Gas (GCG) with the equation of state $p = - \frac{A}{{\rho}^{\alpha}}$ was recently proposed as a candidate for dark energy in the Universe. In this paper we confront the GCG with SNIa data. Specifically we have tested the GCG cosmology in three different classes of models with (1) $\Omega_m= 0.3$, $\Omega_{Ch}= 0.7$; (2) $\Omega_m= 0.05$, $\Omega_{Ch}= 0.95$ and (3) $\Omega_m = 0$, $\Omega_{Ch} = 1$, as well as the model withouth any assumption on $\Omega_m$. The best fitted models are obtained by minimalizing the $\chi^2$ function and $\chi^2$ levels in the $(A_0, \alpha)$ plane. We supplemented our analysis with confidence intervals in the $(A_0, \alpha)$ plane, as well as one-dimensional probability distribution functions for models parameter. The general conclusion is that SNIa data strongly support the Chaplygin gas (with $\alpha = 1$). Extending our analysisby relaxing the flat prior lead to the result that even though the best fitted values of $\Omega_k$ are formally non-zero, still they are close to flat case. It should be viewed as an advantage of the GCG model since in similar analysisof $\Lambda$CDM model high negative value of $\Omega_{k}$ were found to be bestfitted to the data and independent inspiration from CMBR and extragalactic astronomy has been invoked to fix the curvature problem. Our results show clearly that in Generalized Chaplygin Gas cosmology distant $z >1$ supernovae should be brighter than in $\Lambda$CDM model.This prediction seems to be confirmed with new Riess high redshift SNIa sample. Moreover, we argue that with the future SNAP data it would be possible to differentiate between models with various value of $\alpha$ parameter and/or discriminated between GCG, Cardassian and $\Lambda$CDM models

We consider a Friedmann-Robertson-Walker cosmological model with some exotic perfect fluid with spin known as the Weyssenhoff fluid. The possibility that the dark energy may be described in part by the Weyssenhoff fluid is discussed. The observational constraint coming from supernovae type Ia observations is established. This result indicates that, whereas the cosmological constant is still needed to explain current observations, the model with spin fluid is admissible. For high redshifts $z > 1$ the differences between the model with spin fluid and the cold dark matter model with a cosmological constant become detectable observationally for the flat case with $\Omega_{\text{m},0}=0.3$. From the maximum likelihood method we obtain the value of $\Omega_{\text{s},0} = 0.004 \pm 0.016$. This gives us the limit $\Omega_{\text{s},0} > -0.012$ at the $1\sigma$ level. While the model with ``brane effects'' is preferred by the supernovae Ia data, the model with spin fluid is statistically admissible. For comparison, the limit on the spin fluid coming from cosmic microwave background anisotropies is also obtained. The uncertainties in the location of a first peak give the interval $-1.4 \times 10^{-10} < \Omega_{\text{s},0} < -10^{-10}$. From big bang nucleosynthesis we obtain the strongest limit $\Omega_{\text{s},0} \gtrsim -10^{-20}$. The interconnection between the model considered and brane models is also pointed out.

We analyze dynamics of the FRW models with global rotation in terms of dynamical system methods. We reduce dynamics of these models to the FRW models with some fictitious fluid which scales like radiation matter. This fluid mimics dynamically effects of global rotation. The significance of the global rotation of the Universe for the resolution of the acceleration and horizon problems in cosmology is investigated. It is found that dynamics of the Universe can be reduced to the two-dimensional Hamiltonian dynamical system. Then the construction of the Hamiltonian allows for full classification of evolution paths. On the phase portraits we find the domains of cosmic acceleration for the globally rotating universe as well as the trajectories for which the horizon problem is solved. We show that the FRW models with global rotation are structurally stable. This proves that the universe acceleration is due to the global rotation. It is also shown how global rotation gives a natural explanation of the empirical relation between angular momentum for clusters and superclusters of galaxies. The relation $J \sim M^2$ is obtained as a consequence of self similarity invariance of the dynamics of the FRW model with global rotation. In derivation of this relation we use the Lie group of symmetry analysis of differential equation.

It is proposed to use the Lie group theory of symmetries of differential equations to investigate the system of equations describing a static star in a radiative and convective equilibrium. It is shown that the action of an admissible group induces a certain algebraic structure in the set of all solutions, which can be used to find a family of new solutions. We have demonstrated that, in the most general case, the equations admit an infinite parameter group of quasi-homologous transformations. We have found invariants of the symmetries group which correspond to the fundamental relations describing a physical characteristic of the stars such as the Hertzsprung-Russell diagram or the mass-luminosity relation. In this way we can suggest that group invariants have not only purely mathematical sense, but their forms are closely associated with the basic empirical relations.

We show the effectiveness of using the procedure of the potential reconstruction (the V-method) in the estimation of the transition redshift corresponding to switch from deceleration to acceleration phase of the Universe. We investigate the FRW models with dark energy by applying the particle-like description of their dynamics. In this picture the evolution of the universe is represented by motion of a fictitious particle in the one-dimensional potential $V(a)$ or V(z)$ where $a$ is the scale factor and $z$ is the redshift. The V-method solves the inverse problem where we reconstruct the potential function from empirical data. We use the Riess et al.'s gold and gold+silver samples of SN Ia data in this reconstruction. In the same framework we obtain both the estimates of the form of the equation of state for dark energy and the transition redshift. We obtain that transition redshift is $0.38_{-0.07}^{+0.10}$ (gold subset) and $0.37_{-0.07}^{+0.06}$ (gold+silver sample) when the linear model of dark energy ($w(z)=w_{0} + w_{1}z$) and $\Omega_{m,0}=0.3$ are assumed. We compare the estimation of transition redshift with the Riess et al. results. The V-method also alows us to find the value of Hubble function at the moment of transition.

We discuss the equation of motion of the rotating homogenous and isotropic model of the Universe. We show that the model predicts the presence of a minimum in the relation between the mass of an astronomical object and its angular momentum. We show that this relation appears to be universal, and we predict the masses of structures with minimal angular momenta in agreement with observations. In such a manner we suggest the possibility at acquirement of angular momenta of celestial bodies during their formation from the global rotation of the Universe.

We discuss observational constrains coming from supernovae Ia imposed on the behaviour of the Randall-Sundrum models. In the case of dust matter on the brane, the difference between the best-fit Perlmutter model with a $\Lambda$-term and the best-fit brane models becomes detectable for redshifts $z > 1.2$. It is interesting that brane models predict brighter galaxies for such redshifts which is in agreement with the measurement of the $z = 1.7$ supernova. We also demonstrate that the fit to supernovae data can also be obtained, if we admit the "super-negative" dark energy (phantom matter) $p= - (4/3) \varrho$ on the brane, where the dark energy in a way mimics the influence of the cosmological constant. It also appears that the dark energy enlarges the age of the universe which is demanded in cosmology. Finally, we propose to check for dark radiation and brane tension by the application of the angular diameter of galaxies minimum value test. We point out the existence of coincidence problem for the brane tension parameter.

A class of Stephani cosmological models as a prototype of non-homogeneous universe is considered. The non-homogeneity can lead to accelerated evolution which is now observed from the SNIa data. Three samples of type Ia supernovae obtained by Perlmutter at al., Tonry et al. and Knop et al. are taken into account. Different statistical methods (best fits as well as maximum likelihood method) to obtain estimates of the model parameters are used. Stephani model is considered as an alternative to the concordance of Lambda-CDM model in the explanation of the present acceleration of the universe. The model explains the acceleration of the universe at the same level of accuracy as the Lambda-CDM model. From the best fit analysis it follows that the Stephani model is characterized by higher value of nonrelativistic matter density parameter than the Lambda-CDM model. It is also shown that the obtained results are consistent with location of CMB peaks.

We investigate a class of Cardassian scenarios of the universe evolution in notions of the qualitative theory of dynamical systems. This theory allows us to analyze all solutions for all possible initial conditions on the phase plane. In the Cardassian models we find that big-rip singularities are present as a typical behavior in the future if $n<0$. Some exact solutions for the flat Cardassian models as well as a duality relation were found. In turn from the statistical analysis of Knop's SNIa data, without any priors on matter content in the model, we obtain that at the 99% confidence level this big-rip scenario will reach. The potential function for the Hamiltonian description of dynamics is reconstructed from the SNIa data (inverse dynamical problem). We also pointed out the statistical analysis results depend oversensitively on the choice of the model parameter $\mathcal{M}$.

We discuss some observational constraints resulting only from SNIa observations, imposed on the behavior of the original flat Cardassian model, and its extension with the curvature term included. We test the models using the Perlmutter SN Ia data as well as the new Knop and Tonry samples. We estimate the Cardassian model parameters using the best-fitting procedure and the likelihood method. In the fitting procedure we use the density variables for matter, Cardassian fluid and curvature, and include the errors in redshift measurements. For the Perlmutter sample in the non-flat Cardassian model we obtain the high or normal matter density universe ($\Omega_{m,0} \ge 0.3$), while for the flat Cardassian model we have the high density universe. For sample A in the high density universe we also find the negative values of estimates of $n$ which can be interpreted as the phantom fluid effect. For the likelihood method we get that a nearly flat universe is preferred. We show that, if we assume that the matter density is 0.3, then $n \approx 0$ in the flat Cardassian model, which corresponds to the Perlmutter model with the cosmological constant. Testing with the Knop and Tonry SN Ia samples show no significant differences.

We demonstrate a model-independent method of estimating qualitative dynamics of an accelerating universe from observations of distant type Ia supernovae. Our method is based on the luminosity-distance function, optimized to fit observed distances of supernovae, and the Hamiltonian representation of dynamics for the quintessential universe with a general form of equation of state $p=w(a(z))\rho$. Because of the Hamiltonian structure of FRW dynamics with the equation of state $p = w(a(z)) \rho$, the dynamics is uniquelly determined by the potential function $V(a)$ of the system. The effectiveness of this method in discrimination of model parameters of Cardassian evolution scenario is also given. Our main result is the following, restricting to the flat model with the current value of $\Omega_{m,0}=0.3$, the constraints at $2\sigma$ confidence level to the presence of $\rho^{n}$ modification of the FRW models are $-0.50 \lesssim n \lesssim 0.36$.

We discuss a class of phantom ($p < - \varrho$) cosmological models. Except for phantom we admit various forms of standard types of matter and discuss the problem of singularities for these cosmologies. The singularities are different from those of standard matter cosmology since they appear for infinite values of the scale factor. We also find an interesting relation between the phantom models and standard matter models which is like the duality symmetry of string cosmology.

We apply methods of dynamical systems to study the behaviour of universe dominated by the generalized Chaplygin gas. We reduce the dynamics to a 2-dimensional Hamiltonian system and study its behaviour for various ranges of parameters. The dynamics is studied on the phase plane by using methods of qualitative analysis of differential equations. The behaviour of trajectories at infinity is studied in some convenient coordinates introduced on the phase plane. Hence we shown that FRW model with the generalized Chaplygin gas is structurally stable. We clearly find the domains of cosmic acceleration as well as conditions for which the horizon problem is solved. We also define some general class of fluids which generalize the Chaplygin gas. The dynamics of such models in terms of energy conditions is also discussed.

We discuss the problem of universe acceleration driven by global rotation. The redshift-magnitude relation is calculated and discussed in the context of SN Ia observation data. It is shown that the dynamics of considered problem is equivalent to the Friedmann model with additional non-interacting fluid with negative pressure. We demonstrate that the universe acceleration increase is due to the presence of global rotation effects, although the cosmological constant is still required to explain the SN Ia data. We discuss some observational constraints coming from SN Ia imposed on the behaviour of the homogeneous Newtonian universe in which matter rotates relative local gyroscopes. In the Newtonian theory $\Omega_{\text{r},o}$ can be identified with $\Omega_{\omega,0}$ (only dust fluid is admissible) rotation can exist with $\Omega_{\text{r},0}=\Omega_{\omega,0} \le 0$. However, the best-fit flat model is the model without rotation, i.e., $\Omega_{\omega,0}=0$. In the considered case we obtain the limit for $\Omega_{\omega,0}>-0.033$ at 1\sigma$ level$. We are also beyond the model and postulate the existence od additional matter which scales like radiation matter and then analyse how that model fits the SN Ia data. In this case the limits on rotation coming from BBN and CMB anisotropies are also obtained. If we assume that the current estimates are $\Omega_{\text{m},0}\sim 0.3$, $\Omega_{\text{r},0} \sim 10^{-4}$, then the SN Ia data show that $\Omega_{\omega,0} \ge -0.01$ (or $\omega < 2.6 \cdot 10^{-19} \text{rad}/\text{s}$). The statistical analysis gives us that the interval for any matter scaling like radiation is $\Omega_{\text{r},0} in (-0.01, 0.04)$.

We discuss observational constrains coming from supernovae Ia \cite{Perlmutter99} imposed on the behaviour of the Randall-Sundrum models. In the case of dust matter on the brane, the difference between the best-fit general relativistic model with a $\Lambda$-term \cite{Perlmutter99} and the best-fit brane models becomes detectable for redshifts $z > 0.6$. It is interesting that brane models predict brighter galaxies for such redshifts which is in agreement with the measurement of the $z = 1.7$ supernova \cite{Riess01} and with the New Data from the High Z Supernovae Search Team \cite{schmit02}. We also demonstrate that the fit to supernovae data can also be obtained, if we admit the "super-negative" dark energy $p = - (4/3) \varrho$ on the brane, where the dark energy in a way mimics the influence of the cosmological constant. It also appears that the dark energy enlarges the age of the universe which is demanded in cosmology. Finally, we propose to check for dark radiation and brane tension by the application of the angular diameter of galaxies minimum value test.

We apply methods of dynamical systems to study the behaviour of the Randall-Sundrum models. We determine evolutionary paths for all possible initial conditions in a 2-dimensional phase space and we investigate the set of accelerated models. The simplicity of our formulation in comparison to some earlier studies is expressed in the following: our dynamical system is a 2-dimensional Hamiltonian system, and what is more advantageous, it is free from the degeneracy of critical points so that the system is structurally stable. The phase plane analysis of Randall-Sundrum models with isotropic Friedmann geometry clearly shows that qualitatively we deal with the same types of evolution as in general relativity, although quantitatively there are important differences.