The modified supergravity approach is applied to describe a formation of Primordial Black Holes (PBHs) after Starobinsky inflation. Our approach naturally leads to the two-(scalar)-field attractor-type double inflation, whose first stage is driven by scalaron and whose second stage is driven by another scalar field which belongs to a supergravity multiplet. The scalar potential and the kinetic terms are derived, the vacua are studied, and the inflationary dynamics of those two scalars is investigated. We numerically compute the power spectra and we find the ultra-slow-roll regime leading to an enhancement (peak) in the scalar power spectrum. This leads to an efficient formation of PBHs. We estimate the masses of PBHs and we find their density fraction (as part of Dark Matter). We show that our modified supergravity models are in agreement with inflationary observables, while they predict the PBH masses in a range between $$10^{16}$$10 16 g and $$10^{20}$$10 20 g. In this sense, modified supergravity provides a natural top-down approach for explaining and unifying the origin of inflation and the PBHs Dark Matter.

We investigate the production of primordial black holes and their contribution to the presently observed dark matter in a dilaton two-field extension of Starobinsky’s quadratic $f(R)$ model of inflation. The model features a multi-field amplification mechanism which leads to the generation of a sharp peak in the inflationary power spectrum at small wavelengths responsible for the production of primordial black holes. This mechanism is significantly different from single-field models and requires a stochastic treatment during an intermediate phase of the inflationary dynamics. We find that the model leads to a successful phase of effective single-field Starobinsky inflation for wavelengths probed by the cosmic microwave background radiation and explains the observed cold dark matter content in the Universe by the formation of primordial black holes.

We propose a large class of supergravity models in terms of a single chiral matter superfield, leading to (almost) arbitrary single-field inflationary scalar potentials similar to the F-term in rigid supersymmetry. Those scalar potentials are positively semi-definite (in our approximation), and can preserve supersymmetry at the end of inflation. The only scalar superpartner of inflaton is stabilized by a higher-dimensional term in the Kähler potential. We argue that couplings of the inflaton to other sectors of the particle spectrum do not affect the inflationary dynamics, and briefly discuss reheating of the universe by the inflaton decays.