We investigate the initial condition which are simultaneously consistent with perturbation equations in both the radiation-dominated and the mater-dominated epochs. The exact formula for the spectrum transfer-function is derived.

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.

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