Translator Disclaimer
2 August 2016 Proton irradiation results for long-wave HgCdTe infrared detector arrays for Near-Earth Object Camera
Author Affiliations +
HgCdTe detector arrays with a cutoff wavelength of ∼10  μm intended for the Near-Earth Object Camera (NEOCam) space mission were subjected to proton-beam irradiation at the University of California Davis Crocker Nuclear Laboratory. Three arrays were tested—one with 800-μm substrate intact, one with 30-μm substrate, and one completely substrate-removed. The CdZnTe substrate, on which the HgCdTe detector is grown, has been shown to produce luminescence in shorter wave HgCdTe arrays that causes an elevated signal in nonhit pixels when subjected to proton irradiation. This testing was conducted to ascertain whether or not full substrate removal is necessary. At the dark level of the dewar, we detect no luminescence in nonhit pixels during proton testing for both the substrate-removed detector array and the array with 30-μm substrate. The detector array with full 800-μm substrate exhibited substantial photocurrent for a flux of 103  protons/cm2 s at a beam energy of 18.1 MeV (∼750  e−/s) and 34.4 MeV (∼65  e−/s). For the integrated space-like ambient proton flux level measured by the Spitzer Space Telescope, the luminescence would be well below the NEOCam dark current requirement of <200  e−/s, but the pattern of luminescence could be problematic, possibly complicating calibration.



We have been developing sensitive long-wave infrared (LWIR) 10-μm cutoff Hg1xCdxTe detector arrays for use in the proposed Near-Earth Object Camera (NEOCam) Discovery mission [Hg1xCdxTe (mercury cadmium telluride) is a II to VI ternary compound, whose molar cadmium fraction x can be varied to tune to the desired cutoff wavelength].1 NEOCam is a survey mission designed to find, track, and characterize asteroids and comets in our solar system, including most of those greater than 140 m in size that travel close to Earth, the near-Earth objects (NEOs). There is particular emphasis on finding those NEOs with the potential of impacting the Earth. NEOCam will operate at two wavelength ranges 4 to 5  μm and 6 to 10  μm. The shorter wavelength HgCdTe detector arrays have already been designed for the James Webb Space Telescope (JWST),2 and our group at the University of Rochester (UR), in collaboration with NASA Jet Propulsion Laboratory, has been working with Teledyne Imaging Sensors (TIS) to produce the detector arrays that will cover the longer wavelength range. Laboratory testing of several arrays grown on an 800-μm CdZnTe substrate, and hybridized to an H1RG multiplexer, demonstrated that the resulting detector arrays met all NEOCam requirements (see Table 1) for dark current, quantum efficiency, well depth, and noise.3 Short-wave infrared (SWIR) and mid-wave infrared (MWIR) HgCdTe detector arrays utilizing the same multiplexer, or the same family of multiplexers in a larger format, have been or will be employed in other space missions, including the Orbiting Carbon Observatory 2, the Wide-field Infrared Survey Explorer (WISE), the Hubble Space Telescope’s (HST) Wide Field Camera 3, Euclid, and JWST.2,4,57 Detector arrays flown in space must be robust against cosmic ray (CR) hits; therefore, we subjected the arrays from 12 to 63 MeV protons to determine the magnitude of the responses.

Table 1

Minimum NEOCam requirements and detector array characteristics.

NEOCam requirementGoalH1RG-17346H1RG-17354H1RG-16886
Array format1024×10242048×20481024×10241024×10241024×1024
Cutoff wavelength (μm)1010.39.99.910.2
Responsive quantum efficiency (RQE) (%)55 (Non-ARC)55 (Non-ARC)58 (Non-ARC)56 (Non-ARC)58 (Non-ARC)
CDS read noise (e)3630222222
Dark current (e/s)<200<1<1a<1a<1a
Well depth (e)>44,000>65,000>72,400>64,600>57,300
Operability (%)909093.793.596.7
Note: ARC, antireflection coating.

aThe dark current value quoted is the median and histograms are shown in Figs. 1–3.



Waczynski et al.8 showed that SWIR HgCdTe arrays intended for the Wide Field Camera 3 instrument exhibited an elevated signal level in background (nonhit) pixels during 15.7- to 63-MeV proton irradiation. They found evidence that the spatial distribution of the elevated signal level across the array is correlated with the responsivity to flood illumination with 800-nm light and the elevated signal level is proportional to some fraction of the proton energy deposited in the CdZnTe substrate. Energy deposited in the 800-μm-thick substrate can create electron-hole pairs, some of which may radiatively recombine and emit 775-nm radiation, corresponding to the bandgap of the CdZnTe.9

The observed luminescence has been eliminated on SWIR arrays by removing the CdZnTe substrate.8,10,11 MWIR arrays developed for the JWST were also substrate-removed by TIS, utilizing the same process used for SWIR arrays.12 In order to remove the substrate, the volume between the silicon multiplexer and the detector must be epoxy backfilled, because the light-sensitive HgCdTe layer is thin (10  μm). McMurtry et al.3 investigated the effect of epoxy backfill on the substrate-intact LWIR devices and showed that it did not adversely affect the dark current, well depth, quantum efficiency, or noise although an increase in interpixel capacitance was noted.

Fig. 1

Histogram of dark current for pixels with well depth >72,400  e for H1RG-17346, the substrate-removed array.


Fig. 2

Histogram of dark current for pixels with well depth >64,600  e for H1RG-17354, the array with 30-μm substrate.


Fig. 3

Histogram of dark current for pixels with well depth >57,300  e for H1RG-16886, the array with 800-μm substrate.


The first epoxy-backfilled and substrate-removed array produced for this project was also the first substrate-removed LWIR array produced by TIS. This array exhibited excellent dark current. However, 99.7% of pixels had well depth below the NEOCam requirement (second column, Table 1) and the quantum efficiency was extremely low. Following the production of this first substrate-removed array, we received another substrate-removed array, H1RG-17346. A quarter of the array met NEOCam requirements. Subsequently, TIS produced two more substrate-removed arrays that also did not meet NEOCam requirements.

The overall poor performance of the substrate-removed devices tested and the low yield of the LWIR substrate removal process motivated the question as to whether complete substrate removal was necessary. Consequently, we pursued the fabrication of a device with most, but not all, of its substrate removed. H1RG-17354 was delivered with 30-μm substrate remaining. H1RG-17354 exhibited excellent dark current, quantum efficiency, and noise characteristics, and met all NEOCam requirements. However, we were concerned that the remaining 30-μm substrate would produce some degree of “background” luminescence under CR irradiation. To mitigate this concern, we subjected the one successfully substrate-removed device (H1RG-17346), the partially substrate-removed device (H1RG-17354), and a fully substrate-intact device (H1RG-16886) to 12- to 63-MeV proton irradiation at the University of California Davis Crocker Nuclear Laboratory (CNL) cyclotron. Previous tests in our laboratory at UR found no evidence for any substrate luminescence from CR hits, mostly 4-GeV muons, for any of the 10-μm cutoff arrays satisfying the NEOcam requirements (Table 1), including the arrays discussed in this paper.13


Proton Stopping Power

High-energy protons lose energy as they travel through both the CdZnTe substrate and the HgCdTe detector layer. The energy loss can generate charge carriers via ionization. The energy loss for charged heavy particles traveling through a material is calculated to first order by the Bethe–Bloch formula14,15

Eq. (1)

where me is the mass of the electron, c is the speed of light, z is the charge on the particle (proton charge z=+1), Z is the effective atomic number of the material, A is the atomic mass, Tmax is the maximum energy transfer in a single collision, I is the mean excitation energy, and ρ is the density of the material. The constant K is given by

Eq. (2)

where NA is Avogadro’s number and re is the radius of an electron. The Bethe–Bloch formula describes the mean energy loss/unit length for 0.1βγ1000, where β is the ratio of velocity to speed of light in vacuum and γ is the Lorentz factor of the particle. The equation can also include a density correction, not included here, as it is not relevant.

Generally, relatively high-energy particles lose less energy to the medium and have a greater range, whereas relatively low-energy particles may lose enough energy to stop within the medium and thus typically cause more damage or upsets to electronic components. As high-energy protons traverse the detector, the liberated holes are collected by the diode as a current component. The effective ionization potential for long-wavelength cutoff HgCdTe is 1.04  eV/e, corresponding to an efficiency of converting energy into electron-hole pairs of about 10% (the band gap energy is about 0.1 eV).16 The total charge generated by high-energy protons is approximately17

Eq. (3)

chargepathlength*(dE/dx)/ionization energypere.
For this calculation, we assume normal incidence. Luminescence from the CdZnTe substrate is another consequence of ionization by lower energy protons (of the energy range tested). The most damaging energy range of protons is 0.1 to 20 MeV, since <20  MeV protons can stop in the substrate or detector material. The Bragg peak resides within the 30-μm substrate for proton energies <4.8  MeV and within the detector material for energies between 4.8 and 5.1 MeV. For arrays with full 800-μm substrate, the Bragg peak resides just within the substrate for 15.7 MeV protons.8 Although the CR spectrum is composed primarily of protons, alpha particles in the range of 0.1 to 20 MeV have approximately the same flux as protons in that range, and can be more damaging, because they are stopped much more efficiently than protons.18,19 The energy transfer dE/dx (Eq. 1), is directly proportional to z2, where z is the particle charge: therefore, an alpha particle of the same energy as a proton would be stopped in approximately one quarter of the distance.

For a typical spacecraft, some shielding will be present. This shielding tends to stop the relatively low-energy component of the expected CR spectrum, which consists of galactic CRs and solar quiet-time particles, while simultaneously degrading the energy of the remaining higher energy CRs.18 For 3 mm of aluminum shielding in all directions, all protons with energies below 24  MeV will be stopped [see Fig. 4, calculated using the stopping range of ions in matter (SRIM), a selection of software packages that can be downloaded from that calculate many features of the transport of ions in matter. A textbook by the same name describes the methodology and application].20,21 For the same amount of shielding, alpha particles with energy below 95.75 MeV will be stopped. The peak of the alpha particle component of the CR spectrum occurs at the same energy as that of the proton spectrum, but the flux of alpha particles not stopped by the shielding is 10 times lower than that of the protons.

Fig. 4

(a) Simulated proton tracks through 3 mm of aluminum. (b) Energy loss in aluminum as a function of target thickness for 24 MeV protons. Simulation performed using SRIM.


The cosmic protons with initial energies between 24 and 32 MeV will degrade to 0.1 to 18 MeV after passage through 3 mm of aluminum (from a TRIM calculation: TRIM, transport of ions in matter, is a program in the SRIM package that calculates the energy loss of ions moving through matter using a quantum mechanical treatment of atom–atom collisions). These lower energy, damaging CRs are experienced by a spacecraft during the solar storms emanating from the Sun after strong flares lead to coronal mass ejection events, but detectors are exposed to an approximately steady differential flux at energies <18  MeV (14% of the peak galactic CR differential flux at 300 MeV). Higher energy (>28  MeV) protons and CRs will tend to pass directly through the shielding as well as the detector and substrate material and will deposit less energy than lower energy protons in the detector material. The IRAC InSb and Si:As detector arrays on the Spitzer Space Telescope, with similar aluminum shielding, experienced transient rates of pixels hit ranging from 3 to 10  s1 over the course of the cryogenic mission, except during solar flares, which led to much higher rates.22


Experimental Tests and Methods

The primary objectives of the proton radiation tests are to

  • 1. characterize the spatial extent and energy deposited by isolated proton hits;

  • 2. examine any residual effects and recovery time from proton hits;

  • 3. investigate the dark current in nonhit pixels during proton irradiation; and

  • 4. assess array performance following a cumulative life-time dose.

The three detector arrays discussed in this paper are 1024×1024×18-μm-pixel-pitch LWIR HgCdTe infrared arrays on H1RG read-outs and were tested at CNL. The pixel design is very similar to that of the 2.5- and 5-μm HgCdTe TIS arrays, except that it is a 10-μm cutoff. Over several trips, H1RG-17346 and H1RG-17354 were irradiated with 63, 32, and 12 MeV protons at various fluxes. H1RG-16886 was irradiated with 34.4 and 18.1 MeV protons. Intermediate metal foils of various thicknesses were used to attenuate the proton beam to the desired energy. The beam energies quoted were measured at the exit of the cyclotron. Table 2 summarizes the test dates, the array tested, and the nature of the tests, as well as the beam energies utilized. Prior to the CNL proton tests, similar data were obtained in the Ames Laboratory and at CNL without the proton beam irradiating the array under test, in the same test dewar.

Table 2

Summary of data obtained over several radiation runs.

DateDetector arraySubstrate thickness (μm)Measurements
September 26, 2013H1RG-1734607.5 krad(Si) cumulative dose at 63 MeV
August 13, 2014H1RG-1735430Single event data at 63 MeV, 1 krad(Si) and 5 krad(Si) cumulative dose at 63 MeV
September 4, 2014H1RG-173460Single event data for energies 12, 32, and 63 MeV, 5 krad(Si) cumulative dose at 32 MeV
October 2, 2014H1RG-1735430Single event data at 12 and 32 MeV
September 22, 2015H1RG-16886800Single event data at 18.1 and 34.4 MeV

We made use of the NASA Ames array controller and acquisition system, as well as its test dewar, since these have been extensively used for other proton irradiation tests for space experiments, including JWST, Spitzer, WISE, and HST among others. The dewar and experimental setup at CNL has been described in detail by others.8 A 5-mil Kapton window at a side entrance of the dewar passed the beam through three radiation shield windows masked by 1-mil aluminum foil. An external aluminum aperture defined the proton-beam size to 37  mm×37  mm, exceeding the size of the detector array surface being irradiated (<25  mm×25  mm). Normal incidence was used for all tests, with the dewar window situated close to the cyclotron exit beam. The metal windows and masks in the beam path reduced the beam energies incident on the detector array during the first two runs by 0.4 MeV for 63-MeV initial beam energy, by 0.6 MeV for 32-MeV initial beam energy and by 1.3 MeV for 12-MeV initial beam energy (from SRIM simulations). The 18.1- and 34.4-MeV beam energies used in the tests on H1RG-16886 were degraded by 1 and 0.6 MeV, respectively. Also, the tested beam energies had a spread in energy due to the accelerator tuning process, and passage through both our windows and masks; e.g., the 12-MeV beam energy was both degraded and spread in energy such that the incident energies on our detector array under test were actually in the range of 10.57 to 10.93 MeV.

All data were obtained in sample-up-the-ramp (SUTR) mode at two reverse biases (150 and 250 mV). Preliminary dark tests at Ames prior to the trip to CNL, as well as data obtained during proton testing at CNL were obtained primarily in two modes: 4 repetitions with 16 samples per image and 2 repetitions with 64 samples per image. The initial proton tests were obtained at low flux: the minimum calibrated level at CNL is 105  protons/cm2s. CNL Laboratory scientists can reduce the level reliably by a factor of 10 to 104  protons/cm2s, although dosimetry is unavailable. Calibration at lower flux levels than 105  protons/cm2s must be done through analysis of hits on the array itself by reducing the flux until a small number of proton hits are shown in an image. A focal-plane temperature of 35 K was maintained throughout the experiments. We have shown that the detectors operate well between 30 and 42 K, and <40  K is the nominal temperature for the NEOCam 10-μm focal plane (See Fig. 5 in Ref. 3). Pixel voltages read out are converted to electrons by using our laboratory conversion factor, given by the typical nodal capacitance of 42 fF, normalized by the relative gains of the Ames and UR readout electronics.3


Data Analysis

Dark current and well depth maps for all arrays have previously been produced from measurements made in the laboratory at UR. Three sigma (3σ) clipping of the bias- and temperature-specific dark current and well depth maps was used to mask out pixels exhibiting either high dark current or low well depth, or both.

For these SUTR data, differences of consecutive frames were used to find pixels whose SUTR data exhibit a sudden jump in signal greater than 3× the standard deviation of the average signal from each pixel. These pixels are selected as possible proton hits.

For pixels flagged as a proton hit, a 5×5 box around the flagged pixel was masked off to examine the nonhit pixels surrounding the proton strike. We took SUTR data with 5.278 s between samples several times for each flux level and reverse bias level, while resetting the array in between SUTR data sets. Subsequently, the mean and standard deviation of the dark current, the slopes for a steady SUTR, were calculated for nonhit pixels common to all data. The mean “dark current” before irradiation, 0.3  e/s/pixel, is significantly above the upper limit to dark current, 0.04  e/s/pixel, obtained in UR Laboratory.13 We believe there is a slight light leak or glow in the Ames dewar; however, since the NEOCam dark current requirement is <200  e/s, this is an acceptable value.

To investigate the number of pixels affected by a proton hit, pixels previously flagged as a potential hit and a 5×5 region around that pixel were considered. For each transient event, we considered only pixels with a jump in signal 5σ above the median background to be part of the proton hit. After the jump, the pixel continues integrating up the ramp, with the same slope as before the jump.


Results and Discussion


Dark Current Before, During, and After Testing

The procedure outlined in Sec. 5 was used to reduce data obtained while the proton beam was tuned to a low flux level, before irradiation, and after. We provide those data for all three arrays below.



The first detector array we tested was H1RG-17354, an array with 30-μm substrate. Figure 5 shows the dark current and the cumulative dark current immediately before and during irradiation by 12 MeV protons.

Fig. 5

(a) Dark current for nonhit pixels (130,000  pixels), normalized to peak value, during low-dose 12-MeV proton irradiation for H1RG-17354, an array with a 30-μm substrate. The dark current data obtained at CNL with the beam stop in place before radiation testing are overplotted on the radiation data. (b) Cumulative dark current of nonhit pixels during 12-MeV proton irradiation. Cumulative dark current measured in the laboratory at UR, and cumulative dark current measured in the laboratory space outside the beam chamber immediately before irradiation, are overplotted. The applied reverse bias for these data is 250 mV.


Low-fluence data for nonhit pixels do not show an appreciable increase in dark current during irradiation compared with dark current data obtained in the laboratory space outside of the beam chamber before irradiation, within the uncertainty in the measurement. Table 3 summarizes the modal dark current for all beam energies at two different applied biases. The mode was computed using a function in interactive data language that determines the value of an array, where the maximum number of elements of the array are located, given a bin size. The uncertainty in the measurements is included: it is important to note that within the uncertainties in the dark measurement before and during irradiation, there is consistently no increase in modal dark current for all beam energies (the median yields the same result).

Table 3

Modal dark current in the laboratory space outside the beam chamber before irradiation and during irradiation for various beam energies for H1RG-17354, an array with a 30-μm substrate.

Beam energy (MeV)Modal dark current (e−/s/pixel)
150 mV reverse bias250 mV reverse bias
12a (and before irradiation)0.40±0.12 (0.23)0.36±0.09 (0.32)
32a (and before irradiation)0.41±0.12 (0.23)0.46±0.09 (0.32)
63b (and before irradiation)c0.08±0.04 (0.04d)

aData run on October 2, 2014.

bData run on August 13, 2014.

cNo baseline data are available for a 63-MeV beam at applied bias of 150 mV.

dBaseline data for a 63-MeV beam were obtained in the laboratory at NASA Ames in the same dewar.

The cumulative dark current for 250-mV applied bias is plotted in Fig. 5(b) to illustrate the dark current levels of nonhit pixels during low-dose irradiation. All nonhit pixels meet the NEOCam dark current requirement of <200  e/s/pixel during irradiation. In Fig. 5(a), a small dark current range of 0.5 to 2.5  e/s is shown, with a linear scale on the y-axis. A few pixels exhibit higher current up to 150  e/s and are not shown on the plot.

From Eq. (1), the energy loss in the 30-μm-thick substrate for 12 MeV protons is 394 keV, for 32 MeV protons is 176 keV, and for 63 MeV protons is 105 keV, after initially losing energy through the metal windows and masks. Since we observe no extra dark current for this array, we conclude that either the absorbed energy did not lead to substantial luminescence, or the energy was dissipated by other nonradiative processes through the crystal lattice.

In order to further illustrate the signal in the nonhit pixels, a correlated double sample image obtained with a long integration time is shown in Fig. 6. A correlated double sample image is defined as an image obtained immediately after reset subtracted from an image obtained after integrating over some time period. The integration time is, therefore, the time difference between those images.

Fig. 6

Correlated double sample image of H1RG-17354 (an array with a 30-μm substrate), during 12-MeV irradiation integrated for 327 s. The scale on the right side of the image is in electrons. The large multipixel blotches (clusters) result from individual proton hits. A small fraction of isolated hot pixels were identified as inoperable before radiation testing. In between these are the “nonhit” pixels showing dark current 118  e/327  s=0.36  e/s.


Fig. 7

(a) “Dark” current for nonhit pixels during low-dose 12-MeV proton irradiation for 1/4 of the area of H1RG-17346, the substrate-removed array. The “dark” current obtained at CNL with the beam stop in place before radiation testing is overplotted on the radiation data. (b) Cumulative “dark” current of nonhit pixels during 12-MeV proton irradiation. Cumulative dark current of the same pixels, tested in the laboratory at UR following exposure to a cumulative life-time dose of radiation is overplotted, as is the cumulative “dark” current in the laboratory space outside the beam chamber at CNL immediately before irradiation. The applied reverse bias for these data is 250 mV.


Fig. 8

Correlated double sample image of H1RG-17346, the substrate-removed array. This image was obtained by integrating for 74 s, whereas the array was being irradiated with 12 MeV protons. The applied bias for this integration was 250 mV.


Fig. 9

Photocurrent for nonhit pixels during low-dose 18.1-MeV proton irradiation for H1RG-16886. Data were obtained with 250-mV applied bias.


After a total dose of 7.5 krad(Si) (1.5× the life-time dose for NEOCam) was applied on August 13, 2014, the test dewar was returned to Ames, and kept at 35 K while we continued to take data. Although the test dewar was still “hot” (secondaries are detected from the dewar material), we searched for high dark current pixels that were consistent throughout the data frames. Excluding pixels that were inoperable before irradiation and transient high signals that could be mistaken for permanently damaged pixels, we found that for an applied bias of 250 mV, 1.05% newly identified pixels had high dark current and at 150-mV bias, 0.51% newly identified pixels had high dark current. This high dark current is consistent with bias dependent tunneling current due to dislocations induced by bombardment of the protons. The dark current operability of the remainder of pixels was still >90%, and the modal dark current is unchanged within the measurement uncertainty. As the modal dark current did not significantly change after the array was irradiated, we attribute the increase in baseline dark current from 0.04  e/s/pixel obtained on August 13, 2014, to 0.32  e/s/pixel obtained on October 2, 2014 (Table 3), to a small light leak in the Ames dewar. These baseline dark currents are well below the NEOCam dark current specification (Table 1). Although the plots in Fig. 5 are labeled as dark current, the current measured is the photocurrent from the light leak or glow plus the dark current. Hereafter, we refer to this level as “dark” current.



Data obtained with H1RG-17346 (substrate-removed) were reduced using similar methods as those employed for H1RG-17354. We present results for the best region on the array (one quarter of the array) in Fig. 7. The FWHM for a histogram of the dark current data is greater than for the data presented for H1RG-17354 because these data are SUTR-16, rather than SUTR-64, and therefore have a factor of four fewer SUTR samples. Additionally, we have only a baseline measurement before the array was irradiated for an applied bias of 250 mV. Therefore we compare “dark” current before and during irradiation of SUTR-16 data with an applied bias of 250 mV.

The uncertainty in the measurement of “dark” current during this radiation experiment was higher than for the measurements of the previous detector array, and the perceived elevation in the modal current in nonhit pixels in the fully substrate-removed array (Table 4) is within the uncertainties in the measurements. A 74 s correlated double sampled image of H1RG-17346 during 12 MeV proton irradiation is shown in Fig. 8.

Table 4

Summary of the modal “dark” current outside of the beam chamber before irradiation and during irradiation for various beam energies for H1RG-17346, the substrate-removed array.

Beam energy (MeV)Modal dark current (e−/s/pixel)
250 mV reverse bias
12 (and before irradiation)0.33±0.48 (0.14)
32 (and before irradiation)0.14±0.48 (0.14)
63 (and before irradiation)0.38±0.50 (0.14)

Following a cumulative life-time dose of radiation at CNL, H1RG-17346 was warmed up and sent back to UR, where all performance characteristics were retested. For pixels that were operable before particle irradiation there is no change in modal dark current within the measurement uncertainty after the cumulative life-time dose, although the operability decreased by a fraction of 1% for 150-mV applied bias and a few percent for 250-mV applied bias.



The final array tested, H1RG-16886 (full substrate), was irradiated with 18.1 and 34.4 MeV protons that differ from the energies used to test the two previously discussed arrays. The first energy, 18.1 MeV, was chosen to be as close as we were able to tune to 15.7 MeV, the energy at which the Bragg peak resides just within the 800-μm substrate. A beam energy of 34.4 MeV was chosen to be close to the 32-MeV beam energy utilized for the previous arrays tested.

The elevated current in nonhit pixels that we measure, shown in Fig. 9, is two to three orders of magnitude above the dark current plus assumed light leak measured before irradiation. To further illustrate the luminescence in nonhit pixels, a 5.3-s integration time correlated double sampled image obtained with a beam energy of 18.1 MeV, and 250-mV applied bias is shown in Fig. 10.

Fig. 10

Correlated double sample image of H1RG-16886, the array with 800-μm substrate, during 18.1-MeV irradiation integrated for 5.3 s with an applied bias of 250 mV. The scale on the right side of the image is in electrons. The flux level for these data is the same as in Table 5.


Table 5

Modal “dark” current in the laboratory space outside the beam chamber before irradiation and the background photocurrent detected during irradiation for various beam energies for H1RG-16886. The proton flux during irradiation by 18.1 MeV protons is 103±20  protons/cm2-s.

Beam energy (MeV)Modal dark current (e−/s/pixel)
150 mV reverse bias250 mV reverse bias
18.1 (and before irradiation)706±22 (0.3)749±10 (0.4)
34.4 (and before irradiation)81±2 (0.3)66±2 (0.4)

The results in Table 5 are for a proton-beam flux of 103  protons/cm2s. For the 5-μm cutoff wavelength InSb arrays on the Spitzer Space Telescope, which had pixels that were 30  μm in size and an array size of 256×256 (Area=0.59  cm2), the rate of pixels hit was 3 to 10  s1 (Ref. 22; Hora—private communication). Typically, they saw a single CR affecting between 2 and 4 pixels, smaller than what we report in Sec. 6.2, because the pixel size for Spitzer is larger. This translates from 1 to 5  CR/s for the Spitzer arrays. To directly compare with the 103  protons/cm2s from the CNL beam, the total integrated CR rate from Spitzer is 1.7 to 8.5  CR/cm2s, a factor of 12 lower. This reduces the 750  e/s luminescence we measure in the laboratory at CNL to at most 62  e/s in an ambient space-like environment, except when irradiated by solar flares. Since this luminescence current would be spatially and temporally variable, the 800-μm CdZnTe substrate might adversely affect the performance of NEOcam’s detector arrays.

Testing of H1RG-16886 in the UR Laboratory after radiation exposure and after an anneal to room temperature showed that the operability decreased by 1.2% for an applied bias of 150 mV and by 2.7% for an applied bias of 250 mV.


Single Proton Hits

We examined the clusters of pixels corresponding to a single proton hit, and found distinct differences in both the charge distribution and the number of pixels associated with a hit as a function of substrate width, and beam energy. We first examine single proton hit data for the array with 30-μm CdZnTe substrate, H1RG-17354.

These data are background subtracted (SUTR-64) consecutive differences. The total hit is a sum of all pixels in the difference frame. The distribution in charge was bimodal, with separation in bimodality larger at 12-MeV energy than at the higher beam energies. We used an exponential plus Gaussian fit to model the charge distribution data: the exponential is not anticipated to be physical, but likely represents a portion of a second Gaussian [see Fig. 11(a), for a representative plot]. These two distributions are fit to the data above 50,000  e by a least-squares method using the following model:

Eq. (4)

where the parameters A and c describe the shape of an exponential distribution, and the parameters B, xcenter, and w describe the amplitude, center, and width of a Gaussian distribution. A possible explanation as to the origin of the distributions follows.

Fig. 11

(a) Histogram of total charge collected per proton hit in H1RG-17354, an array with a 30-μm substrate (Binsize=104  e). (b) Histogram representing number of pixels affected by a single proton hit. The data shown here were obtained with an applied bias of 250 mV, where the array is irradiated with 12 MeV protons.


In Fig. 11(b), the total number of pixels affected by a proton hit is shown for the data in the Fig. 12(a). SRIM simulations show the beam traveling directly through both the substrate and the detector layer. The spread in beam energies is not expected to produce two distinct distributions, so it is possible that the smaller number of pixels per hit, arising from the “exponential distribution” of charge collected per proton hit, were in fact not caused by protons, but by secondary particles and radiation. When the protons pass through the material surrounding the detector array, a nuclear reaction can occur that releases secondary particles and radiation. An atom can become activated, decay, and release secondaries. These secondaries are also reduced in flux when less beam attenuation is used to produce the higher beam energies. The purple line drawn on the charge distribution in Fig. 11 indicates the charge calculated via the Bethe–Bloch formula assuming the 12-MeV proton-beam energy was attenuated to 10.7  MeV before entering the back surface of the detector. After passage through the 30-μm CdZnTe, in this example, the beam loses 394  keV and travels 12  μm through the HgCdTe to the front surface of the detector, where the junction and depletion region are located.14

Fig. 12

(a) Subarray image of a region on the array at the end of a 169-s integration and (b) the same region of the array immediately after reset.


The mean of the Gaussian charge distribution for all three beam energies is listed in Table 6 for H1RG-17354, corresponding to different numbers of pixels affected by a proton hit.

Table 6

Summary of mean charge collected per proton hit and number of pixels per proton hit cluster at various beam energies for H1RG-17354, an array with a 30-μm substrate.

Beam energy (MeV)Mean total charge per proton hit (e−)Mean number of pixels per clusterMean total charge per proton hit (e−)Mean number of pixels per cluster
150 mV reverse bias250 mV reverse bias

Table 7 shows the estimated charge generated in the HgCdTe detector layer after first losing energy through the metal windows and masks in front of the dewar and then through the 30-μm CdZnTe from Eqs. (1) and (3). These values can be compared to the mean total charge for the observed proton hits (column 4). The two values listed for the proton hits we observe correspond to the two applied biases utilized for testing. Similar results for H1RG-17346 (substrate-removed) and H1RG-16886 (full substrate) are shown in Tables 8Table 9Table 1011.

Table 7

Bethe–Bloch calculation for energy loss in 12-μm-thick HgCdTe after passage through 30-μm CdZnTe.

Beam energy (MeV)HgCdTe dE/dxTotal calculatedMean of total observed proton hits (e−)
(keV/μm)Charge (e−)150 mV bias, 250 mV bias
1216.7193,000217,000, 295,000
327.5988,00067,000, 98,000
634.5552,00064,000, 86,000

Table 8

Summary of mean charge collected per proton hit and number of pixels per proton hit cluster at various beam energies for H1RG-17346, the substrate-removed array. Larger charge events affect more pixels. These cluster sizes as a function of proton energy are similar to those seen for device H1RG-17354 (the array with a 30-μm substrate).

Beam energy (MeV)Mean total charge per proton hit (e−)Mean number of pixels per clusterMean total charge per proton hit (e−)Mean number of pixels per cluster
150 mV reverse bias250 mV reverse bias

Table 9

Charge deposition for H1RG-17346, the substrate-removed array.

Beam energy (MeV)HgCdTe dE/dxTotal calculatedMean of total observed proton hits (e−)
(keV/μm)Charge (e−)150 mV bias, 250 mV bias
1216.7193,000238,000, 238,000
327.6088,00084,000, 85,000
634.5553,00065,000, 65,000

Table 10

Summary of median charge collected per proton hit and number of pixels per proton hit cluster at various beam energies for H1RG-16886, the array with full 800-μm substrate. The smaller proton hit size corresponding to a larger total charge per proton hit is different from the previous observations with the 30-μm substrate and the array with no substrate.

Beam energy (MeV)Mean total charge per proton hit (e−)Mean number of pixels per clusterMean total charge per proton hit (e−)Mean number of pixels per cluster
150 mV reverse bias250 mV reverse bias

Table 11

Bethe–Bloch calculation for H1RG-16886–800-μm CdZnTe and 12-μm-thick HgCdTe detector configuration.

Beam energy (MeV)HgCdTe dE/dxTotal calculatedMean of total observed proton hits (e−)
(keV/μm)Charge (e−)150 mV bias, 250 mV bias
18.111.9138,000184,000, 245,000
34.47.1983,000104,000, 112,000

The proton-beam loses energy as it first passes through the window and metal masks in the beam path, then further as it ionizes material along its path, in both the CdZnTe substrate and the HgCdTe material. For charge spreading across multiple pixels from a single proton hit, the dominant mechanisms are coulomb repulsion and enhanced diffusion within the bulk material. Interpixel capacitance is responsible for only 3  pixels of charge spread.23 In the HgCdTe material, holes associated with a given proton hit from the electron-hole pairs produced at the incident surface will diffuse more than those produced close to a junction, and overall will result in a spread in the charge detected by the diodes. This is similar to the loss of image quality or lower MTF seen as the wavelength of light is changed from longer to shorter, i.e., the long-wavelength photons are typically absorbed near the pn junction after they pass through most of the bulk material, while the short-wavelength photons are absorbed near the backside surface or farthest away from the pn junction, hence allowing for the greatest amount of diffusion spread.

Tables 6, 8, and 10 show that for the lowest energy protons there is a variation in total hit size that is inversely related to the amount of substrate on a given device. The enhanced charge spreading seen in the low-energy proton hits can be attributed to the loss of the charge (hole) due either to radiative or nonradiative, e.g., phonon, processes if that hole is allowed to enter the substrate. A hole (charge carrier) may enter the substrate only if it has more than the bandgap energy of CdZnTe, which is about 1.6 eV.9 The hole mobility in CdZnTe is much lower (<100  cm2/V/s) than HgCdTe. Thus, the hole may be more easily lost (For example, by producing a 775 nm photon, or being dissipated via phonon thermalization) once it enters the CdZnTe. In other words, although all of the charge deposited by a proton will be initially located in a narrow column along that proton track, the charges (holes) will spread laterally due to coulomb repulsion and diffusion. All of those charges (holes) will be able to diffuse into 4π steradians unless there is a boundary such as the pn junction or the HgCdTe to vacuum boundary, in which case the charge (hole) may diffuse into only 2π steradians, i.e., reflect back in the direction of a pn junction. Hence, for the case of an array without a substrate, the holes closest to the HgCdTe to vacuum boundary will be reflected from that surface and produce, on average, a larger amount of lateral diffusion. For an array with a substrate, the holes may be lost if they enter the substrate.

However, for the higher energy proton irradiation, we find almost no change in pixel hit size from one device to the next. This variation in charge spreading versus energy of the incident protons is again related to the total amount of charge initially deposited and thus related to the amount of coulomb repulsion that would occur for the given amount of charge for the same volume.

After the array has been hit with protons while we integrate charge, the array immediately recovers after reset (Fig. 12). Note that the recovered array image shows a typical dark response with the occasional hot and dark pixels. Our results confirm those of Girard et al.,13 who showed that a pixel’s dark current, responsivity, and noise were unaffected 5.5 s after a muon hit.

No significant latent images from the proton hits were observed in the frames obtained immediately after frames exposed to individual hits.



Proton testing on NEOCam developed 10 µm cutoff arrays was conducted in order to assess whether we could detect luminescence, examine any residual effects and recovery time from proton hits, and assess the array performance following a cumulative lifetime dose of irradiation.

The “dark” current operability of nonhit pixels during each frame of low-fluence proton irradiation of H1RG-17354 and H1RG-17346 was 100%, fully meeting the NEOCam requirement of <200  e/s. Proton-induced luminescence is not significant after the removal of 95% of the CdZnTe substrate for our LWIR array, H1RG-17354. The observed “dark” current is unchanged within measurement uncertainties. On the other hand, H1RG-16886, the array with a fully intact CdZnTe substrate did lead to a substantial luminescence for a flux level of 103  protons/cm2s at a beam energy of 18.1 MeV.

Both H1RG-17354 and H1RG-17346 were subjected to a cumulative life-time dose of at least 5 krad(Si). Laboratory testing at UR posttotal life-time dose irradiation conducted on H1RG-17346 (after the array was warmed up and cooled down again) showed that the dark current and well depth operability decreased by a fraction of a percent with an applied bias of 150 mV, and with an applied bias of 250 mV, the operability had decreased a few percent (Table 1). For both applied biases, the detector array still meets NEOCam operability requirements after a life-time dose of radiation.

Measurements of the dark current of one array, H1RG-17354, immediately after a life-time dose of proton irradiation showed that at most 1% of pixels sustained permanent damage that led to substantial dark current. For our LWIR 10-μm cutoff HgCdTe detector arrays, we find that they will still meet the NEOCam dark current and operability requirements after a life-time dose of protons.

For pixels hit by individual protons, the short-term effect is transient. The charge is spread out over 10  pixels for higher energy protons and up to 21 pixels for lower energy protons hitting H1RG-17346 and H1RG-17354, and 10  pixels for low-energy protons hitting H1RG-16886. A typical hit is well below saturation, so data before and after the event are usable.


The research was carried out at UR, the NASA Ames Research Center, the University of California Davis Crocker Nuclear Laboratory, and the Jet Propulsion Laboratory, California Institute of Technology, under contracts with the National Aeronautics and Space Administration. Special thanks to Frank Wolfs for helpful feedback on this paper, and to the referees whose comments improved the paper.



A. Mainzer et al., “Survey simulations of a new near-earth asteroid detection system,” Astron. J., 149 (5), 172 (2015). Google Scholar


B. J. Rauscher et al., “New and better detectors for the JWST near-infrared spectrograph,” Publ. Astron. Soc. Pac., 126 (942), 739 (2014). PASPAU 0004-6280 Google Scholar


C. W. McMurtry et al., “Development of sensitive long-wave infrared detector arrays for passively cooled space missions,” Opt. Eng., 52 (9), 091804 (2013). Google Scholar


J. O. Day et al., “Preflight spectral calibration of the orbiting carbon observatory,” IEEE Trans. Geosci. Remote Sens., 49 (7), 2793 –2801 (2011). IGRSD2 0196-2892 Google Scholar


A. Mainzer et al., “Characterization of flight detector arrays for the wide-field infrared survey explorer,” Proc. SPIE, 7021 70210X (2008). PSISDG 0277-786X Google Scholar


S. M. Baggett et al., “The wide-field camera 3 detectors,” Proc. SPIE, 7021 70211Q (2008). PSISDG 0277-786X Google Scholar


R. Laureijs et al., “Euclid: ESA’s mission to map the geometry of the dark universe,” Proc. SPIE, 8442 84420T (2012). PSISDG 0277-786X Google Scholar


A. Waczynski et al., “Radiation induced luminescence of the CdZnTe substrate in HgCdTe detectors for WFC3,” Proc. SPIE, 5902 59020P (2005). PSISDG 0277-786X Google Scholar


M. A. Quijada and R. Henry, “Temperature evolution of exciton absorptions in Cd1xZnxTe materials,” Proc. SPIE, 6692 669206 (2007). PSISDG 0277-786X Google Scholar


E. C. Piquette et al., “Substrate-removed HgCdTe-based focal-plane arrays for short-wavelength infrared astronomy,” J. Electron. Mater., 37 (9), 1396 –1400 (2008). JECMA5 0361-5235 Google Scholar


R. Smith et al., “Noise and zero point drift in 1.7  μm cutoff detectors for SNAP,” Proc. SPIE, 6276 62760R (2006). PSISDG 0277-786X Google Scholar


J. D. Garnett et al., “2Kx2K molecular beam epitaxy HgCdTe detectors for the James Webb space telescope NIRCam instrument,” Proc. SPIE, 5499 35 (2004). PSISDG 0277-786X Google Scholar


J. J. Girard et al., “Cosmic ray response of megapixel LWIR arrays from TIS,” Proc. SPIE, 9154 91542A (2014). PSISDG 0277-786X Google Scholar


K. Nakamura and Particle Data Group, “Review of particle physics,” J. Phys. G: Nucl. Part. Phys., 37 (7A), 075021 (2010). JPGPED 0954-3899 Google Scholar


H. Bethe, “Zur theorie des durchgangs schneller korpuskularstrahlen durch materie,” Ann. Phys., 397 (3), 325 –400 (1930). ANPYA2 0003-3804 Google Scholar


P. W. Marshall et al., “Proton-induced transients and charge collection measurements in a LWIR HgCdTe focal plane array,” IEEE Trans. Nucl. Sci., 50 (6), 1968 –1973 (2003). IETNAE 0018-9499 Google Scholar


J. C. Pickel et al., “Radiation-induced transient effects in HgCdTe IR focal plane arrays,” Proc. SPIE, 5487 698 (2004). PSISDG 0277-786X Google Scholar


Y. I. Logachev, K. Kecskeméty and M. A. Zeldovich, “Energy spectra of low-flux protons in the inner heliosphere under quiet solar conditions,” Sol. Phys., 208 (1), 141 –166 (2002). SLPHAX 0038-0938 Google Scholar


J. A. Simpson, “Elemental and isotopic composition of the galactic cosmic rays,” Annu. Rev. Nucl. Part. Sci., 33 (1), 323 –382 (1983). ARPSDF 0163-8998 Google Scholar


J. F. Ziegler, “Software documentation and description for the stopping range of ions in matter,” (2016) March ). 2016). Google Scholar


J. F. Ziegler, J. P. Biersack and M. D. Ziegler, SRIM: The Stopping and Range of Ions in Matter, Lulu Press Co., Morrisville, North Carolina (2015). Google Scholar


J. L. Hora et al., “The effects of cosmic rays and solar flares on the IRAC detectors: the first two years of in-flight operation,” Proc. SPIE, 6276 62760J (2006). PSISDG 0277-786X Google Scholar


C. W. McMurtry et al., “James Webb space telescope: characterization of flight candidate NIR InSb arrays,” Proc. SPIE, 5167 144 (2004). PSISDG 0277-786X Google Scholar


Meghan L. Dorn received her BS degree in imaging science from Rochester Institute of Technology and her MS degree in optics from UR in 2012 and 2015, respectively. She is pursuing her PhD in materials science at UR.

Judith L. Pipher has worked on development of infrared detectors and detector arrays for ground-based and space-based astronomy since the mid-1960s. While a PhD student, She flew three far-infrared single detectors on astronomical rocket experiments. She later worked with William Forrest at Rochester on midwave IR arrays for ground-based astronomy, and then midwave IR arrays for the infrared array camera on the Spitzer Space Telescope. Long-wave IR HgCdTe array development began in the late 1990s.

Biographies for the other authors are not available.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Meghan L. Dorn, Judith L. Pipher, Craig W. McMurtry, Spencer Hartman, Amy Mainzer, Mark McKelvey, Robert McMurray, David Chevara, and Joshua Rosser "Proton irradiation results for long-wave HgCdTe infrared detector arrays for Near-Earth Object Camera," Journal of Astronomical Telescopes, Instruments, and Systems 2(3), 036002 (2 August 2016).
Published: 2 August 2016

Back to Top