1 Introduction

Dark matter (DM) is believed to be ubiquitous. It makes up 85% of the mass of the universe [1,2,3,4]. Many candidates including axions, low-mass DM, and weakly interacting massive particles (WIMPs) have been postulated [5,6,7]. WIMPs are expected to generate observable nuclear recoil energy through elastic scattering off nuclei [8]. Despite great efforts made in searching for axions and WIMPs [9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31], DM remains undetected. Recently, low-mass DM in the MeV-scale has become an exciting DM candidate [32,33,34,35]. To directly detect MeV-scale DM, a detector with sensitivity of measuring a single electron–hole (e–h) pair is required, since the energy deposition induced by MeV-scale DM through elastic scattering off electrons or nuclei is in the range of sub-eV to 100 eV [32, 36]. In 2018, Mei et al. proposed to detect MeV-scale DM utilizing germanium internal charge amplification (GeICA) [37] for the charge created by the ionization of impurities [36]. GeICA can potentially achieve a detection energy threshold of \(\sim \) 0.1 eV (100 meV), allowing a large portion of both electronic recoils and nuclear recoils in the range of sub-eV to 100 eV induced by MeV-scale DM to be accessible [36].

GeICA amplifies internal charge through impact ionization, which is a process first observed in Ge diodes a few decades ago [38, 39]. In this process, a charge carrier, electron or hole, with sufficient kinetic energy can knock a bound electron out of its valence state and elevate it to a state in the conduction band, creating an electron–hole pair. Carriers gain sufficient kinetic energy through applying a strong electric field. Impact ionization of Ge atoms requires higher electric field than that of impurities, since the bandgap of Ge is about 0.7 eV (700 meV) while the ionization potential of neutral impurity atoms in Ge is around 0.01 eV (10 meV). When a Ge detector is cooled down to below 10 K, the residual impurities in Ge start to freeze out from the conduction or valence band to localized states [40]. At below 6 K, the localized states become thermally stable and form electric dipole states [40], which are excited neutral impurity states with a binding energy of less than 10 meV. The dipole states can trap charge to form cluster dipole states [40] with even smaller binding energy depending on the operational temperature.

The formation of excited dipole states and cluster dipole states in p-type Ge is depicted in Fig. 1 [40]. The phase space of an immobile negative ion for trapping charge carriers is smaller than that of movable bound holes, whose motion is restricted by the Onsager radius, \(R=\frac{1}{4\pi \varepsilon \varepsilon _{0} K_{B}T}\), where \(\varepsilon \)= 16.2 is the relative permittivity for Ge, \(\varepsilon _{0}\) is the permittivity of free space, \(K_{B}\) is the Boltzmann constant, and T is temperature. Therefore, the probability of forming \(A^{-^{*}}\) states is higher than forming \(A^{+^{*}}\) states in a p-type detector. This is why electrons are trapped more severely than holes in a p-type detector.

Fig. 1
figure 1

Shown are the processes involved in the formation of the excited dipole states and the cluster dipole states in a p-type Ge detector operated at low temperatures, where \(\textbf{p}\) and \(\textbf{q}\) are the corresponding dipole moments

2 The experimental methods and the observed physical phenomenon

The impact ionization of impurities in Ge specimens has been reported for a range of temperatures (\(\sim \) 4.2 K to 298 K) by many authors from 1950s to 1970s [41,42,43,44]. The most recent impact ionization of impurities was reported by Phipps et al. with SuperCDMS-style detectors at 40 milliKelvin (mK) [45, 46], which is similar to this work with a detector made from a USD crystal [47, 48]. The detector was fabricated at Texas A & M University with four channels for charge readout, geometrically similar to SuperCDMS style detectors [49], as shown in Fig. 2.

Fig. 2
figure 2

Shown is the detector studied in this work with four electrodes read-out (\(Q_{1}, Q_{2}, Q_3, and Q_{4}\)). Four \(^{241}Am\) sources were arranged in such a way that one source was over each channel. The opposite side of the detector is grounded with uniform Al electrodes so that the electric field points along the z-axis of the crystal. The size of the detector is 10 cm in diameter and 3.3 cm in thickness with a mass of  \(\sim \) 1.4 kg

The detector was wire-bonded, mounted in a dilution refrigerator and tested at the K100 Detector Testing Facility at the University of Minnesota (UMN).

Fig. 3
figure 3

Shown is the orientation of the Am-241 source mover before it was installed

The detector was cooled down to \(\sim \) 40 mK and then tested in two separate refrigerator runs, Run 67 in 2018 and Run 74 in 2021. In Run 67 four \(^{241}\)Am sources were mounted above each channel on the detector (see Fig. 2). The lead collimators with 0.2 mm holes allow 59.54 keV \(\gamma \) rays through to the detector. Alpha particles from the sources were blocked by the source encapsulation. We observed 59.54 keV peaks in spectra from each channel and performed different measurements over a course of two weeks. In Run 74 a single \(^{241}\)Am source was mounted on a carriage that could be moved by a superconducting stepper motor (see Fig. 3). The source was of a different design with a 0.5 mm collimator hole and a lower rate of \(\gamma \) rays incident on the detector (about 75% as much as the sources used in Run 67).

Fig. 4
figure 4

Shown is the time-dependent charge response from Q4 at positive biases for Run 67. The baseline of 59.54 keV gamma ray from Q4 initially rose linearly for a few minutes, then quickly transitioned to an \(\sim \) exponential fall off over 10 s of minutes

Fig. 5
figure 5

Shown is the time-dependent charge response from Q4 at positive biases for Run 74. The baseline of 59.54 keV gamma ray from Q4 initially rose linearly for a few minutes, then quickly transitioned to an \(\sim \) exponential fall off over 10 s of minutes

3 Experimental data analysis and results

We used the data measured at various bias voltages to characterize the charge collection efficiency (CCE), defined as the ratio of the measured charge energy to 59.54 keV gamma rays, as a function of bias for each channel, and studied how the CCE varies with time under a given bias voltage. This work presents the results of an analysis of the 59.54 keV calibration line behavior, as shown in Figs. 4 and 5. The observation in particular had a rich set of time dependent behaviors, which are the amplitude of the signals from the 59.54 keV line to initially rise with time after bias, then fall to some steady-state value. This phenomenon was only present with positive biases. Since the \(^{241}\)Am sources are positioned at the biased side of the crystal and the mean free path of 59.54 keV gamma rays in Ge is about 0.09 cm, this means that this phenomenon is present for events in which electrons are collected almost immediately and holes are drifted through the full detector thickness. This indicates that the observed behavior has a polarity dependence for all events.

In Run 67 the observed phenomenon was much stronger in the center channel (Q4) and weakest in the outer ring channel (Q1), suggesting a radial dependence. In Run 74, with only one source present with a weaker event rate, the radial dependence is not evident, indicating an overall event rate dependence. In this work, we only study the 59.54 keV line in the center channel (Q4) at positive biases, since this channel outputs the clearest results.

The observed features of the time-dependent charge collection for the 59.54 keV line can be summarized as: (i) time-dependent impact ionization of drifting holes kicks off at a positive bias of 4.5 V (1.36 V/cm) and a mixing chamber (MC) temperature in a range from (30–35) mK during Run 67 and no such a phenomenon is observed with a negatively biased detector when drifting electrons; similarly, we observed the time-dependent impact ionization phenomenon for holes in Run 74 kicking off at a positive bias around 7 V (2.12 V/cm) and a MC temperature in a range from (35–40) mK (ii) the hole impact ionization initially increases linearly with a rate dependent on the bias and overall event rate; and (iii) after reaching a peak the increased charge signal falls exponentially with a fall time dependent on bias. The loss of both types of signal is easily explained as the loss of CCE due to the breakdown of bulk field from trapped charges. The linear increase of the hole impact ionization signal and lack of electron impact ionization indicate that the population of possible hole impact ionization sites starts small and is created by drifting electrons which are captured. This suggests a possible physical mechanism involving the combination of the following three processes: \(e^{-} + A^{0^{*}} \rightarrow A^{-^{*}}\), \(h^{+} + A^{-^{*}} \rightarrow e^{-} + 2\,h^{+} + A^{-}\), and \(h^{+} + A^{-} \rightarrow A^{0^{*}}\). The first process represents the creation of cluster dipole states (\(A^{-^{*}}\)) through drifting electrons across the detector induced by background radiation as shown in Fig. 1; the second stands for hole impact ionization of cluster dipole states; and the third is trapping of \(h^{+}\), which determines the loss of neutralization at the end of series.

There is also a possible process, \(E_{ph} + A^{-^{*}} \rightarrow e^{-} + A^{0{*}}\), which is the impact ionization of cluster dipole states through absorbing Neganov-Luke phonons [50], where \(E_{ph}\) is the energy of phonons created by drifting holes under a given electric field. This process is not a significant contribution in the experiment since the Neganov-Luke phonons have energy smaller than 1 meV [36]. However, the creation of cluster dipole charge states (\(e^{-} + A^{0^{*}} \rightarrow A^{-^{*}}\)) and the impact ionization of cluster dipole states (\(h^{+} + A^{-^{*}} \rightarrow e^{-}\) + 2\(h^{+} + A^{-}\)) inside the detector generate a dynamic electric field, which, in turn, impacts charge transport and charge creation. As a result, in this model the observed time-dependent impact ionization of 59.54 keV \(\gamma \) rays involves the growth of charge states as a function of time, the impact ionization of time-dependent cluster dipole states, and the loss of CCE due to the breakdown of bulk field from trapped charges, which is again a function of time.

Thus, the analysis presented here attempts to quantify the observed behavior in an empirical way. This will help lend insight into likely physical models that dominate the observed behavior. We assume that the detected charge energy, E(t), is related to the input 59.54 keV \(\gamma \) rays through the following equation:

$$\begin{aligned} E(t)= & {} E_{\gamma } \left\{ p_{0}+ p_{1}exp\left[ \frac{p_{2}}{p_{3}}(1-exp(-p_{3} t))\right] \right\} \nonumber \\{} & {} \times exp(-p_{4}t), \end{aligned}$$
(1)

where \(E_{\gamma }\) = 59.54 keV. Other terms in Eq. (1) are explained below:

(1) \(p_{0} + p_{1}\) represents the average CCE (\(\epsilon _{0}\)) at t = 0. If it is greater than 1, it means that charge energy is gained through impact ionization. It is expected that the reaction, \(h^{+} + A^{-^{*}} \rightarrow 2h^{+}+ e^{-} + A^{-}\), dominates the hole impact ionization at t=0 in the range of the applied field [51]. The parameter \(p_{0} + p_{1}\) is the direct measurement of CCE at t = 0 in Eq. (1), as depicted in Figs. 4 and 5.

(2) \(\{p_{0} + p_{1}exp[\frac{p_{2}}{p_{3}}(1-exp(-p_{3} t))]\}exp(-p_{4}t)\) is the average CCE (\(\epsilon _{t})\) for the charge created by the impact ionization at \(t>0\).

We can write \(\epsilon _{t}\) = \(\epsilon _{0}\times M(t)\times d(t)\), where M(t) stands for the charge gained through time-dependent impact ionization while d(t) is a time-dependent charge damping factor that describes the charge trapping due to a dynamic process described in points (3) and (4) which creates more charge states. Note that \(\epsilon _{0}\), d(t), and M(t) are correlated with applied electric field.

(3) \(p_{2}\) represents the rate of creating \(A^{-^{*}}\) states. It depends on the density of the dipole state, the overall event rate, the drift velocity, and the charge trapping cross section. Since both the drift velocity and the charge trapping cross section are field dependent, it is expected \(p_{2}\) has a correlation with applied electric field.

(4) \(p_{3}\) represents the rate of decreasing \(A^{-^{*}}\) states, which is proportional to the applied field and the overall event rate.

(5) \(exp(-p_{4}t)\) represents the loss of neutralization. If we let \(p_{4}=1/\tau \), the factor \(exp(-t/\tau )\) describes the collection of charge created by impact ionization falling exponentially as a function of time depending on bias. The parameter \(\tau \) measures the effective time constant for the loss of charge signal due to the loss of CCE resulting from the breakdown of bulk field from trapped charges.

Fig. 6
figure 6

Shown is the detected charge energy created by impact ionization as a function of applied bias voltage fitted by Eq. (1) using data from Run 67. The error bars represent the vertical width of the band shown in Fig. 4

Fig. 7
figure 7

Shown is the detected charge energy created by impact ionization as a function of applied bias voltage from Run 74 fitted by Eq. (1). The error bars represent the vertical width of the band shown in Fig. 5

Using this empirical model expressed in Eq. (1), we fit the data from Run 67 and Run 74. The observed trends and the fits are shown for different applied bias voltages in Figs. 6 and 7. Since the two runs were taken at two different periods at which the tower temperature of detector housing was slightly different. Thus, we expect that the detector responds slightly different to the impact ionization as it depends on the detector conditions including the thermal effect. For accurate reading in a temperature, the mixing chamber (MC) temperature was taken. The MC temperature in Run 67 was in a range between (30–35) mK for the data series shown in Fig. 6. Similarly in Run 74, the MC temperature was in a range between (35–40) mK for the data series shown in Fig. 7, being higher compared to Run 67.

At a given bias the observed phenomenon is stronger in Run 67 than in Run 74. This can be attributed to differences in the overall event rate in the detector between the two runs. By introducing an additional external source it was observed that the impact ionization effect seen in 59.54 keV events was stronger when the overall event rate in the detector was increased. In Run 74 the overall rate was lower because the \(^{241}\)Am source used produced events at a lower rate, only one source was used rather than four, and a lead shield was erected around the cryostat which reduced the background rate from radioactivity in the lab environment.

We summarize the fitting parameters in Table 1. The fitting to the impact ionization curve using Eq. (1) shows that all the parameters have a good correlation with the electric field. This indicates that the parameters have their own effect under the electric field. The variation of those parameters with the electric field is shown in Fig. 8. Since each fitting parameter is a combination of at least two physical processes, it is difficult to figure out the impact of the applied field on any single process.

Table 1 A summary of the fit parameters using Eq. (1). Note that the values of \(p_2\) and \(p_3\) are omitted for 5.4 V in Run 74 due to the rate of creating and decreasing \(A^{-*}\) states was very small, leading to abnormal values for \(p_2\) and \(p_3\), which do not follow the trend observed at higher bias voltages, where the dependency on the applied bias field is evident
Fig. 8
figure 8

Shown are the fitting parameters, \(p_0\) to \(p_4\) from Eq. (1) as a function of the applied field, E, and fits a linear regression model, which demonstrates the correlation of the parameters with the electric field using data from each run. The fitting functions for those parameters in Run 67 (solid lines fitted to the circles) are: \(p_{0}\)=0.25E+0.69, \(p_{1}\)=0.017E \(-\)0.0087, \(p_{2}\)=2.27E \(-\)2.57, \(p_{3}\)=0.49E \(-\)0.49, and \(p_{4}\)=0.036E \(-\)0.039. Similarly, the fitting functions for those parameters in Run 74 (dotted lines fitted to the squares ) are: \(p_{0}=0.47E-0.11\), \(p_{1}=0.0078E+0.011\), \(p_{2}=2.1E-3.5\), \(p_{3}=0.43E-0.61\), and \(p_{4}=0.022E-0.033\)

4 Discussion and perspective of the physics application

As can be seen from Fig. 8, although the parameters, \(p_{0}\), \(p_{1}\), \(p_{2}\), \(p_{3}\), and \(p_{4}\) obtained from Run 67 and Run 74 have similar tendencies as a function of applied bias field, the values of these parameters differ between the two runs from a minimum of a few percent to a maximum of \(\sim \) 40%, depending on the applied bias field. The cause of this difference is likely due to the overall event rates used in the two runs. Run 67 has a higher event rate compared to Run 74.

For bias 4.5 V and above in Run 67, the value of \(p_{0}+p_{1}\) is greater than 1.0, indicating impact ionization at \(t=0\). In Run 74 the bias required to produce impact ionization at \(t=0\) was greater, 9.4 V. The difference between the runs is apparently attributable to the higher event rate in Run 67 mentioned previously.

For stable operation at which the average CCE is not a function of time for a given electric field, the parameters \(p_{3}\) and \(p_{4}\) from the empirical model (Eq. 1) should be zero. Using the fitted functions obtained from Run 67, \(p_{3}\) = 0 and \(p_{4}\) = 0 at electric fields less than \(\sim \) 1 V/cm. At such a small electric field, the measured CCE, \(p_{0}\) + \(p_{1}\) = 0.9483 and is constant in time. This means that at this electric field, the detector is a normal detector, with no visible impact ionization. The charge collection efficiency is consistent with typical SuperCDMS detectors operated at a similar electric field [46] and suggests that detectors made from USD-grown crystals are suitable as SuperCDMS-style detectors in terms of charge collection.

At a higher electric field, for example at the applied bias of 9 V, the CCE is about 1.4 at \(t = 0\) from Run 67. This means a higher electric field is needed to generate impact ionization. If one assumes the charge collection time per event is less than 2 \(\upmu \)s in this large-size detector, we can predict the CCE as a function of applied electric field using the fitted parameters from Run 67 and the empirical model (Eq. 1). Figure 9 shows that the CCE can be increased by a factor of \(\sim \) 100 when the detector is operated at a field of \(\sim \) 400 V/cm with a charge collection time of 2 \(\upmu \)s. This means that the charge can be amplified by a factor of \(\sim \) 100 at a bias voltage of \(\sim \) 1300 V. Note that this extrapolation was obtained by extending the fitted empirical trend as a function of applied electric field by more than two orders of magnitude beyond the data. Therefore, the discussion above needs to be verified using experimental data in the future.

Fig. 9
figure 9

Shown is the predicted CCE as a function of applied electric field

5 Conclusion

We conclude that the observed time-dependent behavior of charge amplification for the 59.54 keV line is mainly due to the impact ionization of charged states by drifting holes across the detector. We attempted to set out possible mechanisms for the impact ionization, which is mainly due to the creation of cluster dipole states, and impact ionization of cluster dipole states. These two processes dominate the production of additional charge that contributes to the amplification of charge while drifting holes across the detector. Our results suggest that a detector with an appropriate density of charged states, created by controlled radiation from the appropriate density of neutral states, can be developed for searching for MeV-scale dark matter with a detection threshold as low as < 10 meV, which corresponds to recoil energy induced by sub-MeV dark matter particles. The challenge is to test such a detector at a higher bias voltage such as 1300 V without triggering avalanche breakdown. The detector avalanche breakdown occurs at about 12 V in Run 67 and about 15 V in Run 74. This could be increased by improving the fabrication of the electrical contacts. But the overall event rate, from the \(^{241}\)Am source plus background radiation, is also relevant to breakdown. Run 74 has a somewhat lower event rate compared to Run 67, and the avalanche breakdown voltage was somewhat higher. We expect the avalanche breakdown voltage can be much higher if events of only a few electron–hole pairs are generated by a well-controlled optical source with the detector being well-shielded from background radiation. More research and development is needed.