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.

The extremely low flux of ultra-high energy cosmic rays (UHECR) makes their direct observation by orbital experiments practically impossible. For this reason all current and planned UHECR experiments detect cosmic rays indirectly by observing the extensive air showers (EAS) initiated by cosmic ray particles in the atmosphere. The world largest statistics of the ultra-high energy EAS events is recorded by the networks of surface stations. In this paper we consider a novel approach for reconstruction of the arrival direction of the primary particle based on the deep convolutional neural network. The latter is using raw time-resolved signals of the set of the adjacent trigger stations as an input. The Telescope Array Surface Detector is an array of 507 stations, each containing two layers plastic scintillator with an area of $3$ m$^2$. The training of the model is performed with the Monte-Carlo dataset. It is shown that within the Monte-Carlo simulations, the new approach yields better resolution than the traditional reconstruction method based on the fitting of the EAS front. The details of the network architecture and its optimisation for this particular task are discussed.

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 consider the ensemble of very-high-energy gamma-ray sources observed at distances and energies where a significant absorption of gamma rays is expected due to pair production on the extragalactic background light (EBL). Previous studies indicated that spectra of these sources, upon correction for the absorption, exhibit unusual spectral hardenings which happen precisely at the energies where the correction becomes significant. Here, we address this subject with the most recent clean sample of distant gamma-ray blazars, making use of published results of imaging atmospheric Cerenkov telescopes and of the Fermi-LAT Pass 8 data, supplemented by the newest absorption models and individual measurements of sources' redshifts. We perform a search for spectral breaks at energies corresponding to unit optical depth with respect to the absorption on EBL. These energies are different for distant and nearby objects, and consequently, such features may not be related to intrinsic properties of the sources. While in some spectra such breaks are not seen, hardenings at distance-dependent energies are present in many of them, though the overall statistical significance of the effect is lower than reported in previous studies. The dependence of the break strength on the redshift found earlier is not confirmed in the new analysis.

Digital radio arrays are widely used for the low-frequency radio astronomy as well as for detection of air-showers induced by high-energy cosmic rays and neutrinos. Since the radio emission from air-showers forms short broadband pulses with duration of tens nanoseconds, the data acquisition strategies of cosmic-ray and astronomical arrays have significant differences. To perform precise measurement of cosmic rays, the radio array should have absolute amplitude calibration and record the entire electric field on the antenna in the broad frequency range. These requirements are similar to ones defined for the experiments aimed at the detection of weak signal from neutral hydrogen at redshifts of $z$>10, what led us to the application of our experience with Tunka-Rex to this problem. We are developing new experimental setup comprising of four antenna stations, placed on the area of 100 sq.m. Each antenna station consists of two perpendicular loop antennas measuring electric field in the frequency band of 30-80 MHz. The setup records electric fields from all antennas in portions of 50 $\mu$s reaching the spectral resolution of 20 kHz. We expect a flow of redundant data of about 10 GB/day, and plan to exploit this redundancy in order to decrease systematic uncertainty of the measurements by application of digital beam-forming, matched filtering and RFI suppression with neural networks. In the present contribution we describe the design and calibration of the setup, expected performance and data analysis techniques.

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.

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.

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.

It was proposed that the two phenomena, WMAP-Planck haze and Fermi bubbles, may have a common origin. In the present paper we analyze the spatial structure of the haze using the Planck 2018 data release. It is found that the spatial dimensions and locations of WMAP-Planck haze and Fermi bubbles are compatible within the experimental uncertainties. No substructures similar to the Fermi bubbles cocoon are identified in the Planck data. Comparison with the spatial extent of possible synchrotron emission caused by the electron-positron pair emitted by the Galactic center pulsar population and by the decay of dark matter particles in the Galactic center region are performed. Both galactic pulsars and dark matter decay remain viable explanations of the WMAP-Planck haze.

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.

There are multiple techniques to determine the chemical composition of the ultra-high-energy cosmic rays. While most of the methods are primarily sensitive to the average atomic mass, it is challenging to discriminate between the two lightest elements: proton and helium. In this paper, the proton-to-helium ratio in the energy range $10^{18.0} \mbox{eV}$ to $10^{19.3} \mbox{eV}$ is estimated using the tail of the distribution of the depth of the shower maximum $X_{\mathrm{max}}$. Using the exponential decay scale $\Lambda$ measured by the Pierre Auger Observatory and the Telescope Array experiment we derive the 68\%\,CL constraints on the proton-to-helium ratio $\mathrm{p}/\mathrm{He} > 7.3 $ and $\mathrm{p}/\mathrm{He} > 0.43 $ for $10^{18.0} < E < 10^{18.5}$ eV and $10^{18.3} < E < 10^{19.3}$ eV correspondingly. It is shown that the result is conservative with respect to the admixture of heavier elements. We evaluate the impact of the hadronic interaction model uncertainty. The implications for the astrophysical models of the origin of cosmic rays and the safety of the future colliders are discussed.

The Sydney University Giant Air-shower Recorder (SUGAR) measured the energy spectrum of ultra-high-energy cosmic rays reconstructed from muon-detector readings, while the Pierre Auger Observatory, looking at the same Southern sky, uses the calorimetric fluorescence method for the same purpose. Comparison of their two spectra allows us to reconstruct the empirical dependence of the number of muons in the shower on the primary energy for energies between $10^{17}$ and $10^{18.5}$ eV. We compare this dependence with the predictions of hadronic interaction models \mbox{QGSJET-II-04} and \mbox{EPOS-LHC}. The empirically determined number of muons with energies above 0.75 GeV exceeds the simulated one by the factors $\sim$1.67 and $\sim$1.28 for $10^{17}$ eV proton and iron primaries, respectively. The muon excess grows moderately with the primary energy, increasing by an additional factor of $\sim 1.2$ for $10^{18.5}$ eV primaries.

Bursts of gamma ray showers have been observed in coincidence with downward propagating negative leaders in lightning flashes by the Telescope Array Surface Detector (TASD). The TASD is a 700~square kilometer cosmic ray observatory located in southwestern Utah, U.S.A. In data collected between 2014 and 2016, correlated observations showing the structure and temporal development of three shower-producing flashes were obtained with a 3D lightning mapping array, and electric field change measurements were obtained for an additional seven flashes, in both cases co-located with the TASD. National Lightning Detection Network (NLDN) information was also used throughout. The showers arrived in a sequence of 2--5 short-duration ($\le$10~$\mu$s) bursts over time intervals of several hundred microseconds, and originated at an altitude of $\simeq$3--5 kilometers above ground level during the first 1--2 ms of downward negative leader breakdown at the beginning of cloud-to-ground lightning flashes. The shower footprints, associated waveforms and the effect of atmospheric propagation indicate that the showers consist primarily of downward-beamed gamma radiation. This has been supported by GEANT simulation studies, which indicate primary source fluxes of $\simeq$$10^{12}$--$10^{14}$ photons for $16^{\circ}$ half-angle beams. We conclude that the showers are terrestrial gamma-ray flashes (TGFs), similar to those observed by satellites, but that the ground-based observations are more representative of the temporal source activity and are also more sensitive than satellite observations, which detect only the most powerful TGFs.

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.

The present day spectrum of the extragalactic background light (EBL) in UV, optical and IR wavelengths is the integral result of multiple astrophysical processes going on throughout the evolution of the Universe. The relevant processes include star formation, stellar evolution, light absorption and emission by the cosmic dust. The properties of these processes are known with uncertainties which contribute to the EBL spectrum precision. In the present paper we develop a numerical model of the EBL spectrum while maintaining the explicit dependence on the astrophysical parameters involved. We constructed a Markov Chain in the parameter space by using the likelihood function built with the up-to-date upper and lower bounds on the EBL intensity. The posterior distributions built with the Markov Chain Monte Carlo method are used to determine an allowed range of the individual parameters of the model. Consequently, the star formation rate multiplication factor is constrained in the range $1.01 < C_{\mbox{sfr}} < 1.69$ at $68\%$ C.L. The method also results in the bounds on the lifetime, radius, dust particle density and opacity of the molecular clouds that have large ambiguity otherwise. It is shown that there is a reasonable agreement between the model and the intensity bounds while the astrophysical parameters of the best fit model are close to their estimates from literature.

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.

We discuss the physics case for and the concept of a medium-scale axion helioscope with sensitivities in the axion-photon coupling a few times better than CERN Axion Solar Telescope (CAST). Search for an axion-like particle with these couplings is motivated by several persistent astrophysical anomalies. We present early conceptual design, existing infrastructure, projected sensitivity and timeline of such a helioscope (Troitsk Axion Solar Telescope Experiment, TASTE) to be constructed in the Institute for Nuclear Research, Troitsk, Russia. The proposed instrument may be also used for the search of dark-matter halo axions.

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.

Some discrepancies have been reported between observed and simulated muon content of extensive air showers: the number of observed muons exceeded the expectations in HiRes-MIA, Yakutsk and Pierre Auger Observatory data. Here, we analyze the data of the Moscow State University Extensive Air Shower (EAS-MSU) array on E_mu>~10 GeV muons in showers caused by ~(10^17-10^18) eV primary particles and demonstrate that they agree with simulations (QGSJET-II-04 hadronic interaction model) once the primary composition inferred from the surface-detector data is assumed.

Several early Universe scenarios predict a direction-dependent spectrum of primordial curvature perturbations. This translates into the violation of the statistical isotropy of cosmic microwave background radiation. Previous searches for statistical anisotropy mainly focussed on a quadrupolar direction-dependence characterised by a single multipole vector and an overall amplitude $g_*$. Generically, however, the quadrupole has a more complicated geometry described by two multipole vectors and $g_*$. This is the subject of the present work. In particular, we limit the amplitude $g_*$ for different shapes of the quadrupole by making use of Planck 2015 maps. We also constrain certain inflationary scenarios which predict this kind of more general quadrupolar statistical anisotropy.

Results of the search for $\sim (10^{16} - 10^{17.5})$ eV primary cosmic-ray photons with the data of the Moscow State University (MSU) Extensive Air Shower (EAS) array are reported. The full-scale reanalysis of the data with modern simulations of the installation does not confirm previous indications of the excess of gamma-ray candidate events. Upper limits on the corresponding gamma-ray flux are presented. The limits are the most stringent published ones at energies $\sim 10^{17}$ eV.

The Tunka Radio Extension (Tunka-Rex) is a radio detector at the TAIGA facility located in Siberia nearby the southern tip of Lake Baikal. Tunka-Rex measures air-showers induced by high-energy cosmic rays, in particular, the lateral distribution of the radio pulses. The depth of the air-shower maximum, which statistically depends on the mass of the primary particle, is determined from the slope of the lateral distribution function (LDF). Using a model-independent approach, we have studied possible features of the one-dimensional slope method and tried to find improvements for the reconstruction of primary mass. To study the systematic uncertainties given by different primary particles, we have performed simulations using the CONEX and CoREAS software packages of the recently released CORSIKA v7.5 including the modern high-energy hadronic models QGSJet-II.04 and EPOS-LHC. The simulations have shown that the largest systematic uncertainty in the energy deposit is due to the unknown primary particle. Finally, we studied the relation between the polarization and the asymmetry of the LDF.

The Tunka-Radio extension (Tunka-Rex) is a radio detector for air showers in Siberia. From 2012 to 2014, Tunka-Rex operated exclusively together with its host experiment, the air-Cherenkov array Tunka-133, which provided trigger, data acquisition, and an independent air-shower reconstruction. It was shown that the air-shower energy can be reconstructed by Tunka-Rex with a precision of 15\% for events with signal in at least 3 antennas, using the radio amplitude at a distance of 120\,m from the shower axis as an energy estimator. Using the reconstruction from the host experiment Tunka-133 for the air-shower geometry (shower core and direction), the energy estimator can in principle already be obtained with measurements from a single antenna, close to the reference distance. We present a method for event selection and energy reconstruction, requiring only one antenna, and achieving a precision of about 20\%. This method increases the effective detector area and lowers thresholds for zenith angle and energy, resulting in three times more events than in the standard reconstruction.

The Tunka-Rex experiment (Tunka Radio Extension) has been deployed in 2012 at the Tunka Valley (Republic of Buryatia, Russia). Its purpose is to investigate methods for the energy spectrum and the mass composition of high-energy cosmic rays based on the radio emission of air showers. Tunka-Rex is an array of 25 radio antennas distributed over an area of 3 km^2. The most important feature of the Tunka-Rex is that the air-shower radio emission is measured in coincidence with the Tunka-133 installation, which detects the Cherenkov radiation generated by the same atmospheric showers. Joint measurements of the radio emission and the Cherenkov light provide a unique opportunity for cross calibration of both calorimetric detection methods. The main goal of Tunka-Rex is to determine the precision for the reconstruction of air-shower parameters using the radio detection technique. In this article we present the current status of Tunka-Rex and first results, including reconstruction methods for parameters of the primary cosmic rays.

Tunka-Rex, the Tunka Radio extension at the TAIGA facility (Tunka Advanced Instrument for cosmic ray physics and Gamma Astronomy) in Siberia, has recently been expanded to a total number 63 SALLA antennas, most of them distributed on an area of one square kilometer. In the first years of operation, Tunka-Rex was solely triggered by the co-located air-Cherenkov array Tunka-133. The correlation of the measurements by both detectors has provided direct experimental proof that radio arrays can measure the position of the shower maximum. The precision achieved so far is 40 g/cm^2, and several methodical improvements are under study. Moreover, the cross-comparison of Tunka-Rex and Tunka-133 shows that the energy reconstruction of Tunka-Rex is precise to 15 %, with a total accuracy of 20 % including the absolute energy scale. By using exactly the same calibration source for Tunka-Rex and LOPES, the energy scale of their host experiments, Tunka-133 and KASCADE-Grande, respectively, can be compared even more accurately with a remaining uncertainty of about 10 %. The main goal of Tunka-Rex for the next years is a study of the cosmic-ray mass composition in the energy range above 100 PeV: For this purpose, Tunka-Rex now is triggered also during daytime by the particle detector array Tunka-Grande featuring surface and underground scintillators for electron and muon detection.

The Tunka Radio Extension (Tunka-Rex) is an antenna array consisting of 63 antennas at the location of the TAIGA facility (Tunka Advanced Instrument for cosmic ray physics and Gamma Astronomy) in Eastern Siberia, nearby Lake Baikal. Tunka-Rex is triggered by the air-Cherenkov array Tunka-133 during clear and moonless winter nights and by the scintillator array Tunka-Grande during the remaining time. Tunka-Rex measures the radio emission from the same air-showers as Tunka-133 and Tunka-Grande, but with a higher threshold of about 100 PeV. During the first stages of its operation, Tunka-Rex has proven, that sparse radio arrays can measure air-showers with an energy resolution of better than 15\% and the depth of the shower maximum with a resolution of better than 40 g/cm\textsuperscript{2}. To improve and interpret our measurements as well as to study systematic uncertainties due to interaction models, we perform radio simulations with CORSIKA and CoREAS. In this overview we present the setup of Tunka-Rex, discuss the achieved results and the prospects of mass-composition studies with radio arrays.

We discuss the effect of hypothetical violation of Lorentz invariance at high energies on the formation of atmospheric showers by very-high-energy gamma rays. In the scenario where Lorentz invariance violation leads to a decrease of the photon velocity with energy the formation of the showers is suppressed compared to the Lorentz invariant case. Absence of such suppression in the high-energy part of spectrum of the Crab nebula measured independently by HEGRA and H.E.S.S. collaborations is used to set lower bounds on the energy scale of Lorentz invariance violation. These bounds are competitive with the strongest existing constraints obtained from timing of variable astrophysical sources and the absorption of TeV photons on the extragalactic background light. They will be further improved by the next generation of multi-TeV gamma-ray observatories.

Observations by Fermi LAT enabled us to explore the population of non-recycled gamma-ray pulsars with the set of 89 objects. It was recently noted that there are apparent differences in properties of radio-quiet and radio-loud subsets. In particular, average observed radio-loud pulsar is younger than radio-quiet one and is located at smaller galactic latitude. Even so, the analysis based on the full list of pulsars may suffer from selection effects. Namely, most of radio-loud pulsars are first discovered in the radio-band, while radio-quiet ones are found using the gamma-ray data. In this work we perform a blind search for gamma-ray pulsars using the Fermi LAT data alone using all point sources from 3FGL catalog as the candidates. Unlike preceding blind search, the present catalog is constructed with novel semi-coherent method and covers the full range of characteristic ages down to 1 kyr. The search resulted in the catalog of 40 non-recycled pulsars, 26 of which are radio-quiet. There are no statistically significant differences in age and galactic latitude distributions for the radio-loud and radio-quiet pulsars, while the rotation period distributions are marginally different with $2.4\sigma$ pre-trial statistical significance. The fraction of radio-quiet pulsars is estimated as $\epsilon_{RQ}=63\pm 8\%$. The results are in agreement with the predictions of the outer magnitosphere models, while the Polar cap models are disfavored.

The radio technique is a promising method for detection of cosmic-ray air showers of energies around $100\,$PeV and higher with an array of radio antennas. Since the amplitude of the radio signal can be measured absolutely and increases with the shower energy, radio measurements can be used to determine the air-shower energy on an absolute scale. We show that calibrated measurements of radio detectors operated in coincidence with host experiments measuring air showers based on other techniques can be used for comparing the energy scales of these host experiments. Using two approaches, first via direct amplitude measurements, and second via comparison of measurements with air shower simulations, we compare the energy scales of the air-shower experiments Tunka-133 and KASCADE-Grande, using their radio extensions, Tunka-Rex and LOPES, respectively. Due to the consistent amplitude calibration for Tunka-Rex and LOPES achieved by using the same reference source, this comparison reaches an accuracy of approximately $10\,\%$ - limited by some shortcomings of LOPES, which was a prototype experiment for the digital radio technique for air showers. In particular we show that the energy scales of cosmic-ray measurements by the independently calibrated experiments KASCADE-Grande and Tunka-133 are consistent with each other on this level.

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.

One of the possible sources of hadronic cosmic rays (CRs) are newborn pulsars. If this is indeed the case, they should feature diffusive gamma-ray halos produced by interactions of CRs with interstellar gas. In this paper we try to identify extended gamma-ray emission around young pulsars, making use of the 7-year Fermi-LAT data. For this purpose we select and analyze a set of eight pulsars that are most likely to possess detectable gamma-ray halos. We find extended emission that might be interpreted as a gamma-ray halo only in the case of PSR J0007+7303. Its luminosity accords with the total energy of injected cosmic rays $\sim 10^{50}$ erg, although other interpretations of this source are possible. Irrespectively of the nature of this source we put bounds on the luminosity of gamma-ray halos which suggest that pulsars' contribution to the overall energy budget of galactic CRs is subdominant in the GeV-TeV range.

The Moscow State University Extensive Air Shower (EAS-MSU) array studied high-energy cosmic rays with primary energies ~(1-500) PeV in the Northern hemisphere. The EAS-MSU data are being revisited following recently found indications to an excess of muonless showers, which may be interpreted as the first observation of cosmic gamma rays at ~100 PeV. In this paper, we present a complete Monte-Carlo model of the surface detector which results in a good agreement between data and simulations. The model allows us to study the performance of the detector and will be used to obtain physical results in further studies.

TARA (Telescope Array Radar) is a cosmic ray radar detection experiment colocated with Telescope Array, the conventional surface scintillation detector (SD) and fluorescence telescope detector (FD) near Delta, Utah, U.S.A. The TARA detector combines a 40 kW, 54.1 MHz VHF transmitter and high-gain transmitting antenna which broadcasts the radar carrier over the SD array and within the FD field of view, towards a 250 MS/s DAQ receiver. TARA has been collecting data since 2013 with the primary goal of observing the radar signatures of extensive air showers (EAS). Simulations indicate that echoes are expected to be short in duration (~10 microseconds) and exhibit rapidly changing frequency, with rates on the order of 1 MHz/microsecond. The EAS radar cross-section (RCS) is currently unknown although it is the subject of over 70 years of speculation. A novel signal search technique is described in which the expected radar echo of a particular air shower is used as a matched filter template and compared to waveforms obtained by triggering the radar DAQ using the Telescope Array fluorescence detector. No evidence for the scattering of radio frequency radiation by EAS is obtained to date. We report the first quantitative RCS upper limits using EAS that triggered the Telescope Array Fluorescence Detector.

A number of recent studies indicates a significant amount of ionized gas in a form of the hot gas halo around the Milky Way. The halo extends over the region of 100 kpc and may be acountable for the missing baryon mass. In this paper we calculate the contribution of the proposed halo to the dispersion measure (DM) of the pulsars. The Navarro, Frenk & White (NFW), Maller & Bullock (MB) and Feldmann, Hooper & Gnedin (FHG) density distibutions are considered for the gas halo. The data set includes pulsars with the distance known independently from the DM, e.g. pulsars in globular clusters, LMC, SMC and pulsars with known parallax. The results exclude the NFW distribution for the hot gas, while the more realistic MB and FHG models are compatible with the observed dispersion measure.

This paper presents the results of different searches for correlations between very high-energy neutrino candidates detected by IceCube and the highest-energy cosmic rays measured by the Pierre Auger Observatory and the Telescope Array. We first consider samples of cascade neutrino events and of high-energy neutrino-induced muon tracks, which provided evidence for a neutrino flux of astrophysical origin, and study their cross-correlation with the ultrahigh-energy cosmic ray (UHECR) samples as a function of angular separation. We also study their possible directional correlations using a likelihood method stacking the neutrino arrival directions and adopting different assumptions on the size of the UHECR magnetic deflections. Finally, we perform another likelihood analysis stacking the UHECR directions and using a sample of through-going muon tracks optimized for neutrino point-source searches with sub-degree angular resolution. No indications of correlations at discovery level are obtained for any of the searches performed. The smallest of the p-values comes from the search for correlation between UHECRs with IceCube high-energy cascades, a result that should continue to be monitored.

We reconstructed the energy and the position of the shower maximum of air showers with energies $E \gtrsim 100 $PeV applying a method using radio measurements performed with Tunka-Rex. An event-to-event comparison to air-Cherenkov measurements of the same air showers with the Tunka-133 photomultiplier array confirms that the radio reconstruction works reliably. The Tunka-Rex reconstruction methods and absolute scales have been tuned on CoREAS simulations and yield energy and $X_{\mathrm{max}}$ values consistent with the Tunka-133 measurements. The results of two independent measurement seasons agree within statistical uncertainties, which gives additional confidence in the radio reconstruction. The energy precision of Tunka-Rex is comparable to the Tunka-133 precision of $15 %$, and exhibits a $20 %$ uncertainty on the absolute scale dominated by the amplitude calibration of the antennas. For $X_{\mathrm{max}}$, this is the first direct experimental correlation of radio measurements with a different, established method. At the moment, the $X_{\mathrm{max}}$ resolution of Tunka-Rex is approximately $40 $g/cm$^2$. This resolution can probably be improved by deploying additional antennas and by further development of the reconstruction methods, since the present analysis does not yet reveal any principle limitations.

We perform a blind search for the variability of the gamma-ray sky in the energy range E>1 GeV using 308 weeks of the Fermi-LAT data. We use the technique based on the comparison of the weekly photon counts and exposures in sky pixels by means of the Kolmogorov-Smirnov test. We consider the flux variations in the region significant if statistical probability of uniformity is less than $4\times10^{-6}$, which corresponds to 0.05 false detections in the whole set of 12288 pixels. Close inspection of the detected variable regions result in identification of 8 sources without previous known variability. Two of them are included in the second Fermi LAT source catalogue (FBQS J122424.1+243623 and GB6 J0043+3426) and one (3EG J1424+3734) was reported by EGRET and also was included in the First Fermi LAT source catalogue (1FGL), but is missing in the 2FGL. Possible identifications of five other sources are obtained using NED and SIMBAD databases (1RXS J161939.9+765515, PMN J2320-6447, PKS 0226-559, PKS J0030-0211, PMN J0225-2603). These new variable gamma-ray sources demonstrate recurring flaring activity with time scale ~weeks and have hard spectra. Their spectral energy distributions deviate significantly from a simple power-law shape and often peak around ~GeV. These properties of activity are typical for flaring blazars.

The IceCube, Pierre Auger and Telescope Array Collaborations have recently reported results on neutral particles (neutrons, photons and neutrinos) which complement the measurements on charged primary cosmic rays at ultra-high energy. The complementarity between these messengers and between their detections are outlined. The current status of their search is reviewed and a cross-correlation analysis between the available results is performed. The expectations for photon and neutrino detections in the near future are also presented.

We have conducted three searches for correlations between ultra-high energy cosmic rays detected by the Telescope Array and the Pierre Auger Observatory, and high-energy neutrino candidate events from IceCube. Two cross-correlation analyses with UHECRs are done: one with 39 cascades from the IceCube `high-energy starting events' sample and the other with 16 high-energy `track events'. The angular separation between the arrival directions of neutrinos and UHECRs is scanned over. The same events are also used in a separate search using a maximum likelihood approach, after the neutrino arrival directions are stacked. To estimate the significance we assume UHECR magnetic deflections to be inversely proportional to their energy, with values $3^\circ$, $6^\circ$ and $9^\circ$ at 100 EeV to allow for the uncertainties on the magnetic field strength and UHECR charge. A similar analysis is performed on stacked UHECR arrival directions and the IceCube sample of through-going muon track events which were optimized for neutrino point-source searches.

Joint contributions of the Pierre Auger Collaboration and the Telescope Array Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.

Tunka-Rex is a radio detector for cosmic-ray air showers in Siberia, triggered by Tunka-133, a co-located air-Cherenkov detector. The main goal of Tunka-Rex is the cross-calibration of the two detectors by measuring the air-Cherenkov light and the radio signal emitted by the same air showers. This way we can explore the precision of the radio-detection technique, especially for the reconstruction of the primary energy and the depth of the shower maximum. The latter is sensitive to the mass of the primary cosmic-ray particles. In this paper we describe the detector setup and explain how electronics and antennas have been calibrated. The analysis of data of the first season proves the detection of cosmic-ray air showers and therefore, the functionality of the detector. We confirm the expected dependence of the detection threshold on the geomagnetic angle and the correlation between the energy of the primary cosmic-ray particle and the radio amplitude. Furthermore, we compare reconstructed amplitudes of radio pulses with predictions from CoREAS simulations, finding agreement within the uncertainties.

Since its commissioning in autumn 2012, Tunka-Rex, the radio extension of the air-Cherenkov detector Tunka-133, performed three years of air shower measurements. Currently the detector consists of 44 antennas connected to air-Cherenkov and scintillator detectors, respectively, placed in the Tunka valley, Siberia. Triggered by these detectors, Tunka-Rex measures the radio signal up to EeV-scale air-showers. This configuration provides a unique possibility for cross-calibration between air-Cherenkov, radio and particle techniques. We present reconstruction methods for the energy and the shower maximum developed with CoREAS simulations, which allow for a precision competitive with the air-Cherenkov technique. We apply these methods to data acquired by Tunka-Rex in the first year which we use for cross-calibration, and we compare the results with the reconstruction of the energy and the shower maximum by Tunka-133, which provides also a reconstruction for the shower core used for the radio reconstruction. Our methods have shown that the atmospheric depth of the shower maximum ($X_\mathrm{max}$) can be reconstructed with a precision of better than 40 g/cm$^2$ for high quality events, in some cases even when only three antenna stations have signal. The energy precision is comparable with the air-Cherenkov precision of 15%. Soon the results will be checked with the independent data of the second year.

The Tunka observatory is located close to Lake Baikal in Siberia, Russia. Its main detector, Tunka-133, is an array of photomultipliers measuring Cherenkov light of air showers initiated by cosmic rays in the energy range of approximately $10^{16}-10^{18}\,$eV. In the last years, several extensions have been built at the Tunka site, e.g., a scintillator array named Tunka-Grande, a sophisticated air-Cherenkov-detector prototype named HiSCORE, and the radio extension Tunka-Rex. Tunka-Rex started operation in October 2012 and currently features 44 antennas distributed over an area of about $3\,$km$^2$, which measure the radio emission of the same air showers detected by Tunka-133 and Tunka-Grande. Tunka-Rex is a technological demonstrator that the radio technique can provide an economic extension of existing air-shower arrays. The main scientific goal is the cross-calibration with the air-Cherenkov measurements. By this cross-calibration, the precision for the reconstruction of the energy and mass of the primary cosmic-ray particles can be determined. Finally, Tunka-Rex can be used for cosmic-ray physics at energies close to $1\,$EeV, where the standard Tunka-133 analysis is limited by statistics. In contrast to the air-Cherenkov measurements, radio measurements are not limited to dark, clear nights and can provide an order of magnitude larger exposure.

Approximate de Sitter symmetry of inflating Universe is responsible for the approximate flatness of the power spectrum of scalar perturbations. However, this is not the only option. Another symmetry which can explain nearly scale-invariant power spectrum is conformal invariance. We give a short review of models based on conformal symmetry which lead to the scale-invariant spectrum of the scalar perturbations. We discuss also potentially observable features of these models.

In this work we are reporting on the measurement of the proton-air inelastic cross section $\sigma^{\rm inel}_{\rm p-air}$ using the Telescope Array (TA) detector. Based on the measurement of the $\sigma^{\rm inel}_{\rm p-air}$ the proton-proton cross section $\sigma_{\rm p-p}$ value is also determined at $\sqrt{s} = 95_{-8}^{+5}$ TeV. Detecting cosmic ray events at ultra high energies with Telescope Array enables us to study this fundamental parameter that we are otherwise unable to access with particle accelerators. The data used in this report is the hybrid events observed by the Middle Drum fluorescence detector together with the surface array detector collected over five years. The value of the $\sigma^{\rm inel}_{\rm p-air}$ is found to be equal to $567.0 \pm 70.5 [{\rm Stat.}] ^{+29}_{-25} [{\rm Sys.}]$ mb. The total proton-proton cross section is subsequently inferred from Glauber Formalism and Block, Halzen and Stanev QCD inspired fit and is found to be equal to $170_{-44}^{+48} [{\rm Stat.}] _{-17}^{+19} [{\rm Sys.}] $mb.

The Tunka Radio Extension (Tunka-Rex) is an array of 44 radio antenna stations, distributed over 3 km$^{2}$, constituting a radio detector for air showers with an energy threshold around 10$^{17}$ eV. It is an extension to Tunka-133, an air-Cherenkov detector in Siberia, which is used as an external trigger for Tunka-Rex and provides a reliable reconstruction of energy and shower maximum. Each antenna station consists of two perpendicularly aligned active antennas, called SALLAs. An antenna calibration of the SALLA with a commercial reference source enables us to reconstruct the detected radio signal on an absolute scale. Since the same reference source was used for the calibration of LOPES and, in a calibration campaign in 2014, also for LOFAR, these three experiments now have a consistent calibration and, therefore, absolute scale. This was a key ingredient to resolve a longer standing contradiction between measurements of two calibrated experiments. We will present how the calibration was performed and compare radio measurements of air showers from Tunka-Rex to model calculations with the radio simulation code CoREAS, confirming it within the scale uncertainty of the calibration of 18%.

Tunka-Rex is the radio extension of the Tunka cosmic-ray observatory in Siberia close to Lake Baikal. Since October 2012 Tunka-Rex measures the radio signal of air-showers in coincidence with the non-imaging air-Cherenkov array Tunka-133. Furthermore, this year additional antennas will go into operation triggered by the new scintillator array Tunka-Grande measuring the secondary electrons and muons of air showers. Tunka-Rex is a demonstrator for how economic an antenna array can be without losing significant performance: we have decided for simple and robust SALLA antennas, and we share the existing DAQ running in slave mode with the PMT detectors and the scintillators, respectively. This means that Tunka-Rex is triggered externally, and does not need its own infrastructure and DAQ for hybrid measurements. By this, the performance and the added value of the supplementary radio measurements can be studied, in particular, the precision for the reconstructed energy and the shower maximum in the energy range of approximately $10^{17}-10^{18}\,$eV. Here we show first results on the energy reconstruction indicating that radio measurements can compete with air-Cherenkov measurements in precision. Moreover, we discuss future plans for Tunka-Rex.

The angular extension of the gamma-ray image of 3C 279 may be constrained by studying its solar occultations as suggested by Fairbairn et al. (2010). We perform this kind of analysis for seven occultations observed by Fermi-LAT in 2008-2014, using the Fermi-LAT Solar System tools. The results are interpreted in terms of models with extended gamma-ray halo of 3C 279; first constraints on the size and the flux of the halo are reported.

We report on the search for steady point-like sources of neutral particles around 10$^{18}$ eV between 2008 May and 2013 May with the scintillator surface detector of the Telescope Array experiment. We found overall no significant point-like excess above 0.5 EeV in the northern sky. Subsequently, we also searched for coincidence with the Fermi bright Galactic sources. No significant coincidence was found within the statistical uncertainty. Hence, we set an upper limit on the neutron flux that corresponds to an averaged flux of 0.07 km$^{-2}$ yr$^{-1}$ for $E>1$ EeV in the northern sky at the 95% confidence level. This is the most stringent flux upper limit in a northern sky survey assuming point-like sources. The upper limit at the 95% confidence level on the neutron flux from Cygnus X-3 is also set to 0.2 km$^{-2}$ yr$^{-1}$ for $E>0.5$ EeV. This is an order of magnitude lower than previous flux measurements.

It is known that the high energy $(> 100\,\text{MeV})$ emission of gamma-ray bursts is delayed with respect to the low energy emission. However, the dependence of light curves on energy has not been studied for the high energy bands. In this paper we consider the bursts observed by Fermi LAT from 2008 August 4 to 2011 August 1, for which at least $10$ photons were observed with the energy greater than $1\,\text{GeV}$. These include $4$ bursts: GRB 080916C, GRB 090510, GRB 090902B, and GRB 090926A. We use the Kolmogorov-Smirnov test to compare the light curves in the two bands, $100\,\text{MeV} < E < 1\,\text{GeV}$ and $E > 1\,\text{GeV}$. For GRB 080916C and GRB 090510 the light curves in the two bands are statistically compatible. However, for GRB 090926A, the higher-energy light curve is stretched compared to the lower-energy one with a statistical significance of $3.3 \sigma$ and, for GRB 090902B, on the contrary, the lower-energy curve is stretched with $2.3 \sigma$ significance. We argue that the observed diversity of stretching factors may be explained in a simple geometrical model. The model assumes that the jet opening angle depends on the emission energy in a way that the most energetic photons are radiated near the axis of the jet. All the bursts are considered equivalent in their rest frames and the observed light curves differ only due to different redshifts and view directions. The model conforms to the total burst energy constraint and matches the Fermi LAT observations of the fraction of GRBs visible in $100\,\text{MeV} < E < 1\,\text{GeV}$ band, which may be observed at higher energies. The model predicts the distribution of observable stretching factors, which may be tested in the future data. Finally, we propose a method to estimate observer's off-axis angle based on the stretching factor and the fraction of the high-energy photons.

We revisit constraints on the (pseudo)conformal Universe from the non-observation of statistical anisotropy in the Planck data. The quadratic maximal likelihood estimator is applied to the Planck temperature maps at frequencies 143 GHz and 217 GHz as well as their cross-correlation. The strongest constraint is obtained in the scenario of the (pseudo)conformal Universe with a long intermediate evolution after conformal symmetry breaking. In terms of the relevant parameter (coupling constant), the limit is h^2 <0.0013 at 95% C.L. (using the cross-estimator). The analogous limit is much weaker in the scenario without the intermediate stage (h^2 \ln \frac{H_0}{\Lambda}<0.52) allowing the coupling constant to be of order one. In the latter case, the non-Gaussianity in the 4-point function appears to be a more promising signature.

Previous measurements of the composition of Ultra-High Energy Cosmic Rays(UHECRs) made by the High Resolution Fly's Eye(HiRes) and Pierre Auger Observatory(PAO) are seemingly contradictory, but utilize different detection methods, as HiRes was a stereo detector and PAO is a hybrid detector. The five year Telescope Array(TA) Middle Drum hybrid composition measurement is similar in some, but not all, respects in methodology to PAO, and good agreement is evident between data and a light, largely protonic, composition when comparing the measurements to predictions obtained with the QGSJetII-03 and QGSJet-01c models. These models are also in agreement with previous HiRes stereo measurements, confirming the equivalence of the stereo and hybrid methods. The data is incompatible with a pure iron composition, for all models examined, over the available range of energies. The elongation rate and mean values of Xmax are in good agreement with Pierre Auger Observatory data. This analysis is presented using two methods: data cuts using simple geometrical variables and a new pattern recognition technique.

The Fermi Large Area Telescope (LAT) has observed more than a hundred of gamma-ray pulsars, about one third of which are radio-quiet, i.e. not detected at radio frequencies. The most of radio-loud pulsars are detected by Fermi LAT by using the radio timing models, while the radio-quiet ones are discovered in a blind search. The difference in the techniques introduces an observational selection bias and, consequently, the direct comparison of populations is complicated. In order to produce an unbiased sample, we perform a blind search of gamma-ray pulsations using Fermi-LAT data alone. No radio data or observations at optical or X-ray frequencies are involved in the search process. We produce a gamma-ray selected catalog of 25 non-recycled gamma-ray pulsars found in a blind search, including 16 radio-quiet and 9 radio-loud pulsars. This results in the direct measurement of the fraction of radio-quiet pulsars $\varepsilon_{RQ} = 64\pm 10\%$, which is in agreement with the existing estimates from the population modeling in the outer magnetosphere model. The Polar cap models are disfavored due to a lower expected fraction and the prediction of age dependence. The age, gamma-ray energy flux, spin-down luminosity and sky location distributions of the radio-loud and radio-quiet pulsars from the catalog do not demonstrate any statistically significant difference. The results indicate that the radio-quiet and radio-loud pulsars belong to one and the same population. The catalog shows no evidence for the radio beam evolution.

Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage. This can be achieved by combining data from observatories located in both the northern and southern hemispheres. To this end, a joint analysis using data recorded at the Telescope Array and the Pierre Auger Observatory above $10^{19}$ eV is presented in this work. The resulting multipolar expansion of the flux of cosmic rays allows us to perform a series of anisotropy searches, and in particular to report on the angular power spectrum of cosmic rays above $10^{19}$ eV. No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole.

We have searched for intermediate-scale anisotropy in the arrival directions of ultrahigh-energy cosmic rays with energies above 57~EeV in the northern sky using data collected over a 5 year period by the surface detector of the Telescope Array experiment. We report on a cluster of events that we call the hotspot, found by oversampling using 20$^\circ$-radius circles. The hotspot has a Li-Ma statistical significance of 5.1$\sigma$, and is centered at R.A.=146.7$^{\circ}$, Dec.=43.2$^{\circ}$. The position of the hotspot is about 19$^{\circ}$ off of the supergalactic plane. The probability of a cluster of events of 5.1$\sigma$ significance, appearing by chance in an isotropic cosmic-ray sky, is estimated to be 3.7$\times$10$^{-4}$ (3.4$\sigma$).

We point out that violation of Lorentz invariance affects the interaction of high-energy photons with the Earth's atmosphere and magnetic field. In certain parameter region this interaction becomes suppressed and the photons escape observation passing through the atmosphere without producing air showers. We argue that a detection of photon-induced air showers with energies above 10^19 eV, implying the absence of suppression as well as the absence of photon decay, will put tight double-sided limits on Lorentz violation in the sector of quantum electrodynamics. These constraints will be by several orders of magnitude stronger than the existing ones and will be robust against any assumptions about the astrophysical origin of the detected photons.

The expected opacity of the intergalactic space limits the mean free path of TeV gamma rays to dozens of Megaparsecs. However, TeV photons from numerous more distant sources have been detected. This might be interpreted, in each particular case, in terms of hardening of the emitted spectrum caused by presently unknown mechanisms at work in the sources. Here we show that this interpretation is not supported by the analysis of the ensemble of all observed sources. In the frameworks of an infrared-background model with the lowest opacity, we reconstruct the emitted spectra of distant blazars and find that upward spectral breaks appear precisely at those energies where absorption effects are essential. Since these energies are very different for similar sources located at various distances, we conclude that the breaks are artefacts of the incorrect account of absorption and, therefore, the opacity of the Universe for gamma rays is overestimated even in the most conservative model. This implies that some novel physical or astrophysical phenomena should affect long-distance propagation of gamma rays. A scenario in which a part of energetic photons is converted to an inert new particle in the vicinity of the source and reconverts back close to the observer does not contradict to our results. [abridged]

We constrain several models of the early Universe that predict a statistical anisotropy of the cosmic microwave background (CMB) sky. We make use of WMAP9 maps deconvolved with beam asymmetries. As compared to previous releases of WMAP data, they do not exhibit the anomalously large quadrupole of statistical anisotropy. This allows to strengthen the limits on the parameters of models established earlier in the literature. In particular, the amplitude of the special quadrupole is constrained as |g_*|<0.072 at 95% C.L. (-0.046<g_*<0.048 at 68% C.L.) independently of the preferred direction in the sky. The upper limit is obtained on the total number of e-folds in anisotropic inflation with the Maxwellian term nonminimally coupled to the inflaton, namely N_{tot} <N_{CMB} +82 at 95% C.L. (+14 at 68% C.L.) for N_{CMB}=60. We also constrain models of the (pseudo)conformal universe. The strongest constraint is obtained for spectator scenarios involving a long stage of subhorizon evolution after conformal rolling, which reads h^2 < 0.006 at 95% C.L., in terms of the relevant parameter. The analogous constraint is much weaker in dynamical models, e.g., Galilean genesis.

We search for ultra-high energy photons by analyzing geometrical properties of shower fronts of events registered by the Telescope Array surface detector. By making use of an event-by-event statistical method, we derive upper limits on the absolute flux of primary photons with energies above 10^19, 10^19.5 and 10^20 eV based on the first three years of data taken.

Tunka-Rex is the radio extension of Tunka-133 located in Siberia close to Lake Baikal. The latter is a photomultiplier array registering air-Cherenkov light from air showers induced by cosmic-ray particles with initial energies of approximately $10^{16}$ to $10^{18}$ eV. Tunka-Rex extends this detector with 25 antennas spread over an area of 1 km$^2$. It is triggered externally by Tunka-133, and detects the radio emission of the same air showers. The combination of an air-Cherenkov and a radio detector provides a facility for hybrid measurements and cross-calibration between the two techniques. The main goal of Tunka-Rex is to determine the precision of the reconstruction of air-shower parameters using the radio detection technique. It started operation in autumn 2012. We present the overall concept of Tunka-Rex, the current status of the array and first analysis results.

We revisit cosmic microwave background (CMB) constraints on the abundance of millicharged particles based on the Planck data. The stringent limit Omega_{mcp}h^2 < 0.001 (95% CL) may be set using the CMB data alone if millicharged particles participate in the acoustic oscillations of baryon-photon plasma at the recombination epoch. The latter condition is valid for a wide region of charges and masses of the particles. Adding the millicharged component to LCDM shifts prefered scalar spectral index of primordial perturbations to somewhat larger values as compared to minimal model, even approaching Harrison-Zeldovich spectrum under some assumptions.