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).

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.

The method of a determination of the Primary Cosmic Ray mass composition is presented. Data processing is based on the theoretical model representing the integral muon multiplicity spectrum as the superposition of the spectra corresponding to different kinds of primary nuclei. The method consists of two stages. At the first stage, the permissible intervals of primary nuclei fractions f_i are determined on the base of the EAS spectrum vs the total number of muons (E_mu > 235 GeV). At the second stage, the permissible intervals of f_i are narrowed by fitting procedure. We use the experimental data on high multiplicity muon events (n_mu > 114) collected at the Baksan underground scintillation telescope. Within the framework of three components (protons, helium and heavy nuclei), the mass composition in the region 10^15-10^16 eV has been defined: f_p = 0.235 +- 0.02, f_He = 0.290 +- 0.02$, f_H = 0.475 +- 0.03.

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.

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).