We calculate the quantum gravitational corrections to the Mukhanov-Sasaki equation obtained by the canonical quantization of the inflaton-gravity system. Our approach, which is based on the Born-Oppenheimer decomposition of the resulting Wheeler-DeWitt equation, was previously applied to a minimally coupled inflaton. In this article we examine the case of a non minimally coupled inflaton and, in particular, the induced gravity case is also discussed. Finally, the equation governing the quantum evolution of the inflationary perturbations is derived on a de Sitter background. Moreover the problem of the introduction of time is addressed and a generalized method, with respect to that used for the minimal coupling case, is illustrated. Such a generalized method can be applied to the universe wave function when, through the Born-Oppenheimer factorization, we decompose it into a part which contains the minisuperspace degrees of freedom and another which describes the perturbations.

The quasi-isotropic expansion for a simple two-fluid cosmological model, including radiation and string gas is constructed. The first non-trivial order expressions for the metric coefficients, energy densities and velocities are explicitly written down. Their small and large time asymptotics are studied. It is found that the large time asymptotic for the anisotropic component of the metric coefficients grows faster than that of the isotropic (trace-proportional) component.

We study the Wheeler-DeWitt equation for a class of induced gravity models in the minisuperspace approximation. In such models a scalar field nonminimally coupled to gravity determines the effective Newton's constant. For simplicity our analysis is limited to power-law potentials for the scalar field which have exact classical solutions. We show that these models have exact solutions also when quantised. Finally the Einstein Frame form of these solutions is obtained and a classical-quantum correspondence is found. Realistic induced gravity models also must include a symmetry breaking term which is needed in order to obtain a gravitational constant, successful inflation and a subsequent standard cosmological evolution. Nonetheless the potentials considered are important as they may describe the inflationary phase when the symmetry breaking part of the potential is negligible.

We study a mechanism for the amplification of the inflationary scalar perturbation when the inflaton field action is non-canonical, i.e. the inflaton kinetic term has a non-standard form. For such a case the speed of sound of the perturbations generated during inflation is less than one and in general changes with time. Furthermore in such models, even when the scalar field potential is negligible, diverse inflationary attractors may exist. The possible effects of a speed of sound approaching zero during some stage of inflation may lead to a large amplification for the amplitude of the scalar spectrum which, on horizon re-entry during the radiation dominated phase, can collapse and form primordial black holes (PBH) of a mass $M_{\rm BH}\sim 10^{-15}M_{\odot}$ which may constitute a large fraction of the total Dark Matter (DM) today.

We study the cosmological evolution and singularity crossing in the Bianchi-I universe filled with a conformally coupled scalar field and compare them with those of the Bianchi-I universe filled with a minimally coupled scalar field. We also write down the solution for the Bianchi-I Universe in the induced gravity cosmology.

We consider a rather simple method for the description of the Big Bang - Big Crunch cosmological singularity crossing. For the flat Friedmann universe this method gives the same results as more complicated methods, using Weyl symmetry or the transitions between the Jordan and Einstein frames. It is then easily generalized for the case of a Bianchi-I anisotropic universe. We also present early-time and late-time asymptotic solutions for a Bianchi-I universe, filled with a conformally coupled massless scalar field.

We study the relation between the Jordan-Einstein frame transition and the possible description of the crossing of singularities in flat Friedmann universes, using the fact that the regular evolution in one frame can correspond to crossing singularities in the other frame. We show that some interesting effects arise in simple models such as one with a massless scalar field or another wherein the potential is constant in the Einstein frame. The dynamics in these models and in their conformally coupled counterparts are described in detail, and a method for the continuation of such cosmological evolutions beyond the singularity is developed. We compare our approach with some other, recently developed, approaches to the problem of the crossing of singularities.

We illustrate how it is possible to calculate the quantum gravitational effects on the spectra of primordial scalar/tensor perturbations starting from the canonical, Wheeler-De Witt, approach to quantum cosmology. The composite matter-gravity system is analysed through a Born-Oppenheimer approach in which gravitation is associated with the heavy degrees of freedom and matter (here represented by a scalar field) with the light ones. Once the independent degrees of freedom are identified the system is canonically quantised. The differential equation governing the dynamics of the primordial spectra with its quantum-gravitational corrections is then obtained and is applied to diverse inflationary evolutions. Finally, the analytical results are compared to observations through a Monte Carlo Markov Chain technique and an estimate of the free parameters of our approach is finally presented and the results obtained are compared with previous ones.

We study integrable models in the Bianchi I metric case with scalar fields minimally and non-minimally coupled with gravity and the correspondence between their general solutions. Using the model with a minimally coupled scalar field and a constant potential as an example, we demonstrate how to obtain the general solutions of the corresponding models in the Jordan frame.

We study the integrable model with minimally and non-minimally coupled scalar fields and the correspondence of their general solutions. Using the model with a minimally coupled scalar field and a the constant potential as an example we demonstrate the difference between the general solutions of the corresponding models in the Jordan and the Einstein frames.

We consider the precession of a Dirac particle spin in some anisotropic Bianchi universes. This effect is present already in the Bianchi-I universe. We discuss in some detail the geodesics and the spin precession for both the Kasner and the Heckmann-Schucking solutions. In the Bianchi-IX universe the spin precession acquires the chaotic character due to the stochasticity of the oscillatory approach to the cosmological singularity. The related helicity flip of fermions in the very early Universe may produce the sterile particles contributing to dark matter.

We consider the relation between exact solutions of cosmological models having minimally and non-minimally coupled scalar fields. This is done for a particular class of solvable models which, in the Einstein frame, have potentials depending on hyperbolic functions and in the Jordan frame, where the non-minimal coupling is conformal, possess a relatively simple dynamics. We show that a particular model in this class can be generalized to the cases of closed and open Friedmann universes and still exhibits a simple dynamics. Further we illustrate the conditions for the existences of bounces in some sub-classes of the set of integrable models we have considered.

In the framework of a general scalar-tensor theory, we investigate the equivalence between two different parametrizations of fields that are commonly used in cosmology - the so-called Jordan frame and Einstein frame. While it is clear that both parametrizations are mathematically equivalent at the level of the classical action, the question about their mathematical equivalence at the quantum level as well as their physical equivalence is still a matter of debate in cosmology. We analyze whether the mathematical equivalence still holds when the first quantum corrections are taken into account. We explicitly calculate the one-loop divergences in both parametrizations by using the generalized Schwinger-DeWitt algorithm and compare both results. We find that the quantum corrections do not coincide off shell and hence induce an off shell dependence on the parametrization. According to the equivalence theorem, the one-loop divergences should however coincide on shell. For a cosmological background, we show explicitly that the on shell equivalence is indeed realized by a nontrivial cancellation.

The spectrum of primordial perturbations obtained by calculating the quantum gravitational corrections to the dynamics of scalar perturbations is compared with Planck 2013 and BICEP2/{\it Keck Array} public data. The quantum gravitational effects are calculated in the context of a Wheeler-De Witt approach and have quite distinctive features. We constrain the free parameters of the theory by comparison with observations.

We solve a general equation describing the lowest order corrections arising from quantum gravitational effects to the spectrum of cosmological fluctuations. The spectra of scalar and tensor perturbations are calculated to first order in the slow roll approximation and the results are compared with the most recent observations. The slow roll approximation gives qualitatively new quantum gravitational effects with respect to the pure de Sitter case.

We obtain general solutions for some flat Friedmann universes filled with a scalar field in induced gravity models and models including the Hilbert-Einstein curvature term plus a scalar field conformally coupled to gravity. As is well known, these models are connected to minimally coupled models through the combination of a conformal transformation and a transformation of the scalar field. The explicit forms of the self-interaction potentials for six exactly solvable models are presented here. We obtain the general solution for one of the integrable models, namely, the induced gravity model with a power-law potential for the self-interaction of the scalar field. We argue that although being mathematically in a one-to-one correspondence with the solutions in the minimally coupled models, the solutions in the corresponding non-minimally coupled models are physically different. This is because the cosmological evolutions seen by an internal observer connected with the cosmic time can be quite different. The study of a few induced gravity models with particular potentials gives us an explicit example of such a difference.

A general equation, describing the lowest order corrections coming from quantum gravitational effects to the spectrum of cosmological scalar fluctuations is obtained. These corrections are explicitly estimated for the case of a de Sitter evolution.

We investigate particular cosmological models, based either on tachyon fields or on perfect fluids, for which soft future singularities arise in a natural way. Our main result is the description of a smooth crossing of the soft singularity in models with an anti-Chaplygin gas or with a particular tachyon field in the presence of dust. Such a crossing is made possible by certain transformations of matter properties. Some of these cosmological evolutions involving tachyons are compatible with SNIa data. We compute numerically their dynamics involving a first soft singularity crossing, a turning point and a second soft singulatity crossing during recollapse, ending in a Big Crunch singularity.

Both dark energy models and modified gravity theories could lead to cosmological evolutions different from either the recollapse into a Big Crunch or exponential de Sitter expansion. The newly arising singularities may represent true endpoints of the evolution or allow for the extension of geodesics through them. In the latter case the components of the Riemann tensor representing tidal forces diverge. Sudden Future Singularities (SFS) occur at finite time, finite scale factor and finite Hubble parameter, only the deceleration parameter diverges. The energy density of a perfect fluid is regular and its pressure diverges at the SFS. A particular SFS, the Big Brake occurs when the energy density vanishes and the expansion arrives at a full stop at the singularity. Such scenarios are generated by either a particular scalar field (the tachyon field) or the anti-Chaplygin gas. By adding matter to these models, an unwanted feature appears: at the finite scale factor of the SFS the energy density of matter remains finite, implying (for a spatially flat universe) a finite Hubble parameter, hence finite expansion rate, rather then full stop. The universe would further expand through the singularity, this nevertheless seems forbidden as the energy density of the tachyonic field / anti-Chaplygin gas would become ill-defined. This paradox is relieved in the case of the anti-Chaplygin gas by redefining its energy density and pressure in terms of distributions peaked on the singularity. The regular cosmological quantities, continuous across the SFS are the energy density and the square of the Hubble parameter; those allowing for a jump at the SFS are the Hubble parameter and expansion rate (both being mirror-symmetric). The pressure and the decelaration parameter will contain Dirac delta-function contributions peaked on the SFS, however they anyhow would diverge at the singularity.

We employ the superpotential technique for the reconstruction of cosmological models with a non-minimally coupled scalar field evolving on a spatially flat Friedmann-Robertson-Walker background. The key point in this method is that the Hubble parameter is considered as a function of the scalar field and this allows one to reconstruct the scalar field potential and determine the dynamics of the field itself, without a priori fixing the Hubble parameter as a function of time or of the scale factor. The scalar field potentials that lead to de Sitter or asymptotic de Sitter solutions, and those that reproduce the cosmological evolution given by Einstein-Hilbert action plus a barotropic perfect fluid, have been obtained.

We consider the renormalization group improvement in the theory of the Standard Model (SM) Higgs boson playing the role of an inflaton with a strong non-minimal coupling to gravity. At the one-loop level with the running of constants taken into account, it leads to a range of the Higgs mass that is entirely determined by the lower WMAP bound on the cosmic microwave background (CMB) spectral index. We find that the SM phenomenology is sensitive to current cosmological data, which suggests to perform more precise CMB measurements as a SM test complementary to the LHC program. By using the concept of a field-dependent cutoff, we show the naturalness of the gradient and curvature expansion in this model within the conventional perturbation theory range of the SM. We also discuss the relation of these results to two-loop calculations and the limitations of the latter caused by parametrization and gauge dependence problems.

A cosmological model of a flat Friedmann universe filled with a mixture of anti-Chaplygin gas and dust-like matter exhibits a future soft singularity, where the pressure of the anti-Chaplygin gas diverges (while its energy density is finite). Despite infinite tidal forces the geodesics pass through the singularity. Due to the dust component, the Hubble parameter has a non-zero value at the encounter with the singularity, therefore the dust implies further expansion. With continued expansion however, the energy density and the pressure of the anti-Chaplygin gas would become ill-defined, hence from the point of view of the anti-Chaplygin gas only a contraction is allowed. Paradoxically, the universe in this cosmological model would have to expand and contract simultaneously. This obviosly could not happen. We solve the paradox by redefining the anti-Chaplygin gas in a distributional sense. Then a contraction could follow the expansion phase at the singularity at the price of a jump in the Hubble parameter. Although such an abrupt change is not common in any cosmological evolution, we explicitly show that the set of Friedmann, Raychaudhuri and continuity equations are all obeyed both at the singularity and in its vicinity. We also prove that the Israel junction conditions are obeyed through the singular spatial hypersurface. In particular we enounce and prove a more general form of the Lanczos equation.

We study a relation between the cosmological singularities in classical and quantum theory, comparing the classical and quantum dynamics in some models possessing the Big Brake singularity - the model based on a scalar field and two models based on a tachyon-pseudo-tachyon field . It is shown that the effect of quantum avoidance is absent for the soft singularities of the Big Brake type while it is present for the Big Bang and Big Crunch singularities. Thus, there is some kind of a classical - quantum correspondence, because soft singularities are traversable in classical cosmology, while the strong Big Bang and Big Crunch singularities are not traversable.

We develop a technique for the reconstruction of the potential for a scalar field in cosmological models based on induced gravity. The potentials reproducing cosmological evolutions driven by barotropic perfect fluids, a cosmological constant, a Chaplygin gas and a modified Chaplygin gas are constructed explicitly.

We study the cosmological evolution of an induced gravity model with a scale symmetry breaking potential for the scalar field. The radiation to matter transition, following inflation and reheating, influences the dynamics of such a field through its non minimal coupling. We illustrate how, under certain conditions on the potential, such a dynamics can lead to a suitable amount of dark energy explaining the present accelerated expansion.

We expand previous works on the spherical `top-hat' collapse (SC-TH) framework in generalized Chaplygin gas (gCg) dominated universes. Here we allow the collapse in all energetic components within the system. We analyze the non-linear stages of collapse for various choices of parameter \alpha of the gCg model introducing an exact formulation for the so-called effective sound speed, $c_{eff}^2$. We show that, within the SC-TH framework, the growth of the structure becomes faster with increasing values of \alpha.

We investigate a Friedmann universe filled with a tachyon scalar field, which behaves as dustlike matter in the past, while it is able to accelerate the expansion rate of the universe at late times. The comparison with type Ia supernovae (SNIa) data allows for evolutions driving the universe into a Big Brake. Some of the evolutions leading to a Big Brake exhibit a large variation of the equation of state parameter at low redshifts which is potentially observable with future data though hardly detectable with present SNIa data. The soft Big Brake singularity occurs at finite values of the scale factor, vanishing energy density and Hubble parameter, but diverging deceleration and infinite pressure. We show that the geodesics can be continued through the Big Brake and that our model universe will recollapse eventually in a Big Crunch. Although the time to the Big Brake strongly depends on the present values of the tachyonic field and of its time derivative, the time from the Big Brake to the Big Crunch represents a kind of invariant timescale for all field parameters allowed by SNIa.

An exact solution describing the evolution of the type Bang-to-Rip with the phantom divide line crossing is constructed in the Chameleon cosmology model, based on two independent functions of the scalar field.

We investigate the Tolman-Oppenheimer-Volkoff equations for the generalized Chaplygin gas with the aim of extending the findings of V. Gorini, U. Moschella, A. Y. Kamenshchik, V. Pasquier, and A. A. Starobinsky [Phys. Rev. D {\bf 78}, 064064 (2008)]. We study both the standard case, where we reproduce some previous results, and the phantom case. In the phantom case we show that even a superluminal group velocity arising for $\alpha > 1$ cannot prevent the divergence of the pressure at a finite radial distance. Finally, we investigate how a modification of the generalized Chaplygin gas equation of state, required by causality arguments at densities very close to $\Lambda$, affects the results found so far.

We consider the PT symmetric flat Friedmann model of two scalar fields with positive kinetic terms. While the potential of one ("normal") field is taken real, that of the other field is complex. We study a complex classical solution of the system of the two Klein-Gordon equations together with the Friedmann equation. The solution for the normal field is real while the solution for the second field is purely imaginary, realizing classically the "phantom" behavior. The energy density and pressure are real and the corresponding geometry is well-defined. The Lagrangian for the linear perturbations has the correct potential signs for both the fields, so that the problem of stability does not arise. The background dynamics is determined by an effective action including two real fields one normal and one "phantom". Remarkably, the phantom phase in the cosmological evolution is transient and the Big Rip never occurs. Our model is contrasted to well-known quintom models, which also include one normal and one phantom fields.

We investigate whether a tachyonic scalar field, encompassing both dark energy and dark matter-like features will drive our universe towards a Big Brake singularity or a de Sitter expansion. In doing this it is crucial to establish the parameter domain of the model, which is compatible with type Ia supernovae data. We find the 1-sigma contours and evolve the tachyonic sytem into the future. We conclude, that both future evolutions are allowed by observations, Big Brake becoming increasingly likely with the increase of the positive model parameter k.

We consider the renormalization group improvement in the theory of the Standard Model Higgs boson playing the role of an inflaton with a strong non-minimal coupling to gravity. It suggests the range of the Higgs mass $135.6 {\rm GeV} \lesssim M_H\lesssim 184.5 {\rm GeV}$ compatible with the current CMB data (the lower WMAP bound on $n_s$), which is close to the widely accepted range dictated by the electroweak vacuum stability and perturbation theory bounds. We find the phenomenon of asymptotic freedom induced by this non-minimal curvature coupling, which brings the theory to the weak coupling domain everywhere except at the lower and upper boundary of this range. The renormalization group running of the basic quantity ${\boldmath $A_I$}$ -- the anomalous scaling in the non-minimally coupled Standard Model, which analytically determines all characteristics of the CMB spectrum -- brings ${\boldmath $A_I$}$ to small negative values at the inflation scale. This property is crucial for the above results and may also underlie the formation of initial conditions for the inflationary dynamics in quantum cosmology.

We compare the existing observational data on type Ia Supernovae with the evolutions of the universe predicted by a one-parameter family of tachyon models which we have introduced recently in paper \cite{we-tach}. Among the set of the trajectories of the model which are compatible with the data there is a consistent subset for which the universe ends up in a new type of soft cosmological singularity dubbed Big Brake. This opens up yet another scenario for the future history of the universe besides the one predicted by the standard $\Lambda$CDM model.

We study static solutions of the Tolman--Oppenheimer--Volkoff equations for spherically symmetric objects (stars) living in a space filled with the Chaplygin gas. Two cases are considered. In the normal case all solutions (excluding the de Sitter one) realize a three-dimensional spheroidal geometry because the radial coordinate achieves a maximal value (the "equator"). After crossing the equator, three scenarios are possible: a closed spheroid having a Schwarzschild-type singularity with infinite blue-shift at the "south pole", a regular spheroid, and a truncated spheroid having a scalar curvature singularity at a finite value of the radial coordinate. The second case arises when the modulus of the pressure exceeds the energy density (the phantom Chaplygin gas). There is no more equator and all solutions have the geometry of a truncated spheroid with the same type of singularity. We consider also static spherically symmetric configurations existing in a universe filled with the phantom Chaplygin gas only. In this case two classes of solutions exist: truncated spheroids and solutions of the wormhole type with a throat. However, the latter are not asymptotically flat and possess curvature singularities at finite values of the radial coordinate. Thus, they may not be used as models of observable compact astrophysical objects.

We consider a quantum corrected inflation scenario driven by a generic GUT or Standard Model type particle model whose scalar field playing the role of an inflaton has a strong non-minimal coupling to gravity. We show that currently widely accepted bounds on the Higgs mass falsify the suggestion of the paper arXiv:0710.3755 (where the role of radiative corrections was underestimated) that the Standard Model Higgs boson can serve as the inflaton. However, if the Higgs mass could be raised to $\sim 230$ GeV, then the Standard Model could generate an inflationary scenario with the spectral index of the primordial perturbation spectrum $n_s\simeq 0.935$ (barely matching present observational data) and the very low tensor-to-scalar perturbation ratio $r\simeq 0.0006$.

We consider two different toy cosmological models based on two fields (one normal scalar and one phantom) realizing the same evolution of the Bang-to-Rip type. One of the fields (pseudoscalar) interacts with the magnetic field breaking the conformal invariance of the latter. The effects of the amplification of cosmic magnetic fields are studied and it is shown that the presence of such effects can discriminate between different cosmological models realizing the same global evolution of the universe.

We study the procedure of the reconstruction of phantom-scalar field potentials in two-field cosmological models. It is shown that while in the one-field case the chosen cosmological evolution defines uniquely the form of the scalar potential, in the two-field case one has an infinite number of possibilities. The classification of a large class of possible potentials is presented and the dependence of cosmological dynamics on the choice of initial conditions is investigated qualitatively and numerically for two particular models.

We exploit the gauge-invariant formalism to analyse the perturbative behaviour of two cosmological models based on the generalized Chaplygin gas describing both dark matter and dark energy in the present Universe. In the first model we consider the generalized Chaplygin gas alone, while in the second one we add a baryon component to it. We extend our analysis also into the parameter range $\alpha > 1$, where the generalized Chaplygin gas sound velocity can be larger than that of light. In the first model we find that the matter power spectrum is compatible with the observed one only for $\alpha < 10^{-5}$, which makes the generalized Chaplygin gas practically indistinguishable from $\Lambda$CDM. In the second model we study the evolution of inhomogeneities of the baryon component. The theoretical power spectrum is in good agreement with the observed one for almost all values of $\alpha$. However, the growth of inhomogeneities seems to be particularly favoured either for sufficiently small values of $\alpha$ or for $\alpha \gtrsim 3$. Thus, it appears that the viability of the generalized Chaplygin gas as a cosmological model is stronger when its sound velocity is superluminal. We show that in this case the generalized Chaplygin gas equation of state can be changed in an unobservable region in such a way that its equivalent $k$-essence microscopical model has no problems with causality.

We construct a scalar field based cosmological model, possessing a cosmological singularity characterized by a finite value of the cosmological radius and an infinite scalar curvature. Using the methods of the qualitative theory of differential equations, we give a complete description of the cosmological evolutions in the model under consideration. There are four classes of evolutions, two of which have finite lifetimes, while the other two undergo an infinite expansion.

We outline the evaluation of the cosmological constant in the framework of the standard field-theoretical treatment of vacuum energy and discuss the relation between the vacuum energy problem and the gauge-group spontaneous symmetry breaking. We suggest possible extensions of the 't Hooft-Nobbenhuis symmetry, in particular, its complexification till duality symmetry and discuss the compatible implementation on gravity. We propose to use the discrete time-reflection transform to formulate a framework in which one can eliminate the huge contributions of vacuum energy into the effective cosmological constant and suggest that the breaking of time--reflection symmetry could be responsible for a small observable value of this constant.

Motivated by the generalization of quantum theory for the case of non-Hermitian Hamiltonians with PT symmetry, we show how a classical cosmological model describes a smooth transition from ordinary dark energy to the phantom one. The model is based on a classical complex Lagrangian of a scalar field. Specific symmetry properties analogous to PT in non-Hermitian quantum mechanics lead to purely real equation of motion.

We discuss the phenomenon of the smooth dynamical gravity induced crossing of the phantom divide line in a framework of simple cosmological models where it appears to occur rather naturally, provided the potential of the unique scalar field has some kind of cusp. The behavior of cosmological trajectories in the vicinity of the cusp is studied in some detail and a simple mechanical analogy is presented. The phenomenon of certain complementarity between the smoothness of the spacetime geometry and matter equations of motion is elucidated. We introduce a network of cosmological histories and qualitatively describe some of its properties.

We apply the method of comparison equations to study cosmological perturbations during inflation, obtaining the full power spectra of scalar and tensor perturbations to first and to second order in the slow-roll parameters. We compare our results with those derived by means of other methods, in particular the Green's function method and the improved WKB approximation, and find agreement for the slow-roll structure. The method of comparison equations, just as the improved WKB approximation, can however be applied to more general situations where the slow-roll approximation fails.

Inspired by the generalization of quantum theory for the case of non-Hermitian Hamiltonians with CPT symmetry, we construct a simple classical cosmological scalar field based model describing a smooth transition from ordinary dark energy to the phantom one.

Flat FRW perfect fluid cosmologies can be reproduced as particular solutions of suitable field theoretical models. Here we investigate the stability of perfect fluid model trajectories with respect to sets of trajectories of the corresponding field models having generic initial conditions. It is shown that the trajectories of barotropic perfect fluid models and those of the Chaplygin gas model are stable. The total probability to reach the Chaplygin gas regime early enough to achieve a matter dominated stage having a realistic duration is calculated for a scalar field model assuming a natural measure in the space of initial conditions taken near a cosmological singularity. An example is presented of a two-fluid cosmological model where the stability is partially absent.

Simple scalar field cosmological models are considered describing gravity assisted crossing of the phantom divide line. This crossing or (de)-phantomization characterized by the change of the sign of the kinetic term of the scalar field is smooth and driven dynamically by the Einstein equations. Different cosmological scenarios, including the phantom phase of matter are sketched.

We derive a relativistic hodograph equation for a two-dimensional stationary isentropic hydrodynamical motion. For the case of stiff matter, when the velocity of sound coincides with the light speed, the singularity in this equation disappears and the solutions become regular in all hodograph plane.

We review the essential features of the Chaplygin gas cosmological models and provide some examples of appearance of the Chaplygin gas equation of state in modern physics. A possible theoretical basis for the Chaplygin gas in cosmology is discussed. The relation with scalar field and tachyon cosmological models is also considered.

We study the role that tachyon fields may play in cosmology as compared to the well-established use of minimally coupled scalar fields. We first elaborate on a kind of correspondence existing between tachyons and minimally coupled scalar fields; corresponding theories give rise to the same cosmological evolution for a particular choice of the initial conditions but not for any other. This leads us to study a specific one-parameter family of tachyonic models based on a perfect fluid mixed with a positive cosmological constant. For positive values of the parameter one needs to modify Sen's action and use the sigma process of resolution of singularities. The physics described by this model is dramatically different and much richer than that of the corresponding scalar field. For particular choices of the initial conditions the universe, that does mimick for a long time a de Sitter-like expansion, ends up in a finite time in a special type of singularity that we call a "big brake". This singularity is characterized by an infinite deceleration.

One of the most interesting features of braneworld models is the existence of massive gravitational modes in addition to the usual massless one. Mixing of the modes which depends nontrivially on the inter-brane distance can be interpreted as radion-induced gravitational-wave oscillations, a classical analogy to meson and neutrino oscillations. We show that these oscillations arising in M-theory inspired two-brane setups could lead to effects detectable by gravitational-wave interferometers.

We discuss the theory and phenomenology of the interplay between the massless graviton and its massive Kaluza-Klein modes in the Randall-Sundrum two-brane model. The equations of motion of the transverse traceless degrees of freedom are derived by means of a Green function approach as well as from an effective nonlocal action. The second procedure clarifies the extraction of the particle content from the nonlocal action and the issue of its diagonalization. The situation discussed is generic for the treatment of two-brane models if the on-brane fields are used as the dynamical degrees of freedom. The mixing of the effective graviton modes of the localized action can be interpreted as radion-induced gravitational-wave oscillations, a classical analogy to meson and neutrino oscillations. We show that these oscillations arising in M-theory-motivated braneworld setups could lead to effects detectable by gravitational-wave interferometers. The implications of this effect for models with ultra-light gravitons are discussed.

The quasi-isotropic inhomogeneous solution of the Einstein equations near a cosmological singularity in the form of a series expansion in the synchronous system of reference, first found by Lifshitz and Khalatnikov in 1960, is generalized to the case of a two-fluid cosmological model. This solution describes non-decreasing modes of adiabatic and isocurvature scalar perturbations and gravitational waves in the regime when deviations of a space-time metric from the homogeneous isotropic Friedmann-Robertson-Walker (FRW) background are large while locally measurable quantities like Riemann tensor components are still close to their FRW values. The general structure of the perturbation series is presented and the first coefficients of the series expansion for the metric tensor and the fluid energy densities and velocities are calculated explicitly.

In this note two cosmological models representing the flat Friedmann Universe filled with a Chaplygin fluid, with or without dust, are analyzed in terms of the recently proposed "statefinder" parameters. Trajectories of both models in the parameter plane are shown to be significantly different w.r.t. "quiessence" and "tracker" models. The generalized Chaplygin gas model with an equation of state of the form $p = -A/\rho^{\alpha}$ is also analyzed in terms of the statefinder parameters.