We discuss the gravitational creation of superheavy particles $\chi$ in an inflationary scenario with a quartic potential and a non-minimal coupling between the inflaton $\varphi$ and the Ricci curvature: $\xi \varphi^2 R/2$. We show that for large constants $\xi >> 1$, there can be abundant production of particles $\chi$ with masses largely exceeding the inflationary Hubble rate $H_{infl}$, up to $(a~few) \times \xi H_{infl}$, even if they are conformally coupled to gravity. We discuss two scenarios involving these gravitationally produced particles $\chi$. In the first scenario, the inflaton has only gravitational interactions with the matter sector and the particles $\chi$ reheat the Universe. In this picture, the inflaton decays only due to the cosmic expansion, and effectively contributes to dark radiation, which can be of the observable size. The existing limits on dark radiation lead to an upper bound on the reheating temperature. In the second scenario, the particles $\chi$ constitute Dark Matter, if substantially stable. In this case, their typical masses should be in the ballpark of the Grand Unification scale.

We consider Dark Matter composed of an oscillating singlet scalar field. On top of the mass term, the scalar is equipped with a potential spontaneously breaking Z_2-symmetry. This potential dominates at early times and leads to the time-dependent expectation value of the scalar, which decreases in the expanding Universe. As it drops below some critical value, the symmetry gets restored, and the Dark Matter field starts to oscillate around zero. We arrange the spontaneous symmetry breaking through the interaction of the scalar with the Ricci curvature. In that way, superheavy Dark Matter can be produced at very early times. Depending on its mass, the production takes place at inflation (very large masses up to the Grand Unification scale), at preheating, or at radiation-dominated stage (masses 10^{6}-10^{7} Gev).

We discuss the Damour--Esposito-Far\`ese model of gravity, which predicts the spontaneous scalarization of neutron stars in a certain range of parameter space. In the cosmological setup, the scalar field responsible for scalarization is subject to a tachyonic instability during inflation and the matter domination stage, resulting in a large value of the field today. This value feeds into the PPN parameters, which turn out to be in gross conflict with the Solar system measurements. We modify the original Damour--Esposito-Far\`ese model by coupling the scalar to the inflaton field. This coupling acts as an effective mass for the scalar during inflation. For generic couplings that are not extremely small, the scalar (including its perturbations) relaxes to zero with an exponential accuracy by the beginning of the hot stage. While the scalar exhibits growth during the subsequent cosmological stages, the resulting present value remains very small---in a comfortable agreement with the Solar system tests.

In a subclass of scalar-tensor theories, it has been shown that standard general relativity solutions of neutron stars and black holes with trivial scalar field profiles are unstable. Such an instability leads to solutions which are different from those of general relativity and have non-trivial scalar field profiles, in a process called scalarization. In the present work we focus on scalarization due to a non-minimal coupling of the scalar field to the Gauss-Bonnet curvature invariant. The coupling acts as a tachyonic mass for the scalar mode, thus leading to the instability of general relativity solutions. We point out that a similar effect may occur for the scalar modes in a cosmological background, resulting in the instability of cosmological solutions. In particular, we show that a catastrophic instability develops during inflation within a period of time much shorter than the minimum required duration of inflation. As a result, the standard cosmological dynamics is not recovered. This raises the question of the viability of scalar-Gauss-Bonnet theories exhibiting scalarization.

We study in detail the recently proposed mechanism of generating superheavy Dark Matter with the mass larger than the Hubble rate at the end of inflation. A real scalar field constituting Dark Matter linearly couples to the inflaton. As a result of this interaction, the scalar gets displaced from its zero expectation value. This offset feeds into the energy density of Dark Matter. This mechanism is universal and can be implemented in a generic inflationary scenario. Phenomenology of the model is comprised of Dark Matter decay into inflatons, which in turn decay into Standard Model species triggering cascades of high energy particles contributing to the cosmic ray flux. We evaluate the lifetime of Dark Matter and obtain limits on the inflationary scenarios, where this mechanism does not lead to the conflict with the Dark Matter stability considerations/studies of cosmic ray propagation.

We introduce a class of non-supersymmetric models explaining baryogenesis a la Affleck-Dine, which use a decay of two superheavy scalar fields with close masses. These scalars acquire non-zero expectation values during inflation through linear couplings to a function of an inflaton. After the inflaton decay, the model possesses approximate U(1)-invariance, explicitly broken by a small mass splitting. This splitting leads to the baryogenesis in the early Universe. Resulting baryon asymmetry is automatically small for the scalars with the masses about the Grand Unification scale and larger. It is fully determined by the inflaton dynamics and the Lagrangian parameters, i.e., is independent of initial pre-inflationary conditions for the scalars. As a consequence, baryon perturbations are purely adiabatic. We point out a possible origin of the mass splitting: masses of scalars degenerate at some large energy scale may acquire different loop corrections due to the interaction with the inflaton. Compared to electroweak baryogenesis and conventional Affleck-Dine scenarios, our mechanism generically leads to the proton decay suppressed by the powers of the superheavy scalar masses, which makes this scenario potentially testable.

We study the correspondence between models of a self-interacting canonical complex scalar field and P(X)-theories/shift-symmetric k-essence. Both describe the same background cosmological dynamics, provided that the amplitude of the complex scalar is frozen modulo the Hubble drag. We compare perturbations in these two theories on top of a fixed cosmological background. The dispersion relation for the complex scalar has two branches. In the small momentum limit, one of these branches coincides with the dispersion relation of the P(X)-theory. Hence, the low momentum phase velocity agrees with the sound speed in the corresponding P(X)-theory. The behavior of high frequency modes associated with the second branch of the dispersion relation depends on the value of the sound speed. In the subluminal case, the second branch has a mass gap. On the contrary, in the superluminal case, this branch is vulnerable to a tachyonic instability. We also discuss the special case of the P(X)-theories with an imaginary sound speed leading to the catastrophic gradient instability. The complex field models provide with a cutoff on the momenta involved in the instability.

We propose a universal mechanism of producing dark matter and baryon (lepton) charge at the stage of the quasi-de Sitter expansion of the Universe---inflation. The key ingredient of the mechanism is a linear coupling of the field, responsible for generation of dark matter or baryon (lepton) charge, to a function of the inflaton. During inflation this induces almost constant force dragging the corresponding field to the non-zero value. This force explicitly breaks quantum numbers associated with dark matter/baryon abundance at later stages. As a particular realization of the mechanism we introduce a super-heavy complex scalar field with the mass larger than the Hubble rate during the last e-folds of inflation. The global U(1)-symmetry is violated due to the linear coupling of the phase of the complex scalar to the inflaton. The symmetry breaking leads to the generation of a non-zero Noether charge. The latter is directly related to the dark matter abundance, or, alternatively, can be converted into baryon asymmetry, if the complex scalar carries the baryon charge.

We study a particular realization of the cosmological bounce scenario proposed recently by Ijjas and Steinhardt. First, we find that their bouncing solution starts from a divergent sound speed and ends with its vanishing. Thus, the solution connects two strongly coupled configurations. These pathologies are separated from the bouncing regime by only a few Planck times. We then reveal the exact structure of the Lagrangian, which reproduces this bouncing solution. This reconstruction allowed us to consider other cosmological solutions of the theory and analyze the phase space. In particular, we find other bouncing solutions and solutions with superluminal sound speed. These stable superluminal states can be continuously transformed into the solution constructed by Ijjas and Steinhardt. We discuss the consequences of this feature for a possible UV-completion.

Both k-essence and the pressureless perfect fluid develop caustic singularities at finite time. We further explore the connection between the two and show that they belong to the same class of models, which admits the caustic free completion by means of the canonical complex scalar field. Specifically, the free massive/self-interacting complex scalar reproduces dynamics of pressureless perfect fluid/shift-symmetric k-essence under certain initial conditions in the limit of large mass/sharp self-interacting potential. We elucidate a mechanism of resolving caustic singularities in the complete picture. The collapse time is promoted to complex number. Hence, the singularity is not developed in real time. The same conclusion holds for a collection of collisionless particles modelled by means of the Schroedinger equation, or ultra-light axions (generically, coherent oscillations of bosons in the Bose--Einstein condensate state).

Several early Universe scenarios predict a direction-dependent spectrum of primordial curvature perturbations. This translates into the violation of the statistical isotropy of cosmic microwave background radiation. Previous searches for statistical anisotropy mainly focussed on a quadrupolar direction-dependence characterised by a single multipole vector and an overall amplitude $g_*$. Generically, however, the quadrupole has a more complicated geometry described by two multipole vectors and $g_*$. This is the subject of the present work. In particular, we limit the amplitude $g_*$ for different shapes of the quadrupole by making use of Planck 2015 maps. We also constrain certain inflationary scenarios which predict this kind of more general quadrupolar statistical anisotropy.

Motivated by the projectable Horava--Lifshitz model/mimetic matter scenario, we consider a particular modification of standard gravity, which manifests as an imperfect low pressure fluid. While practically indistinguishable from a collection of non-relativistic weakly interacting particles on cosmological scales, it leaves drastically different signatures in the Solar system. The main effect stems from gravitational focusing of the flow of Imperfect Dark Matter passing near the Sun. This entails strong amplification of Imperfect Dark Matter energy density compared to its average value in the surrounding halo. The enhancement is many orders of magnitude larger than in the case of Cold Dark Matter, provoking deviations of the metric in the second order in the Newtonian potential. Effects of gravitational focusing are prominent enough to substantially affect the planetary dynamics. Using the existing bound on the PPN parameter $\beta_{PPN}$, we deduce a stringent constraint on the unique constant of the model.

We consider the branch of the projectable Horava-Lifshitz model which exhibits ghost instabilities in the low energy limit. It turns out that, due to the Lorentz violating structure of the model and to the presence of a finite strong coupling scale, the vacuum decay rate into photons is tiny in a wide range of phenomenologically acceptable parameters. The strong coupling scale, understood as a cutoff on ghosts' spatial momenta, can be raised up to $\Lambda \sim 10$ TeV. At lower momenta, the projectable Horava-Lifshitz gravity is equivalent to General Relativity supplemented by a fluid with a small positive sound speed squared ($10^{-42}\lesssim$) $c^2_s \lesssim 10^{-20}$, that could be a promising candidate for the Dark Matter. Despite these advantages, the unavoidable presence of the strong coupling obscures the implementation of the original Horava's proposal on quantum gravity. Apart from the Horava-Lifshitz model, conclusions of the present work hold also for the mimetic matter scenario, where the analogue of the projectability condition is achieved by a non-invertible conformal transformation of the metric.

We discuss initial conditions for the recently proposed Imperfect Dark Matter (Modified Dust). We show that they are adiabatic under fairly moderate assumptions about the cosmological evolution of the Universe at the relevant times.

At large scales and for sufficiently early times, dark matter is described as a pressureless perfect fluid---dust---non-interacting with Standard Model fields. These features are captured by a simple model with two scalars: a Lagrange multiplier and another playing the role of the velocity potential. That model arises naturally in some gravitational frameworks, e.g., the mimetic dark matter scenario. We consider an extension of the model by means of higher derivative terms, such that the dust solutions are preserved at the background level, but there is a non-zero sound speed at the linear level. We associate this {\it Modified Dust} with dark matter, and study the linear evolution of cosmological perturbations in that picture. The most prominent effect is the suppression of their power spectrum for sufficiently large cosmological momenta. This can be relevant in view of the problems that cold dark matter faces at sub-galactic scales, e.g., the missing satellites problem. At even shorter scales, however, perturbations of Modified Dust are enhanced compared to the predictions of more common particle dark matter scenarios. This is a peculiarity of their evolution in radiation dominated background. We also briefly discuss clustering of Modified Dust. We write the system of equations in the Newtonian limit, and sketch the possible mechanism which could prevent the appearance of caustic singularities. The same mechanism may be relevant in light of the core-cusp problem.

We revisit constraints on the (pseudo)conformal Universe from the non-observation of statistical anisotropy in the Planck data. The quadratic maximal likelihood estimator is applied to the Planck temperature maps at frequencies 143 GHz and 217 GHz as well as their cross-correlation. The strongest constraint is obtained in the scenario of the (pseudo)conformal Universe with a long intermediate evolution after conformal symmetry breaking. In terms of the relevant parameter (coupling constant), the limit is h^2 <0.0013 at 95% C.L. (using the cross-estimator). The analogous limit is much weaker in the scenario without the intermediate stage (h^2 \ln \frac{H_0}{\Lambda}<0.52) allowing the coupling constant to be of order one. In the latter case, the non-Gaussianity in the 4-point function appears to be a more promising signature.

We consider a non-inflationary early Universe scenario in which relevant scalar perturbations get frozen out at some point, but then are defrosted and follow a long nearly Minkowskian evolution before the hot era. This intermediate stage leaves specific imprint on the CMB 3-point function, largely independent of details of microscopic physics. In particular, the CMB bispectrum undergoes oscillations in the multipole l space with roughly constant amplitude. The latter is in contrast to the oscillatory bispectrum enhanced in the flattened triangle limit, as predicted by inflation with non-Bunch--Davies vacuum. Given this and other peculiar features of the bispectrum, stringent constraints imposed by the Planck data may not apply. The CMB 3-point function is suppressed by the inverse duration squared of the Minkowskian evolution, but can be of observable size for relatively short intermediate Minkowskian stage.

We constrain several models of the early Universe that predict a statistical anisotropy of the cosmic microwave background (CMB) sky. We make use of WMAP9 maps deconvolved with beam asymmetries. As compared to previous releases of WMAP data, they do not exhibit the anomalously large quadrupole of statistical anisotropy. This allows to strengthen the limits on the parameters of models established earlier in the literature. In particular, the amplitude of the special quadrupole is constrained as |g_*|<0.072 at 95% C.L. (-0.046<g_*<0.048 at 68% C.L.) independently of the preferred direction in the sky. The upper limit is obtained on the total number of e-folds in anisotropic inflation with the Maxwellian term nonminimally coupled to the inflaton, namely N_{tot} <N_{CMB} +82 at 95% C.L. (+14 at 68% C.L.) for N_{CMB}=60. We also constrain models of the (pseudo)conformal universe. The strongest constraint is obtained for spectator scenarios involving a long stage of subhorizon evolution after conformal rolling, which reads h^2 < 0.006 at 95% C.L., in terms of the relevant parameter. The analogous constraint is much weaker in dynamical models, e.g., Galilean genesis.

The low multipole anomalies of the Cosmic Microwave Background has received much attention during the last few years. It is still not ascertained whether these anomalies are indeed primordial or the result of systematics or foregrounds. An example of a foreground, which could generate some non-Gaussian and statistically anisotropic features at low multipole range, is the very symmetric Kuiper Belt in the outer solar system. In this paper, expanding upon the methods presented by Maris et al. (2011), we investigate the contributions from the Kuiper Belt objects (KBO) to the WMAP ILC 7 map, whereby we can minimize the contrast in power between even and odd multipoles in the CMB, discussed discussed by Kim & Naselsky (2010). We submit our KBO de-correlated CMB signal to several tests, to analyze its validity, and find that incorporation of the KBO emission can decrease the quadrupole-octupole alignment and parity asymmetry problems, provided that the KBO signals has a non-cosmological dipole modulation, associated with the statistical anisotropy of the ILC 7 map. Additionally, we show that the amplitude of the dipole modulation, within a 2 sigma interval, is in agreement with the corresponding amplitudes, discussed by Lew (2008).

Search for the statistical anisotropy in the CMB data is a powerful tool for constraining models of the early Universe. In this paper we focus on the recently proposed cosmological scenario with conformal rolling. We consider two sub-scenarios, one of which involves a long intermediate stage between conformal rolling and conventional hot epoch. Primordial scalar perturbations generated within these sub-scenarios have different direction-dependent power spectra, both characterized by a single parameter h^2. We search for the signatures of this anisotropy in the seven-year WMAP data using quadratic maximum likelihood method, first applied for similar purposes by Hanson and Lewis. We confirm the large quadrupole anisotropy detected in V and W bands, which has been argued to originate from systematic effects rather than from cosmology. We construct an estimator for the parameter h^2. In the case of the sub-scenario with the intermediate stage we set an upper limit h^2 < 0.045 at the 95% confidence level. The constraint on h^2 is much weaker in the case of another sub-scenario, where the intermediate stage is absent.

Scalar cosmological perturbations with nearly flat power spectrum may originate from perturbations of the phase of a scalar field conformally coupled to gravity and rolling down negative quartic potential. We consider a version of this scenario whose specific property is a long intermediate stage between the end of conformal rolling and horizon exit of the phase perturbations. Such a stage is natural, e.g., in cosmologies with ekpyrosis or genesis. Its existence results in small negative scalar tilt, statistical anisotropy of all even multipoles starting from quardupole of general structure (in contrast to the usually discussed single quadrupole of special type) and non-Gaussianity of a peculiar form.