*Article* **Investigations of Muon Flux Variations Detected Using Veto Detectors of the Digital Gamma-rays Spectrometer**

**Krzysztof Gorzkiewicz 1,\*, Jerzy W. Mietelski 1, Zbigniew Ustrnul 2, Piotr Homola 1, Renata Kierepko 1, Ewa Nalichowska <sup>1</sup> and Kamil Brudecki <sup>1</sup>**


**Abstract:** This paper presents the results of cosmic ray muons flux monitoring registered by a digital gamma-ray spectrometer's active shield made of five large plastic scintillators. In traditional, i.e., analogue active shields working in anticoincidence mode with germanium detectors, the generated data are used only as a gating signal and are not stored. However, thanks to digital acquisition applied in designed novel gamma-ray spectrometers enabling offline studies, it has not only become possible to use generated data to reduce the germanium detector background (cosmic rays veto system) but also to initialize long-term monitoring of the muon flux intensity. Furthermore, various analyses methods prove the relevance of the acquired data. Fourier analyses revealed the presence of daily (24 h), near-monthly (27 days) and over bi-monthly (68 days) cycles.

**Keywords:** digital gamma-rays spectrometer; cosmic veto; active shield; muons; muon flux periodicity

#### **1. Introduction**

Low-background gamma-ray spectrometry is commonly used in research studies of materials characterized by trace concentrations of gamma-ray emitting radioisotopes. Hence, it finds applications in various fields of science, such as from neutrino physics to environmental research [1–3]. In such investigations, it is vital to use advanced shielding systems to reduce the background radiation of gamma-ray detectors (mainly HPGe).

One of the main components of background radiation is terrestrial gamma radiation, which spectrometers' passive shield can efficiently reduce. However, passive shield layers and their width must be appropriately selected in order to minimize the impact of internal gamma radiation from traces of radionuclides present in the shield's construction materials and any isotopes produced by interactions of cosmic rays with those materials. In lowbackground detection systems, contributions from air radioactivity, namely radon and its daughter isotopes, are not negligible [2,4].

Another important source of background radiation in gamma-ray spectrometers are particles of secondary cosmic rays. At sea level, the secondary particles flux consists of hadrons, neutrons, gamma quanta, electrons, muons, nucleons and antinucleons [5]. Of these charged particles, muons are the most abundant, with a mean energy of around 4 GeV. The intensity of the muon flux depends on the zenith angle *θ* of the incident particle, which at sea level can be expressed as follows (1):

$$I(\theta) = I(0^{\circ}) \cos^{n(p)}(\theta),\tag{1}$$

where *n*(*p*) is the particle momentum-dependent exponent and *n* ≈ 2 for muons with energies of a few GeV [6]. This relation explains the necessity of shielding gamma-ray

**Citation:** Gorzkiewicz, K.; Mietelski, J.W.; Ustrnul, Z.; Homola, P.; Kierepko, R.; Nalichowska, E.; Brudecki, K. Investigations of Muon Flux Variations Detected Using Veto Detectors of the Digital Gamma-rays Spectrometer. *Appl. Sci.* **2021**, *11*, 7916. https://doi.org/10.3390/ app11177916

Academic Editors: Roberta Sparvoli and Matteo Martucci

Received: 20 July 2021 Accepted: 25 August 2021 Published: 27 August 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

detectors in both vertical and horizontal directions. The most significant contribution to the radiation background is from particles with 0◦ zenith angle (vertical direction) [1].

Secondary cosmic rays passing through the shielding may deposit their energy in the germanium crystal (producing continuous background component) as well as generate neutrons and photons via several processes such as muon-induced hadronic and electromagnetic cascades, muon capture and muon-induced spallation reactions [2,7,8].

The cosmic rays background component can be reduced by using active shields consisting of detectors (plastic scintillators or multiwire Charpak chamber) surrounding a passive shield. In the case of particle detection (in the preset coincidence time window) by the active shield's detector and the germanium detector, the signal from the latter is not stored. This is the principle of the so-called classical (analogue) cosmic ray veto system.

However, the development of digital signal processing systems allowed one to apply digital analyzers (digitizers) as critical components of nuclear spectroscopy electronics and substitute a few discrete electronic devices used in analogue electronics. Such a device provides information about the registration time, energy and pulse shape of each significant signal generated by the detector. Furthermore, using digital analyzer allows all generated data to be stored for later processing (offline), enabling the application of various data exploration techniques. Since 2018, a low-background, gamma-ray spectrometer with an active shielding and digital acquisition system has been operating in the Department of Nuclear Physical Chemistry, at the Institute of Nuclear Physics Polish Academy of Sciences (IFJ PAN), in Krakow, Poland [1].

In this paper, we present the results of analyses of data generated by the spectrometer's active shielding. As mentioned above, these data are used to reduce the germanium detector background. Furthermore, offline data analysis makes it possible to develop a continuous cosmic-ray muons monitoring system. Hence one device, namely a low-background, digital gamma-ray spectrometer, can simultaneously performs experiments from two branches of physics, namely low-background gamma-ray spectrometry and astrophysics.

#### **2. Materials and Methods**

The spectrometer is equipped with a Broad Energy Germanium detector BE5030 (Canberra, USA) with relative efficiency of ≥48% and a composite passive shield in a cubic shape with an internal layer made of lead cast over 2500 years ago. The active shielding consists of five large, five cm thick plastic scintillation detectors EJ-200 with photo-multipliers ET 9900 (Scionix, Nl). These detectors are mounted outside the passive shield. TOP and BOTTOM detectors are placed horizontally while the latter three—FRONT, SIDE and REAR—are placed vertically. The relative positioning of all spectrometer's detectors is depicted in Figure 1. The detectors' preamplifier signals are transmitted directly to the inputs of a digitizer DT5725 (CAEN, Italy), where data acquisition and signal pre-processing are performed.

The digital analyzer DT5725 allows simultaneous acquisition of data generated by up to eight detectors (at this moment, in our setup six inputs are occupied) with a maximum time resolution of 4 ns. Registered data consist of pulse time registration, height (proportional to deposited particle energy) and shape. Raw data generated by the digitizer are stored in a PC as six files (in \*.csv format), which are further processed using purposely written software VETO. Commissioning and optimization processes of the described spectrometer and software development and its properties are discussed in detail in [1].

The data were collected from 1 September 2018 to 30 April 2020. Since the described spectrometer is primarily used to measure low-active gamma-ray emitting samples, the obtained data may be divided into two groups, namely short and long-period data. Short period data consist of data generated by scintillation detectors during single gamma-rays spectrometric measurement, which last up to 6 days. Data acquisition is stopped at the end of measurements in order to replace the sample in the spectrometer's chamber and/or refill the liquid nitrogen dewar. Such breaks last up to 30 min, after which the subsequent

gamma-ray spectrometric measurement along with registration of cosmic ray particles by scintillators is restarted. The long-period data consist of all digitizer output files generated from 1 September 2018 to 30 April 2020. During this time, 256 gamma-spectrometric measurements were carried out, and as a result 1280 output files of scintillation detectors were generated (total size 270 GB). In this paper, data generated by scintillators TOP (horizontal) and FRONT (vertical) were used to analyze the long-term modulations of muon flux.

**Figure 1.** The relative positioning of the low-background, digital gamma-ray spectrometer's detectors. In blue—horizontally (TOP and BOTTOM); in yellow—vertically placed scintillators (FRONT, SIDE and REAR). Each of the scintillators is 5 cm thick. The device is in operation at the Institute of Nuclear Physics Polish Academy of Sciences. Some parts of the passive shield (e.g., the lead walls) are hidden.

For this work, the atmospheric pressure data were obtained thanks to the cooperation with the Institute of Meteorology and Water Management National Research Institute. Data are collected at the Meteorological Station located at the Kraków-Balice Airport. The station is located 6 km west of IFJ PAN, and data were recorded hourly with a high accuracy of 0.1 hPa. The measurements contain pressure values at the station level of 237 m asl (above sea level). It is worth mentioning that the pressure values did not differ from the long-term data series. Only a clear diurnal course was found, which exceeded 1 hPa. The lowest pressure values are recorded in the afternoon, while the highest at night and in the morning. These hours coincide with the strongest and weakest convective movements in the atmosphere, respectively.

In the case of analyses of variations of cosmic rays intensity using detectors located at ground level, atmospheric pressure effects must be taken under consideration [9,10]. The influence of atmospheric pressure on the intensity of the cosmic ray flux is defined as the barometric effect, which includes the following components [11]:


Near the Earth's surface, the dominant factor is absorption one; considering that air density is usually the highest near the ground, the knowledge of the atmospheric pressure at the level at which the detection is performed is sufficient for determining the value of the barometric effect [12].

The normalized deviation of the recorded muon flux from the average depends on the change in atmospheric pressure [13]:

$$\frac{\Delta I}{I} = \beta\_P \Delta P\_\prime \tag{2}$$

where <sup>Δ</sup>*<sup>I</sup> <sup>I</sup>* is the normalized deviation of muon count rate, Δ*P* is the deviation of the atmospheric pressure from the mean and *β<sup>P</sup>* is the so-called barometric coefficient.

The *β<sup>P</sup>* coefficient (expressed in (%/hPa)) can be determined by assuming a linear correlation between changes in the normalized muon flux intensity and variations in atmospheric pressure. Therefore, it is essential to estimate the *β<sup>P</sup>* value only during the most geomagnetically quiet days [14].

Pressure corrected data are analyzed using various statistical tools. In order to investigate any correlation between subsequent registered events, the autocorrelation function (ACF) of a sequence of time intervals between pulses generated by detectors is determined. The autocorrelation function indicates the Pearson correlation coefficient between values of the same series as a function of time lag. Furthermore, Fast Fourier Transform (FFT) of a given time series is performed, allowing the investigation of any periodic components. Prior transformation, mean and linear trends have been subtracted from the sequence.

In the case of long-period data, a threshold normalization procedure was necessary. The spectrometer's configuration and optimization procedures allow a fixed threshold level for the scintillators spectra to 300th ADC channel to be set [1]. This procedure allowed redundant data (mainly registered gamma rays) which did not increase the effectiveness of the cosmic ray veto system to be reduced up to around 24 times However, to limit the influence of any signal threshold level fluctuations or gain changes in the digital acquisition circuit which may have occurred in data collected over a long period and cause uncontrolled changes in the recorded number of counts, normalization of the spectra discrimination level was carried out [14]. The normalization procedure involved cutting off part of the energetic spectra located below the ADC channel containing 30% of the maximum number of counts recorded in one channel in a given spectrum (i.e., the highest point in the energy spectrum, see Figure 2). The prepared output files were used to determine a time series of hourly muon count rates for the whole period considered.

**Figure 2.** Scintillator energy spectra discrimination procedure to eliminate fluctuations in the threshold level and signal gain. The black curve is the scintillation detector spectrum with a fixed threshold level set to the 300th ADC channel; the red dashed line represents 30% of the maximum number of counts registered in one channel in the given spectrum. The blue vertical line indicates the ADC channel defining the normalized discrimination level. The part of the spectrum below this channel is not used in further analyses.

#### **3. Results and Discussion**

#### *3.1. Short Period Data*

Distributions of pulses generated by the scintillators in unit time *t* should correspond to Poisson distributions with an expected value equal to *λt*. This value (for a chosen time interval) depends mainly on the area of the detector and its vertical or horizontal position. Figure 3 presents distributions of the numbers of counts registered in one second by the scintillation detectors TOP, BOTTOM and FRONT. In order to obtain experimental data (the cyan, red and pink lines in Figure 3), Poisson distributions were fitted with expected count rates *λTOP* = 71.93 counts/s (black), *λBOTTOM* = 18.88 counts/s (green) and *λFRONT* = 14.43 counts/s (blue) for TOP, BOTTOM and FRONT detectors, respectively. The horizontal detectors registered more particles in unit time than the FRONT detector, which was placed vertically. Additionally, Figure 3 demonstrates a great amount of data registered by a digital gamma-ray spectrometry system and supports the demand for the optimization process. Such factsare in accordance with previous research results described in [1]. Further data analyses involved the verification of the presence of correlated structures in detector signals. The gamma-ray measurement investigation lasted approximately 119 h.

**Figure 3.** Expected count rate distributions for scintillation detectors: TOP (71.93 counts/s), BOTTOM (18.88 counts/s) and FRONT (14.43 counts/s).

Figure 4a presents values of the ACF in the examined series in cases of up to ten events lag and from 10–1000 events lags. The ACF values oscillate around zero, indicating that there were no correlations in this signal. This means that practically every registered muon can be considered as a single independent event. The same results were achieved for other scintillators, implying that singular scintillators detect non-correlated cosmic-ray muons. Additionally, conducted analyses proved that signals generated by scintillation detectors are stochastic Poisson processes [15].

Figure 4b shows the obtained results from FFT analysis of the sequence of counts registered in one second by the TOP detector. A flat frequency spectrum (and, consequently, flat power spectrum) proves that data generated by the individual scintillation detectors during single gamma-ray spectrometric measurement correspond to white noise, and it is impossible to detect, using single scintillators, correlated muons created in the same air-shower event. This is because the time resolution of an active shield scintillator is mainly determined by two signal shaping parameters: the rise time of the trapezoid generated by the digitizer's energy filter and the length of its flat part [1]. Since the main purpose of the constructed active shield, muon detection, requires correct evaluation of the pulse height, the total shaping time is around 12 μs, much longer than the intervals between any possible registered cosmic-ray particles generated in the same air-showers.

**Figure 4.** (**a**) Autocorrelation function of time intervals between events registered at the TOP scintillator. As shown, there is no correlation in the examined series. (**b**) Discrete Fourier transform analysis of the sequence of counts registered by detector TOP in 1 s. No significant periodicities were found. Similar results were obtained for other detectors: BOTTOM, FRONT, REAR and SIDE.

Bearing in mind the previous results, coincidence events of signals from different scintillators have been studied. Coincidence events from four scintillation detectors, namely TOP, BOTTOM, FRONT and REAR, were found using the VETO software. The detection system registered 471 events of four-fold coincidences during the investigated gammaspectrometric measurement. Figure 5 presents the distribution of the coincidence events found as a function of the maximum time interval between pulses from the same event. Due to the non-linear geometric positioning of the considered scintillators, generated signals must come from at least two correlated muons from the same air-shower. Furthermore, the width of the distribution (σ = 12 ns) is comparable with the digitizer time resolution, and correlated muons time travel through the distance between scintillators. This feature is interesting from the novel CREDO project point of view [16], which focuses on studying cosmic rays and cosmic-ray ensembles. The presented detection system includes five plastic scintillators, which may play the role of a reference detector, that are complementary to other detection systems (e.g., smartphones' cameras) already applied in the CREDO project [16].

#### *3.2. Long-Period Data*

For long-period data, the barometric coefficient *β<sup>P</sup>* was established. During investigations, only data collected on the ten most quiet days of every month were used. The list of those days was acquired from the International Quiet Days (IQD) database (http: //wdc.kugi.kyoto-u.ac.jp/qddays/ accessed on 14 February 2021) [17]. Figure 6a presents relative deviation of the muon count rate <sup>Δ</sup>*<sup>I</sup> <sup>I</sup>* as a function of pressure deviation ΔP. The blue markers represent all data from the period studied, while the red markers represent data recorded only on IQD days. The correlation coefficient between these values and the pressure changes indicates strong negative correlation (*r*1,*TOP* = −0.75). The *β<sup>P</sup>* coefficient was determined by a least squares fit, and its value was *βP,TOP* = −0.168(2) %/hPa. The value was used to correct the data and to eliminate the dependence of the muon flux intensity on the atmospheric pressure using the Equation (2). The obtained results are presented in Figure 6b. The value of the correlation coefficient after data correction *r*2,*TOP* = −0.08 confirms a significant reduction in this dependence. In the case of the FRONT detector, analyses were performed to allow an estimate of *βP*,*FRONT* = −0.153(2) %/hPa, which reduced the data correlation from *r1,FRONT* = −0.72 to *r2,FRONT* = 0.09.

**Figure 5.** Distribution of 471 coincidence events registered by TOP, BOTTOM, FRONT and REAR detectors. To obtain the data, a Gaussian distribution has been fitted (μ = 32 ns; σ = 12 ns).

**Figure 6.** Determination of the barometric coefficient for the data recorded by the TOP detector. (**a**) Relative deviation of muon count rate as a function of pressure deviation for all data (blue markers) and for International Quiet Days (IQD) data—red markers. The correlation coefficient between the IQD data and the pressure variations was *r*<sup>1</sup> = −0.75. Obtained barometric coefficient *β<sup>P</sup>* = −0.168 %/hPa. (**b**) Pressure corrected data; correlation coefficient reduced to *r*<sup>2</sup> = −0.08.

Figure 7a shows the normalized intensity of the muon flux registered by the detector TOP before correction for changes in atmospheric pressure (which are depicted in Figure 7b). The pressure corrected data demonstrated in Figure 7c,d for scintillator TOP and FRONT, respectively, show less variability. It should be noted that seasonal variation is still clearly visible; however, seasonal changes do not affect the final results for the analyzed period (20 months).

**Figure 7.** (**a**) Normalized changes in the muon flux intensity recorded by the TOP detector before correcting for atmospheric pressure, (**b**) changes in atmospheric pressure, (**c**) corrected muon flux intensity Icorr registered by scintillator TOP and (**d**) corrected muon flux intensity Icorr registered by scintillator FRONT.

Fast Fourier transform analysis was conducted on the pressure-corrected data, and the obtained results are presented in Figure 8.

**Figure 8.** Discrete Fourier analysis of long-term muon intensity registered by scintillation detector TOP (**a**) and FRONT (**b**). Periodicities in time series 24 h (along with harmonics 12 h and 8 h), 27 d and 68 d have been found.

In the case of the TOP detector, the FFT algorithm allowed identification of the diurnal periodicity (24 h) with harmonics components (12 h and 8 h) of muon flux caused by the rotation of the Earth and its relative motion in the solar magnetic field, modulated by solar wind [18,19]. Furthermore, a periodic component with a period of around 27 days was identified. This variability is associated with the effect of the Sun's rotation, causing variations in both intensity of solar wind particles and their speed [18,20]. In addition, a periodicity of about 68 days is noticeable in the frequency spectrum but its origin is unclear. Takai et al. (2016) [21] conducted a frequency analysis of the eight-year time series of muon flux recorded during the MARIACHI experiment. The authors identified a signal component, among others, with a period of about 62.5 days (which is the closest

to the result obtained in presented research), but unfortunately its interpretation was not provided. Hence the 68day periodicity needs to be further analyzed.

The results of the FFT analyses of the signal generated by detector FRONT (Figure 8b) confirmed diurnal (with 12 h harmonics). In addition, there are 27 and 68 days periodicities in the horizontal component of cosmic ray flux registered by the vertical scintillator.

#### **4. Conclusions**

This paper presents the capabilities of a low-background digital gamma-ray spectrometer which expand its fields of applications and the analysis techniques to correct and extract information of the data generated by an active shield detector. Usually, such data are used only to reduce the germanium detector radiation background; thus, it is not saved or analyzed. Our approach, involving digital data acquisition and offline analysis, allows both tasks to be fulfilled simultaneously.

In the case of short period time series (i.e., data generated by a single scintillator during a single gamma-spectrometric measurement), research studies indicated that, according to expectations, the number of pulses generated by scintillators in unit time (in our case, 1 s) follows a Poisson distribution, and the time intervals between the pulses are not correlated. Moreover, by using the fast Fourier transform, the absence of periodic structures in these series was demonstrated.

Relatively long period data were defined as the hourly mean count rates of the registered cosmic rays. Investigations demonstrated strong negative correlations between the recorded muon flux and atmospheric pressure at ground level. The estimated barometric coefficients allowed data generated by detectors TOP and FRONT to be corrected for pressure. These corrected data were used in the analysis by using the FFT technique, and periodic components have been identified in both time series, including those related to the rotation of the Earth (with a period of 24 h) and the rotation of the Sun (a period of about 27 days). The periodicity with a 68 days period remains unexplained.

To summarize, by using a digital data acquisition system, it is possible to expand the research potential of the low-background gamma-ray spectrometer by numerous methods to explore collected measurement data and to allow monitoring of cosmic-ray muons flux registered by the active shield's detectors. Moreover, this device may find applications in the CREDO scientific project and other investigations focused on various phenomena correlated with the intensity of muon flux (e.g., earthquakes [22]).

**Author Contributions:** Conceptualization, J.W.M. and K.G.; methodology, K.G., J.W.M. and P.H.; software, K.G.; validation, K.G., J.W.M., Z.U., R.K. and P.H.; formal analysis, K.G.; investigation, K.G., J.W.M., R.K. and Z.U.; data curation, K.G. and Z.U.; writing—original draft preparation, K.G., J.W.M. and R.K.; writing—review and editing, K.G., J.W.M., R.K., Z.U., P.H., E.N. and K.B.; visualization, K.G., E.N. and K.B. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**

