Dynamical systems methods are used to investigate cosmological model with non-minimally coupled scalar field. Existence of an asymptotically unstable de Sitter state distinguishes values of the non-minimal coupling constant parameter $\frac{3}{16}\le\xi<\frac{1}{4}$, which correspond to conformal coupling in higher dimensional theories of gravity.

In dynamical system describing evolution of universe with the flat Friedmann-Robertson-Walker symmetry filled with barotropic dust matter and non-minimally coupled scalar field with a constant potential function an invariant manifold of the de Sitter state is used to obtain exact solutions of the reduced dynamics. Using observational data coming from distant supernovae type Ia, the Hubble function $H(z)$ measurements and information coming from the Alcock-Paczy$\'n$ski test we find cosmological constraints on the non-minimal coupling constant $\xi$ between the scalar curvature and the scalar field. For all investigated models we can exclude negative values of this parameter at the $68\%$ confidence level. We obtain constraints on the non-minimal coupling constant consistent with condition for conformal coupling of the scalar field in higher dimensional theories of gravity.

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

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

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

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.

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.

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.

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