Observing cosmic-ray extensive air showers with a silicon imaging detector

Extensive air showers induced from high-energy cosmic rays provide a window into understanding the most energetic phenomena in the universe. We present a new method for observing these showers using the silicon imaging detector Subaru Hyper Suprime-Cam (HSC). This method has the advantage of being able to measure individual secondary particles. When paired with a surface detector array, silicon imaging detectors like Subaru HSC will be useful for studying the properties of extensive air showers in detail. The following report outlines the first results of observing extensive air showers with Subaru HSC. The potential for reconstructing the incident direction of primary cosmic rays is demonstrated and possible interdisciplinary applications are discussed.

1 Cosmic-ray extensive air showers and surface detector arrays Cosmic rays are energetic particles from the universe discovered by V.F.Hess in 1912 [1].Their origin, acceleration mechanisms and mass composition at the highest energies are still largely unknown.In particular, determining the origin of high-energy cosmic rays is considered to be among the most important problems in modern astrophysics [2,3].The flux of cosmic rays arriving at Earth follows a power-law function of roughly E −3 , where E is the primary energy of the cosmic rays.As such, high-energy cosmic rays are extremely rare and are typically studied indirectly.This is accomplished by measuring the showers of secondary particles produced when primary cosmic rays interact with atmospheric nuclei, called "extensive air showers" (EASs).Several methods are used to observe EASs, such as the detection of fluorescence light [4,5] and radio emission [6][7][8] from the electromagnetic component of an EAS.Another method is the use of surface detector (SD) arrays -collections of ground based particle detectors spread across several hundreds or thousands of square kilometers [5,9,10].Such arrays measure the arrival times and densities of the secondary particles from an EAS at ground level.This information is used to reconstruct the arrival direction, primary energy and mass composition information of the initial cosmic ray.As an example, the Telescope Array experiment in Utah, USA, uses 507 plastic scintillators, deployed with a spacing of 1.2 km, to form a detecting area of ∼700 km 2 to cosmic rays with energies above 10 19 eV [9].By combining measurements from several detectors, the energy and arrival direction of the primary cosmic ray are determined.However, as is typical for surface detectors, these scintillators are unable to measure the properties of individual secondary particles during normal operation.Such information would be useful for identifying the primary particle and studying high-energy particle interactions [11,12].
To this end, this paper details the first results of utilizing the world-class silicon imaging detector Subaru Hyper Suprime-Cam to observe individual tracks of EAS particles with unprecedented statistics and precision.In doing so, we develop a new method of observing EASs.

Observing extensive air showers with Subaru
Hyper Suprime-Cam Subaru Hyper Suprime-Cam (HSC) [13] is the prime-focus wide field camera on Subaru Telescope.Located atop Maunakea in Hawaii, it is the largest silicon imaging detector above an altitude of 4200 m.Subaru HSC typically observes distant stars, galaxies and other interstellar objects in the optical and infrared.This is achieved by employing 112 Charge Coupled Devices (CCDs) each with dimensions of 4176×2048 pixels.Of these 112 CCDs, 104 are used for scientific observation and the remaining 8 are for calibration.The total area of the CCDs which perform observations is 0.196 m 2 .Standard operation of Subaru HSC involves taking images or "exposures" of targets of interest throughout a night.Typical exposure times for individual images are 150 or 200 seconds.During an exposure charged particles from the atmosphere may penetrate into the depletion layer of the CCDs.These particles leave long, thin "tracks" on the final image where the length of a track indicates the particle's angle of entry into the CCD.Normally, these tracks are sparse and randomly directed, and so are treated as unwanted noise which is removed.However, during the standard analysis of the Subaru HSC data, we encountered several images which displayed a far larger number of tracks than usual.Moreover, within each of these images, the majority of the tracks were aligned in a similar direction and were of similar lengths, indicating that the particles were travelling with similar trajectories.In turn we concluded that, for each of these images, the bulk of the detected particles had originated from a single, highenergy cosmic ray induced, EAS.We wish to emphasize the fact that the observation of these EASs was fortuitous and simply occurred during the normal operation of Subaru HSC.An example of one such image is shown in Figure 1.
Figure 1 contains several particle tracks not aligned with the general direction of the shower.These may be deflected particles from the same shower, or randomly directed particles originating from the constant background of low-energy cosmic ray showers.Understanding the background rate of these randomly directed particles is important as, based on how we initially found EAS events in our data, we expect the number of particle tracks in an image where an EAS was present to be significantly greater than the background rate.Thus, by comparing the expected number of "background tracks" in an image to the measured number, we can determine which exposures have observed an EAS.To achieve this, the following section constructs a model to describe the number of background tracks expected to be observed by Subaru HSC in different operating conditions.
Before continuing, it is important to note that with only one silicon imaging detector, as is the case with our current data, it is not possible to determine the energy, arrival direction or mass composition information of the original cosmic rays.However, with additional detectors spaced in a similar fashion to a regular SD array, these properties could in principle be reconstructed.Additionally, we expect that the detailed measurements of individual secondary particles by silicon imaging detectors could be used in conjunction with standard SD arrays as an effective method for studying EASs in detail.Both of these points are discussed further in Section 5.

Modeling the background rate of secondary particles
To effectively search for EASs among Subaru HSC data, it is necessary to formulate a model for the number of randomly directed, background electromagnetic and muonic particles expected to leave tracks on a Subaru HSC image, given various datataking conditions.We label this value N model .The various contributions to N model are described below, where the telescope elevation angle, telescope azimuth angle and sky brightness level are given by ψ tel , ϕ tel , and B sky respectively.Zero degrees azimuth was defined to be south, with azimuth increasing clockwise.

Effective area
We expect the number of detected tracks to be proportional to the effective detection area of the CCDs.This area itself is proportional to the absolute value of the cosine of angular separation between the telescope pointing direction and incoming cosmic ray direction.Given the cosmic ray azimuth angle, ϕ CR , and cosmic ray zenith angle, θ CR , the cosine of angular separation is given by: where α is angular separation.Assuming that the zenith angle dependency of the cosmic ray flux at Earth is cos 2 θ CR , the expected cosmic ray flux at telescope elevation ψ tel is proportional to the solid angle integral Therefore, the term (1 + sin 2 ψ tel )/2 is included in our model, where the normalization factor of 1/2 has been chosen so that the term becomes unity at ψ tel = 90 • (west).

East-west effect
The east-west effect is a suppression of low-energy positive cosmic rays arriving from the east due to the Earth's magnetic field [14][15][16].With our definition of azimuth, this effect can be modelled as proportional to the term sin(ϕ tel ).

Night sky background
We define the night sky background, B sky , to be the median number of ADC counts across the CCD array over the period of an exposure.This value changes with the background photon flux.Empirically, we have observed that higher values of B sky generally result in a larger number of tracks being detected.These additional tracks are believed to be "false" triggers of the track detection algorithm, not actual secondary particles.Since the exact dependence of B sky on the number of detected tracks is unknown, we assume a basic linear relationship in our model.
Adding the above effects, the final parameterisation of N model is where c 0 , c 1 , c 2 and c 3 are fitted parameters.Here, c 0 represents an offset to account for the different responses of the various filters used.
To find these parameters, we fit the above function to data collected by Subaru HSC between March 2014 and January 2020.Subaru labels one measurement of the sky with a specific filter and exposure time as a "visit".For this study, visits with an exposure time of either 150 or 200 seconds were used.Table 1 below summarises the data showing, for each filter, the measurement period, number of visits in each period, exposure time for a single visit, and the corresponding total exposure time (subtotal).Across all filters, this gave approximately 875 hours of data.
For each visit, dedicated software was used to count the number of tracks in the image by searching for different morphologies such as a long-straight, worm-like or wedge-like tracks.The tracks were identified based on a series of conditions, primarily considering the pixel intensity relative to the local background and gradient checks against the point spread functions of the astronomical objects in the image.The procedure is described in depth in Section 4.4 of [17].Pixels suspected to be contaminated by secondary particles were then grouped into contiguous events and further scrutinized in iterative passes.
Once the number of tracks for each visit had been calculated, the telescope's elevation and azimuthal angles, together with the night sky background levels, were extracted.Equation 3 was then fit to the data from each filter using the following method.Initially, the model was fit to all the data, including events containing EASs, using a χ 2 fit.As the objective was to model the number of background tracks we then removed data points with a number of tracks different from the model prediction by 2σ, where σ was the standard deviation of the mean squared errors between the measured data and fitted model.Equation 3 was then fit again to the remaining data and the best fit parameters recorded.We present an example of data and model fitting in Figure 2. The parameterisation results for all filters are shown in Table 2  The dependency on elevation angle (c 1 ) is roughly the same in every filter.Additionally, the azimuthal variation (c 2 ) shows that the number of cosmic ray tracks from the west is greater than that from the east, as expected from the east-west effect.The values for c 3 are comparatively small as typical values for B sky are on the order of 1000 counts.Finally, we note the large variation in c 0 across the different filters.In particular we see that, despite observing the same optical band, HSC-r has a significantly Fig. 2 A simplified example of fitting equation 3 to track numbers from different visits.The horizontal axis shows the telescope elevation angle.The vertical axis shows the number of detected tracks in a visit.The purple squares and green circles show data that was taken with the telescope pointing west and east respectively.There are several points above the main group of events.These points are likely to be EASs and are removed when applying the 2σ threshold (see text).The red line shows the model fit results, with ϕ tel = 0 and B sky = 1000.larger background level than HSC-r2.This may be because B270 contains potassium, a natural radio isotope.

EAS search results
By using the model defined in the previous section, we search for visits in the data-set with a significant excess of tracks compared to our model prediction.For demonstration purposes, we have defined a significant excess to be more than 20σ of each filter.Table 3 shows the significant visits.The telescope elevation angles for these visits, primarily between 50 • -70 • , is simply a consequence of those elevations being the most commonly observed.
To determine whether these visits actually observed an EAS, we checked for clustering in the particle arrival directions.Unfortunately, the true 3D trajectories of each particle are impossible to determine with the CCD image alone.This is because there is no method for determining which end of a track corresponds to the top face/bottom face of the CCD.Furthermore, we cannot tell which end corresponds to the particle entering/exiting the CCD.In other words, for a single track there are four possible trajectories, each of which is shown in Figure 3.If up-going particles are excluded (A and C in Figure 3), then each measured track has two possible directions (B and D in Figure 3).These directions differ by 180 • in azimuth in the frame of the camera   and have similar probabilities.For the sake of demonstration, we choose the trajectory with the smallest zenith angle.Regardless of whether this is the true direction, the observation of clustering will be evidence of an EAS.Choosing one direction in this way limits the effective detection area of the telescope to π sr.For each visit, the directions of the detected tracks were traced back to their position on the sky and plotted as 2D histograms.Two example histograms for visits 073808 and 161642, together with the positions of the raw tracks on the CCDs, are shown in Figure 4.The sky-map histograms have been divided into 12288 bins, with each bin equivalent to ∼ 0.001 sr.A clustering of particle arrival directions is clearly visible for both visits, indicating a common source for these particles i.e. a primary cosmic ray.All significant visits show some level of clustering.The open black circles in Figure 4 represent the telescope pointing direction, with these regions also containing a large number of entries.This is due to both the increased likelihood of detecting particles coming from the telescope pointing detection and an artifact of the analysis procedure which causes the track detection algorithm to misinterpret stars, tracks which cross each other, and coincidental bundles of high signal pixels as short particle tracks.

Discussion and future objectives
For a single detector, the clustering observed in our results is encouraging and is evidence of having observed EASs.With additional silicon imaging detectors, the arrival direction, energy and mass composition of the primary cosmic ray may be able to be reconstructed.The particle density at each detector could be directly calculated based on the number of tracks and related to a primary energy through a lateral distribution fit, whilst the clustering of the secondary particle trajectories from multiple detectors could lead to a reconstruction of the arrival direction.Assuming we are able to distinguish electrons/positrons from muons, which could be accomplished by separating curved (electron/positron) and straight (muon) tracks, mass composition information on an event-by-event basis may also be accessible through calculating the electromagnetic-muonic component ratio of the EAS.
Ultimately, we envision using silicon imaging detectors like Subaru HSC in conjunction with a SD array for highly detailed studies of EASs.With such a setup, the silicon imaging detector's ability to study EAS particle interactions in detail and determine event-by-event mass composition information could be combined with the precise knowledge of a showers energy and arrival detection from the SD array.We intend to test this concept by installing plastic scintillators inside the Subaru HSC building.This will allow us to look for coincident events between the two detectors.
Moving forward, we intend to cross-check the results in Table 3 with other telescopes at Maunakea Observatory, looking for time-coincident excesses evident of an EAS.If present, the particle densities at the different locations could be compared, similar to a SD array as described above, leading to a reconstruction of the primary energy.Possible future interdisciplinary applications include using detailed measurements of the vertex of muon decay inside the CCDs for searches of lepton flavor violating interactions such as µ + → e + +e − +e + , similar to the Mu3e experiment [18], and exotic signal searches for super heavy dark matter [19].

Conclusion
We have reported on the first results of detecting cosmic-ray extensive air showers with Subaru Hyper Suprime-Cam.Visits with a significant excess of secondary particle tracks in the CCDs were found and analysed for clustering of particle arrival directions.suggestions of data analysis, and thank Masato Yamanaka, Kohta Murase, Atsushi Naruko, Nagisa Hiroshima for valuable discussions regarding possible applications.

Fig. 1
Fig. 1 An example of a cosmic-ray extensive air shower recorded by a CCD of Subaru HSC.Dimensions : 4176 × 2048 pixels = 62.6 mm × 30.7 mm.The high density of aligned tracks, indicative of an EAS, is clearly visible in the upper right portion of the image.

Fig. 3 A
Fig. 3 A schematic view of the four possible incident directions of a charged particle penetrating the depletion layer of a CCD.The open arrows indicate the direction of charge transfer inside the depletion layer of the CCD.All four incident directions result in the same track being recorded.The up-going directions (A) and (C) are unlikely for cosmic rays (and hence their associated secondary particles).

Table 1
Summary of the Subaru HSC data used to search for cosmic-ray extensive air showers.

Table 2
together with σ. Results for the parameterisation of N model for each filter.
VisitDate UT Filter ψ tel ϕ tel B sky N track N model N excess

Table 3
Event information of the possible extensive air showers detected by Subaru HSC."HSC-" has been removed from the filter names for aesthetic purposes.