α -event characterization and rejection in point-contact HPGe detectors.

P-type point contact (PPC) HPGe detectors are a leading technology for rare event searches due to their excellent energy resolution, low thresholds, and multi-site event rejection capabilities. We have characterized a PPC detector's response to α particles incident on the sensitive passivated and p + surfaces, a previously poorly-understood source of background. The detector studied is identical to those in the Majorana Demonstrator experiment, a search for neutrinoless double-beta decay ( 0 ν β β ) in 76 Ge. α decays on most of the passivated surface exhibit significant energy loss due to charge trapping, with waveforms exhibiting a delayed charge recovery (DCR) signature caused by the slow collection of a fraction of the trapped charge. The DCR is found to be complementary to existing methods of α identification, reliably identifying α background events on the passivated surface of the detector. We demonstrate effective rejection of all surface α events (to within statistical uncertainty) with a loss of only 0.2% of bulk events by combining the DCR discriminator with previously-used methods. The DCR discriminator has been used to reduce the background rate in the 0 ν β β region of interest window by an order of magnitude in the Majorana Demonstrator  and will be used in the upcoming LEGEND-200 experiment.

Introduction

backgrounds in neutrinoless double beta decay searches

The discovery of neutrinoless double-beta decay () would indicate that neutrinos are Majorana particles and that lepton number is not conserved. It could also provide information on the absolute mass scale of the neutrino, the source of the neutrino’s mass, and the origin of the matter/anti-matter asymmetry of the universe [1-3]. This radioactive decay would occur only rarely; sensitivity limits from current experiments indicate that the half-life is over  years [4, 5]. Detecting it therefore requires large-mass experiments with extremely low background rates and the best possible energy resolution.

P-type point contact (PPC) High-Purity Germanium (HPGe) detectors [6] are a key technology for rare event searches, capable of detection thresholds below 1 keV and full-width at half-maximum (FWHM) energy resolutions of 0.12% (2.5 keV) at the region-of-interest of 2039 keV [7]. They have been used to search for low-energy nuclear recoils from external sources in dark matter searches, [8] and have been proposed for use in coherent neutrino-nuclear scattering experiments [9]. They are also a leading technique for searches, where the detector itself is the source [7, 10]. The experiments using PPC detectors have the lowest background rates of any of the currently-operating experiments [10].

A troubling source of backgrounds for experiments is the decay of radon isotopes and their progeny, particularly Rn, on or near the surface of the detectors. Long-lived -emitters such as Po are the most concerning; the 138 day half-life of this isotope makes the use of timing correlation-based rejection prohibitively inefficient. If its decay is supported by the decay of Pb, the associated 22 year half-life will make the Po rate approximately constant throughout the life of the experiment. Since the particle is emitted with 5.304 MeV and that energy is easily degraded when passing through even a thin layer of inactive material, it can appear in the region of interest (ROI) around the Ge Q-value of 2.039 MeV.

Because the range of a 5.304 MeV particle in germanium is less than 20 m [11], all decays from sources outside the detector are considered surface events in HPGe detectors. There is no evidence for bulk events, which would have to originate from contaminants in the HPGe material. If these events did occur, they would not contribute to the ROI, since their full energy would be deposited inside the detector.

In PPC detectors, the majority of the surface is covered in a 1–2 mm thick lithium-diffused dead layer forming the n contact [12], and is therefore completely insensitive to decays. The remaining sensitive surfaces are those of the central p contact region, formed by boron implantation in the germanium bulk, and the passivated surface, a m layer of amorphous germanium or silicon. The dimensions of these regions depend on the chosen detector geometry; two common designs are the Mirion1 BEGe PPC and ORTEC PPC detectors. The ORTEC PPC detectors used in the Majorana Demonstrator have a passivated surface that covers an entire circular face of the detector, with area of approximately 30 cm per detector.

The response of a PPC detector to events near the passivated surface is difficult to predict, and depends on the detector manufacturing process. Charge trapping has been observed on similar surfaces in other HPGe detector geometries [13], but the charge collection properties can differ depending on the surface treatment and field configuration. In the Majorana Demonstrator detectors, events have been observed in which particles originating on this surface are significantly degraded in energy, appearing in the ROI, as discussed in Ref. [7]. We performed a dedicated study of interactions on this surface, leading to more reliable models of the energy spectrum and the distinctive pulse-shape characteristics of these signals.

The delayed charge recovery effect

Based on the characteristics of interactions, it appears that charge mobility is drastically reduced on or near the passivated surface. Therefore, a fraction of the charge from these interactions is slowly released on the timescale of waveform digitization, leading to a measurable increase in the slope of the tail of the waveform.

An offline digital filter can be used to identify these events, allowing for the efficient rejection of passivated surface events [7]. The goal of such a filter is to detect the presence of slow charge collection occurring after the bulk charge collection has been completed. This filter is used to calculate the Delayed Charge Recovery (DCR) parameter. In a waveform that has been corrected for the electronic response function,this appears as a positive slope of the tail, as seen in Fig. 1.

Fig. 1

Baseline-removed and pole-zero corrected waveforms from PONaMa-1 events, taken with the TUBE scanning system. Both the bulk event (in blue) and the surface event, taken with the source incident at r mm (in red), have 2615 keV of energy. The waveform regions used to calculate the DCR parameter for each waveform, determined as described in Sect. 2.4, are indicated by the shaded boxes, shown in the same color as the waveform they correspond to

The delayed component of the surface signal can be modeled by considering the motion of electrons and electron-holes near the passivated surface. Computational models of charge drift in germanium detectors have shown that on the surface of the detector, the carrier drift velocities may be 10–100 times lower than in the detector bulk [14]. In the “slow surface drift” model, some fraction of the charge carriers are driven to the passivated surface, with the remainder of the carriers being collected normally through the bulk. This behavior could be caused by the presence of a net charge on the passivated surface, or by the self-repulsion of charge carriers in the dense charge cloud created by an interaction; these cases can be distinguished, to some extent, by the differing detector response to such events as a function of interaction position.

Alternatively or in addition to the slow surface drift behavior, some fraction of the charge carriers may be trapped in a few-m-thick region just below the passivated surface and then slowly re-released, dominating the observed slow charge collection behavior. If the detector passivation process causes any damage to the germanium crystal structure (for example, via sputtering causing crystal dislocations), this would create a narrow higher-trapping region near the surface. This would not affect bulk charge collection, but would cause incomplete charge collection in events occurring on or near the surface. In this case, which we term the “near-surface trapping” model, the remainder of the charge carriers can be collected promptly, as they are in bulk events, collected slowly, as in the slow surface drift model, or trapped at the detector surface due to passivated surface charge build-up, contributing negligibly to the signal.

In all of these cases, part of the energy of the event appears as a normal, fast pulse, and the remainder of the charge is collected slowly. Depending on which charge carriers are affected, these models produce different predictions of event energy and DCR as the position of the interaction on the passivated surface changes. A dedicated scan along a radial path of a PPC detector’s passivated surface with a collimated -emitting source can not distinguish between the two causes of delayed charge, but can reveal which charge carriers are being affected. The near-surface trapping and slow surface drift effects may both be present, with the radial behavior governed by the dominant effect.

The expected efficiency of the DCR pulse-shape discriminator can be determined from single-site Th events, as discussed in Ref. [7]. Its leakage (i.e. the fraction of events that are misidentified as signal-like events), however, requires a priori knowledge of the number of events in the detector. Since the DCR varies as a function of interaction position, the leakage can also vary as a function of position, requiring additional knowledge of the spatial event distribution. A dedicated surface scan using a spectroscopic-grade source is one of the few reliable ways to measure the DCR discriminator’s leakage.

Experimental setup and calibration

The TUBE scanner

The TUM (Technische Universität, München) Upside-down BEGe (TUBE) scanner is a custom-built cryostat first made to study the backgrounds in GERDA due to surface interactions on the p electrode and insulating surface of Mirion BEGe PPC detectors [15]. It allows a PPC detector’s passivated surface to be scanned with a collimated source. Key aspects of the scanner design, following minor modifications needed for use with an ORTEC PPC detector, can be seen in Fig. 2.

Fig. 2

A simplified cross-sectional view of the TUBE scanner, showing key dimensions in millimeters. Details of the detector cup, front-end electronics, and cold-finger are removed for clarity. Both the detector holder (unlabeled) and the IR cup are held at the ground potential, with an insulating spacer placed between the detector holder and the n surface of the detector, which is at high voltage

Measurements of the detector’s response were taken with an open 40 kBq Am spectrometry source with a vendor-determined energy peak full-width at half-maximum (FWHM) of 20 keV. The source was mounted in a 1 mm-diameter collimator. The entire 68.9 mm diameter of the detector’s passivated surface was scanned, save for a 6-mm “blind spot” on the detector surface that is occluded by signal electronics components. The p contact is only partially occluded by the signal electronics. The 0 and 180 positions along the scanned diameter are distinguished by assigning them positive and negative radial positions. The source beam had a 1.8 mm-diameter spot size on the detector surface.

When calculated from the background-subtracted data as described in Sect. 3.3, the average event rate (including only source positions where the beam is not occluded) was  mHz. This is consistent with the 17.6 mHz expected rate derived from the collimator geometry and manufacturer-cited source strength. A muon veto system was used to reduce backgrounds from cosmogenic muons at high energy. Events from both the detector and the muon veto system were recorded with a Struck SIS3302 digitizer, sampling at 100 MHz with a trace length of s. Offline analysis was used to reject multi-site events, reducing the rate of background events, as described in Sect. 2.3.

The scanned detector, named PONaMa-1 (PPC from ORTEC made from Natural Material), was produced by ORTEC. Its fabrication was identical and geometry similar to the Ge-enriched detectors used in the Majorana Demonstrator. The hemispherical p contact was made by boron implantation and is m thick. The passivated surface covers nearly an entire circular face of the crystal, and has a radius of 30 mm. The dimensions and other key characteristics of the detector are given in Table 1, and a diagram of the detector is shown in Fig. 3.

Table 1

Dimensions and operating parameters of the PONaMa-1 PPC detector. The dimensions were determined by the detector manufacturer. The n dead layer thickness, capacitance, depletion voltage, leakage current, and resolution were determined by the detector manufacturer and then confirmed with independent measurements conducted as part of the Majorana Demonstrator detector characterization campaign [16]

PONaMa-1 Properties
Diameter	68.9 mm
Height	52.0 mm
n^+ Dead Layer Thickness	1.2 mm
Passivated Surface Diameter	60 mm
p^+ Contact Diameter	3.2 mm
p^+ Contact Depth	2.0 mm
Capacitance	1.8 pF
Depletion Voltage	850 V
Leakage Current	10 pA
Resolution at 1332 keV	2.05 keV

Fig. 3

Diagrams of the commonly-used PPC detector geometries from ORTEC (top) and Mirion (bottom)

See Appendix A for more details about the source, collimator, and scanner.

Scanning and calibration measurements

The detector was operated with a bias voltage of 1050 V, which is 200 V above the depletion voltage of the crystal. Measurements were taken every 1.5 mm along a radius with intermediate 0.75 mm points near the p contact and passivated surface edge. Each measurement lasted 24 h. Several multi-day runs were taken to study the stability of the system, and a subset of scanning positions were repeated non-contiguously to study the long-term stability of the DCR parameters.

The data analyzed herein correspond to two deployments of the detector. The first deployment was a set of scans taken over 217 days of continuous operation. During this time, the detector was kept biased, cooled, and under vacuum. After these 217 days, the detector was warmed and the cryostat was opened. The source position was adjusted to give a higher event rate, and the detector was put back into operation. The data from this second deployment are used only for studies of the detector response stability, and are not included in other analyses.

The high ambient background rate allowed the energy to be calibrated independently for each data set, without the need for dedicated calibration runs. Additional runs with Th and Th sources were also used to confirm the energy estimation performance and calibrate pulse shape discrimination parameters. For a description of the calibration procedure and discussion of the detector energy scale stability, see Appendix B.

Rising-edge-based pulse shape discrimination parameters

In PPC detectors, the charge drift time varies dramatically as a function of the interaction position. This prompt charge drift forms the shape of the rising edge of the waveform. Pulse shape discriminators that are sensitive to this portion of the signal can provide information including the interaction location and whether interactions occurred at multiple sites within the crystal. The rising edge is also sensitive to changes in the electric field within the bulk of the crystal. In the present analysis, one algorithm, A vs. E, is used for multi-site background event rejection; another, A/E, is used to identify near-point-contact events. These algorithms measure the peak value of the current pulse (A), normalized with respect to the event energy (E). They differ in how the normalization is applied, but are otherwise similar. Their implementation is described in Refs. [17, 18]. A third algorithm, which directly measures the drift time of the pulse, is used to track the stability of the detector’s electric field over time.

events are intrinsically single-site, since the particles have short range in germanium. Ambient background events, which are primarily high-energy Compton interactions, are often multi-site. These backgrounds can be reduced via the use of a multi-site event cut. The most common techniques rely on the peak amplitude of the derivative of the waveform (A, which corresponds to the maximum signal current) as a function of energy (E), which is reduced in multi-site events relative to its value in single-site events. Applying a lower bound to this rising-edge-based discrimination parameter preferentially removes multi-site events. In this work, we use the A vs. E algorithm described in Ref. [18].

A similar parameter, A/E [17], can be used to reject (or select) near-point-contact events. The fast drift times of these events increase their peak current values for a given energy. Therefore, an upper-bound cut on A/E preferentially cuts events that occur near the contact. See Refs. [15, 19] for more discussion of this approach.

Figure 4 shows the A vs. E and A/E distributions with respect to energy in a Th calibration data set. A vs. E, which has higher signal efficiency at low energy, is used to reject multi-site events, as in the analysis employed in the Majorana Demonstrator [18]. A/E is used to identify near-point-contact events, as in the GERDA experiment, where it is used to reject and surface events [19]. A/E can also be used to identify multi-site events, but it is not used for that purpose in this analysis.

Fig. 4

The distribution of A vs. E (left) and A/E (right) values with respect to energy is shown for calibration events after the muon veto is applied. The color scale indicates the number of events. A vs. E is used to identify multi-site background events; all events with A vs. E less than -1 (indicated by the red line) are rejected. A/E is used to identify near-p contact events; 99.9% of bulk events with energies between 1 and 2.63 MeV have A/E less than 2 (indicated by the dashed magenta line), and 99.8% have A/E less than 1.5 (indicated by the dashed black line)

A vs. E and A/E were calibrated based on two Th calibrations, one taken at the start of data-taking and the second taken following an observed gain shift. Both are calibrated such that the lower-bound threshold cut on the parameter accepts 90% of events in the Tl double-escape peak. See Refs. [18, 19] for more details.

The upper-bound bulk acceptance values of A/E are calibrated using runs with no or calibration sources present. They are based on the acceptance of events with energies between 1 and 2.63 MeV after the application of a muon veto cut and a basic pile-up cut. No multi-site event rejection cut is applied. The A/E distribution is normalized so that the 99% acceptance value occurs at a value of 1. With this normalization, the 99.9% acceptance value of A/E is found to be . A cut of A/E , which is used to select near-point-contact events in the energy analysis (see Sect. 3.1), is found to accept 99.8% of events. As seen in Fig. 4, employing a similar near-point-contact event cut in A vs. E would require an energy-dependent upper threshold. Therefore, A vs. E is not used to select near-point-contact events in this analysis. The energy-dependent width of the A vs. E distribution is due to charge cloud diffusion; this effect is corrected, to first order, by the energy normalization used for the A/E parameter. Additional discussion of this effect can be found in Ref. [20], and will be included in an upcoming publication from the Majorana Collaboration.

To study the bulk event pulse shape stability in the detector, as described in Appendix E, we also compute the “ to ” drift time of each pulse. The start time of the pulse, , is found by applying a trapezoidal filter with an integration time of 1.5 s and a collection time of 1 s, identifying the time at which the filtered pulse passes a fixed level threshold of 5 ADC counts, and then correcting for the trapezoidal filter timing offset. The relatively high threshold was chosen to give stable results under changing noise conditions. The 50% rise time of the pulse, , is found by linearly interpolating between sampled points and identifying the first time at which the waveform crosses 50% of its maximum value. The difference between the start time and 50% rise of the pulse is taken to be the drift time. The parameters were chosen to provide sensitivity to the sharp initial rise of near-point-contact events without being subject to noise, and the point is used to avoid the degradation of the drift time parameter by the DCR effect itself.

Tail-based pulse shape discrimination parameters

In a PPC detector operated above its depletion voltage, bulk events exhibit a sharp, step-like transition from the rising edge of the waveform to its tail, as seen in the event shown in Fig. 1. All of the charge collection occurs promptly and ends at that transition. Passivated surface events, on the other hand, exhibit additional slow charge collection after the prompt signal has ended. This DCR effect can be measured via pulse shape discrimination parameters that are sensitive to the shape of the waveform tail.

To calculate the DCR parameter, the resting baseline, taken to be the average ADC value of the first s (500 samples) of each waveform, is subtracted from the waveform. Then the shaping effect of the pre-amplifier, which adds a s pole-zero decay constant to the signal, is deconvolved from each waveform. Second-order shaping effects remain, but their effects are minimized by selecting waveform regions for the DCR calculation that are slightly delayed relative to the end of the waveform rise.

The tail slope () of each waveform is found from a two-point slope calculation. This approach minimizes the parameter computation time per waveform. It can be expressed as:where is the average value of the waveform (in ADC) in the s time window starting at . The first point, , is 2 s after the 97% rise point of the waveform, and is set so that the second average is over the final 1 s of the waveform. In this case, 30 s traces are recorded, so is set to be 29 s after the start of the waveform. The delayed charge recovery continues beyond the end of the digitization window. See Fig. 1 (right inset) for an example.

This unnormalized tail slope gives the rate of delayed charge collection. , the energy that is collected as delayed charge during the digitization window, is given by:where a is the linear energy calibration constant and is the length of the tail of the waveform (s, in this case). Since the digitization window is limited to 30 s, this measured delayed charge does not include all of the delayed charge from each event. Given the limitations of the digitizer used and the high ambient background rate, the true total charge recovered cannot be measured in TUBE.

The calculation of the normalized DCR value is based on calibration data sets, which are contiguous sets of Th or Th calibration runs. Following the application of a muon veto and basic pile-up cut, the distribution of values of for single-site calibration events with energies between 1 and 2.38 MeV is fit with a Gaussian peak. The fit range is set to exclude the high- and low-DCR tails. The 99% and 99.9% acceptance values of DCR are calculated from this same event population, including the events in the tails of the distribution.

The normalized DCR value is then given by:where is the centroid of the fit, is the tail slope, calculated as in Eq. (1), and is the 99% acceptance value of the shifted distribution.

After parameter calibration, the mean value of DCR is 0 and 99% of bulk events have DCR . The 99.9% acceptance value is . The DCR parameter was recalibrated during each Th or Th calibration run, with calibration occurring every 2 to 4 weeks. The uncertainty on the 99.9% acceptance cut value captures the fluctuations in the high DCR tail of the distribution over the course of these calibration runs. The normalized DCR distribution for calibration run events is shown in Fig. 5.

Fig. 5

The DCR distribution in a Th calibration data set. Non-muon single-site events with energies between 1 and 2.63 MeV are used to calibrate the DCR parameter. The centroid of a Gaussian fit (shown in red) is used to shift the tail slope , and the distribution is normalized to the 99% acceptance value (in blue). The 99.9% acceptance value (shown as a dashed violet line) is also provided for reference

We note that the delayed charge component of the waveform tail observed is not truly linear; it is more accurately fit by an exponential rise. Fitting an exponential function to each waveform, however, adds computational complexity to the waveform processing without measurably improving event rejection. Therefore the faster first-order DCR algorithm given by Eqs. (1) and (3) is studied here. The full exponential fit can be used to find the time constant of delayed charge collection, as discussed in Sect. 3.2.

Event response

Event energy

By studying the measured energy of fixed-energy events at a variety of positions on the detector surface, we can probe the mechanism of energy degradation as discussed in Sect. 4. These measurements can also be used to create more accurate spectral models of -emitting background sources in low-background experiments such as the Majorana Demonstrator.

Am events distributed in a broad peak are observed for every source position incident on the passivated surface or p contact. When the source beam is partially or entirely incident on the p contact, an peak in the spectrum appears at nearly the full energy of the emitted . For positions incident on the passivated surface, the events are degraded in energy, with the peak energy and width varying with radius.

Across most of the detector surface, muon veto, pile up rejection, and multi-site rejection cuts are applied, and the remaining peak in the energy spectrum is fit with a Gaussian function. For positions with radii smaller than 6 mm, an additional cut requiring A/E is used to identify near-point-contact events, and the peak fitting function contains an additional component accounting for the peak’s low energy tail. The same approach is used for all small-radius ( mm) data sets. The low-energy tail contains a majority of the events in data sets with  mm, so for these data sets, an additional estimated energy range of the observed events is given. This energy range is given in lieu of the mean position of the Gaussian for the smallest-radius data sets ( mm), where the Gaussian+tail model fit fails. See Appendix C for details of the cuts applied in each case and the peak fitting procedure.

All of the peak energies of the fits to the energy spectra are depicted in Fig. 6. For positions with the source incident on the passivated surface, the full width at half maximum (FWHM) of the peaks ranges from 50 keV, at large-radius positions, to 240 keV, at the smallest-radius positions. At the p contact, the width of the peak is narrower, with a FWHM of 21 keV.

Fig. 6

The centroids of the energy peaks in each data set; certain positions were studied in multiple data sets. Negative (positive) radius positions are indicated by blue downward-pointing (red upward-pointing) triangles. For scanning positions with significant low-energy tailing, the black box depicts the estimated full energy range of events. At positions that are partially or completely incident on the point-contact, an additional peak appears at nearly the full incident energy. The vertical error bars (not visible) depict the uncertainty in the peak position from the maximum-likelihood fit of the peaks, and the horizontal errors depict the 0.75 mm estimated uncertainty of the source position. The observed instability of the peak energies is discussed in Appendix E

The peak energies at the positive- and negative-radii scanning positions have a discrepancy of up to 11%. The bulk event energy scale, as measured with dedicated calibration runs and environment background peaks, was stable to within 1 keV (see Appendix B for more details). This apparent asymmetry can be ascribed to the instability of the DCR parameter discussed at length in Appendix E. Subsequent measurements of the same scanning position (with the source unmoved between the measurements) demonstrate an average energy peak drift of per day, to increasing energies. The ordering of the measurements matches the deviation seen in Fig. 6: measurements began with small-magnitude negative radii and moved to ever-larger-magnitude negative radii, followed by measurements moving from large-magnitude positive radii to small-magnitude positive radii. Scans of positive and negative small-radius positions were separated in time by up to three months; these positions also show the largest observed energy deviation. Minimal apparent asymmetry is seen in scans of large-radius positions, which occurred within one month of one another.

Event pulse shape

Depending on the charge collection properties of the detector at their point of incidence, events can exhibit distinctive pulse shapes in a variety of ways. Two complementary discriminators, A/E and DCR, were studied in the TUBE scans.

The A/E values of pulses in PPC detectors depend strongly on the event incidence radius [17], rising as the radius falls to 0. In detectors with small passivated surface radii, events therefore have anomalously high A/E values, a signature used by the GERDA Collaboration to reduce the impact of backgrounds [19]. In these measurements, the DCR effect reliably identifies events at large radii, where A/E loses sensitivity as an event discriminator. We study the value of DCR as a function of position on the passivated surface to determine how the expected event rejection efficiency may vary depending on the geometric distribution of -emitting contaminants and to provide a detailed comparison to rising-edge-based event discriminators. These measurements also allow us to study the mechanism underlying the DCR effect, as discussed in Sect. 4.

In the DCR distribution for each data set with  mm, approximately 90% of events fall in a Gaussian peak. The remaining 10% of events have both degraded energies and DCR values, as seen in Fig. 7. Follow-up investigations in the CAGE (Collimated Alphas, Gammas, and Electrons) scanner, which allows for measurements with a wide range of beam angles of incidence, indicate that these events likely originate from shallower-depth energy depositions [21].2 We find, however, that their DCR value and prompt energy have the same ratio as those events falling in the energy peak, indicating that the mechanism causing their DCR effect is consistent with the mechanism of slow charge release for events in the peak. To simplify the modeling of events, we disregard the events falling outside of the energy peak when investigating the DCR response.

Fig. 7

The DCR and energy distribution measured with the source incident at  mm. The tail of events degraded in energy and DCR is thought to be due to ’s that deposit their energy at shallower depths, as discussed in Sect. 3.2. Events with low values of DCR are caused by pile-up

The DCR distribution for each scan is fit with a Gaussian, and the underlying background events are fit with a step function centered at the mean of the Gaussian. In the narrow fitting range used, the step function accounts for the high-DCR tail of the bulk event distribution. As in the analysis of the peak energy, an additional A/E cut is applied to select near point-contact events when studying small-radius scans. Since the events are better-separated from bulk events in DCR than in energy, this additional cut is only used for data sets with  mm when studying the DCR distribution. Events incident on the point contact itself do not have a distinct peak in DCR. In this case, the peak position is fit using an energy window in which the events dominate the spectrum; a 5 window around the peak energy is used. See Appendix D for details.

The DCR peak positions found in the fits are depicted in Fig. 8, along with their peak widths. At positions with  mm, the DCR values of events are consistently much higher than those of background events, and the DCR parameter is a highly efficient pulse shape discriminator, as shown in Sect. 3.3.

Fig. 8

The centroids of the DCR peaks in each data set; certain positions were studied in multiple data sets. Negative (positive) radius positions are indicated by blue downward-pointing (red upward-pointing) triangles. The vertical error bars give the -width of the DCR distribution peaks, and the horizontal errors depict the 0.75 mm estimated uncertainty of the source position. The dashed (dotted) line indicates the 99% (99.9%) bulk event acceptance DCR cut. The observed instability of the DCR peak position is discussed in Appendix E

The DCR values differ by up to a factor of 2 at the and scanning positions due to instability of the DCR parameter over the course of scanning. As discussed in Appendix E, several indicators of the passivated surface conditions point to positive charge build-up over time, driven by the incomplete event charge collection. This instability had a minor but measurable effect on the observed energy of the events, and a major effect on the observed slow charge recovery rate (i.e. DCR). The instability led to rising values of DCR over time, which would improve the surface rejection efficiency. The surface charge instability did not affect the bulk event charge collection efficiency or energy resolution, and is expected to be minimal in low rate environments like those found in searches.

To further study the charge collection properties at the passivated surface, we perform exponential fits to the tails of a set of pole-zero corrected waveforms from each position included in the model fits. We limit the data used to the first 3000 runs taken (1500 h, taken over the first 78 days of the first 217-day detector deployment), to reduce the effects of the surface charge instability. The average time-constant of delayed charge release, regardless of the position on the passivated surface, is found to be s. If the rate of delayed charge release remains constant over the full course of the event, this implies that time required to achieve full charge collection of an event is approximately  ms.

Since the magnitude of the DCR effect observed in this detector is largest for events with large incident radii on the passivated surface, it is highly complementary to the A/E-based rejection approach. The distributions of events incident at various radii in the A/E vs. DCR parameter space (see Fig. 9) clearly show the expected complementarity. At the smallest radii, the DCR effect is small, and the events are not well-separated from bulk events in the DCR parameter space. These events, however, have larger A/E than 99.9% of bulk events. At larger radii, the reverse is true; the events have large DCR values and A/E values similar to those of many more bulk events. This suggests that the combination of the two parameters should provide effective background rejection with the highest-possible bulk event efficiency.

Fig. 9

The distribution of A/E and DCR values for a range of scanning positions (indicated by the color scale) and a data set with no source shining on the detector surface (in black). All single-site non-muon events with energies between 1 and 6 MeV are included. The dashed blue line indicates the 99.9% acceptance value of A/E and the dotted red line indicates the 99.9% acceptance value of DCR

Survival probability

Though the expected efficiency of the DCR pulse-shape discriminator can be determined using readily-available calibration data (see Ref. [7]), the rate at which background events are misidentified as signal-like can only be determined from this dedicated event study. It is determined here as a function of radial position on the passivated surface using both the DCR-only rejection approach and a combined DCR and A/E-based rejection approach.

This quantitative analysis of the event rejection focuses only on the full-depth events in the Gaussian event energy peak; the survival probability for events in the tail of shallower-depth events is not calculated in this study. The fixed incidence angle of the source leads to a low rate for these events that varies dramatically with the radius of the scanning position. In preliminary simulations, the variation in the rate of this population based on changes to the scanning incidence angle is degenerate with variations due to conditions on the passivated surface (see Sect. 4). Therefore, any efficiency calculated for these events would be both high in uncertainty and of limited applicability to other detectors in other settings. The consistent ratio of the DCR to the prompt energy found in these events and the events in the peak leads us to expect that the rejection will perform similarly well for tail events in the ROI; the integrated energy of the delayed charge is 2-3% of the total prompt signal, corresponding to 40 to 60 keV at . This effect is well above the energy resolution limit of HPGe detectors in searches. This can also be seen from the band separation from the bulk events at the Q energy in Fig. 7. A quantitative study of the rejection efficiency for these events is planned as part of the CAGE measurements [21].

To calculate an survival probability, we calculate the excess counts in the 5 region centered at the mean value of the peak energy before and after applying a cut rejecting events. The expected background rate is determined from the spectral shape measured in a source-free data set. The estimation window determination is described in Appendix F. S and C are the signal regions in the data set taken with and without the source, respectively. The spectra before and after the application of the event rejection cut are indicated by the subscripts u and c, respectively. See Fig. 10 for an example. The rejection efficiency is given by:where is the ratio of the live time in the source run to the live time in the source-free run.

Fig. 10

The energy spectra measured with the TUBE detector without the source incident on the passivated surface, in blue, and with the source incident at r=18.0 mm, in red. Spectra are shown before (the solid lines) and after (the dashed lines) the application of the DCR cut. The hatched energy window is used for the survival probability calculation. The regions indicated here correspond to the variables in Eq. (4)

Event counts (, , etc.) are determined separately in each data set, using the signal energy windows and pulse shape parameters found for that data set. Then all Deployment 1 data sets taken at each particular source incidence position are combined to give a single efficiency evaluation for each position on the detector surface. The measurements in different regions are independent of one another, and given the low statistics, uncertainties on the count rates are assumed to be Poisson-distributed. 90% confidence intervals for the efficiencies are calculated via ROOT’s (see Ref. [22]) TRolke class with the profile likelihood method using the fully frequentist treatment described in Ref. [23]. Where the survival probability is found to be 0 or 1, an upper or lower limit is given.

Due to a slight (1–2) misalignment of the scanning and IR cup axes, the source beam is partially occluded at negative radius scanning positions. This leads to reduced event rates and increased statistical uncertainty at these positions; for some positions, these effects are mitigated by combining multiple data sets, leading to longer run times and higher statistics. For more discussion of factors affecting the uncertainty, see Appendix F.

The survival rates are evaluated with three different event cuts: the DCR cut is tested at bulk acceptance rates of 99 and 99.9%, and a cut using both DCR and A/E pulse shape parameters is also studied. In the latter case, we apply both a 99.9% bulk acceptance DCR cut and a 99.9% bulk acceptance A/E cut eliminating near-p contact events, for a total bulk acceptance of 99.8%.

The survival probability for each cut as a function of passivated surface position is shown in Fig. 11 with upper and lower limits indicated by the arrows. The uncertainties shown are purely statistical. The average survival probability for each cut is also given in Table 2: we consider the positions with  mm separately from those with  mm, since the former have a survival probability central value at or very near 0, and the latter at or very near 1 for cuts using only DCR. Combining these ranges obscures the behavior of the event cut. Any weighting strategy that attempts to account for the true incident event rate at different positions on the detector surface requires an assumption about the event source distribution and the effect of incidence angle. Since the incidence angle effect was not studied in these measurements and the source distribution depends on the experimental setup, all included positions are weighted equally in the average. We average the upper limits or lower limits of the survival probability, as is appropriate for each region. In one case ( mm with the 99% acceptance DCR cut), the 90% confidence interval on the survival probability is large enough to include the entire range from 0 to 100%.

Fig. 11

The survival probability of three possible event cuts, applied at each scanning position. The efficiency of the 99% (99.9%) bulk event acceptance DCR cut is given by the red (blue) points. The black points indicate the efficiency of an event rejection cut using both the DCR and A/E parameters, each applied at the 99.9% bulk event acceptance level. Downward-pointing arrows indicate that the uncertainty given is an upper limit

Table 2

The average survival probability of each cut, found based on the 90% confidence intervals. The average of the upper or lower limits is given, depending on the scanning position range in question. All included positions are weighted equally in the averages

α Event cut	Radius range	α Survival (%)
99.9% bulk acc. DCR	|r|≥3 mm	<0.8
99.9% bulk acc. DCR	|r|<3 mm	>8
99% bulk acc. DCR	|r|≥3 mm	<0.9
99% bulk acc. DCR	|r|<3 mm	>0
99.9% A/E and 99.9% DCR	|r|≥3 mm	<0.8
99.9% A/E and 99.9% DCR	|r|<3 mm	<14

The survival probability in every scan with  mm is consistent with 0%, within 90% confidence, and all scanning positions have a survival probability upper limit of less than 7.5%. Though the 99% bulk acceptance cut does remove more events from each data set than the looser 99.9% bulk acceptance cut, the upper limits found are dominated by the uncertainty due to fluctuations in the background level, leading to a slightly lower upper limit in the looser cut. At the positions closest to the p contact, only the combined DCR and A/E-based identification parameter is effective. Even in the relatively high-noise environment of the TUBE cryostat, this combination of pulse shape discriminators eliminates background events, with a maximum survival probability upper limit of 19%, across the entire detector surface with just 0.2% signal event loss.

Modeling the event response

Using simulations of signal formation in the detector, we can investigate the mechanism underlying the observed energy degradation and delayed charge effect. We first model the potential causes of energy loss to distinguish which charge carrier(s) are responsible for the observed energy degradation. Based on our findings, we then consider which of the affected carriers could be responsible for the observed delayed charge. The two sets of models considered (for energy and DCR) are therefore interrelated but not identical.

The detector may have a true, fully inactive layer on the passivated surface. Such a dead layer would affect electrons and electron-holes in exactly the same way, and the energy loss incurred would not vary by position on the passivated surface. The maximum allowable thickness of this dead layer can be calculated from the maximum observed energy of s on the passivated surface, which is 4.746 MeV. Using the Bethe formula as calculated in the NIST ASTAR database [11], we find that the maximum dead layer that is compatible with our observations is 3.5 m thick, far thicker than anticipated based on our understanding of the passivation process.

In TUBE, the observed energy of interactions varies dramatically depending on the incident radius of the . This indicates that while a thin dead layer may be present, charge loss must also be occurring. The radial dependence of energy indicates that positive and negative charge carrier contributions are affected differently. Therefore, we use a more sophisticated model of signal formation to investigate the detector response to surface events.

As described in Ref. [24], the Shockley–Ramo theorem can be used to calculate the induced signal expected at the point-contact of the detector. It can be written in terms of the weighting potential to give the integrated charge at the point contact:where q is the charge in motion in the detector and Q(t) is the charge induced at the point contact. The value of the weighting potential at a point can be conceptualized as the fraction of the total induced charge resulting from the minority carrier (in this case, electron) motion over the course of the full signal evolution. The value of the weighting potential varies as a function of the initial interaction position within the detector, and can be used to decompose the signal into the portions induced by the positive and negative charges. The weighting potential of the detector studied here is shown in Fig. 12 for , describing interaction positions at the passivated surface.

Fig. 12

The weighting potential at (i.e. events originating at the passivated surface) for the PONaMa-1 detector, calculated using mjd_fieldgen. Note that the electron contribution to the signal is small over most of the surface, and rises dramatically at small radii

Signals are simulated using the mjd_siggen software package [25], which has two parts. From the known impurity gradient and geometry of the detector, the electric field inside the detector is calculated using the stand-alone mjd_fieldgen software. The mjd_siggen software is then used to calculate the expected total energy as a function of radius for events occurring at the passivated surface of this detector. More details on the pulse shape simulation of this detector can be found in Ref. [26].

We limit the data used for fits to siggen-based models of charge collection on the passivated surface to just the first 3000 runs (1500 h, taken over 78 days) of the first 217-day detector deployment period. Using the data from this period allows us to reduce the effect of the surface charge instability while retaining measurements from across a full range of radial source positions. The energy and DCR peak positions found from fits of the 42 included data sets are averaged together by position, with fit uncertainties summed in quadrature. Datasets with  mm are excluded from the fits, since at these scanning positions, the source is partially or fully incident on the p contact.

The siggen charge collection behavior can be modified to reflect various charge loss mechanisms. Several models are considered:Models 2–4 are all potential behaviors exhibited by a detector in which positive charges are present on the passivated surface , with Model 2 being most likely at low surface charge densities, where self-repulsion of the charge cloud dominates the energy loss. Models 3 and 4 reflect the expected behavior at higher positive charge densities, where electrons are effectively pulled to the passivated surface. The instability in the charge drift time seen in these measurements (see Appendix E) suggests that positive charges may be accumulating on the surface over the course of the scanning measurement.

Energy Model 1: Charge collection for both electrons and electron-holes occurs as it does for bulk events. This model is provided for reference.

Energy Model 2: Some fraction of each of the electron and electron-hole signals does not contribute to the prompt energy. The fractions are held constant with respect to radius but allowed to float independently in the fit. This model matches the behavior expected if near-surface trapping and/or slow surface drift of both electrons and electron-holes were contributing to the prompt energy loss, with some measurable fraction of the electron component collected promptly.

Energy Model 3: The electron component does not contribute to the prompt energy, and a fraction of electron-holes (held constant with respect to radius, but allowed to float in the fit) is trapped or slowed instead of being collected promptly. This model matches the behavior expected if near-surface trapping and/or slow surface drift of both electrons and electron-holes were contributing to the prompt energy loss, with none of the electron component collected promptly.

Energy Model 4: The electron component does not contribute to the prompt energy, and electron-holes are collected as they are in bulk events. This model matches the behavior expected if electron surface drift or trapping were the only source of prompt energy loss.

Another category of energy loss models, corresponding to Models 2–4 but with the signs of the charges reversed, is also possible. In these models, electron-holes exhibit surface drift or near-surface trapping, and are partially or entirely missing from the prompt signal. These models, which describe possible behaviors of a detector with a net-negative charge on its passivated surface, lead to energies that fall with increasing incidence radius. Since this is opposite to the behavior we observe, we did not study these models in detail. Under different operating conditions, however, negative charge build-up could occur, and these models would be appropriate descriptions of the observed behavior, as discussed in Ref. [27].

We also study a variation of Models 1–4 that allows for the possibility of a fixed-thickness dead layer over the entire passivated surface. Such a dead layer would degrade incoming particle energies by a fixed amount, regardless of their position. In the case of Models 1 and 4, the dead layer thickness is set to its maximum allowed value of 3.5 m. In the case of Models 2 and 3, it is allowed to float. Table 3 summarizes which parameters are held fixed and which are allowed to float in each model’s fit. The variant of each model with an added dead layer thickness is indicated by the addition of a ‘d’ to the model number.

Table 3

A summary of the energy response models studied. and are the fraction of the electron and electron-hole charge collection signals that contribute to the prompt event energy, and d is the dead layer thickness

E model	Free param.	Fixed param.
1	None	f_e=1,f_h=1,d=0μm
1d	None	f_e=1,f_h=1,d=3.5μm
2	f_h,f_e	d=0μm
2d	f_h,f_e,d	None
3	f_h	f_e=0,d=0μm
3d	f_h,d	f_e=0
4	None	f_h=1,f_e=0,d=0μm
4d	None	f_h=1,f_e=0,d=3.5μm

Using only the centroid values from the passivated surface peaks, we minimize the residual sum of squares between the model and the peak data. Since the systematic uncertainty associated with the charge collection instability is not well-characterized and we cannot assume that this uncertainty is Gaussian-distributed, the absolute value of the reduced (i.e. whether its value is near 1) cannot be used to quantify the goodness-of-fit. We can, however, use it for model comparison, and to find the best-fit value for the electron-hole and electron contributions in Models 2 and 3.

Models 1, 3, 4, and their dead-layer variants are shown in Fig. 13a along with the average peak positions in the data and the energy ranges observed in the smallest-radius data sets. The best-fit is found in Model 3, when 19% of the electron-holes are lost to charge trapping or slow surface drift.

Fig. 13

The results of siggen simulations of energy loss (a) and delayed charge recovery () of events occurring on the passivated surface. The average values of spectral/DCR distribution fits to the data are given by blue circles, with blue boxes indicating the event energy range in cases where the peak position does not adequately capture the distribution shape

The addition of a dead layer does not improve the goodness-of-fit, but it does result in a non-zero best-fit thickness. In this case, the best fit is found when incident particles lose 245 keV upon entry, corresponding to a dead layer that is 1.2 m thick, and 15% of the electron-holes are lost to charge trapping. Though a thin dead layer may be present, these measurements are not able to distinguish its effects from those of incomplete electron-hole component collection. Future measurements of the passivated surface with lower-range incident particles, like low-energy and sources, should allow us to study these effects independently.

The preferred fit in Model 2 was found when the electron component of the signal was 0%. We can conclude that if electrons do contribute to the measured prompt energy, their effect is below the sensitivity threshold of these measurements. Determining this contribution would require measurements to be conducted using a setup with a smaller source spot size on the detector surface. This would allow a finer-grained investigation of ’s incident at small-radius positions, where the weighting potential is both high and rapidly varying.

Model 3 matches the observed energy most closely, but does not match the data exactly. A possible source of the deviation is the instability, which is described in more detail in Appendix E. As discussed there, the prompt energy of the passivated surface events increased over time. This corresponds to a decrease in the electron-hole loss fraction over time or an increase in the fraction of trapped electron-holes that are released into the bulk within the first 3 s of the signal. This is the behavior that would be observed if the passivated surface were becoming more positively charged over time. Variation in the trapping layer properties as a function of radius could also be responsible for some of the deviation. Both of these possibilities will be explored in future measurements.

Given this near-match, we conclude that all or nearly all of the negative charge is being trapped and/or slowed for interactions near the passivated surface, and that a fraction of the positive charge is also being trapped and/or slowed. These behaviors occur for events on the passivated surface regardless of the radial position of the interaction. Events incident on the p contact, on the other hand, do not show indications of significant charge trapping. The average energy loss observed at these positions is consistent with the loss seen in scans of the point contact of BEGe-type PPC detectors [15], and corresponds to a dead layer of 0.3 to 0.5 m.

Since we can conclude that both electrons and electron-holes are trapped and/or slowed in passivated surface events, the slow release of charges from either component could be responsible for the observed DCR effect. We consider the effects of both possibilities, and compare them to total energy collected as delayed charge at each scanning position. In this study, we disregard the effects of a thin fully-dead layer at the passivated surface, which would be minimal. The delayed energy is calculated as given in Eq. (2), which gives the amount of energy collected as delayed charge in the first 18 s of waveform digitization following the end of the prompt rise.

Three models of the DCR effect are considered:In Models B and C, the amount of delayed energy recovered in the first 18 s is a constant fraction of the total energy induced by the relevant charge carrier at each position. This is the case regardless of whether surface drift or near-surface trapping is responsible for the DCR: the surface charge drift velocity is constant for a particular drift direction relative to the crystal axis in the former case, and the charge trapping and subsequent re-release is assumed to affect a constant fraction of the charge carriers in the latter case.

DCR Model A: Both charge carrier species are allowed to contribute to the energy released as delayed charge. Each contributes a constant fraction of the component’s total prompt energy, as described below.

DCR Model B: The DCR effect is caused by the electron component of the signal, whether from electron drift on the passivated surface or trapped charge re-release.

DCR Model C: The DCR effect is caused by the electron-hole component of the signal, whether from electron-hole drift on the passivated surface or trapped charge re-release.

As in the case of the energy models, these correspond to the potential behaviors exhibited by a detector in which a net-positive charge is on the passivated surface. A different set of models, with the signs of the charges reversed, would be relevant if the behavior of the energies indicated that a net-negative charge was present on the passivated surface.

As seen in Fig. 13b, siggen simulations of Models B and C show opposite behavior as a function of the radial position on the passivated surface. Table 4 summarizes the free and fixed parameters in each model. The best-fit is found in Model C, when 10.2% of the missing electron-hole component (i.e. 1.9% of the total electron-hole component of the signal) is collected during waveform digitization. Adding the effect of a passivated surface dead layer increases these fractions slightly, but does not change the qualitative behavior. As in the energy degradation model, the electron component contribution to the delayed energy was driven to 0 in Model A. Slow electron transport or charge re-release may also be occurring, but it contributes negligibly to the signal shape in the 18 s window studied here.

Table 4

A summary of the delayed charge collection models studied. and are the fraction of the electron and electron-hole charge collection signals that contribute to the delayed event energy

DCR model	Free param.	Fixed param.
A	g_e,g_h	None
B	g_e	g_h=0
C	g_h	g_e=0

Model C correctly predicts that low DCR values will be observed at the smallest radii. Though the value of DCR at each position changes over time in these measurements (as discussed in Appendix E), DCR remains lower at small radii than at large radii. This indicates that though the electron-hole trapping and/or release rate may be unstable, under the conditions found in the TUBE scanner, the electron-holes are always responsible for the DCR effect.

Though siggen-based pulse shape simulations allow us to determine that both charge species are being lost or delayed, and that the electron-holes are responsible for the DCR effect, more detailed simulations that model the behavior of the full charge cloud are needed to determine whether surface charge drift effects are sufficient to explain the observed pulse shapes. These simulations are currently being developed, and will be discussed in a future publication.

Backgrounds in searches

The complementarity of rising-edge-based and delayed-charge-based background identification strategies motivates the differing approaches taken by the Majorana and GERDA collaborations. For detectors in which the -sensitive surface is small in radius, like those used in the GERDA experiment, a drift-time-based pulse shape parameter like A/E can be used as the sole strategy to reject events with high signal efficiency. The enriched detectors used in the Majorana Demonstrator, however, have much larger-radius -sensitive surfaces; the signal efficiency of such a cut would be unacceptably low. In this case, the DCR parameter is a far more effective alternative.

In the Gerda Experiment’s point-contact detectors, 234 events at high energy are identified in 5.8 kg year of exposure [19], indicating an rate of at least 0.110 counts/(kg day). Since Mirion’s BEGe geometry [28] (which has only a small annular passivated ditch region) is used for the point-contact detectors, a cut on high-A/E events can serve as the sole identifier of events [29]. This rising-edge-based cut incurs a fiducial volume loss of % in the Th double-escape (DEP) peak, which is used as a proxy for events [30].

Given the larger-radius passivated region of the Majorana Demonstrator enriched detectors, a similar approach would incur a fiducial volume loss of over 10%. The DCR method of event identification, on the other hand, has a signal sacrifice of between 0.8 and 3.1%, depending on the noise conditions of the data set [7]. Using the DCR discriminator, an rate of at least 0.276 counts/(kg day) is identified between 1 and 5.5 MeV. Unlike the high-A/E discriminator, this method is relatively insensitive to near-p interactions; if all the events between 2.7 and 5.5 MeV that are retained after the application of all analysis cuts [7] are assumed to be due to interactions, the unidentified ’s contribute an additional  counts/(kg day). Based on the results shown above (see Fig. 9), the implementation of an additional rising-edge-based estimator tuned to sacrifice 0.1% of bulk events would be expected to eliminate this remaining background contribution. Future analyses of the Majorana Demonstrator data will employ both DCR and high-A vs. E discriminators to identify events. The cut levels of both discriminators will be optimized based on in situ measurements.

The DCR cut is most effective at identifying highly-degraded passivated surface events, and is less effective for events where a large fraction of the incident energy is collected promptly. This trend can be seen in Fig. 14, where the survival fraction of events with energies above 2615 keV increases with increasing energy. Given this behavior, events that lie in the ROI and survive the DCR cut would have to be significantly degraded in energy before reaching the detector, without additionally exhibiting delayed charge energy loss. Such events can still be identified using a high A/E or A vs. E discriminator over much of the detector surface, and their rate will drive the optimal cut levels determined for the DCR and rising-edge-based discriminators.

Fig. 14

Energy spectra from the Majorana Demonstrator’s enriched detectors, including all data from Datasets 0 through 6 [7]. The red lines indicate the events kept by all other analyses (data cleaning, environmental, muon veto, multiplicity, and A vs. E multi-site cuts). The blue line shows the events remaining after all cuts, including the DCR cut. The inset shows the effect of the DCR cut in the background estimation window. The regions excluded due to backgrounds are shaded in green and the 10 keV window centered on is shaded in blue

The spectrum of identified events shown in Fig. 14 also demonstrates structure from 1 to 3 MeV, a feature that can be used along with charge transport and collection simulations of the passivated surface to constrain the spatial distribution of the emitters when building the full Majorana Demonstrator background model.

Background projections based on the Majorana assay program and simulations predict a flat background between 1950 and 2350 keV after rejecting possible peaks within that energy range. We exclude  keV ranges centered at 2103 keV (Tl single escape peak), 2118 and 2204 keV (Bi), and 2039 keV () from the background estimation window. In the remaining 360 keV window, the use of the DCR discriminator reduces the background rate from  counts/(keV kg year) to  counts/(keV kg year), over an order of magnitude reduction [7]. Suspected events do survive at energies between 2.7 and 5 MeV, indicating that some events may also survive in the region of interest. These events may be occurring on or near the contact, with the incoming being degraded in energy due to interaction in inactive materials (i.e. in the source of the background events). These events would be effectively suppressed by the implementation of a high-A vs. E cut, which will be used in future Majorana analyses.

The upcoming LEGEND experiment [31] will combine techniques from both the GERDA and Majorana Demonstrator experiments, including the use of both rising-edge and DCR-based rejection. LEGEND-200, the first stage of this project, will use existing PPC detectors of varying geometries, as well as new inverted-coaxial detectors [32]. Our experience with PPC detectors demonstrates that the effect of a large passivated surface has to be taken into account when rejecting events in analysis, and that the unique charge collection properties at this surface can help reduce the impact of such events in low-background experiments.

We are also exploring the possibility of detector designs that minimize or avoid passivation entirely. Further study of these and existing BEGe-type detectors will focus on the effectiveness of the DCR technique in these alternative geometries. LEGEND, like GERDA, will use a liquid argon shield, potentially leading to K background events from the decay of Ar. Future studies of charge collection near the passivated surface will determine whether these events exhibit a DCR component, allowing them to be identified.

Conclusion

This paper presented data taken with a p-type point contact detector produced by ORTEC that reproduces the design of the Majorana Demonstrator enriched detectors. The design of the TUBE cryostat and detector mount allows the passivated surface and p contact of the detector to be scanned using a collimated source. This enables a detailed study of the detector response to backgrounds from surface deposition of Rn progeny on parts and the detectors themselves. Due to degradation of the particle energy in inactive materials and poor charge collection in the detector itself, these events can be a significant background source in experiments like the search for in Ge. The work shown here, however, demonstrates that these events can be effectively rejected based on their pulse shape characteristics, while maintaining high signal event efficiency.

We find that ’s incident on the p contact do not exhibit the DCR effect and are observed at nearly their full energy. ’s incident on the passivated surface, on the other hand, exhibit significant charge trapping and can be highly degraded in energy, depending on their radius of incidence. Modeling the process of signal-formation in this detector, we find that the observed event energies are most consistent with a total loss of the electron-driven component of the signal, and a loss of between 15 and 19% of the electron-hole-driven component. The energy recovered as delayed charge in the first 18 s following the prompt signal corresponds to a recovery of 10.2% of that trapped electron-hole charge. The trapped electron-holes appear to be released with a 0.9 ms time constant. The sharply reduced energies and DCR effects measured at these near-p positions suggest that electrons do not play a significant role in the signal formation in the timescale of digitization. The formation mechanism of the DCR signal makes the DCR-based pulse shape discriminator highly effective at identifying events on the passivated surface of the detector, particularly for ’s incident at positions far from the p contact.

Subsequent measurements of this same detector in the event characterization test-stand described in Ref. [27] have found opposite behavior from that observed with the TUBE scanner, with the DCR component of events decreasing with increasing event incidence radius. In this case, the DCR parameter can still be used to identify events in the energy range relevant for searches. These measurements seem to indicate that the passivated surface becomes negatively charged, instead of the positive charge build-up we observe. This indicates that environmental conditions may impact the observed event response. In this case the vacuum conditions are believed to play a major role in the differing observed behavior. Additional measurements that study the effect of detector bias voltage, the position of the ground plane, the incidence angle, and the detector temperature are underway.

Detailed pulse-shape simulations that correctly model the effects of self-repulsion and diffusion in the dense charge cloud created by interactions have been developed in response to the behavior observed in this study. These simulations will be used to study the effect of varying surface charge densities and polarity on the pulse shape and energy observed in interactions. Along with additional detector surface characterization measurements that are being performed as part of the LEGEND experiment, these studies will inform the background model being developed for the Majorana Demonstrator.

Based on the presence of the DCR component and existing rising-edge-based discrimination techniques, we are able to reliably identify events on all of the detector surface with 100% efficiency in this test-stand measurement, while incurring 0.2% bulk event sacrifice. Both discriminators are expected to perform even more efficiently in the low-noise environment of the Majorana Demonstrator. Using the combined discriminator technique, we achieve a full order-of-magnitude improvement in the bulk event acceptance rate over the approach using only the rising-edge discriminator [19].

Applying the DCR technique in the Majorana Demonstrator, we find that the background rate in the ROI is reduced by an order of magnitude. Given the observed position dependence of the reconstructed energy, the spectrum of rejected events should allow us to identify the source of Rn backgrounds with appropriate simulations in an upcoming full background model of the Majorana Demonstrator. The Majorana Collaboration is also developing an improved version of the DCR discriminator that corrects for bulk-trapping effects, which is beyond the scope of this paper. Efforts have also begun to study whether deep learning techniques can be used to identify the pulse shapes associated with events. These new methods, along with the DCR technique, will be used to identify passivated surface events in other detector geometries, like the inverted-coaxial p-type point contact detectors planned for use in the upcoming LEGEND experiment.