We address the question of regular primordial black holes with de Sitter interior, their remnants and gravitational vacuum solitons G-lumps as heavy dark matter candidates providing signatures for inhomogeneity of early universe, which is severely constrained by the condition that the contribution of these objects in the modern density does not exceed the total density of dark matter. Primordial black holes and their remnants seem to be most elusive among dark matter candidates. However, we reveal a nontrivial property of compact objects with de Sitter interior to induce proton decay or decay of neutrons in neutron stars. The point is that they can form graviatoms, binding electrically charged particles. Their observational signatures as dark matter candidates provide also signatures for inhomogeneity of the early universe. In graviatoms, the cross-section of the induced proton decay is strongly enhanced, what provides the possibility of their experimental searches. We predict proton decay paths induced by graviatoms in the matter as an observational signature for heavy dark matter searches at the IceCUBE experiment.

Source-free equations of nonlinear electrodynamics minimally coupled to gravity (NED-GR) admit regular axially symmetric asymptotically Kerr-Newman solutions, which describe electrically charged rotating black holes and spinning solitons. Asymptotic analysis of solutions shows the existence of de Sitter vacuum interior which has the properties of a perfect conductor and an ideal diamagnetic. The Kerr ring singularity (a naked singularity in the case without horizons) is replaced with a superconducting current, which serves as a nondissipative source of the Kerr-Newman fields and can be responsible for an unlimited life time of a spinning object.

In nonlinear electrodynamics coupled to gravity, regular spherically symmetric electrically charged solutions satisfy the weak energy condition and have obligatory de Sitter centre. By the G\"urses-G\"ursey algorithm they are transformed to spinning electrically charged solutions asymptotically Kerr-Newman for a distant observer. Rotation transforms de Sitter center into de Sitter vacuum surface which contains equatorial disk $r=0$ as a bridge. We present general analysis of the horizons, ergoregions and de Sitter surfaces, as well as the conditions of the existence of regular solutions to the field equations. We find asymptotic solutions and show that de Sitter vacuum surfaces have properties of a perfect conductor and ideal diamagnetic, violation of the weak energy condition is prevented by the basic requirement of electrodynamics of continued media, and the Kerr ring singularity is replaced with the superconducting current.

We present a family of spherically symmetric multi-horizon spacetimes with a vacuum dark fluid, associated with a time-dependent and spatially inhomogeneous cosmological term. The vacuum dark fluid is defined in a model-independent way by the symmetry of its stress-energy tensor, i.e., its invariance under Lorentz boosts in a distinguished spatial direction ($p_r=-\rho$ for spherical symmetry), which makes the dark fluid essentially anisotropic and allows its density to evolve. The related cosmological models belong to the Lemaitre class of models with anisotropic fluids and describe a universe with several scales of vacuum energy related to phase transitions during its evolution. The typical behavior of solutions and the number of spacetime horizons are determined by the number of vacuum scales. We study in detail a model with three vacuum scales: GUT, QCD and that responsible for the present accelerated expansion. The model parameters are fixed by the observational data and by analyticity and causality conditions. We find that our Universe has three horizons. During the first inflation the Universe enters a T-region which makes the expansion irreversible. After the second phase transition at the QCD scale the Universe enters an R-region, where for a long time its geometry remains almost pseudo-Euclidean. After crossing the third horizon related to the present vacuum density, the Universe should enter the next T-region with inevitable expansion.

We consider the process of decay of symmetric vacuum state as evaporation of a Bose condensate of physical Higgs particles, defined over asymmetric vacuum state. Energy density of their selfinteraction is identified with cosmological constant $\Lambda$ in the Einstein equation. $\Lambda$ decay then provides dynamical realization of spontaneous symmetry breaking. The effective mechanism is found for damping of coherent oscillations of a scalar field, leading to slow evaporation regime as the effective mechanism for $\Lambda$ decay responsible for inflation without special fine-tuning of the microphysical parameters. This mechanism is able to incorporate reheating, generation of proper primordial fluctuations, and nonzero cosmological constant today.