Cosmological evolution and particle creation in $R^2$-modified gravity are considered for the case of the dominant decay of the scalaron into a pair of gauge bosons due to conformal anomaly. It is shown that in the process of thermalization superheavy dark matter with the coupling strength typical for the GUT SUSY can be created. Such dark matter would have the proper cosmological density if the particle mass is close to $10^{12}$ GeV.

The review is devoted to consideration of possible observational consequences of modified gravity theories, suggested for explanation of the contemporary accelerated expansion of the universe. The major attention is paid to F(R)-models. It is shown that in systems with rising energy density high frequency and large amplitude oscillations of the curvature scalar, R(t), are induced. These oscillations lead to the production of elementary particles, which may be observed in the spectra of energetic cosmic rays. In the background of such oscillating solutions gravitation repulsion between finite-size objects becomes possible. Since the Lagrangian is a non-linear function of curvature, equations of motion become of higher (4th) order and exhibit very rich pattern of new physical effects. In particular, the evolution of density perturbations is strongly different from that in General Relativity, amplified due to both parametric resonance and anti-friction phenomena.

Production of massive stable relics in $R^2$-modified gravity is considered. It is shown that the cosmological evolution and kinetics of massive species differs significantly from those in the conventional cosmology. The results are applied to the lightest supersymmetric particles and it is argued that they are viable candidates for the constituents of dark matter, if their mass is about 1000 TeV.

Universe history in $R^2$-gravity is studied from "beginning" up to the present epoch. It is assumed that initially the curvature scalar $R$ was sufficiently large to induce the proper duration of inflation. Gravitational particle production by the oscillating $R(t)$ led to a graceful exit from inflation, but the cosmological evolution in the early universe was drastically different from the standard one till the universe age reached the value of the order of the inverse decay rate of the oscillating curvature $R(t)$. This deviation from the standard cosmology might have a noticeable impact on the formation of primordial black holes and baryogenesis. At later time, after exponential decay of the curvature oscillations, cosmology may return to normality.

Some problems of spontaneous and gravitational baryogenesis are discussed. Gravity modification due to the curvature dependent term in gravitational baryogensis scenario is considered. It is shown that the interaction of baryonic fields with the curvature scalar leads to strong instability of the gravitational equations of motion and as a result to noticeable distortion of the standard cosmology.

General features of spontaneous baryogenesis are studied. The relation between the time derivative of the (pseudo)goldstone field and the baryonic chemical potential is revisited. It is shown that this relation essentially depends upon the representation chosen for the fermionic fields with non-zero baryonic number (quarks). The calculations of the cosmological baryon asymmetry are based on the kinetic equation generalized to the case of non-stationary background. The effects of the finite interval of the integration over time are also included into consideration. All these effects combined lead to a noticeable deviation of the magnitude of the baryon asymmetry from the canonical results.

Generation of the cosmological baryon asymmetry in frameworks of spontaneous baryogenesis is studied in detail. It is shown that the relation between baryonic chemical potential and the time derivative of the (pseudo)Goldstone field essentially depends upon the representation chosen for the fermionic fields with non-zero baryonic number (quarks). Kinetic equation is modified and numerically solved in equilibrium for the case of time dependent external background or finite integration time to be applicable to the case when energy conservation law is formally violated.

Evolution of density and metric perturbations in the background of high frequency oscillations of curvature in F(R) gravity is considered. In addition to the usual Jeans-like instability new effects of amplification of perturbations, associated with parametric resonance and antifriction phenomena, are found.

Gravitational instability in classical Jeans theory, General Relativity, and modified gravity is considered. The background density increase leads to a faster growth of perturbations in comparison with the standard theory. The transition to the Newtonian gauge in the case of coordinate dependent background metric functions is studied. For modified gravity a new high frequency stable solution is found.

We study particle production in infrared-modified gravitational theories in the contemporary universe. It is shown that in astronomical systems with rising mass density, the curvature scalar may oscillate with very high frequency. These oscillations lead to efficient particle production and in an interesting range of model parameters could be a source of energetic cosmic rays. This effect either excludes some models of modified gravity or suggests a new mechanism of cosmic ray production.

Spherically symmetric solutions in F(R) theories in astronomical systems with rising energy density are studied. The range of parameters is established for which the flat space-time approximation for the background metric is valid. For the solutions in which the curvature scalar oscillates with large amplitude and high frequency, found in our previous papers, it is shown that the analysis of the Jeans instability is strongly modified. It is discovered that for large astronomical objects modified gravity can be repulsive, so such objects shrink forming relatively thin shells instead of quasi uniform bodies. This may explain the formation of cosmic voids.

It is shown that F(R)-modified gravitational theories lead to curvature oscillations in astrophysical systems with rising energy density. The frequency and the amplitude of such oscillations could be very high and would lead to noticeable production of energetic cosmic ray particles.

Some models of modified gravity and their observational manifestations are considered. It is shown, that gravitating systems with mass density rising with time evolve to a singular state with infinite curvature scalar. The universe evolution during the radiation dominated epoch is studied in $R^2$-extended gravity theory. Particle production rate by the oscillating curvature is calculated. Possible implications of the model for cosmological creation of non-thermal dark matter are discussed.

The Universe evolution during the radiation-dominated epoch in the R^2-extended gravity theory is considered. The equations of motion for R and H are solved analytically and numerically. The particle production rate by the oscillating curvature is calculated in one-loop approximation and the back reaction of particle production on the evolution of R is taken into account. Possible implications of the model for cosmological creation of non-thermal dark matter is discussed.

Observational manifestations of some models of modified gravity, which have been suggested to explain the accelerated cosmological expansion, are analyzed for gravitating systems with time dependent mass density. It is shown that if the mass density rises with time, the system evolves to the singular state with infinite curvature scalar. The corresponding characteristic time is typically much shorter than the cosmological time.