We analyze the global behaviour of the growth index of cosmic inhomogeneities in an isotropic homogeneous universe filled by cold non-relativistic matter and dark energy (DE) with an arbitrary equation of state. Using a dynamical system approach, we find the critical points of the system. That unique trajectory for which the growth index $\gamma$ is finite from the asymptotic past to the asymptotic future is identified as the so-called heteroclinic orbit connecting the critical points $(\Omega_m=0,~\gamma_{\infty})$ in the future and $(\Omega_m=1,~\gamma_{-\infty})$ in the past. The first is an attractor while the second is a saddle point, confirming our earlier results. Further, in the case when a fraction of matter (or DE tracking matter) $\varepsilon \Omega^{\rm tot}_m$ remains unclustered, we find that the limit of the growth index in the past $\gamma_{-\infty}^{\varepsilon}$ does not depend on the equation of state of DE, in sharp contrast with the case $\varepsilon=0$ (for which $\gamma_{-\infty}$ is obtained). We show indeed that there is a mathematical discontinuity: one cannot obtain $\gamma_{-\infty}$ by taking $\lim_{\varepsilon \to 0} \gamma^{\varepsilon}_{-\infty}$ (i.e. the limits $\varepsilon\to 0$ and $\Omega^{\rm tot}_m\to 1$ do not commute). We recover in our analysis that the value $\gamma_{-\infty}^{\varepsilon}$ corresponds to tracking DE in the asymptotic past with constant $\gamma=\gamma_{-\infty}^{\varepsilon}$ found earlier.

In this paper we will study $R^2$-like inflation in a non-local modification of gravity which contains quadratic in Ricci scalar and Weyl tensor terms with analytic infinite derivative form-factors in the action. It is known that the inflationary solution of the local $R+R^2$ gravity remains a particular exact solution in this model. It was shown earlier that the power spectrum of scalar perturbations generated during inflation in the non-local setup remains the same as in the local $R+R^2$ inflation, whereas the power spectrum of tensor perturbations gets modified due to the non-local Weyl tensor squared term. In the present paper we go beyond 2-point correlators and compute the non-Gaussian parameter $f_{NL}$ related to 3-point correlations generated during inflation, which we found to be different from those in the original local inflationary model and scenarios alike based on a local gravity. We evaluate non-local corrections to the scalar bi-spectrum which give non-zero contributions to squeezed, equilateral and orthogonal configurations. We show that $f_{NL}\sim O(1)$ with an arbitrary sign is achievable in this model based on the choice of form-factors and the scale of non-locality. We present the predictions for the tensor-to-scalar ratio, $r$, and the tensor tilt, $n_t$. In contrast to standard inflation in a local gravity, here the possibility $n_t$>0 is not excluded. Thus, future CMB data can probe non-local behaviour of gravity at high space-time curvatures.

We study a metric cubic gravity theory considering odd-parity modes of linear inhomogeneous perturbations on a spatially homogeneous Bianchi type I manifold close to the isotropic de Sitter spacetime. We show that in the regime of small anisotropy, the theory possesses new degrees of freedom compared to General Relativity, whose kinetic energy vanishes in the limit of exact isotropy. From the mass dispersion relation we show that such theory always possesses at least one ghost mode as well as a very short-time-scale (compared to the Hubble time) classical tachyonic (or ghost-tachyonic) instability. In order to confirm our analytic analysis, we also solve the equations of motion numerically and we find that this instability is developed well before a single e-fold of the scale factor. This shows that this gravity theory, as it is, cannot be used to construct viable cosmological models.

Using redshift space distortion data, we perform model-independent reconstructions of the growth history of matter inhomogeneity in the expanding Universe using two methods: crossing statistics and Gaussian processes. We then reconstruct the corresponding history of the Universe background expansion and fit it to type Ia supernovae data, putting constraints on $(\Omega_\mathrm{m,0},\sigma_{8,0})$. The results obtained are consistent with the concordance flat-$\Lambda$CDM model and General Relativity as the gravity theory given the current quality of the inhomogeneity growth data.

With Planck CMB observations, we established the spectral amplitude and tilt of the primordial power spectrum. Evidence of a red spectral tilt ($n_\mathrm{s}=0.96$) at $8\sigma$ provides strong support for the inflationary mechanism, especially the slow-roll of the effective scalar field in its nearly flat potential as the generator of scalar primordial perturbations. With the next generation of Large Scale Structure surveys, we expect to probe primordial physics beyond the overall shape and amplitude of the main, smooth and slowly-changing part of the inflaton potential. Using the specifications for the upcoming Euclid survey, we investigate to what extent we can constrain the inflation potential beyond its established slow-roll behaviour. We provide robust forecasts with Euclid and Planck mock data from 9 fiducial power spectra that contain suppression and wiggles at different cosmological scales, using the Wiggly Whipped Inflation framework to generate these features in the primordial spectrum. We include both Euclid cosmic shear and galaxy clustering, with a conservative cut-off for non-linear scales. Using MCMC simulations, we obtain an improvement in constraints in the WWI potential, as well an improvement for the background cosmology parameters. We find that apart from improving the constraints on the overall scale of the inflationary potential by 40-50%, we can also identify oscillations in the primordial spectrum that are present within intermediate to small scales ($k\sim0.01-0.2\,\mathrm{Mpc^{-1}}$).

Emergence of $R^2$ inflation which is the best fit framework for CMB observations till date comes from the attempts to attack the problem of quantization of gravity which in turn have resulted in the trace anomaly discovery. Further developments in trace anomaly and different frameworks aiming to construct quantum gravity indicate an inevitability of non-locality in fundamental physics at small time and length scales. A natural question would be to employ the $R^2$ inflation as a probe for signatures of non-locality in the early Universe physics. Recent advances of embedding $R^2$ inflation in a string theory inspired non-local gravity modification provides very promising theoretical predictions connecting the non-local physics in the early Universe and the forthcoming CMB observations.

Primordial black holes (PBHs) generated by gravitational collapse of large primordial overdensities can be a fraction of the observed dark matter. In this paper, we introduce a mechanism to produce a large peak in the primordial power spectrum (PPS) in two-field inflationary models characterized by two stages of inflation based on a large non-canonical kinetic coupling. This mechanism is generic to several two-field inflationary models, due to a temporary tachyonic instability of the isocurvature perturbations at the transition between the two stages of inflation. We numerically compute the primordial perturbations from largest scales to the small scales corresponding to that of PBHs using an extension of BINGO (BI-spectra and Non-Gaussianity Operator). Moreover we numerically compute the stochastic background of gravitational waves (SBGW) produced by second order scalar perturbations within frequencies ranging from nano-Hz to KHz that covers the observational scales corresponding to Pulsar Timing Arrays, Square Kilometer Array to that of Einstein telescope. We discuss the prospect of its detection by these proposed and upcoming gravitational waves experiments.

We use $Planck$ 2018 data to constrain the simplest models of scalar-tensor theories characterized by a coupling to the Ricci scalar of the type $F(\sigma) R$ with $F(\sigma) = N_{pl}^2 + \xi \sigma^2$. We update our results with previous $Planck$ and BAO data releases obtaining the tightest constraints to date on the coupling parameters, that is $\xi < 5.5 \times 10^{-4}$ for $N_{pl}=0$ (induced gravity or equivalently extended Jordan-Brans-Dicke) and $(N_{pl} \sqrt{8 \pi G})-1 < 1.8 \times 10^{-5}$ for $\xi = -1/6$ (conformal coupling), both at 95% CL. Because of a modified expansion history after radiation-matter equality compared to the $\Lambda$CDM model, all these dynamical models accommodate a higher value for $H_0$ and therefore alleviate the tension between $Planck$/BAO and distance-ladder measurement from SNe Ia data from $4.4\sigma$ at best to $2.3\sigma$. We show that all these results are robust to changes in the neutrino physics. In comparison to the $\Lambda$CDM model, partial degeneracies between neutrino physics and the coupling to the Ricci scalar allow for smaller values $N_{\rm eff} \sim 2.8$, $1\sigma$ lower compared to the standard $N_{\rm eff} = 3.046$, and relax the upper limit on the neutrino mass up to 40%.

We generalize the notion of constant-roll inflation earlier introduced in General Relativity (GR) and $f(R)$ gravity to inflationary models in more general scalar-tensor gravity. A number of novel exact analytic solutions for a FLRW spatially flat cosmological background is found for this case. All forms of the scalar field potential and its coupling to gravity producing the exact de Sitter solution, while the scalar field is varying, are presented. In the particular cases of induced gravity and GR with a non-minimally coupled scalar field, all constant-roll inflationary solutions are found. In the former case they represent power-law inflation, while in the latter case the solution is novel and more complicated. Comparison of scalar perturbations generated during such inflation in induced gravity with observational data shows that the constant-roll parameter should be small, similar to constant-roll inflation in GR and $f(R)$ gravity. Then the solution reduces to the standard slow-roll one with small corrections.

It has been argued that the mixed Higgs-$R^2$ model acts as the UV extension of the Higgs inflation, pushing up its cut-off scale in the vacuum close up to the Planck scale. In this letter, we study the inflaton oscillation stage after inflation, focusing on the effective mass of the phase direction of the Higgs field, which can cause a violent preheating process. We find that the "spikes" in the effective mass of the phase direction observed in the Higgs inflation still appear in the mixed Higgs-$R^2$ model. While the spikes appear above the cut-off scale in the Higgs-only case, they appear below the cut-off scale when the model is extended with $R^2$ term though reheating cannot be completed in the violent particle production regime since the spikes get milder.

We study the cosmological constraints on the variation of the Newton's constant and on post-Newtonian parameters for simple models of scalar-tensor theory of gravity beyond the extended Jordan-Brans-Dicke theory. We restrict ourselves to an effectively massless scalar field with a potential $V \propto F^2$, where $F(\sigma)=N_{pl}^2+\xi\sigma^2$ is the coupling to the Ricci scalar considered. We derive the theoretical predictions for cosmic microwave background (CMB) anisotropies and matter power spectra by requiring that the effective gravitational strength at present is compatible with the one measured in a Cavendish-like experiment and by assuming adiabatic initial condition for scalar fluctuations. When comparing these models with $Planck$ 2015 and a compilation of baryonic acoustic oscilation (BAO) data, all these models accomodate a marginalized value for $H_0$ higher than in $\Lambda$CDM. We find no evidence for a statistically significant deviation from Einstein's general relativity. We find $\xi < 0.064$ ($|\xi| < 0.011$) at 95 % CL for $\xi > 0$ (for $\xi < 0$, $\xi \ne -1/6$). In terms of post-Newtonian parameters, we find $0.995 < \gamma_{\rm PN} < 1$ and $0.99987 < \beta_{\rm PN} < 1$ ($0.997 < \gamma_{\rm PN} < 1$ and $1 < \beta_{\rm PN} < 1.000011$) for $\xi >0$ (for $\xi < 0$). For the particular case of the conformal coupling, i.e. $\xi=-1/6$, we find constraints on the post-Newtonian parameters of similar precision to those within the Solar System.

We study the energy-momentum tensor and helicity of gauge fields coupled through $g \phi F \tilde{F}/4$ to a pseudo-scalar field $\phi$ driving inflation. Under the assumption of a constant time derivative of the background inflaton, we compute analitically divergent and finite terms of the energy density and helicity of gauge fields for any value of the coupling $g$. We introduce a suitable adiabatic expansion for mode functions of physical states of the gauge fields which correctly reproduces ultraviolet divergences in average quantities and identify corresponding counterterms. Our calculations shed light on the accuracy and the range of validity of approximated analytic estimates of the energy density and helicity terms previously existed in the literature in the strongly coupled regime only, i.e. for $g \dot \phi/(2H) \gg 1$. We discuss the implications of our analytic calculations for the backreaction of quantum fluctuations onto the inflaton evolution.

We investigate constraints on some key cosmological parameters by confronting metastable dark energy models with different combinations of the most recent cosmological observations. Along with the standard $\Lambda$CDM model, two phenomenological metastable dark energy models are considered: (\romannumeral1) DE decays exponentially, (\romannumeral2) DE decays into dark matter. We find that: (1) when considering the most recent supernovae and BAO data, and assuming a fiducial $\Lambda$CDM model, the inconsistency in the estimated value of the $\Omega_{\rm{m,0}}h^2$ parameter obtained by either including or excluding Planck CMB data becomes very much substantial and points to a clear tension~\citep{sahni2014model,zhao2017dynamical}; (2) although the two metastable dark energy models that we study provide greater flexibility in fitting the data, and they indeed fit the SNe Ia+BAO data substantially better than $\Lambda$CDM, they are not able to alleviate this tension significantly when CMB data are included; (3) while local measurements of the Hubble constant are significantly higher relative to the estimated value of $H_0$ in our models (obtained by fitting to SNe Ia and BAO data), the situation seems to be rather complicated with hints of inconsistency among different observational data sets (CMB, SNe Ia+BAO and local $H_0$ measurements). Our results indicate that we might not be able to remove the current tensions among different cosmological observations by considering simple modifications of the standard model or by introducing minimal dark energy models. A complicated form of expansion history, different systematics in different data and/or a non-conventional model of the early Universe might be responsible for these tensions.

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 perform here a global analysis of the growth index $\gamma$ behaviour from deep in the matter era till the far future. For a given cosmological model in GR or in modified gravity, the value of $\gamma(\Omega_{m})$ is unique when the decaying mode of scalar perturbations is negligible. However, $\gamma_{\infty}$, the value of $\gamma$ in the asymptotic future, is unique even in the presence of a nonnegligible decaying mode today. Moreover $\gamma$ becomes arbitrarily large deep in the matter era. Only in the limit of a vanishing decaying mode do we get a finite $\gamma$, from the past to the future in this case. We find further a condition for $\gamma(\Omega_{m})$ to be monotonically decreasing (or increasing). This condition can be violated inside general relativity (GR) for varying $w_{DE}$ though generically $\gamma(\Omega_{m})$ will be monotonically decreasing (like $\Lambda$CDM), except in the far future and past. A bump or a dip in $G_{\rm eff}$ can also lead to a significant and rapid change in the slope $\frac{d\gamma}{d\Omega_{m}}$. On a $\Lambda$CDM background, a $\gamma$ substantially lower (higher) than $0.55$ with a negative (positive) slope reflects the opposite evolution of $G_{\rm eff}$. In DGP models, $\gamma(\Omega_{m})$ is monotonically increasing except in the far future. While DGP gravity becomes weaker than GR in the future and $w^{DGP}\to -1$, we still get $\gamma_{\infty}^{DGP}= \gamma_{\infty}^{\Lambda CDM}=\frac{2}{3}$. In contrast, despite $G^{DGP}_{\rm eff}\to G$ in the past, $\gamma$ does not tend to its value in GR because $\frac{dG^{DGP}_{\rm eff}}{d\Omega_{m}}\Big|_{-\infty}\ne 0$.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

We combine model-independent reconstructions of the expansion history from the latest Pantheon supernovae distance modulus compilation and measurements from baryon acoustic oscillation to test some important aspects of the concordance model of cosmology namely the FLRW metric and flatness of spatial curvature. We then use the reconstructed expansion histories to fit growth measurement from redshift-space distortion and obtain strong constraints on $(\Omega_\mathrm{m},\gamma,\sigma_8)$ in a model independent manner. Our results show consistency with a spatially flat FLRW Universe with general relativity to govern the perturbation in the structure formation and the cosmological constant as dark energy. However, we can also see some hints of tension among different observations within the context of the concordance model related to high redshift observations ($z > 1$) of the expansion history. This supports earlier findings of Sahni et al. (2014) & Zhao et al. (2017) and highlights the importance of precise measurement of expansion history and growth of structure at high redshifts.

The form of the primordial power spectrum (PPS) of cosmological scalar (matter density) perturbations is not yet constrained satisfactorily in spite of the tremendous amount of information from the Cosmic Microwave Background (CMB) data. While a smooth power-law-like form of the PPS is consistent with the CMB data, some PPS with small non-smooth features at large scales can also fit the CMB temperature and polarization data with similar statistical evidence. Future CMB surveys cannot help distinguish all such models due to the cosmic variance at large angular scales. In this paper, we study how well we can differentiate be- tween such featured forms of the PPS not otherwise distinguishable using CMB data. We ran 15 N-body DESI-like simulations of these models to explore this approach. Showing that statistics such as the halo mass function and the two-point correlation function are not able to distinguish these models in a DESI-like survey, we advocate to avoid reducing the dimensionality of the problem by demonstrating that the use of a simple three-dimensional count-in-cell density field can be much more effective for the purpose of model distinction.

We study the quantum backreaction from inflationary fluctuations of a very light, non-minimally coupled spectator scalar and show that it is a viable candiate for dark energy. The problem is solved by suitably adapting the formalism of stochastic inflation. This allows us to self-consistently account for the backreaction on the background expansion rate of the Universe where its effects are large. This framework is equivalent to that of semiclassical gravity in which matter vacuum fluctuations are included at the one loop level, but purely quantum gravitational fluctuations are neglected. Our results show that dark energy in our model can be characterized by a distinct effective equation of state parameter (as a function of redshift) which allows for testing of the model at the level of the background.

We consider an inflationary scenario where the rate of inflaton roll defined by $\ddot\phi/H\dot \phi$ remains constant. The rate of roll is small for slow-roll inflation, while a generic rate of roll leads to the interesting case of `constant-roll' inflation. We find a general exact solution for the inflaton potential required for such inflaton behaviour. In this model, due to non-slow evolution of background, the would-be decaying mode of linear scalar (curvature) perturbations may not be neglected. It can even grow for some values of the model parameter, while the other mode always remains constant. However, this always occurs for unstable solutions which are not attractors for the given potential. The most interesting particular cases of constant-roll inflation remaining viable with the most recent observational data are quadratic hilltop inflation (with cutoff) and natural inflation (with an additional negative cosmological constant). In these cases even-order slow-roll parameters approach non-negligible constants while the odd ones are asymptotically vanishing in the quasi-de Sitter regime.

We consider the prospects of probing features in the primordial power spectrum with future Cosmic Microwave Background (CMB) polarization measurements. In the scope of the inflationary scenario, such features in the spectrum can be produced by local non-smooth pieces in an inflaton potential (smooth and quasi-flat in general) which in turn may originate from fast phase transitions during inflation in other quantum fields interacting with the inflaton. They can fit some outliers in the CMB temperature power spectrum which are unaddressed within the standard inflationary ${\mathrm{\Lambda}}$CDM model. We consider Wiggly Whipped Inflation (WWI) as a theoretical framework leading to improvements in the fit to the Planck 2015 temperature and polarization data in comparison with the standard inflationary models, although not at a statistically significant level. We show that some type of features in the potential within the WWI models, leading to oscillations in the primordial power spectrum that extend to intermediate and small scales can be constrained with high confidence (at 3$\sigma$ or higher confidence level) by an instrument as the Cosmic ORigins Explorer (CORE). In order to investigate the possible confusion between inflationary features and footprints from the reionization era, we consider an extended reionization history with monotonic increase of free electrons with decrease in redshift. We discuss the present constraints on this model of extended reionization and future predictions with CORE. We also project, to what extent, this extended reionization can create confusion in identifying inflationary features in the data.

We analyze a two-field inflationary model consisting of the Ricci scalar squared ($R^2$) term and the standard Higgs field non-minimally coupled to gravity in addition to the Einstein $R$ term. Detailed analysis of the power spectrum of this model with mass hierarchy is presented, and we find that one can describe this model as an effective single-field model in the slow-roll regime with a modified sound speed. The scalar spectral index predicted by this model coincides with those given by the $R^2$ inflation and the Higgs inflation implying that there is a close relation between this model and the $R^2$ inflation already in the original (Jordan) frame. For a typical value of the self-coupling of the standard Higgs field at the high energy scale of inflation, the role of the Higgs field in parameter space involved is to modify the scalaron mass, so that the original mass parameter in the $R^2$ inflation can deviate from its standard value when non-minimal coupling between the Ricci scalar and the Higgs field is large enough.

We propose a new class of metastable dark energy (DE) phenomenological models in which the DE decay rate does not depend on external parameters such as the scale factor or the curvature of the Universe. Instead, the DE decay rate is assumed to be a constant depending only on intrinsic properties of DE and the type of a decay channel, similar to case of the radioactive decay of unstable particles and nuclei. As a consequence, the DE energy density becomes a function of the proper time elapsed since its formation, presumably in the very early Universe. Such a natural type of DE decay can profoundly affect the expansion history of the Universe and its age. Metastable DE can decay in three distinct ways: (i) exponentially, (ii) into dark matter, (iii) into dark radiation. Testing metastable DE models with observational data we find that the decay half-life must be many times larger than the age of the Universe. Models in which dark energy decays into dark matter lead to lower values of the Hubble parameter at large redshifts relative to $\Lambda$CDM. Consequently these models provide a better fit to cosmological BAO data (especially data from from high redshift quasars) than concordance ($\Lambda$CDM) cosmology.

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.

It is natural to expect a consistent inflationary model of the very early Universe to be an effective theory of quantum gravity, at least at energies much less than the Planck one. For the moment, $R+R^2$, or shortly $R^2$, inflation is the most successful in accounting for the latest CMB data from the PLANCK satellite and other experiments. Moreover, recently it was shown to be ultra-violet (UV) complete via an embedding into an analytic infinite derivative (AID) non-local gravity. In this paper, we derive a most general theory of gravity that contributes to perturbed linear equations of motion around maximally symmetric space-times. We show that such a theory is quadratic in the Ricci scalar and the Weyl tensor with AID operators along with the Einstein-Hilbert term and possibly a cosmological constant. We explicitly demonstrate that introduction of the Ricci tensor squared term is redundant. Working in this quadratic AID gravity framework without a cosmological term we prove that for a specified class of space homogeneous space-times, a space of solutions to the equations of motion is identical to the space of backgrounds in a local $R^2$ model. We further compute the full second order perturbed action around any background belonging to that class. We proceed by extracting the key inflationary parameters of our model such as a spectral index ($n_s$), a tensor-to-scalar ratio ($r$) and a tensor tilt ($n_t$). It appears that $n_s$ remains the same as in the local $R^2$ inflation in the leading slow-roll approximation, while $r$ and $n_t$ get modified due to modification of the tensor power spectrum. This class of models allows for any value of $r<0.07$ with a modified consistency relation which can be fixed by future observations of primordial $B$-modes of the CMB polarization. This makes the UV complete $R^2$ gravity a natural target for future CMB probes.

In this paper we investigate the past evolution of an anisotropic Bianchi I universe in $R+R^2$ gravity. Using the dynamical system approach we show that there exists a new two-parameters set of solutions that includes both an isotropic "false radiation" solution and an anisotropic generalized Kasner solution, which is stable. We derive the analytic behaviour of the shear from a specific property of $f(R)$ gravity and the analytic asymptotic form of the Ricci scalar when approaching the initial singularity. Finally we numerically checked our results.

The previously introduced class of two-parametric phenomenological inflationary models in General Relativity in which the slow-roll assumption is replaced by the more general, constant-roll condition is generalized to the case of $f(R)$ gravity. A simple constant-roll condition is defined in the original Jordan frame, and exact expressions for a scalaron potential in the Einstein frame, for a function $f(R)$ (in the parametric form) and for inflationary dynamics are obtained. The region of the model parameters permitted by the latest observational constraints on the scalar spectral index and the tensor-to-scalar ratio of primordial metric perturbations generated during inflation is determined.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

We study quintessential inflation using a generalized exponential potential $V(\phi)\propto exp(-\lambda \phi^n/Mpl^n), n>1$, the model admits slow-roll inflation at early times and leads to close-to-scaling behaviour in the post inflationary era with an exit to dark energy at late times. We present detailed investigations of the inflationary stage in the light of the Planck 2015 results, study post-inflationary dynamics and analytically confirm the existence of an approximately scaling solution. Additionally, assuming that standard massive neutrinos are non-minimally coupled, makes the field $\phi$ dominant once again at late times giving rise to present accelerated expansion of the Universe. We derive observational constraints on the field and time-dependent neutrino masses. In particular, for $n=6 (8)$, the parameter $\lambda$ is constrained to be,$\log \lambda > -7.29 (-11.7)$; the model produces the spectral index of the power spectrum of primordial scalar (matter density) perturbations as $ n_s = 0.959 \pm 0.001 (0.961 \pm 0.001)$ and tiny tensor-to-scalar ratio, $r<1.72 \times 10^{-2} (2.32 \times 10^{-2})$ respectively. Consequently, the upper bound on possible values of the sum of neutrino masses $\Sigma m_{\nu} \lesssim 2.5$ eV significantly enhances compared to that in the standard $\Lambda$CDM model.

We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 \mu$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations.

The previously proposed class of phenomenological inflationary models in which the assumption of inflaton slow-roll is replaced by the more general, constant-roll condition is compared with the most recent cosmological observational data, mainly the Planck ones. Models in this two-parametric class which remain viable appear to be close to the slow-roll ones, and their inflaton potentials are close to (but still different from) that of the natural inflation model. Permitted regions for the two model parameters are presented.

The growth index $\gamma$ is an interesting tool to assess the phenomenology of dark energy (DE) models, in particular of those beyond general relativity (GR). We investigate the possibility for DE models to allow for a constant $\gamma$ during the entire matter and DE dominated stages. It is shown that if DE is described by quintessence (a scalar field minimally coupled to gravity), this behaviour of $\gamma$ is excluded either because it would require a transition to a phantom behaviour at some finite moment of time, or, in the case of tracking DE at the matter dominated stage, because the relative matter density $\Omega_m$ appears to be too small. An infinite number of solutions, with $\Omega_m$ and $\gamma$ both constant, are found with $w_{DE}=0$ corresponding to Einstein-de Sitter universes. For all modified gravity DE models satisfying $G_{\rm eff}\ge G$, among them the $f(R)$ DE models suggested in the literature, the condition to have a constant $w_{DE}$ is strongly violated at the present epoch. In contrast, DE tracking dust-like matter deep in the matter era, but with $\Omega_m <1$, requires $G_{\rm eff} > G$ and an example is given using scalar-tensor gravity for a range of admissible values of $\gamma$. For constant $w_{DE}$ inside GR, departure from a quasi-constant value is limited until today. Even a large variation of $w_{DE}$ may not result in a clear signature in the change of $\gamma$. The change however is substantial in the future and the asymptotic value of $\gamma$ is found while its slope with respect to $\Omega_m$ (and with respect to $z$) diverges and tends to $-\infty$.

Combining the stochastic and $\delta N$ formalisms, we derive non perturbative analytical expressions for all correlation functions of scalar perturbations in single-field, slow-roll inflation. The standard, classical formulas are recovered as saddle-point limits of the full results. This yields a classicality criterion that shows that stochastic effects are small only if the potential is sub-Planckian and not too flat. The saddle-point approximation also provides an expansion scheme for calculating stochastic corrections to observable quantities perturbatively in this regime. In the opposite regime, we show that a strong suppression in the power spectrum is generically obtained, and comment on the physical implications of this effect.

The $R+R^2$, shortly named "$R^2$" ("Starobinsky") inflationary model, represents a fully consistent example of a one-parameter inflationary scenario. This model has a "graceful exit" from inflation and provides a mechanism for subsequent creation and final thermalization of the standard matter. Moreover, it produces a very good fit of the observed spectrum of primordial perturbations. In the present paper we show explicitly that the $R^2$ inflationary spacetime is an exact solution of a range of weakly non-local (quasi-polynomial) gravitational theories, which provide an ultraviolet completion of the $R^2$ theory. These theories are ghost-free, super-renormalizable or finite at quantum level, and perturbatively unitary. Their spectrum consists of the graviton and the scalaron that is responsible for driving the inflation. Notably, any further extension of the spectrum leads to propagating ghost degrees of freedom. We are aimed at presenting a detailed construction of such theories in the so called Weyl basis. Further, we give a special account to the cosmological implications of this theory by considering perturbations during inflation. The highlight of the non-local model is the prediction of a modified, in comparison to a local $R^2$ model, value for the ratio of tensor and scalar power spectra $r$, depending on the parameters of the theory. The relevant parameters are under control to be successfully confronted with existing observational data. Furthermore, the modified $r$ can surely meet future observational constraints.

Quantum effects derived through conformal anomaly lead to an inflationary model that can be either stable or unstable. The unstable version requires a large dimensionless coefficient of about $5\times 10^8$ in front of the $R^2$ term that results in the inflationary regime in the $R+R^2$ ("Starobinsky") model being a generic intermediate attractor. In this case the non-local terms in the effective action are practically irrelevant, and there is a 'graceful exit' to a low curvature matter-like dominated stage driven by high-frequency oscillations of $R$ -- scalarons, which later decay to pairs of all particles and antiparticles, with the amount of primordial scalar (density) perturbations required by observations. The stable version is a genuine generic attractor, so there is no exit from it. We discuss a possible transition from stable to unstable phases of inflation. It is shown that this transition is automatic if the sharp cut-off approximation is assumed for quantum corrections in the period of transition. Furthermore, we describe two different quantum mechanisms that may provide a required large $\,R^2$-term in the transition period.

With the Planck 2015 Cosmic Microwave Background (CMB) temperature and polarization data, we search for possible features in the primordial power spectrum (PPS). We revisit the Wiggly Whipped Inflation (WWI) framework and demonstrate how generation of some particular primordial features can improve the fit to Planck data. WWI potential allows the scalar field to transit from a steeper potential to a nearly flat potential through a discontinuity either in potential or in its derivatives. WWI offers the inflaton potential parametrizations that generate a wide variety of features in the primordial power spectra incorporating most of the localized and non-local inflationary features that are obtained upon reconstruction from temperature and polarization angular power spectrum. At the same time, in a single framework it allows us to have a background parameter estimation with a nearly free-form primordial spectrum. Using Planck 2015 data, we constrain the primordial features in the context of Wiggly Whipped Inflation and present the features that are supported both by temperature and polarization. WWI model provides more than $13$ improvement in $\chi^2$ fit to the data with respect to the best fit power law model considering combined temperature and polarization data from Planck and B-mode polarization data from BICEP and Planck dust map. We use 2-4 extra parameters in the WWI model compared to the featureless strict slow roll inflaton potential. We find that the differences between the temperature and polarization data in constraining background cosmological parameters such as baryon density, cold dark matter density are reduced to a good extent if we use primordial power spectra from WWI. We also discuss the extent of bispectra obtained from the best potentials in arbitrary triangular configurations using the BI-spectra and Non-Gaussianity Operator (BINGO).

We present a systematic analysis of homogeneous and isotropic cosmologies in a particular Horndeski model with Galileon shift symmetry, containing also a $\Lambda$-term and a matter. The model, sometimes called Fab Five, admits a rich spectrum of solutions. Some of them describe the standard late time cosmological dynamic dominated by the $\Lambda$-term and matter, while at the early times the universe expands with a constant Hubble rate determined by the value of the scalar kinetic coupling. For other solutions the $\Lambda$-term and matter are screened at all times but there are nevertheless the early and late accelerating phases. The model also admits bounces, as well as peculiar solutions describing "the emergence of time". Most of these solutions contain ghosts in the scalar and tensor sectors. However, a careful analysis reveals three different branches of ghost-free solutions, all showing a late time acceleration phase. We analyze the dynamical stability of these solutions and find that all of them are stable in the future, since all their perturbations stay bounded at late times. However, they all turn out to be unstable in the past, as their perturbations grow violently when one approaches the initial spacetime singularity. We therefore conclude that the model has no viable solutions describing the whole of the cosmological history, although it may describe the current acceleration phase. We also check that the flat space solution is ghost-free in the model, but it may acquire ghost in more general versions of the Horndeski theory.

We consider the possibility to produce a bouncing universe in the framework of scalar-tensor gravity models in which the scalar field potential may be negative, and even unbounded from below. We find a set of viable solutions with nonzero measure in the space of initial conditions passing a bounce, even in the presence of a radiation component, and approaching a constant gravitational coupling afterwards. Hence we have a model with a minimal modification of gravity in order to produce a bounce in the early universe with gravity tending dynamically to general relativity (GR) after the bounce.

We address the cosmological role of an additional ${\cal O}(1)$ eV sterile neutrino in modified gravity models. We confront the present cosmological data with predictions of the FLRW cosmological model based on a variant of $f(R)$ modified gravity proposed by one of the authors previously. This viable cosmological model which deviation from general relativity with a cosmological constant $\Lambda$ decreases as $R^{-2n}$ for large, but not too large values of the Ricci scalar $R$ provides an alternative explanation of present dark energy and the accelerated expansion of the Universe. Various up-to-date cosmological data sets exploited include Planck CMB anisotropy, CMB lensing potential, BAO, cluster mass function and Hubble constant measurements. We find that the CMB+BAO constraints strongly the sum of neutrino masses from above. This excludes values $\lambda\sim 1$ for which distinctive cosmological features of the model are mostly pronounced as compared to the $\Lambda$CDM model, since then free streaming damping of perturbations due to neutrino rest masses is not sufficient to compensate their extra growth occurring in $f(R)$ gravity. Thus, we obtain $\lambda>8.2$ ($2\sigma$) with cluster systematics and $\lambda>9.4$ ($2\sigma$) without that. In the latter case we find for the sterile neutrino mass $0.47\,\,\rm{eV}$$\,<\,$$m_{\nu,\,\rm{sterile}}$$\,<\,$$1\,\,\rm{eV}$ ($2\sigma$) assuming the active neutrinos are massless, not significantly larger than in the standard $\Lambda$CDM with the same data set: $0.45\,\,\rm{eV}$$\,<\,$$m_{\nu,\,\rm{sterile}}$$\,<\,$$0.92\,\,\rm{eV}$ ($2\sigma$). However, a possible discovery of a sterile neutrino with the mass $m_{\nu,\,\rm{sterile}} \approx 1.5\,$eV motivated by various anomalies in neutrino oscillation experiments would favor cosmology based on $f(R)$ gravity rather than the $\Lambda$CDM model.

Baryon Acoustic Oscillations (BAO) allow us to determine the expansion history of the Universe, thereby shedding light on the nature of dark energy. Recent observations of BAO's in the SDSS DR9 and DR11 have provided us with statistically independent measurements of $H(z)$ at redshifts of 0.57 and 2.34, respectively. We show that these measurements can be used to test the cosmological constant hypothesis in a model independent manner by means of an improved version of the $Om$ diagnostic. Our results indicate that the SDSS DR11 measurement of $H(z) = 222 \pm 7$ km/sec/Mpc at $z = 2.34$, when taken in tandem with measurements of $H(z)$ at lower redshifts, imply considerable tension with the standard $\Lambda$CDM model. Our estimation of the new diagnostic $Omh^2$ from SDSS DR9 and DR11 data, namely $Omh^2 \approx 0.122 \pm 0.01$, which is equivalent to $\Omega_{0m}h^2$ for the spatially flat $\Lambda$CDM model, is in tension with the value $\Omega_{0m}h^2 = 0.1426 \pm 0.0025$ determined for $\Lambda$CDM from Planck+WP. This tension is alleviated in models in which the cosmological constant was dynamically screened (compensated) in the past. Such evolving dark energy models display a pole in the effective equation of state of dark energy at high redshifts, which emerges as a smoking gun test for these theories.

Motivated by the idea that inflation occurs at the GUT symmetry breaking scale, in this paper we construct a new class of large field inflaton potentials where the inflaton starts with a power law potential; after initial period of relative fast roll that lasts until after a few e-folds inside the horizon, it transits to the attractor of the slow roll part of the potential with a lower power. Due to the initial fast roll stages of inflation, we find a suppression in scalar primordial power at large scales and at the same time the choice of the potential can provide us a tensor primordial spectrum with high amplitude. This suppression in scalar power with a large tensor-to-scalar ratio helps us to reconcile the Planck and BICEP2 data in a single framework. We find that a transition from a cubic to quadratic form of inflaton potential generates an appropriate suppression in power of scalar primordial spectrum that provides significant improvement in fit compared to power law model when compared with Planck and BICEP2 data together. We calculate the extent of non-Gaussianity, specifically, the bispectrum for the best fit potential and show that it is consistent with Planck bispectrum constraints.

Motivated by BICEP2 results on the CMB polarization B-mode which imply primordial gravitational waves are produced when the Universe has the expansion rate of about $H \approx 10^{14}$ GeV, and by deviations from a smooth power-law behaviour for multipoles $\ell <50$ in the CMB temperature anisotropy power spectrum found in the WMAP and Planck experiments, we have expanded our class of large field inflationary models that fit both the BICEP2 and Planck CMB observations consistently. These best-fitted large field models are found to have a transition from a faster roll to the slow roll $V(\phi)=m^2 \phi^2/2$ inflation at a field value around 14.6~${\rm M_{Pl}}$ and thus a potential energy of $V(\phi) \sim (10^{16}\,{\rm GeV})^4$. In general this transition with sharp features in the inflaton potential produces not only suppression of scalars relative to tensor modes at small $k$ but also introduces wiggles in the primordial perturbation spectrum. These wiggles are shown to be useful to explain some localized features in the CMB angular power spectrum and can also have other observational consequences. Thus, primordial GW can be used now to make a tomography of inflation determining its fine structure. The resulting Wiggly Whipped Inflation scenario is described in details and the anticipated perturbation power spectra, CMB power spectra, non-Gaussianity and other observational consequences are calculated and compared to existing and forthcoming observations.

Combining Planck CMB temperature [1] and BICEP2 B-mode polarization data [2,3] we show qualitatively that, assuming inflationary consistency relation, the power-law form of the scalar primordial spectrum is ruled out at more than $3\sigma$ CL. This is an important finding, since the power-law form of the scalar primordial spectrum is one of the main assumptions of concordance model of cosmology and also a direct prediction of many inflationary scenarios. We show that a break or step in the form of the primordial scalar perturbation spectrum, similar to what we studied recently analyzing Planck data [4] can address both Planck and BICEP2 results simultaneously. Our findings also indicate that the data may require more flexibilities than what running of scalar spectral index can provide. Finally we show that an inflaton potential, originally appeared in [5] can generate both the step and the break model of scalar primordial spectrum in two different limits. The discussed potential is found to be favored by Planck data but marginally disfavored by BICEP2 results as it produces slightly lower amplitude of tensor primordial spectrum. Hence, if the tensor-to-scalar ratio ($r$) quoted by BICEP2 persists, it is of importance that we generate inflationary models with large $r$ and at the same time provide suppression in scalar primordial spectrum at large scales.

We study effects on the luminosity distance of a local inhomogeneity seeded by primordial curvature perturbations of the type predicted by the inflationary scenario and constrained by the cosmic microwave background radiation. We find that a local underdensity originated from a one, two or three standard deviations peaks of the primordial curvature perturbations field can induce corrections to the value of a cosmological constant of the order of $0.6%,1%,1.5%$ respectively. These effects cannot be neglected in the precision cosmology era in which we are entering. Our results can be considered an upper bound for the effect of the monopole component of the local non linear structure which can arise from primordial curvature perturbations and requires a fully non perturbative relativistic treatment.

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in May 2013 as a large-class mission for investigating within the framework of the ESA Cosmic Vision program a set of important scientific questions that require high resolution, high sensitivity, full-sky observations of the sky emission at wavelengths ranging from millimeter-wave to the far-infrared. PRISM's main objective is to explore the distant universe, probing cosmic history from very early times until now as well as the structures, distribution of matter, and velocity flows throughout our Hubble volume. PRISM will survey the full sky in a large number of frequency bands in both intensity and polarization and will measure the absolute spectrum of sky emission more than three orders of magnitude better than COBE FIRAS. The aim of this Extended White Paper is to provide a more detailed overview of the highlights of the new science that will be made possible by PRISM

We compare two different probes of the expansion history of the universe, namely, luminosity distances from type Ia supernovae and angular diameter distances from galaxy clusters, using the Bayesian interpretation of Crossing statistic [1, 2] in conjunction with the assumption of cosmic duality relation. Our analysis is conducted independently of any a-priori assumptions about the nature of dark energy. The model independent method which we invoke searches for inconsistencies between SNIa and galaxy cluster data sets. If detected such an inconsistency would imply the presence of systematics in either of the two data sets. Simulating observations based on expected WFIRST supernovae data and X-ray eROSITA + SZ Planck cluster data, we show that our method allows one to detect systematics with high precision and without advancing any hypothesis about the nature of dark energy.

It is shown that the tension between recent neutrino oscillation experiments, favoring sterile neutrinos with masses of the order of 1 eV, and cosmological data which impose stringent constraints on neutrino masses from the free streaming suppression of density fluctuations, can be resolved in models of the present accelerated expansion of the Universe based on f(R) gravity.

Baryon Acoustic Oscillations (BAO) provide an important standard ruler which can be used to probe the recent expansion history of our universe. We show how a simple extension of the Om diagnostic, which we call Om3, can combine standard ruler information from BAO with standard candle information from type Ia supernovae (SNIa) to yield a powerful novel null diagnostic of the cosmological constant hypothesis. A unique feature of Om3 is that it requires minimal cosmological assumptions since its determination does not rely upon prior knowledge of either the current value of the matter density and the Hubble constant, or the distance to the last scattering surface. Observational uncertainties in these quantities therefore do not affect the reconstruction of Om3. We reconstruct Om3 using the Union 2.1 SNIa data set and BAO data from SDSS, WiggleZ and 6dFGS. Our results are consistent with dark energy being the cosmological constant. We show how Om and Om3 can be used to obtain accurate model independent constraints on the properties of dark energy from future data sets such as BigBOSS.

Using the exact Lemaitre-Bondi-Tolman solution with a non-vanishing cosmological constant $\Lambda$, we investigate how the presence of a local spherically-symmetric inhomogeneity can affect apparent cosmological observables, such as the deceleration parameter or the effective equation of state of dark energy (DE), derived from the luminosity distance under the assumption that the real space-time is exactly homogeneous and isotropic. The presence of a local underdensity is found to produce apparent phantom behavior of DE, while a locally overdense region leads to apparent quintessence behavior. We consider relatively small large scale inhomogeneities which today are not linear and could be seeded by primordial curvature perturbations compatible with CMB bounds. Our study shows how observations in an inhomogeneous $\Lambda$CDM universe with initial conditions compatible with the inflationary beginning, if interpreted under the wrong assumption of homogeneity, can lead to the wrong conclusion about the presence of "fake" evolving dark energy instead of $\Lambda$.

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.

We search for hidden mirror symmetries at large angular scales in the WMAP 7 year Internal Linear Combination map of CMB temperature anisotropies using global pixel based estimators introduced for this aim. Two different axes are found for which the CMB intensity pattern is anomalously symmetric (or anti-symmetric) under reflection with respect to orthogonal planes at the 99.84(99.96)% CL (confidence level), if compared to a result for an arbitrary axis in simulations without the symmetry. We have verified that our results are robust to the introduction of the galactic mask. The direction of such axes is close to the CMB kinematic dipole and nearly orthogonal to the ecliptic plane, respectively. If instead the real data are compared to those in simulations taken with respect to planes for which the maximal mirror symmetry is generated by chance, the confidence level decreases to 92.39 (76.65)%. But when the effect in question translates into the anomalous alignment between normals to planes of maximal mirror (anti)-symmetry and these natural axes mentioned. We also introduce the representation of the above estimators in the harmonic domain, confirming the results obtained in the pixel one. The symmetry anomaly is shown to be almost entirely due to low multipoles, so it may have a cosmological and even primordial origin. Contrary, the anti-symmetry one is mainly due to intermediate multipoles that probably suggests its non-fundamental nature. We have demonstrated that these anomalies are not connected to the known issue of the low variance in WMAP observations and we have checked that axially symmetric parts of these anomalies are small, so that the axes are not the symmetry ones.

Our goal is to see how density waves of different scale combine to form voids between galaxy systems of various scale. We perform numerical simulations of structure formation in cubes of size 100 and 256 Mpc/h, with resolutions 256^3 and 512^3 particles and cells. To understand the role of density perturbations of various scale we cut power spectra at scales from 8 to 128 Mpc/h, using in all series identical initial random realisations. We find that small haloes and short filaments form all over the simulation box, if perturbations only up to scale 8 Mpc/h are present. The phenomenon of large multi-scale voids in the cosmic web requires the presence of an extended spectrum of primordial density perturbations. The void phenomenon is due to the action of two processes: the synchronisation of density perturbations of medium and large scales, and the suppression of galaxy formation in low-density regions by the combined action of negative sections of medium- and large-scale density perturbations, so that their densities are less than the mean density, and thus during the evolution their densities decrease.

We follow the evolution of galaxy systems in numerical simulation. Our goal is to understand the role of density perturbations of various scales in the formation and evolution of the cosmic web. We perform numerical simulations with the full power spectrum of perturbations, and with spectrum cut at long wavelengths. Additionally, we have one model, where we cut the intermediate waves. We analyze the density field and study the void sizes and density field clusters in different models. Our analysis shows that the fine structure (groups and clusters of galaxies) is created by small-scale density perturbations of scale $\leq 8$ \Mpc. Filaments of galaxies and clusters are created by perturbations of intermediate scale from $\sim 8$ to $\sim 32$ \Mpc, superclusters of galaxies by larger perturbations. We conclude that the scale of the pattern of the cosmic web is determined by density perturbations of scale up to $\sim 100$ \Mpc. Larger perturbations do not change the pattern of the web, but modulate the richness of galaxy systems, and make voids emptier. The stop of the increase of the scale of the pattern of the cosmic web with increasing scale of density perturbations can probably be explained as the freezing of the web at redshift $z\simeq 0.7$.

We address the issue why a cosmological constant (dark energy) possesses a small positive value instead of being zero. Motivated by the cosmic landscape picture, we mimic the dark energy by a scalar field with potential wells and show that other degrees of freedom interacting with it can localize this field by decoherence in one of the wells. Dark energy can then acquire a small positive value. We also show that the additional degrees of freedom enhance the tunneling rate between the wells. The consideration is performed in detail for the case of two wells and then extended to a large number of wells.

It is known that scalar-tensor theory of gravity admits regular crossing of the phantom divide line $w_{\DE}=-1$ for dark energy, and existing viable models of present dark energy for its particular case -- $f(R)$ gravity -- possess one such crossing in the recent past, after the end of the matter dominated stage. It was recently noted that during the future evolution of these models the dark energy equation of state $w_{\DE}$ may oscillate with an arbitrary number of phantom divide crossings. In this paper we prove that the number of crossings can be infinite, present an analytical condition for the existence of this effect and investigate it numerically. With the increase of the present mass of the scalaron (a scalar particle appearing in $f(R)$ gravity) beyond the boundary of the appearance of such oscillations, their amplitude is shown to decrease very fast. As a result, the effect quickly becomes small and its beginning is shifted to the remote future.

We study $F(R)$ modified gravity models which are capable of driving the accelerating epoch of the Universe at the present time whilst not destroying the standard Big Bang and inflationary cosmology. Recent studies have shown that a weak curvature singularity with $|R|\to\infty$ can arise generically in viable $F(R)$ models of present dark energy (DE) signaling an internal incompleteness of these models. In this work we study how this problem is cured by adding a quadratic correction with a sufficiently small coefficient to the $F(R)$ function at large curvatures. At the same time, this correction eliminates two more serious problems of previously constructed viable $F(R)$ DE models: unboundedness of the mass of a scalar particle (scalaron) arising in $F(R)$ gravity and the scalaron overabundance problem. Such carefully constructed models can also yield both an early time inflationary epoch and a late time de Sitter phase with vastly different values of $R$. The reheating epoch in these combined models of primordial and present dark energy is completely different from that of the old $R + R^{2}/6M^{2}$ inflationary model, mainly due to the fact that values of the effective gravitational constant at low and intermediate curvatures are different for positive and negative $R$. This changes the number of e-folds during the observable part of inflation that results in a different value of the primordial power spectrum index.

We find the natural embedding of the (R+R^2)-inflationary model into the recently constructed N=1 F(\cal R)-supergravity. It gives a simple and viable realization of chaotic inflation in supergravity. The only requirement for a slow-roll inflation is the existence of the (\cal R)^3-term with an anomalously large coefficient in Taylor expansion of the F(\cal R) function, where \cal R is the covariantly-chiral scalar supercurvature superfield.

The standard field-theoretical approach to the slow-roll inflation is introduced. We then show as, in order to calculate the mean square of the canonical gauge invariant quantum fluctuations associated to a generic field, the logarithm of the scale factor has to be used as the time variable in the Fokker-Planck equation in the stochastic approach. Then we compute the growth of different test fields with a small effective mass during slow-roll inflationary models, comparing the results with the one for the gauge invariant canonical fluctuation associated to the inflaton, the Mukhanov variable. We find that in most of the single fields inflationary models such fluctuation grows faster than any test field with a non-negative effective mass, with the exception of hybrid models.

We have investigated the evolution of a homogeneous isotropic background of the Universe and inhomogeneous subhorizon matter density perturbations in viable $f(R)$ models of present dark energy and cosmic acceleration analytically and numerically. It is found that viable $f(R)$ models generically exhibit recent crossing of the phantom boundary $w_{\rm DE}=-1$. Furthermore, it is shown that the growth index of perturbations depends both on time and wavenumber. This anomalous growth may explain properties of the observational matter power spectrum from the SDSS data and can also partially counteract the spectrum suppression by massive neutrinos making larger values of the total sum of neutrino rest masses possible.

The effect of massive neutrinos on matter power spectrum is discussed in the context of $f(R)$ gravity. It is shown that the anomalous growth of density fluctuations on small scales due to the scalaron force can be compensated by free streaming of neutrinos. As a result, models which predict observable deviation of the equation-of-state parameter $w_\DE$ from $w_\DE=-1$ can be reconciled with observations of matter clustering if the total neutrino mass is $O(0.5 \eV)$.

We compute the growth of the mean square of quantum fluctuations of test fields with small effective mass during a slowly changing, nearly de Sitter stage which took place in different inflationary models. We consider a minimally coupled scalar with a small mass, a modulus with an effective mass $ \propto H^2$ (with $H$ as the Hubble parameter) and a massless non-minimally coupled scalar in the test field approximation and compare the growth of their relative mean square with the one of gauge-invariant inflaton fluctuations. We find that in most of the single field inflationary models the mean square gauge invariant inflaton fluctuation grows {\em faster} than any test field with a non-negative effective mass. Hybrid inflationary models can be an exception: the mean square of a test field can dominate over the gauge invariant inflaton fluctuation one on suitably choosing parameters. We also compute the stochastic growth of quantum fluctuation of a second field, relaxing the assumption of its zero homogeneous value, in a generic inflationary model; as a main result, we obtain that the equation of motion of a gauge invariant variable associated, order by order, with a generic quantum scalar fluctuation during inflation can be obtained only if we use the number of e-folds as the time variable in the corresponding Langevin and Fokker-Planck equations for the stochastic approach. We employ this approach to derive some bounds in the case of a model with two massive fields.

Some recent papers have claimed the existence of static, spherically symmetric wormhole solutions to gravitational field equations in the absence of ghost (or phantom) degrees of freedom. We show that in some such cases the solutions in question are actually not of wormhole nature while in cases where a wormhole is obtained, the effective gravitational constant G_eff is negative in some region of space, i.e., the graviton becomes a ghost. In particular, it is confirmed that there are no vacuum wormhole solutions of the Brans-Dicke theory with zero potential and the coupling constant \omega > -3/2, except for the case \omega = 0; in the latter case, G_eff < 0 in the region beyond the throat. The same is true for wormhole solutions of F(R) gravity: special wormhole solutions are only possible if F(R) contains an extremum at which G_eff changes its sign.

Evolution of a background space-time metric and sub-horizon matter density perturbations in the Universe is numerically analyzed in viable $f(R)$ models of present dark energy and cosmic acceleration. It is found that viable models generically exhibit recent crossing of the phantom boundary $w_{\rm DE}=-1$. Furthermore, it is shown that, as a consequence of the anomalous growth of density perturbations during the end of the matter-dominated stage, their growth index evolves non-monotonically with time and may even become negative temporarily.

For a class of viable cosmological models in $f(R)$ gravity which deviation from the Einstein gravity decreases as a inverse power law of the Ricci scalar $R$ for large $R$, an analytic solution for density perturbations in the matter component during the matter dominated stage is obtained in terms of hypergeometric functions. An analytical expression for the matter transfer function at scales much less than the present Hubble scale is also obtained.

We present numerical calculation of the evolution of a background space-time metric and sub-horizon matter density perturbations in viable $f(R)$ gravity models of present dark energy and cosmic acceleration. We found that viable models generically exhibit recent crossing of the phantom boundary $w_{DE}=-1$. Moreover, as a consequence of the anomalous growth of density perturbations during the end of the matter-dominated stage, their growth index evolves non-monotonically with time and may even become negative.

There exist consistent low energy effective field theories describing gravity in the Higgs phase that allow the coexistence of massive gravitons and the conventional 1/r potential of gravity. In an effort to constrain the value of the graviton mass in these theories, we study the tensor contribution to the CMB temperature anisotropy and polarization spectra in the presence of a non-vanishing graviton mass. We find that the observation of a B-mode signal consistent with the spectrum predicted by inflationary models would provide the strongest limit yet on the mass of an elementary particle -- a graviton -- at a level of m\lesssim 10^(-30) eV\approx(10 Mpc)^(-1). We also find that a graviton mass in the range between (10 Mpc)^(-1) and (10 kpc)^(-1) leads to interesting modifications of the polarization spectrum. The characteristic signature of a graviton mass in this range would be a plateau in the B-mode spectrum up to angular multipoles of l\sim 100. For even larger values of the graviton mass the tensor contribution to the CMB spectra becomes strongly suppressed.

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 investigate the course of cosmic expansion in its `recent past' using the Constitution SN Ia sample (which includes CfA data at low redshifts), jointly with signatures of baryon acoustic oscillations (BAO) in the galaxy distribution and fluctuations in the cosmic microwave background (CMB). Earlier SN Ia data sets could not address this issue because of a paucity of data at low redshifts. Allowing the equation of state of dark energy (DE) to vary, we find that a coasting model of the universe (q_0=0) fits the data about as well as LCDM. This effect, which is most clearly seen using the recently introduced `Om' diagnostic, corresponds to an increase of Om(z) and q(z) at redshifts z \lleq 0.3. In geometrical terms, this suggests that cosmic acceleration may have already peaked and that we are currently witnessing its slowing down. The case for evolving DE strengthens if a subsample of the Constitution set consisting of SNLS+ESSENCE+CfA SN Ia data is analysed in combination with BAO+CMB using the same statistical methods. The effect we observe could correspond to DE decaying into dark matter (or something else). A toy model which mimics this process agrees well with the combined SN Ia+BAO+CMB data.

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.

For a large class of dark energy (DE) models, for which the effective gravitational constant is a constant and there is no direct exchange of energy between DE and dark matter (DM), knowledge of the expansion history suffices to reconstruct the growth factor of linearized density perturbations in the non-relativistic matter component on scales much smaller than the Hubble distance. In this paper we develop a non-parametric method for extracting information about the perturbative growth factor from data pertaining to the luminosity or angular size distances. A comparison of the reconstructed density contrast with observations of large scale structure and gravitational lensing can help distinguish DE models such as the cosmological constant and quintessence from models based on modified gravity theories as well as models in which DE and DM are either unified, or interact directly. We show that for current SNe data, the linear growth factor at z = 0.3 can be constrained to 5%, and the linear growth rate to 6%. With future SNe data, such as expected from the JDEM mission, we may be able to constrain the growth factor to 2-3% and the growth rate to 3-4% at z = 0.3 with this unbiased, model-independent reconstruction method. For future BAO data which would deliver measurements of both the angular diameter distance and Hubble parameter, it should be possible to constrain the growth factor at z = 2.5 to 9%. These constraints grow tighter with the errors on the datasets. With a large quantity of data expected in the next few years, this method can emerge as a competitive tool for distinguishing between different models of dark energy.

It is proved that all thin-shell wormholes built from two identical regions of vacuum static, spherically symmetric space-times have a negative shell surface energy density in any scalar-tensor theory of gravity with a non-ghost massless scalar field and a non-ghost graviton.

A sudden small change in the second derivative of the inflaton potential can result in a universal local feature in the spectrum of primordial perturbations generated during inflation. The exact solution describing this feature \cite{minu} is characterized by a step in the spectral index modulated by characteristic oscillations and results in a large running of the spectral index localized over a few e-folds of scale. In this paper we confront this step-like feature with the 5 year WMAP results and demonstrate that it provides a better fit to this data than a featureless initial spectrum. If such a feature exists at all, then it should lie at sufficiently large scales $k_0 \lesssim 0.003 {\rm Mpc}^{-1}$ corresponding to $l \lesssim 40$. The sign of the effect is shown to correspond to the negative running of $n_s$ localized near this scale. This feature could arise as a result of a `mini-waterfall'-type fast second order phase transition experienced by an auxiliary heavy field during inflation, in a model similar to hybrid inflation (though for a different choice of parameters). If this is the case, then the auxiliary field should be positively coupled to the inflaton.

A new and powerful probe of the origin and evolution of structures in the Universe has emerged and been actively developed over the last decade. In the coming decade, non-Gaussianity, i.e., the study of non-Gaussian contributions to the correlations of cosmological fluctuations, will become an important probe of both the early and the late Universe. Specifically, it will play a leading role in furthering our understanding of two fundamental aspects of cosmology and astrophysics: (i) the physics of the very early universe that created the primordial seeds for large-scale structures, and (ii) the subsequent growth of structures via gravitational instability and gas physics at later times. To date, observations of fluctuations in the Cosmic Microwave Background (CMB) and the Large-Scale Structure of the Universe (LSS) have focused largely on the Gaussian contribution as measured by the two-point correlations (or the power spectrum) of density fluctuations. However, an even greater amount of information is contained in non-Gaussianity and a large discovery space therefore still remains to be explored. Many observational probes can be used to measure non-Gaussianity, including CMB, LSS, gravitational lensing, Lyman-alpha forest, 21-cm fluctuations, and the abundance of rare objects such as clusters of galaxies and high-redshift galaxies. Not only does the study of non-Gaussianity maximize the science return from a plethora of present and future cosmological experiments and observations, but it also carries great potential for important discoveries in the coming decade.

We prove that the stochastic and standard field-theoretical approaches produce exactly the same results for the amount of light massive scalar field fluctuations generated during inflation in the leading order of the slow-roll approximation. This is true both in the case for which this field is a test one and inflation is driven by another field, and the case for which the field plays the role of inflaton itself. In the latter case, in order to calculate the average of the mean square of the gauge-invariant inflaton fluctuation, the logarithm of the scale factor $a$ has to be used as the time variable in the Fokker-Planck equation in the stochastic approach. The implications of particle production during inflation for the second stage of inflation and for the moduli problem are also discussed. The case of a massless self-interacting test scalar field in a de Sitter background with a zero initial renormalized mean square is also considered in order to show how the stochastic approach can easily produce results corresponding to diagrams with an arbitrary number of scalar field loops in the field-theoretical approach (explicit results up to 4 loops inclusive are presented).

We introduce two new diagnostics of dark energy (DE). The first, Om, is a combination of the Hubble parameter and the cosmological redshift and provides a "null test" of dark energy being a cosmological constant. Namely, if the value of Om(z) is the same at different redshifts, then DE is exactly cosmological constant. The slope of Om(z) can differentiate between different models of dark energy even if the value of the matter density is not accurately known. For DE with an unevolving equation of state, a positive slope of Om(z) is suggestive of Phantom (w < -1) while a negative slope indicates Quintessence (w > -1). The second diagnostic, "acceleration probe"(q-probe), is the mean value of the deceleration parameter over a small redshift range. It can be used to determine the cosmological redshift at which the universe began to accelerate, again without reference to the current value of the matter density. We apply the "Om" and "q-probe" diagnostics to the Union data set of type Ia supernovae combined with recent data from the cosmic microwave background (WMAP5) and baryon acoustic oscillations.

We continue investigation of the hidden plane-mirror symmetry in the distribution of excursion sets in cosmic microwave background (CMB) temperature anisotropy maps, previously noticed in the three-year data of the Wilkinson microwave anisotropy probe (WMAP), using the WMAP 5 years maps. The symmetry is shown to have higher significance, \chi^2 < 1.7, for low multipoles \ell < 5, while disappearing at larger multipoles, \chi^2 > 3.5 for \ell > 10. The study of the sum and difference maps of temperature inhomogeneity regions, along with simulated maps, confirms its existence.The properties of these mirroring symmetries are compatible with those produced by the Sachs-Wolfe effect in the presence of an anomalously large component of horizon-size density perturbations, independent of one of the spatial coordinates, and/or a slab-like spatial topology of the Universe.

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

A model is developed in which the inflaton potential experiences a sudden small change in its second derivative (the effective mass of the inflaton). An exact treatment demonstrates that the resulting density perturbation has a quasi-flat power spectrum with a break in its slope (a step in n_s). The step in the spectral index is modulated by characteristic oscillations and results in large running of the spectral index localised over a few e-folds of scales. A field-theoretic model giving rise to such behaviour of the inflationary potential is based on a fast phase transition experienced by a second scalar field weakly coupled to the inflaton. Such a transition is similar to that which terminates inflation in the hybrid inflationary scenario. This scenario suggests that the observed running of the spectral index in the WMAP data may be caused by a fast second order phase transition which occurred during inflation.

We study inhomogeneities in the distribution of the excursion sets in the Cosmic Microwave Background (CMB) temperature maps obtained by the three years survey of the Wilkinson Microwave Anisotropy Probe (WMAP). At temperature thresholds |T|<90 \mu K, the distributions of the excursion sets with over 200 pixels are concentrated in two regions, nearly at the antipodes, with galactic coordinates l= 94^\circ.7, b= 34^\circ.4 and l= 279^\circ.8, b= -29^\circ.2. The centers of these two regions drift towards the equator when the temperature threshold is increased. The centers are located close to one of the vectors of \ell =3 multipole. The two patterns of the substructures in the distribution of the excursion sets are mirrored, with \chi^2=0.7-1.5. There is no obvious origin of this effect in the noise structure of WMAP, and there is no evidence for a dependence on the galactic cut. Would this effect be cosmological, it could be an indication of an anomalously large component of horizon-size density perturbations, independent of one of the spatial coordinates, and/or of a non-trivial slab-like spatial topology of the Universe.

We derive new constraints on the Hubble function H(phi) and subsequently on the inflationary potential V(phi) from WMAP 3-year data combined with the Sloan Luminous Red Galaxy survey (SDSS-LRG), using a new methodology which appears to be more generic, conservative and model-independent than in most of the recent literature, since it depends neither on the slow-roll approximation, nor on any extrapolation scheme for the potential beyond the observable e-fold range, nor on additional assumptions about initial conditions for the inflaton velocity. This last feature represents the main improvement of this work, and is made possible by the reconstruction of H(phi) prior to V(phi). Our results only rely on the assumption that within the observable range, corresponding to ~ 10 e-folds, inflation is not interrupted and the function H(phi) is smooth enough for being Taylor-expanded at order one, two or three. We conclude that the variety of potentials allowed by the data is still large. However, it is clear that the first two slow-roll parameters are really small while the validity of the slow-roll expansion beyond them is not established.

For higher-derivative f(R) gravity where R is the Ricci scalar, a class of models is proposed which produce viable cosmology different from the LambdaCDM one at recent times and satisfy cosmological, Solar system and laboratory tests. These models have both flat and de Sitter space-times as particular solutions in the absence of matter. Thus, a cosmological constant is zero in flat space-time, but appears effectively in a curved one for sufficiently large R. A 'smoking gun' for these models would be small discrepancy in values of the slope of the primordial perturbation power spectrum determined from galaxy surveys and CMB fluctuations. On the other hand, a new problem for dark energy models based on f(R) gravity is pointed which is connected with possible overproduction of new massive scalar particles (scalarons) arising in this theory in the very early Universe.

Two previously proposed conjectures--gravitational trans-Planckian particle creation in the expanding universe, and the existence of ultra-heavy stable particles with masses up to the Planck scale (wimpzillas)--are combined in a proposal for trans-Planckian particle creation of wimpzillas. This new scenario leads to a huge enhancement in their production compared to mechanisms put forward earlier. As a result, it requires the trans-Planckian particle creation parameter to be rather small to avoid overproduction of such particles, much less than that is required for observable effects in the primordial perturbation spectrum. This ensures also that wimpzillas are mainly created at the end of primordial inflation. Conditions under which trans-Planckian wimpzillas can constitute the present dark matter are determined.

From photometric observations of elliptical galaxies, among which are both radio galaxies and radio-quiet objects, an investigation was carried out of the relationship `redshift -- age of the stellar system' $(\Delta z/\Delta t)$. By means of this relationship cosmological parameters $H(z)$ and $\Omega_\Lambda$ are estimated. Ages of stellar systems are determined within the framework of evolution models of synthetic spectra PEGASE and GISSEL. This approach can be considered as time study of objects of the early Universe independent of other cosmological models. Construction of a pooled sample is described, containing 220 objects from different populations of elliptical galaxies, for which an analysis of the upper limit of the age of formation of a stellar system was performed. These data were used to estimate the boundaries of determination of the cosmological parameters $H_0$ and $\Lambda$--term: $H_0=72\pm10$ and $\Omega_\Lambda=0.8\pm0.1$ in the model GISSEL and $H_0=53\pm10$, and $\Omega_\Lambda=0.8\pm0.1$ in the model PEGASE.

We propose a non-parametric method of smoothing supernova data over redshift using a Gaussian kernel in order to reconstruct important cosmological quantities including H(z) and w(z) in a model independent manner. This method is shown to be successful in discriminating between different models of dark energy when the quality of data is commensurate with that expected from the future SuperNova Acceleration Probe (SNAP). We find that the Hubble parameter is especially well-determined and useful for this purpose. The look back time of the universe may also be determined to a very high degree of accuracy (\lleq 0.2 %) in this method. By refining the method, it is also possible to obtain reasonable bounds on the equation of state of dark energy. We explore a new diagnostic of dark energy-- the `w-probe'-- which can be calculated from the first derivative of the data. We find that this diagnostic is reconstructed extremely accurately for different reconstruction methods even if \Omega_m is marginalized over. The w-probe can be used to successfully distinguish between $\Lambda$CDM and other models of dark energy to a high degree of accuracy.

Primordial fluctuations in inflationary cosmology acquire classical properties through decoherence when their wavelengths become larger than the Hubble scale. Although decoherence is effective, it is not complete, so a significant part of primordial correlations remains up to the present moment. We address the issue of the pointer states which provide a classical basis for the fluctuations with respect to the influence by an environment (other fields). Applying methods from the quantum theory of open systems (the Lindblad equation), we show that this basis is given by narrow Gaussians that approximate eigenstates of field amplitudes. We calculate both the von Neumann and linear entropy of the fluctuations. Their ratio to the maximal entropy per field mode defines a degree of partial decoherence in the entropy sense. We also determine the time of partial decoherence making the Wigner function positive everywhere which, for super-Hubble modes during inflation, is virtually independent of coupling to the environment and is only slightly larger than the Hubble time. On the other hand, assuming a representative environment (a photon bath), the decoherence time for sub-Hubble modes is finite only if some real dissipation exists.

We present here some recent results concerning scalar-tensor Dark Energy models. These models are very interesting in many respects: they allow for a consistent phantom phase, the growth of matter perturbations is modified. Using a systematic expansion of the theory at low redshifts, we relate the possibility to have phantom like DE to solar system constraints.

We reconstruct dark energy properties from two complementary supernova datasets -- the newly released Gold+HST sample and SNLS. The results obtained are consistent with standard $\Lambda$CDM model within $2\sigma$ error bars although the Gold+HST data favour evolving dark energy slightly more than SNLS. Using complementary data from baryon acoustic oscillations and the cosmic microwave background to constrain dark energy, we find that our results in this case are strongly dependent on the present value of the matter density $\Omega_m$. Consequently, no firm conclusions regarding constancy or variability of dark energy density can be drawn from these data alone unless the value of $\Omega_m$ is known to an accuracy of a few percent. However, possible variability is significantly restricted if this data is used in conjunction with supernova data.

It is proved that no wormholes can be formed in viable scalar-tensor models of dark energy admitting its phantom-like ($w < -1$) behaviour in cosmology, even in the presence of electric or magnetic fields, if the non-minimal coupling function $f(\Phi)$ is everywhere positive and the scalar field $\Phi$ itself is not a ghost. Some special static, spherically symmetric wormhole solutions may exist if $f(\Phi)$ is allowed to reach zero or to become negative, so that the effective gravitational constant becomes negative in some region making the graviton a ghost. If $f$ remains non-negative, such solutions require severe fine tuning and a very peculiar kind of model. If $f < 0$ is allowed, it is argued (and confirmed by previous investigations) that such solutions are generically unstable under non-static perturbations, the instability appearing right near transition surfaces to negative $f$.

We consider the viability of dark energy (DE) models in the framework of the scalar-tensor theory of gravity, including the possibility to have a phantom DE at small redshifts $z$ as admitted by supernova luminosity-distance data. For small $z$, the generic solution for these models is constructed in the form of a power series in $z$ without any approximation. Necessary constraints for DE to be phantom today and to cross the phantom divide line $p=-\rho$ at small $z$ are presented. Considering the Solar System constraints, we find for the post-Newtonian parameters that $\gamma_{PN}<1$ and $\gamma_{PN,0}\approx 1$ for the model to be viable, and $\beta_{PN,0}>1$ (but very close to 1) if the model has a significantly phantom DE today. However, prospects to establish the phantom behaviour of DE are much better with cosmological data than with Solar System experiments. Earlier obtained results for a $\Lambda$-dominated universe with the vanishing scalar field potential are extended to a more general DE equation of state confirming that the cosmological evolution of these models rule them out. Models of currently fantom DE which are viable for small $z$ can be easily constructed with a constant potential; however, they generically become singular at some higher $z$. With a growing potential, viable models exist up to an arbitrary high redshift.

This review summarizes recent attempts to reconstruct the expansion history of the Universe and to probe the nature of dark energy. Reconstruction methods can be broadly classified into parametric and non-parametric approaches. It is encouraging that, even with the limited observational data currently available, different approaches give consistent results for the reconstruction of the Hubble parameter $H(z)$ and the effective equation of state $w(z)$ of dark energy. Model independent reconstruction using current data allows for modest evolution of dark energy density with redshift. However, a cosmological constant (= dark energy with a constant energy density) remains an excellent fit to the data. Some pitfalls to be guarded against during cosmological reconstruction are summarized and future directions for the model independent reconstruction of dark energy are explored.

Archeops is a balloon--borne experiment widely inspired by the Planck satellite and by its High Frequency Instrument (HFI). It is mainly dedicated to measure the Cosmic Microwave Background (CMB) temperature anisotropies at high angular resolution (about 12 arcminutes) over a large fraction of the sky (around 30 %) in the millimetre and submillimetre range at 143, 217, 353 and 545 GHz. Further, the Archeops 353 GHz channel consists of three pairs of polarized sensitive bolometers designed to detect the polarized diffuse emission of Galactic dust. We present in this paper the update of the instrumental setup as well as the inflight performance for the last Archeops flight campaign in February 2002 from Kiruna (Sweden). We also describe the processing and analysis of the Archeops time ordered data for that campaign which lead to the measurement of the CMB anisotropies power spectrum in the multipole range l=10-700 (Benoit et al. 2003a, Tristram et al. 2005) and to the first measurement of the dust polarized emission at large angular scales and its polarized. We present maps of 30 % of the sky of the Galactic emission, including the Galactic plane, in the four Archeops channels at 143, 217, 353 and 545 GHz and maps of the CMB anisotropies at 143 and 217 GHz. These are the firstever available sub--degree resolution maps in the millimetre and submillimetre range of the large angular-scales Galactic dust diffuse emission and CMB temperature anisotropies respectively.

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.

Presented is the exact solution of the problem of finding a potential of an inflaton scalar field for which adiabatic perturbations generated during a de Sitter (inflationary) stage in the early Universe have the exactly flat (or, the Harrison-Zeldovich) initial spectrum. This solution lies outside the scope of the slow-roll approximation and higher-order corrections to it. The potential found depends on two arbitrary physical constants, one of those determines an amplitude of perturbations. For small (zero) values of the other constant, a long (infinite) inflationary stage with slow rolling of the inflaton field exists.

We reconstruct the equation of state $w(z)$ of dark energy (DE) using a recently released data set containing 172 type Ia supernovae without assuming the prior $w(z) \geq -1$ (in contrast to previous studies). We find that dark energy evolves rapidly and metamorphoses from dust-like behaviour at high $z$ ($w \simeq 0$ at $z \sim 1$) to a strongly negative equation of state at present ($w \lleq -1$ at $z \simeq 0$). Dark energy metamorphosis appears to be a robust phenomenon which manifests for a large variety of SNe data samples provided one does not invoke the weak energy prior $\rho + p \geq 0$. Invoking this prior considerably weakens the rate of growth of $w(z)$. These results demonstrate that dark energy with an evolving equation of state provides a compelling alternative to a cosmological constant if data are analysed in a prior-free manner and the weak energy condition is not imposed by hand.

In a recent paper (astro-ph/0311364) Alam et al argued that the SNe data of Tonry et al 2003 and Barris et al 2003 appear to favour DE which evolves in time, provided no other priors are invoked. (The effect of invoking priors such as the age of the Universe, the values of $H_0$ and $\Omega_m$ and CMB/LSS observations could modify this conclusion, as demonstrated in astro-ph/0403687 and other recent papers.) The approach adopted by Alam et al to reconstruct the properties of DE was severely (and, as we shall show below -- unfairly) criticized by J\"onsson et al in astro-ph/0404468. In this paper we re-examine the parametrisation used in astro-ph/0311364 and show that, contrary to the claims of J\"onsson et al, the results obtained from this reconstruction are robust and therefore representative of the true nature of dark energy.

We investigate the behaviour of dark energy using the recently released supernova data of Riess et al ~(2004) and a model independent parameterization for dark energy (DE). We find that, if no priors are imposed on $\Omega_{0m}$ and $h$, DE which evolves with time provides a better fit to the SNe data than $\Lambda$CDM. This is also true if we include results from the WMAP CMB data. From a joint analysis of SNe+CMB, the best-fit DE model has $w_0 < -1$ at the present epoch and the transition from deceleration to acceleration occurs at $z_T = 0.39 \pm 0.03$. However, DE evolution becomes weaker if the $\Lambda$CDM based CMB results $\Omega_{0m} = 0.27 \pm 0.04$, $h = 0.71 \pm 0.06$ are incorporated in the analysis. In this case, $z_T = 0.57 \pm 0.07$. Our results also show that the extent of DE evolution is sensitive to the manner in which the supernova data is sampled.

We investigate four-dimensional effective theories which are obtained by dimensional reduction of multidimensional cosmological models with factorizable geometry and consider the interaction between conformal excitations of the internal space (geometrical moduli excitations) and Abelian gauge fields. It is assumed that the internal space background can be stabilized by minima of an effective potential. The conformal excitations over such a background have the form of massive scalar fields (gravitational excitons) propagating in the external spacetime. We discuss cosmological and astrophysical implications of the interaction between gravexcitons and four-dimensional photons as well as constraints arising on multidimensional models of the type considered in our paper. In particular, we show that due to the experimental bounds on the variation of the fine structure constant, gravexcitons should decay before nucleosynthesis starts. For a successful nucleosynthesis the masses of the decaying gravexcitons should be m>10^4 GeV. Furthermore, we discuss the possible contribution of gravexcitons to UHECR. It is shown that, at energies of about 10^{20}eV, the decay length of gravexcitons with masses m>10^4 GeV is very small, but that for m <10^2 GeV it becomes much larger than the Greisen-Zatsepin-Kuzmin cut-off distance. Finally, we investigate the possibility for gravexciton-photon oscillations in strong magnetic fields of astrophysical objects. The corresponding estimates indicate that even the high magnetic field strengths of magnetars are not sufficient for an efficient and copious production of gravexcitons.

The coming few years are likely to witness a dramatic increase in high quality Sn data as current surveys add more high redshift supernovae to their inventory and as newer and deeper supernova experiments become operational. Given the current variety in dark energy models and the expected improvement in observational data, an accurate and versatile diagnostic of dark energy is the need of the hour. This paper examines the Statefinder diagnostic in the light of the proposed SNAP satellite which is expected to observe about 2000 supernovae per year. We show that the Statefinder is versatile enough to differentiate between dark energy models as varied as the cosmological constant on the one hand, and quintessence, the Chaplygin gas and braneworld models, on the other. Using SNAP data, the Statefinder can distinguish a cosmological constant ($w=-1$) from quintessence models with $w \geq -0.9$ and Chaplygin gas models with $\kappa \leq 15$ at the $3\sigma$ level if the value of $\om$ is known exactly. The Statefinder gives reasonable results even when the value of $\om$ is known to only $\sim 20%$ accuracy. In this case, marginalizing over $\om$ and assuming a fiducial LCDM model allows us to rule out quintessence with $w \geq -0.85$ and the Chaplygin gas with $\kappa \leq 7$ (both at $3\sigma$). These constraints can be made even tighter if we use the Statefinders in conjunction with the deceleration parameter. The Statefinder is very sensitive to the total pressure exerted by all forms of matter and radiation in the universe. It can therefore differentiate between dark energy models at moderately high redshifts of $z \lleq 10$.

We present the first determination of the Galactic polarized emission at 353 GHz by Archeops. The data were taken during the Arctic night of February 7, 2002 after the balloon--borne instrument was launched by CNES from the Swedish Esrange base near Kiruna. In addition to the 143 GHz and 217 GHz frequency bands dedicated to CMB studies, Archeops had one 545 GHz and six 353 GHz bolometers mounted in three polarization sensitive pairs that were used for Galactic foreground studies. We present maps of the I, Q, U Stokes parameters over 17% of the sky and with a 13 arcmin resolution at 353 GHz (850 microns). They show a significant Galactic large scale polarized emission coherent on the longitude ranges [100, 120] and [180, 200] deg. with a degree of polarization at the level of 4-5%, in agreement with expectations from starlight polarization measurements. Some regions in the Galactic plane (Gem OB1, Cassiopeia) show an even stronger degree of polarization in the range 10-20%. Those findings provide strong evidence for a powerful grain alignment mechanism throughout the interstellar medium and a coherent magnetic field coplanar to the Galactic plane. This magnetic field pervades even some dense clouds. Extrapolated to high Galactic latitude, these results indicate that interstellar dust polarized emission is the major foreground for PLANCK-HFI CMB polarization measurement.

We explore the fate of the universe given the possibility that the density associated with `dark energy' may decay slowly with time. Decaying dark energy is modeled by a homogeneous scalar field which couples minimally to gravity and whose potential has {\em at least one} local quadratic maximum. Dark energy decays as the scalar field rolls down its potential, consequently the current acceleration epoch is a transient. We examine two models of decaying dark energy. In the first, the dark energy potential is modeled by an analytical form which is generic close to the potential maximum. The second potential is the cosine, which can become negative as the field evolves, ensuring that a spatially flat universe collapses in the future. We examine the feasibility of both models using observations of high redshift type Ia supernovae. A maximum likelihood analysis is used to find allowed regions in the $\lbrace m, \phi_0\rbrace$ plane ($m$ is the tachyon mass modulus and $\phi_0$ the initial scalar field value; $m\sim H_0$ and $\phi_0\sim M_P$ by order of magnitude). For the first model, the time for the potential to drop to half its maximum value is larger than $\sim 8$ Gyrs. In the case of the cosine potential, the time left until the universe collapses is always greater than $\sim 18$ Gyrs (both estimates are presented for $\om = 0.3$, $m/H_0 \sim 1$, $H_0 \simeq 70$ km/sec/Mpc, and at the 95.4% confidence level).

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.

We introduce a new cosmological diagnostic pair $\lbrace r,s\rbrace$ called Statefinder. The Statefinder is dimensionless and, like the Hubble and deceleration parameters $H(z)$ and $q(z)$, is constructed from the scale factor of the Universe and its derivatives only. The parameter $r(z)$ forms the next step in the hierarchy of geometrical cosmological parameters used to study the Universe after $H$ and $q$, while the parameter $s(z)$ is a linear combination of $q$ and $r$ chosen in such a way that it does not depend upon the dark energy density $\Omega_X(z)$. The Statefinder pair $\lbrace r,s\rbrace$ is algebraically related to the the dark energy pressure-to-energy ratio $w=p/\epsilon$ and its time derivative, and sheds light on the nature of dark energy/quintessence. Its properties allow to usefully differentiate between different forms of dark energy with constant and variable $w$, including a cosmological constant ($w = -1$). The Statefinder pair can be determined to very good accuracy from a SNAP type experiment.

We analyze the cosmological constraints that Archeops places on adiabatic cold dark matter models with passive power-law initial fluctuations. Because its angular power spectrum has small bins in l and large l coverage down to COBE scales, Archeops provides a precise determination of the first acoustic peak in terms of position at multipole l_peak=220 +- 6, height and width. An analysis of Archeops data in combination with other CMB datasets constrains the baryon content of the Universe, Omega(b)h^2 = 0.022 (+0.003,-0.004), compatible with Big-Bang nucleosynthesis and with a similar accuracy. Using cosmological priors obtainedfrom recent non-CMB data leads to yet tighter constraints on the total density, e.g. Omega(tot)=1.00 (+0.03,-0.02) using the HST determination of the Hubble constant. An excellent absolute calibration consistency is found between Archeops and other CMB experiments, as well as with the previously quoted best fit model.The spectral index n is measured to be 1.04 (+0.10,-0.12) when the optical depth to reionization, tau, is allowed to vary as a free parameter, and 0.96 (+0.03,-0.04) when tau is fixed to zero, both in good agreement with inflation.

We present a determination by the Archeops experiment of the angular power spectrum of the cosmic microwave background anisotropy in 16 bins over the multipole range l=15-350. Archeops was conceived as a precursor of the Planck HFI instrument by using the same optical design and the same technology for the detectors and their cooling. Archeops is a balloon-borne instrument consisting of a 1.5 m aperture diameter telescope and an array of 21 photometers maintained at ~100 mK that are operating in 4 frequency bands centered at 143, 217, 353 and 545 GHz. The data were taken during the Arctic night of February 7, 2002 after the instrument was launched by CNES from Esrange base (Sweden). The entire data cover ~ 30% of the sky.This first analysis was obtained with a small subset of the dataset using the most sensitive photometer in each CMB band (143 and 217 GHz) and 12.6% of the sky at galactic latitudes above 30 degrees where the foreground contamination is measured to be negligible. The large sky coverage and medium resolution (better than 15 arcminutes) provide for the first time a high signal-to-noise ratio determination of the power spectrum over angular scales that include both the first acoustic peak and scales probed by COBE/DMR. With a binning of Delta(l)=7 to 25 the error bars are dominated by sample variance for l below 200. A companion paper details the cosmological implications.

We consider observational constraints on creation of particles induced by hypothetical trans-Planckian effects during the current stage of the Universe expansion. We show that compatibility with the diffuse gamma-ray background measured by the EGRET experiment strongly restricts this creation. In particular, it rules out the possibility to detect signatures of such short distance effects in anisotropies of the cosmic microwave background radiation. On the other hand, a possibility that some part of the ultra-high energy cosmic rays originates from new trans-Planckian physics remains open.

We consider scalar Born-Infeld type theories with arbitrary potentials V(T) of a scalar field T. We find that for models with runaway potentials V(T) the generic inhomogeneous solutions after a short transient stage can be very well approximated by the solutions of a Hamilton-Jacobi equation that describes free streaming wave front propagation. The analytic solution for this wave propagation shows the formation of caustics with multi-valued regions beyond them. We verified that these caustics appear in numerical solutions of the original scalar BI non-linear equations. Our results include the scalar BI model with an exponential potential, which was recently proposed as an effective action for the string theory tachyon in the approximation where high-order spacetime derivatives of T are truncated. Since the actual string tachyon dynamics contain derivatives of all orders, the tachyon BI model with an exponential potential becomes inadequate when the caustics develop because high order spatial derivatives of T become divergent. BI type tachyon theory with a potential decreasing at large T could have interesting cosmological applications because the tachyon field rolling towards its ground state at infinity acts as pressureless dark matter. We find that inhomogeneous cosmological tachyon fluctuations rapidly grow and develop multiple caustics. Any considerations of the role of the tachyon field in cosmology will have to involve finding a way to predict the behavior of the field at and beyond these caustics.

We consider the evolution of FRW cosmological models and linear perturbations of tachyon matter rolling towards a minimum of its potential. The tachyon coupled to gravity is described by an effective 4d field theory of string theory tachyon. In the model where a tachyon potential $V(T)$ has a quadratic minimum at finite value of the tachyon field $T_0$ and $V(T_0)=0$, the tachyon condensate oscillates around its minimum with a decreasing amplitude. It is shown that its effective equation of state is $p=-\epsilon/3$. However, linear inhomogeneous tachyon fluctuations coupled to the oscillating background condensate are exponentially unstable due to the effect of parametric resonance. In another interesting model, where tachyon potential exponentially approaches zero at infinity of $T$, rolling tachyon condensate in an expanding universe behaves as pressureless fluid. Its linear fluctuations coupled with small metric perturbations evolve similar to these in the pressureless fluid. However, this linear stage changes to a strongly non-linear one very early, so that the usual quasi-linear stage observed at sufficiently large scales in the present Universe may not be realized in the absence of the usual particle-like cold dark matter.

We discuss the structure of parametric resonance which occurs in the process of reheating after inflation with two interacting scalar fields. It is found that, for the case of a not too large coupling constant, a quasi-homogeneous part of the second, initially subdominant scalar field may not be neglected due to its stochastic growth during inflation. This fact has strong consequences for the reheating stage: dynamics of background fields becomes chaotic after inflation -- it consists of subsequent chaotic and regular eras, and a tachyonic instability for inhomogeneous perturbations arises in the quantum reheating problem. This instability may pose a problem for the standard reheating scenario. In order to avoid it, the coupling constant should be either sufficiently small, or very large.

We show that the amplitude of the second acoustic peak in the newly released BOOMERANG-98 and MAXIMA-I data is compatible with the standard primordial nucleosynthesis and with the locally broken-scale-invariant matter power spectrum suggested by recent measurements of the power spectrum in the range $20 - 200$ Mpc. If the slope of matter density perturbations on large scales is $n \approx 1$, the Hubble constant is $0.5< h < 0.75$, and r.m.s. mass fluctuations at 8 Mpc are $0.65 \le \sigma_8 \le 0.75$, then for a Universe approximately 14 Gyr old our best fit within the nucleosynthesis bound $\Omega_Bh^2= 0.019 \pm 0.0024$ requires $0.3 \le \Omega_m \le 0.5$. Cluster abundances further constraint the matter density to be $\Omega_m \approx 0.3$. The CMB data alone are not able to determine the detailed form of the matter power spectrum in the range $0.03<k<0.06$ \hmpc where deviations from the scale-invariant spectrum are expected to be most significant, but they do not contradict the existence of the previously claimed peak at $k\sim 0.05$ \hmpc, and a depression at $k\sim 0.035$ \hmpc.

We study cosmological perturbations generated from quantum fluctuations in multi-field inflationary scenarios in generalized Einstein theories, taking both adiabatic and isocurvature modes into account. In the slow-roll approximation, explicit closed-form long-wave solutions for field and metric perturbations are obtained by the analysis in the Einstein frame. Since the evolution of fluctuations depends on specific gravity theories, we make detailed investigations based on analytic and numerical approaches in four generalized Einstein theories: the Jordan-Brans-Dicke (JBD) theory, the Einstein gravity with a non-minimally coupled scalar field, the higher-dimensional Kaluza-Klein theory, and the $R+R^2$ theory with a non-minimally coupled scalar field. We find that solutions obtained in the slow-roll approximation show good agreement with full numerical results except around the end of inflation. Due to the presence of isocurvature perturbations, the gravitational potential $\Phi$ and the curvature perturbations ${\cal R}$ and $\zeta$ do not remain constant on super-horizon scales. In particular, we find that negative non-minimal coupling can lead to strong enhancement of ${\cal R}$ in both the Einstein and higher derivative gravity, in which case it is difficult to unambiguously decompose scalar perturbations into adiabatic and isocurvature modes during the whole stage of inflation.

It is investigated if predictions of the inflationary scenario regarding spectra of scalar and tensor perturbations generated from quantum vacuum fluctuations are robust with respect to a modification of the dispersion law for frequencies beyond the Planck scale. For a large class of such modifications of special and general relativity, for which the WKB condition is not violated at super-high frequencies, the predictions remain unchanged. The opposite possibility is excluded by the absence of large amount of created particles due to the present Universe expansion. Creation of particles in the quantum state minimizing the energy density of a given mode at the moment of Planck boundary crossing is prohibited by the latter argument, too (contrary to creation in the adiabatic vacuum state which is very small now).

Observations of high-redshift supernovae indicate that the universe is accelerating. Here we present a {\em model-independent} method for estimating the form of the potential $V(\phi)$ of the scalar field driving this acceleration, and the associated equation of state $w_\phi$. Our method is based on a versatile analytical form for the luminosity distance $D_L$, optimized to fit observed distances to distant supernovae and differentiated to yield $V(\phi)$ and $w_\phi$. Our results favor $w_\phi\simeq -1$ at the present epoch, steadily increasing with redshift. A cosmological constant is consistent with our results.

The appearance of peaks in various primordial fluctuation Fourier power spectra is a generic prediction of the inflationary scenario. We investigate whether future experiments, in particular the satellite experiment PLANCK, will be able to detect the possible appearance of these peaks in the B-mode polarization multipole power spectrum. This would yield a conclusive proof of the presence of a primordial background of gravitational waves.

Recent observations of Type 1a supernovae indicating an accelerating universe have once more drawn attention to the possible existence, at the present epoch, of a small positive Lambda-term (cosmological constant). In this paper we review both observational and theoretical aspects of a small cosmological Lambda-term. We discuss the current observational situation focusing on cosmological tests of Lambda including the age of the universe, high redshift supernovae, gravitational lensing, galaxy clustering and the cosmic microwave background. We also review the theoretical debate surrounding Lambda: the generation of Lambda in models with spontaneous symmetry breaking and through quantum vacuum polarization effects -- mechanisms which are known to give rise to alarge value of Lambda hence leading to the `cosmological constant problem'. More recent attempts to generate a small cosmological constant at the present epoch using either field theoretic techniques, or by modeling a dynamical Lambda-term by scalar fields are also extensively discussed. Anthropic arguments favouring a small cosmological constant are briefly reviewed. A comprehensive bibliography of recent work on Lambda is provided.

Gravitons produced from quantum vacuum fluctuations during an inflationary stage in the early Universe have zero entropy as far as they reflect the time evolution (squeezing) of a pure state, their large occupation number notwithstanding. A non-zero entropy of the gravitons (classical gravitational waves (GW) after decoherence) can be obtained through coarse graining. The latter has to be physically justified {\it and} should not contradict observational constraints. We propose two ways of coarse graining for which the fixed temporal phase of each Fourier mode of the GW background still remains observable: one based on quantum entanglement, and another one following from the presence of a secondary GW background. The proposals are shown to be mutually consistent. They lead to the result that the entropy of the primordial GW background is significantly smaller than it was thought earlier. The difference can be ascribed to the information about the regular (inflationary) initial state of the Universe which is stored in this background and which reveals itself, in particular, in the appearance of primordial peaks (acoustic peaks in the case of scalar perturbations) in the multipole spectra of the CMB temperature anisotropy and polarization.

The present acceleration of the Universe strongly indicated by recent observational data can be modeled in the scope of a scalar-tensor theory of gravity. We show that it is possible to determine the structure of this theory (the scalar field potential and the functional form of the scalar-gravity coupling) along with the present density of dustlike matter from the following two observable cosmological functions: the luminosity distance and the linear density perturbation in the dustlike matter component as functions of redshift. Explicit results are presented in the first order in the small inverse Brans-Dicke parameter 1/omega.

The existence of a positive and possibly varying Lambda-term opens a much wider field of possibilities for the future of our Universe than it was usually thought before. Definite predictions may be made for finite (though very large) intervals of time only, as well as in other branches of science. In particular, our Universe will continue to expand as far as the Lambda-term remains positive and does not decay to other forms of matter, even if the Universe is closed. Two new effects due to the presence of a constant Lambda-term are discussed: reversal of a sign of the redshift change with time for sufficiently close objects and inaccessibility of sufficiently distant objects in the Universe for us. A number of more distant and speculative possibilities for the future evolution of the Universe is listed including hitting a space-time singularity during an expansion phase. Finally, in fantastically remote future, a part of our Universe surrounding us can become supercurved and superdense due to various quantum-gravitational effects. This returns us to the past, to the origin of our Universe from a superdense state about 14 Gy ago. According to the inflationary scenario, this state was almost maximally symmetric (de Sitter-like). Though this scenario seems to be sufficient for the explanation of observable properties of the present Universe, and its predictions have been confirmed by observations, the question of the origin of the initial de Sitter (inflationary) state itself remains open. A number of conjectures regarding the very origin of our Universe, ranging from "creation from nothing" to "creation from anything", are discussed.

The distribution of Abell clusters of galaxies is analysed to study the regularity of the supercluster-void network. A new geometric method sensitive to the regularity of the location of clusters is applied. We find that the supercluster-void network resembles a cubical lattice over the whole space investigated. The distribution of rich superclusters is not isotropic: along the main axis of the network it is periodic with a step of length about 130 Mpc (for Hubble constant h=1), whereas along the diagonal of the network the period is larger. This large-scale inhomogeneity is compatible with recent CMB data.

We consider the primordial gravitational wave (GW) background in a class of spatially-flat inflationary cosmological models with cold dark matter (CDM), a cosmological constant, and a broken-scale-invariant (BSI) steplike primordial (initial) spectrum of adiabatic perturbations produced in an exactly solvable inflationary model where the inflaton potential has a rapid change of its first derivative at some point. In contrast to inflationary models with a scale-free initial spectrum, these models may have a GW power spectrum whose amplitude (though not its shape) is arbitrary for fixed amplitude and shape of the adiabatic perturbations power spectrum. In the presence of a positive cosmological constant, the models investigated here possess the striking property that a significant part of the large-angle CMB temperature anisotropy observed in the COBE experiment is due to primordial GW. Confronting them with existing observational data on CMB angular temperature fluctuations, galaxy clustering and peculiar velocities of galaxies, we find that for the best parameter values Omega_Lambda=0.7 and h=0.7, the GW contribution to the CMB anisotropy can be as large as that of the scalar fluctuations.

We calculate the mean power spectrum of galaxies using published power spectra of galaxies and clusters of galaxies. On small scales we use the power spectrum derived from the 2-dimensional distribution of APM galaxies, on large scales we use power spectra derived from 3-dimensional data for galaxy and cluster samples. Spectra are reduced to real space and to the amplitude of the power spectrum of APM galaxies. Available data indicate the presence of two different populations in the nearby Universe. Clusters of galaxies sample a relatively large region in the Universe where rich, medium and poor superclusters are well represented. Their mean power spectrum has a spike on scale 120 h^{-1}Mpc, followed by an approximate power-law spectrum of index n = -1.9 towards small scales. The power spectrum found from LCRS and IRAS 1.2 Jy surveys is flatter around the maximum, which may represent regions of the Universe with medium-rich and poor superclusters.

Observed power spectrum of matter found in Papers I and II is compared with analytical power spectra. Spatially flat cold and mixed dark matter models with cosmological constant and open models are considered. The primordial power spectrum of matter is determined using the power spectrum of matter and the transfer functions of analytical models. The primordial power spectrum has a break in amplitude. We conclude that a scale-free primordial power spectrum is excluded if presently available data on the distribution of clusters and galaxies represent the true mass distribution of the Universe.

A new method to determine the bias parameter of galaxies relative to matter is suggested. The method is based on the assumption that gravity is the dominating force which determines the formation of the structure in the Universe. Due to gravitational instability the galaxy formation is a threshold process: in low-density environments galaxies do not form and matter remains in primordial form. We investigate the influence of the presence of void and clustered populations to the power spectrum of matter and galaxies. The power spectrum of galaxies is similar to the power spectrum of matter; the fraction of total matter in the clustered population determines the difference between amplitudes of fluctuations of matter and galaxies, i.e. the bias factor. To determine the fraction of matter in voids and clustered population we perform numerical simulations. The fraction of matter in galaxies at the present epoch is found using a calibration through the sigma_8 parameter.

Though predictions of the simplest inflationary cosmological models with cold dark matter, flat space and approximately flat initial spectrum of adiabatic perturbations are remarkably close to observational data, we have to go beyond them and to introduce new physics not yet discovered in laboratories to account for all data. Two extensions of these models which seem to be the most actual at present time are discussed. The first one is the possibility that we are living at the beginning of a new inflation-like era. Then classical cosmological tests, like the luminosity distance or the angular size of distant objects as functions of redshift, as well as the behaviour of density perturbations in a dustlike matter component including baryons as a function of redshift, can provide information sufficient for the unambiguous determination of an effective potential of a corresponding present inflaton scalar field. The second, unrelated extension is a possibility of broken-scale-invariant cosmological models which have localized steps or spikes in the primordial perturbation spectrum. These features can be produced by fast phase transitions in physical fields other than an inflaton field in the early Universe during inflation and not far from the end of it. At present, it seems that the only scale in the spectrum around which we might see something of this type is $k=0.05 h Mpc^{-1}$.

It is shown that if a variable cosmological term in the present Universe is described by a scalar field with minimal coupling to gravity and with some phenomenological self-interaction potential $V(\varphi)$, then this potential can be unambiguously determined from the following observational data: either from the behaviour of density perturbations in dustlike matter component as a function of redshift (given the Hubble constant additionally), or from the luminosity distance as a function of redshift (given the present density of dustlike matter in terms of the critical one).

We consider the instability of fluctuations in an oscillating scalar field which obeys the Sine-Gordon equation. We present simple closed-form analytic solutions describing the parametric resonance in the Sine-Gordon model. The structure of the resonance differs from that obtained with the Mathieu equation which is usually derived with the small angle approximation to the equation for fluctuations. The results are applied to axion cosmology, where the oscillations of the classical axion field, with a Sine-Gordon self-interaction potential, constitute the cold dark matter of the universe. When the axion misalignment angle at the QCD epoch, $\theta_0$, is small, the parametric resonance of the axion fluctuations is not significant. However, in regions of larger $\theta_0$ where axion miniclusters would form, the resonance may be important. As a result, axion miniclusters may disintegrate into finer, denser clumps. We also apply the theory of Sine-Gordon parametric resonance to reheating in the Natural Inflation scenario. The decay of the inflaton field due to the self-interaction alone is ineffective, but a coupling to other bosons can lead to preheating in the broad resonance regime. Together with the preheating of fermions, this can alter the reheating scenario for Natural Inflation.

Two extensions of ideas lying in the basis of the inflationary scenario of the early Universe and their effect on the large scale structure of the present-day Universe are discussed. The first of them is the possibility of fast phase transitions in physical fields other than an inflaton scalar field during inflation and nor far from the end of it. This results in the appearance of specific features in the inflaton effective potential which, in turn, lead to the generation of localized spikes in the primordial perturbation spectrum. At present, there seems to exist one scale in the spectrum, $k=0.05h$ Mpc$^{-1}$, around which we might see something of this type. The second one is the possibility that we are living at the beginning of a new inflation-like era now. Then observations of clustering of high-redshift objects can provide information sufficient for the unambiguous determination of the effective potential of a corresponding present inflaton scalar field.

The behaviour of quantum metric perturbations produced during inflation is considered at the stage after the second Hubble radius crossing. It is shown that the classical correlation between amplitude and momentum of a perturbation mode, previously shown to emerge in the course of an effective quantum-to-classical transition, is maintained for a sufficiently long time, and we present the explicit form in which it takes place using the Wigner function. We further show with a simple diffraction experiment that quantum interference, non-expressible in terms of a classical stochastic description of the perturbations, is essentially suppressed. Rescattering of the perturbations leads to a comparatively slow decay of this correlation and to a complete stochastization of the system.

We show how to determine the typical phase space volume $\Gamma$ for primordial gravitational waves produced during an inflationary stage, which is invariant under squeezing. An expression for $\Gamma$ is found in the long wavelength regime. The quasi-classical entropy of a pure vacuum initial state defined as the logarithm of $\Gamma$ modulo a constant remains zero in spite of the generation of fluctuations (creation of real gravitons).

According to the inflationary scenario for the very early Universe, all inhomogeneities in the Universe are of genuine quantum origin. On the other hand, looking at these inhomogeneities and measuring them, clearly no specific quantum mechanical properties are observed. We show how the transition from their inherent quantum gravitational nature to classical behaviour comes about -- a transition whereby none of the successful quantitative predictions of the inflationary scenario for the present-day universe is changed. This is made possible by two properties. First, the quantum state for the spacetime metric perturbations produced by quantum gravitational effects in the early Universe becomes very special (highly squeezed) as a result of the expansion of the Universe (as long as the wavelength of the perturbations exceeds the Hubble radius). Second, decoherence through the environment distinguishes the field amplitude basis as being the pointer basis. This renders the perturbations presently indistinguishable from stochastic classical inhomogeneities.

A class of spatially-flat models with cold dark matter (CDM), a cosmological constant and a broken-scale-invariant (BSI) steplike primordial (initial) spectrum of adiabatic perturbations, generated in an exactly solvable inflationary model where the inflaton potential has a rapid change of its first derivative at some point, is confronted with existing observational data on angular fluctuations of the CMB temperature, galaxy clustering and peculiar velocities of galaxies. If we locate the step in the initial spectrum at k=0.05 h/Mpc, where some feature in the spectrum of Abell clusters of galaxies was found that could reflect a property of the initial spectrum, and if the large scales flat plateau of the spectrum is normalized according to the COBE data, the only remaining parameter of the spectrum is p - the ratio of amplitudes of the metric perturbations between the small scales and large scales flat plateaux. Allowed regions in the plane of parameters (Omega=1-Omega_Lambda, H_0) satisfying all data have been found for p lying in the region (0.8-1.7). Especially good agreement of the form of the present power spectrum in this model with the form of the cluster power spectrum is obtained for the inverted step (p<1, p=0.7-0.8), when the initial spectrum has slightly more power on small scales.

We calculate temperature anisotropies of the cosmic microwave background (CMB) for several initial power spectra of density perturbations with a built-in scale suggested by recent optical data on the spatial distribution of rich clusters of galaxies. Using cosmological models with different values of spectral index, baryon fraction, Hubble constant and cosmological constant, we compare the calculated radiation power spectrum with the CMB temperature anisotropies measured by the Saskatoon experiment. We show that spectra with a sharp peak at 120 h^{-1} Mpc are in agreement with the Saskatoon data. The combined evidence from cluster and CMB data favours the presence of a peak and a subsequent break in the initial matter power spectrum. Such feature is similar to the prediction of an inflationary model where an inflaton field is evolving through a kink in the potential.

Reheating after inflation occurs due to particle production by the oscillating inflaton field. In this paper we describe the perturbative approach to reheating, and then concentrate on effects beyond the perturbation theory. They are related to the stage of parametric resonance called preheating. It may occur in an expanding universe if the initial amplitude of oscillations of the inflaton field is large enough. We investigate a simple model of a massive inflaton field coupled to another scalar field X. Parametric resonance in this model is very broad. It occurs in a very unusual stochastic manner, which is different from the parametric resonance in the case when the expansion of the universe is neglected. Quantum fields interacting with the oscillating inflaton field experience a series of kicks which occur with phases uncorrelated to each other. We call this process stochastic resonance. We develop the theory of preheating taking into account the expansion of the universe and backreaction of produced particles, including the effects of rescattering. The process of preheating can be divided into several distinct stages. At the first stage the backreaction of created particles is not important. At the second stage backreaction increases the frequency of oscillations of the inflaton field, which makes the process even more efficient than before. Then the effects related to scattering of X-particles terminate the resonance. We calculate the density of X-particles and their quantum fluctuations with all backreaction effects taken into account. This allows us to find the range of masses and coupling constants for which one has efficient preheating. In particular, under certain conditions this process may produce particles with a mass much greater than the mass of the inflaton field.

We consider preheating in the theory $1/4 \lambda \phi^4 + 1/2 g^2\phi^2\chi^2 $, where the classical oscillating inflaton field $\phi$ decays into $\chi$-particles and $\phi$-particles. The parametric resonance which leads to particle production in this conformally invariant theory is described by the Lame equation. It significantly differs from the resonance in the theory with a quadratic potential. The structure of the resonance depends in a rather nontrivial way on the parameter $g^2/\lambda$. We construct the stability/instability chart in this theory for arbitrary $g^2/\lambda$. We give simple analytic solutions describing the resonance in the limiting cases $g^2/\lambda\ll 1$ and $g^2/\lambda \gg 1$, and in the theory with $g^2=3\lambda$, and with $g^2 =\lambda$. From the point of view of parametric resonance for $\chi$, the theories with $g^2=3\lambda$ and with $g^2 =\lambda$ have the same structure, respectively, as the theory $1/4 \lambda \phi^4$, and the theory $\lambda /(4 N) (\phi^2_i)^2$ of an N-component scalar field $\phi_i$ in the limit $N \to \infty$. We show that in some of the conformally invariant theories such as the simplest model $1/4 \lambda\phi^4$, the resonance can be terminated by the backreaction of produced particles long before $<\chi^2>$ or $<\phi^2 >$ become of the order $\phi^2$. We analyze the changes in the theory of reheating in this model which appear if the inflaton field has a small mass.

We present a new data analysis method to study rectangular T^3 ``small universes" with one or two of its dimensions significantly smaller than the present horizon (which we refer to as T^1- and T^2-models, respectively). We find that the 4 year COBE/DMR data set a lower limit on the smallest cell size for T^1- and T^2-models of 3000h^{-1} Mpc, at 95% confidence, for a scale invariant power spectrum (n=1).

Although the cubic T^3 "small universe" has been ruled out by COBE/DMR results as an interesting cosmological model, we still have the possibility of living in a universe with a more anisotropic topology such as a rectangular T^3 "small universe" with one or two of its dimensions significantly smaller than the present horizon (which we refer to as T^1- and T^2-models, respectively). In order to rule out these anisotropic topologies as well, we apply a new data analysis method that searches for the specific kind of symmetries that these models should produce. We find that the 4 year COBE/DMR data set a lower limit on the smallest cell size for T^1- and T^2-models of 3000h^{-1} Mpc, at 95% confidence, for a scale invariant power spectrum (n=1). These results imply that ALL toroidal universes (cubes and rectangles) are ruled out as interesting cosmological models.

We investigate properties of the correlation function of clusters of galaxies using geometrical models. On small scales the correlation function depends on the shape and the size of superclusters. On large scales it describes the geometry of the distribution of superclusters. If superclusters are distributed randomly then the correlation function on large scales is featureless. If superclusters and voids have a tendency to form a regular lattice then the correlation function on large scales has quasi-regularly spaced maxima and minima of decaying amplitude; i.e., it is oscillating. The period of oscillations is equal to the step size of the grid of the lattice. We calculate the power spectrum for our models and compare the geometrical information of the correlation function with other statistics. We find that geometric properties (the regularity of the distribution of clusters on large scales) are better quantified by the correlation function. We also analyse errors in the correlation function and the power spectrum by generating random realizations of models and finding the scatter of these realizations.

We use rich clusters of galaxies in the Northern and Southern Galactic hemispheres up to a redshift z=0.12 to determine the cluster correlation function. We show that superclusters of galaxies and voids between them form a moderately regular network. As a result the correlation function determined for clusters located in rich superclusters oscillates: it has a series of regularly spaced secondary maxima and minima. The scale of the supercluster-void network, determined from the period of oscillations, is P=115 Mpc. The amplitude of oscillations is about 3 times larger than the estimated error. We compare the observed cluster correlation function with similar functions derived for popular models of structure formation. We find that the production of the observed cluster correlation function in any model with a smooth transition of the power spectrum from a Harrison-Zeldovich regime with positive spectral index on long wavelengths to a negative spectral index on short wavelengths is highly unlikely. The power spectrum must have an extra peak located at the wavelength equal to the period of oscillations of the correlation function.

Using a new compilation of available data on galaxy clusters and superclusters we present evidence for a quasiregular three-dimensional network of rich superclusters and voids, with the regions of high density separated by about 120 Mpc. We calculate the power spectrum for clusters of galaxies; it has a peak on the wavelength equal to the step of the network; the excess in the amplitude of the spectrum over that of the cold dark matter model is by a factor of 1.4. The probability that the spectrum can be formed within the framework of the standard cosmogony is very small. If the cluster distribution reflects the distribution of all matter (luminous and dark), then there must exists some hithero unknown process that produces regular structure on large scales.

Transition from quantum to semiclassical behaviour and loss of quantum coherence for inhomogeneous perturbations generated from a non-vacuum initial state in the early Universe is considered in the Heisenberg and the Schr\"odinger representations, as well as using the Wigner function. We show explicitly that these three approaches lead to the same prediction in the limit of large squeezing (i.e. when the squeezing parameter $|r_k|\to \infty$): each two-modes quantum state (k, -k) of these perturbations is equivalent to a classical perturbation that has a stochastic amplitude, obeying a non-gaussian statistics which depends on the initial state, and that belongs to the quasi-isotropic mode (i.e. it possesses a fixed phase). The Wigner function is not everywhere positive for any finite $r_k$, hence its interpretation as a classical distribution function in phase space is impossible without some coarse graining procedure. However, this does not affect the transition to semiclassical behaviour since the Wigner function becomes concentrated near a classical trajectory in phase space when $|r_k|\to \infty$ even without coarse graining. Deviations of the statistics of the perturbations in real space from a Gaussian one lie below the cosmic variance level for the N-particles initial states with N=N(|k|) but may be observable for other initial states without statistical isotropy or with correlations between different k modes. As a way to look for this effect, it is proposed to measure the kurtosis of the angular fluctuations of the cosmic microwave background temperature.

New systematic classification of cosmological models of the present Universe is introduced. After making the comparison of these models with all existing observational data three viable models remain: the cold dark matter model with the cosmological constant (which becomes the most reasonable one if the Hubble constant $H_0>60$ km/s/Mpc); the mixed cold-hot dark matter model (models with two and especially three types of neutrinos with equal masses are in a slightly better agreement with observational data than the model with one massive neutrino); the pure cold dark matter model with a step-like initial spectrum of perturbations. The two latter models require $H_0\leq 60$ km/s/Mpc.

Spectrum of initial scalar and tensor perturbations created during an inflationary stage producing a closed or open FRW universe now is discussed. In the closed case, the CMB temperature anisotropy ${\Delta T/T}$ generated by scalar perturbations is enhanced for low multipoles. It is argued that in the open case there is no suppression of low multipoles. A possibility of the existence of a preferred space direction in the open case is noted.

Transition to the semiclassical behaviour and the decoherence process for inhomogeneous perturbations generated from the vacuum state during an inflationary stage in the early Universe are considered both in the Heisenberg and the Schr\"odinger representations to show explicitly that both approaches lead to the same prediction: the equivalence of these quantum perturbations to classical perturbations having stochastic Gaussian amplitudes and belonging to the quasi-isotropic mode. This equivalence and the decoherence are achieved once the exponentially small (in terms of the squeezing parameter $r_k$) decaying mode is neglected. In the quasi-classical limit $|r_k|\to \infty$, the perturbation mode functions can be made real by a time-independent phase rotation, this is shown to be equivalent to a fixed relation between squeezing angle and phase for all modes in the squeezed-state formalism. Though the present state of the gravitational wave background is not a squeezed quantum state in the rigid sense and the squeezing parameters loose their direct meaning due to interaction with the environment and other processes, the standard predictions for the rms values of the perturbations generated during inflation are not affected by these mechanisms (at least, for scales of interest in cosmological applications). This stochastic background still occupies a small part of phase space.

At the first stage of reheating after inflation, parametric resonance may rapidly transfer most of the energy of an inflaton field $\phi$ to the energy of other bosons. We show that quantum fluctuations of scalar and vector fields produced at this stage are much greater than they would be in a state of thermal equilibrium. This leads to cosmological phase transitions of a new type, which may result in a copious production of topological defects and in a secondary stage of inflation after reheating.

Skewness of temperature fluctuations of the cosmic microwave background (CMB) produced by initially Gaussian adiabatic perturbations with the flat (Harrison-Zeldovich) spectrum, which arises due to non-linear corrections to a gravitational potential at the matter-dominated stage, is calculated quantitatively. For the standard CDM model, the effect appears to be smaller than expected previously and lies below the cosmic variance limit even for small angles. The sign of the skewness is opposite to that of the skewness of density perturbations.

The gravitational waves (GW) background generated in a double inflationary model, with two scalar fields mutually interacting through gravity only, and its relative contribution $T/S$ to large-angle temperature fluctuations of the relic microwave background are investigated in detail. The relation between $T/S$ and the slope of the GW spectrum $n_T$ is shown to be a discriminative test between a slow-roll inflation driven by one scalar field and more complicated models. It is found that the GW amplitude is not exactly zero in minima of spectral oscillations, this property being an observational, in principle, manifestation of GW being in a squeezed vacuum state during inflation.

Mixed cold-hot dark matter cosmological models (CHDM) with $\Omega_{tot}=1$, approximately flat initial spectrum of adiabatic perturbations and 1, 2 or 3 types of massive neutrinos are compared and tested using recent observational data. The models with 2 or 3 neutrino types of equal mass permit as the best fit larger values of both the Hubble constant ($H_0\le 60$ for 2 types, $H_0\le 65$ for 3 types) and the total $\Omega_{\nu}$ (up to 0.3 for 3 types) than the model with 1 massive type. Also, they have less problems with abundances of early compact objects including $Ly-\alpha$ clouds.

Spectrum of density perturbations in the Universe generated from quantum-gravitational fluctuations in slow-roll-over inflationary scenarios with the Brans-Dicke gravity is calculated. It is shown that after inflation the isocurvature mode of perturbations may be neglected as compared to the adiabatic mode, and that an amplitude of the latter mode is not significantly different from that in the Einstein gravity. However, the account of the isocurvature mode is necessary to obtain the quantitatively correct spectrum of adiabatic perturbations.

The theory of reheating of the Universe after inflation is developed. The transition from inflation to the hot Universe turns out to be strongly model-dependent and typically consists of several stages. Immediately after inflation the field $\phi$ begins rapidly rolling towards the minimum of its effective potential. Contrary to some earlier expectations, particle production during this stage does not lead to the appearance of an extra friction term $\Gamma\dot\phi$ in the equation of motion of the field $\phi$. Reheating becomes efficient only at the next stage, when the field $\phi$ rapidly oscillates near the minimum of its effective potential. We have found that typically in the beginning of this stage the classical inflaton field $\phi$ very rapidly (explosively) decays into $\phi$-particles or into other bosons due to broad parametric resonance. This stage cannot be described by the standard elementary approach to reheating based on perturbation theory. The bosons produced at this stage, as well as some part of the classical field $\phi$ which survives the stage of explosive reheating, should further decay into other particles, which eventually become thermalized. The last stages of decay can be described in terms of perturbation theory. Complete reheating is possible only in those theories where a single massive $\phi$-particle can decay into other particles. This imposes strong constraints on the structure of inflationary models. On the other hand, this means that a scalar field can be a cold dark matter candidate even if it is strongly coupled to other fields.

The mixed cold-hot dark matter cosmological model (CHDM) with $\Omega_{tot}=1$ and a falling power-law initial spectrum of Gaussian adiabatic perturbations ($n>1$) is tested using recent obserbational data. It is shown that its fit to the data becomes worse with the growth of $n-1$, and may be considered as unreasonable for $n>1.1$ for all possible values of the Hubble constant. Thus, the CHDM model with a falling initial spectrum is worse than the same model with the approximately flat $(|n-1|<0.1)$ spectrum. On the other hand, the CHDM model provides a rather good fit to the data if $n$ lies in the range $(0.9-1.0)$, the Hubble constant $H_0 < 60$ km/s/Mpc ($H_0 < 55$ for $n=1$) and the neutrino energy density $\Omega_{\nu}< 0.25$. So, the CHDM model provides the best possibility for the realization of the simplest variants of the inflationary scenario having the effective slope $n\approx (0.95-0.97)$ between galaxy and horizon scales, including a modest contribution of primordial gravitational wave background to large-angle $\Delta T/T$ fluctuations of the cosmic microwave background (resulting in the increase of their total {\it rms} amplitude by $(5-10)%$) expected in some variants. A classification of cosmological models according to the number of fundamental parameters used to fit observational data is presented, too.

Behaviour of a weekly self-interacting scalar field with a small mass in the de Sitter background is investigated using the stochastic approach (including the case of a double-well interaction potential). Existence of the de Sitter invariant equilibrium quantum state of the scalar field in the presence of interaction is shown for any sign of the mass term. The stochastic approach is further developed to produce a method of calculation of an arbitrary anomalously large correlation function of the scalar field in the de Sitter background, and expressions for the two-point correlation function in the equilibrium state, correlation time and spatial physical correlation radius are presented. The latter does not depend on time that implies that the characteristic size of domains with positive and negative values of the scalar field remains the same on average in the equilibrium state in spite of the expansion of the $t=const$ hypersurface of the de Sitter space-time.

This is an old article which has never been posted for public use and which only appeard in a relatively hard-to-get Proceedings of the Sakharov Conference in Moscow (May, 1991). The subject of this article has received a lot of attention in last two years. Although we acknowledge that there has been a very big progress in the field of relating the spectrum of fluctuations to scalar field potentials of inflation models, which have outdated some of our results (particularly because we never took into account the tensor fluctuations which have been shown to be crucial for reconstructing the inflaton potential) we think that some of the other ideas and approaches which we presented in this article maybe still of interest. We did not rewrite the article to include any of the newer approaches and we did not add any of the newer references. The reader is advised to refer to recent papers (for example, see the paper astro-ph/9403001 by Edward Kolb and Sharon Vadas) for state of the art review and up-to-date references. ==The original abstract of the article==: We prove that, in stochastic approach, there exists an equivalence relation between different inflationary models under some redefinition of field and time variables. The postinflationary physics is insensitive to it and one can say that related theories are in fact indistinguishable from the viewpoint

Dynamics of long-wave isocurvature perturbations during an inflationary stage in multiple (multi-component) inflationary models is calculated analytically for the case where scalar fields producing this stage interact between themselves through gravity only. This enables to determine correct amplitudes of such perturbations produced by vacuum quantum fluctuations of the scalar fields during the multiple inflationary stage. Exact matching to a post-inflationary evolution that gives the amplitude of isocurvature perturbations in the cold dark matter model with radiation is performed in the case where a massive inflaton field remains uncoupled from usual matter up to the present time. For this model, isocurvature perturbations are smaller than adiabatic ones in the region of the break in the perturbation spectrum which arises due to a transition between the two phases of inflation, but they may be much bigger and have a maximum at much shorter scales. The case of an inflaton with a quartic coupling which remains uncoupled after inflation is considered, too.

We consider double-inflationary models with two noninteracting scalar fields, a light scalar field $\phi_l$ with potential $\frac{1}{2}m_l^2\phi_l^2$ and a heavy scalar field $\phi_h$ with potential $\frac{\lambda}{n}\phi_h^n$ with $n=2,4$. CDM with the initial spectrum of adiabatic perturbations produced in these models is compared with observations. These models contain two more free parameters than the standard CDM model with an initial scale-invariant spectrum. We normalize our spectra to COBE DMR and compare the predictions with observations on the biasing factor, large-scale peculiar velocities, quasar and galaxy formation and the Stromlo-APM counts-in-cells analysis. The model with $n=4$ is excluded by the data while for the $n=2$ model, taking cosmic variance into account, a small window of parameters compatible with observations is found.

We compare different nonlinear approximations to gravitational clustering in the weakly nonlinear regime, using as a comparative statistic the evolution of non-Gaussianity which can be characterised by a set of numbers $S_p$ describing connected moments of the density field at the lowest order in $<\delta^2>$: $<\delta^n>_c \simeq S_n<\delta^2>^{n-1}$. Generalizing earlier work by Bernardeau (1992) we develop an ansatz to evaluate all $S_p$ in a given approximation by means of a generating function which can be shown to satisfy the equations of motion of a homogeneous spherical density enhancement in that approximation. On the basis of the values of we show that approximations formulated in Lagrangian space (such as the Zeldovich approximation and its extensions) are considerably more accurate than those formulated in Eulerian space such as the Frozen Flow and Linear Potential approximations. In particular we find that the $n$th order Lagrangian perturbation approximation correctly reproduces the first $n+1$ parameters $S_n$. We also evaluate the density probability distribution function for the different approximations in the quasi-linear regime and compare our results with an exact analytic treatment in the case of the Zeldovich approximation.

Nonlinear approximation methods such as the Zeldovich approximation, and more recently the frozen flow and linear potential approximations, are sometimes used to simulate nonlinear gravitational instability in the expanding Universe. We investigate the relative accuracy of these approximations by comparing them with the exact solution using second order perturbation theory. We evaluate the density and velocity fields in these approximations to second order, and also determine the skewness parameter $S_3 = \langle \delta^3 \rangle / \langle (\delta^{(1)})^2\rangle ^2$ for each of the approximations again in second order. We find that $S_3 = 4, ~3, ~3.4$ for the Zeldovich approximation, the frozen flow and the linear potential approximations respectively as compared to $S_3 = {34/7}$ for the exact solution. Our results show that, of all the approximations considered, the Zeldovich approximation is the most accurate in describing the weakly nonlinear effects of gravity. Moreover, the Zeldovich approximation is much closer to the exact results for matter and velocity distributions than the other approximations if the slope of the power spectrum of density perturbations is $-3< n \le -1$.

CDM models with non-scale-free step-like spectra of adiabatic perturbations produced in a realistic double inflationary model are compared with recent observational data. The model contains two additional free parameters relatively to the standard CDM model with the flat ($n=1$) initial spectrum. Results of the COBE experiment are used for the determination of a free overall spectrum normalization. Then predictions for the galaxy biasing parameter, the variance for "counts in cells", the galaxy angular correlation function, bulk flow peculiar velocities and the Mach number test are obtained. Also considered are conditions for galaxy and quasar formation. Observational data strongly restricts allowed values for the two remaining model parameters. However, a non-empty region for them satisfying all considered tests is found.

If the Universe has the topology of a 3-torus ($T^3$), then it follows from recent data on large-scale temperature fluctuations $\Delta T/T$ of the cosmic microwave background that either the minimal size of the torus is at least of the order of the present cosmological horizon, or the large-scale $\Delta T/T$ pattern should have a symmetry plane (in case of the effective $T^1$ topology) or a symmetry axis (in case of the effective $T^2$ topology), the latter possibility being probably excluded by the data.