Purposeful investigation of radiation fluxes strongly delayed in relation to the main particles front of extensive air shower (EAS) was undertaken at the Tien Shan Mountain Cosmic Ray Station. It was found that the passage of the EAS can be accompanied by the delayed thermal neutrons and by the soft $(30-50)$ keV gamma rays, mostly concentrated within a region of about $(5-10)$ m around shower axis, where the integral radiation fluence can vary in the limits of $(10^{-4}-1)$ cm$^{-2}$ for neutrons, and of $(0.1-1000)$ cm$^{-2}$ for gamma rays. The dependence of signal multiplicity on the shower size $N_e$ has a power shape both for the neutron and gamma ray components, with a sharp increase of its power index around the value of $N_e\approx 10^6$, which corresponds to the position of the $3\cdot10^{15}$ eV knee in the primary cosmic ray spectrum. Total duration of detectable radiation signal after the EAS passage can be of some tens of milliseconds in the case of neutron component, and up to a few whole seconds for gamma rays. The delayed accompaniment of low-energy radiation particles can be an effective probe to study the interaction of the hadronic component of EAS.

We present a description of the new complex installation for the study of extensive air showers which was created at the Tien Shan mountain cosmic ray station, as well as the results of the test measurements made there in 2014-2016. At present, the system for registration of electromagnetic shower component consists of $\sim$100 detector points built on the basis of plastic scintillator plates with the sensitive area of 0.25m$^2$ and 1m$^2$, spread equidistantly over $\sim$10$^4$m$^2$ space. The dynamic range of scintillation amplitude measurements is currently about $(3-7)\cdot 10^4$, and there is a prospect of it being extended up to $\sim$10$^6$. The direction of shower arrival is defined by signal delays from a number of the scintillators placed cross-wise at the periphery of the detector system. For the investigation of nuclear active shower components there was created a multi-tier 55m$^2$ ionization-neutron calorimeter with a sum absorber thickness of $\sim$1000g/cm$^2$, typical spatial resolution of the order of 10cm, and dynamic range of ionization measurement channel about $\sim$10$^5$. Also, the use of saturation-free neutron detectors is anticipated for registration of the high- and low-energy hadron components in the region of shower core. A complex of underground detectors is designed for the study of muonic and penetrative nuclear-active components of the shower. The full stack of data acquisition, detector calibration, and shower parameters restoration procedures are now completed, and the newly obtained shower size spectrum and lateral distribution of shower particles occur in agreement with conventional data. Future studies in the field of $10^{14}-10^{17}$eV cosmic ray physics to be held at the new shower installation are discussed.

The events of multiple neutron production under 2000g/cm$^2$ thick rock absorber were studied at the Tien~Shan mountain cosmic ray station, at the altitude of 3340m above the sea level. From comparison of the experimental and Geant4 simulated neutron multiplicity spectra it follows that the great bulk of these events can be explained by interaction of cosmic ray muons with internal material of the neutron detector. In synchronous operation of the underground neutron monitor with the Tien~Shan shower detector system it was found that the characteristics of the muonic component of extensive air showers which is seemingly responsible for generation of the neutron events underground do change noticeably within the energy range of the knee of primary cosmic ray spectrum. Some peculiar shower events were detected when the neutron signal reveals itself only $\sim$(100--1000)\,$\mu$s after the passage of the shower particles front which probably means an existence of corresponding delay of the muon flux in such events.

We reanalyze the dataset collected during the years 1998--2003 by the deep underwater neutrino telescope NT200 in the lake Baikal with the low energy threshold (10 GeV) in searches for neutrino signal from dark matter annihilations near the center of the Milky Way. Two different approaches are used in the present analysis: counting events in the cones around the direction towards the Galactic Center and the maximum likelihood method. We assume that the dark matter particles annihilate dominantly over one of the annihilation channels $b\bar{b}$, $W^+W^-$, $\tau^+\tau^-$, $\mu^+\mu^-$ or $\nu\bar{\nu}$. No significant excess of events towards the Galactic Center over expected neutrino background of atmospheric origin is found and we derive 90% CL upper limits on the annihilation cross section of dark matter.

We analyse sensitivity of the gigaton volume telescope Baikal-GVD for detection of neutrino signal from dark matter annihilations or decays in the Galactic Center. Expected bounds on dark matter annihilation cross section and its lifetime are found for several annihilation/decay channels.

We have analyzed a data set taken over 2.76 years live time with the Baikal neutrino telescope NT200. The goal of the analysis is to search for neutrinos from dark matter annihilation in the center of the Sun. Apart from the conventional annihilation channels $b\bar{b}$, $W^+W^-$ and $\tau^+\tau^-$ we consider also the annihilation of dark matter particles into monochromatic neutrinos. From the absence of any excess of events from the direction of the Sun over the expected background, we derive 90% upper limits on the fluxes of muons and muon neutrinos from the Sun, as well as on the elastic cross sections of dark matter scattering on protons.

The Prototyping phase of the BAIKAL-GVD project has been started in April 2011 with the deployment of a three string engineering array which comprises all basic elements and systems of the Gigaton Volume Detector (GVD) in Lake Baikal. In April 2012 the version of engineering array which comprises the first full-scale string of the GVD demonstration cluster has been deployed and operated during 2012. The first stage of the GVD demonstration cluster which consists of three strings is deployed in April 2013. We review the Prototyping phase of the BAIKAL-GVD project and describe the configuration and design of the 2013 engineering array.

The pure leptonic decay B_s -> mu mu is strongly suppressed in the Standard Model (SM), but can have large enhancements in Supersymmetry, especially at large values of tanbe. New limits on this decay channel from recent LHC data have been used to claim that these limits restrict the SUSY parameter space even more than the direct searches. However, direct searches are hardly dependent on tanbe, while BR(B_s -> mu mu) is proportional to tanbe^6. The relic density constraint requires large tanbe in a large region of the parameter space, which can lead to large values of B_s -> mu mu. Nevertheless, the experimental upper limit on BR(B_s -> mu mu) is not constraining the parameter space of the CMSSM more than the direct searches and the present Higgs limits, if combined with the relic density. We also observe SUSY parameter regions with negative interferences, where the B_s -> mu mu value is up to a factor three below the SM expectation, even at large values of tanbe.

We report theoretical and experimental results of on-going feasibility studies to detect cosmic neutrinos acoustically in Lake Baikal. In order to examine ambient noise conditions and to develop respective pulse detection techniques a prototype device was created. The device is operating at a depth of 150 m at the site of the Baikal Neutrino Telescope and is capable to detect and classify acoustic signals with different shapes, as well as signals from neutrino-induced showers.

The excess of diffuse galactic gamma rays above 1 GeV, as observed by the EGRET telescope on the NASA Compton Gamma Ray Observatory, shows all the key features from Dark Matter (DM) annihilation: (i) the energy spectrum of the excess is the same in all sky directions and is consistent with the gamma rays expected for the annihilation of WIMPs with a mass between 50-100 GeV; (ii) the intensity distribution of the excess in the sky is used to determine the halo profile, which was found to correspond to the usual profile from N-body simulations with additional substructure in the form of two doughnut-shaped structures at radii of 4 and 13 kpc; (iii) recent N-body simulations of the tidal disruption of the Canis Major dwarf galaxy show that it is a perfect progenitor of the ringlike Monoceros tidal stream of stars at 13 kpc with ring parameters in agreement with the EGRET data; (iiii) the mass of the outer ring is so large, that its gravitational effects influence both the gas flaring and the rotation curve of the Milky Way. Both effects are clearly observed in agreement with the DMA interpretation of the EGRET excess.

Some part of the relic Dark Matter is distributed in small-scale clumps which survived structure formation in inflation cosmological scenario. The annihilation of DM inside these clumps is a strong source of stable charged particles which can have a substantial density near the clump core. The streaming of the annihilation products from the clump can enhance irregularities in the galactic magnetic field. This can produce small scale variations in diffusion coefficient affecting propagation of Cosmic Rays.

The extragalactic background (EGB) of diffuse gamma rays can be determined by subtracting the Galactic contribution from the data. This requires a Galactic model (GM) and we include for the first time the contribution of dark matter annihilation (DMA), which was previously proposed as an explanation for the EGRET excess of diffuse Galactic gamma rays above 1 GeV. In this paper it is shown that the newly determined EGB shows a characteristic high energy bump on top of a steeply falling soft contribution. The bump is shown to be compatible with a contribution from an extragalactic DMA signal from weakly interacting massive particles (WIMPs) with a mass between 50 and 100 GeV in agreement with the EGRET excess of the Galactic diffuse gamma rays and in disagreement with earlier analysis. The remaining soft contribution of the EGB is shown to resemble the spectra of the observed point sources in our Galaxy.

The diffuse galactic EGRET gamma ray data show a clear excess for energies above 1 GeV in comparison with the expectations from conventional galactic models. This excess shows all the features expected from Dark Matter WIMP Annihilation: a)it is present and has the same spectrum in all sky directions, not just in the galactic plane, as expected for WIMP annihilation b) it shows an interesting substructure in the form of a doughnut shaped ring at 14 kpc from the centre of the galaxy, where a ring of stars indicated the probable infall of a dwarf galaxy. From the spectral shape of the excess the WIMP mass is estimated to be between 50 and 100 GeV, while from the intensity the halo profile is reconstructed, which is shown to explain the peculiar change of slope in the rotation curve at about 11 kpc (due to the ring of DM at 14 kpc). Recently it was claimed by Bergstrom et al. that the DM interpretation of the EGRET gamma ray excess is excluded by the antiproton fluxes, since in their propagation model with isotropic diffusion the flux of antiprotons would be far beyond the observed flux. However, the propagation can be largely anisotropic, because of the convection of particles perpendicular to the disc and inhomogeneities in the local environment. It is shown that anisotropic propagation can reduce the antiproton yield by an order of magnitude, while still being consistent with the B/C ratio. Therefore it is hard to use antiprotons to search for {\it light} DM particles, which yield a similar antiproton spectrum as the background, but the antiprotons are a perfect means to tune the many degenerate parameters in the propagation models.

Recently it was shown that the excess of diffuse Galactic gamma rays above 1 GeV traces the Dark Matter halo, as proven by reconstructing the peculiar shape of the rotation curve of our Galaxy from the gamma ray excess. This can be interpreted as a Dark Matter annihilation signal. In this paper we investigate if this interpretation is consistent with Supersymmetry. It is found that the EGRET excess combined with all electroweak constraints is fully consistent with the minimal mSUGRA model for scalars in the TeV range and gauginos below 500 GeV.

Elsaesser and Mannheim fit a contribution of Dark Matter Annihilation (DMA) to the extragalactic contribution of the galactic diffuse gamma ray flux, as deduced from the EGRET data by Strong, Moskalenko and Reimer.They find a WIMP mass of 515{+110}{-75} GeV and quote a systematic error of 30\%. However, they do not include large systematic uncertainties from the fact that the determination of the extragalactic flux (EGF) requires a model for the subtraction of the Galactic flux from the data.The data used were obtained with a model without Galactic DM, so one expects additional uncertainty in the region where DMA contributes. Including a Galactic DMA contribution reduces the significance and the WIMP mass. The latter then becomes compatible with the Galactic excess of diffuse gamma rays, which posseses all the properties of DMA with a much higher significance than the extragalactic excess.

The public data from the EGRET space telescope on diffuse Galactic gamma rays in the energy range from 0.1 to 10 GeV are reanalyzed with the purpose of searching for signals of Dark Matter annihilation (DMA). The analysis confirms the previously observed excess for energies above 1 GeV in comparison with the expectations from conventional Galactic models. In addition, the excess was found to show all the key features of a signal from Dark Matter Annihilation (DMA): a) the excess is observable in all sky directions and has the same shape everywhere, thus pointing to a common source; b) the shape corresponds to the expected spectrum of the annihilation of non-relativistic massive particles into - among others - neutral $\pi^0$ mesons, which decay into photons. From the energy spectrum of the excess we deduce a WIMP mass between 50 and 100 GeV, while from the intensity of the excess in all sky directions the shape of the halo could be reconstructed. The DM halo is consistent with an almost spherical isothermal profile with substructure in the Galactic plane in the form of toroidal rings at 4 and 14 kpc from the center. These rings lead to a peculiar shape of the rotation curve, in agreement with the data, which proves that the EGRET excess traces the Dark Matter.

The diffuse galactic EGRET gamma ray data show a clear excess for energies above 1 GeV in comparison with the expectations from conventional galactic models. The excess is seen with the same spectrum in all sky directions, as expected for Dark Matter (DM) annihilation. This hypothesis is investigated in detail. The energy spectrum of the excess is used to limit the WIMP mass to the 50-100 GeV range, while the sky maps are used to determine the halo structure, which is consistent with a triaxial isothermal halo with additional enhancement of Dark Matter in the disc. The latter is strongly correlated with the ring of stars around our galaxy at a distance of 14 kpc, thought to originate from the tidal disruption of a dwarf galaxy. It is shown that this ring of DM with a mass of $\approx 2\cdot 10^{11} M_\odot$ causes the mysterious change of slope in the rotation curve at $R=1.1R_0$ and the large local surface density of the disc. The total mass of the halo is determined to be $3\cdot 10^{12} M_\odot$. A cuspy profile is definitely excluded to describe the gamma ray data. These signals of Dark Matter Annihilation are compatible with Supersymmetry for boost factors of 20 upwards and have a statistical significance of more than $10\sigma$ in comparison with the conventional galactic model. The latter combined with all features mentioned above provides an intriguing hint that the EGRET excess is indeed a signal from Dark Matter Annihilation.

The process of sunspot formation is considered with the account of heat effects. According to the Le Chatelier principle, a local overheating must precede to the cooling of solar surface in the places of sunspot formation. The sunspot dynamics is a process close to the surface nucleate-free boiling in a thin layer with formation of bubbles (or craters), so we focus on the analogy between these two processes. Solar spots and surface nucleate-free boiling in a thin layer have similarities in formation conditions, results of impact on the surface were they have been formed, periodicity, and their place in the hierarchy of self-organization in complex systems. The difference is in the working medium and method of channelling of extra energy from the overheated surface -for boiling process, the energy is forwarded to generation of vapor, and in sunspots the solar energy is consumed to formation of a strong magnetic field. This analogy explains the problem of a steady brightness (temperature) of a spot that is independent of the spot size and other characteristics.