Characterization of inverted coaxial 76 Ge detectors in GERDA for future double- β decay experiments.

Neutrinoless double- β decay of 76 Ge is searched for with germanium detectors where source and detector of the decay are identical. For the success of future experiments it is important to increase the mass of the detectors. We report here on the characterization and testing of five prototype detectors manufactured in inverted coaxial (IC) geometry from material enriched to 88% in 76 Ge. IC detectors combine the large mass of the traditional semi-coaxial Ge detectors with the superior resolution and pulse shape discrimination power of point contact detectors which exhibited so far much lower mass. Their performance has been found to be satisfactory both when operated in vacuum cryostat and bare in liquid argon within the Gerda setup. The measured resolutions at the Q-value for double- β decay of 76 Ge ( Q β β  = 2039 keV) are about 2.1 keV full width at half maximum in vacuum cryostat. After 18 months of operation within the ultra-low background environment of the GERmanium Detector Array (Gerda) experiment and an accumulated exposure of 8.5 kg · year, the background index after analysis cuts is measured to be 4 . 9 - 3.4 + 7.3 × 10 - 4 counts / ( keV · kg · year ) around Q β β . This work confirms the feasibility of IC detectors for the next-generation experiment Legend.

Introduction

As of now the hypothetical neutrinoless double beta () decay eluded detection in any of the candidate isotopes (e.g. Ge  [1, 2], Te  [3], Xe  [4, 5]). A discovery would have far reaching consequences in particle physics, revealing for the first time lepton number violation and hence bringing an additional input to the unresolved matter-antimatter asymmetry creation scenarios [6]. In absence of a signal at the level of  year half-life, next generation experiments are under preparation with significant increase of the isotope mass and further suppression of the background. Latest Ge-based experiments, GERmanium Detector Array (Gerda) and Majorana Demonstrator, benefited from the high background rejection capability of point contact Ge detectors. They also feature an excellent energy resolution at the Ge double- decay Q-value (2039 keV), referred to as hereafter, as compared to other detector types. Among the many challenges inherent to these experiments, one is to scale up the isotope mass by one order of magnitude while lowering the background, which requires a change of germanium detector design.

Previous studies have demonstrated the excellent performance of inverted coaxial (IC) point contact Ge detectors [7, 8]. Their new design allows to increase the detector mass of Broad Energy Ge (BEGe) [9] and P-type, point contact (PPC) [10] detectors from 0.7–1.0 kg up to 2.0–3.0 kg. As a result, the amount of surrounding radioactive components like cables and mechanical supports are lowered per detector mass as well as the surface-to-volume ratio. The latter lowers the rate of surface background events per mass unit.

The background rejection is based on the analysis of the time profile of the induced current signal, called pulse shape discrimination (PSD). Ge double- decays in solid germanium feature a single-site energy deposition within  mm (single-site event, SSE). Background events from Compton scattered rays or surface events result in a different pulse shape [11], as will be discussed below.

Gerda operated the germanium detectors in a 64 m liquid argon (LAr) cryostat [12]. In May 2018, five enriched IC Ge detectors, produced in collaboration with Mirion Technologies [13], have been added to the setup. The total mass of enriched Ge detectors has increased by 9.6 to 44.6 kg. Overall they contribute with 8.5 kgyear exposure, i.e. 8% of the total exposure of 103.7 kgyear. This has been the first time IC detectors have been operated in LAr and important experience has been gained. PSD together with the veto based on argon scintillation light resulted in the lowest background around , thus allowing Gerda to remain background-free over its entire exposure [1, 14].

The IC detector design has been selected by the Legend collaboration, which aims to deploy 200 kg of Ge detectors in the modified Gerda setup by 2021 and ultimately 1 ton in a subsequent stage in order to overpass the year half-life sensitivity [15].

In the following, the characterization of the five detectors in vendor vacuum cryostats, performed in March 2018, is described together with a summary of the performance. Final results within the Gerda LAr cryostat after 18 months of operation are shown in the last section.

Inverted coaxial design and detector production

In 2017, 20 kg of germanium, enriched at 87.7(5)% in Ge, have been purchased in form of GeO powder. The reduction and purification process to convert the GeO to 50 cm electronics grade Ge material took place at the PPM Pure Metals company in Goslar, Germany. About 95% of the material had the required quality and 1.6% was lost during etching and cutting of the bars. The rest of the material was too little for another purification process (zone-refinement) and has been stored for future processing.

The Ge material bars were shipped to Mirion Technologies – Oak Ridge, US in order to further proceed with crystal pulling. Three large crystals were grown and cut into 5 slices with a total mass of 9.826 kg. A loss of 1.5 kg from etching, cutting and grinding has been reported at this stage. Another 4 kg of Ge material with unsatisfactory purity or remaining of cutting and grinding processes will enter a future production cycle after chemical purification. The final detector fabrication stage took place in February 2018 for about two weeks at Mirion Technologies – Olen, Belgium. A total of 9.618 kg were converted into 5 working IC detectors. Overall, the estimated production mass yield was 51%, slightly below the achieved 53.3% during the previous BEGe detectors production [9].

During production, care was taken to minimize the detector exposure to cosmic rays. It meant storing detectors in underground places near manufacturer sites and transporting them across ocean in a shielded container [9]. A total of 36 days-equivalent time above ground over 1 year has been estimated.

In Fig. 1(left), the geometry of an IC detector is sketched. The detector is cylindrically symmetric. The important novelty of the IC design is the well opposite to the p+ contact. The electric field (see Fig. 1(right)) created by the bulk net impurities and the externally applied high voltage (Fig. 2(left)) pushes holes towards the p+ contact and electrons towards the n+ contact. The well is also important to deplete the detector at the usual operational voltage of about 4000 V. Detectors without the well (e.g. BEGes) are limited to less than 1 kg of mass while IC detectors might have a mass above 3 kg. Typical charge collection times of such detectors lie within 1–1.5 s (see Fig. 3) which is 2–3 times larger than that of BEGes.

Fig. 1

Left: Main IC detector features. Middle: ADL calculation of the weighting potential. Right: Electric field strength in kV/cm. The minimum required electric field is 200 V/cm (dark blue). The black dashed lines show electron drift paths ending at the n+ contact while white solid lines are the hole drift paths reaching the signal contact

Fig. 2

Configurations used for detector characterization. Left: Setup for depletion voltage estimation with a mixed 1.5 kBq source of Co–Cs–Am; bias and readout circuits are indicated. Middle: Setup for PSD studies with a flood top and side 13 kBq Th source and with a collimated 250 kBq Th source for lateral scans. The vacuum cryostat (cyan) and detector holder (orange), both made of aluminum, are added here for illustration. Right: Setup for scans with the collimated 4.3 MBq Am source: lateral at 3 azimuthal angles (dashed lines), 2 orthogonal directions on top (solid lines) and a circular one (dotted lines)

Fig. 3

Waveform examples of SSEs (top) and MSEs (bottom) of detector 50A after applying a moving window average. The amplitude of the maximum current A and of energy E are explicitly shown. The boundaries of the rise time are estimated at 0.5 and 90% of the maximum charge amplitude

Given an impurity concentration profile of the bare Ge crystal from the manufacturer after the growth, the geometrical dimensions were optimized via the electrostatic simulation package ADL [16] by Gerda in order to ensure a minimal electric field of 200 V/cm and a maximum depletion voltage of 4000 V. Based on this optimization, the detectors were finally machined by the manufacturer.

The current signal I(t) induced by the drifting charges q as a function of time t can be calculated with the weighting potential according to the Shockley–Ramo theorem [17]:with being the drift velocity of the charges. Figure 1(middle) shows the weighting potential. I(t) is maximal in the vicinity of the p+ contact. Since all holes follow the same path at the end of their drift, close to the p+ contact (funnel effect [18]), the maximum A of the induced current is independent of the starting point. As a consequence A/E, with E being the deposited energy, is constant for all SSEs, while it will be reduced for multiple well-separated depositions (multi-site events, MSEs – see Fig. 3) and reduced or enhanced for surface events. Latter events, originating in Gerda from K decays on the n+ contact and particles on the p+ contact [19], are not studied in this paper due to the inadequate setup. The PSD for such point contact detectors is therefore primarily based on the ratio A/E [11, 20].

In order to understand the response of the new large detectors and evaluate the effect of this well on the detector performance, an extensive characterization with radioactive sources was needed.

Detector characterization in vacuum cryostat

The characterization took place in the HADES (High Activity Disposal Experimental Site) underground laboratory in Mol, Belgium [21], benefiting from a 500 m water equivalent overburden with negligible hadronic background and low Ge activation rate from atmospheric muons. A description of the HEROICA platform that was used can be found in [22]. Similar measurements were performed and the same analysis protocols were applied as in [23]. A summary of the general detector information is given in Table 1.

Table 1

Diameter (D), height (H), well depth (W) and mass (M) of the five inverted coaxial detectors. Their well diameter is 10.5 mm. The dimensions were provided by the manufacturer while masses were measured by Gerda

Det. ID	48A	48B	50A	50B	74A
D±0.2 (mm)	74.6	72.6	74.2	72.6	76.6
H±0.3 (mm)	80.4	80.5	80.4	85.4	82.3
W±0.5 (mm)	47.4	56.0	40.0	53.9	52.4
M±0.5 (g)	1918.9	1815.8	1881.1	1928.7	2072.9

Experimental setup

Each detector was mounted in a vacuum cryostat and cooled down to 95 K by a cold finger immersed in liquid nitrogen. In total, three test benches were used, two static tables with lead and copper shields and one unshielded scanning table (see Fig. 2), with the following goals.A Gerda-like data acquisition system [24] was used for these measurements. The signal was first amplified before being digitized at 100 MHz by a VME-Struck FlashADC (FADC). For each triggered event, a 10 s long trace and a 25 MHz down-sampled trace of 160 s length were stored to disk for offline processing. The former trace was used to study the shape of the current pulse, i.e. A and signal rise time, while the latter was used for the energy reconstruction. The standard Gerda analysis software, Gelatio [25] has been used for this study. For the energy reconstruction, a trapezoidal digital filter has been selected with the following parameters: 4 s rise and fall time and 2 s flat top. Due to occasional FADC failures, part of the data have been recorded by a Multi-Channel Analyser (MCA). This acquisition mode provided an energy spectrum measured at 10s shaping time without pulse shape information.

Determine the nominal bias voltage. The analysis of the 1333 keV line properties of Co was performed for various applied bias voltages.

Probe the geometrical detector response. The homogeneity of the detectors’ surface response was finely scanned by embedding a highly collimated Am source on a 3D movable arm [22]. The source was typically moved by steps of 1 mm every 5 min. The 3 cm thick collimator with an aperture of 1 mm diameter produces a spot for which 95% of the 60 keV rays are absorbed in Ge within 3 mm.

Estimate the best achievable energy resolution. Construct the resolution curve from all available lines full width at half maximum (FWHM) as measured in optimal electronic noise data taking condition.

Evaluate the pulse shape discrimination performance. The Th energy spectrum and pulse shape parameters were recorded for typically a few hours using a 12 kBq non-collimated source to probe the overall detector bulk response. Additionally, due to the large size of the IC detectors, these parameters have been coarsely scanned at four different heights with a 250 kBq Th source mounted behind a 20 cm long collimator of 10 cm diameter with an aperture of 2.5 mm diameter. This special measurement was performed at the Technische Universität München underground laboratory.

Nominal bias voltage

The energy resolution of the detectors has been studied under various applied bias voltages in order to retrieve the nominal bias voltage. It is defined here as the voltage where the energy resolution reaches its minimal and stable value. The bias voltage was varied from 2000 V up to 4700 V in steps of 100 V and data taken with a mixed 1.5 kBq source of Co–Cs–Am (see Fig. 4). The results of the five detectors are illustrated for the 1333 keV Co line in Fig. 5 and the retained nominal bias voltage for each detector is listed in Table 2. The peak position and peak integral variations were also studied. For each detector, the deduced nominal bias voltage agreed within 100 V with the estimate from the energy resolution measurements. All nominal bias voltages were found to be lower than 4000 V, fulfilling the required specification.

Fig. 4

Spectrum taken with detector 50A and the mixed source of Co, Cs, Am for the determination of the nominal bias voltage. The inset shows the fit to the Co 1333 keV line. The MCA module was used for this measurement using a Gaussian energy filter, thus explaining the significant tail of the line from ballistic deficit

Fig. 5

Energy resolution FWHM of the Co 1333 keV line as a function of the applied bias voltage. The arrows show depletion voltages and resolutions reported by the manufacturer. The statistical uncertainties are less than the widths of the markers

Table 2

Nominal bias voltages from scans compared to simulations and manufacturer depletion voltages. Also shown are the operational voltages of the measurements at HADES and Gerda

Det. ID	48A	48B	50A	50B	74A
Nominal volt. (V)	2900	3200	2900	3200	3800
ADL simu. (V)	3010	3335	2140	3320	3620
Manufacturer (V)	3500	3500	3000	3500	4000
Set in HADES (V)	4000	4000	4000	4000	4500
Set in Gerda (V)	3300	3200	3700	3800	4400

Theoretical field calculations were performed with the ADL software [16] for gradually increasing bias voltages. Full depletion was assumed when a minimal electric field of 10 V/cm was reached everywhere in the detector bulk. The deduced depletion voltages agreed with the data at the level of 100 V, except for detector 50A (see below).

The disagreement can mainly be attributed to the uncertainty on the net impurity concentration profile (amplitude and gradient) provided by the manufacturer. As a general rule, the operational voltage in the Gerda experiment was set to 400–600 V above the quoted nominal bias voltage (see Table 2). Detector 50A exhibits a peculiar behaviour: different from the other detectors, its energy resolution remains rather constant between 2000 and 2800 V, and the ADL simulation predicts a nominal bias voltage that is lower by 760 V than measured. The ADL simulation does explain these observations showing a pinch-off (see discussion in [23]) at bias voltages within 2100–2700 V, i.e. a persisting very low electric field region close to the p+ contact. Hence a much higher bias voltage is needed to establish sufficient field strength at the p+ contact.

As to the recommended operational voltage, we verified that the norm of the electric field fulfilled the minimal electric field specification of 200 V/cm (see Sect. 2).

Detector 48A could not be operated beyond the nominal bias voltage in Gerda because the noise from the strongly increased leakage current led to unacceptable data quality.

For characterization measurements in vacuum cryostats, described in the following, the voltages recommended by the manufacturer, which are typically 500–1000 V higher than the depletion voltages, were applied (see Table 2).

Energy resolution

The energy resolution of all available prominent lines from the full data collection (cf. Fig. 2) is shown for each detector in Fig. 6. The FWHM resolution is estimated from a Gaussian fit to the data after subtraction of the background observed in the side bands. Results compatible with the specification from the manufacturer were found. A global fit () to the data is used to interpolate the energy resolution at which ranges from 2.04 to 2.19 keV FWHM in vacuum cryostat.

Fig. 6

Energy resolution FWHM as a function of ray energy; data are taken at the bias voltage recommended by the manufacturer. The dashed lines show fits to the data, performed for each detector separately. The statistical uncertainties are less than the widths of the markers

Surface response homogeneity

In this section, detector response parameters of interest for the data analysis in Gerda are detailed as a function of the interaction position near the detector surface. We investigated the charge collection efficiency, the [0.5–90%] signal rise time and the pulse shape discrimination parameter by means of scans with a highly collimated low energy ray source (60 keV line from a 4.3 MBq Am source). In total, about 1500 measurements were performed with an exposure of 5 min each. For the sake of conciseness, we show in Fig. 8 and comment below a selection of relevant observables for detector 50A and events originating from the prominent 60 keV line. The number of counts and the peak position are obtained after Compton background and tail subtraction (see Fig. 7). Events found in the tail are typically from rays with partial charge collection due to energy loss at the detector surface [23, 26].

Fig. 8

Results of the lateral (left) and top radial (right) scans for the 60 keV Am line obtained with detector 50A. Statistical uncertainties are less than the widths of the markers. The dashed gray lines on the left(right) show the expected holder ditch (well) position (see Fig. 2)

Fig. 7

Example of Am FEP data recorded with the highly collimated 4.3 MBq source, positioned above the upper surface of detector 50A. The model (red line) is shown together with its decomposition into a Gaussian, a tail and a shoulder functions

Firstly, the lateral scans (Fig. 8a) feature a specific count rate profile mainly affected by the non homogeneous width of the detector holder (cf. two 9 mm ditches in Fig. 2, middle). The profile is well reproduced at different angles. This specific detector has a slight taper on the upper part at heights above 40 mm, leading to a thinner full charge collection depth (FCCD) in this region hence explaining the visible higher count rate starting at  mm. The charge collection efficiency of the 60 keV line, quantified via the peak position (Fig. 8b), decreases as the energy deposition occurs further away from the p+ contact. This behavior, observed at all angles and for all detectors, is understood as charge trapping along the drift path. The effect is lower ( %) than the energy resolution at (0.13%) but sizable. Thanks to its correlation with the charge collection time, it can effectively be corrected for. The reported rise time [0.5–90%] is the average value of all the events in the peak. Its mean distribution (Fig. 8c) increases linearly in the lower part of the detector before reaching a plateau. This plateau is best explained by the charge drift paths (see Fig. 1 right), all converging at a height of 5 cm toward the same low electric field region. Finally, the normalized pulse shape parameter is also scrutinized. is the position of the Gaussian peak of the A/E distribution of 60 keV events. As seen in Fig. 8d) it features two plateaus that are understood via the simulation of charge cloud dynamics. This dynamics comprises diffusion, acceleration and repulsion of the charges that, when summed, increase the spatial distribution of the charge clouds as they start drifting towards the p+ contact. As a consequence of a larger spatial distribution, a lower current amplitude A for the upper scan positions is observed. Electrostatic simulations, performed with the Siggen software [27], reproduce this feature both qualitatively and quantitatively, as detailed in [28].

Secondly, the top radial scans shown in Fig. 8, right, also have an interesting pattern due to the presence of the well. While the average count rate is in agreement with the upper lateral measurement, an increase at the well location is explained by both an increased acceptance and a thinner FCCD originating from the fabrication process. Same observations apply to the charge collection. The mean rise time profile is stable at 1240 ns but peaks off to 820 ns in the well, close to the p+ contact. The shift observed at the well level confirms the interpretation given for the lateral scan since no difference in between well and top events is expected from pure electrostatic simulation, without charge cloud dynamics.1 The small angular shifts apart from the well can originate from a slight misalignment of the well or the groove ( mm) relative to the central detector axis that would create an asymmetric weighting field.

The outcome of these exhaustive scans is summarized for each detector in Table 3 in form of the observed maximal variation of each parameter between the various source positions. From this table, one can appreciate that the shift of the 60 keV line for the entire data collection lies within the same order of magnitude of the energy resolution at , 0.1 and 0.15% in vacuum cryostat and in LAr (see Sect. 4), respectively. As a result, despite the large drift time, a sufficiently good homogeneity of the charge collection over the whole detector volume is ensured. Concerning , shifts of its average value vary from 0.2% up to 1.0%. This is of the same order of magnitude as the typical A/E resolution found at higher energies, close to in vacuum cryostat (see Fig. 10). However, the shift is much lower than what was achievable in Gerda2 [11]. The minimal and maximal rise times are taken from the lateral measurement, i.e. the most relevant values for benchmarking Th source data (see Sect. 3.5). For detector 50B, only MCA data were taken, hence no pulse shape parameters are available and the peak position is less reliable.

Table 3

Maximal variation of peak position (PP) and at the 60 keV line. The minimum and maximum rise time (RT) of the 60 keV line is reported for comparison with Th data

Det. ID	48A	48B	50A	50B	74A
ΔPP (%)	0.06	0.20	0.08	0.33	0.13
ΔA/E (%)	0.21	1.05	0.43	–	0.75
min RT (ns)	970	925	980	–	920
max RT (ns)	1190	1200	1240	–	1150

Fig. 10

Correlation of A/E with the signal rise time [0.5–90%] of detector 50A for indicated heights H of the lateral collimator. The p+ contact is at  mm. The gray dashed lines show the A/E cut position. These four datasets correspond to Tl DEP events from the highly collimated Th source

Pulse shape discrimination performance

Pulse-shape analysis provides a simple but powerful technique to discriminate signal-like SSEs from background-like MSEs which is crucial for realizing the background-free condition [14]. In the following analysis, the Tl double escape peak (DEP) has been used as a proxy for SSEs while MSEs have been considered from various Compton scattering regions, including the Compton continuum at  keV, referred to as CC @ . During the characterization campaign, three types of data were recorded with a Th source: two flood measurements from top and lateral positions and collimated lateral measurements at four different heights (see Figs. 2 and 9) with the aim to resolve the spatial dependency of the A/E parameter.

Fig. 9

Th spectrum taken with detector 50A and the source in lateral position. The main lines, Tl double escape (DEP), Tl single escape (SEP) and Tl and Bi full energy (FEP) peaks, used in the pulse analysis, are emphasized together with the  keV Compton continuum region. The inset shows a fit to the 2615 keV line

Correlation between A/E and rise time

From Sect. 3.4, the A/E ratio of SSEs is expected to be distributed around two values from energy depositions that happen close to and far away from the p+ contact. These two values are close but could be resolved in the low electronic noise environment of the vacuum cryostat. On Fig. 10, the lateral scan data show the two SSE A/E populations from the Tl DEP line of detector 50A: short rise time with high A/E and long rise time with low A/E. At the closest position to the p+ contact ( mm), the dominant SSE population has a rise time between 800 ns and 950 ns. This is about 50 ns shorter than what is observed in the lateral Am source scan (Fig. 8c).

This observation can be explained in part by the fact that the 2615 keV rays penetrate deeper in the detector and hence probe a region on average closer to the p+ contact. Additionally, Am source lateral scans are performed along a single detector orientation as opposed to flood measurements that probes a large fraction of the detector bulk. As a result, scans produce pulses with a specific charge velocity, sensitive to the closest crystal axis, that is not representative of the average bulk velocity. At higher positions up to the top ( mm), the main population is moving towards longer rise times of about 1200 ns. The shift between the two A/E populations is found to be 0.4% for this detector which is compatible with the observed variation in Am source measurements.

Background rejection

The rejection of MSEs for the two rise time populations (Fig. 11) was studied by considering the lateral Th source flood measurement (Fig. 2). The effect of charge cloud diffusion is here most enhanced due to well distributed event statistics across the detector height. The PSD cut is set on the low A/E distribution side to accept 90% of the Tl DEP events. Various peaks of interest featuring variable proportions of SSEs and MSEs have been investigated. Their survival fractions after applying the cut, for four different analysis configurations, are summarized in Table 4 for detector 50A. The lateral analysis dataset combines all data without any correction or selection. In addition, the corrected analysis shows the PSD performance when the linear rise time dependence of A/E (cf. Fig. 10) is corrected for. This energy dependent linear correction aims to cancel the position dependence of the SSE A/E by empirically minimizing the width of the SSE band A/E distribution. This is done by (1) fitting its rise time dependence at many Compton continuum energy regions above 1 MeV:and, (2) applying this correction to all events:

Fig. 11

Left (right): Rise time (A/E) distribution of Tl DEP events obtained from the lateral Th source flood measurement with detector 50A. The lateral, short and long analysis datasets selection are emphasized

Table 4

Survival fractions (in %) and their statistical uncertainty obtained with detector 50A and the low-sided PSD cut at indicated line regions and the Compton continuum around (CC@)

Event class	Lateral	Corrected	Short RT	Long RT
^208Tl DEP	90(2)	90(2)	90(3)	90(3)
^212Bi FEP	9.3(4)	8.3(4)	10.2(7)	5.9(5)
CC @ Q_ββ	33.4(6)	33.8(6)	37.1(9)	30.6(8)
^208Tl SEP	6.3(4)	6.2(4)	8.3(6)	4.3(5)
^208Tl FEP	8.05(8)	8.47(8)	11.0(1)	6.8(1)

The short RT and long RT analysis datasets comprise only the events that fall into the short and long rise time population, respectively (cf. Fig. 11). Despite a double A/E peak for SSEs, the lateral analysis shows excellent rejection of MSEs (32% acceptance of Compton continuum at ) as compared to the Gerda BEGe detectors (40%) [23]. Similar performance is observed if the linear rise time correction is applied to A/E. However, this position dependence correction comes at the expense of a better understanding of its energy dependence and will require additional studies on other classes of events found in Gerda [19]. The short RT and long RT analysis quantify the ability of our PSD method to discriminate SSEs and MSEs in the lower and upper parts of the detector, respectively. As expected, it does perform significantly better for the upper part as the separation in time between the hits is on average larger. Care must be taken not to overinterpret this qualitative observation as MSEs do feature, by definition, longer rise time.

The background rejection performance of all five IC detectors for the top Th source measurements, without correction, is summarized in Table 5. The values are comparable to the ones obtained with BEGe and PPC detectors [20, 23].

Table 5

Survival fractions (in %) and their statistical uncertainty obtained with top Th source position and the low-sided PSD cut at indicated energy regions

Event class	48A	48B	50A	50B	74A
^208Tl DEP	90(2)	90(1)	90(1)	90(1)	90(3)
^212Bi FEP	6.9(9)	8.0(5)	8.6(4)	7.9(4)	12.1(9)
CC @ Q_ββ	37.1(8)	35.8(4)	35.7(6)	36.6(5)	34.5(8)
^208Tl SEP	6.4(7)	6.6(3)	6.5(3)	5.6(3)	7.6(6)
^208Tl FEP	9.8(1)	8.70(7)	9.28(8)	8.51(7)	–

Long-term performance in the GERDA cryostat

From July 2018 to November 2019, the five IC detectors have been operated within the Gerda LAr cryostat. Among these detectors, detector 48B could not be fully depleted due to high leakage current over the entire data taking and has been biased with [3000-3200] V only. As a consequence it was used in anti-coincidence mode only and excluded from the high level decay analysis. The averaged FWHM energy resolution at of the four others, observed during 18 months of operation, was 2.9(1) keV. During this data taking period, the background rate of particles was studied to validate our detector handling procedure [19], from detector production to integration in LAr. We observed a constant rate of 4.5 particles per month. This is compatible with the one from BEGe detectors before the upgrade if normalized by the total p+ contact surface.3 It demonstrates the successful handling of the detectors during their fabrication and integration. The weekly Th source calibrations were used to study the PSD performance of these detectors in a similar manner as described in [11]. The much higher electronics noise did not allow to resolve the two SSE populations shown in Fig. 11. While accepting 90% of the Tl DEP events with the low-side A/E cut, the survival fractions of Tl single escape peak (SEP) events was 7.4(7)% and 35(2)% for CC @  on average. This is, respectively, 2 and 7% better than the achieved rejection performance with BEGe detectors [29]. None of the high energy particles survived the high A/E cut, set to +3 standard deviations of the A/E SSE band distribution.

Figure 12 shows the combined physics spectrum of the four working IC detectors at the cumulated exposure of 8.5 kgyear, before and after LAr veto and PSD cuts. For a detailed description of the various features (2 decay continuum, lines from K and K and high energy particles) we refer to [1, 19]. An additional line was found on top of the usual energy spectrum at 1125 keV before analysis cuts. This peak originates from Zn which produces a 1115.5 keV EC ray accompanied by a 10 keV deexcitation X-ray. This cascade is however highly suppressed by the PSD cut. Such a detector bulk contamination has already been reported for some of the enriched Gerda coaxial detectors in [30] at the beginning of the Heidelberg-Moscow experiment. With a half-life of 244 days, this line had not been observed in previous Gerda datasets since coaxial and BEGe detectors had been stored at least 2 years underground before being deployed.

Fig. 12

Spectrum measured with the IC detectors at the exposure of 8.5 kgyear in Gerda prior to and after indicated analysis cuts. The inset shows a zoom in the background analysis window. The only surviving event at 2058.9 keV was recorded with detector 74A on October 9, 2018, 01:09:14 (UTC)

The final Gerda analysis led to a single event surviving the cuts in the [1930-2190] keV region of interest (see inset of Fig. 12) which, with the omission of the two lines at 2104 ± 5 keV from Tl SEP and 2119 ± 5 keV from Bi, was used to deduce the background index [1]. For the total exposure of 8.5 kg.year the value is .

Conclusion

Five enriched inverted coaxial Ge detectors have been characterized in terms of surface response, energy resolution and background rejection efficiency. High performance level was found for these last two critical items despite a peculiar distribution of the PSD parameter A/E across the detector bulk. After the Gerda upgrade in spring 2018, an exposure of 8.5 kgyear has been validated until fall 2019. The standard analysis of Gerda reported a background rejection compatible with that of BEGe detectors for the background and slightly improved for multi-site events at , while the active mass of the IC detectors is up to a factor 3 larger. This is the first successful and stable long-term operation of IC detectors in LAr. These very promising results serve as a proxy for estimating the PSD performance in ultra-low background environment for the upcoming Legend-200 experiment, currently under preparation at the INFN Laboratori Nazionali del Gran Sasso.