Evidence for a large-scale supergalactic cosmic ray multiplet (arrival directions correlated with energy) structure is reported for ultra-high energy cosmic ray (UHECR) energies above 10$^{19}$ eV using seven years of data from the Telescope Array (TA) surface detector and updated to 10 years. Previous energy-position correlation studies have made assumptions regarding magnetic field shapes and strength, and UHECR composition. Here the assumption tested is that, since the supergalactic plane is a fit to the average matter density of the local Large Scale Structure (LSS), UHECR sources and intervening extragalactic magnetic fields are correlated with this plane. This supergalactic deflection hypothesis is tested by the entire field-of-view (FOV) behavior of the strength of intermediate-scale energy-angle correlations. These multiplets are measured in spherical cap section bins (wedges) of the FOV to account for coherent and random magnetic fields. The structure found is consistent with supergalactic deflection, the previously published energy spectrum anisotropy results of TA (the hotspot and coldspot), and toy-model simulations of a supergalactic magnetic sheet. The seven year data post-trial significance of this supergalactic structure of multiplets appearing by chance, on an isotropic sky, is found by Monte Carlo simulation to be 4.2$\sigma$. The ten years of data post-trial significance is 4.1$\sigma$. Furthermore, the starburst galaxy M82 is shown to be a possible source of the TA Hotspot, and an estimate of the supergalactic magnetic field using UHECR measurements is presented.

Motivated by the detection of a significant dipole structure in the arrival directions of ultrahigh-energy cosmic rays above 8 EeV reported by the Pierre Auger Observatory (Auger), we search for a large-scale anisotropy using data collected with the surface detector array of the Telescope Array Experiment (TA). With 11 years of TA data, a dipole structure in a projection of the right ascension is fitted with an amplitude of 3.3+- 1.9% and a phase of 131 +- 33 degrees. The corresponding 99% confidence-level upper limit on the amplitude is 7.3%. At the current level of statistics, the fitted result is compatible with both an isotropic distribution and the dipole structure reported by Auger.

Ultra high energy cosmic rays provide the highest known energy source in the universe to measure proton cross sections. Though conditions for collecting such data are less controlled than an accelerator environment, current generation cosmic ray observatories have large enough exposures to collect significant statistics for a reliable measurement for energies above what can be attained in the lab. Cosmic ray measurements of cross section use atmospheric calorimetry to measure depth of air shower maximum ($X_{\mathrm{max}}$), which is related to the primary particle's energy and mass. The tail of the $X_{\mathrm{max}}$ distribution is assumed to be dominated by showers generated by protons, allowing measurement of the inelastic proton-air cross section. In this work the proton-air inelastic cross section measurement, $\sigma^{\mathrm{inel}}_{\mathrm{p-air}}$, using data observed by Telescope Array's Black Rock Mesa and Long Ridge fluorescence detectors and surface detector array in hybrid mode is presented. $\sigma^{\mathrm{inel}}_{\mathrm{p-air}}$ is observed to be $520.1 \pm 35.8$[Stat.] $^{+25.0}_{-40}$[Sys.]~mb at $\sqrt{s} = 73$ TeV. The total proton-proton cross section is subsequently inferred from Glauber formalism and is found to be $\sigma^{\mathrm{tot}}_{\mathrm{pp}} = 139.4 ^{+23.4}_{-21.3}$ [Stat.]$ ^{+15.0}_{-24.0}$[Sys.]~mb.

We present an upper limit on the flux of ultra-high-energy down-going neutrinos for $E > 10^{18}\ \mbox{eV}$ derived with the nine years of data collected by the Telescope Array surface detector (05-11-2008 -- 05-10-2017). The method is based on the multivariate analysis technique, so-called Boosted Decision Trees (BDT). Proton-neutrino classifier is built upon 16 observables related to both the properties of the shower front and the lateral distribution function.

The surface detector (SD) of the Telescope Array (TA) experiment allows one to indirectly detect photons with energies of order $10^{18}$ eV and higher and to separate photons from the cosmic-ray background. In this paper we present the results of a blind search for point sources of ultra-high energy (UHE) photons in the Northern sky using the TA SD data. The photon-induced extensive air showers (EAS) are separated from the hadron-induced EAS background by means of a multivariate classifier based upon 16 parameters that characterize the air shower events. No significant evidence for the photon point sources is found. The upper limits are set on the flux of photons from each particular direction in the sky within the TA field of view, according to the experiment's angular resolution for photons. Average 95% C.L. upper limits for the point-source flux of photons with energies greater than $10^{18}$, $10^{18.5}$, $10^{19}$, $10^{19.5}$ and $10^{20}$ eV are $0.094$, $0.029$, $0.010$, $0.0073$ and $0.0058$ km$^{-2}$yr$^{-1}$ respectively. For the energies higher than $10^{18.5}$ eV, the photon point-source limits are set for the first time. Numerical results for each given direction in each energy range are provided as a supplement to this paper.

We present the results of the search for ultra-high-energy photons with nine years of data from the Telescope Array surface detector. A multivariate classifier is built upon 16 reconstructed parameters of the extensive air shower. These parameters are related to the curvature and the width of the shower front, the steepness of the lateral distribution function, and the timing parameters of the waveforms sensitive to the shower muon content. A total number of two photon candidates found in the search is fully compatible with the expected background. The $95\%\,$CL limits on the diffuse flux of the photons with energies greater than $10^{18.0}$, $10^{18.5}$, $10^{19.0}$, $10^{19.5}$ and $10^{20.0}$ eV are set at the level of $0.067$, $0.012$, $0.0036$, $0.0013$, $0.0013~\mbox{km}^{-2}\mbox{yr}^{-1}\mbox{sr}^{-1}$ correspondingly.

The paper is a script of a lecture given at the ISAPP-Baikal summer school in 2018. The lecture gives an overview of the Tunka Advanced Instrument for cosmic rays and Gamma Astronomy (TAIGA) facility including historical introduction, description of existing and future setups, and outreach and open data activities.

The results on ultra-high-energy cosmic rays (UHECR) mass composition obtained with the Telescope Array surface detector are presented. The analysis employs the boosted decision tree (BDT) multivariate analysis built upon 14 observables related to both the properties of the shower front and the lateral distribution function. The multivariate classifier is trained with Monte-Carlo sets of events induced by the primary protons and iron. An average atomic mass of UHECR is presented for energies $10^{18.0}-10^{20.0}\ \mbox{eV}$. The average atomic mass of primary particles shows no significant energy dependence and corresponds to $\langle \ln A \rangle = 2.0 \pm 0.1 (stat.) \pm 0.44 (syst.)$. The result is compared to the mass composition obtained by the Telescope Array with $\mbox{X}_{\mbox{max}}$ technique along with the results of other experiments. Possible systematic errors of the method are discussed.

The Pierre Auger Collaboration (Auger) recently reported a correlation between the arrival directions of cosmic rays with energies above 39 EeV and the flux pattern of 23 nearby starburst galaxies (SBGs). In this Letter, we tested the same hypothesis using cosmic rays detected by the Telescope Array experiment (TA) in the 9-year period from May 2008 to May 2017. Unlike the Auger analysis, we did not optimize the parameter values but kept them fixed to the best-fit values found by Auger, namely 9.7% for the anisotropic fraction of cosmic rays assumed to originate from the SBGs in the list and 12.9{\deg} for the angular scale of the correlations. The energy threshold we adopted is 43 EeV, corresponding to 39 EeV in Auger when taking into account the energy-scale difference between two experiments. We find that the TA data is compatible with isotropy to within 1.1{\sigma} and with the Auger result to within 1.4{\sigma}, meaning that it is not capable to discriminate between these two hypotheses.

The Telescope Array observatory utilizes fluorescence detectors and surface detectors to observe air showers produced by ultra high energy cosmic rays in the Earth's atmosphere. Cosmic ray events observed in this way are termed hybrid data. The depth of air shower maximum is related to the mass of the primary particle that generates the shower. This paper reports on shower maxima data collected over 8.5 years using the Black Rock Mesa and Long Ridge fluorescence detectors in conjunction with the array of surface detectors. We compare the means and standard deviations of the observed $X_{\mathrm{max}}$ distributions with Monte Carlo $X_{\mathrm{max}}$ distributions of unmixed protons, helium, nitrogen, and iron, all generated using the QGSJet~II-04 hadronic model. We also perform an unbinned maximum likelihood test of the observed data, which is subjected to variable systematic shifting of the data $X_{\mathrm{max}}$ distributions to allow us to test the full distributions, and compare them to the Monte Carlo to see which elements are not compatible with the observed data. For all energy bins, QGSJet~II-04 protons are found to be compatible with Telescope Array hybrid data at the 95% confidence level after some systematic $X_{\mathrm{max}}$ shifting of the data. Three other QGSJet~II-04 elements are found to be compatible using the same test procedure in an energy range limited to the highest energies where data statistics are sparse.

One of the uncertainties in interpretation of ultra-high energy cosmic ray (UHECR) data comes from the hadronic interaction models used for air shower Monte Carlo (MC) simulations. The number of muons observed at the ground from UHECR-induced air showers is expected to depend upon the hadronic interaction model. One may therefore test the hadronic interaction models by comparing the measured number of muons with the MC prediction. In this paper, we present the results of studies of muon densities in UHE extensive air showers obtained by analyzing the signal of surface detector stations which should have high $\it{muon \, purity}$. The muon purity of a station will depend on both the inclination of the shower and the relative position of the station. In 7 years' data from the Telescope Array experiment, we find that the number of particles observed for signals with an expected muon purity of $\sim$65% at a lateral distance of 2000 m from the shower core is $1.72 \pm 0.10{\rm (stat.)} \pm 0.37 {\rm (syst.)}$ times larger than the MC prediction value using the QGSJET II-03 model for proton-induced showers. A similar effect is also seen in comparisons with other hadronic models such as QGSJET II-04, which shows a $1.67 \pm 0.10 \pm 0.36$ excess. We also studied the dependence of these excesses on lateral distances and found a slower decrease of the lateral distribution of muons in the data as compared to the MC, causing larger discrepancy at larger lateral distances.

We report on a measurement of the cosmic ray energy spectrum by the Telescope Array Low-Energy Extension (TALE) air fluorescence detector. The TALE air fluorescence detector is also sensitive to the Cherenkov light produced by shower particles. Low energy cosmic rays, in the PeV energy range, are detectable by TALE as "Cherenkov Events". Using these events, we measure the energy spectrum from a low energy of $\sim 2$ PeV to an energy greater than 100 PeV. Above 100 PeV TALE can detect cosmic rays using air fluorescence. This allows for the extension of the measurement to energies greater than a few EeV. In this paper, we will describe the detector, explain the technique, and present results from a measurement of the spectrum using $\sim 1000$ hours of observation. The observed spectrum shows a clear steepening near $10^{17.1}$ eV, along with an ankle-like structure at $10^{16.2}$ eV. These features present important constraints on galactic cosmic rays origin and propagation models. The feature at $10^{17.1}$ eV may also mark the end of the galactic cosmic rays flux and the start of the transition to extra-galactic sources.

An intermediate-scale energy spectrum anisotropy has been found in the arrival directions of ultra-high energy cosmic rays of energies above $10^{19.2}$ eV in the northern hemisphere, using 7 years of data from the Telescope Array surface detector. A relative energy distribution test is done comparing events inside oversampled spherical caps of equal exposure, to those outside, using the Poisson likelihood ratio. The center of maximum significance is at $9^h$$16^m$, $45^{\circ}$. and has a deficit of events with energies $10^{19.2}$$\leq$$E$$<$$10^{19.75}$ eV and an excess for $E$$\geq$$10^{19.75}$ eV. The post-trial probability of this energy anisotropy, appearing by chance anywhere on an isotropic sky, is found by Monte Carlo simulation to be $9$$\times$$10^{-5}$ ($3.74$$\sigma_{global}$).

The energy of the ultrahigh energy spectral cutoff was measured, integrating over the northern hemisphere sky, by the Telescope Array (TA) collaboration, to be $10^{19.78 \pm 0.06}$ eV, in agreement with the High Resolution Fly's Eye (HiRes) experiment, whereas the Pierre Auger experiment, integrating over the southern hemisphere sky, measured the cutoff to be at 10$^{19.62 \pm 0.02}$ eV. An 11% energy scale difference between the TA and Auger does not account for this difference. However, in comparing the spectra of the Telescope Array and Pierre Auger experiments in the band of declination common to both experiments ($-15.7^{\circ} < \delta < 24.8^{\circ}$) we have found agreement in the energy of the spectral cutoff. While the Auger result is essentially unchanged, the TA cutoff energy has changed to $10^{19.59 \pm 0.06}$ eV. In this paper we argue that this is an astrophysical effect.

The Telescope Array (TA) experiment is located in the western desert of Utah, USA, and observes ultra high energy cosmic rays (UHECRs) in the Northern hemisphere. At the highest energies, $E>10$~EeV, the shape of cosmic ray energy spectrum may carry an imprint of the source density distribution along the line of sight different in different directions of the sky. In this study, we search for such directional variations in the shape of the energy spectrum using events observed with the Telescope Array's surface detector. We divide the TA field of view into two nearly equal-exposure regions: the "on-source" region which we define as $\pm 30^\circ$ of the supergalactic plane containing mostly nearby structures, and the complementary "off-source" region where the sources are further away on average. We compare the UHECR spectra in these regions by fitting them to the broken power law and comparing the resulting parameters. We find that the off-source spectrum has an earlier break at highest energies. The chance probability to obtain such or larger difference in statistically equivalent distributions is estimated as $6.2\pm1.1\times10^{-4}$ ($3.2\sigma$) by a Monte-Carlo simulation. The observed difference in spectra is in a reasonable quantitative agreement with a simplified model that assumes that the UHECR sources trace the galaxy distribution from the 2MRS catalogue, primary particles are protons and the magnetic deflections can be neglected.

Cosmic rays in the energy range $10^{18.0}$ - $10^{18.5}$ eV are thought to have a light, probably protonic, composition. To study their origin one can search for anisotropy in their arrival directions. Extragalactic cosmic rays should be isotropic, but galactic cosmic rays of this type should be seen mostly along the galactic plane, and there should be a shortage of events coming from directions near the galactic anticenter. This is due to the fact that, under the influence of the galactic magnetic field, the transition from ballistic to diffusive behavior is well advanced, and this qualitative picture persists over the whole energy range. Guided by models of the galactic magnetic field that indicate that the enhancement along the galactic plane should have a standard deviation of about 20$^\circ$ in galactic latitude, and the deficit in the galactic anticenter direction should have a standard deviation of about 50$^\circ$ in galactic longitude, we use the data of the Telescope Array surface detector in $10^{18.0}$ to $10^{18.5}$ eV energy range to search for these effects. The data are isotropic. Neither an enhancement along the galactic plane nor a deficit in the galactic anticenter direction is found. Using these data we place an upper limit on the fraction of EeV cosmic rays of galactic origin at 1.3% at 95% confidence level.

We propose the liquid-scintillator detector LENA (Low Energy Neutrino Astronomy) as a next-generation neutrino observatory on the scale of 50 kt. The outstanding successes of the Borexino and KamLAND experiments demonstrate the large potential of liquid-scintillator detectors in low-energy neutrino physics. LENA's physics objectives comprise the observation of astrophysical and terrestrial neutrino sources as well as the investigation of neutrino oscillations. In the GeV energy range, the search for proton decay and long-baseline neutrino oscillation experiments complement the low-energy program. Based on the considerable expertise present in European and international research groups, the technical design is sufficiently mature to allow for an early start of detector realization.

A new EAS Cherenkov light array, Tunka-133, with ~1 km^2 geometrical area has been installed at the Tunka Valley (50 km from Lake Baikal) in 2009. The array permits a detailed study of cosmic ray energy spectrum and mass composition in the energy range 10^16 - 10^18 eV with a uniform method. We describe the array construction, DAQ and methods of the array calibration.The method of energy reconstruction and absolute calibration of measurements are discussed. The analysis of spatial and time structure of EAS Cherenkov light allows to estimate the depth of the EAS maximum X_max. The results on the all particles energy spectrum and the mean depth of the EAS maximum X_max vs. primary energy derived from the data of two winter seasons (2009 -- 2011), are presented. Preliminary results of joint operation of the Cherenkov array with antennas for detection of EAS radio signals are shown. Plans for future upgrades -- deployment of remote clusters, radioantennas and a scintillator detector network and a prototype of the HiSCORE gamma-telescope -- are discussed.

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 new EAS Cherenkov array Tunka-133 with about 1 km**2 geometric acceptance area is installed in the Tunka Valley (50 km from Lake Baikal). The array willpermit a detailed study of cosmic ray energy spectrum and mass composition in the energy range of 10**15 - 10**18 eV with uniform method. The array consistsof 19 clusters, each composed of 7 optical detectors with 20 cm PMTs. Since November 2008, the first 12 clusters are in operation, commissioning of the whole array is planned for September 2009 (At the time of submission of this paperto electronic arXiv(February 2010) the completed Tunka-133 array is already taking data). We describe the array construction and DAQ, preliminary results and plans for the future development: deployment of radio-antennas and muon detectors network.

The final analysis of the Extensive Air Shower (EAS) maximum X_max depth distribution derived from the data of Tunka-25 atmospheric Cherenkov light array in the energy range 3.10^15 - 3.10^16 eV is presented. The perspectives of X_max studies with the new Cherenkov light array Tunka-133 of 1 km^2 area, extending the measurements up to 10^18 eV, are discussed.

The new EAS Cherenkov array TUNKA-133, with about 1 km**2 sensitive area, is being installed in the Tunka Valley. The investigated energy range is 10**15-10**18 eV. It will consist of 133 optical detectors based on EMI9350 PMTs. Optical detectors are grouped into 19 clusters with 7 detectors each. The detectors are connected to the cluster box with RG-58 cables. Every PMT signal is digitized in the cluster box with 200 MHz FADC. The cluster boxes are connected to the data acquisition center with a 1 Gb/s optical link. A detailed description of the data acquisition system (DAQ) is presented.

The new EAS Cherenkov array Tunka-133, with about 1 km**2 sensitive area, is being installed in the Tunka Valley since the end of 2005. This array will permit a detailed study of the cosmic ray energy spectrum and the mass composition in the energy range of 10**15-10**18 eV with a unique method. The array will consist of 19 clusters, each composed of 7 optical detectors. The first cluster started operation in October 2006. We describe the data acquisition system and present preliminary results from data taken with the first cluster.

In this paper we present the results of extensive studies of scintillators for hybrid phototubes with luminescent screens. The results of the developments of such phototubes with a variety of scintillators are presented. New scintillator materials for such kind of application are discussed. The requirements for scintillators to use in such hybrid phototubes are formulated. It is shown that very fast and highly efficient inorganic scintillators like ZnO:Ga will be ideal scintillators for such kind of application.

In this article we describe results of a photoelectron backscattering effect in vacuum phototubes: classical photomultipliers (PMT) and hybrid phototubes (PH). Late pulses occurring in PMTs are attributed to the photoelectron backscattering and distinguished from pulses due to an anode glow effect. The late pulses are measured in a number of PMTs and HPs with various photocathode sizes covering 1-50 cm range and different types of the first dynode materials and construction designs. It is shown that the late pulses are a generic feature of all vacuum photodetectors - PMTs and PHs and they don't deteriorate dramatically amplitude and timing responses of vacuum phototubes.

In this very short note we review some historical aspects of photomultiplier tube invention. It is our tribute to the memory of great Soviet-Russian physicist and engineer Leonid Aleksandrovitch Kubetsky whose life and scientific achievements are described briefly. Particular efforts are made to shed light on a controversial issue of who invented the first photomultiplier tube. It is asserted that if to recognize L.A.Kubetsky's priority on the photomultiplier tube invention the last Beaune Conference would be held on the eve of the 75th Anniversary of that great event.

We present results of an improved analysis of the experimental data of the EAS Cherenkov array Tunka-25. A new function to fit the Cherenkov light lateral distribution LDF at core distances from 0 to 350 m has been developed on the base of CORSIKA simulations and applied to the analysis of Tunka data. Two methods to estimate the EAS maximum position have been used. The one is based on the pulse FWHM, the other on the light LDF. We present the primary energy spectrum in the energy range 10^15 - 10^17 eV. The use of the depth of the EAS maximum to determine the mean mass composition is discussed.

The project of an EAS Cherenkov array in the Tunka valley/Siberia with an area of about 1 km2 is presented. The new array will have a ten times bigger area than the existing Tunka-25 array and will permit a detailed study of the cosmic ray energy spectrum and the mass composition in the energy range from 10^15 to 10^18 eV.

I present here some reflections and very speculative remarks on the detection of relativistic magnetic monopoles by currently operating deep underwater/ice neutrino telescopes.

A new primary mass independent method of energy measurement has been developed by exploiting: a) the joint analysis of the shower size, obtained by EAS-TOP, and of the EAS atmospheric Cherenkov light lateral distribution, obtained by the QUEST array, and b) simulations based on the CORSIKA/QGSJET code. The method is based on the correlation between the size/energy ratio and the steepness of Cherenkov light lateral distribution and has been compared with a "classical" one based on the Cherenkov light flux at a fixed distance (175 m) from the EAS core. An absolute energy calibration of the EAS atmospheric Cherenkov light flux has been obtained.

The project of an EAS Cherenkov array in the Tunka valley/Siberia with an area of about 1 km^2 is presented. The new array will have a ten times bigger area than the existing Tunka-25 array and will permit a detailed study of the cosmic ray energy spectrum and the mass composition in the energy range from 10^15 to 10^18 eV.

New results of 300 hours of operation of the Tunka array are presented. An improved parametrization of the Cherenkov light lateral distribution function (LDF), based on CORSIKA Monte Carlo simulations and the experiment QUEST, has been used for the reconstruction of EAS parameters. The corrected energy spectrum in the knee region is obtained. The mean depth of the EAS maximum has been derived both from the analysis of LDF steepness and the FWHM of Cerenkov light pulse. The mean mass composition around the knee is estimated.

We present the results of studies of a nanosecond light pulser built following J.S.Kapustinsky et al original design and using bright InGaN/GaN ultraviolet and blue LEDs produced by NICHIA CHEMICAL. It is shown how timing characteristics of the pulser depend on the type of LED and the value of power supply voltage.

We present the results of stability studies of nanosecond light sources based on single quantum well (SQW) InGaN/GaN ultra bright blue LEDs. It is shown that the light yield of such light sources and their timing characteristics don't deteriorate even after 10^10 total pulses. The longterm stability of the sources light yield is better than 1%.

The results of direct measurements of group velocity of light in the lake Baikal water at the depth of 1100 m are presented. The lake Baikal water dispersion has been measured at three wavelengths: 370 nm, 470 nm and 525 nm. The results are in a rather good agreement with theoretical predictions.