Domains @physik.uzh.ch
Gator: a low-background counting facility at the Gran Sasso Underground Laboratory
Gator: a low-background counting facility at the
Gran Sasso Underground Laboratory
L. Baudisa, A.D. Ferellaa, A. Askina, J. Angleb, E. Aprilec, T. Brucha, A.Kisha, M. Laubensteind, A. Manalaysaya, T. Marrodán Undagoitiaa and M. Schumanna
aPhysics Institute,
University of Zürich, Winterthurerstrasse 190, CH-8057, Zürich, Switzerland
bDepartment of Physics,
University of Florida, Gainesville, FL 32611, USA
cDepartment of Physics,
Columbia University, New York, NY 10027, USA
cGran Sasso National Laboratory,
Assergi, L'Aquila, 67010, Italy
E-mail: laura.baudis@physik.uzh.ch Corresponding author.
A low-background germanium spectrometer has been installed and is being operated in an ultra-low background shield (the Gator facility) at the Gran Sasso underground laboratory in Italy (LNGS). With an integrated rate of ∼0.16 events/min in the energy range between 100-2700 keV, the background is comparable to those of the world's most sensitive germanium detectors. After a detailed description of the facility, its background sources as well as the calibration and efficiency measurements are introduced. Two independent analysis methods are described and compared using examples from selected sample measurements. The Gator facility is used to screen materials for XENON, GERDA, and in the context of next-generation astroparticle physics facilities such as DARWIN.
Preprint typeset in JINST style - HYPER VERSION
L. Baudisa††thanks: Corresponding author., A.D. Ferellaa, A. Askina, J. Angleb, E. Aprilec, T. Brucha, A.Kisha, M. Laubensteind, A. Manalaysaya, T. Marrodán Undagoitiaa and M. Schumanna
E-mail: laura.baudis@physik.uzh.ch
\abstract@cs
Keywords: HPGe spectrometers; low level γ-ray spectrometry.
Gamma-ray spectroscopy offers a standard method for material screening and selection for rare-event searches, such as the direct detection of Weakly Interacting Massive Particles (WIMPs)  [2] or the search for the neutrinoless double beta decay  [3] . An ultra-low level germanium spectrometer in a dedicated low-background shield (the Gator facility) has been built and installed at the Laboratori Nazionali del Gran Sasso (LNGS), Italy. While the facility is being operated mainly in the context of the XENON program  [4] [5] [6] [7] , recently it also has been used for the GERDA experiment  [8] as well as for R&D purposes within low-background particle astrophysics.
This paper is organized as follows: in the next section, the Gator facility is described in detail. In the third section, calibration measurements, efficiency determinations for the different sample geometries, and cross-checks with standard, calibrated samples are described. In the fourth section, two different data analysis methods are introduced. The fifth section describes the determination of the main background sources of the facility as well as selected results for sample measurements. In the final section, a summary and discussion of future plans are presented.
2 The Gator facility
The design of the facility has been inspired by the layout of the world's most sensitive germanium spectrometers, operated at LNGS in connection with the Borexino and GERDA experiments  [9] [10] . The core consists of a high-purity, p-type coaxial germanium (HPGe) detector with 2.2 kg of sensitive mass with a relative efficiency of 100.5%111The quoted efficiency is defined relative to a 7.62 cm × 7.62 cm NaI(Tl) crystal, for the 1.33 MeV 60Co photo-peak, at a source-detector distance of 25 cm  [11] .. The detector construction has been performed in close cooperation with Canberra semiconductors  [12] , using only materials with ultra-low intrinsic radioactive contamination. The cryostat is made of ultra-low activity, oxygen-free copper with the cooling provided by a copper dipstick ('cold finger') in thermal contact with a liquid nitrogen bath. The cryostat is of U-type, with the cooling rod shaped in a right angle below the cryostat to avoid direct line-of-sight to the outside (see Figure 1). While the field-effect transistor (FET) used for the charge readout is cooled and close to the detector, the preamplifier is placed outside the low-background shield, since it contains more radioactive components.
The shield of the detector has been designed to provide a large sample capacity, an ultra-low background and easy access to the germanium spectrometer itself. The sample chamber, with a dimension of 25×25×33 cm3, is surrounded by 5 cm (7 cm for the base plate) of oxygen-free, radio-pure copper from Norddeutsche Affinerie  [13] . Residual surface contaminations of the copper plates were removed by treating them with diluted sulfuric and citric acid solutions, followed by cleaning with deionized water. All steps were performed under clean-room conditions. The copper is surrounded by four layers of lead from Plombum  [14] , each 5 cm thick. The inner 5 cm lead layer has a nominal 210Pb activity of 3 Bq/kg, while the outer 3 layers were built from lead with a higher 210Pb activity of 75 Bq/kg. All lead bricks were cleaned with ethanol before being installed in the shield. Their arrangement is such that no direct line-of-sight to the HPGe detector is possible. Large copper plates, which close the sample cavity and carry the upper lead bricks, were placed on sliding rails allowing for easy opening and closing of the shield: The full sample chamber can be accessed when the top shield is open. The lead is surrounded by 5 cm of polyethylene as a shield against ambient neutrons, and the entire shield is enclosed in an airtight aluminum box in order to prevent radon entering the system.
Figure 1: Schematic figure of the Gator facility at LNGS. The HPGe detector (a) with its cold finger (b) and dewar (c), the open sample chamber, the copper (d) and lead (e) shield with the sliding door, the glove box (f) and polyethylene shield (g) can be seen.
Figure 1 shows a schematic view of the detector and shield configuration. A sample handling and glove box made of plexiglass, including an airlock system, is placed on top of the aluminum housing. The entire system is continuously flushed at slight over-pressure with boil-off nitrogen gas to suppress radon diffusion into the shield. A 5 mm inner diameter PTFE tube allows sealed calibration sources to be brought close to the germanium detector. The samples to be screened are first cleaned in an ultra-sonic bath of ethanol (where applicable), then enclosed in a sealed plastic bag for transportation to the underground site, after which they are stored for a few days under nitrogen atmosphere above the closed chamber. This allows trapped radon and plate-out 220Rn and 222Rn progenies to decay before the actual counting starts.
To remotely check the stability of the detector with time, a monitor system has been installed. The liquid nitrogen level, the flow of nitrogen gas, the leakage current across the germanium diode and the high-voltage applied to the diode are read every 2 minutes, while the overall trigger rate is read every 6 hours. The liquid nitrogen level is measured with a level-meter consisting of two 40 cm long concentric metal tubes acting as a capacitor whose capacitance is read out with a universal transducer interface board. The nitrogen gas flow into the shield is monitored and regulated with an electronic flowmeter and the leakage current is measured as the voltage drop across the feed-back resistor and read out with an analog-to-digital converter. The relevant parameters are plotted versus time and can be monitored on a web interface which is refreshed every 10 minutes. In case pre-defined thresholds for these parameters are not met, email and SMS alarms are being issued. The data acquisition consists of a high-voltage unit, a spectroscopy amplifier, and a a self triggering, 16 k channel multichannel analyzer. The spectra, acquired and cleared every 6 hours, are saved in ascii files and analyzed offline.
Prior to its installation at LNGS, the germanium spectrometer was operated in the Soudan underground laboratory in northern Minnesota, within the SOLO facility  [15] . At the Soudan laboratory, the detector was used for the screening of XENON10  [4] materials, and several background runs were acquired  [16] . During its water- and ground-based transportation to LNGS, the detector was exposed for several months to the cosmic ray flux at the Earth's surface, leading to cosmic activation of the crystal and the surrounding copper of the cryostat. The shield was improved in October 2008, and it was cleaned once again in February 2009.
Table 1 shows the integral background in the 100-2700 keV region in the Soudan configuration, and for three measurements taken at LNGS. The integral counting rate at LNGS has been constantly decreasing due the decay of cosmogenic radio-nuclides such as 54Mn, 57Co, 58Co and 65Zn, with typical half-lives around one year on the one hand, and due to an improved shield and overall sealing of the system on the other hand. This is reflected in the decrease of prominent lines from 214Bi and 214Pb, as shown in Table 2. It gives the integral counting rates in the ±3σ-regions for the main primordial γ-lines, as well as for the 137Cs, 60Co and 40K lines, along with a comparison with the GeMPI detector, which is one of the world's most sensitive low-background spectrometers  [9] .
Lifetime [days]
Rate [events/min]
Gator at Soudan
22.96 0.842 ± 0.005
Gator at LNGS (09-2007)
Table 1: Integral background counting rates for Gator as measured at Soudan and at LNGS in three different runs. The integral is evaluated in the energy range [100, 2700] keV.
Chain/nuclide
Peak integral background rate [counts/day]
[keV]
Gator (Soudan)
Gator (LNGS, 09-2007)
GeMPI [9]
239 232Th/212Pb 1.1 ± 0.7 0.7 ± 0.1 0.13 ± 0.08 <0.5 NA
911 232Th/228Ac 0.9 ± 0.3 0.4 ± 0.2 0.4 ± 0.1 <0.5 <0.2
352 238U/214Pb 4.9 ± 0.7 4.3 ± 0.7 1.1 ± 0.2 0.7 ± 0.3 <0.5
609 238U/214Bi 4.5 ± 0.5 4.0 ± 0.5 1.1 ± 0.2 0.6 ± 0.2 0.50±0.45
1120 238U/214Bi 1.6 ± 0.3 2.7 ± 0.4 1.3 ± 0.2 0.3 ± 0.1 NA
1765 238U/214Bi 1.3 ± 0.2 1.5 ± 0.3 0.2 ± 0.1 0.08 ± 0.06 NA
662 137Cs 2.9 ± 0.4 0.5 ± 0.3 <0.5 0.3 ± 0.1 NA
1173 60Co 0.5 ± 0.1 0.5 ± 0.2 0.5 ± 0.2 0.5 ± 0.1 0.6±0.4
1461 40K 5.8 ± 0.4 0.5 ± 0.2 0.4 ± 0.1 0.5 ± 0.1 0.6±0.4
2615 208Tl 0.7 ± 0.1 0.2 ± 0.1 0.2 ± 0.1 0.2 ± 0.1 NA
Table 2: Background counting rates (in events/day) in the ±3σ-regions for the main primordial and the gamma lines of 137Cs, 60Co and 40K.
Figure 2 (top) shows the comparison between the latest background spectrum acquired at LNGS (2010), a spectrum taken in the SOLO facility at Soudan (2007) and the background of the GeMPI detector [9] . It also shows (bottom) the Gator spectrum underground at LNGS prior to its installation in the shield, inside the shield, and inside the shield with the radon protection system on. The background decrease is more than four orders of magnitude, and shows that a careful shielding is needed even when the detector is operated deep underground.
Figure 2: (Top) Background spectra of Gator at Soudan (red), at LNGS (black) and the spectrum of the GeMPI detector [9] (blue). (Bottom) Gator background spectrum at LNGS: outside the shield (black), inside the shield (blue) and inside the shield with the radon protection system on (red), clearly showing the suppression of the main gamma lines associated with radon decays.
3 Calibration measurements and efficiency determination
The HPGe detector is calibrated regularly with radioactive sources such as 109Cd, 133Ba, 137Cs, 60Co, 57Co, 22Na, 54Mn and 228Th. In Figure 3 (left) the comparison of the spectrum obtained from a 60Co calibration with the one from a Monte Carlo simulation of the source-detector geometry is shown. The FWHM of the two 60Co lines at 1173 keV and 1332 keV are 2.5 keV and 3.0 keV, respectively. The energy resolution of the detector, defined here as the ratio of the σ to the mean energy of the gamma line, is shown in Figure 3 (right).
Figure 3: (Left) Comparison of the Gator spectrum from a 60Co calibration (black, data points) with the one from a Monte Carlo simulation of the source-detector geometry (red, solid). (Right) Energy resolution (defined here as the ratio of σ/energy) as a function of energy. The solid curves represent a fit using the function σ2(E)=E2(2.35×10−7)+E(7.70×10−4)+(4.43×10−1).
The calculation of the detection efficiencies of the various gamma lines used in the analysis of the experimental data (see Section 4) is based on Monte Carlo simulations using Geant4  [17] . For each measured sample, a detailed geometry is included into a Geant4 model of the facility. The efficiency ϵ of a specific γ-line is defined as the ratio between the number of events detected in the line to the number of gammas of that energy emitted by the source. In order to simulate each decay chain, the G4RadioactiveDecay class, which takes into account the branching ratios for the different gamma lines in one decay, is used.
To cross-check the efficiency determination, a measurement of two extended sources and a comparison with the certified values for their activities has been performed. The sources used for this measurements, which had similar dimensions and weights, are CANMET-STSD2 (from the Canada Centre for Mineral and Energy Technology) and IAEA-Soil6 (from the International Atomic Energy Agency). Both sources are soils from different places on Earth, which have been thermally treated and sieved several times in order to destroy any remaining organic matter. The homogeneity of the material is certified by the provider. Tables 3 shows the results Gator's screening of these two sources and a comparison with the certified values. Within the uncertainties of the measurements, we find a good agreement. We also show in Figure 4 the comparison of the efficiencies as determined by the Monte Carlo method, with those from the data, as a function of energy, as well as the relative difference among these.
These results indicate that the measurements performed with our spectrometer provide a reliable value for the activity of a given sample. Further cross-check were done by using well-calibrated, commercially available point sources.
Nuclide
STSD2 Activity [Bq/kg]
Gator results Certified values
75±4 70±5
226Ra
230±30 230±10
Soil6 Activity [Bq/kg]
Table 3: Results from screening the CANMET-STSD2 (493 g) and IAEA Soil6 (530 g) sources and comparison with certified values provided by two agencies (for details, see text).
Figure 4: Upper panel: Efficiencies as a function of energy for the two certified sources, as determined from data (filled symbols) and Monte Carlo (open symbols). Lower panel: relative difference between simulated (MC) and measured (D) efficiencies.
4 Data analysis methods
Two methods are used to determine the concentrations of radioactive nuclides in a given sample. In the first method, the most prominent γ-lines are analyzed, using efficiencies as determined by a full Monte Carlo simulation. In the case of 238U, the γ-lines from the daughters of 226Ra (214Pb and 214Bi) and in case of 232Th the γ-lines from 228Ac and from 212Pb, 212Bi and 208Tl are used. In the second method, the overall data spectrum is compared to the one obtained from a Monte Carlo simulation, after subtracting the measured background spectrum, and the activities are determined from the best fit. The results from the two methods agree within the statistical errors, as shown in Section 5. Both methods are explained in more detail below.
4.1 Analysis of gamma-lines
The first method is based on counting events at the location of the most prominent lines, after subtracting the background spectrum closest in time. The Compton background, estimated from the regions left and right of a peak, is subtracted as well. The decision on whether the signal exceeds the background is based on comparing the net signal number of counts
Snet=S−B⋅tS/tB−BC (4.0)
with the so-called detection limit Ld (the level of a true net signal that can be detected with a given probability) as defined in [19] for a ∼95% C.L.:
\specialhtml:<Aname="ref−ld">\specialhtml:</A>Ld=2.86+4.78√BC+B⋅tStB+1.36. (4.0)
S is the number of counts in the ±3σ-region around a peak, B and BC are the number of background and Compton-background counts in the same region, and tS, tB are the measuring times for signal and background, respectively. For each peak, three cases are considered  [19] :
Snet<0: the upper limit is set to Ld (no net contribution from a signal)
0<Snet<Ld: the upper limit is set to Snet+Ld (there is an indication of a signal, but it can not be confirmed for the existing background level and sample exposure)
Snet>Ld: the detection limit is exceeded (clear indication for a signal at 95% C.L.)
For the third case, the specific activity and its 1 σ error is calculated as
\specialhtml:<Aname="ref−eq:activity">\specialhtml:</A>A[Bq/kg]=Snetr⋅ϵ⋅m⋅twithΔAA=ΔSnetSnet. (4.0)
with the peak efficiencies ϵ as determined by Monte Carlo simulations, the branching ratio r for the specific line, the mass m of the sample (in kg), and the measuring time t (in seconds). For the case in which an upper limit is reported, Snet is replaced by Ld or by (Snet+Ld) in equation (4.1). As concrete examples, we show the above quantities in Table 4, together with the determined specific activities or upper limits for a copper and a stainless steel sample, using different gamma lines from their measured spectrum.
Used line [keV]
B⋅tS/tB
Activity [mBq/kg]
239 0 0 93 49 Snet<Ld <0.33
352 -5 19 71 48 Snet<Ld <0.36
1173 42 6 19 27 Snet>Ld 0.24±0.06
352 66 7 58 42 Snet>Ld 4.3±0.9
1173 236 2 19 25 Snet>Ld 7.2±0.9
1461 -3 3 10 21 Snet<Ld <5.7
Table 4: Examples of upper limit or specific activities calculation. Details are given in the text.
4.2 Fit of the data to a simulated spectrum
The individual activities are also determined by a global fit based on a χ2-minimization of the simulated spectrum to the experimental data. The aim is to model →y, the measured number of counts in each energy bin, as a function of →x, the energy bin value from the Monte Carlo simulation, with a functional dependence of the form
→y=M∑k=0ak⋅fk(→x),\specialhtml:<Aname="ref−eq:chi2">\specialhtml:</A> (4.0)
where M is the number of simulated radioactive isotopes. ak≥0 are the scaling factors for each isotope, which are kept as free parameters, and fk(→x) are M Monte Carlo spectra which are already smeared with the energy resolution of the detector. For uncorrelated statistical errors in each energy bin i, the χ2 value is defined as
χ2=N∑i=0[yi−∑Mk=0ak⋅fk(xi)]2σ2i,\specialhtml:<Aname="ref−eq:chi2">\specialhtml:</A> (4.0)
where σi=√yi is the variance in the observed number of counts in each bin and N is the number of energy bins over which the fit is performed. The statistical uncertainty in the Monte Carlo component is negligible. Before the analysis, the measured background spectrum is subtracted from the sample spectrum. The outcome of this procedure are the scaling factors ak±Δak for every decay chain/isotope k. In the subsequent analysis, these factors are assumed to describe a normalized Gaussian probability distribution function (PDF) with mean a0=ak and σ=Δak.
The method applied to decide between a real detection and an upper limit is based on [18] . If the lower of the two symmetric limits providing 95% statistical coverage is positive, a peak detection is claimed and the activity is calculated from ak, taking into account the measuring time t in seconds, the sample mass m in kg, and the number of simulated events nsim. The relative ±1σ error is given by Δak/ak:
\specialhtml:<Aname="ref−eq:limitschi2">\specialhtml:</A>Ak[Bq/kg]=nsim akm twithΔAkAk=Δakak. (4.0)
If the lower limit is equal or less than zero, no detection can be claimed at 95% C.L. and an upper limit is given. Its value aup is determined by the 95% quantile of the positive part of the Gaussian PDF defined by a0=ak and σ=Δak. The limit on the activity (95% C.L.) is calculated using aup in equation (4.2).
As an example, Figure 5 shows a measurement of a stainless steel sample: the data spectrum is compared with the Monte Carlo sum spectrum, and the individual contributions from 238U, 232Th, 40K, and 60Co, as given by the best-fit, are shown. The derived activities and upper limits for this sample are given in Table 6.
Figure 5: Best-fit of the Monte Carlo simulations to the measured data for a stainless steel sample. Data is shown with error bars (red), the full spectrum as determined by the simulation is shown as the grey solid curve. The individual contributions from the decays of 60Co (blue), 238U (green), 232Th (magenta), 54Mn (cyan) and 58Co (brown) are also shown, along with the measured background spectrum (black).
In this section, we first outline results obtained from a detailed study of the background of the facility, using the analysis method described in Section 4.2. We then present screening results for a few selected samples, the goal being to compare the outcome of the two analysis methods introduced in the previous section. Results from a much larger selection of screened samples are given and discussed in detail in [21] .
5.1 Background Analysis
To model the residual background of the Gator facility, the detailed geometry of the crystal, cryostat system and shields has been simulated with Geant4. The following potential contributions to the background have been simulated: the natural decay chains of 238U, 232Th and 40K decays in the copper of the shield and of the cryostat, the decays of the cosmogenic radio-nuclides 54Mn, 65Zn,58Co, 57Co, 60Co in the Ge crystal, the Cu of the cryostat and of the shield, 210Pb decays in the innermost Pb layer, and 222Rn decays inside the shield.222A background run up to an energy of ∼8 MeV was acquired for a few weeks and revealed no U/Th contaminations close to the crystal, which would be visible as alpha peaks above a few MeV. Other materials in the cryostat (PTFE, mylar and kapton) were neglected in the simulation due to their much lower mass compared to the copper (a few tens of grams and milligrams, respectively).
Figure 6 (top) shows the comparison of the latest measured Gator background to the sum spectrum obtained from the Monte Carlo simulations. On the bottom, the individual contributions of the main simulated background components are shown. The resulting activities for the individual chains and radio-nuclides, using the best fit of the data to the Monte Carlo simulations, are presented in Table 5. The background of the detector above ∼500 keV is thus currently dominated by the residual 238U, 232Th, 40K and 60Co activities in the Cu of the shield, and, to a lower extent, in the copper of the cryostat. Below ∼500 keV, the background is dominated by 210Bi bremsstrahlung from 210Pb decays in the innermost lead shield. The obtained 210Pb activity is with 5.7±0.5 Bq/kg higher than 3 Bq/kg, which is the value specified by the provider.
Specific activity [μBq/kg] Activity [μBq/m3] Activity [Bq/kg]
(chain)
Cu (shield) Ge crystal Cu (cryostat) Sample chamber Inner Pb shield
226Ra (238U)
56 ± 11 8 ± 5
228Th (232Th)
27 ± 7 4 ± 2
32 ± 13 <1.30 11 ± 6
60Co
8 ± 4 <0.80 1.3 ± 0.4
<0.27 <0.11 <0.22
222Rn (238U)
210Pb
Table 5: Activities of residual and cosmogenic radionuclides inside the copper shield, the Ge crystal and the cryostat and of 222Rn inside the shield.
Figure 6: (Top) Comparison of the sum of the simulated spectra from natural and cosmogenic radionuclides in the detector and shield materials (black) with the observed background spectrum (red data points). (Bottom) The individual, best-fit contributions to the observed spectrum are shown: natural radioactivity in Cu (blue), cosmogenic radio-nuclides in Ge and Cu (green), 222Rn decays inside the shield (magenta) and 210Pb decays in the Pb shield (yellow).
5.2 Sample Analysis
Prior to the screening of a specific sample, an estimate of the minimal measuring time tmin, based on the mass, shape and the targeted activity Atarget is performed:
tmin=LdAtarget⋅r⋅m⋅ϵ, (5.0)
where m is the mass of the sample, r is the branching ratio, ϵ is the detection efficiency and the detection limit Ld is defined in equation (4.1). Typical numbers for tmin are 1-3 weeks in order to achieve a sensitivity below a few mBq/kg. Figure 7 shows the sensitivity versus counting time for different sample masses using a specific gamma line (right) and for various radio-isotopes for a specific sample (left).
Figure 7: (Left) Sensitivity versus time for 40K (black), the 238U chain (using the 214Bi 609 keV line, red), the 232Th chain (using the 212Pb 239 keV line, green) and 60Co (blue) for a 13.5 kg sample of PTFE. (Right) Sensitivity versus counting time for different PTFE sample masses: 4 kg (black), 8 kg (red), 12 kg (green) and 16 kg (blue) using the 609 keV gamma line from 214Bi.
The activities of the screened materials are calculated using the analysis methods described in Section 4.1 and 4.2. As evident from Table 6, where a few results from screened samples are provided, an excellent agreement between the results of the two methods was achieved. It provides a cross-check on the reliability of the used methods and on the obtained activities (or upper limits) for the samples that were screened. A large set of measurements is described and discussed in detail in  [21] .
Method Activity 226Ra 228Th 40K 60Co
(amount, time)
χ2 mBq/kg <0.5 <0.5 <0.9 0.24 ±0.08
(18 kg, 20 d)
γ-lines mBq/kg <0.3 <0.3 <1.3 0.24 ±0.06
χ2 mBq/kg 4.1±0.6 1.4±0.4 <4.6 7.4±1.1
(6.6 kg, 6.8 d)
γ-lines mBq/kg 4.3±0.9 <1.8 <5.7 7.2±0.9
χ2 Bq/kg 17±2 3.1±0.7 53±4 <0.6
(35 g, 0.7 d)
γ-lines Bq/kg 15±2 3.8±0.8 42±6 <0.7
χ2 mBq/PMT <0.2 0.20±0.09 8±1 0.6±0.1
(22 pieces, 5.5 d)
γ-lines mBq/PMT <0.2 0.18±0.06 11±2 0.6±0.1
Table 6: Comparison of results obtained by the two different analysis methods using four different screening samples.
The Gator screening facility at LNGS, described in detail in this paper, includes one of the world's highest sensitivity HPGe spectrometers. It allows to measure large samples and to reach a sensitivity on specific activities below ∼mBq/kg when using typical sample masses of a few kg and measuring times extending from one to several weeks. The integral counting rate of the detector in the 100-2700 keV energy region has decreased by more than a factor of 5 from (0.842±0.005) counts/min at Soudan to (0.157±0.001) counts/min at LNGS. This decrease is mainly due to the improved shield and radon protection system, but also because of the decay of cosmogenic isotopes in the crystal and surrounding copper. The background of the facility, which was modeled using a detailed geometry of the spectrometer and its shield, and Monte Carlo simulations with Geant4, is dominated by U/Th/K decays in the copper of the cryostat and the shield, and from residual 210Pb in the innermost lead shield. This results from a comparison of background measurements with the simulations, by choosing the best fit of the MC spectra to the measured data. Two different data analysis methods show an excellent agreement when applied on a variety of samples. Detailed screening results for the XENON100 [6] experiment are presented and discussed in  [21] . At present, the Gator facility is operated under stable conditions at LNGS. It is used to screen materials for the construction of the XENON1T experiment, for the next phase of the GERDA project, as well as for a future noble liquid dark matter search facility, DARWIN  [20] .
We thank James Beaty from the University of Minnesota for help with the detector operation at the Soudan laboratory and the machine shop crew at the RWTH Aachen, in particular Dipl. Ing. Michael Wlochal for the collaboration in the design and construction of the shield. We thank Giuseppina Mosca and Stefano Nisi from the chemistry lab at LNGS for the assistance in cleaning the samples, the LNGS staff, in particular Ing. Piergiorgio Aprili, for their continuous support and the XENON collaboration for the help in maintaining the facility. This work is supported by the Volkswagen Foundation, by the University of Zurich, and by the Swiss National Foundation Grant No. 20-118119 and No. 20-126993.
[2] G. Jungman, M. Kamionkowski, and K. Griest, Phys. Rept. 267, 195 (1996), hep-ph/9506380.
[3] G. Heusser, Ann. Rev. Nucl. Part. Sci. 45, 543 (1995).
[4] J. Angle et al. (XENON10 Collaboration), Phys. Rev. Lett. 100, 021303 (2008), 0706.0039.
[5] E. Aprile and L. Baudis, (XENON100 Collaboration), PoS IDM2008, 018 (2008), 0902.4253.
[6] E. Aprile et al. (XENON100 Collaboration), Phys. Rev. Lett. 105, 131302 (2010), 1005.0380.
[7] E. Aprile et al. (XENON100 Collaboration) Phys. Rev. D 83, 082001 (2011) 1101.3866.
[8] Gerda experiment, http://www.mpi-hd.mpg.de/gerda/home.html.
[9] M. L. H. Neder, G. Heusser, Applied Radiation and Isotopes 53, 191 (2000).
[10] G. Heusser, M. Laubenstein, N. Neder, Proc. of Intern. Conf. Isotop. Environm. (2004).
[11] G. Knoll, Radiation Detection and Measurement (John Wiley and Sons, New York, 2000).
[12] Canberra, http://www.canberra.com/.
[13] Aurubis, http://www.aurubis.com/en/home/.
[14] Plombum, http://www.plombum.republika.pl/.
[15] SOLO facility, http://particleastro.brown.edu/SOLO/.
[16] J. Angle, Ph.D. thesis, University of Florida, Florida (2008).
[17] S. Agostinelli et al. (GEANT4), Nucl. Instrum. Meth. A506, 250 (2003).
[18] M. Heisel, F. Kaether, and H. Simgen, Applied Radiation and Isotopes 67, 741 (2009), ISSN 0969-8043.
[19] C. Hurtgen, S. Jerome, and M. Woods, Applied Radiation and Isotopes 53, 45 (2000), ISSN 0969-8043.
[20] L. Baudis (DARWIN consortium) (2010), PoS IDM2010, 122 (2010), 1012.4764.
[21] E. Aprile et al. (XENON100 Collaboration), Astroparticle Physics 35, 43-49 (2011).