The gravitational-wave event S190521g -- a likely binary black hole merger in the accretion disk of an active galactic nucleus -- was accompanied by an optical counterpart. Such dense environments around luminous energy release regions are favourable for high energy $\gamma$-ray production. We report on a search for high energy $\gamma$-rays from the direction of the candidate electromagnetic counterpart to the S190521g event using publicly-available data of the Fermi-LAT space $\gamma$-ray telescope. No significant signal was found. We present upper limits on the spectral energy distribution of the source in the 100 MeV - 300 GeV energy range. We discuss the importance of studying S190521g-like transients in the context of cosmic ray acceleration, $\gamma$-ray and neutrino production in such sources.

Primary very high energy $\gamma$-rays from $\gamma$-ray bursts (GRBs) are partially absorbed on extragalactic background light (EBL) photons with subsequent formation of intergalactic electromagnetic cascades. Characteristics of the observable cascade $\gamma$-ray signal are sensitive to the strength and structure of the extragalactic magnetic field (EGMF). GRB 190114C was recently detected with the MAGIC imaging atmospheric Cherenkov telescopes, for the first time allowing to estimate the observable cascade intensity. We inquire whether any constraints on the EGMF strength and structure could be obtained from publicly-available $\gamma$-ray data on GRB 190114C. We present detailed calculations of the observable cascade signal for various EGMF configurations. We show that the sensitivity of the Fermi-LAT space $\gamma$-ray telescope is not sufficient to obtain such constraints on the EGMF parameters. However, next-generation space $\gamma$-ray observatories such as MAST would be able to detect pair echoes from GRBs similar to GRB 190114C for the EGMF strength below 10^{-17}--10^{-18} G.

The SPHERE-2 balloon-borne detector designed for extensive air shower (EAS) observations using EAS optical Vavilov-Cherenkov radiation (``Cherenkov light''), reflected from the snow-covered surface of Lake Baikal is described. We briefly discuss the concept behind the reflected Cherenkov light method, characterize the conditions at the experimental site and overview the construction of the tethered balloon used to lift the SPHERE-2 telescope above the surface. This paper is mainly dedicated to a detailed technical description of the detector, including its optical system, sensitive elements, electronics, and data acquisition system (DAQ). The results of some laboratory and field tests of the optical system are presented.

Context. Ultrahigh energy protons and nuclei from extragalactic cosmic ray sources produce secondary $\gamma$-rays, electrons and positrons on cosmic microwave background (CMB) and extragalactic background light (EBL) photons. These secondary particles initiate intergalactic electromagnetic (EM) cascades, eventually resulting in observable fluxes of $\gamma$-rays in the GeV-TeV energy domain. The total spectrum of such cascade $\gamma$-rays of hadronic nature is significantly harder than the one usually expected from blazars. The spectra of some sources known as "extreme TeV blazars" were claimed to be well-described by this "intergalactic hadronic cascade model". Aims. We calculate the shape of the observable spectrum inside the point spread function (PSF) of a typical imaging atmospheric Cherenkov telescope (IACT), as well as the observable angular distibution of $\gamma$-rays, for the first time taking into account the effect of primary proton deflection in a realistic extragalactic magnetic field (EGMF). We investigate if these deflections significantly change the shape of the observable spectrum and whether the observable angular distibution could be detected with currently operating or projected $\gamma$-ray telescopes. Methods. We estimate the width of the observable $\gamma$-ray angular distribution from simple geometrical considerations. In addition, we employ a hybrid code previously delevoped by us, supplementing it with new routines to compute the observable angle for detectable $\gamma$-rays. Results. The observable point-like spectrum at milti-TeV energies is much softer than the one averaged over all values of the observable angle. The observable angular distribution is found to be sufficiently broad to enable a robust identification of the extended emission around the sources with next-generation $\gamma$-ray telescopes. \abridged

We explore the concept of liquid Argon time projection chamber (TPC) for gamma-ray astronomy in the 100 MeV --- 1 TeV energy range. We propose a basic layout for such a telescope called MAST. Using a last-generation rocket such as Falcon Heavy, it is possible to launch a detector with the effective area and the differential sensitivity about one order of magnitude better than the Fermi-LAT ones. At the same time, the MAST concept allows for an excellent angular resolution, 3-10 times better than the Fermi-LAT one depending on the energy, and good energy resolution ($\approx$ 20 % at 100 MeV and 6-10 % for the 10 GeV --- 1 TeV energy range). We show that such a telescope would be instrumental in a broad range of long-standing astrophysical problems.

A compact device lifted over the ground surface might be used to observe optical radiation of extensive air showers (EAS). Here we consider spatial and temporal characteristics of Vavilov-Cherenkov radiation ("Cherenkov light") reflected from the snow surface of Lake Baikal, as registered by the SPHERE-2 detector. We perform detailed full direct Monte Carlo simulations of EAS development and present a dedicated highly modular code intended for detector response simulations. Detector response properties are illustrated by example of several model EAS events. The instrumental acceptance of the SPHERE-2 detector was calculated for a range of observation conditions. We introduce the concept of "composite model quantities", calculated for detector responses averaged over photoelectron count fluctuations, but retaining EAS development fluctuations. The distortions of EAS Cherenkov light lateral distribution function (LDF) introduced by the SPHERE-2 telescope are understood by comparing composite model LDF with the corresponding function as would be recorded by an ideal detector situated at the ground surface. We show that the uncertainty of snow optical properties does not change our conclusions, and, moreover, that the expected performance of the SPHERE experiment in the task of cosmic ray mass composition study in the energy region $\sim$10 PeV is comparable with other contemporary experiments. Finally, we compare the reflected Cherenkov light method with other experimental techniques and briefly discuss its prospects.

Preliminary results on the development of a separation method for Cerenkov (CL) and fluorescence (FL) light from EAS are shown. The results are based on the measurement of attenuation coefficients of CL and FL for different filters. A total of six optical filters were investigated: filters from optical glass UFS-1, UFS-5, FS6 (analogue BG3) and interference filters SL 360-50, SL 280-380, FF01-375/110. The measurements were performed using silicon photomultipliers (SiPM). To improve existing fluorescent light detectors, a segment of 7 SiPM was developed, which would be able to separate both components of the light flux from EAS at the level of primary data processing.

We review extragalactic $\gamma$-ray propagation models with emphasis on the electromagnetic (EM) cascade process in the magnetized expanding Universe. We consider cascades initiated by primary protons of ultra-high energy accelerated by blazars and show that the observable spectrum is similar to the universal spectrum of a purely EM cascade. We also present a detailed calculation of the observable angular distribution for the case of EM cascades developing from relatively nearby (<20 Mpc) sources. Finally, we calculate the point-like source differential sensitivity of a novel liquid Argon time projection chamber $\gamma$-ray telescope and show that its sensitivity is several times better than the Fermi LAT sensitivity in the 100 MeV -- 100 GeV energy range.

We review the physics of intergalactic electromagnetic cascades in the presence of the extragalactic magnetic field (EGMF). Various regimes of intergalactic electromagnetic cascades are considered depending on the number of cascade generations, the value of the cascade electron deflection angle, and the relations between the EGMF coherence length, typical cascade $\gamma$-ray mean free path, and electron energy loss length. We also review contemporary constraints on the EGMF parameters and explore the sensitivity of various $\gamma$-ray instruments to the EGMF parameters.

Secondary {\gamma}-rays from intergalactic cascades may contribute to observable spectra of blazars, also modifying observable angular and temporal distributions. In this paper we briefly review basic features of intergalactic electromagnetic cascade physics, suggest a new approximation for {\gamma}-ray mean free path, consider angular patterns of magnetically broadened cascade emission, and present an example of a fit to the observable blazar spectrum.

Recent progress in very high energy (VHE, E >100 GeV) $\gamma$-ray observations, together with advances in the extragalactic background light (EBL) modelling, allows to search for new phenomena such as $\gamma$-axion-like particle ($\gamma \rightarrow$ ALP) oscillation and to explore the extragalactic magnetic field (EGMF) strength and structure. These studies are usually performed by searching for some deviation from the so-called absorption-only model, that accounts for only primary photon absorption on the EBL and adiabatic losses. In fact, there exist several indications that the absorption-only model is incomplete. We present and discuss the intergalactic electromagnetic cascade model (IECM) --- the simplest model that allows to coherently explain all known anomalies. This model has a number of robust signatures that could be searched for with present and future instruments. The IECM model may serve as a new background template, allowing to make future searches for $\gamma \rightarrow$ ALP oscillation more robust. A detailed account of our calculations is available in astro-ph/1609.01013v2 (A\&A,in print).

Context. Most of the studies on extragalactic {\gamma}-ray propagation performed up to now only accounted for primary gamma-ray absorption and adiabatic losses ("absorption-only model"). However, there is growing evidence that this model is oversimplified and must be modified in some way. In particular, it was found that the intensity extrapolated from the optically-thin energy range of some blazar spectra is insufficient to explain the optically-thick part of these spectra. This effect was interpreted as an indication for {\gamma}-axion-like particle (ALP) oscillation. On the other hand, there are many hints that a secondary component from electromagnetic cascades initiated by primary {\gamma}-rays or nuclei may be observed in the spectra of some blazars. Aims. We study the impact of electromagnetic cascades from primary {\gamma}-rays or protons on the physical interpretation of blazar spectra obtained with imaging Cherenkov telescopes. Methods. We use the publicly-available code ELMAG to compute observable spectra of electromagnetic cascades from primary {\gamma}-rays. For the case of primary proton, we develop a simple, fast and reasonably accurate hybrid method to calculate the observable spectrum. We perform the fitting of the observed spectral energy distributions (SEDs) with various physical models: the absorption-only model, the "electromagnetic cascade model" (for the case of primary {\gamma}-rays), and several versions of the hadronic cascade model (for the case of primary proton). We distinguish the following species of hadronic cascade models: 1) "basic hadronic model", where it is assumed that the proton beam travels undisturbed by extragalactic magnetic field and that all observable {\gamma}-rays are produced by primary protons through photohadronic processes with subsequent development of electromagnetic cascades /abridged

Most of the recent research on extragalactic $\gamma$-ray propagation focused on the study of the $\gamma\gamma \rightarrow e^{+}e^{-}$ absorption process ("absorption-only model"). Starting from a possible anomaly at very high energies (VHE, E >100 GeV), we briefly review several existing deviations from this model. The exotic interpretation of the VHE anomaly is not supported by the recent works. On the other hand, the process of intergalactic electromagnetic cascade development naturally explains these effects. We discuss phenomenology of intergalactic cascades and the main spectral signatures of the electromagnetic cascade model. We also briefly consider the hadronic cascade model; it also may explain the data, but requires low strength of magnetic field around the source of primary protons or nuclei.

We briefly review contemporary extragalactic {\gamma}-ray propagation models. It is shown that the Extragalactic Magnetic Field (EGMF) strength and structure are poorly known. Strict lower limits on the EGMF strength in voids are of order 10^{-17}--10^{-20} G, thus allowing a substantial contribution of a secondary component generated by electromagnetic cascades to the observable spectrum. We show that this "electromagnetic cascade model" is supported by data from imaging Cherenkov telescopes and the Fermi LAT detector.

For some blazars, the gamma-ray absorption features due to pair-production on the Extragalactic Background Light (EBL) are fainter than expected. The present work reviews the main models that could explain this paradox, with emphasis on conservative ones, that do not include any new physics. The models that are intrinsic to the source, do allow a very hard primary spectrum, but fail to explain a regular redshift dependence of the anomaly starting energy. The model that includes a contribution from secondary photons produced by cosmic rays (CR) near the Earth seems to require a well collimated CR beam, what is hard to achieve. Finally, the model with secondary photons produced in electromagnetic (EM) cascades initiated by primary gamma-rays is considered. In principle, it allows to decrease the statistical significance of the anomaly and, while requiring quite low EGMF strength B, does not contradict to most contemporary constraints on the B value. Additionally, it is shown that the recently observed correlation between directions to hard gamma-ray sources and voids in the Large Scale Structure is a natural feature of the EM cascade model.

We present a very brief overview of some recent gamma-ray observations of selected blazars to reveal an indication for a considerable or even dominant contribution of secondary gamma-rays from electromagnetic cascades to the observable spectra in the 1-500 GeV energy range.

Nowadays there is compelling evidence for the existence of dark matter in the Universe. A general consensus has been expressed on the need for a directional sensitive detector to confirm, with a complementary approach, the candidates found in conventional searches and to finally extend their sensitivity beyond the limit of neutrino-induced background. We propose here the use of a detector based on nuclear emulsions to measure the direction of WIMP-induced nuclear recoils. The production of nuclear emulsion films with nanometric grains is established. Several measurement campaigns have demonstrated the capability of detecting sub-micrometric tracks left by low energy ions in such emulsion films. Innovative analysis technologies with fully automated optical microscopes have made it possible to achieve the track reconstruction for path lengths down to one hundred nanometers and there are good prospects to further exceed this limit. The detector concept we propose foresees the use of a bulk of nuclear emulsion films surrounded by a shield from environmental radioactivity, to be placed on an equatorial telescope in order to cancel out the effect of the Earth rotation, thus keeping the detector at a fixed orientation toward the expected direction of galactic WIMPs. We report the schedule and cost estimate for a one-kilogram mass pilot experiment, aiming at delivering the first results on the time scale of six years.

We present an event-by-event study of cosmic ray (CR) composition with the reflected Cherenkov light method. The fraction of CR light component above 5 PeV was reconstructed using the 2013 run data of the SPHERE experiment which observed optical Vavilov-Cherenkov radiation of extensive air showers, reflected from snow surface of Lake Baikal. Additionally, we discuss a possibility to improve the elemental groups separability by means of multidimensional criteria.

We give an overview of the SPHERE experiment based on detection of reflected Vavilov-Cherenkov radiation (Cherenkov light) from extensive air showers in the energy region E>10^{15} eV. A brief history of the reflected Cherenkov light technique is given; the observations carried out with the SPHERE-2 detector are summarized; the methods of the experimental datasample analysis are described. The first results on the primary cosmic ray all-nuclei energy spectrum and mass composition are presented. Finally, the prospects of the SPHERE experiment and the reflected Cherenkov light technique are given.

It is well known that the effect of gamma-ray absorption on extragalactic background light (EBL) is weakly expressed in the spectra of some blazars. It is shown that a secondary component generated by electromagnetic cascades might considerably decrease the statistical significance of this anomaly. Observational results indicate the existence of the cascade component in the spectra of extragalactic gamma-ray sources, thus supporting the proposed model.

The present paper highlights the data analysis status of the OPERA experiment. The experiment was designed to perform the neutrino interactions analysis on event-by-event basis, and optimized to search for $\nu_{\mu} \rightarrow \nu_{\tau}$ oscillation in appearance mode, also allowing to perform a $\nu_{e}$ appearance search. We review the data simulation and the analysis chains implemented to search for $\nu_{\tau}$ interactions. The main kinematical parameters sensitive to the neutrino flavor are discussed, the uncertainties of the event parameter estimation are reviewed, and the main sources of background for the $\nu_{\mu} \rightarrow \nu_{\tau}$ oscillation search are examined. The topologies of the two first $\nu_{\tau}$ candidate events are presented. Finally, we review the status of the $\nu_{e}$ appearance search and present the constraints set by the OPERA experiment on the mixing angle $\theta_{13}$ and on the LSND/MiniBooNE anomaly.

First preliminary results of the balloon-borne experiment SPHERE-2 on the all-nuclei primary cosmic rays (PCR) spectrum and primary composition are presented. The primary spectrum in the energy range $10^{16}$--$5\cdot10^{17}$ eV was reconstructed using characteristics of Vavilov-Cherenkov radiation of extensive air showers (EAS), reflected from a snow surface. Several sources of systematic uncertainties of the spectrum were analysed. A method for separation of the primary nuclei' groups based on the lateral distribution function' (LDF) steepness parameter is presented. Preliminary estimate of the mean light nuclei' fraction $f_{30-150}$ at energies $3\cdot10^{16}$--$1.5\cdot10^{17}$ eV was performed and yielded $f_{30-150}$= (21$\pm$11)%.

Here is presented the current state of the SPHERE-2 balloon-borne experiment. The detector is elevated up to 1 km above the snow surface and registers the reflected Vavilov-Cherenkov radiation from extensive air showers. This method has good sensitivity to the mass-composition of the primary cosmic rays due to its high resolution near the shower axis. The detector consists of a 1500 mm spherical mirror with a 109 PMT cluster in its focus. The electronics record a signal pulse profile in each PMT. In the last 2 years the detector was upgraded: time resolution of pulse registration was enhanced up to 12.5 ns, channel sensitivity was increased by a factor of 3, a new LED-based relative PMT calibration method was introduced, and new hardware and etc. was installed.