The Starobinsky model of modified gravity suggested to explain dark energy may be also considered in the astrophysical context. Recently it has been pointed out that in contracting regions curvature oscillations around the GR value may lead to the production of high energy particles which contribute to the cosmic ray flux. We revisit these calculations in the Einstein frame and show that the continuous approximation for the matter density used in the original calculations is not valid. We show that this problem is generic in $F(R)$-gravity models introduced to describe the dark energy. We go beyond the approximation and find the rate of particle production to be negligible.

String scenarios generically predict that we live in a so called axiverse: the Universe with about a hundred of light axion species which are decoupled from the Standard Model particles. However, the axions can couple to the inflaton which leads to their production after inflation. Then, these axions remain in the expanding Universe contributing to the dark radiation component, which is severely bounded from present cosmological data. We place a general constraint on the axion production rate and apply it to several variants of reasonable inflaton-to-axion couplings. The limit merely constrains the number of ultralight axions and the relative strength of inflaton-to-axion coupling. It is valid in both large and small field inflationary models irrespectively of the axion energy scales and masses. Thus, the limit is complementary to those associated with the Universe overclosure and axion isocurvature fluctuations. In particular, a hundred of axions is forbidden if inflaton universally couples to all the fields at reheating. In the case of gravitational sector being responsible for the reheating of the Universe (which is a natural option in all inflationary models with modified gravity), the axion production can be efficient. We find that in the Starobinsky $R^2$-inflation even a single axion (e.g. the standard QCD-axion) is in tension with the Planck data, making the model inconsistent with the axiverse. The general conclusion is that an inflation with inefficient reheating mechanism and low reheating temperature may be in tension with the presence of light scalars.

We study the creation and evolution of cosmological perturbations in renormalizable quadratic gravity with a Weyl term. We adopt a prescription that implies the stability of the vacuum at the price of introducing a massive spin-two ghost state, leading to the loss of unitarity. The theory may still be predictive regardless the interpretation of non-unitary processes provided that their rate is negligible compared to the Universe expansion rate. This implies that the ghost is effectively stable. In such a setup, there are two scalar degrees of freedom excited during inflation. The first one is the usual curvature perturbation whose power spectrum appears to coincide with that of single-field inflation. The second one is a scalar component of the ghost encoded in the shift vector of the metric in the uniform inflaton gauge. The amplitudes of primordial tensor and vector perturbations are strongly suppressed. After inflation the ghost field starts to oscillate and its energy density shortly becomes dominant in the Universe. For all ghost masses allowed by laboratory constraints ghosts should have "overclosed" the Universe at temperatures higher than that of primordial nucleosynthesis. Thus, the model with the light Weyl ghost is ruled out.

We study the creation of high energy SM particles in the Starobinsky model of dark energy (a variant of F(R)-gravity) inside the regions contracting due to the Jeans instability. In this modification of gravity the additional degree of freedom -- scalaron -- behaves as a particle with mass depending on matter density. So when the mass changes light scalarons could be created at a non-adiabatic stage. Later scalaron mass grows and could reach large values, even the value of 10^{13} GeV, favored by early-time inflation. Heavy scalarons decay contributing to the cosmic ray flux. We analytically calculated the number density of created particles for the case of exponential (Jeans) contraction and found it negligibly small provided the phenomenologically viable and cosmologically interesting range of model parameters. We expect similar results for a generic model of F(R)-gravity mimicking cosmological constant.

Recent cosmological data favour R^2-inflation and some amount of non-standard dark radiation in the Universe. We show that a framework of high energy scale invariance can explain these data. The spontaneous breaking of this symmetry provides gravity with the Planck mass and particle physics with the electroweak scale. We found that the corresponding massless Nambu--Goldstone bosons -- dilatons -- are produced at reheating by the inflaton decay right at the amount needed to explain primordial abundances of light chemical elements and anisotropy of the cosmic microwave background. Then we extended the discussion on the interplay with Higgs-inflation and on general class of inflationary models where dilatons are allowed and may form the dark radiation. As a result we put a lower limit on the reheating temperature in a general scale invariant model of inflation.

We introduce conformal coupling of the Standard Model Higgs field to gravity and discuss the subsequent modification of $R^2$-inflation. The main observation is a lower temperature of reheating which happens mostly through scalaron decays into gluons due to the conformal (trace) anomaly. This modifies all predictions of the original $R^2$-inflation. To the next-to-leading order in slow roll parameters we calculate amplitudes and indices of scalar and tensor perturbations produced at inflation. The results are compared to the next-to-leading order predictions of $R^2$-inflation with minimally coupled Higgs field and of Higgs-inflation. We discuss additional features in gravity wave signal that may help to distinguish the proposed variant of $R^2$-inflation. Remarkably, the features are expected in the region available for study at future experiments like BBO and DECIGO. Finally, we check that (meta)stability of electroweak vacuum in the cosmological model is consistent with recent results of searches for the Higgs boson at LHC.