Measurement of fast neutron background in the Gran Sasso underground laboratory using a geyser-concept bubble-chamber

A measurement of the fast neutron background in the Gran Sasso underground laboratory has been performed by the MOSCAB thermodynamic bubble chamber, a neutron detector which identifies fast neutrons passing through the metastable fluid by scattering ions. The result of this measurement confirms the results previously obtained using detectors based on totally different techniques, thus corroborating the previous results in terms of both integral intensity and spectral characteristics, at least in the energy region where the current measurement technique is sensitive.

Experiments dedicated to the search for very rare events require very low background rates. The main background source above ground at almost every energy is the cosmic radiation, which consists of electrons, muons, gammas and neutrons produced by the interaction of primary cosmic rays with the Earth atmosphere. Cosmic radiation can be attenuated or completely absorbed by the rock overburden of a deep underground laboratory. In this case, main external sources of background for a hosted detector will be cosmic ray muons surviving up to the laboratory depth, and the products of local radioactivity. Additionally, high energy muons can interact with the rock surrounding the laboratory giving rise to showers of secondaries and, in some extreme cases, atmospheric and solar neutrinos can also be considered as background sources. a e-mail: massimo.corcione@uniroma1.it b e-mail: walter.fulgione@lngs.infn.it (corresponding author) c e-mail: alfonso.santagata@enea.it An accurate knowledge of the background characteristics of an underground laboratory, in terms of both energy spectra and related intensities, is essential to design the appropriate shielding of any experimental equipment to be hosted; moreover, a continuous monitoring of the time variations of the intensity of the background sources during data taking is recommended, see e.g. [1].
Among all the sources of background, neutrons play an exceptional role. In fact, due to their capability to mimic the searched signals, neutrons are the most important background source for a wide variety of underground experiments, such as those aimed to the search for Weakly Interacting Massive Particles (WIMPs) in which case the seasonal variation of the neutron rate could also simulate the expected modulation of the dark matter flux.
In general, the neutron flux in underground laboratories is orders of magnitude lower than that above ground [2]. Consequently, its measurement requires the use of detectors with high efficiency, high signal discrimination capability and long-term stability. The MOSCAB geyser-concept bubble-chamber, which has been in operation in Hall C of the INFN Gran Sasso National Laboratory (LNGS) during the period 2018-2021, has proven to have these characteristics [3].
The INFN Gran Sasso National Laboratory, located at the average depth of 3800 m w.e. under the Gran Sasso massif, a limestone mountain in the center of Italy, can ensure particularly low background conditions to the hosted experiments. Notably, the neutron flux at the LNGS underground laboratory has been measured in the past by several research groups, in different locations and employing different tech-niques [2,[4][5][6][7][8][9][10][11]. A discussion on the low energy component of the neutron flux at the Gran Sasso National Laboratory can be found in [12], according to which the flux is dominated by neutrons produced in the concrete layer covering the cavity surfaces, thus implying that, at least for fast neutrons, it is practically the same in the three halls of the laboratory, their differences being mainly due to differences in the environmental humidity. Fission and (α, n) reactions contribute more or less equally to the total production rates both in the rock and in the concrete, but, while the spontaneous fission of 238 U, which dominates with respect to 235 U and 232 Th, mainly gives rise to neutrons with energies below 4 MeV, (α, n) reactions are mainly responsible of the production of neutrons with higher energies, up to about 8 MeV. At higher energies the contribution of muon-induced neutrons becomes dominant, even if their flux is definitely lower.
A set of measurements carrying information about the neutron energy spectrum have been executed at different sites inside the LNGS underground laboratory as reported in Fig.  1. Within the estimated uncertainties, in the energy interval 1-20 MeV the different data sets are in fair agreement, identifying a single energy spectrum.
The differential neutron flux dφ n /d E may be written as the product between the total neutron flux Φ n and the normalised energy distribution function g(E) which describes the shape of the spectrum: With this notation, the experimental data were fitted by the function where S 3,k (E) is the cubic spline with k breaks (or knots) [13] obtained minimizing the Residual Sum of Squares between the experimental values of the flux (F i ) and the corresponding estimated values in the n energy intervals The residual sum of squares tends to decrease with increasing the number of breaks, which, on the other hand, leads to an increase of the number of parameters (complexity) in the model. The trade-off between the goodness of fit and the simplicity of the model was performed applying the Akaike information criterion [14] which resulted in the optimal value of 3 knots. The uncertainties on the best fit function were calculated using the Monte Carlo method [15]. For each experimental value a log-normal probability density function was  Fig. 1.
In this paper we present the comparison between the expected event rate induced in the MOSCAB detector by the neutron energy spectrum shown in Fig. 1 and the experimental results obtained during a campaign of background measurement performed in Hall C of the Gran Sasso underground laboratory between 2020 and 2021.
The MOSCAB geyser-concept bubble-chamber, exhaustively described in a previous study [16], consists of a bubble chamber operated just controlling the temperatures of C 3 F 8 both in the metastable liquid phase and in the vapour phase. Its main strength is that it can be operated in such thermodynamic conditions to be insensitive to electron recoils, minimum ionizing particles and α particles, maintaining its sensitivity to nuclear recoils. In fact, the production of a vapour bubble in a metastable liquid requires that a given amount of energy is released in a well defined volume, which results in two conditions that must be simultaneously met: an energy threshold represented by the critical energy, E c , and a threshold in ionization power.
MOSCAB detector has been in operation in the LNGS underground laboratory from 2018 to 2021. During this period, thanks to the low background conditions existing inside the laboratory, we performed several calibration runs to validate the detector response functions computed for the detector operated in different conditions. Among them, the response function G(E c , E) of the detector irradiated by an isotropic flux of neutrons, which we will use here, are reported in Fig. 2.
The bubble event rate due to the neutron background in the Gran Sasso LaboratoryṄ ev may be computed by integrating the product between the response function G(E c , E) and the differential neutron flux dφ n /d E: where the function G B (E c ) evaluates the fraction of the neutron energy spectrum that can generate a bubble in the metastable liquid at the thermodynamic conditions corresponding to the critical energy E c . The results of this calculation, conducted for E c ranging from 75 to 550 keV, are reported in Fig. 3 with the error band produced by the propagation of the errors of the neutron spectrum of Fig. 1 and the response functions of Fig. 2. In the same figure, the experimental results obtained during 16 runs of data taking, for a total of 2586 h, with the detector operated at critical energies ranging from 96 to 514 keV, are also reported for comparison purposes. The related details are listed in Table 1. Vertical error bars in Fig. 3 represent the standard deviation of the distributions or, when larger, the Poisson fluctuations. Horizontal bar visualizes the standard deviation of the corresponding critical energy distributions.
It can be seen that less than 100 events have been detected during more than 100 days of live time of the detector, as in the cited range of critical energy the detector is sensitive only to neutron-induced recoils. The intrinsic background of the detector, which is dominated by the products of the Fig. 3 The background count rate of the MOSCAB detector placed in Hall C of the LNGS underground laboratory plotted versus the critical energy, E c and compared with the count rate expected from the environmental fast neutron flux Table 1 Results of the background campaign in Hall C of the LNGS underground laboratory in the period August 2020 -October 2021 The experimental data in the range 100 keV ≤ E c ≤ 150 keV shown by Fig. 3 support this interpretation and then, since we are aware that the background runs acquired at critical energies lower than about 150 keV are partly contaminated by signals due to the 222 Rn decay chain, in the analysis we introduced an additional energy cut at E c = 150 keV, which corresponds to an energy threshold E th n = 1.9 MeV for neutrons scattering with Carbon and to an energy threshold E th n = 2.7 MeV for neutrons scattering with Fluorine [3].
The remaining 26 events collected during more than 95 days of life time of the detector operated at 11 different metastability conditions have allowed us to obtain a measurement of the neutron flux intensity. In this regard, by keeping unchanged the shape of the neutron energy spectrum described by the energy distribution function g(E) defined by Eq. 1, we computed the optimal value of the normalization parameter Φ n as: where O(E c,i ) is the event rate of the i-th experimental run, According to our measurements, the neutron flux intensity integrated between 1 and 20 MeV turns out to be: 20 MeV 1 MeV dφ n d E d E = (0.38 ± 0.14) · 10 −6 n · cm −2 · s −1 which is in excellent agreement with the results reported in [12], i.e., 0.33 ± 0.07 ≤ Φ n · 10 6 ≤ 0.58 ± 0.13 n · cm −2 · s −1 . Moreover, the integration between 2.5 and 20 MeV results in a total neutron flux: 20 MeV which does not differ much from the neutron flux above 2.5 MeV of (0.9 ± 0.6) · 10 −7 n · cm −2 · s −1 measured by the GALLEX Collaboration [8] with a completely different technique.
Data availability statement This manuscript has no associated data or the data will not be deposited. [Authors' comment:The experimental data discussed in this paper are given in Table 1.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 . SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.