We present the earliest astronomical observation of a high energy neutrino error box in which its variability was discovered after high-energy neutrinos detection. The one robotic telescope of the MASTER global international network (Lipunov et al. 2010) automatically imaged the error box of the very high-energy neutrino event IceCube-170922A. Observations were carried out in minute after the IceCube-170922A neutrino event was detected by the IceCube observatory at the South Pole. MASTER found the blazar TXS 0506+056 to be in the off-state after one minute and then switched to the on-state no later than two hours after the event. The effect is observed at a 50-sigma significance level. Also we present own unique 16-years light curve of blazar TXS 0506+056 (518 data set).

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.

Carpet is an air-shower array at Baksan, Russia, equipped with a large-area muon detector, which makes it possible to separate primary photons from hadrons. We report first results of the search for primary photons with energies E>100 TeV. The experiment's ongoing upgrade and future sensitivity are also discussed.

Using data of the Baksan Underground Scintillation Telescope (BUST) we have made a search for muon neutrinos and antineutrinos with energies above 1 GeV coinciding with the gravitational wave event GW170817 that was recorded on August 17, 2017 by the Advanced LIGO and Advanced Virgo observatories. This is a first detection of the new type of events occurring as a result of a merger of two neutron stars in a binary system. A short gamma-ray burst GRB170817A accompanying this event is an evidence of particle acceleration in the source whose precise position was determined by detection of the subsequent optical signal. No neutrino signals were found with the BUST in the interval $\pm 500$ s around the moment of the gravitational wave event GW170817, as well as during the next 14 days. The upper limits on integral fluxes of muon neutrino and antineutrino from the source are derived.

The experiment on recording neutrino bursts operates since the mid-1980. As the target, we use two parts of the facility with the total mass of 240 tons. The current status of the experiment and some results related to the investigation of background events and the stability of facility operation are presented. Over the period of June 30, 1980 to December 31, 2018, the actual observational time is 33.02 years. No candidate for the stellar core collapse has been detected during the observation period. An upper bound of the mean frequency of core collapse supernovae in our Galaxy is 0.070 year$^{-1}$ (90\% CL).

Early results of the search for E_\gamma>1 PeV cosmic photons from point sources with the data of Carpet-2, an air-shower array equipped with a 175 m^2 muon detector, are presented. They include 95% CL upper limits on PeV photon fluxes from stacked directions of high-energy IceCube neutrino events and from four predefined sources, Crab, Cyg X-3, Mrk 421 and Mrk 501. An insignificant excess of events from Mrk 421 will be further monitored. Prospects of the use of the upgraded installation, Carpet-3 (410 m^2 muon detector), scheduled to start data taking in 2019, for searches of E_\gamma>100 TeV photons, are briefly discussed.

Carpet-2 is an air-shower array at Baksan Valley, Russia, equipped with a large-area (175 m^2) muon detector, which makes it possible to separate primary photons from hadrons. We report the first results of the search for primary photons with energies E_\gamma>1 PeV, directionally associated with IceCube high-energy neutrino events, in the data obtained in 3080 days of Carpet-2 live time.

The matrices formed of silicon photomultipliers (SiPMs) are very promising photosensors for the scintillation detectors. The use of SiPM matrices with appropriate optical collector gives, in principle, a possibility to do a snapshot of glowing track of charged particle traversing a scintillator. The prototype of such scintillation detector is under development now in INR RAS. The preliminary results of characterization study of the matrix ArrayC-60035-64P-PCB (SensL company) for the prototype of such detector are presented.

The present-day status of the problem of searching for primary cosmic gamma rays at energies above 100 TeV is discussed, as well as a proposal for a new experiment in this field. It is shown that an increase of the area of the muon detector of the Carpet-2 air shower array up to 410 square meters, to be realized in 2016, will make this array quite competitive with past and existing experiments, especially at modest energies. Some preliminary results of measurements made with smaller area of the muon detector are presented together with estimates of expected results to be obtained with the coming large-area muon detector.

Observing a high-statistics neutrino signal from the supernova explosions in the Galaxy is a major goal of low-energy neutrino astronomy. The prospects for detecting all flavors of neutrinos and antineutrinos from the core-collapse supernova (ccSN) in operating and forthcoming large liquid scintillation detectors (LLSD) are widely discussed now. One of proposed LLSD is Baksan Large Volume Scintillation Detector (BLVSD). This detector will be installed at the Baksan Neutrino Observatory (BNO) of the Institute for Nuclear Research, Russian Academy of Sciences, at a depth of 4800 m.w.e. Low-energy neutrino astronomy is one of the main lines of research of the BLVSD.

According to the Baksan underground scintillation telescope data the barometric coefficients for six intervals of zenith angles of particle registration were obtained for ten years (2003-2012). The obtained barometric coefficients exceed an order of magnitude the theoretical values. Nevertheless, these results are in good agreement with the data obtained for the barometric effect of high-energy muons in a number papers of other authors. And there is a noticeable increase of the absolute values of the measured barometric coefficients with increase of the threshold energy of muons.

Differential and integral energy spectra of cosmic ray muons in the energy range from several TeV to ~ 1 PeV obtained by means of the analysis of multiple interactions of muons (pair meter technique) in the Baksan underground scintillation telescope (BUST) are presented. The results are compared with preceding BUST data on muon energy spectrum based on electromagnetic cascade shower measurements and depth-intensity curve analysis, with calculations for different muon spectrum models, and also with data of other experiments.

The multichannel data acquisition system is intended to be used in the EMMA experiment studying cosmic rays. The array will be in the Pihasalmi mine (central Finland) at a depth of about 85 m. The scintillator counters (SC-1) of the array are cast plastic scintillators with a wavelength of each SC-1 detector is 12.2{\times}12.2{\times}3.0 cm^3. 16 SC-1 detectors are placed in the metal case of 50.0{\times}50.0{\times}13.0 cm^3 dimension. Each case, called SC-16 detector, contains electronics of preliminary processing of signals and operating mode stabilization. The whole of the array will contain 96 detectors SC-16. It will make 1536 channels placed in three planes (48+24+24 detectors). The array will allow us to measure the time of flight of particles between SC-16 detectors and the coordinates of SC-1 fired detectors. This paper presents the function diagram of data acquisition system that includes electronic of detectors, the hodoscope pulse channels, the trigger block and VME blocks.

The main goal of the Andyrchy-BUST experiment is to study the primary cosmic rays spectrum and composition around the knee. The experimental data on the knee, as observed in the electromagnetic and high energy muon components, are presented. The electromagnetic component in our experiment is measured using the "Andyrchy" EAS array. High energy muon component (with 230 GeV threshold energy of muons) is measured using the Baksan Underground Scintillation Telescope (BUST). The location of the "Andyrchy" right above the BUST gives us a possibility for simultaneous measurements of both EAS components.

Current status and results of the experiment on recording neutrino bursts are presented. The observation time (since 1980) is 24.7 years. The upper bound of collapse frequency in our Galaxy is 0.093 $y^{-1}$ (90% CL).

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.

The lateral distribution function of high energy muons in EAS around the knee ($5.9\le \lg N_e \le 7.1$) has been measured for near vertical showers ($\theta \le 20^{\circ}$, effective muon threshold energy is 230 GeV). The measurements have been performed at the Baksan Underground Scintillation Telescope (BUST). The electromagnetic component is measured by the "Andyrchy" EAS array, located above the BUST. The knee in EAS size spectrum is found to be at $\lg N_e \approx 6.3$. The experimental results are compared with Monte Carlo simulations.

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.

22 Sep 2005 Swift-BAT triggered and located GRB050922C. The light curve shows the intense broad peak with $T_{90}$ of $(5 \pm 1)$ s. The Nordic Optical Telescope has obtained spectra of the afterglow with several absorption features corresponding to a redshift of $z = 2.17 \pm 0.03$. Observation of optical transient of GRB050922C was carried out in the R-band with the 60-cm telescope equipped with a CCD on Peak Terskol (North Caucasus). The OT magnitude was fading from R $\approx 16$ to $\approx 17.5$. Detection of an oscillatory phenomenon in the R post-burst light curve is described in this work. Analysis of the R data reveals coherent harmonic with a period of $0.0050 \pm 0.0003$ days (7.2 min) during observing run of about 0.05 days ($\sim 70$ min). Amplitude of oscillations is about 0.05 magnitude. The simplest model suggests that GRB050922C may result from tidal disruption of a white dwarf star by a black hole of about one thousand solar mass. The periodicity in the light curve can be identified with relativistic precession of an accretion disc.

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.

A new type of cosmic-ray experiment is under construction in the Pyh\"asalmi mine in the underground laboratory of the University of Oulu, Finland. It aims to study the composition of cosmic rays at and above the knee region. The experiment, called EMMA, will cover approximately 150 square-metres of detector area. The array is capable of measuring the multiplicity and the lateral distribution of underground muons, and the arrival direction of the air shower. The full-size detector is expected to run by the end of 2007.

The results of measuring diurnal variations of extensive air showers with energy above 100 TeV detected by the Andyrchi air shower array in 1995-2002 are presented. It is shown that solar diurnal wave of intensity is formed by pressure variations and Compton-Getting effect (due to the orbital motion of the Earth). The data on cosmic ray anisotropy of the Baksan underground scintillation telescope (BUST) are also presented for energies above 2.5 TeV covering the period 1982-1998. The diurnal sidereal wave measured by both Andyrchi and BUST is due to anisotropy of the cosmic rays.

The estimation of light velocity dispertion caused by quantum fluctuations influence on electromagnetic wave propogation in four-dimensional space-time is presented. Analytical cosmological solutions for the flat and open Universe for the case of cosmic vacuum model are obtained. The time delay of cosmic gamma-ray burst radiation is calculated. The delay is explained by propogation of dispersive electromagnetic wave in the expanding Universe. It is shown that the delay value depends on a model of the Universe. We conclude that to discern a model of the Universe measurement accuracy of parameter \Delta t / \Delta E_{\gamma} should be better than 10^{-5} s / MeV.

We present the experimental data on the knee, as observed in the electromagnatic and high energy muon (E\mu \ge 230 Gev) components.