The current Universe is composed by a mixture of relativistic species, baryonic matter, dark matter and dark energy which evolve in a non-trivial way at perturbative level. An advanced description of the cosmological dynamics should include non-standard features beyond the simplistic approach idealized by the standard cosmology in which cosmic components do not interact, are adiabatic and dissipationless. We promote a full perturbative analysis of linear scalar perturbations of a cosmological model containing baryons, dark matter and a scalar field allowing for the presence of relative entropic perturbations between the three fluids. By including the latter feature we establish a new set of equations for the scalar cosmological perturbations. As a consequence of this new approach, we show that tiny departures from a constant scalar field equation of state $w_{S}=-1$ damage structure formation in a non-acceptable manner. Hence, by strongly constraining $w_S$ our results provide compelling evidence in favor of the standard cosmological model and rule out a large class of dynamical dark energy models.

Models with non-gravitational interactions between the dark matter and dark energy components are an alternative to the standard cosmological scenario. These models are characterized by an interaction term, and a frequently used parameterization is $Q = 3\xi H \rho_{x}$, where $H$ is the Hubble parameter and $\rho_{x}$ is the dark energy density. Although current observations support such a model for negative values of the interaction parameter $\xi$, we show here that this interval of values of $\xi$ leads the model to predict a violation of the Weak Energy Condition (WEC) for the dark matter density, regardless of the value of the equation-of-state parameter of the dark energy component. This violation is accompanied by unphysical instabilities of matter perturbations.

On the basis of a previously established scalar-tensor extension of the $\Lambda$CDM model we develop an effective fluid approach for the matter growth function. This extended $\Lambda$CDM (henceforth $e_{\Phi}\Lambda$CDM) cosmology takes into account deviations from the standard model both via a modified background expansion and by the inclusion of geometric anisotropic stresses as well as of perturbations of the geometric dark-energy equivalent. The background dynamics is governed by an explicit analytic expression for the Hubble rate in which modifications of the standard model are given in terms of a single constant parameter [1]. To close the system of fluid-dynamical perturbation equations we introduce two phenomenological parameters through which the anisotropic stress is related both to the total energy density perturbation of the cosmic substratum and to relative perturbations in the effective two-component system. We quantify the impact of deviations from the standard background, of anisotropic stresses and of non-vanishing perturbations of the effective dark-energy component on the matter growth rate function $f \sigma_8$ and confront the results with recent redshift-space distortion (RSD) measurements.

Extensions of Einstein's General Relativity (GR) can formally be given a GR structure in which additional geometric degrees of freedom are mapped on an effective energy-momentum tensor. The corresponding effective cosmic medium can then be modeled as an imperfect fluid within GR. The imperfect fluid structure allows us to include, on a phenomenological basis, anisotropic stresses and energy fluxes which are considered as potential signatures for deviations from the cosmological standard $\Lambda$-cold-dark-matter ($\Lambda$CDM) model. As an example, we consider the dynamics of a scalar-tensor extension of the standard model, the $e_{\Phi}\Lambda$CDM model. We constrain the magnitudes of anisotropic pressure and energy flux with the help of redshift-space distortion (RSD) data for the matter growth function $f \sigma_8$.

We reconsider the dynamics of the Universe in the presence of interactions in the cosmological dark sector. A class of interacting models is introduced via a real function $f\left(r\right)$ of the ratio $r$ between the energy densities of the (pressureless) cold dark matter (CDM) and dark energy (DE). The subclass of models for which the ratio $r$ depends only on the scale factor is shown to be equivalent to unified models of the dark sector, i.e. models for which the CDM and DE components can be combined in order to form a unified dark fluid. For specific choices of the function $f\left(r\right)$ we recover several models already studied in the literature. We analyse various special cases of this type of interacting models using a suitably modified version of the CLASS code combined with MontePython in order to constrain the parameter space with the data from supernova of type SNe Ia (JLA), the Hubble constant $H_{0}$, cosmic chronometers (CC), baryon acoustic oscilations (BAO) and data from the Planck satellite (Planck TT). Our analysis shows that even if data from the late Universe ($H_{0}$, SNe Ia and CC) indicate an interaction in the dark sector, the data related to the early Universe (BAO and Planck TT) constrain this interaction substantially, in particular for cases in which the background dynamics is strongly affected.

Both scalar fields and (generalized) Chaplygin gases have been widely used separately to characterize the dark sector of the Universe. Here we investigate the cosmological background dynamics for a mixture of both these components and quantify the fractional abundances that are admitted by observational data from supernovae of type Ia and from the evolution of the Hubble rate. Moreover, we study how the growth rate of (baryonic) matter perturbations is affected by the dark-sector perturbations.

Yes, but only for a parameter value that makes it almost coincide with the standard model. We reconsider the cosmological dynamics of a generalized Chaplygin gas (gCg) which is split into a cold dark matter (CDM) part and a dark energy (DE) component with constant equation of state. This model, which implies a specific interaction between CDM and DE, has a $\Lambda$CDM limit and provides the basis for studying deviations from the latter. Including matter and radiation, we use the (modified) CLASS code \cite{class} to construct the CMB and matter power spectra in order to search for a gCg-based concordance model that is in agreement with the SNIa data from the JLA sample and with recent Planck data. The results reveal that the gCg parameter $\alpha$ is restricted to $|\alpha|\lesssim 0.05$, i.e., to values very close to the $\Lambda$CDM limit $\alpha =0$. This excludes, in particular, models in which DE decays linearly with the Hubble rate.

An interaction between dark matter and dark energy, proportional to the product of their energy densities, results in a scaling behavior of the ratio of these densities with respect to the scale factor of the Robertson-Walker metric. This gives rise to a class of cosmological models which deviate from the standard model in an analytically tractable way. In particular, it becomes possible to quantify the role of potential dark-energy perturbations. We investigate the impact of this interaction on the structure formation process. Using the (modified) CAMB code we obtain the CMB spectrum as well as the linear matter power spectrum. It is shown that the strong degeneracy in the parameter space present in the background analysis is considerably reduced by considering \textit{Planck} data. Our analysis is compatible with the $\Lambda$CDM model at the $2\sigma$ confidence level with a slightly preferred direction of the energy flow from dark matter to dark energy.

We consider cosmological backreaction effects in Buchert's averaging formalism on the basis of an explicit solution of the Lema\^itre-Tolman-Bondi (LTB) dynamics which is linear in the LTB curvature parameter and has an inhomogeneous bang time. The volume Hubble rate is found in terms of the volume scale factor which represents a derivation of the simplest phenomenological solution of Buchert's equations in which the fractional densities corresponding to average curvature and kinematic backreaction are explicitly determined by the parameters of the underlying LTB solution at the boundary of the averaging volume. This configuration represents an exactly solvable toy model but it does not adequately describe our "real" Universe.

We construct a cosmological scalar-tensor-theory model in which the Brans-Dicke type scalar $\Phi$ enters the effective (Jordan-frame) Hubble rate as a simple modification of the Hubble rate of the $\Lambda$CDM model. This allows us to quantify differences between the background dynamics of scalar-tensor theories and general relativity (GR) in a transparent and observationally testable manner in terms of one single parameter. Problems of the mapping of the scalar-field degrees of freedom on an effective fluid description in a GR context are discused. Data from supernovae, the differential age of old galaxies and baryon acoustic oscillations are shown to strongly limit potential deviations from the standard model.

Based on the Lema\^itre-Tolman-Bondi (LTB) metric we consider two flat inhomogeneous big-bang models. We aim at clarifying, as far as possible analytically, basic features of the dynamics of the simplest inhomogeneous models and to point out the potential usefulness of exact inhomogeneous solutions as generalizations of the homogeneous configurations of the cosmological standard model. We discuss explicitly partial successes but also potential pitfalls of these simplest models. Although primarily seen as toy models, the relevant free parameters are fixed by best-fit values using the Joint Light-curve Analysis (JLA)-sample data. On the basis of a likelihood analysis we find that a local hump with an extension of almost 2 Gpc provides a better description of the observations than a local void for which we obtain a best-fit scale of about 30 Mpc. Future redshift-drift measurements are discussed as a promising tool to discriminate between inhomogeneous configurations and the $\Lambda$CDM model.

We present a new approximation to include fully general relativistic pressure and velocity in Newtonian hydrodynamics. The energy conservation, momentum conservation and two Poisson's equations are consistently derived from Einstein's gravity in the zero-shear gauge assuming weak gravity and action-at-a-distance limit. The equations show proper special relativity limit in the absence of gravity. Our approximation is complementary to the post-Newtonian approximation and the equations are valid in fully nonlinear situations.

We consider an effective viscous pressure as the result of a backreaction of inhomogeneities within Buchert's formalism. The use of an effective metric with a time-dependent curvature radius allows us to calculate the luminosity distance of the backreaction model. This quantity is different from its counterpart for a "conventional" spatially flat bulk viscous fluid universe. Both expressions are tested against the SNIa data of the Union2.1 sample with only marginally different results for the distance-redshift relation and in accordance with the $\Lambda$CDM model. Future observations are expected to be able to discriminate among these models on the basis of indirect measurements of the curvature evolution.

A suitable nonlinear interaction between dark matter with an energy density $\rho_{M}$ and dark energy with an energy density $\rho_{X}$ is known to give rise to a non-canonical scaling $\rho_{M} \propto \rho_{X}a^{-\xi}$ where $\xi$ is a parameter which generally deviates from $\xi =3$. Here we present a covariant generalization of this class of models and investigate the corresponding perturbation dynamics. The resulting matter power spectrum for the special case of a time-varying Lambda model is compared with data from the SDSS DR9 catalogue. We find a best-fit value of $\xi = 3.25$ which corresponds to a decay of dark matter into the cosmological term. Our results are compatible with the $\Lambda$CDM model at the 2$\sigma$ confidence level.

We discuss the possibility to implement a viscous cosmological model, attributing to the dark matter component a behaviour described by bulk viscosity. Since bulk viscosity implies negative pressure, this rises the possibility to unify the dark sector. At the same time, the presence of dissipative effects may alleviate the so called small scale problems in the $\Lambda$CDM model. While the unified viscous description for the dark sector does not lead to consistent results, the non-linear behaviour indeed improves the situation with respect to the standard cosmological model.

We investigate a Boltzmann gas at equilibrium with its center of mass moving on a circular geodesics in the Schwarzschild field. As a consequence of Tolman's law we find that a central comoving observer measures oscillations of the temperature and of other thermodynamic quantities with twice the frequencies that are known from test-particle motion. We apply this scheme to the gas dynamics in the gravitational fields of the planets of the Solar System as well as to strong-field configurations of neutron stars and black holes.

The observational fact that the present values of the densities of dark energy and dark matter are of the same order of magnitude, $\rho_{de0}/\rho_{dm0} \sim \mathcal{O}(1)$, seems to indicate that we are currently living in a very special period of the cosmic history. Within the standard model, a density ratio of the order of one just at the present epoch can be seen as coincidental since it requires very special initial conditions in the early Universe. The corresponding "why now" question constitutes the cosmological "coincidence problem". According to the standard model the equality $\rho_{de} = \rho_{dm}$ took place "recently" at a redshift $z \approx 0.55$. The meaning of "recently" is, however, parameter dependent. In terms of the cosmic time the situation looks different. We discuss several aspects of the "coincidence problem", also in its relation to the cosmological constant problem, to issues of structure formation and to cosmic age considerations.

We consider the perturbation dynamics for the cosmic baryon fluid and determine the corresponding power spectrum for a $\Lambda(t)$CDM model in which a cosmological term decays into dark matter linearly with the Hubble rate. The model is tested by a joint analysis of data from supernovae of type Ia (SNIa) (Constitution and Union 2.1), baryonic acoustic oscillation (BAO), the position of the first peak of the anisotropy spectrum of the cosmic microwave background (CMB) and large-scale-structure (LSS) data (SDSS DR7). While the homogeneous and isotropic background dynamics is only marginally influenced by the baryons, there are modifications on the perturbative level if a separately conserved baryon fluid is included. Considering the present baryon fraction as a free parameter, we reproduce the observed abundance of the order of $5\%$ independently of the dark-matter abundance which is of the order of $32\%$ for this model. Generally, the concordance between background and perturbation dynamics is improved if baryons are explicitly taken into account.

Interactions inside the cosmological dark sector influence the cosmological dynamics. As a consequence, the future evolution of the Universe may be different from that predicted by the $\Lambda$CDM model. We review main features of several recently studied models with nongravitational couplings between dark matter and dark energy.

We consider the dynamics of a cosmological substratum of pressureless matter and holographic dark energy with a cutoff length proportional to the Ricci scale. Stability requirements for the matter perturbations are shown to single out a model with a fixed relation between the present matter fraction $\Omega_{m0}$ and the present value $\omega_{0}$ of the equation-of-state parameter of the dark energy. This model has the same number of free parameters as the $\Lambda$CDM model but it has no $\Lambda$CDM limit. We discuss the consistency between background observations and the mentioned stability-guaranteeing parameter combination.

The Reduced Relativistic Gas (RRG) model was introduced by A. Sakharov in 1965 for deriving the cosmic microwave background (CMB) spectrum. It was recently reinvented by some of us to achieve an interpolation between the radiation and dust epochs in the evolution of the Universe. This model circumvents the complicated structure of the Boltzmann-Einstein system of equations and admits a transparent description of warm-dark-matter effects. It is extended here to include, on a phenomenological basis, an out-of-equilibrium interaction between radiation and baryons which is supposed to account for relevant aspects of pre-recombination physics in a simplified manner. Furthermore, we use the tight-coupling approximation to explore the influence of both this interaction and of the RRG warmness parameter on the anisotropy spectrum of the CMB. The predictions of the model are very similar to those of the {\Lambda}CDM model if both the interaction and the dark-matter warmness parameters are of the order of $10^{-4}$ or smaller. As far as the warmness parameter is concerned, this is in good agreement with previous estimations on the basis of results from structure formation.

We assume cold dark matter to possess a small bulk-viscous pressure which typically attenuates the growth of inhomogeneities. Explicit calculations, based on Eckart's theory of dissipative processes, reveal that for viscous cold dark matter the usual Newtonian approximation for perturbation scales smaller than the Hubble scale is no longer valid. We advocate the use of a neo-Newtonian approach which consistently incorporates pressure effects into the fluid dynamics and correctly reproduces the general relativistic dynamics. This result is of interest for numerical simulations of nonlinear structure formation involving nonstandard dark-matter fluids. We obtain upper limits on the magnitude of the viscous pressure by requiring that relevant perturbation amplitudes should grow sufficiently to enter the nonlinear stage.

The split of a generalised Chaplygin gas with an equation of state p = -A/\rho^{\alpha} into an interacting mixture of pressureless matter and a dark-energy component with equation of state p_{\Lambda} = - \rho_{\Lambda} implies the existence of non-adiabatic pressure perturbations. We demonstrate that the square of the effective (non-adiabatic) sound speed c_s of the medium is proportional to the ratio of the perturbations of the dark energy to those of the dark matter. Since, as demonstrated explicitly for the particular case \alpha = -1/2, dark-energy perturbations are negligible compared with dark-matter perturbations on scales that are relevant for structure formation, we find |c_s^2| << 1. Consequently, there are no oscillations or instabilities which have plagued previous adiabatic Chaplygin-gas models.

We assume the cosmological dark sector to consist of pressureless matter and holographic dark energy with a cutoff length proportional to the Ricci scale. The requirement of separate energy-momentum conservation of the components is shown to establish a relation between the matter fraction and the (necessarily time-dependent) equation-of-state parameter of the dark energy. Focusing on intrinsically adiabatic pressure perturbations of the dark-energy component, the matter perturbations are found as linear combinations of the total energy-density perturbations of the cosmic medium and the relative (nonadiabatic) perturbations of the components. The resulting background dynamics is consistent with observations from supernovae of type Ia, baryonic acoustic oscillations and the differential age of old objects. The perturbation dynamics, on the other hand, is plagued by instabilities which excludes any phantom-type equation of state. The only stable configuration is singled out by a fixed relation between the present matter fraction $\Omega_{m0}$ and the present value $\omega_{0}$ of the equation-of-state parameter of the dark energy. However, this instability-avoiding configuration is only marginally consistent with the observationally preferred background values of the mentioned parameters.

We argue that interactions in the dark sector may have a crucial impact on the cosmological dynamics. In particular, the future cosmic evolution may be very different from that predicted by the LCDM model. An example is a scenario in which the currently observed accelerated expansion is an interaction-induced transient phenomenon. We discuss such type of behavior on the basis of a two-fluid toy model.

We consider a set of non-linear interactions between dark matter and dark energy which comprises couplings proportional to products of (powers of) the energy densities of both dark components and of the total energy. We demonstrate that under such conditions the final state of the Universe may differ substantially from that of the standard LCDM model. In particular, the ratio of the energy densities of dark matter and dark energy may approach a stable finite value or oscillate about such a value. Stationary solutions of this type require a phantom-type "bare" equation of state of the dark energy which, however, does not lead to a big-rip singularity. A corresponding analytic solution for a particular case is shown to be consistent with the supernova type Ia (SNIa) data from the Union2 set. For a broader class of interactions, a dynamical system analysis classifies stationary points with emphasis on their potential relevance for the coincidence problem.

We demonstrate that creation of dark-matter particles at a constant rate implies the existence of a cosmological term that decays linearly with the Hubble rate. We discuss the cosmological model that arises in this context and test it against observations of the first acoustic peak in the cosmic microwave background (CMB) anisotropy spectrum, the Hubble diagram for supernovas of type Ia (SNIa), the distance scale of baryonic acoustic oscillations (BAO) and the distribution of large scale structures (LSS). We show that a good concordance is obtained, albeit with a higher value of the present matter abundance than in the \Lambda CDM model. We also comment on general features of the CMB anisotropy spectrum and on the cosmic coincidence problem.

Interactions between dark matter and dark energy which result in a power-law behavior (with respect to the cosmic scale factor) of the ratio between the energy densities of the dark components (thus generalizing the LCDM model) have been considered as an attempt to alleviate the cosmic coincidence problem phenomenologically. We generalize this approach by allowing for a variable equation of state for the dark energy within the CPL-parametrization. Based on analytic solutions for the Hubble rate and using the Constitution and Union2 SNIa sets, we present a statistical analysis and classify different interacting and non-interacting models according to the Akaike (AIC) and the Bayesian (BIC) information criteria. We do not find noticeable evidence for an alleviation of the coincidence problem with the mentioned type of interaction.

Any non-gravitational coupling between dark matter and dark energy modifies the cosmological dynamics. Interactions in the dark sector are considered to be relevant to address the coincidence problem. Moreover, in various models the observed accelerated expansion of the Universe is a pure interaction phenomenon. Here we review recent approaches in which a coupling between both dark components is crucial for the evolution of the Universe.

According to the standard $\Lambda$CDM model, the accelerated expansion of the Universe will go on forever. Motivated by recent observational results, we explore the possibility of a finite phase of acceleration which asymptotically approaches another period of decelerated expansion. Extending an earlier study on a corresponding homogeneous and isotropic dynamics, in which interactions between dark matter and dark energy are crucial, the present paper also investigates the dynamics of the matter perturbations both on the Newtonian and GR levels and quantifies the potential relevance of perturbations of the dark-energy component. In the background, the model is tested against the Supernova type Ia (SNIa) data of the Constitution set and on the perturbative level against growth rate data and the data of the 2dFGRS project. Our results indicate, that a transient phase of accelerated expansion is not excluded by current observations.

The Universe is modeled as consisting of pressureless baryonic matter and a bulk viscous fluid which is supposed to represent a unified description of the dark sector. In the homogeneous and isotropic background the \textit{total} energy density of this mixture behaves as a generalized Chaplygin gas. The perturbations of this energy density are intrinsically nonadiabatic and source relative entropy perturbations. The resulting baryonic matter power spectrum is shown to be compatible with the 2dFGRS and SDSS (DR7) data. A joint statistical analysis, using also Hubble-function and supernovae Ia data, shows that, different from other studies, there exists a maximum in the probability distribution for a negative present value of the deceleration parameter. Moreover, the unified model presented here favors a matter content that is of the order of the baryonic matter abundance suggested by big-bang nucleosynthesis. A problem of simple bulk viscous models, however, is the behavior of the gravitational potential and the reproduction of the CMB power spectrum.

Different holographic dark-energy models are studied from a unifying point of view. We compare models for which the Hubble scale, the future event horizon or a quantity proportional to the Ricci scale are taken as the infrared cutoff length. We demonstrate that the mere definition of the holographic dark-energy density generally implies an interaction with the dark-matter component. We discuss the relation between the equation-of-state parameter and the energy density ratio of both components for each of the choices, as well as the possibility of non-interacting and scaling solutions. Parameter estimations for all three cutoff options are performed with the help of a Bayesian statistical analysis, using data from supernovae type Ia and the history of the Hubble parameter. The $\Lambda$CDM model is the clear winner of the analysis. According to the Bayesian Information Criterion ($BIC$), all holographic models should be considered as ruled out, since the difference $\Delta BIC$ to the corresponding $\Lambda$CDM value is $> 10$. According to the Akaike Information Criterion ($AIC$), however, we find $\Delta AIC$ $< 2$ for models with Hubble-scale and Ricci-scale cutoffs, indicating, that they may still be competitive. As we show for the example of the Ricci-scale case, also the use of certain priors, reducing the number of free parameters to that of the $\Lambda$CDM model, may result in a competitive holographic model.

We investigate a spatially flat Friedmann-Lema\^itre-Robertson-Walker cosmology in which a decaying vacuum term causes matter production at late times. Assuming a decay proportional to the Hubble rate, the ratio of the background energy densities of dark matter and dark energy changes with the cosmic scale factor as $a^{-3/2}$. The intrinsically non-adiabatic two-component perturbation dynamics of this model is reduced to a single second-order equation. Perturbations of the vacuum term are shown to be negligible on scales that are relevant for structure formation. On larger scales, dark-energy perturbations give a somewhat higher contribution but remain always smaller than the dark-matter perturbations.

We consider cosmological scenarios originating from a single imperfect fluid with bulk viscosity and apply Eckart's and both the full and the truncated M\"uller-Israel-Stewart's theories as descriptions of the non-equilibrium processes. Our principal objective is to investigate if the dynamical properties of Dark Matter and Dark Energy can be described by a single viscous fluid and how such description changes when a causal theory (M\"uller-Israel-Stewart's, both in its full and truncated forms) is taken into account instead of Eckart's non-causal theory. To this purpose, we find numerical solutions for the gravitational potential and compare its behaviour with the corresponding LambdaCDM case. Eckart's and the full causal theory seem to be disfavoured, whereas the truncated theory leads to results similar to those of the LambdaCDM model for a bulk viscous speed in the interval 10^{-11} << c_b^2 < 10^{-8}. Tentatively relating such value to a square propagation velocity of the order of T/m of perturbations in a non-relativistic gas of particles with mass m at the epoch of matter-radiation equality, this may be compatible with a mass range 0.1 GeV < m << 100 GeV.

We model the cosmological substratum by a viscous fluid that is supposed to provide a unified description of the dark sector and pressureless baryonic matter. In the homogeneous and isotropic background the \textit{total} energy density of this mixture behaves as a generalized Chaplygin gas. The perturbations of this energy density are intrinsically non-adiabatic and source relative entropy perturbations. The resulting baryonic matter power spectrum is shown to be compatible with the 2dFGRS and SDSS (DR7) data. A joint statistical analysis, using also Hubble-function and supernovae Ia data, shows that, different from other studies, there exists a maximum in the probability distribution for a negative present value $q_{0} \approx - 0.53$ of the deceleration parameter. Moreover, different from other approaches, the unified model presented here favors a matter content that is of the order of the baryonic matter abundance suggested by big-bang nucleosynthesis.

We investigate the cosmological perturbation dynamics for a universe consisting of pressureless baryonic matter and a viscous fluid, the latter representing a unified model of the dark sector. In the homogeneous and isotropic background the \textit{total} energy density of this mixture behaves as a generalized Chaplygin gas. The perturbations of this energy density are intrinsically non-adiabatic and source relative entropy perturbations. The resulting baryonic matter power spectrum is shown to be compatible with the 2dFGRS and SDSS (DR7) data. A joint statistical analysis, using also Hubble-function and supernovae Ia data, shows that, different from other studies, there exists a maximum in the probability distribution for a negative present value $q_0 \approx - 0.53$ of the deceleration parameter. Moreover, while previous descriptions on the basis of generalized Chaplygin gas models were incompatible with the matter power spectrum data since they required a much too large amount of pressureless matter, the unified model presented here favors a matter content that is of the order of the baryonic matter abundance suggested by big-bang nucleosynthesis.

We constrain an interacting, holographic dark energy model, first proposed by two of us in [1], with observational data from supernovae, CMB shift, baryon acoustic oscillations, x-rays, and the Hubble rate. The growth function for this model is also studied. The model fits the data reasonably well but still the conventional $\Lambda$CDM model fares better. Nevertheless, the holographic model greatly alleviates the coincidence problem and shows compatibility at $1\sigma$ confidence level with the age of the old quasar APM 08279+5255.

We find exact power-law solutions for scalar-tensor theories and clarify the conditions under which they can account for an accelerated expansion of the Universe. These solutions have the property that the signs of both the Hubble rate and the deceleration parameter in the Jordan frame may be different from the signs of their Einstein-frame counterparts. For special parameter combinations we identify these solutions with asymptotic attractors that have been obtained in the literature through dynamical-system analysis. We establish an effective general-relativistic description for which the geometrical equivalent of dark energy is associated with a time dependent equation of state. The present value of the latter is consistent with the observed cosmological ``constant". We demonstrate that this type of power-law solutions for accelerated expansion cannot be realized in f(R) theories.

Recent investigations seem to favor a cosmological dynamics according to which the accelerated expansion of the Universe may have already peaked and is now slowing down again \cite{sastaro}. As a consequence, the cosmic acceleration may be a transient phenomenon. We investigate a toy model that reproduces such a background behavior as the result of a time-dependent coupling in the dark sector which implies a cancelation of the "bare" cosmological constant. With the help of a statistical analysis of Supernova Type Ia (SNIa) data we demonstrate that for a certain parameter combination a transient accelerating phase emerges as a pure interaction effect.

The dynamics of TeVeS in a homogeneous and isotropic universe is shown to be equivalent to the dynamics of an interacting two-component system, consisting of a scalar field and a "fluid", related to the matter part, with explicitly given coupling term. Scaling solutions (solutions with a constant ratio of the energy densities of both components) in the "Einstein frame" are found which are exponentially expanding or contracting with no remaining freedom for Bekenstein's $F$ function. In the "physical frame" these solutions are of the power-law type. An equivalent General Relativity (GR) picture of the dynamics suggests that it is the scalar field which plays the role of dark matter, while the "matter" has to mimic (phantom-type) dark energy.

The generalized Chaplygin gas (GCG) model is the prototype of a unified model of dark energy (DE) and dark matter (DM). It is characterized by equation-of-state (EoS) parameters $A$ and $\alpha$. We use a statistical analysis of the 2dFGRS data to constrain these parameters. In particular, we find that very small (close to zero) and very large values ($\alpha\gg 1$) of the equation-of-state parameter $\alpha$ are preferred. To test the validity of this type of unification of the dark sector we admit the existence of a separate DM component in addition to the Chaplygin gas and calculate the probability distribution for the fractional contributions of both components to the total energy density. This analysis favors a model for which the Universe is nearly entirely made up of the separate DM component with an almost negligible Chaplygin gas part. This confirms the results of a previous Newtonian analysis.

We model the dark sector of the cosmic substratum by a viscous fluid with an equation of state $p=-\zeta \Theta$, where $\Theta$ is the fluid-expansion scalar and $\zeta$ is the coefficient of bulk viscosity for which we assume a dependence $\zeta \propto \rho^{\nu}$ on the energy density $\rho$. The homogeneous and isotropic background dynamics coincides with that of a generalized Chaplygin gas with equation of state $p = - A/\rho^{\alpha}$. The perturbation dynamics of the viscous model, however, is intrinsically non-adiabatic and qualitatively different from the Chaplygin-gas case. In particular, it avoids short-scale instabilities and/or oscillations which apparently have ruled out unified models of the Chaplygin-gas type. We calculate the matter power spectrum and demonstrate that the non-adiabatic model is compatible with the data from the 2dFGRS and the SDSS surveys. A $\chi^{2}$-analysis shows, that for certain parameter combinations the viscous-dark-fluid (VDF) model is well competitive with the $\Lambda$CDM model. These results indicate that \textit{non-adiabatic} unified models can be seen as potential contenders for a General-Relativity-based description of the cosmic substratum.

Interactions between dark matter and dark energy with a given equation of state are known to modify the cosmic dynamics. On the other hand, the strength of these interactions is subject to strong observational constraints. Here we discuss a model in which the transition from decelerated to accelerated expansion of the Universe arises as a pure interaction phenomenon. Various cosmological scenarios that describe a present stage of accelerated expansion, like the LCDM model or a (generalized) Chaplygin gas, follow as special cases for different interaction rates. This unifying view on the homogeneous and isotropic background level is accompanied by a non-adiabatic perturbation dynamics which can be seen as a consequence of a fluctuating interaction rate.

We model the cosmic medium as the mixture of a generalized Chaplygin gas and a pressureless matter component. Within a neo-Newtonian approach we compute the matter power spectrum. The 2dFGRS data are used to discriminate between unified models of the dark sector and different models, for which there is separate dark matter, in addition to that accounted for by the generalized Chaplygin gas. Leaving the corresponding density parameters free, we find that the unified models are strongly disfavored. On the other hand, using unified model priors, the observational data are also well described, in particular for small and large values of the generalized Chaplygin gas parameter $\alpha$.

We model the cosmic medium as the mixture of a generalized Chaplygin gas and a pressureless matter component. Within a neo-Newtonian approach (in which, different from standard Newtonian cosmology, the pressure enters the homogeneous and isotropic background dynamics) we compute the matter power spectrum. The 2dFGRS data are used to discriminate between unified models of the dark sector (a purely baryonic matter component of roughly 5 percent of the total energy content and roughly 95 percent generalized Chaplygin gas) and different models, for which there is separate dark matter, in addition to that accounted for by the generalized Chaplygin gas. Leaving the corresponding density parameters free, we find that the unified models are strongly disfavored. On the other hand, using unified model priors, the observational data are also well described, in particular for small and large values of the generalized Chaplygin gas parameter $\alpha$. The latter result is in agreement with a recent, more qualitative but fully relativistic, perturbation analysis in Gorini et al.

A duality transformation that interrelates expanding and contracting cosmological models is shown to single out a duality invariant, interacting two-component description of any irrotational, geodesic and shearfree cosmic medium with vanishing three curvature scalar. We apply this feature to a system of matter and radiation, to a mixture of dark matter and dark energy, to minimal and conformal scalar fields, and to an enlarged Chaplygin gas model of the cosmic substratum. We extend the concept of duality transformations to cosmological perturbations and demonstrate the invariance of adiabatic pressure perturbations under these transformations.

We propose a scenario in which the dark components of the Universe are manifestations of a single bulk viscous fluid. Using dynamical system methods, a qualitative study of the homogeneous, isotropic background scenario is performed in order to determine the phase space of all possible solutions. The specific model which we investigate shares similarities with a generalized Chaplygin gas in the background but is characterized by non-adiabatic pressure perturbations. This model is tested against supernova type Ia and matter power spectrum data. Different from other unified descriptions of dark matter and dark energy, the matter power spectrum is well behaved, i.e., there are no instabilities or oscillations on small perturbation scales. The model is competitive in comparison with the currently most popular proposals for the description of the cosmological dark sector.

We demonstrate that a transition from decelerated to accelerated cosmic expansion arises as a pure interaction phenomenon if pressureless dark matter is coupled to holographic dark energy whose infrared cutoff scale is set by the Hubble length. In a spatially flat universe the ratio of the energy densities of both components remains constant through this transition, while it is subject to slow variations for non-zero spatial curvature. The coincidence problem is dynamized and reformulated in terms of the interaction rate. An early matter era is recovered since for negligible interaction at high redshifts the dark energy itself behaves as matter. A simple model for this dynamics is shown to fit the SN Ia data. The constant background energy density ratio simplifies the perturbation analysis which is characterized by non-adiabatic features.

Different models of the cosmic substratum which pretend to describe the present stage of accelerated expansion of the Universe like the $\Lambda$CDM model or a Chaplygin gas, can be seen as special realizations of a holographic dark energy cosmology if the option of an interaction between pressurless dark matter and dark energy is taken seriously. The corresponding interaction strength parameter plays the role of a cosmological constant. Differences occur at the perturbative level. In particular, the pressure perturbations are intrinsically non-adiabatic.

Perturbations in a Chaplygin gas, characterized by an equation of state $p = -A/\rho$, may acquire non-adiabatic contributions if spatial variations of the parameter $A$ are admitted. This feature is shown to be related to a specific internal structure of the Chaplygin gas. We investigate how perturbations of this type modify the adiabatic sound speed. A reduction of the effective sound speed compared with the adiabatic value is expected to suppress oscillations in the matter power spectrum. This text is an abridged version of the following reference: W. Zimdahl and J.C. Fabris, Classical and Quantum Gravity, {\bf 22}, 4311(2005).

We review the notion of holographic dark energy and assess its significance in the light of the well documented cosmic acceleration at the present time. We next propose a model of holographic dark energy in which the infrared cutoff is set by the Hubble scale. The model accounts for the aforesaid acceleration and, by construction, is free of the cosmic coincidence problem.

Interactions within the cosmic medium modify its equation of state. We discuss implications of interacting dark energy models both for the spatially homogenous background and for the perturbation dynamics.

In this Letter we demonstrate that any interaction of pressureless dark matter with holographic dark energy, whose infrared cutoff is set by the Hubble scale, implies a constant ratio of the energy densities of both components thus solving the coincidence problem. The equation of state parameter is obtained as a function of the interaction strength. For a variable degree of saturation of the holographic bound the energy density ratio becomes time dependent which is compatible with a transition from decelerated to accelerated expansion.

Perturbations in a Chaplygin gas, characterized by an equation of state $p = -A/\rho$, may acquire non-adiabatic contributions if spatial variations of the parameter $A$ are admitted. This feature is shown to be related to a specific internal structure of the Chaplygin gas. We investigate how perturbations of this type modify the adiabatic sound speed and influence the time dependence of the gravitational potential which gives rise to the Integrated Sachs-Wolfe effect in the anisotropy spectrum of the cosmic microwave background.

Motivated by the structure of one-loop vacuum polarization effects in curved spacetime we discuss a non-minimal extension of the Einstein-Maxwell equations. This formalism is applied to Bianchi I models with magnetic field. We obtain several exact solutions of the non-minimal system including those which describe an isotropization process. We show that there are inflationary solutions in which the cosmological constant is determined by the non-minimal coupling parameters. Furthermore, we find an isotropic de Sitter solution characterized by a "screening" of the magnetic field as a consequence of the non-minimal coupling.

Using relativistic kinetic theory, we study spherically symmetric, static equilibrium configurations of a collisionless Maxwell-Boltzmann gas with non-standard self-interactions, modelled by an effective one--particle force. The resulting set of equilibrium conditions represents a generalization of the classical Tolman-Oppenheimer-Volkov equations. We specify these conditions for two types of Lorentz--like forces: one coupled to the 4-acceleration and the 4--velocity and the other one coupled to the Riemann tensor. We investigate the weak field limits in each case and show that they lead to various Newtonian type configurations that are different from the usual isothermal sphere characterizing the conventional Newtonian Maxwell--Boltzmann gas. These configurations could provide a plausible phenomenological and theoretical description of galactic dark matter halo structures. We show how the self--interaction may act phenomenologically as an effective cosmological constant and discuss possible connections with Modified Newtonian Dynamics(MOND).

We investigate the role of a possible coupling of dark matter and dark energy. In particular, we explore the consequences of such an interaction for the coincidence problem, i.e., for the question, why the energy densities of dark matter and dark energy are of the same order just at the present epoch. We demonstrate, that, with the help of a suitable coupling, it is possible to reproduce any scaling solution $\rho_X \propto \rho_M a^\xi$, where $a$ is the scale factor of the Robertson-Walker metric and $\xi$ is a constant parameter. $\rho_X$ and $\rho_M$ are the densities of dark energy and dark matter, respectively. Furthermore, we show that an interaction between dark matter and dark energy can drive the transition from an early matter dominated era to a phase of accelerated expansion with a stable, stationary ratio of the energy densities of both components.

We find exact solutions leading to power law accelerated expansion for a homogeneous, isotropic and spatially flat universe, dominated by an interacting mixture of cold dark matter and a tachyonic field such that the ratio of the energy densities of both components at late times is constant and no coincidence problem arises.

Non-canonical scaling of the form $\rho_{x} \propto \rho_{m} a^{\xi}$, where $\rho_{x}$ and $\rho_{m}$ are the energy densities of dark energy and dark matter, respectively, provides a natural way to address the coincidence problem -why the two densities are comparable today. This non-canonical scaling is achieved by a suitable interaction between the dark components. We study the observational constraints imposed on such models from the cosmic microwave background (CMB) anisotropy spectrum. Using the recent WMAP results for the location of the CMB peaks, we determine the admissible parameter space and find that interacting models are well consistent with current observational bounds.

We argue that the recently introduced "statefinder parameters" (Sahni et al., JETP Lett. 77, 201 (2003)), that include the third derivative of the cosmic scale factor, are useful tools to characterize interacting quitessence models. We specify the statefinder parameters for two classes of models that solve, or at least alleviate, the coincidence problem.

We consider the homogeneous and isotropic cosmological fluid dynamics which is compatible with a homothetic, timelike motion, equivalent to an equation of state $\rho + 3P = 0$. By splitting the total pressure $P$ into the sum of an equilibrium part $p$ and a non-equilibrium part $\Pi$, we find that on thermodynamical grounds this split is necessarily given by $p = \rho$ and $\Pi = - (4/3)\rho$, corresponding to a dissipative stiff (Zel'dovich) fluid.

A curvature self-interaction of the cosmic gas is shown to mimic a cosmological constant or other forms of dark energy, such as a rolling tachyon condensate or a Chaplygin gas. Any given Hubble rate and deceleration parameter can be traced back to the action of an effective curvature force on the gas particles. This force self-consistently reacts back on the cosmological dynamics. The links between an imperfect fluid description, a kinetic description with effective antifriction forces, and curvature forces, which represent a non-minimal coupling of gravity to matter, are established.

We show that a suitable interaction between a scalar field and a matter fluid in a spatially homogeneous and isotropic spacetime can drive the transition from a matter dominated era to an accelerated expansion phase and simultaneously solve the coincidence problem of our present Universe. For this purpose we study the evolution of the energy density ratio of these two components. We demonstrate that a stationary attractor solution is compatible with an accelerated expansion of the Universe. We extend this study to account for dissipation effects due to interactions in the dark matter fluid. Finally, Type Ia supernovae and primordial nucleosynthesis data are used to constrain the parameters of the model.

We show that with the help of a suitable coupling between dark energy and cold dark matter it is possible to reproduce any scaling solution $\rho _{X}\propto \rho_{M}a^{\xi}$, where $\rho_{X}$ and $\rho_{M}$ are the densities of dark energy and dark matter, respectively. We demonstrate how the case $\xi = 1$ alleviates the coincidence problem. Future observations of supernovae at high redshift as well as quasar pairs which are planned to discriminate between different cosmological models will also provide direct constraints on the coupling between dark matter and dark energy.

Specific internal self-interactions in gaseous cosmic fluids are shown to give rise to effective negative pressures which may violate the strong energy condition. On this basis we discuss the transition from an initial de Sitter phase to a subsequent Friedmann-Lema\^{\i}tre-Robertson-Walker period as a non-equilibrium configuration of an ultrarelativistic gas. An accelerated expansion of the present universe is obtained as the consequence of a negative internal friction force which is self-consistently exerted on the non-relativistic cold dark matter (CDM) type microscopic constituents of the cosmic gas.

We demonstrate that a suitable coupling between a quintessence scalar field and a pressureless cold dark matter (CDM) fluid leads to a constant ratio of the energy densities of both components which is compatible with an accelerated expansion of the Universe.

Accelerated expansion of the Universe may result from an anti-frictional force that is self-consistently exerted on cold dark matter (CDM). Cosmic anti-friction is shown to give rise to an effective negative pressure of the cosmic medium. While other models introduce a component of dark energy besides ``standard'' CDM, we resort to a phenomenological one-component model of CDM with internal self-interactions. We demonstrate how the dynamics of the LambdaCDM model may be recovered as a special case of cosmic anti-friction.

Introducing an effective refraction index of an isotropic cosmic medium, we investigate the cosmological fluid dynamics which is consistent with a conformal, timelike symmetry of a corresponding ``optical'' metric. We demonstrate that this kind of symmetry is compatible with the existence of a negative viscous pressure and, consequently, with cosmological entropy production. We establish an exactly solvable model according to which the viscous pressure is a consequence of a self-interacting one-particle force which is self-consistently exerted on the microscopic particles of a relativistic gas. Furthermore, we show that a sufficiently high decay rate of the refraction index of an ultrarelativistic cosmic medium results in an inflationary expansion of the universe.

We explain an accelerated expansion of the present universe, suggested from observations of supernovae of type Ia at high redshift, by introducing an anti-frictional force that is self-consistently exerted on the particles of the cosmic substratum. Cosmic anti-friction, which is intimately related to ``particle production'', is shown to give rise to an effective negative pressure of the cosmic medium. While other explanations for an accelerated expansion (cosmological constant, quintessence) introduce a component of dark energy besides ``standard'' cold dark matter (CDM) we resort to a phenomenological one-component model of CDM with internal self-interactions. We demonstrate how the dynamics of the LambdaCDM model may be recovered as a special case of cosmic anti-friction. We discuss the connection with two-component models and obtain an attractor behavior for the ratio of the energy densities of both components which provides a possible phenomenological solution to the coincidence problem.

We describe the ``deflationary'' evolution from an initial de Sitter phase to a subsequent Friedmann-Lema\^{\i}tre-Robertson-Walker (FLRW) period as a specific non-equilibrium configuration of a self-interacting gas. The transition dynamics corresponds to a conformal, timelike symmetry of an ``optical'' metric, characterized by a refraction index of the cosmic medium which continously decreases from a very large initial value to unity in the FLRW phase.

We establish a general thermodynamic scheme for cosmic fluids with internal self-interactions and discuss equilibrium and non-equilibrium aspects of such systems in connection with (generalized) symmetry properties of the cosmological dynamics. As an example we construct an exactly solvable gas dynamical model of a ``deflationary'' transition from an initial de Sitter phase to a subsequent Friedmann-Lema\^{\i}tre-Robertson-Walker (FLRW) period. We demonstrate that this dynamics represents a manifestation of a conformal symmetry of an ``optical'' metric, characterized by a specific effective refraction index of the cosmic medium.

We reveal unifying thermodynamic aspects of so different phenomena as the cosmological electron-positron annihilation, the evaporation of primordial black holes with a narrow mass range, and the ``deflationary'' transition from an initial de Sitter phase to a subsequent standard Friedmann-Lema\^{\i}tre-Robertson-Walker begin (FLRW) behavior.

We discuss two different realizations of a ``deflationary'' scenario of the early universe which imply a smooth transition from an initial de Sitter stage to a subsequent radiation-dominated FLRW period. Thermodynamically, this transition is a non-equilibrium process which may also be interpreted as the decay of an initial cosmic vacuum into relativistic matter.

Combining the equivalence between cosmological particle creation and an effective viscous fluid pressure with the fact that the latter represents a dynamical degree of freedom within the second-order Israel-Stewart theory for imperfect fluids, we reconsider the possibility of accelerated expansion in fluid cosmology. We find an inherent self-limitation for the magnitude of an effective bulk pressure which is due to adiabatic (isentropic) particle production. For a production rate which depends quadratically on the Hubble rate we confirm the existence of solutions which describe a smooth transition from inflationary to noninflationary behavior and discuss their interpretation within the model of a decaying vacuum energy density. An alternative formulation of the effective imperfect fluid dynamics in terms of a minimally coupled scalar field is given. The corresponding potential is discussed and an entropy equivalent for the scalar field is found.

Combining the second-order entropy flow vector of the causal Israel-Stewart theory with the conformal Killing-vector property of $u_{i}/T$, where $u_{i}$ is the four-velocity of the medium and T its equilibrium temperature, we investigate generalized equilibrium states for cosmological fluids with nonconserved particle number. We calculate the corresponding equilibrium particle production rate and show that this quantity is reduced compared with the results of the previously studied first-order theory. Generalized equilibrium for massive particles turns out to be compatible with a dependence $\rho \propto a ^{-2}$ of the fluid energy density $\rho$ on the scale factor a of the Robertson-Walker metric and may be regarded as a realization of so-called K-matter.

We show that a relativistic gas may be at ``global'' equilibrium in the expanding universe for any equation of state $0 < p \leq \rho /3$, provided that the gas particles move under the influence of a self-interacting, effective one-particle force in between elastic binary collisions. In the force-free limit we recover the equilibrium conditions for ultrarelativistic matter which imply the existence of a conformal timelike Killing vector.

Primordial black holes with a narrow mass range are regarded as a nonrelativistic fluid component with an equation of state for dust. The impact of the black hole evaporation on the dynamics of the early universe is studied by resorting to a two-fluid model. We find periods of intense radiation reheating in the initial and final stages of the evaporation.

We show that a de Sitter spacetime is a solution of Einstein's field equations with the energy momentum tensor of a self-interacting, classical Maxwell-Boltzmann gas in collisional equilibrium. The self-interaction is described by a four-force which is quadratic in the (spatially projected) particle four-momenta. This force does not preserve the particle number and gives rise to an exponential increase in the comoving entropy of the universe while the temperature of the latter remains constant. These properties of a gas universe are related to the existence of a ``projector-conformal'' timelike Killing vector representing a symmetry which is ``in between'' the symmetries characterized by a Killing vector and those characterized by a conformal Killing vector.

The particles of a classical relativistic gas are supposed to move under the influence of a quasilinear (in the particle four-momenta), self-interacting force inbetween elastic, binary collisions. This force which is completely fixed by the equilibrium conditions of the gas, gives rise to an effective viscous pressure on the fluid phenomenological level. Earlier results concerning the possibility of accelerated expansion of the universe due to cosmological particle production are reinterpreted. A phenomenon such as power law inflation may be traced back to specific self-interacting forces keeping the particles of a gas universe in states of generalized equilibrium.

Using the existence of a covariant conserved quantity on large perturbation scales in a spatially flat perfect fluid or scalar field universe, we present a general formula for gauge-invariantly defined comoving energy density perturbations which encodes the entire linear perturbation dynamics in a closed time integral. On this basis we discuss perturbation modes in different cosmological epochs.

With the help of a conformal, timelike Killing-vector we define generalized equilibrium states for cosmological fluids with particle production. For massless particles the generalized equilibrium conditions require the production rate to vanish and the well known ``global'' equilibrium of standard relativistic thermodynamics is recovered as a limiting case. The equivalence between the creation rate for particles with nonzero mass and an effective viscous fluid pressure follows as a consequence of the generalized equilibrium properties. The implications of this equivalence for the cosmological dynamics are discussed, including the possibility of a power-law inflationary behaviour. For a simple gas a microscopic derivation for such kind of equilibrium is given on the basis of relativistic kinetic theory.

Previously defined covariant and gauge-invariant perturbation variables, representing, e.g., the fractional spatial energy density gradient on hypersurfaces of constant expansion, are used to simplify the linear perturbation analysis of a classical scalar field. With the help of conserved quantities on large scales we establish an exact first-order relation between comoving fluid energy density perturbations at `reentry' into the horizon and corresponding scalar field energy density perturbations at the first Hubble scale crossing during an early de Sitter phase of a standard inflationary scenario.

The linear cosmological perturbation theory of an almost homogeneous and isotropic perfect fluid universe is reconsidered and formally simplified by introducing new covariant and gauge-invariant variables with physical interpretations on hypersurfaces of constant expansion, constant curvature or constant energy density. The existence of conserved perturbation quantities on scales larger than the Hubble scale is discussed. The quantity which is conserved on large scales in a flat background universe may be expressed in terms of the fractional, spatial gradient of the energy density on constant expansion hypersurfaces or, alternatively, with the help of expansion or curvature perturbation variables on hypersurfaces of constant energy density. For nonvanishing background curvature the perturbation dynamics is most suitably described in terms of energy density perturbations on hypersurfaces of constant curvature.

It is shown that two reacting cosmological fluids, each of them perfect on its own, which exchange energy and momentum without preserving particle numbers, give rise to an entropy producing `reactive' bulk stress of the system as a whole, as soon as the detailed balance between decay and inverse decay processes is perturbed. This demonstrates explicitly that particle generation is dynamically equivalent to an effective bulk pressure. We derive a semiquantitative formula for the corresponding new kinetic coefficient and evaluate the latter for the out-of-equilibrium decay of heavy, nonrelativistic particles into radiation. It turns out that the associated reactive bulk viscosity may be more than one order of magnitude larger than the conventional bulk viscosity, calculated, e.g., in radiative hydrodynamics.

The reheating process in inflationary universe models is considered as an out-of-equilibrium mixture of two interacting and reacting fluids, and studied within the framework of causal, irreversible thermodynamics. The evolution of the temperature and the decay rate as determined by causal thermodynamics are estimated at different stages of the process. A simple model is also used to find the perturbations of the expansion rate, including the possibility of damped oscillations.

Particle production processes in the expanding universe are described within a simple kinetic model. The equilibrium conditions for a Maxwell-Boltzmann gas with variable particle number are investigated. We find that radiation and nonrelativistic matter may be in equilibrium at the same temperature provided the matter particles are created at a rate that is half the expansion rate. Using the fact that the creation of particles is dynamically equivalent to a nonvanishing bulk pressure we calculate the backreaction of this process on the cosmological dynamics. It turns out that the `adiabatic' creation of massive particles with an equilibrium distribution for the latter necessarily implies power-law inflation. Exponential inflation in this context is shown to become inconsistent with the second law of thermodynamics after a time interval of the order of the Hubble time.

Different cooling rates of interacting fluid components of the expanding cosmic medium give rise to entropy producing bulk stresses. With the help of effective bulk pressures, on the other hand, one may phenomenologically take into account cosmological particle production processes. We report on recent progress in bulk viscous cosmology both on the phenomenological level and on the level of relativistic kinetic theory.

A universe consisting of two interacting perfect fluids with the same 4-velocity is considered. A heuristic mean free time argument is used to show that the system as a whole cannot be perfect as well but neccessarily implies a nonvanishing bulk viscosity. A new formula for the latter is derived and compared with corresponding results of radiative hydrodynamics.

Recently, the phenomenological description of cosmological particle production processes in terms of effective viscous pressures has attracted some attention. Using a simple creation rate model we discuss the question to what extent this approach is compatible with the kinetic theory of a relativistic gas. We find the effective viscous pressure approach to be consistent with this model for homogeneous spacetimes but not for inhomogeneous ones.

The full causal M\"uller-Israel-Stewart (MIS) theory of dissipative processes in relativistic fluids is applied to a flat, homogeneous and isotropic universe with bulk viscosity. It is clarified in which sense the so called truncated version is a reasonable limiting case of the full theory. The possibility of bulk viscosity driven inflationary solutions of the full theory is discussed. As long as the particle number is conserved almost all these solutions exhibit an exponential increase of the temperature. Assuming that the bulk viscous pressure of the MIS theory may also be interpreted as an effective description for particle production processes, the thermodynamical behaviour of the Universe changes considerably. In the latter case the temperature increases at a lower rate or may remain constant during a hypothetical de Sitter stage, accompanied by a substantial growth of the comoving entropy.