Primordial black holes (PBHs) are black holes which may form in the early Universe through the gravitational collapse of primordial cosmological density fluctuations. Due to Hawking radiation these PBHs are supposed to evaporate by emitting particles. Recent developments in the experimental searching for evaporating PBHs in the local Universe are reviewed. The multimessenger techniques of searching for signals from evaporating PBHs are discussed.

We consider the inflationary model in which the inflaton $\phi$ couples to another scalar field $\chi$ via the interaction $g^2(\phi-\phi_0)^2\chi^2$ with a small coupling constant $g$ ($g^2 \sim 10^{-7}$). We assume that there is a sequence of "trapping points" $\phi_{0i}$ along the inflationary trajectory where particles of $\chi$-field become massless and are rather effectively produced. We calculate the power spectrum of inflaton field fluctuations originated from a backreaction of $\chi$-particles produced, using the Schwinger's "in-in" formalism. We show that the primary curvature power spectrum produced by these backreaction effects is blue, which leads to a strong overproduction of primordial black holes (PBHs) in subsequent radiation era.

We study the process of primordial black hole (PBH) formation at the beginning of radiation era for the cosmological scenario in which the inflaton is a pseudo-Nambu-Goldstone boson (axion) and there is a coupling of the inflaton with some gauge field. In this model inflation is accompanied by the gauge quanta production and a strong rise of the curvature power spectrum amplitude at small scales (along with non-Gaussianity) is predicted. We show that data on PBH searches can be used for a derivation of essential constraints on the model parameters in such an axion inflation scenario. We compare our numerical results with the similar results published earlier, in the work by Linde et al.

We consider the mixed inflaton-curvaton scenario in which quantum fluctuations of the curvaton field during inflation lead to a relatively large curvature perturbation spectrum at small scales. We use the model of chaotic inflation with quadratic potential including supergravity corrections leading to a large positive tilt in the power spectrum of the curvaton field. The model is characterized by the strongly inhomogeneous curvaton field in the Universe and large non-Gaussianity of curvature perturbations at small scales. We obtained the constraints on the model parameters considering the process of primordial black hole (PBH) production in radiation era.

We consider the early Universe scenario which allows for production of non-Gaussian curvature perturbations at small scales. We study the peculiarities of a formation of primordial black holes (PBHs) connected with the non-Gaussianity. In particular, we show that PBH constraints on the values of curvature perturbation power spectrum amplitude are strongly dependent on the shape of perturbations and can significantly (by two orders of magnitude) deviate from the usual Gaussian limit ${\cal P}_\zeta \lesssim 10^{-2}$. We give examples of PBH mass spectra calculations and PBH constraints for the particular case of the curvaton model.

We consider several early Universe models that allow for production of large curvature perturbations at small scales. As is well known, such perturbations can lead to formation of primordial black holes (PBHs). We briefly review the today's situation with PBH constraints and then focus on two models in which strongly non-Gaussian curvature perturbations are predicted: the hybrid inflation waterfall model and the curvaton model. We show that PBH constraints on the values of curvature perturbation power spectrum amplitude are strongly dependent on the shape of perturbations and can significantly (by two orders of magnitude) deviate from the usual Gaussian limit ${\cal P}_\zeta \lesssim 10^{-2}$. We give examples of PBH mass spectra calculations for both inflationary models.

We consider the process of primordial black hole (PBH) formation originated from primordial curvature perturbations produced during waterfall transition (with tachyonic instability), at the end of hybrid inflation. It is known that in such inflation models, rather large values of curvature perturbation amplitudes can be reached, which can potentially cause a significant PBH production in the early Universe. The probability distributions of density perturbation amplitudes in this case can be strongly non-Gaussian, which requires a special treatment. We calculated PBH abundances and PBH mass spectra for the model, and analyzed their dependence on model parameters. We obtained the constraints on the parameters of the inflationary potential, using the available limits on $\beta_{PBH}$.

We carried out numerical calculations of a contribution of the waterfall field to the primordial curvature perturbation (on uniform density hypersurfaces) $\zeta$, which is produced during waterfall transition in hybrid inflation scenario. The calculation is performed for a broad interval of values of the model parameters. We show that there is a strong growth of amplitudes of the curvature perturbation spectrum in the limit when the bare mass-squared of the waterfall field becomes comparable with the square of Hubble parameter. We show that in this limit the primordial black hole constraints on the curvature perturbations must be taken into account. It is shown that, in the same limit, peak values of the curvature perturbation spectra are far beyond horizon, and the spectra are strongly non-Gaussian.

We perform a consistent calculation of primordial black hole (PBH) mass spectrum and second-order induced gravitational wave (GW) background produced from primordial scalar perturbations in radiation era of the early Universe. It is shown that the maximal amplitudes of the second order GW spectrum that can be approached without conflicting with the PBH data do not depend significantly on the shape of primordial perturbation spectrum. The constraints on the GW background obtained in previous works are extended to a wider GW frequency range. We discuss the applicability of the currently available pulsar timing limits for obtaining the constraints on scalar power spectrum and PBH abundance and show that they can be used for strongly constraining the PBH number density in the PBH mass range $\sim (0.03 - 10) M_{\odot}$.

We study, by numerical methods, the time evolution of scalar perturbations in radiation era of Randall-Sundrum braneworld cosmology. Our results confirm an existence of the enhancement of perturbation amplitudes (near horizon crossing), discovered recently. We suggest the approximate solution of equations of the perturbation theory in the high energy regime, which predicts that the enhancement factor is asymptotically constant, as a function of scale. We discuss the application of this result for the problem of primordial black hole production in braneworld cosmology.

We calculate the frequency dependence of gravitational wave background arising at second order of cosmological perturbation theory due to mixing of tensor and scalar modes. The calculation of the induced gravitational background is performed for two special cases: for the power spectrum of scalar perturbations which has a peak at some scale and for the scalar spectrum predicted by the inflationary model with the running mass potential. We show that the amplitudes of the induced gravitational background, in the frequency region $\sim 10^{-3} - 10^3 $Hz, are effectively constrained by results of studies of primordial black hole production in early universe. We argue that though today's LIGO bound on ${\cal P}_{\cal R}(k)$ is weaker than the PBH one, Advanced LIGO will be able to set a stronger bound, and in future the ground-based interferometers of LIGO type will be suitable for obtaining constraints on PBH number density in the mass range $\sim 10^{11} - 10^{15}$ g.

The today's energy density of the induced (second order) gravitational wave background in the frequency region $\sim 10^{-3} - 10^3$ Hz is constrained using the existing limits on primordial black hole production in the early Universe. It is shown, in particular, that at frequencies near $\sim 40$ Hz (which is the region explored by LIGO detector), the value of the induced part of $\Omega_{GW}$ cannot exceed $(1-3)\times 10^{-7}$. The spread of values of the bound is caused by the uncertainty in parameters of the gravitational collapse of black holes.

The technical procedure of a search of bursts of very high energy gamma rays from evaporation of primordial black holes on air-shower array "Andyrchy" of Baksan Neutrino Observatory of Institute for Nuclear Research is described. The theoretical model used in the present work assumes that the chromosphere around the evaporating black hole does not form. For minimization of the cosmic ray background the method of multidimensional analysis of modelled as well as experimentally detected events is applied. The new upper limit on the concentration of evaporating primordial black holes in the local region of Galaxy is obtained. The comparison of the results of different experiments is given.

Experimental data of arrays "Andyrchy" and "Carpet-2" of Baksan Neutrino Observatory (Institute for Nuclear Research), obtained in the regime of a detection of the single cosmic-ray component, are used for a search of the bursts of cosmic gamma rays from evaporating primordial black holes. Different theoretical models of the evaporation process are used for the analysis. Distributions of the counting rate fluctuations on both arrays agree with the expectations from the cosmic ray background. The new constraints on the concentration of evaporating primordial black holes in the local region of Galaxy are obtained. The comparison of the results of different experiments is given.

We calculate the primordial black hole (PBH) mass spectrum produced from a collapse of the primordial density fluctuations in the early Universe using, as an input, several theoretical models giving the curvature perturbation power spectra with large (~ 0.01 - 0.1) values at some scale of comoving wave numbers k. In the calculation we take into account the explicit dependence of gravitational (Bardeen) potential on time. Using the PBH mass spectra, we further calculate the neutrino and photon energy spectra in extragalactic space from evaporation of light PBHs, and the energy density fraction contained in PBHs today (for heavier PBHs). We obtain the constraints on the model parameters using available experimental data (including data on neutrino and photon cosmic backgrounds). We briefly discuss the possibility that the observed 511 keV line from the Galactic center is produced by annihilation of positrons evaporated by PBHs.

We consider the examples of single-field inflation models predicting large amplitudes of the curvature perturbation power spectrum at relatively small scales. It is shown that in models with an inflationary potential of double-well type the peaks in the power spectrum, having, in maximum, the amplitude ~0.1, can exist (if parameters of the potential are chosen appropriately). It is shown also that the spectrum amplitude of the same magnitude (at large k values) is predicted in the model with the running mass potential, if the positive running, n', exists and is about 0.005 at cosmological scales. Estimates of the quantum diffusion effects during inflation in models with the running mass potential are given.

Temporal and energy characteristics of the very-high-energy gamma-ray bursts from evaporating primordial black holes have been calculated by assuming that the photospheric and chromospheric effects are negligible. The technique of searching for such bursts on shower arrays is described. We show that the burst time profile and the array dead time should be taken into account to interpret experimental data. Based on data from the Andyrchy array of the Baksan Neutrino Observatory (Institute for Nuclear Research, Russian Academy of Sciences), we have obtained an upper limit on the number density of evaporating primordial black holes in a local region of space with a scale size of ~10^{-3} pc. Comparison with the results of previous experiments is made.

Cosmic gamma-ray bursts with primary-photon energies > 10 GeV are sought in the data from the Andyrchy array obtained in the mode of detection of a single cosmic-ray component during a net observation period of 2005.4 d. The distribution of fluctuations of the detector counting rate agrees with the expected cosmic-ray background, the only exception being an event with a deviation of 7.9 sigma. Constraints on the number density of evaporating primordial black holes in a local region of the Galaxy are obtained for the chromospheric evaporation models.

Possibilities of an experimental search for gamma-ray bursts from primordial black hole (PBH) evaporations in space are reconsidered. It is argued that the corresponding constraints which can be obtained in experiments with cosmic ray detectors strongly depend on theoretical approach used for a description of the PBH evaporation process. Predictions of several theoretical models for gamma-ray spectra from final stages of PBH life (integrated over time) are given.

We calculate neutrino and photon energy spectra in extragalactic space from evaporation of primordial black holes, assuming that the power spectrum of primordial density fluctuations has a strong bump in the region of small scales. The constraints on the parameters of this bump based on neutrino and photon cosmic background data are obtained.

We consider the modification of extragalactic cosmic ray spectrum caused by cosmic ray interactions with infrared background photons which are present in the extragalactic space together with relic photons. It is assumed that cosmic ray spectrum at superhigh energies has extragalactic origin and is proton dominated.

Diffuse background of high energy neutrinos arising from interactions of cosmic ray protons with far infrared radiation background in extragalactic space is calculated. It is assumed that cosmic ray spectrum at superhigh energies has extragalactic origin and is proton dominated.The cosmological evolution of extragalactic sources of cosmic ray protons as well as infrared-luminous galaxies is taken into account in the calculation.

A review of the recent achievements in high energy neutrino physics and, partly, neutrino astrophysics is presented. It is argued that experiments with high energy neutrinos of natural origin can be used for a search of new physics effects beyond the electroweak scale.

Extragalactic background of high energy neutrinos arising from the interactions of cosmic ray protons with far-infrared extragalactic background radiation is calculated. The main assumption is that the cosmic ray spectrum at energies higher than $10^8 GeV$ has extragalactic origin and, therefore, is proton dominated. All calculations are performed with taking into account the possible cosmological evolution of extragalactic sources of cosmic ray protons as well as infrared-luminous galaxies.

We obtain the constraints on the diffuse neutrino background from primordial black hole evaporations using the recent data of KamLAND experiment setting an upper limit on the antineutrino flux from the Sun.

We present results of the detailed Monte Carlo calculation of the rates of double-bang events in 1 km$^3$ underwater neutrino telescope with taking into account the effects of $\tau$-neutrino propagation through the Earth. As an input, the moderately optimistic theoretical predictions for diffuse neutrino spectra of AGN jets are used.

If neutrinos produced in astrophysics sources oscillate, neutrino telescopes can have a chance to detect signals from tau-neutrinos. We present an estimate of the 'double bang' event rate produced by tau neutrinos in cubic kilometer scale detectors.

The formula for the initial mass spectrum of primordial black holes (PBHs), which can be used for a general case of the scale dependent spectral index, and for a wide class of models of the gravitational collapse, is derived. The derivation is based on the Press and Schechter formalism. The comparative analysis of different types of initial mass spectra used in concrete calculations is carried out. It is shown that densities of background radiations ($\nu$, $\gamma$) from PBH evaporations depend rather strongly on a type of the gravitational collapse and on a taking into account the spread of horizon masses at which PBHs can form. Constraints on parameters of the primordial density perturbation amplitudes based on PBH evaporation processes and on atmospheric and solar neutrino data are obtained.

Comparative analysis of different types of initial mass spectra of primordial black holes (PBHs) is carried out. It is assumed that primordial density fluctuations have power law spectrum and the spectral index is constant throughout all scales. It is shown that densities of background radiations (neutrinos, gamma-quanta) from PBH evaporations strongly depend on the type of gravitational collapse and on the way of taking into account the spread of horizon masses at which PBHs can form. Constraints on spectral index values based on PBH evaporation process and on atmospheric and solar neutrino data are obtained.

The process of baryogenesis through the evaporation of black holes formed at the end of inflation phase is considered. The increase of black hole mass due to accretion from the surrounding radiation after the reheating is taken into account. It is shown that the influence of the accretion on the baryogenesis is important only in the case when the initial values of black hole mass are larger than $\sim 10^{4}{g}$. The behavior of calculated baryon asymmetry, as a function of model parameters, is studied.

We calculated the energy spectra and the fluxes of electron neutrinos in extragalactic space emitted in the process of the evaporation of primordial black holes (PBHs) in the early universe. It was assumed that PBHs are formed by a blue power-law spectrum of primordial density fluctuations. In the calculations of neutrino spectra the spectral index of density fluctuations and the reheating temperature were used as free parameters. The absorption of neutrinos during propagation in the space was taken into account. We obtained the bounds on the spectral index assuming validity of the standard picture of gravitational collapse and using the available data of several experiments with atmospheric and solar neutrinos. The comparison of our results with the previous constraints (which had been obtained using diffuse photon background data) shows that such bounds are quite sensitive to an assumed form of the initial PBH mass function.

We calculated the energy spectra and the fluxes of electron neutrino emitted in the process of evaporation of primordial black holes (PBHs) in the early universe. It was assumed that PBHs are formed by a blue power-law spectrum of primordial density fluctuations. We obtained the bounds on the spectral index of density fluctuations assuming validity of the standard picture of gravitational collapse and using the available data of several experiments with atmospheric and solar neutrinos. The comparison of our results with the previous constraints (which had been obtained using diffuse photon background data) shows that such bounds are quite sensitive to an assumed form of the initial PBH mass function.

The vertical sea-level muon spectrum at energies above 1 GeV and the underground/underwater muon intensities at depths up to 18 km w.e. are calculated. The results are particularly collated with a great body of the ground-level, underground, and underwater muon data. In the hadron-cascade calculations, the growth with energy of inelastic cross sections and pion, kaon, and nucleon generation in pion-nucleus collisions are taken into account. For evaluating the prompt muon contribution to the muon flux, we apply two phenomenological approaches to the charm production problem: the recombination quark-parton model and the quark-gluon string model. To solve the muon transport equation at large depths of homogeneous medium, a semi-analytical method is used. The simple fitting formulas describing our numerical results are given. Our analysis shows that, at depths up to 6-7 km w. e., essentially all underground data on the muon intensity correlate with each other and with predicted depth-intensity relation for conventional muons to within 10%. However, the high-energy sea-level data as well as the data at large depths are contradictory and cannot be quantitatively decribed by a single nuclear-cascade model.

Supersymmetric models allow for stable non-topological solitons, Q-balls, which can be produced in the early Universe and contribute to dark matter. Experimental signature of electrically neutral Q-balls is, in fact, the same as is expected for superheavy magnetic monopoles catalyzing baryon decay. Here we use the upper limits on monopole flux obtained with deep underwater Cherenkov array "Gyrlyanda" which operated in the Baikal lake in 1984-90 with 267 days of live time to obtain the limit on Q-ball flux. The last has been found to be equal to 3.9 x 10^{-16} cm^{-2} sr^{-1} s^{-1} (90% CL). This result is discussed and compared with other restrictions.