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.

Strong explosion of a compact star surrounded by a thick stellar wind drives a fast ($>0.1c$) radiation mediated shock (RMS) that propagates in the wind, and ultimately breaks out gradually once photons start escaping from the shock transition layer. In exceptionally strong or aspherical explosions the shock velocity may even be relativistic. The properties of the breakout signal depend on the dynamics and structure of the shock during the breakout phase. Here we present, for the first time, spectra and lightcurves of the breakout emission of fast Newtonian and mildly relativistic shocks, that were calculated using self-consistent Monte-Carlo simulations of finite RMS with radiative losses. We find a strong dependence of the $\nu F_\nu$ peak on shock velocity, ranging from $\sim 1$ keV for $v_s/c=0.1$ to $\sim 100$ keV for $v_s/c=0.5$, with a shift to lower energies as losses increase. For all cases studied the spectrum below the peak exhibits a nearly flat component ($F_\nu \sim \nu^0$) that extends down to the break frequency below which absorption becomes important. This implies much bright optical/UV emission than hitherto expected. The computed lightcurves show a gradual rise over tens to hundreds of seconds for representative conditions. The application to SN 2008D/XRT 080109 and the detectability limits are also discussed. We predict a detection rate of about one per year with eROSITA.

Chemical extraction using a molecular recognition resin named "Empore Radium Rad Disk" was developed to improve sensitivity for the low concentration of radium (Ra). Compared with the previous method, the extraction process speed was improved by a factor of three and the recovery rate for $^{226}$Ra was also improved from 81$\pm$4% to $>$99.9%. The sensitivity on the 10$^{-1}$ mBq level was achieved using a high purity germanium detector. This improved method was applied to determine $^{226}$Ra in Gd$_2$(SO$_4$)$_3{\cdot}$8H$_2$O which will be used in the Super-Kamiokande Gadolinium project. The improvement and measurement results are reported in this paper.

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.

Radiation mediated shocks (RMS) play a key role in shaping the early emission observed in many transients. In most cases, e.g., shock breakout in supernovae, llGRBs and neutron star mergers, the upstream plasma is devoid of radiation, and the photons that ultimately reach the observer are generated predominantly inside and downstream of the shock. Predicting the observed spectrum requires detailed calculations of the shock structure and thermodynamic state that account properly for the shock microphysics. We present results of self-consistent Monte-Carlo simulations of photon-starved RMS, that yield the shock structure and emission for a broad range of shock velocities, from sub-relativistic ($\beta_{sh} = 0.1$) to highly relativistic ($\Gamma_{sh} = 20$). Our simulations confirm that in relativistic RMS the immediate downstream temperature is regulated by exponential pair creation, ranging from $50$ keV at $\beta_{sh}=0.5$ to $200$ keV at $\Gamma_{sh}=20$. At lower velocities the temperature becomes sensitive to the shock velocity, with $kT\sim 0.5$ keV at $\beta_{sh}=0.1$. We also confirm that in relativistic shocks the opacity is completely dominated by newly created pairs, which has important implications for the breakout physics. We find the transition to pair dominance to occur at $\beta_{sh}=0.5$ roughly. In all cases examined, the spectrum below the $\nu F_\nu$ peak has been found to be substantially softer than the Planck distribution. This has important implications for the optical emission in fast and relativistic breakouts, and their detection. The applications to GRB 060218 and GRB 170817A are discussed.

We report a search for astronomical neutrinos in the energy region from several GeV to TeV in the direction of the blazar TXS0506+056 using the Super-Kamiokande detector following the detection of a 100 TeV neutrino from the same location by the IceCube collaboration. Using Super-Kamiokande neutrino data across several data samples observed from April 1996 to February 2018 we have searched for both a total excess above known backgrounds across the entire period as well as localized excesses on smaller time scales in that interval. No significant excess nor significant variation in the observed event rate are found in the blazar direction. Upper limits are placed on the electron and muon neutrino fluxes at 90\% confidence level as $6.03 \times 10^{-7}$ and $4.52 \times 10^{-7}$ to $9.26 \times 10^{-10}$ [${\rm erg}/{\rm cm}^2/{\rm s}$], respectively.

An important issue for rare-event-search experiments, such as the search for dark matter or neutrinoless double beta decay, is to reduce radioactivity of the detector materials and the experimental environment. The selection of materials with low radioactive impurities, such as isotopes of the uranium and thorium chains, requires a precise measurement of surface and bulk radioactivity. Focused on the first one, an alpha-particle detector has been developed based on a gaseous micro-time-projection chamber. A low-alpha mu-PIC with reduced alpha-emission background was installed in the detector. The detector offers the advantage of position sensitivity, which allows the alpha-particle contamination of the sample to be imaged and the background to be measured at the same time. The detector performance was measured by using an alpha-particle source. The measurement with a sample was also demonstrated and the sensitivity is discussed.

Supernova detection is a major objective of the Super-Kamiokande (SK) experiment. In the next stage of SK (SK-Gd), gadolinium (Gd) sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons. A core-collapse supernova will be preceded by an increasing flux of neutrinos and anti-neutrinos, from thermal and weak nuclear processes in the star, over a timescale of hours; some of which may be detected at SK-Gd. This could provide an early warning of an imminent core-collapse supernova, hours earlier than the detection of the neutrinos from core collapse. Electron anti-neutrino detection will rely on inverse beta decay events below the usual analysis energy threshold of SK, so Gd loading is vital to reduce backgrounds while maximising detection efficiency. Assuming normal neutrino mass ordering, more than 200 events could be detected in the final 12 hours before core collapse for a 15-25 solar mass star at around 200 pc, which is representative of the nearest red supergiant to Earth, $\mathrm{\alpha}$Ori (Betelgeuse). At a statistical false alarm rate of 1 per century, detection could be up to 10 hours before core collapse, and a pre-supernova star could be detected by SK-Gd up to 600 pc away. A pre-supernova alert could be provided to the astrophysics community following gadolinium loading.

Long duration gamma-ray bursts (GRBs), the brightest events since the Big Bang itself, are believed to originate in an ultra-relativistic jet breaking out from a massive stellar envelope. Despite decades of study, there is still no consensus on their emission mechanism. One unresolved question is the origin of the tight correlation between the spectral peak energy Ep and peak luminosity Lp discovered in observations. This Yonetoku relation is the tightest correlation found in the properties of the prompt phase of GRB emission, providing the best diagnostic for the radiation mechanism. Here we present 3D hydrodynamical simulations, and post-process radiation transfer calculations, of photospheric emission from a relativistic jet. Our simulations reproduce the Yonetoku relation as a natural consequence of viewing angle. Although jet dynamics depend sensitively on luminosity, the Ep-Lp correlation holds regardless. This result strongly suggests that photospheric emission is the dominant component in the prompt phase of GRBs.

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

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.

In the standard synchrotron afterglow model, a power law of electrons is responsible for all aspects of photon production and absorption. Recent numerical work has shown that the vast majority of particles in the downstream medium are actually "thermal" particles, which were shock-heated but did not enter the Fermi acceleration process (the name stands in contrast to the nonthermal high-energy tail, rather than connoting a Maxwellian distribution). There are substantial differences at optical and higher energies when these thermal electrons participate in the afterglow, but early work along these lines ignored the radio end of the electromagnetic spectrum. We report here on an extension of previous Monte Carlo simulations of gamma-ray burst afterglows. The model now includes the synchrotron self-absorption (SSA) process and so can simulate afterglows across the entire EM spectrum, and several orders of magnitude in time. In keeping with earlier work, inclusion of the thermal electrons increases the SSA frequency by a factor of 30, and the radio intensity by a factor of 100. Furthermore, these changes happen with no modification to the late optical or X-ray afterglow. Our results provide very strong evidence that thermal electrons must be considered in any multiwavelength model for afterglows.

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.

We explore the physics of relativistic radiation mediated shocks (RRMSs) in the regime where photon advection dominates over photon generation. For this purpose, a novel iterative method for deriving a self-consistent steady-state structure of RRMS is developed, based on a Monte-Carlo code that solves the transfer of photons subject to Compton scattering and pair production/annihilation. Systematic study is performed by imposing various upstream conditions which are characterized by the following three parameters: the photon-to-baryon inertia ratio $\xi_{u *}$, the photon-to-baryon number ratio $\tilde{n}$, and the shock Lorentz factor $\gamma_u$. We find that the properties of RRMSs vary considerably with these parameters. In particular, while a smooth decline in the velocity, accompanied by a gradual temperature increase is seen for $\xi_{u*} \gg 1$, an efficient bulk Comptonization, that leads to a heating precursor, is found for $\xi_{u*} \lesssim 1$. As a consequence, although particle acceleration is highly inefficient in these shocks, a broad non-thermal spectrum is produced in the latter case. The generation of high energy photons through bulk Comptonization leads, in certain cases, to a copious production of pairs that provide the dominant opacity for Compton scattering. We also find that for certain upstream conditions a weak subshock appears within the flow. For a choice of parameters suitable to gamma-ray bursts, the radiation spectrum within the shock is found to be compatible with that of the prompt emission, suggesting that subphotospheric shocks may give rise to the observed non-thermal features despite the absence of accelerated particles.

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.

We developed a three-dimensional radiative transfer code for an ultra-relativistic background flow-field by using the Monte Carlo (MC) method in the context of gamma-ray burst (GRB) emission. For obtaining reliable simulation results in the coupled computation of MC radiation transport with relativistic hydrodynamics which can reproduce GRB emission, we validated radiative transfer computation in the ultra-relativistic regime and assessed the appropriate simulation conditions. The radiative transfer code was validated through two test calculations: (1) computing in different inertial frames and (2) computing in flow-fields with discontinuous and smeared shock fronts. The simulation results of the angular distribution and spectrum were compared among three different inertial frames and in good agreement with each other. If the time duration for updating the flow-field was sufficiently small to resolve a mean free path of a photon into ten steps, the results were thoroughly converged. The spectrum computed in the flow-field with a discontinuous shock front obeyed a power-law in frequency whose index was positive in the range from 1 to 10 MeV. The number of photons in the high-energy side decreased with the smeared shock front because the photons were less scattered immediately behind the shock wave due to the small electron number density. The large optical depth near the shock front was needed for obtaining high-energy photons through bulk Compton scattering. Even one-dimensional structure of the shock wave could affect the results of radiation transport computation.

We explore physical properties of the shocked external medium (i.e., a shell) in 3C 84 associated with the recurrent radio lobe born around 1960. In the previous work of Ito et al., we investigated a dynamical and radiative evolution of such a shell after the central engine stops the jet launching and we found that a fossil shell emission overwhelms that of the rapidly fading radio lobe. We apply this model to 3C 84 and find the followings: (i) The fossil shell made of shocked diffuse ambient matter with the number density of 0.3 cm$^{-3}$ radiates bright Inverse-Compton (IC) emission with the seed photons of the radio emission from the central compact region and the IC emission is above the sensitivity threshold of the Cherenkov Telescope Array (CTA). (ii) When the fossil shell is produced in a geometrically thick ionized plasma with the number density of $10^{3}$ cm$^{-3}$ and the field strength in the shell may reach about 17 mG in the presence of magnetic fields amplification and the radio emission becomes comparable to the sensitivity of deep imaging VLBI observations. A possible production of ultra high energy cosmic-rays (UHECRs) in the dense shocked plasma is also argued.

We analyzed the Suzaku data of M86 and its adjacent regions to study the extended emission around it. The M86 core, the plume, and the tail extending toward the northwest were clearly detected, as well as the extended halo around them. From the position angle $\sim45\deg$ to $\sim275\deg$, the surface brightness distribution of the core and the extended halo was represented relatively well with a single $\beta$-model of $\beta\sim0.5$ up to 15'-20'. The X-ray spectra of the core was represented with a two-temperature model of $kT\sim 0.9$ keV and $\sim0.6$ keV. The temperatures of the core and the halo have a positive gradient in the center, and reach the maximum of $kT\sim1.0$ keV at $r\sim7'$, indicating that the halo gas is located in a larger scale potential structure than that of the galaxy. The temperatures of the plume and the tail were $0.86\pm0.01$ keV and $1.00\pm0.01$ keV. We succeeded in determining the abundances of $\alpha$-element separately, for the core, the plume, the tail, and the halo for the first time. Abundance ratios with respect to Fe were consistent with the solar ratios everywhere, except for Ne. The abundance of Fe was $\sim0.7$ in the core and in the plume, while that in the tail was $\sim1.0$, but the difference was not significant considering the uncertainties of the ICM. The abundance of the halo was almost the same up to $r\sim10'$, and then it becomes significantly smaller (0.2-0.3) at $r\gtrsim10'$, indicating the gas with low metal abundance still remains in the outer halo. From the surface brightness distribution, we estimated the gas mass ($\sim3\times10^{10}M_\odot$) and the dynamical mass ($\sim3\times10^{12}M_\odot$) in $r<100$~kpc. The gas mass to the dynamical mass ratio was $10^{-3}$-$10^{-2}$, suggesting a significant fraction of the halo gas has been stripped.

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.

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.

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 analyse the environment of the supermassive black hole (SMBH) in the centre of a massive elliptical galaxy NGC 1275 in the Perseus cluster, hosting the radio source 3C 84. We focus on the young radio lobe observed inside the estimated Bondi accretion radius. We discuss the momentum balance between the jet associated with the lobe and the surrounding gas. The results are compared with the proper motion of the radio lobe obtained with the very long baseline interferometry. We find that under assumption of a high-density environment >~ 100 cm^-3), the jet power must be comparable to the Eddington luminosity --- this is clearly inconsistent with the current moderate activity of 3C 84, which indicates instead that the jet is expanding in a very low density region (<~1 cm^-3), along the rotation axis of the accretion flow. The power required for the jet to expand in the low-density environment is comparable to the past average jet power estimated from the X-ray observations. We estimate the classical Bondi accretion rate, assuming that (1) gas accretion is spherically symmetric, (2) accretion is associated with the jet environment, and (3) the medium surrounding the jet is representative of the properties of the dominant accreting gas. We find that Bondi accretion is inconsistent with the estimated jet power. This means that either accretion of the cold gas in the NGC 1275 is more efficient than that of the hot gas, or the jets are powered by the SMBH spin.

We explore the photospheric emission from a relativistic jet breaking out from a massive stellar envelope based on relativistic hydrodynamical simulations and post-process radiation transfer calculations in three dimensions. To investigate the impact of three-dimensional (3D) dynamics on the emission, two models of injection conditions are considered for the jet at the center of the progenitor star: one with periodic precession and another without precession. We show that structures developed within the jet due to the interaction with the stellar envelope, as well as due to the precession, have a significant imprint on the resulting emission. Particularly, we find that the signature of precession activity by the central engine is not smeared out and can be directly observed in the light curve as a periodic signal. We also show non-thermal features, which can account for observations of gamma-ray bursts, are produced in the resulting spectra, even though only thermal photons are injected initially and the effect of non-thermal particles is not considered.

We investigate shell emission associated with dying radio loud AGNs. First, based on our recent work by Ito et al. (2015), we describe the dynamical and spectral evolutions of shells after stopping the jet energy injection. We find that the shell emission overwhelms that of the radio lobes soon after stopping the jet energy injection because fresh electrons are continuously supplied into the shell via the forward shock while the radio lobes rapidly fade out without jet energy injection. We find that such fossil shells can be a new class of target sources for SKA telescope. Next, we apply the model to the nearby radio source 3C84. Then, we find that the fossil shell emission in 3C84 is less luminous in radio band while it is bright in TeV gamma-ray band and it can be detectable by CTA.

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.

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.

Matter mixing is one important topic in the study of core-collapse supernova (CCSN) explosions. In this paper, we perform two-dimensional hydrodynamic simulations to reproduce the high velocity $^{56}$Ni clumps observed in SN 1987A. This is the first time that large density perturbation is proposed in the CCSN progenitor to generate Rayleigh-Taylor (RT) instability and make the effective matter mixing. In the case of a spherical explosion, RT instability is efficient at both C+O/He and He/H interfaces of the SN progenitor. Radial coherent structures shown in perturbation patterns are important for obtaining high velocity $^{56}$Ni clumps. We can also obtain matter mixing features and high velocity $^{56}$Ni clumps in some cases of aspherical explosion. We find that one of the most favorable models in our work has a combination of bipolar and equatorially asymmetric explosions in which at least 25\% of density perturbation is introduced at different composition interfaces of the CCSN progenitor. These simulation results are comparable to the observational findings of SN 1987A.

We examine the fate of a dead radio source in which jet injection from the central engine has stopped at an early stage of its evolution ($t = t_j \lesssim 10^5$ yr). To this aim, we theoretically evaluate the evolution of the emission from both the lobe and the shell, which are composed of shocked jet matter and a shocked ambient medium, respectively. Based on a simple dynamical model of expanding lobe and shell, we clarify how the broadband spectrum of each component evolves before and after the cessation of the jet activity. It is shown that the spectrum is strongly dominated by the lobe emission while the jet is active ($t \leq t_j$). On the other hand, once the jet activity has ceased ($t > t_j$), the lobe emission fades out rapidly, since fresh electrons are no longer supplied from the jet. Meanwhile, shell emission only shows a gradual decrease, since accelerated electrons are continuously supplied from the bow shock that is propagating into the ambient medium. As a result, overall emission from the shell overwhelms that from the lobe at wide range of frequencies from radio up to gamma-ray soon after the jet activity has ceased. Our result predicts a new class of dead radio sources that are dominated by shell emission. We suggest that the emission from the shell can be probed in particular at a radio wavelengths with the Square Kilometer Array (SKA) phase 1.

We calculate the electron acceleration in random superluminal strong waves (SLSWs) and radiation from them by using numerical methods in the context of the termination shock of the pulsar wind nebulae. We pursue the electrons by solving the equation of motion in the analytically expressed electromagnetic turbulences. These consist of primary SLSW and isotropically distributed secondary electromagnetic waves. Under the dominance of the secondary waves, all electrons gain nearly equal energy. On the other hand, when the primary wave is dominant, selective acceleration occurs. The phase of the primary wave felt by the electrons moving nearly along the wavevector changes very slowly compared to the oscillation of the wave, which is called "phase locked", and such electrons are continuously accelerated. This acceleration by SLSWs may play a crucial role in the pre-acceleration for the shock acceleration. In general, the radiation from the phase-locked population is different from the synchro-Compton radiation. However, when the amplitude of the secondary waves is not extremely weaker than that of the primary wave, the typical frequency can be estimated from the synchro-Compton theory by using the secondary waves. The primary wave does not contribute to the radiation, because the SLSW accelerates electrons almost linearly. This radiation can be observed as a radio knot at the upstream of the termination shock of the pulsar wind nebulae without counter parts in higher frequency range.

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.

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.

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

Thin and large area NaI(Tl) scintillator to search for WIMPs dark matter was developed. The performance of thin and wide area NaI(Tl) showed good enough to search for dark matter. The energy threshold was as low as 2keV and the energy resolution was about 24% in FWHM at 60keV.

We explore the spectral and polarization properties of photospheric emissions from stratified jets in which multiple components, separated by a sharp velocity shear regions, are distributed in lateral direction. Propagation of thermal photons injected at high optical depth region are calculated until they escape from the photosphere. It is found that presence of the lateral structure within the jet leads to non-thermal feature of the spectra and significant polarization signal in the resulting emission. The deviation from thermal spectra as well as the polarization degree tends to be enhanced as the velocity gradient in the shear region increases. In particular, we show that emissions from multi-component jet can reproduce the typical observed spectra of gamma-ray bursts (GRBs) irrespective to the position of the observer when a velocity shear region is closely spaced in various lateral ($\theta$) positions. The degree of polarization associated in the emission is significant (> few%) at wide range of observer angles and can be higher than 30%.

Joint contributions of the Pierre Auger and Telescope Array Collaborations to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013: cross-calibration of the fluorescence telescopes, large scale anisotropies and mass composition.

We explore photospheric emissions from stratified two-component jets, wherein a highly relativistic spine outflow is surrounded by a wider and less relativistic sheath outflow. Thermal photons are injected in regions of high optical depth and propagated until they escape at the photosphere. Due to the presence of shear in velocity (Lorentz factor) at the boundary of the spine and sheath region, a fraction of the injected photons are accelerated via a Fermi-like acceleration mechanism such that a high energy power-law tail is formed in the resultant spectrum. We show, in particular, that if a velocity shear with a considerable variance in the bulk Lorentz factor is present, the high energy part of observed Gamma-ray Bursts (GRBs) photon spectrum can be explained by this photon acceleration mechanism. We also show that the accelerated photons may also account for the origin of the extra hard power-law component above the bump of the thermal-like peak seen in some peculiar bursts (e.g., GRB 090510, 090902B, 090926A). It is demonstrated that time-integrated spectra can also reproduce the low energy spectrum of GRBs consistently due to a multi-temperature effect when time evolution of the outflow is considered. Finally, we show that the empirical Ep-Lp relation can be explained by differences in the outflow properties of individual sources.

We search for correlations between positions of extragalactic objects and arrival directions of Ultra-High Energy Cosmic Rays (UHECRs) with primary energy $E \ge 40$ EeV as observed by the surface detector array of the Telescope Array (TA) experiment during the first 40 months of operation. We examined several public astronomical object catalogs, including the Veron-Cetty and Veron catalog of active galactic nuclei. We counted the number of TA events correlated with objects in each catalog as a function of three parameters: the maximum angular separation between a TA event and an object, the minimum energy of the events, and the maximum redshift of the objects. We determine combinations of these parameters which maximize the correlations, and calculate the chance probabilities of having the same levels of correlations from an isotropic distribution of UHECR arrival directions. No statistically significant correlations are found when penalties for scanning over the above parameters and for searching in several catalogs are taken into account.

We perform two-dimensional axisymmetric hydrodynamic simulations of matter mixing in aspherical core-collapse supernova explosions of a 16.3 $M_{\odot}$ star with a compact hydrogen envelope. Observations of SN 1987A have provided evidence that $^{56}$Ni synthesized by explosive nucleosynthesis is mixed into fast moving matter ($\gtrsim$ 3,500 km s$^{-1}$) in the exploding star. In order to clarify the key conditions for reproducing such high velocity of $^{56}$Ni, we revisit matter mixing in aspherical core-collapse supernova explosions. Explosions are initiated artificially by injecting thermal and kinetic energies around the interface between the iron core and the silicon-rich layer. Perturbations of 5% or 30% amplitude in the radial velocities are introduced at several points in time. We found that no high velocity $^{56}$Ni can be obtained if we consider bipolar explosions with perturbations (5% amplitude) of pre-supernova origins. If large perturbations (30% amplitude) are introduced or exist due to some unknown mechanism in a later phase just before the shock wave reaches the hydrogen envelope, $^{56}$Ni with a velocity of 3,000 km s$^{-1}$ can be obtained. Aspherical explosions that are asymmetric across the equatorial plane with clumpy structures in the initial shock waves are investigated. We found that the clump sizes affect the penetration of $^{56}$Ni. Finally, we report that an aspherical explosion model that is asymmetric across the equatorial plane with multiple perturbations of pre-supernova origins can cause the penetration of $^{56}$Ni clumps into fast moving matter of 3,000 km s$^{-1}$. We show that both aspherical explosion with clumpy structures and perturbations of pre-supernova origins may be necessary to reproduce the observed high velocity of $^{56}$Ni. To confirm this, more robust three-dimensional simulations are required.

We present detailed chemical abundances for the bright carbon-enhanced metal-poor (CEMP) star BD+44 493, previously reported on by Ito et al. Our measurements confirm that BD+44 493 is an extremely metal-poor ([Fe/H]=-3.8) subgiant star with excesses of carbon and oxygen. No significant excesses are found for nitrogen and neutron-capture elements (the latter of which place it in the CEMP-no class of stars). Other elements that we measure exhibit abundance patterns that are typical for non-CEMP extremely metal-poor stars. No evidence for variations of radial velocity have been found for this star. These results strongly suggest that the carbon enhancement in BD+44 493 is unlikely to have been produced by a companion asymptotic giant-branch star and transferred to the presently observed star, nor by pollution of its natal molecular cloud by rapidly-rotating, massive, mega metal-poor ([Fe/H] < -6.0) stars. A more likely possibility is that this star formed from gas polluted by the elements produced in a "faint" supernova, which underwent mixing and fallback, and only ejected small amounts of elements of metals beyond the lighter elements. The Li abundance of BD+44 493 (A(Li)=log(Li/H)+12=1.0) is lower than the Spite plateau value, as found in other metal-poor subgiants. The upper limit on Be abundance (A(Be)=log(Be/H)+12<-1.8) is as low as those found for stars with similarly extremely-low metallicity, indicating that the progenitors of carbon- (and oxygen-) enhanced stars are not significant sources of Be, or that Be is depleted in metal-poor subgiants with effective temperatures of ~5400K.

We measure the spectrum of cosmic rays with energies greater than $10^{18.2}$ eV with the Fluorescence Detectors (FDs) and the Surface Detectors (SDs) of the Telescope Array Experiment using the data taken in our first 2.3-year observation from May 27 2008 to September 7 2010. A hybrid air shower reconstruction technique is employed to improve accuracies in determination of arrival directions and primary energies of cosmic rays using both FD and SD data. The energy spectrum presented here is in agreement with our previously published spectra and the HiRes results.