The Telescope Array (TA) observatory utilizes fluorescence detectors and surface detectors (SDs) to observe air showers produced by ultra high energy cosmic rays in 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 yr using the Black Rock Mesa and Long Ridge fluorescence detectors in conjunction with the array of SDs. We compare the means and standard deviations of the observed X max distributions with Monte Carlo X 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 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 TA hybrid data at the 95% confidence level after some systematic X 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 this paper we report the first close, high‐resolution observations of downward‐directed terrestrial gamma‐ray flashes (TGFs) detected by the large‐area Telescope Array cosmic ray observatory, obtained in conjunction with broadband VHF interferometer and fast electric field change measurements of the parent discharge. The results show that the TGFs occur during strong initial breakdown pulses (IBPs) in the first few milliseconds of negative cloud‐to‐ground and low‐altitude intracloud flashes and that the IBPs are produced by a newly identified streamer‐based discharge process called fast negative breakdown. The observations indicate the relativistic runaway electron avalanches (RREAs) responsible for producing the TGFs are initiated by embedded spark‐like transient conducting events (TCEs) within the fast streamer system and potentially also by individual fast streamers themselves. The TCEs are inferred to be the cause of impulsive sub‐pulses that are characteristic features of classic IBP sferics. Additional development of the avalanches would be facilitated by the enhanced electric field ahead of the advancing front of the fast negative breakdown. In addition to showing the nature of IBPs and their enigmatic sub‐pulses, the observations also provide a possible explanation for the unsolved question of how the streamer to leader transition occurs during the initial negative breakdown, namely, as a result of strong currents flowing in the final stage of successive IBPs, extending backward through both the IBP itself and the negative streamer breakdown preceding the IBP.

Long γ-ray bursts (GRBs) are the most dramatic examples of massive stellar deaths, often associated with supernovae. They release ultra-relativistic jets, which produce non-thermal emission through synchrotron radiation as they interact with the surrounding medium. Here we report observations of the unusual GRB 101225A. Its γ-ray emission was exceptionally long-lived and was followed by a bright X-ray transient with a hot thermal component and an unusual optical counterpart. During the first 10 days, the optical emission evolved as an expanding, cooling black body, after which an additional component, consistent with a faint supernova, emerged. We estimate its redshift to be z = 0.33 by fitting the spectral-energy distribution and light curve of the optical emission with a GRB-supernova template. Deep optical observations may have revealed a faint, unresolved host galaxy. Our proposed progenitor is a merger of a helium star with a neutron star that underwent a common envelope phase, expelling its hydrogen envelope. The resulting explosion created a GRB-like jet which became thermalized by interacting with the dense, previously ejected material, thus creating the observed black body, until finally the emission from the supernova dominated. An alternative explanation is a minor body falling onto a neutron star in the Galaxy.

We report on a measurement of the cosmic ray composition by the Telescope Array Low-Energy Extension (TALE) air fluorescence detector (FD). By making use of the Cherenkov light signal in addition to air fluorescence light from cosmic ray (CR) induced extensive air showers, the TALE FD can measure the properties of the cosmic rays with energies as low as ∼ 2 PeV and exceeding 1 EeV. In this paper, we present results on the measurement of X max distributions of showers observed over this energy range. Data collected over a period of ∼ 4 years was analyzed for this study. The resulting X max distributions are compared to the Monte Carlo (MC) simulated data distributions for primary cosmic rays with varying composition and a 4-component fit is performed. The comparison and fit are performed for energy bins, of width 0.1 or 0.2 in log 10 (E/eV), spanning the full range of the measured energies. We also examine the mean X max value as a function of energy for cosmic rays with energies greater than 10 15.8 eV. Below 10 17.3 eV, the slope of the mean X max as a function of energy (the elongation rate) for the data is significantly smaller than that of all elements in the models, indicating that the composition is becoming heavier with energy in this energy range. This is consistent with a rigidity-dependent cutoff of events from galactic sources. Finally, an increase in the X max elongation rate is observed at energies just above 10 17 eV indicating another change in the cosmic rays composition.

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 yr period from 2008 May to 2017 May. 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 for the angular scale of the correlations. The energy threshold that 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σ and with the Auger result to within 1.4σ, meaning that it is not capable to discriminate between these two hypotheses.