1 Introduction

Liquid noble gas detectors play an important role in rare-event search experiments looking for dark matter interactions or neutrinoless double beta decay [1]. One of their key features is ease of scalability. With larger target masses, external radioactivity can be better shielded through fiducialization. This is not the case for internal backgrounds that are intrinsic to the liquid gas target. First, these are medium- to long-lived radioisotopes of the target itself. For instance, liquid argon detectors need to take special care to avoid [2]. In the case of liquid xenon detectors, the two-neutrino double-beta emitter becomes relevant at the multi-ton scale [3]. Second, and more relevant for xenon detectors, are the radionuclides from radon and krypton. Both elements are inert gases that cannot be removed by established purification techniques based on hot gas purifiers commonly used in the field [4,5,6,7]. Radon and krypton dissolve in the liquid xenon target and cannot be excluded by standard fiducialization techniques which otherwise allow the rejection of background [5]. In this work we describe how krypton and radon backgrounds are assessed in the XENON100 experiment. Dark matter data from the three main XENON100 science runs (SRs), with exposure times of 101, 223 and 153 days each, are examined. The runs themselves and the corresponding detector conditions are outlined in Table 1 of [8].

2 The XENON100 detector

The XENON100 detector [5] is a cylindrical dual-phase time projection chamber (TPC) of \(30.5\,\hbox {cm}\) height and \(30.6\,\hbox {cm}\) diameter. It is located at the Laboratori Nazionali del Gran Sasso (LNGS) and uses about \(62\,\hbox {kg}\) of liquid xenon (LXe) as a target which is monitored by two arrays of photomultiplier tubes (PMTs), with one being at the top and one at the bottom of the TPC. Its primary goal is to search for dark matter in the form of weakly interacting massive particles (WIMPs).

Incoming particles are detected via their interactions with the LXe, generating xenon scintillation photons (S1 signal) as well as ionization electrons. The electrons then drift towards the top of the TPC due to a homogeneous drift field applied across the LXe volume. At the top of the TPC, the electrons are accelerated into a region of gaseous xenon (GXe) by an extraction field. Due to the moving electrons interacting with the GXe, proportional scintillation photons are created (S2 signal) [9]. Both S1 and S2 signals, measured in photoelectrons (PE), are detected by the PMT arrays. The delay between the S1 and the S2 signals of an interaction, combined with the hit pattern of the S2 signal on the top array, allows the reconstruction of all 3 coordinates of the interaction vertex. The X and Y coordinates are defined relative to the TPC’s central axis (where \(X = Y = 0\)) and are determined with a resolution of \(\sigma _{X/Y} < 3\,\hbox {mm}\). The Z coordinate has a resolution of \(\sigma _{Z} < 0.3\,\hbox {mm}\) and is defined with respect to the liquid gas interface at the top (\(Z = 0\)) and the cathode electrode, which is used to create the drift field, at the bottom of the TPC (\(Z = -30.5\,\hbox {cm}\)).

In addition, the S2/S1 ratio allows discrimination between nuclear recoils (NRs), which WIMPs are expected to induce, and electronic recoils (ERs), produced by \(\gamma \)-rays and \(\beta \)-particles. For example, in XENON100 WIMP analyses, \(99.75\,\%\) of ERs can be rejected at the price of an energy-dependent NR acceptance of 30–50% by utilizing this feature [10]. As tails of the ER distribution contaminate the NR region, it is of paramount importance to measure the abundance of and and to estimate their impact on the detector’s sensitivity to WIMP interactions.

3 Radon and thoron

The decays of (radon) and (thoron), as well as their daughters, are illustrated in Fig. 1 with the half-lives, \(\mathscr {Q}\)-values and branching ratios used throughout this work [11]. Radon and thoron are produced in the decay chains of the primordial nuclides and , respectively. Both of these nuclides are present, at least at trace level, in all materials, making it necessary to carefully screen and select all detector components [12]. The radon concentration of the air underground at LNGS, emanating from the surrounding rock, has been found to be of \(\mathscr {O}(100\,\hbox {Bq}/\hbox {m}^{3})\) [10]. For this reason, the inner cavity of the detector’s shield is continuously flushed with boil-off nitrogen [5], minimizing the amount of ambient radon and thoron that could potentially enter.

Levels of radon and thoron inside the LXe target of XENON100 are determined by the emanation of either isotope from surfaces inside the detector and the xenon purification system. Additionally, one period before SR1 and two periods during SR2 and SR3, were identified where air leaks of \(\mathscr {O}(10^{-3}\,\hbox {mbar 1}/\hbox {s})\) and \(\mathscr {O}(10^{-5}\,\hbox {mbar 1}/\hbox {s})\), respectively, developed at the purification system’s diaphragm pump, leading to a variation of the radon background over time (see [13] and Sec. 3.3).

Fig. 1
figure 1

Illustrations of the (radon) and (thoron) chains (ignoring decay modes with a branching ratio \(\le 0.1\,\%\)). Half-lives, branching ratios and \(\mathscr {Q}\)-values are taken from [11]. Solid boxes mark the isotopes that are quantified in this work

After entering the LXe, radon and thoron are able to reach the fiducial volume used for the WIMP search via diffusion and convection [14]. As a consequence, \(\beta \)-decaying daughter nuclides of both isotopes can contribute to the low-energy ER background. Contributions from \(\alpha \)-decays are not relevant because the involved \(\alpha \)-particle energies are two orders of magnitude larger than the energies expected from WIMP-induced NRs, which are at \(\mathscr {O}(10\,\hbox {keV})\) [10].

Some of the progeny of the chains’ nuclei have short half-lives compared to the event window of XENON100, which has a length of \(400\,{\upmu \,s}\) and is centered on the triggering signal [5]. This aspect results in two decays being recorded within the same event (delayed coincidence signature). An example for this are decays of (radon chain) and (thoron chain) which are followed by the decays of their polonium daughters (BiPo coincidence). This causes multiple S1 and S2 signals to be present in an event, making it possible to identify and reject them (see Sec. 3.2).

Of the \(\beta \)-decaying nuclides in either chain, many have a significant likelihood of decaying under prompt emission of \(\gamma \)-rays, which gives the same kind of signature as mentioned above. However, certain \(\beta \)-decaying nuclides from either chain are able to decay without \(\gamma \)-ray emission. If they are well-separated in time from accompanying decays, they are very likely to elude identification. Such nuclides are , (radon chain) and (thoron chain).

Radon and thoron concentrations in the LXe target can be inferred by selecting and counting events from decays of their chains. Especially suited for this task are \(\alpha \)-decays, as they produce large, monoenergetic signals resulting in a distinct event signature. Another explicit signature is the BiPo delayed coincidence as outlined above. This coincidence has already been successfully utilized by, for instance, the Borexino and SuperNEMO collaborations for estimating radioactive background levels inside their detectors [15, 16] and by the XENON collaboration for assessing the suitability of a thoron source for calibrating tonne-scale LXe detectors [14]. The focus of this section is to describe the selection of \(\alpha \)-decays and BiPo events. Results thereof are presented in Sec. 3.3.

3.1 Alpha event selection

For the analysis of \(\alpha \)-decaying nuclides in XENON100, and are used because they are the cleanest \(\alpha \) populations available as explained further below. is covered in Sec. 3.2 as part of the BiPo coincidence.

To select \(\alpha \)-decays, a set of cuts, based on criteria described in [10, 17], is applied to the data. We require at least one S1 signal with a minimum of two PMTs in coincidence. Any secondary S1 signal has to be below \(1600\,\hbox {PE}\) (motivated in [17]) to avoid decay pileup and multi-scatters. In addition, at least one S2 signal must be present with at least \(25\,\%\) of its area observed by the top PMT array to reject mis-identified signals. Due to the large signal sizes of \(\alpha \)-decays, which are found to be of \(\mathscr {O}(10^4\,\hbox {PE})\) for S1s and \(\mathscr {O}(10^5\,\hbox {PE})\) for S2s, acceptance losses due to the above mentioned cuts are considered to be negligible.

Finally, the detector volume used for this analysis is restricted to \(R = \sqrt{X^2+Y^2} < 135\,\hbox {mm}\) and \(10\,\hbox {mm}< Z < 260\,\hbox {mm}\) (\(\hat{=}\ 40.5\,\hbox {kg}\) LXe or \(65.3\,\%\) of the active volume). This excludes regions close to the TPC walls, which suffer from reduced light and charge collection efficiencies, and regions with insufficient separation between and .

A potential population close to the PTFE wall enclosing the TPC is also selected for further studies. Selection criteria are the same as above, with the following differences: \(R \ge 135\,\hbox {mm}\) is required, and the largest S2 signal must be smaller than \(8\times 10^4\,\hbox {PE}\). The latter criterion is motivated by the observation of S2 signals well below those seen from and for wall population events. Reduced S2 signals correspond to charge losses which can result from, for example, decays happening close to or within the walls. In the latter case, decay products can still enter the TPC, but lose energy in the process as they need to traverse the wall material. The S1 signal cannot be used as the only parameter for nuclide discrimination in this case, as the detector’s energy resolution is insufficient to separate the peaks of and in the S1 spectrum (see Fig. 2).

Fig. 2
figure 2

Corrected S1\(_{\alpha }\) spectra (SR1). The left peak of the regular \(\alpha \) selection spectrum (red, circular markers) is attributed to , the right one to , with their bounds (\(3\sigma \) intervals) indicated by dashed lines. The tailing peak in the wall selection spectrum (green, triangular markers) is considered to belong to and has its bounds indicated by dotted lines

The largest S1 and S2 signals are interpreted as belonging to an \(\alpha \)-decay and are correspondingly named S1\(_{\alpha }\) and S2\(_{\alpha }\). A position-dependent area correction for \(\alpha \)-decays [17] is applied to S1\(_{\alpha }\) in order to account for PMT saturation, which affects both the observed signal size and position reconstruction. Looking at the corrected S1\(_{\alpha }\) spectrum for events happening at \(R < 135\,\hbox {mm}\) (Fig. 2, red circular markers), we identify two peaks, which are attributed to the \(\alpha \)-decays of and , respectively.

To determine peak positions and extents while accounting for slight tailing, we fit a sum of two Crystal Ball functions (defined in equation (F-1) of [18]) and a constant to the peaks. Events are classified as containing a or decay if the area of their S1\(_{\alpha }\) is within \(3\sigma \) of the respective peak mean (bounds as shown in Fig. 2). This choice is valid as the peaks are, in good approximation, symmetric. Fits are done separately for each SR due to changes of detector parameters affecting positions and widths of both peaks [8]. In all SRs, the peak bounds determined according to this method do either not overlap, or overlap negligibly. In the latter case, the bound separating both peaks is determined by the arithmetic mean of the overlapping bounds in order to ensure events to be attributed to a single peak only. Leakage of the peaks beyond the boundaries assigned to them are estimated to be \(< 1\,\%\) and are thus considered negligible. For consistency, the same procedure is utilized for the wall population, as it also shows a peak in the S1 spectrum (Fig. 2, green triangular markers), using a single Crystal Ball function plus a constant for fitting.

In the thoron chain, the number of \(\alpha \)-decays from and can, in principle, be inferred from the S1\(_{\alpha }\) spectrum via peak fitting (see [19]). However, in the XENON100 background, the peak overlaps the peak region due to insufficient energy resolution. In addition, there is no indication of being present in the S1\(_{\alpha }\) spectrum of the fiducial volume used, while, at the same time, it is negligible in the rest of the sensitive volume compared to the wall population. As thus no direct evidence of them exists in the fiducial volume, and are not taken into account in this analysis, even though they are present in the detector as demonstrated by being measured, which belongs to the nuclei discussed in detail in the following section.

3.2 BiPo event selection

The decays of and are often recorded within the same event as the decays of their daughter nuclei, and . This is due to the short half-life of the polonium isotopes compared to the event window recorded by the XENON100 data acquisition system. The S1 signals generated by the \(\beta \) decays of the bismuth isotopes (S1\(_{\beta }\)) are smaller than those generated by the \(\alpha \)-decays of the polonium daughters (S1\(_{\alpha }\)), because the \(\beta \)-decay \(\mathscr {Q}\)-values (see Fig. 1) and ionization densities are lower than those of the \(\alpha \)-decays [20].

The result is a delayed coincidence signature of one S1 signal being followed by a larger one. For selecting such events, we require at least two S1 signals with at least twofold PMT coincidence and the correct time order (S1\(_{\beta }\) before S1\(_{\alpha }\)). Both signals need to be larger than \(200\,\hbox {PE}\) and S1\(_{\alpha }\) has to pass a data quality cut on the fraction of its area observed by the top PMT array to reject signals seen almost exclusively by the bottom array. Such a signal topology is virtually impossible to occur for an \(\alpha \)-decay happening inside the TPC due to the large amounts of scintillation photons generated (see Sec. 3.1).

In order to distinguish delayed coincidences from Bi and Po decays (called BiPos in the following) from either chain additional constraints are applied exploiting the fact that has a much shorter half-life than (\(T_{1/2} = 300\,\hbox {ns}\) vs. \(T_{1/2} = 162\,{\upmu }\hbox {s}\)). are selected by requiring S1\(_{\alpha }\) to occur at least \(7\,{\upmu }\hbox {s}\) after S1\(_{\beta }\), which removes more than \(99.99\,\%\) of events. For , the time difference has to be between \(0.5\,{\upmu }\hbox {s}\) and \(2\,{\upmu }\hbox {s}\). The lower bound ensures that both signals are individually identified with \(\sim 100\,\%\) efficiency by the data processor, while the upper bound removes about \(99\,\%\) of events.

Due to the possibility of \(\gamma \)-radiation accompanying the Bi-decays, the S2\(_{\alpha }\) signal falling outside the event window, and signal losses because of the spatial distribution of events as detailed in Sec. 3.3, no constraints are required on the number of S2 signals and their parameters. In fact, as the number of S2 signals is expected to vary and event reconstruction is not optimized for pairing S1 and S2 signals when multiple physical interactions are present, signal matching has to be done separately. A match requires the absolute time difference between a pair of S2 signals to be within \(\sim 1\,{\upmu }\hbox {s}\) of the one between the S1 signals (detailed in [21]). The S2 which occurs earlier is assigned to S1\(_{\beta }\). If no match is found, the largest S2 is assigned to S1\(_{\beta }\). We then recalculate positions and signal corrections (for S1\(_{\alpha }\) and those mentioned in [5]) for each event, as both depend on pairing S1 with S2 signals. Events without any S2 signal are not rejected, but are assumed to have occurred in a charge-insensitive region such as below the cathode, with a set of default coordinates assigned to them (\(R = 0\,\hbox {cm}\), \(Z = -30.5\,\hbox {cm}\)).

The data processor does not search for S1 signals occurring after a sufficiently large S2 signal within the same event [5]. This behavior is intentional as the processor has been developed for the analysis of single interaction events. However, this reduces the acceptance of the BiPo event selection because the S1\(_{\alpha }\) signal might occur after the first S2 signal of the Bi-decay which happens after S1\(_{\beta }\) within the maximum drift time of \(176\,{\upmu }\hbox {s}\) [5]. This loss in acceptance as well as the one resulting from the finite size of the event window is accounted for by summing up weights \(\varepsilon \) for each BiPo event, defined by

$$\begin{aligned} \varepsilon ^{-1} = \exp \left( -\lambda \ \varDelta t_{ min }\right) -\exp \left( -\lambda \ \varDelta t\right) . \end{aligned}$$
(1)

\(\lambda \) is the decay constant of the corresponding polonium isotope, \(\varDelta t\) is the time difference between S1\(_{\beta }\) and the first S2 peak (or the end of the event window, if no S2 is present), and \(\varDelta t_{ min }\) is the minimum time difference between S1\(_{\alpha }\) and S1\(_{\beta }\) allowed by the selection criteria. Thus, the right side of equation (1) is the probability of a polonium decay to occur within the given constraints in time.

Acceptance losses caused by the S1\(_{\beta }\) size criterion, however, cannot be reliably predicted without depending on measured and rates. This is caused by events that happen on the cathode, whose S1 signals are shadowed by the cathode grid which results in a modification of the expected S1\(_{\beta }\) spectrum. The relevance of cathode events is explained in the following section. Losses induced by other quality cuts are negligible.

3.3 Results and discussion

The rate evolution of each decay is shown in Fig. 3. To verify that the selected event populations represent the correct nuclei, an exponential decay plus a constant offset is fitted to the rate decrease of that is observed during the two months of SR1 (Fig. 3, top left). A small air leak was closed before this period, resulting in the decay of the excess radon which is visible in the rate evolution.

Fig. 3
figure 3

Radon event rates during the three SRs. Before starting SR1, entered the detector via an air leak which was subsequently closed. The decay of this additional radon contribution is clearly visible in the rates, and the decay constant is compatible with expectations from literature values (see text for more details). Note the rate correlation among , (selected by \(\alpha \)-counting) and (selected by delayed coincidence analysis) which are shown in the top row. The rates of (delayed coincidence analysis) and (wall event selection), shown in the bottom rows, do not correlate with them. Regions affected by a leak are shaded in blue. Mismatches between rates of the chain nuclides are explained in the text

Fig. 4
figure 4

Spatial distributions of the radon populations identified in XENON100. Dashed lines indicate the TPC’s radius. For and , the main S2 signal is required to be larger than \(80000\,\hbox {PE}\) to separate both from at \(R > 135\,\hbox {mm}\) (losses introduced by this cut are negligible at smaller radii). (Top) XY distributions. Due to vertex reconstruction artifacts caused by PMT saturation, the structure of the top PMT array is visible in the and data. (Bottom) \(R^2Z\) distributions. Only events whose reconstructed Z is within the TPC’s actual height are shown. The same reconstruction artifacts as for the top plots are visible in the and distributions. The number of events seems to diminish toward the top of the TPC, but this is actually an artifact induced by the peak finder peculiarities described in Sec. 3.2 and by multi-scatters induced by \(\gamma \)-rays which accompany the decay [17]

The half-life given by the fit is \(T_{1/2} = (3.81\pm 0.12)\,\hbox {d}\), which is in perfect agreement with the literature value for (\(T_{1/2} = 3.82\,\hbox {d}\)). In addition, the relative positions of the peaks in the S1\(_{\alpha }\) spectrum match the expectations given by the \(\mathscr {Q}\)-values of the individual decays, with a constant light yield of \(\sim 3.7\,\hbox {PE}/\hbox {keV}\) observed for all nuclides. Furthermore, the rates assigned to , and (radon chain) correlate with each other, while no correlation with the rates from (thoron chain) and (radon chain) can be seen. While the latter is also part of the radon chain, secular equilibrium is broken due to the long half-life of (\(T_{1/2} = 22.2\,\hbox {y}\)). While thoron can also enter the detector via leaks, it has a much shorter half-life (\(T_{1/2} = 55.8\,\hbox {s}\)) than , which results in a large suppression as it is more likely that it decays before reaching the TPC [14].

Table 1 Average specific rates for all SRs (statistical errors only). The leak period of SR1 ends on February 7, 2010, and the leak period of SR3 lasts from June 27, 2013, to December 1, 2013. Note that the rate concentration is large compared to the other nuclides because it is concentrated at the PTFE wall enclosing the TPC

The condition on the S2 peak size introduced to select decays originating from the TPC’s PTFE walls does not specifically select . However, considering that its rate is not correlated with the remainder of the radon chain, and taking into account similar observations made by the LUX experiment [22], we conclude that the wall population indeed consists of . The spatial distribution that includes it (see Fig. 4) shows that it is located almost exclusively at \(R^2 > 200\,\hbox {cm}^{2}\), while the largest fiducial volume used for XENON100 WIMP analyses requires \(R^2 < 200\,\hbox {cm}^{2}\) among other constraints [23]. For \(R^2 < 180\,\hbox {mm}^{2}\), a small number of events, likely caused by and leakage, can be seen. However, it is evident that these events are negligible compared to those happening at \(R^2 \ge 180\,\hbox {mm}^{2}\) as well as to those belonging to and .

Computing the average specific rates (Table 1) yields \((48.0\pm 0.4)\,{\upmu }\hbox {Bq}/\hbox {kg}\), \((64.3\pm 0.4)\,{\upmu }\hbox {Bq}/\hbox {kg}\) and \((68.3\pm 0.4)\,{\upmu }\hbox {Bq}/\hbox {kg}\) for in SR1 to SR3 respectively. Periods of increased average rates and fluctuations are observed in SR2 and SR3. These increases are caused by tiny air leaks in the diaphragm pump used in the xenon purification system, leading to a correlation of the rates inside and outside of the detector [13]. Restricting the rate average to periods not affected by a leak gives \((38.3\pm 0.4)\,{\upmu }\hbox {Bq}/\hbox {kg}\) (SR1) and \((41.8\pm 0.9)\,{\upmu }\hbox {Bq}/\hbox {kg}\) (SR3). We thus conclude that constant emanation of radon from detector materials results in a base rate of, on average, \(40\,{\upmu }\hbox {Bq}/\hbox {kg}\).

A direct measurement of the emanation at room temperature by means of miniaturized proportional counters was performed in summer 2012 between SR2 and SR3 [24]. It resulted in \((9.3\pm 1.0)\,\hbox {mBq}\) and \((2.6\pm 0.5)\,\hbox {mBq}\) being measured for the XENON100 detector and gas system, respectively, leading to an expected specific rate of \((74 \pm 7)\,{\upmu }\hbox {Bq}/\hbox {kg}\) assuming homogeneous mixing of in the full LXe inventory. Inside the TPC, the assumption of homogeneous mixing is valid, with the exception of (Fig. 4). The apparent decrease of events towards the top of the TPC is caused by losses induced by the peak finding algorithm as explained in Sec. 3.2 and by \(\gamma \)-rays, which accompany the decay, scattering off the LXe at a different position than the original decay. Measurements with an external source suggest, that the homogeneous admixing of radon throughout the entire LXe inventory takes place within a few hours [19]. The environmental conditions of the direct measurements with proportional counters differed from the standard operation conditions, as the detector and gas system were at different temperatures and exposed to nitrogen or helium, respectively. Both the increased temperature and reduced stopping power are known to impact the emanation rate of (for example, see [25, 26]) and we consider the direct measurement to be a weak confirmation of our results.

A priori, we expect the radon chain from to to be in secular equilibrium, as the longest-lived daughter nuclide in this part of the chain, , has a half-life of \(26.9\,\min \) (Fig. 1). This is short compared to both the time scales of the SRs, which lasted for several months (Fig. 3), and the time scale of the target purification, which is about \(5\) days per revolution. However, we observe only about \(50\,\%\) of the expected amount of events and about \(86\,\%\) of events (Table 1). Acceptance losses due to cuts are negligible for the \(\alpha \)-events from because of their high-energy signature. Thus, we have to consider additional causes for this mismatch. The most appealing one is radon daughters plating out onto the cathode due to convection and drift in the electric field. Radon daughters which remain ionized were, for example, observed in the EXO-200 TPC [27], and plating of radon progeny onto the cathode of a LXe TPC has already been reported by the ZEPLIN-III collaboration [28]. The precise motivation for the plate-out hypothesis is the observation of a surplus of events in the cathode region for and (Fig. 4), which is visible even when rejecting events without a proper S2 signal (which we assign to \(Z = -30.5\,\hbox {cm}\), the height of the cathode, by default). In addition, it has been observed in calibration data, that the drift field affects the motion of daughters inside the detector [14]. While nuclide velocities inside the TPC are dominated by convection, which contributes up to \(\sim 5\,\hbox {mm}/\hbox {s}\) to up-/downward motion along the Z axis, a constant contribution of \(\sim 1\,\hbox {mm}/\hbox {s}\) towards the cathode is observed which is attributed to the drift field (500–533\(\hbox {V}\) \(/\) \(\hbox {c}\) \(\hbox {m}\) depending on the SR [8]). As a consequence of the plate-out, decays happening on the cathode are shadowed, leading to losses in the S1 signal and thus a lower acceptance of BiPo and events.

No cathode accumulation is seen in the distribution. Such an effect could be hidden due to the reduced discrimination power at the bottom of the TPC between and (see Sec. 3.1). In addition, the effect on is assumed to be enhanced because of the repeated chance of collecting ionized daughters with every decay. A larger fraction of remaining ionized compared to , as suggested in [27], might also play a role.

rates are larger compared to those of other radon chain nuclides by a factor of \(\sim 4.2\) in the outermost part of the detector in periods not affected by a leak. However, one has to take into account that the volume within which is selected is by a factor of \(\sim 4.5\) smaller than a volume without any requirement on R (Sec. 3.1). Averaging the activity without constraining R gives, for instance, \((38.6\pm 0.3)\,{\upmu }\hbox {Bq/kg}\) in SR1. Because this rate is still larger than the one observed for and does not correlate with rates of preceding chain decays, we assume surface contamination of the PTFE walls due to air exposure during TPC assembly to be the origin of the population (analogous to observations made in [29]). Under this assumption, we find a activity per unit area of PTFE in the range from 0.6–0.9 \({\upmu }\) \(\hbox {Bq}\) \(/\) \(\hbox {cm}\) \(^{2}\).

4 Krypton

Natural krypton is present at the parts-per-billion (ppb) level in commercially available xenon produced in air separation plants. It contains the radioisotope , which is an almost pure beta emitter and has a relatively long half-life of 10.76 years. Krypton spreads throughout the liquid xenon target where it can induce low-energy events that may leak into the WIMP search region. To mitigate the -induced background, the xenon target typically is purified by means of adsorption or distillation before starting a measurement [30,31,32,33]. However, re-contamination due to even tiny air leaks readily increases the concentration of .

Natural sources result in a constant equilibrium content of 0.09 PBq in the atmosphere [34]. In addition, is produced alongside plutonium in spent nuclear fuel and irradiated breeding targets. The noble gas remains therein until it is released by nuclear fuel reprocessing facilities during the extraction of plutonium. These anthropogenic sources increase the atmospheric concentration of by orders of magnitude [35]. In present-day northern atmosphere, the activity of is approximately \(1.4\,{\hbox {Bq/m}}^{3}\) [36, 37]. This number roughly corresponds to a relative isotopic abundance of / = \(2\times 10^{-11}\,\hbox {mol/mol}\) [38]. However, the concentration varies across both time and space due to location and duty cycles of reprocessing plants as well as region-specific meteorological conditions [35]. To our knowledge, no atmospheric monitoring data is publicly available for the region around LNGS and for the relevant period of time. A single measurement (October 1, 2009) using miniaturized proportional counters of an air sample drawn underground close to the XENON100 detector exists, resulting in \((1.33 \pm 0.16)\,{\hbox {Bq/m}}^{3}\) [39], or / = \((2.11 \pm 0.25)\times 10^{-11}\,\hbox {mol/mol}\) in agreement with the expected average value.

In the following section we will discuss an in situ analysis technique to uniquely identify decays, quantify the krypton abundance during the investigated SRs, and compare the results to external measurements using a gas chromatographic system and a rare gas mass spectrometer (RGMS) [40].

4.1 Delayed coincidence analysis

disintegrates by \(\beta ^{-}\) emission to the ground state or, in \(0.438\,\%\) of all cases, to its second excited level. The half-life of the latter is \(1.015\,{\upmu }\hbox {s}\). This decay mode offers a unique feature for identification. In more than \(99\,\%\) of all cases the prompt \(\beta ^{-}\) emission with endpoint energy of \(173\,\hbox { keV}\) is followed by a single \(514\,\hbox { keV}\) gamma. This clear, delayed coincidence signature allows for an in situ analysis of krypton concentrations in the XENON100 detector despite the tiny branching ratio and low statistics. Energy levels, branching ratios and half lives are taken from [11].

A set of basic cuts is applied in order to reject electronic noise and to ensure data quality, closely following the procedure outlined in [10]. We require a twofold PMT coincidence level for both the largest and next-to-largest S1 signals, as well as a minimum width of both the S1 waveforms. In addition, no light must be seen by the PMTs observing the LXe volume outside of the TPC (veto volume) in coincidence with the two S1 signals. Background events due to increased electronic noise in SR2 and SR3 are also removed. Finally, we require that at least one S2 signal be identified in each recorded event trace.

delayed coincidence events are selected by requiring that the largest S1 (S1\(_\gamma \)) follows the next-to-largest S1 (S1\(_\beta \)) within a time window of 0.5–4.9 \({\upmu }\) \(\hbox {s}\). The acceptance of this criterion is \(67.5\,\%\). In addition, we demand that the reconstructed S1\(_\gamma \) and S1\(_\beta \) energies fall within amply defined energy ranges: for the gamma interaction this is three times the detector resolution (taken from [5]) around the expected value, i.e., from 330–698 keV. The maximal accepted S1\(_\beta \) energy is \(219\,\hbox { keV}\), i.e., the decay’s endpoint energy of \(173\,\hbox { keV}\) plus twice the detector’s resolution.

Detector-specific acceptance losses for small energy S1\(_\beta \) deposits are avoided by requiring signals to exceed \(14\,\hbox {PE}\), corresponding to \(5.8\,\hbox { keV}\) [41]. The acceptance of the latter condition is computed to be \(91.7\,\%\), using the \(\beta \)-Fermi-Function and the GEANT4 implementation thereof [42,43,44]. Poisson-like fluctuations in S1\(_\beta \) that affect the transformation from energy to S1 are negligible compared to the remaining uncertainties and ignored in the following. events originating from the decay chain (see Sec. 3) close to or on the PTFE wall that encloses the TPC constitute the background. These events are successfully removed by requiring that the sum of all identified S2 signals, i.e., the ones from the \(\beta \)-particle and the \(\gamma \)-particle, fall within the expected energy region of 514–687 keV. Conservatively, the region is enlarged by five times the S2 energy resolution.

We use the energy and interaction-type dependent S1 light yield and S2 gain to convert the energy ranges into S1 and S2 light signals. In SR1, the statistics in events is sufficiently high in order to determine both the S1 light yield at \(514\,\hbox { keV}_\gamma \) and the S2 gain at \(514\,\hbox { keV}_\gamma + 48\,\hbox { keV}_\beta \). The latter corresponds to the S2 sum signal of the monoenergetic \(\gamma \)-particle and the \(\beta \)-electron with an average energy of \(48\,\hbox { keV}\). For our purpose, we can assume the S1 light yield and S2 gain are constant throughout the three SRs investigated [8]. Finally, we use NEST [45], evaluated at \(0.5\,\hbox {kV}/\hbox {cm}\) similar to the XENON100 drift field, and the measured light yield at \(122\,\hbox { keV}\), to convert the upper bound of our acceptance window for \(\beta \) particles into an S1 value.

To convert the number of identified delayed coincidence events into a krypton concentration (always given in \(\hbox {mol}\) \(/\) \(\hbox {mol}\)), we have to account for the lifetime, the amount of xenon, the cut acceptances of \((61.7\pm 2.0)\,\%\) in total, and the relative isotopic abundance of . For simplicity, we assume \(2\times 10^{-11}\,\hbox {mol}/\hbox {mol}\) for the latter and resume the discussion in the context of induced uncertainties at the end of this section.

Fig. 5
figure 5

Observed event distributions for all three SRs. Displayed in black are events passing all cuts. Shown in red are those events that pass all selection cuts except for the condition on the S2 sum. See text for details. The red (34 kg) and blue (48 kg) contours indicate the two fiducial volumes used in the publications of the three SRs [8, 23, 47]. Dashed lines indicate the TPC radius

4.2 Results and discussion

Figure 5 shows the event distributions inside the TPC for the three SRs. Drawn in black are events passing all data selection criteria. Plotted in red are the events that pass all criteria except for the condition on the S2 sum. They are clearly clustered close to the PTFE wall enclosing the TPC, while the former (black events) are distributed throughout the TPC. For large radii, however, a reduced acceptance for events passing all selection criteria becomes obvious. We attribute this to a \(514\,\hbox { keV}\) \(\gamma \)-ray’s mean free path of roughly \(2\,\hbox {cm}\) in liquid xenon [46]. Close to the wall these \(\gamma \)-rays can exit the TPC undetected and we lose the characteristic pattern of the delayed coincidence. Fitting an exponential decay to the time delay between S1\(_\beta \) and S1\(_\gamma \), we find \(T_{1/2} = 1.08(_{-0.14}^{+0.18}) \, {\upmu }\hbox {s}\) and \(T_{1/2} = (0.32 \pm 0.04)\,{\upmu }\hbox {s}\) for the events passing and failing the S2 sum condition, respectively. This supports the hypothesis that the former (black events) indeed are due to , while the latter (red events) are caused by the delayed coincidence.

Table 2 Overview of the concentration measurements that were performed during the three SRs using mass spectrometry
Table 3 Result of the delayed coincidence study for the three SRs and considering three different fiducial volumes (FV). The number of tagged events is converted into a krypton concentration (DC). The concentrations can be compared to the corresponding SR-averaged off-line measurements (RGMS). See text for details
Table 4 Estimates for the average ER background induced by the chain and (below \(100\,\hbox {keV}\), before applying ER/NR discrimination). Values are inferred from different measurements. Only delayed coincidence values for the \(34\,\hbox {kg}\) fiducial volume are used, as they are affected the least by acceptance losses as outlined in Sec. 4.2

Table 2 lists all measurements that were performed off-line with the RGMS setup using gas samples drawn from the purification loop of XENON100. In this loop, about \(5\,\hbox {slpm}\) of xenon are continuously evaporated from the liquid xenon phase. Due to this large mass flow, we assume the Kr concentrations of these gaseous samples to represent the liquid xenon target. Employing the model describing the time evolution of the krypton concentration from [13] and using the available RGMS measurements, we compute the run-averaged krypton concentrations of \((340\pm 60)\,\hbox {ppt}\), \((11.0\pm 1.7)\,\hbox {ppt}\) and \((6.3\pm 1.0)\,\hbox {ppt}\) for SR1 to SR3, respectively. Along with the RGMS-derived concentrations, Table 3 lists the number of identified delayed coincidence events and the resulting krypton concentrations ( /Xe DC) found in the three SRs in the full TPC and two smaller fiducial volumes. As discussed above, the chance to miss a delayed coincidence event increases with radius due to \(514\,\hbox { keV}\) gammas escaping the TPC undetected. To account for this, we consider the innermost fiducial volume (\(34\,\hbox {kg}\)) only, where we find \(370_{-50}^{+60}\) ppt (SR1), \(15_{-7}^{+10}\) ppt (SR2) and \(18_{-9}^{+14}\) ppt (SR3), in good agreement with the average concentrations derived from the RGMS measurements. Limited statistics prevails in the uncertainties of the delayed coincidence method. For example, in SR2, we select only 5 events in \(7.6\,\hbox {td}\) of exposure. Comparing the measurements of SR1 where statistics is most favorable, we find that the gas samples drawn from the liquid in fact represent the entire xenon target. However, there is a small indication of higher concentrations from the delayed coincidence analysis in SR2 and SR3, i.e., we compute probabilities of \(0.058\) (\(0.30\)) for finding 4 (5) or more events in SR3 (SR2) based on the corresponding RGMS-estimated concentrations. This hints at a background gaining more importance with reduced krypton concentration and increased exposure times by, for example, random coincidences due to altered noise conditions [8]. Alternatively, underestimating the abundance of that entered the detector through the air leaks in SR2 and SR3 could cause the surplus in events observed in the delayed coincidence analysis. In contrast to SR1, where the krypton was introduced in a short period of time for which we have a direct measurement of the / ratio, plumes arriving from the two nearest reprocessing plants La Hague, France, and Sellafield, England, could alter the abundance of for SR2 and SR3. We account for such an effect by averaging the activity concentration in ambient air monitored by the German Federal Office for Radiation Protection (BfS) [48] at Mount Schauinsland close to Freiburg, Germany, during the relevant periods of time. We find correction factors of +10, +30 and +50% with respect to our initial assumption of / = \(2\times 10^{-11}\,\hbox {mol}/\hbox {mol}\) for SR1 to SR3, respectively. The resulting krypton concentrations are \(340_{-50}^{+60}\) ppt, \(12_{-5}^{+8}\) ppt and \(12_{-6}^{+9}\) ppt, increasing the probabilities for the observed number of delayed coincidences to 0.17 (0.50) in SR3 (SR2). We assume this estimate to serve as a conservative upper limit only. The distance from the dominant sources La Hague and Sellafield to the monitoring station at Mount Schauinsland is only half of the distance to the underground laboratories. Increased concentrations due to reprocessing cycles are supposed to be reduced at LNGS. In fact, simulations suggest variations in central Italy to be only on the order of \(0.5\,\hbox {Bq}/\hbox {m}^{3}\) [49]. Yet, for future experiments a local monitoring station is desirable to reduce this large systematic uncertainty.

5 Summary and conclusions

In this work, we presented techniques for selecting decays of the radon ( ) and thoron ( ) chains and those of . These methods allow us to estimate the contributions of the involved nuclei to the ER background and to study the distribution of background sources within the LXe target. Furthermore, they provide complementary values to those gained via direct measurements of the emanation rate and the concentration of natural krypton in the xenon target.

Background rates are given in units of differential rate (mDRU = \(1\times 10^{-3}\,\hbox {events}/(\hbox {kg}\,\hbox {day} \ \hbox {keV})\)). Monte Carlo studies [50], combined with the assumption of / = \(2\times 10^{-11}\,\hbox {mol}/\hbox {mol}\), yield conversion factors of \(0.029\,\hbox {mDRU}/({\upmu }\hbox {Bq/kg})\) and to relate and concentrations to ER rates, respectively. The chain has not been simulated due to and its daughters being more abundant as observed in data. The contribution of the chain can be estimated by using the observed rates. In the XENON100 science runs covered by this work, the xenon purity was affected by three air leaks of different leak rates. Consequently, we divide the radon-induced ER background into a constant pedestal driven by emanation and a variable offset due to leaking into the detector. In most parts of SR1, we do not observe a variable component due to external radon and infer the pedestal ER rate to be \(1.4\,\hbox {m\, DRU}\). In SR2 (SR3) we find the variable offset to account for \(0.5\,\hbox {m\, DRU}\) (\(0.6\,\hbox {m\, DRU}\)) on average. This corresponds to 35–40% of the total -induced ER background.

However, as we explained in Sec. 3.3, this overestimates the induced ER background because of plate-out effects. With being the most relevant \(\beta \)-emitter and ER background source of the chain [3], the actual ER background is smaller, assuming that the discrepancies between the decay rates arise mostly due to plate-out. We find the background index reduced to \(86\,\%\) (\(50\,\%\)) if we take ( ) rates to assess the effective activity concentration. Table 4 lists the radon contribution to the ER background using the different assumptions.

ER background from can be estimated via both RGMS and delayed coincidence measurements. While the RGMS measurements are more precise than the delayed coincidence measurements, we have to account for systematic uncertainties in the / ratio. Delayed coincidence measurements suffer from limited statistics, especially in SR2 and SR3, but constitute a direct measurement of the concentration which does not rely on any assumption for the krypton ratio.

Based on the analysis procedures detailed in this work, we can quantify the amount of air that entered through the leaks by means of the two tracers radon and krypton. In all cases, we find the -based estimate to be a factor of approximately two lower than the one from . From measurements with a spiked source, we know that within only two hours radon homogeneously admixes throughout the entire LXe inventory [19]. The agreement between in situ delayed coincidence and external RGMS measurements suggests that we do not miss a significant fraction of krypton in the liquid xenon target. To resolve the apparent tension, we conclude that krypton is enriched in the gaseous part of the detector beyond the expected value of \(\sim 10\) [51].

The contribution of in SR1 is \((14\pm 2)\,\hbox {m\, DRU}\) – one order of magnitude larger than . Krypton removal by cryogenic distillation results in being dominant in SR2 and SR3. For instance in SR2, and contribute \(29\,\%\) and \(11\,\%\) to the total ER background, respectively. This emphasizes the necessity to understand and control the target-intrinsic backgrounds radon and krypton and, in particular, it outlines the importance of radon screening and dedicated material selection campaigns complemented by online radon removal techniques for current and future liquid noble gas experiments like nEXO [52], DarkSide-20k [53], LZ [54], XENON1T, XENONnT and DARWIN [55].